页岩气压裂缝复杂裂缝网络动态特性分析

页岩气压裂缝复杂裂缝网络动态特性分析

传统的压力裂裂缝模拟通常不考虑自然裂缝的存在，并在井眼两侧形成垂直于最小主应力方向的对称双侧形裂缝。th。但从直接的矿井压裂测试［２］到间接的微地震监测［３］,都表明在天然裂缝发育的页岩气储层,会形成不对称、不规则的复杂裂缝网络。开启和扩展天然裂缝是形成缝网的基础。众多学者通过实验和理论分析,研究了人工裂缝与天然缝相交后的物理现象,并建立了裂缝延伸的判定准则［４－５］;Warpinski［６］采用线摩擦理论考虑作用在裂缝面的剪应力引起的剪切滑移破坏,并以摩尔库伦准则分析作用在裂缝面的正应力引起的张性破坏。Beugelsdijk等［７］通过实验分析了水平应力差、排量和黏度对分支缝延伸的影响。研究表明,与天然缝相交后形成的裂缝网络复杂程度不仅与地应力相关,而且受到岩石力学参数、天然缝参数、压裂施工参数以及工作液物性所影响［８－１０］。压裂过程中人工缝与天然缝相交后,分支缝的起裂与延伸机理对形成复杂裂缝网络的几何尺寸和复杂程度有重要影响。笔者以断裂力学理论为基础,引入起裂与延伸判定计算准则;在此基础上,考虑附加应力场影响,建立页岩气压裂复杂裂缝网络模型,在求解过程中以等压力梯度值划分裂缝延伸方向上的计算节点,将节点压力作为关键变量显式求解裂缝网络几何尺寸参数;并将起裂与延伸准则作为分支缝扩展依据,避免常规拟三维裂缝模拟中对缝长延伸的复杂拟合过程。与线网模型［１１］,该模型能够充分考虑储层中任意分布的天然缝对整个裂缝网络结构的影响,井眼两侧的缝网整体形态是非对称且内部的分支缝非均匀间距,更符合微地震测试结果。与非常规裂缝模型相比［１２］,建立的起裂与延伸判断准则时考虑了剪应力与正应力共同作用下的裂缝尖端周向应力,能够合理解释延伸过程中裂缝转向的物理现象并能计算分支缝扩展延伸中的角度变化;并且通过改进数值求解方法能有效提高计算速度。1裂缝形成和延长条件分析1.1第一、二代椭圆完整积分在地应力作用下椭圆形天然缝表面的剪应力所产生的Ⅱ型应力强度因子为:式中K(κ)和E(κ)分别为第一、二类椭圆完整积分;τ为剪应力,MPa;υ为泊松比;ω为剪应力方向与椭圆长轴夹角;θ为椭圆裂缝边缘任意点与长轴夹角;a和b分别为裂缝长、短半轴长度,m。1.2应力强度因子在水力裂缝延伸过程中,不断扩展的裂缝长度和缝内压力重新分布将导致由正应力引起的裂缝尖端应力强度因子动态变化。根据初始起裂位置沿着裂缝延伸方向上按照等压力梯度划分节点,再使用叠加原理得到在不同节点压力作用下的裂缝尖端应力强度因子的值。压力节点划分如图1所示。由不同节点压力所形成的应力强度因子为［１３］:式中l（ｉ）为第i个压力节点与初始位置的间距;l为人工裂缝长度;pｎｅｔ（ｉ）为第i个节点处净压力。在此基础上可得到正应力作用下的Ⅰ型应力强度因子为:式中n为划分的压力节点总数。1.3应变弹性影响系数多条分支裂缝在延伸的同时,相互之间会受到作用在临近裂缝的正、剪应力影响,而形成附加应力场。按照Crouch等人提出的一种解析计算公式,并经Ol-son［１４］改进后,其表达式为:式中Cｎｓｉｊ、Cｎｎｉｊ分别为平面应变弹性影响系数;Dｓｊ和Dｎｊ分别是各个微段上剪应力和正应力引起的不连续位移(即基于流固耦合计算而得到的裂缝宽度);Gｉｊ为Ol-son提出的三维修正系数,引入缝高和裂缝间距对附加应力场影响。在数值计算中的每一时步都需计算该附加应力场,然后将其作用在地应力场上,用于进行延伸和起裂压力计算公式和迭代求解不同节点处的裂缝几何参数。2裂缝起裂判定当压裂液进入天然缝,只要缝内压力大于裂缝面正应力及岩石抗张强度,将会开启天然缝并形成具有一定宽度的分支裂缝。根据裂缝面上正应力计算公式:式中T０为岩石抗张强度,MPa;σ为裂缝面正应力,MPa;σ２为垂向应力,MPa;σ１、σ３分别为水平应力,MPa;α、β、γ分别为裂缝面法线与主应力方向夹角。可得到不同节点处的净压力为:将不同节点处的净压力带入式(3)可得到裂缝起裂判定计算方程:式中m为裂缝开启时对应的节点数,m≤n。基于Erdogan提出的最大周向应力理论［１５］,可推导Ⅰ-Ⅱ复合型裂缝延伸角度方程:以及Ⅰ-Ⅱ复合型裂缝延伸判定计算方程:3缝网接收点流量的重新平均分配模型假设条件:①天然裂缝都是垂直分布;②在无天然裂缝区域设定地层是均质、各向同性的连续线弹性体;③储层厚度较大,不发生穿层现象;④裂缝在垂直方向上的延伸速度小于在长度方向上的延伸速度。复杂裂缝网络模型包含以下控制方程。1)起裂与延伸判定公式:见式(7)、(9)。2)缝宽方程假设产层足够厚,缝内净压pｎｅｔ=pｆ(x,t)-σ′,σ′=σ+σｎｉ,得到裂缝横截面上任意z位置的缝宽方程为:3)压降方程根据Nolte平板间压降公式,并引入管道形状ue788ｊ(n)因子的方法,可得到缝内压降方程:4)高度方程5)连续性方程模型初、边值条件与拟三维模型类似,但在裂缝相交后,交点处流量需重新平均分配。式中xｊ为裂缝长度,m;wｊ(x,z,t)为缝宽,m;hｊ(x,t)为裂缝高度,m;pｎｅｔ(x,t)为缝内压力,MPa;qｊ(x,t)为流量,m３/min;E为杨氏模量,GPa;n为流态指数,无因次;K为稠度系数,Pa·sｎ;Cｔ为滤失系数,m·min－０．５;Aｊ为裂缝横切面积,m２;t为施工时间,min;τ(x)为t时刻液体到达缝内x处的时间,min;下标j代表缝网中某条分支缝。改进的复杂裂缝模型求解方法如下所述。1)先以二维PKN模型的缝长公式计算在不同时步下的初始裂缝长度:,并根据井底或交点压力作为初始条件在缝长方向上以等压降值划分计算节点:,将得到的节点压力带入起裂判定方程计算实际能够开启的裂缝长度。2)根据节点处的压力值,带入缝高方程求出不同节点处的缝高。3)将得到的缝高和该节点处的压力带入缝宽方程可求解不同节点处的缝宽。4)将前面求得的几何尺寸和初始流量分布带入压降方程,可得到节点处的新的压力值,此时需要对比前面等压降节点处的差值,并调整和拟合初始节点压力,使其满足|p′ｎｅｔ(i)-pｎｅｔ(i)|≤ε。5)进一步将一维流动的连续性方程按照不同节点对应的参数值带入相应的缝宽、缝高值,以求出不同节点处的流量分布,将该流量值与初始假设的节点流量分布比较。若不一致,则改变初始流量分布关系式,重复第1步至第5步,最终实现流量拟合|q′（ｉ）-q（ｉ）|≤η。6)在时步更新之前,需要考虑裂缝延伸判定条件,用该时刻的裂缝尖端方位、缝内压力值和附加应力场为已知条件判定裂缝是否会在下一个时刻继续扩展,并计算动态变化的延伸角。若满足条件则进入下一时步,重复步骤第1步至第6步,直至裂缝停止扩展。针对页岩气压裂实际过程提出裂缝网络模拟程序流程如图2所示。4模拟结果分析4.1缝内压力增大,部位变化根据延伸角计算公式,可分析水平地应力差、缝内净压和夹角对分支缝延伸轨迹的影响(图3)。当水平应力差不变、缝内净压力的增大,延伸角度减小,分支缝延伸轨迹的曲率半径增大,不仅能有效扩大缝网宽度,也易于在空间上与其他缝相遇。而当净压力不变、水平应力差逐步增大时,延伸角变化越大,延伸轨迹的曲率半径减小,裂缝面在较小的范围内迅速垂直于最小主应力方向,使得缝内压力递减迅速,延伸或起裂难度加大。为提高压裂改造效果,应选取天然缝较为发育或者水平地应力差较小的适宜体积压裂的有利区块或层位,并通过优化施工参数和改进压裂材料,提高缝内净压力。4.2岩石化系统裂缝的发育根据裂缝延伸判断准则,定量分析水平应力差、裂缝半径、夹角、岩石断裂韧性对延伸所需压力的影响(图4)。模拟结果表明:决定天然缝能否延伸的关键因素是自身长度和方位角,而岩石断裂韧性的变化对裂缝延伸压力的影响不明显。因此对于页岩气长水平井段压裂应选择天然缝发育程度较高,密度较大的区域射孔;在布井时不一定完全按照最小主应力方向,要充分考虑天然缝的走向对整个裂网形态影响。4.3岩相色谱-井底密度特征根据页岩气压裂的“大液量、大排量、低黏度、低滤失”等特点,结合某海相页岩气井地质、施工参数,模拟了复杂裂缝网络形成过程和各分支裂缝几何尺寸,具体输入参数如下:水平最大主应力为37.5MPa,水平最小主应力为34.5MPa,杨氏模量为23.37GPa,泊松比为0.25,岩石抗张强度为3MPa,压裂液黏度为30mPa·s(施工/模拟可控关键变量),压裂液密度为1.02g/cm３,综合滤失系数为0.0009m·min－０．５(施工/模拟可控关键变量),岩石断裂韧性为1.21MPa·m０．５,平均施工排量为8.5m３(施工/模拟可控关键变量),井底压力为55MPa(施工/模拟可控关键变量),计算时间步长为2min。天然缝参数位置随机分布,手动输入,与最小水平主应力方向夹角为±0°~40°,倾角为90°。将整个缝网与微地震解释结果比较,几何尺寸能够实现与微震数据点分布范围的较好拟合,而且能够以裂缝复杂程度的差异形态解释数据点分布密度的变化(图5)。对比验证了模型计算结果,并能初步用于压裂设计优化。5模拟结果验证1)通过引入起裂与延伸判断准则、建立复杂裂缝网络模型、改进数值计算方法,实现了对含有大量天然缝的页岩气藏压裂裂缝网络的模拟研究,并结合微地震测试结果验证模型的准确性。2)