CN113109162B - 一种基于热流固耦合的岩石起裂压力计算方法 - Google Patents
一种基于热流固耦合的岩石起裂压力计算方法 Download PDFInfo
- Publication number
- CN113109162B CN113109162B CN202110394846.8A CN202110394846A CN113109162B CN 113109162 B CN113109162 B CN 113109162B CN 202110394846 A CN202110394846 A CN 202110394846A CN 113109162 B CN113109162 B CN 113109162B
- Authority
- CN
- China
- Prior art keywords
- fluid
- pressure
- mpa
- equation
- stress
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
- G01N3/10—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces generated by pneumatic or hydraulic pressure
- G01N3/12—Pressure testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/003—Generation of the force
- G01N2203/0042—Pneumatic or hydraulic means
- G01N2203/0048—Hydraulic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0064—Initiation of crack
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Earth Drilling (AREA)
Abstract
本发明公开一种基于热流固耦合的岩石起裂压力计算方法,包括根据测井资料、室内三轴压缩试验和巴西劈裂实验确定目标地层的基本参数;根据力平衡方程、本构方程、传输方程、质量守恒方程、热量守恒方程、边界条件以及最大拉应力破坏准则建立基于热流固耦合的Laplace空间的岩石起裂压力计算模型;根据目标地层的基本参数、基于热流固耦合的Laplace空间的岩石起裂压力计算模型以及Stehfest数值反演方法确定岩石起裂压力。本发明充分考虑了井筒中流体渗流、热扩散与岩石变形之间的耦合关系;能够更准确的计算地层起裂压力,为压裂施工设计提供高精度的参数,满足储层评价和压裂改造指导需要。
Description
技术领域
本发明涉及一种基于热流固耦合的岩石起裂压力计算方法,属于油气储层增产技术领域。
背景技术
压裂是油气井增产、注水井增注的一项重要技术措施。从确定施工规模、正确选择压裂设备到预测施工效果都要涉及到起裂压力,准确地预测压裂井的起裂压力是水力压裂成功实施的关键步骤,直接影响储层改造效果。因此,地层起裂压力是压裂设计和施工工艺的一项重要技术参数。
目前,传统的地层起裂压力计算方法有马修斯和凯利法、伊顿法、安德森法和黄氏计算法等,这些主要方法在二十世纪六十年代至八十年代形成。其中大多计算方法缺乏严格的理论依据,考虑因素不全,与实际结果相差较大,只适用于特定地区,并且忽略了温度变化和流体渗流对岩石破裂的影响,更未考虑井筒周围应力和孔隙压力随时间的变化。
因此,亟需建立一种基于热流固耦合的岩石起裂压力计算方法,综合考虑井筒中流体渗流、热扩散与岩石变形之间的耦合关系,确定井筒周围的最大应力集中,从而得到岩石起裂压力,为压裂施工设计提供高精度的参数,对实际钻井完井过程具有重要的指导意义。
发明内容
为了克服现有技术中的问题,本发明提供一种基于热流固耦合的岩石起裂压力计算方法。
本发明解决上述技术问题所提供的技术方案是:一种基于热流固耦合的岩石起裂压力计算方法,包括:
根据测井资料、室内三轴压缩试验和巴西劈裂实验确定目标地层的基本参数;
根据力平衡方程、本构方程、传输方程、质量守恒方程、热量守恒方程、边界条件以及最大拉应力破坏准则建立热流固耦合的Laplace空间的岩石起裂压力计算模型;
根据目标地层的基本参数、热流固耦合的Laplace空间的岩石起裂压力计算模型以及Stehfest数值反演方法确定岩石起裂压力。
进一步的技术方案是,所述力平衡方程:
σij,j+fi=0
式中:σij,j为应力张量,无量纲;fi为力的指标符号,无量纲;
本构方程:
其中:
式中:εij为固体应变张量,无量纲;σij为固体应力张量,无量纲;σkk为总主应力,MPa;p为孔隙压力,MPa;ζ为孔隙体积变化,无量纲;T为温度,K;δij为克罗内克函数;G为体积剪切模量,MPa;α为Biot系数,无量纲;v、vu分别为压力波及区泊松比和压力未波及区泊松比;αs为固体基质的体积热膨胀系数,1/K;αf为孔隙流体的体积热膨胀系数,1/K;
传输方程:
傅里叶定律:hi=-kTT,i
式中:qi为流体流量,m3/s;hi为热通量,W/m2;k为固体基质渗透率,D;μ为流体黏度,Pa·s;p,i为压力张量,无量纲;kT为饱和流体岩石的热传导系数,W/(m·K);T,i为温度张量,无量纲;
质量守恒方程:
式中:ρf为流体密度,kg/m3;φρf为流体和固体的混合密度,kg/m3;vf为流体速度m/s;
热量守恒方程:
式中:为饱和流体的熵,J/(molK);Q为流体质量通量,kg/(m2·s);ef为单位质量流体中孔隙流体的内能,J;p为孔隙压力,MPa;ηf为单位质量流体中孔隙流体的熵,J/(molK);为每单位体积的流体质量,Kg;qh为热通量,W/m2;
边界条件:
井壁处的边界条件为:
无穷远处的边界条件为:
式中:σrr为径向应力,MPa;σrθ为剪应力,MPa;pm为泥浆压力,MPa;Tm为泥浆温度,K;psh为井壁处压力,MPa;Tsh为井壁处温度,K;
最大拉应力破坏准则:
当井壁有效周向应力超过岩石的抗拉强度后,即认为岩石发生破坏,判别式如下:
σθ-αpp=-σ′t
式中:α为Biot系数,无量纲;σ′t为岩石原始抗张强度,MPa;pp为孔隙压力,MPa。
进一步的技术方案是,所述根据力平衡方程、本构方程、传输方程、质量守恒方程、热量守恒方程、边界条件以及最大拉应力破坏准则建立基于热流固耦合的Laplace空间的岩石起裂压力计算模型包括:
根据力平衡方程、本构方程、传输方程、质量守恒、热量守恒方程建立场方程;
根据场方程、边界条件建立Laplace空间井壁周围地层的孔隙压力方程和总周向应力方程;
将最大拉应力破坏准则进行Laplace变换;
根据Laplace空间井壁周围地层的总周向应力方程和孔隙压力方程、Laplace变换后的最大拉应力破坏准则建立基于热流固耦合的Laplace空间的岩石起裂压力计算模型。
进一步的技术方案是,所述场方程为:
式中:uij,j为固体基质的位移张量,无量纲;p,i为压力张量,无量纲;T,i为温度张量,无量纲;G为剪切模量,MPA;v为压力波及区泊松比;α为Biot系数,无量纲;αs为固体基质的体积热膨胀系数,1/K。
进一步的技术方案是,所述Laplace空间井壁周围地层的孔隙压力方程:
Laplace空间井壁周围地层的总周向应力方程:
其中:
D1=2(vu-v)K1(β)
D2=β(1-v)K2(β)
式中:s为与时间t相关的Laplace变换量,无量纲;Pf为地层孔隙压力,MPa;a为井眼半径,m;r为井眼轴线到地层中一点的半径,m;θ为任意径向与x轴的极角,°;ω为井斜方位与水平最大主应力的夹角,°;γ为井斜角,°;P0为各向同性压缩应力,MPa;η为常量;cf为流体扩散系数,m2/s;cT为热扩散系数,m2/s;cfT为热-流体耦合压力系数,MPa/K;Tm为泥浆温度,K;Tf为地层温度,K;pm为钻孔中的泥浆压力,MPa;pf为地层孔隙压力,MPa;αm为固体热膨胀系数,K-1;α为Biot系数,无量纲;K0、K1、K2分别为第二类零阶、一阶、二阶修正Bessel函数;ξ、ξT、β、βT、C1、C2、C3为变量,无量纲;v、vu分别为压力波及区和压力未波及区泊松比;B为Skempton孔隙压力系数,无量纲;S0为偏应力,MPa;σH、σh、σν分别为最大、最小水平主应力和垂向应力,MPa。
进一步的技术方案是,所述热流固耦合的Laplace空间的岩石起裂压力计算模型:
其中:
D1=2(vu-v)K1(β)
D2=β(1-v)K2(β)
式中:为Laplace空间下基于热流固耦合的起裂压力,MPa;Pf为地层孔隙压力,MPa;a为井眼半径,m;r为井眼轴线到地层中一点的半径,m;θ为任意径向与x轴的极角,°;P0为各向同性压缩应力,MPa;η为常量;cf为流体扩散系数,m2/s;cT为热扩散系数,m2/s;cfT为热-流体耦合压力系数,MPa/K;Tm为泥浆温度,K;Tf为地层温度,K;pm为钻孔中的泥浆压力,MPa;pf为地层孔隙压力,MPa;αm为固体热膨胀系数,K-1;α为Biot系数,无量纲;K0、K1、K2分别为第二类零阶、一阶、二阶修正Bessel函数;ξ、ξT、β、βT、C1、C2、C3为变量,无量纲;v、vu分别为压力波及区和压力未波及区泊松比;B为Skempton孔隙压力系数,无量纲;S0为偏应力,MPa;σH、σh、σν分别为最大、最小水平主应力和垂向应力,MPa。
进一步的技术方案是,所述根据目标地层的基本参数、热流固耦合的Laplace空间的岩石起裂压力计算模型以及Stehfest数值反演方法确定岩石起裂压力包括:基于热流固耦合的Laplace空间的岩石起裂压力计算模型,利用Stehfest数值反演方法,将Laplace空间下的岩石起裂压力计算模型反演到实空间,从而利用基本参数求解出岩石起裂压力。
进一步的技术方案是,所述Stehfest数值反演方法中的反演公式为:
进一步的技术方案是,所述反演公式中的N=8。
本发明具有以下有益效果:本发明充分考虑了井筒中流体渗流、热扩散与岩石变形之间的耦合关系,基于本构和传输方程与力平衡、质量和热量守恒方程,获得场方程,从而得到与时间相关的井壁周向应力和孔隙压力,最后根据最大拉应力破坏准则,得Laplace空间下起裂压力计算模型,再反演得起裂压力;能够更准确的计算地层起裂压力,为压裂施工设计提供高精度的参数,满足储层评价和压裂改造指导需要。
附图说明
图1是本发明的流程框图;
图2是实施例中的不同井起裂压力柱状图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明的一种基于热流固耦合的岩石起裂压力计算方法;包括以下步骤:
步骤S1、基于测井资料、室内三轴压缩试验和巴西劈裂实验,获得目标地层的基本参数,所述基本参数包括最大水平主应力、最小水平主应力、上覆岩层应力、岩石有效应力系数、岩石孔隙度、原始孔隙压力、原始地层温度、岩石体积热膨胀系数、岩石杨氏弹性模量、泊松比、岩石原始抗张强度等参数;
步骤S2、基于力平衡方程、本构方程、传输方程、质量守恒方程和热量守恒方程,将这五种方程结合形成场方程,并且通过流体和热量传输方程推出流体和热扩散方程;
(1)力平衡方程
σij,j+fi=0 (1)
式中:σij,j为应力张量,无量纲。fi为力的指标符号,无量纲。
(2)非等温导热条件下热孔弹性介质耦合本构方程
其中:
式中:εij为固体应变张量,无量纲;σij为固体应力张量,无量纲;σkk为总主应力,MPa;p为孔隙压力,MPa;ζ为孔隙体积变化,无量纲;T为温度,K;δij为克罗内克函数;G为体积剪切模量,MPa;α为Biot系数,无量纲;v、vu分别为压力波及区泊松比和压力未波及区泊松比;αs为固体基质的体积热膨胀系数,1/K;αf为孔隙流体的体积热膨胀系数,1/K。
(3)传输方程
忽略热对流效应和热渗流作用,可分别得到流体和热量传输方程:
傅里叶定律:hi=-kTT,i (6)
式中:qi为流体流量,m3/s;hi为热通量,W/m2;k为固体基质渗透率,D;μ为流体黏度,Pa·s;p,i为压力张量,无量纲;kT为饱和流体岩石的热传导系数,W/(m·K);T,i为温度张量,无量纲。
(4)流体在饱和多孔介质中的质量守恒方程
式中:ρf为流体密度,kg/m3;φρf为流体和固体的混合密度,kg/m3;vf为流体速度m/s。
(5)热量守恒方程
式中:为饱和流体的熵,J/(molK);Q为流体质量通量,kg/(m2·s);ef为单位质量流体中孔隙流体的内能,J;p为孔隙压力,MPa;ηf为单位质量流体中孔隙流体的熵,J/(molK);为每单位体积的流体质量,Kg;qh为热通量,W/m2。
通过(5)式和(7)式可推得流体扩散方程为:
其中:
其中:
式中:cf为流体扩散系数,m2/s;cfT为热—流体耦合压力系数,MPa/K;p,kk为应力张量,无量纲;k为固体基质的渗透率,%;G为剪切模量,MPa;B为Skempton孔隙压力系数,无量纲;v、vu分别为压力波及区和压力未波及区泊松比;φ为孔隙度,%。
通过(6)式和(8)式可推得热扩散方程为:
其中:
式中:cT为热扩散系数,m2/s;T,kk为温度张量,无量纲;ρm为总密度(流体和固体),kg/m3;cm为总比热容,J/(kg·K)。
将(1)、(2)、(3)、(5)、(6)、(7)、(8)结合得到场方程为:
式中:uij,j为固体基质的位移张量,无量纲;p,i为压力张量,无量纲;T,i为温度张量,无量纲;
步骤S3、基于步骤S2得到的场方程,将边界条件带入场方程中,得到Laplace空间的总周向应力和孔隙压力。
将边界条件带入场方程,得到Laplace空间的孔隙压力和总周向应力。
井壁承受不同的压力和温度,适用以下的边界条件。
井壁处(r=a)的边界条件为:
无穷远处(r→∞)的边界条件为:
式中:σrr为径向应力,MPa;σrθ为剪应力,MPa;pm为泥浆压力,MPa;Tm为泥浆温度,K;psh为井壁处压力,MPa;Tsh为井壁处温度,K。
将(17)、(18)式带入(16)式中,得到的Laplace空间的孔隙压力为:
Laplace空间的总周向应力为:
其中:
D1=2(vu-v)K1(β) (32)
D2=β(1-v)K2(β) (33)
式中:s为与时间t相关的Laplace变换量,无量纲;Pf为地层孔隙压力,MPa;a为井眼半径,m;r为井眼轴线到地层中一点的半径,m;θ为任意径向与x轴的极角,°;ω为井斜方位与水平最大主应力的夹角,°;γ为井斜角,°;P0为各向同性压缩应力,MPa;η为常量;cf为流体扩散系数,m2/s;cT为热扩散系数,m2/s;cfT为热-流体耦合压力系数,MPa/K;Tm为泥浆温度,K;Tf为地层温度,K;pm为钻孔中的泥浆压力,MPa;pf为地层孔隙压力,MPa;αm为固体热膨胀系数,K-1;α为Biot系数,无量纲;K0、K1、K2分别为第二类零阶、一阶、二阶修正Bessel函数;ξ、ξT、β、βT、C1、C2、C3为变量,无量纲;ν、νu分别为压力波及区和压力未波及区泊松比;B为Skempton孔隙压力系数,无量纲;S0为偏应力,MPa;σH、σh、σν分别为最大、最小水平主应力和垂向应力,MPa。
根据应力坐标变换,上述公式可应用到任何斜度的井。斜井中地应力(σv,σH,σh)被转换成新的坐标,如下所示:
式中:σxx、σyy、σzz分别为坐标系(x,y,z)中正应力分量,MPa;σxy、σxz、σyz分别为坐标系(x,y,z)中剪应力分量,MPa。
步骤S4、基于步骤S3得到的Laplace空间的总周向应力和孔隙压力,将其带入最大拉应力破坏准则,得到基于热流固耦合的Laplace空间的岩石起裂压力计算模型;
破裂发生在井壁周向应力最小处,即θ=0°或θ=180°。根据最大拉应力破坏准则,当井壁有效周向应力超过岩石的抗拉强度后,即认为岩石发生破坏,判别式如下:
σθ-αpp=-σ′t (35)
式中:α为Biot系数,无量纲;σt为岩石原始抗张强度,MPa。pp为孔隙压力,MPa。
将(35)式进行Laplace变换得:
当θ=0°或θ=180°,r=a时,将(19)式和(20)式带入(36)式中,即可得到的Laplace空间下基于热流固耦合的岩石起裂压力计算模型为:
步骤S5、基于步骤S4得到的基于热流固耦合的Laplace空间的岩石起裂压力计算模型,利用Stehfest数值反演方法,将Laplace空间下的岩石起裂压力计算模式反演到实空间,从而求解起裂压力。
反演公式为(38)式~(40)式:
根据(37)式、(38)式、(39)式、(40)式可以求解得实空间的岩石起裂压力。
本实例实施中,剪切模量G=760MPa,Biot系数=0.966,压力波及区泊松比v=0.219,压力不波及区泊松比vu=0.461,渗透率k=1×10-20m2,孔隙度φ=20%,固体热膨胀系数αm=1.8×10-5K-1,流体热膨胀系数αf=3.0×10-4K-1,热扩散系数cT=1.6×10-4K-1,Skempton孔隙压力系数B=0.915,流体黏度μ=3.0×10-4mPa·s,流体扩散系数cf=6.0×10-8m2/s,计算时间为120s。
根据不同地层岩石参数及测试数据,起裂压力计算结果如表1所示。
表1
根据表1画出不同井起裂压力柱状图,如图2所示。
为更好地说明预测结果的准确性,将本发明与黄氏模型、Stephen模型和Eaton模型进行对比,如图2所示,结果显示:本发明计算结果与实际地层起裂压力较接近,其他方法的预测结果则相对偏低。统计得出,黄氏模型的误差为38%,Stephen模型的误差为28%,Eaton模型的误差为28.7%,本发明模型的误差为4.29%。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (7)
1.一种基于热流固耦合的岩石起裂压力计算方法,其特征在于,包括:
根据测井资料、室内三轴压缩试验和巴西劈裂实验确定目标地层的基本参数;
根据力平衡方程、本构方程、传输方程、质量守恒方程、热量守恒方程、边界条件以及最大拉应力破坏准则建立基于热流固耦合的Laplace空间的岩石起裂压力计算模型;
建立基于热流固耦合的Laplace空间的岩石起裂压力计算模型包括:
根据力平衡方程、本构方程、传输方程、质量守恒、热量守恒方程建立场方程;
根据场方程、边界条件建立Laplace空间井壁周围地层的孔隙压力方程和总周向应力方程;
所述Laplace空间井壁周围地层的孔隙压力方程:
Laplace空间井壁周围地层的总周向应力方程:
其中:
D1=2(vu-v)K1(β)
D2=β(1-v)K2(β)
式中:s为与时间t相关的Laplace变换量,无量纲;Pf为地层孔隙压力,MPa;a为井眼半径,m;r为井眼轴线到地层中一点的半径,m;θ为任意径向与x轴的极角,°;ω为井斜方位与水平最大主应力的夹角,°;γ为井斜角,°;P0为各向同性压缩应力,MPa;η为常量;cf为流体扩散系数,m2/s;cT为热扩散系数,m2/s;cfT为热-流体耦合压力系数,MPa/K;Tm为泥浆温度,K;Tf为地层温度,K;pm为钻孔中的泥浆压力,MPa;pf为地层孔隙压力,MPa;αm为固体热膨胀系数,K-1;α为Biot系数,无量纲;K0、K1、K2分别为第二类零阶、一阶、二阶修正Bessel函数;ξ、ξT、β、βT、C1、C2、C3为变量,无量纲;v、vu分别为压力波及区和压力未波及区泊松比;B为Skempton孔隙压力系数,无量纲;S0为偏应力,MPa;σH、σh、σν分别为最大、最小水平主应力和垂向应力,MPa;
将最大拉应力破坏准则进行Laplace变换;
根据Laplace空间井壁周围地层的总周向应力方程和孔隙压力方程、Laplace变换后的最大拉应力破坏准则建立基于热流固耦合的Laplace空间的岩石起裂压力计算模型;
根据目标地层的基本参数、基于热流固耦合的Laplace空间的岩石起裂压力计算模型以及Stehfest数值反演方法确定岩石起裂压力。
2.根据权利要求1所述的一种基于热流固耦合的岩石起裂压力计算方法,其特征在于,所述力平衡方程:
σij,j+fi=0
式中:σij,j为应力张量,无量纲;fi为力的指标符号,无量纲;
本构方程:
其中:
式中:εij为固体应变张量,无量纲;σij为固体应力张量,无量纲;σkk为总主应力,MPa;p为孔隙压力,MPa;ζ为孔隙体积变化,无量纲;T为温度,K;δij为克罗内克函数;G为体积剪切模量,MPa;α为Biot系数,无量纲;v、vu分别为压力波及区泊松比和压力未波及区泊松比;αs为固体基质的体积热膨胀系数,1/K;αf为孔隙流体的体积热膨胀系数,1/K;
传输方程:
傅里叶定律:hi=-kTT,i
式中:qi为流体流量,m3/s;hi为热通量,W/m2;k为固体基质渗透率,D;μ为流体黏度,Pa·s;p,i为压力张量,无量纲;kT为饱和流体岩石的热传导系数,W/(m·K);T,i为温度张量,无量纲;
质量守恒方程:
式中:ρf为流体密度,kg/m3;φρf为流体和固体的混合密度,kg/m3;vf为流体速度m/s;
热量守恒方程:
式中:为饱和流体的熵,J/(molK);Q为流体质量通量,kg/(m2·s);ef为单位质量流体中孔隙流体的内能,J;p为孔隙压力,MPa;ηf为单位质量流体中孔隙流体的熵,J/(mol K);为每单位体积的流体质量,Kg;qh为热通量,W/m2;
边界条件:
井壁处的边界条件为:
无穷远处的边界条件为:
式中:σrr为径向应力,MPa;σrθ为剪应力,MPa;pm为泥浆压力,MPa;Tm为泥浆温度,K;psh为井壁处压力,MPa;Tsh为井壁处温度,K;
最大拉应力破坏准则:
当井壁有效周向应力超过岩石的抗拉强度后,即认为岩石发生破坏,判别式如下:
σθ-αpp=-σ′t
式中:α为Biot系数,无量纲;σ′t为岩石原始抗张强度,MPa;pp为孔隙压力,MPa。
4.根据权利要求1所述的一种基于热流固耦合的岩石起裂压力计算方法,其特征在于,所述基于热流固耦合的Laplace空间的岩石起裂压力计算模型:
其中:
D1=2(vu-v)K1(β)
D2=β(1-v)K2(β)
式中:为Laplace空间下基于热流固耦合的起裂压力,MPa;Pf为地层孔隙压力,MPa;a为井眼半径,m;r为井眼轴线到地层中一点的半径,m;θ为任意径向与x轴的极角,°;P0为各向同性压缩应力,MPa;η为常量;cf为流体扩散系数,m2/s;cT为热扩散系数,m2/s;cfT为热-流体耦合压力系数,MPa/K;Tm为泥浆温度,K;Tf为地层温度,K;pm为钻孔中的泥浆压力,MPa;pf为地层孔隙压力,MPa;αm为固体热膨胀系数,K-1;α为Biot系数,无量纲;K0、K1、K2分别为第二类零阶、一阶、二阶修正Bessel函数;ξ、ξT、β、βT、C1、C2、C3为变量,无量纲;v、vu分别为压力波及区和压力未波及区泊松比;B为Skempton孔隙压力系数,无量纲;S0为偏应力,MPa;σH、σh、σν分别为最大、最小水平主应力和垂向应力,MPa。
5.根据权利要求1所述的一种基于热流固耦合的岩石起裂压力计算方法,其特征在于,所述根据目标地层的基本参数、热流固耦合的Laplace空间的岩石起裂压力计算模型以及Stehfest数值反演方法确定岩石起裂压力包括:基于热流固耦合的Laplace空间的岩石起裂压力计算模型,利用Stehfest数值反演方法,将Laplace空间下的岩石起裂压力计算模型反演到实空间,从而利用基本参数求解出岩石起裂压力。
7.根据权利要求6所述的一种基于热流固耦合的岩石起裂压力计算方法,其特征在于,所述反演公式中的N=8。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110394846.8A CN113109162B (zh) | 2021-04-13 | 2021-04-13 | 一种基于热流固耦合的岩石起裂压力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110394846.8A CN113109162B (zh) | 2021-04-13 | 2021-04-13 | 一种基于热流固耦合的岩石起裂压力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113109162A CN113109162A (zh) | 2021-07-13 |
CN113109162B true CN113109162B (zh) | 2022-02-01 |
Family
ID=76716422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110394846.8A Active CN113109162B (zh) | 2021-04-13 | 2021-04-13 | 一种基于热流固耦合的岩石起裂压力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113109162B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113935092B (zh) * | 2021-10-23 | 2022-08-23 | 西南石油大学 | 一种裂缝性地层井壁稳定双重孔隙热弹性耦合评价方法 |
CN116151139B (zh) * | 2022-09-14 | 2024-01-30 | 清华大学 | 一种深海钻井时井筒稳定性的预测方法、装置及设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102135478A (zh) * | 2011-01-07 | 2011-07-27 | 中国海洋石油总公司 | 一种测定天然气水合物沉积物体变的三轴试验装置 |
CN103344539A (zh) * | 2013-06-27 | 2013-10-09 | 中国科学院武汉岩土力学研究所 | 渗透压力和相分比例可控的岩石长期多场耦合实验装置与测试方法 |
CN105352996A (zh) * | 2015-10-20 | 2016-02-24 | 中国矿业大学 | 一种测试地下煤田火区覆岩温度变化的模型试验方法 |
CN108830020A (zh) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | 一种基于热流固耦合理论的模拟海上油田微压裂增注裂缝扩展的方法 |
CN109975119A (zh) * | 2019-04-15 | 2019-07-05 | 四川大学 | 一种岩石双轴压缩***设计方法 |
CN110005389A (zh) * | 2019-03-07 | 2019-07-12 | 西南石油大学 | 一种基于热流固耦合作用的超深层砂岩缝网改造评价方法 |
CN111488703A (zh) * | 2019-01-28 | 2020-08-04 | 西南石油大学 | 基于扩展有限元模拟层状岩石断裂韧性的方法 |
CN112177598A (zh) * | 2019-07-01 | 2021-01-05 | 中国石油化工股份有限公司 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9080946B2 (en) * | 2012-06-15 | 2015-07-14 | Ingrain, Inc. | Digital rock analysis systems and methods with multiphase flow REV determination |
CN102735549B (zh) * | 2012-07-05 | 2014-07-02 | 重庆大学 | 多功能真三轴流固耦合压力室 |
CN107609258B (zh) * | 2017-09-07 | 2018-07-27 | 西南石油大学 | 一种页岩重复压裂转向裂缝起裂压力的计算方法 |
CN109063383A (zh) * | 2018-09-19 | 2018-12-21 | 西南石油大学 | 基于微尺度重建模型的热-流-固多场耦合模拟方法 |
-
2021
- 2021-04-13 CN CN202110394846.8A patent/CN113109162B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102135478A (zh) * | 2011-01-07 | 2011-07-27 | 中国海洋石油总公司 | 一种测定天然气水合物沉积物体变的三轴试验装置 |
CN103344539A (zh) * | 2013-06-27 | 2013-10-09 | 中国科学院武汉岩土力学研究所 | 渗透压力和相分比例可控的岩石长期多场耦合实验装置与测试方法 |
CN105352996A (zh) * | 2015-10-20 | 2016-02-24 | 中国矿业大学 | 一种测试地下煤田火区覆岩温度变化的模型试验方法 |
CN108830020A (zh) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | 一种基于热流固耦合理论的模拟海上油田微压裂增注裂缝扩展的方法 |
CN111488703A (zh) * | 2019-01-28 | 2020-08-04 | 西南石油大学 | 基于扩展有限元模拟层状岩石断裂韧性的方法 |
CN110005389A (zh) * | 2019-03-07 | 2019-07-12 | 西南石油大学 | 一种基于热流固耦合作用的超深层砂岩缝网改造评价方法 |
CN109975119A (zh) * | 2019-04-15 | 2019-07-05 | 四川大学 | 一种岩石双轴压缩***设计方法 |
CN112177598A (zh) * | 2019-07-01 | 2021-01-05 | 中国石油化工股份有限公司 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
Non-Patent Citations (4)
Title |
---|
A lined hole in a viscoelastic rock under biaxial far-field stress;sofia;《international journal of rock mechanics and mining sciences》;20181231;第106卷;350-363 * |
关于Laplace 数值反演Stehfest算法适用性的一点思考;何家欢;《油气井测试》;20150831;第18-20页 * |
多孔岩体热流固耦合有限元模拟;张强林;《三峡大学学报》;20080630;第30卷(第3期);第18-21页 * |
岩石损伤过程中的热-流-力耦合模型及其引用初探;朱万成;《岩土力学》;20091231;第3851-3857页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113109162A (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113109162B (zh) | 一种基于热流固耦合的岩石起裂压力计算方法 | |
CN108952700B (zh) | 一种各向异性地层井壁破裂压力确定方法 | |
CN103206206B (zh) | 用于计算超深井地层破裂压力的方法 | |
CN109356567B (zh) | 深水浅部地层井壁稳定性预测方法 | |
CN103912269B (zh) | 页岩气储层地层破裂压力梯度录井确定方法 | |
US10928282B2 (en) | Method and device for determining elasticity of cement stone utilized in well cementing of oil-gas well | |
CN102182453B (zh) | 井壁坍塌分析方法 | |
CN107169248B (zh) | 一种特殊地层安全泥浆密度窗口确定方法 | |
CN108843303B (zh) | 一种基于泥岩蠕变模型的油水井套管损坏预测方法 | |
CN111353205B (zh) | 用于致密气藏产水气井地层压力和动态产能的计算方法 | |
CN109033504B (zh) | 一种油水井套管损坏预测方法 | |
CN114372352A (zh) | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 | |
CN111927417A (zh) | 一种页岩气分段压裂水平井组储量动用状况评价方法 | |
CN112302607B (zh) | 一种致密气藏压裂水平井人工裂缝参数解释方法 | |
CN114595504B (zh) | 一种页岩气藏压裂水平井产出剖面反演解释方法 | |
CN111963164A (zh) | 一种针对多裂缝发育储层的井壁坍塌压力评价方法 | |
CN113189647B (zh) | 一种横观各向同性页岩地层脆性指数预测方法 | |
CN116205163A (zh) | 一种天然气水合物藏数值模拟中的地层稳定性分析方法 | |
CN115293066B (zh) | 考虑地层渗流传热效应的气井温度场计算方法 | |
CN116717224A (zh) | 一种针对低渗致密储层复杂缝网的压裂产能预测方法 | |
CN109522579B (zh) | 水平井压裂施工破裂压力预测方法 | |
CN113530534B (zh) | 储层改造方法及装置 | |
CN112651113A (zh) | 适用于东海压力衰竭地层的安全泥浆密度窗口判定方法 | |
CN116341174A (zh) | 一种裂缝性砂岩储层临界出砂生产压差预测方法及*** | |
CN115186524A (zh) | 低渗透长缝压裂井组co2驱临界注气参数的评测方法 |
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 |