Задачи моделирования трехмерных геоэлектромагнитных полей в неоднородных средах относятся к одним из наиболее вычислительно трудоемких. Уменьшение вычислительных затрат и повышение точности решения этих задач, несомненно, является актуальной задачей, поскольку от этого напрямую зависит качество, а часто и сама возможность адекватного восстановления геоэлектрического строения среды.
Основными методами, применяемыми для численного моделирования трехмерных геоэлектромагнитных полей, являются метод интегральных уравнений (МИУ), метод конечных разностей (МКР) и метод конечных элементов (МКЭ). При этом известно, что МИУ, в разработку которого существенный вклад внесли [5, 6, 17, 18, 49, 70, 76, 77, 108, 109, 131 и др.], является вычислительно эффективным при расчетах полей в средах с небольшим количеством трехмерных объектов, с ростом же количества трехмерных неоднородностей (что более соответствует реальным условиям) этот метод начинает существенно уступать по вычислительной эффективности таким сеточным методам, как МКР (для задач геоэлектрики он разрабатывался [19, 39, 48, 68, 72, 107 и др.]) и МКЭ, в развитие которого для решения задач геоэлектромагнетизма значимый вклад внесли [14, 23, 42, 45, 64, 118, 128, 130 и др.]. МКР и МКЭ основаны на аппроксимациях дифференциальных уравнений в частных производных. При этом по вычислительной эффективности МКР, в свою очередь, довольно сильно уступает МКЭ из-за необходимости использования в МКР прямоугольных регулярных сеток — для получения решения требуемой точности это приводит к дискретным аналогам гораздо большей размерности.
В принципе, МКР и МКЭ уже достаточно давно применяются для решения трехмерных задач электроразведки. Основным преимуществом этих методов перед МИУ является разреженность матриц систем уравнений, получаемых в результате аппроксимаций соответствующих трехмерных задач. Но вместе с тем МКР и МКЭ требуют включения в расчетную область помимо самих трехмерных объектов довольно большого пространства вокруг них, и при этом для достижения хорошей точности необходимы достаточно подробные сетки, причем мелкие ячейки нужны не только в ЗБ-объектах, но и в окружающем их пространстве, особенно вблизи источников поля и приемников. В результате, при решении задач с контролируемыми источниками, если ЗБ-объекты в них дают относительно слабые отклики, МКР и МКЭ требуют очень высоких вычислительных затрат на получение решения с приемлемой точностью.
Однако разрабатываемые в последние годы модификации МКЭ с включением в него возможности выделения поля простой структуры (что является стандартным для МИГУ и не использовалось до самого последнего времени в МКЭ) — так называемого нормального поля или поля вмещающей среды — существенно расширили его возможности при решении трехмерных задач геоэлектрики. Присущие МКЭ возможности учета геометрии расчетной области любого уровня сложности и его гибкость при построении сеток в сочетании с технологиями, использующими выделение нормального поля, позволяют сделать этот метод мощным и удобным инструментом при решении трехмерных задач геоэлектромагнетизма.
При использовании МКЭ для моделирования трехмерных геоэлектромагнитных полей основная вычислительная нагрузка, как правило, приходится на решение систем линейных алгебраических уравнений (СЛАУ), получаемых в результате конечноэлементной аппроксимации исходной краевой задачи.
При решении трехмерных задач, особенно векторных, конечноэлемент-ные аппроксимации приводят к СЛАУ с разреженными матрицами очень большого размера с довольно большой шириной ленты, поэтому для решения таких СЛАУ применяют, в основном, итерационные методы, которые рассмотрены, например, в работах [7, 13, 55, 67, 88, 112, 115 и др.]. При этом скорость сходимости различных итерационных методов может сильно (в десятки, а иногда и в сотни раз) различаться, причем при решении конкретной практической задачи, как правило, заранее неизвестно, какой из методов решения может оказаться наиболее быстрым. Поэтому исследование различных методов решения СЛАУ, получаемых в результате конечноэлементной аппроксимации трехмерных задач, также является актуальной задачей.
При моделировании трехмерных векторных полей в методе конечных элементов наиболее широко распространенными являются следующие два подхода. Первый из них (узловой или скалярный МКЭ) состоит в том, что каждая из компонент искомого поля ищется в виде линейной комбинации скалярных базисных функций (которые обычно связывают с узлами конечноэлементной сетки). Второй подход (векторный МКЭ) позволяет строить математические модели непосредственно в терминах векторных переменных (используются векторные базисные функции). В условиях, когда токами смещения можно пренебречь (моделирование электромагнитных полей на не очень высоких частотах), можно использовать постановки, как для векторных, так и для скалярных базисных функций, и в зависимости от решаемой задачи преимущество может иметь как одна, так и другая постановка.
Отметим, что даже при использовании технологии выделения поля вычислительные затраты, требуемые для расчета трехмерных полей в сложных средах, содержащих много трехмерных включений, могут оставаться довольно большими. Дополнительный эффект сокращения вычислительных затрат может быть получен за счет оптимизации сеток с обоснованным изменением шагов и устранением так называемых «лишних» узлов из регулярных сеток. В этом случае в сетках могут появляться терминальные узлы, и необходима разработка специальных подходов к построению дискретных аналогов. Кроме того, вычислительные затраты при решении задач геоэлектромагнетизма с использованием МКЭ довольно сильно зависят от размера расчетной области. Ее границы должны быть удалены от мест регистрации сигнала настолько, чтобы краевые условия практически не влияли на получаемые в приемниках значения, поэтому важной является задача определения удаленных границ расчетной области (границ так называемого большого объема или «бака»). Реализация возможности определения минимально допустимого размера «бака» позволяет дополнительно сократить размер соответствующих дискретных аналогов и, соответственно, уменьшить вычислительные затраты.
В целом, все рассмотренные проблемы, связанные со снижением вычислительных затрат при конечноэлементном моделировании трехмерных геоэлектромагнитных полей, решению которых посвящена данная диссертационная работа, являются отдельными аспектами основной научной проблемы создания эффективных методов численного моделирования трехмерных гармонических по времени электромагнитных полей.
Цели и задачи исследования.
1. Построение эффективных конечноэлементных аппроксимаций для решения задачи моделирования трехмерных гармонических по времени электромагнитных полей. Для достижения данной цели сформулированы и решены следующие задачи:
— разработка и программная реализация конечноэлементных вычислительных схем моделирования гармонических по времени электромагнитных полей на трехмерных несогласованных сетках на основе узлового и векторного МКЭ с использованием выделения основной части поля;
— построение дискретных аналогов при использовании метода деревьев-кодеревьев и совместного использования скалярного и векторного потенциалов в векторном МКЭ на трехмерных сеткахреализация разработанных вычислительных схем в программном комплексе.
2. Разработка и реализация эффективных процедур решения конечноэлементных СЛАУ, возникающих при конечноэлементной аппроксимации трехмерных гармонических по времени электромагнитных полей. Для достижения данной цели сформулированы и решены следующие задачи:
— экспериментальное исследование существующих методов решения СЛАУ, возникающих при использовании различных конечноэлементных постановок для моделирования гармонических электромагнитных полей;
— разработка методов решения СЛАУ, получаемых при использовании узлового МКЭ для моделирования трехмерных гармонических электромагнитных полей;
— программная реализация методов решения конечноэлементных СЛАУ, получаемых в результате аппроксимации трехмерных гармонических по времени электромагнитных полейсравнение производительности разработанных программ с существующими математическими библиотеками.
На защиту выносятся.
1. Математические постановки и схемы конечноэлементного моделирования трехмерных гармонических по времени электромагнитных полей с использованием узлового и векторного МКЭ с выделением части поля на сетках с терминальными узлами.
2. Методы решения СЛАУ, возникающих при моделировании трехмерных электромагнитных полей на высоких частотах с использованием векторного МКЭ и на низких частотах с использованием узлового МКЭ.
3. Программные модули конечноэлементного моделирования трехмерных гармонических геоэлектромагнитных полей и решения СЛАУ в подсистеме ЭР-МТЗ программного комплекса ОеоЕМ.
4. Результаты использования разработанных методов и программного комплекса для решения задач геоэлектромагнетизма.
Достоверность результатов.
Вычислительные схемы конечноэлементного моделирования трехмерных электромагнитных полей протестированы путем сравнения с решением задач меньшей размерности, с решениями, полученными другими авторами, а также путем сравнения решений одних и тех же задач с использованием представлений электромагнитного поля через различные потенциалы в векторном и в узловом МКЭ.
Научная новизна.
1. Разработана и реализована вычислительная схема для моделирования трехмерных магнитотеллурических полей в математической постановке для векторного МКЭ с использованием векторного и скалярного потенциалов, технологии выделения поля вмещающей горизонтально-слоистой среды и нерегулярных сеток с терминальными узлами.
2. Разработана и реализована вычислительная схема для моделирования трехмерных гармонических геоэлектромагнитных полей с использованием метода деревьев-кодеревьев и технологии выделения части поля.
3. Разработан новый метод решения СЛАУ, возникающих при моделировании трехмерных гармонических по времени электромагнитных полей на низких частотах с использованием узлового МКЭ — метод минимизации невязки с использованием блочной релаксации (МНБР). Практическая значимость.
Предлагаемые в данной работе способы построения дискретных аналогов и методы решения СЛАУ реализованы в подсистемах ЭР-МТЗ и ЭРИАН программного комплекса ОеоЕМ, который использовался при проектировании поисково-разведочных работ и интерпретации практических данных. Публикации.
По результатам выполненных исследований опубликовано 29 работ, в том числе 7 статей в журналах, входящих в перечень изданий, рекомендуемых Высшей аттестационной комиссией для защиты диссертаций, 2 статьи в сборнике научных трудов, 18 публикаций в материалах конференций, 2 свидетельства о государственной регистрации программы для ЭВМ. Апробация работы.
Основные результаты работы были представлены и докладывались на: международной научно-практической конференции «Геобайкал 2012» (Иркутск, 2012) — всероссийской школе-семинаре им. М. Н. Бердичевского и Л. Л. Ваньяна по электромагнитным зондированиям Земли, ЭМЗ-2011 (Санкт-Петербург, 2011) — рабочем семинаре «Наукоемкое программное обеспечение» в рамках Ершовской конференции по информатике (РБГП) (Новосибирск, Академгородок, 2011) — IX и X международной конференции «Актуальные проблемы электронного приборостроения» АПЭП-2008, АПЭП-2010 (Новосибирск, 2008,.
2010) — всероссийской конференции по вычислительной математике КВМ-2011 (Новосибирск, Академгородок, 2011) — научной конференции Лаборатории НГУ-Интел «Технологии высокопроизводительных вычислений» (Новосибирск, Академгородок, 2009, 2010) — IV международной научно-технической конференции молодых специалистов, аспирантов и студентов «Математическое и компьютерное моделирование естественнонаучных и социальных проблем» (Пенза, 2010) — всероссийской научной конференции молодых ученых «Наука. Технологии. Инновации» (Новосибирск, 2008, 2009, 2010) — российской научно-технической конференции «Информатика и проблемы телекоммуникаций» (Новосибирск, 2007, 2008, 2010, 2011) — VI Всероссийской научно-практической конференции студентов, аспирантов и молодых ученых «Технологии Microsoft в теории и практике программирования» (Томск, 2009).
Личный вклад.
Разработан метод решения СЛАУ, возникающих при моделировании гармонических по времени трехмерных электромагнитных полей на низких частотах с использованием узлового МКЭ. Разработаны, реализованы и апробированы программные модули конечноэлементного моделирования и решения СЛАУ в программном комплексе GeoEM. Разработан и реализован алгоритм подбора размера расчетной области для решения трехмерных задач магнитотеллурики. Выполнена компьютерная реализация методов решения конечноэлементных СЛАУ в многоядерных системах с общей памятью и реализация запуска задач магнитотеллурических зондирований по частотам в системах с распределенной памятью. Разработаны и исследованы методы построения дискретных аналогов при моделировании трехмерных электромагнитных полей с использованием метода деревьев-кодеревьев.
В совместных публикациях автору принадлежат следующие результаты. В работах [94, 117, 119] автором проведены расчеты трехмерных электромагнитных полей, в работах [92, 98, 120, 122, 132, 133] автором разработаны методы решения конечноэлементных СЛАУ, в работах [116, 123] автором разработаны и реализованы модули распараллеливания расчетов трехмерных электромагнитных полей и выполнены вычислительные эксперименты, в [90, 97, 113] проведены сравнения итерационных методов решения СЛАУ.
Работа выполнялась при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009;2013 годы (автор принимал участие в 9 проектах, в том числе в 6 проектах являлся руководителем) и ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007;2013 годы» (в 3 проектах). Работа выполнялась в рамках государственных заданий высшим учебным заведениям в части проведения НИР (шифр заявки 8.874.2011), а также была поддержана грантом Фонда содействия отечественной науке в рамках программы «Лучшие аспиранты РАН» в 2010 г.
Структура работы.
Диссертационная работа состоит из введения, четырех глав, заключения, списка использованных источников (133 наименования), приложения и содержит 40 рисунков и 52 таблицы. Общий объем работы -191 страница.
ЗАКЛЮЧЕНИЕ
.
