Магнитотеллурическое зондирование измеряет естественные электромагнитные поля земли, и его главная слабость в том, что тот же источник поля, грозовая активность планеты, работает и как полезный сигнал, и как помеха. Атмосферники, импульсные излучения грозовых разрядов, приходят в точку измерения одновременно во всех компонентах поля и искажают передаточную функцию, из которой считается сопротивление породы. Подавление этих помех и есть основная работа обработки полевых данных, и от неё зависит достоверность результата зондирования. В статье разобрано, как работает метод, что делают с данными атмосферники и какие приёмы возвращают достоверность.
Магнитотеллурическое зондирование и его сигналы
Метод измеряет естественные вариации электрического и магнитного поля земли и извлекает из их отношения сопротивление породы. Передаточная функция считается как отношение электрической компоненты к магнитной, и сопротивление породы вычисляется по формуле ρa = |Z|²/(μ0ω), где Z импеданс среды, ω угловая частота. Для практического примера импеданс 1 Ом на частоте 10 Гц даёт сопротивление ρa = 1/(1,2566e-6 × 62,83) ≈ 1,27e4 Ом·м, типичное значение для сопротивлительных пород платформ. Рабочий диапазон метода охватывает частоты от 1e-4 до 1e3 Гц, и чем ниже частота, тем глубже проникает поле, и глубина зондирования растёт с периодом. Электрическая компонента измеряется заземлёнными линиями длиной в десятки и сотни метров, магнитная индукционными датчиками, и обе компоненты снимаются синхронно с точностью времени, которая для когерентной обработки критична. Естественный источник поля работает в узком окне по уровню, и сигналы в спокойные дни слабы, что делает метод зависимым от уровня естественной активности, и грозовая активность планеты работает одновременно и как источник, и как помеха. Уровень естественного поля меняется с широтой и временем суток, и практические окна зондирования строятся по активности источника. В спокойные дни сигналы слабы, и накопление требует часов, в возмущённые дни активность выше и накопление сокращается, но вместе с активностью растут и помехи от возмущений. Сезонная динамика добавляет свою картину, грозовой сезон приносит импульсные помехи, зимний спокойный источник при малой активности, и практические окна зондирования строятся по пересечению обеих динамик. Для планирования работ используют прогнозы геомагнитной активности, и запись планируется на окна с активностью источника, достаточной для короткого накопления.
Атмосферники как источник помех
Грозовая активность планеты создаёт около 40-50 разрядов в секунду, и каждый сильный разряд излучает импульс, который распространяется в волноводе земля-ионосфера на тысячи километров. Атмосферник приходит в точку измерения как импульсный выброс во всех компонентах, и его уровень на 10-100 дБ выше среднего уровня естественного поля в полосе измерения. Практическая тонкость метода в том, что атмосферник квазисинхронен во всех компонентах, он приходит одновременно в электрическую и магнитную компоненты, и его вклад в передаточную функцию не сокращается простым отношением, поскольку импульс искажает обе компоненты по-разному в зависимости от геометрии прихода. Одиночный сильный атмосферник смещает оценку передаточной функции на десятки процентов, и серия атмосфериков в записи искажает сопротивление по всему спектру. Отдельный класс помех, резонансы волновода земля-ионосфера на частотах 7,8, 14 и 20 Гц, которые создают устойчивые пики в спектре поля внутри рабочего диапазона метода, и они, в отличие от импульсных атмосфериков, квазистационарны. Суммарная картина помех объясняет, почему зондирование в грозовой сезон требует иных приёмов обработки, чем в спокойный, и почему одна и та же точка в разные условия даёт разные сопротивления при плохой обработке. Структура атмосферника заслуживает отдельного взгляда. Импульс разряда длительностью в миллисекунды приходит с серией повторных приходов, отражённых от ионосферы, и серия эхо растягивается на секунды. В частотной области серия эхо создаёт интерференционную структуру с провалами и пиками, и её масштаб в спектре поля виден как быстрое колебание уровня между соседними частотами. Для обработки это означает, что атмосферник искажает не один отсчёт спектра, а интерференционную полосу частот, и его исключение работает по окнам во времени, где серия эхо живёт. Расстояние до разряда работает на структуру, ближние разряды дают резкие выбросы, дальние серии эхо, и обработка закрывает оба класса окнами времени.
Дистанционная опорная станция в шумных условиях
Главный приём подавления локальных помех, вторая станция в десятках километров от рабочей. Дистанционная опорная станция измеряет те же компоненты поля, и её данные используются в кросс-спектрах, где передаточная функция считается не от собственной магнитной компоненты рабочей станции, а от магнитной компоненты опорной. Логика приёма простая, атмосферник, пришедший в опорную станцию, приходит и в рабочую с задержкой, и его вклад в кросс-спектры сокращается, поскольку помеха коррелирована между станциями слабее, чем полезный сигнал. Локальная помеха, работающая только на рабочей станции, сокращается ещё сильнее, она некоррелирована с опорной вовсе. Числовой пример показывает масштаб приёма, оценка передаточной функции с погрешностью 20-30 процентов по собственной магнитной компоненте в шумных условиях падает до менее чем 5 процентов с дистанционной опорой, и разница между этими оценками определяет достоверность результата. Цена приёма тоже известна, вторая станция с датчиками и синхронизацией удваивает полевую работу, и расстановка опорной станции требует её выноса за пределы локальных помех, что в практических работах означает дистанции 5-10 км и больше. Практика зондирования сделала дистанционную опору стандартом в шумных условиях, и без неё результаты в грозовой сезон не принимаются. Расстановка опорной станции заслуживает отдельного взгляда. Вынос опорной за пределы локальных помех противоречит её работе на тот же источник поля, и практический компромисс, дистанция, при которой помехи некоррелированы, а сигнал ещё коррелирован. Источник поля крупномасштабен, и его корреляция на дистанциях в десятки километров сохраняется, локальные помехи коррелируют на сотнях метров, и окно дистанций 5-10 км закрывает обе задачи. В практических работах используют и две опорные станции, и их кросс-спектры дают дополнительное подавление, и каждая опорная добавляет полевую работу, но и качество вместе с ней.
Робастные методы оценивания передаточной функции
Вторая линия обороны работает на статистике оценивания. Классическая оценка передаточной функции по методу наименьших квадратов чувствительна к выбросам, и один атмосферник с уровнем в 100 раз выше сигнала тянет оценку за собой пропорционально квадрату своего уровня. Робастные методы, M-оценки с весовой функцией, снижают вес выбросов, и оценка становится нечувствительной к отдельным атмосферикам. Весовая функция назначает каждому измерению вес по его отклонению от текущей модели, выбросы получают малый вес, и оценка сходится к значению, определяемому основной массой данных. Числовой пример показывает масштаб, выброс в 100 раз выше сигнала при обычном усреднении тянет оценку пропорционально своей величине, при робастном весе 0,01 его вклад падает до уровня обычного измерения, и серия выбросов перестаёт искажать результат. Дополнительный приём, медианные оценки вместо средних, медиана устойчива к выбросам аналитически, и для серий кросс-спектров медианное накопление работает там, где среднее проваливается. Скорость сходимости робастных методов ниже классических, и практические реализации работают в два прохода, сначала классическая оценка для определения выбросов, затем робастная с весами для финального значения. Выбор весовой функции заслуживает отдельного взгляда. Мягкие веса дают плавное понижение вклада выбросов, но не убирают их полностью, жёсткие веса отсекают выбросы полностью, но рискуют отсечь и сильный полезный сигнал. Практические реализации работают с жёсткими весами по порогу, вычисленному из дисперсии данных, и порог выбирается так, чтобы его превышало менее одного процента измерений. Проверка выбора идёт по повторным оценкам с разными порогами, устойчивая оценка при изменении порога в разы считается достоверной, плавающая требует пересмотра.
Накопление и режекция выбросов
Третья линия, работа с самой записью. Накопление по длинным записям повышает отношение сигнала к шуму пропорционально корню из времени, и запись в часы даёт выигрыш в десятки раз против минутных записей. Прямое накопление атмосфериков, которые приходят в обе станции, работает не всегда, и перед накоплением запись чистится от выбросов, где окна с обнаруженными атмосфериками исключаются из обработки. Обнаружение выбросов работает по порогу на уровне поля, и окна с превышением порога в 5-10 раз над средним исключаются. Числовой пример показывает масштаб, атмосферник длительностью в миллисекунды на частоте опроса в сотни герц занимает единичные отсчёты, и исключение их окон теряет доли процента данных при сокращении искажений на десятки процентов. Режекция на гармониках сети работает отдельным контуром, гармоники 50, 150 и 250 Гц режутся узкополосными фильтрами, и промышленный фон сокращается на десятки децибел. Для резонансов волновода на 7,8 и 14 Гц работают не фильтры, а накопление по длинным записям, где стационарные резонансы дают устойчивые пики в спектре, и их вклад в передаточную функцию оценивается и исключается аналитически. Сумма приёмов, режекция выбросов, накопление и робастная оценка, и определяет достоверность финального результата.
Расчёт кажущегося сопротивления и его погрешности
Финальный расчёт сопротивления ведётся из передаточной функции по формуле ρa = |Z|²/(μ0ω), и погрешность входной функции удваивается в сопротивлении. Погрешность передаточной функции считается из когерентности сигналов, и при когерентности 0,9 погрешность импеданса составляет единицы процентов, а сопротивления вдвое больше. Для практического примера импеданс 1 Ом на 10 Гц с погрешностью 5 процентов даёт сопротивление 1,27e4 Ом·м с погрешностью 10 процентов, и фазовая погрешность добавляет искажение в инверсию профиля. Фазовая компонента импеданса работает на глубинность профиля не меньше, чем модуль, и её погрешность искажает границы слоёв при инверсии. Практические критерии качества данных включают когерентность выше 0,8, погрешность сопротивления ниже 10 процентов и согласованность повторных измерений, и данные, не отвечающие критериям, отбрасываются. Кривая сопротивления по частотам проверяется на гладкость, резкие скачки между соседними частотами говорят о помехах, и участки с скачками либо пересчитываются, либо исключаются. Такая дисциплина качества превращает кривую зондирования в достоверный профиль, и её нарушение превращает красивую кривую в артефакт.
Порядок обработки полевых данных
Обработка полевых данных укладывается в пять шагов:
- синхронизация записей рабочей и опорной станций по точному времени;
- обнаружение и исключение окон с атмосферниками и гармониками сети;
- расчёт кросс-спектров с дистанционной опорой и робастное оценивание передаточной функции;
- расчёт сопротивления и фазы по частотам с контролем когерентности и погрешностей;
- сглаживание кривой, проверка повторными измерениями и передача в инверсию.
Такой порядок превращает зашумлённую запись в достоверный профиль, и его дисциплина определяет, что увидит инверсия. Практика зондирования подтверждает масштаб приёмов, станции с дистанционной опорой и робастной обработкой дают повторяющиеся кривые в условиях, где классическая обработка расходится на десятки процентов, и достоверность результата здесь не роскошь, а условие работы метода. Подавление атмосфериков остаётся той частью зондирования, где физика помех известна, и работа с ней, от опорной станции до робастного оценивания, определяет, достоверна ли кривая сопротивления, переданная в инверсию, или это произведение грозовой активности, а не породы.