CN112257230A - 一种水平井套管交变压力作用下的套管力学分析计算方法 - Google Patents

一种水平井套管交变压力作用下的套管力学分析计算方法 Download PDF

Info

Publication number
CN112257230A
CN112257230A CN202010996834.8A CN202010996834A CN112257230A CN 112257230 A CN112257230 A CN 112257230A CN 202010996834 A CN202010996834 A CN 202010996834A CN 112257230 A CN112257230 A CN 112257230A
Authority
CN
China
Prior art keywords
casing
sleeve
stress
plastic
elastic
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.)
Pending
Application number
CN202010996834.8A
Other languages
English (en)
Inventor
段志锋
欧阳勇
吴学升
来轩昂
黄占盈
高云文
蒙华军
陈春宇
赵巍
张燕娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010996834.8A priority Critical patent/CN112257230A/zh
Publication of CN112257230A publication Critical patent/CN112257230A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开一种水平井套管交变压力作用下的套管力学分析计算方法,包括以下步骤:将受到交变压力的套管分为塑性区域和弹性区域,套管内交变压力f(pi)大于套管产生塑性变形的最大应力时,则会在套管内壁处开始产生塑性变形,随着套管内压力的加大,弹塑性边界也会慢慢的向套管外壁移动,用rρ1来表示套管弹塑性边界:当r1≤r≤rρ1时,分析得出套管塑性区内的应力分量;当r≥rρ1时,给出弹性区域的应力情况和套管弹性区外径的径向位移的计算方法。本发明用于解决现有的水平井压裂过程中,当套管内存在交变压力时套管力学分析不明确、计算方法不确定等问题。为水平井进行交变压力的压裂时套管受力分析提供参考,对保持井筒完整性具有重要意义。

Description

一种水平井套管交变压力作用下的套管力学分析计算方法
技术领域
本发明涉及石油与天然气工程技术领域,特别涉及一种水平井套管交变压力作用下的套管力学分析计算方法。
背景技术
随着石油天然气勘探开发工作的不断深入,页岩气、致密油气以及储气库建设工作的开展,油井面临的井下的地质条件也越来越昔刻;井身结构变得越来越复杂,更容易发生水泥环密封失效的问题。而套管内压的变化是导致水泥环与套管脱胶的最主要问题之一。如何能够准确地计算出在进行交变压力施工作业时水泥环与套管第一胶结面压力大小是目前很难攻克的一个问题。针对此方面的研究,研究者们在理论和数值模型方面做了大量的工作,并取得了一定的结果。但是在针对水平井段套管在交变压力问题下,压裂变化导致的第一胶结面受力的研究依然还存在着很多的欠缺,对水平井段第一胶结面受力分析还需要进一步的提高。因此,能够针对不同套管内压,进行套管受力分析计算就具有重要的研究意义。
在开展水平井套管交变压力作用下的套管力学分析的过程中,不仅需要考虑施工作业时的交变压力,还要分析套管的受力状态。为此,本发明内容重点阐述与保护的内容是在进行交变压力作用过程中,设计一种水平井套管力学分析的计算方法。
发明内容
本发明的目的是给出一种水平井套管交变压力作用下的套管力学分析计算方法,并且能通过分析得到套管力学分布和径向位移。它主要用于解决现有的水平井压裂过程中,当套管内存在交变压力时套管力学分析不明确、计算方法不确定等问题。为水平井进行交变压力的压裂时套管受力分析提供参考,对保持井筒完整性具有重要意义。
为达到上述目的,本发明采用如下技术方案:
一种水平井套管交变压力作用下的套管力学分析计算方法,包括以下步骤:
将受到交变压力的套管分为塑性区域和弹性区域,套管内交变压力f(pi)大于套管产生塑性变形的最大应力时,则会在套管内壁处开始产生塑性变形,随着套管内压力的加大,弹塑性边界也会慢慢的向套管外壁移动,用rρ1来表示套管弹塑性边界:当r1≤r≤rρ1时,分析得出套管塑性区内的应力分量;当r≥rρ1时,给出弹性区域的应力情况和套管弹性区外径的径向位移的计算方法。
作为本发明的进一步改进,塑性区域和弹性区域的应力解由平衡方程与应力方程为依据推出,套管弹性区外径的径向位移需满足边界位置Tresca准则。
作为本发明的进一步改进,套管塑性区域的受力分析为:
当r1≤r≤rρ1时,套管处于塑性状态,平衡方程变为:
Figure BDA0002692812780000021
利用边界条件来解式(1)的微分方程得到:
Figure BDA0002692812780000022
式(2)即为套管塑性区应力分量。
作为本发明的进一步改进,套管弹性区域的受力分析如下:
当r≥rρ1时,此时套管处在弹性区域中,应用弹性应力与应变的关系得:
Figure BDA0002692812780000031
式中:A和B为任意常数;
在当弹性区域的应力状态必须满足r=r2时的边界条件,式(2)在r=rρ1处的连续条件为:
Figure BDA0002692812780000032
式中:上标(a)为塑性区内的应力分量,上标(b)为弹性区内的应力分量;
运用弹性区域的边界条件
Figure BDA0002692812780000033
公式(4)中的两个常数以确定为如下形式:
Figure BDA0002692812780000034
弹性区域在交变压力f(pi)的作用下在弹塑性的边界位置同样满足Tresca准则即
Figure BDA0002692812780000035
将上式简化为:
Figure BDA0002692812780000036
式(6)给出了弹性区域的应力解,当r=r2的时候得到:
Figure BDA0002692812780000037
式中,q2为套管塑性半径上的应力;σ0为原始水平地应力;rρ1为套管产生塑性变形的半径;rρ2为水泥环产生塑性变形的半径;q1为套管外壁与水泥环内径的层间压力。
作为本发明的进一步改进,已知弹性区域的应力情况,求得套管弹性区外径的径向位移:
Figure BDA0002692812780000041
式中,
Figure BDA0002692812780000042
为套管弹性区外径的径向位移;μ1为套管泊松比;E1为套管弹性模量;r2为套管外半径。
作为本发明的进一步改进,计算方法中的计算过程中进行了以下假设:
1)套管为弹塑性各向同性材料;
2)原场地应力为水平均布应力;
3)套管、水泥环以及地层三者是紧密连接的,无滑动,满足应力与位移的连续条件;
4)井眼垂直且套管居中。
与现有技术相比,本发明具有以下优点:
本发明在交变压力作用下,针对水平井套管力学问题而设计的一种力学分析计算方法主要由套管的塑性区域和弹性区域的受力分析两部分组成。计算表达式主要由套管塑性区内的应力分量、弹性区域的应力解、套管弹性区外径的径向位移组成。解决了现阶段套管施加交变压力时,套管受力不明确的问题。
附图说明
图1是本发明组合系统弹塑性分析图。
图中:套管内半径为r1,m;套管外半径(即水泥环内半径)为r2,m;井眼半径(即水泥环外半径)为r3,m;地层岩石圈半径(实际为无限大)为r4,m;套管产生塑性变形的半径为rρ1,m;水泥环产生塑性变形的半径为rρ2,m;套管内压为P1,MPa;套管弹性模量为E1,GPa;套管泊松比为μ1;水泥环弹性模量为E2,GPa;水泥环泊松比为μ2;地层弹性模量为E3,GPa;地层泊松比为μ3;图中的q1表示为套管外壁与水泥环内径的层间压力,MPa;q4为水泥环外壁与井眼内径的层间压力,MPa;q2为套管塑性半径上的应力,MPa;q3为水泥环塑性半径上的应力,MPa;q5为地层塑性半径上的应力,MPa。
图2是本发明套管弹塑性区域示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
一种水平井套管交变压力作用下的套管力学分析计算方法,包括以下步骤:
将受到交变压力的套管分为塑性区域和弹性区域,套管内交变压力f(pi)大于套管产生塑性变形的最大应力时,则会在套管内壁处开始产生塑性变形,随着套管内压力的加大,弹塑性边界也会慢慢的向套管外壁移动,用rρ1来表示套管弹塑性边界:当r1≤r≤rρ1时,分析得出套管塑性区内的应力分量;当r≥rρ1时,给出弹性区域的应力情况和套管弹性区外径的径向位移的计算方法。
图1为地层-水泥环-套管组合系统弹塑性分析图;图2为从组合体中单独将套管部分取出来的套管弹塑性区域示意图。
根据材料的自身特性可以知道,由于套管的弹性模量远远大于水泥环以及地层的弹性模量,所以在套管试压过程中套管首先承担大部分内压力,水泥环及地层承担着远远小于套管的内压力。根据这一过程我们需要对该系统组合做出如下的假设:
1)套管为弹塑性各向同性材料;
2)原场地应力为水平均布应力;
3)套管、水泥环以及地层三者是紧密连接的,无滑动,满足应力与位移的连续条件;
4)井眼垂直且套管居中良好。
当套管收到交变压力时,套管的套管弹塑性边界用rρ1来表示。当r1≤r≤rρ1时,分析得出套管塑性区内的应力分量;当r≥rρ1时,给出弹性区域的应力情况和套管弹性区外径的径向位移的计算方法。
实施例
套管受力问题是一个轴对称问题,所以所有的剪应力和剪应变均为零,同时可以知道σz=0。平衡方程和应变位移关系用下式表示为:
Figure BDA0002692812780000061
这些控制方程要利用下面的边界条件求解:
Figure BDA0002692812780000062
由于考虑的是套管内压力不断增高的情况,所以此时σr≤0且σθ>0,可以观察得到σθ>σz=0≥σr。因此,在这种情况下的Tresca准则由如下形式表达:
σθr=σ0 (3)
套管内交变压力f(pi)大于套管产生塑性变形的最大应力时,则会在套管内壁处开始产生塑性变形,随着套管内压力的加大,弹塑性边界也会慢慢的向套管外壁移动。用rρ1来表示套管弹塑性边界。
(1)套管塑性区域受力分析
当r1≤r≤rρ1时,套管处于塑性状态。平衡方程变为:
Figure BDA0002692812780000071
利用式(2)所示边界条件来解式(4)的微分方程得到:
Figure BDA0002692812780000072
式(5)即为套管塑性区内的应力分量。
(2)套管弹性区域受力分析
当r≥rρ1时,此时套管处在弹性区域中,应用弹性应力与应变的关系可得:
Figure BDA0002692812780000073
式中:X和Y为任意常数。
在当弹性区域的应力状态必须满足r=r2时的边界条件,式(5)在r=rρ1处的连续条件为:
Figure BDA0002692812780000074
式中:上标(a)为塑性区内的应力分量,上标(b)为弹性区内的应力分量。
运用弹性区域的边界条件
Figure BDA0002692812780000075
公式(6)中的两个常数可以确定为如下形式:
Figure BDA0002692812780000076
弹性区域在压力pi的作用下在弹塑性的边界位置同样满足Tresca准则即
Figure BDA0002692812780000081
可将上式简化为:
Figure BDA0002692812780000082
式(9)给出了弹性区域的应力解,当r=r2的时候得到:
Figure BDA0002692812780000083
式中,q2为套管塑性半径上的应力,MPa;σ0为原始水平地应力,MPa;rρ1为套管产生塑性变形的半径,m;rρ2为水泥环产生塑性变形的半径,m;q1为套管外壁与水泥环内径的层间压力,MPa。
由于知道了弹性区域的应力情况,可以求得套管弹性区外径的径向位移:
Figure BDA0002692812780000084
式(8)中,
Figure BDA0002692812780000085
为套管弹性区外径的径向位移;μ1为套管泊松比;E1为套管弹性模量;r2为套管外半径(即水泥环内半径)。
以上的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,包括以下步骤:
将受到交变压力的套管分为塑性区域和弹性区域,套管内交变压力f(pi)大于套管产生塑性变形的最大应力时,则会在套管内壁处开始产生塑性变形,随着套管内压力的加大,弹塑性边界向套管外壁移动,用rρ1来表示套管弹塑性边界:当r1≤r≤rρ1时,分析得出套管塑性区内的应力分量;当r≥rρ1时,给出弹性区域的应力情况和套管弹性区外径的径向位移的计算方法。
2.根据权利要求1所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,塑性区域和弹性区域的应力解由平衡方程与应力方程为依据推出,套管弹性区外径的径向位移需满足边界位置Tresca准则。
3.根据权利要求2所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,
套管塑性区域的受力分析为:
当r1≤r≤rρ1时,套管处于塑性状态,平衡方程变为:
Figure FDA0002692812770000011
利用边界条件来解式(1)的微分方程得到:
Figure FDA0002692812770000012
式(2)即为套管塑性区应力分量。
4.根据权利要求2所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,
套管弹性区域的受力分析如下:
当r≥rρ1时,此时套管处在弹性区域中,应用弹性应力与应变的关系得:
Figure FDA0002692812770000021
式中:A和B为任意常数;
在当弹性区域的应力状态必须满足r=r2时的边界条件,式(2)在r=rρ1处的连续条件为:
Figure FDA0002692812770000022
式中:上标(a)为塑性区内的应力分量,上标(b)为弹性区内的应力分量;
运用弹性区域的边界条件
Figure FDA0002692812770000023
公式(4)中的两个常数以确定为如下形式:
Figure FDA0002692812770000024
弹性区域在交变压力f(pi)的作用下在弹塑性的边界位置同样满足Tresca准则即
Figure FDA0002692812770000025
将上式简化为:
Figure FDA0002692812770000026
式(6)给出了弹性区域的应力解,当r=r2的时候得到:
Figure FDA0002692812770000027
式中,q2为套管塑性半径上的应力;σ0为原始水平地应力;rρ1为套管产生塑性变形的半径;rρ2为水泥环产生塑性变形的半径;q1为套管外壁与水泥环内径的层间压力。
5.根据权利要求2所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,
已知弹性区域的应力情况,求得套管弹性区外径的径向位移:
Figure FDA0002692812770000031
式中,
Figure FDA0002692812770000032
为套管弹性区外径的径向位移;μ1为套管泊松比;E1为套管弹性模量;r2为套管外半径。
6.根据权利要求2所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,计算方法中的计算过程中进行了以下假设:
1)套管为弹塑性各向同性材料;
2)原场地应力为水平均布应力;
3)套管、水泥环以及地层三者是紧密连接的,无滑动,满足应力与位移的连续条件;
4)井眼垂直且套管居中。
CN202010996834.8A 2020-09-21 2020-09-21 一种水平井套管交变压力作用下的套管力学分析计算方法 Pending CN112257230A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010996834.8A CN112257230A (zh) 2020-09-21 2020-09-21 一种水平井套管交变压力作用下的套管力学分析计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010996834.8A CN112257230A (zh) 2020-09-21 2020-09-21 一种水平井套管交变压力作用下的套管力学分析计算方法

Publications (1)

Publication Number Publication Date
CN112257230A true CN112257230A (zh) 2021-01-22

Family

ID=74232850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010996834.8A Pending CN112257230A (zh) 2020-09-21 2020-09-21 一种水平井套管交变压力作用下的套管力学分析计算方法

Country Status (1)

Country Link
CN (1) CN112257230A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115012861A (zh) * 2022-06-28 2022-09-06 湖南科技大学 适用于海底钻机保压取心技术的压力损失程度预测方法
CN117235839A (zh) * 2022-09-16 2023-12-15 中国葛洲坝集团三峡建设工程有限公司 基于自平衡应力场分析的井筒水泥环安全载荷计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109165434A (zh) * 2018-08-13 2019-01-08 中国科学院武汉岩土力学研究所 高应力地下洞室流变性软岩稳定性分析的解析计算方法
CN109522588A (zh) * 2018-09-27 2019-03-26 天津大学 一种油井水泥环的水-热-力学多场耦合模型的建立及计算方法
CN109657358A (zh) * 2018-12-21 2019-04-19 中国地质大学(北京) 考虑围岩与支护结构相互作用的圆形隧道力学计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109165434A (zh) * 2018-08-13 2019-01-08 中国科学院武汉岩土力学研究所 高应力地下洞室流变性软岩稳定性分析的解析计算方法
CN109522588A (zh) * 2018-09-27 2019-03-26 天津大学 一种油井水泥环的水-热-力学多场耦合模型的建立及计算方法
CN109657358A (zh) * 2018-12-21 2019-04-19 中国地质大学(北京) 考虑围岩与支护结构相互作用的圆形隧道力学计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
吴健: "固井水泥环密封完整性研究", 中国优秀硕士学位论文全文数据库工程科技Ⅰ辑, no. 02, pages 019 - 60 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115012861A (zh) * 2022-06-28 2022-09-06 湖南科技大学 适用于海底钻机保压取心技术的压力损失程度预测方法
CN115012861B (zh) * 2022-06-28 2023-08-15 湖南科技大学 适用于海底钻机保压取心技术的压力损失程度预测方法
CN117235839A (zh) * 2022-09-16 2023-12-15 中国葛洲坝集团三峡建设工程有限公司 基于自平衡应力场分析的井筒水泥环安全载荷计算方法
CN117235839B (zh) * 2022-09-16 2024-04-19 中国葛洲坝集团三峡建设工程有限公司 基于自平衡应力场分析的井筒水泥环安全载荷计算方法

Similar Documents

Publication Publication Date Title
CN108868748B (zh) 一种页岩气水平井重复压裂裂缝开启压力的计算方法
Liu et al. Analysis on integrity of cement sheath in the vertical section of wells during hydraulic fracturing
Liu et al. Effect of local loads on shale gas well integrity during hydraulic fracturing process
Lecampion et al. Interface debonding as a controlling mechanism for loss of well integrity: Importance for CO2 injector wells
CN109538185B (zh) 一种耦合温度场条件下的多层套管井井筒完整性分析模型
US20030163257A1 (en) Method for selection of cementing composition
EP2583214B1 (en) Systems and methods for wellbore optimization
WO2017121112A1 (zh) 一种多级水力喷射压裂管柱性能评价方法
CN102182453B (zh) 井壁坍塌分析方法
Zhao et al. Experimental investigation on the basic law of directional hydraulic fracturing controlled by dense linear multi-hole drilling
CN112257230A (zh) 一种水平井套管交变压力作用下的套管力学分析计算方法
Meng et al. Numerical investigation of casing shear deformation due to fracture/fault slip during hydraulic fracturing
Tian et al. Problems in the wellbore integrity of a shale gas horizontal well and corresponding countermeasures
CN110543745A (zh) 一种热采井水泥环封固完整性的评价方法
Yang et al. Effects of in-situ stress and elastic parameters of cement sheath in salt rock formation of underground gas storage on seal integrity of cement sheath
CN114896775A (zh) 一种评价固井水泥环在压裂工况下层间封隔的方法
Zhang et al. Practice and understanding of sidetracking horizontal drilling in old wells in Sulige Gas Field, NW China
Liu et al. Effect of expandable cement on increasing sealing ability of cement sheath in shale gas wells
Liu et al. Study on optimization design of permanent packer slip structure
Huang et al. Local-integral coupling model of tubular strings with connectors and its application in periodic sticking analyses
Wei et al. Numerical investigation of the factors affecting the cement sheath integrity in hydraulically fractured wells
Zhao et al. Anti-channeling cementing technology for long horizontal sections of shale gas wells
Fan et al. Study on the sealing integrity of cement sheath during volume fracturing of shale gas of horizontal well
CN111859701A (zh) 一种固井水泥浆塑性体积收缩与孔隙压降评价方法
Shi et al. Research on casing damage mechanism and prevention and control under creep behavior of salt-gypsum layer

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