非常规油气藏压裂水平井动态缝网模拟方法及应用.docx

《非常规油气藏压裂水平井动态缝网模拟方法及应用.docx》由会员分享，可在线阅读，更多相关《非常规油气藏压裂水平井动态缝网模拟方法及应用.docx（12页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、非常规油气藏压裂水平井动态缝网模拟方法及应用非常规油气藏压裂水平井进入生产期，压裂缝网形态并不是 静态不变的，而是随着开发过程的进行，缝内压力发生变化, 表现出动态缝网特征l-4o 2022年，Vermilye等5采 用生产期微地震监测手段，获取一口压裂水平井投产后由于 应力和流体扰动所产生的微地震事件点，证实了缝网的动态 变化特征，并提出了动态生产体积APV (activeproductionvolume)概念。监测结果说明储层中的 压裂改造范围在生产期迅速缩小，2年后的APV约为初始改 造体积的一半，3年后的APV几乎消失。这与目前非常规油 气藏生产井初期产量相对较高但随后快速递减并维持长

2、期 低产的生产规律相吻合6-7。压裂缝网的动态变化对生产期储层的供给能力和井的产量 大小具有重要影响，是开发调整方案制定和实施的重要依据 8-10L但目前尚未建立针对实际储层中缝网动态变化过 程的数值模拟方法，即目前用于油藏数值模拟的缝网模型保 持为压裂后的初始缝网形态，不随储层压力发生变化。本研 究在构建压裂缝网形变数学模型和数值模型的基础上，采用 传导率修正技术将裂缝形变模型与离散裂缝数值模拟器进 行耦合，研发了具有模拟动态缝网功能的数值模拟器，并开 展实例应用。图9数值模拟不同生产时期的动态缝网模型 a.第 7d； b.第 22d； c.第 44d； d.第 80d； e.第 146d；

3、 f.第226dFig. 9Dynamicfracturenetworkmodelsatdifferentproduct ionperiodsinthenumericalsimulator3.2实际井动态缝网模拟 为研究油田实际井的动态缝网变化特征及其对井产量的影 响，选取鄂尔多斯盆地某致密油藏EN125井开展动态缝网构 建和油藏数值模拟。EN125井水平段长780m,压裂段数7段, 泵注总液量1482nl3,总砂(陶粒)量176nl3,排量3m3/min, 井周围局部发育有天然裂缝。采用压裂模拟方法预测了该井 的压裂缝网模型，如图10a所示。图10图10EN125井压裂缝网模型(a)及对应的

4、油藏数值模拟网格模型(b)Fig. lOFracturenetworkmodel( a )andcorrespondingreservoirnumericalsimulationmeshmod el (b) forWellEN125图10b为EN125井离散裂缝数值模拟网格模型，基质渗透率 0. 36X10-3 um2,基质孔隙度8%,裂缝初始宽度范围110mll1, 原始油藏压力18MPa,流体粘度3. 2mPa* s,模型考虑了致 密油藏启动压力梯度，实验测试结果为0. 18MPa/mo拟合后 的压裂裂缝形变参数：剩余裂缝宽度系数Ores=0.1,裂缝 刚度中二2MPa,极限有效应力o2=

5、12MPa,油藏总正应力在 各网格间均匀分布on=18MPa。图11为模拟1年得到的EN125井在不同时期的动态缝网模 型。可见，裂缝动态闭合主要发生在生产初期的12个月内， 到后期逐渐放缓。整个过程中动态缝网的变化规律为：主裂 缝沟通的天然裂缝首先闭合，水力裂缝的宽度变小；而后水 力裂缝由两端向中间逐渐闭合。图11说明，压裂形成的初 始缝网在1年后会发生显著变化，缝网规模明显变小，与文 献中的实际井监测结果相符5o图11图11EN125井不同生产时期的动态缝网模型a.第1时间步；b.第3时间步；c.第15时间步；d.第31时 同步；e.第41时间步；f.第49时间步Fig. 1IDynami

6、cfracturenetworkmodelsatdifferentproduc tionperiodsforWellEN125为更直观研究动态缝网对压裂水平井产量的影响，采用离散裂缝数值模拟器对EN125井开展了两种情况下的模拟。一种 情况是考虑裂缝闭合，另一种情况是不考虑裂缝闭合，其他 参数均一致。将两种情况下的模拟产液量与该井实际产液量 进行拟合，结果如图12所示。图12 图12EN125井裂缝闭合对产量的影响Fig.12EffectoffractureclosureonWellEN125production图12说明，当不考虑裂缝闭合时，模拟产液量明显高于井 的实际产液量，特别是初期产量

7、吻合率较低，因为此期间储 层中裂缝闭合对产量的影响并未在模拟器中表达。而在相同 物性和流体参数条件下，当在模拟器中考虑裂缝闭合时，模 拟产液量与实际产液量差异较小，特别是初期产量较容易拟 合上。因此，动态缝网对致密油藏压裂水平井的产量影响显 著，在数值模拟过程中考虑动态缝网更加符合致密油藏生产 实际，开发方案设计精度也将更高。4结论1）压裂缝网是非常规油气藏中的主要渗流通道，油气井投 产后，由于压裂缝内压力发生变化，裂缝发生动态形变，表 现出动态缝网特征，对油气井产能具有重要影响。2）建立了表征裂缝形变的双曲数学模型及描述裂缝宽度非 均匀分布的二次函数模型，可同时模拟裂缝宽度和长度的动 态变化

8、。3）采用传导率修正方式，将裂缝形变数学模型与离散裂缝 数值模拟模型耦合，形成了具有模拟裂缝形变功能的离散裂 缝数值模拟器。4）建立了一口压裂水平井典型模型和一口致密油藏实际压 裂水平井模型，应用模拟器进行模拟，构建了不同生产时期的动态缝网，验证了动态缝网模拟功能和可靠性。1裂缝形变数学模型非常规油气藏水平井压裂形成缝网后，随着开采进行，裂缝 中流体压力会产生较大波动，加之裂缝面上岩石应力的作用, 导致开采过程中裂缝属性（宽度和渗透率）发生变化。实验 说明:11-13,压裂裂缝宽度随裂缝面上有效应力的变化遵 循双曲模型（图1）。曲线的变化分为3个阶段：随着裂缝面 上有效应力的增大（缝内压力降低

9、），初始阶段裂缝宽度呈 非线性急剧下降；裂缝空间压缩到一定程度后，支撑剂的抗 压强度发挥作用，此时继续增大有效应力那么使裂缝宽度呈缓 慢拟线性下降；当裂缝空间压缩逼近极限时，裂缝宽度对有 效应力变化不敏感，此时宽度变化进入平缓段。3个阶段在 横坐标上的分界点为。1 （有效应力，MPa）和。2,其中。2 为裂缝宽度下降至极限点的极限有效应力，MPa；对应的宽 度称为剩余裂缝宽度wr, m；当有效应力大于。2后，所建 立的数学模型规定裂缝宽度不随有效应力的变化而变化，保 持为wro图1图1裂缝宽度随有效应力变化的双曲模型示意图Fig. ISchematicsketchofhyperbolicmod

10、elshowingvariati onoffracturewidthwitheffectivestress将图1的双曲模型用如下公式进行定量表征： Dwf=wr（。n 三。2） wf=wmax 1 - B （onon+中）（0 o n o 2） wf=wmax（。nWO） wf=wr o nN。2wf=wmaxl- Bon。n+W Oono 2wf=wmax o nWO （ 1） 式中：wf为裂缝宽度，m； wmax为初始裂缝宽度，m； B为 与裂缝力学性能相关的无因次特征参数；中为裂缝刚度，MPa； on为裂缝面上有效应力，MPa,其计算公式为：。n二。normal - pf。n=。norm

11、al-pf (2)式中：0 nomial为裂缝面上的总正应力，MPa； pf为裂缝内 的流体压力，MPa。特征参数B的计算表达式为：3=(1- 0 res) (a 2+W。2) P=l-0resa2+W。2 (3)式中：Ores为无因次剩余裂缝宽度系数，即剩余裂缝宽度 与最大裂缝宽度的比值：。res=wrwmax。res=wrwmax (4)裂缝刚度甲与双曲模型曲线的曲率密切相关，中越大，曲率越小，意味着裂缝宽度下降的速度越小，曲线越接近线性, 如图2所示，图中计算参数分别为。2=20MPa, 0 res=O. lo图2图2裂缝变形双曲模型曲线形态与裂缝刚度之间的关系Fig. 2Hyperbo

12、liccurvesforfracturedeformationunderdif ferentfracturestiffness 在实际储层中，沿裂缝走向的裂缝宽度分布并不是一个定值, 而是呈“两头窄中间宽”的分布形态。为实现对裂缝变形更 加真实的模拟，先对裂缝进行线剖分，然后采用二次函数对 各个线单元的初始裂缝宽度wmax进行赋值（图3）。图3图3实际裂缝宽度的非均匀分布及其表征方法a.实际压裂裂缝形态；b.模拟裂缝形态（对裂缝进行线剖分）Fig. 3Non- uniformdistributionoffracturewidthalongar ealfractureandthecharacte

13、rizationmethodofthefractur ewidthdistribution将描述裂缝宽度随应力变化的双曲模型与描述裂缝宽度非 均匀分布的二次函数模型进行结合，那么不仅可模拟生产过程 中裂缝宽度的变化，还可以模拟裂缝长度的变化，使裂缝变 形模拟更贴近实际（图4）。图4图4模拟裂缝长度和宽度动态变化示意图Fig. 4Schematicsketchofthedynamicfracturewidthandlen gthmodeledinthesimulator2裂缝形变数值模拟方法采用基于非结构化网格的离散裂缝模型（DFM）开展裂缝形 变模拟，离散裂缝模型在模拟过程中保持了裂缝原始分布

14、和 形态，对裂缝形变过程的模拟具有较好适应性。技术思路为：通过传导率修正的方式，将裂缝形变数学模型与离散裂缝数 值模拟模型进行耦合，形成具有模拟裂缝形变功能的离散裂 缝数值模拟器。1.1 离散裂缝模型传导率计算在DFM模型中，任意两个相邻网格i和j之间的传导率可采 用以下公式计算14：Tij=TiTjTi+TjTij=TiTjTi+Tj (5)式中：Tij为i和j之间的全传导率，m3； Ti和Tj分别为 网格i和网格j与网格交界面之间的半传导率，m3o上式即为DFM模型中常采用的两点流动近似(TPFA)方法 该方法分开计算不同类型介质(基质、裂缝)的 传导率，更易于观察不同类型介质中流动机理的

15、不同。单个网格与网格交界面之间半传导率的计算公式为：Ti=AijKidTi=AijKid (6),mo式中：Aij为网格i与j的网格交界面面积，m2； Ki为网格 i渗透率，m2； d为网格交界面到网格中心点(基质网格为三角形形心，裂缝网格为四边形对角线交点0)垂直距离 对于裂缝网格: d=wf2d=wf2 (7)如图5所示，DFM模型的网格类型有两种，分别为基质网格(m)和裂缝网格(f),对应着两种网格半传导率Tm和Tf。网格的邻接类型有3种，分别为基质-基质()、基质-裂 缝（）和裂缝-裂缝（），对应着3种全传导率Tm-m, Tm-f 和Tf-fo对裂缝形变进行模拟的核心方法就是利用裂缝形

16、变 数学模型对裂缝半传导率进行修正，进而修正与裂缝网格相 关的全传导率，到达裂缝形变与流体流动耦合模拟的目的。图5模拟器中2种网格和3种传导率示意图.典型非结构化离散裂缝网格模型；b.离散裂缝模型中的网格连接和传导率（，为网格交界面类型编号。）Fig. 5Schematicsketchoftwokindsofgridsandthreekindsoftransmissibilityinthesimulator2. 2裂缝形变数值表征将裂缝形变模型与裂缝网格半传导率耦合，即将公式(1)代入公式(6)： Tf (p)=Tf01 e res(o o 2) Tf (p) =Tf011 - B(onon+

17、中)(0 o n o 2) Tf (p) =Tf 0(。nO) Tfp=Tf01 0 res。n三。2Tfp二TfOll-B onon+e Oono 2Tfp=TfO。nWO (8) 式中：TfO为裂缝网格初始半传导率，川3,计算式为:TfO=2Ai,fKfwmaxTfO=2Ai, fKfwmax (9) 式中：Ai, f表示网格i与裂缝网格f交界面面积，m2； Kf为裂缝网格渗透率，m2o考虑裂缝形变效应的裂缝半传导率 具有时变特征，公式(8)中on与裂缝中的流体压力密切 相关。考虑裂缝变形效应的跨尺度耦合流动全传导率计算仍 基于TPFA方法公式(5),此时，与裂缝网格相关的全传 导率(基质

18、-裂缝，裂缝-裂缝)也具备时变特征。即，基质 网格(m)-裂缝网格(f)流动全传导率计算公式为：Tm- f(p)=TmTf(p)Tm+Tf(p)Tm-fp=TmTfpTm+Tfp (10)裂缝网格(fl)-裂缝网格(f2)流动全传导率计算公式为：Tfl- f2(p)=Tfl(p)Tf2(p)Tf1(p)+Tf2(p)Tfl-f2p=TflpTf2 pTflp+Tf2p (11)在模拟计算过程中，首先计算出只与网格几何参数有关的裂 缝网格初始半传导率矩阵，而后在每一步模拟过程中对半传 导率和全传导率进行修正。3动态缝网模拟实例应用3. 1典型模型动态缝网模拟建立致密油藏中一 口多段压裂水平井典型

19、模型，如图6所示,以直观研究裂缝形变效应对油藏开发的影响。典型模型的尺 寸为400mX200m,水平井长度设计为300m,压裂段数为7段，半缝长为60mo 图6图6用于裂缝形变效应模拟的多段压裂水平井网格模型Fig. 6Meshmodelcontainingmulti- stagefracturesforthe simulationoffracturedeformation表1总结了涉及到的模拟参数，模拟时间为365d。裂缝应力 参数设置为：剩余裂缝宽度系数9res=0. 1,极限有效应力 。2=12MPa ,油藏总正应力在各网格间均匀分布 。normal=20MPao为模拟水力压裂缝的缝宽在

20、长度方向上的 非均匀分布现象，对每一条裂缝的不同位置赋了不同的缝宽 初值，每条裂缝的宽度呈现“两头窄中间宽”的分布规律, 缝宽范围为110mmo表1多段压裂水平井模型数值模拟参数TableISimulationparametersofthemulti- stagefracture dwellmodel参数取值参数取值油藏压力/MPa20原油压缩系数 /MPa-1 IX 10-3基质渗透率/ (10-3 um2) 0.01地层原油粘 度/ (mPa s) 5基质孔隙度/%8地层原油密度/ (kg m-3) 850裂缝渗透率/ (10-3uni2) 5原油体积系数1. 13岩石压 缩系数/MPa-

21、11X10-5生产时定井底流压/MPalO模型中设置的裂缝形变程度分为无裂缝形变、弱裂缝形变 (里二10MPa)和强裂缝形变(中=2MPa)。利用模拟器计算得 到的水平井产量和累计产量变化曲线如图7所示。由图可知, 裂缝形变显著降低了油井初产，其后一直保持低产水平，主 要因为在衰竭式开采条件下，开井初期裂缝中的压力快速下 降，导致裂缝快速闭合，储层供给能力急剧降低，从而极大 降低了水平井产能，如图8所示不同裂缝刚度对应的不同模 拟时间点压力分布图。图8说明，考虑裂缝形变效应后，裂 缝的导流能力显著降低，井周围的压降幅度明显缩小。此外, 比照可看出，中二2MPa时，裂缝导流能力的减弱首先发生在

22、裂缝的两端，且减弱程度最大。当裂缝中间处的导流能力处 于相对较高的水平时，裂缝两端的压降程度仍然很小，这主 要是因为裂缝两端的宽度小，因此在压力下降过程中优先闭 合。图7图7不同裂缝刚度对应的水平井生产动态曲线Fig. 7Productiondynamiccurvesobtainedunderdifferentf racturestiffness图8图8不同裂缝刚度对应的压力分布变化Fig. 8Reservoirpressurevariationunderdifferentfractu restiffness图9为不同时期缝网裂缝宽度分布图，即动态缝网模型。可 看出，裂缝形变的过程由裂缝两端向中间延伸。在生产初期, 裂缝形态变化较剧烈，闭合速度较快；到了生产中后期，闭 合速度变慢；模拟期末，裂缝两端的缝宽基本接近极限宽度, 这与图1所示的裂缝宽度的双曲变化规律相吻合。图9