Применение индукционного метода в нефтегазовой геофизике
Введение…………………………………………………………
1 Физические основы метода индукционного метода……………………..6
1.1 Приближенная теория низкочастотных индукционных методов……..7
1.2 Пространственный фактор элементарного кольца……………………..
2 Техника и методика работ
Обычный низкочастотный
индукционный метод с
2.1Зонды обычного низкочастотного индукционного метода……………..
3 Интерпретация диаграмм метода индукционного каротажа……………..
3.1Формы кривых и определение границ пластов…………………………..
3.2 Задачи, решаемые методом индукционного каротажа………………….
3.3 Кривая ИМ………………………………………………………
4 Применение индукционного метода в нефтегазовой геофизике…………
Введение
Геофизические методы исследования скважин - один из разделов
прикладной геофизики. Они применяются
для решения геологических и технических
задач, связанных с поисками, разведкой
и разработкой месторождений полезных
ископаемых, а также с изучением гидрогеологических
и других особенностей исследуемых районов.
Наиболее широкое применение геофизические методы получили при изучении нефтяных и газовых скважин в процессе их бурения, опробования и эксплуатации, т. е. в промысловой геофизике.
Исследование скважин
Изучение геологических разрезов скважин - наиболее важное направление. При этом используются электрические, магнитные, радиоактивные, термические, акустические, механические, геохимические и другие методы. Применение их основано на изучении физических, естественных и искусственных полей различной природы. Интенсивность того или иного поля определяется разнообразными факторами, в первую очередь физическими свойствами горных пород - электрическим удельным сопротивлением диэлектрической и магнитной проницаемостью, электрохимической активностью, радиоактивностью, и т. д. При геофизических исследованиях скважин регистрируются диаграммы или производятся точечные измерения различных физических параметров: кажущегося электрического сопротивления, потенциалов собственной и вызванной поляризации пород, силы тока, сопротивления заземления, электродных потенциалов, интенсивности гамма-излучения, плотности тепловых и надтепловых нейтронов, температуры, напряженности магнитного поля, скорости и времени распространения упругих колебаний, продолжительности бурения и др. Характер изменения указанных параметров по стволу скважины, отраженный на регистрируемых диаграммах, зависит от многих факторов и главным образом от физических свойств пород. Последние в свою очередь находятся в тесной связи с литолого-петрографической, коллекторской, продуктивной и другими характеристиками. Разработанная теория геофизических методов и петрофизические зависимости позволяют производить научно обоснованную интерпретацию результатов исследований. При изучении геологических разрезов скважин на основе интерпретации комплекса данных геологической и геофизической документации решаются следующие задачи: 1) геофизическое расчленение разрезов и выявление геофизических коррелятивов (реперов); 2) определение пород, слагающих разрезы скважин; 3) выявление коллекторов и изучение их свойств (пористости, проницаемости, глинистости и др.); 4) выявление и определение местоположения различных полезных ископаемых (нефти, газа, каменного угля, каменной соли, руды, термальных, минеральных и пресных вод и др.); 5) количественная оценка нефтегазонасыщения, а в некоторых случаях угленасыщения, оруденения, а также минерализации пластовых вод.
Обобщающая интерпретация
Изучение технического состояния скважин проводится с помощью комплекса различных методов геофизики. В этом случае выполняются следующие основные операции: 1) определение искривления скважин инклинометрами - инклинометрия; 2) установление фактического диаметра скважин с помощью каверномеров - кавернометрия; 3) определение профиля сечения скважины и обсадных колонн - профилеметрия; 4) определение высоты подъема, характера распределения и степени сцепления цемента в затрубном пространстве термическим, радиоактивным, акустическим методами - цементометрия; 5) выявление мест притоков и затрубной циркуляции вод в скважинах электрическим, термическим и радиоактивным методами - притокометрия; 6) определение водопоглощающих горизонтов и контролирование гидравлического разрыва пласта термическим и радиоактивным методами; 7) определение уровней жидкости, местоположения башмаков обсадных колонн и металлических предметов, оставленных в скважинах при авариях, глубин расположения забоев скважин и решение многих других важных нефтепромысловых задач.
Контроль за разработкой месторождений нефти и газа предусматривает следующие определения: 1) динамика водонефтяных, водогазовых и газонефтяных контактов; 2) дебита и состава флюидов в скважинах; 3) профилей отдачи и приемистости пластов - дебитометрия и расходометрия; 4) интервалов прорыва нагнетаемых вод; 5) нефтеотдачи пластов.
Прострелочно-взрывные и другие работы в скважинах включают перфорацию обсадных труб для сообщения скважины с пластом, отбор образцов пород из стенок пробуренных скважин для уточнения геологического разреза и торпедирование, производимое с разными целями. Эти работы выполняются с учетом данных геофизических методов исследования скважин с помощью стреляющих перфораторов, боковых грунтоносов и торпед, спускаемых в скважины и действующих в комплекте с соответствующим геофизическим оборудованием.
Геофизические методы изучения скважин являются важнейшим и неотъемлемым звеном в геологических, буровых и эксплуатационных работах, производимых на нефтяных и газовых промыслах, угольных и рудных месторождениях, в гидрогеологических и инженерно-геологических изысканиях.
1 Физические основы индукционных методов
Индукционные методы, основанные на изучении в скважинах переменного электромагнитного поля низкой и высокой частоты, разработаны достаточно детально. Низкочастотными индукционными методами изучают электромагнитное переменное поле ультразвуковой частоты 20 - 60 кГц, высокочастотными – переменные поля частотой 1- 10 МГц.
В группу низкочастотных индукционных методов входят обычный индукционный метод с продольным датчиком, индукционный метод с поперечным датчиком, индукционный метод переходных процессов, частотный индукционный метод и др. К высокочастотным индукционным методам относятся обычный высокочастотный индукционный метод (ВИМ) (амплитудный метод), волновой метод проводимости (ВМП) (фазовый метод), высокочастотное индукционное изопараметрическое зондирование. Индукционные методы применяются для исследования вторичного электромагнитного поля среды, ЭДС которого прямо пропорциональна электропроводности горных пород. Вторичное электромагнитное поле возникает в окружающей среде за счет вихревых токов, которые индуцированы катушкой, питающейся от помещенного в скважину генератора переменного тока. Индукционные методы принципиально отличаются от всех методов стационарного и квазистационарного электрических полей, прежде µ всего, тем, что для создания вторичного электромагнитного поля в горных породах не требуется непосредственного (гальванического) контакта зондовой установки с окружающей средой. Если в методах КС, СЗ, ТМ и ВП электрический ток распространяется в горные породы от токовых электродов через слой проводящей жидкости (промывочной), то в индукционных методах электроды как таковые не используются, и вторичное электромагнитное поле формируется в горных породах за счет индуктивной связи первичного электромагнитного поля со средой, окружающей зонд. Следовательно, индукционные методы позволяют изучать разрезы скважин, пробуренных с обычной промывочной жидкостью и с промывочной жидкостью на нефтяной или другой основе, плохо проводящей электрический ток. Индукционные методы отличаются также характером распределения вторичных токов, индуцированных генераторной катушкой в горных породах: их токовые линии лежат в плоскостях, перпендикулярных к оси генераторной катушки. В однородной среде линии вихревых токов представляют собой окружности с центрами на оси прибора. При таком распределении токовых линий можно более точно определить истинное удельное сопротивление пластов, а влияние электропроводности вмещающих пород на показания индукционных методов существенно уменьшается.
Простейший зонд индукционного метода может быть составлен из двух катушек (генераторной и измерительной), опущенных в скважину. Расстояние между серединами генераторной и измерительной катушек есть длина LH индукционного зонда. Генераторная катушка зонда подключена к генератору переменного тока ультразвуковой частоты 20–60 кГц и питается стабилизированным по частоте и амплитуде током. Измерительная катушка зонда через усилитель и фазочувствительный элемент подключена посредством кабеля к регистрирующему прибору, расположенному на поверхности. Переменный ток, протекающий по генераторной катушке, создает переменное магнитное поле (прямое и первичное), которое, в свою очередь, индуцирует в среде, окружающей зонд, вихревые токи, формирующие вторичное переменное магнитное поле той же частоты, что и первичное. Первичное и вторичное переменные магнитные поля индуцируют ЭДС в измерительной катушке. Непосредственное воздействие первичного поля на приемную катушку не связано с горными породами, поэтому ЭДС, индуцированная прямым полем, компенсируется встречной ЭДС, равной первой по величине и противоположной по фазе, с помощью дополнительных катушек или специальных электронных устройств.
Электродвижущая сила, генерируемая вторичным полем в измерительной катушке, состоит из двух составляющих – активной и реактивной. Регистрирующим прибором фиксируется сигнал активной составляющей ЭДС, наиболее тесно связанной с электропроводностью окружающей среды. В случае низкой проводимости среды ЭДС активной составляющей прямо пропорциональна ее электропроводности. С ростом электропроводности среды ЭДС активного сигнала увеличивается медленнее и по более сложному закону. Нарушение пропорциональности между активным сигналом и электропроводностью среды связано с взаимодействием вихревых токов. Это явление называется скин-эффектом. Чем выше частота тока и электропроводность среды, тем значительнее взаимодействие вихревых токов и, следовательно, существеннее влияние скин-эффекта на показания индукционного метода. Активный сигнал фиксируется на поверхности измерительным устройством в виде кривой, отражающей изменение электропроводности пород по разрезу скважины. Точкой записи кривой является середина расстояния между центрами генераторной и измерительной катушек. Единицей измерения электропроводности пород является сименс на метр (См/м) – величина, обратная ом-метру (Ом∙м). На практике используют тысячную долю сименса на метр – миллисименс на метр (мСм/м).
В индукционных методах измеряется эффективная удельная электропроводимость σэф, зависящей от проводимостей пласта, промывочной жидкости, вмещающих пород, от диаметра скважины, мощности пласта, а также размера и конструкции зонда. В связи с этим эффективная электропроводность в общем случае отличается от истинной удельной электропроводности изучаемого пласта σпл.
1.1 Приближенная теория низкочастотных индукционных методов
Решение прямой задачи индукционных методов состоит в получении зависимости измеряемой ЭДС электромагнитного поля от электропроводности среды, ее геометрии и параметров зонда. Приближенная теория низкочастотного индукционного метода впервые была разработана X. Г. Доллем в 1949 г, затем получила развитие в работах С. М. Аксельрода, Ю. Н. Антонова и других исследователей. В основе приближенной теории Долля лежат два допущения: 1) все индуцированные в окружающей среде вихревые токи имеют одну и ту же фазу, сдвинутую относительно фазы тока в генераторной катушке на 90°; это означает, что взаимодействие вихревых токов отсутствует, т. е. при решении прямой задачи явление скин-эффекта не принимается во внимание; 2) амплитуда плотности тока в любой точке среды рассчитывается по упрощенной формуле и определяется только пространственным фактором и удельной электропроводностью участка среды. Эти допущения справедливы лишь тогда, когда частота тока питания и электропроводность среды сравнительно невелики. При высокой частоте тока или большой электропроводности пород явление скин-эффекта существенно изменяет характер распространения электромагнитного поля, и в этом случае при решении прямой задачи необходимо использовать строгую теорию. Приближенная теория низкочастотных индукционных методов позволяет сравнительно просто установить зависимость эффективной удельной электропроводности изучаемого пространства от электропроводности отдельных сред, их размеров и положения относительно индукционного зонда, а также наглядно представить физическую сущность этого метода. С помощью приближенной теории можно путем элементарных расчетов решить прямую задачу описываемых методов для плоских и цилиндрических поверхностей раздела между средами. Формулы, полученные на основании этой теории, следует рассматривать как асимптотические, справедливые при электромагнитных волнах, очень длинных по сравнению с радиусом скважины, мощностью пласта и диаметром, зоны проникновения фильтрата промывочной жидкости. И так, имеем однородную изотропную среду удельной электропроводности σ абсолютной диэлектрической проницаемости µа и магнитной проницаемости εа. На оси скважины расположен двухкатушечный индукционный зонд. Диаметр скважины dc→0. Центры генераторной и измерительной катушек расположены на общей оси, на расстоянии LH одна от другой. Ось генераторной катушки совпадает с осью скважины (рисунок 1). Генераторная и измерительная катушки имеют соответственно высоты lг и lп; число витков пг и пп, радиусы которых rг и rп; площади каждого витка ; общие площади витков sгnг =Sг и sпnп =Sп. Условимся считать, что размеры катушек значительно меньше расстояния Lи, т. е.(lг, rг, lп, rп) << Lп. Это допущение позволяет рассматривать катушки как точечные. Генераторную катушку можно заменить для упрощения расчетов магнитным диполем с переменным моментом.
Генераторная катушка питается переменным током с амплитудой Ia, мгновенное значение которого
I =Iаехр(-i ωτ) (1)
где τ – время; ω=2πf – угловая частота; f – циклическая частота.
Упомянутый выше магнитный диполь, ось которого совпадает с осью генераторной катушки, создает в окружающем пространстве электромагнитное поле. Необходимо определить величину ЭДС, которая наводится вихревыми токами в изучаемой среде, и установить связь между наведенной ЭДС и удельной электропроводностью однородной среды, частотой поля и параметрами зонда.
Для решения задачи введем цилиндрическую систему координат rzψ, начало которой расположим в точке O, являющейся серединой расстояния между центрами генераторной и измерительной катушек. Разобьем все изучаемое пространство на элементарные тороиды, представляющие собой участки породы с горизонтальными поверхностями. Единичный тороид – это горизонтальное кольцо радиуса rк с центром на оси скважины (см. рисунок 1). Он назван X. Г. Доллем элементарным кольцом.
С учетом (1) магнитный момент магнитного диполя
М = nг sг I= nг sг Ia eхр(-i ωτ.) (2)
Магнитный диполь создает в окружающем пространстве магнитное поле напряженностью
Н = М/2πLг3.
или, принимая во внимание (2),
(4)
где Lг – расстояние от центра генераторной катушки до рассматриваемой точки среды.
Как известно, величина магнитного потока, пронизывающего замкнутую поверхность,
(5)
где В=μ0μН= μаН – вектор магнитной индукции; φ – угол между нормалью к элементарной площадке и силовыми линиями магнитного поля; ds – площадь сечения элементарной площадки.
Для элементарного кольца с радиусом rk, площадь которого магнитные силовые силовые линии пересекают под углом 90○( φ=0), формула (5) имеет вид
(6)
(Нz – вертикальная компонента напряженности магнитного диполя, поскольку силовые линии магнитного поля направлены вверх вдоль оси скважины).
Подставив (4) в (6) и проведя интегрирование, получим
(7)
ПК-генераторная и приёмная катушки; Т-элементарное кольцо(торонд); У-уси-
литель.
Рисунок 1- Схема индукционного зонда с элементарным витком ТГК
Изменение магнитного потока во времени создает ЭДС электромагнитной индукции е в элементарном кольце. Согласно закону электромагнитной индукции Фарадея величина
е= –dФ/dτ
Знак минус в формуле (8) соответствует закону Ленца, согласно которому индукционный ток, взаимодействуя с магнитным полем, вызывает силу, направленную противоположно действию магнитного диполя.
Взяв производную в соответстви
(9)
Под действием этой ЭДС в единичном элементарном кольце возникает вихревой ток силой
(10)
где Rк = 2πrк/ds∙σ – сопротивление единичного элементарного кольца; ds — площадь его сечения, равная единице.
Ток создает в окружающем пространстве магнитное поле. В соответствии с (3) в точке, где расположена измерительная катушка, напряженность этого поля
(11)
где Lп — расстояние от точек элементарного кольца до центра измерительной катушки; Sк – площадь элементарного кольца; Мк – магнитный момент кольца.
Величина вторичного магнитного потока, пронизывающего витки измерительной катушки,
Фвт=Вsпnп=μаHsпnп.
Подставив в (12) выражение (11), получим
(13)
Изменение вторичного магнитного потока во времени создает в измерительной катушке ЭДС
e’= –dФвт/dτ=iωФВт.
Окончательная формула для ЭДС в измерительной катушке с учетом (13) и соотношения ω=2πƒ имеет вид
(15)
Одновременно с ЭДС изучаемой среды в измерительной катушке генерируется ЭДС прямого поля генераторной катушки. Поскольку величина ЭДС прямого поля не связана с параметрами среды, то она исключается с помощью специального устройства.
Учтя соотношение μa = μ0μ и умножив правую часть равенства (15) на Lи/2 и 2/Lи, представим выражение для е' в следующем виде:
где
- коэффициент индукционного зонда;
-пространственный (геометрический) фактор элементарного кольца.
Введение множителя Lи/2 в (18) позволяет свести интегральное значение пространственного фактора всего однородного пространства к безразмерной величине, равной единице.
Коэффициент зонда Ки согласно формуле (17) определяется его конструкцией – длиной, числом и площадью витков генераторной и измерительной катушек; силой и частотой тока, питающего генераторную катушку. Знак минус в (17) свидетельствует о том, что активная составляющая э.д.с, индуцируемая вторичным магнитным полем в измерительной катушке, находится в противофазе с током, питающим генераторную катушку.
Из формулы (16) следует, что величина э.д.с, создаваемая элементарным кольцом породы при Ки=const, прямо пропорциональна удельной электропроводности среды σ, квадрату относительной магнитной проницаемости μ2 и пространственному фактору Вк.
1.2 Пространственный фактор элементарного кольца
В теории Долля введение понятия пространственного фактора элементарного кольца занимает одно из центральных мест. Рассмотрим его физический и геометрический смысл.
Выразим расстояния от центров генераторной и измерительной катушек до оси элементарного кольца через цилиндрические координаты r и z.Из рисунке 2 следует
(19)
Выражение (18) с учетом соотношений (19) перепишем в виде
(20)
Из (20) следует, что значение пространственного фактора элементарного кольца определяется его вертикальным расположением относительно катушек и горизонтальным расстоянием от оси зонда при фиксированной длине зонда.
Выясним геометрический смысл пространственного фактора Вк.
Для треугольника ABC на рисунке 2 на основании теоремы синусов запишем
sin γ/Lи = sinα/Lп.
Из прямоугольного треугольника ADC следует
Подставив (22) в (21), найдем
На основании (18) и (23)
Вк=sin3γ/2Lи2
Из (24) следует, что пространственный фактор полностью определяется величиной угла γ между сторонами Lг и Lп, под которыми виден зонд из точек элементарного кольца. В этом состоит геометрический смысл пространственного фактора. Следовательно, все элементарные кольца, из точек которых виден зонд под одним и тем же углом, имеют одинаковые пространственные факторы. Геометрическим местом сечений элементарных колец с равным пространственными факторами являются окружности, проходящие через центры генераторной и измерительной катушек. Значения пространственного фактора изменяются от 0 до 1. Максимальную величину он имеет в случаях, когда sinγ=1, т.е. угол γ = 90°. Сечения этих колец лежат на окружности с диаметром, равным длине зонда (рисунок 2).
Представление пространственного фактора элементарным кольцом позволяет получить сигнал от различных участков среды. В случае однородной среды величины э.д.с, создаваемые ее участками, зависят только от их пространственных факторов. Физический смысл пространственного фактора элементарного кольца заключается в том, что он определяет ту долю сигнала на выходе измерительной катушки зонда, которую вносят различные участки изучаемой среды. Выражения для ЭДС вторичного магнитного поля в однородном и неоднородном пространствах. Выражение (16) определяет величину сигнала, создаваемую в измерительной катушке одним элементарным кольцом породы. Электродвижущая сила всего исследуемого пространства Е будет равна сумме ЭДС элементарных колец породы, т.е.
(25)
где ds – площадь сечения кольца.
1,2,3–сечения колец, из которых виден зонд под углом γ, π-γ и π/2 соответственно.
Рисунок 2- Геометрические места сечений элементарных колец с Одинаковыми пространственными факторами при | sinγ |=| sin(π-γ) | (а) и sinπ/2 =1 (б).
Удельная электропроводность в общем случае может быть непрерывной функцией координат точек пространства.
В случае однородной немагнитной среды (µ =1)
(26)
ибо радиусы элементарных колец изменяются от 0 до ∞, а координата z – от –∞ до +∞.
Существует строгое доказательство соотношения
т. е. суммарный пространственный фактор однородной изотропной среды равен единице. В этом случае отношение Е/КИ равно на основании (92) истинно удельной электропроводности среды
(27)
Для неоднородной немагнитной среды с коаксиально-цилиндрическими и плоско-параллельными горизонтальными поверхностями раздела, имеющей участки с электропроводностями пласта σпл, зоны проникновения σзп скважины σр и вмещающих пород, полная ЭДС, будет равна сумме сигналов от каждого участка среды. Вклад каждого участка в суммарный сигнал равен произведению его электропроводности на пространственный фактор. Полная ЭДС в такой неоднородной среде
Е' = Ки(σплВи + σзпВзп+σрВр + σвмВвм), (28)
где
– пространственные факторы соответственно пласта,
зоны проникновения, промывочной жидкости и вмещающих пород.
Пространственные факторы каждого участка среды представляют собой сумму пространственных факторов Вк,п Вк,зп Вк,р, Вк,вм элементарных колец по их площадям сечения sп, sр, sвм.
По аналогии с (27) и на основании (28) отношение E'/Kп в случае неоднородной среды есть ее эффективная удельная электропроводность
Измеряемая эффективная электропроводность в общем случае отличается от истинной удельной электропроводности пласта вследствие искажающего влияния на величину сигнала параметров скважины, зоны проникновения, вмещающих пород и соседних пластов. Это отличие тем существеннее, чем неоднороднее исследуемая среда.
Из (29) следует, что удельная электропроводность участка среды и его пространственный фактор равнозначны при формировании полного сигнала. Участок среды с высокой электропроводностью, но малым пространственным фактором может создать такой же сигнал, как и участок среды с низкой электропроводностью, но большим пространственным фактором. Чем выше электропроводность отдельных участков среды, тем значительнее величина сигнала от них; при весьма низкой электропроводности какого-либо участка среды (σ→0) его вклада в полный сигнал практически не будет.
Из сказанного выше следует, что кривая индукционного метода наиболее дифференцирована против пород высокой электропроводности и сглажена против пород с низкой электропроводностью.
2.Техника и методика работ
Обычный низкочастотный индукционный метод с продольным датчиком
Обычный низкочастотный индукционный метод основан на изучении электромагнитного поля продольного (вертикального) датчика, ось которого совпадает с осью скважины. В этом случае вихревые токи, индуцированные первичным полем, расположены в плоскостях, перпендикулярных к оси скважины, и не пересекают поверхностей раздела горизонтальных слоев.
На результаты измерении простейшим двухкатушечным зондом с целью определения истинного удельного сопротивления пласта значительное искажающее влияние оказывают скважина, зона проникновения и вмещающие породы, а также прямой сигнал от генераторной катушки. Для уменьшения влияния указанных факторов и исключении прямого сигнала. Х. Г. Долль предложил многокатушечные фокусирующие зонды.
2.1 Зонды обычного низкочастотного индукционного метода
Многокатушечный зонд представляет собой систему катушек, укрепленные на одном изоляционном стержне (рисунок 1). Генераторная ГК и измерительная ПК катушки являются основными, остальные называются компенсационными К и фокусирующими Ф(ФМ) в генераторной ФГ и измерительной ФП цепях. Компенсационные катушки служат для исключения в в измерительной катушке ЭДС прямого поля, индуцируемого генераторной катушкой.
Фокусирующие катушки предназначены для уменьшения влияния неполезных сигналов путем создания в приемной цепи э,д.с., обратных по знаку э.д.с , вызванных вихревыми токами, которые циркулируют в скважине, зоне проникновения и вмещающих породах. Фокусирующее действие катушек достигается путем подбора числа их витков, расположения и включения их относительно главных катушек. Число дополнительных катушек, их взаимное положение и число витков должны быть такими, чтобы в значительной мере было исключено влияние скважины, зоны проникновения и вмещающих пород, а измеряемая эффективная электропроводность была как можно ближе к истинному значению электропроводности пласта. Компенсационные и фокусирующие катушки включаются последовательно с главными, но их витки намотаны противоположно виткам генераторной и измерительной катушек.