李勇明，陈曦宇，赵金洲*，申峰，乔红军

水平井分段多簇压裂缝间干扰研究

李勇明1，陈曦宇1，赵金洲1*，申峰2，乔红军2

1.“油气藏地质及开发工程”国家重点实验室·西南石油大学，四川成都610500
2.陕西延长石油（集团）有限责任公司，陕西西安710075

水平井分段多簇压裂在现场得到了广泛运用，其压裂过程中普遍存在缝间干扰现象。缝间干扰有助于形成复杂裂缝网络以提高储层导流能力，但是也会导致起裂困难，甚至形成砂堵。因此有必要对分段多簇压裂的缝间干扰问题进行研究。对此，基于弹性力学建立了分析多簇裂缝诱导应力的数学模型，从起裂压力、裂缝宽度、簇间距等多方面研究了缝间干扰对水平井分段多簇压裂施工的影响。模拟结果显示，诱导应力会导致起裂压力升高、裂缝变窄，严重时将造成压裂施工失败。通过进行分析，给出了起裂过程及延伸过程中缝间干扰的影响关系。分析认为，利用缝间干扰提高改造体积时应当控制簇间距防止对压裂施工造成负面影响。研究结论对优化水平井分段多簇压裂设计具有指导意义。

缝间干扰；诱导应力场；起裂压力；裂缝宽度；多簇压裂

引言

随着近年来水平井钻井技术及大型水力压裂技术的不断发展，低渗透非常规油气资源成为新一轮开发热点［1-4］。Mayerhofer M等指出［5］，压裂改造体积是影响低渗透储层增产效果至关重要的一项参数。对于低渗透非常规油气藏，为了追求更大的压裂改造体积，压裂段、簇数近年来明显增加，井网也在不断加密以获得更高产量。然而，这些措施同时也加剧了缝间干扰、井间干扰问题［6］。缝间干扰有助于形成复杂裂缝网络以提高储层导流能力。但是缝间干扰也可能造成裂缝起裂延伸困难，甚至形成砂堵，从而影响施工效果。随着水平井分段多簇压裂施工日益增多，国外学者不断研究其缝间干扰问题。Cipolla C等［7］进行理论研究指出，当压裂主裂缝导流能力达到一定程度后，继续增加射孔簇数或者减小簇间距对最终增产效果影响甚微，因此应当适当控制射孔簇数。Olson J［8］采用位移不连续法（DDM）分析多簇裂缝应力场，研究缝间干扰以优化簇间距。Ajani A［9］等统计上百口水平压裂井进行分析，系统研究了页岩压裂中存在的干扰问题。本文通过建立分析多簇裂缝诱导应力的数学模型，研究了缝间干扰对水平井分段多簇压裂施工的影响，对优化分段多簇压裂设计具有指导意义。

1 起裂过程应力干扰分析

进行水平井分段多簇压裂时，压裂液注入使得射孔憋压起裂。由于地层岩石存在非均质性，压裂段内各簇裂缝很难同时起裂。部分射孔簇起裂较晚，会受到先破裂延伸裂缝的诱导应力干扰，导致破裂压力升高而加大破裂难度（图1）。

图1 水平井分段多簇压裂示意图Fig.1 The diagram of multi-stage horizontal fracturing

忽略多簇压裂造成的影响可能会低估破裂压力，而错误的预估可能造成压裂施工失败。因此，需要修正应力干扰下的破裂压力计算模型进行预测。假设已破裂裂缝面垂直，其缝高为2c，缝内净压力为pnet，处于均质弹性地层中，根据图2所示，建立垂直裂缝周围应力场模型。

图2 垂直裂缝诱导应力场模型Fig.2 The induced stress field of vertical fracture

假定压应力为正，拉应力为负。通过傅里叶变换可以求解平面应变假设下垂直裂缝周围应力场分布解。但由于破裂初期缝长较短，考虑平面应变假设误差较大。因此借鉴修正因子方法进行修正，近似获得有限缝长的诱导应力

其中：

为了计算已破裂延伸裂缝的诱导应力，需要先计算其裂缝半长l。在此采用二维延伸模型计算其裂缝长度。裂缝延伸半长计算公式为［10］

为求得时间t下裂缝i内的压裂液流量qi，需要求解裂缝缝口净压力

根据物质平衡，总流量等于流入每簇的压裂液流量之和

基于式（2）～式（4）采用牛顿-拉夫逊迭代算法进行求解，获得起裂后延伸t时刻的裂缝半长l并代入式（1）。将已先破裂延伸的裂缝诱导应力场作用到未起裂处原应力场上，获得尚未起裂点处局部地应力场

修正后的局部应力分布式可计算诱导应力干扰下的破裂压力，分析应力干扰对起裂过程的影响。

2 延伸过程应力干扰分析

裂缝延伸过程中，充满压裂液的水力裂缝在周围局部区域产生压应力，附近其他裂缝延伸受到抑制并且宽度下降。这个现象被称为应力阴影，已被实验研究所证实［10-11］。裂缝宽度是水力压裂施工最重要的参数之一。只有形成足够的裂缝宽度，才能在施工过程中确保足够的支撑剂进入到储层深处，同时避免发生砂堵事故。由于近年来压裂段、簇数提高，忽视应力阴影而对每条裂缝单独进行压裂设计容易高估部分裂缝宽度，导致优化设计失效。

实施水平井分段多簇压裂时，裂缝形态一方面受到储层岩石性质、地应力分布等地质因素影响，另一方面受到射孔簇设置，压裂液性质及排量等工程因素影响。因此，模拟分析裂缝延伸过程，需要同时解决裂缝延伸过程中的弹性力学问题以及压裂液在缝内流动的流体润滑问题。

对于多簇裂缝延伸中的力学问题，采用解析式求解变得过于困难。选用边界元法中的位移不连续法［12-13］，将裂缝离散成单元，数值求解多裂纹力学问题。把边界条件（裂缝面所受流体压力）代入控制方程，求出边界上的位移不连续量（缝宽及裂缝面滑移量），最后获得整个压裂区域的应力场与位移场。这种方法不但具有边界元的降维优势，同时能简单通过增加边界单元模拟裂纹扩展，在水力裂缝延伸模拟中具有极大优势。

裂缝周围区域的应力场可以通过位移不连续量Dx、Dy表示为

函数f（x，y）为

位移不连续法基于二维平面应变假设，适合于缝高远远大于缝长的情况。然而，实际压裂裂缝缝高有限且通常小于缝长。因此采用式（6）求解误差较大。Olosn J E［14］提出三维因子用于修正平面应变假设带来的误差。三维因子Gf可由式（8）得到

因而，修正后的式（6）变为

由于受到诱导应力干扰，裂缝延伸时路径会发生偏转。复合裂纹断裂理论主要用于研究裂纹的延伸方向及延伸条件。其中应用最广泛的方法是最大正应力理论及能量密度因子理论。本文选用最大正应力理论计算裂缝的延伸方向及条件，其理由在于大量实验证明最大正应力理论在脆性材料上判断裂缝延伸方向上准确，其理论的假设要点是：

（1）裂纹扩展方向：周向正应力最大的方向为裂纹延伸方向。

（2）裂缝扩展条件：应力强度因子超过I型临界值。

根据最大正应力理论，要使得周向正应力最大，水力裂缝扩展方向θc应当满足

裂缝尖端I型、II型应力强度因子KI、KII可以根据线弹性力学用尖端单元的Dx、Dy计算获取

由于常单元位移不连续法中假设一个单元内的位移不连续量即缝宽为常数，导致数值求解的缝宽略大于解析解值，尤其在裂缝宽度下降梯度较大的尖端区域误差较大。因此为降低应力强度因子计算误差，式（12）中经验常数C取值0.806。

假设裂缝延伸过程中，压裂液完全充满裂缝。实施分段多簇压裂作业时，压裂液泵入压裂段中，段内各簇裂缝流入压裂液流量之和等于压裂液总排量

假定压裂液不可压缩，其在裂缝内的流动满足泊肃叶公式

当压裂液表现为幂律流体特性时，黏度µ可以由式（15）计算

压裂过程中压裂液满足物质平衡式

式（8）～（9）组成延伸过程的力学方程组，式（13）～（16）组成延伸过程的压裂液流动方程组。针对这个非线性方程组问题，通过假设一组初始缝宽和净压力，迭代求解，最终获得收敛的一组数值解。

3 实例计算分析

采用某油田的低渗透油藏区块，计算分析缝间干扰对水平井分段多簇压裂所造成的影响，油藏区域基本参数见表1。

表1 油藏基本参数Tab.1 The parameter of reservoir

图3显示了水平井中单条破裂裂缝的诱导应力。

图3裂缝无因次诱导应力分析Fig.3 The dimensionless analysis of fracture induced stress

图3 中横坐标为远离裂缝的距离与裂缝高度之比，纵坐标为裂缝产生的诱导应力与裂缝本身净压力之比。蓝色、绿色、红色曲线分别对应x、y、z方向的诱导应力。裂缝诱导应力随距离增大其迅速减小。距离裂缝0.5倍缝高以内位置（d＜c）受到诱导应力作用明显，而远离裂缝达裂缝缝高2倍以上时（d＞4c），其诱导应力几乎为零。诱导应力作用范围受到裂缝高度所控制。

裂缝延伸早期，由于裂缝缝长较短，此时缝长对诱导应力作用影响同样显著。图3中同时对比了不同半缝长情况下的诱导应力。如图中箭头方向所示，随裂缝半长l从c/3增长到8c，诱导应力逐渐增大。这说明当某簇裂缝先行破裂延伸后，随着裂缝缝长增加，其他未破裂簇所受诱导应力干扰将会加剧。但是，随着裂缝长度的增长，诱导应力增长趋势减缓，最终趋于平面应变假设下的情况。

假设邻簇裂缝已经破裂延伸（破裂时间设定为0），尚未起裂处所需破裂压力会随时间变化，针对诱导应力对起裂压力的影响，采用最大张应力破裂准则对区块数据（表1）进行计算分析（图4）。分析显示，随着时间增加，所需破裂压力会短时间陡升，之后缓慢上升。破裂压力的这种上升趋势由两方面因素所影响：一方面随着已破裂缝缝长快速增加，诱导应力会迅速上升并趋于稳定；另一方面已破裂裂缝净压力也缓慢上升，从而缓慢增大应力干扰。同样，图4显示簇间距对裂缝起裂压力影响也很明显。较低的簇间距会使得破裂压力显著提升；而随着簇间距增大，起裂受到的干扰逐渐下降，所需破裂压力趋近于单簇压裂时的情况。

图4 簇间距对破裂压力影响Fig.4 The effect of fracture space on initiation pressure

综上所述，在低簇间距下未能快速起裂的射孔簇，随着时间增加可能会更加难以起裂，最终导致压裂簇失效。这种簇间距较小时部分射孔簇未能压开的现象，也在一些现场施工也得到了证实。因此，在通过设置较小簇间距增大改造效果时，也必须综合考虑其带来的负面影响。

同样采用表1数据，以目标区块段内压3簇为例，限制缝高为储层高度，采用建立的模型模拟裂缝延伸时簇间距对缝间干扰的影响，结果见图5。随着簇间距减小，一方面，应力干扰使得裂缝倾向于扭曲延伸；另一方面，裂缝延伸时受到邻缝压应力影响而变得更为窄小。

图5 应力干扰对裂缝宽度的影响Fig.5 The impact of stress shadow on widths

为了定量分析多簇裂缝缝间干扰对裂缝宽度的影响。可以粗略采用式（17）计算裂缝所受干扰［15-16］，分析结果见图6。

图6 不同簇间距下裂缝宽度无因次分析Fig.6 The dimensionless analysis of fracture width

图6中横坐标为簇间距与裂缝高度之比，纵坐标为裂缝平均宽度与其未受应力干扰的宽度之比。当簇间距减少到缝高的0.5倍时，裂缝宽度减少明显。随着段内簇数N增加，缝间干扰影响增大，缝宽也在逐渐下降，但下降的趋势逐渐减缓。当簇间距大于两倍裂缝高度时，诱导应力对裂缝宽度的影响较小，此时可以根据油气分布实际情况进一步减少簇间距或者增加簇数以提高储层改造体积。

图7显示不同簇位置下应力干扰对缝宽的影响。裂缝随着簇间距缩小而宽度变窄时，外边缘簇裂缝明显宽于中部簇裂缝。显然，这是因为中部裂缝由于位置原因受到更强的诱导应力所致。随压裂时间推进，中部裂缝所受抑制将会更加明显。因此进行多簇压裂施工时，应当重点保证中部簇有效破裂及延伸，防止砂堵。

压裂液黏度同样显著影响裂缝形态，图8模拟了裂缝延伸时压裂液黏度对裂缝尺寸以及缝间干扰的影响。液体黏度减小使得流体在缝内所受阻力下降，高黏度压裂液能够形成宽裂缝，相比之下低黏度流体使得流体波及更远范围，形成长而窄的裂缝。此外，压裂液黏度增大虽然使得裂缝净压力增大，导致诱导压应力增大，但对裂缝的抑制程度却影响不大。影响裂缝抑制程度的主控因素是裂缝高度及簇间距。

图7 不同位置裂缝宽度变化Fig.7 The change of widths with fracture positions

图8 压裂液黏度对裂缝尺寸的影响Fig.8 The impact of viscosity on fractures

4 结语

（1）进行多簇压裂时，由于段内各簇裂缝难以同时起裂，先破裂裂缝会造成邻近起裂簇所需起裂压力上升。在低簇间距下未能快速起裂的射孔簇，将会随着时间增加更加难以起裂，影响压裂施工。对此，必须充分重视应力干扰对起裂压力所造成的影响。

（2）随着簇间距缩小，裂缝受到应力干扰更为严重。一方面，裂缝倾向于扭曲延伸；另一方面，裂缝受到邻缝压应力而变窄小。随着簇间距缩小，其裂缝缝宽不断减小；段内簇数增加，缝间干扰影响也逐渐增大，同样导致缝宽下降。施工时，应当控制簇间距避免缝宽降低过多造成裂缝延伸困难甚至砂堵事故。

（3）进行分段多簇压裂施工时，中部簇缝宽严重受到压抑。施工设计时应当加大排量保证中部簇有效破裂延伸，防止砂堵。

（4）压裂液黏度显著影响裂缝形态，但不会明显影响应力阴影现象对裂缝的压抑程度。影响应力阴影较大的因素为裂缝高度及簇间距。

（5）水平井分段多簇压裂近年来使用广泛，应力阴影的影响非常明显。对此，不能忽视裂缝被压抑的负面效果，而需要有效控制缝间干扰影响程度，提高改造体积。

符号说明

c-半缝高，m；

［1］吴奇，胥云，王腾飞，等.增产改造理念的重大变革体积改造技术概论［J］.天然气工业，2011，31（4）：7-12. WU Qi，XU Yun，WANG Tengfei，et al.The revolution of reservoir stimulation：An introduction of volume fracturing［J］.Natural Gas Industry，2011，31（4）：7-12.

［2］唐颖，张金川，张琴，等.页岩气井水力压裂技术及其应用分析［J］.天然气工业，2010，30（10）：33-38. TANG Ying，ZHANG Jinchuan，ZHANG Qin，et al.An analysis of hydraulic fracturing technology in shale gas wells and i ts application［J］.Natural Gas Industry，2010，30（10）：33-38.

［3］李勇明，彭瑀，王中泽.页岩气压裂增产机理与施工技术分析［J］.西南石油大学学报（自然科学版），2013，35（2）：90-96. LI Yongming，PENG Yu，Wang Zhongze，et al.Analysis of shale gas fracture stimulation mechanism and operating techniques［J］.Journal of Southwest Petroleum University（Science&TechnologyEdition），2013，35（2）：90-96.

［4］KING G E.Thirty years of gas shale fracturing：What have we learned？［C］//SPE Annual Technical Conference and Exhibition，Society of Petroleum Engineers，2010.

［5］MAYERHOFER M，LOLON E，WARPINSKI N，et al. What is stimulated rock volume？［C］.SPE 119890，2008.

［6］LAWAL H，JACKSON G，ABOLO N，et al.A novel approach to modeling and forecasting frac hits in shale gas wells［C］.SPE 164898，2013.

［7］CIPOLLAC，LOLONE，ERDLEJ，etal.Reservoirmodeling in shale-gas reservoirs［J］.SPE Reservoir Evaluation&Engineering.2010，13（4）：638-653.

［8］OLSON J，DAHI TALEGHANI A.Modeling simultaneous growth of multiple hydraulic fractures and their interaction with natural fractures［C］.SPE 119739，2009.

［9］AJANI A，KELKAR M.Interference study in shale plays［C］.SPE 151045，2012.

［10］VALK P，ECONOMIDES M J.Hydraulic fracture mechanics［M］.Wiley，New York，1995.

［11］GERMANOVICH L N，ASTAKHOV D K.Fracture closure in extension and mechanical interaction of parallel joints［J］.Journal of Geophysical Research（Solid Earth），2004，109（B2）：B02208-1.

［12］CROUCH S L.Solution of plane elasticity problems by thedisplacementdiscontinuitymethod.I.Infinitebodysolution［J］.International Journal for Numerical Methods in Engineering，1976，10（2）：301-343.

［13］CROUCH S L，STARFIELD A M.Boundary element methods in solid mechanics：With applications in rock mechanics and geological engineering：Winchester，Massachusetts［M］.Allen and Unwin Ltc，1983.

［14］OLSON J E.Predicting fracture swarms—The influence of subcritical crack growth and the crack-tip process zone on joint spacing in rock［J］.Geological Society，2004，231（1）：73-88.

［15］MEYER B，BAZAN L.A discrete fracture networkmodel for hydraulically induced fractures-theory［C］.SPE 140514，2011.

［16］DAHI-TALEGHANI A，OLSON J E.Numerical modeling of multistranded-hydraulic-fracture propagation：Accounting for the interaction between induced and natural fractures［J］.SPE Journal，2011，16（3）：575-581.

乔红军，1970年生，男，汉族，河南新乡人，高级工程师，主要从事油气田开发方面工作。E-mail：hongjunqiao2006@126.com

编辑：王旭东

编辑部网址：http：//zk.swpuxb.com

The Effects of Crack Interaction in Multi-stage Horizontal Fracturing

LI Yongming1，CHEN Xiyu1，ZHAO Jinzhou1*，SHEN Feng2，QIAO Hongjun2
1.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation，Southwest Petroleum University，Chengdu，Sichuan 610500，China
2.Shaanxi Yanchang Petroleum（Group）Co.Ltd.，Xian，Shaanxi 710075，China

Multi-stage horizontal fracturing has been widely used in the field with the inevitable crack interaction.Crack interaction is conducive to initiation and extension of fractures，which forms complex fracture network and improves reservoir conductivity.However，crack interaction also has adverse effects on hydraulic fracturing resulting in poor stimulation.Based on the theory of elasticity，mathematical models are established to study the effects of crack interaction in multiple cluster fractures，from the aspects of initiation pressure，fracture width and the position of fracture.The study shows that it is hard to get successful stimulation with serious crack interaction because the phenomenon such as change of fracture width resulting from crack interaction.To get the relationship of cluster space，fracture height and crack interaction，a sensitivity analysis is conducted.The result indicates that the cluster space should be reasonable to avoid the adverse effects of crack interaction.The study result contributes to the optimization of the stimulation of multiple stages horizontal well.

crack interaction；induced stress field；initiation pressure；fracture width；multiple cluster stimulation

李勇明，1974年生，男，汉族，四川营山人，教授，博士，主要从事油气藏压裂酸化工程技术方面的基础理论及应用技术研究。E-mail：swpifrac@163.com

陈曦宇，1989年生，男，汉族，重庆人，博士研究生，主要从事油气藏压裂酸化工程技术方面的基础理论及应用技术研究。E-mail：cxyswpu@gmail.com

赵金洲，1962年生，男，汉族，湖北仙桃人，教授，博士生导师，主要从事油气田开采和增产新技术新理论的研究工作。E-mail：zhaojz@swpu.edu.cn

申峰，1970年生，男，汉族，河南新乡人，高级工程师，主要从事油气田开发方面工作。E-mail：sf88899614@163.com

10.11885/j.issn.1674-5086.2013.11.02.01

1674-5086（2016）01-0076-08

TE357

A

http：//www.cnki.net/kcms/detail/51.1718.TE.20160108.0916.012.html

李勇明，陈曦宇，赵金洲，等．水平井分段多簇压裂缝间干扰研究［J］．西南石油大学学报（自然科学版），2016，38（1）：76-83．

LI Yongming，CHEN Xiyu，ZHAO Jinzhou，et al．The Effects of Crack Interaction in Multi-stage Horizontal Fracturing［J］．Journal of Southwest Petroleum University（Science&Technology Edition），2016，38（1）：76-83．*

2013-11-02网络出版时间：2016-01-08

赵金洲，E-mail：zhaojz@swpu.edu.cn

国家重点基础研究发展计划（973计划）（2013CB228004）；国家自然科学基金（51344005）；西南石油大学研究生创新基金（CXJJ2015043）。