天然气地球科学, 2021, 32(1): 109-118 doi: 10.11764/j.issn.1672-1926.2020.11.006

天然气开发

页岩气储层水力压裂扩展有限元模拟方法及应用

张瑛堃,1, 陈尚斌,1,2, 李学元1, 王慧军1

1.中国矿业大学资源与地球科学学院,江苏 徐州 221008

2.煤层气资源与成藏过程教育部重点实验室,江苏 徐州 221008

Hydraulic fracturing simulation technology of shale gas reservoir and application of extended finite element method

ZHANG Ying-kun,1, CHEN Shang-bin,1,2, LI Xue-yuan1, WANG Hui-jun1

1.School of Resources and Geosciences,China University of Mining and Technology,Xuzhou 221008,China

2.Key Laboratory of the Ministry of Education on Coalbed Methane Resources and Accumulation Process,Xuzhou 221008,China

通讯作者: 陈尚斌(1983-),男,甘肃通渭人,教授,博士生导师,主要从事页岩气地质研究.E-mail: shangbinchen@163.com.

收稿日期: 2020-08-20   修回日期: 2020-10-29   网络出版日期: 2020-12-31

基金资助: 国家自然科学基金.  41772141
江苏省自然科学基金.  BK20181362

Received: 2020-08-20   Revised: 2020-10-29   Online: 2020-12-31

作者简介 About authors

张瑛堃(1997-),女,河北唐山人,硕士研究生,主要从事页岩气地质与开发研究.E-mail:m15512011623@163.com. , E-mail:m15512011623@163.com

摘要

页岩气储层低孔低渗,需用水力压裂等方法进行储层改造方可获得经济产能。储层改造中裂缝的形态和分布对体积改造效果至关重要。为了研究压裂裂缝的模拟方法,系统调研和对比了储层水力压裂模拟常用方法,开展了扩展有限元模拟,研究表明:①水力压裂物理模拟实验能够直观观测裂缝的形态及展布特征,但因试样尺寸等问题难以代表实际储层压裂情形;②常用的数值模拟方法有边界元法、非常规裂缝模型、离散化缝网模型和扩展有限元法等,这些方法各有优缺点,需做有针对性的改进才能更好地模拟真实页岩储层压裂情况;③应用扩展有限元法模拟水力压裂和分段顺序压裂过程中裂缝的延伸情况,得到射孔方向与最大水平主应力之间夹角和诱导应力对压裂裂缝的影响,夹角越大,裂缝偏转角度越小,偏转距离越大,初始破裂压力越高,裂缝稳定延伸的压力也越大,而诱导应力的存在会抑制压裂裂缝的延伸。对实际压裂工程中射孔方向的选择和分段压裂射孔间距的设计具有指导意义。

关键词: 页岩储层 ; 水力压裂 ; 物理实验 ; 数值模拟 ; 扩展有限元

Abstract

Low porosity and low permeability of shale gas reservoirs require hydraulic fracturing and other methods to achieve economic productivity. The shape and distribution of fractures are very important to the volume transformation. In order to study the simulation methods of hydraulic fracture, the common methods of reservoir hydraulic fracture simulation are systematically investigated and compared, and the extended finite element simulation is carried out. The results show that: (1) The physical experiment of hydraulic fracturing can visually observe the fracture morphology and distribution characteristics, but it is difficult to represent the actual fracturing situation of the reservoir due to the sample size and other problems. (2) The commonly used numerical simulation methods include boundary element method (BEM), unconventional fracture model (UFM), discrete fracture network (DFN) and extended finite element method (XFEM). These methods have their own advantages and disadvantages, which need to be improved to better simulate the real shale reservoir fracturing. (3) The extended finite element method is used to simulate the fracture extension of hydraulic fracturing and staged sequential fracturing. The influence of the angle between the perforation direction and the direction of the maximum horizontal principal stress and the induced stress on the fracturing pressure is obtained. The larger the angle is, the smaller the fracture deflection angle is, while the larger the fracture deflection distance is. The larger the angle is, the higher the initial fracturing pressure is, and the greater the pressure of stable fracture extension is. The induced stress will hinder the fracture extension. The simulation results have guiding significance for the selection of perforation direction and the design of perforation spacing in staged fracturing in actual engineering.

Keywords: Shale reservoir ; Hydraulic fracturing ; Physical experiment ; Numerical simulation ; Extended finite element

PDF (2068KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张瑛堃, 陈尚斌, 李学元, 王慧军. 页岩气储层水力压裂扩展有限元模拟方法及应用. 天然气地球科学[J], 2021, 32(1): 109-118 doi:10.11764/j.issn.1672-1926.2020.11.006

ZHANG Ying-kun, CHEN Shang-bin, LI Xue-yuan, WANG Hui-jun. Hydraulic fracturing simulation technology of shale gas reservoir and application of extended finite element method. Natural Gas Geoscience[J], 2021, 32(1): 109-118 doi:10.11764/j.issn.1672-1926.2020.11.006

0 引言

近年来,我国页岩气勘探开发取得重要进展,已实现初步商业化开发1-2。尽管我国页岩气资源潜力巨大,但页岩气储层渗透率极低,要实现长期稳产,需依赖于理想的储层改造技术。目前应用最广泛的储层改造方法为水平井分段水力压裂技术3-4。水力压裂的效果取决于裂缝的形态及延伸展布特征5。我国页岩储层埋藏深,物性及所受地应力复杂,直接监测体积改造效果困难6。物理实验和数值模拟等方法可以作为直接监测的有效替代方法,在储层压裂研究中被广泛应用7-12。这些方法可以观察分析压裂裂缝的起裂压力、形态及其与天然裂缝的联系,评价水力压裂的效果,进行压裂设计方案优化。已有学者13-15研究了射孔方向与地应力之间的关系对压裂裂缝起裂的影响,但关于不同射孔方向下单裂缝和分段顺序压裂裂缝的扩展形态及规律,尚未有较为系统的研究。本文系统调研了水力压裂裂缝扩展模拟方法,对比分析了其优点和局限性,并采用扩展有限元法对压裂裂缝的扩展进行了数值模拟,与物理试验结果对比,模拟了不同射孔方向和分段顺序压裂下裂缝的延伸情况。

1 页岩储层水力压裂物理实验与模拟方法

1.1 物理实验方法

页岩储层水力压裂物理实验是指用天然露头或者人工浇铸试样,加载三向应力,用井筒注入液体,模拟水力压裂过程,直到裂缝延伸至试样边界16。国外页岩储层水力压裂物理模拟实验的文献较少17-19。国内学者设计了大尺寸真三轴模拟压裂试验装置,模拟水力压裂过程,分析压裂裂缝的延伸机理。侯振坤等6采用大尺寸真三轴试验系统,研究了水力压裂裂缝的形态、形成机制和水力裂缝与天然裂缝的相互作用方式(表1);解经宇等20采用真三轴压裂改造模拟器开展水平井水力压裂物理模拟实验,研究了射孔参数(射孔类型、射孔相位角和射孔轴向间距)和页岩物理力学性质对压裂裂缝形态的影响;考佳玮等21采用真三轴水力压裂系统研究了不同注液方式(连续注液、变排量注液和交替注液)对裂缝形态的影响。大尺寸真三轴实验用压裂液泵注系统向模拟井眼注入流体,采用声发射监测仪器监测压裂过程中裂缝的扩展情况,试验后,可采用工业CT扫描或者将试样剖开的方法观察裂缝分布和形态,分析压裂裂缝的延伸机理30


1.2 数值模拟方法

页岩储层压裂过程中形成的裂缝形态复杂,甚至呈网状,应用常规的模型不能模拟页岩储层中压裂裂缝的延伸情况。近年来,国内外学者对常规模型进行改进,提出了一些适用于模拟页岩储层压裂的数值模拟方法,这些方法的基本原理、优势和局限性对比列入表1。

(1)边界元法(Boundary Element Method,BEM)。边界元法是在有限差分和有限元的影响下发展起来的数值模拟方法31。该方法基于控制微分方程的基本解建立边界积分方程来表达问题,通过对边界分元插值离散,化为代数方程组求解。

(2)非常规裂缝模型(Unconventional Fracture Model,UFM)。在拟三维模型的基础上,KRESSE等23-24考虑储层的应力场以及相邻裂缝间的应力干扰,建立了非常规裂缝模型。UFM模型的控制方程包括支撑剂的沉降方程、G—W裂缝延伸准则和裂缝变形方程。利用该模型可以计算出裂缝的宽度等参数。

(3)离散化缝网模型(Discrete Fracture Network,DFN)。离散化缝网模型最早由MEYER等32提出,是目前模拟页岩气水力压裂复杂缝网较为成熟的模型。该模型假设缝网为拟椭圆状球体,主裂缝垂直于最小水平主应力σh,并与次生裂缝相交形成正交裂缝网络25

(4)扩展有限元法(Extend Finite Element Method,XFEM)。扩展有限元法是在经典有限元法的基础上增加广义节点自由度,改进插值形函数提出的33-34,是目前解决不连续问题最准确的方法,在国内外得到广泛应用。

2 扩展有限元法模型建立及应用实例

2.1 模型建立

2.1.1 初始岩体建模

基于前人35物理实验数据,应用ABAQUS2017数值模拟软件建立水力压裂裂缝延伸的扩展有限元二维模型。模型尺寸为0.3 m×0.3 m,模拟射孔长度为0.03 m,选择射孔方向与最大水平主应力夹角为60°,其他模型参数如表2所示。

表2   模型参数[35]

Table 2  Model parameters [35]

弹性模量E/GPa泊松比单轴抗压强度/MPa抗拉强度/MPa渗透率(K)/(10-3 μm2)
8.4020.2328.342.590.1
孔隙度(φ)/%最大水平主应力(σH)/MPa最小水平主应力(σh)/MPa试验排量/(m3/s)压裂液表观黏度/(mPa·s)
1.85612.1×10-973

新窗口打开| 下载CSV


2.1.2 初始裂缝模拟

在水力压裂过程中,流体流经初始裂缝时,初始裂缝不发生偏转,因此可以用初始裂缝近似表示射孔长度。本文定义一条直线为初始裂缝。

扩展有限元模型中裂缝独立于网格存在,因此对岩石骨架模型进行网格划分时不需要考虑裂缝的状态,采用简单的四边形网格即可。本文应用的数值模拟软件ABAQUS自2017版本以后,要对初始裂缝进行网格划分,在相互作用中指定初始裂缝位置。

2.1.3 断裂起始准则

当应力或应变满足定义的起裂准则时,判断裂缝开始起裂。扩展有限元方法中常用的断裂起始准则有最大主应力准则和最大主应变准则36

最大主应力准则可以表示为:f=σmaxσmaxa (1)

式中:σmax为最大主应力,MPa;σmaxa为最大临界主应力,MPa;σmax=0, σmax<0  σmax,σmax0

最大主应变准则可以表示为:f=εmaxεmaxa (2)

式中:εmax为最大主应变;εmaxa为最大临界主应变。以上两式中,f为断裂标准,当1.0f1.0+ftol时发生断裂,ftol为公差,默认为0.05。

为了提高计算结果的收敛性,选择最大主应力准则作为判断裂缝起裂的依据。

2.1.4 应力平衡分析

水力压裂过程中,岩石受到的地应力平衡,故模型中建立地应力平衡分析。忽略岩石中流体的黏性,用虚功原理表示应力平衡方程如下37- 38VδεTDepdεdtdV+VδεTDepmso+poε3KsdpodtdV-VδεTmso+poξdpodtdV=VδuTdfdtdV+SδuTdτdtdS (3)

式中:δε为虚位移,m;Dep为弹塑性矩阵;m=[1,1,1,0,0,0];Ks为固体颗粒的压缩模量,Pa;so为固相含液饱和度;po为孔隙液体压力,Pa;f为岩石体力,Pa;τ为岩体的面力,Pa;ξ为毛细压力与饱和度关系的参数;δu为虚应变。

2.1.5 边界条件施加

(1)位移边界条件。因为本文建立的是水力压裂二维模型,所以要约束x、y方向的边界位移,即U1=0,U2=0(ABAQUS中“1”代表x方向,“2”代表y方向)。

(2)孔压边界条件。定义孔压边界条件po=pob,即边界孔压为pob。若为静压边界条件,则定义po=0。

(3)初始地应力条件。定义场条件中的初始地应力,设置水平最大主应力、水平最小主应力等。

(4)初始孔隙比。定义场条件的初始孔隙比,模型中默认Void Ratio=φφ为模型孔隙度,%。

2.1.6 流体注入模拟

忽略流体的可压缩性,选择注入点后,修改模型输入文件,通过子命令*cflow设置流体注入。假设流体在岩体中流动遵循达西渗流规律,而在裂缝内流体为牛顿流体,其流动遵循如下规律39

Q=-ktp

式中:Q为试验排量, m3/s;kt为流体流动系数;p为流体流动压力,Pa。

2.2 模拟结果验证

模拟结果如图1(a)所示,当射孔方向与最大水平主应力方向的夹角为60°时,裂缝在初始裂缝两侧地层中均发生偏转,转向方向沿垂直最小水平主应力方向,转向后沿平行最大水平主应力方向延伸,且初始裂缝两侧的裂缝并不对称,右侧转向半径略大于左侧。姜浒等35物理试验结果如图1(b)所示,可见模拟结果与试验结果十分接近。

图1

图1   水力压裂裂缝延伸形态

Fig.1   Fracture extension pattern of hydraulic fracturing


2.3 扩展有限元法应用
2.3.1 不同射孔方向水力裂缝的延伸形态

用2.1.1中参数建立不同射孔角度的扩展有限元模型(图2),模拟水力压裂过程中压裂裂缝的延伸情况。设置射孔方向与最大水平主应力方向的夹角α分别为0°、15°、30°、45°、60°和75°。

图2

图2   模型射孔方向与最大水平主应力夹角示意

Fig.2   Angle between model perforation direction and the maximum horizontal principal stress


模型计算结果如图3,当射孔方向与最大水平主应力的夹角α=0°时,裂缝沿最小水平主应力方向(即阻力最小的方向)裂开,之后沿着最大水平主应力方向延伸,不发生偏转;当α不为0°时,由于最大、最小水平主应力存在应力差,导致裂缝发生不同程度的偏转。由图3可知,射孔方向与最大水平主应力的夹角α越大,裂缝偏转角度越小,偏转距离越大。

图3

图3   不同射孔方向扩展有限元模型计算结果

Fig.3   Calculation results of finite element model for different perforation directions


模型计算完成后进行后处理即可得到初始断裂单元的压力随时间的变化曲线(图4),不同射孔角度的压力随时间变化趋势相似:第一阶段为压力急剧增加阶段,直到峰值。随着压裂液注入,初始断裂单元压力逐渐上升,达到最大值时初始断裂发生,此时的压力称为破裂压力。第二阶段为台阶式锯齿状下降阶段。初始断裂发生后,压力下降,待压裂液充满初始裂缝后,压力又上升达到下一个峰值,此时压裂液突破主裂缝的某一点,形成次级裂缝,如此往复。第三阶段为持续稳定阶段。在该阶段,裂缝稳定延伸,压力曲线趋于平稳。

图4

图4   不同射孔方向下压力随时间的变化曲线

Fig.4   Curve of pressure under different perforation directions changing with time


图4可知,射孔方向与最大水平主应力方向的夹角α不同,破裂压力不同。当α=0°时,破裂压力最小;当α不为0°时,破裂压力增大。采用不同射孔方向,初始断裂单元达到破裂压力的注入时间不同,裂缝稳定延伸的压力也不同,除45°外,射孔方向与最大水平主应力方向的夹角越大,裂缝稳定延伸的压力也越大;夹角为45°时,裂缝延伸的压力最小。

提取图4中不同射孔方向的破裂压力得到破裂压力随射孔方向变化关系曲线(图5),破裂压力随射孔方向与最大水平主应力之间夹角α的增大呈增大趋势。当α=0°时,即射孔方向与最大水平主应力方向一致,此时裂缝沿着最大水平主应力方向延伸不发生偏转,裂缝受到最小水平主应力作用而起裂,破裂压力最小;当α不为0°时,裂缝朝最大主应力方向发生明显的偏转,最终沿着平行于最大水平主应力方向延伸。裂缝起裂受到最大最小水平主应力的共同作用,且裂缝转向需要消耗更多的能量,破裂压力增大。随射孔方向与最大水平主应力之间夹角α从0°增加到75°,破裂压力从5.98 MPa增大至17.36 MPa。上述趋势说明当射孔方向与最大水平主应力方向存在夹角时,由于水平应力差的存在,导致裂缝受到的阻力增大,使破裂压力遵循最小能量消耗原理变化,而裂缝也朝能量消耗最小的方向延伸。

图5

图5   破裂压力随射孔方向变化关系曲线

Fig.5   Relation curve of fracturing pressure with perforation direction


综合图4图5可知,当射孔方向与最大水平主应力的夹角小于45°时,裂缝较早达到破裂压力,破裂压力较小且裂缝稳定延伸时的压力也较小。实际压裂工程中,选择射孔方向时,可优先考虑与最大水平主应力夹角小于45°的方向。沿着该方向射孔,可以适当减小注液压力和压裂液的排量,有利于形成稳定延伸的压裂裂缝。

影响水力压裂裂缝延伸的因素还有注入压裂液排量、水平应力差等。实际压裂工程中,可以根据模拟结果选择合适的射孔方向,以达到较好的压裂效果。

2.3.2 分段顺序压裂水力裂缝的扩展形态

以射孔方向与最大水平主应力夹角45°为例,创建分段压裂扩展有限元模型,用2条直线表示2次压裂产生的初始裂缝(图6),模型参数同2.1节。在单裂缝扩展模型基础上,增加泄压分析和二次压裂分析,模拟分段顺序水力压裂的裂缝扩展形态。

图6

图6   分段压裂初始模型

Fig.6   Initial model of staged fracturing


模型计算结果(图7)显示,二次压裂压裂液注入初期,初始裂缝靠近一次压裂裂缝尖端的一侧先开裂,并沿垂直最小水平主应力方向发生偏转,偏转后存在靠近一次压裂裂缝延伸的趋势;且在二次压裂初期,一次压裂裂缝尖端与二次压裂初始裂缝之间形成一条新的短裂缝[图7(a)]。随着压裂液不断注入,初始裂缝两侧均开裂并延伸:靠近一次压裂裂缝尖端一侧,主裂缝偏向一次压裂裂缝延伸,初期形成的短裂缝垂直最小水平主应力发生偏转,最终二者相交并停止扩展;远离一次压裂裂缝一侧,主裂缝的延伸情况与一次压裂裂缝的延伸情况一致,垂直最小水平主应力偏转后沿平行最大水平主应力方向延伸[图7(b)]。

图7

图7   分段压裂模型计算结果

Fig.7   Calculation results of staged fracturing model


由此可知,一次压裂的水力裂缝周围产生了诱导应力,改变了岩体原始应力分布。一次压裂产生的水力裂缝尖端出现应力集中,因此二次压裂过程中靠近一次压裂裂缝一侧受到较大的地应力作用首先开裂,偏转后偏向应力集中方向延伸。由于存在一次压裂裂缝和二次压裂主裂缝尖端的应力集中现象,导致二者之间的应力值明显大于附近的地应力,因此在压裂初期形成一个较小的次生裂缝。次生裂缝与主裂缝相交后,因受到诱导应力影响,导致应力差不足以满足裂缝延伸所需,该侧裂缝停止扩展。远离一次压裂裂缝一侧受到诱导应力的影响很小,可忽略不计,因此这一侧主裂缝的延伸符合单次压裂裂缝延伸的规律,即初期裂缝起裂后垂直最小水平主应力偏转,最终沿平行最大水平主应力方向延伸。

综上可知,分段压裂过程中,先成裂缝产生的诱导应力会使后续压裂裂缝的延伸方向发生改变,并对其延伸存在一定的阻碍作用,且这种影响随距离的增加而减小。因此在实际工程中,可适当增大压裂段之间的距离以使裂缝充分延伸,增加裂缝的导流能力,进而改善增产效果。

3 结论

(1)压裂裂缝的扩展及分布是影响水力压裂效果的关键,常采用物理实验和数值模拟开展水力压裂模拟。物理实验模拟过程较准确,得到的结果更直观,但代表性有限,经济性不强;数值模拟主要有边界元法、非常规裂缝模型、离散化缝网模型和扩展有限元法等,各有优势和不足,需进一步改进才能更加真实地模拟页岩储层水力压裂过程。

(2)扩展有限元模型适用范围广,求解速度快,精度高,是页岩储层水力压裂模拟的主要发展方向。应用扩展有限元方法模拟水力压裂裂缝延伸情况与物理实验结果十分接近,表明该方法有效可靠。

(3)应用扩展有限元法模拟不同射孔方向下压裂裂缝的延伸情况,得到初始断裂单元压力随时间的变化曲线,表明压裂液注入到结束压力变化分3个阶段:急剧增加阶段、锯齿状下降阶段和持续稳定阶段。

(4)射孔方向与最大水平主应力的夹角越大,裂缝偏转角度越小,偏转距离越大,岩石破裂压力也越大;夹角小于45°时,具有较小的破裂压力和裂缝稳定延伸压力,是实际工程中较为理想的选择。

(5)分段顺序压裂过程中,先成裂缝产生的诱导应力对后续压裂裂缝的延伸存在一定的不利影响,可适当增大压裂段间距以弱化其影响,对分段压裂设计具有一定指导意义。

参考文献

邹才能,杨智,朱如凯,等. 中国非常规油气勘探开发与理论技术进展[J]. 地质学报,2015,89(6):979-1007.

[本文引用: 1]

ZOU C N,YANG Z,ZHU R K,et al. Progress in China's unconventional oil & gas exploration and development and theoretical technologies[J]. Acta Geologica Sinica,2015,89(6):979-1007.

[本文引用: 1]

董大忠,王玉满,李新景,等. 中国页岩气勘探开发新突破及发展前景思考[J]. 天然气工业,2016,36(1):19-32.

[本文引用: 1]

DONG D Z,WANG Y M,LI X J,et al. Breakthrough and prospect of shale gas exploration and development in China[J]. Natural Gas Industry,2016,36(1):19-32.

[本文引用: 1]

郑哲敏. 关于中国页岩气持续开发工程科学研究的一点认识[J]. 科学通报,2016,61(1):34-35.

[本文引用: 1]

ZHENG Z M. Some knowledge on engineering scientific research of shale gas sustainable development in China[J]. Science Bulletin,2016,61(1):34-35.

[本文引用: 1]

续震,卢鹏. 国内水力压裂技术现状[J]. 石化技术,2017,24(5):280.

[本文引用: 1]

XU Z,LU P. Current situation of domestic hydraulic fracturing technology[J]. Petrochemical Industry Technology,2017,24 (5):280.

[本文引用: 1]

蒋廷学,贾长贵,王海涛,等. 页岩气网络压裂设计方法研究[J]. 石油钻探技术,2011,39(3):36-40.

[本文引用: 1]

JIANG T X,JIA C G,WANG H T,et al. Study on network fracturing design method in shale gas[J]. Petroleum Drilling Techniques,2011,39(3):36-40.

[本文引用: 1]

侯振坤,杨春和,王磊, 等. 大尺寸真三轴页岩水平井水力压裂物理模拟试验与裂缝延伸规律分析[J]. 岩土力学,2016,37(2):407-414.

[本文引用: 2]

HOU Z K,YANG C H,WANG L,et al. Hydraulic fracture propagation of shale horizontal well by large-scale true triaxial physical simulation test[J]. Rock and Soil Mechanics,2016,37(2):407-414.

[本文引用: 2]

周治东,程万,魏子俊,等. 基于BEM的水力裂缝起裂与扩展数值模拟[J]. 地球物理学进展,2020,35(2):807-814.

[本文引用: 1]

ZHOU Z D,CHENG W,WEI Z J,et al. Numerical simulation of hydraulic fracture initiation and propagation based on BEM[J]. Progress in Geophysics,2020,35(2): 807-814.

[本文引用: 1]

梁天成,付海峰,刘云志,等. 水力压裂裂缝扩展声发射破裂机制判定方法研究[J]. 实验力学,2019,34(2):358-364.

LIANG T C,FU H F,LIU Y Z,et al. On the determination method of rupture mechanism in acoustic emission used in hydraulic fracturing fracture propagation[J]. Journal of Experimental Mechanics,2019,34(2):358-364.

梁冰,岳鹭飞,孙维吉. 页岩矿物组分对裂缝扩展影响的数值模拟分析[J]. 海相油气地质,2019,24(4):97-101.

LIANG B,YUR L F,SUN W J. The influence of shale mineral composition on crack growth:Numerical simulation[J]. Marine Origin Petroleum Geology,2019,24(4):97-101.

侯冰,陈勉,谭鹏,等. 页岩气藏缝网压裂物理模拟的声发射监测初探[J]. 中国石油大学学报:自然科学版,2015,39(1):66-71.

HOU B,CHEN M,TAN P,et al. Monitoring of hydraulic fracture network by acoustic emission method in simulated tri-axial fracturing system of shale gas reservoirs[J]. Journal of China University of Petroleum:Edition of Natural Science,2015,39(1):66-71.

曾义金,周俊,王海涛,等. 深层页岩真三轴变排量水力压裂物理模拟研究[J]. 岩石力学与工程学报,2019,38(9):1758-1766.

ZENG Y J,ZHOU J,WANG H T,et al. Research on true triaxial hydraulic fracturing in deep shale with varying pumping rates[J]. Chinese Journal of Rock Mechanics and Engineering,2019, 38(9):1758-1766.

LLANOS E M,JEFFREY R G,HILLIS R,et al. Hydraulic fracture propagation through an orthogonal discontinuity: A laboratory,analytical and numerical study[J]. Rock Mechanics and Rock Engineering,2017,50(8):2101-2118.

[本文引用: 1]

孙峰,唐梅荣,张翔,等. 水平井定面射孔近井筒的破裂形态[J]. 天然气工业,2019,39(4):62-68.

[本文引用: 1]

SUN F,TANG M R,ZHANG X,et al. Fracture geometry near the wellbore of a horizontal well with in-plane perforation[J]. Natural Gas Industry,2019,39(4):62-68.

[本文引用: 1]

刘奉银,王少宾,徐楠. 定面射孔水力压裂裂缝扩展研究[C]∥中国力学学会,北京理工大学. 中国力学大会,2017暨庆祝中国力学学会成立60周年大会论文集.北京:中国力学学会,2017:1791-1800.

LIU F Y,WANG S B,XU N. Study on fracture expansion of hydraulic fracturing with fixed perforation[C]∥Chinese Society of Theoretical and Applied Mechanics, Beijing Institute of Technology. Proceedings of Proceedings of China mechanical congress,2017 and celebration of the 60th anniversary of the founding of the Chinese society of Mechanics. Beijing: Chinese Society of Theoretical and Applied Mechanics, 2017:1791-1800.

刘合,兰中孝,王素玲,等. 水平井定面射孔条件下水力裂缝起裂机理[J]. 石油勘探与开发,2015,42(6):794-800.

[本文引用: 1]

LIU H,LAN Z X,WANG S L,et al. Hydraulic fracture initiation mechanism in the definite plane perforating technology of horizontal well[J]. Petroleum Exploration and Development,2015,42(6):794-800.

[本文引用: 1]

陈勉,庞飞,金衍. 大尺寸真三轴水力压裂模拟与分析[J]. 岩石力学与工程学报,2000,19(S1):868-872.

[本文引用: 1]

CHEN M,PANG F,JIN Y. Experiments and analysis on hydraulic fracturing by a large-size triaxial simulator[J]. Chinese Journal of Rock Mechanics and Engineering,2000,19(S1):868-872.

[本文引用: 1]

STANCHITS S,SURDI A,GATHOGO P. Onset of hydraulic fracture initiation monitored by acoustic emission and volumetric deformation measurements[J].Rock Mechanics and Rock Engineering,2014,47(5):1521-1532.

[本文引用: 1]

LECAMPION B,DESROCHES J,JEFFREY R G,et al. Experiments versus theory for the initiation and propagation of radial hydraulic fractures in low-permeability materials[J]. Journal of Geophysical Research: Solid Earth,2017,122(2):1239-1263.

DEHGHAN A N,GOSHTASBI K,AHANGARI K,et al. Experimental investigation of hydraulic fracture propagation in fractured blocks[J]. Bulletin of Engineering Geology and the Environment,2015,74(3):887-895.

[本文引用: 1]

解经宇,蒋国盛,王荣璟,等. 射孔对页岩水力裂缝形态影响的物理模拟实验[J]. 煤炭学报,2018,43(3):776-783.

[本文引用: 1]

JIE J Y,JIANG G S,WANG R J,et al. Experimental investigation on the influence of perforation on the hydraulic fracture geometry in shale[J]. Journal of China Coal Society,2018,43(3):776-783.

[本文引用: 1]

考佳玮,金衍,付卫能,等. 注液方式对深层页岩裂缝形态影响实验研究[J]. 地下空间与工程学报,2019,15(1):32-38.

[本文引用: 1]

KAO J W,JIN Y,FU W N,et al. Experiments research on the influence of injection method on hydraulic fracture morphology of deep shale[J]. Chinese Journal of Underground Space and Engineering,2019,15(1):32-38.

[本文引用: 1]

陈铭,张士诚,胥云,等. 基于互补算法的水力裂缝与天然裂缝相互作用边界元模型[J]. 岩石力学与工程学报,2018,37(S2):3947-3957.

CHEN M,ZHANG S C,XU Y,et al. Boundary element model for mechanical interaction between hydraulic fracture and natural fracture based on complementary algorithm[J].Chinese Jou-rnal of Rock Mechanics and Engineering,2018,37(S2):3947-3957.

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

[本文引用: 1]

GU H,COHEN C,WENG X,et al. Modeling of hydraulic-fracture-network propagation in a naturally fractured formation[J]. SPE Production & Operations,2011,26(4):368-380.

[本文引用: 1]

卞晓冰,侯磊,蒋廷学,等. 深层页岩裂缝形态影响因素[J]. 岩性油气藏,2019,31(6):161-168.

[本文引用: 1]

BIAN X B,HOU L,JIANG T X,et al. Influencing factors of fracture geometry in deep shale gas wells[J]. Lithologic Reservoirs,2019,31(6):161-168.

[本文引用: 1]

赵猛,范锡彦. 页岩气缝网优化的数值模拟[J]. 大庆石油地质与开发,2019,38(4):167-174.

ZHAO M,FAN X Y. Optimized numerical simulation of the fracture network in the shale gas[J]. Petroleum Geology & Oilfield Development in Daqing,2019,38(4):167-174.

欧阳伟平,孙贺东,韩红旭.致密气藏水平井多段体积压裂复杂裂缝网络试井解释新模型[J].天然气工业,2020,40(3):74-81.

OUYANG W P,SUN H D,HAN H X. A new well test interpretation model for complex fracture networks in horizontal wells with multi-stage volume fracturing in tight gas reservoirs[J]. Natural Gas Industry,2020,40(3):74-81.

徐加祥,丁云宏,杨立峰,等. 基于扩展有限元的水力压裂缝间干扰及裂缝形态分析[J]. 天然气地球科学,2018,29(9):1356-1363.

XU J X,DING Y H,YANG L F,et al. Analysis of stress interference and geometry of hydraulic fractures based on the extended finite element method[J].Natural Gas Geoscience,2018,29(9):1356-1363.

曾青冬,姚军. 水平井多裂缝同步扩展数值模拟[J]. 石油学报,2015,36(12):1571-1579.

ZENG Q D, YAO J. Numerical simulation of multiple fra-ctures simultaneous propagation in horizontal wells[J]. Acta Petrolei Sinica,2015,36(12):1571-1579.

衡帅,杨春和,曾义金,等. 页岩水力压裂裂缝形态的试验研究[J]. 岩土工程学报,2014,36(7):1243-1251.

[本文引用: 1]

HENG S,YANG C H,ZENG Y J,et al. Experimental study on hydraulic fracture geometry of shale[J]. Chinese Journal of Geotechnical Engineering,2014,36(7):1243-1251.

[本文引用: 1]

HOSSAIN M M,RAHMAN M K. Numerical simulation of complex fracture growth during tight reservoir stimulation by hydraulic fracturing[J]. Journal of Petroleum Science & Engineering,2008,60(2):86-104.

[本文引用: 1]

MEYER B R,BAZAN L W,JACOT R H. Optimization of Multiple Transverse Hydraulic Fractures in Horizontal Wellbores[C]// SPE Unconventional Gas Conference. SPE 131732. Richardson: Society of Petroleum Engineers,2010.

[本文引用: 1]

MOES N,DOLBOW J,BELYTSCHKO T. A finite element method for crack growth without remeshing[J]. International Journal for Numerical Methods in Engineering,1999,46(1):131-150.

[本文引用: 1]

BELYTSCHKO T,SUKUMAR N,DAUX C,et al. Arbitrary branched and intersecting cracks with the extended finite element method[J]. International Journal for Numerical Methods in Engineering,2000,48(12):1741-1760.

[本文引用: 1]

姜浒,陈勉,张广清,等. 定向射孔对水力裂缝起裂与延伸的影响[J]. 岩石力学与工程学报,2009,28(7):1321-1326.

[本文引用: 4]

JIANG H,CHEN M,ZHANG G Q,et al. Impact of oriented perforation on hydraulic fracture initiation and propagation[J]. Chinese Journal of Rock Mechanics and Engineering,2009,28(7):1321-1326.

[本文引用: 4]

高红,郑颖人,冯夏庭. 岩土材料最大主剪应变破坏准则的推导[J]. 岩石力学与工程学报,2007,26(3):518-524.

[本文引用: 1]

GAO H,ZHENG Y R,FENG X T. Deduction of failure criterion for geomaterials based on maximum principal shear strain[J]. Chinese Journal of Rock Mechanics and Engineering,2007,26(3):518-524.

[本文引用: 1]

张晓咏,戴自航. 应用ABAQUS程序进行渗流作用下边坡稳定分析[J]. 岩石力学与工程学报,2010,29(S1):2927-2934.

[本文引用: 1]

ZHANG X Y,DAI Z H. Analysis of slope stability under seepage by using ABAQUS program[J]. Chinese Journal of Rock Mechanics and Engineering,2010,29(S1):2927-2934.

[本文引用: 1]

龚迪光,曲占庆,李建雄,等. 基于ABAQUS平台的水力裂缝扩展有限元模拟研究[J]. 岩土力学,2016,37(5):1512-1520.

[本文引用: 1]

GONG D G,QU Z Q,LI J X,et al. Extended finite element simulation of hydraulic fracture based on ABAQUS platform[J]. Rock and Soil Mechanics,2016,37(5):1512-1520.

[本文引用: 1]

薛炳,张广明,吴恒安,等. 油井水力压裂的三维数值模拟[J].中国科学技术大学学报,2008,38(11):1322-1325,1347.

[本文引用: 1]

XUE B,ZHANG G M,WU H A,et al. Three-dimensional numerical simulation of hydraulic fracture in oil wells[J]. Journal of University of Science and Technology of China,2008,38(11):1322-1325,1347.

[本文引用: 1]

/