СПОСОБ ИНДУКЦИОННОГО КАРОТАЖА
INDUCTION WELL LOGGING

Предлагаемый способ относится к области геофизических исследований переменным электромагнитным полем, возбуждаемым индуктивным способом в геологоразведочных скважинах, и предназначен для определения электрического сопротивления вмещающих горных пород.
Область преимущественного применения - изучение электрических свойств горных пород в геологоразведочных скважинах.
Известен способ индукционного каротажа скважин, предложенный Г. Доллем в 1946 г. [1] и развитый в работах Аксельрода СМ., Кауфмана А.А. и Плюснина М.И. [2, 3, 4] и реализующая этот способ аппаратура различных российских и зарубежных фирм.
Этот способ индукционного каротажа основан на возбуждении индукционных токов во вмещающих горных породах гармоническим магнитным полем от тока в генераторной катушке с моментом, направленным вдоль оси скважины. Сопротивление вмещающих пород определяют по величине индуцируемого токами вторичного магнитного поля, измеряемого приемной катушкой, соосной генераторной [1-4].
Развитием индукционного каротажа с определением удельного сопротивления пород являются разработка способов высокочастотного изопараметрического каротажного индукционного зондирования скважин (ВИКИЗ) [5] и высокочастотного электромагнитного метода каротажа с определением диэлектрической проницаемости пород [6].
Однако известным способам присущи существенные недостатки. В способах индукционного каротажа, использующих соосные генераторные и измерительные катушки, с моментами вдоль оси скважины, определение удельного сопротивления пород возможно до величин не более 100 Ом·м [3, 4].
В работе [7] предложено выполнять индукционный каротаж с магнитными передатчиками и приемниками (индукционные катушки), ориентированными в поперечном направлении относительно продольной оси скважинного прибора, снимающими это ограничение.
Дальнейшим развитием индукционного каротажа является разработка многозондовых установок, содержащих системы трех ортогональных генераторных и измерительных зондов и измерений магнитного поля на разных частотах [8-10]. К недостаткам данного способа относятся сложность ее технического исполнения и необходимость применения мощных вычислительных систем для интерпретации измеренных данных [8].
Наиболее близким техническим решением является способ [11], взятый нами в качестве способа-прототипа.
В способе-прототипе, который относится к методам аэроэлектроразведки, электромагнитное поле в окружающем пространстве возбуждают гармоническими токами в двух взаимоперпендикулярных генераторных контурах, размещенных на одном самолете, сила тока в которых сдвинута по фазе на 90° относительно друг друга. Равные по величине магнитные моменты контуров направлены ортогонально направлению движения самолета. Результирующий вектор переменного магнитного поля вращается вокруг оси направления движения с частотой, равной частоте возбуждающего тока в генераторных контурах.
Измерения магнитного поля выполняют на некотором расстоянии от генераторных контуров с помощью второго самолета, буксирующего гондолу с двумя совмещенными измерительными контурами с равными моментами, ориентированными ортогонально друг другу и направлению движения самолетов. Ось гондолы с измерительными контурами находится на продолжении оси первого самолета. В однородной среде наводимые в измерительных контурах величины ЭДС, пропорциональные составляющим магнитного поля, равны по амплитуде, но сдвинуты по фазе на 90°.
Однако способ-прототип имеет недостаток, связанный с тем, что измерения магнитного поля проводят в условиях сбалансированного поля в однородной среде. Это достигается тем, что величина ЭДС в одном из измерительных контуров, сдвинутая по фазе на 90°, вычитается из ЭДС во втором измерительном контуре и общий сигнал в однородной среде равен нулю. Наличие проводящего тела около установки нарушает балансировку, что и отмечается регистратором, размещенным на борту второго самолета [11]. Таким образом, данный способ наблюдений пригоден для обнаружения локальных проводящих объектов и применяется для поисков проводящих руд. Этот способ не может быть использован для определения величины удельного сопротивления вмещающей среды, поскольку ее влияние на измеряемое магнитное поле является скомпенсированным.
Цель предлагаемого технического решения - увеличение диапазона определения удельного сопротивления горных пород и увеличение дальности зоны исследований при проведении геофизических исследований в геологоразведочных скважинах.
x
x
x
x
y
y
y
y
х
y
x
y
Поставленная цель достигается тем, что в предлагаемом способе индукционного каротажа индукционные токи во вмещающих горных породах возбуждают гармоническим магнитным полем от токов, равных по амплитуде, но находящихся в противофазе с частотой ω, в системе двух генераторных катушек с равными моментами М=JNS=JNS=Μ, ориентированными ортогонально друг другу и оси скважины Ζ, и перемещаемых по исследуемой скважине, а по фазочувствительным измерениям квадратур Re и Jm горизонтальных составляющих магнитной индукции В и В на различных расстояниях L и L по оси скважины определяют удельное электрическое сопротивление вмещающих горных пород.
Величина радиального момента общей системы двух генераторных катушек в течение одного периода остается постоянной, а направление момента изменяется по азимуту, описывая окружность с круговой частотой ω.
Сущность изобретения иллюстрируется следующими чертежами.
x
X
O
x
X
X
x
x
X
y
y
Y
2
На фиг. 1 изображена структурная схема устройства, с помощью которого реализуется предлагаемый способ. На фиг. 2 представлена зависимость Reh составляющей магнитного числа, создаваемого магнитным диполем с моментом Μ, от волнового параметра поля ρ=ωμL/ρ. На фиг. 3 представлена зависимость Jmh составляющей магнитного числа, создаваемого магнитным диполем с моментом Μ, от волнового параметра поля р. На фиг. 4 представлена зависимость фазового угла  создаваемого гармоническим магнитным диполем с моментом Μ, от волнового параметра поля р. На фиг. 5 представлены кривые ReB и JmB составляющей магнитной индукции, создаваемой магнитным диполем с моментом Μ при разносе Lx=2 м, по оси Z, пересекающей горизонтальный контакт пород с различными удельными сопротивлениями. На фиг. 6 представлены кривые ReB и JmB составляющей магнитной индукции, создаваемой магнитным диполем с моментом Μ при разносе Ly=3 м, по оси Z, пересекающей горизонтальный контакт пород с различными удельными сопротивлениями.
X
Y
Устройство (фиг. 1) содержит: силовой генератор 1, датчики фаз тока 2 и 3, совмещенные генераторные катушки 4 и 5 с горизонтальными моментами Μ и Μ, измерительные катушки 6 и 7 горизонтальных составляющих магнитной индукции Вх и By, усилители сигналов 8 и 9, аналого-цифровые фазочувствительные измерители 10 и 11, регистрирующее устройство 12.
Предлагаемый способ реализуется следующим образом.
Электромагнитное поле в окружающем пространстве создают гармоническими токами, протекающими в возбуждающих катушках 4 и 5 с равными ортогональными моментами, с помощью силового генераторного устройства 1, сила тока в которых изменяется в противофазе с круговой частотой ω. Выходные напряжения с измерительных катушек 6 и 7, пропорциональные величинам горизонтальных составляющих вектора магнитной индукции, поступают через усилители сигналов 8 и 9 в аналого-цифровые фазочувствительные преобразователи 10 и 11, в которых с помощью сигналов с датчиков фаз токов 2 и 3 осуществляется определение значений параметров квадратур Re и Jm величин составляющих Вх и By магнитной индукции, регистрируемые в устройстве 12.
X
Переменный магнитный диполь с моментом Μ=JSN в безграничной среде с удельным электрическим сопротивлением ρ создает в экваториальной плоскости магнитное поле с единственной составляющей индукции, совпадающей с моментом [12]
- волновое число среды.
2
2
2
O
x
В случае x=y=0, R=z≡|L| - разнос между диполем и точкой измерений, обозначив kL=jωµL/ρ=jp, где p - вещественный безразмерный волновой параметр поля, безразмерное магнитное число h имеет вид:
y
Y
Аналогично описывается h - составляющая магнитного числа гармонического магнитного диполя с моментом Μ.
O
2
Зависимость горизонтальной составляющей магнитного числа магнитного диполя от волнового параметра поля p=ωµL/ρ для известных величин частоты ω и разноса L используют для определения удельного сопротивления вмещающей среды.
x
X
На фиг. 2-4 представлены зависимости квадратур Re и Jm составляющей магнитного числа h и фазового угла  создаваемого гармоническим магнитным диполем с моментом Μ, от волнового параметра поля р.
x
Областями однозначного определения волнового параметра поля p в низкочастотной области спектра по измерениям квадратур Rehx, Jmhx и фазового угла φh являются:
x
- для квадратуры Re h: р<3.2;
x
- для квадратуры Jm h: р<0.7;
x
- для фазового угла φh: p<0.6.
Наилучшим разрешением для определения волнового параметра поля p по измерениям магнитного поля обладает квадратура Rehx, амплитуда которой в диапазоне 0.01<р<3.2 изменяется в 1.4 раза относительно величины поля на постоянном токе (р<0.01). В то же время квадратура Jmhx обладает большей чувствительностью к изменениям волнового параметра поля р, поскольку измеряется относительно нулевой величины при p<0.01.
O
X
Y
2
Для однозначного определения волнового параметра поля и удельного сопротивления среды ρ=ωμL/p измерения двух горизонтальных составляющих Вх и By выполняют на различных разносах Lx и Ly в фазе с соответствующим им максимальным амплитудам моментов магнитных диполей Μи Μ.
Величины разносов L и оптимальные частоты ω измерений магнитного поля определяются по экстремальным значениям удельного электрического сопротивления вмещающей среды по всему разрезу скважины в области однозначных значений определения волнового параметра поля р.
Как видно из фиг. 2-4, предлагаемый способ позволяет определять удельное электрическое сопротивление вмещающих горных пород по значениям измеряемых величин горизонтальных составляющих магнитной индукции, возникающих от токов во вмещающих горных породах, индуцируемых магнитным полем от гармонических токов, равных по амплитуде и находящихся в противофазе в системе двух генераторных катушек с равными моментами, ориентированных ортогонально друг другу и к оси скважины и перемещаемых по исследуемой скважине.
Предлагаемый способ обладает большим радиусом исследования пород в окрестности скважины, поскольку магнитное поле магнитного диполя в квазистационарном случае на полярной оси вдвое превосходит по величине магнитное поле в экваториальной плоскости [12], используемое в известных способах индукционного каротажа [1, 5, 6].
Для определения горизонтальных границ сред, обладающих различными удельными электрическими сопротивлениями, с помощью предложенного способа выполнено математическое моделирование для подошвенной установки (точка измерений расположена ниже источника электромагнитного поля) с магнитным диполем, момент которого параллелен границе раздела сред.
x
x
X
1
2
2
6
На фиг. 5 представлены кривые ReB и JmB по оси Ζ составляющей магнитной индукции, создаваемой магнитным диполем с моментом Μ=1 А*м, пересекающей горизонтальный контакт пород с удельными сопротивлениями ρ=500 Ом·м, ρ=100 Ом·м при разносе Lx=2 м на частоте f=10Гц.
y
y
Y
1
2
2
6
На фиг. 6 представлены кривые ReB и JmB по оси Ζ составляющей магнитной индукции, создаваемой магнитным диполем с моментом Μ=1 А*м, пересекающей горизонтальный контакт пород с удельными сопротивлениями ρ=500 Ом·м, ρ=100 Ом·м при разносе Ly=3 м на частоте f=10Гц.
1
2
1
2
х
у
При наличии горизонтальной границы раздела сред с различными удельными сопротивлениями перемещение установки по скважине из среды с удельным сопротивлением ρ в среду с удельным сопротивлением ρ приводит к изменению величин составляющих магнитного поля вблизи границы сред в зависимости от контрастностей сред по удельному электрическому сопротивлению ρ/ρ, а также характеристик установки (моментов генераторных катушек М и М, частоты тока f, разносов между генераторными и измерительными катушками Lx и Ly).
x
x
y
y
Из графиков, представленных на фиг. 5 и 6, следует, что контакт сред с различными удельными электрическими сопротивлениями можно определить по положениям характерных точек локальных экстремумов на кривых ReB, JmB, ReB и JmB с учетом разносов Lx и Ly.
Как видно из фиг. 5 и 6, предлагаемый способ с использованием питающей установки, создающей вращающееся гармоническое магнитное поле, позволяет определять горизонтальные границы горных пород, обладающие различными удельными электрическими сопротивлениями, по поведению горизонтальных составляющих магнитной индукции, измеряемых вдоль оси скважины.
Сущность заявляемого изобретения выражается в совокупности существенных признаков, достаточных для достижения технического результата, который выражается в увеличении диапазона определения удельного электрического сопротивления горных пород в скважинах и большем радиусе исследуемой зоны.
Заявленная совокупность существенных признаков находится в прямой причинно-следственной связи с достигаемым результатом. Анализ современного уровня техники показал, что предлагаемое техническое решение соответствует критериям "новизна" и "изобретательский уровень" и может быть промышленно реализовано при использовании существующих технических средств.
Источники информации
1. Долль Г. Теория индукционного метода исследования разрезов скважин и его применение в скважинах, пробуренных с глинистым раствором на нефти. // Вопросы промысловой геофизики. - М.: Гостоптехиздат, 1957. С. 252-274.
2. Аксельрод С.М. Способ индукционного каротажа. А.с №133536. БИ №22, 1960.
3. Кауфман А.А. Теория индукционного каротажа. - Новосибирск: Наука, 1965. 236 с.
4. Плюснин М.И. Индукционный каротаж. - М.: Недра, 1968. 142 с.
5. Антонов Ю.Н. Изопараметрическое каротажное зондирование (обоснование ВИКИЗ). // Геология и геофизика, 1980. №6. С. 81-91.
6. Даев Д.С. Способ диэлектрического индуктивного каротажа. А.с. №206735. Изобретения, пром. образцы, товарные знаки. №18, 1969.
7. Табаровский Л.А., Эпов М.И. Радиальные характеристики индукционных фокусирующих зондов с поперечными датчиками в анизотропной среде// Геология и геофизика, №7, 1979, с. 69-77.
8. Андерсон Б., Бербер Т., Леверидж Р. Трехмерный индукционный каротаж: старые измерения под новым углом. Нефтегазовое обозрение. 2008. URL: ia/Files/resources/oilfield_review/russia08/sum08/05_triaxialinduction.pdf.
9. Dzevat Omeragic, Keli Sun, Qiming Li, Tarek Habashy. Multilevel workflow method for extracting resistivity anisotropy data from 3D induction measurements. EP №2320251 A2. 2011.
10. Gupta P., Kriegshauser B.F., Fanini O.N. Conductivity anisotropy estimation method for inversion processing of measurements made by a transverse electromagnetic induction logging instrument. Патент US 5999883. 1998. URL: http://www.google.com/patents/US5999883.html.
11. Шауб Ю.Б. Основы аэроэлектроразведки методом вращающегося магнитного поля. - Л.: Гостоптехиздат, 1963. 228 с.
12. Бурсиан В.Р. Теория электромагнитных полей, применяемых в электроразведке. - Л.: Недра, 1972. 368 с.