CN112257230A - 一种水平井套管交变压力作用下的套管力学分析计算方法 - Google Patents
一种水平井套管交变压力作用下的套管力学分析计算方法 Download PDFInfo
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 26
- 238000004364 calculation method Methods 0.000 title claims abstract description 21
- 230000009471 action Effects 0.000 title claims abstract description 16
- 238000006073 displacement reaction Methods 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 19
- 230000008569 process Effects 0.000 claims abstract description 8
- 239000004568 cement Substances 0.000 claims description 24
- 239000011229 interlayer Substances 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 2
- 230000006872 improvement Effects 0.000 description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 239000003209 petroleum derivative Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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
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时,套管处于塑性状态,平衡方程变为:
利用边界条件来解式(1)的微分方程得到:
式(2)即为套管塑性区应力分量。
作为本发明的进一步改进,套管弹性区域的受力分析如下:
当r≥rρ1时,此时套管处在弹性区域中,应用弹性应力与应变的关系得:
式中:A和B为任意常数;
在当弹性区域的应力状态必须满足r=r2时的边界条件,式(2)在r=rρ1处的连续条件为:
式中:上标(a)为塑性区内的应力分量,上标(b)为弹性区内的应力分量;
式(6)给出了弹性区域的应力解,当r=r2的时候得到:
式中,q2为套管塑性半径上的应力;σ0为原始水平地应力;rρ1为套管产生塑性变形的半径;rρ2为水泥环产生塑性变形的半径;q1为套管外壁与水泥环内径的层间压力。
作为本发明的进一步改进,已知弹性区域的应力情况,求得套管弹性区外径的径向位移:
作为本发明的进一步改进,计算方法中的计算过程中进行了以下假设:
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。平衡方程和应变位移关系用下式表示为:
这些控制方程要利用下面的边界条件求解:
由于考虑的是套管内压力不断增高的情况,所以此时σr≤0且σθ>0,可以观察得到σθ>σz=0≥σr。因此,在这种情况下的Tresca准则由如下形式表达:
σθ-σr=σ0 (3)
套管内交变压力f(pi)大于套管产生塑性变形的最大应力时,则会在套管内壁处开始产生塑性变形,随着套管内压力的加大,弹塑性边界也会慢慢的向套管外壁移动。用rρ1来表示套管弹塑性边界。
(1)套管塑性区域受力分析
当r1≤r≤rρ1时,套管处于塑性状态。平衡方程变为:
利用式(2)所示边界条件来解式(4)的微分方程得到:
式(5)即为套管塑性区内的应力分量。
(2)套管弹性区域受力分析
当r≥rρ1时,此时套管处在弹性区域中,应用弹性应力与应变的关系可得:
式中:X和Y为任意常数。
在当弹性区域的应力状态必须满足r=r2时的边界条件,式(5)在r=rρ1处的连续条件为:
式中:上标(a)为塑性区内的应力分量,上标(b)为弹性区内的应力分量。
式(9)给出了弹性区域的应力解,当r=r2的时候得到:
式中,q2为套管塑性半径上的应力,MPa;σ0为原始水平地应力,MPa;rρ1为套管产生塑性变形的半径,m;rρ2为水泥环产生塑性变形的半径,m;q1为套管外壁与水泥环内径的层间压力,MPa。
由于知道了弹性区域的应力情况,可以求得套管弹性区外径的径向位移:
以上的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,包括以下步骤:
将受到交变压力的套管分为塑性区域和弹性区域,套管内交变压力f(pi)大于套管产生塑性变形的最大应力时,则会在套管内壁处开始产生塑性变形,随着套管内压力的加大,弹塑性边界向套管外壁移动,用rρ1来表示套管弹塑性边界:当r1≤r≤rρ1时,分析得出套管塑性区内的应力分量;当r≥rρ1时,给出弹性区域的应力情况和套管弹性区外径的径向位移的计算方法。
2.根据权利要求1所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,塑性区域和弹性区域的应力解由平衡方程与应力方程为依据推出,套管弹性区外径的径向位移需满足边界位置Tresca准则。
4.根据权利要求2所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,
套管弹性区域的受力分析如下:
当r≥rρ1时,此时套管处在弹性区域中,应用弹性应力与应变的关系得:
式中:A和B为任意常数;
在当弹性区域的应力状态必须满足r=r2时的边界条件,式(2)在r=rρ1处的连续条件为:
式中:上标(a)为塑性区内的应力分量,上标(b)为弹性区内的应力分量;
式(6)给出了弹性区域的应力解,当r=r2的时候得到:
式中,q2为套管塑性半径上的应力;σ0为原始水平地应力;rρ1为套管产生塑性变形的半径;rρ2为水泥环产生塑性变形的半径;q1为套管外壁与水泥环内径的层间压力。
6.根据权利要求2所述的一种水平井套管交变压力作用下的套管力学分析计算方法,其特征在于,计算方法中的计算过程中进行了以下假设:
1)套管为弹塑性各向同性材料;
2)原场地应力为水平均布应力;
3)套管、水泥环以及地层三者是紧密连接的,无滑动,满足应力与位移的连续条件;
4)井眼垂直且套管居中。
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)
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)
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 | 中国地质大学(北京) | 考虑围岩与支护结构相互作用的圆形隧道力学计算方法 |
-
2020
- 2020-09-21 CN CN202010996834.8A patent/CN112257230A/zh active Pending
Patent Citations (3)
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)
Title |
---|
吴健: "固井水泥环密封完整性研究", 中国优秀硕士学位论文全文数据库工程科技Ⅰ辑, no. 02, pages 019 - 60 * |
Cited By (4)
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 |