唐 洋,何 胤,李 旺,舒將軍,何玉發(fā),孫 鵬,姚佳鑫
(1.西南石油大學(xué) 機(jī)電工程學(xué)院,四川 成都 610500; 2.南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(湛江),廣東 湛江 524088;3.中海油研究總院有限責(zé)任公司,北京 100027)
在進(jìn)行深水油氣測試時,當(dāng)遇到臺風(fēng)、海嘯等緊急情況及其他安全隱患時,可能導(dǎo)致浮動式鉆井平臺及鉆井船被摧毀,下部管柱受到破壞,造成極大的人員傷亡和環(huán)境污染[1-2]。所以遇到特殊海況時必須采用一個關(guān)鍵性工具——水下測試樹(Subsea Test Tree,SSTT)來封堵管柱內(nèi)高壓油氣,切斷連續(xù)油管或鋼絲繩,并迅速斷開作業(yè)管柱,使得鉆井平臺可以撤離,實(shí)現(xiàn)對深水測試作業(yè)的安全保護(hù)[3-4]。
球閥剪切模塊是實(shí)現(xiàn)管柱斷開的保障,球閥剪斷連續(xù)油管及電纜等,使得絞車可以上提連續(xù)油管,測試管柱內(nèi)的防噴球閥等部件能夠關(guān)閉,防止被連續(xù)油管卡住。所以球閥能否成功剪切連續(xù)油管及剪切性能的好壞直接影響到深水測試過程作業(yè)安全。
目前國內(nèi)外對于井下管柱剪切的分析如下:John等[5]評估了API油管剪切要求,并利用有限元方法分析了閘板防噴器剪切油管中各個參數(shù)對剪切力影響;Amy 等[6]采用FEA有限元分析和CFD流體分析結(jié)合的方式分析了閘板防噴器剪切鉆桿過程及相關(guān)分析方法討論;劉冰等[7]分析了閘板防噴器剪切鉆桿過程中的力學(xué)性能;孟祥喻等[8]通過理論計(jì)算與工程應(yīng)用參數(shù)相結(jié)合提出了深水水下剪切閘板防噴器的剪切能力評估方法,并通過實(shí)例對該評估方法進(jìn)行了應(yīng)用,應(yīng)用結(jié)果顯示該方法所計(jì)算的數(shù)值與廠家提供的數(shù)據(jù)誤差在 3% 之內(nèi);劉冬等[9]采用有限元軟件模擬油管變形過程,表明油管受力過程中最易出現(xiàn)應(yīng)力集中和裂紋的部位是水平外徑和上下對徑4個塑形鉸處。
雖然目前國內(nèi)外專家對于油管的剪切分析很多,但是對于球閥剪切連續(xù)油管的分析在國內(nèi)外卻非常少,目前僅國外斯倫貝謝公司Seranton等[10]曾做過球閥對不同直徑連續(xù)油管的剪切實(shí)驗(yàn),但是由于實(shí)驗(yàn)條件限制,只得到不同管厚油管剪切斷口形狀。就國內(nèi)來說,水下測試樹的研究尚處于發(fā)展初期,對于球閥剪切能力分析處于空白,水下測試樹剪切功能雖然僅在危險情況下使用,但對其研究卻是不可或缺的。
本文基于有限元理論建立了球閥剪切連續(xù)油管的數(shù)值仿真模型,得到了球閥剪切油管剪切力矩變化規(guī)律,分析了球閥大小、球閥剪切速度、連續(xù)油管承載的拉力等工況對球閥剪切力的影響,并與實(shí)驗(yàn)結(jié)果進(jìn)行比較,研究工作為水下測試樹球閥結(jié)構(gòu)設(shè)計(jì)和選用提供理論依據(jù)。
水下測試樹剪切球閥部分結(jié)構(gòu)組成如圖1所示。其剪切原理為:由多功能控制模塊發(fā)出控制信號,并輸入液壓油;當(dāng)液壓油從液壓流道1流入,液動活塞會被向上推動,液動活塞上移的過程中會帶動驅(qū)動滑塊上移,驅(qū)動滑塊推動球閥旋轉(zhuǎn)使得球閥剪切連續(xù)油管至最終關(guān)閉。
圖1 水下測試樹原理Fig.1 Principle of subsea test tree
金屬剪切過程中的真實(shí)應(yīng)力和應(yīng)變特性十分復(fù)雜,涉及彈塑性力學(xué)、金屬物理、損傷與斷裂、計(jì)算方法等多學(xué)科的內(nèi)容。所以,選擇用于表征連續(xù)油管的材料強(qiáng)度和破壞模型的選擇是數(shù)值模擬的一個重要方面,在金屬剪切有限元分析中常用到Johnson-Cook塑性本構(gòu)模型[11-12],該方法適用于大變形分析,其表達(dá)式見式(1),而本文分析忽略溫度影響,理論分析式如式(2)所示:
(1)
(2)
經(jīng)典塑形理論描述的是一種均勻的、連續(xù)的介質(zhì)在外力作用下產(chǎn)生的不可恢復(fù)的位移或者滑移現(xiàn)象[13-14]。金屬在從變形到斷裂的過程中,會因?yàn)檩d荷或者溫度等外部條件發(fā)生變化,進(jìn)而材料內(nèi)部產(chǎn)生微觀或者宏觀的缺陷,即材料將受到損傷,損傷的拓展會導(dǎo)致材料最終斷裂。本文所分析的球閥剪切油管斷裂就是因?yàn)榇笞冃握T發(fā)的材料內(nèi)部空洞和微裂紋的產(chǎn)生和拓展,被稱為“韌性損傷”[15]。Shear damage斷裂模型是基于單元積分點(diǎn)的等效塑性應(yīng)變,能夠預(yù)測由于金屬內(nèi)部空隙成核、成長、集結(jié)引起的損傷萌發(fā)。當(dāng)材料失效參數(shù)ω超過1時,則假定為材料失效。失效參數(shù)定義見式(3):
(3)
對剪切過程進(jìn)行簡化如圖2所示,只分析油管被剪切階段。本文所選用的球閥外徑D1為270 mm,球閥通孔直徑d1為165 mm,連續(xù)油管外徑D2為50.8 mm。連續(xù)油管損傷計(jì)算模型使用Joshnon-cook損傷準(zhǔn)則和本構(gòu)方程,采用Shear damage斷裂準(zhǔn)則,在分析油管時,連續(xù)油管網(wǎng)格采用六面體網(wǎng)格并細(xì)化剪切處連續(xù)油管網(wǎng)格,將球閥設(shè)置為剛體,其網(wǎng)格劃分如圖3所示。
圖2 球閥剪切模型Fig.2 Shearing model of ball valve
圖3 網(wǎng)格劃分Fig.3 Meshing
為了提高有限元模型計(jì)算的準(zhǔn)確性,進(jìn)行了一系列的參數(shù)設(shè)置比較,以確定更符合實(shí)際情況的仿真環(huán)境,相關(guān)參數(shù)設(shè)置以及材料模型設(shè)置參數(shù)見表1~2。
在球閥剪切過程中,隨著球閥旋轉(zhuǎn),球閥推動油管到閥座上,油管被壓扁,然后被剪斷,先后經(jīng)歷彈性變形、塑形變形、裂紋產(chǎn)生與擴(kuò)展、材料斷裂分離等階段。不同剪切階段的應(yīng)力云圖如圖4所示。在無外界因素影響下球閥剪切油管所需要的力矩隨剪切時間的變化曲線如圖5所示。從圖4(a)可知,油管在被球閥剪切過程中,油管先被球閥壓扁,此時剪切力矩增大,如圖5中的O-A段;壓扁后的油管分為2層,2層之間有空隙和連接,在擠壓油管的過程中,到達(dá)時間點(diǎn)A時,第1層油管開始被剪切,此時從非剪切階段過渡到剪切段,金屬從塑性變形到塑形斷裂,剪切力矩在A處會有一定下降;其后隨著剪切深度的增加,剪切力矩持續(xù)增加直到B點(diǎn),剪切如圖4(b)所示;經(jīng)過B點(diǎn)之后,第1層油管完全被剪斷,此時球閥刃口到達(dá)2層管之間連接部分,如圖4(c)所示,剪切面積大幅度減小,剪切力矩有一個快速的下降,如圖5中B-C段所示;隨后,球閥開始對第2層管剪切,剪切力矩的變化如C-D段所示,其剪切如圖4(d)所示。
表1 仿真模型參數(shù)Table 1 Parameters of simulation model
表2 球閥及CT90連續(xù)油管材料參數(shù)Table 2 Materials parameters of ball valve and CT90 coiled tubing
根據(jù)以上分析可以得到球閥剪切連續(xù)油管中剪切力矩的變化規(guī)律:油管在被剪斷的過程中,球閥的剪切力矩一直發(fā)生變化,油管被壓扁后分為2層,剪切過程出現(xiàn)2次峰值,故在后續(xù)對剪切油管的分析中,我們只需要提取出每條曲線的2個峰值B,D點(diǎn)剪切力矩大小進(jìn)行比較,就可以知道不同剪切條件下,剪切油管所需要的最大剪切力矩變化情況。
圖4 剪切過程中應(yīng)力云圖Fig.4 Stress nephogram during shearing process
圖5 標(biāo)準(zhǔn)剪切過程中剪切力矩隨時間變化Fig.5 Change of shearing torque with time during standard shearing process
連續(xù)油管由于作業(yè)情況不同,其拉力載荷會發(fā)生變化,為了比較不同拉力載荷對球閥剪切力矩的影響,在連續(xù)油管下端添加75,100,125 kN拉力,圖6為不同拉力載荷作用下球閥剪切力矩隨著剪切時間的變化。由圖6可以看出,當(dāng)油管承受拉力時,無論是第1層油管剪切過程還是第2層油管剪切過程,所需要的剪切力矩都減小,尤其是在剪切第2層油管過程中由于拉力作用,剪切力矩下降約30%。
圖6 不同拉力下剪切力矩比較Fig.6 Comparison of shearing torque under different tensions
分析可知,當(dāng)連續(xù)油管在井下由于自身重力或其他因素受到拉力時,球閥在剪切過程中需要的剪切力矩會更小,可以適當(dāng)減小設(shè)計(jì)安全系數(shù)。另外,由于剪切力矩大幅度減小是在第2層管剪切過程中,而第1層管剪切所需力矩依然很大,所以在設(shè)計(jì)球閥傳動機(jī)構(gòu)時,必須保證球閥剪切過程中前1/2弧度處于最大力矩傳動。
由于球閥驅(qū)動機(jī)構(gòu)不同,球閥會有大小不同的剪切速度,給球閥剪切施加不同的剪切速度進(jìn)行仿真計(jì)算。不同剪切速度下球閥剪切力矩隨時間變化曲線如圖7所示。由圖7可以明顯看出,在剪切第1層管的過程的,剪切速度越大,剪切力矩越大,但是在第2段剪切過程中即剪切第2層管時,剪切速度越大,剪切力矩反而越小。
圖7 不同剪切速度下球閥剪切力矩變化Fig.7 Change of shearing torque of ball valve under different shearing velocities
由分析可以得出,在保證剪切時間滿足緊急情況要求的同時,盡量減小球閥關(guān)閉速度,避免產(chǎn)生過大的剪切力矩。
水下測試樹型號不同,會有不同尺寸的球閥,為研究球閥大小對球閥剪切力的影響,選用了外徑240,255,270,285,300 mm的球閥尺寸進(jìn)行分析。得到剪切力矩變化如圖8所示。從圖8可以看出,球閥直徑越大,切斷連續(xù)油管所需要的最大剪切力越大,在剪切第1層管和第2層管中所需要的最大剪切力差別不大。
圖8 不同球閥直徑下剪切力矩變化情況Fig.8 Change of shearing torque under different diameters of ball valve
由分析可知,在相同內(nèi)徑下,球閥外徑越大,剪切連續(xù)油管所需要的剪切力矩就越大,所以在設(shè)計(jì)時應(yīng)盡量選用較小的球閥外徑。
為了驗(yàn)證剪切球閥剪切連續(xù)油管有限元模型的可靠性,同Seranton等[10]所做的球閥剪切連續(xù)油管試驗(yàn)進(jìn)行對比驗(yàn)證。
保證分析環(huán)境與試驗(yàn)環(huán)境一致,剪切球閥及連續(xù)油管主要參數(shù)如表3所示。通過使用同試驗(yàn)一致的4種厚度連續(xù)油管進(jìn)行有限元分析,得到不同的油管斷口,并測量每個斷口的長徑和短徑,得到不同厚度剪切分析下的長短徑比,并與試驗(yàn)進(jìn)行對比,計(jì)算其相對偏差δ,進(jìn)而判斷分析模型可靠性。
表3 剪切球閥及連續(xù)油管結(jié)構(gòu)參數(shù)Table 3 Structure parameters of shearing ball valve and coiled tubing
4組連續(xù)油管剪切之后的斷口形狀和仿真得到的斷口形狀如圖9所示,得到的CT90連續(xù)油管被剪斷時,連續(xù)油管斷口為橢圓狀,其長徑b和短徑a比值大小如表4所示。
所選用的有限元計(jì)算模型與剪切試驗(yàn)獲得的油管剪切變形程度數(shù)據(jù)相對偏差在5% 內(nèi),符合數(shù)據(jù)偏差要求。試驗(yàn)結(jié)果驗(yàn)證了球閥剪切連續(xù)油管分析模型的可靠性。
表4 不同管厚斷口仿真和試驗(yàn)結(jié)果長短徑之比Table 4 Simulation and test results of ratios of long radius to short radius for fracture with different pipe thicknesses
1)本文建立了球閥剪切連續(xù)油管模型,得到了球閥剪切連續(xù)油管所需剪切力變化規(guī)律:連續(xù)油管被壓扁后分為2層及中間連接,剪切過程中產(chǎn)生2個剪切力峰值,通過分析這2個峰值大小變化情況可以得出油管在不同剪切情況下的剪切性能,有效驗(yàn)證了球閥剪切連續(xù)油管的機(jī)理。
2)對連續(xù)油管承受拉力、球閥剪切速度、球閥大小3個因素進(jìn)行了分析,得出結(jié)果:承受拉力越大,球閥直徑越小、球閥剪切速度越小,球閥切斷連續(xù)油管所需要的剪切力越小。同時根據(jù)實(shí)際剪切過程,提出相應(yīng)的解決辦法,用以提高剪切性能,為深水測試樹的國產(chǎn)化研制提供參考依據(jù),也為其現(xiàn)有產(chǎn)品在現(xiàn)場作業(yè)提供重要的指導(dǎo)數(shù)據(jù)。
3)將球閥剪切連續(xù)油管有限元模型計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行比較,兩者的偏差較小,表明所建立計(jì)算模型的可靠性,因此,在水下測試樹產(chǎn)品剪切性能測評和深水油氣測試作業(yè)施工設(shè)計(jì)等過程中,可以直接采用本文所建立的數(shù)值仿真模型和方法進(jìn)行分析與計(jì)算,從而可有效降低實(shí)驗(yàn)成本和提高施工設(shè)計(jì)效率。