CN106529092A - 钻井套管的剪切变形力有限元计算方法 - Google Patents
钻井套管的剪切变形力有限元计算方法 Download PDFInfo
- Publication number
- CN106529092A CN106529092A CN201611154853.6A CN201611154853A CN106529092A CN 106529092 A CN106529092 A CN 106529092A CN 201611154853 A CN201611154853 A CN 201611154853A CN 106529092 A CN106529092 A CN 106529092A
- Authority
- CN
- China
- Prior art keywords
- finite element
- shearing
- well casing
- sleeve pipe
- casing
- 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
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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种钻井套管的剪切变形力有限元计算方法，该方法是基于有限元软件模拟由于地层的错动导致钻井套管发生剪切变形的计算，该方法包括有11个步骤。本发明的效果是通过钻井套管的剪切变形力有限元计算方法，能够精确的模拟出套管在井下受到地层错动力的作用时管体的形貌，以及管体所能承受的最大剪切错动力的大小。通过有限元计算可大幅度降低进行钻井固井套管剪切试验的成本。本发明计算的剪切力结果与实物试验进行了对比，误差在5％以内，本发明的计算结果对套管管串设计有重要的指导意义。
Description
技术领域
本发明涉及钻井和采油技术领域，涉及一种钻井套管的剪切变形力有限元计算方法。
背景技术
在钻井抽采过程中，油层出砂后会造成地层疏松，形成空洞。导致地层支撑的损失和滑移，将对钻井固井套管产生横向的剪切作用，很容易导致套管的剪切破坏甚至完全错断，会使钻井被迫停止抽采，造成巨大的经济损失。
中国专利文献号为ZL201610333516.7公开了一种钻井固井套管剪切试验系统，可以对不同管径的固井套管施加剪切力，用来模拟岩层的相对位移对固井套管的剪切作用。但是该文献只能通过试验的方法得到套管的剪切变形力，试验耗费大量的物力财力，并且试验设备的最大剪切载荷有一定的限制，对所进行的试验规格也有一定的局限性。
发明内容
本发明的目的是提供一种钻井套管的剪切变形力有限元计算方法，可以定量分析套管的侧向位移变化和套管能承受的最大剪切力，能为预防套管损坏和屈曲变形提供了理论依据，确定套管在井下的剪切破坏力的大小，从而更好的为油田选择合适的套管。为钻井工程提供技术支持。
为实现上述目的，本发明提供的钻井套管的剪切变形力有限元计算方法包括以下步骤：
1)通过有限元软件作图建立钻井套管的三维薄壁圆筒几何模型，该薄壁圆筒的外径是钻井套管的外径D，壁厚为套管的壁厚t，该薄壁圆筒长度L为1000mm；
2)通过有限元软件在上述薄壁圆筒几何模型的管体总长对称位置建立两个剪切环的三维薄壁圆筒几何模型，剪切环的内径dj与管体外径D相同，剪切环的外径Dj大于管体外径D，剪切环的长度s为20～60mm，两个剪切环的距离a为127mm；
3)所述套管为碳钢或低合金钢材质，在有限元软件中根据材料的特性设置所述套管的弹性模量、应力应变曲线、泊松比、材料密度的性能；
4)设置两个剪切环为刚体，即在任何情况下都不发生任何变形；
5)在有限元软件中，在计算前需要设置的加载过程为单个或多个步骤即分析步，设定钻井套管的剪切变形为准静态，设置分析步的形式为动态隐式分析，计算时间为1秒，迭代时长0.01秒；
6)对所建立的套管三维薄壁圆筒的几何模型进行网格划分，网格划分的方式是先设置种子点，后进行有限元软件的自动划分网格，在壁厚方向上至少划分3层网格，网格的形式为六面体单元；
7)约束该套管两端面的径向位移为0；
8)设置套管的外表面与两个剪切环的内表面之间的接触形式为面-面接触；
9)对两个剪切环中的一个施加径向位移，另一个质心施加大小相等方向相反的位移，位移值至少为4mm；
10)通过有限元软件的求解器实现步骤1～9所完成的前处理步骤的运算求解，计算的结果将保存在后处理文件中；
11)从步骤10得到的计算结果中提取出管体所受的剪切力的值，并作出随着剪切环的径向位移增加的管体的剪切力值曲线，对所述曲线进行分析，剪切力的值随着剪切环的位移增加而增速减缓，之后开始有突增和突降的突变，在突变之前的最大值为所述套管能承受的剪切力最大值。
本发明的效果是，通过钻井套管的剪切变形力有限元计算方法，能够精确的模拟出套管在井下受到地层错动力的作用时管体的形貌，以及管体所能承受的最大剪切错动力的大小。通过有限元计算可大幅度降低进行钻井固井套管剪切试验的成本。本发明计算的剪切力结果与实物试验进行了对比，误差在5％以内，本发明的计算结果对套管管串设计有重要的指导意义。
附图说明
图1是本发明的套管的剪切变形力有限元计算方法的三维有限元模型示意图；
图2是本发明的剪切环示意图
图3是本发明的套管网格划分示意图；
图4是本发明的套管的剪切变形力有限元计算方法提取出的套管所承受的剪切力与剪切环错动位移的曲线图。
图5-1、5-2是本发明的套管的剪切变形力有限元计算方法的有限元计算结果示意图。
图中：
A：管体；B：剪切环；C：受到剪切力后变形的管体；D：管体外径；t：管体壁厚；L：管体长度；s：剪切环的长度度；a：两个剪切环之间的距离；Dj：剪切环外径；dj：剪切环内径
具体实施方式
结合附图对本发明的钻井套管的剪切变形力有限元计算方法加以说明。
本发明提供的钻井套管的剪切变形力有限元计算方法包括以下步骤：
1)通过有限元软件作图建立钻井套管的三维薄壁圆筒几何模型，该薄壁圆筒的外径是钻井套管的外径D，壁厚为套管的壁厚t，该薄壁圆筒长度L为1000mm，如图1的A所示。
2)通过有限元软件在上述薄壁圆筒几何模型的管体总长对称位置建立两个剪切环的三维薄壁圆筒几何模型，剪切环的内径dj与管体外径D相同，剪切环的外径Dj大于管体外径D，剪切环的长度s为20～60mm，两个剪切环的距离a为127mm，如图1的B和图2所示。
3)所述套管为碳钢或低合金钢材质，在有限元软件中根据材料的特性设置所述套管的弹性模量、应力应变曲线、泊松比、材料密度的性能。
4)设置两个剪切环为刚体，即在任何情况下都不发生任何变形。
5)在有限元软件中，在计算前需要设置的加载过程为单个或多个步骤即分析步，设定钻井套管的剪切变形为准静态，设置分析步的形式为动态隐式分析，计算时间为1秒，迭代时长0.01秒。
6)对所建立的套管三维薄壁圆筒的几何模型进行网格划分，网格划分的方式是先设置种子点，后进行有限元软件的自动划分网格，在壁厚方向上至少划分3层网格，网格的形式为六面体单元，如图3所示。
7)约束该套管两端面的径向位移为0。
8)设置套管的外表面与两个剪切环的内表面之间的接触形式为面-面接触；
9)对两个剪切环中的一个施加径向位移，另一个质心施加大小相等方向相反的位移，位移值至少为4mm。
10)本步骤对1～9步建立的有限元模型进行求解。在求解时，将方程在空间上采用有限元方法离散后，变成常微分方程：
F＝M(u)+C(u)+K(u)
通过动态隐式分析所用的NewMark法求解该方程，任一时刻的位移、速度、加速度都相互关联，采用迭代和求解联立方程实现求解。计算的结果将保存在后处理文件中。
11)从步骤10)得到的计算结果中提取出管体所受的剪切力的值，并作出随着剪切环的径向位移增加的管体的剪切力值曲线，对所述曲线进行分析，剪切力的值随着剪切环的位移增加而增速减缓，之后开始有突增和突降的突变，在突变之前的最大值为所述套管能承受的剪切力最大值。如图4所示，在突变之前的最大值为所述套管能承受的剪切力最大值。图5-1、5-2为剪切力导致的管体发生变形示意图。
下面给出两个具体计算实例。
实施例1
套管规格外径D＝139.7mm，壁厚t＝12.7mm，钢级Q125，其弹性模量E＝206GPa、泊松比μ＝0.3、屈服强度σS＝862～1034MPa、抗拉强度σb>931MPa。根据上述有限元计算法进行建模，管体长度1000mm，两个剪切环宽度s＝40mm，相距a＝127mm。图4为通过有限元软件提取的剪切力与剪切环的移动的距离的关系图，可以得到剪切力在突变前的最大值1625KN，所以该值为该套管能承受的最大剪切力。
实施例2
套管规格外径D＝146.7mm，壁厚t＝15.9mm，钢级TP140V，其弹性模量E＝206GPa、泊松比μ＝0.3、屈服强度σS＝985～1080MPa、抗拉强度σb>1034MPa。根据上述有限元计算法进行建模，管体长度1000mm，两个剪切环宽度s＝40mm，相距a＝127mm。通过有限元软件可以得到剪切力在突变前的最大值2405KN，所以该值为该套管能承受的最大剪切力。
实施例3
套管规格外径D＝177.8mm，壁厚t＝10.36mm，钢级P110，其弹性模量E＝206GPa、泊松比μ＝0.3、屈服强度σS＝758～965MPa、抗拉强度σb>862MPa。根据上述有限元计算法进行建模，管体长度1000mm，两个剪切环宽度s＝50mm，相距a＝127mm。通过有限元软件可以得到剪切力在突变前的最大值1454KN，所以该值为该套管能承受的最大剪切力。
Claims (1)
1.一种钻井套管的剪切变形力有限元计算方法，该方法是基于有限元软件模拟由于地层的错动导致钻井套管发生剪切变形的计算，其特征是：该方法包括以下步骤：
1)通过有限元软件作图建立钻井套管的三维薄壁圆筒几何模型，该薄壁圆筒的外径是钻井套管的外径D，壁厚为套管的壁厚t，该薄壁圆筒长度L为1000mm；
2)通过有限元软件在上述薄壁圆筒几何模型的管体总长对称位置建立两个剪切环的三维薄壁圆筒几何模型，剪切环的内径dj与管体外径D相同，剪切环的外径Dj大于管体外径D，剪切环的长度s为20～60mm，两个剪切环的距离a为127mm；
3)所述套管为碳钢或低合金钢材质，在有限元软件中根据材料的特性设置所述套管的弹性模量、应力应变曲线、泊松比、材料密度的性能；
4)设置两个剪切环为刚体，即在任何情况下都不发生任何变形；
5)在有限元软件中，在计算前需要设置的加载过程为单个或多个步骤即分析步，设定钻井套管的剪切变形为准静态，设置分析步的形式为动态隐式分析，计算时间为1秒，迭代时长0.01秒；
6)对所建立的套管三维薄壁圆筒的几何模型进行网格划分，网格划分的方式是先设置种子点，后进行有限元软件的自动划分网格，在壁厚方向上至少划分3层网格，网格的形式为六面体单元；
7)约束所述钻井套管两端面的径向位移为0；
8)设置所述钻井套管的外表面与两个剪切环的内表面之间的接触形式为面-面接触；
9)对所述两个剪切环中的一个施加径向位移，另一个质心施加大小相等方向相反的位移，位移量至少为4mm；
10)通过有限元软件的求解器实现步骤1～9所完成的前处理步骤的运算求解，计算的结果将保存在后处理文件中；
11)从步骤10)得到的计算结果中提取出管体所受的剪切力的值，并作出随着剪切环的径向位移增加的管体的剪切力值曲线，对所述曲线进行分析，剪切力的值随着剪切环的位移增加而增速减缓，之后开始有突增和突降的突变，在突变之前的最大值为所述套管能承受的剪切力最大值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611154853.6A CN106529092B (zh) | 2016-12-14 | 2016-12-14 | 钻井套管的剪切变形力有限元计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611154853.6A CN106529092B (zh) | 2016-12-14 | 2016-12-14 | 钻井套管的剪切变形力有限元计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106529092A true CN106529092A (zh) | 2017-03-22 |
| CN106529092B CN106529092B (zh) | 2020-01-21 |
Family
ID=58339550
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611154853.6A Active CN106529092B (zh) | 2016-12-14 | 2016-12-14 | 钻井套管的剪切变形力有限元计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106529092B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106991235A (zh) * | 2017-04-01 | 2017-07-28 | 中国石油天然气集团公司 | 固井水泥环完整性分析预测方法及装置 |
| CN107169240A (zh) * | 2017-06-22 | 2017-09-15 | 贵州财经大学 | 基于阶梯孔的回拖计算方法和装置 |
| CN108279173A (zh) * | 2018-01-02 | 2018-07-13 | 中国石油天然气集团公司 | 一种套管抗剪切性能评价方法 |
| CN110717284A (zh) * | 2019-08-16 | 2020-01-21 | 中国石油天然气集团有限公司 | 一种套管抗剪切承载能力分析及试验评价方法 |
| CN111625979A (zh) * | 2020-05-28 | 2020-09-04 | 中国船舶工业集团公司第七0八研究所 | 一种桩腿插销孔的强度校核加载方法 |
| CN113673126A (zh) * | 2021-07-28 | 2021-11-19 | 中国石油大学(北京) | 一种钻井用多层套管环空带压有限元计算方法和装置 |
| CN113824070A (zh) * | 2021-09-22 | 2021-12-21 | 深圳市骏鼎达新材料股份有限公司 | 一种具有自定位功能的保护套管 |
| CN113836656A (zh) * | 2021-09-14 | 2021-12-24 | 临海伟星新型建材有限公司 | 采用有限元算法的氟硅改性pert阻隔内衬管缩径量计算方法 |
| WO2022227488A1 (zh) * | 2021-04-28 | 2022-11-03 | 天津钢管制造有限公司 | 双层套管抗外压挤毁能力的获取方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105760564A (zh) * | 2014-12-19 | 2016-07-13 | 中国石油天然气股份有限公司 | 分析油层套管失效的方法和装置 |
| CN106018125A (zh) * | 2016-05-19 | 2016-10-12 | 中国矿业大学 | 一种钻井固井套管剪切实验系统 |
-
2016
- 2016-12-14 CN CN201611154853.6A patent/CN106529092B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105760564A (zh) * | 2014-12-19 | 2016-07-13 | 中国石油天然气股份有限公司 | 分析油层套管失效的方法和装置 |
| CN106018125A (zh) * | 2016-05-19 | 2016-10-12 | 中国矿业大学 | 一种钻井固井套管剪切实验系统 |
Non-Patent Citations (2)
| Title |
|---|
| 于浩 等: "页岩气体积压裂过程中套管失效机理研究", 《中国安全生产科学技术》 * |
| 蒋可 等: "页岩气水平井固井质量对套管损坏的影响", 《天然气工业》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106991235B (zh) * | 2017-04-01 | 2020-03-10 | 中国石油天然气集团公司 | 固井水泥环完整性分析预测方法及装置 |
| CN106991235A (zh) * | 2017-04-01 | 2017-07-28 | 中国石油天然气集团公司 | 固井水泥环完整性分析预测方法及装置 |
| CN107169240A (zh) * | 2017-06-22 | 2017-09-15 | 贵州财经大学 | 基于阶梯孔的回拖计算方法和装置 |
| CN107169240B (zh) * | 2017-06-22 | 2020-09-11 | 贵州财经大学 | 基于阶梯孔的回拖计算方法和装置 |
| CN108279173A (zh) * | 2018-01-02 | 2018-07-13 | 中国石油天然气集团公司 | 一种套管抗剪切性能评价方法 |
| CN108279173B (zh) * | 2018-01-02 | 2020-08-07 | 中国石油天然气集团公司 | 一种套管抗剪切性能评价方法 |
| CN110717284A (zh) * | 2019-08-16 | 2020-01-21 | 中国石油天然气集团有限公司 | 一种套管抗剪切承载能力分析及试验评价方法 |
| CN111625979A (zh) * | 2020-05-28 | 2020-09-04 | 中国船舶工业集团公司第七0八研究所 | 一种桩腿插销孔的强度校核加载方法 |
| CN111625979B (zh) * | 2020-05-28 | 2023-08-15 | 中国船舶工业集团公司第七0八研究所 | 一种桩腿插销孔的强度校核加载方法 |
| WO2022227488A1 (zh) * | 2021-04-28 | 2022-11-03 | 天津钢管制造有限公司 | 双层套管抗外压挤毁能力的获取方法 |
| CN113673126A (zh) * | 2021-07-28 | 2021-11-19 | 中国石油大学(北京) | 一种钻井用多层套管环空带压有限元计算方法和装置 |
| CN113673126B (zh) * | 2021-07-28 | 2024-02-13 | 中国石油大学(北京) | 一种钻井用多层套管环空带压有限元计算方法和装置 |
| CN113836656A (zh) * | 2021-09-14 | 2021-12-24 | 临海伟星新型建材有限公司 | 采用有限元算法的氟硅改性pert阻隔内衬管缩径量计算方法 |
| WO2023039932A1 (zh) * | 2021-09-14 | 2023-03-23 | 临海伟星新型建材有限公司 | 采用有限元算法的氟硅改性pert阻隔内衬管缩径量计算方法 |
| CN113824070B (zh) * | 2021-09-22 | 2022-11-18 | 深圳市骏鼎达新材料股份有限公司 | 一种具有自定位功能的保护套管 |
| CN113824070A (zh) * | 2021-09-22 | 2021-12-21 | 深圳市骏鼎达新材料股份有限公司 | 一种具有自定位功能的保护套管 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106529092B (zh) | 2020-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106529092A (zh) | 钻井套管的剪切变形力有限元计算方法 | |
| CN105319603B (zh) | 致密砂岩储层复杂网状裂缝的预测方法 | |
| CN108868748B (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
| CN109271738B (zh) | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 | |
| CN104361211B (zh) | 岩石统计损伤本构模型的构建和应用方法 | |
| CN112036098A (zh) | 一种深层油气藏水力裂缝扩展数值模拟的方法 | |
| CN110005407B (zh) | 盐穴储气库稳定性评价方法 | |
| US9836561B2 (en) | Optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics | |
| CN113820750B (zh) | 基于弹塑性力学定量预测泥岩构造裂缝的方法 | |
| CN104005750A (zh) | 基于钻柱-钻头-岩石-井筒系统全尺寸钻进模拟的钻井提速评价方法 | |
| CN104963657A (zh) | 基于高含水期监测数据约束的油藏数值模拟方法及装置 | |
| CN109033504B (zh) | 一种油水井套管损坏预测方法 | |
| CN105089595A (zh) | 水平压裂裂缝导流作用下的油藏数值模拟方法及装置 | |
| CN108952657B (zh) | 一种水平井平台压裂裂缝长度确定方法及装置 | |
| CN106372297A (zh) | 砂土岩溶地层中盾构与溶洞安全垂直距离的确定方法 | |
| CN106555575A (zh) | 一种深层页岩的地层起裂压力的计算方法 | |
| CN107991188B (zh) | 一种基于岩心残余应力水平预测水力裂缝复杂性的方法 | |
| Xi et al. | Research on the influence of strike-slip fault slippage on production casing and control methods and engineering application during multistage fracturing in deep shale gas wells | |
| Yu et al. | On how asymmetric stimulated rock volume in shales may impact casing integrity | |
| Zhang et al. | Strata movement and stress evolution when mining two overlapping panels affected by hard stratum | |
| Yan et al. | A multiscale analysis of adjacent fault dislocation mechanism induced by tunnel excavation based on continuous-discrete coupling method | |
| CN115324557A (zh) | 基于多因素分析预测压裂诱发套管变形风险程度的方法 | |
| CN115081294B (zh) | 高水头钢筋混凝土衬砌压力隧洞充水过程分析方法 | |
| CN114662316A (zh) | 一种计算空间三维天然裂缝壁面内摩擦系数的方法 | |
| CN106372295A (zh) | 砂土岩溶地层中盾构与溶洞安全水平距离的确定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190628 Address after: No. 396 Jintang Highway, Dongli District, Tianjin 300301 Applicant after: Tianjin Steel Tube Manufacturing Co., Ltd. Address before: No. 396 Jintang Highway, Dongli District, Tianjin 300301 Applicant before: Tianjin Steel Pipe Group Co., Ltd. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |