Способ электромагнитного каротажного изопараметрического зондирования

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

 

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

Известен способ электромагнитного каротажного зондирования, основанный на измерении характеристик переменного гармонического магнитного поля комплексом геометрически подобных трехкатушечных зондов, состоящих из генераторной и пары разнесенных друг от друга измерительных катушек. Генераторные катушки питают токами синхронных гармонических частот, значения которых обратно пропорциональны длинам зондов, и измеряют разности фаз между векторными значениями ЭДС синхронных колебаний. Устройства, реализующие способ, отличаются количеством зондов и тем же количеством измеряемых характеристик электромагнитного поля. Устройства могут отличаться видом измеряемой характеристики магнитного поля (амплитуды, фазы и различные функции этих величин), формируемых на основе измеряемых откликов поля от среды.

Известный способ реализован в устройстве для каротажного электромагнитного зондирования согласно авторскому свидетельству СССР №1004940, МПК G01V 3/18, 1983, а также в устройстве для электромагнитного каротажа скважин согласно патенту РФ №206305392875, МПК 6 G01V 3/18. Опубликован 27.06.96. Бюллетень №18.

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

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

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

Поставленная цель достигается тем, что в способе, при котором измерения проводят набором из нескольких зондов, включающих генераторные и измерительные катушки, предлагается зондирование выполнять набором из N-го количества двухкатушечных зондов различной длины: Z₁<Z₂<Z₃<…<Z_N-1<Z_N, состоящих их одной генераторной и одной измерительной катушки, генераторные катушки зондов питать током различных частот: ω₁>ω₂>ω₃>…>ω_N-1>ω_N при равных величинах произведений длины зонда в квадрате на частоту и измерять каждой смежной парой зондов Z₁ и Z₂, Z₂ и Z₃, …, Z_N-1 и Z_N пары фаз φ₁ и φ₂, φ₂ и φ₃, …, φ_N-1 и φ_N несинхронных колебаний ЭДС гармонических частот, по которым вычислять для каждой пары смежных зондов приведенную разность фаз путем умножения вычитаемой фазы на величину отношения частот большую единицы и(или) путем умножения уменьшаемой фазы на величину отношения частот меньшую единицы по которым определять удельную электропроводность среды с большей и(или) меньшей разрешающей способностью.

Величины отношения частот для вычитаемых фаз в приведенных разностях фаз могут быть определены в пределах

Величины отношение частот для уменьшаемых фаз в приведенных разностях фаз могут быть определены в пределах

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

На фиг.1 представлено схематическое изображение системы двухкатушечных зондов для устройства, реализующего предлагаемый способ, основанный на применении несинхронных гармонических колебаний. Зонды пронумерованы в порядке возрастания их длины от 1 до N. Каждый из зондов состоит из генераторной катушки 1.1-1.N и измерительной катушки 2.1-2.N.

На фиг.2 показаны результаты сравнительных расчетов для приведенной разности фаз Δψ (ось ординат) от обобщенного параметра р (ось абсцисс) для двух пар коэффициентов приведения (отношения частот) по предлагаемому способу и данные прототипа ВИКИЗ. [Технология исследования нефтегазовых скважин на основе ВИКИЗ. Методическое руководство. Ред. Эпов М.И., Антонов Ю.Н. Новосибирск: НИЦ ОИГГМ СО РАН, Издательство СО РАН, 2000.]

Предлагаемый способ электромагнитного каротажного изопараметрического зондирования основан на теории поля магнитного диполя в сферических и цилиндрических координатных системах [А.А.Кауфман. Введение в теорию геофизических методов. Часть 2. Электромагнитные поля. Перевод с англ. - М., Недра. 2000].

В сферической системе координат уравнения для электрической Е_ϕ и магнитных H_R, Н_θ компонент поля магнитного диполя, момент которого ориентирован по оси θ=0, определены для однородной изотропной среды

где R - расстояние от центра (О) диполя до точки наблюдения (Р);

ϕ - угол между полуплоскостями ZOX и ZOP;

θ - угол ZOP между лучами OZ и ОР;

ω - круговая частота колебаний магнитного поля;

µ - абсолютная магнитная проницаемость среды, равная 4π·10^-7 [Гн/м];

М - момент переменного магнитного диполя, которым может служить многовитковая катушка с током I=I₀e^-iωt, где I₀ - амплитуда тока;

k - волновое число (приведено ниже).

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

С учетом этого компоненты поля магнитного диполя е_ϕ, h_R, h_θ в проводящей поляризующейся среде в относительных единицах запишутся в виде

В цилиндрических координатах относительная компонента магнитного поля h_z на оси магнитного диполя (в соответствии с (5 и 8)), описывается следующим выражением:

здесь Z - расстояние от центра генераторного диполя по оси z;

- вертикальная компонента поля на оси диполя в воздухе.

Волновое число определяется электрическими параметрами среды и частотой колебания магнитного поля диполя:

Здесь

и

Здесь а и b - соответственно реальная и мнимая компоненты волнового числа k,

ε - абсолютная диэлектрическая проницаемость среды - ε=ε_отнε₀,

ε₀ - диэлектрическая проницаемость вакуума, равная 8,85·10^-12 [Гн/м];

σ - удельная проводимость среды [См/м].

Как видно из выражений (12) и (13), влияние параметров среды на поле в сильной а мере зависит от отношения, что адекватно соотношению токов проводимости к токам смещения.

Реальная компонента а волнового числа является фазовой постоянной и определяет длину волны в среде - (расстояние, на котором фаза плоской волны меняется на 2π. Мнимая компонента b носит название коэффициента поглощения и определяет глубину проникновения поля - , называемую часто скин-слоем. Глубина b проникновения поля соответствует расстоянию в [м], на котором поле плоской волны убывает в е раз.

Произведя в уравнении (10) разделение на реальную и мнимую компоненты поля, получаем

Относительная амплитуда поля на оси диполя запишется в виде

Фаза гармонического колебания на осевом расстоянии Z от магнитного диполя

В известном способе электромагнитного индукционного зондирования измеряют разность фаз между двумя синхронными колебаниями. В точках размещения измерительных катушек на расстояниях Z₁ и Z₂ (Z₁>Z₂) от генераторной катушки разность фаз запишется, в соответствии с (17), в следующем виде:

Соответственно, разность амплитуд, с учетом (14, 15) и (16), представляется в виде

При регистрации амплитуд на расстояниях Z₁ и Z₂ необходимо обеспечивать равенство прямых полей в соответствии с уравнением (10) - .

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

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

в уравнениях для компонент волнового числа (12) и (13).

При этом реальная и мнимая компоненты волнового числа

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

Выражение для фазы гармонического колебания поля магнитного диполя, в соответствии с условием (20), определяется следующим соотношением:

Здесь - безразмерный параметр; i=1 и 2, где 1 индекс длинного зонда, 2 - индекс короткого зонда в составе трехкатушечного зонда.

Численная величина р_i зависит от длины зонда Z_i и частоты колебаний электромагнитного поля - ω. Произведение этих параметров определяет основные характеристики трехкатушечных зондов с N-м их количеством в устройстве, реализующем способ-прототип.

В уравнениях (22) - const₁>const₂. Этими константами определяются численные величины электродинамического изопараметра с размерностью .

Примером устройства, реализующего известный способ электромагнитного каротажного изопараметрического зондирования, является комплекс ВИКИЗ. [Технология исследования нефтегазовых скважин на основе ВИКИЗ. Методическое руководство. Ред. Эпов М.И., Антонов Ю.Н. - Новосибирск: НИЦ ОИГГМ СО РАН, Издательство СО РАН, 2000.]

Комплекс ВИКИЗ состоит из N-го количества трехкатушечных зондов, которые пронумерованы цифрами в порядке увеличения длины 7, 2, 3, …, N-1, N. Генераторные катушки зондов питают гармоническими частотами ω₁>ω₂>ω₃>…>ω_N-1>ω_N, которые уменьшаются по мере увеличения длины зондов.

Все возможные отношения длин между одноименными катушками смежных зондов являются величинами постоянными, что соответствует их геометрическому подобию. Характеристической величиной трехкатушечных зондов является относительная длина между измерительными катушками в долях максимальной длины трехкатушечного зонда: ΔZ_N=(Z₁-Z₂)/Z₁=const₃. Эта величина является геометрическим изопараметром трехкатушечных зондов.

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

Выражение для разности фаз изопараметрических зондов, в соответствии с уравнением (18) и условия (20), имеет вид:

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

Геометрический изопараметр трехкатушечных зондов не зависит от частоты:

Здесь ΔZ_N - величина постоянная для всех трехкатушечных зондов с синхронными гармоническими колебаниями.

Таким образом, в известном способе - прототипе, реализованном в устройстве ВИКИЗ, измеряют разности фаз между парами синхронных гармонических колебаний S₁ и S₂, которые имеют вид:

Векторная разность амплитуд является разностью фаз или сдвигом фаз двух синхронных колебаний

Разность фаз двух синхронных колебаний не зависит от выбора начала счета времени.

Действительно, пусть t - время, отсчитываемое по первому секундомеру вступление сигнала в ближнюю измерительную катушку, t* - время вступления в дальнюю катушку зонда по второму секундомеру. Оба секундомера имеют одинаковый ход. Очевидно, что сигнал во второй катушке появится с запозданием на время τ. Следовательно, мгновенный отсчет времени в разнесенных измерительных катушках зонда будет иным, а именно:

Подставляя этот отсчет в уравнения для S₁ и S₂, получаем уравнения для обоих колебаний при счете времени вступления сигнала во вторую, более удаленную, катушку зонда:

где

Таким образом, фазы φ₁* и φ₂* отличны от φ₁ и φ₂ при отсчете времени по второй измерительной катушке, более удаленной от источника поля. Однако, как следует из приведенного уравнения (30)

разность фаз между фазами двух синхронных колебаний не зависит от выбора начала счета времени.

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

В предлагаемом способе измерения разности фаз выполняют между асинхронными (несинхронными) гармоническими колебаниями. Базами измерения являются пары смежных двухкатушечных зондов, генераторные катушки которых питаются токами различных гармонических частот. Длины двухкатушечных зондов и их различные частоты определяются следующим уравнением:

Согласно уравнению (31) каждая пара смежных зондов, например, по мере уменьшения длины является измерительной базой несинхронных гармонических колебаний электромагнитных полей. Такими измерительными базами являются пары зондов со следующими номерами в соответствии с фиг.1: (N и N-1), …, (3 и 2), (2 и 1). Указанными парами зондов измеряют фазы несинхронных колебаний: и , …, , и , и . Разность фаз между несинхронными колебаниями вычисляют приведением одной из фаз к единому счету времени со второй фазой. Таким путем получают приведенную разность фаз несинхронных колебаний.

Сущность приведенной разности фаз несинхронных колебаний рассмотрим на примере двух различных гармонических колебаний двух смежных зондов с номерами 1 и 2 в соответствии с уравнением (31):

Тогда, при другом начале счета времени, имеем согласно (27)

где

откуда разность

Из (36) следует, что разность фаз (φ₁-φ₂) зависит от начала отсчета времени.

Однако от начала отсчета времени не зависят величины

определяемые как «приведенные разности фаз» несинхронных колебаний: [Г.С.Горелик. Колебания и волны. Введение в акустику, радиофизику и оптику. Государственное издательство физико-математической литературы (издание второе). М., 1959. С.29].

Соотношения (37) нетрудно проверить путем приравнивания уравнений (34) и (35) относительно времени τ:

Умножая левую и правую части уравнения (38) на ω₁ или ω₂(ω₁>ω₂), получаем два значения приведенных разностей фаз, которые отличаются частотной зависимостью:

1) фаза, соответствующая меньшей частоте колебания, умножается на отношение частот - , что соответствует максимальной частотной зависимости приведенной разности фаз

2) фаза, соответствующая большей частоте колебания, умножается на отношение частот - , что соответствует минимальной частотной зависимости приведенной разности фаз

В соответствии со схемой на фиг.1 и уравнений (39) и (40) для набора зондов приведенные разности фаз Δψ будут выражены двумя рядами.

1. Первый ряд разностей фаз приводят к фазе более высокой частоты

2. Второй ряд разностей фаз приводят к фазе низкой частоты

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

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

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

В качестве примера на фиг.2 приведены результаты расчета приведенных разностей фаз - Δψ (ось ординат в градусах) в зависимости от обобщенного параметра - p_N (ось абсцисс) для модели однородной изотропной среды. Шифры кривых соответствуют предельным величинам отношения частот. Одна пара предельных отношений большей частоты к меньшей частоте находится в диапазоне от 1,2 до 1,4. Другая пара предельных отношений частот заключена между значениями от до или - от 0,714 до 0,833. Для сопоставления приведен график тех же зависимостей для зондового комплекса ВИКИЗ.

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

1. Способ электромагнитного изопараметрического каротажного зондирования, в котором измерения проводят набором зондов, состоящих из генераторных и измерительных катушек, отличающийся тем, что зондирование выполняют набором из N-го количества 2-х катушечных зондов различной длины: Z₁<Z₂<Z₃<…<Z_N-1<Z_N, состоящих из одной генераторной и одной измерительной катушки, генераторные катушки зондов питают током различных гармонических частот: ω₁>ω₂>ω₃>…>ω_N-1>ω_N при равных величинах произведений длины зонда в квадрате на частоту: и измеряют каждой смежной парой зондов Z₁ и Z₂, Z₂ и Z₃, …, Z_N-1 и Z_N пары фаз φ₁ и φ₂ ф₂ и φ₃, …, φ_N-1 и φ_N несинхронных колебаний ЭДС гармонических частот, по которым вычисляют для каждой пары смежных зондов приведенную разность фаз путем умножения вычитаемой фазы на величину отношения частот большую единицы ; ; …; и (или) путем умножения уменьшаемой фазы на величину отношения частот меньшую единицы ; ; …; , по которым определяют удельную электропроводность среды с большей и (или) меньшей разрешающей способностью.

2. Способ по п.1, отличающийся тем, что величины отношения частот для вычитаемой фазы определяют в пределах

3. Способ по п.1 или 2, отличающийся тем, что величину отношения частот для уменьшаемой фазы определяют в пределах