![页岩气压裂缝复杂裂缝网络动态特性分析_第1页](http://file4.renrendoc.com/view/6744e0d22f608edc816590fdea88f527/6744e0d22f608edc816590fdea88f5271.gif)
![页岩气压裂缝复杂裂缝网络动态特性分析_第2页](http://file4.renrendoc.com/view/6744e0d22f608edc816590fdea88f527/6744e0d22f608edc816590fdea88f5272.gif)
![页岩气压裂缝复杂裂缝网络动态特性分析_第3页](http://file4.renrendoc.com/view/6744e0d22f608edc816590fdea88f527/6744e0d22f608edc816590fdea88f5273.gif)
![页岩气压裂缝复杂裂缝网络动态特性分析_第4页](http://file4.renrendoc.com/view/6744e0d22f608edc816590fdea88f527/6744e0d22f608edc816590fdea88f5274.gif)
![页岩气压裂缝复杂裂缝网络动态特性分析_第5页](http://file4.renrendoc.com/view/6744e0d22f608edc816590fdea88f527/6744e0d22f608edc816590fdea88f5275.gif)
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
页岩气压裂缝复杂裂缝网络动态特性分析
传统的压力裂裂缝模拟通常不考虑自然裂缝的存在,并在井眼两侧形成垂直于最小主应力方向的对称双侧形裂缝。th。但从直接的矿井压裂测试[2]到间接的微地震监测[3],都表明在天然裂缝发育的页岩气储层,会形成不对称、不规则的复杂裂缝网络。开启和扩展天然裂缝是形成缝网的基础。众多学者通过实验和理论分析,研究了人工裂缝与天然缝相交后的物理现象,并建立了裂缝延伸的判定准则[4-5];Warpinski[6]采用线摩擦理论考虑作用在裂缝面的剪应力引起的剪切滑移破坏,并以摩尔库伦准则分析作用在裂缝面的正应力引起的张性破坏。Beugelsdijk等[7]通过实验分析了水平应力差、排量和黏度对分支缝延伸的影响。研究表明,与天然缝相交后形成的裂缝网络复杂程度不仅与地应力相关,而且受到岩石力学参数、天然缝参数、压裂施工参数以及工作液物性所影响[8-10]。压裂过程中人工缝与天然缝相交后,分支缝的起裂与延伸机理对形成复杂裂缝网络的几何尺寸和复杂程度有重要影响。笔者以断裂力学理论为基础,引入起裂与延伸判定计算准则;在此基础上,考虑附加应力场影响,建立页岩气压裂复杂裂缝网络模型,在求解过程中以等压力梯度值划分裂缝延伸方向上的计算节点,将节点压力作为关键变量显式求解裂缝网络几何尺寸参数;并将起裂与延伸准则作为分支缝扩展依据,避免常规拟三维裂缝模拟中对缝长延伸的复杂拟合过程。与线网模型[11],该模型能够充分考虑储层中任意分布的天然缝对整个裂缝网络结构的影响,井眼两侧的缝网整体形态是非对称且内部的分支缝非均匀间距,更符合微地震测试结果。与非常规裂缝模型相比[12],建立的起裂与延伸判断准则时考虑了剪应力与正应力共同作用下的裂缝尖端周向应力,能够合理解释延伸过程中裂缝转向的物理现象并能计算分支缝扩展延伸中的角度变化;并且通过改进数值求解方法能有效提高计算速度。1裂缝形成和延长条件分析1.1第一、二代椭圆完整积分在地应力作用下椭圆形天然缝表面的剪应力所产生的Ⅱ型应力强度因子为:式中K(κ)和E(κ)分别为第一、二类椭圆完整积分;τ为剪应力,MPa;υ为泊松比;ω为剪应力方向与椭圆长轴夹角;θ为椭圆裂缝边缘任意点与长轴夹角;a和b分别为裂缝长、短半轴长度,m。1.2应力强度因子在水力裂缝延伸过程中,不断扩展的裂缝长度和缝内压力重新分布将导致由正应力引起的裂缝尖端应力强度因子动态变化。根据初始起裂位置沿着裂缝延伸方向上按照等压力梯度划分节点,再使用叠加原理得到在不同节点压力作用下的裂缝尖端应力强度因子的值。压力节点划分如图1所示。由不同节点压力所形成的应力强度因子为[13]:式中l(i)为第i个压力节点与初始位置的间距;l为人工裂缝长度;pnet(i)为第i个节点处净压力。在此基础上可得到正应力作用下的Ⅰ型应力强度因子为:式中n为划分的压力节点总数。1.3应变弹性影响系数多条分支裂缝在延伸的同时,相互之间会受到作用在临近裂缝的正、剪应力影响,而形成附加应力场。按照Crouch等人提出的一种解析计算公式,并经Ol-son[14]改进后,其表达式为:式中Cnsij、Cnnij分别为平面应变弹性影响系数;Dsj和Dnj分别是各个微段上剪应力和正应力引起的不连续位移(即基于流固耦合计算而得到的裂缝宽度);Gij为Ol-son提出的三维修正系数,引入缝高和裂缝间距对附加应力场影响。在数值计算中的每一时步都需计算该附加应力场,然后将其作用在地应力场上,用于进行延伸和起裂压力计算公式和迭代求解不同节点处的裂缝几何参数。2裂缝起裂判定当压裂液进入天然缝,只要缝内压力大于裂缝面正应力及岩石抗张强度,将会开启天然缝并形成具有一定宽度的分支裂缝。根据裂缝面上正应力计算公式:式中T0为岩石抗张强度,MPa;σ为裂缝面正应力,MPa;σ2为垂向应力,MPa;σ1、σ3分别为水平应力,MPa;α、β、γ分别为裂缝面法线与主应力方向夹角。可得到不同节点处的净压力为:将不同节点处的净压力带入式(3)可得到裂缝起裂判定计算方程:式中m为裂缝开启时对应的节点数,m≤n。基于Erdogan提出的最大周向应力理论[15],可推导Ⅰ-Ⅱ复合型裂缝延伸角度方程:以及Ⅰ-Ⅱ复合型裂缝延伸判定计算方程:3缝网接收点流量的重新平均分配模型假设条件:①天然裂缝都是垂直分布;②在无天然裂缝区域设定地层是均质、各向同性的连续线弹性体;③储层厚度较大,不发生穿层现象;④裂缝在垂直方向上的延伸速度小于在长度方向上的延伸速度。复杂裂缝网络模型包含以下控制方程。1)起裂与延伸判定公式:见式(7)、(9)。2)缝宽方程假设产层足够厚,缝内净压pnet=pf(x,t)-σ′,σ′=σ+σni,得到裂缝横截面上任意z位置的缝宽方程为:3)压降方程根据Nolte平板间压降公式,并引入管道形状ue788j(n)因子的方法,可得到缝内压降方程:4)高度方程5)连续性方程模型初、边值条件与拟三维模型类似,但在裂缝相交后,交点处流量需重新平均分配。式中xj为裂缝长度,m;wj(x,z,t)为缝宽,m;hj(x,t)为裂缝高度,m;pnet(x,t)为缝内压力,MPa;qj(x,t)为流量,m3/min;E为杨氏模量,GPa;n为流态指数,无因次;K为稠度系数,Pa·sn;Ct为滤失系数,m·min-0.5;Aj为裂缝横切面积,m2;t为施工时间,min;τ(x)为t时刻液体到达缝内x处的时间,min;下标j代表缝网中某条分支缝。改进的复杂裂缝模型求解方法如下所述。1)先以二维PKN模型的缝长公式计算在不同时步下的初始裂缝长度:,并根据井底或交点压力作为初始条件在缝长方向上以等压降值划分计算节点:,将得到的节点压力带入起裂判定方程计算实际能够开启的裂缝长度。2)根据节点处的压力值,带入缝高方程求出不同节点处的缝高。3)将得到的缝高和该节点处的压力带入缝宽方程可求解不同节点处的缝宽。4)将前面求得的几何尺寸和初始流量分布带入压降方程,可得到节点处的新的压力值,此时需要对比前面等压降节点处的差值,并调整和拟合初始节点压力,使其满足|p′net(i)-pnet(i)|≤ε。5)进一步将一维流动的连续性方程按照不同节点对应的参数值带入相应的缝宽、缝高值,以求出不同节点处的流量分布,将该流量值与初始假设的节点流量分布比较。若不一致,则改变初始流量分布关系式,重复第1步至第5步,最终实现流量拟合|q′(i)-q(i)|≤η。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/cm3,综合滤失系数为0.0009m·min-0.5(施工/模拟可控关键变量),岩石断裂韧性为1.21MPa·m0.5,平均施工排量为8.5m3(施工/模拟可控关键变量),井底压力为55MPa(施工/模拟可控关键变量),计算时间步长为2min。天然缝参数位置随机分布,手动输入,与最小水平主应力方向夹角为±0°~40°,倾角为90°。将整个缝网与微地震解释结果比较,几何尺寸能够实现与微震数据点分布范围的较好拟合,而且能够以裂缝复杂程度的差异形态解释数据点分布密度的变化(图5)。对比验证了模型计算结果,并能初步用于压裂设计优化。5模拟结果验证1)通过引入起裂与延伸判断准则、建立复杂裂缝网络模型、改进数值计算方法,实现了对含有大量天然缝的页岩气藏压裂裂缝网络的模拟研究,并结合微地震测试结果验证模型的准确性。2)
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论