, 卢义玉1,2, 葛兆龙1,2, 仲建宇1,21. 重庆大学 煤矿灾害动力学与控制国家重点实验室, 重庆 400044;
2. 重庆大学 复杂煤气层瓦斯抽采国家地方联合工程实验室, 重庆 400044
收稿日期：2016-01-06
基金项目：“十二五”国家科技重大专项资金资助项目(2011ZX05065)；国家自然科学基金资助项目(51374258)；****和创新团队计划项目(IRT13043)。
作者简介：程玉刚(1991-)，男，河南安阳人，重庆大学博士研究生;
卢义玉(1972-)，男，湖北京山人，重庆大学教授，博士生导师, 教育部“****奖励计划”特聘教授。
摘要：基于孔隙水压力梯度对孔壁及裂缝尖端应力作用, 提出了煤层导向压裂扩展方法.通过研究有效应力变化, 建立了导向压裂起裂压力及裂缝扩展压力计算模型.采用RFPA3D-Flow流固耦合软件建立了孔隙压力梯度场下的导向压裂数值模型.对7组不同数值模型计算分析, 研究了导向孔的布置方式、距离、控制水压三个主要因素对导向压裂作用效果.研究结果表明:孔隙压力梯度能有效降低压裂孔起裂压力及裂缝扩展压力.起裂压力及裂缝扩展压力降低幅度与导向孔控制水压呈正相关关系, 与导向孔距离呈负相关关系.数值计算得到的起裂压力变化规律与理论计算预测相符合, 论证了导向压裂起裂压力计算模型的正确性以及导向压裂裂缝控制方法的可行性.
关键词：水力压裂导向压裂起裂压力孔隙水压力裂缝扩展
Influence of Pore Water Pressure Gradient on Guiding Hydraulic Fracturing in Underground Coal Mine
CHENG Yu-gang1,2
, LU Yi-yu1,2, GE Zhao-long1,2, ZHONG Jian-yu1,21. State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, China;
2. National & Local Joint Engineering Laboratory of Gas Drainage in Complex Coal Seam, Chongqing University, Chongqing 400044, China
Corresponding author: CHENG Yu-gang, E-mail: cyg19910213@163.com
Abstract: Based on the influence of pore water pressure gradient on the stress of hole wall and crack tip, a new crack-guiding method with guiding drilling is proposed. According to research on changing of the effective stress, the calculation model of the fracturing pressure and fracture propagation pressure is established. By using RFPA3D-Flow fluid solid coupling software, the numerical model of oriented fracturing under pore water pressure gradient field is established. 7 groups of different guiding drilling models are calculated to study the guide fracturing effect of three main factors, viz. arrangement of guide borehole, distance, and water pressure. The results show that the pore pressure gradient can effectively reduce the crack initiation pressure and crack propagation pressure. The pressure of initiation and crack propagation is positively correlated with the control pressure of the guiding hole and inverse correlation between the distance from the guide hole. The variation law of the initial crack pressure calculated by the numerical calculation is consistent with the theoretical prediction. The study demonstrates that the guiding fracturing pressure calculation model is correct and the control method of the fracture orientation is feasible.
Key Words: hydraulic fracturingguiding fracturinginitiation pressurepore water pressurecrack propagation
近年来, 水力压裂技术凭借其快速、大范围改善煤层透气性的特点成为了煤矿井下水治瓦斯技术的研究热点[1].随着矿井的开采深度增加地应力越来越大, 水力压裂时的起裂压力升高, 导致煤体内某些区域产生高应力集中, 这些应力集中区则成为了煤与瓦斯突出的隐患.另外, 由于常规水力压裂形成的裂缝面往往是一条主裂缝, 在煤层中沿某单一方向增透, 导致裂缝两侧形成增透空白带.煤矿井下水力压裂的目的是在减少对顶底板损伤的情况下, 尽可能实现均衡增透, 因此水压裂缝的起裂压力及裂缝数量、扩展方向是判断压裂效果优劣的直接因素.
陈勉、刘勇等研究发现, 射孔技术形成的微裂隙环会降低起裂压力, 且对裂缝扩展有一定导向能力[2-3], 但是射孔技术主要应用于油气系统, 煤层相对于油气储层较薄, 若将射孔压裂技术直接应用于煤矿井下, 难度较高, 操作复杂.目前煤矿井下已有的导向水力压裂技术主要有预置缝槽定向压裂[4]和定向孔控制压裂[5].前者在裂缝脱离预置缝槽的控制范围后仍会沿最大地应力方向发生转向；后者需要大量卸压孔来进行导向, 施工量较大.可以发现, 现有导向压裂技术虽然能在一定程度上提高压裂效果, 但也存在自身方法的不完善性.Brnuo等[6]通过对含双缺口的岩样进行孔隙水压力试验, 发现有效的孔隙水压会对裂缝的扩展路径有影响, 冷雪峰等[7]通过数值分析对该试验进行了验证.以上试验研究均表明孔隙水压力是渗流损伤作用机制中的重要影响参数, 但是关于煤矿井下起裂压力及裂缝扩展压力的研究成果中, 一般忽略了孔隙水压力影响, 因此在水力压裂相关领域中, 采用孔隙水压力梯度来控制煤层导向压裂还有待研究.
基于以上分析, 本文提出了煤矿井下导向孔诱导起裂及裂缝扩展方法, 即在压裂孔周围布置一定数量的导向孔, 导向孔在压裂前提前注入高压水.本文拟通过分析煤层内部形成孔隙水压力梯度对压裂孔围岩应力影响, 建立导向压裂起裂压力及裂缝扩展压力计算模型, 采用RFPA3D-Flow对导向起裂压力变化规律及裂缝诱导作用进行验证分析.
1 导向压裂起裂及扩展机理1.1 孔隙水压力梯度对起裂压力的影响假设穿层钻孔与煤层面垂直, 则钻孔受力示意图如图 1所示, 其中σv为竖直方向主应力, σh, max, σh, min分别为水平方向最大、最小主应力, 煤层倾角为α, 煤层倾向与最大水平主应力夹角为β.对于远场地应力转化为坐标系(X, Y, Z)下压裂孔围岩应力状态可参考文献[8], 本文不再赘述.
图 1(Fig. 1)
 | 图 1 压裂孔受力示意图Fig.1 Stress state of borehole |
通过压裂孔围岩应力计算得到钻孔坐标系(r, θ, Z)中孔壁处应力状态:
 | (1) |
在图 1沿A-A面将煤层剖切, 如图 2所示.
图 2(Fig. 2)
 | 图 2 A-A剖面Fig.2 A-A profile |
根据厚壁平面径向流理论, 假设在无限大地层中, 导向孔保持稳定注入高压水, 其影响范围为R, 则距离导向孔d处的压力分布方程及边界条件可写成:
 | (2) |
解得导向孔附近压力分布为
 | (3) |
 | (4) |
 | (5) |
 | (6) |
 | (7) |
 | (8) |
根据式(8) 表达式, 起裂压力pf降低幅度与导向孔注水压力pw呈正相关关系, 与压裂孔距导向孔之间的距离d呈负相关关系.当pw越大, d值越小, 起裂压力降低幅度越大; 反之, 当pw越小, d值越大, 起裂压力降低幅度越小.
1.2 孔隙压力梯度对扩展压力的影响煤层内可能会存在天然裂隙或节理, 当水压裂缝沟通这些天然裂隙时, 网络化裂缝范围更广, 以下分析裂缝在受孔隙水压力梯度场影响时的扩展压力.假设忽略裂缝宽度, 裂缝受力情况如图 3所示.
图 3(Fig. 3)
 | 图 3 裂缝受力示意图Fig.3 Stress state of a stressed crack |
将图 3中的坐标系(X, O, Y)转化为裂缝坐标系(X′, O, Y′), 得到裂缝所受远场应力状态为
 | (9) |
 | (10) |
 | (11) |
 | (12) |
 | (13) |
由式(13) 可知, 孔隙水压力梯度会降低裂缝扩展压力, 能够诱导裂缝向导向孔扩展.导向孔的水压与距离影响效果表现出与压裂孔起裂压力相同的规律.由于裂缝扩展的自组织行为, 为了减小扩展所需要的能量, 裂缝会沿扩展压力较低方向扩展.
2 数值计算分析2.1 模型模拟条件的确定为了验证以上理论分析, 本文采用RFPA3D-Flow来研究煤岩破裂问题, 模拟复杂三维应力下的流固耦合难题[8], 定性地研究在不同孔隙压力梯度条件下水力压裂的变化规律, 得到水压裂缝的起裂、诱导扩展特征.
设计如图 4所示压裂模型，模型在煤层走向和倾向的长度均为20m, 煤层及顶底板均厚5m.压裂钻孔位于模型的对称中心轴线, 长度贯穿整个煤层, 压裂孔与导向孔之间距离为d, 其钻孔半径均为0.2m.压裂孔初始水头为500m(每100 m水头代表 1MPa水压), 每一计算步增加20m水头, 导向孔保持水压为p.由于岩层均质度比煤层均质度高, 设置顶底板均质度m=10, 煤层均质度m=5.根据中国大陆浅层地壳实测地应力分布规律研究[9], 设定模型在埋深500m左右时X方向加载最大水平主应力20MPa, Y方向加载最小水平主应力10MPa, Z方向加载中间主应力13.5MPa.模型中涉及的煤岩层材料具体参数设置如表 1所示.
图 4(Fig. 4)
 | 图 4 数值模型(弹性模量)Fig.4 Numerical model (elastic modulus) |
表 1(Table 1)
 表 1 材料的物理力学参数Table 1 Physico-mechanical parameters of material
| 表 1 材料的物理力学参数 Table 1 Physico-mechanical parameters of material |
2.2 方案设计根据理论研究, 首先针对导向孔是否可以真正对裂缝起诱导作用进行计算分析.分别设计了常规压裂与右侧含有一个导向孔(d=4m, p=5MPa)的2组模型进行空白对照, 裂缝扩展形态如图 5所示.可以看出常规压裂裂缝面基本沿最大地应力方向(X方向)扩展, 由于材料的非均质性, 裂缝会因为细观强度的差别发生些许偏转, 符合以往学者的研究.而当有导向孔时, 裂缝变为向导向孔位置扩展, 在孔隙水压力梯度作用下, 裂缝扩展方向得到明显改变.
图 5(Fig. 5)
 | 图 5 常规压裂与导向压裂裂缝对比Fig.5 Crack propagation of ordinary fracturing and guiding fracturing (a)—无导向孔；(b)—有导向孔. |
由于现场实施压裂条件有限, 导向孔数量过多、控制距离过远对于压裂泵的注水流量和注水时间要求太高, 结合其他学者研究结果[7], 设计如表 2所示7组计算模型.其中, A, B, C为导向孔布置方式对照组; B, D, E为导向孔控制距离对照组; B, F, G为控制水压对照组.每个模型中导向孔的距离、水压、数量如表 2所示.
表 2(Table 2)
表 2 导向孔参数Table 2 Control drilling parameters
| 表 2 导向孔参数 Table 2 Control drilling parameters |
2.3 结果分析2.3.1 导向孔数量分析为方便观察裂缝形态, 沿模型煤层中部(Z=7.5m)做剖切, 从+Z方向向-Z方向俯视, 观察不同导向孔布置方式下的裂缝扩展形态, 如图 6所示.
图 6(Fig. 6)
 | 图 6 导向孔布置方式对裂缝扩展影响Fig.6 Influence of controlling drilling arrangement on crack propagation (a)—模型A；(b)—模型B; (c)—模型C. |
模型A中主要表现为一条裂缝, 扩展方向与诱导方向一致, 增透面积有限.模型B中4个导向孔在压裂孔附近呈正四边形分布, 形成的分支裂缝较多, 增透范围均匀.结合模型A和模型B, 模型C设置6个导向孔, 对比模型B发现, 导向孔数量增加后裂缝数量反而减少, 且诱导作用没有模型B效果好.A, B, C 3组对比结果发现导向孔并非数量越多诱导越明显, 因此现场实施导向压裂时根据增透需求制定不同的导向孔布置数量(方式).
2.3.2 导向孔控制距离、控制水压分析根据各模型计算结果, 得到不同导向孔条件下压裂孔起裂压力, 见表 3.结果显示, 随着导向孔距离减小或导向孔控制水压升高, 压裂孔的起裂压力逐渐降低, 这与计算模型预测结果相符合.
表 3(Table 3)
表 3 不同导向孔模型下压裂孔起裂压力Table 3 Initiation pressure of different controlling-drilling models
| 表 3 不同导向孔模型下压裂孔起裂压力 Table 3 Initiation pressure of different controlling-drilling models |
裂缝的扩展形态是考察导向压裂的主要标准.图 7为不同模型在裂缝扩展趋于稳定时的裂缝形态.结合表 3发现, 虽然模型D起裂压力最低, 但是由于导向孔控制距离过小, 水压裂缝在越过导向孔后扩展逐渐受最大主应力方向控制, 表现出常规水压裂缝扩展的特征, 裂缝网控制范围较小.而模型E由于导向孔距离较远, 裂缝在前期扩展时主要受地应力影响, 沿X方向扩展, 导致主裂缝数量较少, 裂缝在后期接近导向孔时才向导向孔偏转, 在导向孔的控制区域内裂缝扩展数量没有模型B多, 且起裂压力相对较高.模型F, 模型G裂缝扩展形态类似, 均由于导向孔控制水压偏低, 部分导向孔失效, 裂缝在导向孔的控制区域内形成的裂缝网分布不均匀.
图 7(Fig. 7)
 | 图 7 不同导向孔模型下水压裂缝扩展结果Fig.7 Crack propagation in different controlling-drilling models (a)—模型B；(b)—模型D; (c)—模型E；(d)—模型F; (e)—模型G. |
对比5组模型计算结果可知, 导向孔只有保持适当的导向孔控制距离和水压, 才能促进地应力与孔隙水压力梯度对裂缝扩展的竞争作用, 更有利于裂缝衍生出分支裂缝.
3 结论1) 通过分析导向孔形成的孔隙水压力梯度对孔壁及裂缝尖端应力影响, 建立了导向压裂起裂压力及裂缝扩展压力计算模型.
2) 孔隙水压力梯度能有效降低起裂压力及裂缝扩展压力, 形成诱导作用.起裂压力及裂缝扩展压力降低幅度与导向孔注水压力呈正相关关系, 与导向孔距离呈负相关关系.当导向孔控制水压越大, 距离越近, 起裂压力降低幅度越大; 反之, 起裂压力降低幅度越小.
3) 采用RFPA3D-FLOW建立了三维地应力条件下的渗流-应变导向压裂数值模型, 通过7组不同导向孔条件, 分析了导向孔布置方式、水压、距离对导向压裂的影响效果.数值计算得到的起裂压力变化规律及诱导作用与理论分析一致, 为煤矿井下降低起裂压力及控制水压裂缝扩展方向提供了一定参考价值.
参考文献
| [1] | Tian K Y, Zheng J Y. The application of hydraulic fracturing outburst prevention measures[J].Procedia Engineering, 2011, 26: 495–500.DOI:10.1016/j.proeng.2011.11.2197 |
| [2] | 陈勉, 金衍, 张广清. 石油工程岩石力学[M]. 北京: 科学出版社, 2008: 56-282. ( Chen Mian, Jin Yan, Zhang Guang-qing. Oil engineering rock mechanics[M]. Beijing: Science Press, 2008: 56-282.) |
| [3] | 刘勇, 卢义玉, 魏建平, 等. 降低井下煤层压裂起裂压力方法研究[J].中国安全科学学报, 2013, 23(9): 98–100. ( Liu Yong, Lu Yi-yu, Wei Jian-ping, et al. A novel method for decreasing coal seam fracture initiation pressure[J].China Safety Science Journal, 2013, 23(9): 98–100.) |
| [4] | Yan F Z, Lin B Q, Zhu C J, et al. A novel ECBM extraction technology based on the integration of hydraulic slotting and hydraulic fracturing[J].Journal of Natural Gas Science and Engineering, 2015, 22: 571–579.DOI:10.1016/j.jngse.2015.01.008 |
| [5] | 李全贵, 翟成, 林柏泉, 等. 定向水力压裂技术研究与应用[J].西安科技大学学报, 2011, 31(6): 735–739. ( Li Quan-gui, Zhai Cheng, Lin Bai-quan, et al. Research and application of directional hydraulic fracturing technology[J].Journal of Xi'an University of Science and Technology, 2011, 31(6): 735–739.) |
| [6] | Bruno M S, Nakagawa F M. Pore pressure influence on tensile fracture propagation in sedimentary rock[J].International Journal of Rock Mechanics and Mining Sciences, 1991, 28(4): 261–273.DOI:10.1016/0148-9062(91)90593-B |
| [7] | 冷雪峰, 唐春安, 李连崇, 等. 非均匀孔隙压力下水压致裂的数值试验[J].东北大学学报(自然科学版), 2003, 24(3): 288–291. ( Leng Xue-feng, Tang Chun-an, Li Lian-chong, et al. Numerical simulation of hydrofracturing in anisotropic pore pressure fields[J].Journal of Northeastern University(Natural Science), 2003, 24(3): 288–291.) |
| [8] | Li L C, Tang C A, Fu Y F. Influence of heterogeneity on fracture behavior in multi-layered materials subjected to thermo-mechanical loading[J].Computational Materials Science, 2009, 46(3): 667–671.DOI:10.1016/j.commatsci.2009.03.041 |
| [9] | 景峰, 盛谦, 张勇慧, 等. 中国大陆浅层地壳实测地应力分布规律研究[J].岩石力学与工程学报, 2007, 26(10): 2056–2062. ( Jing Feng, Sheng Qian, Zhang Yong-hui, et al. Research on distribution rule of shallow crustal geostress in China mainland[J].Chinese Journal of Rock Mechanics and Engineering, 2007, 26(10): 2056–2062.DOI:10.3321/j.issn:1000-6915.2007.10.014) |
