Способ морской геоэлектроразведки с фокусировкой электрического тока
Изобретение относится к области геофизических исследований, а более конкретно к способам морской геоэлектроразведки с использованием контролируемых искусственных источников электромагнитного поля. Сущность: при помощи дипольного источника возбуждают электромагнитное поле в толще исследуемой среды, посылая в нее прямоугольные импульсы тока с паузами между ними. Вдоль заданного профиля проводят в течение импульса тока геометрическое зондирование и в течение паузы зондирование на переходных процессах. Измерения проводят при помощи зафиксированной на морском дне измерительной установки, состоящей из пяти электродов: центрального и расположенных равноудаленно вокруг него четырех по вершинам квадрата, две противоположные стороны которого параллельны оси профиля. Измеряют первые разности электрических потенциалов между внешними электродами и центральным. При прохождении дипольного источника через разные точки обеспечивают эквипотенциальность замкнутой линии, проходящей через четыре внешних электрода измерительной установки и, тем самым, исключают внутри области, ограниченной контуром этой линии, горизонтальную составляющую плотности тока в каждой точке зондирования. Из значений измеренных первых разностей электрических потенциалов рассчитывают два множества нормированных интерпретируемых электрических параметров, которые не подвержены боковому влиянию трехмерных геологических неоднородностей, находящихся вне точки зондирования. С использованием полученных параметров находят модель среды и строят временные разрезы этой модели по электропроводности элементов среды, коэффициенту вызванной поляризации и постоянной времени спада разности потенциалов вызванной поляризации. Технический результат: возможность полного исключения в точке зондирования горизонтальной составляющей плотности тока и, тем самым, исключения бокового влияния, что дает возможность при больших глубинах моря находить и с высокой контрастностью оконтуривать скопления углеводородов при сравнительно простой электрической схеме измерителя. 3 з.п. ф-лы, 6 ил.
Изобретение относится к области геофизических исследований, а более конкретно к способам морской геоэлектроразведки с использованием контролируемых искусственных источников электромагнитного поля, и предназначено для поисков и оконтуривания скоплений углеводородов (УВ) на континентальном шельфе на основе фокусировки электрического тока и раздельного определения и картирования свойственных каждому из элементов (горизонтов) толщи осадочных отложений горных пород и обладающих аномальными значениями в зонах скопления УВ следующих трех, необходимых для решения поставленной задачи, электрофизических параметров: удельной электропроводности, вызванной поляризации и постоянной времени спада разности потенциалов вызванной поляризации. В совокупности эти три электрофизические параметра позволяют отличить залежь УВ от вмещающих пород.
Способы геоэлектроразведки, в том числе и морской, с контролируемым возбуждением исследуемой среды электрическим током известны (методы сопротивлений на постоянном и переменном токе, в том числе на переходных процессах, на основе дипольно-осевой установки ABMN). Но они предназначены для определения только одного электрофизического параметра из перечисленных выше трех, а именно кажущегося электрического сопротивления, и не конкретного желаемого участка исследуемого пространства, а всего пространства, куда проникает распространяющийся по закону диффузии электрический ток источника. Этого далеко недостаточно для поисков и оконтуривания ныне предлагаемых для геологической разведки скоплений УВ, залегающих на глубине свыше 1000 м ниже уровня морского дна.
Отметим, что согласно теории и практике геоэлектроразведки на основе дипольно-осевой установки ABMN можно выделить на глубинах в пределах около 1000 м (но не оконтурить) высококонтрастные залежи УВ (толщиной не менее 50-100 м и с электрическим сопротивлением, превышающим сопротивление окружающих пород в 50 и более раз). Но такие залежи встречаются редко и к настоящему времени в большинстве своем уже разведаны.
Попытки зондирования на основе дипольно-осевой установки ABMN с контролируемым возбуждением тока с целью поисков скоплений УВ за последние 90 лет многократно предпринимались в различных вариантах, но из-за их низкой эффективности для поисков скоплений УВ не нашли широкого применения. Дипольное электрическое зондирование обладает высоким разрешением на глубинах до 500 м и на практике успешно используется для поисков залегающих воды и руд, а также в инженерных изысканиях и археологии.
По результатам зондирования на основе дипольно-осевой установки ABMN с контролируемым возбуждением тока вычисляется кажущееся электрическое сопротивление ρ с помощью универсальной формулы:
где
I - измеряемый скачок силы тока в дипольном электрическом источнике;
ΔU - измеряемая разность потенциалов на концах приемных электродов MN;
К - геометрический коэффициент зондирующей установки («Электроразведка», Справочник геофизика. Ред. А.Г.Тархов М.: Недра, 1980, с.237 и с.422-406) [1].
При таком подходе, который применяется обычно при всех традиционных способах определения электрического сопротивления в геоэлектроразведке с контролируемым источником тока, вследствие распространения тока по закону диффузии, получают лишь суммарные сведения о всех элементах строения исследуемой среды, в которой развивается электрическое поле, так как в ней распределение в пространстве измеряемого тока I источника ничем не контролируется. И информации об этом распределении в реально существующих трехмерно-неоднородных средах нет.
В числе этих способов сопротивлений в морской геоэлектроразведке используется донное геометрическое зондирование низкочастотным переменным током на основе дипольно-осевой установки ABMN, названное Норвежской компанией EMGS «(SBL - sea bed logging)» (S.E.Johansen, H.E.F.Amundsen, T.Rosten, S.Ellingsrud, T.Eidesmo and A.H.Bhuyian. Subsurface hydrocarbons detected by electromagnetic sounding. First Break, 23, March 2005, 31-36) [2]. Он также известен под названием «controlled soutce electromagnetic sounding (CSEM)» или «offshore hydrocarbons mapping (OHM)».
В публикации [2] описано донное профилирование установкой SBL с разносом 6,5 км между АВ и MN через одно из крупнейших в мире морское газовое месторождение «Troll West Gas Province (TWGP)». Месторождение залегает на глубине примерно от 1070 м до 1230 м относительно поверхности морского дна и простирается по профилю от 10,5 км до 17,5 км (средняя глубина моря - 330 м). Оно отличается аномально высоким электрическим сопротивлением (250-1000 Ом·м) по сравнению с сопротивлением вмещающих пород (1-2 Ом·м) и большой толщиной - 160 м. (Отметим, что подавляющее большинство залежей УВ обладают электрическим сопротивлением, величина которого по сравнению с сопротивлением законтурных отложений отличается в пределах от пяти до тридцати раз, а их толщины находятся в пределах от нескольких метров до нескольких десятков метров.)
Тем не менее, из описания [2] и приложенных к нему рисунков следует, что кривая профилирования способом SBL искажается влиянием этого геометрически крупного и высококонтрастного по электрическому сопротивлению месторождения уже на шестом километре при подходе к нему измерительной установки. А находящиеся справа от него четыре структуры между 17,5 км и 24 км задавлены его боковым влиянием, и задача определения наличия или отсутствия в них углеводородов не решается. О выделении границ контура залежей здесь говорить не приходится, так как переход от законтурной зоны до основной залежи на кривой сопротивлений (реально измеренной и теоретической) растягивается более чем на 3 км. Для такой яркой залежи УВ по ее геометрическим и электрофизическим параметрам это крайне скудная информация.
Кроме этого при малых глубинах моря способ CSEM (OHM, SBL) не дает даже такой скудной информации о залежи УВ.
Таким образом, способы сопротивлений, в том числе и геометрическое зондирование без фокусировки электрического тока, малопригодны для поисков и оконтуривания большинства встречающихся в геологических недрах залежей УВ на глубинах свыше 1000 м относительно уровня морского дна, по крайней мере, по трем причинам: первая - регистрируется только один из необходимых для этой цели трех электрофизических параметров исследуемой среды (кажущееся сопротивление), что далеко не всегда достаточно для выявления скоплений углеводородов в толще геологических пород; вторая - регистрируемый параметр (кажущееся сопротивление) для той же цели слишком грубый, так как им из-за отсутствия вертикальной фокусировки электрического тока регистрируется сопротивление объема всех геологических объектов исследуемой среды, в которой развивается электрическое поле источника тока, то есть результаты измерений существенно искажены влиянием боковых геологических неоднородностей; третья - не регистрируется параметр вызванной поляризации η, обладающий также наряду с электрическим сопротивлением аномальным значением в зоне скоплений УВ.
В геоэлектроразведке влияние боковых геологических неоднородностей называют боковым влиянием, которое при поисках скоплений углеводородов на глубинах их залегания свыше одного километра сказывается на результаты измерений по латерали за несколько километров от точки зондирования, и эти скопления становятся практически малозаметными.
Искажающее боковое влияние вызвано тем, что электрический ток в пространстве не распространяется, например, в виде акустического луча, как в сейсморазведке, а растекается по закону диффузии по направлению наименьшего электрического сопротивления. И если дипольное электрическое зондирование на малых глубинах исследований обладает высоким разрешением, то с увеличением глубины исследования, вследствие рассеяния тока источника, оно теряет это разрешение.
Морское ДЭЗ из-за низкого электрического сопротивления морской воды можно сравнить с каротажным зондированием градиент-зондами AMN в условиях заполнения скважины соленым буровым раствором. В заполненных соленым раствором скважинах промысловые геофизики отказались от такого зондирования более пятидесяти лет тому назад, а удельное электрическое сопротивление пластов горных пород определяют при помощи основанного на принципе радиальной фокусировки тока бокового каротажа (Laterolog. Doll H.G.).
Задача морской электроразведки состоит в том, чтобы определить наличие углеводородов в традиционно встречающихся сейсмических структурах, которые при наличии в них УВ отличаются по электрическому сопротивлению от окружающих их пород не более чем в 20-30 раз. Таких структур в сотни раз больше, чем структур типа TWGP. Например, все месторождения УВ в Западной Сибири, в Обской губе и на шельфе Карского моря, в том числе и крупнейшие по своим геометрическим размерам, отличаются по электрическому сопротивлению от вмещающих пород по каротажным данным всего лишь в пределах от трех до тридцати раз. И если после проведения сейсморазведочных исследований такие структуры, как TWGP в Норвежском море, не нуждаются в дополнительной проверке их на наличие углеводородов, то из-за высокой стоимости морского бурения большинство структур без проверки на наличие в них углеводородов разбуривать рискованно. Так, например, Established UK oil and gas consultancy Hannon Westwood (см. «Upstream boom likely to frustrate North Sea investment opportunities». First Break, 25 Januar 2007, стр.22-24) не рекомендует вкладывать инвестиции в бурение оставшихся неразбуренными средних и малых сейсмических структур в Северном море, так как успешность бурения в настоящее время на этих структурах составляет всего лишь около 20% (одна удачная скважина из пяти пробуренных).
Неблагоприятная для сейсморазведки ситуация сложилась и на Сахалинском шельфе. Там на выявленных сейсморазведкой крупных структурах Западно-Шмидтовского блока («Сахалин-4») компанией British Petroleum пробурено две разведочные скважины по данным прессы стоимостью в $103 миллиона, в которых углеводородов не обнаружено (см. газета «Ведомости» от 06.03.2008). Также оказалась сухой скважина, пробуренная на одной из крупнейших по геометрическим размерам Адмиралтейской структуре в Баренцевом море (см. Обметко В.В. и др. Перспективы нефтегазоносности Адмиралтейского мегавала. Международная конференция «Нефть и газ Арктического шельфа - 2008». Мурманск. 12-14 ноября 2008 г.). Это связано с тем, что скопления углеводородов не всегда контролируются только структурным фактором.
Известен способ морской геоэлектроразведки с фокусировкой электрического тока, который благодаря поддержанию равенства нулю первой осевой или ортогональной разности электрических потенциалов в находящейся в морском слое измерительной точке зондирования и определению нескольких электрофизических параметров элементов геологической среды мало подвержен влиянию боковых неоднородностей и с вполне достаточной степенью вероятности позволяет определять наличие углеводородов в выявленных сейсморазведочными исследованиями структурах и их оконтуривать (Рыхлинский Н.И., Давыдычева С.Н., Лисин А.С. "Способ морской геоэлектроразведки с фокусировкой электрического тока". Патент РФ №2884555 от 01 июня 2005 г., Бюл. №27, 2006) [3].
Недостатком этого способа является то, что им осуществляется неполная вертикальная фокусировка тока: в первом его варианте исключается в точке зондирования только ортогональная горизонтальная составляющая плотности тока, а во втором - только осевая горизонтальная составляющая плотности тока.
Известен способ морской геоэлектроразведки с круговой фокусировкой электрического тока (Рыхлинская Е.Н. "Способ морской геоэлектроразведки с фокусировкой электрического тока (варианты)". Патент РФ №2351958 от 12 ноября 2007 г., Бюл. №10, 2009) [4]. Этот способ хотя и исключает горизонтальную компоненту плотности электрического тока jxy, но он предназначен для решения других не менее важных задач. Он в силу специальной конструкции измерительных датчиков позволяет измерять только изменение электрических параметров геологического разреза вдоль измерительного профиля, что важно для точного определения границ контура залежей УВ. Тем не менее, в нем не заявлено, что токовый электрический диполь может размещаться на поверхности моря. Кроме этого для его осуществления требуется относительно сложная электрическая схема измерителя.
В предложенном способе при сравнительно простой электрической схеме измерителя решается задача обнаружения и оконтуривания с четко выраженной границей контура нефтегазовых залежей на континентальном шельфе как на больших и малых глубинах моря, так и на больших глубинах залегания этих залежей.
Технический результат, позволяющий решать данную задачу, заключается в возможности полного исключения в точке зондирования горизонтальной составляющей плотности тока jxy как ортогональной jy, так и осевой jx, что полностью исключает боковое влияние и по этой причине дает возможность при больших глубинах моря, по крайней мере до четырех километров, находить и с высокой контрастностью оконтуривать скопления углеводородов в толще геологической среды.
Этот технический результат достигается тем, что в способе морской геоэлектроразведки, при котором возбуждают электромагнитное поле в толще исследуемой среды, пропуская через нее периодические прямоугольные импульсы тока с паузами после каждого из них при помощи горизонтальных дипольных электрических источников; в каждой точке зондирования в период импульса тока и на протяжении каждой паузы после выключения тока измеряют с постоянным интервалом времени Δt последовательность мгновенных значений первых разностей электрических потенциалов прямого поля и поля переходных процессов, при этом обеспечивают условие равенства нулю результирующих первых разностей электрических потенциалов между внешними измерительными электродами; формируют интерпретируемые параметры и, используя их и дифференциальное уравнение математической физики для напряженности электрического поля дипольного источника в электрохимически поляризующейся проводящей среде
где
∇× - оператор rot;
- напряженность электрического поля дипольного источника, выраженная в уравнении для случая гармонического изменения величины электрического поля по времени, V/m;
µ - магнитная проницаемость - величина, постоянная для немагнитных сред, к числу которых относятся осадочные геологические породы, и равна 4π·10-7 Гн/м;
σ(iωσoητ) - частотно-зависимая электропроводность элементов среды,
σo - электропроводность элементов среды без учета влияния вызванной поляризации, S/m,
η - коэффициент их вызванной поляризации,
τ - постоянная времени спада разности потенциалов вызванной поляризации;
решают математическую обратную задачу, определяя присущие каждому элементу среды три электрофизических параметра: удельную электропроводность σо, вызванную поляризацию η и постоянную времени спада разности потенциалов вызванной поляризации τ;
и строят три временных разреза по этим параметрам;
согласно изобретению прокладывают на дне моря измерительный профиль, который проходит через зафиксированные на морском дне измерительные установки, каждая из которых состоит из пяти электродов: центрального и расположенных равноудаленно вокруг него четырех по вершинам квадрата, две противоположные стороны которого параллельны оси профиля;
в пределах каждого периода «импульс-пауза» проводят геометрическое зондирование при включенном токе и зондирование на переходных процессах на протяжении паузы после выключения тока, измеряя первые разности электрических потенциалов между каждым из четырех внешних электродов измерительной установки и центральным;
при этом измерения в каждой фиксированной на измерительном профиле измерительной установке осуществляют при прохождении горизонтального дипольного источника по токовому профилю, проложенному параллельно измерительному и сдвинутому относительно него в плане в ортогональном направлении по оси у на расстояние (у=-b) и по высоте на толщину h водного слоя, посылая токовые импульсы в исследуемую среду при всех положениях дипольного источника от точки с координатами [(x=-L), (y=-b), (z=+h)] до точки с координатами [(x=+L), (y=-b), (z=+h)], затем дипольный источник разворачивают и переводят на другой параллельный токовый профиль, проложенный и сдвинутый в плане относительно измерительного в противоположную сторону по оси у на расстояние (у=+b) и по высоте на толщину h, и продолжают измерение при его движении в обратном направлении от точки с координатами [(x=+L), (y=+b), (z=+h)] до точки с координатами [(x=-L), (y=+b), (z=+h)],
используя измеренные значения первых разностей электрических потенциалов, определяют два независимых от силы тока источника и горизонтальной компоненты плотности тока (jx=0 и jy=о) в точке зондирования интерпретируемых параметра - один Rz(to) на основе геометрического зондирования при всех положениях дипольного источника, вычисляемый по формуле:
Rz(to)=[ΔUM1N(IB1A1, to)-ΔUM3N(IB1A1, to)]:
:{[ΔUM1N(IB1A1, to)+ΔUM2N(IB1A1, to)+ΔUM3N(IB1A1, to)+ΔUM4N(IB1A1, to)]+k1(to)·[ΔUM1N(IA2B2, to)+ΔUM2N(IA2B2, to)+ΔUM3N(IA2B2, to)+ΔUM4N(IA2B2, to)]+k2(to)·[ΔUM1N(IA3B3, to)+ΔUM2N(IA3B3, to)+ΔUM3N(IA3B3, to)+ΔUM4N(IA3B3, to)]+k3(to)·[ΔUM1N(IB4A4, to)+ΔUM2N(IB4A4, to)+ΔUM3N(IB4A4, to)+ΔUM4N(IB4A4, to)]},
и другой Rz(ti) на основе зондирования на переходных процессах при четырех, выбранных методом итераций, наиболее информативных разносах с координатами дипольного источника [(х=-a), (y=-b), (z=+h)], [(х=+а), (у=-b), (z=+h)], [(x=+a), (y=+b), (z=+h)] и [(х=-а), (у=+b), (z=+h)] из всех прозондированных, вычисляемый по формуле:
Rz(ti)=[ΔUM1N(IB1A1, ti)-ΔUM3N(IB1A1, ti)]:
:{[ΔUM1N(IB1A1, ti)+ΔUM2N(IB1A1, ti)+ΔUM3N(IB1A1, ti)+ΔUM4N(IB1A1, ti)]+k1(ti)·[ΔUM1N(IA2B2, ti)+ΔUM2N(IA2B2, ti)+ΔUM3N(IA2B2, ti)+ΔUM4N(IA2B2, ti)]+k2(ti)·[ΔUM1N(IA3B3, ti)+ΔUM2N(IA3B3, ti)+ΔUM3N(IA3B3, ti)+ΔUM4N(IA3B3, ti)]+k3(ti)·[ΔUM1N(IB4A4, ti)+ΔUM2N(IB4A4, ti)+ΔUM3N(IB4A4, ti)+ΔUM4N(IB4A4, ti)]},
где
k1(to), k2(to), k3(to) - коэффициенты фокусировки при геометрическом зондировании, обеспечивающие эквипотенциальность замкнутой линии, проходящей через четыре внешних электрода измерительной установки, и, тем самым, исключение внутри области, ограниченной контуром этой линии, горизонтальной составляющей плотности тока jxy в каждой точке зондирования на всех геометрических разносах в период импульса тока, которые определяют из системы трех уравнений:
,
,
,
k1(ti), k2(ti), k3(ti) - коэффициенты фокусировки при зондировании на переходных процессах в каждой точке зондирования в паузе тока на всех временах переходных процессов, определяемые из системы трех уравнений и обеспечивающие эквипотенциальность замкнутой линии, проходящей через четыре внешних электрода измерительной установки:
,
,
,
t0 - момент времени при пропускании токового импульса, когда электрическое поле переходных процессов не отличается от своего установившегося значения, соответствующего постоянному току;
ti - моменты времени, при которых измеряют сигналы переходных процессов через равные интервалы времени Δt на протяжении всей паузы после выключения тока;
ΔUM1N(IB1A1, to), ΔUM2N(IB1A1, to), ΔUM3N(IB1A1, to), ΔUM4N(IB1A1, to),
ΔUM1N(IA2B2, to), ΔUM2N(IA2B2, to), ΔUM3N(IA2B2, to), ΔUM4N(IA2B2, to),
ΔUM1N(IA3B3, to), ΔUM2N(IA3B3, to), ΔUM3N(IA3B3, to), ΔUM4N(IA3B3, to),
ΔUM1N(IB4A4, to), ΔUM2N(IB4A4, to), ΔUM3N(IB4A4, to), ΔUM4N(IB4A4, to)
- мгновенные значения первых разностей электрических потенциалов между каждым из четырех внешних электродов М1М2М3М4 измерительной установки и центральным N, измеренные при времени tо пропускания тока в дипольном источнике при его прохождении по обоим токовым профилям через равноудаленные от центрального электрода N измерительной установки четыре точки соответственно с координатами [-х, (y=-b), (z=+h)], [+x, (y=-b), (z=+h)], [+х, (у=+b), (z=+h) и [-х, (у=+b), (z=+h)];
ΔUM1N(IB1A1, ti), ΔUM2N(IB1A1, ti), ΔUM3N(IB1A1, ti), ΔUM4N(IB1A1, ti),
ΔUM1N(IA2B2, ti), ΔUM2N(IA2B2, ti), ΔUM3N(IA2B2, ti), ΔUM4N(IA2B2, ti),
ΔUM1N(IA3B3, ti), ΔUM2N(IA3B3, ti), ΔUM3N(IA3B3, ti), ΔUM4N(IA3B3, ti),
ΔUM1N(IB4A4, ti), ΔUM2N(IB4A4, ti), ΔUM3N(IB4A4, ti), ΔUM4N(IB4A4, ti)
- мгновенные значения первых разностей электрических потенциалов между каждым из четырех внешних электродов М1М2М3М4 измерительной установки и центральным N, измеренные при времени ti переходных процессов в паузе тока в дипольном источнике при его прохождении по обоим токовым профилям через равноудаленные от центрального электрода N измерительной установки четыре точки соответственно с координатами [(x=-a), (y=-b), (z=+h)], [(x=+a), (y=-b), (z=+h)], [(х=+а), (y=+b), (z=+h)] и [(x=-a), (y=+b), (z=+h)].
Также технический результат достигается тем, что в способе морской геоэлектроразведки согласно изобретению токовые профили для прохождения горизонтального дипольного источника прокладывают на заданной глубине моря или на дневной поверхности моря.
Также технический результат достигается тем, что в способе морской геоэлектроразведки согласно изобретению расстояние L равно шести километрам и более, расстояния а и b равны одному километру и более, глубина погружения h дипольного источника устанавливается в пределах от двадцати метров выше уровня морского дна до дневной поверхности.
Также технический результат достигается тем, что в способе морской геоэлектроразведки согласно изобретению измерительные установки размещают на измерительном профиле с равным шагом в 200÷1000 метров.
Сущность изобретения поясняется чертежами.
На фиг.1 дана блок-схема устройства для реализации предложенного способа.
На фиг.2 дана схема расположения в плане на морском дне группы измерительных установок на измерительном профиле зондирования и два профиля, параллельных измерительному, вдоль которых перемещают дипольный источник в толще водного слоя моря выше морского дна на расстояние h или по дневной поверхности моря.
На фиг.3 дана схема движения дипольного источника по первому из двух параллельных профилей относительно измерительного.
На фиг.4 дана схема движения дипольного источника по второму параллельному профилю относительно измерительного.
На фиг.5 показана форма электрического поля с фокусировкой электрического тока в направлении вертикальной координаты Z путем поддержания равенства нулю первых разностей электрических потенциалов между каждой парой внешних измерительных электродов от результирующего действия четырех равноудаленных от центрального электрода N измерительной установки дипольных источников. При равенстве нулю разностей электрических потенциалов между четырьмя внешними электродами M1M2M3M4 токи от четырех токовых источников АВ, проникшие внутрь образованного этими электродами квадрата, фокусируются и в дальнейшем могут распространяться только в вертикальном направлении вдоль координаты Z. Кроме этого вертикальное направление электрического тока при условии его фокусировки поддерживается на основе изопериметрической вариации согласно функции Лагранжа (принцип наименьшего действия) (Ландау Л.В., Лифшиц Е.М. Теория поля. 2 изд. М. - Л. 1948). В плоскости M1M2M3M4 N потенциал электрического поля приобретает форму эллипсоида вращения с минимумом в точке N.
На фиг.6 показаны формы одиночных импульсов в функции времени t: а) - форма одиночного прямоугольного импульса тока I в цепи дипольного источника АВ; б) - форма импульсов первой и второй разностей электрических потенциалов.
Устройство для осуществления способа (фиг.1) содержит измерительную установку для измерения первых разностей электрических потенциалов. Эта установка расположена на измерительном профиле - 1, проложенном по морскому дну.
Здесь 2 - N, 3 - M1, 4 - M2, 5 - М3, 6 - М4 - электроды измерительной установки: внешние M1M2M3M4 и центральный N (таких измерительных пятиэлектродных датчиков на измерительном профиле 1 может быть расставлено любое количество (фиг.2)); 7 - измеритель для измерения первой разности потенциалов электрического поля ΔUM1N между внешним электродом M1 измерительной установки и центральным N; 8 - измеритель для измерения первой разности потенциалов электрического поля ΔUM2N между внешним электродом M2 измерительной установки и центральным N; 9 - измеритель для измерения первой разности потенциалов электрического поля ΔUM3N между внешним электродом М3 измерительной установки и центральным N; 10 - измеритель для измерения первой разности потенциалов электрического поля ΔUM4N между внешним электродом М4 измерительной установки и центральным N.
Для обеспечения фокусировки электрического поля в точке зондирования N параллельно измерительному профилю 1 справа и слева от него размещают в толще водного слоя моря выше морского дна на расстояние h или на поверхности моря два профиля: первый 11 и второй 12, вдоль которых перемещают горизонтальный токовый дипольный источник, состоящий из двух электродов 13-А и 14-В, питаемых генератором 15 прямоугольных импульсов тока с паузами между каждым из импульсов.
На фиг.3 показано продвижение дипольного источника АВ по первому профилю 11, расположенному в толще водного слоя моря выше морского дна на расстоянии h или на поверхности моря параллельно измерительному 1 и сдвинутому от него в плане на расстояние (у=-b). Продвижение дипольного источника АВ по профилю 11 осуществляют от его начала [-х, (y=-b),(z=+h)] до конца [+х, (y=-b),
(z=+h)]. Затем дипольный источник АВ (фиг.4) перебрасывают на второй профиль 12, параллельный измерительному и сдвинутый относительно него в плане на расстояние (у=+b) и по высоте на толщину h. Продвижение дипольного источника АВ по профилю 12 осуществляют от точки с координатами [+х, (y=+b), (z=+h)] до точки с координатами [-х, (y=+b), (z=+h)].
Такое сложное профилирование необходимо для измерения четырех первых разностей электрических потенциалов между каждым из четырех внешних электродов измерительной установки M1 M2 М3 М4 и центральным N для достижения равенства потенциалов между внешними электродами M1 M2 М3 М4 от суммарного действия дипольного источника при его прохождении через четыре точки с координатами [-х, (у=-b), (z=+h)], [+x, (y=-b), (z=+h)], [+x, (у=+b, (z=+h)] и [-x, (y=+b), (z=+h)]. Это видно, например, из измеренных четырех первых разностей электрических потенциалов: ΔUM1N-ΔUM2N=ΔUM1M2=UM1-UM2=0, ΔUM2N-ΔUM3N=ΔUM2M3=UM2-UM3=0, ΔUM3N-ΔUM4N=ΔUM3M4=UM3-UM4=0. Откуда следует, что UM1=UM2=UM3=UM4.
Равенство потенциалов между четырьмя внешними электродами измерительной установки обеспечивает эквипотенциальность замкнутой линии, проходящей через эти четыре внешних электрода, и, тем самым, исключение внутри контура этой линии горизонтальной компоненты плотности тока jxy в каждой точке зондирования, что обеспечивает фокусировку тока вдоль вертикальной координаты Z внутри этой замкнутой эквипотенциальной линии.
На фиг.6 (а) показана форма одиночного прямоугольного импульса тока I в цепи дипольного источника АВ в функции времени t. Здесь Т - период (импульс тока - пауза).
На фиг.6 (б) показана форма одного из импульсов ΔnU, где n=1 и 2. Также показано одно из значений ΔnU(t0) в период пропускания тока в дипольный источник АВ и одно из значений ΔnU(ti) в паузе тока.
Рассмотрим основы предложенного способа, его осуществление и новые возможности морской геоэлектроразведки.
Предложенный способ морской геоэлектроразведки позволяет при помощи трех профилей (фиг.1) исключать боковое влияние путем фокусировки тока. Это дает большую глубинность зондирования и высокое разрешение по латерали.
Известно, что электромагнитное поле в плохопроводящей физической среде распространяется по времени t согласно вытекающему из первого и второго уравнений Максвелла дифференциальному затухающему волновому уравнению математической физики для напряженности электрического поля, в том числе и в случае ее импульсного изменения:
где
∇× - оператор rot;
Е - напряженность электрического поля, В/м;
µ - магнитная проницаемость - величина, постоянная для немагнитных сред, к числу которых относятся осадочные геологические породы, она равна 4π·10-7 Гн/м;
σо - электропроводность неполяризующейся среды, См/м;
ε - диэлектрическая проницаемость, Ф/м. (В.А.Говорков. Электрические и магнитные поля. М.: Госэнергоиздат, 1960, с.257-263) [5].
В случае высокопроводящей среды, к которой относятся осадочные отложения горных пород, в связи с тем что σо численно многократно больше ε, второй член в правой части уравнения (2) мал по сравнению с первым, и его отбрасывают (Л.Л.Ваньян. Основы электромагнитных зондирований. М.: Недра, 1965, с.28-30) [6]. Физически это означает, что токами смещения в проводящих средах пренебрегают ввиду их малости по сравнению с токами проводимости. Тогда уравнение (2) принимает вид
Это уравнение в геоэлектроразведке имеет аналитическое решение лишь для одномерных осесимметричных сред, в частности для сред с горизонтальными неограниченно простирающимися границами раздела. При этом отметим, что в реальности геологическая среда всегда трехмерно-неоднородна, так как в ней, во-первых, присутствуют приповерхностные локальные неоднородности, во-вторых, в целом геологическая среда вдоль профиля исследований постоянно меняет свои электрофизические параметры. Однако уравнение (3), как было сказано выше, аналитически решено только для одномерных осесимметричных сред, в том числе для среды с горизонтальными плоскопараллельными границами раздела.
Поэтому использование решения уравнения (3) в обратных задачах геоэлектроразведки для поисков и оконтуривания нефтегазовых залежей допустимо лишь в том случае, когда при полевых измерениях осуществляется фокусировка электрического тока дипольного источника электромагнитного поля, так как в этом случае в точке N (фиг.5) форма распространения поля от результирующего действия четырех дипольных источников B1A1, A2B2, A3B3, B4A4 при поддержании нулю результирующих первых разностей электрических потенциалов между каждой парой внешних измерительных электродов практически всегда одинаковая как в трехмерно-неоднородной среде, так и в одномерной с плоскопараллельными горизонтальными границами раздела.
Из данных ниже формул (4), (5), (6) и (7) вытекает, что согласно закону Ома суммарная составляющая горизонтальной плотности тока jxy в точке зондирования N согласно принципу суперпозиции равна нулю.
Таким образом, в точке измерения происходит фокусировка тока, приводящая к исключению горизонтальной составляющей плотности тока jxy и к сохранению только вертикальной компоненты плотности тока jz, что позволяет исключать искажающее боковое влияние на результаты зондирования и позволяет корректно решать обратную задачу в трехмерно-неоднородных средах, используя известное аналитическое решение одномерного уравнения (14) для слоистых сред с плоскопараллельными границами раздела.
Чтобы обеспечить исключение горизонтальной компоненты плотности jxy в точке зондирования N с координатами [(х=0), (у=0), (z=0)], в предложенном способе выведены соответствующие формулы измеряемых электрических параметров: один Rz(to) на основе геометрического зондировании при всех положениях дипольного источника, вычисляемый по формуле:
и другой Rz(ti) на основе зондирования на переходных процессах при четырех, выбранных методом итераций, наиболее информативных разносах с координатами дипольного источника [(х=-а), (y=-b), (z=+h)], [(х=+а), (y=-b), (z=+h)], [(х=+а), (y=+b), (z=+h)] и [(х=-а), (у=+b), (z=+h)] из всех прозондированных, вычисляемый по формуле:
где
k1(to), k2(to), k3(to) - коэффициенты фокусировки при геометрическом зондировании, обеспечивающие эквипотенциальность замкнутой линии, проходящей через четыре внешних электрода измерительной установки, и, тем самым, исключение внутри области, ограниченной контуром этой линии, горизонтальной составляющей плотности тока jxy в каждой точке зондирования на всех геометрических разносах в период импульса тока, которые определяют из системы трех уравнений:
,
,
k1(ti), k2(ti), k3(ti) - коэффициенты фокусировки при зондировании на переходных процессах в каждой точке зондирования в паузе тока на всех временах переходных процессов, определяемые из системы трех уравнений и обеспечивающие эквипотенциальность замкнутой линии, проходящей через четыре внешних электрода измерительной установки:
,
,
t0 - момент времени при пропускании токового импульса, когда электрическое поле переходных процессов не отличается от своего установившегося значения, соответствующего постоянному току;
ti - моменты времени, при которых измеряют сигналы переходных процессов через равные интервалы времени Δt на протяжении всей паузы после выключения тока;
ΔUM1N(IB1A1, to), ΔUM2N(IB1A1, to), ΔUM3N(IB1A1, to), ΔUM4N(IB1A1, to),
ΔUM1N(IA2B2, to), ΔUM2N(IA2B2, to), ΔUM3N(IA2B2, to), ΔUM4N(IA2B2, to),
ΔUM1N(IA3B3, to), ΔUM2N(IA3B3, to), ΔUM3N(IA3B3, to), ΔUM4N(IA3B3, to),
ΔUM1N(IB4A4, to), ΔUM2N(IB4A4, to), ΔUM3N(IB4A4, to), ΔUM4N(IB4A4, to)
- мгновенные значения первых разностей электрических потенциалов между каждым из четырех внешних электродов М1М2М3М4 измерительной установки и центральным N, измеренные при времени t0 пропускания тока в дипольном источнике при его прохождении по обоим токовым профилям через равноудаленные от центрального электрода N измерительной установки четыре точки соответственно с координатами [-х, (y=-b), (z=+h)], [+x, (y=-b), (z=+h)], [+x, (y=+b), (z=+h)] и [-х, (у=+b), (z=+h)];
ΔUM1N(IB1A1, ti), ΔUM2N(IB1A1, ti), ΔUM3N(IB1A1, ti), ΔUM4N(IB1A1, ti),
ΔUM1N(IA2B2, ti), ΔUM2N(IA2B2, ti), ΔUM3N(IA2B2, ti), ΔUM4N(IA2B2, ti),
ΔUM1N(IA3B3, ti), ΔUM2N(IA3B3, ti), ΔUM3N(IA3B3, ti), ΔUM4N(IA3B3, ti),
ΔUM1N(IB4A4, ti), ΔUM2N(IB4A4, ti), ΔUM3N(IB4A4, ti), ΔUM4N(IB4A4, ti),
- мгновенные значения первых разностей электрических потенциалов между каждым из четырех внешних электродов М1М2М3М4 измерительной установки и центральным N, измеренные при времени ti переходных процессов в паузе тока в дипольном источнике при его прохождении по обоим токовым профилям через равноудаленные от центрального электрода N измерительной установки четыре точки соответственно с координатами [(x=-a), (y=-b), (z=+h)], [(x=+a), (y=-b), (z=+h)], [(x=+a), (у=+b), (z=+h)] и [(х=-а), (у=+b), (z=+h)].
Формулы (4) и (5) независимо от величины тока в диполе АВ в любой трехмерно-неоднородной среде в каждой точке пространства в направлении вертикальной координаты z в точке зондирования N в плоскости, проходящей через измерительные электроды на любых расстояниях х между дипольным источником АВ и точкой зондирования N при геометрическом зондировании и на всем протяжении времени ti при зондировании на переходных процессах, обеспечивают равенство нулю горизонтальной компоненты плотности тока jxy. Это происходит всегда независимо от того, что в процессе изменения расстояния х между дипольным источником АВ и точкой измерения N или в процессе изменения времени t переходных процессов множители k1(to), k2(to), k3(to) и k1(ti), k2(ti), k3(ti) меняются. Благодаря этому при решении обратной задачи исключается эффект бокового влияния, т.е. электрическое поле в трехмерно-неоднородной среде в точке зондирования N, описанное формулами (4) и (5), всегда практически совпадет с полем в одномерной горизонтально-слоистой среде с неограниченными границами раздела.
Это позволяет решать обратную задачу в точке измерения N для трехмерно-неоднородной среды, используя известное аналитическое решение для одномерной среды с горизонтально-слоистыми границами раздела.
Отметим также, что уравнение (3) является уравнением распространения по времени электромагнитного поля в проводящей неполяризующейся среде, которое совпадает с известным в математической физике уравнением теплопроводности или диффузии и которое в геофизике в способах сопротивлений обычно используют для изучения распространения переменного электромагнитного поля вглубь толщи исследуемых геологических пород.
При этом считают, что электропроводность σо того или иного геологического горизонта является главным и практически единственным определяющим его электрические свойства параметром, обладает своим постоянным значением для каждого горизонта и не зависит от частоты возбуждения электромагнитного поля. Однако геологическим осадочным породам при их возбуждении применяемым в геофизике переменным низкочастотным электрическим током свойственна вызванная им поляризация η. Вызванная поляризация есть безразмерная величина, зависящая от электрохимической активности осадочных горных пород. Она определяется как отношение разностей потенциалов, измеренных на образце исследуемой породы после выключения импульсов тока через 0,5 сек (ΔUВП) и до выключения (ΔU). Это отношение обычно выражают в процентах:
Вызванная поляризация осадочных геологических пород имеет уникальную среди физических параметров стабильность и практически не зависит от состава пород и их температуры. Она для ионопроводящих (осадочных) пород зависит от многих факторов: влажности и пористости, состава и концентрации раствора в порах породы, структуры и размера пор, содержания глинистых минералов и т.д. (В.А.Комаров. Электроразведка методом поляризации. Л.: Наука, 1980, с.392) [7]. И, самое главное, вызванная поляризация несет основную информацию о присутствии в геологической среде обладающих высокой степенью этой поляризации нефтегазовых залежей.
Установлено (W.H.Pelton, S.H.Ward, P.G.Hallof, W.R.Sill and P.H.Nelson. Mineral discrimination and removal of inductive coupling with multi-frequency JP, Geophysics 43, 1978, c.588-603) [8], что электропроводность осадочных горных пород не постоянна, а зависит от вызванной поляризации и от частоты возбуждения электрического поля по предложенной, в частности, K.S.Cole и R.H.Cole в форме гармонического его изменения по времени эмпирической формуле:
в которой эта электропроводность зависит от ω, σо, η и τ,
где:
η - вызванная поляризация пород, безразмерная величина, обычно выражающаяся в процентах;
τ - постоянная времени, определяющая скорость спада разности потенциалов, связанной с вызванной поляризацией, сек;
ω - гармоническая частота электрического возбуждения, Гц;
с - безразмерный показатель степени, который хотя и не является физическим параметром горных пород, но и от него зависит σ(iωσоητ).
Вызванная поляризация η на низких частотах электрического возбуждения, в отличие от диэлектрической проницаемости ε, численно не столь мала по сравнению с электропроводностью σo для осадочных геологических пород, измеренной, например, при токах высоких частот (ω→∞), когда, как это видно из формулы (9), вызванная поляризация не проявляется. Следовательно, вызванной поляризацией при изучении с целью поисков и оконтуривания нефтегазовых залежей геоэлектрических параметров осадочных геологических пород на низкочастотном переменном токе уже пренебрегать нельзя.
Известно (Электроразведка. Справочник геофизика. Ред. В.К.Хмелевский и др. М.: Недра, 1989, книга вторая, с.99-102) [9], что для тех или иных осадочных геологических пород через 0,5 сек после выключения импульса возбуждающего тока величина вызванной поляризацией разности потенциалов несмотря на ее интенсивный спад еще сохраняет уровни, численные значения которых составляют от 0,2% до 10% от численных значений разностей потенциалов прямого поля, связанных с электропроводностью σo, измеренной, как отмечалось выше, при токах высоких частот, когда вызванная поляризация не проявляется. Чтобы сохранить по форме формулу (9), тепловое уравнение (3) запишем для случая гармонического изменения величины электромагнитного поля по времени, имея ввиду, что:
и учитывая то, что:
и
Тогда уравнение (3) для проводящей неполяризующейся среды с учетом преобразования (10) примет вид:
Но поскольку электропроводность осадочных горных пород не постоянна, а зависит от вызванной поляризации и от частоты возбуждения по формуле (9), то уравнение (12) с учетом этой формулы приобретает уже три определяющих свойства поляризующейся среды параметра -σо, η, τ вместо одного σо и для случая гармоничного изменения величины электромагнитного поля по времени принимает вид:
а в общем виде с учетом (9) -
Замена частотно-независимой электропроводности σо, которая присутствует в уравнении (12), на частотно-зависимую σ(iω), присутствующую в уравнении (14), математически корректна и теоретически доказана (А.К.Куликов, Е.А.Шемякин. Электроразведка фазовым методом вызванной поляризации. М.: Недра, 1978 г., стр.24-26) [10] и (Дж.Р.Уайт. Гео-электро-магнетизм. М.: Недра, 1987 г., стр.61-62) [11].
Для предложенного способа задача обнаружения нефтегазовых залежей в исследуемой толще горных пород как математическая обратная задача решается согласно уравнению (14) в функции расстояния х между дипольным источником АВ и точкой измерения N, в функции времени t переходных процессов и, как следствие этого, в функции, зависящей от расстояния х и времени t глубины проникновения электромагнитного поля, по трем независимым друг от друга параметрам среды: электропроводности σo, вызванной поляризации η и постоянной времени τ спада разности электрических потенциалов вызванной поляризации.
Эта задача как обратная математическая задача решается для предложенного способа путем использования всего массива определяемых этим способом, по крайней мере, двух независимых от силы тока источников нормированных электрических параметров, описанных формулами (4) (5), (6) и (7).
Следует отметить, что при зондировании предложенным способом в каждой измерительной точке профиля на всем его протяжении получают цифровую информацию по первым разностям электрических потенциалов с шагом дискретности через равные интервалы времени Δt как при каждом включении тока в дипольном источнике в процессе зондирования (геометрическое зондирование), так и в паузе тока (зондирование на переходных процессах).
При этом, если для геометрического зондирования получение информации об отклике среды от каждого определяющего разнос зондирующей установки импульса тока необходимо, то для зондирования на переходных процессах достаточно информации от одной зондирующей установки, но с таким разносом, который несет наибольшую информацию о зондируемой геологической среде (отметим, что при больших разносах сигналы переходных процессов слабые и их трудно выделить на фоне помех, а при малых разносах зондирование на переходных процессах обладает малой глубинностью). Поэтому в решении обратной задачи, в частности, поставлена также задача найти методом итераций наиболее информативный разнос зондирующей установки среди всех прозондированных.
В конечном результате решения обратной задачи с учетом итераций находят модель среды, наиболее близкую по геометрическому строению и электрическим параметрам к исследуемой, строят геоэлектрические разрезы полученных параметров σо, η и τ и выделяют на них участки с аномальными значениями, соответствующими в плане положению нефтегазовых залежей.
Пример конкретного выполнения
На фиг.1 представлена блок-схема устройства для реализации предложенного способа. На блок-схеме показан помещенный в толщу морской воды на высоту h относительно морского дна или на дневную поверхность дипольный источник АВ (13 и 14), питаемый генератором 15 прямоугольных импульсов тока. Отметим, что при помещении дипольного источника на дневную поверхность упрощается конструкция кабельной косы дипольного источника и создается возможность повысить величину электрического тока в цепи дипольного источника. Помещение дипольного источника на дневную поверхность позволяет вести исследования предлагаемым способом при малых глубинах моря (до 300 м). При больших глубинах моря (свыше 300 м) из-за поглощения сигналов переходных процессов в толще высокопроводящей морской воды горизонтальный дипольный источник погружают на заданную глубину и его прохождение вдоль двух токовых профилей осуществляют на этой глубине.
Устройство (фиг.1) содержит измерительную установку для измерения разностей электрических потенциалов. Эта установка расположена на измерительном профиле 1, проложенном по морскому дну. Здесь 2 - N, 3 - M1, 4 - M2, 5 - М3, 6 - М4 - электроды измерительной установки.
Устройство осуществляет измерение четырех первых разностей электрических потенциалов, где: 7 - измеритель для измерения первой разности потенциалов электрического поля ΔUM1N между внешним электродом M1 измерительной установки и центральным N; 8 - измеритель для измерения первой разности потенциалов электрического поля ΔUM2N между внешним электродом M2 измерительной установки и центральным N; 9 - измеритель для измерения первой разности потенциалов электрического поля ΔUM3N между внешним электродом М3 измерительной установки и центральным N; 10 - измеритель для измерения первой разности потенциалов электрического поля ΔUM4N между внешним электродом М4 измерительной установки и центральным N.
Измеренные первые разности потенциалов электрического поля обрабатывают по формулам (4), (5), (6) и (7) при помощи компьютера для получения численных значений Rz(to) и Rz(ti), используемых для решения обратной задачи.
Отметим, что по данным математического моделирования расстояние L равно шести километрам и более, расстояния а и b равны одному километру и более, глубина погружения h дипольного источника устанавливается в пределах от двадцати метров выше уровня морского дна до дневной поверхности.
Следует также отметить, что этот способ может быть использован и для исследований на суше, где измерительные датчики и токовые электрические диполи размещают на поверхности земли.
1. Способ морской геоэлектроразведки, при котором возбуждают электромагнитное поле в толще исследуемой среды, пропуская через нее периодические прямоугольные импульсы тока с паузами после каждого из них при помощи горизонтальных дипольных электрических источников; в каждой точке зондирования в период импульса тока и на протяжении каждой паузы после выключения тока измеряют с постоянным интервалом времени At последовательность мгновенных значений первых разностей электрических потенциалов прямого поля и поля переходных процессов, при этом обеспечивают условие равенства нулю результирующих первых разностей электрических потенциалов между внешними измерительными электродами;
формируют интерпретируемые параметры и, используя их и дифференциальное уравнение математической физики для напряженности электрического поля дипольного источника в электрохимически поляризующейся проводящей среде
,
где ∇x - оператор rot;
- напряженность электрического поля дипольного источника, выраженная в уравнении для случая гармонического изменения величины электрического поля по времени, V/m;
µ - магнитная проницаемость - величина постоянная для немагнитных сред, к числу которых относятся осадочные геологические породы, и равна 4π·10-7 Гн/м;
σ(iωσ0ητ) - частотно-зависимая электропроводность элементов среды;
σ0 - электропроводность элементов среды без учета влияния вызванной поляризации, S/m;
η - коэффициент их вызванной поляризации;
τ - постоянная времени спада разности потенциалов вызванной поляризации;
решают математическую обратную задачу, определяя присущие каждому элементу среды три электрофизических параметра: удельную электропроводность σо, вызванную поляризацию η и постоянную времени спада разности потенциалов вызванной поляризации τ;
и строят три временных разреза по этим параметрам;
отличающийся тем, что прокладывают на дне моря измерительный профиль, который проходит через зафиксированные на морском дне измерительные установки, каждая из которых состоит из пяти электродов: центрального и расположенных равноудаленно вокруг него четырех по вершинам квадрата, две противоположные стороны которого параллельны оси профиля;
в пределах каждого периода «импульс-пауза» проводят геометрическое зондирование при включенном токе и зондирование на переходных процессах на протяжении паузы после выключения тока, измеряя первые разности электрических потенциалов между каждым из четырех внешних электродов измерительной установки и центральным;
при этом измерения в каждой фиксированной на измерительном профиле измерительной установке осуществляют при прохождении горизонтального дипольного источника по токовому профилю, проложенному параллельно измерительному и сдвинутому относительно него в плане в ортогональном направлении по оси у на расстояние (у=-b) и по высоте на толщину h водного слоя, посылая токовые импульсы в исследуемую среду при всех положениях дипольного источника от точки с координатами [(x=-L), (y=-b), (z=+h)] до точки с координатами [(x=+L), (y=-b), (z=+h)], затем дипольный источник разворачивают и переводят на другой параллельный токовый профиль, проложенный и сдвинутый в плане относительно измерительного в противоположную сторону по оси у на расстояние (y=+b) и по высоте на толщину h, и продолжают измерение при его движении в обратном направлении от точки с координатами [(x=+L), (y=+b), (z=+h)] до точки с координатами [(x=-L), (у=+b), (z=+h)];
используя измеренные значения первых разностей электрических потенциалов определяют два независимых от силы тока источника и горизонтальной компоненты плотности тока (jx=0 и jy=0) в точке зондирования интерпретируемых параметра - один Rz(t0) на основе геометрического зондирования при всех положениях дипольного источника, вычисляемый по формуле
Rz(t0)=[ΔUM1N(IB1A1, t0)-ΔUM3N(IB1A1, t0)]:
:{[ΔUM1N(IB1A1, t0)+ΔUM2N(IB1A1, t0)+ΔUM3N(IB1A1, t0)+ΔUM4N(IB1A1, t0)]+
+k1(t0)·[ΔUM1N(IA2B2, t0)+ΔUM2N(IA2B2, t0)+ΔUM3N(IA2B2, t0)+ΔUM4N(IA2B2, t0)]+
+k2(t0)·[ΔUM1N(IA3B3, t0)+ΔUM2N(IA3B3, t0)+ΔUM3N(IA3B3, t0)+ΔUM4N(IA3B3, t0)]+
+k3(t0)·[ΔUM1N(IB4A4, t0)+ΔUM2N(IB4A4, t0)+ΔUM3N(IB4A4, t0)+ΔUM4N(IB4A4, t0)]},
и другой Rz(ti) на основе зондирования на переходных процессах при четырех, выбранных методом итераций, наиболее информативных разносах с координатами дипольного источника [(х=-а), (y=-b), (z=+h)], [(х=+а), (y=-b), (z=+h)], [(x=+a), (y=+b), [(z=+h)] и [(х=-а), (y=+b), (z=+h)] из всех прозондированных, вычисляемый по формуле
Rz(ti)=[ΔUM1N(IB1A1, ti)-ΔUM3N(IB1A1, ti)]:
:{[ΔUM1N(IB1A1, ti)+ΔUM2N(IB1A1, ti)+ΔUM3N(IB1A1, ti)+ΔUM4N(IB1A1, ti)]+
+k1(ti)·[ΔUM1N(IA2B2, ti)+ΔUM2N(IA2B2, ti)+ΔUM3N(IA2B2, ti)+ΔUM4N(IA2B2, t0)]+
+k2(t0)·[ΔUM1N(IA3B3, ti)+ΔUM2N(IA3B3, ti)+ΔUM3N(IA3B3, t0)+ΔUM4N(IA3B3, ti)]+
+k3(ti)·[ΔUM1N(IB4A4, ti)+ΔUM2N(IB4A4, ti)+ΔUM3N(IB4A4, ti)+ΔUM4N(IB4A4, ti)]},
где
ki(t0), k2(t0), k3(t0), - коэффициенты фокусировки при геометрическом зондировании, обеспечивающие эквипотенциальность замкнутой линии, проходящей через четыре внешних электрода измерительной установки и, тем самым, исключение внутри области, ограниченной контуром этой линии горизонтальной составляющей плотности тока jxy в каждой точке зондирования на всех геометрических разносах в период импульса тока, которые определяют из системы трех уравнений:
k1(ti), k2(ti), k3(ti), - коэффициенты фокусировки при зондировании на переходных процессах в каждой точке зондирования в паузе тока на всех временах переходных процессов, определяемые из системы трех уравнений, и обеспечивающие эквипотенциальность замкнутой линии, проходящей через четыре внешних электрода измерительной установки:
t0 - момент времени при пропускании токового импульса, когда электрическое поле переходных процессов не отличается от своего установившегося значения, соответствующего постоянному току;
ti - моменты времени, при которых измеряют сигналы переходных процессов через равные интервалы времени Δt на протяжении всей паузы после выключения тока;
ΔUM1N(IB1A1, t0), ΔUM2N(IB1A1, t0), ΔUM3N(IB1A1, t0), ΔUM4N(IB1A1, t0),
ΔUM1N(IA2B2, t0), ΔUM2N(IA2B2, t0), ΔUM3N(IA2B2, t0), ΔUM4N(IA2B2, t0),
ΔUM1N(IA3B3, t0), ΔUM2N(IA3B3, t0), ΔUM3N(IA3B3, t0), ΔUM4N(IA3B3, t0),
ΔUM1N(IB4A4, t0), ΔUM2N(IB4A4, t0), ΔUM3N(IB4A4, t0), ΔUM4N(IB4A4, t0),
- мгновенные значения первых разностей электрических потенциалов между каждым из четырех внешних электродов М1М2М3М4 измерительной установки и центральным N, измеренные при времени t0 пропускания тока в дипольном источнике при его прохождении по обоим токовым профилям через равноудаленные от центрального электрода N измерительной установки четыре точки, соответственно, с координатами [-х, (y=-b), (z=+h)], [+х, (y=-b), (z=+h)], [+х, (y=+b), (z=+h] и [-х, (y=+b), (z=+h)];
ΔUM1N(IB1A1, ti), ΔUM2N(IB1A1, ti), ΔUM3N(IB1A1, ti), ΔUM4N(IB1A1, ti),
ΔUM1N(IA2B2, ti), ΔUM2N(IA2B2, ti), ΔUM3N(IA2B2, ti), ΔUM4N(IA2B2, ti),
ΔUM1N(IA3B3, ti), ΔUM2N(IA3B3, ti), ΔUM3N(IA3B3, ti), ΔUM4N(IA3B3, ti),
ΔUM1N(IB4A4, ti), ΔUM2N(IB4A4, ti), ΔUM3N(IB4A4, ti), ΔUM4N(IB4A4, ti),
- мгновенные значения первых разностей электрических потенциалов между каждым из четырех внешних электродов M1M2M3M4 измерительной установки и центральным N, измеренные при времени ti переходных процессов в паузе тока в дипольном источнике при его прохождении по обоим токовым профилям через равноудаленные от центрального электрода N измерительной установки четыре точки, соответственно с координатами [(х=-а), (y=-b), (z=+h)], [(x=+a), (y=-b), (z=+h)], [(x=+a), (y=+b), (z-+h)] и [{х=-а), (y=+b), {z-+h)].
2. Способ морской геоэлектроразведки по п.1, отличающийся тем, что токовые профили для прохождения горизонтального дипольного источника прокладывают на заданной глубине моря или на дневной поверхности моря.
3. Способ морской геоэлектроразведки по п.1, отличающийся тем, что расстояние L равно шести километрам и более, расстояния а и b равны одному километру и более, глубина погружения h дипольного источника устанавливается в пределах от двадцати метров выше уровня морского дна до дневной поверхности.
4. Способ морской геоэлектроразведки по п.1, отличающийся тем, что измерительные установки размещают на измерительном профиле с равным шагом в 200-1000 м.