Использование методов численного дифференцирования для обработки геофизических данных электромагнитных зондирований с целью повышения разрешающей способности

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

Рубрика Геология, гидрология и геодезия
Вид статья
Язык русский
Дата добавления 02.02.2019
Размер файла 75,1 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://allbest.ru

Московский государственный университет имени М.В. Ломоносова

Факультет вычислительной математики и кибернетики

УДК 550.837

Использование методов численного дифференцирования для обработки геофизических данных электромагнитных зондирований с целью повышения разрешающей способности

Ж.Г. Ингтем, В.И. Дмитриев

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

Интегральный сплайн позволяет, как сглаживать геофизические данные, полученные с погрешностью, так и одновременно вычислять производную исследуемого процесса. Для того чтобы выделить с хорошей точностью области, где расположена неоднородность, достаточно исследовать поведение интегрального сплайна построенного по полученным геофизическим данным, а также поведение его первой и второй производных. Построение интегрального сплайна заключается в определении сплайна -го порядка через производную -го порядка , которая считается кусочно-постоянной функцией. Исходя из того, что необходимо исследовать поведение производных до второго порядка включительно, достаточно построить сплайн второго порядка.

Для построения интегрального сплайна решается следующая задача

(1)

где - измеренные значения в точках , - длина отрезка на котором строится сплайн, - количество измеренных значений, является кусочно-постоянной функцией на сетке построения сплайна , - шаг сетки, ; при , ; .

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

(где параметр регуляризации, зависящий от погрешности измерения ) позволяет устойчиво определять вектор параметров , и, следовательно, устойчиво построить интегральный сплайн .

Первая производная интегрального сплайна является кусочно-линейной, а вторая кусочно-постоянной. Для того, чтобы иметь гладкое представление первой и второй производных необходимо построить интегральный сплайн по значениям этих производных взятых в середине отрезков , .

Вычислительный эксперимент на модели МТ-полей в двумерных неоднородных средах проводился следующим образом:

- Расчет кажущегося сопротивления для неоднородной среды при фиксированной частоте в зависимости от на земной поверхности;

- Внесение погрешностей в кажущееся сопротивление случайным образом;

- Вычисление и с помощью интегрального сплайна, и их анализ.

На рис. 1 приведены графики профильных кривых , и при для случая двух неоднородностей находящихся в двухслойной среде ; ; . Проводимость неоднородностей , их глубина залегания , неоднородности расположены при и , т.е. расстояние между неоднородностями . Погрешность в равна от среднего аномалии . Из рис.1а видно, что по различить аномалии от каждого тела практически невозможно (точками показаны неточные данные). Первая производная (рис. 1б) уже разделяет неоднородности, а вторая производная (рис.1в) даёт чёткое положение каждой неоднородности (жирные горизонтальные отрезки обозначают расположение неоднородностей по оси, но никак не связанны с глубиной из-за масштаба). интегральный сплайн геофизический зондирование

а) б) в)

Рис. 1. Пример разделения аномалий от двух тел

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

Библиографический список

1. Дмитриев В.И., Мерщикова Н.А. Использование дифференциальных характеристик магнитного поля при решении обратной задачи электромагнитного зондирования двумерно-неоднородных сред// Физика Земли 1991 № 6 с. 59-64.

2. Dmitriev V.I., Dmitrieva I.V., Ingtem J.G. Integral form of the spline function in approximation problems // Computational Mathematics and Modeling, Consultants Bureau (United States) 2013, том 24, № 4, с. 488-497.

Аннотация

УДК 550.837

Использование методов численного дифференцирования для обработки геофизических данных электромагнитных зондирований с целью повышения разрешающей способности. Ж.Г. Ингтем, В.И. Дмитриев. Московский государственный университет имени М.В. Ломоносова факультет Вычислительной математики и кибернетики.

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

Ключевые слова: Численное дифференцирование, Интегральный сплайн, магнитотеллурическое зондирование.

Abstract

Using methods of numerical differentiation for processing geophysical data electromagnetic soundings for the purpose of raising the resolution. J.G. Ingtem, V.I. Dmitriev

The application of differential methods in electrical prospecting of mineral resources allows to increase the resolving power for the interpretation. This work is focused on the application of the integral form of the spline function in the treatment of electromagnetic sounding data.

Keywords: Numerical differentiation, Integral spline, magnetotelluric sounding.

Размещено на Allbest.ru

...

Подобные документы

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.