«Ембімұнайгаз» АҚ кен орындарында көлденең ұңғымалары құрылысында геомеханиканы қолдану

Толық мәтін

Введение

Геомеханическое моделирование для оптимизации бурения горизонтальных скважин (далее – ГС) в ТОО «КМГ Инжиниринг» начали применять с 2019 г. на месторождениях АО «Эмбамунайгаз» (далее – ЭМГ) после начала возникновения осложнений при переходе к бурению ГС.

На сегодняшний день проведены работы по 1D геомеханическому моделированию для 12 ГС, входящих в план сопровождения скважин ЭМГ на месторождениях Гран, Актобе, Вост. Молдабек, Жанаталап, С. Нуржанов, С. Балгимбаев, Юго-Зап. Камышитовое, Забурунье. Выполнена разработка рекомендаций по планированию безаварийного бурения, включающего следующие виды работ:

• построение опорных 1D геомеханических моделей;
• сокращение рисков при бурении;
• выработка рекомендаций.
• Основные задачи проведенных работ:
• разработка рекомендаций по планированию безаварийного бурения;
• сбор, контроль качества и анализ имеющихся геолого-геофизических и технологических данных;
• формирование рабочего цифрового проекта;
• построение модели механических свойств геологической среды (опорных и проектной);
• расчет устойчивости ствола скважин (опорных и проектной);
• предоставление необходимого набора параметров для дальнейшего моделирования дизайна гидроразрыва пластов (далее – ГРП) при необходимости.

Для создания качественной и реалистичной модели механических свойств необходим достаточно большой объем данных по геологическому строению месторождения, сейсморазведке, геофизическим исследованиям скважин, анализам керна, гидродинамическим исследованиям и др. Рассматриваемые месторождения отмечаются хорошей изученностью и наличием необходимой информации, за исключением геомеханических анализов керна и широкополосного акустического каротажа (далее – АКШ).

Построение моделей опорных скважин

Анализ бурения скважин

Процесс построения геомеханической модели традиционно начинается с проведения аудита имеющихся геолого-геофизических данных и материалов по бурению опорных скважин. Для калибровки геомеханической модели необходимо проведение анализа данных геолого-технологических исследований (далее – ГТИ) и инцидентов при бурении. Наиболее важными являются следующие параметры:

1. плотность бурового раствора на входе и выходе;
2. расход бурового раствора на входе и выходе;
3. общий объем емкостной системы;
4. скорость спускоподъемных операций;
5. механическая скорость проходки;
6. газопоказания;
7. вес на крюке.

В рамках аудита были проанализированы суточные буровые рапорты и финальные отчеты по бурению на предмет наличия инцидентов, возникших при бурении. На следующем планшете в качестве примера представлена сводная информация о конструкции скв. 108 м. Актобе, стратиграфической разбивке, используемом типе бурового раствора и плотности, а также график «глубина – день» с отображенными осложнениями как по опорным скважинам, так и по ГС (рис. 1).

 

Рисунок 1. Анализ бурения скв. 108, м. Актобе

 

Анализ данных геофизических исследований, расчет горного давления и упруго-прочностных характеристик

Анализ данных геофизических исследований является неотъемлемой частью построения 1D геомеханических моделей. На большинстве изучаемых месторождений ЭМГ присутствуют запись продольной акустической волны и каротаж объёмной плотности пород, однако отсутствует запись поперечной акустической волны, необходимой для расчета упруго-прочностных параметров.

Для восстановления поперечной волны (далее – DTS) была использована зависимость, полученная по м. С. Нуржанов, которое относится к Эмбенско-Прорвинской группе. В анализе участвовало 3 скв. (337, 538, 540). После восстановления DTS проводился расчет упругих параметров. Калибровка динамических параметров в статические и расчет прочностных характеристик породы производились по результатам керновых исследований. Зависимость продольной акустической волны от поперечной представлена на рис. 2.

 

Рисунок 2. Зависимость продольной акустической волны от поперечной

 

Литостатическое (горное, вертикальное) давление среды создается весом вышележащих пород. Оно рассчитано на основе данных плотностного каротажа в интервалах его наличия и синтетической кривой плотности, соответствующей тренду нормального уплотнения в интервалах отсутствия значений. На рис. 3 представлены примеры расчета горного давления на скв. 116 и 120 м. Актобе.

 

Рисунок 3. Расчет горного давления, м. Актобе

 

Как было отмечено ранее, для расчета упруго-прочностных свойств, таких как коэффициент Пуассона, модуль Юнга, прочность при одноосном сжатии, прочность при растяжении, угол внутреннего трения, используются каротажные, а именно: кривые акустических скоростей пробега, кривая объёмной плотности и гамма-каротаж. На рис. 4 отражен пример полученных упруго-прочностных свойств по исследуемым скважинам.

 

Рисунок 4. Упруго-прочностные свойства, м. Актобе

 

Калибровка минимального горизонтального напряжения

Основную сложность при построении геомеханических моделей представляло отсутствие прямых замеров напряжения, таких как гидравлические тесты extended leak-off test (далее – ELOT) и ГРП в исследуемых интервалах и пластах.

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

• анализ данных фрак-листа по месторождениям ЭМГ;
• отбраковка нерелевантных значений давления смыкания трещины ГРП и пластового давления;
• корректировка пластового давления на основании базы данных гидродинамических исследований скважин (далее – ГДИС);
• построение зависимостей давления смыкания трещины ГРП от пластового давления по всем месторождениям;
• оценка общего тренда;
• калибровка тектонических констант в триасовом и нижнеюрском комплексе;
• использование полученных коэффициентов для расчета минимального и максимального напряжений по рассматриваемым месторождениям.

Рис. 5 отражает полученную общую зависимость по месторождениям ЭМГ. Как видно, все результаты анализа укладываются в единый тренд, что свидетельствует об общей тектонической нагрузке региона.

 

Рисунок 5. Общая зависимость давления смыкания трещины ГРП от пластового давления на месторождениях ЭМГ

 

Для калибровки тектонических констант была выбрана скв. 538 м. С. Нуржанов. По данной скважине имеются запись широкополосной акустики и замеры ГДИС, и проведено ГРП с оценкой смыкания трещины. Также данная скважина лежит на полученном тренде по Прорвинской группе (рис. 6). Оцененное пластовое давление составило 242,3 атм на абсолютной глубине 3080 м, забойное давление смыкания трещины ГРП – 410 атм и 389 атм по проведенному замещению и миниГРП соответственно.

 

Рисунок 6. Зависимость давления смыкания трещины ГРП от пластового давления для Прорвинской группы месторождений

 

По результатам калибровки (рис. 7) получены следующие тектонические коэффициенты: коэффициент минимального горизонтального напряжения εℎ = 0,0001, коэффициент максимального горизонтального напряжения εΗ = 0,00023, коэффициент Биот = 0,9. Стоит отметить, что полученные коэффициенты совпадают с тектоническими коэффициентами на м. Кожасай, что свидетельствует об общем тектоническом региональном тренде.

 

Рисунок 7. Калибровка минимального горизонтального напряжения по скв. 538, м. С. Нуржанов

 

Расчет устойчивости скважин

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

Ниже представлены примеры расчета геомеханических моделей для вертикальных и ГС. На скв. 120 м. Актобе при проработке интервала 2277–2361 м отмечались осыпи, представленные плотными глинистыми породами. После обработки бурового раствора и увеличения плотности с 1,25 г/см³ (красная пунктирная линия) до 1,30 г/см³ (зеленая линия) бурение транспортной секции закончено без осложнений. Данный факт отлично согласуется с полученной геомеханической моделью и расчетным градиентом обвалообразования, который увеличивается с 1,20 г/см³ до 1,30–1,35 г/см³ с глубины 2277 м. При расчете использована прочность породы, полученная по керновым зависимостям.

Таким образом, нестабильность ствола скв. 120 м. Актобе связана с механической нестабильностью глинистой покрышки. Расчетный градиент обвалообразования составил до 1,30–1,35 г/см³ (рис. 8).

 

Рисунок 8. Расчет стабильности ствола скв. 120, м. Актобе

 

Ниже приведен расчет устойчивости скв. 2675 м. Вост. Молдабек (рис. 9). При спуске эксплуатационной колонны происходили многочисленные посадки, которые могут быть связаны с геомеханической нестабильностью пород. Расчетный градиент начала поглощения составил от 1,60 г/см³ до 1,70 г/см³, градиент обвалообразования – от 0,60 г/см³ до 1,25–1,30 г/см³. При бурении данной скважины использовалась плотность бурового раствора 1,22 г/см³, анализ каверномера указывает на отсутствие проблем со стабильностью ствола скважины, что подтверждает расчетный градиент обвалообразования.

 

Рисунок 9. Расчет стабильности ствола скв. 2675, м. Вост. Молдабек

 

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

1. Происходит сужение ствола скважины в интервалах с высокими гамма-показаниями. Данные интервалы соответствуют глинистым отложениям.
2. Происходит увеличение диаметра ствола скважины в интервалах с более низкими гамма-показаниями. Данные интервалы соответствуют песчаникам.
3. Полученная тенденция может отражать размыв слабосцементированных песчаников и образование каверн. Напротив, в интервалах более плотных глинистых отложений или глинистых песчаниках происходит образование корки или выпучивание глин.
4. В разрезе присутствуют плотные карбонатные пропластки (низкие гамма- показания) с номинальным диаметром ствола. Данные пропластки выступают своего рода ступеньками между кавернозными интервалами и могут препятствовать нормальному спуску колонн. Также возможны осложнения с хождением инструмента и спуском колонн в интервалах сужения ствола скважины.

 

Рисунок 10. Кросс-плот зависимости гамма-каротажа от диаметра скважины, скв. 2675, м. Вост. Молдабек

 

Анализ прихвата ГС 718 м. С. Нуржанов

При бурении горизонтальной секции с диаметром Ø152,4 мм ГС 718 м. С. Нуржанов произошел прихват на глубине 2060 м. По сообщению с буровой, буровой инструмент потерял подвижность, циркуляция полная, яс не активируется, возможно, прихвачен. Данные признаки указывают на дифференциальный тип прихвата, для ликвидации которого необходимо снижение статической плотности бурового раствора. Однако при снижении плотности бурового раствора повышаются риски нестабильности ствола скважины.

Расчет рисков обвалообразования при снижении плотности бурового раствора проведен на основании обновленной 1D геомеханической модели по опорной скв. 704 с учетом бурения транспортной секции ГС 718. Оценка рисков производилась с учетом глубины обрушения ствола скважины. Обрушение до 5% является безопасным, от 5% до 10% – приемлемым, свыше 10% существует высокий риск механического прихвата в силу катастрофического обрушения ствола скважины (рис. 11).

 

Рисунок 11. Пример глубины обрушения ствола скважины

 

Наибольший риск обвалообразования в горизонтальной секции соответствует глинистому интервалу 1877–1884 м. Расчетный градиент обвалообразования при нулевом обрушении составил 1,35 г/см³, при обрушении 5% – 1,18 г/см³, при обрушении 10% – 1,12 г/см³. Расчеты согласуются с тем, что плотности 1,25 г/см³ было достаточно для поддержания стабильности ствола скважины, обвального шлама при бурении скв. 718 не наблюдалось (рис. 12).

 

Рисунок 12. Обновление 1D геомеханической модели по скв. 704

 

Плотности бурового раствора 1,18 г/см³ достаточно для поддержания стабильности ствола скважины в горизонтальной секции, за исключением глинистых интервалов.

Рекомендации

Рекомендуется снизить плотность до 1,18 г/см³ для снижения репрессии на пласт и повышения эффективности ликвидации дифференциального прихвата.

Согласно геомеханическим расчетам возможно снижение плотности от 1,18 до 1,12 г/см³ с приемлемым обрушением от 5% до 10%.

Снижение ниже 1,12 г/см³ нежелательно в силу повышения рисков катастрофического обрушения ствола скважины и механического прихвата.

По результатам предоставленных рекомендаций было проведено снижение плотности бурового раствора до 1,18 г/см³, прихват ликвидирован.

Выводы

Ниже представлены основные выводы и рекомендации по результатам проведения исследований и одномерного геомеханического моделирования:

1. Благодаря качественно выполненным керновым геомеханическим исследованиям скв. 704 м. С. Нуржанов, полученные результаты используются для построения моделей опорных скважин месторождений ЭМГ.
2. По данным проведённого анализа ГРП построена общая зависимость давления смыкания трещины от пластового давления, что необходимо для калибровки минимального горизонтального напряжения. Также в результате анализа подобраны тектонические коэффициенты по скв. 538 м. С. Нуржанов, необходимые для расчетов при моделировании.
3. Наиболее вероятное направление максимального горизонтального напряжения на месторождениях ЭМГ составляет 100–150°. Наиболее благоприятный азимут бурения ГС и расположение скважин с многостадийным ГРП (далее – МГРП) 35 – 215°.
4. Построены геомеханические модели и проведен детальный анализ по скв. 116 и 120, м. Актобе, скв. 264, м. Жанаталап, скв. 2675, м. Вост. Молдабек, скв. 718, м. С. Нуржанов. Произведен расчет горного давления, упруго-прочностных свойств, построены концептуальные модели по скв. 74 м. Гран, скв. 240 м. С. Балгимбаев и др.
5. Основными причинами осложнений при бурении ГС, согласно геомеханическим моделям, являются:

• • нестабильность ствола скв. 120 м. Актобе связана с механической нестабильностью глинистой покрышки. Расчетный градиент обвалообразования составил до 1,30–1,35 г/см³;
  • причина поглощений на ГС м. Жанаталап – превышение градиента порового давления в трещиноватой карбонатной пачке К2. До глубины 200–250 м бурение происходит на депрессии из-за альтитуды скважин ниже уровня моря;
  • на м. Молдабек в ГС происходит размыв слабосцементированных песчаников и образование каверн. Напротив, в интервалах более плотных глинистых отложений или глинистых песчаниках происходит образование корки или выпучивание глин;
  • наибольший риск обвалообразования ГС 718 м. С. Нуржанов в горизонтальной секции соответствует глинистому интервалу 1877–1884 м. Расчетный градиент обвалообразования при нулевом обрушении составил 1,35 г/см³, при обрушении 5% – 1,18 г/см³, при обрушении 10% – 1,12 г/см³;
  • в соответствии с выявленными причинами выданы следующие рекомендации по предотвращению осложнений:

1. 1. 1. Рекомендуется проведение керновых исследований для верхнеюрских пластов м. Вост. Молдабек, Жанаталап, С. Балгимбаев, Досмухамбетовское.
      2. Рекомендуется провести микроимиджеровые исследования или запись многорычажного каверномера на месторождениях с планируемыми операциями МГРП для качественного заложения азимута скважин и увеличения объёма стимуляции коллектора.
      3. Для более точных расчетов прочностных характеристик пород рекомендуется на месторождениях провести исследования АКШ как минимум в одной опорной вертикальной скважине.
      4. Рекомендуется проведение тестов ELOT. Это необходимо для калибровки градиента начала поглощения, определения градиента ГРП, расчета глубины спуска колонн.
      5. Рекомендуется производить предбуровое 1D геомеханическое моделирование ГС для оценки рисков нестабильности ствола скважины и снижения аварийности на м. С. Нуржанов, Актобе, Досмухамбетовское и их аналогах, а также скважин с ГРП.

Авторлар туралы

Павел Сергеевич Романенко

«ҚМГ Инжиниринг» ЖШС

Email: p.romanenko@niikmg.kz

ведущий инженер службы онлайн бурения

Қазақстан, Нұр-Сұлтан қ-сы

Илья Георгиевич Алехин

DeGolyer and MacNaughton Corp

Email: I.Alekhin@niikmg.kz

эксперт геомеханик

Қазақстан, Нұр-Сұлтан қ-сы

Канат Берикханович Ашимов

«ҚМГ Инжиниринг» ЖШС

Email: K.Ashimov@niikmg.kz

эксперт службы онлайн бурения

Қазақстан, Нұр-Сұлтан қ-сы

Сарсенбай Абилович Губашев

«ҚМГ Инжиниринг» ЖШС Атыраулық филиалы

Хат алмасуға жауапты Автор.
Email: Gubashev.S@kmge.kz

директор департамента проектирования бурения

Қазақстан, Атырау қ-сы

Әдебиет тізімі

Қосымша файлдар