CN105019877B - 一种页岩水平井压裂破碎区体积的计算方法 - Google Patents
一种页岩水平井压裂破碎区体积的计算方法 Download PDFInfo
- Publication number
- CN105019877B CN105019877B CN201510352975.5A CN201510352975A CN105019877B CN 105019877 B CN105019877 B CN 105019877B CN 201510352975 A CN201510352975 A CN 201510352975A CN 105019877 B CN105019877 B CN 105019877B
- Authority
- CN
- China
- Prior art keywords
- fracture
- mpa
- pressure
- formula
- rock
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种页岩水平井压裂破碎区体积计算方法,包括以下步骤:A、可以运用裂缝扩展理论建立页岩压裂过程中的水力裂缝延伸模型;B、可以运用岩石力学理论建立岩压裂过程中的应力场变化模型;C、可以运用渗流力学理论建立岩压裂过程中的压力场变化模型;D、可以运用、岩石破坏理论建立岩压裂过程中的天然裂缝破坏准则模型;E、结合上述4个步骤通过计算流程即可对压裂破碎区体积进行计算和表征。有益效果是解决了现有页岩水平井压裂时对破碎区体积计算过程费时费力费钱且导致压裂破碎区体积评价不能广泛推广应用问题。
Description
技术领域
本发明涉及页岩气开发领域,尤其是涉及一种页岩水平井压裂破碎区体积计算方法。
背景技术
目前,水平井缝网压裂技术是全球页岩气藏开发过程中广泛使用的增产措施,该技术能够显著提高页岩气开发的经济性和高效性。体积压裂的基本特点是采用大液量、大排量、大砂量、小粒径、低砂比的施工参数进行压裂施工,目的在于利用人工裂缝最大化地破坏、激活和沟通储层中的天然裂缝,形成较大的压裂破碎区体积,从而提高水平井附近区域内的有效渗透性。
对于页岩气藏水平井分段分簇的体积压裂前期设计和后期评估来说,准确表征压裂破碎区的形态和体积具有非常重要的指导性作用。目前,压裂破碎区评价方法主要有微地震监测法、倾斜仪测量法和理论模型计算法。其中,理论模型计算主要包括连续介质模型和离散裂缝网络模型。以上方法各有优劣,微地震监测法较为可靠,但成本高昂;倾斜仪测量法虽成本较低,但精度不理想,应用受限;连续介质模型理论简单,方法易用,但其研究主要集中于宏观层面,未能深入细观乃至微观层面;离散裂缝网络模型的理论基础较完备,研究尺度精细至细观层面,但建模复杂,求解所需已知参数较多,通常需要借助高级数值计算工具,如有限元(FEM)、扩展有限元(XFEM)、边界元(BEM)等,运算量较大,运算耗时较长,通常局限于二维模型计算,不太适合现场快速应用。此外,一部分理论模型没有将人工裂缝之间的应力干扰因素进行考虑。然而,页岩气藏水平井分段压裂时,簇间距通常为20~50m,各簇人工裂缝产生的诱导应力会显著改变原始地应力分布,从而影响压裂破碎区体积的形态。
因此,对页岩水平井压裂破碎区体积进行表征和计算非常有必要,面对目前各类已有方法的局限性,需要探索一种能够快速、准确、经济的页岩水平井压裂破碎区体积评价方法。
发明内容
本发明提供一种解决了现有页岩水平井压裂时对破碎区体积计算过程费时费力费钱且导致压裂破碎区体积评价不能广泛推广应用问题的页岩水平井压裂破碎区体积计算方法。
本发明的技术方案是:一种页岩水平井压裂破碎区体积计算方法:包括以下步骤:
A、可以运用裂缝扩展理论建立页岩压裂过程中的水力裂缝延伸模型;
B、可以运用岩石力学理论建立岩压裂过程中的应力场变化模型;
C、可以运用渗流力学理论建立岩压裂过程中的压力场变化模型;
D、可以运用岩石破坏理论建立岩压裂过程中的天然裂缝破坏准则模型;
E、结合上述4个步骤通过计算流程即可对压裂破碎区体积进行计算和表征。
进一步的是,所述步骤A中水力裂缝延伸模型包括以下计算公式:
裂缝缝宽方程:
缝内流体流动方程:
物质平衡方程:
求解条件:
W(L,t)=0在y≥=L(t) (4)
式中:W(y)—延伸方向上y处的最大缝宽,mm;L(t)—t时刻水力裂缝长度,m;E—岩石杨氏模量,MPa;p(y)—缝内y处压力,MPa;σn(y)—延伸方向上y处裂缝面受到的正应力,MPa;hf—裂缝缝高,m;ν—岩石泊松比,无因次;μ—压裂液粘度,mPa·s;cL—滤失系数,m/s1/2;q—缝内流量,m3/s;Qi—第i条水力裂缝内压裂液分流量,m3/s;Qpumping—压裂施工泵注总排量,m3/s;M—水裂缝条数,即射孔簇数。
结合求解条件式(4)~式(6),采用隐式有限差分法对式(1)和式(3)进行求解,得到页岩压裂时各条水力裂缝的几何参数;
进一步的是,所述步骤B中应力场变化模型包括以下计算公式:
利用弹性力学理论,计算水力裂缝在附近区域产生的诱导应力场:
Δσyy=ν(Δσx+Δσz) (8)
各几何参数表达式如下:
式中:Δσxx、Δσyy、Δσzz—x、y、z方向上法向诱导应力,MPa;Δσxz—x面上z方向切向诱导应力,MPa;pnet—人工裂缝内净压力,MPa;ν—岩石泊松比,无量纲;c—裂缝高度hf的1/2,m;r、r1、r2—几何长度参数,m;θ、θ1、θ2—几何角度参数,°;
利用叠加原理计算出多条水力裂缝存在时的储层应力场:
式中:σxx、σyy、σzz、σxy、σyz、σxz—初始地应力值分量,MPa;ΣΔσxx、ΣΔσyy、ΣΔσzz、ΣΔσxz—多条裂缝产生的诱导应力值各分量之和,MPa;σ′xx、σ′yy、σ′zz、σ′xy、σ′yz、σ′xz—当前地应力值分量,MPa。
进一步的是,所述步骤C中压力场变化模型包括以下计算公式:
利用双重介质理论来计算缝网压裂,在计算Laplace域内各向异性页岩储层压力增量为:
上式各中间参数和无因次量表达式为:
u=sf(s) (27)
式中:域内压力场,MPa;L—参考长度,可取水平井长度,m;hf—裂缝高度,m;hr—油藏厚度,m;hrD—无因次油藏厚度,无量纲;kmx—基质***x方向上渗透率,D;kmy—基质***y方向上渗透率,D;kmz—基质***z方向上渗透率,D;km—基质***等效渗透率,D;—单位面积流量,m/s;Q—注入排量,m3/s;μ—液体粘度,mPa·s;s—Laplace变量;Lf—裂缝半长,m;K0—0阶Bessel函数;u—自定义函数;ω—弹性储能,无量纲;λ—无因次窜流系数,无量纲;α—岩块形状因子,无量纲;φ—孔隙度,无量纲;ct—总压缩系数,MPa-1;kfx—裂缝***x方向渗透率,D;kfy—裂缝***y方向渗透率,D;kfz—裂缝***z方向渗透率,D;kf—裂缝***等效渗透率,D;xw—井底x坐标,m;yw—井底y坐标,m;zw—井底z坐标,m;xD—无因次x坐标,无量纲;yD—无因次y坐标,无量纲;xwD—井底无因次x坐标,无量纲;ywD—井底无因次y坐标,无量纲;下标m表示基质***;下标m表示天然裂缝***;
求解出Laplace域中的储层压力场后,利用Stehfest数值反演即可得到储层实际压力场Δp(x,y,z,t):
其中:
式中:t—时间,s;N—求和项数,值越大越精确,一般取6~18之间的偶整数;
利用叠加原理计算出多条水力裂缝存在时的储层压力场:
p′(x,y,z,t)=pi(x,y,z,0)+ΣΔp(x,y,z,t) (33)
式中:pi—初始储层压力场,MPa;ΣΔp—储层压力增量叠加场,MPa;p′—当前储层压力场,即天然裂缝***压力场,MPa。
进一步的是,所述步骤D中天然裂缝破坏准则模型包括以下计算公式:
天然裂缝的张性破坏判别式为:
pnf>pn+St (34)
剪切破坏判别式为:
其中:
pn=σxxnxnx+σxynynx+σxznznx+σxynxny+σyynyny (36)
+σyznzny+σxznxnz+σyznynz+σzznznz
pτ=[(σxxnx+σxyny+σxznz)2
+(σxynx+σyyny+σyznz)2+(σxznx+σyzny+σzznz)2 (37)
-(σxxnxnx+σxynynx+σxznznx+σxynxny+σyynyny
+σyznzny+σxznxnz+σyznynz+σzznznz)2]1/2
式中:σij—储层应力张量分量,i,j指标取x,y,z,MPa;ni—裂缝法向矢量分量,i指标取x,y,z,MPa;α—天然裂缝与水平主应力方向夹角(逼近角),°;—天然裂缝与水平面夹角(倾角),°;pn—天然裂缝壁面所受正应力值,MPa;pτ—天然裂缝壁面所受切应力值,MPa;Kf—天然裂缝摩擦系数,无量纲;pnf—天然裂缝内压力,等于当前储层压力p′,MPa;St—天然裂缝抗张强度,MPa;τ0—天然裂缝内聚力,MPa。
进一步的是,所述步骤E中的计算流程包括以下步骤:
1)根据水平井筒射孔参数确定分段压裂的簇间距和簇数,建立相应的三维笛卡尔直角坐标系;
2)利用上述公式(1)、(3)、(4)~(6)计算出裂缝扩展空间几何尺寸;
3)利用上述公式(7)~(11)、(17)求解出当前地应力场σ′ij;
4)利用式(18)、(31)、(32)求解出当前储层压力场p′;
5)利用式(34)、(35),带入σ′ij和p′判断储层内任一点的天然裂缝是否发生破坏,以及破坏的类型;
6)根据天然裂缝破坏判断结果,利用空间数值积分方法,分别计算张性破碎区体积和剪切破碎区体积,并将两者的空间并集算作总体破碎区体积。
本发明的有益效果是,与现有技术相比,该方法考虑页岩水平井分段分簇压裂过程中流体扩散渗流场和裂缝诱导应力场同时改变对页岩体天然裂缝的触发破坏机制,建立了一种适合页岩水平井分段分簇压裂破碎区体积计算新模型。由于该模型的计算方法与页岩压裂破碎区实际的物理演化过程相一致,为此,可对压裂破碎区进行准确的定量评价,包括剪切破碎区体积和张性破碎区体积,并获得压裂破碎区体积的三维形态,实现页岩水平井压裂评价的经济化、快速化、准确化,解决了现有页岩水平井压裂时对破碎区体积计算过程费时费力费钱,导致压裂破碎区体积评价不能广泛推广应用的问题。
附图说明
图1是实施例本发明方法的计算流程图;
图2是水力裂缝延伸模型示意图;
图3是应力场变化模型示意图;
图4是压力场变化模型示意图;
图5是天然裂缝破坏准则模型示意图;
图6是实施例地层x-z截面上σx分布图;
图7是实施例地层x-z截面上σy分布图;
图8是实施例地层x-z截面上σz分布图;
图9是实施例储层压力分布图;
图10是实施例的总体压裂破碎区体积空间形态展布图。
具体实施方式
下面结合附图和实施例对本发明做进一步的详细说明,但不构成对发明的任何限制。
实施例1:
一种页岩水平井压裂破碎区体积计算方法,包括以下步骤:
已知某一典型页岩气藏地质参数和水平井压裂施工参数如下表:
表1实例1实际参数
首先运用裂缝扩展理论建立页岩压裂过程中的水力裂缝延伸模型且该模型的示意图如图2所示,又结合以下公式计算实施例1中三条水力裂缝延伸情况:
裂缝缝宽方程:
缝内流体流动方程:
物质平衡方程:
求解条件:
W(L,t)=0在y≥=L(t) (42)
式中:W(y)—延伸方向上y处的最大缝宽,mm;L(t)—t时刻水力裂缝长度,m;E—岩石杨氏模量,MPa;p(y)—缝内y处压力,MPa;σn(y)—延伸方向上y处裂缝面受到的正应力,MPa;hf—裂缝缝高,m;ν—岩石泊松比,无因次;μ—压裂液粘度,mPa·s;cL—滤失系数,m/s1/2;q—缝内流量,m3/s;Qi—第i条水力裂缝内压裂液分流量,m3/s;Qpumping—压裂施工泵注总排量,m3/s。
结合求解条件式(42)~式(44),采用隐式有限差分法对式(39)和式(41)进行求解,得到页岩压裂时各条水力裂缝的几何参数;
然后运用岩石力学理论建立岩压裂过程中的应力场变化模型且该模型的示意图如图3所示,结合以下公式计算单条水力裂缝在储层中产生的诱导应力场:
Δσyy=ν(Δσx+Δσz) (46)
各几何参数表达式如下:
式中:Δσxx、Δσyy、Δσzz—x、y、z方向上法向诱导应力,MPa;Δσxz—x面上z方向切向诱导应力,MPa;pnet—人工裂缝内净压力,MPa;ν—岩石泊松比,无量纲;c—裂缝高度hf的1/2,m;r、r1、r2—几何长度参数,m;θ、θ1、θ2—几何角度参数,°。
计算多条水力裂缝在储层中产生的诱导应力场:
式中:σxx、σyy、σzz、σxy、σyz、σxz—初始地应力值分量,MPa;ΣΔσxx、ΣΔσyy、ΣΔσzz、ΣΔσxz—多条裂缝产生的诱导应力值各分量之和,MPa;σ′xx、σ′yy、σ′zz、σ′xy、σ′yz、σ′xz—当前地应力值分量,MPa。
再运用渗流力学理论建立岩压裂过程中的压力场变化模型且该模型的示意图如图4所示;结合以下公式计算Laplace域内各向异性页岩储层压力增量:
上式各中间参数和无因次量表达式为:
u=sf(s) (65)
式中:域内压力场,MPa;L—参考长度,可取水平井长度,m;hf—裂缝高度,m;hr—油藏厚度,m;hrD—无因次油藏厚度,无量纲;kmx—基质***x方向上渗透率,D;kmy—基质***y方向上渗透率,D;kmz—基质***z方向上渗透率,D;km—基质***等效渗透率,D;—单位面积流量,m/s;Q—注入排量,m3/s;μ—液体粘度,mPa·s;s—Laplace变量;Lf—裂缝半长,m;K0—0阶Bessel函数;u—自定义函数;ω—弹性储能,无量纲;λ—无因次窜流系数,无量纲;α—岩块形状因子,无量纲;φ—孔隙度,无量纲;ct—总压缩系数,MPa-1;kfx—裂缝***x方向渗透率,D;kfy—裂缝***y方向渗透率,D;kfz—裂缝***z方向渗透率,D;kf—裂缝***等效渗透率,D;xw—井底x坐标,m;yw—井底y坐标,m;zw—井底z坐标,m;xD—无因次x坐标,无量纲;yD—无因次y坐标,无量纲;xwD—井底无因次x坐标,无量纲;ywD—井底无因次y坐标,无量纲;下标m表示基质***;下标m表示天然裂缝***。
利用Stehfest数值反演得到储层实际压力场Δp(x,y,z,t):
其中:
式中:t—时间,s;N—求和项数,值越大越精确,一般取6~18之间的偶整数。
计算出多条水力裂缝存在时的储层压力场:
p′(x,y,z,t)=pi(x,y,z,0)+ΣΔp(x,y,z,t) (71)
式中:pi—初始储层压力场,MPa;ΣΔp—储层压力增量叠加场,MPa;p′—当前储层压力场,即天然裂缝***压力场,MPa。
最后运用岩石破坏理论建立岩压裂过程中的天然裂缝破坏准则模型且该模型的示意图如图5所示;然后将以上计算所得的当前应力场和压力场带入天然裂缝破坏判别式中进行储层内任一点的破坏状态进行判断:
天然裂缝的张性破坏判别式为:
pnf>pn+St (72)
剪切破坏判别式为:
其中:
pn=σxxnxnx+σxynynx+σxznznx+σxynxny+σyynyny (74)
+σyznzny+σxznxnz+σyznynz+σzznznz
pτ=[(σxxnx+σxyny+σxznz)2
+(σxynx+σyyny+σyznz)2+(σxznx+σyzny+σzznz)2 (75)
-(σxxnxnx+σxynynx+σxznznx+σxynxny+σyynyny
+σyznzny+σxznxnz+σyznynz+σzznznz)2]1/2
式中:σij—储层应力张量分量,i,j指标取x,y,z,MPa;ni—裂缝法向矢量分量,i指标取x,y,z,MPa;α—天然裂缝与水平主应力方向夹角(逼近角),°;—天然裂缝与水平面夹角(倾角),°;pn—天然裂缝壁面所受正应力值,MPa;pτ—天然裂缝壁面所受切应力值,MPa;
Kf—天然裂缝摩擦系数,无量纲;pnf—天然裂缝内压力,等于当前储层压力p′,MPa;St—天然裂缝抗张强度,MPa;τ0—天然裂缝内聚力,MPa。
最后计算压裂破碎区体积,根据计算储层内任意点的破坏状态,对压裂破碎区进行三维表征,结合空间数值积分方法,分别计算张性破碎区体积和剪切破碎区体积,并将两者的空间并集算作总体破碎区体积。
计算过程中求得储层应力场(图6、7、8)、储层压力场(图9)和总体破碎区体积三维空间展布形态(图10);并计算出压裂总体破碎区体积为3.39×106m3。
该方法考虑页岩水平井分段分簇压裂过程中流体扩散渗流场和裂缝诱导应力场同时改变对页岩体天然裂缝的触发破坏机制,建立了一种适合页岩水平井分段分簇压裂破碎区体积计算新模型。由于该模型的计算方法与页岩压裂破碎区实际的物理演化过程相一致,为此,可对压裂破碎区进行准确的定量评价,包括剪切破碎区体积和张性破碎区体积,并获得压裂破碎区体积的三维形态,实现页岩水平井压裂评价的经济化、快速化、准确化,解决了现有页岩水平井压裂时对破碎区体积计算过程费时费力费钱,导致压裂破碎区体积评价不能广泛推广应用的问题。
Claims (5)
1.一种页岩水平井压裂破碎区体积的计算方法,其特征在于,包括以下步骤:A、运用裂缝扩展理论建立页岩压裂过程中的水力裂缝延伸模型;
B、运用岩石力学理论建立岩压裂过程中的应力场变化模型;
C、运用渗流力学理论建立岩压裂过程中的压力场变化模型;
D、运用岩石破坏理论建立岩压裂过程中的天然裂缝破坏准则模型;
E、结合上述4个步骤通过计算流程即可对压裂破碎区体积进行计算和表征;
所述步骤A中水力裂缝延伸模型包括以下计算公式:
裂缝缝宽方程:
缝内流体流动方程:
物质平衡方程:
求解条件:
W(L,t)=0在y≥=L(t) (4)
结合求解条件式(4)~式(6),采用隐式有限差分法对式(1)和式(3)进行求解,得到页岩压裂时各条水力裂缝的几何参数;
式中:W(y)—延伸方向上y处的最大缝宽,mm;L(t)—t时刻水力裂缝长度,m;E—岩石杨氏模量,MPa;p(y)—缝内y处压力,MPa;σn(y)—延伸方向上y处裂缝面受到的正应力,MPa;hf—裂缝缝高,m;ν—岩石泊松比,无因次;μ—压裂液粘度,mPa·s;cL—滤失系数,m/s1/2;q—缝内流量,m3/s;Qi—第i条水力裂缝内压裂液分流量,m3/s;Qpumping—压裂施工泵注总排量,m3/s;M—水裂缝条数,即射孔簇数。
2.如权利要求1所述的一种页岩水平井压裂破碎区体积的计算方法,其特征在于,所述步骤B中应力场变化模型包括以下计算公式:
利用弹性力学理论,计算水力裂缝在附近区域产生的诱导应力场:
Δσyy=ν(Δσx+Δσz) (8)
各几何参数表达式如下:
式中:Δσxx、Δσyy、Δσzz—x、y、z方向上法向诱导应力,MPa;Δσxz—x面上z方向切向诱导应力,MPa;pnet—人工裂缝内净压力,MPa;ν—岩石泊松比,无量纲;c—裂缝高度hf的1/2,m;r、r1、r2—几何长度参数,m;θ、θ1、θ2—几何角度参数,°;
利用叠加原理计算出多条水力裂缝存在时的储层应力场:
式中:σxx、σyy、σzz、σxy、σyz、σxz—初始地应力值分量,MPa;∑Δσxx、∑Δσyy、∑Δσzz、∑Δσxz—多条裂缝产生的诱导应力值各分量之和,MPa;σ′xx、σ′yy、σ′zz、σ′xy、σ′yz、σ′xz—当前地应力值分量,MPa。
3.如权利要求2所述的一种页岩水平井压裂破碎区体积的计算方法,其特征在于,所述步骤C中压力场变化模型包括以下计算公式:
利用双重介质理论来计算缝网压裂,在计算Laplace域内各向异性页岩储层压力增量为:
上式各中间参数和无因次量表达式为:
u=sf(s) (27)
式中:
—Laplace域内压力场,MPa;L—参考长度,可取水平井长度,m;hf—裂缝高度,m;hr—油藏厚度,m;hrD—无因次油藏厚度,无量纲;kmx—基质***x方向上渗透率,D;kmy—基质***y方向上渗透率,D;kmz—基质***z方向上渗透率,D;km—基质***等效渗透率,D;—单位面积流量,m/s;Q—注入排量,m3/s;μ—液体粘度,mPa·s;s—Laplace变量;Lf—裂缝半长,m;K0—0阶Bessel函数;u—自定义函数;ω—弹性储能,无量纲;λ—无因次窜流系数,无量纲;α—岩块形状因子,无量纲;φ—孔隙度,无量纲;ct—总压缩系数,MPa-1;—裂缝***x方向渗透率,D;—裂缝***y方向渗透率,D;—裂缝***z方向渗透率,D;kf—裂缝***等效渗透率,D;xw—井底x坐标,m;yw—井底y坐标,m;zw—井底z坐标,m;xD—无因次x坐标,无量纲;yD—无因次y坐标,无量纲;xwD—井底无因次x坐标,无量纲;ywD—井底无因次y坐标,无量纲;下标m表示基质***;下标m表示天然裂缝***;
求解出Laplace域中的储层压力场后,利用Stehfest数值反演即可得到储层实际压力场Δp(x,y,z,t):
其中:
式中:t—时间,s;N—求和项数,值越大越精确,取6~18之间的偶整数;
利用叠加原理计算出多条水力裂缝存在时的储层压力场:
p′(x,y,z,t)=pi(x,y,z,0)+∑Δp(x,y,z,t) (33)
式中:pi—初始储层压力场,MPa;∑Δp—储层压力增量叠加场,MPa;p′—当前储层压力场,即天然裂缝***压力场,MPa。
4.如权利要求3所述的一种页岩水平井压裂破碎区体积的计算方法,其特征在于,所述步骤D中天然裂缝破坏准则模型包括以下计算公式:
天然裂缝的张性破坏判别式为:
pnf>pn+St (34)
剪切破坏判别式为:
其中:
式中:σij—储层应力张量分量,i,j指标取x,y,z,MPa;ni—裂缝法向矢量分量,i指标取x,y,z,MPa;α—天然裂缝与水平主应力方向夹角(逼近角),°;—天然裂缝与水平面夹角(倾角),°;pn—天然裂缝壁面所受正应力值,MPa;pτ—天然裂缝壁面所受切应力值,MPa;Kf—天然裂缝摩擦系数,无量纲;pnf—天然裂缝内压力,等于当前储层压力p′,MPa;St—天然裂缝抗张强度,MPa;τ0—天然裂缝内聚力,MPa。
5.如权利要求4所述的一种页岩水平井压裂破碎区体积的计算方法,其特征在于,所述步骤E中的计算流程包括以下步骤:
1)根据水平井筒射孔参数确定分段压裂的簇间距和簇数,建立相应的三维笛卡尔直角坐标系;
2)利用上述公式(1)、(3)、(4)~(6)计算出裂缝扩展空间几何尺寸;
3)利用上述公式(7)~(11)、(17)求解出当前地应力场σ′ij;
4)利用式(18)、(31)、(32)求解出当前储层压力场p′;
5)利用式(34)、(35),带入σ′ij和p′判断储层内任一点的天然裂缝是否发生破坏,以及破坏的类型;
6)根据天然裂缝破坏判断结果,利用空间数值积分方法,分别计算张性破碎区体积和剪切破碎区体积,并将两者的空间并集算作总体破碎区体积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510352975.5A CN105019877B (zh) | 2015-06-24 | 2015-06-24 | 一种页岩水平井压裂破碎区体积的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510352975.5A CN105019877B (zh) | 2015-06-24 | 2015-06-24 | 一种页岩水平井压裂破碎区体积的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105019877A CN105019877A (zh) | 2015-11-04 |
CN105019877B true CN105019877B (zh) | 2017-08-08 |
Family
ID=54410120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510352975.5A Active CN105019877B (zh) | 2015-06-24 | 2015-06-24 | 一种页岩水平井压裂破碎区体积的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105019877B (zh) |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105298464B (zh) * | 2015-11-26 | 2018-02-16 | 东北石油大学 | 不含天然裂缝的致密储层体积压裂树状随机缝网描述方法 |
CN105550410B (zh) * | 2015-12-07 | 2018-08-31 | 西南石油大学 | 计算页岩储层水力压裂倾斜裂缝诱导应力的方法 |
CN105547994B (zh) * | 2015-12-31 | 2018-03-13 | 中国石油天然气股份有限公司 | 一种岩石摩擦系数的测试方法 |
CN105735960B (zh) * | 2016-03-22 | 2017-05-17 | 西南石油大学 | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 |
CN105913155B (zh) * | 2016-04-21 | 2019-05-07 | 中国石油天然气股份有限公司 | 考虑应力干扰及压裂液滤失的致密油产能预测方法及*** |
CN105804736B (zh) * | 2016-05-16 | 2019-03-15 | 西安科技大学 | 一种近水平煤层基本顶岩层极限跨距的计算方法 |
CN107524432B (zh) * | 2016-06-21 | 2019-07-02 | 中国石油化工股份有限公司 | 一种压裂施工中的实时调参方法 |
CN107780913B (zh) * | 2016-08-26 | 2020-03-27 | 中国石油化工股份有限公司 | 一种水平井筒穿越多层的常压页岩气体积压裂方法 |
CN106326591B (zh) * | 2016-08-31 | 2019-09-06 | 西南石油大学 | 水力压裂过程中裂缝内压裂液的压力场获取方法及装置 |
CN106650100B (zh) * | 2016-12-23 | 2020-01-10 | 西南石油大学 | 一种实验页岩藏储层水平井交替体积压裂方法 |
CN106844909B (zh) * | 2017-01-05 | 2018-02-06 | 西南石油大学 | 一种致密储层压裂改造体积区的计算方法 |
CN106761647B (zh) * | 2017-01-13 | 2020-08-14 | 中国石油化工股份有限公司 | 一种估算页岩储层压后平面改造面积的方法 |
CN108343433B (zh) * | 2018-02-28 | 2019-11-05 | 西南石油大学 | 用于计算页岩微裂缝缝宽变化下气体质量传输参数的方法 |
CN108468538B (zh) * | 2018-03-23 | 2020-04-21 | 西南石油大学 | 一种页岩水压裂缝扩展预测方法 |
CN108533242A (zh) * | 2018-04-17 | 2018-09-14 | 成都理工大学 | 一种煤层气压裂井排采压力安全窗分析方法 |
CN109236255B (zh) * | 2018-09-03 | 2020-07-10 | 中国石油大学(华东) | 一种水平井压裂潜力评价方法和装置 |
CN110175723B (zh) * | 2019-05-31 | 2020-04-21 | 西南石油大学 | 一种页岩气井重复压裂裂缝起裂方向与延伸路径预测方法 |
CN110516407B (zh) * | 2019-09-16 | 2022-02-11 | 西南石油大学 | 一种裂缝性储层水平井段内多簇压裂裂缝复杂度计算方法 |
CN111428425B (zh) * | 2020-03-19 | 2020-12-29 | 西南石油大学 | 一种页岩油藏变裂缝渗透率分段压裂水平井产能计算方法 |
WO2021212442A1 (zh) * | 2020-04-24 | 2021-10-28 | 中国矿业大学(北京) | 模拟储层压裂应力场演化的实验方法及*** |
CN111859260B (zh) * | 2020-07-03 | 2024-07-16 | 中国科学技术大学 | 一种利用压裂施工数据获得srv区域的方法 |
CN111720104B (zh) * | 2020-08-04 | 2022-03-11 | 西南石油大学 | 一种预测裂缝性储层多级压裂裂缝形态的方法 |
CN112949129B (zh) * | 2021-03-03 | 2023-03-03 | 西南石油大学 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
CN113846870B (zh) * | 2021-10-12 | 2022-12-06 | 上海建工集团股份有限公司 | 以既有裂缝为导向的混凝土支撑静力破碎方法 |
CN117633409B (zh) * | 2024-01-25 | 2024-04-09 | 中国科学院地质与地球物理研究所 | 计算页岩油气储层压裂缝网渗流参数的方法、***及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103760602A (zh) * | 2014-01-27 | 2014-04-30 | 中国石油天然气集团公司 | 一种三维水力压裂及微地震事件正演模拟的方法及*** |
CN103773356A (zh) * | 2014-01-21 | 2014-05-07 | 成都新柯力化工科技有限公司 | 一种页岩气开采用颗粒料及其制备方法 |
CN104278980A (zh) * | 2013-07-09 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种采用改造体积优化致密油水平井缝网参数的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9618652B2 (en) * | 2011-11-04 | 2017-04-11 | Schlumberger Technology Corporation | Method of calibrating fracture geometry to microseismic events |
-
2015
- 2015-06-24 CN CN201510352975.5A patent/CN105019877B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104278980A (zh) * | 2013-07-09 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种采用改造体积优化致密油水平井缝网参数的方法 |
CN103773356A (zh) * | 2014-01-21 | 2014-05-07 | 成都新柯力化工科技有限公司 | 一种页岩气开采用颗粒料及其制备方法 |
CN103760602A (zh) * | 2014-01-27 | 2014-04-30 | 中国石油天然气集团公司 | 一种三维水力压裂及微地震事件正演模拟的方法及*** |
Non-Patent Citations (2)
Title |
---|
"水平井分段多簇压裂缝间干扰影响分析";赵金洲等;《天然气地球科学》;20150331 * |
"重新界定破碎性地层概念及其稳定性分析";唐一元等;《内蒙古石油化工》;20081231 * |
Also Published As
Publication number | Publication date |
---|---|
CN105019877A (zh) | 2015-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105019877B (zh) | 一种页岩水平井压裂破碎区体积的计算方法 | |
Wu et al. | Physical simulation study on the hydraulic fracture propagation of coalbed methane well | |
Zhu et al. | Coupled flow-stress-damage simulation of deviated-wellbore fracturing in hard-rock | |
CN106285646A (zh) | 基于多信息融合的钻井漏失层位识别方法 | |
CN105484710B (zh) | 一种选择斜井压裂射孔方位的方法 | |
CN110750930B (zh) | 一种基于裂缝连续体模型预测裂缝性储层应力演化的方法 | |
CN105074125A (zh) | 针对微震事件对裂缝几何形状进行校准的方法 | |
GB2563241A (en) | Hydraulic fracturing simulation | |
CN102575510A (zh) | 控制地质力学油藏***中的出砂的计算机实施的***和方法 | |
Zhou et al. | Numerical modeling and investigation of fracture propagation with arbitrary orientation through fluid injection in tight gas reservoirs with combined XFEM and FVM | |
CN112682016B (zh) | 油气田薄互层储层穿层压裂参数的确定方法及压裂方法 | |
Wu et al. | Upper limit analysis of stability of the water-resistant rock mass of a Karst tunnel face considering the seepage force | |
Kong et al. | Seepage problems on fractured rock accompanying with mass loss during excavation in coal mines with karst collapse columns | |
Vishal et al. | Three dimensional modelling of propagation of hydraulic fractures in shale at different injection pressures | |
Zhao et al. | Numerical study of the effect of natural fractures on shale hydraulic fracturing based on the continuum approach | |
Tian et al. | Shale gas production from reservoirs with hierarchical multiscale structural heterogeneities | |
Ou et al. | Numerical analysis of seepage flow characteristic of collapse column under the influence of mining | |
Ravaji et al. | Introducing optimized validated meshing system for wellbore stability analysis using 3D finite element method | |
CN106321082A (zh) | 用于求取盐膏岩地层最大套管外载荷的方法及*** | |
Zhou et al. | Fluid flow through porous sandstone with overprinting and intersecting geological structures of various types | |
CN104712299B (zh) | 适合气井控水增气压裂的设计方法 | |
CN110850470A (zh) | 一种利用地震资料计算地应力差异系数的方法 | |
Nadimi | State-based peridynamics simulation of hydraulic fracture phenomenon in geological media | |
CN104392131A (zh) | 一种水驱沙过程中破碎岩石渗流场计算方法 | |
CN107066744A (zh) | 富水地层中地铁车站的抗震分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |