О новом подходе к геометрической 3D-инверсии данных электроразведки с восстановлением параметров сопротивления и поляризуемости во временной области в проводящих средах
Описание подхода к геометрической 3D-инверсии с одновременным разделением сигналов становления поля и вызванной поляризации в проводящих средах и восстановлением координат границ трехмерных неоднородностей, их удельного сопротивления и поляризуемости.
| Рубрика | Физика и энергетика |
| Вид | статья |
| Язык | русский |
| Дата добавления | 02.02.2019 |
| Размер файла | 2,2 M |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Новосибирский государственный технический университет
о новом ПОДХОДе К геометрической 3D-инверсии данных электроразведки с восстановлением параметров сопротивления и поляризуемости во временной области в проводящих средах
М.Г. Персова, Ю.Г. Соловейчик,
Д.В. Вагин, Ю.И. Кошкина, О.С. Трубачева
Аннотация
инверсия поляризация сигнал разделение
В работе рассматривается подход к геометрической 3D-инверсии с одновременным разделением сигналов становления поля и вызванной поляризации в проводящих средах и восстановлением координат границ трехмерных неоднородностей, их удельного сопротивления и поляризуемости. Демонстрируется работоспособность предлагаемого подхода на синтетических данных.
Ключевые слова: нелинейная 3D-инверсия, электроразведка становлением поля, вызванная поляризация.
Annotation
Marina G. Persova, Yuri G. Soloveichik, Denis V. Vagin, Yulia I. Koshkina, Olga S. Trubacheva
A new approach to geometrical 3D-inversion of geoelectric prospecting data with recovering parameters of resistance and polarizability of a conducting medium in time domain
The article considers an approach to the geometrical 3D-inversion with the simultaneous division of field formation and induced polarization in conducting media and the recovery of boundary coordinates of three-dimensional inhomogeneities, their resistivity and polarizability. The performance of the proposed approach is demonstrated on synthetic data.
Keywords: nonlinear 3D-inversion, geoelectric prospecting with field formation, induced polarization
Основная часть
Известно, что важным поисковым признаком некоторых месторождений полезных ископаемых является наличие зон повышенной поляризации (ВП). При этом одновременный подбор параметров сопротивления и поляризуемости при работе в проводящих средах зачастую приводит к сильной неоднозначности решения соответствующих обратных задач.
В данной работе предлагается подход к геометрической 3D-инверсии с восстановлением координат границ трехмерных неоднородностей, их удельного сопротивления и поляризуемости, позволяющий за счет разделения сигналов становления поля и вызванной поляризации снизить области эквивалентности. При этом 3D-инверсию предлагается выполнять в два этапа. На первом этапе осуществляется восстановление трехмерной структуры проводимости среды с одновременным разделением сигналов становления поля и вызванной поляризации. Эта задача решается путем минимизации функционала вида
, (1)
где - ошибки (невязки) в сигналах, - сигналы, зарегистрированные в -м приемнике в моменты времени , - соответствующие теоретические сигналы, полученные в результате конечноэлементного решения прямой трехмерной задачи [1]; - включающий в себя координаты границ трехмерных объектов и их удельные сопротивления вектор искомых параметров , описывающих трехмерную геоэлектрическую модель; - вектор параметров , полученных на предыдущей итерации процедуры нелинейной 3D-инверсии; ; и - параметры регуляризации; - некоторые веса, отражающие уровень погрешности при приеме сигнала в -м приемнике и масштаб изменения принимаемого сигнала по времени; - дополнительные параметры, позволяющие выделить из измеренных сигналов их составляющие, соответствующие процессам ВП; , - функции спада, - количество слоев (зон) в геоэлектрической модели с различными функциями спада.
При выполнении 3D-инверсии не только теоретические сигналы , но и производные по параметрам трехмерной модели рассчитываются путем конечно-элементного решения соответствующих трехмерных задач в специальной постановке [2].
На втором этапе осуществляется подбор параметров поляризуемости по остаточным сигналам, получившимся после вычитания из практических данных сигналов становления поля, рассчитанных для полученной на первом этапе трехмерной модели проводимости. Соответствующий математический аппарат представлен в работе [1]. Заметим, что поскольку ВП имеет более локальное влияние (при латеральном удалении приемников от зон поляризации сигналы от них очень быстро затухают), а так же, как правило, поляризующиеся объекты ищутся на определенной глубине, на которую сфокусирована установка, то поиск параметров ВП может быть достаточно успешно осуществлен с использованием ячеистых структур.
Продемонстрируем работоспособность предлагаемых методов на синтетических данных морского варианта электроразведки ДНМЭ [3,4]. Геоэлектрическая модель, для которой посредством 3D-моделирования синтезировались данные, представлена на рис. 1. Электромагнитное поле возбуждалось погруженной на 20 м горизонтальной электрической линией длиной 700 м, перемещаемой по профилю (43 положения), показанному на рис. 1. Измерения выполнялись семиэлектродной установкой (в диапазоне времен от 0.1 мс до 4 с), с помощью которой измеряются как первая, так и вторая разности потенциалов на трех электрических линиях.
Рис. 1 "Истинная" геоэлектрическая модель
На рис. 2 представлен разрез удельного сопротивления "истинной" геоэлектрической модели, а на рис. 3 - стартовая модель для первого этапа 3D-инверсии, которая представляет собой "слои" пробных объектов (объектов-"зародышей"), удельное сопротивление которых совпадает с удельным сопротивлением выбранной тем или иным способом (например, подобранной с помощью 1D-инверсии по совокупности всех данных) вмещающей среды (в том месте, где расположен пробный объект). Искомыми параметрами в нелинейной 3D-инверсии (выполняемой на первом этапе) были взяты удельные сопротивления объектов-"зародышей" верхних трех слоев (расположенных до слабопроводящего фундамента) и X-координаты латеральных границ между объектами (поскольку в рассматриваемом примере измерения выполнялись вдоль одного профиля, границы по Y для объектов-"зародышей" фиксировались). Четвертый слой объектов-"зародышей" был задан на границе с фундаментом, при этом сопротивление этих объектов было зафиксировано равным сопротивлению фундамента, а в качестве поисковых параметров для этих объектов-"зародышей" были взяты координаты двух границ по X и верхняя граница по Z. Разрезы удельного электрического сопротивления для геоэлектрических моделей, полученных после 2-й, 4-й и 11-й итераций представлены на рис. 4-6. Из сравнения рис. 2 и рис. 6 видно, что полученная в результате 3D-инверсии геоэлектрическая модель достаточно хорошо соответствует истинной. Изменение функционала невязки по итерациям представлено в табл. 1.
Таблица 1
Зависимость значения функционала невязки от номера итерации 3D-инверсии
|
№ ит. |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
|
|
0.08 |
0.04 |
0.021 |
0.012 |
0.0086 |
0.0057 |
0.0049 |
0.0035 |
0.0033 |
||
|
№ ит. |
9 |
10 |
11 |
|||||||
|
0.0028 |
0.00276 |
0.00274 |
Рис. 2 Разрез удельного сопротивления "истинной" геоэлектрической модели
Рис. 3 Разрез удельного сопротивления стартовой геоэлектрической модели
Рис. 4 Разрез удельного сопротивления геоэлектрической модели, полученной на 2-й итерации 3D-инверсии
Рис. 5 Разрез удельного электрического сопротивления геоэлектрической модели, полученной на 4-й итерации 3D-инверсии
Рис. 6 Разрез удельного электрического сопротивления геоэлектрической модели, полученной на 11-й итерации 3D-инверсии
Для полученной на этом этапе геоэлектрической модели были рассчитаны сигналы становления поля и затем вычтены из "практических" (рассчитанных с учетом ВП для истинной модели) данных. Остаточные сигналы и полученная в результате первого этапа 3D-инверсии модель сопротивления были поданы на вход процедуры, реализующей 3D-инверсию данных ВП [2]. 3D-инверсия данных ВП выполнялась в ячеистой структуре (с поиском значений поляризуемости в каждой ячейке) с учетом полученной в результате первого этапа инверсии модели сопротивления. Разрезы истинного и полученного в результате второго этапа 3D-инверсии распределения поляризуемости приведены на рис. 7. Представленные результаты свидетельствуют о том, что полученная модель поляризуемости также достаточно хорошо соответствует истинной модели.
а
б
Рис. 7 Разрез "истинного" (а) и полученного (б) в результате второго этапа 3D-инверсии распределения поляризуемости
Следует отметить, что использованная для создания синтетических данных геоэлектрическая модель была построена на основе результатов обработки практических данных, поэтому процентное соотношение влияния процессов ВП и становления поля в синтезированных кривых, которые подавались на вход разработанной процедуры 3D-инверсии, соответствовало реальным условиям морской электроразведки (в первой разности потенциалов не превышало 20-30% в области поздних времен, что было сопоставимо с влиянием в сигнале трехмерных неоднородностей сопротивления). Этот факт подтверждает корректность представленной верификации разработанной процедуры 3D-инверсии, а также работоспособность этой процедуры при ее использовании в практических работах.
Заметим также, что на практике трехмерные неоднородности сопротивления (особенно соляные купола и изменение рельефа дна) будут оказывать влияние и по третьей координате, поэтому необходимо использование площадных систем наблюдения. При этом разработанное программное обеспечение может быть использовано без изменений (естественно с добавлением объектов-"зародышей" в стартовую модель).
Библиографический список
1. М.Г. Персова, Ю.Г. Соловейчик, Г.М. Тригубович. Компьютерное моделирование геоэлектромагнитных полей в трехмерных средах методом конечных элементов. // Физика Земли -2011. vol. 47. № 2. С. 3-14
2. М.Г. Персова, Ю.Г. Соловейчик, Г.М. Тригубович, М.Г. Токарева. Методы и алгоритмы восстановления трехмерной структуры проводимости и поляризуемости среды по данным электромагнитных зондирований на основе конечноэлементного 3D-моделирования. Физика Земли. 2013. №3. C. 30-45.
3. Veeken P.C.H., Legeydo P., Davidenko Yu., Kudryavceva E., Ivanov S. and A. Chuvaev. Benefits of the induced polarization geoelectrical method to hydrocarbon exploration. //GEOPHYSICS, VOL. 74, NO. 2, MARCH-APRIL 2009; рр. B47-B59.
4. Flekkoy E., Legeydo P. Using DNME technology in Nord Sea: Case history. // The 2nd International CSEM Conference. CSEM in hydrocarbon exploration and exploitation, Oslo, May 14 - 15, 2013.
Размещено на Allbest.ru
...Подобные документы
Сущность геофизического электроразведочного метода вызванной поляризации. Аппаратура и схемы измерительных установок. Методика проведения полевых работ. Определение значений кажущихся поляризуемости и сопротивления. Интерпретация результатов измерения.
курсовая работа [2,4 M], добавлен 19.06.2012Основные принципы геометрической оптики. Изучение законов распространения световой энергии в прозрачных средах на основе представления о световом луче. Астрономические и лабораторные методы измерения скорости света, рассмотрение законов его преломления.
презентация [1,5 M], добавлен 07.05.2012Ток и плотность тока проводимости. Закон Ома в дифференциальной форме. Стороннее электрическое поле. Законы Кирхгофа в дифференциальной форме. Уравнение Лапласа для электрического поля в проводящей среде. Дифференциальная форма закона Джоуля-Ленца.
презентация [512,3 K], добавлен 13.08.2013Создание технических средств метрологического обеспечения контроля качества полупроводниковых материалов. Анализ установки по измерению удельного электрического сопротивления четырехзондовым методом. Измерение сопротивления кремния монокристаллического.
дипломная работа [1,2 M], добавлен 24.07.2012Исследование электрического поля методом зонда. Температурная зависимость сопротивления проводников и полупроводников. Определение удельного заряда электрона. Магнитное поле кругового тока и измерение горизонтальной составляющей магнитного поля Земли.
учебное пособие [4,6 M], добавлен 24.11.2012Определение эквивалентной емкости схемы и энергии, запасенной ею. Расчет эквивалентного сопротивления и токов. Описание основных характеристик магнитного поля. Расчет тока в электрической лампочке и сопротивления ее нити накала, при подключении сеть 220В.
контрольная работа [32,4 K], добавлен 17.10.2013Инструменты и методы создания объектов в среде Elcut, решение задачи и визуализации результатов расчета. Распределение токов в проводящей среде. Создание геометрической модели, состоящей из электродов, один из которых имеет потенциал "+1В", другой "-1В".
лабораторная работа [175,6 K], добавлен 26.06.2015Нетепловые процессы ЭМ полей. Основы электродинамики нетепловых процессов в материальных средах. О физическом смысле поля электромагнитного векторного потенциала. Электродинамические аспекты теории нетеплового действия электрического тока в металлах.
реферат [139,7 K], добавлен 20.01.2008Общая характеристика сопротивления материалов. Анализ прочности, жесткости, устойчивости. Сущность схематизации геометрии реального объекта. Брус, оболочка, пластина, массив как отдельные тела простой геометрической формы. Особенности напряжения.
презентация [263,5 K], добавлен 22.11.2012Составление схемы замещения электропередачи и определение ее параметров. Определение волнового сопротивления. Определение радиуса расщепления фазы. Отыскание границ области по ограничениям на радиус провода. Расчеты режима работы электропередачи.
курсовая работа [5,1 M], добавлен 31.08.2011Начальные параметры ударной волны, образующейся движением пластины. Параметры воздуха на фронте ударной волны в момент подхода волны к преграде. Расчет параметров продуктов детонации в начальный момент отражения от жесткой стенки и металлической пластины.
курсовая работа [434,5 K], добавлен 20.09.2011Регулирование скорости тягового электродвигателя при изменении магнитного поля. Пересчет характеристик при изменении магнитного поля и смешанном возбуждении. Особенности магнитного потока при шунтировании сопротивления и изменением числа витков обмотки.
презентация [321,9 K], добавлен 14.08.2013Выявление характера зависимостей составляющих основного удельного сопротивления движению при перемещении под током и без него. Использование метода имитационного моделирования. Анализ снижения аэродинамического коэффициента при уменьшении отпора хода.
отчет по практике [91,3 K], добавлен 15.07.2017Силовые линии напряженности электрического поля для однородного электрического поля и точечных зарядов. Поток вектора напряженности. Закон Гаусса в интегральной форме, его применение для полей, созданных телами, обладающими геометрической симметрией.
презентация [342,6 K], добавлен 19.03.2013Понятие электрического тока. Поведение потока электронов в разных средах. Принципы работы вакуумно-электронной лучевой трубки. Электрический ток в жидкостях, в металлах, полупроводниках. Понятие и виды проводимости. Явление электронно-дырочного перехода.
презентация [2,3 M], добавлен 05.11.2014Макроскопическое электромагнитное поле в сплошных неподвижных средах. Уравнения Максвелла в дифференциальной форме. Энергия электромагнитного поля и теорема Пойнтинга. Применение метода комплексных амплитуд. Волновой характер электромагнитного поля.
реферат [272,7 K], добавлен 19.01.2011Виды геометрической симметрии источников магнитного поля. Двойственность локальной идеализации токового источника. Опытное обнаружение безвихревого вида электромагнитной индукции. Магнито-термический эффект.
статья [57,7 K], добавлен 02.09.2007Изучение свойств рассеяния оптического излучения в конденсированных средах в результате его взаимодействия собственными упругими колебаниями. Уравнения полей и гидродинамики в жидкостях. Решение укороченных уравнений с учетом стрикционной нелинейности.
курсовая работа [2,7 M], добавлен 24.06.2015Теория диэлектрических волноводов. Анализ распространения волн в плоском оптическом волноводе с геометрической точки зрения и с точки зрения электромагнитной теории. Распределение электромагнитного поля и зависимость свойств волновода от его параметров.
курсовая работа [5,4 M], добавлен 07.05.2012Аанализ характеристик распространения электромагнитного поля с векторными компонентами электрической и магнитной напряженности, как составляющих единого электродинамического поля в виде плоских волн в однородных изотропных материальных средах.
реферат [121,1 K], добавлен 16.02.2008
