CN115659736B - 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 - Google Patents
一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 Download PDFInfo
- Publication number
- CN115659736B CN115659736B CN202211280274.1A CN202211280274A CN115659736B CN 115659736 B CN115659736 B CN 115659736B CN 202211280274 A CN202211280274 A CN 202211280274A CN 115659736 B CN115659736 B CN 115659736B
- Authority
- CN
- China
- Prior art keywords
- fracture
- fracturing
- stress
- hydraulic fracture
- hydraulic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 96
- 230000008569 process Effects 0.000 title claims abstract description 53
- 230000035699 permeability Effects 0.000 claims abstract description 45
- 230000008859 change Effects 0.000 claims abstract description 19
- 238000012512 characterization method Methods 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 10
- 238000010008 shearing Methods 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims description 38
- 239000011435 rock Substances 0.000 claims description 36
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000006073 displacement reaction Methods 0.000 claims description 23
- 238000011065 in-situ storage Methods 0.000 claims description 15
- 239000003795 chemical substances by application Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000004913 activation Effects 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 8
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 230000000704 physical effect Effects 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000001939 inductive effect Effects 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000000877 morphologic effect Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000007704 transition Effects 0.000 abstract description 3
- 239000000243 solution Substances 0.000 description 10
- 230000008878 coupling Effects 0.000 description 8
- 238000010168 coupling process Methods 0.000 description 8
- 238000005859 coupling reaction Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 5
- 230000000903 blocking effect Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001976 improved effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，包括计算注入暂堵剂之前水力裂缝延伸路径和几何尺寸；计算注入暂堵剂之后水力裂缝内压力变化；计算压裂过程中地层应力场变化、储层压力场变化和储层渗透率场变化；计算深层页岩水平井缝内转向压裂过程中张性缝网扩展区体积与剪切缝网扩展区体积；输出压裂过程中压裂缝网扩展区体积数据，绘制压裂过程中水力裂缝延伸半长、水力裂缝缝内压力、缝网扩展区体积随时间变化曲线，以及压裂结束时水力裂缝延伸形态图与压裂缝网扩展区展布图。本发明解决了深层页岩气水平井通常采用缝内转向压裂工艺提高缝网建造效率，但难以对缝内转向压裂缝网体积进行准确表征的问题。
Description
技术领域
本发明涉及一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，属于页岩气开发技术领域。
背景技术
随着国内外油气田勘探开发的不断深入，以页岩气为代表的非常规资源将接替未来油气资源。近年来，通过实施水平井分段多簇缝网压裂，我国中浅层页岩气开发已经取得了阶段性成果。然而，我国有超过60％的页岩气资源埋存于3500m以上的深层，页岩气开发正逐步向潜力巨大的深层进军。与中浅层页岩相比，深层页岩脆性较弱，水平地应力差较大，天然裂缝发育情况更为复杂，复杂缝网建造难度更大。为此，深层页岩气水平井通常采用缝内转向压裂工艺提高缝网建造效率。
目前，缝内转向压裂工艺已广泛应用于重庆涪陵与川南盆地等国内主要深层页岩气产区(习传学等人2018；曹学军等人，2019；沈骋等人，2022)。曲宝龙(2019)、杨恒林等人(2022)通过真三轴大物模实验，结合三维扫描技术，发现缝内转向压裂过程中，暂堵剂对水力裂缝封堵后，缝内压力有一定增加，可以促进天然裂缝激活，形成复杂缝网。刘建升等人(2016)、董志刚等人(2017)在缝内转向压裂矿场施工过程中，通过微地震监测，证实了缝内转向工艺对深层页岩气压裂缝网扩展有一定促进作用。
综上所述，目前所有深层页岩气水平井缝内转向压裂工艺研究中，主要通过室内实验与现场微地震监测等方法，对缝内转向压裂效果进行定性评价，缺乏对缝内转向压裂缝网建造情况进行定量表征。因此，亟需建立一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，充分考虑深层页岩气水平井缝内转向压裂过程中，暂堵后水力裂缝内压力上升，促进天然裂缝破坏激活，提高缝网复杂度的力学机制，进而对缝内转向压裂缝网扩展行为进行准确模拟计算，为深层页岩气水平井缝内转向压裂工艺参数优化设计奠定理论基础。
发明内容
为了克服现有技术中的深深层页岩气水平井通常采用缝内转向压裂工艺提高缝网建造效率，但难以对缝内转向压裂缝网体积进行准确表征的问题，本发明提供一种深层页岩气水平井缝内转向压裂缝网扩展计算方法。
本发明解决上述技术问题所提供的技术方案是：一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，包括以下步骤：
根据深层页岩气目标储层地质工程情况，获取储层物性参数、三轴地应力参数、岩石力学参数、天然裂缝参数以及压裂工程参数；
根据裂缝延伸理论，建立深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型，利用岩石力学参数、三轴地应力参数和压裂工程参数计算深层页岩水平井缝内转向压裂过程中注入暂堵剂之前水力裂缝延伸路径和几何尺寸；
根据流体力学理论，建立深层页岩气水平井缝内转向压裂暂堵后水力裂缝增压模型，利用岩石力学参数和压裂工程参数计算深层页岩水平井缝内转向压裂过程中注入暂堵剂之后水力裂缝内压力变化；
根据岩石力学理论，建立深层页岩气水平井缝内转向压裂地层应力模型，利用三轴地应力参数和岩石力学参数计算深层页岩水平井压裂过程中地层应力场变化；
根据渗流力学理论，建立深层页岩气水平井缝内转向压裂储层压力模型，利用储层物性参数计算深层页岩水平井压裂过程中储层压力场变化；
根据断裂力学理论，建立深层页岩气水平井缝内转向压裂储层渗透率模型，利用天然裂缝参数和三轴地应力参数计算深层页岩水平井压裂过程中储层渗透率场变化；
根据空间数值积分方法，建立深层页岩气水平井缝内转向压裂缝网扩展表征模型，利用天然裂缝破坏点坐标数据，计算深层页岩水平井缝内转向压裂过程中张性缝网扩展区体积与剪切缝网扩展区体积，并将两者空间并集计为缝网扩展区体积；
根据深层页岩气水平井缝内转向压裂缝网扩展计算结果，输出压裂过程中压裂缝网扩展区体积数据，绘制压裂过程中水力裂缝延伸半长随时间变化曲线、压裂过程中水力裂缝缝内压力随时间变化曲线、压裂过程中缝网扩展区体积随时间变化曲线、压裂结束时水力裂缝延伸形态图、压裂结束时压裂缝网扩展区展布图。
进一步的技术方案是，所述储层物性参数包括基质孔隙度、基质渗透率、原始孔隙压力；所述三轴地应力参数包括水平最小主应力、水平最大主应力、垂向应力；所述岩石力学参数包括杨氏模量、泊松比、断裂韧性；所述天然裂缝参数包括天然裂缝倾角、逼近角、内聚力、摩擦系数、抗张强度、切向刚度、法向刚度；所述压裂工程参数包括压裂排量、压裂液量、压裂液粘度、压裂液密度、射孔簇数、射孔簇间距、射孔孔眼数量、射孔孔眼直径、暂堵剂注入时间。
进一步的技术方案是，所述深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型包括：
水力裂缝物质平衡方程：
式中：qf为水力裂缝内流量，m3/s；hf为水力裂缝高度，m；wf为水力裂缝宽度，m；s为水力裂缝长度方向坐标，m；t为时间，s；vleakoff为压裂液滤失速度，m/s；
水力裂缝压降方程：
式中：pf为水力裂缝缝内压力，Pa；μ为压裂液黏度，Pa·s；q为裂缝内流量，m3/s；
水力裂缝偏转方程：
式中：KⅠ、KⅡ为水力裂缝尖端单元第一、第二类应力强度因子，Pa·m1/2；E为岩石杨氏模量，Pa；Dn、Ds为水力裂缝尖端单元法向、切向应变，m；ν为岩石泊松比，无量纲；a为裂缝尖端单元长度的一半，m；θHF为水力裂缝延伸转向角，°；
水力裂缝闭合应力方程：
σc＝σhcos2θHF+σHsin2θHF
式中：σH为最大水平主应力，Pa；σh为最小水平主应力，Pa，σc为水力裂缝闭合应力，Pa；
水力裂缝高度方程：
式中：KIC为岩石断裂韧性，Pa/m0.5；
水力裂缝流量分配方程：
式中：pheel为水平井跟端压力，Pa；pfi,i为第i条裂缝首个单元内压力，Pa；Δppf,i为第i条裂缝射孔孔眼处的摩阻压降，Pa；Δpw,j为第j水平井段沿程压降，Pa；npf为射孔孔眼数量，个；dpf为射孔孔眼直径，m；α为孔眼流量系数，取0.8～0.85，无量纲；ρ为压裂液密度，kg/m3；μ为压裂液黏度，Pa·s；Lw,j为第j段水平井长度，m；qw,j为第j段水平井流量，m3/s；qT为压裂液总流量，m3/s；qi为第i条裂缝内总流量，m3/s；dw为水平井筒直径，m；i表示裂缝编号；j表示水平井段编号；
水力裂缝延伸初始和边界条件方程：
式中：Lf为水力裂缝长度，m。
进一步的技术方案是，所述深层页岩气水平井缝内转向压裂暂堵后水力裂缝增压模型包括：
暂堵后水力裂缝宽度增值方程：
式中：ΔWf为暂堵后水力裂缝宽度变化，m；qf为暂堵后水力裂缝流量，m3/s；Tdivert为暂堵后裂缝延伸总时间，m；Hf,divert为暂堵时水力裂缝高度，m；Lf,divert为暂堵时水力裂缝长度，m；Vleakoff为裂缝内液体滤失体积，m3；s为裂缝长度坐标，m；h为裂缝高度坐标，m；t为时间，s；
暂堵后水力裂缝压力增值方程：
式中：Δpf为暂堵后水力裂缝压力变化，Pa；E为岩石杨氏模量，Pa；ν为岩石泊松比，无量纲；
暂堵后水力裂缝宽度方程：
Wf(s,t0+t)＝Wf(s,t0)+ΔWf
式中：t0为暂堵时间，s；
暂堵后水力裂缝压力方程：
pf(s,t0+t)＝pf(s,t0)+Δpf
式中：Δpf为暂堵后水力裂缝压力变化，Pa；t0为暂堵时间，s。
进一步的技术方案是，所述深层页岩气水平井缝内转向压裂地层应力模型包括：
水力裂缝离散单元应力-应变平衡方程：
式中：(Sf)j为水力裂缝离散单元j的切向位移，m；(Wf)j为水力裂缝离散单元j的法向位移，即裂缝宽度，m；pf为水力裂缝内压力，Pa；σc为水力裂缝壁面闭合应力，Pa；(Att)ij、(Ant)ij、(Atn)ij、(Ann)ij为水力裂缝离散单元j的切向位移和法向位移分别在水力裂缝离散单元i上引起的切向应力分量和法相应力分量，Pa/m；N为水力裂缝离散单元总数，无量纲；
水力裂缝诱导应力方程：
Δσzz＝ν(Δσxx+Δσyy)
式中：Δσxx、Δσyy、Δσzz、Δσxy为地层诱导应力分量，Pa；n为水力裂缝离散单元全局坐标z轴与局部坐标ζ轴夹角余弦值，无量纲；l为水力裂缝离散单元全局坐标x轴与局部坐标ξ轴夹角余弦值，无量纲；ζ、ξ为水力裂缝离散单元局部坐标，m；Fk为Papkovitch函数偏导方程，k∈{3～6}；
水力压裂地层应力方程：
式中：为原始地应力值分量，Pa；σxx、σyy、σzz、σxy、σyz、σxz为水力压裂过程中地层应力分量，Pa；
水力压裂三轴地应力方程：
式中：V为地应力张量矩阵特征值矩阵，Pa；D为地应力张量矩阵特征向量矩阵，m；eig为特征值与特征向量算子；σ为地应力张量矩阵，m；σh、σH、σv为三轴地应力(最小水平主应力、最大水平主应力、垂向应力)，Pa；x1、y1、z1、x2、y2、z2、x3、y3、z3为三轴地应力方向向量分量值，m。
进一步的技术方案是，所述深层页岩气水平井缝内转向压裂储层压力模型包括：
储层内压裂液流动连续性方程：
式中：K为渗透率张量，D；μ为压裂液粘度，Pa·s；qsc为在地面标准情况下，储层单位体积内的压裂液流入点源流量，s-1；Ct为地层综合压缩系数，Pa-1；pr为储层压力，Pa；t为时间，s；
水力裂缝-储层压裂液滤失方程：
式中：vleakoff为压裂液滤失速度，m/s；nfs为垂直于裂缝壁面的单位向量。
进一步的技术方案是，所述深层页岩气水平井缝内转向压裂储层渗透率模型包括：
天然裂缝壁面受力方程：
其中：
式中：pn为天然裂缝壁面所受正应力值，Pa；pτ为天然裂缝壁面所受切应力值，Pa；σij为应力张量分量，Pa；ni为天然裂缝单位法线向量分量，无量纲；θ为天然裂缝逼近角，即与最大水平主应力方向夹角，°；为天然裂缝倾角，°；i、j、k为坐标指标，取值x、y、z；
天然裂缝张性破坏方程：
pnf＞pn+St
式中：pnf为天然裂缝内流体压力，即储层压力，Pa；St为天然裂缝抗张强度，Pa；
天然裂缝剪切破坏方程：
pτ＞τ0+Kf·(pn-pnf)
式中：Kf为天然裂缝摩擦系数，无量纲；τ0为天然裂缝内聚力，Pa；
天然裂缝破坏激活渗透率方程：
式中：knf为天然裂缝发生破坏激活后的渗透率，D；Kn为天然裂缝法向刚度，Pa/m；pnf为天然裂缝内的流体压力，Pa；pn为天然裂缝壁面所受正应力值，Pa；Ks为天然裂缝切向刚度，Pa/m；Δpτ为天然裂缝壁面所受切应力值，Pa；为天然裂缝剪切膨胀角，°；
天然裂缝激活储层渗透率方程：
其中：
式中：K为储层渗透率张量，D；Kij为储层渗透率张量分量，D；knf为天然裂缝发生破坏激活后的渗透率，D；为储层基质渗透率张量分量，D；i、j为坐标指标，取值x、y、z。
进一步的技术方案是，所述深层页岩气水平井缝内转向压裂缝网扩展表征模型包括：
缝网扩展体积表征方程：
式中：Vtensile-SRV、Vshear-SRV、Vtotal-SRV为张性缝网扩展区体积、剪切缝网扩展区体积、缝网扩展区体积，m3；ε为储层网格单元；εtensile、εshear为含有天然裂缝张性破坏点和剪切破坏点坐标的储层网格单元集合；Δx(ε)、Δy(ε)、Δz(ε)为ε网格单元x、y、z方向网格边长，m。
本发明具有以下有益效果：
与现有技术相比，该方法专门针对深层页岩气水平井缝内转向压裂缝网扩展力学行为，建立了深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型、暂堵后水力裂缝增压模型、地层应力模型、储层压力模型、储层渗透率模型、缝网扩展表征模型，从而提出了一种深层页岩气水平井缝内转向压裂缝网扩展计算方法；由于该方法充分考虑了暂堵后水力裂缝内压力上升，促进天然裂缝破坏激活，促进复杂缝网扩展的力学机制，可对缝内转向压裂缝网扩展行为进行准确模拟计算；解决了深层页岩气水平井通常采用缝内转向压裂工艺提高缝网建造效率，但难以对缝内转向压裂缝网体积进行准确表征的问题。
附图说明
图1是本发明方法的计算流程框图；
图2是实施例目标井压裂过程中水力裂缝延伸半长随时间变化曲线；
图3是实施例目标井压裂过程中水力裂缝缝内压力随时间变化曲线；
图4是实施例目标井压裂过程中缝网扩展区体积随时间变化曲线；
图5是实施例目标井压裂结束时水力裂缝延伸形态图；
图6是实施例目标井压裂结束时压裂缝网扩展区展布图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明的一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，包括以下步骤：
步骤1、根据深层页岩气目标储层地质工程情况，获取储层物性参数、三轴地应力参数、岩石力学参数、天然裂缝参数以及压裂工程参数；
其中所述储层物性参数包括基质孔隙度、基质渗透率、原始孔隙压力；所述三轴地应力参数包括水平最小主应力、水平最大主应力、垂向应力；所述岩石力学参数包括杨氏模量、泊松比、断裂韧性；所述天然裂缝参数包括天然裂缝倾角、逼近角、内聚力、摩擦系数、抗张强度、切向刚度、法向刚度；所述压裂工程参数包括压裂排量、压裂液量、压裂液粘度、压裂液密度、射孔簇数、射孔簇间距、射孔孔眼数量、射孔孔眼直径、暂堵剂注入时间；
步骤2、根据裂缝延伸理论，建立深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型；
水力裂缝物质平衡方程：
式中：qf为水力裂缝内流量，m3/s；hf为水力裂缝高度，m；wf为水力裂缝宽度，m；s为水力裂缝长度方向坐标，m；t为时间，s；vleakoff为压裂液滤失速度，m/s；
水力裂缝压降方程：
式中：pf为水力裂缝缝内压力，Pa；μ为压裂液黏度，Pa·s；q为裂缝内流量，m3/s；
水力裂缝偏转方程：
式中：KⅠ、KⅡ为水力裂缝尖端单元第一、第二类应力强度因子，Pa·m1/2；E为岩石杨氏模量，Pa；Dn、Ds为水力裂缝尖端单元法向、切向应变，m；ν为岩石泊松比，无量纲；a为裂缝尖端单元长度的一半，m；θHF为水力裂缝延伸转向角，°；
水力裂缝闭合应力方程：
σc＝σhcos2θHF+σHsin2θHF (4)
式中：σH为最大水平主应力，Pa；σh为最小水平主应力，Pa，σc为水力裂缝闭合应力，Pa；
水力裂缝高度方程：
式中：KIC为岩石断裂韧性，Pa/m0.5；
水力裂缝流量分配方程：
式中：pheel为水平井跟端压力，Pa；pfi,i为第i条裂缝首个单元内压力，Pa；Δppf,i为第i条裂缝射孔孔眼处的摩阻压降，Pa；Δpw,j为第j水平井段沿程压降，Pa；npf为射孔孔眼数量，个；dpf为射孔孔眼直径，m；α为孔眼流量系数，取0.8～0.85，无量纲；ρ为压裂液密度，kg/m3；μ为压裂液黏度，Pa·s；Lw,j为第j段水平井长度，m；qw,j为第j段水平井流量，m3/s；qT为压裂液总流量，m3/s；qi为第i条裂缝内总流量，m3/s；dw为水平井筒直径，m；i表示裂缝编号；j表示水平井段编号；
水力裂缝延伸初始和边界条件方程：
式中：Lf为水力裂缝长度，m；
步骤3、联立方程式(1)～(6)，结合初始和边界条件方程(7)，通过三点有限差分方法与牛顿迭代法，计算深层页岩水平井缝内转向压裂过程中注入暂堵剂之前水力裂缝延伸路径和几何尺寸；
具体包括水力裂缝的长度、高度、宽度、缝内压力、以及偏转角度；
步骤4、根据流体力学理论，建立深层页岩气水平井缝内转向压裂暂堵后水力裂缝增压模型；
暂堵后水力裂缝宽度增值方程：
式中：ΔWf为暂堵后水力裂缝宽度变化，m；qf为暂堵后水力裂缝流量，m3/s；Tdivert为暂堵后裂缝延伸总时间，m；Hf,divert为暂堵时水力裂缝高度，m；Lf,divert为暂堵时水力裂缝长度，m；Vleakoff为裂缝内液体滤失体积，m3；s为裂缝长度坐标，m；h为裂缝高度坐标，m；t为时间，s；
暂堵后水力裂缝压力增值方程：
式中：Δpf为暂堵后水力裂缝压力变化，Pa；E为岩石杨氏模量，Pa；ν为岩石泊松比，无量纲。
暂堵后水力裂缝宽度方程：
Wf(s,t0+t)＝Wf(s,t0)+ΔWf (10)
式中：t0为暂堵时间，s；
暂堵后水力裂缝压力方程：
pf(s,t0+t)＝pf(s,t0)+Δpf (11)
步骤5、利用方程式(8)、(9)，分别计算暂堵后水力裂缝宽度增值和压力增值，再利用方程式(10)、(11)，分别计算暂堵后水力裂缝宽度和压力；
步骤6、根据岩石力学理论，建立深层页岩气水平井缝内转向压裂地层应力模型；
水力裂缝离散单元应力-应变平衡方程：
式中：(Sf)j为水力裂缝离散单元j的切向位移，m；(Wf)j为水力裂缝离散单元j的法向位移，即裂缝宽度，m；pf为水力裂缝内压力，Pa；σc为水力裂缝壁面闭合应力，Pa；(Att)ij、(Ant)ij、(Atn)ij、(Ann)ij为水力裂缝离散单元j的切向位移和法向位移分别在水力裂缝离散单元i上引起的切向应力分量和法相应力分量，Pa/m；N为水力裂缝离散单元总数，无量纲；
水力裂缝诱导应力方程：
/>
Δσzz＝ν(Δσxx+Δσyy) (16)
式中：Δσxx、Δσyy、Δσzz、Δσxy为地层诱导应力分量，Pa；n为水力裂缝离散单元全局坐标z轴与局部坐标ζ轴夹角余弦值，无量纲；l为水力裂缝离散单元全局坐标x轴与局部坐标ξ轴夹角余弦值，无量纲；ζ、ξ为水力裂缝离散单元局部坐标，m；Fk为Papkovitch函数偏导方程，k∈{3～6}；
水力压裂地层应力方程：
式中：为原始地应力值分量，Pa；σxx、σyy、σzz、σxy、σyz、σxz为水力压裂过程中地层应力分量，Pa；
水力压裂三轴地应力方程：
式中：V为地应力张量矩阵特征值矩阵，Pa；D为地应力张量矩阵特征向量矩阵，m；eig为特征值与特征向量算子；σ为地应力张量矩阵，m；σh、σH、σv为三轴地应力(最小水平主应力、最大水平主应力、垂向应力)，Pa；x1、y1、z1、x2、y2、z2、x3、y3、z3为三轴地应力方向向量分量值，m；
步骤7、利用方程式(12)，计算水力裂缝各离散单元的切向位移和法向位移，即裂缝宽度，代入方程式(1)、(2)中耦合求解；利用方程式(13)～(16)，计算水力裂缝诱导应力分量；利用方程式(17)，计算水力压裂过程中地层应力分量；利用方程式(18)，计算水力压裂过程中三轴地应力，代入方程式(4)中耦合求解；
步骤8、根据渗流力学理论，建立深层页岩气水平井缝内转向压裂储层压力模型；
储层内压裂液流动连续性方程：
式中：K为渗透率张量，D；μ为压裂液粘度，Pa·s；qsc为在地面标准情况下，储层单位体积内的压裂液流入点源流量，s-1；Ct为地层综合压缩系数，Pa-1；p r 为储层压力，Pa；t为时间，s；
水力裂缝-储层压裂液滤失方程：
式中：vleakoff为压裂液滤失速度，m/s；nfs为垂直于裂缝壁面的单位向量；
步骤9、利用方程式(19)，计算压裂过程中储层压力；利用方程式(20)，计算压裂液滤失速度，代入方程式(1)、(8)中耦合求解；
步骤10、根据断裂力学理论，建立深层页岩气水平井缝内转向压裂储层渗透率模型；
天然裂缝壁面受力方程：
其中：
式中：pn为天然裂缝壁面所受正应力值，Pa；pτ为天然裂缝壁面所受切应力值，Pa；σij为应力张量分量，Pa；ni为天然裂缝单位法线向量分量，无量纲；θ为天然裂缝逼近角，即与最大水平主应力方向夹角，°；为天然裂缝倾角，°；i、j、k为坐标指标，取值x、y、z。/>
天然裂缝张性破坏方程：
pnf＞pn+St (23)
式中：pnf为天然裂缝内流体压力，即储层压力，Pa；St为天然裂缝抗张强度，Pa。
天然裂缝剪切破坏方程：
pτ＞τ0+Kf·(pn-pnf) (24)
式中：Kf为天然裂缝摩擦系数，无量纲；τ0为天然裂缝内聚力，Pa。
天然裂缝破坏激活渗透率方程：
式中：knf为天然裂缝发生破坏激活后的渗透率，D；Kn为天然裂缝法向刚度，Pa/m；pnf为天然裂缝内的流体压力，Pa；pn为天然裂缝壁面所受正应力值，Pa；Ks为天然裂缝切向刚度，Pa/m；Δpτ为天然裂缝壁面所受切应力值，Pa；为天然裂缝剪切膨胀角，°；
天然裂缝激活储层渗透率方程：
其中：
式中：K为储层渗透率张量，D；Kij为储层渗透率张量分量，D；knf为天然裂缝发生破坏激活后的渗透率，D；为储层基质渗透率张量分量，D；i、j为坐标指标，取值x、y、z；
步骤11、利用方程式(21)、(22)，计算然裂缝壁面所受正应力值和切应力值；利用方程式(23)、(24)，判断天然裂缝破坏状态；利用方程式(25)，计算天然裂缝发生破坏激活后的渗透率；利用方程式(26)、(27)，计算天然裂缝发激活后的储层渗透率，代入方程式(19)、(20)中耦合求解；
步骤12、根据空间数值积分方法，建立深层页岩气水平井缝内转向压裂缝网扩展表征模型；
缝网扩展体积表征方程：
式中：Vtensile-SRV、Vshear-SRV、Vtotal-SRV为张性缝网扩展区体积、剪切缝网扩展区体积、缝网扩展区体积，m3；ε为储层网格单元；εtensile、εshear为含有天然裂缝张性破坏点和剪切破坏点坐标的储层网格单元集合；Δx(ε)、Δy(ε)、Δz(ε)为ε网格单元x、y、z方向网格边长，m；
步骤13、利用方程式(28)，计算深层页岩水平井缝内转向压裂过程中张性缝网扩展区体积与剪切缝网扩展区体积，并将两者空间并集计为缝网扩展区体积；
步骤14、根据深层页岩气水平井缝内转向压裂缝网扩展计算结果，输出压裂过程中压裂缝网扩展区体积数据，绘制压裂过程中水力裂缝延伸半长随时间变化曲线、压裂过程中水力裂缝缝内压力随时间变化曲线、压裂过程中缝网扩展区体积随时间变化曲线、压裂结束时水力裂缝延伸形态图、压裂结束时压裂缝网扩展区展布图。
实施例1：
一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，包括以下步骤：
已知某深层页岩气区块储层物性参数、三轴地应力参数、岩石力学参数、天然裂缝参数、压裂工程参数如表1所示。
表1 某深层页岩气区块相关数据表
/>
首先，根据裂缝延伸理论，建立深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型，又结合公式计算实施例1中注入暂堵剂之前水力裂缝延伸路径和几何尺寸：
具体步骤为：①输入岩石力学参数、三轴地应力参数和压裂工程参数；②联立方程式(1)～(6)，结合初始和边界条件方程(7)，通过三点有限差分方法与牛顿迭代法，计算深层页岩水平井缝内转向压裂过程中注入暂堵剂之前水力裂缝延伸路径和几何尺寸；
随后，根据岩石力学理论，建立深层页岩气水平井缝内转向压裂地层应力模型，又结合公式计算实施例1中注入暂堵剂之后水力裂缝内压力变化：
具体步骤为：①输入岩石力学参数和压裂工程参数；②利用方程式(8)、(9)，分别计算暂堵后水力裂缝宽度增值和压力增值；③利用方程式(10)、(11)，分别计算暂堵后水力裂缝宽度和压力。
随后，根据岩石力学理论，建立深层页岩气水平井缝内转向压裂地层应力模型，又结合公式计算实施例1中地层应力场变化：
具体步骤为：①输入三轴地应力参数和岩石力学参数；②利用方程式(12)，计算水力裂缝各离散单元的切向位移和法向位移，即裂缝宽度，代入方程式(1)、(2)中耦合求解；③利用方程式(13)～(16)，计算水力裂缝诱导应力分量；④利用方程式(17)，计算水力压裂过程中地层应力分量；⑤利用方程式(18)，计算水力压裂过程中三轴地应力，代入方程式(4)中耦合求解；
随后，根据渗流力学理论，建立深层页岩气水平井缝内转向压裂储层压力模型，又结合以下公式计算实施例1中储层压力场变化：
具体步骤为：①输入储层物性参数；②利用方程式(19)，计算压裂过程中储层压力；③利用方程式(20)，计算压裂液滤失速度，代入方程式(1)、(8)中耦合求解。
随后，根据断裂力学理论，建立深层页岩气水平井缝内转向压裂储层渗透率模型，又结合公式计算实施例1中储层渗透率场变化：
具体步骤为：①输入天然裂缝参数和三轴地应力参数；②利用方程式(21)、(22)，计算然裂缝壁面所受正应力值和切应力值；③利用方程式(23)、(24)，判断天然裂缝破坏状态；④利用方程式(25)，计算天然裂缝发生破坏激活后的渗透率；⑤利用方程式(26)、(27)，计算天然裂缝发激活后的储层渗透率，代入方程式(19)、(20)中耦合求解。
随后，根据空间数值积分方法，建立深层页岩气水平井缝内转向压裂缝网扩展表征模型，又结合以下公式计算实施例1中张性缝网扩展区体积与剪切缝网扩展区体积，并将两者空间并集计为缝网扩展区体积：
具体步骤为：①提取天然裂缝破坏点坐标数据；②利用方程式(28)，计算深层页岩水平井缝内转向压裂过程中张性缝网扩展区体积与剪切缝网扩展区体积，并将两者空间并集计为缝网扩展区体积。
最后，运用如图1所示的本发明方法的数值计算流程框图开展实例计算，并根据计算结果，输出压裂过程中压裂缝网扩展区体积数据(如表1所示)，并绘制压裂过程中水力裂缝延伸半长随时间变化曲线(如图2所示)、压裂过程中水力裂缝缝内压力随时间变化曲线(如图3所示)、压裂过程中缝网扩展区体积随时间变化曲线(如图4所示)、压裂结束时水力裂缝延伸形态图(如图5所示)、压裂结束时压裂缝网扩展区展布图(如图6所示)。
表2 压裂过程中压裂缝网扩展区体积数据表
该方法专门针对深层页岩气水平井缝内转向压裂缝网扩展力学行为，建立了深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型、暂堵后水力裂缝增压模型、地层应力模型、储层压力模型、储层渗透率模型、缝网扩展表征模型，从而提出了一种深层页岩气水平井缝内转向压裂缝网扩展计算方法；由于该方法充分考虑了暂堵后水力裂缝内压力上升，促进天然裂缝破坏激活，促进复杂缝网扩展的力学机制，可对缝内转向压裂缝网扩展行为进行准确模拟计算；解决了深层页岩气水平井通常采用缝内转向压裂工艺提高缝网建造效率，但难以对缝内转向压裂缝网体积进行准确表征的问题。
以上所述，并非对本发明作任何形式上的限制，虽然本发明已通过上述实施例揭示，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。
Claims (2)
1.一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，其特征在于，包括以下步骤：
根据深层页岩气目标储层地质工程情况，获取储层物性参数、三轴地应力参数、岩石力学参数、天然裂缝参数以及压裂工程参数；
根据裂缝延伸理论，建立深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型，利用岩石力学参数、三轴地应力参数和压裂工程参数计算深层页岩水平井缝内转向压裂过程中注入暂堵剂之前水力裂缝延伸路径和几何尺寸；
根据流体力学理论，建立深层页岩气水平井缝内转向压裂暂堵后水力裂缝增压模型，利用岩石力学参数和压裂工程参数计算深层页岩水平井缝内转向压裂过程中注入暂堵剂之后水力裂缝内压力变化；
根据岩石力学理论，建立深层页岩气水平井缝内转向压裂地层应力模型，利用三轴地应力参数和岩石力学参数计算深层页岩水平井压裂过程中地层应力场变化；
根据渗流力学理论，建立深层页岩气水平井缝内转向压裂储层压力模型，利用储层物性参数计算深层页岩水平井压裂过程中储层压力场变化；
根据断裂力学理论，建立深层页岩气水平井缝内转向压裂储层渗透率模型，利用天然裂缝参数和三轴地应力参数计算深层页岩水平井压裂过程中储层渗透率场变化；
根据空间数值积分方法，建立深层页岩气水平井缝内转向压裂缝网扩展表征模型，利用天然裂缝破坏点坐标数据，计算深层页岩水平井缝内转向压裂过程中张性缝网扩展区体积与剪切缝网扩展区体积，并将两者空间并集计为缝网扩展区体积；
根据深层页岩气水平井缝内转向压裂缝网扩展计算结果，输出压裂过程中压裂缝网扩展区体积数据，绘制压裂过程中水力裂缝延伸半长随时间变化曲线、压裂过程中水力裂缝缝内压力随时间变化曲线、压裂过程中缝网扩展区体积随时间变化曲线、压裂结束时水力裂缝延伸形态图、压裂结束时压裂缝网扩展区展布图；
所述深层页岩气水平井缝内转向压裂暂堵前水力裂缝延伸模型包括：
水力裂缝物质平衡方程：
式中：qf为水力裂缝内流量，m3/s；hf为水力裂缝高度，m；wf为水力裂缝宽度，m；s为水力裂缝长度方向坐标，m；t为时间，s；vleakoff为压裂液滤失速度，m/s；
水力裂缝压降方程：
式中：pf为水力裂缝缝内压力，Pa；μ为压裂液黏度，Pa·s；qf为水力裂缝内流量，m3/s；
水力裂缝偏转方程：
式中：KⅠ、KⅡ为水力裂缝尖端单元第一、第二类应力强度因子，Pa·m1/2；E为岩石杨氏模量，Pa；Dn、Ds为水力裂缝尖端单元法向、切向应变，m；ν为岩石泊松比，无量纲；a为裂缝尖端单元长度的一半，m；θHF为水力裂缝延伸转向角，°；
水力裂缝闭合应力方程：
σc＝σh cos2θHF+σH sin2θHF
式中：σH为最大水平主应力，Pa；σh为最小水平主应力，Pa，σc为水力裂缝闭合应力，Pa；
水力裂缝高度方程：
式中：KIC为岩石断裂韧性，Pa/m0.5；
水力裂缝流量分配方程：
式中：pheel为水平井跟端压力，Pa；pfi,i为第i条裂缝首个单元内压力，Pa；Δppf,i为第i条裂缝射孔孔眼处的摩阻压降，Pa；Δpw,j为第j水平井段沿程压降，Pa；npf为射孔孔眼数量，个；dpf为射孔孔眼直径，m；α为孔眼流量系数，取0.8～0.85，无量纲；ρ为压裂液密度，kg/m3；μ为压裂液黏度，Pa·s；Lw,j为第j段水平井长度，m；qw,j为第j段水平井流量，m3/s；qT为压裂液总流量，m3/s；qi为第i条裂缝内总流量，m3/s；dw为水平井筒直径，m；i表示裂缝编号；j表示水平井段编号；
水力裂缝延伸初始和边界条件方程：
式中：Lf为水力裂缝长度，m；
所述深层页岩气水平井缝内转向压裂暂堵后水力裂缝增压模型包括：
暂堵后水力裂缝宽度增值方程：
式中：ΔWf为暂堵后水力裂缝宽度变化，m；qf为暂堵后水力裂缝流量，m3/s；Tdivert为暂堵后裂缝延伸总时间，m；Hf,divert为暂堵时水力裂缝高度，m；Lf,divert为暂堵时水力裂缝长度，m；Vleakoff为裂缝内液体滤失体积，m3；s为裂缝长度坐标，m；h为裂缝高度坐标，m；t为时间，s；
暂堵后水力裂缝压力增值方程：
式中：Δpf为暂堵后水力裂缝压力变化，Pa；E为岩石杨氏模量，Pa；ν为岩石泊松比，无量纲；
暂堵后水力裂缝宽度方程：
Wf(s,t0+t)＝Wf(s,t0)+ΔWf
式中：t0为暂堵时间，s；
暂堵后水力裂缝压力方程：
pf(s,t0+t)＝pf(s,t0)+Δpf
式中：Δpf为暂堵后水力裂缝压力变化，Pa；t0为暂堵时间，s；
所述深层页岩气水平井缝内转向压裂地层应力模型包括：
水力裂缝离散单元应力-应变平衡方程：
式中：(Sf)j为水力裂缝离散单元j的切向位移，m；(Wf)j为水力裂缝离散单元j的法向位移，即裂缝宽度，m；pf为水力裂缝内压力，Pa；σc为水力裂缝壁面闭合应力，Pa；(Att)ij、(Ant)ij、(Atn)ij、(Ann)ij为水力裂缝离散单元j的切向位移和法向位移分别在水力裂缝离散单元i上引起的切向应力分量和法相应力分量，Pa/m；N为水力裂缝离散单元总数，无量纲；
水力裂缝诱导应力方程：
Δσzz＝ν(Δσxx+Δσyy)
式中：Δσxx、Δσyy、Δσzz、Δσxy为地层诱导应力分量，Pa；n为水力裂缝离散单元全局坐标z轴与局部坐标ζ轴夹角余弦值，无量纲；l为水力裂缝离散单元全局坐标x轴与局部坐标ξ轴夹角余弦值，无量纲；ζ、ξ为水力裂缝离散单元局部坐标，m；Fk为Papkovitch函数偏导方程，k∈{3～6}；
水力压裂地层应力方程：
式中：为原始地应力值分量，Pa；σxx、σyy、σzz、σxy、σyz、σxz为水力压裂过程中地层应力分量，Pa；
水力压裂三轴地应力方程：
式中：V为地应力张量矩阵特征值矩阵，Pa；D为地应力张量矩阵特征向量矩阵，m；eig为特征值与特征向量算子；σ为地应力张量矩阵，m；σh、σH、σv为三轴地应力(最小水平主应力、最大水平主应力、垂向应力)，Pa；x1、y1、z1、x2、y2、z2、x3、y3、z3为三轴地应力方向向量分量值，m；
所述深层页岩气水平井缝内转向压裂储层压力模型包括：
储层内压裂液流动连续性方程：
式中：K为渗透率张量，D；μ为压裂液粘度，Pa·s；qsc为在地面标准情况下，储层单位体积内的压裂液流入点源流量，s-1；Ct为地层综合压缩系数，Pa-1；p r 为储层压力，Pa；t为时间，s；
水力裂缝-储层压裂液滤失方程：
式中：vleakoff为压裂液滤失速度，m/s；nfs为垂直于裂缝壁面的单位向量；
所述深层页岩气水平井缝内转向压裂储层渗透率模型包括：
天然裂缝壁面受力方程：
其中：
式中：pn为天然裂缝壁面所受正应力值，Pa；pτ为天然裂缝壁面所受切应力值，Pa；σij为应力张量分量，Pa；ni为天然裂缝单位法线向量分量，无量纲；θ为天然裂缝逼近角，即与最大水平主应力方向夹角，°；为天然裂缝倾角，°；i、j、k为坐标指标，取值x、y、z；
天然裂缝张性破坏方程：
pnf>pn+St
式中：pnf为天然裂缝内流体压力，即储层压力，Pa；St为天然裂缝抗张强度，Pa；
天然裂缝剪切破坏方程：
pτ>τ0+Kf·(pn-pnf)
式中：Kf为天然裂缝摩擦系数，无量纲；τ0为天然裂缝内聚力，Pa；
天然裂缝破坏激活渗透率方程：
式中：knf为天然裂缝发生破坏激活后的渗透率，D；Kn为天然裂缝法向刚度，Pa/m；pnf为天然裂缝内的流体压力，Pa；pn为天然裂缝壁面所受正应力值，Pa；Ks为天然裂缝切向刚度，Pa/m；Δpτ为天然裂缝壁面所受切应力值，Pa；为天然裂缝剪切膨胀角，°；
天然裂缝激活储层渗透率方程：
其中：
式中：K为储层渗透率张量，D；Kij为储层渗透率张量分量，D；knf为天然裂缝发生破坏激活后的渗透率，D；为储层基质渗透率张量分量，D；i、j为坐标指标，取值x、y、z；
所述深层页岩气水平井缝内转向压裂缝网扩展表征模型包括：
缝网扩展体积表征方程：
式中：Vtensile-SRV、Vshear-SRV、Vtotal-SRV为张性缝网扩展区体积、剪切缝网扩展区体积、缝网扩展区体积，m3；ε为储层网格单元；εtensile、εshear为含有天然裂缝张性破坏点和剪切破坏点坐标的储层网格单元集合；Δx(ε)、Δy(ε)、Δz(ε)为ε网格单元x、y、z方向网格边长，m。
2.根据权利要求1所述的一种深层页岩气水平井缝内转向压裂缝网扩展计算方法，其特征在于，所述储层物性参数包括基质孔隙度、基质渗透率、原始孔隙压力；所述三轴地应力参数包括水平最小主应力、水平最大主应力、垂向应力；所述岩石力学参数包括杨氏模量、泊松比、断裂韧性；所述天然裂缝参数包括天然裂缝倾角、逼近角、内聚力、摩擦系数、抗张强度、切向刚度、法向刚度；所述压裂工程参数包括压裂排量、压裂液量、压裂液粘度、压裂液密度、射孔簇数、射孔簇间距、射孔孔眼数量、射孔孔眼直径、暂堵剂注入时间。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211280274.1A CN115659736B (zh) | 2022-10-19 | 2022-10-19 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211280274.1A CN115659736B (zh) | 2022-10-19 | 2022-10-19 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115659736A CN115659736A (zh) | 2023-01-31 |
| CN115659736B true CN115659736B (zh) | 2023-11-03 |
Family
ID=84989353
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211280274.1A Active CN115659736B (zh) | 2022-10-19 | 2022-10-19 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115659736B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117408168B (zh) * | 2023-12-15 | 2024-03-15 | 中国石油大学(华东) | 一种融合物理条件约束的压裂裂缝扩展快速预测方法 |
| CN117633409B (zh) * | 2024-01-25 | 2024-04-09 | 中国科学院地质与地球物理研究所 | 计算页岩油气储层压裂缝网渗流参数的方法、系统及设备 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110134984A (zh) * | 2019-03-06 | 2019-08-16 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩压裂过程中复杂裂缝扩展影响因素的分析方法 |
| CN110469303A (zh) * | 2019-07-04 | 2019-11-19 | 西南石油大学 | 一种基于四类改造体积的体积压裂参数优化设计方法 |
| US10689972B1 (en) * | 2019-05-31 | 2020-06-23 | Southwest Petroleum University | Automatic diagnosis method for wellhead pressure curve of hydraulic fracturing in shale gas horizontal well |
| CN111322050A (zh) * | 2020-04-24 | 2020-06-23 | 西南石油大学 | 一种页岩水平井段内密切割暂堵压裂施工优化方法 |
| CN112878982A (zh) * | 2020-12-31 | 2021-06-01 | 西南石油大学 | 一种考虑裂缝长期导流能力的深层页岩气产能预测方法 |
| CN112949129A (zh) * | 2021-03-03 | 2021-06-11 | 西南石油大学 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
| CN113971378A (zh) * | 2021-10-27 | 2022-01-25 | 西南石油大学 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
| WO2022252591A1 (zh) * | 2021-05-31 | 2022-12-08 | 中国矿业大学 | 一种水力压裂与甲烷原位燃爆协同致裂增透方法 |
-
2022
- 2022-10-19 CN CN202211280274.1A patent/CN115659736B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110134984A (zh) * | 2019-03-06 | 2019-08-16 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩压裂过程中复杂裂缝扩展影响因素的分析方法 |
| US10689972B1 (en) * | 2019-05-31 | 2020-06-23 | Southwest Petroleum University | Automatic diagnosis method for wellhead pressure curve of hydraulic fracturing in shale gas horizontal well |
| CN110469303A (zh) * | 2019-07-04 | 2019-11-19 | 西南石油大学 | 一种基于四类改造体积的体积压裂参数优化设计方法 |
| CN111322050A (zh) * | 2020-04-24 | 2020-06-23 | 西南石油大学 | 一种页岩水平井段内密切割暂堵压裂施工优化方法 |
| CN112878982A (zh) * | 2020-12-31 | 2021-06-01 | 西南石油大学 | 一种考虑裂缝长期导流能力的深层页岩气产能预测方法 |
| CN112949129A (zh) * | 2021-03-03 | 2021-06-11 | 西南石油大学 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
| WO2022252591A1 (zh) * | 2021-05-31 | 2022-12-08 | 中国矿业大学 | 一种水力压裂与甲烷原位燃爆协同致裂增透方法 |
| CN113971378A (zh) * | 2021-10-27 | 2022-01-25 | 西南石油大学 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
Non-Patent Citations (1)
| Title |
|---|
| "深层超深层页岩气水平井缝口暂堵压裂的裂缝调控模拟";胡东风，等;《天然气工业》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115659736A (zh) | 2023-01-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115659736B (zh) | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 | |
| Jiang et al. | Experimental and numerical study on hydraulic fracture propagation in coalbed methane reservoir | |
| CN108825196B (zh) | 一种筛管大斜度井的暂堵体积压裂工艺方法 | |
| Tan et al. | Experiments and analysis on hydraulic sand fracturing by an improved true tri-axial cell | |
| CN108468538A (zh) | 一种页岩水压裂缝扩展预测方法 | |
| CN105019877A (zh) | 一种页岩水平井压裂破碎区体积的计算方法 | |
| CN105735960A (zh) | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 | |
| CN112682016B (zh) | 油气田薄互层储层穿层压裂参数的确定方法及压裂方法 | |
| Gong et al. | Excavation collapse of Hangzhou subway station in soft clay and numerical investigation based on orthogonal experiment method | |
| CN112253187B (zh) | 基于粘土基浆液超前注浆改性底板硬岩抑制采动破坏深度的方法 | |
| LI et al. | Research on random propagation method of hydraulic fracture based on zero-thickness cohesive element | |
| Lu et al. | Numerical analysis of complex fracture propagation under temporary plugging conditions in a naturally fractured reservoir | |
| CN113987965B (zh) | 一种暂堵转向裂缝的预测方法及装置 | |
| CN108612518A (zh) | 一种煤层气井径向微井眼钻孔、水力压裂参数的确定方法 | |
| Rachez et al. | 3D-hydromechanical behavior of a stimulated fractured rock mass | |
| Zhao et al. | Study on catastrophe mechanism of hydraulic fracturing fracks initiation in radial boreholes | |
| CN113107363B (zh) | 一种促进裂隙错动提升自支撑能力的水平井施工方法 | |
| Liu et al. | Numerical simulation study on the distribution law of deviatoric stress of floor under the influence of mining | |
| CN104712299B (zh) | 适合气井控水增气压裂的设计方法 | |
| Lv | Impact of hydraulic fracturing and borehole spacing on gas drainage along a coal seam | |
| Liu et al. | Redevelopment of fractures and permeability changes after multi-seam mining of shallow closely spaced coal seams | |
| Yang et al. | Visualization and characterization of experimental hydraulic fractures interacting with karst fracture-cavity distributions | |
| Dong et al. | Study on the dosage of temporary plugging Agent for steering fracturing | |
| CN115217456B (zh) | 一种提高页岩气水平井产量的渗流体积压裂方法 | |
| CN117328861B (zh) | 一种模拟水力裂缝穿越岩层的计算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |