RU2730576C1 - Способ прогнозирования гидроразрыва пласта, способ гидроразрыва пласта, способы прогнозирования рисков гидроразрыва пласта - Google Patents
Способ прогнозирования гидроразрыва пласта, способ гидроразрыва пласта, способы прогнозирования рисков гидроразрыва пласта Download PDFInfo
- Publication number
- RU2730576C1 RU2730576C1 RU2019134378A RU2019134378A RU2730576C1 RU 2730576 C1 RU2730576 C1 RU 2730576C1 RU 2019134378 A RU2019134378 A RU 2019134378A RU 2019134378 A RU2019134378 A RU 2019134378A RU 2730576 C1 RU2730576 C1 RU 2730576C1
- Authority
- RU
- Russia
- Prior art keywords
- fracture
- hydraulic
- hydraulic fracturing
- fracturing
- proppant
- Prior art date
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 49
- 239000000463 material Substances 0.000 claims abstract description 76
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 230000021615 conjugation Effects 0.000 claims abstract description 5
- 239000012530 fluid Substances 0.000 claims description 100
- 238000002347 injection Methods 0.000 claims description 36
- 239000007924 injection Substances 0.000 claims description 36
- 238000013479 data entry Methods 0.000 claims description 12
- 239000000835 fiber Substances 0.000 claims description 12
- 239000000654 additive Substances 0.000 claims description 11
- 229920000642 polymer Polymers 0.000 claims description 7
- -1 proppant Substances 0.000 claims description 3
- 239000007788 liquid Substances 0.000 abstract description 10
- 230000000694 effects Effects 0.000 abstract description 9
- 238000013459 approach Methods 0.000 abstract description 7
- 230000000638 stimulation Effects 0.000 abstract description 7
- 239000000126 substance Substances 0.000 abstract description 4
- 238000005755 formation reaction Methods 0.000 description 50
- 238000009826 distribution Methods 0.000 description 45
- 238000012546 transfer Methods 0.000 description 26
- 238000004364 calculation method Methods 0.000 description 19
- 239000000243 solution Substances 0.000 description 19
- 239000000725 suspension Substances 0.000 description 13
- 238000011109 contamination Methods 0.000 description 11
- 238000011010 flushing procedure Methods 0.000 description 9
- 239000011435 rock Substances 0.000 description 9
- 238000005086 pumping Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000000518 rheometry Methods 0.000 description 6
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 description 5
- 238000013461 design Methods 0.000 description 5
- 239000003349 gelling agent Substances 0.000 description 5
- 239000002245 particle Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 230000000977 initiatory effect Effects 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 230000001902 propagating effect Effects 0.000 description 3
- 238000012502 risk assessment Methods 0.000 description 3
- 230000004936 stimulating effect Effects 0.000 description 3
- 229920003169 water-soluble polymer Polymers 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 239000003638 chemical reducing agent Substances 0.000 description 2
- 230000001268 conjugating effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000012065 filter cake Substances 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000008719 thickening Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 239000004971 Cross linker Substances 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- HJUFTIJOISQSKQ-UHFFFAOYSA-N fenoxycarb Chemical compound C1=CC(OCCNC(=O)OCC)=CC=C1OC1=CC=CC=C1 HJUFTIJOISQSKQ-UHFFFAOYSA-N 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000009828 non-uniform distribution Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 239000004094 surface-active agent Substances 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
- E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F9/00—Arrangements for program control, e.g. control units
- G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
- G06F9/44—Arrangements for executing specific programs
- G06F9/455—Emulation; Interpretation; Software simulation, e.g. virtualisation or emulation of application or operating system execution engines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Environmental & Geological Engineering (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Operations Research (AREA)
Abstract
Изобретение относится к нефтегазовой промышленности и может найти применение при стимулировании подземного пласта с помощью операции гидравлического разрыва (ГРП) пласта, в частности, при использовании методов математического моделирования, которые позволяют делать прогноз геометрии трещины ГРП и размещения в ней жидкостей, расклинивающих агентов (проппанта), волокон и других материалов. Настоящее решение является способом прогнозирования ГРП, который расширяет возможности существующего псевдотрехмерного (P3D) подхода. Способ гидроразрыва пласта содержит ввод данных о расписании гидроразрыва, ввод данных для подмодели распространения трещины гидроразрыва в пласте. Ввод данных для подмодели транспорта материалов гидроразрыва в трещине. Моделирование гидроразрыва пласта с помощью подмодели распространения трещины гидроразрыва и подмодели транспорта материалов гидроразрыва. При этом сопряжение подмоделей проводят по гидродинамическому сопротивлению для ячеек моделирования. Вывод данных моделирования на момент окончания расписания гидроразрыва и/или закрытия трещины. 12 з.п. ф-лы, 11 ил., 17 табл.
Description
ОБЛАСТЬ ТЕХНИКИ
Решение относится к нефтегазовой промышленности и может найти применение при стимулировании подземного пласта с помощью операции гидравлического разрыва (ГРП) пласта, в частности, при использовании методов математического моделирования, которые позволяют делать прогноз геометрии трещины ГРП и размещения в ней жидкостей, расклинивающих агентов (проппанта), волокон и других материалов.
ПРЕДШЕСТВУЮЩИЙ УРОВЕНЬ ТЕХНИКИ
Гидроразрыв пласта (ГРП) является признанным способом повышения продуктивности скважины. Постоянное развитие новых технологий гидроразрыва пласта требует создания более точных инструментов прогнозирования ГРП. Например, в технологии HiWAY™ (торговая марка компании Шлюмберже) проппант размещают в трещине неоднородно, такое размещение в виде островков проппанта достигается путем закачки в скважину пульсов жидкости ГРП с проппантом («грязный пульс») и жидкости без проппанта («чистый пульс»). Прохождение последовательных пульсов жидкости ГПП через кластеры перфорации создает в трещине ГРП островки с высокой концентрацией проппанта. Такие островки проппанта выполняют функцию зон расклинивания, а узкие зоны без проппанта образуют сеть каналов в трещине ГРП, что повышает проводимость всей трещины ГРП и повышает эффект стимулирования нефтеносного пласта. При этом расчет (моделирование) расположения расклиненных и нерасклиненных зон в трещине (неравномерное размещение проппанта) позволяет прогнозировать последующий дебит нефтегазоносного пласта после ГРП-стимуляции. Настоящее решение является способом прогнозирования ГРП, который расширяет возможности существующего псевдотрехмерного (P3D) подхода для технологии HiWAY и других современных методов гидроразрыва, обеспечивающих неравномерное распределение проводимости внутри трещины ГРП.
Ниже приведен обзор существующих методов прогнозирования гидроразрыва пласта. Более подробный обзор также доступен в статье J. Adachi et al. Computer Simulation of Hydraulic Fractures, Int. J. of Rock Mechanics and Mining Sciences, 44 (2007) 739-757.
На сегодня в нефтегазовой индустрии при проведении операции ГРП проводят моделирование этого сложного и динамического процесса. Если давление в скважине выше давления гидроразрыва пласта, происходит распространение трещины с одновременной фильтрацией части жидкости в пористый пласт. Хотя операции стимулирования с помощью ГРП широко распространены, из-за осложняющих факторов сохраняется риск проведения таких операций ГРП, которые не достигают поставленных целей по увеличению притока нефтегазовой жидкости через трещину ГРП в ствол скважины. Метод проб и ошибок для идентификации рисков является слишком затратным при проведении операций ГРП.
Во время операций ГРП (закачивание жидкости в ствол скважины при давлении выше давления разрыва пласта) в объеме трещины ГРП имеет место массоперенос материала ГРП (проппанта), а также перенос тепла (обычно жидкость ГРП холоднее температуры пласта). Динамическая модель для процесса гидроразрыва требует совмещения моделей, которые описывают эти взаимосвязанные процессы в трещине ГРП: распространение трещины в пласте (геомеханическая модель) и перенос материалов ГРП (гидродинамическая модель переноса в пространстве трещины жидкой среды, содержащей несущую (вязкую) жидкость, частицы проппанта, волокна и различные добавки). Несущая жидкость должна иметь достаточно высокую вязкость, чтобы обеспечить эффективный перенос проппанта (материал ГРП - или «расклинивающий агент») вглубь трещины. Обычно несущая жидкость - это неньютоновская жидкость в виде водорастворимого полимера (загущающий полимер) или ионно-сшитый водорастворимый полимер (сшитый гель). Реология суспензии проппанта в геле (неньютоновская жидкость) описывается уравнениями степенного закона вязкости (power law rheology). Для двух подмоделей (геомеханической для трещины и гидродинамической для переноса материалов ГРП по объему трещины) необходимо провести сопряжение по одному из физических параметров, чтобы решение по двум подмоделям было физически согласованным.
Далее приведены публикации, описывающие построение моделей ГРП различного уровня сложности.
Модели PKN [Т.K. Perkins, L.R. Kern. Widths of Hydraulic Fractures. Journal of Petroleum Technology, 13 (1961), 9, SPE-89-A] и KGD [A Rapid Method of Predicting Width and Extent of Hydraulically Induced Fractures, 21 (1969), 12, SPE-2458-РА].
Распространение трещины. Высота предполагается постоянной, деформации горной породы происходят в одной плоскости трещины ГРП.
Транспорт (перенос) материалов ГРП. Модель переноса материала внутри трещины не рассчитывается.
Радиальная модель.
Распространение трещины. Предполагается осесимметричная деформация горной породы относительно оси скважины (механически изотропная порода).
Транспорт материалов ГРП. Модель переноса материла внутри трещины не рассчитывается (равномерная жидкость ГРП).
Псевдотрехмерные модели (например, Pseudo3D из FracCADE™, торговая марка компании Шлюмберже).
Распространение трещины. Плоская вертикальная трещина переменной высоты, представленная в виде ячеек, подобных модели PKN. Широко используется в нефтегазовой индустрии. Применима, когда полудлина трещины намного больше ее высоты.
Транспорт материалов ГРП. Рассчитывается по одномерной модели с тем же разрешением, что и у подмодели распространения трещины.
Planar 3D (например, US 6876959, GOHFER™, торговая марка компаний Stim-Lab и Marathon Oil Company, TerraFrac™, торговая марка компании TerraTek, StimPlan™, торговая марка компании NSI).
Распространение трещины. Плоская вертикальная трещина переменной высоты с любым соотношением длины и высоты. Обладает улучшенной точностью по сравнению с псевдотрехмерными моделями в случае, когда пласт-коллектор состоит из нескольких слоев с существенно различающимися свойствами: минимальным горизонтальным сжимающим напряжением, модулем Юнга и другими.
Транспорт материалов ГРП. Рассчитывается по двумерной подмодели, имеющей то же разрешение, что и подмодель распространения трещины (одна мелкая сетка для двух подмоделей).
Модель Planar 3D требует большего времени расчета, чем псевдотрехмерные модели и ограничена для полевого применения (не подходит для быстрой оценки рисков при проведении ГРП).
Трехмерные модели (например, SPE-151585-PA).
Распространение трещины. Криволинейная траектория трещины с произвольными направлениями распространения, не лежащими в одной плоскости.
Транспорт материалов ГРП. Модель переноса материла внутри трещины не рассчитывается.
Модели ГРП в пластах с естественной трещиноватостью (например, модель UFM™, торговая марка компании Шлюмберже, патент US 8412500, международная заявка WO 2008093264).
Распространение трещины. Основано на псевдотрехмерном подходе и подмодели взаимодействия с естественными трещинами. Созданная трещина ГРП встречается с существующей естественной трещиной и дальнейший рост сети трещин описывается специально построенной геомеханической моделью (геомеханика для пересекающихся трещин).
Транспорт материалов ГРП. Рассчитывается только по одномерной модели.
Существующие модели ГРП не позволяют прогнозировать реальные двумерные распределения материалов ГРП (проппанта, волокон, концентрации геля). Настоящее решение направлено на преодоление этих недостатков, а также направлено на прогноз рисков, сопровождающих операции ГРП и на выбор такого расписания ГРП и скважинных операций, которые снижают риски неэффективного ГРП.
КРАТКОЕ ОПИСАНИЕ СУЩНОСТИ
Настоящее решение направлено на создание способа прогнозирования гидроразрыва пласта с учетом процессов массопереноса и теплопереноса. Решение содержит следующие операции:
1. Псевдотрехмерная подмодель распространения трещины (известная из текущего уровня техники). Для расчета в подмодели распространения трещины вводят данные о механических свойствах пласта и данные о реологических параметрах жидкости ГРП (суспензии), которые считаются постоянными в каждой ячейке используемой сетки. Сетка подмодели распространения трещины более крупная, чем для подмодели транспорта материалов ГРП. Ячейку подмодели распространения трещины получают объединением более мелких ячеек для подмодели транспорта материалов ГРП.
2. Двумерная подмодель транспорта материалов ГРП, утечек в пласт и теплообмена с породой, имеющая более высокое разрешение по пространству и времени по сравнению с моделью распространения трещины из п. 1.
3. Метод сопряжения подмоделей из пунктов 1 и 2. В нем обеспечивается равенство гидродинамических сопротивлений для ячейки модели распространения трещины для случаев: (а) постоянных реологических параметров, используемых в п. 1 и (б) распределений реологических параметров, рассчитанных в подмодели транспорта материалов ГРП в п. 2.
Настоящее решение расширяет возможности псевдотрехмерного подхода по прогнозированию ГРП. Во-первых, способ прогнозирования позволяет рассчитать распределение материалов гидроразрыва пласта в трещине с более высоким разрешением, чем в существующих моделях. При этом время расчета остается сопоставимым и достаточным для проведения прогноза ГРП на месте. Во-вторых, алгоритм сопряжения из п. 3 дает возможность учесть влияние распределений материалов ГРП и температуры, рассчитанных с высоким разрешением, на распространение трещины.
Также предложен способ гидроразрыва пласта.
Также предложен способ прогнозирования продуктивности скважины.
Настоящее решение также направлено на прогнозирование рисков при проведении ГРП: явление перепродавки жидкости ГРП, загрязнения трещины водорастворимым гелирующим полимером и частичное закупоривание трещины проппантом (бриджинг). Моделирование с помощью сопряженных подмоделей позволяет осуществить выбор скважинных операций, снижающих риски неэффективного ГРП.
Способ также позволяет прогнозировать параметры ГРП для канального ГРП (неоднородное размещение проппанта в трещине ГРП с образованием каналов).
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Фиг. 1. Изображение трещины ГРП с геометрией, рассчитанной по псевдотрехмерной подмодели распространения трещины (показано одно полукрыло трещины ГРП, моделирование второго крыла идентично).
Фиг. 2. Схема сопряжения подмодели распространения трещины и подмодели переноса материалов ГРП.
Фиг. 3. Результаты расчета закачивания жидкости промывки в трещину по двумерной подмодели переноса материалов ГРП: концентрация проппанта (а) и проводимость трещины (б). Обозначено явление перепродавки проппанта вблизи отверстий перфораций.
Фиг. 4. Результаты расчета закачивания жидкости промывки в трещину по одномерной подмодели переноса материалов ГРП: концентрация проппанта (а) и проводимость трещины (б).
Фиг. 5. Результаты расчета пульсирующего закачивания проппанта в трещину по двумерной подмодели переноса материалов ГРП: концентрация проппанта (а) и проводимость трещины (б).
Фиг. 6. Прогноз кумулятивной добычи для скважины после ГРП, выполненного по Фиг. 5.
Фиг. 7. Распределение объемной концентрации геля (а) и проводимости (б) после закрытия трещины.
Фиг. 8. Распределение температуры несущей жидкости в трещине на момент остановки закачивания (а) и распределение проводимости (б). Расчет сделан по двумерной модели переноса материалов ГРП.
Фиг. 9. Распределение температуры несущей жидкости в трещине на момент остановки закачивания (а) и распределение проводимости (б). Расчет сделан по одномерной модели переноса материалов ГРП.
Фиг. 10. Результаты расчета бриджинга проппанта в узкой трещине по двумерной подмодели роста трещины: концентрация проппанта (а) и проводимость трещины (б).
Фиг. 11. Результаты расчета бриджинга проппанта в узкой трещине по одномерной подмодели роста трещины: концентрация проппанта (а) и проводимость трещины (б).
Используемые термины:
Гидроразрыв пласта (ГРП) - метод интенсификации работы нефтяных или газовых скважин, заключающийся в создании трещины в пласте для обеспечения притока добываемого флюида (газ, вода, конденсат, нефть, либо их смесь). Трещина создается путем закачки материалов ГРП в пласт.
Материалы ГРП - несущие жидкости, расклинивающие агенты, волокна, а также добавки (например, понизители трения, разрушители, сшиватели, стабилизаторы, поверхностно-активные вещества).
Проппант - расклинивающий агент (прочный калиброванный песок или его аналоги). Размеры проппанта позволяют закачивать суспензию проппанта в загущенной жидкости (гелированной жидкости) в ствол скважины, затем через перфорационные отверстия в образованную трещину. Размер проппанта меньше ширины трещины ГРП, что обеспечивает перенос проппанта вглубь трещины.
Свойства проппанта - размер проппанта (сито США), средний диаметр и удельная плотность проппанта позволяют смоделировать перенос частиц проппанта вглубь трещины и оседание проппанта (стоксовское оседание крупных частиц в жидкости).
Проппированная область трещины - область, расклиненная проппантом. Проппированная область имеет проницаемость выше проницаемости окружающей породы.
Распределение проводимости трещины ГРП - функция, принимающая значения локальной проводимости трещины в каждой точке (x, y).
Локальная проводимость трещины - это произведение толщины трещины после ее закрытия и проницаемости упаковки проппанта и каналов, образовавшихся в результате ГРП. Параметры средней проницаемости проппантной упаковки [м2] измеряют в лаборатории в соответствии со стандартом ISO 13503-5.
Риски проведения ГРП - факторы, действия которых снижают общую гидравлическую проводимость трещины ГРП и снижают эффективность стимулирования пласта ниже запланированного уровня.
Частичное закупоривание или бриджинг (риск закупоривания) - остановка движения проппанта (или проппанто-волоконной смеси) в узкой трещине из-за взаимодействия между собой, с жидкостью и со стенками трещины. При проведении ГРП возникновение бриджинга проявляется в виде резкого роста давления в забое скважины и сигнализирует о риске выполнения ГРП не по графику. При моделировании потока жидкости ГРП (суспензии проппанта в несущей жидкости) бриджинг проявляется в виде увеличения гидродинамического сопротивления потоку и остановки транспорта проппанта в направлении потока (нулевая скорость проппанта в моделируемом потоке частиц).
Разрешение (модели, расчета) - количество ячеек сетки на единицу длины (пространственное разрешение), количество шагов по времени на единицу времени (временное разрешение).
Свойства жидкости:
Реологические параметры - параметры K и n степенного закона вязкости жидкости, согласно которому
τ=K(∂u/∂ξ)n, μ=K(∂u/∂ξ)n-1,
где τ - это напряжение сдвига, K - коэффициент густоты потока (consistency index), n - показатель поведения жидкости (behavior index), μ - кажущаяся вязкость жидкости (apparent viscosity), ∂u/∂ξ - скорость сдвига (градиент скорости вдоль оси ξ, перпендикулярной к скорости сдвига слоев жидкости). Параметры K и n несущих жидкостей измеряются в лаборатории в соответствии со стандартами ISO 13503-1 и API 13М.
Параметры скорости фильтрации жидкости:
Коэффициент утечек, м/с0.5 и коэффициент мгновенных потерь, м3/м2 - это параметры модели Картера утечки жидкости в пласт, согласно которой
где VL - это объем утечки жидкости в пласт за время t (время, прошедшее с начала процесса утечек в пласт), AL - площадь поверхности, через которую происходит утечка. Подробно модель утечки описана в книге M.J. Economides, K.G. Nolte, Reservoir Stimulation, 3rd Edition, p.2-17.
Тип реологии несущей жидкости:
Ньютоновская жидкость - жидкость со степенным законом вязкости, в котором n=1. Для такой жидкости реологические параметры включают в себя только вязкость μ=K. Напряжение сдвига определяется формулой τ=K ∂u/∂ξ.
Неньютоновская жидкость - жидкость, вязкость которой зависит от скорости сдвига. В случае степенного закона вязкости это жидкость, для которой n≠1.
Механические свойства пласта - набор параметров, описывающих породу в зоне ГРП и позволяющих оценить геометрию трещины ГРП. Кроме обычных прочностных характеристик породы (модуля Юнга [ГПа], коэффициента Пуассона) необходимо при моделировании трещины учитывать следующие параметры:
Вязкость разрушения KIC (или критический коэффициент интенсивности напряжения) - механический параметр имеет размерность [кПа⋅м1/2] и задает порог механического напряжения, при котором происходит дальнейшее распространение трещины в твердой породе.
Напряжение закрытия (МПа) - если давление внутри трещины ниже этого порога, то происходит закрытие (смыкание) стенок трещины. При этом процесс ГРП прекращается.
ПОДРОБНОЕ ОПИСАНИЕ
В одном или нескольких вариантах осуществления в соответствии с настоящим раскрытием способы могут быть направлены на прогнозирование гидроразрыва пласта и содержать выполнение следующих операций:
Проектирование гидроразрыва пласта
Проектирование ГРП содержит этапы, перечисленные ниже.
1. Получение механических и проводящих свойств пласта из данных каротажа и других средств измерений.
2. Сбор информации о скважине: траектория, заканчивание (место перфораций обсадной трубы) и другие параметры.
3. Получение данных лабораторных измерений свойств материалов гидроразрыва пласта: жидкостей (реологические параметры, плотность), свойств проппанта (плотность и размеры), волокон (например, влияние на скорость осаждения проппанта).
4. Выбор расписания закачки ГРП: выбор материала закачки (вязкие жидкости, проппант и добавки), расход и объем закачки материалов ГРП, концентрации проппанта, волокон и добавок, максимальное давление во время ГРП (превышает давление гидроразрыва).
5. Прогнозирование операции гидроразрыва пласта: моделирование распространения трещины и транспорта материалов ГРП в ней, расчет распределения проводимости трещины.
6. Вывод данных моделирования на момент окончания расписания ГРП (локальная концентрация проппанта, концентрация остаточного геля, распределение проводимости для расклиненной трещины).
7. Дополнительно рассчитывают прогноз продуктивности скважины для полученного в пп. 5-6 распределения проводимости трещины ГРП. Если требуемые параметры скважины достигнуты, то переход к п. 8, или переход к п. 4 и изменение расписания закачки ГРП (изменение объемов материалов ГРП или их пропорций).
8. Проведение гидроразрыва на основе выбранного расписания ГРП (по пункту 4).
Для прогноза продуктивности скважины (этап 7) для заданного распределения проводимости трещины (пласт простимулирован с помощью ГРП) могут быть использованы, например, пакеты программ ECLIPSE™ и INTERSECT™, торговые марки компании Шлюмберже. Эти пакеты программ позволяют работать со сложными моделями трещины ГРП и давать прогноз дебита нефтегазовой жидкости по времени для производительного пласта.
Прогнозирование гидроразрыва пласта
Данный раздел содержит описание предлагаемого способа прогнозирования ГРП из п. 5 процесса проектирования гидроразрыва пласта.
В модели ГРП время закачки разбивается на последовательные шаги t1, t2, …, tN, где t1 - это начало закачки материалов гидроразрыва пласта в скважину, tN - конечный рассматриваемый момент времени. Расчет производится последовательно для шагов от tk до tk+1, k=1, …, N-1. На каждом из шагов по времени необходимо выполнить приведенные ниже этапы.
1. Подмодель распространения трещины
На данном этапе производится расчет распространения трещины по псевдотрехмерной модели. В ней область трещины покрывается одномерной расчетной сеткой, показанной на Фиг. 1. В каждой ячейке этой сетки задаются постоянные реологические параметры суспензии, полученные для предыдущего временного шага t=tk. На первом временном шаге берутся свойства жидкостей, закачиваемых в скважину первыми.
2. Подмодель транспорта материалов ГРП
В данной подмодели область трещины покрывается двумерной расчетной сеткой, которая имеет более высокое разрешение, чем сетка из п. 1. В ячейках этой двумерной сетки для всех материалов ГРП (проппанта, волокон, несущих жидкостей, а также добавок) рассчитываются компоненты скорости движения в трещине. Затем рассчитывается перенос материалов и баланс тепловой энергии на шаге от tk до tk+1.
3. Сопряжение подмоделей распространения трещины и транспорта материалов ГРП.
При сопряжении в каждой ячейке сетки из п. 1 проводится расчет постоянных эффективных реологических параметров K* и n* суспензии, обеспечивающих такое же гидродинамическое сопротивление ячейки, что и распределения K(х,у) и n(х,у), полученные в подмодели транспорта в п. 2 при t=tk+1.
Ниже приведены разделы с детальным описанием каждого этапа по пп. 1-3.
Подмодель распространения трещины
Для описания распространения трещины используется псевдотрехмерная (P3D) модель, в которой вводится одномерная сетка из ячеек, показанная на Фиг. 1. Каждая ячейка в любой момент времени является прямоугольником в плоскости трещины со сторонами параллельными вертикальной и горизонтальным осям ориентации трещины. Положения ячеек изменяются на каждом временном шаге. Поэтому такая сетка называется подвижной.
Для расчета на каждом шаге по времени необходимо задать значения реологических параметров суспензии во всех ячейках сетки (параметры K и n). При этом параметры внутри каждой ячейки считаются постоянными.
Полное описание модели распространения трещины приведено в главе 6 из книги M.J. Economides, K.G. Nolte, Reservoir Stimulation, pp. 6-16 - 6-23, 3d edition, John Willey, 2000.
Подмодель транспорта материалов ГРП
Для расчета распределения давления в трещине и скорости перемещения суспензии при ГРП, используются уравнения теории смазки и переноса. Например, в случае ньютоновской жидкости и проппанта они имеют вид:
где
w - это толщина трещины;
uƒ и up - усредненные по массе и поперек трещины поля скорости жидкости (fluid) и проппанта соответственно;
μs(c) - кажущаяся вязкость суспензии, зависящая от концентрации проппанта с;
g - ускорение свободного падения.
Подробно эта модель описана, например, в публикации [S. Boronin, A. Osiptsov, Two-Continua Model of Suspension Flow in a Hydraulic Fracture, Doklady Physics, 2010, Vol.55, No. 4]. В данном решении предлагается использовать этот подход для случая произвольного количества материалов ГРП (проппант, волокна, добавки).
В конце каждого временного шага подмодель транспорта позволяет рассчитать концентрации всех материалов ГРП в узлах сетки (Фиг. 2, схема слева), имеющей более высокое разрешение по сравнению с грубой сеткой (для расчета распространения трещины, Фиг. 2, схема справа). Затем по концентрациям материалов рассчитывают распределение реологических параметров суспензии (степенной закон для реологии). В каждой ячейке сетки подмодели распространения трещины в результате, как правило, получается неравномерное распределение этих реологических свойств.
Сопряжение подмоделей распространения трещины и транспорта материалов ГРП
Задачей сопряжения подмоделей распространения трещины и транспорта материалов ГРП является обеспечение равенства гидродинамических сопротивлений ячейки модели распространения трещины для случаев: (а) постоянных эффективных параметров K* и n* и (б) распределений K(х,у) и n(х,у), рассчитанных в подмодели транспорта материалов ГРП. Эквивалентный набор ячеек моделирования для двух указанных подмоделей изображен на Фиг. 2 (слева - для подмодели транспорта материалов ГРП, справа - для подмодели распространения трещины).
Согласно варианту раскрытия, крупную ячейку для подмодели распространения трещины разбивают на мелкие ячейки сетки, где рассчитывают подмодель транспорта материалов ГРП (Фиг. 2). Это позволяет упростить и ускорить моделирование геометрии трещины ГРП без снижения точности расчета.
Решение задачи сопряжения производится с помощью метода расчета обобщенных параметров K* и n*, описанного ниже. Рассмотрим его на примере одной прямоугольной ячейки где координаты задают положения ее левой, правой, верхней и нижней границ, соответственно.
1. Рассчитать течения через ячейку С, возникающие вследствие перепадов давления Δр1=0.9 Δр и Δр2=1.1 Δр между левой и правой границами, соответственно. При этом использовать распределения реологических параметров K(х,у) и n(x,у), рассчитанные в подмодели транспорта, и известное распределение толщины трещины w(x,y). Здесь Δр - это перепад давления между левой и правой границами ячейки С, полученный в подмодели транспорта на последнем временном шаге.
2. Вычислить суммарные потоки q1 и q2 в направлении оси х через ячейку С, соответствующие перепадам давлений Δр1 и Δр2 из п. 1.
3. Вычислить такие значения K* и n*, чтобы для течения через ячейку C с этими постоянными реологическими параметрами и толщиной w(x,y) из п. 1 суммарные потоки совпадают с q1 и q2 при перепадах давления Δр1 и Δр2, соответственно. Такое условие дает два уравнения для определения искомых постоянных K* и n* (для крупной ячейки С).
В одном или нескольких вариантах осуществления в соответствии с настоящим раскрытием способы могут быть направлены на прогнозирование продуктивности скважины после гидроразрыва пласта и содержать выполнения: ввода данных о расписании ГРП; ввода данных для подмодели распространения трещины ГРП в пласте; ввода данных для подмодели транспорта материалов ГРП в трещине ГРП; моделирования ГРП с помощью подмодели распространения трещины ГРП и подмодели транспорта материалов ГРП, при этом сопряжение подмоделей проводят по гидродинамическому сопротивлению для ячеек моделирования; вывода данных моделирования на момент окончании расписания ГРП и/или закрытия трещины. После получения прогноза о распределении проводимости трещины ГРП, используя программы для моделирования притока из производительного коллектора, получают прогноз продуктивности скважины ГРП по времени (в виде накопленной продуктивности или текущей продуктивности скважины, пересекающей пласт).
В одном или нескольких вариантах осуществления в соответствии с настоящим раскрытием способы могут быть направлены на прогнозирование различных ситуаций, связанных с риском проведения гидроразрыва пласта.
Так в некоторых вариантах осуществления способы могут выполняться для прогнозирования перепродавки в трещине ГРП и содержать выполнения: ввода данных о расписании ГРП; ввода данных для подмодели распространения трещины ГРП в пласте; ввода данных для подмодели транспорта материалов ГРП в трещине ГРП; моделирования ГРП с помощью подмодели распространения трещины ГРП и подмодели транспорта материалов ГРП, при этом сопряжение подмоделей проводят по гидродинамическому сопротивлению для ячеек моделирования; моделирования ГРП для этапа закачивания порции промывочной жидкости в область трещины ГРП; вывода данных моделирования на момент окончании расписания ГРП и/или закрытия трещины. Наконец, получения данных о проводимости трещины ГРП вблизи перфорационных отверстий. При необходимости, изменяют расписание ГРП для снижения эффекта перепродавки проппанта при закачивании порции для очистки ствола скважины от проппанта (продавочная жидкость).
В некоторых вариантах осуществления способы могут выполняться для прогнозирования пульсирующей закачки проппанта и содержать выполнения: ввода данных о расписании ГРП; ввода данных для подмодели распространения трещины ГРП в пласте; ввода данных для подмодели транспорта материалов ГРП в трещине ГРП; моделирования ГРП с помощью подмодели распространения трещины ГРП и подмодели транспорта материалов ГРП, при этом сопряжение подмоделей проводят по гидродинамическому сопротивлению для ячеек моделирования; моделирования ГРП для последовательности закачки импульсов жидкости без проппанта и жидкости с проппантом; вывода данных моделирования на момент окончании расписания ГРП и/или закрытия трещины. По выведенным данным получают и оценивают данные о проводимости трещины ГРП.
В некоторых вариантах осуществления способы могут выполняться для прогнозирования загрязнения трещины ГРП гелирующим полимером и содержать выполнения: ввода данных о расписании ГРП; ввода данных для подмодели распространения трещины ГРП в пласте; ввода данных для подмодели транспорта материалов ГРП в трещине ГРП; моделирования ГРП с помощью подмодели распространения трещины ГРП и подмодели транспорта материалов ГРП, при этом сопряжение подмоделей проводят по гидродинамическому сопротивлению для ячеек моделирования; моделирования ГРП для последовательности закачки импульсов жидкости без проппанта и жидкости с проппантом; вывода данных моделирования на момент окончании расписания ГРП и/или закрытия трещины. Далее определяют зоны с высокой концентрацией гелирующего полимера (то есть зоны загрязнения трещины). При высоких рисках загрязнения пересматривают расписание ГРП (концентрацию гелирующего агента в жидкости ГРП).
В некоторых вариантах осуществления способы могут выполняться для прогнозирования закупоривания трещины ГРП и содержать выполнения: ввода данных о расписании ГРП; ввода данных для подмодели распространения трещины ГРП в пласте; ввода данных для подмодели транспорта материалов ГРП в трещине ГРП; моделирования ГРП с помощью подмодели распространения трещины ГРП и подмодели транспорта материалов ГРП, при этом сопряжение подмоделей проводят по гидродинамическому сопротивлению для ячеек моделирования; вывода данных моделирования на момент окончании расписания ГРП и/или закрытия трещины. Далее определяют зоны низкой проводимости внутри трещины ГРП или на входе в трещину ГРП. При необходимости пересматривают расписание ГРП.
ПРИМЕРЫ
В данном разделе представлены примеры, демонстрирующие возможности предлагаемого решения (моделирования ГРП) по сравнению с доступным уровнем техники. Например, в сравнении с моделью, объединяющей псевдотрехмерную (P3D) подмодель роста трещины и одномерную (1D) подмодель переноса материалов ГРП. Далее этот сравнительный вариант моделирования ГРП обозначен как «известный уровень техники».
Пример 1. Прогнозирование перепродавки жидкости промывки в трещину ГРП
Ситуация
Для промывки скважины после ГРП в нее закачивается низковязкая жидкость без проппанта. При провелении завершающего этапа ГРП существует риск перепродавки -проникновения в призабойную зону жидкости промывки и вымывание проппанта из этой зоны. Это приводит к выносу проппанта из участков, примыкающих к перфорациям. Избыточная перепродавка (наличие расклинненных зон в трещине ГРП) снижает проводящие свойства трещины, что, в свою очередь, снижает продуктивность скважины. С другой стороны, недостаточная промывка скважины приводит к избыточному количеству проппанта, который остался вне трещины ГРП и требует дополнительных скважинных операций и создает проблемы при размещении скважинного оборудования. Для оценки рисков, связанных с перепродавкой (или «недопродавкой»), проектирование ГРП нуждается в прогнозе реального размещения проппаната после промывки чистой жидкостью под высоким давлением (выше давления гидроразрыва) и в идентификации нерасклиненных зон внутри трещины ГРП.
Решение
Проведены расчеты с помощью предлагаемого решения. Для примера рассматривали трещину, распространяющуюся в пластах, описанных в таблице 1. Параметры закачки и свойства материалов ГРП приведены в таблицах 2, 3, 4. Утечки в пласт нулевые (предположение о низкопроницаемом пласте). Закачка жидкости ГРП (неньютоновская жидкость) производится в вертикальную скважину через перфорационный интервал от 4091 м до 4092 м (истинная глубина). Инициация трещины ГРП происходит в пласте 2. Промывка на финальной стадии ГРП проведена водной жидкостью (ньютоновская жидкость с низкой явзкостью).
На Фиг. 3 приведены результаты расчета для завершающей стадии ГРП -моделирование потока жидкостей на пространственной сетке с числом ячеек 257*257. Результаты моделирования закачивания жидкости промывки 2 (водной жидкости на завершающей стадии ГРП) представлены в виде двух характеристик: концентрация проппанта (Фиг. 3а) и проводимость трещины (Фиг. 3б). Отмечено явление перепродавки проппанта вблизи отверстий перфораций (создает риск локального смыкания стенок трепщине в зонах обедненных проппантом). Отверстия перфорации (а также «кластер перфорации» - совокупность отверстий) изображены в виде небольших треугольников между стволом скважины и объемом трещины ГРП. На Фиг. 3 видно, что возле перфораций при закачке жидкости 2 (жидкость промывки) образуется поток с неусточивостью Саффмана-Тейлора (неравномерность течения двух жидкостей из-за большого контраста вязкости). Эта неравномерность потока влияет на окончательное распределение (вымывание) проппанта и проводимости проппантой упаковки вблизи перфорационного кластера. В одномерной модели транспорта материалов ГРП этот эффект не учитывается (см. Фиг. 4) и риск локального снижения проводимости трещины вблизи перфораций не моделируется и не учитывается.
Вывод
Предложенное решение позволяет точнее прогнозировать распределение проводимости трещины после перепродавки (по сравнению с псевдотрехмерной моделью и одномерной моделью траспорта материалов ГРП). На основе этого прогноза оператор выбирает (напрмер, уменьшает) объем жидкости промывки и повторяет прогноз для нового расписания ГРП, что снижает риск появления нерасклиненных зон в трещине. С другой стороны, маломасшабная перепродавка проппанта вблизи кластера перфораций создает преимущества для дальнейших скважинных операций (ствол скважины свободен от остатков проппанта).
Пример 2. Моделирование пульсирующей закачки проппанта (неравномерное размещение проппанта в трещине ГРП)
Ситуация
Необходимо сделать прогноз распределения материалов ГРП и проводимости трещины при чередующейся закачке проппанта короткими пульсами переменно с пульсами чистой жидкости без проппанта.
Решение
Проведены расчеты с помощью предлагаемого решения. Для примера рассматривали трещину, полученную после закачки расписания, указанного в таблице 5. Остальные параметры жидкости и проппанта совпадают с теми, что указаны выше в примере 1.
На Фиг. 5 показаны результаты расчета, полученные с помощью двумерной модели транспорта материалов ГРП. На графиках видно неравномерное распределение проппанта, полученное с помощью подачи пульсирующей закачки. На Фиг. 5(б) видно, что проводимость трещины ГРП зависит от полученного распределения проппанта и каналов (для каналов характерна низкая концентрация проппанта вокруг проппантных агрегатов, непроппированные зоны). На Фиг. 6 приведен прогноз накопленной продуктивности скважины для полученного распределения проводимости трещины. Пунктирной линией отображается целевой уровень продуктивности для периода 120 дней. Видно, что выбранное расписание закачки позволяет достичь этого уровня, то есть этот вариант ГРП достигает поставленной задачи.
Если уровень продуктивности не достигается, то оператор пересматривает параметры закачки или стратегию перфорации обсадной колонны.
Вывод
Предложенное решение позволяет прогнозировать неравеномерное распределение проводимости трещины и материалов ГРП при пульсирующей закачке проппанта. Полученные данные позволяют прогнозировать продуктивность скважины. На основе этого прогноза оператор принимает решение, например, о выборе скважин, где проведение ГРП с пульсирующей закачкой создает запланированную продуктивность.
Пример 3. Моделирование загрязнения трещины гелирующим агентом
Ситуация
Основная цель ГРП - повышение продуктивности скважины. Эта цель достигается созданием в породе трещины, которая увеличивает поверхность контакта между скважиной и пластом и обеспечивает транспорт углеводородов к перфорациям. При проведении ГРП низкая концентрация гелирующего агента снижает транспортные способности несущей жидкости, а избыточная концентрация гелирующего агента вызывает загрязнение породы и образование фильтрационной корки. Выбор между этими крайностями найден при помощи моделирования распределения концентрации геля на момент окончания ГРП (прогноз зон загрязнения). Технология ГРП подразумевает использование веществ, которые снижают результативность стимуляции, создавая низкопроводящую фильтрационную корку на стенках трещины или закупоривая пространство между частицами проппанта (загрязнение упаковки проппанта). Прогноз загрязнения трещины необходим для выбора концентрации геля, положения перфораций и, возможно, стратегии очистки трещины, которые дали бы желаемый прирост продуктивности скважины.
Основным средством прогнозирования загрязнения трещины является математическое моделирование оседания геля в результате фильтрации жидкости гидроразрыва в пласт. Однако распределение геля, полученное при использовании одномерного (упрощенного) приближения для течения в трещине не может служить основанием для выбора расписания закачки ГРП.
Решение
Проведены расчеты с помощью предлагаемого решения. Для примера рассматривали трещину, распространяющуюся в пластах, описанных в Таблице 6. Свойства жидкости и проппанта приведены в таблицах 7, 8. Таблица 9 содержит параметры закачки. Закачка производится в вертикальную скважину через перфорационный интервал от 2720 м до 2725 м (истинная глубина). Утечки жидкости в пласт описывали моделью Картера, представленной, например, в статье A. Settari, General Model of Fluid Flow (Leakoff) From Fractures Induced in Injection Operations, SPE-18197-MS. В модели для ГРП задавали параметры: коэффициент утечки (leakoff coefficient) и коэффициент мгновенных потерь (spurt loss).
На Фиг. 7 представлено распределение концентрации геля в трещине после ее закрытия. Трещина инициируется в пласте 4 рядом с перфорированным участком скважины и с течением времени достигает пласта 2 через промежуточный пласт 3 с другими свойствами.
Концентрация геля в среднем пласте 3 максимальна. Это значит, что при выводе трещины на добывающий режим, гелевая пробка на этом участке затруднит транспорт углеводородов из верхней части трещины к перфорациям, как это можно видеть из распределения проводимости на Фиг. 7(б). Область высокой проводимости в нижней части трещины рядом с перфорациями, обозначенная прямоугольником на Фиг. 7(а), сообщается с верхней частью трещины через узкий проводящий участок, которого недостаточно для эффективного транспорта добываемой нефтегазовой жидкости. Чтобы оптимизировать дизайн, необходимо пересмотреть перфорационную стратегию, например, дополнительно проперфорировать участок скважины в пласте 2.
В одномерной модели транспорта материалов ГРП указанный эффект не учитывается, так как не рассчитывается распределение концентрации геля по вертикали. В результате нельзя сделать вывод о контакте нефтенасыщенных пластов, пересеченных трещиной, с перфорациями. Следовательно, риски потери притока из этих пластов не могут быть предсказаны.
Вывод
Предложенный способ позволяет моделировать загрязнение трещины, возникающее в процессе ГРП, и принимать информированное решение о возможных способах оптимизации графика закачки ГРП и заканчивания скважины. Например, при высокой концентрации геля (gel damage) оператор принимает решение о проведении дополнительного этапа - очистка трещины (fracture cleanup).
Пример 4. Расчет распределения температуры жидкости ГРП
Ситуация
При проведении гидроразрыва химические добавки к жидкости подбираются с учетом температуры пласта таким образом, чтобы при нагревании избежать деградации загущающего полимера, то есть сохранить вязкость и несущие способности жидкости с момента входа в трещину и до ее закрытия. Для разрушения несущей жидкости (например, сшитого геля) добавляют понизителя вязкости (разрушители геля). Активация разрушителей геля обычно производится при температуре выше критической для данного геля, и деструкция геля нежелательна до окончания ГРП и окончательного размещения проппанта для максимальной проводимости трещины. Различные деградируемые материалы, применяемые для улучшения траспорта проппанта (и для отклонения потока от одной трещины ГРП к другой трещине ГРП) также имеют ограничения по температуре. Исходя из вышесказанного, предсказание времени, за которое закачиваемые порции жидкости прогреваются от поверхностной температуры (в устье скважины) до температуры пласта актуально при составлении расписания закачки. Преждевременый прогрев несущей жидкости ГРП (температурное разрушение вязкости) приводит к риску осаждения проппанта (оседание проппанта на дне трещины) и снижению проводимости трещины (неполное расклинивание трещины).
Решение
Проведено моделирование с помощью предлагаемого решения (две сопряженные подмодели). Для примера рассматривали ситуацию, где таблица 10 представляет параметры пластов. Таблица 11 содержит реологию закачиваемой жидкости. Отметим, что при температуре 38°С вязкость жидкости >100 сР при 170 с-1, а при нагревании до 79°С жидкость имеет вязкость 1 сР, т.е. практически полностью теряет несущие свойства. Коэффициент утечек в пласт Ct=3.9Е-6 м/с0.5. Температура пласта 80°С. Температура жидкости на по верности 27°С. Таблица 12 содержит свойства проппанта. Таблица 13 представляет параметры закачки (расписание ГРП). Закачка стадий производится в вертикальную скважину через перфорационный кластер, созданный в интервале от 2716 м до 2746 м (истинная глубина). Инициация трещины происходит в пласте номер 2 (зона перфорации).
Таблица 11 показывает, что для температурного интервала 38-66°С реология вязкого геля (колонка 4) меняется незначительно даже при экспозиции 1-2 часа (характерное время проведения ГРП). При более высоких температурах (>66°С) происходит разрушение геля (достигается вязкость, равная вязкости воды). Такое изменение вязкости с повышением температуры требует учета температурного фактора при прогнозировании операции ГРП.
В этом примере исследовали влияние двумерной и одномерной моделей переноса тепла в трещине ГРП на распределение проппанта в ней. Распределение температуры влияет на вязкость и скорость оседания проппанта, которые, в свою очередь, влияют на конечные распределения концентарции проппанта и проводимости трещины ГРП.
На Фиг. 8(а) и Фиг. 9(а) представлены распределения температур жидкости в трещине на момент окончания закачки. Видно, что форма охлажденной области для рассмотренных двух вариантов различается. При одномерной модели область трещины состоит из вертикальных полосок с одинаковой температурой, а в двумерной модели область состоит из дуг. Упрощенная (одномерная) модель показывает, что верхний и нижний кончики трещины контактируют с охлажденной жидкостью, в то время как двумерная модель предсказывает нагретую жидкость у всех кончиков трещины. Отличие в площади охлажденной области объясняется различными полями скоростей и объемами утечек в пласт. Область распространения трещины, обозначенная сплошной линией, составила 265 м на Фиг. 8(б) и 240 м на Фиг. 9(б). Высоты трещин в обоих моделях совпадают, т.к. были искуственно ограничены с помощью непроницаемых барьеров.
На Фиг. 8(6) и Фиг. 9(б) представлены рассчитанные распределения проводимости в трещине. На Фиг. 8(б) проппант осел в основном на дне трещины, что привело к пониженной проводимости в центральной части трещины по сравнению с распределением проводимости из Фиг. 9(б). Это вызвано тем, что температура жидкости в этой области на Фиг. 9(а) ниже, чем на Фиг. 8(а) и жидкость при одномерном расчете имеет большую вязкость в соответствии с таблицей 11. Прогноз распределения температуры и проводимости трещины показал наличие риска неэффективного ГРП.
Вывод
Предложенное решение позволяет моделировать влияние двумерного распределения температуры на реологию несущей жидкости, учитывать ее влияние на геометрию трещины, перенос и оседание проппанта. На основе этих прогнозов оператор может принять решение, например, о выборе типа несущей жидкости, соответствующей полученным температурам в трещине. Так, для трещины ГРП с более высокой средней температурой потребутся гелирующий агент, приспособленный к более высокой температуре.
Пример 5. Моделирование бриджинга проппанта и оценка риска проведения ГРП
Ситуация
Технология осуществления ГРП включает в себя закачку в скважину жидкости ГРП (гель на основе водорастворимого полимера, в некоторых случаях водная жидкость с добавками) при давлениях выше давления разрыва нефтеносного пласта. Для поддержания трещины в открытом состоянии используют проппант. В процессе закачки проппанта возможно возникновение частичного, либо полного закупоривания трещины (screen out). Это происходит из-за того, что проппант, переносимый несущей жидкостью, создает пробку в отверстиях перфорации или внутри трещины. Это вызывает значительное сопротивление течению жидкости и регистрируемый резкий рост давления на устье скважины. Рост трещины ГРП в длину прекращается, работа по стимуляции пласта останавливается.
На практике дизайн работы (расписание ГРП) создается таким образом, чтобы избежать неконтролируемого закупоривания трещины. Для оценки рисков, связанных с бриджингом при ГРП необходимо иметь возможность прогнозировать размещение проппанта в трещине с высокой точностью. Двумерное моделирование транспорта материалов ГРП позволяет спрогнозировать размещение проппанта, в то время как одномерный подход может приводить к значительным погрешностям и неправильной оценке рисков.
Решение
Проведены расчеты с помощью предлагаемого решения. Для примера рассматривали трещину, распространяющуюся в пластах, описанных в таблице 14. Инициация происходит в пласте 6. Параметры закачки и свойства жидкостей приведены в таблицах 15, 16, 17. Закачка жидкости ГРП производится в вертикальную скважину через 3 перфорационных интервала (истинная глубина): 1) от 3430 м до 3437 м; 2) от 3439 м до 3446 м; 3) от 3449 м до 3460 м.
На Фиг. 10 и Фиг. 11 представлены распределения концентрации проппанта в трещине на момент окончания закачки для одномерного и двумерного подходов моделирования транспорта материалов ГРП. В одномерной модели бриджинг проппанта произошел на расстоянии 75 м от скважины после закачки 7 м3 жидкости на 4 шаге расписания. Основная причина данного эффекта заключается в том, что в одномерной подмодели транспорта бриджинг в одной ячейке полностью останавливает рост трещины в длину. В представленном примере полученная конечная длина трещины равна 84 м (Фиг. 11).
Сопряжение двумерной модели транспорта с псевдотрехмерной моделью распространения трещины позволяет рассчитать с высоким разрешением размещение проппанта с учетом бриджинга в локальных областях. На Фиг. 10 видна зона бриджинга (обозначена белой пунктирной линией) в интервале по оси х от 30 м до 50 м. При этом из-за гетерогенного распределения проппанта трещина продолжает расти и заполняться суспензией. В представленном примере конечная длина трещины равна 115 м (Фиг. 10(б)).
Вывод
Предложенное решение позволяет реалистично моделировать обтекание суспензией зон, в которых произошел бриджинг проппанта, что оказывает существенное влияние на геометрию трещины и перенос проппанта в целом. При использовании одномерной подмодели траспорта материалов ГРП (предыдущий уровень техники) получается прогноз ГРП, в котором длина трещины меньше, чем при использовании двумерной подмодели. На основе двумерного расчета переноса материалов оператор принимает информированное решение, например, о выборе концентрации проппанта или скорости закачивания суспензии, что направлено на уменьшение бриджинга внутри трещины.
Очевидно, что описанные выше варианты осуществления не должны рассматриваться в качестве ограничения объема патентных притязаний настоящего раскрытия. Для любого специалиста в данной области техники понятно, что есть возможность внести множество изменений в описанные выше способы и, без отхода от принципов раскрытия, заявленные в формуле.
Claims (25)
1. Способ гидроразрыва пласта, содержащий:
а) ввод данных о расписании гидроразрыва;
б) ввод данных для подмодели распространения трещины гидроразрыва в пласте;
в) ввод данных для подмодели транспорта материалов гидроразрыва в трещине;
г) моделирование гидроразрыва пласта с помощью подмодели распространения трещины гидроразрыва и подмодели транспорта материалов гидроразрыва,
при этом сопряжение подмоделей проводят по гидродинамическому сопротивлению для ячеек моделирования;
д) вывод данных моделирования на момент окончания расписания гидроразрыва и/или закрытия трещины.
2. Способ по п. 1, в котором выполняют операцию гидроразрыва на основе расписания гидроразрыва.
3. Способ по п. 1, в котором расписание гидроразрыва содержит данные по концентрации проппанта, данные по концентрации волокон и добавок, данные по свойствам жидкости гидроразрыва.
4. Способ по п. 1, в котором подмодель распространения трещины гидроразрыва в пласте описывают с помощью уравнения баланса для давления и напряжения на стенках трещины.
5. Способ по п. 1, в котором подмодель транспорта материалов гидроразрыва описывают с помощью уравнений баланса массы для жидкости гидроразрыва, проппанта, волокон и добавок, а также уравнений баланса импульса и энергии.
6. Способ по п. 1, в котором вывод данных моделирования содержит вывод данных на момент окончания расписания гидроразрыва или закрытия трещины для распределения концентрации проппанта, проводимости трещины, температуры жидкости гидроразрыва.
7. Способ по п. 1, в котором материалы для гидроразрыва содержат проппант, волокна, добавки.
8. Способ по п. 1, дополнительно содержащий:
получение данных о продуктивности скважины после гидроразрыва пласта на основе данных моделирования.
9. Способ по п. 1, в котором
вводят данные для расписания гидроразрыва, описывающие последовательную закачку пульсов жидкости гидроразрыва с проппантом и без проппанта;
получают прогноз проводимости трещины; и
выполняют гидроразрыв пласта на основе данных о проводимости трещины.
10. Способ по п. 1, содержащий:
(а) вычисление данных о зонах в смоделированной трещине гидроразрыва;
(б) изменение расписания гидроразрыва на основе моделирования трещины на момент окончания расписания гидроразрыва и данных на этапе (а).
11. Способ по п. 10, в котором зоны в смоделированной трещине гидроразрыва имеют низкую гидравлическую проводимость вблизи отверстий перфорации, что приводит к ситуации перепродавки проппанта.
12. Способ по п. 10, в котором зоны в смоделированной трещине гидроразрыва имеют высокую концентрацию гелирующего полимера, что приводит к загрязнению трещины гелирующим полимером.
13. Способ по п. 10, в котором зоны в смоделированной трещине гидроразрыва имеют низкую локальную гидравлическую проводимость, что приводит к закупориванию трещины.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/RU2017/000282 WO2018203765A1 (ru) | 2017-05-02 | 2017-05-02 | Способ прогнозирования рисков гидроразрыва пласта |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2730576C1 true RU2730576C1 (ru) | 2020-08-24 |
Family
ID=64016966
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2019134378A RU2730576C1 (ru) | 2017-05-02 | 2017-05-02 | Способ прогнозирования гидроразрыва пласта, способ гидроразрыва пласта, способы прогнозирования рисков гидроразрыва пласта |
Country Status (4)
Country | Link |
---|---|
US (1) | US11499406B2 (ru) |
CA (1) | CA3062854A1 (ru) |
RU (1) | RU2730576C1 (ru) |
WO (1) | WO2018203765A1 (ru) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2745142C1 (ru) * | 2020-09-17 | 2021-03-22 | Общество с ограниченной ответственностью «Газпромнефть Научно-Технический Центр» | Способ и система моделирования трещин гидроразрыва пласта бесконечно-конечной проводимости и поперечно-продольного расположения относительно горизонтального ствола скважины |
RU2775034C1 (ru) * | 2021-07-28 | 2022-06-27 | Публичное акционерное общество "Газпром нефть" | Способ выбора оптимального дизайна гидроразрыва пласта на основе интеллектуального анализа полевых данных для увеличения добычи углеводородного сырья |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113833448A (zh) * | 2020-06-23 | 2021-12-24 | 中国石油化工股份有限公司 | 一种高通道压裂方法及系统 |
US11933165B2 (en) | 2021-03-15 | 2024-03-19 | Saudi Arabian Oil Company | Hydraulic fracture conductivity modeling |
US11921250B2 (en) * | 2022-03-09 | 2024-03-05 | Saudi Arabian Oil Company | Geo-mechanical based determination of sweet spot intervals for hydraulic fracturing stimulation |
US11867047B2 (en) | 2022-06-08 | 2024-01-09 | Saudi Arabian Oil Company | Workflow to evaluate the time-dependent proppant embedment induced by fracturing fluid penetration |
CN115034636A (zh) * | 2022-06-21 | 2022-09-09 | 中海石油(中国)有限公司 | 基于储层双甜点可压裂性三维展布的压裂风险评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2561114C2 (ru) * | 2010-12-30 | 2015-08-20 | Шлюмбергер Текнолоджи Б.В. | Система и способ для выполнения операций интенсификации добычи в скважине |
EA022370B1 (ru) * | 2010-12-10 | 2015-12-30 | Бейкер Хьюз Инкорпорейтед | Способ усовершенствования модели коллектора и повышения отдачи трещиноватых пластов |
RU2575947C2 (ru) * | 2011-11-04 | 2016-02-27 | Шлюмбергер Текнолоджи Б.В. | Моделирование взаимодействия трещин гидравлического разрыва в системах сложных трещин |
WO2016032489A1 (en) * | 2014-08-28 | 2016-03-03 | Landmark Graphics Corporation | Optimizing multistage hydraulic fracturing design based on three-dimensional (3d) continuum damage mechanics |
WO2016153953A1 (en) * | 2015-03-20 | 2016-09-29 | Schlumberger Technology Corporation | Stacked height growth fracture modeling |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5228510A (en) * | 1992-05-20 | 1993-07-20 | Mobil Oil Corporation | Method for enhancement of sequential hydraulic fracturing using control pulse fracturing |
US6876959B1 (en) | 1999-04-29 | 2005-04-05 | Schlumberger Technology Corporation | Method and apparatus for hydraulic fractioning analysis and design |
US20030205376A1 (en) * | 2002-04-19 | 2003-11-06 | Schlumberger Technology Corporation | Means and Method for Assessing the Geometry of a Subterranean Fracture During or After a Hydraulic Fracturing Treatment |
US8412500B2 (en) | 2007-01-29 | 2013-04-02 | Schlumberger Technology Corporation | Simulations for hydraulic fracturing treatments and methods of fracturing naturally fractured formation |
EA025473B1 (ru) * | 2011-02-23 | 2016-12-30 | Лэндмарк Графикс Корпорейшн | Способ и система для определения подходящих вариантов гидроразрыва формации |
-
2017
- 2017-05-02 WO PCT/RU2017/000282 patent/WO2018203765A1/ru active Application Filing
- 2017-05-02 CA CA3062854A patent/CA3062854A1/en active Pending
- 2017-05-02 RU RU2019134378A patent/RU2730576C1/ru active
- 2017-05-02 US US16/610,231 patent/US11499406B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EA022370B1 (ru) * | 2010-12-10 | 2015-12-30 | Бейкер Хьюз Инкорпорейтед | Способ усовершенствования модели коллектора и повышения отдачи трещиноватых пластов |
RU2561114C2 (ru) * | 2010-12-30 | 2015-08-20 | Шлюмбергер Текнолоджи Б.В. | Система и способ для выполнения операций интенсификации добычи в скважине |
RU2575947C2 (ru) * | 2011-11-04 | 2016-02-27 | Шлюмбергер Текнолоджи Б.В. | Моделирование взаимодействия трещин гидравлического разрыва в системах сложных трещин |
WO2016032489A1 (en) * | 2014-08-28 | 2016-03-03 | Landmark Graphics Corporation | Optimizing multistage hydraulic fracturing design based on three-dimensional (3d) continuum damage mechanics |
WO2016153953A1 (en) * | 2015-03-20 | 2016-09-29 | Schlumberger Technology Corporation | Stacked height growth fracture modeling |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2745142C1 (ru) * | 2020-09-17 | 2021-03-22 | Общество с ограниченной ответственностью «Газпромнефть Научно-Технический Центр» | Способ и система моделирования трещин гидроразрыва пласта бесконечно-конечной проводимости и поперечно-продольного расположения относительно горизонтального ствола скважины |
RU2775034C1 (ru) * | 2021-07-28 | 2022-06-27 | Публичное акционерное общество "Газпром нефть" | Способ выбора оптимального дизайна гидроразрыва пласта на основе интеллектуального анализа полевых данных для увеличения добычи углеводородного сырья |
RU2786303C1 (ru) * | 2022-07-01 | 2022-12-19 | Общество с ограниченной ответственностью "Газпромнефть-Технологические партнерства" (ООО "Газпромнефть-Технологические партнерства") | Способ определения геометрии трещин ГРП с учетом разделения и взаимодействия потоков жидкости ГРП между трещинами |
RU2796158C1 (ru) * | 2022-08-18 | 2023-05-17 | Шлюмберже Текнолоджи Б.В. | Способ определения концентрации расклинивающего агента в жидкости гидроразрыва и способ выполнения гидроразрыва пласта |
Also Published As
Publication number | Publication date |
---|---|
CA3062854A1 (en) | 2019-11-29 |
WO2018203765A1 (ru) | 2018-11-08 |
US20200056460A1 (en) | 2020-02-20 |
US11499406B2 (en) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2730576C1 (ru) | Способ прогнозирования гидроразрыва пласта, способ гидроразрыва пласта, способы прогнозирования рисков гидроразрыва пласта | |
Xu et al. | Modeling dynamic behaviors of complex fractures in conventional reservoir simulators | |
Velikanov et al. | New fracture hydrodynamics and in-situ kinetics model supports comprehensive hydraulic fracture simulation | |
US11560776B2 (en) | Methods and systems of modeling fluid diversion treatment operations | |
US9175555B2 (en) | Fluid injection completion techniques | |
WO2020247621A1 (en) | Fracturing-fluid formula workflow | |
Jabbari et al. | Hydraulic fracturing design for horizontal wells in the Bakken formation | |
US11341298B2 (en) | Method for reservoir stimulation analysis and design based on lagrangian approach | |
Hwang et al. | A 3-dimensional fracture propagation model for long-term water injection | |
Abdrazakov et al. | Proppant transport capabilities of high-viscosity friction reducers at relatively low rates: Advanced modeling and field validation | |
Sharma et al. | Hydraulic fracturing design and 3D modeling: a case study from Cambay Shale and Eagleford Shale | |
McClure et al. | Applying a Combined Hydraulic Fracturing, Reservoir, and Wellbore Simulator: Staged Field Experiment# 3, Cluster Spacing, and Stacked Parent/Child Frac Hits | |
Chaplygin et al. | Case Studies of Re-Fracturing Achimov Reservoirs with High-Viscous Friction Reducer on Salym Group of Oilfields | |
Suppachoknirun | Evaluation of multi-stage hydraulic fracturing techniques to optimize production in naturally fractured reservoirs using a DFN-based numerical technique | |
Molina et al. | A computational fluid dynamics approach to predict pressure drop and flow behavior in the near wellbore region of a frac-packed gas well | |
Bajestani et al. | Effect of hydraulic fracture and natural fractures interaction in fracture propagation | |
Zheng | Development of a fully integrated equation of state compositional hydraulic fracturing and reservoir simulator | |
Acharya | Hydraulic-fracture-treatment design simulation | |
US20230034507A1 (en) | Formation stimulation with acid etching model | |
Almarri et al. | Numerical feasibility of near-wellbore cooling as a novel method for reducing breakdown pressure in hydraulic fracturing | |
Shin | Simultaneous propagation of multiple fractures in a horizontal well | |
Soomro et al. | Numerical Simulation on Controlling and Optimizing of Fracture Behaviors for Hydraulic Fracturing in Low Permeability Formation | |
Sepehrnoori et al. | An extension of the embedded discrete fracture model for modeling dynamic behaviors of complex fractures | |
Nguyen Huu | Fracture Treatment Design in the Lower Miocene Reservoir, Offshore Viet Nam | |
Hanna Menhem | The Impact Stress on the Fracture Conductivity and the Productivity of the Marcellus Shale Horizontal Well |