Численные методы решения интегральных уравнений в задачах электромагнитного зондирования неоднородных сред тема диссертации и автореферата по ВАК РФ 01.01.07, кандидат физико-математических наук Кругляков, Михаил Сергеевич

• Кругляков, Михаил Сергеевич
• кандидат физико-математических наук
• 2011

Оглавление диссертации кандидат физико-математических наук Кругляков, Михаил Сергеевич

Введение

Глава 1. Модифицированный метод интегральных токов в электродинамике неоднородных сред.

1.1. Вывод интегрального уравнения по объему.

1.2. Численные методы решения интегральных уравнений электродинамики

1.3. Модифицированный метод интегральных токов в уравнениях электродинамики

Глава 2. Квазитрехмерная задача электромагнитного зондирования

2.1. Постановка квазитрехмерной задачи.

2.2. Квазитрехмерный метод.

2.3. Вычислительные аспекты алгебаризации.

2.4. Численное исследование берегового эффекта.

Глава 3. Осесимметрическая задача электромагнитного зондирования

3.1. Построение интегрального уравнения.

3.2. Асимптотическое поведение поля.

3.3. Вычислительные аспекты решения интегрального уравнения в осесимметрическом случае.

3.4. Численное исследование островного эффекта

Введение диссертации (часть автореферата) на тему «Численные методы решения интегральных уравнений в задачах электромагнитного зондирования неоднородных сред»

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

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

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

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

• Растекание тока по суше вблизи источника. Точное моделирование такого растекания в настоящее время не представляется возможным, однако для морских измерений это растекание необходимо учитывать, например, заменой реального источника неким эквивалентным диполем с неизвестным моментом.

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

• В случае источника, расположенного на острове, возникает так называемый островной эффект, когда ток замыкается вокруг острова, образуя гигантскую токовую петлю.

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

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

До недавнего времени метод интегральных уравнений практически не использовался для решения крупномасштабных задач в силу технических ограничений на размер матрицы системы линейных уравнений. Однако с появлением суперкомпьютеров, позволяющих обращать полные матрицы порядок которых составляет несколько миллионов, использование этого метода становится вполне реальным. В настоящее время метод интегральных уравнений является основным для решения прямых задач электромагнитного зондирования в следующих коллективах: группа под руководством д.ф.-м.н профессора В. И. Дмитриева в МГУ, группа под руководством академика РАН Е. П. Велихова в НИЦ "Курчатовский Институт" (при взаимодейстии с группой сотрудников ИФЗ под руководством д.ф.-м.н С. М. Коротаева).

Традиционный численный метод решения интегральных уравнений электродинамики — метод коллокаций [28, 50] и его развитие — метод сжимающих интегральных уравнений [32], являются эффективными при размерах аномалии порядка нескольких длин волн в этой аномалии. Однако для крупномасштабного моделирования — когда размеры аномалии составляют несколько десятков длин волн, использование этих методов приводит к необходимости решать системы линейных уравнений высоких (106-108) порядков.

Другим подходом к численному решению интегральных уравнений электродинамики является метод интегральных токов [6, 28]. Если требуется рассчитать электромагнитное поле в точках, удаленных от высокопроводящей аномалии, расположенной в изолирующей среде, то этот метод является существенно менее требовательным к вычислительным ресурсам, нежели метод ко л локаций. Это обстоятельство позволяет использовать данный метод для ряда задач крупномасштабного моделирования, но при расчетах поля внутри аномалии метод интегральных токов предъявляет примерно те же требования к вычислительным ресурсам, что и метод коллокаций.

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

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

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

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

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

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

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

Заключение диссертации по теме «Вычислительная математика», Кругляков, Михаил Сергеевич

Заключение

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

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

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

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

Список литературы диссертационного исследования кандидат физико-математических наук Кругляков, Михаил Сергеевич, 2011 год

1. Gupta Pravin K, Niwas Sri, Chaudhary. Neeta. Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function // J. of Earth Sys.Sci. 1979. - Vol. 115, no. 3. - Pp. 1287-1305.

2. H.Takashi, M.Mori Double exponential formaulas for numerical integration // Publ.RIMS, Kyoto Univ. 1974. - Vol. 9. - Pp. 721-741.

3. J.A. Stratton. Electromagnetic Theory. — NY: MGH, 1941. — 631 pp.

4. M.S. Kruglyakov. The impedance method in the remote sounding problems // Electromagnetic Studies: Theory and Applications.Proceedings of the International Workshop in memory of M.N.Berdichevsky and P.Weidelt.- Moscow: 2010. P. 64.

5. M.S. Zhdanov, V.I. Dmitriev, A.V. Gribenko. Integral Electric Current Method in 3-D Electromagnetic Modeling for Large Conductivity Contrast // IEEE Transaction on Geoscience and remote sensing. — 2007. — Vol. 45, no. 5. Pp. 1282-1290.

6. T. Ooura. A Double Exponential Formula for the Fourier Transforms // Publ.RIMS, Kyoto Univ. 2005. - Vol. 41. - P. 971-977.

7. T.Ooura, M.Mori. The double exponential formula for oscillatory functions over the half infinite interval // J. Comput. Appl.Math. — 1991. — Vol. 38. Pp. 353-360.

8. T. Ooura, M.Mori. A robust double exponential formular for Fourier type integrals // J. Comput. Appl.Math. — 1999. — Vol. 112. Pp. 229-241.

9. V.I. Dmitriev, M.S. Kruglyakov. Fast computation of the field of a vertical magnetic dipole above a layered medium // Computational Mathematics and Modeling. 2008. - Vol. 19, no. 3. - Pp. 263-270.

10. W.L. Anderson. Numerical integration of related Hankel transforms of order 0 and 1 by adaptive digital filtering // Geophysics. — 1979. — Vol. 44. — Pp. 1287-1305.

11. Zhang Shanjie, Jin Jianming. Computation of special functions. — NY.: Wiley, 1996. 717 pp.

12. Акивис M.A., Гольберг В.В. Тензорное исчисление. — М.: Наука. Гл.ред.физ.-мат.лит., 1969. — 352 с.

13. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Модели и методы магни-тотеллурики. — М: Бином, 2001. — 630 с.

14. Бердический М.Н., Дмитриев В. И. Модели и методы магнитотеллури-ки. — М: Научный мир, 2009. — 680 с.

15. Борисенко А.И., Тарапов И.Е. Векторный анализ и начала тензорного исчисления. — 3-е. изд. — М.: Высшая Школа, 1966. — 251 с.

16. Бэйтман Г., Эрдели А. Высшие трансцендентные функции: В 3 т. Пер с англ. — М.: Наука. Гл.ред.физ.-мат.лит., 1974. — Т. 2.Функции Бесселя, функции параболического цилиндра, ортогональные многочлены. — 296 с.

17. Васильева А.Б., Тихонов А.Н. Интегральные уравнения. — 2-е, стреотип. изд. М.: ФИЗМАТЛИТ, 2002. - 360 с.

18. Ватсон Г.Н. Теория Бесселевых функций. Часть первая: Пер с анл. — М.: Изд. Иностр. Лит., 1949. — 795 с.

19. Геоэлектрические исследования с мощным источником тока на Балтийском щите:. — М: Наука, 1989. — 272 с.

20. Градштейн И. С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. — 4-е, перер изд. — М.: Физматгиз, 1963. — 1100 с.

21. Гюнтер Н. М. Теория потенциала и её применение к основным задачам математической физики. — М.: Гостехиздат, 1953. — 415 с.

22. Даугавет И.К. Теория приближенных методов. Линейные уравнения. — 2-е стер изд. СПБ: БХВ-Петербург, 2006. - 288 с.

23. Дмитриев В.И. Электромагнитные поля в неоднородных средах. — М.: Изд-во Моск.унт-та, 1969. — 131 с.

24. Дмитриев В.И., Захаров Е.В. Интегральные уравнения в краевых задачах электродинамики.Учебное пособие. — М.: Изд-во Моск.унт-та, 1967.- 316 с.

25. Дмитриев В.И., Захаров Е.В. Метод интегральных уравнений в вычислительной электродинамике. — М.: МАКС Пресс, 2008. — 316 с.

26. Дмитриев В. И, Кругляков М.С. Математическое моделирование поля вертикального магнитного диполя в двумерной неоднородной среде. // Прикладная математика и информатика. — 2010. — № 36. — С. 56-66.

27. Дмитриев В.И, Мерщикова H.A. Математическое моделирование поля кабеля в двумерной неоднородной среде // Прикладная математика и информатика. — 2010. — № 35. — С. -.

28. Дородницын A.A. Численные методы решения дифференциальных и интегральных уравнений и квадратурные формулы.Сборник статей. — М: Наука, 1964. — 354 с.

29. Жданов М.С. Теория обратных задач и реглуляризации в геофизике: Пер.с англ. — М.: Научный мир, 2007. — 712 с.

30. Изучение возможностей магнитотеллурического зондирования в Северном Ледовитом океане с помощью численного моделирования / С.М. Ко-ротаев, М.С. Жданов, М.С. Кругляков и др. // Физика Земли. — 2010.- № 9. С. 35-47.

31. Ильин В.А., Позняк Э.Г. Основы математического анализа: В 2-х т. — 4-е стер изд. М.: ФИЗМАТЛИТ, 2002. - Т. II. - 464 с.

32. Колмогоров А.Н., Фомин С. В. Элементы теории функций и функционального анализа. 7-е изд. - М.: ФИЗМАТЛИТ, 2004. - 572 с.

33. Кругляков М. С. Модифицированный метод интегральных токов в электродинамике неоднородных сред. // Прикладная математика и информатика. — 2010. — № 35. С. 25-35.

34. Кругляков М.С. О вычислительных методах в осесимметрических задачах электродинамики / / Вестник Московского Университета. Вычислительная Математика и Кибернетика. — 2010. — № 4.- С. 10-16.

35. Лифанов И. К. Метод синуглярных интегральных уравнений и вычислительный эксперемент. — М.: ТОО Янус, 1995. — 520 с.

36. M.S. Kruglyakov. Estimating the influence zone of a vertical magnetic dipole in aerial prospecting // Computational Mathematics and Modeling. ~ 2010.- Vol. 21, no. 1. Pp. 30-40.

37. Михлин С.Г. Курс математической физики. — М.: Физматгиз, 1968. — 576 с.

38. Михлин С.Г., Смолицкий X.JI. Приближенные методы решения дифференциальных и интегральных уравнений. — М.: Наука.Гл.ред.физ.-мат.лит., 1965. — 384 с.

39. Михлин С. Г. Многомерные сингулярные интегралы и интегральные уравнения. — JL: Физматгиз, 1962. — 256 с.

40. Модели морских электромагнитных зондирований / И.С. Барашков, В.И. Дмитриев, М.С. Кругляков, H.A. Мерщикова // Инновационные электромагнитные методы геофизики. — М: Научный Мир, 2011 (в печати) .

41. Мусхелишвили Н.И. Сингулярные интегральные уравнения. — 3-е, испр и дополн. изд. — М.: Наука.Гл.ред.физ.-мат.лит., 1968. — 512 с.

42. Перспективы магнитотеллурических зондирований на некоторых крупных структурах Северного Ледовитого океана / С.М. Коротаев, М.С. Жданов, М.С. Кругляков и др. // Физика Земли. — 2010. — № 9.- С. 48-54.

43. Перспективы магнитотеллурических зондирований на некоторых крупных структурах Северного Ледовитого океана / Е.П. Велихов, С.М. Коротаев, М.С. Кругляков и др. // Нефтегазопромысловый инжиниринг.- 2010. № 4. - С. 16-17.

44. Понтрягин Л. С. Обыкновенные дифференциальные уравнения. — 4-е изд. — М.: Наука.Гл.ред.физ.-мат.лит., 1974. — 331 с.

45. Садовничий В.А. Теория операторов. — 2-е изд. — М.: Изд-во Моск.унт-та, 1986. — 368 с.

46. Самохин A.B. Интегральные уравнения и итерационные методы в электромагнитном рассеянии. — М.: Радио и свзяь, 1998. — 160 с.

47. Тихонов А.Ii. Об ассимтотическом поведении интегралов, содержащих Бесселевы функции // ДАН СССР. 1959. - Т. 125, № 5. - С. 982-985.

48. Тихонов А.Н., Самарский А.А. Уравнения математической физики. -М.: Наука, 1977. 735 с.

49. Треногий В.А. Функциональный анализ. — М: Наука, 1980. — 496 с.

50. Трикоми Ф. Интегральные уравнения: Пер с англ. — М.: Изд. Иностр Лит., 1960. 301 с.

51. Хенл X., Мауэ А., Вестпфалъ К. Теория дифракции: Пер с англ. — М. МИР, 1964. 428 с.