автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Конечноэлементное моделирование гармонических электромагнитных полей от тороидальной катушки в осесимметричных и трехмерных средах

На правах рукописи

Волкова Алла Владимировна

КОНЕЧНОЭЛЕМЕНТНОЕ МОДЕЛИРОВАНИЕ ГАРМОНИЧЕСКИХ ЭЛЕКТРОМАГНИТНЫХ ПОЛЕЙ ОТ ТОРОИДАЛЬНОЙ КАТУШКИ В ОСЕСИММЕТРИЧНЫХ И ТРЕХМЕРНЫХ СРЕДАХ

Специальность 05.13.18 — Математическое моделирование, численные методы и комплексы программ

- 8 ОНТ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Новосибирск - 2009

003479017

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет»

Научный руководитель: кандидат технических наук, доцент

Персова Марина Геннадьевна

Официальные оппоненты: доктор технических наук

Могилатов Владимир Сергеевич

доктор физико-математических наук, профессор Селезнев Вадим Александрович

Ведущая организация: Сибирский научно-исследовательский институт

геологии, геофизики и минерального сырья, г.Новосибирск

Защита состоится «15» октября 2009 г. в 16 часов на заседании диссертационного совета Д 212.173.06 при Государственном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет» по адресу: 630092, г.Новосибирск, пр. К. Маркса, 20.

С диссертацией можно ознакомиться в библиотеке Новосибирского государственного технического университета.

Автореферат разослан «/"» сентября 2009 года.

Ученый секретарь диссертационного совета

Чубич В.М.

Общая характеристика работы

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

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

Довольно часто для вычисления характеристик гармонических электромагнитных полей применяют аналитические и полуаналитические методы, использующие некоторые упрощения математической модели поля (В.С.Могилатов, М.И.Эпов, Л.А Табаровский и др.). Основное достоинство этих методов состоит в том, что базирующиеся на их использовании алгоритмы позволяют достаточно быстро, с небольшими вычислительными затратами получать характеристики изучаемого поля, а к их недостаткам можно отнести невозможность учета всех особенностей решаемой задачи, в том числе и полностью адекватного учета структуры среды. Методы численного моделирования являются с этой точки зрения более универсальными, хотя в тех частных случаях, когда аналитические и полуаналитические методы применимы, могут проигрывать им из-за гораздо более высокой сложности получения конечного результата.

В настоящее время при решении дифференциально-краевых задач, описывающих различные электромагнитные поля, широко используется метод конечных элементов (МКЭ) как один из наиболее эффективных численных методов решения задач математической физики (А. ВоБэауИ, В.\УЫ1е, В.П. Ильин и др.).

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

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

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

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

Цели и задачи исследования

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

- разработаны конечноэлементные схемы моделирования осесимметричных электромагнитных полей, возбуждаемых тороидальным гармоническим источником, для двух постановок: скалярной - для (р -компоненты напряженности магнитного поля и векторной - для г- и г-компоненты напряженности электрического поля;

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

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

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

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

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

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

Научная новизна

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

нитных полей, основанные на использовании узлового (для 11^) и векторного

(для Ег, Ег) МКЭ.

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

3. Проведено сравнение тороидального и петлевого источников для некоторых задач электрокаротажа.

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

На защиту выносятся:

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

(для ЕГ,ЕЖ).

2. Математическая постановка и конечноэлементные схемы моделирования трехмерных гармонических электромагнитных полей, основанные на выделении осесимметричной части поля.

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

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

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

- источник смещается с оси скважины к ее границе;

- скважина параллельна границе раздела сред;

- источник перемещается в скважине, которая наклонена к пересекаемому ею разделу сред или слою.

Достоверность результатов

Адекватность математических моделей и разработанных конечноэлементных схем и вычислительных процедур подтверждены следующими экспериментами:

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

2. Верификация трехмерных расчетов проводилась на осесимметричных моделях путем сравнения с решениями, получаемыми на основе двумерных постановок для Нр и для (Ег, Ег).

3. Одна и та же трехмерная задача решалась в трех постановках: для вектор-потенциала А и скалярного потенциала V (узловой МКЭ), для напряжен-

ности электрического поля Е и для напряженности магнитного поля Н (векторный МКЭ).

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

Практическая значимость работы и реализация результатов

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

Личный вклад

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

В совместных публикациях автору принадлежат

- в работе [2] — построение конечноэлементных сеток (с шестигранными ячейками) и проведение расчетов электромагнитных полей;

- в работах [1,4,8,9,10] - построение конечноэлементных аппроксимаций, разработка алгоритмов построения сеток, проведение конечноэлементных расчетов и анализ результатов.

Апробация работы

Основные результаты работы были представлены на VII международном научном симпозиуме имени академика М.А.Усова 2004 «Проблемы геологии и освоения недр» (Томск, 2004 г.); Новосибирской научно-технической конференции имени Попова «Информатика и проблемы телекоммуникаций» (Новосибирск 2005 г., 2008 г.); Всероссийской научной конференции молодых ученых «Наука. Технологии. Инновации» (Новосибирск 2005 г., 2008 г.); IX международной научно-технической конференции «Актуальные проблемы электронного приборостроения» (Новосибирск 2008 г.).

Публикации

По результатам выполненных исследований опубликовано 10 работ, из них 2 статьи в изданиях, рекомендованных ВАК, 4 статьи в сборниках научных трудов и 4 статьи в материалах конференций.

Структура работы

Диссертационная работа изложена на 143 страницах и состоит из введения, четырех глав, заключения, списка использованных источников (115 наименований), приложения и содержит 66 рисунков и 18 таблиц.

Основное содержание работы

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

векторной для г- и г-компоненты напряженности электрического поля (Ег, Ez). Рассмотрены две процедуры, позволяющие решить осесимметричную задачу как скалярную для Н' и как векторную для (Er, Ez). Для обоснования

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

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

rot

— rotH

+ fiiuiH = —fj,iu> H°.

(О

Комплекснозначная вектор-функция Н° в правой части уравнения (1) определяется соотношениями:

Ни

(о,Я°,о),

1тН"

вЯ»,

2лт

R еЯ° = 0.

О вне ¿Г,

Фактически уравнение (1) является скалярным уравнением для комплекс-нозначной функции Н^ (г, 2). Оно может быть записано также в виде

-div

ill

а

+ =

(2)

igrad +

a err2 rdr{

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

rot irotE + агшЕ = icjJ , (3)

в котором для осесимметричного случая Е = (Ег,0,Ег) - комплекснозначная вектор-функция, а Ег и Ег - комплекснозначные функции координат гиг.

Плотность сторонних токов Л = (7г,0,72) описывает возбуждающий поле электрический ток в тороидальной катушке, который определяется вещественными составляющими компонент Зт и Зг.

Значение ЭДС в приемной тороидальной катушке может быть вычислено по формуле

£ =-к»шН^(глр,гпр), (4)

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

По значениям ЭДС могут быть получены значения амплитуды зы Ф принимаемых сигналов:

= ^¡(Ъе£)2 +(1т£)2,

180°

Ф =-агЫд^ 11е 1т £ ).

и фа-

(5)

(6)

Компоненты напряженности электрического поля при использовании постановки (1) могут быть вычислены по формуле

Е,

1 дНш

а дг

+

дг г

(7)

При использовании модели (3) значение Н1р {г ,гпр) для вычисления амплитуды и фазы £ может быть получено как

Нч> = — го^Е

(8)

либо вычислено по формуле

При существенных токах смещения уравнения (1) и (3) принимают вид: 1

го!

rot

а + кие

-пЛЕ

-гоШ

+ шрН. = —гшрН0,

+ [сгш — еш2 )Е = шЗ,

(9)

где е - диэлектрическая проницаемость среды.

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

На рис. 1 для тороидального источника приведены результаты расчетов полей в цилиндрически-слоистых средах, полученные с использованием разработанных методов и с использованием полуаналитического метода вычисления электромагнитного поля, который был описан в работе B.C. Могилатова.

Рис. 1. Сравнение результатов расчета предложенным численными методом (а) и полуаналитическим методом (б)

Было также проведено сравнение постановок (2) и (3) для различных по контрастности геоэлектрических моделей. Полученные результаты показывают, что при расчетах поля в слабоконтрастных по проводимости средах постановки для и для (Ег, Ег) примерно одинаковы по эффективности (т.е. соотношению точности и вычислительных затрат). Скалярная постановка в этом случае имеет некоторое преимущество в том плане, что рассмотренные в работе решатели СЛАУ (комплексный метод сопряженных градиентов СОСй и комплексный метод сопряженных невязок СОСЯ, метод бисопряженных градиентов ВСО, вМЯЕБ, локально оптимальная схема) работают для нее намного надежнее. При расчете же полей с большими контрастами в проводимости преимущество по точности имеет векторная постановка для ( Ет, Е2).

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

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

где i - мнимая единица, - магнитная проницаемость вакуума, ш - круговая частота (ш = 2т>, и - частота тока в источнике), а - удельная проводимость среды, Е' - комплекснозначная функция координат г тл. г, 3(р - комплексно-значная плотность тока, вещественная часть которой определяется через <5-функцию как /5(г — г0,г — (т.е. ток с амплитудой I сосредоточен на окружности {г = г„,2 = 20}), а мнимая часть равна нулю.

Для петлевого источника ЭДС в приемной петле определяется в виде

где гпр - радиус приемной петли, г - ее 2-координата.

Проведенные расчеты показали следующее. При радиусе зоны проникновения гзп = 0.6 м относительная аномалия в принимаемом сигнале Ф (фазе ЭДС) для тороидальной катушки и петли примерно одинаковая. Однако характер поведения фазы в зависимости от положения зонда различается довольно сильно: если Ф(г) для петли при любой из рассмотренных толщин аномального слоя имеет один экстремум в окрестности середины этого слоя, то для тороидальной катушки отклонение Ф (2) от нормального фона чаще всего имеет два экстремума, которые расходятся с увеличением толщины слоя И.

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

Для тороидальной катушки большее (по сравнению с петлей) изменение в Ф(г) дает наличие в среде менее проводящего слоя при увеличенной до 1 м зоне проникновения. Особенно ярко это проявляется при малом разносе.

В качестве примера на рис. 2 приведены графики регистрируемых в приемниках значений фазы Ф (см. соотношение (6)) при перемещении зонда с источником в виде петли и тороидальной катушки (и соответствующими приемниками) в зависимости от значений г -координаты приемника (источник расположен выше, его г -координата равна г + й) для двух вариантов радиуса зоны проникновения (гзп = 0.6 м и гт = 1 м) для геоэлектрической модели со слоем толщиной 2 м и удельной проводимостью 0.01 См/м во вмещающей среде 0.1 м при разносе между приемником и источником й = 0.5 м. В этих расчетах радиус скважины 0.1 м, удельная проводимость в ней 1 См/м, радиус зоны про-

— + (хгшЕ1П =

(И)

(12)

никновения гзп = 0.6 м и гзп = 1 м, удельная проводимость в ней 0.3 См/м; середина слоя совпадает с плоскостью г = 0.

4 г,м

1202 кГц

800 кГц

адоТсгц

-2

0.0 б)

4 г,м

1200 кГц ю

-2 0.0

1_800 кГц

400 кГц

4 г,м

_11200'кГц

ЮОцГц

400'кГц

0.0

4 Ъм

в) г)

Рис. 2. Графики Ф(г) при гш = 0.6 м для петли (а) и для тора (б) и при гзп = 1 м для петли (в) и для тора (г)

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

Рис. 3. Возможная конструкция прибора (а) и графики зависимости значений ЭДС и ее фазы (б-г) при профилировании без учета прибора (треугольнички) и с учетом прибора с двумя размерами ободка (прямоугольнички - больший ободок, кружки - меньший ободок) для моделей без слоя и со слоем

Было рассмотрено, какое влияние может оказывать конструкция каротажного прибора на сигналы в приемнике. Для этого было выполнено профилирование для среды со слоем, верхней границей которого является плоскость г = 0.5 м, а нижней - плоскость 2 = —0.5 м. Значения удельной проводимости следующие: в скважине 1 См/м, в зоне проникновения 0.3 См/м, в слое 0.01 См/м, во вмещающей среде 0.1 См/м. Радиус скважины гскв = 0.1 м, радиус зоны проникновения гзп = 0.6 м, внутренний и внешний радиусы тороидальных катушек 0.04 и 0.05 м соответственно, радиус прибора 0.05 м. Приемная катушка удалена от генераторной на 1 м.

На рис. 3 приведены графики значений реальной и мнимой компонент ЭДС и ее фазы, полученные при профилировании без прибора и с прибором, у которого проводящие металлические проводки распространяются от углублений для приемных тороидальных катушек в обе стороны (по г) на 0.02 м и 0.01 м. Как видно из графиков, прибор и конструктивные изменения в нем ока-

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

В работе рассмотрена также ситуация, когда весь прибор (за исключением, естественно, генераторной и приемной катушек) окружен металлической оболочкой, контактирующей со средой в скважине. В этом случае резко увеличивается уровень электромагнитного поля (и, соответственно, значения сигналов в приемниках), но при это реакция в фазе ЭДС и в разности фаз (для двух приемников) на слой с измененной проводимостью существенно уменьшается. Например, при профилировании вдоль скважины, находящейся в среде с удельной проводимостью 0.1 См/м со слоем с проводимостью 0.01 См/м и толщиной 1 м (г = 0.1 м, г = 0.6 м, сг = 1 См/м, <7 = 0.3 См/м), измене' СКВ 311 СКВ 311 '

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

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

Сначала рассмотрим ситуацию, когда для расчета трехмерного поля от тороидального источника применяется узловой МКЭ. Будем считать, что среда однородна по магнитной проницаемости, т.е. /i = fiQ - магнитная проницаемость вакуума. В этом случае система уравнений для расчета добавочного поля, определяемого отличием удельного сопротивления а трехмерной среды от удельного сопротивления <т° осесимметричной среды, имеет вид:

-—ДА + голтА + agrad V = (<т-а0)Е°, (13)

-div(<TgradF)-iwdiv(oA) = -div((tx - <7° )е°), (14)

Характеристики суммарного поля могут быть определены по формулам:

В = rot А - — rotE0, (15)

и

Е = -шА - grad V + Е°. (16)

В системе уравнений (13)—(14) нормальная составляющая Е° вектора напряженности электрического поля считается известной комплекснозначной

вектор-функцией, найденной при моделировании электромагнитного поля от тороидального источника тока в соответствующей <т° осесимметричной среде.

С учетом выражений (7), значение Е° в декартовой системе координат может быть получено из осесимметричного поля, описываемого осесимметричной задачей (1) или (3), по следующим формулам (г = yjx2 + у ):

= = (17)

4 ' г v ' га oz

E°{X,y,Z) = y.Er(r,z) = -Z^dJ^, (18)

"4 г га oz

Эквивалентная вариационная постановка для модели (13)—(14) имеет вид:

—JgradA • gгadФdfi + шJстАФйП + JcrgradVVdil = а п п

= f(a-a°)E0Vdfl, УФ е/, (20)

п

~~ J crgrad V • gradф<ifi -f iui J aAgгadФбíf2 =,

Q SI

= /(<?- a°)E°grad$dn, УФ e (21)

n

где Ф - вещественная пробная функция, i//1 - гильбертово пространство функций с суммируемыми с квадратом первыми производными, все компоненты А и V принадлежат ¡¿Л.

При использовании векторного МКЭ с учетом разделения искомого поля на осесимметричную и трехмерную составляющие уравнение для расчета добавочного поля Е, определяемого влиянием трехмерных неоднородностей удельного сопротивления, будет иметь следующий вид:

rot rotE j + гысгЕ = -iu: ( а - <т°) Е°, (22)

где значения напряженности нормального электрического поля Е° могут быть получены пересчетом решения #Дг,2:) осесимметричной задачи (1) или (3) в

декартовую систему координат с помощью соотношений (17)-(19).

Эквивалентная вариационная формулировка для уравнения (22) имеет

вид

J — rotE ■ rot4dil + iuif aEYdd = -iujJ"(ct-ct° )e(23) n ^ n n

Решение Е (его вещественная и мнимая компоненты) и пробные вектор-функции *Р должны принадлежать пространству ,//mt, содержащему все векторные функции 0, такие что J" 0 • QdQ < оо и J rot© • rotQdfi < оо.

n а

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

1_ J_

о

rot

—rotH I + цгсиН = —rot а

rotH^ , (24)

2D

--— lOLlJ

2D

где Н - вектор напряженности магнитного поля, а Н и - это вычисленное из решения уравнения (1) или (3) (с учетом (8)) осесимметричное поле, пересчитанное в декартовую систему координат.

Эквивалентная вариационная формулировка для уравнения (24) имеет

вид

J-rotH ■ rotVdQ + шJnWVdn =

= - j i- - rotH2D • rotT^n. (25)

nv- a

Конечноэлементные аппроксимации строились с использованием квазирегулярных сеток с ячейками в виде параллелепипедов и шестигранников. При вычислении вкладов в правые части конечноэлементных СЛАУ значения Е° и Н2Й (см. уравнения (20), (21), (23), (25)) принимались постоянными на каждом конечном элементе и равными значениям этих величин в центре элемента.

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

В качестве одного из примеров рассмотрим следующую модель. Во вмещающую среду с удельной проводимостью 0.1 См/м поместим скважину радиусом 0.1 м с удельной проводимостью 1 См/м, вокруг которой зададим зону проникновения с удельной проводимостью 0.3 См/м, а также осесимметричный (аномальный) объект в виде полого цилиндра с внутренним радиусом 0.6 м, внешним радиусом 5 м. Вдоль оси г этот объект лежит в интервале от -1 м до -1.5 м. Его удельную проводимость примем равной 0.01 См/м.

Будем полагать, что генераторная катушка содержит 100 витков, которые намотаны на стенку полого цилиндра высотой 0.01 м с внутренним радиусом 0.05 м и внешним радиусом 0.06 м.

Сравним поведение вещественной и мнимой составляющих ЭДС, полученной в результате решения соответствующих осесимметричных задач (суммарной и нормальной) и трехмерной задачи, вдоль ствола скважины при г = 0.055 м (что соответствует положению оси приемной катушки). Суммарную ЭДС обозначим через £s, а нормальную (т.е. соответствующую модели, в которой вместо аномального объекта задана вмещающая среда) - через £п. Че-

рез Е3-30 обозначим сумму £" и £а-3°, где £а-'ю - это значение ЭДС, вычисленное путем решения соответствующей трехмерной задачи.

На рис. 4. приведены графики зависимости мнимых и вещественных составляющих ЭДС £", £* и £3-зв и относительное отклонение £"-зп от £'.

-2 -1.5 -1 -0.5 7,м

-2 -1.5 -1 -0.5 г,и

в) г)

Рис. 4. Графики зависимости мнимых (а) и вещественных (б) составляющих £п, Е" и £3-30 и относительное отклонение £3-30 от Е3 (е, г)

В главе 3 также приводится сравнение вычислительных затрат при расчете двумя способами: с включением источника в трехмерную задачу и по технологии с выделением поля. С учетом осесимметричности рассматриваемой модели, расчет трехмерных полей выполнялся в четверти полной расчетной области с соответствующими условиями симметрии. В проводимых исследованиях сетки для обеих технологий были оптимизированы в сечении Х2 (в зависимости от технологии), а количество интервалов по координате (р было взято равным 4, 6 и 8. Оценка погрешности получаемых решений проводилась относительно решения осесимметричной задачи в двумерной постановке, погрешность которого не превышала 0.1 %.

В результате расчетов было получено, что для получения решения, погрешность которого не превышает 1.5%, при использовании технологии с

включением источника в трехмерную задачу количество узлов по <р должно быть не менее девяти, а сечение в плоскости ХТ должно содержать порядка 3000 узлов. При использовании же технологии с выделением поля для получения решения с такой же точности по координате ф оказалось достаточно пяти узлов, а в сечении плоскостью - около 500 узлов. Вычислительные затраты для обеих технологий приведены в табл. 1 (время счета приведено для компьютера 1те1(1*)Соге(ТМ) 2^иаё СРи@2.40СНг).

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

Таблица 1

Сравнение вычислительных затрат при использовании технологии с включением источника в трехмерную задачу и при использовании технологии с выделением поля

Решение трехмерной задачи Число узлов в сетке Число итераций решателя Время решения задачи

По технологии без включения источника в трехмерную задачу 1881 2993 6 сек

По технологии с включением источника в трехмерную задачу (4 шага по ф) 15752 12207 4мин 2сек

По технологии с включением источника в трехмерную задачу (6 шагов по ф) 22352 13452 бмин ЗОсек

По технологии с включением источника в трехмерную задачу (8 шагов по ф) 29128 17520 13мин 20сек

Также проводились сравнения вычислительных затрат при использовании модели с векторным МКЭ и модели с узловым МКЭ в диапазоне частот от 800 кГц и ниже. Было получено, что узловой МКЭ имеет очевидные преимущества от частот 200 кГц и ниже. Для частот от 800 кГц и выше более эффективным становится использование векторного МКЭ и его преимущества становятся еще более существенными при расчете полей в высококонтрастных средах.

В четвертой главе рассматриваются методы построения конечноэле-ментных сеток для моделирования трехмерных полей в нескольких характерных для задач электрокаротажа ситуациях - при смещении источника с оси скважины к ее границе, для параллельной скважине границе раздела сред (или слоя) и при пересечении скважиной раздела двух сред под некоторым углом и наклонного слоя. Приводятся примеры решения соответствующих трехмерных задач. На рис. 5 приведены примеры продольных сечений (плоскостью у = 0) конечноэлементных сеток для задач пересечения скважиной раздела двух сред и слоя, наклоненного на 10 и 65 градусов.

На рис. 6 приведены примеры поперечных сечений (плоскостью 2 = 0) конечноэлементных сеток для задач смещения источника с оси скважины к ее границе.

а) б) в)

Рис. 5. Примеры продольных сечений конечноэлементных сеток для задач пересечения скважиной раздела двух сред (а) и слоя, наклоненного на 10° (б) и 65° (в)

а)

Рис. 6. Примеры поперечных сечений конечноэлементных сеток для задач со смещением источника с оси скважины к ее границе: а) - при относительно небольшом смещении оси прибора с оси скважины;

б) - при касании прибора границы скважины

При решении задачи исследования влияния на измеряемый по стволу скважины сигнал положения границы раздела сред, параллельной оси скважины, использовались конечноэлементные сетки, содержащие порядка 20000 узлов. При этом время счета с использованием модели для векторного МКЭ в зависимости от положения границы варьировалось в интервале от 2.5 мин. до 1 мин. Для этой задачи для верификации был проведен также расчет с использованием модели узлового МКЭ. Различие в результатах не превысило 0.5 %.

При решении задачи со скважиной, которая пересекает раздел сред, каротажный зонд перемещался вдоль скважины (по оси Z) из среды с удельной проводимостью =0.2 См/м в среду с удельной проводимостью сг2 =0.05 См/м. Были проведены расчеты для трех значений частоты v тороидального и петлевого источников (400 кГц, 800 кГц и 1200 кГц) при трех значениях угла а между осью скважины и плоскостью раздела сред (10°, 35° и 90°). Результаты расчетов представлены в виде графиков зависимости фазы от z -координаты положения приемника (место пересечения оси скважины с разделом сред z=0) на рис.7. Из этих результатов видно, что для петлевого источника участок перехода значений фазы от среды с удельной проводимостью (за зоной проникновения) Cj =0.2 См/м к среде с удельной проводимостью сг2 =0.05 См/м существенно меньше при всех рассмотренных углах наклона.

Приведем еще пример решения задачи со скважиной, пересекающей под углом 10° слой с удельной электрической проводимостью 0.2 См/м (см. рис. 5,6). Удельная электрическая проводимость вмещающей среды равна 0.05 См/м. Используемые конечноэлементные сетки для рассматриваемой задачи содержали порядка 50000-60000 узлов, а время счета каждого положения источника составило порядка 10 мин. на компьютере класса Intel(R)Xeon(R) CPU Х5355 @ 2.66GHz. На рис. 8-9 приведены результаты расчетов в виде графиков фазы, реальной и мнимой частей ЭДС для петлевого и тороидального источника тока при профилировании вдоль оси скважины, пересекающей наклонный на 10° слой (приемник удален от источника на 1 м, место пересечения слоя и оси скважины z = —13 м). Для оценки влияния слоя на этих же графиках приведены значения нормального поля - для геоэлектрической модели, не содержащей слой (незакрашенные значки).

При решении задачи оценки влияния смещения зонда внутри скважины относительно ее оси параметры скважины и радиус зоны проникновения были взяты следующие: гскв = 0.1 м, (тскв = 1 См/м, гзп = 0.6 м. Удельная проводимость зоны проникновения была взята равной 0.3 См/м, а среды за зоной проникновения - 0.1 См/м. Для рассматриваемой задачи количество узлов для большинства положений зонда составило порядка 60000, а для последнего — когда зонд касается стенки скважины - примерно в два раза больше. Поскольку в данной задаче в качестве аномалии при использовании технологии выделения поля задавалась скважина и потребовались столь подробные сетки, то время счета составило от 30 мин. до 1.5 часов. В результате расчетов было получено,

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

а) б)

Рис. 7. Графики зависимости фазы ЭДС Фот значений г -координаты приемника для трех значений частоты и тороидального (а) и петлевого (б) источников и для трех значений угла а между осью скважины и плоскостью раздела сред.

Обозначения:

1 - а=10°, у=1200 кГц; 2 - а=35°, у=1200 кГц; 3 - а=90°, у=1200 кГц; 4 - а=10°, у=800 кГц; 5 - а=35°, у=800 кГц; 6 - а=90°, у=800 кГц; 7 - а=10°, у=400 кГц; 8 - а=35°, у=400 кГц; 9 - сс=90°, у=400 кГц

а) б)

Рис. 8. Графики фазы при профилировании вдоль оси скважины в среде с наклонным слоем для тороидального источника (а) и петлевого источника (б)

-10 -5 0.0 г.м -25 -20 -15 -10 -5 0.0 г,м

а) 6)

Рис. 9. Графики реальной (треугольнички) и мнимой (прямоугольнички) частей ЭДС при профилировании вдоль оси скважины в среде с наклонным слоем для тороидального источника (а) и петлевого источника (б)

Заключение

Основные результаты проведенных в диссертационной работе исследований состоят в следующем.

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

2. Разработаны и программно реализованы алгоритмы автоматического перестроения двумерных сеток при перемещении источника вдоль оси скважины. Проведено сравнение откликов от различных слоев для тороидального и петлевого источников. Реализована возможность учета влияния каротажного прибора.

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

4. Построены аппроксимации трехмерных задач на базе трех различных моделей электромагнитного поля - системы уравнений для потенциалов (А,К), векторного уравнения для напряженности электрического поля Е и

векторного уравнения для напряженности магнитного поля Н.

5. На основе метода выделения осесимметричной части поля разработаны вычислительные схемы конечноэлементного моделирования гармонических трехмерных электромагнитных полей от индукционных источников. Построенные численные процедуры протестированы на осесимметричных задачах и путем решения одной и той же трехмерной задачи в различных постановках: для (А, У), для Н и для Е.

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

Публикации автора по теме диссертации:

1. Волкова A.B. Программно-математическое обеспечение моделирования осесимметричных полей в задачах электрокаротажа для двух источников индукционного типа / A.B. Волкова, М.Г. Персова // Научный вестник НГТУ. - 2009. - № 2(35). - С. 3-12.

2. Персова М.Г. «Проблемы и возможности электромагнитных площадных геофизических зондирований при использовании телеметрических систем регистрации данных» / М.Г. Персова, Ю.Г. Соловейчик, Г.М. Тригубович, М.В. Абрамов, A.B. Зинченко (A.B. Волкова) // Автометрия. - 2007. -№ 2. - Т.43 - С. 45-54.

3. Волкова A.B. Математические модели для конечноэлементных расчетов трехмерных гармонических полей, вызываемых тороидальным током / A.B. Волкова // Сб. науч. тр. НГТУ. - 2009. - № 2(56). - С. 33-36.

4. Волкова A.B. Программный комплекс моделирования осесимметричных электромагнитных полей от тороидальной катушки в задачах электрокаротажа Тор-ЭК / A.B. Волкова, М.Г. Персова // Инновации в науке и образовании (Телеграф отраслевого фонда алгоритмов и программ). - 2008. -№9(44).-С. 61.

5. Волкова A.B. Сравнительная оценка разрешающей способности двух источников индукционного типа в осесимметричной среде / A.B. Волкова //

Материалы Российской научно-технической конференции «Информатика и проблемы телекоммуникаций». - Новосибирск: СибГУТИ, 2008 г. -Том 1.-С. 40-43.

6. Волкова A.B. Учет объектов сложной геометрии при моделировании гармонических электромагнитных полей от тороидального источника тока / A.B. Волкова // Материалы IX международной конференции актуальные проблемы элеюронного приборостроения, Новосибирск: НГТУ, 2008 г. - Т. 6. - С. 85-88.

7. Зинченко A.B. (Волкова A.B.) Моделирование трехмерного гармонического электромагнитного поля, возбуждаемого тороидальным источником / A.B. Зинченко (A.B. Волкова) // Сб. науч. тр. НГТУ. - 2006. -№ 1(43).-С. 61-66.

8. Персова М.Г. «Математическая модель и конечноэлементная аппроксимация осесимметричного поля, возбуждаемого тороидальным током» / М.Г. Персова, A.B. Зинченко (A.B. Волкова) // Сб. науч. тр. НГТУ. - 2005. -№3(41).-С. 33-38.

9. Персова М.Г. Моделирование гармонического электромагнитного поля для решения задач электрокаротажа / М.Г. Персова, A.B. Зинченко (A.B. Волкова) // Материалы Российской научно-технической конференции «Информатика и проблемы телекоммуникаций». - Новосибирск: СибГУТИ, 2005 г.-Том 1.-С. 137-139.

10-Персова М.Г. «Расчет осесимметричного электромагнитного поля, вызванного током в тороидальной катушке» / М.Г. Персова, A.B. Зинченко (A.B. Волкова) // Труды VII Международного научного симпозиума имени академика М.А.Усова «Проблемы геологии и освоения недр». - Томск: ТПУ, 2004 г.-С. 400-401.

Отпечатано в типографии Новосибирского государственного технического университета 630092, г. Новосибирск, пр. К. Маркса, 20, тел./факс (383) 346-08-57 формат 60 X 84/16 объем 1.5 пл., тираж 110 экз.. заказ № 1250 подписано в печать 08.09.09 г.

ОГЛАВЛЕНИЕ.

ВВЕДЕНИЕ.

ГЛАВА 1. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ОСЕСИММЕТРИЧНЫХ ГАРМОНИЧЕСКИХ ЭЛЕКТРОМАГНИТНЫХ ПОЛЕЙ ДЛЯ ТОРОИДАЛЬНОГО ИСТОЧНИКА.

1.1 Математические модели.

1.2 вариационные постановки.:.

1.3 Построение конечноэлементных аппроксимаций.

1.4 Методы решения конечноэлементных СЛАУ.

1.5 Сравнение с полуаналитическим подходом.

1.6 Сравнение результатов, полученных с использованием различных постановок.

1.7 Выводы.

ГЛАВА 2. ПРИМЕРЫ РЕШЕНИЯ ЗАДАЧ ЭЛЕКТРОКАРОТАЖА В ОСЕСИММЕТРИЧНЫХ ОБЛАСТЯХ.

2.1 Построение двумерных сеток.

2.2 Сравнение тороидального и петлевого источников при профилировании в осесимметричной среде.

2.3 Учет возможной конструкции прибора.

2.4 Выводы.

ГЛАВА 3. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ТРЕХМЕРНЫХ ГАРМОНИЧЕСКИХ ЭЛЕКТРОМАГНИТНЫХ ПОЛЕЙ ДЛЯ ТОРОИДАЛЬНОГО ИСТОЧНИКА.

3.1 Математические модели.

3.2 Верификация решения трехмерной задачи.

3.3 Сравнение вычислительных затрат на получение решения по технологии с выделением поля и по технологии с включением источника в трехмерную задачу.

3 А Сравнение вычислительных затрат векторного и узлового МКЭ при расчете низких частот.

3.5 Выводы.

ГЛАВА 4. КОНЕЧНОЭЛЕМЕНТНЫЕ ТРЕХМЕРНЫЕ СЕТКИ В ХАРАКТЕРНЫХ ДЛЯ ЗАДАЧ ЭЛЕКТРОКАРОТАЖА СИТУАЦИЯХ.

4.1 Построение трехмерных конечноэлементных сеток для сред с наклонными слоями.

4.2 Построение сетки в околоскважинной зоне.

4.3 Результаты решения трехмерных задач.

4.3.1 Задача со слоем, границы которого параллельны оси скважины

4.3.2 Задача со скважиной, которая пересекает раздел двух сред.

4.3.3 Задача со скважиной, которая пересекает наклонный слой.

4.3.4 Задача о влиянии смещения зонда к стенке скважины.

4.4 Выводы.

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

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

Довольно часто для вычисления характеристик гармонических электромагнитных полей применяют аналитические и полуаналитические методы, использующие некоторые упрощения математической модели поля. Эти методы разрабатывались в трудах таких ученых, как В.С.Могилатов, М.И.Эпов и др. [1,26,53,83,98,110,112,114]. Основное достоинство этих методов состоит в том, что базирующиеся на их использовании алгоритмы позволяют достаточно быстро, с небольшими вычислительными затратами получать характеристики изучаемого поля, а к их недостаткам можно отнести невозможность учета всех особенностей решаемой задачи, в том числе и полностью адекватного учета структуры среды. Методы численного моделирования являются с этой точки зрения более универсальными, хотя в тех частных случаях, когда аналитические и полуаналитические методы применимы, могут проигрывать им из-за гораздо более высокой сложности получения конечного результата.

В настоящее время при решении дифференциально-краевых задач, описывающих различные электромагнитные поля, широко используется метод конечных элементов (МКЭ) как один из наиболее эффективных численных методов решения задач математической физики [16,17,24,28,52,58,73].

Однако применение стандартных вычислительных конечноэлементных схем для моделирования трехмерных электромагнитных полей требует очень больших вычислительных затрат. Это приводит либо к большим погрешностям решения из-за использования недостаточно подробной сетки, либо к чрезмерно большому времени счета при решении важных практических задач. Предложенный в [34,35,44,58,59,60,61,62,64] подход к конечноэлементному моделированию с выделением двумерной части поля позволяет при решении многих трехмерных задач во много раз снизить вычислительные затраты и тем самым сделать эти задачи доступными для решения с необходимой точностью при относительно небольших вычислительных затратах. В предлагаемой диссертационной работе будут построены основанные на таком подходе конечноэлемент-ные схемы моделирования электромагнитных полей с гармоническими источниками типа тороидальная катушка и токовая петля.

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

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

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

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

- разработаны конечноэлементные схемы моделирования осесиммет-ричных электромагнитных полей, возбуждаемых тороидальным гармоническим источником, для двух постановок: скалярной - для </?-компоненты напряженности магнитного поля и векторной - для г - и х -компоненты напряженности электрического поля;

- разработаны алгоритмы автоматического перестроения конечноэле-ментных сеток при перемещении источника вдоль оси скважины, на основе чего реализованы методы автоматического расчета соответствующих гармонических полей с выдачей значений в одном или нескольких приемниках; выполнена программная реализация разработанных конечноэлемент-ных схем моделирования осесимметричных полей, вызываемых гармоническим током в тороидальной катушке, на основе скалярной (для Н(р ) и векторной (для Ег, Е2) постановок. Моделирование трехмерных электромагнитных полей от гармонических индукционных источников в виде тороидальной катушки и круговой петли. Для этого были решены следующие задачи:

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

Научная новизна

1. Предложены конечноэлементные схемы моделирования возбуждаемых гармоническим током в тороидальной катушке осесимметричных электромагнитных полей, основанные на использовании узлового (для Н ) и векторного (для Ег, Е2) МКЭ.

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

3. Проведено сравнение тороидального и петлевого источников для некоторых задач электрокаротажа.

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

На защиту выносятся:

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

2. Математическая постановка и конечноэлементные схемы моделирования трехмерных гармонических электромагнитных полей, основанные на выделении осесимметричной части поля.

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

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

Достоверность результатов

Адекватность математических моделей и разработанных конечноэле-ментных схем и вычислительных процедур подтверждены следующими экспериментами:

1. Решение гармонической осесимметричной задачи в цилиндрически-слоистых средах сравнивалось с результатами, полученными другими авторами с использованием полуаналитических методов [26].

2. Верификация трехмерных расчетов проводилась на осесимметричных моделях путем сравнения с решениями, получаемыми на основе двумерных постановок для Н .и для {Ег, Ег ).

3. Одна и та же трехмерная задача решалась в трех постановках: для вектор-потенциала А и скалярного потенциала V (узловой МКЭ), для напряженности электрического поля Е и напряженности магнитного поля Н (векторный МКЭ).

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

Практическая значимость работы и реализация результатов

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

Личный вклад

Автором лично разработаны и программно реализованы конечноэле-ментные схемы моделирования осесимметричных гармонических электромагнитных полей для тороидального источника, вычислительные схемы решения трехмерных гармонических электромагнитных задач на основе выделения осесимметричной части поля, выполнены расчеты осесимметричных и трехмерных полей для источников в виде тороидальной катушки и круговой токовой петли. В совместных публикациях автору принадлежат в работе [40] - построение конечноэлементных сеток (с шестигранными ячейками) и проведение расчетов электромагнитных полей. в работах [9,10,31,33] — построение конечноэлементных аппроксимаций, разработка алгоритмов построения сеток, проведение конечно-элементных расчетов и анализ результатов.

Апробация работы

Основные результаты работы были представлены на VII международном научном симпозиуме имени академика М.А.Усова 2004 «Проблемы геологии и освоения недр» (Томск, 2004 г.); Новосибирской научно-технической конференции имени Попова «Информатика и проблемы телекоммуникаций» (Новосибирск 2005 г., 2008 г.); Всероссийской научной конференции молодых ученых «Наука. Технологии. Инновации» (Новосибирск 2005 г., 2008 г.); IX международная научно-техническая конференция «Актуальные проблемы электронного приборостроения» (Новосибирск 2008 г.).

Публикации

По результатам выполненных исследований опубликовано 10 работ, из них 2 статьи в изданиях, рекомендованных ВАК, 4 статьи в сборниках научных трудов и 4 статьи в материалах конференций.

Структура работы

Диссертационная работа изложена на 143 страницах и состоит из введения, четырех глав, заключения, списка использованных источников (115 наименований), приложения и содержит 66 рисунка и 18 таблиц.

Основные результаты проведенных в диссертационной работе исследований состоят в следующем.

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

2. Разработаны и программно реализованы алгоритмы автоматического перестроения двумерных сеток при перемещении источника вдоль оси скважины. Проведено сравнение откликов от различных слоев для тороидального и петлевого источников. Реализована возможность учета влияния каротажного прибора.

3. На основе метода выделения осесимметричной части поля разработаны вычислительные схемы конечноэлементного моделирования гармонических трехмерных электромагнитных полей от индукционных источников. Построенные численные процедуры протестированы на осесимметричных задачах и путем решения одной и той же трехмерной задачи в различных постановках (для Н и для Е). Проведено сравнение этих постановок по точности и вычислительным затратам.

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

5. Построены аппроксимации трехмерных задач на базе трех различных моделей электромагнитного поля - системы уравнений для потенциалов (А, У), векторного уравнения для напряженности электрического поля Е и векторного уравнения для напряженности магнитного поля Н .

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

ЗАКЛЮЧЕНИЕ

1. Альпин JIM. Теория полей, применяемых в разведочной геофизике / JIM. Альпин, Д.С. Даев, А.Д. Каринский. М.: Недра, - 1985. - 407 с.

2. Антонов Ю.Н. ВИКИЗ обоснование и интерпретация / Ю.Н. Антонов, М.И. Эпов // Состояние и пути развития высокочастотного электромагнитного каротажа: Материалы научно-практической конференции, Новосибирск, НИЦ ОИГГиМ СО РАН, 1998. - С. 6-14.

3. Антонов Ю.Н. Эффективность и перспективы метода ВИКИЗ/ Ю.Н.Антонов, М.И.Эпов, К.Н.Каюров // Каротажник. 1998. - № 53. -С. 30-51.

4. Антонов Ю.Н. Интерпретация данных ВИКИЗ в комплексе актуальных задач нефтяной геологии / Ю.Н.Антонов, М.И.Эпов, Ю.Н. Карогодин, Д.П. Юшин, Н.К Глебочева// Каротажник. 2002. -№ 100. - С. 172-193.

5. Антонов Ю.Н. Интерпретация диаграмм ВИКИЗ / Ю.Н.Антонов // КароIтажник. 1999. - № 64. - С. 3-43.

6. Волкова A.B. Математические модели для конечноэлементных расчетов трехмерных гармонических полей, вызываемых тороидальным током / A.B. Волкова // Сб. науч. тр. НГТУ. 2009. - № 2(56). - С. 33-36.

7. Волкова A.B. Программно-математическое обеспечение моделирования осесимметричных полей в задачах электрокаротажа для двух источников индукционного типа / A.B. Волкова, М.Г. Персова // Научный вестник НГТУ. 2009. - № 2(35). - С. 3-7.

8. Ельцов И.Н. Технология интерпретации в системе МФС ВИКИЗ+ / И.Н. Ельцов // Состояние и пути развития высокочастотного электромагнитного каротажа, Новосибирск, НИЦ ОИГГМ СО РАН, 1998. С. 40-41.

9. Иванов В.Т. Методы решения прямых и обратных задач электрокаротажа / В.Т. Иванов, М.С. Масютина. М. Наука, 1983. - 143 с.

10. Ильин В.П. Методы неполной факторизации для решения алгебраических систем / В.П.Ильин. М.:Физматлит, 1995. - 288 с.

11. Ильин В.П. Методы и технологии конечных элементов / В.П. Ильин. — Новосибирск: Изд. ИВМиМГ, 2007. 371 с.

12. Кулон Ж.-Л. -К. САПР в электротехнике / Ж.-Л. Кулон, Ж.-К. Сабоннадьер; Перевод с фр. В. А. Соколова; Под ред. Э. К. Стрельбицкого. — М.: Мир, 1988.-204 с.

13. Мартаков C.B. Прямые задачи электромагнитного каротажа / C.B. Марта-ков, М.И. Эпов // Геология и геофизика. 1999. - Т.40. - №2. - С. 249-254.

14. Митчел Э. Метод конечных элементов для уравнений с частными производными / Э.Митчел, Р.Уэйт. М.: Мир, 1981. - 216 с.

15. Могилатов B.C. Возбуждение слоистых геоэлектрических сред гармоническим током / B.C. Могилатов, Г.А.Борисов // Сибирский журнал индустриальной математики. 2003. - Т.6. -№1(13). - С. 77 - 98.

16. Персова М.Г. «Математическая модель и конечноэлементная аппроксимация осесимметричного поля, возбуждаемого тороидальным током» / М.Г.Персова, А.В.Зинченко (A.B. Волкова) // Сб. науч. тр. НГТУ. 2005. -№3(41).-С. 33-38.

17. Персова М.Г. «Математическое моделирование постоянного электрического поля при заряде и измерениях в обсаженных скважинах» / М.Г.Персова, Ю.Г.Соловейчик // Сибирский журнал индустриальной математики.-2006.-№ 1(25)-С. 116-125.

18. Персова М.Г. «Моделирование трехмерных стационарных магнитных полей вертикальной электрической линии» / М.Г. Персова // Научный вестник НГТУ. Новосибирск, НГТУ, 2006. - №1(22) - С. 113-122.

19. Персова М.Г. Зондирование становлением поля трехмерных сред и проблемы интерпретации / М.Г. Персова // Сибирский журнал индустриальной математики. 2009. - № 2(38) - С. 84-96.

20. Персова М.Г. «Применение шестигранных конечных элементов для решения задач геоэлектрики» / М.Г. Персова, Д.В.Шилак // Сб. науч. тр. НГТУ. 2005. - № 2(40). - С. 9-14.

21. Писсанецки С. Технология разреженных матриц / С.Писсанецки. — М.: Мир, 1988.-410 с.

22. Рояк М.Э. Алгоритмы построения нерегулярных треугольных и тетраэдральных сеток / М.Э.Рояк, Ю.Г.Соловейчик // Сб. науч. тр. НГТУ. Новосибирск: НГТУ, 1996. -№2(4). - С. 39-46.

23. Рояк С.Х. Конечноэлементное моделирование гармонических электромагнитных полей: дис. . канд. тех. наук : 05.13.18 / С.Х.Рояк — Новосибирск, 2000. -206 с.-Библиогр.: с. 162-171.

24. Сабоннадьер Ж.-К. Метод конечных элементов и САПР: Пер. с франц. / Ж.-К. Сабоннадьер, Ж.-Л. Кулон. -М.: Мир, 1989. 192 с.

25. Светов Б.С. Основы геоэлектрики / Б.С.Светов. — М.: Издательство ЛКИ, 2008.-256 с.

26. Эпов М.И. Оценка электрической анизотропии по диаграммам ВИКИЗ / М.И. Эпов, К.В. Сухорукова, М.Н. Никитенко // НТВ "Каротажник". -1999.-№54.-С. 17-29.

27. Соловейчик Ю.Г. Метод конечных элементов для решения скалярных и векторных задач / Ю.Г. Соловейчик, М.Э. Рояк, М.Г. Персова // Учебное пособие. Сер. «Учебники НГТУ». Новосибирск: НГТУ, 2007. — 899 с.

28. Соловейчик Ю.Г. Вычислительные схемы МКЭ-моделирования трехмерных электромагнитных и тепловых полей в сложных областях: Дис. . д-ра техн. наук : 05.13.16 / Ю.Г. Соловейчик. — Новосибирск, 1997. 333 с. -Библиогр.: с. 307-333.

29. Соловейчик Ю.Г. Математическое моделирование на базе метода конечных элементов трехмерных электрических полей в задачах электроразведки / Ю.Г.Соловейчик, М.Э.Рояк, В.С.Моисеев, A.B. Васильев // Физика Земли. -1997.-№ 9. С. 67-71.

30. Соловейчик Ю.Г. Моделирование нестационарных электромагнитных полей в трехмерных средах методом конечных элементов / Ю.Г. Соловейчик, М.Э. Рояк, B.C. Моисеев, Г.М.Тригубович // Физика Земли. 1998. - № 10. -С. 78-84.

31. Соловейчик Ю.Г. Математическое моделирование процесса становления осесимметричного поля вертикальной электрической линии / Ю.Г.Соловейчик, М.Г. Персова, Г.М.Тригубович // Сибирский журнал индустриальной математики. 2003. - Т.6. - №2(14) - с. 107-125.

32. Соловейчик Ю.Г. Использование векторного МКЭ для расчёта становления осесимметричного поля вертикальной электрической линии / Ю.Г.Соловейчик, М.Г.Персова, Ю.В.Тракимус // Доклады АН ВШ январь-июнь №1(2), 2004 С. 76-86.

33. Соловейчик Ю.Г. Решение трехмерных стационарных задач электроразведки на нерегулярных параллелепипеидальных сетках / Ю.Г. Соловейчик, М.Г.Токарева, М.Г. Персова // Вестник ИрГТУ. Иркутск. — 2004. — № 1.-С. 45-60.

34. Соловейчик Ю.Г. Разработка системы интерпретации электромагнитных полей в задачах индукционной геоэлектроразведки/ Ю.Г. Соловейчик,

35. М.Э. Рояк, Г.М.Тригубович, A.B. Чернышев // Доклады СО АН ВШ январь-июнь №1(5), 2002-С. 105-114.

36. Стренг Г. Теория метода конечных элементов / Г.Стренг, Дж.Фикс М.: Мир, 1977.-350 с.

37. Тракимус Ю.В. Исследование точности расчета векторным МКЭ нестационарного осесимметричного электромагнитного поля / Ю.В.Тракимус // Сб. науч. тр. НГТУ. Новосибирск, 2005. - № 4(42). - С. 57-62.

38. Токарева М.Г. Алгоритм оптимизации прямоугольных сеток для решения задач электроразведки / М.Г.Токарева, М.Г.Персова, А.Г. Задорожный // Сб. науч. тр. НГТУ. Новосибирск, 2002. - № 2(28). - С. 41-48.

39. Воронежский государственный университет, 2003. 742 с.) - г. Мирный — 2003.-С. 301-311.

40. Уэйт Дж. Р. Геоэлектромагнетизм / Дж. Р.Уэйт. — М.: Недра, 1987. 235 с.

41. Эпов М.И. Электромагнитные методы исследования скважин / М.И. Эпов -Новосибирск: Наука, 1979. 104 с.

42. Эпов М.И. Система одномерной интерпретации данных высокочастотных индукционных каротажных зондирований / М.И. Эпов, М.Н. Никитенко // Геология и геофизика. 1993.-Т. 34. -№2. - С. 124-130.

43. Эпов М.И. Особенности высокочастотных индукционных зондирований в скважинах с горизонтальным завершением / М.И. Эпов, К.В. Сухорукова, М.Н. Никитенко, Ю.Н. Антонов // Геология и геофизика. 1998. - Т. 39. — №5. - С. 649-656.

44. Эпов М.И. Оценка электрической анизотропии по диаграммам ВИКИЗ / М.И. Эпов, К.В. Сухорукова, М.Н. Никитенко // НТВ "Каротажник" -1999.-№54, С. 17-29.

45. Эпов М.И. Характеристики пространственного разрешения зондов ВИКИЗ / М.И. Эпов, В.Н. Ульянов // Электр, и электромагн. методы исследования в нефтегазовых скважинах — Новосибирск, Изд-во СО РАН, НИЦ ОИГГГМ, 1999.-С. 196-211.

46. Эпов М.И. Интегрированная интерпретация данных ВИКИЗ по площади / М.И. Эпов, A.M. Пестерев, В.Н. Ульянов // Электр, и электромагн. методы исследования в нефтегазовых скважинах Новосибирск, Изд-во СО РАН, НИЦ ОИГТТМ, 1999. - С. 235-239.

47. Albanese R. Analysis of three-dimensional electromagnetic fields using edge elements / R.Albanese, G. Rubinacci // J.Comput.Phys. 1993. - Vol. 108. -P. 236-245.

48. Bossavit A. Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements / A. Bossavit Academic Press (Boston), 1998.

49. Bossavit A. Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism / A. Bossavit // IEE Proc., 135, Pt.A, 1988. — P.493-500.

50. Cendes Z.J. Vector finite elements for electromagnetic field calculations / Z.J.Cendes // IEEE Trans. Magn. 1991. - Vol. 27. № 5 - P. 3958-3966.

51. Cingoski V. An Improved Method for Magnetic Flux Density Visualization using Three-Dimensional Edge Finite Elements / V.Cingoski, H. Yamashita // J. Applied Phys. 1994. - Vol.75. - № 10. - P. 6042-6044.

52. Dyczij-Edlinger R. A fast vector potencial method using tangentially continuous vector finite elements / R. Dyczij-Edlinger, G. J.-F. Peng, Lee // IEEE Trans. Microwave Theory & Tech. 1998. Vol. 46. P. 863-868.

53. Fletcher R. Conjugate gradient methods for indefinite systems / R. Fletcher // In G. A. Watson, editor, Proceedings of the Dundee Biennal Conference on Numerical Analysis 1974. Springer Verlag, New York, 1975. - P. 73-89.

54. Habashy T.M. Beyond the Born and Rytov approximation: a nonlinear approach to electromagnetic scattering / T.M Habashy, R.W. Groom // Spies B.R.J.Geophys. Res. 1993. 98. B2. P. 1759-1775

55. Lanczos C. Solution of systems of linear equations by minimized iterations. / C. Lanczos // Journal of Research of the National Bureau of Standards, 1952-№49.-P. 33-53.

56. Nedelec J.C. A new family of mixed finite elements in R3 / J.C.Nedelec // Numer. Math. №50, 1986 .-P. 57-81.

57. Nedelec J.C. Mixed finite elements in R3 / J.C. Nedelec // Numer. Math. №35, 1980.-P. 315-341.

58. Rodrigue G.A vector finite element time-domain method for solving Maxwell's equations on unstructured hexahedral grids / G.Rodrigue, D.White // SIAM J. Sci. Comput. 2001. - Vol. 23; № 3. - P. 683-706.

59. Rieben R.N. A high order mixed vector finite element method for solving the time dependent Maxwell equations on unstructured grids / R.N. Rieben, G.H. Rodrigue, D.A. White.// Journal of Computational Physics vol.204, 2005. -P. 490-519.

60. Ren Z. High Order Differential Form-Based Elements for the for the Computation of Electromagnetic Field / Z. Ren , N.Ida // IEEE Transactions On Magnetics, VOL. 36, №. 4, 2000. P. 1472-1478.

61. Saad Y. GMRES: a generalized minimal residual algorithm for solving non-symmetric linear systems / Y. Saad, M. Schultz // SIAM J. Sci. Comput. 1986. -vol. 7.-P. 856-869.

62. Sogabe T. A COCR method for solving complex symmetric linear systems / T. Sogabe, S.-L. Zhang. // Journal of Computational and Applied Mathematics, 199(2007), P. 297-303.

63. Soloveichik Y.G. Iterative method for solving finite element systems of algebraic equations / Y.G. Soloveichik // Computers Math. Applic — Vol.33, 1996 — P. 87-90.

64. Torres-Verdin C. Rapid 2.5 — dimensional forward modeling and inversion via a new scattering approximation / C. Torres-Verdin, T.M. Habashy // Radio Sci. 1994.-29.-№4.-P. 1051-1079.

65. H.A. van der Vorst, A Petrov-Galerkin type method for solving Ax=b, where A is symmetric complex / H.A. van der Vorst, J.B.M. Melissen // IEEE Transaction on Magnetics, Vol.26. № 2(1990). - P. 706-708.

66. Zhdanov M.S. Quasi-linear approximation in 3-D EM modeling/ M.S.Zhdanov, S. Fang // Geophysics. 1996a. 61 P. 646-665.

67. Zhdanov M.S. 3-D quasi-linear electromsgnetic modeling and inversion/ M.S.Zhdanov, S.Fang // Three Dimensional Electromsgnetics, Seg Monograph. -Tulsa, OK, 1999-P. 233-255.

68. Zhdanov M.S. Quasi-analitical approximations and series in electromagnetic modeling /M.S. Zhdanov, V.I.Dmitriev, S.Fang, G Hursan // Geophysics. 2000. 65. P. 1746-1757.

69. Zhou L. Residual smoothing techniques for iterative methods / L. Zhou, H. Walker.// SIAM J. Sci. Computing. vol. 15. - № 2 - 1994 - P. 297-312.