CN112949129A - 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 - Google Patents
一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 Download PDFInfo
- Publication number
- CN112949129A CN112949129A CN202110236276.XA CN202110236276A CN112949129A CN 112949129 A CN112949129 A CN 112949129A CN 202110236276 A CN202110236276 A CN 202110236276A CN 112949129 A CN112949129 A CN 112949129A
- Authority
- CN
- China
- Prior art keywords
- fracturing
- cluster
- fracture
- horizontal well
- crack
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明公开一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,包括计算压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力;计算各簇裂缝尖端转向角度;计算压裂过程中各簇裂缝分得压裂液流量;继续进行下一时步内的裂缝延伸计算,直到压裂时间结束,确定压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度;绘制各簇裂缝长度随时间变化曲线、压裂结束时各簇裂缝延伸空间展布图。本发明考虑了深层页岩水平井各簇射孔位置处地层应力非均质性变化对裂缝起裂延伸的影响,可对各簇裂缝的起裂时间与延伸路径进行准确预测;解决了深层页岩地层应力分布不均匀,水力压裂时多簇裂缝起裂延伸不同步,缺乏快速准确的多簇裂缝异步起裂延伸计算方法的问题。
Description
技术领域
本发明涉及一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,属于页岩气开发技术领域。
背景技术
近年来,随着水平井分段分簇压裂技术的广泛应用,国内外已逐步实现了中浅层页岩气的商业化高效开发,并正在向深层页岩气领域进军。相比于中浅层页岩气,深层页岩气地层应力高、压力高、温度高,并且地质构造相对复杂,地层应力变化较大。因此,深层页岩气水平井段不同位置处地层应力条件存在明显差异,导致压裂过程中不同射孔簇位置处的水力裂缝起裂时间有所不同。由此可见,深层页岩地层应力非均匀特征,可显著影响水平井压裂多簇裂缝起裂延伸行为。
中浅层页岩气地层地质构造运动相对平缓,地层应力在一定区域内变化不大,水平井压裂过程中各簇水力裂缝同步起裂延伸。2010年,Meyer等人基于有限差分方法,模拟计算了页岩气水平井压裂裂缝同步起裂正交平面延伸行为。压裂过程中,由于多簇裂缝相互之间存在应力干扰效应,可能出现非均匀、非平面延伸现象。2015年,郭建春等人基于有限元方法,模拟计算了页岩气水平井压裂多簇裂缝同步起裂非均匀平面延伸行为。随后,Wu等人(2015),Kresse和Weng(2018)、Liu等人(2019)、Zheng等人(2020)分别基于边界元、离散元等方法,模拟计算了页岩气水平井压裂多簇裂缝同步起裂非均匀非平面延伸行为。
综上所述,目前所有页岩气水平井压裂裂缝起裂延伸模拟研究中,多簇裂缝均为同步起裂,仅适用于中浅层页岩气地层。而深层页岩气水平井压裂过程中,由于地层应力分布不均匀,多簇裂缝通常存在异步起裂延伸的现象。因此,亟需建立一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,充分考虑深层页岩水平井各簇射孔位置处地层应力非均质性变化对裂缝起裂延伸的影响,进而对各簇裂缝的起裂时间与延伸路径进行准确预测,为深层页岩压裂优化设计奠定理论基础。
发明内容
为了克服现有技术中的深层页岩地层应力分布不均匀,水力压裂时多簇裂缝起裂延伸不同步,缺乏快速准确的多簇裂缝异步起裂延伸计算方法的问题,本发明提供一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法。
本发明解决上述技术问题所提供的技术方案是:一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,包括:
根据目标深层页岩水平井的地层参数、射孔参数、压裂参数以及深层页岩水平井压裂水力裂缝起裂延伸模型确定压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力;
根据目标深层页岩水平井的地层参数、射孔参数以及深层页岩水平井压裂水力裂缝转向延伸模型确定压裂过程中各簇裂缝尖端转向角度;
根据目标深层页岩水平井的地层参数、射孔参数、压裂参数、缝内压力以及深层页岩水平井压裂多簇裂缝流量分配模型确定压裂过程中各簇裂缝分得压裂液流量;
根据压裂过程中各簇裂缝尖端转向角度、各簇裂缝分得压裂液流量继续进行下一时步内的裂缝延伸计算,直到压裂时间结束,确定压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度数据;
根据压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度数据绘制压裂过程中各簇裂缝长度随时间变化曲线、压裂结束时各簇裂缝延伸空间展布图。
进一步的技术方案是,所述深层页岩水平井压裂水力裂缝起裂延伸模型包括:
物质平衡方程:
式中:q为裂缝内流量,m3/s;hf为裂缝高度,m;wf为裂缝开度,m;s为裂缝长度方向坐标,m;t为时间,s;qL为压裂液滤失速度,m/s;
缝内流体压降方程:
式中:p为裂缝内压力,Pa;s为裂缝长度方向坐标,m;μ为液体粘度,Pa·s;
压裂液体滤失方程:
式中:CL为地层滤失系数,m/s0.5;τ为开始滤失时间,s;
裂缝宽度方程:
wf(s)=(Dn)j (4)
式中:(Dn)j为裂缝缝长s处对应的裂缝j单元法向位移量,m;
裂缝高度方程:
式中:KIc为页岩断裂韧性,Pa·m0.5;
裂缝闭合应力方程:
σc=σhmincos2(θ)+σHmaxsin2(θ) (6)
式中:θ为裂缝逼近角,即与地层最小水平主应力方向夹角,°;σhmin为地层最小水平主应力,Pa;σHmax为地层最大水平主应力,Pa;
裂缝起裂延伸边界条件与初始条件方程:
式中:Q为水力压裂泵注排量,m3/min;Lf为裂缝半长,m。
进一步的技术方案是,根据目标深层页岩水平井的地层参数、射孔参数、压裂参数以及深层页岩水平井压裂水力裂缝起裂延伸模型确定压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力包括:
收集目标深层页岩水平井地层参数、射孔参数、压裂参数;
联立方程式(1)~(6),结合边界条件与初始条件方程式(7),通过有限差分方法与Picard迭代法,分别计算压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力。
进一步的技术方案是,所述深层页岩水平井压裂水力裂缝转向延伸模型包括:
裂缝离散单元坐标转换方程:
式中:ζij、ξij为局部坐标值,m;xi、xj为全局坐标下裂缝i、j单元中心的x轴坐标值,m;yi、yj为全局坐标下裂缝i、j单元中心的y轴坐标值,m;
裂缝离散单元应力-应变平衡方程:
其中:
式中:(σt)i、(σn)i为裂缝i单元在局部坐标系内所受切应力和正应力,Pa;(Dn)j、(Ds)j为裂缝j单元的法向位移量与切向位移量,m;(Att)ij、(Ant)ij、(Atn)ij、(Ann)ij为裂缝j单元切向位移和法向位移不连续量分别在i单元上引起的切向应力分量和法相应力分量;E为岩石杨氏模量,Pa;ν为地层泊松比,无量纲;nj为全局坐标y轴与j单元局部坐标ζ轴夹角余弦值,无量纲;lj为全局坐标x轴与j单元局部坐标ξ轴夹角余弦值,无量纲;F3~F6为Papkovitch函数偏导方程;
裂缝尖端单元转向角度方程:
(Dn)tipsinθtip+(Ds)tip(3cosθtip-1)=0 (14)
式中:(Dn)tip、(Ds)tip为裂缝尖端单元法向、切向位移量,m;θtip为裂缝尖端转向角度,°。
进一步的技术方案是,根据目标深层页岩水平井的地层参数、射孔参数以及深层页岩水平井压裂水力裂缝转向延伸模型确定压裂过程中各簇裂缝尖端转向角度包括:
将多簇裂缝划分为离散单元,利用方程式(8)计算全局坐标下所有裂缝离散单元中心坐标值;
利用方程式(9)~(13)计算所有裂缝离散单元的法向位移量与切向位移量;
利用方程式(14)计算各簇裂缝尖端转向角度。
进一步的技术方案是,所述深层页岩水平井压裂多簇裂缝流量分配模型包括:
流体沿程压降方程:
其中:
式中:pheel为水平井跟端压力,Pa;pfi,i为第i簇裂缝缝口处压力,Pa;Δppf,i为第i簇射孔孔眼处的摩阻压降,Pa;Δpw,j为第j段水平井段内流体流动压降,Pa;pnet,i为第i簇裂缝缝口净压力,Pa;σhmin,i为第i簇裂缝射孔位置地层最小水平主应力,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表示各段水平井段编号。
进一步的技术方案是,根据目标深层页岩水平井的地层参数、射孔参数、压裂参数、缝内压力以及深层页岩水平井压裂多簇裂缝流量分配模型确定压裂过程中各簇裂缝分得压裂液流量包括:
从目标深层页岩水平井地层最小水平主应力数据中,提取各簇射孔位置处地层最小水平主应力,带入方程式(16)中;
联立方程式(15)~(17),通过牛顿迭代法,构造Jacobian系数矩阵,计算压裂过程中各簇裂缝分得压裂液流量。
本发明具有以下有益效果:与现有技术相比,该方法专门针对深层页岩水平井压裂过程中多簇裂缝起裂延伸力学行为复杂,建立了深层页岩水平井压裂水力裂缝起裂延伸模型、水力裂缝转向延伸模型、多簇裂缝流量分配模型,从而提出了一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法;由于该方法充分考虑了深层页岩水平井各簇射孔位置处地层应力非均质性变化对裂缝起裂延伸的影响,可对各簇裂缝的起裂时间与延伸路径进行准确预测;解决了深层页岩地层应力分布不均匀,水力压裂时多簇裂缝起裂延伸不同步,缺乏快速准确的多簇裂缝异步起裂延伸计算方法的问题。
附图说明
图1是本发明方法的计算流程框图;
图2是实施例目标井压裂过程中各簇裂缝长度随时间变化曲线;
图3是实施例目标井压裂段各簇水力裂缝延伸空间展布图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明的一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,包括以下步骤:
步骤1、收集目标深层页岩水平井地层参数(地层最大水平主应力、地层最小水平主应力、岩石杨氏模量、岩石泊松比、储层垂深、储层厚度、储层岩石断裂韧性、地层滤失系数)、射孔参数(簇间距、簇数、射孔孔眼数量、射孔孔眼直径、孔眼流量系数)、压裂参数(压裂排量、压裂液量、泵注时间、压裂管柱内径、压裂液粘度);
步骤2、运用流固耦合理论,建立深层页岩水平井压裂水力裂缝起裂延伸模型;
物质平衡方程:
式中:q为裂缝内流量,m3/s;hf为裂缝高度,m;wf为裂缝开度,m;s为裂缝长度方向坐标,m;t为时间,s;qL为压裂液滤失速度,m/s;
缝内流体压降方程:
式中:p为裂缝内压力,Pa;s为裂缝长度方向坐标,m;μ为液体粘度,Pa·s;
压裂液体滤失方程:
式中:CL为地层滤失系数,m/s0.5;τ为开始滤失时间,s;
裂缝宽度方程:
wf(s)=(Dn)j (4)
式中:(Dn)j为裂缝缝长s处对应的裂缝j单元法向位移量,m;
裂缝高度方程:
式中:KIc为页岩断裂韧性,Pa·m0.5;
裂缝闭合应力方程:
σc=σhmincos2(θ)+σHmaxsin2(θ) (6)
式中:θ为裂缝逼近角,即与地层最小水平主应力方向夹角,°;σhmin为地层最小水平主应力,Pa;σHmax为地层最大水平主应力,Pa;
裂缝起裂延伸边界条件与初始条件方程:
式中:Q为水力压裂泵注排量,m3/min;Lf为裂缝半长,m。
步骤3、联立方程式(1)~(6),结合边界条件与初始条件方程式(7),通过有限差分方法与Picard迭代法,分别计算压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力;
步骤4、运用岩石力学理论,建立深层页岩水平井压裂水力裂缝转向延伸模型;
裂缝离散单元坐标转换方程:
式中:ζij、ξij为局部坐标值,m;xi、xj为全局坐标下裂缝i、j单元中心的x轴坐标值,m;yi、yj为全局坐标下裂缝i、j单元中心的y轴坐标值,m;
裂缝离散单元应力-应变平衡方程:
其中:
式中:(σt)i、(σn)i为裂缝i单元在局部坐标系内所受切应力和正应力,Pa;(Dn)j、(Ds)j为裂缝j单元的法向位移量与切向位移量,m;(Att)ij、(Ant)ij、(Atn)ij、(Ann)ij为裂缝j单元切向位移和法向位移不连续量分别在i单元上引起的切向应力分量和法相应力分量;E为岩石杨氏模量,Pa;ν为地层泊松比,无量纲;nj为全局坐标y轴与j单元局部坐标ζ轴夹角余弦值,无量纲;lj为全局坐标x轴与j单元局部坐标ξ轴夹角余弦值,无量纲;F3~F6为Papkovitch函数偏导方程;
裂缝尖端单元转向角度方程:
(Dn)tipsinθtip+(Ds)tip(3cosθtip-1)=0 (14)
式中:(Dn)tip、(Ds)tip为裂缝尖端单元法向、切向位移量,m;θtip为裂缝尖端转向角度,°;
步骤5、将多簇裂缝划分为离散单元,利用方程式(8)计算全局坐标下所有裂缝离散单元中心坐标值;再利用方程式(9)~(13)计算所有裂缝离散单元的法向位移量与切向位移量;最后利用方程式(14)计算各簇裂缝尖端转向角度;
步骤6、运用流体力学理论,建立深层页岩水平井压裂多簇裂缝流量分配模型;
流体沿程压降方程:
其中:
式中:pheel为水平井跟端压力,Pa;pfi,i为第i簇裂缝缝口处压力,Pa;Δppf,i为第i簇射孔孔眼处的摩阻压降,Pa;Δpw,j为第j段水平井段内流体流动压降,Pa;pnet,i为第i簇裂缝缝口净压力,Pa;σhmin,i为第i簇裂缝射孔位置地层最小水平主应力,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表示各段水平井段编号;
步骤7、从目标深层页岩水平井地层最小水平主应力数据中,提取各簇裂缝射孔位置处地层最小水平主应力,带入方程式(16)中;再联立方程式(15)~(17),通过牛顿迭代法,构造Jacobian系数矩阵,计算压裂过程中各簇裂缝分得压裂液流量;
步骤8、返回步骤2,继续进行下一时步内的裂缝延伸计算,直到压裂时间结束;
其中,下一时步内的各簇裂缝尖端转向角度由上一时步的步骤5内获得;下一时步内的各簇裂缝分得压裂液流量由上一时步的步骤7内获得。
步骤9、数据输出:输出压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度数据;
步骤10、图像绘制:绘制压裂过程中各簇裂缝长度随时间变化曲线、压裂结束时各簇裂缝延伸空间展布图。
实施例1:
一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,包括以下步骤:
已知某深层页岩气藏区域内目标水平井地质、射孔、压裂数据如表1所示。
表1某深层页岩气藏区域内目标水平井相关数据表
首先,运用流固耦合理论,建立深层页岩水平井压裂水力裂缝起裂延伸模型,又结合以下公式计算实施例1中六簇水力裂缝起裂延伸情况:
具体步骤为:①输入目标深层页岩水平井地层参数、射孔参数、压裂参数;②联立方程式(1)~(6),结合边界条件与初始条件方程式(7),通过有限差分方法与Picard迭代法,分别计算压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力。
随后,运用岩石力学理论建立深层页岩水平井压裂水力裂缝转向延伸模型,又结合以下公式计算实施例1中六簇水力裂缝转向延伸情况:
具体步骤为:①将多簇裂缝划分为离散单元,利用方程式(8)计算全局坐标下所有裂缝离散单元中心坐标值;②利用方程式(9)~(13)计算所有裂缝离散单元的法向位移量与切向位移量;③利用方程式(14)计算各簇裂缝尖端转向角度。
随后,运用流体力学理论建立深层页岩水平井压裂多簇裂缝流量分配模型,又结合以下公式计算实施例1中六簇水力裂缝流量分配情况:
具体步骤为:①从目标深层页岩水平井地层最小水平主应力数据中,提取各簇裂缝射孔位置处地层最小水平主应力,带入方程式(16)中;②联立方程式(15)~(17),通过牛顿迭代法,构造Jacobian系数矩阵,计算压裂过程中各簇裂缝分得压裂液流量。
最后,运用如图1所示的本发明方法的数值计算流程框图开展实例计算,并根据计算结果,输出压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度数据(如表2所示),分别绘制压裂过程中各簇裂缝长度随时间变化曲线(如图2所示)、压裂结束时各簇裂缝延伸空间展布图(如图3所示)。
表2压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度数据表
簇号 | 水力裂缝半长(m) | 水力裂缝宽度(mm) | 水力裂缝高度(m) |
1 | 21 | 5.97 | 46.33 |
2 | 164 | 5.36 | 41.59 |
3 | 139 | 7.04 | 54.64 |
4 | 136 | 6.98 | 54.17 |
5 | 76 | 5.36 | 41.59 |
6 | 218 | 6.40 | 49.68 |
该方法专门针对深层页岩水平井压裂过程中多簇裂缝起裂延伸力学行为复杂,建立了深层页岩水平井压裂水力裂缝起裂延伸模型、水力裂缝转向延伸模型、多簇裂缝流量分配模型,从而提出了一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法。由于该方法充分考虑了深层页岩水平井各簇射孔位置处地层应力非均质性变化对裂缝起裂延伸的影响,可对各簇裂缝的起裂时间与延伸路径进行准确预测。解决了深层页岩地层应力分布不均匀,水力压裂时多簇裂缝起裂延伸不同步,缺乏快速准确的多簇裂缝异步起裂延伸计算方法的问题。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (7)
1.一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,其特征在于,包括:
根据目标深层页岩水平井的地层参数、射孔参数、压裂参数以及深层页岩水平井压裂水力裂缝起裂延伸模型确定压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力;
根据目标深层页岩水平井的地层参数、射孔参数以及深层页岩水平井压裂水力裂缝转向延伸模型确定压裂过程中各簇裂缝尖端转向角度;
根据目标深层页岩水平井的地层参数、射孔参数、压裂参数、缝内压力以及深层页岩水平井压裂多簇裂缝流量分配模型确定压裂过程中各簇裂缝分得压裂液流量;
根据压裂过程中各簇裂缝尖端转向角度、各簇裂缝分得压裂液流量继续进行下一时步内的裂缝延伸计算,直到压裂时间结束,确定压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度数据;
根据压裂结束时各簇裂缝长度、裂缝高度、裂缝宽度数据绘制压裂过程中各簇裂缝长度随时间变化曲线、压裂结束时各簇裂缝延伸空间展布图。
2.根据权利要求1所述的一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,其特征在于,所述深层页岩水平井压裂水力裂缝起裂延伸模型包括:
物质平衡方程:
式中:q为裂缝内流量,m3/s;hf为裂缝高度,m;wf为裂缝开度,m;s为裂缝长度方向坐标,m;t为时间,s;qL为压裂液滤失速度,m/s;
缝内流体压降方程:
式中:p为裂缝内压力,Pa;s为裂缝长度方向坐标,m;μ为液体粘度,Pa·s;
压裂液体滤失方程:
式中:CL为地层滤失系数,m/s0.5;τ为开始滤失时间,s;
裂缝宽度方程:
wf(s)=(Dn)j (4)
式中:(Dn)j为裂缝缝长s处对应的裂缝j单元法向位移量,m;
裂缝高度方程:
式中:KIc为页岩断裂韧性,Pa·m0.5;
裂缝闭合应力方程:
σc=σhmincos2(θ)+σHmaxsin2(θ) (6)
式中:θ为裂缝逼近角,即与地层最小水平主应力方向夹角,°;σhmin为地层最小水平主应力,Pa;σHmax为地层最大水平主应力,Pa;
裂缝起裂延伸边界条件与初始条件方程:
式中:Q为水力压裂泵注排量,m3/min;Lf为裂缝半长,m。
3.根据权利要求2所述的一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,其特征在于,根据目标深层页岩水平井的地层参数、射孔参数、压裂参数以及深层页岩水平井压裂水力裂缝起裂延伸模型确定压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力包括:
收集目标深层页岩水平井地层参数、射孔参数、压裂参数;
联立方程式(1)~(6),结合边界条件与初始条件方程式(7),通过有限差分方法与Picard迭代法,分别计算压裂过程中各簇裂缝的长度、高度、宽度,以及缝内压力。
4.根据权利要求1所述的一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,其特征在于,所述深层页岩水平井压裂水力裂缝转向延伸模型包括:
裂缝离散单元坐标转换方程:
式中:ζij、ξij为局部坐标值,m;xi、xj为全局坐标下裂缝i、j单元中心的x轴坐标值,m;yi、yj为全局坐标下裂缝i、j单元中心的y轴坐标值,m;
裂缝离散单元应力-应变平衡方程:
其中:
式中:(σt)i、(σn)i为裂缝i单元在局部坐标系内所受切应力和正应力,Pa;(Dn)j、(Ds)j为裂缝j单元的法向位移量与切向位移量,m;(Att)ij、(Ant)ij、(Atn)ij、(Ann)ij为裂缝j单元切向位移和法向位移不连续量分别在i单元上引起的切向应力分量和法相应力分量;E为岩石杨氏模量,Pa;ν为地层泊松比,无量纲;nj为全局坐标y轴与j单元局部坐标ζ轴夹角余弦值,无量纲;lj为全局坐标x轴与j单元局部坐标ξ轴夹角余弦值,无量纲;F3~F6为Papkovitch函数偏导方程;
裂缝尖端单元转向角度方程:
(Dn)tipsinθtip+(Ds)tip(3cosθtip-1)=0 (14)
式中:(Dn)tip、(Ds)tip为裂缝尖端单元法向、切向位移量,m;θtip为裂缝尖端转向角度,°。
5.根据权利要求4所述的一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,其特征在于,根据目标深层页岩水平井的地层参数、射孔参数以及深层页岩水平井压裂水力裂缝转向延伸模型确定压裂过程中各簇裂缝尖端转向角度包括:
将多簇裂缝划分为离散单元,利用方程式(8)计算全局坐标下所有裂缝离散单元中心坐标值;
利用方程式(9)~(13)计算所有裂缝离散单元的法向位移量与切向位移量;
利用方程式(14)计算各簇裂缝尖端转向角度。
6.根据权利要求1所述的一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,其特征在于,所述深层页岩水平井压裂多簇裂缝流量分配模型包括:
流体沿程压降方程:
其中:
式中:pheel为水平井跟端压力,Pa;pfi,i为第i簇裂缝缝口处压力,Pa;Δppf,i为第i簇射孔孔眼处的摩阻压降,Pa;Δpw,j为第j段水平井段内流体流动压降,Pa;pnet,i为第i簇裂缝缝口净压力,Pa;σhmin,i为第i簇裂缝射孔位置地层最小水平主应力,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表示各段水平井段编号。
7.根据权利要求6所述的一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法,其特征在于,根据目标深层页岩水平井的地层参数、射孔参数、压裂参数、缝内压力以及深层页岩水平井压裂多簇裂缝流量分配模型确定压裂过程中各簇裂缝分得压裂液流量包括:
从目标深层页岩水平井地层最小水平主应力数据中,提取各簇射孔位置处地层最小水平主应力,带入方程式(16)中;
联立方程式(15)~(17),通过牛顿迭代法,构造Jacobian系数矩阵,计算压裂过程中各簇裂缝分得压裂液流量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110236276.XA CN112949129B (zh) | 2021-03-03 | 2021-03-03 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110236276.XA CN112949129B (zh) | 2021-03-03 | 2021-03-03 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112949129A true CN112949129A (zh) | 2021-06-11 |
CN112949129B CN112949129B (zh) | 2023-03-03 |
Family
ID=76247453
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110236276.XA Active CN112949129B (zh) | 2021-03-03 | 2021-03-03 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112949129B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113469832A (zh) * | 2021-08-09 | 2021-10-01 | 中国石油天然气股份有限公司 | 一种致密砂岩加砂压裂用液强度设计方法 |
CN113971378A (zh) * | 2021-10-27 | 2022-01-25 | 西南石油大学 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
CN113982568A (zh) * | 2021-11-19 | 2022-01-28 | 中国石油天然气股份有限公司 | 致密油气多裂缝竞争起裂井底压力预测方法 |
CN114441337A (zh) * | 2022-01-14 | 2022-05-06 | 中国石油大学(北京) | 一种金属断裂韧性裂纹尖端张开位移的检测方法 |
CN115659736A (zh) * | 2022-10-19 | 2023-01-31 | 西南石油大学 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
CN116029167A (zh) * | 2023-03-31 | 2023-04-28 | 中国石油大学(华东) | 一种水平井多簇压裂非均匀射孔的定量设计方法 |
CN117408168A (zh) * | 2023-12-15 | 2024-01-16 | 中国石油大学(华东) | 一种融合物理条件约束的压裂裂缝扩展快速预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105019877A (zh) * | 2015-06-24 | 2015-11-04 | 西南石油大学 | 一种页岩水平井压裂破碎区体积的计算方法 |
CN106650100A (zh) * | 2016-12-23 | 2017-05-10 | 西南石油大学 | 一种实验页岩藏储层水平井交替体积压裂方法 |
CN108442911A (zh) * | 2018-02-28 | 2018-08-24 | 西南石油大学 | 一种页岩气水平井重复压裂水力裂缝参数优化设计方法 |
CN110175723A (zh) * | 2019-05-31 | 2019-08-27 | 西南石油大学 | 一种页岩气井重复压裂裂缝起裂方向与延伸路径预测方法 |
CN111322050A (zh) * | 2020-04-24 | 2020-06-23 | 西南石油大学 | 一种页岩水平井段内密切割暂堵压裂施工优化方法 |
-
2021
- 2021-03-03 CN CN202110236276.XA patent/CN112949129B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105019877A (zh) * | 2015-06-24 | 2015-11-04 | 西南石油大学 | 一种页岩水平井压裂破碎区体积的计算方法 |
CN106650100A (zh) * | 2016-12-23 | 2017-05-10 | 西南石油大学 | 一种实验页岩藏储层水平井交替体积压裂方法 |
CN108442911A (zh) * | 2018-02-28 | 2018-08-24 | 西南石油大学 | 一种页岩气水平井重复压裂水力裂缝参数优化设计方法 |
CN110175723A (zh) * | 2019-05-31 | 2019-08-27 | 西南石油大学 | 一种页岩气井重复压裂裂缝起裂方向与延伸路径预测方法 |
CN111322050A (zh) * | 2020-04-24 | 2020-06-23 | 西南石油大学 | 一种页岩水平井段内密切割暂堵压裂施工优化方法 |
Non-Patent Citations (3)
Title |
---|
LAN REN 等: "Stimulated Reservoir Volume Estimation for Shale Gas Fracturing:Mechanism and Modeling Approach", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
任岚 等: "页岩压裂水平井增产改造体积的动态演化模型", 《应用数学和力学》 * |
林然: "页岩压裂水平井SRV动态模拟的理论研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113469832A (zh) * | 2021-08-09 | 2021-10-01 | 中国石油天然气股份有限公司 | 一种致密砂岩加砂压裂用液强度设计方法 |
CN113971378A (zh) * | 2021-10-27 | 2022-01-25 | 西南石油大学 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
CN113971378B (zh) * | 2021-10-27 | 2022-08-02 | 西南石油大学 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
CN113982568A (zh) * | 2021-11-19 | 2022-01-28 | 中国石油天然气股份有限公司 | 致密油气多裂缝竞争起裂井底压力预测方法 |
CN113982568B (zh) * | 2021-11-19 | 2024-04-16 | 中国石油天然气股份有限公司 | 致密油气多裂缝竞争起裂井底压力预测方法 |
CN114441337A (zh) * | 2022-01-14 | 2022-05-06 | 中国石油大学(北京) | 一种金属断裂韧性裂纹尖端张开位移的检测方法 |
CN115659736A (zh) * | 2022-10-19 | 2023-01-31 | 西南石油大学 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
CN115659736B (zh) * | 2022-10-19 | 2023-11-03 | 西南石油大学 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
CN116029167A (zh) * | 2023-03-31 | 2023-04-28 | 中国石油大学(华东) | 一种水平井多簇压裂非均匀射孔的定量设计方法 |
CN117408168A (zh) * | 2023-12-15 | 2024-01-16 | 中国石油大学(华东) | 一种融合物理条件约束的压裂裂缝扩展快速预测方法 |
CN117408168B (zh) * | 2023-12-15 | 2024-03-15 | 中国石油大学(华东) | 一种融合物理条件约束的压裂裂缝扩展快速预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112949129B (zh) | 2023-03-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112949129B (zh) | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 | |
CN110210144B (zh) | 一种暂堵剂促进水平井压裂裂缝均匀扩展的优化设计方法 | |
CN110222477B (zh) | 维持水平井分段压裂裂缝均衡扩展的射孔参数优化方法 | |
CN113389534B (zh) | 一种水平井密切割压裂裂缝扩展预测与设计参数优化方法 | |
CN106874544B (zh) | 一种页岩储层改造体积的地质表征方法 | |
CN110175723B (zh) | 一种页岩气井重复压裂裂缝起裂方向与延伸路径预测方法 | |
CN113850029B (zh) | 一种页岩气水平井密切割压裂射孔参数优化设计方法 | |
CN113076676B (zh) | 非常规油气藏水平井压裂缝网扩展与生产动态耦合方法 | |
CN113971378B (zh) | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 | |
CN116401897B (zh) | 一种基于近似解和能量方程的半解析裂缝扩展模拟方法 | |
CN115618759A (zh) | 页岩气地层压裂施工参数优化方法 | |
CN111734380A (zh) | 一种水平井多级压裂裂缝扩展形态的快速预测方法 | |
CN112761609B (zh) | 一种水力压裂作业中高效铺置支撑剂的优化方法 | |
CN112878977B (zh) | 一种致密储层水平井多簇压裂射孔簇间距优化方法 | |
CN111271042A (zh) | 一种针对天然裂缝发育储层水力压裂形成复杂缝网的排量优化方法 | |
CN112096362B (zh) | 非常规储层多簇射孔竞争起裂与扩展模拟方法及装置 | |
CN111950209A (zh) | 一种水力压裂裂缝延伸轨迹的控制方法 | |
CN110439519A (zh) | 一种基于极限限流设计的压裂方法及系统 | |
CN106777663B (zh) | 一种考虑天然裂缝的压裂液滤失速度计算方法 | |
Dang et al. | Simulation of effective fracture length of prepad acid fracturing considering multiple leak-off effect | |
CN115659736A (zh) | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 | |
CN114372428B (zh) | 砂砾岩储层水平井段内多簇压裂裂缝延伸跨尺度模拟方法 | |
CN114233261A (zh) | 一种低成本分段压裂实现油气井均匀改造的方法 | |
CN111734394B (zh) | 一种确定致密油藏压裂井不定常流井底压力的方法 | |
CN115374681B (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 |