CN109657334A - 一种超静定结构的受力特性分析方法 - Google Patents

一种超静定结构的受力特性分析方法 Download PDF

Info

Publication number
CN109657334A
CN109657334A CN201811536615.0A CN201811536615A CN109657334A CN 109657334 A CN109657334 A CN 109657334A CN 201811536615 A CN201811536615 A CN 201811536615A CN 109657334 A CN109657334 A CN 109657334A
Authority
CN
China
Prior art keywords
tubing string
packer
force
statically
statically indeterminate
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
Application number
CN201811536615.0A
Other languages
English (en)
Other versions
CN109657334B (zh
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.)
Xian Shiyou University
Original Assignee
Xian Shiyou University
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 Xian Shiyou University filed Critical Xian Shiyou University
Priority to CN201811536615.0A priority Critical patent/CN109657334B/zh
Publication of CN109657334A publication Critical patent/CN109657334A/zh
Application granted granted Critical
Publication of CN109657334B publication Critical patent/CN109657334B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

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)
  • Earth Drilling (AREA)

Abstract

一种超静定结构的受力特性分析方法,将封隔器坐封后的四封隔器管柱简化,管柱被四个固定端固定,每个固定端能提供三个反力,但是通常井中的管柱只承受轴向力,因此每个固定端只提供水平反力,从而此结构简化为3次超静定结构,如果去掉超静定结构中的三个多余约束,成为力法的基本结构;若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3,可以得到含有多余未知力的静定结构,求解后各封隔器管柱段轴向力和轴向变形的目的,本发明克服了传统分析方法中分段逐段求解各封隔器之间约束的局限性,可同时计算封隔器之间的每段管柱轴向力,所建立模型更符合多封隔器管柱坐封后的实际受力情况,计算结果也更精确。

Description

一种超静定结构的受力特性分析方法
技术领域
本发明涉及油气井开采技术领域,特别涉及一种超静定结构的受力特性分析方法。
背景技术
油气井开采过程中,为提高油气采收率,增大储层接触面积,常采用直井分层和水平井分段技术。为实现这两项技术,需要采用多封隔器(即多于三个封隔器)。封隔器的作用为定位和密封,从力学角度来说,多封隔器的存在引入了多余约束的问题,即原来的油井管结构变成了超静定结构。
从上世纪60年代开始至今,有很多国内外专家学者对油气井管柱力学进行了研究。1980年,Hammerlindl以双封隔器单一管柱为例,分析了管柱在四种基本效应下的轴向力及变形;1993年,冯建华等人建立了双封隔器复合管柱的轴向力计算模型,并介绍了封隔器受力及其承受压差的计算方法;2014年,吕占国分析了压裂过程中水平井多封隔器管柱的残余轴向力及强度安全性;2015年,练章华等人建立了多封隔器坐封等三种工况下的水平井压裂管柱力学计算模型,并对多级压裂管柱进行了安全性评价;2016年,张智等人建立了四种基本效应下的多封隔器复合管柱受力及变形计算模型;2017年,胡志强等人建立了双封隔器密闭环空压力预测模型,对不同工况下的封隔器管柱进行了安全校核。
现有的技术缺陷在于双封隔器结构中没有考虑超静定的约束条件,仅建立静力平衡方程进行求解,而多封隔器结构从变形协调条件建立超静定结构的变形方程,分段逐段计算管柱的受力和变形。
发明内容
为了克服上述现有技术的不足,本发明的目的在于提供一种超静定结构的受力特性分析方法,具有精确计算各封隔器管柱段轴向力和轴向变形的目的,特点在于充分考虑了前一个封隔器的坐封状态对后一个封隔器的影响,可同时计算封隔器间每段管柱的轴向力,而不是像传统的计算公式分段逐段进行计算管柱的受力和变形。
为了实现上述目的,本发明采用的技术方案是:
一种超静定结构的受力特性分析方法,包括以下步骤;
步骤一:模型简化;
封隔器坐封后的四封隔器管柱简化,管柱被四个固定端固定,其中从左往右第一个固定端为封隔器1对管柱的固定约束,第二个固定端为封隔器2对管柱的固定约束,第三个固定端为封隔器3对管柱的固定约束,第四个固定端为封隔器4对管柱的固定约束;
步骤二:超静定次数的判定;
每个固定端能提供水平X方向、竖直Y方向和弯矩M三个反力,由以上简化模型可知此结构为9次超静定结构,但是通常井中的管柱只承受轴向力,因此每个固定端只提供X方向的水平反力,从而此结构简化为3次超静定结构,其中三个多余约束分别为封隔器2对管柱的轴向约束、封隔器3对管柱的轴向约束以及封隔器4对管柱的轴向约束;
如果去掉超静定结构中的三个多余约束,则原结构变为静定结构,成为力法的基本结构;若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3,可以得到含有多余未知力的静定结构,此静定结构称为力法的基本体系;
步骤三:基本方程;
由于原超静定结构在任何情况下的变形量均为零,因此基本体系沿多余未知力X1、X2和X3方向的位移Δ1、Δ2和Δ3应与原结构相同,即:
由叠加原理可得此三次超静定结构的力法基本方程为:
式中:系数δi1在数值上等于基本结构在单位力X1=1单独作用下沿Xi方向产生的位移;系数δi2在数值上等于基本结构在单位力X2=1单独作用下沿Xi方向产生的位移;系数δi3在数值上等于基本结构在单位力X3=1单独作用下沿Xi方向产生的位移;自由项ΔiP为基本结构在外界因素(温度、压力、外载荷等)作用下沿Xi方向产生的位移;i=1,2,3;
方程(2)中各个系数及自由项的计算公式如下:
i,j=1,2,3
式中:是基本结构由于Xi=1作用而产生的内力,N;L为管柱的长度,m;E为管柱的弹性模量,Pa;A为管柱的横截面积,m2;εi为外界因素(温度、压力、外载荷等)使管柱产生的线应变;
步骤四:基本方程的求解
在基本体系中,当X1=1单独作用时,由平衡条件可得杆1的内力为1,杆2的内力为0;当X2=1单独作用时,由平衡条件可得杆1和杆2的内力均为1;当X3=1单独作用时,由平衡条件可得杆1、杆2和杆3的内力均为1,则可得方程中各个系数及自由项的大小;
将各个系数及自由项的计算所得值代入基本方程(2),则有:
超静定结构中各杆的内力可以根据叠加原理用下式计算:
解方程组(6),可得三段管柱的轴向力分别为:
由于水平段管柱温度改变量和压力改变量处处相同,即:
ΔT1=ΔT2=ΔT3=ΔT4=ΔT5=ΔT
Δpi1=Δpi2=Δpi3=Δpi4=Δpi5=Δpi
Δpo1=Δpo2=Δpo3=Δpo4=Δpo5=Δpo
因此,由温度效应和鼓胀效应引起的水平段管柱线应变也处处相等,即:
ε温1=ε温2=ε温3=ε温4=ε温5=βΔT
从而可知,由温度效应和鼓胀效应引起的水平段管柱轴向力处处相等,即:
F温1=F温2=F温3=-EAβΔT
本发明的有益效果:
本发明克服了传统分析方法中分段逐段求解各封隔器之间约束的局限性,可同时计算封隔器之间的每段管柱轴向力,所建立模型更符合多封隔器管柱坐封后的实际受力情况,计算结果也更精确。
附图说明
图1为水平井四封隔器管柱简图。
图2为四封隔器管柱简化模型。
图3为基本力系示意图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
以水平井四封隔器单一尺寸管柱为研究对象(如图1所示),从左往右分别为封隔器1、封隔器2、封隔器3、封隔器4;井口到水平段管柱的垂深为h;封隔器1、2之间的管柱长度为l1,为第一段管柱;封隔器2、3之间的管柱长度为l2,为第二段管柱;封隔器3、4之间的管柱长度为l3,为第三段管柱。
(1)模型简化
封隔器坐封后的四封隔器管柱简化模型如图2所示,由于本专利主要研究四个封隔器之间的三个管柱段所受轴向力的大小,所以封隔器1以前的管柱段以及封隔器4以后的管柱段不再考虑,略去。从图中可以看出,管柱被四个固定端固定,其中从左往右第一个固定端为封隔器1对管柱的固定约束,第二个固定端为封隔器2对管柱的固定约束,第三个固定端为封隔器3对管柱的固定约束,第四个固定端为封隔器4对管柱的固定约束。
(2)超静定次数的判定
因为每个固定端能提供水平X方向、竖直Y方向和弯矩M三个反力,而对于图2所示结构只能列出三个平衡方程,所以由以上简化模型可知此结构为9次超静定结构。但是通常井中的管柱只承受轴向力,因此每个固定端只提供X方向的水平反力,从而此结构简化为3次超静定结构,其中三个多余约束分别为封隔器2对管柱的轴向约束、封隔器3对管柱的轴向约束以及封隔器4对管柱的轴向约束。
如果去掉超静定结构中的三个多余约束,则原结构变为静定结构,成为力法的基本结构;若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3,可以得到含有多余未知力的静定结构,此静定结构称为力法的基本体系,如图3所示。
(3)基本方程
由于原超静定结构在任何情况下的变形量均为零,因此基本体系沿多余未知力X1、X2和X3方向的位移Δ1、Δ2和Δ3应与原结构相同,即:
由叠加原理可得此三次超静定结构的力法基本方程为:
式中:系数δi1在数值上等于基本结构在单位力X1=1单独作用下沿Xi方向产生的位移;系数δi2在数值上等于基本结构在单位力X2=1单独作用下沿Xi方向产生的位移;系数δi3在数值上等于基本结构在单位力X3=1单独作用下沿Xi方向产生的位移;自由项ΔiP为基本结构在外界因素(温度、压力、外载荷等)作用下沿Xi方向产生的位移;i=1,2,3。
方程(2)中各个系数及自由项的计算公式如下:
i,j=1,2,3
式中:是基本结构由于Xi=1作用而产生的内力,N;L为管柱的长度,m;E为管柱的弹性模量,Pa;A为管柱的横截面积,m2;εi为外界因素(温度、压力、外载荷等)使管柱产生的线应变。
(4)基本方程的求解
由图3可知,在基本体系中,当X1=1单独作用时,由平衡条件可得杆1的内力为1,杆2的内力为0;当X2=1单独作用时,由平衡条件可得杆1和杆2的内力均为1;当X3=1单独作用时,由平衡条件可得杆1、杆2和杆3的内力均为1。则可得方程中各个系数及自由项的大小如表1和表2所示。
表1系数的计算
表2自由项的计算
将表中各个系数及自由项的计算所得值代入基本方程(2),则有:
超静定结构中各杆的内力可以根据叠加原理用下式计算:
解方程组(6),可得三段管柱的轴向力分别为:
由于水平段管柱温度改变量和压力改变量处处相同,即:
ΔT1=ΔT2=ΔT3=ΔT4=ΔT5=ΔT
Δpi1=Δpi2=Δpi3=Δpi4=Δpi5=Δpi
Δpo1=Δpo2=Δpo3=Δpo4=Δpo5=Δpo
因此,由温度效应和鼓胀效应引起的水平段管柱线应变也处处相等,即:
ε温1=ε温2=ε温3=ε温4=ε温5=βΔT
从而可知,由温度效应和鼓胀效应引起的水平段管柱轴向力处处相等,即:
F温1=F温2=F温3=-EAβΔT
由所推导的公式,可以定量求解多封隔器管柱这种超静定结构在不同工况下的温度效应和鼓胀效应,及在这两种效应下各个管柱段间的力学特性,进而分析管柱的强度安全性,为安全施工提供参数。
本专利运用结构力学中的力法对以超静定结构为特点的多封隔器管柱轴向力进行分析,可同时计算封隔器之间的每段管柱轴向力,更符合实际工况下多封隔器管柱的受力情况,计算结果也更精确。
本专利以直井多封隔器单一尺寸管柱为研究对象,考虑温度效应和鼓胀效应对管柱的影响,运用结构力学中的力法对以超静定结构为特点的多封隔器管柱轴向力进行分析,建立用于同时计算封隔器之间的每段管柱轴向力的计算模型,最后以四封隔器为例推导两种效应下的管柱轴向力计算公式,并和传统的计算方法进行比较,计算两者之间的差值。

Claims (1)

1.一种超静定结构的受力特性分析方法,其特征在于,包括以下步骤;
步骤一:模型简化;
封隔器坐封后的四封隔器管柱简化,管柱被四个固定端固定,其中从左往右第一个固定端为封隔器1对管柱的固定约束,第二个固定端为封隔器2对管柱的固定约束,第三个固定端为封隔器3对管柱的固定约束,第四个固定端为封隔器4对管柱的固定约束;
步骤二:超静定次数的判定;
每个固定端能提供水平X方向、竖直Y方向和弯矩M三个反力,由以上简化模型可知此结构为9次超静定结构,但是通常井中的管柱只承受轴向力,因此每个固定端只提供X方向的水平反力,从而此结构简化为3次超静定结构,其中三个多余约束分别为封隔器2对管柱的轴向约束、封隔器3对管柱的轴向约束以及封隔器4对管柱的轴向约束;
如果去掉超静定结构中的三个多余约束,则原结构变为静定结构,成为力法的基本结构;若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3,可以得到含有多余未知力的静定结构,此静定结构称为力法的基本体系;
步骤三:基本方程;
由于原超静定结构在任何情况下的变形量均为零,因此基本体系沿多余未知力X1、X2和X3方向的位移Δ1、Δ2和Δ3应与原结构相同,即:
由叠加原理可得此三次超静定结构的力法基本方程为:
式中:系数δi1在数值上等于基本结构在单位力X1=1单独作用下沿Xi方向产生的位移;系数δi2在数值上等于基本结构在单位力X2=1单独作用下沿Xi方向产生的位移;系数δi3在数值上等于基本结构在单位力X3=1单独作用下沿Xi方向产生的位移;自由项ΔiP为基本结构在外界因素作用下沿Xi方向产生的位移;i=1,2,3;
方程(2)中各个系数及自由项的计算公式如下:
i,j=1,2,3
式中:是基本结构由于Xi=1作用而产生的内力,N;L为管柱的长度,m;E为管柱的弹性模量,Pa;A为管柱的横截面积,m2;εi为外界因素使管柱产生的线应变;
步骤四:基本方程的求解
在基本体系中,当X1=1单独作用时,由平衡条件可得杆1的内力为1,杆2的内力为0;当X2=1单独作用时,由平衡条件可得杆1和杆2的内力均为1;当X3=1单独作用时,由平衡条件可得杆1、杆2和杆3的内力均为1,则可得方程中各个系数及自由项的大小。
CN201811536615.0A 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法 Active CN109657334B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811536615.0A CN109657334B (zh) 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811536615.0A CN109657334B (zh) 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法

Publications (2)

Publication Number Publication Date
CN109657334A true CN109657334A (zh) 2019-04-19
CN109657334B CN109657334B (zh) 2023-03-31

Family

ID=66114274

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811536615.0A Active CN109657334B (zh) 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法

Country Status (1)

Country Link
CN (1) CN109657334B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111651922A (zh) * 2020-06-03 2020-09-11 重庆大学 力法计算的模拟手算分析方法、系统、装置及存储介质
CN112507503A (zh) * 2019-09-16 2021-03-16 深圳市建筑设计研究总院有限公司 一种次生自反力结构的加载及预反力计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020166374A1 (en) * 2001-02-08 2002-11-14 Demarchos Andronikos S. Method for analysing a completion system
CN102720185A (zh) * 2012-02-03 2012-10-10 招商局重庆交通科研设计院有限公司 h型抗滑桩桩身结构段分式计算方法
CN104166792A (zh) * 2014-08-06 2014-11-26 中国科学院工程热物理研究所 一种预应力混凝土连续刚构桥温度作用有限元分析方法
CN106548030A (zh) * 2016-11-03 2017-03-29 中国石油天然气股份有限公司 一种多封隔器复合管柱受力分析计算方法
WO2017121112A1 (zh) * 2016-01-14 2017-07-20 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 一种多级水力喷射压裂管柱性能评价方法
CN107563000A (zh) * 2017-08-01 2018-01-09 哈尔滨工程大学 一种内燃机曲轴系统纵扭耦合刚度计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020166374A1 (en) * 2001-02-08 2002-11-14 Demarchos Andronikos S. Method for analysing a completion system
CN102720185A (zh) * 2012-02-03 2012-10-10 招商局重庆交通科研设计院有限公司 h型抗滑桩桩身结构段分式计算方法
CN104166792A (zh) * 2014-08-06 2014-11-26 中国科学院工程热物理研究所 一种预应力混凝土连续刚构桥温度作用有限元分析方法
WO2017121112A1 (zh) * 2016-01-14 2017-07-20 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 一种多级水力喷射压裂管柱性能评价方法
CN106548030A (zh) * 2016-11-03 2017-03-29 中国石油天然气股份有限公司 一种多封隔器复合管柱受力分析计算方法
CN107563000A (zh) * 2017-08-01 2018-01-09 哈尔滨工程大学 一种内燃机曲轴系统纵扭耦合刚度计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
丁瑞娜等: "超静定结构的分析与应用", 《城市地理》 *
沈蒲生等: "超静定结构徐变效应的力法分析方法", 《铁道科学与工程学报》 *
王仁华: ""结构力学"课程中力法与位移法的对偶关系剖析", 《中国建设教育》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112507503A (zh) * 2019-09-16 2021-03-16 深圳市建筑设计研究总院有限公司 一种次生自反力结构的加载及预反力计算方法
CN111651922A (zh) * 2020-06-03 2020-09-11 重庆大学 力法计算的模拟手算分析方法、系统、装置及存储介质

Also Published As

Publication number Publication date
CN109657334B (zh) 2023-03-31

Similar Documents

Publication Publication Date Title
Li et al. The self-propelled force model of a multi-orifice nozzle for radial jet drilling
CN102680176B (zh) 一种管材实物应力腐蚀试验机
Fan et al. Experimental investigation on fire resistance of stainless steel columns with square hollow section
CN109657334A (zh) 一种超静定结构的受力特性分析方法
CN103410516B (zh) 小直径隧洞全断面开挖围岩变形预警工法
CN104373053B (zh) 一种地下储气库注采管柱设计方法
CN107084821A (zh) 一种高温高压水平井封隔器模拟试验装置
CN102322989B (zh) 一种井下抽油杆径向压力测量装置
CN103411771A (zh) 管线球阀试验用的大型环境模拟试验装置
CN105510148A (zh) 一种高温下封隔器胶筒接触应力测试装置及方法
Yu et al. Interference effects between two high-rise buildings on wind-induced torsion
CN111795788A (zh) 一种机械加载式封隔器胶筒模拟加载密封实验装置及方法
CN203432794U (zh) 管线球阀试验用的大型环境模拟试验装置
Qian-bei et al. The research of post-buckling about slender rod string in wellbore based on energy method and experiment
Guo et al. Nonlinear flow-induced vibration response characteristics of a tubing string in HPHT oil&gas well
CN105651606A (zh) 一种高温高压焊接构件结构完整性的评定方法
CN105259041A (zh) 一种模拟注蒸汽井套管热强度的测试方法及装置
CN104318054A (zh) 延迟焦化焦炭塔在线监测及安全评估系统
CN104453850B (zh) 多级油管柱参数预测方法及装置
CN201060097Y (zh) 一种连续管多轴复合载荷低周疲劳试验装置
CN204203007U (zh) 一种高温下封隔器胶筒接触应力测试装置
CN204920904U (zh) 剪销式旁通阀
CN203519239U (zh) 一种活塞式菱形结构光纤Bragg光栅渗压传感器
CN103670372B (zh) 一种热采井套管柱应变判断系统及其方法
CN106767356B (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