УДК 550.837
ВОЗБУЖДЕНИЕ ПЕРЕХОДНЫХ ПРОЦЕССОВ ВЫЗВАННОЙ ПОЛЯРИЗАЦИИ С ПОМОЩЬЮ ИСТОЧНИКА НАПРЯЖЕНИЯ
Николай Олегович Кожевников
Институт нефтегазовой геологии и геофизики им. А. А. Трофимука СО РАН, 630090, Россия, г. Новосибирск, пр. Академика Коптюга, 3, доктор геолого-минералогических наук, профессор, главный научный сотрудник лаборатории геоэлектрики, тел. (383)333-28-16, e-mail: KozhevnikovNO@ipgg. sbras.ru
Евгений Юрьевич Антонов
Институт нефтегазовой геологии и геофизики им. А. А. Трофимука СО РАН, 630090, Россия, г. Новосибирск, пр. Академика Коптюга, 3, доктор физико-математических наук, доцент, главный научный сотрудник лаборатории геоэлектрики, тел. (383)333-28-16, e-mail: [email protected]
В традиционном методе вызванной электрической поляризации (ВП) воздействие на среду осуществляется посредством прямоугольных импульсов тока, а измеряемой величиной является напряжение. В статье на примере лабораторных измерений рассматривается альтернативный способ возбуждения с помощью прямоугольных импульсов напряжения. В этом случае откликом является переходный процесс тока, который зависит от свойств образца. Это позволяет определить поляризационные параметры образца, измерив магнитное поле тока.
Ключевые слова: вызванная поляризация, источник напряжения, источник тока.
TRANSIENT INDUCED POLARIZATION RESPONSE TO A RECTANGULAR VOLTAGE PULSE
Nikolay O. Kozhevnikov
Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, 630090, Russia, Novosibirsk, 3 Koptyug Prospect, D. Sc., Professor, Principal Scientist, tel. (383)333-28-16, e-mail: [email protected]
Evgeniy Yu. Antonov
Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, 630090,
Russia, Novosibirsk, 3 Koptyug Prospect, D. Sc., Principal Scientist, tel. (383)333-28-16, e-mail: [email protected]
In the IP method, one usually excites the earth using rectangular current pulses and records the transient voltage response. By the example of laboratory studies, we discuss energizing a geologic material with a rectangular voltage pulse. In this case, a measured quantity is the transient IP current. We show that the IP current response and the associated magnetic field depend on the sample electrical properties. This allows IP studies to be done by recording the transmitter line current or the associated magnetic field.
Key words: induced polarization, voltage source, current source.
При изучении и описании переходных процессов ВП и характеризующих их параметров обычно предполагают, что среда возбуждается прямоугольными
импульсами тока [1, 2]. Отклик среды представляет собой устанавливающееся напряжение и называется «фундаментальной» переходной характеристикой (basic transient response) [3].
Рассмотренный способ изучения поляризующихся сред не является единственным. Можно воздействовать на среду прямоугольными импульсами напряжения, а в качестве отклика регистрировать переходный процесс тока. По поводу вызванной поляризации геологических сред опубликовано множество работ, в которых обсуждаются природа этого явления и возможности его использования при решении разнообразных геологических задач. Вместе с тем в литературе, посвященной методу ВП, способ возбуждения с помощью источника напряжения обычно не рассматривается. В этой связи целесообразно рассмотреть этот способ и сравнить его с традиционным. В данной статье такое рассмотрение выполнено применительно к лабораторным измерениям ВП на образцах с использованием идеального источника напряжения Us(t). По определению внутреннее сопротивление такого источника равно нулю, поэтому выходное напряжение не зависит от параметров нагрузки.
Предположим, что к образцу (длина l, сечение S, удельная электропроводность а) подключен идеальный источник, который в момент t=0 мгновенно создает и затем поддерживает напряжение U0 (рис. 1). Такой импульс можно представить как U(t)=U0l(t), где l(t) - единичная функция Хевисайда.
Bit)
Us(t) Us(t)
е-
Рис. 1. Образец с подключенным к нему источником напряжения Us(t)
Если образец поляризующийся, его удельная электропроводность является комплексной и зависит от частоты:
а
(-) =
а
1 +
гОат)
1 + (jar)c (1 -ri)
(1)
где ю - круговая частота; у'=(-1)ш; а0 - удельная электропроводность на нулевой частоте; ц - поляризуемость; т - постоянная времени; с - показатель степени. Диапазон значений параметра ц - от нуля до единицы, т - от нуля
до бесконечности. Параметр с может принимать значения в пределах от 1 (одно время релаксации) до 0 (бесконечно широкое равномерное распределение времен релаксации).
Найдем переходную характеристику «по току» при с=1 (дебаевская релаксация). Это позволит без потери общности выводов обойтись простой математикой. Очевидно, при с=1
о («) =
О
1 +
/атГ
1 + /ат(1 -г)
(2)
В частотной области удельная электропроводность образца, ток I и напряжение и связаны соотношением (закон Ома):
I (» = Ц/7 = и (») у о
Я(/а) I
1 +
/аГт
1 + /ат(1 -г)
(3)
Заменим в (3)/а на комплексную переменную я. С учетом того, что изображение ступени напряжения равно Ц0/я,
I (у) =
Ц, у,
у I
■О
1+
1 + ят(1 -г)
(4)
Обратное преобразование Лапласа дает:
I(1) =
Ц^
_
I
О
1 + г ет(1-г)
1 -г
(5)
Очевидно, при Ц0=1В это выражение определяет переходную характеристику тока ^(1) при воздействии на образец единичной ступени напряжения. Прямоугольный импульс напряжения амплитудой Ц0 (напряжение включается в момент 1=11, выключается в момент 1=12) можно записать в виде: Цу(^=ио[1(Н) - 1(1 - 12)].
На рис. 2 показан прямоугольный импульс напряжения длительностью 1 с и переходный процесс тока ВП.
Первое слагаемое в (5) соответствует току проводимости, второе - экспоненциально убывающему току поляризации. С течением времени ток поляризации спадает до пренебрежимо малой величины и остается лишь ток проводимости. Во время действия импульса напряжения (рис. 2, а) ток в цепи равен сумме токов проводимости и поляризации (рис. 2, б). При снятии напряжения в момент 1=12 ток проводимости становится равным нулю, а направление тока поляризации изменяется на противоположное по отношению к тому, которое наблюдалось во время действия импульса напряжения.
Рассматривая поляризацию образца импульсом напряжения, необходимо отметить некоторые особенности такого способа возбуждения. Согласно (5), посто-
янная времени, характеризующая скорость спада тока (11), определяется выражением:
= т(1 ■
Т = т
(6)
100
т
(Л
20
10
-10
0.5
1.5
2.5
I С
- V (б)
1 1 1 1
0.5
1.5
2.5
I, С
Рис. 2. Импульс напряжения (а) и реакция поляризующегося образца (б). Параметры образца: а0=0.01 См/м, 1=0.1 с, ц=0.5, 5=10 см2, /=10 см
Таким образом, ц<т. Когда т ^ 0. На различие между т и 11 впервые обратил внимание Б.С. Светов с соавторами [4], сравнивая формулы Коул-Коул для удельного электрического сопротивления и удельной электропроводности. В приведенном примере (рис. 2) т=0.1 с, ц=0.5, Т1=0.05 с. Кроме этого, при начальное значение тока поляризации неограниченно возрастает.
При воздействии на образец прямоугольным импульсом тока откликом на него является переходный процесс напряжения ВП. При этом ток и создаваемое им магнитное поле Б^) не зависят от свойств образца, поскольку полностью контролируются источником. Однако при использовании источника напряжения ток, протекающий по цепи и, соответственно, через образец, зависит от а0, т и ц. Таким образом, эти параметры можно определить, измерив ток в цепи либо его магнитное поле.
Практическая реализация этого способа осложнена тем, что в переходный процесс тока вносит вклад вызванная поляризация не только среды, но и питающих электродов. Очевидно, погрешность определения поляризационных параметров среды должна быть тем меньше, чем меньше - в относительном вы-
ражении - поляризация электродов по сравнению с поляризацией среды. В работе [5] показано, что эту погрешность можно снизить путем повышения плотности тока, протекающего через электрод во время действия импульса напряжения. Однако это не решает проблему, поскольку поляризация электродов является неконтролируемым процессом и ее вклад в суммарную переходную характеристику ВП трудно оценить, особенно в полевых исследованиях. При изучении ранних стадий вызванной поляризации дополнительным источником погрешности, особенно при высоких сопротивлениях заземлений, является емкость линии [5]. Такая же проблема возникает в лаборатории, хотя в этом случае имеется больше возможностей для контроля и снижения влияния поляризации электродов.
В этой связи представляет интерес индуктивный способ возбуждения. В этом случае на среду воздействует электрическое поле, что в терминах теории электрических цепей эквивалентно возбуждению с помощью источника напряжения [3]. Откликом на воздействие является электрический ток. При этом проблема, связанная с поляризацией питающих электродов, не возникает. Заземленная линия является источником смешанного типа [6], поэтому для индуктивного возбуждения лучше использовать незаземленную петлю, связанную с подстилающей средой исключительно посредством электромагнитной индукции. При использовании незаземленной петли или рамки в качестве приемника снимается и проблема поляризации приемных электродов. Дополнительными преимуществами индуктивного источника является нечувствительность к трансверсальной анизотропии и возможность изучать объекты, перекрытые непроводящими экранами.
Общий подход к интерпретации данных импульсной индуктивной электроразведки с учетом индукционно-вызванной поляризации (ВПИ) заключается в расширении модельной базы путем использования при расчетах комплексной, частотно-зависимой удельной электропроводности. В настоящее время инверсия индукционных переходных характеристик в рамках модели горизонтально-слоистой среды с комплексной, частотно-зависимой электропроводностью представляет рутинную процедуру [7]. Уже появились публикации, в которых представлены результаты инверсии с учетом ВПИ на основе трехмерных моделей распределения поляризуемости и удельной электропроводности [8].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Комаров В.А. Электроразведка методом вызванной поляризации. - Л.: Недра, 1980. -
391 с.
2. Sumner J.S. Principles of induced polarization for geophysical exploration. - Amsterdam: Elsevier, 1976 - 227 p.
3. Wait J.R. A simple view of the I.P. influence in an inductive E.M. prospecting system // IEEE Trans. Geosc. and Remote Sens. - 1983. - Vol. 21 (24). - P. 505-506.
4. Светов Б.С., Агеев В.В., Лебедева Н.А. Поляризуемость горных пород и феномен выскоразрешающей электроразведки // Геофизика. - 1996. - №4. - С. 42-52.
5. Шестернев Д.М., Карасев А.П., Оленченко В.А. Исследование криолитозоны методом РСВП. - Новосибирск: Издательство СО РАН, 2003. - 238 с.
6. Куликов А.В., Шемякин Е.А. Электроразведка фазовым методом вызванной поляризации. - М.: Недра, 1978. - 157 с.
7. Антонов Е.Ю., Кожевников Н.О., Корсаков М.А. Автоматизированная система для интерпретации данных индукционных импульсных электромагнитных зондирований с учетом индукционно-вызванной поляризации // Геология и геофизика. - 2014. - Т. 55, № 8. -С. 1282-1293.
8. Marchant D., Haber E., Oldenburg D.W. Three-dimensional modeling of IP effects in time-domain electromagnetic data // Geophysics. - 2014. - E303-E314.
© Н. О. Кожевников, Е. Ю. Антонов, 2017