正交节理与应力比值系数对水力裂缝扩展影响规律研究

2020-09-25张伯虎周昌满郑永香刘建军

油气藏评价与开发 2020年5期

张伯虎，周昌满，郑永香，刘建军,2

（1.西南石油大学地球科学与技术学院，四川成都610500；2.中国科学院武汉岩土力学研究所，湖北武汉430071）

水力压裂是油田增产和致密油气开发的关键技术，水力裂缝可以有效提高储层的渗透能力[1-3]。水力裂缝扩展机理是预测压裂后最终裂缝网络形态的基础[4]。前人针对水力裂缝扩展规律开展了一系列实验研究和数值模拟[5-6]。杨文波等[7]借鉴北美致密储层经验，在东胜气田开展了体积压裂实验。郑永香等[8]基于断裂力学分析了裂缝内注水对周围地应力场的干扰。GUO 等[9]研究表明泵压曲线波动是由于压裂液严重漏失至天然裂缝或断层。TAN 等[10]研究多因素水平层理页岩垂直裂缝的萌生和破裂，并提出5 种起裂和扩展模式。刘海龙等[11]研究了定向射孔水力压裂的起裂压力。XU等[12]和SUN等[13]研究层理对水压裂缝的影响，表明水压裂缝由层理和地应力共同作用。WANG 等[14]考虑页岩各向异性的数值模型，认为裂缝扩展是由原岩应力和各向异性特性综合决定。侯冰等[15]采用位移不连续方程，研究了水力裂缝在随机分布的天然裂缝页岩地层影响下形成的复杂缝网形态。ZHANG等[16]采用二维边界元法对水力裂缝在层状沉积岩中层理面处偏转进行数值研究，发现低黏度的流体产生的水力裂缝在和层理面的相交处会受到抑制，转而流入层理面。

国内外学者研究表明，水力裂缝的扩展容易受到地应力状态、节理、断层以及随机天然裂缝的影响，而基于正交节理（水平节理和垂直节理正交）和应力状态对水力压裂裂缝扩展行为的研究却较少。本文采用块体离散元，对三维正交裂缝的扩展规律进行模拟，分析了不同应力条件对裂缝缝内压力和裂缝宽度的影响，初步讨论了产生缝高限制的原因，得出剪应力是节理相交处产生阻碍作用的主要因素，并对裂缝边缘应力状态进行了分析，借助真三轴水力压裂实验，验证了控制应力比值系数可以使水力裂缝穿过水平节理的结论。本文提出的应力比值系数与正交节理等对水力裂缝扩展时是否穿层与转向、是否形成复杂缝网等储层改造效果有较大的影响。

1 数值模型与方案设计

1.1 块体离散元技术

块体离散元技术是一种可以用于三维不连续建模的先进工具，其处理非线性变形和破坏都集中在节理面处的问题非常有优势，其模拟不受静态或者动态荷载下的影响，可以真实地反映不连续介质的力学响应。利用岩体中的不连续面把岩体分割成许多完整的小岩块，把这些岩块视为刚体或可变形体，根据相关力学法则来计算节理岩体中对应的节理面的剪应力和正应力，并自动判别各块体之间的接触点，将其作为应力边界条件来模拟岩体的分离、转动以及大变形等情况[17-18]。因此，可以用来模拟水力压裂的液体注入与裂缝扩展过程。近年来已经有很多学者采用离散元模拟水力压裂：MCLENNAN 等[19]采用块体离散元软件3DEC 模拟了正交裂缝网络条件下裂缝的扩展规律；HAMIDI 等[20]采用3DEC 研究了不同流体性质、流体注入速率、地应力以及岩石力学性质对水力裂缝扩展的影响；NAGEL等[21]采用3DEC模拟了不同工况影响下天然裂缝剪切和张开的破坏规律；ZHENG 等[22]采用块体离散元对水力裂缝在含有胶结裂缝的地层进行了数值模拟。尽管模拟水力压裂的注入过程时，初始水力裂缝和水力裂缝面需要预置[23]，且水力裂缝只能按照预置的裂缝面进行扩展，但与其他有着许多假设条件的商业化软件相比，其仍然是一个比较好的模拟裂隙流体方法。

1.2 模型的建立

采用地层面的弱化来模拟节理，从而建立含正交节理的数值模型，模型尺寸为20 m×20 m×10 m，如图1a所示。考虑2个水平层理面和1个竖向节理面，并进行相交，将注入点设在竖向节理面的中心，如图1b所示。应力边界条件如图1c所示。

图1 模型示意图Fig.1 Model schematic

地应力大小与埋深有关，具体数值如式（1）所示。

式中：σx，σy，σz为在x，y，z方向地应力，MPa；ρ为岩体密度，取2 600 kg/m3；g为重力加速度，取9.81 m/s2；h为注入点处的地层深度，取3 000 m[24]；k为应力比值系数，为天然水平地应力与竖向地应力的比值。

为了分析水力裂缝在节理面内的扩展情况，通过改变地应力状态来研究遇到相交节理时水力裂缝的扩展行为。在竖向节理面上布置注入点和监测点，如图2所示，点1为注入点，点2～9为监测点。监测点处的压力值突然上升，说明压裂液到达该点，该点水力裂缝开始扩展。通过测定任意两点间压裂液到达的时间差，可以确定裂缝的扩展速度。其他参数取值如表1所示，相关属性借鉴川南龙马溪组页岩参数[25-26]。

图2 监测点Fig.2 Monitoring points

表1 模型参数Table1 Model parameter

1.3 方案设计

为了探究地应力状态对水力裂缝遭遇正交节理的影响规律，设置3 组应力状态，应力比值系数k分别取0.8、1.0、1.2，注入速度均取0.03 m3/s，形成3个计算方案。

2 结果分析

2.1 阻碍作用与三维裂缝延伸方向分析

从数值模拟结果发现当竖向水力裂缝扩展到正交节理相交处时，裂缝会受到阻碍作用而扩展缓慢，并可能沿着缝宽方向延伸。但当裂缝扩展达到一定位置时，裂缝会突破阻碍而继续扩展，其延伸方向受k值影响。当k为0.8 和1.0 时，竖向水力裂缝会贯穿水平节理继续沿竖向延伸；当k为1.2时，裂缝转向水平节理扩展。

如图3所示，当k取0.8 和1.0 时，水力裂缝直接穿过水平节理；当k取1.2时（为了更好展示裂缝的转向，将上部水平节理面隐藏，其结果与下部水平节理面对称），水力裂缝未能穿过水平节理，而是转向水平节理面扩展。根据文献[27-30]研究结果，当水力裂缝与节理面的夹角逼近90°时，水力裂缝将直接穿过，并不会发生转向。相关结论假设原始水力裂缝垂直于最小主应力而得出，但在实际地层中存在不垂直于最小主应力的弱结构面，且该弱结构面强度较低，水力裂缝将会优先在弱结构面内扩展，因此，水力裂缝与天然弱面夹角接近90°时仍可能转向。综上，在判断不同逼近角下的裂缝扩展状态时，应综合考虑裂缝面所处的应力状态，单独以逼近角判断裂缝形态是不准确的。

图3 不同应力场状态下的裂缝扩展形态Fig.3 Crack propagation morphology under different stress field conditions

与二维状态下裂缝只在平面空间内延伸不同，三维状态下，当水力裂缝到达水平节理与垂直节理相交处时，裂缝并没有直接穿过或转向，而是沿着缝宽方向扩展，说明相交节理会对裂缝扩展造成阻碍作用，导致裂缝沿着缝宽方向扩展。这与实际压裂过程中形成裂缝的长度远大于其高度有关。当裂缝沿缝宽扩展达到一定程度后会突破其阻碍作用，继续竖向延伸或转向水平扩展。因此，三维条件下的节理面内裂缝扩展与二维裂缝扩展有较大的区别，某一节理面内裂缝的扩展维度将从一维线扩展转为二维面扩展。

2.2 裂缝宽度演化与分布规律

在不同应力条件下，当水力裂缝克服节理相交处的阻碍作用后，水力裂缝将表现为直接穿过或者转向。分别选取k=1.0 和k=1.2 两种应力情况，并选择4 个扩展阶段分析裂缝最大宽度的变化，如图4a所示。在压裂初期，水力裂缝受到阻碍作用，裂缝宽度均逐渐变大。但当水力裂缝克服阻碍作用后，相交面阻力消失，裂缝将穿过相交面，裂缝宽度随应力比值系数变化而有所变化。当k=1.0时，裂缝的最大宽度逐渐趋于稳定，并在后期有所减小，但由于裂缝面上的正应力未发生变化，所以裂缝宽度基本保持不变。当k=1.2 时，裂缝突破相交面的阻碍后，裂缝转向进入水平节理面，此时裂缝最大宽度下降，这是因为水平节理面上的正应力低，裂缝扩展压力也较竖向节理面低，因此缝内压力会降低，导致水力裂缝的水平宽度降低。

2.3 裂缝压力演化与分布规律

由于受到阻碍作用，水力裂缝中的缝内压力会产生一定幅度的波动，但水力裂缝的尖端远离起裂点后的缝内压力容易受到k值影响。

图4b为不同k值时注入压力变化规律。在注入初期，注入压力快速上升，当达到破裂压力时裂缝开始扩展，注入压力基本保持稳定。在扩展初期注入压力具有较大波动，是因为裂缝扩展使裂缝空间突然增大，流体压力降低；随着压裂液的补充，压力迅速得到恢复，后期基本平衡。当k=1.2时，裂缝在扩展后期会发生转向，随后注入压力降低，其值与k=1.0的结果趋于一致。这说明注入压力会受到裂缝扩展的影响，可以用来表征裂缝当前的扩展状态。因此，可以根据注入压力大小来判定裂缝的扩展难易程度。

图4c、图4d 反映不同监测点的缝内压力变化规律。当监测点处压力上升并与注入点相同时，说明该处开裂，开裂后缝内压力不再变化。当k=1.0时，在起裂初期，由于附近裂缝的开启，注入压力会有一定幅度的波动，当裂缝尖端远离起裂点后，缝内压力保持稳定。当k=1.2时，裂缝先沿竖向裂缝发展，此时缝内压力保持稳定；但当水力裂缝遇到水平层理后，裂缝转而沿着水平方向发展，缝内压力产生突降，随后压力会基本保持不变，使水平向裂缝持续发展。

3 相交裂缝的阻碍作用机制

上述结果表明，节理面相交处对水力裂缝的扩展具有阻碍作用。前人对此也做过解释，起初认为节理面上下岩层间的弹性模量和应力差异是影响缝高的关键因素，但这只适用于不同性质岩层内的裂缝扩展。后来发现弹性模量和泊松比并不是缝高限制的必要因素，在岩性一致的岩层内也会出现缝高限制现象。而影响缝高扩展的主要因素在于节理面上的剪应力，其产生的阻碍作用限制了裂缝的扩展。图5为裂缝发生转向前后水平与竖向节理面上的剪应力分布规律。由图5可知，裂缝在未突破阻碍作用前，水平节理面上有以相交线为对称轴的剪应力；当裂缝穿过后剪应力值减小，并且由于水力裂缝沿水平节理面扩展，故在竖向节理面上也产生了剪应力。这说明剪应力是导致节理相交处产生阻碍作用的主要因素。相对而言，由于节理上的剪应力较大，导致水平节理面和竖向节理面滑动；水力裂缝扩展到此处时会产生不连续滑动，进而使水力裂缝尖端钝化，造成裂纹尖端的应力释放传递到节理表面阻碍水力裂缝继续扩展。只有当缝内压力足够大的时候，裂纹尖端释放的应力才会使节理破坏，进而使水力裂缝突破阻碍作用而向前扩展。

图4 地层应力状态对裂缝扩展的影响规律Fig.4 Effect of formation stress state on crack propagation

图5 裂缝发生转向前后的节理面剪应力分布Fig.5 Shear stress distribution of joint surface before and after crack passing through

4 水力裂缝克服阻碍作用的机制

在均质理想模型下，竖向节理面内的水力裂缝扩展受到阻碍作用后会沿着水平方向扩展，其扩展形状如图6a所示。当裂缝的长度和宽度达到一定比例后，水力裂缝克服阻碍作用向外扩展。探索在同等的缝内压力条件下，水力裂缝最初与交界处相遇时不能克服阻力，但最终又能克服阻碍作用的原因。基于此，建立如图6a的理想模型，在模型上加载水平应力σx、σy和竖向应力σz，计算裂缝边缘的法向应力Sz分布规律。图中a为裂缝长度，b为裂缝宽度，分别取a=2 m，b=1 m，缝内液体压力p取82 MPa（k=1，注入速度为0.01 m3/s），取σx=σy=σz=75 MPa。裂缝边缘的法向应力Sz如图6b所示。在同等缝内压力条件下，长边的法向应力大于短边的法向应力，故水力裂缝会突破阻碍作用继续扩展。

为分析裂缝边缘法向应力受缝内压力p和竖向应力σz的影响，设置5组缝内压力值和7组竖向应力值，其他参数不变。分析方案如表2所示。计算结果如图6c和图6d所示。

表2 裂缝边缘应力分析计算参数Table2 Crack edge stress analysis programMPa

在竖向应力σz一定而缝内压力p大于75 MPa时，p增大则长短边之间的法向应力差值增大。同理，在p一定而σz小于82 MPa 时，随着σz增加，长边与短边的差值减小。这说明p大于σz情况下，p与σz的差值越大，长短边的应力差值越大。综上，裂缝的起裂和扩展压力与裂缝的几何形态有关，单独通过注入压力来定义起裂和扩展的相关压力不够准确。

5 实验验证及工程意义

以真三轴水力压裂实验来验证通过控制地应力状态实现水压裂缝穿透层理的目标。实验中压裂液添加绿色颜料便于判断压裂液走向。相关参数见表3，其中，σv表示竖向主应力，σH表示水平最大主应力，σh表示水平最小主应力。

表3 真三轴水力压裂物理模拟实验参数Table3 Physics simulation experimental parameters of true triaxial hydraulic fracturing

压裂前各面天然裂缝情况如图7a 所示，岩样存在水平节理。图7b为注入绿色压裂液形成的主裂缝面展开图，发现两侧主裂缝面均有绿色。图7c 为右侧主裂缝面局部展示图，绿色压裂液直接穿过层理处裂缝。图7d 为打开层理处裂缝展开图，可以发现水力裂缝没有沿着层理面扩展，而是直接穿过层理面。综上，可以发现水力裂缝沿着垂直于最小水平主应力方向扩展，受竖向主应力影响较大，水力裂缝扩展至层理时，直接贯穿层理，而非沿着层理面扩展，进而实现了通过控制地应力状态实现水压裂缝穿透层理的目标。

图7 压裂前后试样Fig.7 Samples before and after fracturing

水力裂缝的缝高控制对于实际压裂工程至关重要。一般在页岩气开采过程中，必须使裂缝控制在产气区域，使裂缝尖端扩展到地层弱面时停止延伸，以控制裂缝高度。因为延伸的裂缝若突破地层弱面会导致含水层流入，严重影响天然气的开采。当然，穿层压裂也是实际压裂工程的常用手段，比如夹层砂岩和泥岩中的水力压裂，裂缝必须穿过泥岩地层，垂直连接孤立的含气层，通过这些裂缝进入井中，达到采集油气的目的[31-33]。因此，裂缝扩展受应力比值和阻碍等影响的研究成果无疑会给这些压裂目标和手段提供理论指导。