摘要：深层页岩气主要通过水平井分段多簇密切割压裂进行增产，但由于深层地应力变化较快以及缝间应力的干扰导致水力裂缝长短不一，亟需优化调控。针对深层页岩压裂多簇裂缝起裂时间不同、延伸不均匀的现象，建立了水力裂缝异步起裂延伸模型，同时考虑了各簇射孔位置地应力差异和缝间应力干扰对水力裂缝异步起裂、非均匀延伸的影响。模拟结果表明：采用常规均匀射孔工艺，射孔簇位置地层破裂压力越高，该簇裂缝起裂时间越晚，外侧裂缝延伸更快，裂缝延伸非均匀程度较大;采用非均匀限流射孔工艺，通过合理调整各簇射孔孔眼数量，可使各簇裂缝起裂时间与延伸速度差异明显减小，保证压裂结束时各簇裂缝长度基本一致，实现“抑长促短，均匀延伸”的目的。现场监测各簇进液比例与模型计算结果较为一致，验证了模型的可靠性。研究成果可有效降低目前非均匀限流射孔工艺参数设计的盲目性，对深层页岩气压裂优化设计具有重要的理论价值和指导作用。
关键词：深层页岩;密切割压裂;非均匀射孔;数值模拟;矿场应用
0引言
中国页岩气资源储量丰富，开发潜力巨大[1]，但必须通过水平井分段多簇压裂技术实现有效开发[2⁃3]。有超过60%的页岩气资源埋藏于3500m以下的深层[4]，开发难度较大。深层页岩气压裂通常采用密切割方式，即通过增加射孔簇数、缩小簇间距，增强水力裂缝之间干扰效应，促进缝网形成，提高增产效果[5]。
目前，国内外学者已围绕页岩气水平井分段多簇压裂水力裂缝延伸开展了大量研究。何青等[6]、衡峰等[7]研究发现多簇水力裂缝起裂和延伸存在竞争干扰机制，中间裂缝受到挤压而开度受限、延伸不足。史吉辉等[8]基于有限元方法模拟了多簇水力裂缝的平面扩展的过程，发现簇间距是影响缝间干扰的最主要因素。薛明宇等[9]结合黏聚单元方法与扩展有限单元法，架构了多簇水力裂缝内的二维流固耦合模型，分析了压裂施工方案对裂缝扩展的影响规律。曾庆磊等[10]基于扩展有限元法模拟岩石中任意簇裂缝沿任意路径转向延伸，并考虑各簇裂缝间流量动态分配。Y.Cheng[11]利用边界元模型研究发现平行裂缝缝间应力集中会导致不对称延伸。谭超[12]运用位移不连续方法，并通过多物理场耦合迭代方法模拟多裂缝干扰延伸。R.Lin等[13]、陈曦宇[14]基于位移不连续法、有限差分法与有限体积法，建立完全流固耦合的多簇裂缝延伸数值模型。时贤等[15]基于多层压裂流量动态分配思想，建立了多簇裂缝同步扩展数学模型。
深层页岩地层构造复杂，地应力变化快，各射孔簇破裂压力通常存在差异[16]，水力裂缝无法同时起裂，甚至无法起裂[17⁃18]。此外，密切割压裂缝间应力干扰效应更强，裂缝起裂后延伸速度差异较大，非均匀程度更大[19]。目前裂缝延伸研究主要针对常规中浅层页岩，而针对深层页岩的研究相对较少，同时在多缝延伸模拟研究过程中较少考虑地应力的变化[20⁃22]。本文通过建立深层页岩气水平井压裂水力裂缝异步起裂延伸模型，在多簇裂缝流量分配方程中引入起裂压力，考虑各簇射孔破裂压力差异与缝间应力干扰对水力裂缝异步起裂、非均线延伸的影响，突破了传统裂缝延伸模型的局限性，对深层页岩气压裂优化设计起到一定的指导作用。
1数学模型
深层页岩气水平井分段多簇密切割缝网压裂过程中，多条裂缝从不同射孔簇起裂并延伸，而深层地应力变化较快，各射孔簇位置的岩石破裂压力可能存在差异，因此各条水力裂缝起裂时间有所不同，存在有先有后的“异步起裂”现象。此外，由于缝间应力干扰效应，裂缝起裂后的延伸速度也各不相同，延伸路径可能发生一定程度的偏转，最终导致各条水力裂缝非均匀延伸(图1)。
基于岩石力学、断裂力学、流体力学等基础理论，建立了深层页岩气水平井压裂水力裂缝异步起裂延伸模型，包括物质平衡方程、缝内流动方程、裂缝开度方程、裂缝高度方程、裂缝转向方程，以及各簇裂缝流量分配方程。模型的主要基本假设条件包括：地层岩石为线弹性均匀介质;缝内流体为牛顿不可压缩液体;每簇射孔仅起裂并延伸一条水力裂缝。该模型将重点考虑各簇射孔位置地层破裂压力差异与缝间应力干扰对水力裂缝异步起裂延伸行为的影响。
1.7计算方法与流程
联立方程式(1)~式(10)，结合各簇水力裂缝延伸边界条件与初始条件方程式(15)，利用有限差分方法求解物质平衡方程与缝内流体压降方程;利用边界元方法求解裂缝离散单元应力—应变平衡方程，求取裂缝开度;利用牛顿迭代方法构造雅各布矩阵，求解各簇水力裂缝流量分配方程。即可计算得到各簇水力裂缝延伸过程中几何参数，包括半长、宽度、高度，以及各簇水力裂缝分配所得流量，并可通过压裂结束后各簇水力裂缝延伸半长变异系数定量表征裂缝延伸非均匀程度。
2矿场应用与分析
利用本文建立的水力裂缝异步起裂延伸数值表征模型，在川南深层页岩气区块开展了矿场应用，模拟了XH1⁃3井第24段压裂过程中水力裂缝的非均匀延伸行为。
2.1目标井概况
XH1⁃3井位于渝西区块西山构造内，地处川东南坳褶带，构造西北翼末端为威远—龙女寺构造群，构造东南翼隔向斜与东山和西温泉构造平行，构造西南段与新店子构造呈低鞍状相对。XH1⁃3井为水平井，完钻井深6370m，垂深4178m，目标层位为龙马溪组，地层岩石弹性模量为47.22GPa，泊松比0.23，地层温度118.8℃，储层压力75.21MPa，压力系数1.72。XH1⁃3井使用外径139.7mm套管完井，水平段1516m，采用分段多簇压裂工艺进行增产改造。水平段压裂共分为25段，每段射孔3―6簇，簇间距8~13m，每簇射孔孔眼8―12个。
相关知识推荐：大庆测井人员的石油论文怎样发表
由测井数据解释成果可知，该井水平段地层应力变化较快。其中，第24号压裂段内的地层水平最小主应力变化幅度超过2MPa，可能造成该段压裂过程中各簇水力裂缝延伸长度严重不均。因此，本文以XH1⁃3井第24号压裂段为例，进行水力裂缝异步起裂延伸数值模拟。该压裂段地应力与压裂工程参数见表1。
2.2模拟结果与分析
基于该井段地质与工程参数，利用本文建立的深层页岩气水平井压裂水力裂缝异步起裂延伸模型对均匀射孔条件下的各簇水力裂缝起裂延伸行为进行模拟。首先，忽略深层页岩特有的地应力变化特征，假设每簇射孔位置的地层最小水平主应力相同。通过数值模拟计算，得到均匀地层应力下压裂结束时水力裂缝三维空间分布情况(图2)，以及压裂期间各条水力裂缝延伸半长和时间的关系(图3)。
由图2、图3的模拟计算结果可知，即使忽略地应力变化，由于缝间应力干扰效应的影响，各簇裂缝也呈现出了非均匀延伸的特征，总体来看，内侧的2、3、4号簇裂缝延伸受限，长度较短，而外侧的1、5号簇裂缝延伸较长。但因为各簇射孔位置地应力相同，各簇裂缝均从压裂开始时立刻起裂延伸，不存在先后依次异步起裂的现象。
考虑深层页岩地应力变化特征，根据表1中的各簇射孔地应力实际数据取值，通过数值模拟计算，得到非均匀地层应力下压裂结束时水力裂缝三维空间分布(图4)，以及压裂期间各条水力裂缝延伸半长和时间关系(图5)。
由图4和图5的模拟计算结果可知，若考虑深层页岩地应力非均匀特征，该段压裂各簇裂缝半长平均121.4m，但长短各异，呈现出更为显著的非均匀特征。其中，由于1号簇位置的地层最小水平主应力最高，比相邻的2号簇高1.48MPa，导致1号簇裂缝破裂压力更高，起裂时间最晚，直到压裂开始后78min时，该簇裂缝才开始起裂延伸，最终延伸半长也最短，仅有55.3m;此外，2、5号簇地层最小水平主应力略高于3、4号簇，起裂时间也较晚，分别为压裂开始后的36min、14min。但在各簇裂缝延伸过程中，由于缝间应力干扰效应，位于内侧的3、4号簇裂缝延伸受到外侧2、5号簇裂缝的限制，虽然起裂时间较早，但延伸速度较慢，最终延伸半长仅分别为110.5m和120.1m，而2、5号簇裂缝最终延伸半长分别为153.3m和167.6m。压裂结束时，各簇水力裂缝半长的变异系数(标准差/平均值)为0.36，非均匀程度较大。
整体来看，由于目标深层页岩气地层应力变化较快，各簇裂缝起裂时间各不相同，起裂时间与破裂压力呈正相关关系，即破裂压力越高，裂缝起裂越晚。此外，由于缝间应力干扰效应，裂缝延伸速度也各不相同，内侧裂缝延伸较快，外侧裂缝延伸较慢，各簇裂缝呈现出异步起裂、非均匀延伸的现象。
为有效降低各簇裂缝非均匀程度，该段压裂可采用非均匀限流射孔工艺，通过调整各射孔簇的孔眼数量，从而改变各簇孔眼摩阻大小，优化各簇裂缝流量分配，促进裂缝均匀延伸。为此，以均匀射孔工艺下的裂缝延伸模拟结果为基础，将该压裂段1―5簇射孔数量由统一的8孔，调整为20、10、14、13、8孔。利用本文模型，对上述孔眼数量情况下各簇水力裂缝起裂延伸行为进行模拟，得到非均匀地层应力下非均匀射孔压裂结束时水力裂缝三维图如图6所示，压裂期间各条水力裂缝延伸半长随时间变化曲线如图7所示。
由图6和图7的模拟计算结果可知，采用非均匀限流射孔工艺后，该段压裂各簇裂缝半长平均130.6m，延伸更加均匀。其中，破裂压力最高的1号簇裂缝由于孔眼数量增多，摩阻下降，分得流量增多，压裂开始后39min即起裂延伸，比均匀射孔压裂提前了39min，且最终延伸半长达到了101.4m。此外，尽管2、5号簇裂缝起裂后延伸速度仍然较快，但由于孔眼数量比3、4号簇更少，起裂时间明显推迟，最终延伸半长有所缩短，与3、4号簇裂缝延伸半长基本相同。压裂结束时，各簇水力裂缝半长的变异系数为0.13，非均匀程度显著降低。但由于1号簇破裂压力过高，裂缝延伸长度仍然明显小于其他簇裂缝，仍然有进一步优化空间，如适当扩大该簇射孔孔径，或采用段内暂堵转向压裂工艺。
考虑地层应力变化裂缝异步起裂延伸的情况下，均匀射孔与非均匀射孔压裂各簇裂缝延伸半长对比如图8所示。此外，利用模型计算了该井段压裂期间各簇射孔进液比例，并与现场井下分布式光纤监测得到的压裂期间各簇射孔进液比例进行对比，如图9所示。
由图9可以看出，均匀射孔压裂1号簇进液比例明显较小，仅占总注入液量的3.6%，而5号簇进液比例过大，达到了39.7%。模型预测调整各簇孔眼数量之后，1号簇进液比例显著提高，达到了12.8%，其余各簇进液比例差距也明显减小。通过现场井下分布式光纤监测，各射孔簇进液比例变异系数从0.68显著降低至0.32，表明非均匀限流射孔工艺实现了XH1⁃3井第24段压裂各簇裂缝“抑长促短，均匀延伸”的目的。此外，现场井下分布式光纤监测结果与模型预测结果平均相对误差仅为7.1%，一致性较好，证实了本文模型的准确性。
3结论
(1)由于深层页岩地层应力变化快与缝间应力干扰效应，若采用常规的均匀射孔工艺进行压裂，各簇裂缝起裂时间与延伸速度各不相同。模拟发现，射孔簇位置地层破裂压力越高，该簇裂缝起裂时间越晚，且外侧裂缝延伸更快，裂缝延伸长段不一，非均匀程度显著，影响压裂增产效果。
(2)深层页岩压裂可采用非均匀限流射孔工艺，调整各簇射孔孔眼数量——适当增加高破裂压力簇与内侧簇射孔数量，从而有效降低各簇裂缝非均匀程度。模拟结果表明，各簇裂缝起裂时间与延伸速度差异明显减小，压裂结束时各簇裂缝半缝长变异系数仅约为0.1，延伸较为均匀一致。
(3)通过现场井下分布式光纤监测，各射孔簇进液比例变异系数从0.68显著降低至0.32，实现了“抑长促短，均匀延伸”的目的，且矿场实际监测的各簇进液比例与本文模型计算结果较为一致，验证了本文模型的可靠性。——论文作者：宋毅1林然2黄浩勇1任岚2岳文翰1孙映2
