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

页岩气压裂缝复杂裂缝网络动态特性分析 传统的压力裂裂缝模拟通常不考虑自然裂缝的存在，并在井眼两侧形成垂直于最小主应力方向的对称双侧形裂缝。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) 缝