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

Обложка

Цитировать

Полный текст

Аннотация

Исследованы некоторые эффекты вызванной поляризации (ВП) в рудных телах и вмещающих породах при частотах выше 1 000 000 Гц. Инструментально изучены фазы реакции комплекса сопротивлений аппаратурой системы «Геозонд». Надежность высокочастотных ВП-измерений была подтверждена независимыми заверочными горными работами. Полученные результаты подтверждают наличие ВП-эффектов в рудных телах и вмещающих породах и предполагают необходимость их учета в интерпретации электромагнитных данных, в частности в данных индукционного каротажа. Целью математического компьютерного моделирования в настоящей работе было изучение распространения высокочастотного электромагнитного поля от линейного источника тока в материальной среде путем решения уравнений Максвелла методом конечно-разностных приближений (итераций) во временной области. Обработка результатов измерений велась по известной схеме - методом итеративного подбора в автоматическом режиме, что позволяет использовать указанный метод для сквозного расчета электромагнитных полей в сложнопостроенных трехмерно неоднородных структурах.

Полный текст

Введение[**] Математическое компьютерное моделирование структуры электромагнитного поля не потеряло своей актуальности, несмотря на большое количество работ в данной области. Причиной этого является разнообразие физических сочетаний состояния вещества горных пород в реальных геологических обстановках и необходимость поиска твердых полезных ископаемых новых геолого-промышленных типов. Сочетания петрофизических свойств горных пород создают среды самых различных форм и особенностей - нелинейные, анизотропные, неоднородные, слоистые, складчатые и прочие. Целью математического компьютерного моделирования в настоящей работе было изучение распространения высокочастотного электромагнитного поля от линейного источника тока в материальной среде путем решения уравнений Максвелла методом конечно-разностных приближений (итераций) во временной области. Подобные модельные исследования ранее проводились различными авторами [1-9]. Отличие данной работы от перечисленных выше в попытке решить обозначенную проблему для высокочастотного электромагнитного поля с изучением эффектов вызванной поляризации. В таком решении видится возможность практического применения этой работы при поисках рудных месторождений (как коренных, так и россыпных). Результаты исследования Строгое решение системы уравнений Максвелла в общем виде пока неизвестно. Исключение составляют несколько специальных случаев. Преодолению проблемы препятствует удовлетворение решения граничным условиям. Однако решение задач электродинамики сейчас интенсифицируется [10-16]. Чаще всего при решениях нестационарных уравнений Максвелла применяется конечно-разностный метод во временной области (К-Р ВО). Этот метод был дан в работе Kane S. Yee в 1966 г. [17]. Решая задачу о рассеянии электромагнитного импульса на проводящей поверхности, он предложил перемежать компоненты поля в трехмерной матрице пространства так, чтобы электрические компоненты поля (Е) располагались в середине ребер матрицы, а магнитные компоненты поля (Н) располагались в серединах граней этой матрицы (рис. 1 и 2). Это позволяет иметь систему конечно-разностных уравнений, во многом удовлетворяющую граничным условиям, включающим проводящие и полупроводящие поверхности. Рис. 1. Расположение составляющих электромагнитного поля (Е и Н) в элементарной ячейке трехмерной матрицы расчетной области [Figure 1. Location of the components of the electromagnetic field (E and H) in the unit cell of the three-dimensional matrix of the calculated area] Рис. 2. Направления векторов электромагнитного поля (Е и Н) относительно элементарной ячейки трехмерной матрицы расчетной области [Figure 2. Directions of the electromagnetic field vectors (E and H) relative to the unit cell of the three-dimensional matrix of the calculated area] Поперечная электромагнитная волна [Secondary electromagnetic wave] Продольная электромагнитная волна [Primary electromagnetic wave] Рис. 3. Распределение векторов и величин амплитуд электрических (красная линия графиков) и магнитных (синяя линия графиков) составляющих переменного электромагнитного поля в продольной и поперечной электромагнитных волнах [18] [Figure 3. Distribution of vectors and amplitudes of electric (red line of graphs) and magnetic (blue line of graphs) components of the alternating electromagnetic field in the primary and the secondary electromagnetic waves [18]] В представленной на рис. 3 модели автоматически удовлетворяются граничные условия и для поверхностей раздела сред (благодаря смене положения в пространстве электрической и магнитной составляющих поля). Это позволяет использовать данный метод для сквозного расчета электромагнитных полей в сложнопостроенных трехмерно неоднородных структурах. В этой модели каждая компонента вектора магнитного поля находится между четырьмя компонентами электрического поля, и наоборот. Такое положение компонент полей дает возможность реализовать и центрально-разностный аналог пространственных производных вихревых уравнений Максвелла, и естественную геометрию воспроизведения интегрального вида закона Фарадея в простой пространственной ячейке. Трехмерные конечно-разностные уравнения во временной области для прямоугольной системы координат записывают так, как показано ниже. Уравнения Максвелла: Расположение компонент поля: i = 1,…, Nx + 1, j = 1,…, Ny, k = 1,…, Nz, i = 1,…, Nx + 1, j = 1,…, Ny, k = 1,…, Nz; i = 1,…, Nx + 1, j = 1,…, Ny, k = 1,…, Nz, i = 1,…, Nx, j = 1,…, Ny + 1, k = 1,…, Nz + 1; i = 1,…, Nx + 1, j = 1,…, Ny, k = 1,…, Nz + 1; i = 1,…, Nx + 1, j = 1,…, Ny + 1, k = 1,…, Nz. Конечно-разностные уравнения: где а ( ) берутся в местах расположения H-компонент, и где а ( ) берутся в местах расположения E-компонент. В свободном пространстве Электрические и магнитные составляющие поля вычисляются попеременно через полушаги во времени - сначала магнитная составляющая, затем электрическая. Следует отметить, что правильность вычислений сильно зависит от соотношения пространственной и временной детальности (дискретизации) трехмерных матриц с результатами моделирования. Для этого размеры ячеек трехмерных матриц необходимо выбирать кратными в 2, 3, 4 и более раз по отношению к размеру объекта, который необходимо выделить по результатам моделирования. Численное моделирование экспериментального электромагнитного рассеяния волн на объектах сложной формы, выполненное в работе [11], показало, что время счета, то есть число полных периодов падающей синусоидальной волны, необходимое для установления синусоидального режима в расчетной области, зависит от размера объекта и его диэлектрической и магнитной однородности. Это условие позволяет рассматривать поведение электромагнитных волн в горных породах по аналогии с сейсмическими волнами (рис. 3). Такой подход дает возможность вычленить и отфильтровать некоторые явления (ранее трактовавшиеся как помехи, внутренние шумы аппаратуры и пр.), возникающие при измерениях вызванной поляризации на высоких частотах от объектов сложной формы. Фильтрация ведется на основе анализа фазовых величин магнитных компонент электромагнитного поля (рис. 4). Из рис. 3 и 4 видно, что изменение амплитуды магнитной составляющей поля в продольной волне отстает от поперечной на четверть периода колебаний несущей частоты. Дальнейшая обработка результатов измерений велась по известной схеме, методом итеративного подбора в автоматическом режиме. Обычно число итераций составляло не менее 10 и не более 20 с поправочным коэффициентом 0,01-0,02. Большее число итераций или более высокая величина поправочного коэффициента не имеют смысла и часто приводят к искажению результата при визуализации (проверялось после заверочных буровых и горных работ). Рис. 4. Поляризация плоской электромагнитной волны [19] [Figure 4. The polarization of a plane electromagnetic wave [19]] Первые эксперименты по измерению эффектов вызванной поляризации электромагнитного поля на высоких частотах были начаты нами в 1995 г. на аппаратуре «Геозонд-РЧ». Недостатки аппаратуры (супергетеродинная схема приемника, узкий частотный спектр и пр.), не позволявшие в должной мере обеспечить чистоту эксперимента, привели к необходимости выбора новой аппаратурной схемы, ее реализации и отладки. Для визуализации результатов измерений сначала использовался цифровой осциллограф Visual DMM 740. Позже визуализация и запись результатов измерений производилась на виртуальный осциллограф Zetlab, инсталлируемый в ноутбук или планшет. Эксперименты проводились в полевых условиях на Южном Урале, в Монголии, Казахстане и Центральной России. Заверочные работы по проверке результатов геологической интерпретации полученных измерений выполнялись во всех случаях и показали хорошую сходимость с результатами экспериментов. Заключение Ограниченный размер статьи не позволяет привести результаты всех экспериментов. На рис. 5 представлен геофизический разрез, полученный на участке Дзун-Чулут в Монголии по одному из профилей через межгорную впадину. зона вероятного развития слабой кварц-сульфидной минерализации zone of probable development of weak quartz-sulfide mineralization Рис. 5. Результаты измерения электрических сопротивлений в разрезе по бортам долины ручья Узун-Чулун [Figure 5. Results of measuring electrical resistances in the section along the sides of the valley of the Uzun-Chulun stream] Отчетливо выделяются участки пониженных электрических сопротивлений в долине, соответствующие обводненным породам в самой нижней части разреза. В бортах долины выделяются участки с повышенными сопротивлениями, соответствующие выходам магматических пород, и участки пониженных сопротивлений, соответствующие выходам углефицированных сланцев.

×

Об авторах

Владимир Юрьевич Абрамов

Российский университет дружбы народов

Автор, ответственный за переписку.
Email: geophy-rudn@mail.ru

доцент департамента недропользования и нефтегазового дела Инженерной академии РУДН; кандидат геолого-минералогических наук

Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, д. 6

Список литературы

  1. Дмитриев В.И., Серебренников Н.Н. Численный расчет электрического поля точечного источника в слоистой среде с осесимметричным включением // Известия вузов. Геология и разведка. 1987. № 2. С. 112-117.
  2. Иванов В.Т., Кризский В.Н. Решение некоторых задач электроразведки методом граничных интегральных уравнений // Известия вузов. Геология и разведка. 1993. № 4. С. 112-117.
  3. Манаева Е.Н., Постников Е.Б. Электромагнитные процессы вблизи фронта ступенчатого импульса при электромагнитом зондировании Земли в случае слабонеоднородной проводящей среды // Известия вузов. Геология и разведка. 2002. № 3. С. 104-108.
  4. Математическое моделирование электромагнитных полей: материалы международного проекта COMMEMI. М.: Наука, 1992. 198 с.
  5. Пашкова Т.И., Постников Е.Б., Соболев С.В. О взаимодействии электромагнитных импульсов ступенчатой и прямоугольной форм с поверхностью Земли // Известия вузов. Геология и разведка. 1999. № 5. С. 115-119.
  6. Редозубов А.А. Об ориентировке полярных диаграмм ηk и ρk в анизотропной среде // Известия вузов. Геология и разведка. 1994. № 2. С. 119-122.
  7. Талалов А.Д. Модель электрических свойств водонасыщенных горных пород для частотного диапазона 103-109 Гц // Известия вузов. Геология и разведка. 2003. № 1. С. 75-80.
  8. Талалов А.Д., Даев Д.С. Лабораторные исследования частотой зависимости электрических свойств глинистых пород в диапазоне 20-109 Гц // Известия вузов. Геология и разведка. 1997. № 6. С. 123-129.
  9. Burtman V., Endo M., Zhdanov M.S., Ingeman-Nielsen T. High-frequency induced polarization measurements of hydrocarbon-bearing rocks // SEG Technical Program Expanded Abstracts. 2011. Pp. 677-681. https://doi.org/10.1190/1.3628168
  10. Жданов М.С., Спичак В.В. Математическое моделирование электромагнитных полей в трехмернонеоднородных средах. М.: Наука, 1992. 188 с.
  11. Тафлав А., Умашанкар К.Р. Численное моделирование рассеяния электромагнитных волн и вычисление эффективной площади отражения целей конечно-разностным методом во временной области // ТИИЭР. 1989. Т. 77. № 5. С. 57-76.
  12. Zhou Pei-bai. Numerical analysis of electromagnetic fields. Berlin: Springer, 1993. 406 p.
  13. Binns K.J., Lawrenson P.J., Trowbridge C.W. The analytical and numerical solution of electric and magnetic fields. John Wiley & Sons, 1992. 486 p.
  14. Booton R.C. Computation methods for electromagnetic and microwaves. John Wiley & Sons, 1992. 192 p.
  15. Paul C.R. Introduction to electromagnetic compatibility. John Wiley & Sons, 1992. 784 p.
  16. Сильвестр П., Феррари Р. Метод конечных элементов для радиоинженеров и инженеров-электриков. М.: Мир, 1986. 229 с.
  17. Yee K.S. Numerical solution of initial boundary value problems involving Maxwell’s equations in isotopic media // IEEE Trans. and Antennas and Prop. 1966. Vol. AP-14. No. 3. Pp. 302-307.
  18. Кроуфорд Ф. Волны: учебное руководство. 3-е изд., испр. М.: Наука, 1984. (Берклеевский курс физики. Т. 3).
  19. Селунский А.Б., Кузьмин А.В., Комарова Н.Ю. Тепловое поглощение плоской электромагнитной волны произвольной поляризации на синусоидальной границе водной поверхности. М.: ИКИ РАН, 2013.

© Абрамов В.Ю., 2020

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution 4.0 International License.

Данный сайт использует cookie-файлы

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

О куки-файлах