Вестник МГТУ, 2026, Т. 29, № 2.

Обработка данных в режиме МТЗ была выполнена с помощью программы VMTU-1d для обработки и одномерной интерпретации данных станции VMTU-10. В ней реализована стандартная процедура обработки магнитотеллурических данных, предложенная Семеновым в 1985 г. (Семенов, 1985), которая актуальна и сегодня. Для каждой измеренной компоненты поля рассчитываются усредненные спектры и корреляционные соотношения между сопряженными и антисопряженными компонентами электромагнитного поля. Далее по частотам, удовлетворяющим критериям отбора, определяется отклик среды Zxx, Zxy, Zyx, Zyy, по которому определяется кажущееся удельное сопротивление: РXV f = 2 ^ ■ IZxyf 2 , Фxy ( f ) = argiZxy ( f ) ) , (1) РУ* (Л> = 2 ^ ^ IZyxf 2 , Фyx (f ) = argiZyx ( f )), (2) где д0 - магнитная проницаемость вакуума, f - частота. С помощью вращения матрицы тензора импеданса по формуле Z '(а) = R(a) ■ ( Z Z Zxx Z x Z Z V yx yy У TX 4 - 1 ^ \ I C0S(a) Sin(a) I R(a) ,R(a) = \ (3) - sin(a) cos(a) У определены направления главных компонент поля относительно магнитного меридиана, значения кажущегося удельного сопротивления pxy, pyx (рис. 2, а) по этим направлениям (ТЕ и ТМ), фазы импеданса ф^, 9 yx (рис. 2, б) по формулам (1, 2), а также строятся полярные диаграммы (рис. 2) главных и дополнительных компонент тензора импеданса (Бердичевский, 2009). Дополнительно рассчитаны параметры асимметрии Свифта skews и Бара skewb, определяющие размерность нижнего полупространства: | z „ + Z l skew =]------------ [, (4) s Z - Z xy yx\ skew b= |lm(Zxy■z ; +Zxx■O j (5) Z - Z xy yx\ где Zxy, Zyx- главные, Zxx, Zyy - дополнительные компоненты тензора импеданса. Для дальнейшего рассмотрения выбирались данные, соответствующие двумерному приближению (Бердичевский, 2009). Полярные диаграммы главной компоненты в случае одномерного и двумерного приближения имеют форму окружности, эллипса, или "расширенной" восьмерки, а диаграммы дополнительной компоненты - форму четырехлистника значительно меньшей амплитуды. В случае трехмерного, ассиметричного приближения диаграммы обеих компонент имеют сложные, близкие по амплитуде формы, а параметры Свифта и Бара в одномерном и двумерном приближении стремятся к нулю. На рис. 2 показаны вычисленные по формуле (3) полярные диаграммы главной компоненты (ху), обозначенные синим цветом, и дополнительной (хх), обозначенные красным. Как видно из рисунка, имеется согласованность между формой полярных диаграмм и значением инвариантов skew. Таким образом, произведен отбор частот, соответствующий двумерному приближению, при этом часть точек была исключена полностью. В результате для трех профилей было отобрано 28 точек МТЗ. В каждой точке определены направления главных компонент поля относительно магнитного меридиана, значения кажущегося удельного сопротивления pxy, pyx (рис. 2, а), фазы импеданса ^ xy, qyx (рис. 3) по формулам (1, 2), отобранные для двумерного приближения. На рис. 3, 4, 5 приведены кривые кажущегося удельного сопротивления и фазы импеданса, соответственно для профиля 1, 2 и 3. В ряде точек видны узкие частотные диапазоны и довольно большой разброс по погрешности, это связано как с присутствием высоких промышленных помех, так и с наличием отклика сигнала от трехмерных неоднородностей, которые были исключены с помощью параметров Свифта и Бара, а также полярных диаграмм. За исключением нескольких точек наблюдений, общей чертой для амплитудных кривых является нисходящий вид на периодах больше 1 с, максимумы амплитудных кривых находятся на уровне 104 Ом-м и выше. Фазовые кривые в пределах погрешности, в основном, практически одинаковы для обеих поляризаций, что подтверждает корректность выбора критериев отбора данных.