CC BY
43
_ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА_
2017 Геология Том 16, № 2
ГЕОФИЗИКА, ГЕОФИЗИЧЕСКИЕ МЕТОДЫ ПОИСКОВ ПОЛЕЗНЫХ ИСКОПАЕМЫХ
УДК 550.832.75
Разработка технологии индукционного каротажа в нестационарном режиме при бурении субгоризонтальных участков стволов нефтяных и газовых скважин
В.К. Теплухина, А.Н. Ратушнякь, В.И. Костицынс, Ван Сяолуна
аООО «ПетроТул-НБ», 452606, Республика Башкортостан, Октябрьский, ул. Садовое Кольцо, 4 А. E-mail: [email protected] ьИнститут геофизики УрО РАН, 620016, Екатеринбург, ул. Амундсена, 100. E-mail: [email protected]
сПермский государственный национальный исследовательский университет, 614990, Пермь, ул. Букирева, 15. E-mail: [email protected] компания L.L.C. «Leader», 215211, Китай, Сучжоу, район Уцзян, проспект Фэньху, 558.E-mail: [email protected]
(Статья поступила в редакцию 20 мая 2017 г. )
Рассмотрены результаты разработки технологии модуля индукционного комплекса (МИР), предназначенного для измерения сопротивления пород в процессе бурения горизонтальных и наклонно-направленных скважин на нефть и газ. Представлены материалы производственно-технологического опробования при проведении эксплуатационного бурения промысловых скважин на территории Республики Башкортостан. Выполнен сопоставительный анализ данных каротажа сопротивления, полученных после завершения бурения, и результатов, зафиксированных комплексом МИР в составе LWD с телесистемой «ТАРГЕТ» в процессе бурения. Приведены технические характеристики параметров комплекса МИР.
Ключевые слова: сопротивление пород, горизонтальные и наклонные скважины, месторождения нефти и газа, модуль индукционного комплекса (МИР), технология, LWD, телесистема «ТАРГЕТ». DOI: 10.17072/psu.geol.16.2.122
Введение
Одним из наиболее эффективных методов формирования оптимальной системы разработки является разбуривание нефтяных и газовых месторождений горизонтальными и многоствольными наклонно-направленными скважинами. Это приводит к увеличению площади фильтрации и в значительной степени по-
вышает эффективность разработки низкопроницаемых коллекторов.
Исследования скважин в процессе бурения LWD (LoggingWhileDrilling) в значительной степени позволяют оптимизировать время на анализ геологической информации в связи с существенным уменьшением зоны проникновения фильтрата бурового раствора в структуру нефтяного или газового коллектора, что
© Теплухин В.К., Ратушняк А.Н., Костицын В.И., Сяолун Ван , 2017
122Щ1Ц|Н|
позволяет сократить время его освоения и, это особенно актуально при разработке пластов малой мощности, вести процесс геонавигации траектории ствола скважины в соответствии с морфологией пласта.
В компании «ПетроТул» в 2015-2016 гг. разработана специализированная технология модуля индукционного резистивиметра МИР, позволяющая проводить детальные исследования удельного сопротивления пород в процессе бурения в режиме on-line, основанная на изучении характеристик комплекса электромагнитных полей в нестационарном режиме.
Теоретическое обоснование
Применение технологии импульсного изучения электромагнитного поля вместо изучения гармонического сигнала (Антонов, 1980; Бурсиан, 1972; Ваньян, 1965; Кауфман, Соколов, 1972 и др.) продиктовано результатами детальных физико-теоретических исследований авторов статьи и большими потенциальными возможностями нестационарного электромагнитного поля в прикладном применении, указанными еще в ранней работе М.И. Плюснина (1968).
Для измерений переменного магнитного поля применяют многовитковые катушки, ориентированные по оси скважины, в которых индуктивно наводимая электродвижущая сила (ЭДС) связана с магнитным полем законом индукции Фа-радея.
Определение ЭДС в контуре через циркуляцию напряженности электрического поля S (w) = <f Etdl по формуле связи через вектор-потенциал E = rot^A* в
цилиндрической системе координат (R, ф, z) дает выражение, учитывающее размер измерительного контура радиуса а:
S(w) = f Etdl = -
Mza sin0(1 + kR)e-kR, 2R
sin# = a =
а
R Л
a2 + L2
(1)
где Ег - напряженность электрического поля, Мz - магнитный момент, А - вектор-потенциал магнитного типа; н - циклическая (круговая) частота тока (2л/), 1/с; /Л0 - магнитная постоянная = 4я,10"7 Гн/м; к - волновое число среды; Ь - размер установки (расстояние между генераторным и приемным контурами по оси 7).
Определение ЭДС через изменение потока (йФ) индукции магнитного во времени (й) приводит к аналогичному выражению, если измерительный контур имеет малый радиус, величина которого значительно меньше разноса: а < Ь, бШ^ « a/L,
dФ
£(w)=--ф = -ja/doSH Hz =
dt
M
(2)
= -jaMoS„-MzT (1 + kL)e-kL
2 лL
где 8н - эффективная площадь соленоида приемной катушки; Hz - вертикальная составляющая напряженности магнитного поля. Откуда следует, что квадратурные компоненты R £ = а -а- J H и
e m z
Jm £ = Ь -а- ЯН, а, Ь - сопб1.
т е z 5 '
На основании зависимости (2) произведены расчеты реальной компоненты поля в зависимости от значений удельного сопротивления среды (рис. 1). Представленные графики наглядно демонстрируют частотные зависимости ЭДС в проводящей среде, где поле является практически реальной компонентой.
В непроводящей среде ЭДС является мнимой величиной
£(м;)=-]аМоМ4г = -jаJK,
где K = д
S Г SИ
2лЙ
2nL
- коэффициент взаим-
ной индукции измерительной и генераторной катушек; и 8н - площади генераторного и измерительного контуров, например, для контуров а = г = 5 см при Ь = 1 м коэффициент К « 10-8 Гн.
Чтобы избежать частотной зависимости измеряемой ЭДС, в измерительной схеме применяют преобразователь «ток-напряжение» с индуктивностью (измери-
тельной катушкой) в цепи обратной связи. Это позволяет регистрировать величину ЭДС, прямо пропорциональную напряженности магнитного поля с учетом сдвига фазы ЭДС относительно магнитного поля на угол л/2.
Рис. 1. Частотные зависимости Яе£ в
контуре малого радиуса. Шифр кривых -сопротивление среды (Омм)
Взаимоиндукцию измерительной и генераторной катушек устраняют при фазировании измерительного канала перед скважинными измерениями магнитного поля, добиваясь величины Jm£ = 0 в непроводящей среде (в воздухе).
При импульсном режиме изменения тока выражение для измеряемой нестационарной величины ЭДС в контуре малого радиуса а получим с помощью обратного преобразования Лапласа-Карсона ф^кт, Ргиёшкоу, 1965).
Формула для нестационарной ЭДС примет вид
e(t) = -р sin в • u se-272
422ЛR4
(3)
=- pj
S r S И
ж
42ПR5
• s(u ),
где р - удельное электрическое сопротивление, I - сила тока, питающего генераторный контур, u - безразмерный параметр становления, Sr и Sh - эффективная площадь генераторного и измерительного
контуров соответственно, s{u) = u5e~"2'2.
На полярной оси диполя при размерах измерительного контура а < L и R « L, sin 6 ~ a / L в проводящей среде (уста-
новка для каротажа скважин) S(t) = -р- J • dк •£(u),
S rS и-
где d =
(4)
коэффициент уста-
новки или постоянная прибора (Doll 1949; Anderson, Berber, Leveridge. 2008).
Временная характеристика ЭДС определяется только поведением безразмерной функции - электрического
числа s{u). Ее анализ = 0 ^ показывает, что максимум электрического числа
£(и ) = 4,58869 соответствует параметру u = л/5 = 2,235 (или параметра
и2 осч
становления поля p = — = 2,5 ). При этом £ (и )• u = 10,2606 = ^nst. На
max max
рис. 2 представлена временная зависимость ЭДС S(t) от диполя с моментом
M* = 1 Ам2 на разносе L = 1 м в витке радиусом а = 0,01 м.
Рис. 2. Временная зависимость нестацио-
/V
нарнои ЭДС S(t) = pJdKS(u) от р среды в измерительном контуре малого размера (шифр кривых, Омм)
Величину удельного сопротивления среды р по данным индукционного каротажа можно определять двумя способами:
- по времени t, соответствующему экстремуму измеряемой ЭДС,
p = ß0L2/10/1 = 125,66 • L2/1; (5)
- по величине экстремума измеряемой ЭДС
£ (t)
-ext^ = к £ (t),
J • d • 4,58869 extryj
где К =
1
J - - 4,58869
зонда,
¿„ =
2жЬ5
коэффициент коэффициент
установки.
Важнейшим направлением при исследованиях стандартным индукционным каротажом является изучение влияния значений проводимости среды на величину измерения ЭДС, оценка разрешающей способности технологии по значениям проводимости и определение максимального значения сопротивления среды, при котором возможны измерения.
Модели исследования с модулем индукционного резистивиметра (МИР)
Основными моделями для оценки эффективности зондовой системы резисти-виметра МИР являлись:
1.Стенд с концентрически расположенными проводниками радиусом до 0,9 м и с возможностью изменения сопротивления каждого контура в диапазоне 1-400 Омм;
2. Стенд с концентрически расположенными проводниками радиусом до 2,5 м и с изменением сопротивления каждого контура в диапазоне от 1 до 400 Омм;
3. Объемная «большая» модель емкостью 5 м3 с изменением сопротивления электролита в диапазоне от 7 до 200 Омм и измерением на оси модели;
4. Объемная «малая» модель емкостью 0,75 м3 с изменением сопротивления электролита в диапазоне 7-200 Ом м и измерением на оси модели;
5. Модельная (тестировочная) скважина в четвертичных отложениях;
6. Рабочие (производственные) наклонно-направленные скважины.
Типовые зависимости значения ЭДС на приемном зонде, полученные при наличии одного замкнутого контура с изменением сопротивления в диапазоне 0-400 Омм при условии, что контур находится в безграничной непроводящей
среде (отсутствуют дополнительные осложняющие факторы), представлен на
рис. 3.
Рис. 3. Зависимость ЭДС от значения сопротивления проводящего контура (шифр кривых) на расстоянии 0,925 м в проводящей среде (14,97 Омм). Длина зонда Ь = 1,8 м
Исследования с резистивиметром МИР могут проводиться в диапазоне рабочих температур от 0 до 120°С при изменении рабочего давления от 0 до 60 мПа.
В табл. 1 приведены сведения о точности измерений удельного электрического сопротивления при различных радиусах исследования в диапазоне от 0 до 300 Ом.м. В табл. 2 представлены данные о глубинности исследований (расстоянию от оси скважины) и вертикальном разрешении при удельном электрическом сопротивлении от 1 до 10 Омм и в диапазоне от 10 до 50 Омм.
С модулем индукционного резистиви-метра МИР проведены испытания на системах двухмерных и трехмерных (объемных) моделей. Были исследованы возможности модуля для решения задач геонавигации, приближения и удаления относительно границ со средами с разной электрической проводимостью.
Результаты исследований представлены на рис. 4. Как видим, форма основного сигнала сохраняется при перемещении прибора от центра сред к границе. При соосном положении осей прибора и емкостей новый экстремум практически не
проявляется (вариант 1). При смещении оси прибора в сторону границы сред на величину 1/2 радиуса малой емкости новый экстремум становится сопоставим по амплитуде с основным (вариант 2), а при дальнейшем приближении существенно превышает основной максимум (вариант 3).
Таблица 1. Сведения о точности измерений удельного электрического сопротивления при различных радиусах исследований_
Радиус исследования Диапазон Точность, %
0,75-1,0 м 0-200 Омм 5
1,5-2,0 м 0-200 Омм 5
Таблица 2. Глубина исследования и верти-
кальное разрешение
Пара- Глубина исследова- Верти-
метр ния, см каль-
измере- малый боль- ное
ния радиус шой ра- разре-
диус шение,
см
Удельное сопротивление среды 1- 10 Омм
Ампли- 50-75 150-200 1
туда
Удельное сопротивление среды 10-50Омм
Ампли- 60-100 160-220 1
туда
Отсюда можно сделать вывод, что при приближении к вертикальной границе сред (варианты 2 и 3) на времени 0,2-0,4 микросекунды появляется новый экстремум, амплитуда которого растет с приближением к субвертикальной радиальной границе сред с различными значениями сопротивления. Указанное обстоятельство открывает перспективное направление для проведения геонавигации - управление траекторией бурения в маломощных коллекторах.
Опытно-методические работы с рези-стивиметром МИР
Бурение горизонтального участка в скв. 3125 гс (Бузовьязовское месторождение, Башкортостан) было проведено с
глубины 2 316 м до глубины 2 412 м (проектный забой) за 2 суток с регистрацией необходимых технологических параметров и данных измерений удельного электрического сопротивления в радиусах 0,75-1,0 м и в диапазоне глубин 2 268-2 412 м c измерением параметров в радиусах 1,5-2,0 м.
Бурение аналогичного участка в скв. 210 гс1 (Шкаповское месторождение, Башкортостан) проведено с глубины 2 152,2 м до 2 343 м по стволу также с регистрацией всех технологических параметров, в том числе и данных электрического сопротивления со стробированием по глубине 0,2 м.
Компоновка низа бурильной колонны включала следующее оборудование: долото, винтовой забойный двигатель, модули гамма-каротажа и инклинометра, разделитель, модуль резистивиметра МИР, комплекс телесистемы с электромагнитным каналом связи «TARGET» (ООО «ПетроТул-НБ»). При проводке скв. 210 гс1 в компоновку бурильной колонны был дополнительно включен модуль резистивиметра АРБ(Канада) для контроля при сопоставлении измеренных данных сопротивления.
Проведение исследований скважин-ным резистивиметром МИР непосредственно в процессе бурения позволило получить значения сопротивления в условиях незначительного проникновения фильтрата бурового раствора. По данным бурения и резистивиметра МИР горизонтальные участки стволов скважины вскрыли относительно однородные коллекторы, представленные тонкозернистыми светло-серыми и коричнево-серыми известняками.
При сопоставлении данных МИР с результатами каротажа ГИС (ОАО «Баш-нефтегеофизика») по скв. 3125 гс и данными резистивиметра APS по скв. 210 гс1 установлено, что показания резистиви-метров APS и МИР практически совпадают (рис. 5). Расхождение данных находится в пределах допустимой погрешности измерений (5 %).
/ \ Условия проведения эксперимента О о • Емкость 5 м3 Емкость 0.75 м3 Прибор
Вариант 1 [ 0 ) Емкости и прибор отцентрированы \ / ---------- 11111111111111111 IN 1111111111 и 1111 ЭДС, MB
i /
Жо 2.0 V i nil nil ми IiiiiIii i
ХТХ Вариант 2 I 1 \ Емкости отцентрированы, прибор смещен 1 [ 0 ] —— ...................'"i»»(""j'r ЭДС, MB i T, мкс
......Цо'""...... ' —- % r
Вариант 3 / \ Емкости отцентрированы, прибор смещен / [ 1 на радиус к стенке малой емкости 1 1 -■ к II ii ii 1111 II ii 11 II ■ 11111 > 111111111 ii f;s................ iFTT T, мкс
mi'IiimIiiiiIiiiiIi'iiIiiiiIiiiiI'I" -2.0 IV .............'"" niliiiiliiiiliiiiiiii
Рис. 4. Вид в плане и результаты эксперимента по изучению степени удаления от границ сред с различными значениями удельного электрического сопротивления
Рис. 5. Сопоставлении данных МИР («ПетроТул-НБ») и резистивиметра APS (Канада) по скв. 210 гс1 (Шкапово, Башкортостан)
Выводы
1. Модуль резистивиметра МИР технологически совместим с телесистемой «TARGET», подключен к стандартной шине передачи данных RS-485. Конструкция модуля МИР допускает свободный допуск к блокам памяти, позволяет скачивать информацию по значениям сопротивления с необходимой детальностью стробирования по глубине.
2. Резистивиметр МИР может работать во всех типах бурового раствора, включая соленасыщенные и растворы на нефти (Ратушняк, Байдиков, Теплухин, 2016). Значения удельного сопротивления пород для управления траекторией бурения предоставляются в режиме реального времени.
3. Резистивиметр МИР практически может быть адаптирован к работе в составе телесистем с гидравлическим каналом связи.
4. Результаты резистивиметрии МИР и дальнейшей интерпретации доступны специалистам еще во время бурения в режиме реального времени, что позволяет оперативно реагировать на изменение геологической обстановки, уточнять в комплексе с гамма-каротажом структурные элементы пласта коллектора и оперативно проводить целевое ориентированное бурение (геонавигация) в процессе проводки ствола скважины.
Библиографический список
Антонов Ю.Н. Изопараметрическое каротажное зондирование (обоснование ВИКИЗ) //
Геология и геофизика, 1980. № 6. С. 81-91.
Бурсиан В.Р. Теория электромагнитных полей, применяемых в электроразведке. Л.: Недра, 1972. 245 с.
Ваньян Л.Л. Основы электромагнитных зондирований. М.: Недра, 1965. 109 с.
Кауфман А.А., Соколов В.П.Теория индукционного каротажа методом переходных процессов. Новосибирск: Наука, 1972. 108 с.
Кудрявцев Ю.И. Некоторые вопросы теории индукционного каротажа // Прикладная геофизика, 1960. Вып. 28. С. 101-115.
Никитина В.Н. Общее решение осесиммет-ричной задачи теории индукционного каротажа // Известия АН СССР. Серия геофизическая. 1960. № 4. С. 607-616.
Плюснин М.И. Индукционный каротаж. М.: Недра, 1968. 142 с.
Плюснин М.И., Вильге Б.И. Обоснование индукционного каротажа методом переходных процессов // Изв. вузов. Геология и разведка. 1969. № 5. С. 158-165.
Диткин В.А., Прудников А.П. Справочник по операционному исчислению. М.: Высшая школа, 1965. 466 с.
Ратушняк А.Н., Байдиков С.В., Теплухин В.К. Индукционный каротаж скважин с учетом влияния проводящей промывочной жидкости // Уральский геофизический вестник. 2016. № 1 (27). С. 40-47.
Doll H. G. Introduction to Induction Loggingand Application to Logging of Wells Drilled With Oil-Based Mud // Journal of Petroleum Technology. 1949. Vol. 1, № 6. P. 148-162.
Андерсон Б., Бербер Т., Леверидж Р. Трехмерный индукционный каротаж: старые измерения под новым углом // Нефтегазовое обозрение. 2008.
URL : http: //www .slb.com/~/media/File s/resour-ces/oilfield_review/russia08/sum08/05_triaxi alinduction.pdf.
Development of Induction Logging Technology for Non-Stationary Conditions While Drilling Subhorozontal Sections of Oil and Gas Wells
V.K. Teplukhina, A.N. Ratushnyakb, V.I. Kostitsync, Wang Xiaolongd
aL.L.C. «PetroTul-NB», 4A Sadovoye Koltso, Oktyabrskiy 452606, Bashkortostan Republic, Russia. E-mail: [email protected]
bInstitute of Geophysics, Ural Branch of Russian Academy of Science, 100 Amundsen Str, Ekaterinburg 620016, Russia. E-mail: [email protected] cPerm State University, 15 Bukireva Str, Perm 614990, Russia. E-mail: [email protected]
dL.L.C. Leader, 558 Avenue Fenhu, Wujiang District, Builing, Suzhou 215211, China
The article presents the results of development of the induction logging module designed for resistivity measurements while drilling horizontal and directional wells for oil and gas. Materials of the tool testing during the drilling of wells in the oil field in the Republic of Bashkortostan are shown. We introduce the comparative analysis of resistivity data obtained after completion of the drilling, and the results obtained with developed complex as part of the LWD telemetry system TARGET during the drilling. The technical characteristics of complex MIR are presented.
Key words: rock resistivity; drilling; horizontal and directional wells; oil and gas fields; induction logging complex (MIR); LWD; telemetric system TARGET.
References
Antonov Y.N. 1980. Izoparametricheskoe karota-zhnoe zondirovanie (obosnovanie VIKIZ) [Isoparametric induction logging]. Geologiya i geofizika. 6:81-91. (in Russian) Bursian V.R. 1972. Teoriya elektromagnitnykh poley, primenyaemykh v elektrorazvedke [Theory of electromagnetic fields used in electrical exploration methods]. Leningrad, Nedra, p. 245. (in Russian) Vanyan L.L. 1965. Osnovy elektromagnitnykh zondirovaniy [Basics of electromagnetic sounding]. Moskva, Nedra, p. 109. (in Russian)
Kaufman A.A., Sokolov V.P. 1972. Teoriya in-duktsionnogo karotazha metodom perekhod-nykh processov [Theory of transient induction logging]. Novosibirsk, Nauka, p. 108. (in Russian)
Kudriavzev Y.I. 1960. Nekotorie voprosy teorii induktsionnogo karotazha [Some aspects of theory of induction logging]. Priklandnaya geofizika. 28:101-115. (in Russian) Nikitina V.N. 1960. Obshchee reshenie osesi-metrichnoy zadachi teorii induktsionnogo ka-rotazha [General solution for axisymmetric problem of theory of induction logging]. Izvestiya AN SSSR. Seriya geofizicheskaya. 4:607-616. (in Russian)
Plyusnin M.I. 1968. Induktsionnyy karotazh [Induction logging]. Moskva, Nedra, p. 142. (in Russian)
Pliusnin M.I., Vilge B.I. 1969. Obosnovanie in-duktsionnogo karotazha metodom perekhod-nykh protsessov [Verification of transient induction logging]. Izv. vuzov. Geologiya i razvedka. 5:158-165. (in Russian) Ditkin V.A., Prudnikov A.P. 1965. Spravochnik po operatsionnomu ischisleniyu [Reference book for operator calculus]. Moskva, Vysshaya shkola, p. 466. (in Russian) Ratushnyak A.N., Baydikov S.V., Teplukhin V.K. 2016. Induktsionnyy karotazh skvazhin s uchetom vliyaniya provodyashchey promyvochnoy zhidkosti [Induction logging with consideration of influence of the conductive mud]. Uralskiy geofizicheskiy vest-nik. 1(27):40-47. (in Russian) Doll H.G. 1949. Introduction to Induction Logging and Application to Logging of Wells Drilled With Oil-Based Mud. Journal of Petroleum Technology. 1(6): 148-162. Anderson B., Berber T., Leveridge R. 2008. Trekhmernyy induktsionnyy karotazh: starye izmereniya pod novym uglom [3D Induction logging: old measurements by new point of view]. Neftegazovoe obozrenie. URL: http://www.slb.com/~7media/Files/resour-ces/oilfield_review/russia08/sum08/05_triaxi alinduction.pdf. (in Russian)
