石油工程

苏里格砂泥薄互储层缝控压裂造缝机制及穿层判别准则

  • 谢锦阳 , 1, 2 ,
  • 侯冰 , 1, 2, 3 ,
  • 何明舫 4 ,
  • 刘欣佳 4 ,
  • 魏靖依 5
展开
  • 1 中国石油大学(北京)油气资源与工程全国重点实验室,北京 102249
  • 2 中国石油大学(北京)石油工程教育部重点实验室,北京 102249
  • 3 中国石油大学(北京)克拉玛依校区石油学院,新疆克拉玛依 834000
  • 4 中国石油长庆油田公司油气工艺研究院,西安 710021
  • 5 西南石油大学,成都 610500
侯冰(1979-),男,辽宁北镇人,博士,中国石油大学(北京)克拉玛依校区石油学院教授,目前主要从事石油工程岩石力学、水力压裂等方面的研究工作。地址:新疆克拉玛依市克拉玛依区安定路355号,中国石油大学(北京)克拉玛依校区石油学院,邮政编码:834000。E-mail:

谢锦阳(1999-),男,辽宁葫芦岛人,中国石油大学(北京)在读博士生,目前主要从事多岩性薄互层穿层压裂及页岩分布式光纤监测水力压裂等方面的研究工作。地址:北京市昌平区府学路18号,中国石油大学(北京)石油工程学院,邮政编码:102249。E-mail:

Copy editor: 唐俊伟

收稿日期: 2024-06-03

  修回日期: 2024-09-25

  网络出版日期: 2024-10-15

基金资助

国家重点研发计划-中美政府间国际合作项目(2022YFE0129800)

中国石油天然气集团有限公司-中国石油大学(北京)战略合作科技专项(ZLZX2020-02)

Fracture-controlled fracturing mechanism and penetration discrimination criteria for thin sand-mud interbedded reservoirs in Sulige gas field, Ordos Basin, China

  • XIE Jinyang , 1, 2 ,
  • HOU Bing , 1, 2, 3 ,
  • HE Mingfang 4 ,
  • LIU Xinjia 4 ,
  • WEI Jingyi 5
Expand
  • 1 National Key Laboratory of Petroleum Resources and Engineering, China University of Petroleum, Beijing 102249, China
  • 2 MOE Key Laboratory of Petroleum Engineering, China University of Petroleum, Beijing 102249, China
  • 3 College of Petroleum Engineering, China University of Petroleum-Beijing at Karamay, Karamay 834000, China
  • 4 Oil& Gas Technology Research Institute, PetroChina Changqing Oilfield Company, Xi’an 710021, China
  • 5 Southwest Petroleum University, Chengdu 610500, China

Received date: 2024-06-03

  Revised date: 2024-09-25

  Online published: 2024-10-15

摘要

针对苏里格气田二叠系盒8段砂泥薄互储层压裂施工中裂缝穿层扩展判别与压裂效果评价存在的问题,建立考虑层间非均质性的砂泥薄互储层地质力学模型,开展水力裂缝穿层扩展实验,揭示砂泥薄互储层中人工裂缝的起裂-延伸-交互-穿层机制。基于非常规裂缝模型的数值模拟,阐明了裂缝在薄互层中的垂向起裂与扩展特征,构建了裂缝穿层的判别准则与压裂效果评价方法。研究表明:层间应力差是直接影响穿层压裂裂缝形态的主要地质因素,苏里格气田层间应力差异系数(ξIL)小于0.40时,裂缝可突破隔层在目标砂岩层内扩展;0.40≤ξIL<0.45时,裂缝突破隔层但无法在目标砂层内有效扩展;ξIL≥0.45时,裂缝仅在射孔储层扩展,无法实现穿层扩展。提高压裂液黏度与排量能弥补能量损耗并突破隔层限制,苏里格气田压裂施工参数控制区域内采用高黏度(50~100 mPa·s)压裂液、高排量(12~18 m3/min)施工条件有利于裂缝穿层。

本文引用格式

谢锦阳 , 侯冰 , 何明舫 , 刘欣佳 , 魏靖依 . 苏里格砂泥薄互储层缝控压裂造缝机制及穿层判别准则[J]. 石油勘探与开发, 2024 , 51(5) : 1150 -1159 . DOI: 10.11698/PED.20240366

Abstract

Considering the problems in the discrimination of fracture penetration and the evaluation of fracturing performance in the stimulation of thin sand-mud interbedded reservoirs in the eighth member of Shihezi Formation of Permian (He-8 Member) in the Sulige gas field, a geomechanical model of thin sand-mud interbedded reservoirs considering interlayer heterogeneity was established. The experiment of hydraulic fracture penetration was performed to reveal the mechanism of initiation-extension-interaction-penetration of hydraulic fractures in the thin sand-mud interbedded reservoirs. The unconventional fracture model was used to clarify the vertical initiation and extension characteristics of fractures in thin interbedded reservoirs through numerical simulation. The fracture penetration discrimination criterion and the fracturing performance evaluation method were developed. The results show that the interlayer stress difference is the main geological factor that directly affects the fracture morphology during hydraulic fracturing. When the interlayer stress difference coefficient is less than 0.4 in the Sulige gas field, the fractures can penetrate the barrier and extend in the target sandstone layer. When the interlayer stress difference coefficient is not less than 0.4 and less than 0.45, the factures can penetrate the barrier but cannot extend in the target sandstone layers. When the interlayer stress difference coefficient is greater than 0.45, the fractures only extend in the perforated reservoir, but not penetrate the layers. Increasing the viscosity and pump rates of the fracturing fluid can compensate for the energy loss and break through the barrier limit. The injection of high viscosity (50-100 mPa·s) fracturing fluid at high pump rates (12-18 m3/min) is conducive to fracture penetration in the thin sand-mud interbedded reservoirs in the Sulige gas field.

0 引言

鄂尔多斯盆地苏里格气田二叠系盒8段为典型的致密砂岩储层,随着开发逐渐转向砂泥薄层交互严重的边缘储层,其纵向小层高效动用难度大,水平井钻遇率低,单井产能差,人工裂缝起裂与延伸困难。实现水力裂缝纵向可控沟通目的层,是水平井穿层压裂设计优化的关键问题。目前穿层压裂技术主要应用于煤层气与页岩油等非常规油气储层的开发,其同样具有纵向多产层、厚度不均、层内与层间强非均质的特征[1]。国内外学者对页岩油储层水力裂缝的穿层扩展特征、界面强度、压裂施工参数等方面进行了大量物理模拟与数值模拟研究[2-4],在理论与应用方面均取得了很好的效果。苏里格气田砂泥薄互储层具有纵向上砂层发育分散、泥岩夹层薄等特征,地质条件复杂,如果将页岩油穿层压裂基础理论应用于该气田的开发,必须明确苏里格气田裂缝垂向扩展的影响规律,对穿层效果进行量化评价。
针对苏里格气田盒8段的砂泥多薄互叠置储层,本文将大尺寸真三轴水力压裂模拟实验与基于非常规裂缝模型的数值模拟(简称非常规裂缝数值模拟)相结合,充分考虑实际砂泥薄互储层组合特点,以地质力学参数为基础,通过优化关键压裂施工参数,揭示了裂缝在多层系储层中的垂向扩展规律,实现了穿层效果的量化评价,绘制了苏里格气田水平井穿层压裂施工参数优化设计图版,实现了一次改造提升产量的目标,可为现场压裂施工提供技术支持。

1 砂泥薄互储层地质力学模型

苏里格气田主力气层为盒8段致密砂岩,上段为棕褐色泥岩夹浅灰色细砂岩、浅灰色泥质砂岩;下段以浅灰色、灰白色中砂岩为主,夹浅灰色泥质砂岩、深绿灰色泥岩。测井解释视孔隙度最大为13.04%,最小为1.38%;视渗透率最大为1.67×10−3 μm2,最小为0.01×10−3 μm2,物性差异非常大,纵向上难以有效沟通多个储层,严重影响穿层压裂改造效果。
采集井下岩心6块开展压缩破坏实验,其主要层位为盒8段及山西组,深度2 911.30~3 860.53 m,经测试岩石力学参数如表1所示。图1为致密砂岩与泥岩岩心压缩破坏实验前后对比照片。单轴压缩条件下,由于砂岩本身脆性强,1号砂岩岩心在不施加围压时,破坏形态同时呈现拉张破坏与剪切破坏;2号砂岩岩心施加围压20 MPa时,破坏形态过渡为剪切破坏;3号与4号砂岩岩心分别施加围压25,30 MPa时,剪切破坏特征愈发明显;5号泥岩岩心不施加围压时呈现典型拉张破坏现象;6号泥岩岩心施加围压20 MPa时,破坏形态过渡为剪切破坏。
表1 岩石力学参数测试结果
编号 岩性 围压/
MPa
弹性模量/
GPa
泊松比 抗压强度/
MPa
黏聚力/
MPa
内摩擦角/
(°)
1 砂岩 0 23.17 0.35 70.02 19.5 43.0
2 20 36.96 0.30 88.61
3 25 56.77 0.28 157.50
4 30 60.51 0.26 209.88
5 泥岩 0 56.04 0.25 75.92 35.5 20.3
6 20 56.98 0.24 101.36
图1 致密砂岩与泥岩试样变形破坏特征
建立层间强非均质性测井解释模型,需明确动态岩石力学参数。根据岩石力学基础理论,采用现场声波测井资料与实验室静态弹性模量、静态泊松比数据,进行动静参数的转换与拟合[5],可得动态与静态弹性模量间、动态与静态泊松比间关系式。依据苏里格气田盒8段岩石力学参数测试结果,采用最小二乘法原理拟合得动、静态弹性模量的关系式为(相关系数0.895 3):
${{E}_{\text{s}}}=0.653\ 1{{E}_{\text{d}}}-3.233\ 9$
动、静态泊松比的关系式为(相关系数0.880 2):
${{\nu }_{\text{s}}}=0.163\ 4+0.312\ 4{{\nu }_{\text{d}}}$
苏里格气田盒8段不同岩性地层叠置,采用传统地应力计算模型进行计算误差较大,难以准确体现层间应力的差异性。因此,充分考虑盒8段砂泥岩层间差异以及沉积特点,引入泥质含量作为影响因素,重新构建层间强非均质性地质力学模型。
泥质含量较正公式为:
$f\left( GR \right)=\left\{ \begin{array}{*{35}{l}} 0.8\ \ \ \ GR50 \\ 0.9\ \ \ \ 50\le GR100 \\ 1.0\ \ \ \ 100\le GR150 \\ 1.2\ \ \ \ 150\le GR200 \\ \end{array} \right.$
强非均质性地质力学模型为:
$\left\{ \begin{align} & {{\sigma }_{\text{H}}}=\left[ M+\frac{1}{2\left( 1-\nu _{\text{d}}^{\text{2}} \right)}\left( K+L \right) \right]f\left( GR \right) \\ & {{\sigma }_{\text{h}}}=\left[ M+\frac{1}{2\left( 1-\nu _{\text{d}}^{\text{2}} \right)}\left( K-L \right) \right]f\left( GR \right) \\ \end{align} \right.$
其中
$\begin{align} & K=0.013{{E}_{\text{d}}}\left( 1+{{\nu }_{\text{d}}} \right) \\ & M\text{=}\frac{{{\nu }_{\text{d}}}}{1-{{\nu }_{\text{d}}}}\left( {{\sigma }_{\text{v}}}-\alpha {{p}_{\text{p}}} \right)+\alpha {{p}_{\text{p}}} \\ & L=0.01{{E}_{\text{d}}}\left( 1-{{\nu }_{\text{d}}} \right) \\ \end{align}$
按照层间应力差大小划分出易穿层区域(0~3 MPa)、压裂施工参数控制区域(3~6 MPa)和难穿层区域(大于6 MPa)。提取苏里格气田某区块盒8段砂-泥岩层界(深度3 300 m)的层间水平最小地应力差,利用高斯-克里金插值法重构该区块层间应力差分布图(见图2)。其中压裂施工参数控制区域占80.36%,难穿层区域占10.28%,易穿层区域占9.36%,由此可见穿透多层扩展难度大,80%以上区域需要优选压裂施工参数。因此,充分提高纵向改造程度,必须厘清盒8段多岩性叠置储层的压裂造缝机制,明确多薄互层裂缝垂向起裂与扩展形态特征与层间应力差的关系。
图2 盒8段3 300 m处砂-泥岩井区层间应力差分布图

2 单向薄夹层缝控压裂造缝机制

2.1 真三轴水力压裂裂缝形态物理模拟实验

为揭示砂泥薄互储层水力裂缝起裂-延伸-交互-穿层机制,采用真三轴水力压裂物理模拟系统开展原位环境下砂泥岩薄互层水力裂缝穿层扩展实验。采集苏里格气田二叠系盒8段砂岩露头、泥岩露头各1块加工成砂岩和泥岩岩板并按图3所示尺寸制作叠置试样1块(记为1#试样),同时采集现场砂泥叠置露头岩样1块(记为2#试样),分别开展高层间应力差(1#试样)与低层间应力差(2#试样)条件下砂泥薄互储层水力压裂物理模拟实验,试样尺寸为300 mm×300 mm×300 mm。实验步骤为:①用砂纸打磨加工好的4块长方体薄岩板面,然后按图3a方式拼接成试样。②在图3a所示砂岩储层岩板中心位置钻长15.0 cm、直径2.5 cm的沉孔,孔底部用食盐填充以模拟裸眼完井。③将直径2.5 cm的橡胶塞固定在沉孔中心,沿橡胶塞外壁面缓慢灌入混凝土,待凝固后取掉胶塞形成井眼。④将长13.0 cm、直径1.5 cm的井筒用固井胶封固在井眼内,模拟单簇射孔水力裂缝穿透泥岩层扩展至目标砂岩层的压裂实验。⑤将试样如图3b所示放置于真三轴原位应力测试围压釜后,加载三向应力,注入压裂液开展压裂实验。⑥实验结束后拆卸设备,取出并打开试样,使用3D扫描设备获取裂缝面信息,实验结束。
图3 砂泥薄层穿层压裂物理模拟试样及测试装置
实验参数采用三维水力压裂模型相似准则[6]确定,考虑几何条件、介质条件、边界条件等设置相似准则系数。经计算,设定实验排量36 mL/min;1#试样垂向主应力20 MPa,水平最大主应力15 MPa,水平最小主应力10 MPa,水平最小主应力仅加载于泥岩层,以模拟现场层间应力差10 MPa的条件;2#试样垂向主应力20 MPa,水平最大主应力15 MPa,水平最小主应力5 MPa,水平最小主应力仅加载于泥岩层,以模拟现场层间应力差5 MPa的条件。实验过程中压裂液加入绿色荧光粉作为示踪剂便于直接观察裂缝激活区域。

2.2 砂泥岩试样水力压裂裂缝穿层扩展特征

采用三维扫描成像设备对裂缝面进行扫描,并对三维裂缝面数据进行数值处理分析,得压裂裂缝面三维重构图,可直观地观察层间应力差影响下的穿层压裂裂缝形态,明确压裂造缝机制。
图4为高层间应力差条件下1#试样压裂测试后试样展开图。可以看到,裂缝在砂岩板射孔孔眼起裂并垂向扩展,遇层界面后受高层间应力差影响,仅能垂直于井筒方向延伸。如②号展开图所示,岩样边缘见大量绿色荧光,表明当裂缝延伸至砂岩层边界后发生压裂液漏失;而黄色圆圈标注的裂缝位置未见绿色荧光,表明裂缝并非是穿层扩展至此处。
图4 1#试样在10 MPa层间应力差下穿层裂缝形态
图5为1#试样压裂裂缝面照片与三维重构图,图5a中红色虚线区域未见绿色荧光表明水力裂缝未能穿透层界面;图5b中水力裂缝在层界面上发生钝化现象,形成“一”字钝化缝,裂缝形态简单。1#试样水力压裂测试泵压曲线为典型的两翼型水平主裂缝起裂-扩展特征(见图6)。实验结果表明高应力差条件下,裂缝难以穿层形成有效导流通道,因此钻井过程中应尽量避免钻遇高层间应力差段或水力压裂施时采用其他方法实现穿层。
图5 1#试样压裂裂缝面照片与三维重构图
图6 1#试样压裂泵压曲线
图7为低层间应力差条件下2#试样压裂测试后试样展开图。可以看到,砂泥交互叠置露头表面存在大量天然裂缝,水力裂缝在射孔孔眼处起裂,裂缝纵向扩展至层界面处受阻,随后在近井筒段发生拉伸破坏形成水平主缝(见图8a中Ⅱ号裂缝),水力裂缝向上扩展过程中激活上部层间缝后,裂缝又发生转向扩展。经短暂憋压后,裂缝穿透泥岩隔层并激活层界面发生转向,最终其下部与Ⅲ号天然裂缝沟通(见图8b),形成了复杂缝网。重构裂缝面后发现,受层间应力差影响,垂直缝扩展过程中呈现非均匀扩展现象,同时激活多个层间界面形成粗糙、形态复杂的裂缝网络。
图7 2#试样在5 MPa层间应力差下穿层裂缝形态
图8 2#试样压裂裂缝面照片与三维重构图
2#试样水力压裂测试泵压曲线(见图9)整体呈憋压后水平主缝起裂、憋压后近井层界面起裂、憋压后远井层界面起裂“三段式”模式,破裂压力按12.5,8.0,5.0 MPa逐级降低,主裂缝明显呈“工”字型裂缝组合特征。
图9 2#试样压裂泵压曲线
实验结果表明低应力差条件下裂缝垂向扩展能力有效提升,但受层界面影响裂缝会发生转向,限制裂缝进一步垂向扩展。建议在水力压裂施工过程中先采用高黏度压裂液优先保证主裂缝有效扩展。而后采用低黏度压裂液沟通天然裂缝,有效提高远井地带裂缝复杂程度,增加改造体积,实现“控近扩远”储层改造[7]

3 缝控压裂非常规裂缝数值模拟

3.1 非常规裂缝模型数值计算方法

为揭示层间应力差条件下裂缝在多薄砂泥互层中垂向起裂与扩展形态特征,基于非常规裂缝模型,采用Kinetix方法[8]建立流固耦合模型。该方法充分考虑了储层非均质及应力各向异性,可以更加精细地模拟复杂裂缝形态。

3.1.1 水力压裂基本控制方程

假设流体沿裂缝扩展方向流动,满足基于局部质量守恒定律的不可压缩流体连续性方程:
$\frac{\partial q}{\partial s}+\frac{\partial \left( {{H}_{\text{f}}}\bar{w} \right)}{\partial t}+{{q}_{\text{l}}}=0$
其中压裂液渗流速度可根据卡特滤失模型获得:
${{q}_{\text{l}}}=2{{h}_{\text{l}}}{{u}_{\text{l}}}$
压裂泵入流体总体积等于裂缝中的流体体积与从裂缝中漏失的流体体积之和,可得总流体的体积平衡方程为:
$\int_{0}^{t}{Q\left( t \right)}\text{d}t=\int_{0}^{L(t)}{h}\bar{w}\text{d}s+\int_{{{H}_{\text{l}}}}{\int_{0}^{t}{\int_{0}^{L(t)}{2}}}{{q}_{\text{l}}}\text{d}s\text{d}t\text{d}{{h}_{\text{l}}}$
裂缝内幂律流体满足泊肃叶定律动量方程[9]
$\frac{\partial p}{\partial s}=-{{\alpha }_{0}}\frac{1}{{{{\bar{w}}}^{2{n}'+1}}}\frac{q}{{{H}_{\text{f}}}}{{\left| \frac{q}{{{H}_{\text{f}}}} \right|}^{{n}'-1}}$
其中
${{\alpha }_{0}}=\frac{2{K}'}{\phi {{\left( {{n}'} \right)}^{{{n}'}}}}{{\left( \frac{4{n}'+2}{{{n}'}} \right)}^{{{n}'}}}$
$\phi \left( {{n}'} \right)=\frac{1}{{{H}_{\text{f}}}}\int_{{{H}_{\text{f}}}}{{{\left[ \frac{w(z)}{{\bar{w}}} \right]}^{\frac{2{n}'+1}{{{n}'}}}}}\text{d}z$
为了提高工程应用计算效率,假设缝宽为二维平面应变,由弹性方程与流体压力的相关关系可得:
$w\left( x,y,z \right)=w\left[ p\left( x,y \right),h,z \right]$
同时在裂缝尖端需要满足如下边界条件:
$\left\{ \begin{array}{*{35}{l}} p_{\text{n,t }}^{{}}\left( z \right)=0 \\ q=0 \\ w=0 \\ \end{array} \right.$

3.1.2 基于均衡方程的多岩性叠置储层裂缝扩展模型

裂缝高度可以通过对比裂缝尖端处的应力强度因子与断裂尖端处地层的断裂韧性来确定。若任一尖端的应力强度因子大于地层的断裂韧性,则缝高继续扩展,反之缝高则停止扩展。传统拟三维模型(见图10a)仅考虑了上、下裂缝尖端附近垂直方向上流体压力梯度的影响,导致缝高被高估。叠置式裂缝扩展模型[10]将每个单元分割为垂直连接的元素单独计算,在层位应力差较大条件下,当水力裂缝越过岩性界面时产生新的元素,会将邻层中心设为新起裂点,将相邻元素作为新的独立单元扩展计算(见图10b),更能精确解释层间应力非均质性的影响。
图10 拟三维模型与叠置式裂缝扩展模型水力裂缝垂向扩展示意图
基于水力能量动态守恒,考虑岩石性质对缝高扩展的阻碍作用[11],保证裂缝扩展时能量供给速率不小于能量消耗速率,得水力裂缝穿过非均质邻层所需最小流体压力为:
$\Delta {{p}_{\text{e}}}=\frac{{{E}_{\text{a}}}}{{{h}_{\text{a}}}S}\left[ \frac{K_{\text{Ic}\left( \text{b} \right)}^{2}}{{{E}_{\text{b}}}}-\frac{K_{\text{Ic}\left( \text{a} \right)}^{2}}{{{E}_{\text{a}}}} \right]$
其中
$S=\frac{1}{h}\int_{0}^{h}{\sqrt{\frac{h-z}{z}}}p_{\text{n,t }}^{{}}\left( z \right)\text{d}z$
叠置式模型顶、底尖端应力强度因子计算公式为:
${{K}_{\text{Iu}}}=\sqrt{\frac{2}{\text{ }\!\!\pi\!\!\text{ }h}}\sum\limits_{j=1}^{{{N}_{\text{s }}}}{\int_{{{H}_{j-1}}}^{{{H}_{j}}}{\left[ {{p}_{\text{d},j}}+{{\rho }_{\text{f},j}}g\left( {{h}_{\text{d},j}}-a \right)-\sigma \left( a \right) \right]}}\sqrt{\frac{a}{h-a}}\text{d}a$
${{K}_{\text{Il}}}=\sqrt{\frac{2}{\text{ }\!\!\pi\!\!\text{ }h}}\sum\limits_{j=1}^{{{N}_{\text{s }}}}{\int_{{{H}_{j-1}}}^{{{H}_{j}}}{\left[ {{p}_{\text{d},j}}+{{\rho }_{\text{f},j}}g\left( {{h}_{\text{d},j}}-a \right)-\sigma \left( a \right) \right]}}\sqrt{\frac{h-a}{a}}\text{d}a$
由于叠置元素的顶部或底部不一定是断裂尖端,因此需要考虑相邻元素应用叠加原理计算裂缝宽度:
$w\left( h,z \right)=\frac{4\left( 1-{{v}^{2}} \right)}{\text{ }\!\!\pi\!\!\text{ }E}\sum\limits_{j=1}^{{{N}_{\text{s}}}}{\int_{{{H}_{j-1}}}^{{{H}_{j}}}{\left[ {{p}_{\text{d},j}}+{{\rho }_{\text{f},j}} \right.g\left( {{h}_{\text{d},j}}-a \right)}}-\left. \sigma \left( a \right) \right]{{\cosh }^{-1}}\frac{z\left( \frac{h-2a}{h} \right)+a}{|z-a|}\text{d}a$

3.2 缝控压裂数值模型及模拟参数设置

为更好地定量表征水力裂缝在砂泥岩互层条件下的穿层形态特征,采用砂泥岩单向薄互层堆叠模型(见图11)进行模拟。结合现场测井数据以及岩石力学参数测试结果,设置模型参数如表2所示。
图11 非常规裂缝模型建模示意图
表2 模型储层应力参数及岩石力学参数设置表
储层岩性 垂向主应力/
MPa
水平最大主应力/
MPa
水平最小主应力/
MPa
弹性模量/
GPa
泊松比 抗拉强度/
MPa
断裂韧性/
(Pa•m0.5
孔隙度/
%
渗透率/
10−3 μm2
砂岩 60 45~51 35~41 36 0.30 9.0 4 6.8 0.47
泥岩 60 53 43 56 0.25 7.5 3 5.1 0.01

3.3 层间应力差影响下的穿层压裂造缝机制

苏里格气田砂泥岩互层段纵向上水平应力差异为1~8 MPa,通过对数值模拟结果进行分析发现层间应力差对水力裂缝穿层扩展形态影响显著,在不同层间应力差下,主要呈现4种裂缝形态组合。为进一步量化这种影响,定义层间应力差异系数为:
${{\xi }_{\text{IL}}}=\frac{{{\sigma }_{\text{mh}}}-{{\sigma }_{\text{sh}}}}{{{\sigma }_{\text{v}}}-{{\sigma }_{\text{sh}}}}$
缝宽变化对支撑剂运移铺置影响显著,缝宽减小会使支撑剂在粗糙壁面裂缝内堆积,有效裂缝长度较小[12]。因此本文将裂缝宽度大于1 mm作为有效裂缝扩展区的边界以贴合实际工程问题。

3.3.1 “葫芦”型裂缝

以上述数值模型为基础,设定层间应力差异系数ξIL<0.2,模拟压裂时的裂缝扩展形态(见图12)。可以看到,由于层间应力差较小且垂向主应力为最大主应力,压裂液容易穿透夹层,激活区域边缘平滑,未发生非连续性扩展现象。由于砂岩孔隙度与渗透率高于泥岩,水力裂缝更趋向于在砂岩层段内延伸;同时受射孔位置影响,射孔砂岩层张开裂缝缝长高于目标砂岩层,下部穿层段中,在水力能量优势作用下裂缝仍能继续扩展穿层,最终形成非对称水力裂缝激活区域,裂缝整体呈现“葫芦”型。
图12 ξIL<0.2时裂缝扩展形态

3.3.2 “工”字缝

设定层间应力差异系数0.20≤ξIL<0.40,模拟压裂时的裂缝扩展形态(见图13)。由于随着层间应力差的逐渐增大,压裂液穿透夹层所需能量逐渐增高,表现为穿层缝在泥岩夹层的缝宽减小。由于水力裂缝的扩展遵循低流动阻力、低能量耗散原则,在应力差的作用下,裂缝在泥岩层向目标砂岩层扩展所需能量小于水平扩展,最终发生非连续性扩展,上部穿层裂缝与射孔段裂缝整体形态呈“工”字缝。
图13 0.20≤ξIL<0.40时裂缝扩展形态

3.3.3 “T”型缝

设定层间应力差异系数0.40≤ξIL<0.45,模拟压裂时裂缝的扩展形态(见图14)。在此阶段层间应力差较大,因此水平主裂缝沿最大水平主应力方向延伸更明显。当裂缝穿透泥岩夹隔层进入目标砂岩层后,裂缝发生转向后水平延伸,但延伸距离较短,最终在闭合压力作用下,泥岩隔层中的穿层缝大部分闭合,缝宽减少,导致裂缝形态整体呈“T”型,此类穿层缝相比之下难以有效沟通目标砂岩层。
图14 0.40≤ξIL<0.45时裂缝扩展形态

3.3.4 “一”字钝化缝

同样设定层间应力差异系数ξIL≥0.45,模拟压裂时裂缝的扩展形态(见图15)。由于高层间应力差作用较强,水平主裂缝主要沿最大水平主应力方向延伸,所形成裂缝的长度最长,同时裂缝无法穿透隔层,在层界面产生钝化效应,最终呈现“一”字钝化缝。
图15 ξIL≥0.45时裂缝扩展形态
综上可知,当层间应力差作用较强时,裂缝一般无法穿透隔层或穿透隔层的能力较差,应考虑采用如径向井等措施尝试突破隔层,实现裂缝沟通。

4 穿层压裂缝高扩展判别准则及图版

为实现压裂施工参数控制区域(见图2)裂缝有效穿层扩展,优化压裂施工参数,以层间应力差6.0 MPa单向砂泥薄互层段为研究对象,设置施工排量为6~18 m3/min、压裂液黏度为10~100 mPa·s,基本可以涵盖目前压裂施工的参数范围。为方便裂缝穿层判别,将水力裂缝穿层施工参数阈值(以下简称“施工阈值”)定义为裂缝穿透泥岩隔层顶部并开始在目标砂岩储层扩展的最低参数。

4.1 缝控压裂缝高连通率穿层判别准则

由于单向砂泥薄互层中,目标砂岩储层位于射孔储层上部,故研究重点为上半缝区域。定义缝控压裂缝高连通率公式为(负值计为0):
${{H}_{\text{c}}}=\frac{{{H}_{\text{a}}}-{{H}_{\text{m}}}}{{{H}_{\text{p}}}+{{H}_{\text{t}}}}\times 100\%$
带入各参数可得该模型条件下其穿层判别准则为缝控压裂缝高连通率大于50%。以上述设定排量、压裂液黏度施工参数为基础,模拟缝控压裂缝高连通率的变化情况,绘制成三维图版(见图16)。图中缝控压裂缝高连通率为0~50%的区域表示水力裂缝在泥岩隔层穿层扩展,形成“T”型裂缝,但未能有效实现穿层;缝控压裂缝高连通率高于50%的红色区域表示水力裂缝已经穿透泥岩隔层进入目标砂岩储层,形成“工”字裂缝。
图16 缝控压裂缝高连通率穿层判别三维图版

4.2 储层改造面积有效率穿层判别准则

侯冰等[13]基于页岩水力压裂实验提出裂缝沟通面积(SRA)作为水力压裂效果的评价指标。储层改造后,SRA值越大,增产稳产效果越好。为更好评价压裂效果,进一步定义储层改造面积有效率为有效裂缝面积与总改造面积之比:
${{S}_{\text{e}}}=\frac{{{S}_{\text{s}}}}{{{S}_{\text{s}}}+{{S}_{\text{m}}}}\times 100\%$
同样以上述设定施工参数为基础,模拟储层改造面积有效率的变化情况并绘制成三维图版(见图17)。由图可以看到,图左侧储层改造面积有效率为70%~82%的红色区域表示裂缝难以穿过泥岩隔层,裂缝仅在射孔储层内扩展并沟通产层;随着压裂液黏度及施工排量的增加,裂缝在泥岩隔层中扩展,但该部分为非目标产层,因此储层改造面积有效率快速下降;压裂液黏度及施工排量进一步增加,裂缝穿过泥岩隔层进入目标砂岩层中扩展并沟通产层,故储层改造面积有效率缓慢上升。可见穿层演化过程中储层改造面积有效率存在拐点,该拐点即为穿层阈值点。压裂裂缝是否穿层,判别准则可定义为储层改造面积有效率由降转增的时间节点。
图17 储层改造面积有效率穿层判别三维图版

4.3 储层改造体积有效率穿层效果评价

图17显示随着压裂液黏度及施工排量继续增加,由于沟通了目标储层,储层改造面积有效率出现逐渐回升的趋势,因此在给定大于穿层阈值的任意储层改造面积有效率值的情况下无法判别其处于未穿层或是穿层状态。同时由于未考虑采用高黏压裂液压裂时裂缝扩展引起的缝宽变化,该方法仅能作为一种穿层判别准则,还亟需一种新的穿层压裂效果评价方法。因此定义一种考虑缝宽影响的储层改造体积有效率作为穿层压裂的效果评价方法:
${{V}_{\text{e}}}=\frac{{{S}_{\text{s}}}}{{{S}_{\text{s}}}+{{S}_{\text{m}}}}\frac{{{{\bar{w}}}_{n}}}{{{{\bar{w}}}_{\max }}}\times 100\%$
该方法考虑了裂缝宽度系数影响,以此绘制可实现裂缝穿层压裂效果评价的图版(见图18)。由图可知在相同施工排量下,储层改造体积有效率随压裂液黏度的增加而增加;图中从黄色过渡至红色的区域储层改造体积有效率为50%~64%,为裂缝有效穿层区域,关键施工参数可在此范围进行优化。
图18 储层改造体积有效率穿层压裂效果评价图版
对比缝控压裂缝高连通率与储层改造面积有效率穿层判别准则发现两者的关键施工参数阈值相同,适应性较好,均可以作为苏里格气田盒8段砂泥薄互层的穿层判别准则。缝控压裂缝高连通率穿层判别准则可与水力压裂电磁监测有效结合,描述裂缝几何形状和空间展布,可量化穿层压裂效果,实现可控穿层。储层改造面积有效率穿层判别准则及储层改造体积有效率穿层效果评价方法可与微地震监测结果有效结合,适用于多岩性叠置的非常规油气储层水力压裂改造效果监测,实现穿层判别与压裂施工参数优化。同时综合上述穿层判别准则与穿层效果评价图版可以发现,黏度越高(50~100 mPa·s)、排量越大(12~18 m3/min),改造储层的效果越好,穿层上半缝高最高可达13 m,储层改造体积有效率最高可达64%。

5 穿层压裂效果现场验证

苏里格气田苏A井水平段井眼轨迹位于盒8x2砂岩层,埋藏深度3 300~3 312 m,砂体厚度12 m,孔隙度6%,渗透率0.3×10-3 μm2,含气饱和度30%,储层品质较差,属于差气层。为提高单井产量,设计将该层上方厚度5 m的气层同步改造,增大单井控制区域。两个气层间间隔了5 m的泥岩夹层,泥岩GR值为120~160 API,利用式(3)、式(4)及式(18)计算层间应力差异系数约为0.28,具有穿层压裂可行性。
为了形成“工”字型裂缝动用纵向资源,采用单段单簇180°垂向定向射孔工艺诱导裂缝起裂方向。考虑实际情况将关键施工参数优化为排量16 m3/min、压裂液黏度80 mPa·s。压裂施工过程中,采用广域电磁裂缝监测技术实时监测裂缝扩展动态,电磁监测压裂液三维波及范围监测结果图19所示。压裂液最大波及高度为23 m,缝控压裂缝高连通率约为63%,高于其穿层阈值(55%);结合微地震事件点解释测算,压裂液波及总面积3 334 m2,压裂过程中储层改造面积有效率随时间推移出现明显拐点,最终有效率约为58%,两种方法均证明裂缝成功穿透上侧厚度5 m的泥岩夹层与上部气层有效沟通,实现了人工裂缝的穿层改造目标。压裂液波及总体积32 866 m3,其储层改造体积有效率约为64%,评价为裂缝有效穿层。试采结果与储层品质相近的邻井对比,苏A井平均日产气3.2×104 m3,较邻井提产12%,穿层压裂效果良好,同时证明研究成果具有较好的实用性。
图19 电磁监测压裂液三维波及图

6 结论

层间应力差是直接影响穿层压裂裂缝形态的主要地质因素,层间应力差异系数小于0.40时,裂缝可突破隔层在目标砂岩层内扩展,可形成“葫芦”型或“工”字型穿层裂缝形态;层间应力差异系数大于0.40且小于0.45时,裂缝突破隔层但缝宽较小,未在目标储层内充分扩展,裂缝整体形态呈“T”型裂缝,难以形成高效导流通道;层间应力差异系数大于0.45时,裂缝仅在射孔储层扩展,形成钝化缝,无法实现穿层扩展。
控制压裂施工中压裂液黏度和施工排量,可有效实现缝控穿层。在苏里格气田二叠系盒8段采用高黏度(50~100 mPa·s)压裂液,高排量(12~18 m3/min)施工能提供充足的裂缝扩展能量,有利于裂缝穿层,进一步有效改造油藏。
经现场应用检验证实缝控压裂缝高连通率穿层判别准则、储层改造面积有效率穿层判别准则与储层改造体积有效率穿层效果评价方法具有较好的实用性,可为现场压裂施工提供技术支撑。
符号注释:
a——元素底部到元素中某个点的距离,m;E——弹性模量,Pa;EaEb——射孔层段、隔层段弹性模量,Pa;EdEs——动、静态弹性模量,Pa;fGR)——泥质含量影响因子,无因次;GR——自然伽马,API;g——重力加速度,m/s2h——裂缝高度,m;ha——射孔层段厚度,m;hcp——裂缝底部到射孔中心的裂缝高度,m;hcp1hcp2——裂缝底部到1号、2号元素起裂中心的高度,m;hd,j——第j个元素的裂缝底部到产油层中心高度,m;hl——流体滤失高度,m;H1H2——裂缝底部到元素1、元素2顶部的高度,m;Ha——水力裂缝上半缝高,m;Hc——缝控压裂缝高连通率,%;Hf——裂缝内流体高度,m;Hl——总漏失区域高度,m;Hm——上层泥岩厚度,m;Hp——射孔中心到储隔层界面的距离,m;Ht——目标砂岩层的厚度,m;i——层位编号;j——元素编号;Klc1Klc2——第1,2层的断裂韧性,Pa·m0.5K°——稠度系数,Pa·snKIuKIl——叠置式裂缝扩展模型裂缝顶、底部端点应力强度因子,Pa·m0.5KIc(a)KIc(b)——射孔层、隔层段断裂韧性,Pa·m0.5L(t)——裂缝长度,m;MKL——中间变量,Pa;n°——流体幂律指数,无因次;Ns——模拟裂缝截面的元素个数;p(xy)——位置s的流体压力,Pa;p——流体压力,Pa;pd,j——第j个元素起裂中心的流体压力,Pa;Δpe——穿层所需流体压力,Pa;pn,t(z)——裂逢尖端净压力,Pa;pp——孔隙压力,Pa;q——缝长方向体积流量,m3/s;ql——单位长度裂缝中流体滤失体积速度,m2/s;Q(t)——总体积流量,m3/s;s——裂缝中某点到裂缝端点的距离,m;S——中间变量;Se——储层改造面积有效率,%;Sm——泥岩隔层裂缝沟通面积,m2Ss——砂岩储层裂缝沟通面积,m2t——泵注时间,s;ul——滤失速度,m/s;Ve——储层改造体积有效率,%;w——裂缝尖端宽度,m;w(z)——位置s处的裂缝缝宽,m;w(hz)——考虑最小水平应力及流体压力条件下裂缝剖面的解析解,m;w(xyz)——空间某一点的裂缝缝宽,m;$\bar{w}$——位置s的平均缝宽,m;${{\bar{w}}_{n}}$——第n组压裂施工参数条件下水力裂缝平均宽度,m;${{\bar{w}}_{\max }}$——水力裂缝最大平均宽度,m;z——裂缝中的某点位置,m;α——有效应力系数,无因次;α0——中间变量,Pa·sn$\nu $——泊松比,无因次;${{\nu }_{\text{d}}}$${{\nu }_{\text{s}}}$——动、静态泊松比,无因次;ρf,j——裂缝第j个元素中流体密度,kg/m3σHσhσv——水平最大主应力、最小主应力及上覆岩层压力,Pa;σi——第i层的原地应力,Pa;σmh——泥岩隔层的水平最小主应力,Pa;$\sigma \left( a \right)$——在位置a处的地层应力,Pa;σsh——砂岩储层的水平最小主应力,Pa;$\phi \left( {{n}'} \right)$——中间变量;ξIL——层间应力差异系数,无因次。
[1]
雷群, 杨战伟, 翁定为, 等. 超深裂缝性致密储集层提高缝控改造体积技术: 以库车山前碎屑岩储集层为例[J]. 石油勘探与开发, 2022, 49(5): 1012-1024.

DOI

LEI Qun, YANG Zhanwei, WENG Dingwei, et al. Techniques for improving fracture-controlled stimulated reservoir volume in ultra-deep fractured tight reservoirs: A case study of Kuqa piedmont clastic reservoirs, Tarim Basin, NW China[J]. Petroleum Exploration and Development, 2022, 49(5): 1012-1024.

[2]
TAN P, JIN Y, YUAN L, et al. Understanding hydraulic fracture propagation behavior in tight sandstone-coal interbedded formations: An experimental investigation[J]. Petroleum Science, 2019, 16(1): 148-160.

[3]
HOU B, CHANG Z, FU W, et al. Fracture initiation and propagation in a deep shale gas reservoir subject to an alternating-fluid-injection hydraulic-fracturing treatment[J]. SPE Journal, 2019, 24(4): 1-17.

[4]
TAN P, CHEN Z W, FU S H, et al. Experimental investigation on fracture growth for integrated hydraulic fracturing in multiple gas bearing formations[J]. Geoenergy Science and Engineering, 2023, 231(Part A): 212316.

[5]
楼一珊, 庄锦江, 黄荣樽. 岩石动、静弹性参数相关性的研究及其应用[J]. 江汉石油学院学报, 1989, 11(2): 62-69.

LOU Yishan, ZHUANG Jinjiang, HUANG Rongzun. Study on the interrelation of dynamic and static elastic parameters of rocks and its application[J]. Journal of Jianghan Petroleum Institute, 1989, 11(2): 62-69.

[6]
柳贡慧, 庞飞, 陈治喜. 水力压裂模拟实验中的相似准则[J]. 石油大学学报(自然科学版), 2000, 24(5): 45-48.

LIU Gonghui, PANG Fei, CHEN Zhixi. Development of scaling laws for hydraulic fracture simulation tests[J]. Journal of the University of Petroleum, China, 2000, 24(5): 45-48.

[7]
刘合, 黄有泉, 蔡萌, 等. 松辽盆地古龙页岩油储集层压裂改造工艺实践与发展建议[J]. 石油勘探与开发, 2023, 50(3): 603-612.

DOI

LIU He, HUANG Youquan, CAI Meng, et al. Practice and development suggestions of hydraulic fracturing technology in the Gulong shale oil reservoirs of Songliao Basin, NE China[J]. Petroleum Exploration and Development, 2023, 50(3): 603-612.

[8]
WENG X W, KRESSE O, CHUPRAKOV D, et al. Applying complex fracture model and integrated workflow in unconventional reservoirs[J]. Journal of Petroleum Science and Engineering, 2014, 124: 468-483.

[9]
KRESSE O, WENG X W, GU H R, et al. Numerical modeling of hydraulic fractures interaction in complex naturally fractured formations[J]. Rock Mechanics and Rock Engineering, 2013, 46(3): 555-568.

[10]
COHEN C E, KRESSE O, WENG X W. Stacked height model to improve fracture height growth prediction, and simulate interactions with multi-layer DFNs and ledges at weak zone interfaces[R]. SPE 184876-MS, 2017.

[11]
XU W Y, PRIOUL R, BERARD T, et al. Barriers to hydraulic fracture height growth: A new model for sliding interfaces[R]. SPE 194327-MS, 2019.

[12]
黄海, 郑永, 王毅, 等. 粗糙壁面压裂裂缝内支撑剂运移铺置特征[J]. 石油勘探与开发, 2024, 51(2): 399-408.

DOI

HUANG Hai, ZHENG Yong, WANG Yi, et al. Characteristics of proppant transport and placement within rough hydraulic fractures[J]. Petroleum Exploration and Development, 2024, 51(2): 399-408.

[13]
侯冰, 陈勉, 李志猛, 等. 页岩储集层水力裂缝网络扩展规模评价方法[J]. 石油勘探与开发, 2014, 41(6): 763-768.

HOU Bing, CHEN Mian, LI Zhimeng, et al. Propagation area evaluation of hydraulic fracture networks in shale gas reservoirs[J]. Petroleum Exploration and Development, 2014, 41(6): 763-768.

文章导航

/