Сглаживающие кубические сплайны

МГУ им. М.В. Ломоносова, Геологический Факультет

Кафедра геофизических методов исследования земной коры

Практикум по электроразведке (5 курс)

Построение и S - интерпретация кривых МТЗ

Практическое задание

В рамках настоящей работы нужно выполнить сокращенный вариант обработки и интерпретации данных МТЗ. Исходными являются записи МТЗ, относящиеся к четырем точкам регионального профиля (шаг по профилю 20 км) и априорная геолого-геофизическая информация. Результатом выполнения задачи должен стать геоэлектрический разрез по этому профилю. При выполнении данной работы необходимо сделать следующее :

1. Для первой точки МТЗ с помощью программы СORRECT выполнить второй этап обработки данных МТЗ - построить сглаженные частотные зависимости компонент тензора импеданса (4 амплитудных и 4 фазовых кривых, а также эффективные кривые) (см. раздел 4 данной разработки).

2. С помощью программы MTS-PROFILE выполнить анализ полученных кривых Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru и определить, какие из них наименее искажены влиянием приповерхностных неоднородностей (см. раздел 5).

3. По программе MTS-SI провести контролируемую трансформацию неискаженных кривых с целью получения кривых Сглаживающие кубические сплайны - student2.ru , отражающих нарастание интегральной проводимости с глубиной (см. раздел 6).

4. Для каждой точки зондирования по полученной кривой Сглаживающие кубические сплайны - student2.ru с учетом априорной информации определить параметры геоэлектрического разреза (сопротивления и мощности слоев) (см. раздел 7). По результатам интерпретации кривых построить геоэлектрический разрез по профилю.

Программы для работы с данными находятся в каталоге C:\STUDENTS, сами файлы данных и файл с априорной информацией read_me - в каталоге C:\STUDENTS\N, где N - номер вашего варианта.

Тензор импеданса

Возникновение метода магнитотеллурического зондирования (МТЗ) связано с именем А.Н. Тихонова. В 1950 году им было показано, что отношение ортогональных горизонтальных электрических и магнитных компонент магнитотеллурического поля, получившее название импеданса Сглаживающие кубические сплайны - student2.ru , зависит от частоты и свойств среды, но не зависит от источника.

Первоначально использовалась горизонтально-слоистая модель Земли (модель Тихонова - Каньяра), в которой импеданс Сглаживающие кубические сплайны - student2.ru может быть определен на основе полевых данных как по формуле Сглаживающие кубические сплайны - student2.ru , так и по формуле Сглаживающие кубические сплайны - student2.ru , где Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru и Сглаживающие кубические сплайны - student2.ru - комплексные амплитуды соответствующих компонент электромагнитного поля. Однако практика очень скоро показала, что в районах с горизонтально-неоднородным строением среды импеданс Сглаживающие кубические сплайны - student2.ru зависит от ориентации установки. М.Н. Бердичевским и Т. Кантуэллом было предложено перейти от скалярного импеданса к тензору импеданса

Сглаживающие кубические сплайны - student2.ru

Тогда горизонтальные компоненты электрического и магнитного полей связаны соотношениями :

Сглаживающие кубические сплайны - student2.ru (1.1)

Компоненты тензора Сглаживающие кубические сплайны - student2.ru и Сглаживающие кубические сплайны - student2.ru называются главными, и именно по ним строятся кривые МТЗ. Компоненты Сглаживающие кубические сплайны - student2.ru и Сглаживающие кубические сплайны - student2.ru называются дополнительными, для горизонтально-слоистой среды они равны нулю. Кроме того, при интерпретации данных МТЗ часто используется такая величина, как эффективный импеданс

Сглаживающие кубические сплайны - student2.ru

Замечательным свойством эффективного импеданса является то, что он является инвариантной величиной, т.е. не меняется при повороте системы координат.

Для последующего анализа удобно пересчитать кривые компонент тензора импеданса в соответствующие кривые кажущегося сопротивления по формулам

Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru (1.2)

Кривые кажущегося сопротивления можно дополнить фазовыми кривыми

Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru (1.3)
2. Первый этап обработки данных МТЗ

При проведении полевых работ с помощью цифровых электроразведочных станций (ЦЭС-2, ЦЭС-М, ЦАИС и др.), как правило, регистрируются вариации пяти компонент МТ-поля - двух горизонтальны компонент электрического поля (Ех, Еу) и трех компонент магнитного поля (Hх, Hу, Hz). Компонента Hz при определении тензора импеданса не используется, однако дает полезную информацию о горизонтальной неоднородности разреза.

Особенностью метода МТЗ по отношению к другим методам зондирования на постоянном и переменном токе является трудоемкий процесс обработки полевых материалов, предшествующий построению кривых зондирования. Целью обработки является разделение естественного электромагнитного поля Земли, имеющего случайный характер, на гармонические составляющие, по которым находятся компоненты тензора импеданса для набора соответствующих частот.

Сглаживающие кубические сплайны - student2.ru Рис. 1. Результаты первого этапа обработки МТ-данных (компонента Сглаживающие кубические сплайны - student2.ru ).

Выделение гармонической составляющей сигнала (узкополосную фильтрацию) в частотной области можно представить как умножение спектра сигнала на частотную характеристику фильтра. Однако удобнее осуществлять узкополосную фильтрацию не в частотной, а непосредственно во временной области. При этом, согласно теореме о свертке, умножение в частотной заменяется во временной облас-ти сверткой сигнала с импульсной харак-теристикой фильтра.

В результате свертки мы находим значения гармони-ческих составляю-щих каждой компо-ненты электромаг-нитного поля для различных промежут-ков времени. Под-ставляя эти значения в импедансное соот-ношение (1.1), полу-чаем переопреде-ленную систему уравнений относительно компонент тензора импеданса, которая решается методом наименьших квадратов. Таким образом, для каждого участка записи МТ-поля получаем оценки значений компонент тензора импеданса. Разбивая всю запись на несколько участков, мы получаем оценки компонент тензора импеданса на каждом периоде (рис. 1). По этим оценкам можно определить среднее значение и среднее квадратическое отклонение компонент Сглаживающие кубические сплайны - student2.ru для данного периода.

Большинство оценок тензора импеданса группируется вблизи среднего значения. В то же время некоторые из них, отвечающие промежуткам времени с высоким уровнем помех, могут сильно отличаться от среднего значения.

Сглаживающие кубические сплайны

Задача о построении гладкой зависимости по результатам эксперимента является общей для многих разделов геофизики и других естественных наук. Она представляет собой обратную задачу, при решении которой применяется метод регуляризации. В этой задаче требуется провести кривую (сплайн) через облако полевых данных так, чтобы эта кривая, с одной стороны, проходила как можно ближе к средним для каждого периода значениям, а с другой стороны, удовлетворяла некоторому требованию (априорной информации). В качестве априорной информации принимается условие максимальной гладкости проводимой кривой частотной зависимости.

Сглаживающий кубический сплайн - это гладкая функция, на каждом отрезке представляющая собой многочлен третьей степени :

Сглаживающие кубические сплайны - student2.ru , где Сглаживающие кубические сплайны - student2.ru

В узлах сплайна (при Сглаживающие кубические сплайны - student2.ru ) функция Сглаживающие кубические сплайны - student2.ru удовлетворяет граничным условиям, которые заключаются в непрерывности самой функции Сглаживающие кубические сплайны - student2.ru , а также ее первой Сглаживающие кубические сплайны - student2.ru и второй Сглаживающие кубические сплайны - student2.ru производных. Коэффициенты многочленов Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru и Сглаживающие кубические сплайны - student2.ru находятся из условия минимума функционала

Сглаживающие кубические сплайны - student2.ru (3.1)

где Сглаживающие кубические сплайны - student2.ru - параметр сглаживания, Сглаживающие кубические сплайны - student2.ru - среднее значение сглаживаемой величины на периоде Сглаживающие кубические сплайны - student2.ru , Сглаживающие кубические сплайны - student2.ru - среднее квадратическое отклонение сглаживаемой величины на периоде Сглаживающие кубические сплайны - student2.ru , а Сглаживающие кубические сплайны - student2.ru - значение сглаживающего сплайна на периоде Сглаживающие кубические сплайны - student2.ru . Величина

Сглаживающие кубические сплайны - student2.ru

является взвешенной невязкой экспериментальных данных и сплайна, а величина

Сглаживающие кубические сплайны - student2.ru Рис. 2. Подбор оптимального параметра сглаживания.

Сглаживающие кубические сплайны - student2.ru

определяет гладкость искомого сплайна.

Таким образом, при минимизации (3.1) получается решение, в некоторой степени удовлетворяющее двум противоречащим требованиям :

a) максимальной гладкости

b) минимальной невязке сплайна и экспериментальных данных.

Критерием, определяющим то, какому из требований отдается предпочтение, является точность экспериментальных данных. Чем ниже точность, тем более гладким должен быть сплайн. Технически соотношение между этими требованиями определяет параметр Сглаживающие кубические сплайны - student2.ru , изменяющийся от 0 до 1 и играющий роль параметра регуляризации (рис. 2). При Сглаживающие кубические сплайны - student2.ru априорная информация о гладкости кривой не учитывается и аппроксимирующая кривая проходит по средним значениям экспериментальных данных. При уменьшении Сглаживающие кубические сплайны - student2.ru гладкость кривой возрастает, но начинает увеличиваться невязка сплайна и экспериментальных данных. При Сглаживающие кубические сплайны - student2.ru сплайн превращается в прямую.

Выбирая оптимальное значение параметра Сглаживающие кубические сплайны - student2.ru , лучше начать с его крайнего значения, например 1. Невязка исходной (экспериментальной) и сглаженной кривых составит при этом 0. Далее следует постепенно уменьшать Сглаживающие кубические сплайны - student2.ru до тех пор, пока невязка исходной и сглаженной кривых не станет равна дисперсии исходных данных. Соответствующее этой ситуации значение Сглаживающие кубические сплайны - student2.ru и будет оптимальным.

Наши рекомендации