徐冠軍雷金晶馮定后鑫
1.中國石化江漢油田分公司采油工藝研究院機(jī)械采油研究所 湖北武漢 430000 2.中原油田普光分公司開發(fā)管理部 四川宣漢 636150 3.荊州市長江大學(xué)機(jī)械工程學(xué)院 湖北荊州 434023 4.中原油田普光分公司采氣廠 四川宣漢 636150
壓縮式封隔器膠筒坐封過程數(shù)值分析
徐冠軍1雷金晶2馮定3后鑫4
1.中國石化江漢油田分公司采油工藝研究院機(jī)械采油研究所 湖北武漢 430000 2.中原油田普光分公司開發(fā)管理部 四川宣漢 636150 3.荊州市長江大學(xué)機(jī)械工程學(xué)院 湖北荊州 434023 4.中原油田普光分公司采氣廠 四川宣漢 636150
本文介紹了壓縮式封隔器密封的工作原理,總結(jié)國內(nèi)井下封隔工具的應(yīng)用現(xiàn)狀,指出封隔器研究的發(fā)展方向;在此基礎(chǔ)上研究了有限元分析的基本理論和模型參數(shù)的確定方法,建立膠筒的有限元分析模型,用ANSYS軟件對(duì)壓縮式封隔器膠筒坐封過程進(jìn)行有限元分析,分析在150℃溫度時(shí)膠筒的軸向壓縮位移與坐封載荷的關(guān)系,以及膠筒接觸壓力和坐封載荷的關(guān)系,分析結(jié)果與實(shí)際情況基本吻合,為壓縮式裸眼封隔器的在裸眼井的應(yīng)用提供了的理論指導(dǎo)。
壓縮式封隔器;擴(kuò)張式封隔器;封隔器膠筒;有限元分析
壓縮式封隔器的工作原理:當(dāng)膠筒承受軸向載荷時(shí),封隔器膠筒將產(chǎn)生徑向大變形,使膠筒與套管之間產(chǎn)生接觸壓力,藉此封隔環(huán)空,隔絕產(chǎn)層,以控制產(chǎn)(注)液,保護(hù)套管。膠筒所用材料為橡膠,橡膠最顯著的特性就是它的超彈性,即在很小的力的作用下就能產(chǎn)生很大的變形,正是這一性質(zhì)使得橡膠材料被制成密封元件并在工程中得到廣泛應(yīng)用。
但是,國內(nèi)還沒有專門用于深水井下作業(yè)的相關(guān)封隔工具,因此,研究封隔器膠筒的材料以及研究用于惡劣環(huán)境下的井下措施工具以及技術(shù),是國內(nèi)各大研究單位以及相關(guān)高校的當(dāng)務(wù)之急。
橡膠材料在較短時(shí)間內(nèi)及恒定的環(huán)境溫度下通常被視為各向同性不可壓縮材料,其應(yīng)變能密度函數(shù)W是變形張量不變量的函數(shù),即:
式中:λ1、λ2和λ3是3個(gè)主伸長率。
Rivlin用唯象法從純數(shù)學(xué)角度出發(fā)導(dǎo)出應(yīng)變能密度函數(shù)的最一般的形式為:
式中:Cijk為Rivlin材料常數(shù)。
橡膠材料可以近似認(rèn)為是不可壓縮的,因此I3=1,式(5)可簡化為:
取式(6)的前兩項(xiàng)得到Mooney-Rivlin模型:
其中C10和C01為Rivlin材料常數(shù),均為正定常數(shù)。
2.2 Mooney-Rivlin模型材料常數(shù)的確定
指視野中垂直范圍內(nèi)所能清晰觀察到的物象界限。當(dāng)不用微調(diào)就能看清楚一個(gè)物像的平面及上下結(jié)構(gòu),這個(gè)能夠看清的厚度即焦點(diǎn)深度。
對(duì)兩參數(shù)Mooney-Rivlin模型,Mooney-Rivlin材料模型的應(yīng)力-應(yīng)變方程可表示為:
這個(gè)方程可用σ/2(λ-λ-2)對(duì)1/λ作圖。在1/λ=1時(shí),相應(yīng)值為C10+C01并且直線的斜率為C01,初始剪切模量與材料常數(shù)的關(guān)系是G=2(C10+C01)。如果材料被假定為不可壓縮的,那么初始拉伸模量是E=6(C10+C01),G=2(C10+C01)和楊氏模量的值近似等于3倍的剪切模量。
利用橡膠材料的硬度與C10和C01存在一定的關(guān)系,根據(jù)橡膠材料硬度與彈性模量的實(shí)驗(yàn)數(shù)據(jù),經(jīng)擬合可得:
利用硬度計(jì)可測得橡膠試樣的硬度HA,通過式(9)可算出橡膠彈性模量E,根據(jù)公式E=6(C10+C01)可以得到C10+C01的值,通過公式d=(1-2μ)/(C10+C01)可以測出模型參數(shù)d的值。利用經(jīng)驗(yàn)公式C01=0.5C10求解下面的方程就可得到C10。
由實(shí)驗(yàn)測得在常溫25℃和高溫150℃下膠筒的邵氏A硬度分別為94和91,通過公式(8)、(9)即可計(jì)算出的膠筒的Mooney-Rivlin模型材料參數(shù)。
封隔器采用兩膠筒形式,上、下膠筒、中心管和套管結(jié)構(gòu)。上下膠筒力學(xué)參數(shù)和幾何尺寸完全相同,兩膠筒由隔環(huán)隔開。由于封隔器中心管、膠筒、套管以及膠筒所受載荷均為軸對(duì)稱分布,故取過軸線的剖面建立有限元計(jì)算模型, 膠筒采用四節(jié)點(diǎn)四邊形超彈單元PLANE182,中心管、套管和隔環(huán)采用四節(jié)點(diǎn)四邊形等參數(shù)單元PLANE42,用面面接觸單元CONTA171和TARGE169模擬膠筒與套管、中心管和剛性隔環(huán)間的接觸,并考慮接觸時(shí)的摩擦作用,通過實(shí)驗(yàn)測定膠筒與鋼材之間的摩擦因數(shù)為0.3,鋼與鋼之間的摩擦因數(shù)為0.1,中心管和套管上下兩端采用固定約束,橡膠筒在上端隔環(huán)的作用下壓縮膨脹,與套管接觸,隨著坐封壓力不斷增加,膠筒與套管的接觸壓力也隨之增大,從而起到密封套管內(nèi)上下環(huán)空壓差的作用。
建立膠筒有限元模型,對(duì)膠筒進(jìn)行非線性接觸分析。膠筒具有不可壓縮性,因此取膠筒泊松比為0.4996。
4.1 坐封過程中接觸壓力分布規(guī)律
為模擬封隔器坐封過程中膠筒與套管內(nèi)壁接觸情況,穩(wěn)定溫度150℃不變,分步向膠筒上端隔環(huán)施加軸向荷載,觀察膠筒的變形情況。
由于膠筒與套管內(nèi)壁接觸壓力和摩擦阻力的存在,2個(gè)膠筒所受到的軸向壓縮載荷由上至下依次降低,所以盡管上下膠筒的材料和幾何形狀完全相同,但是下膠筒的變形、與套管的接觸壓力均小于上膠筒,坐封完成以后,上膠筒的接觸壓力比較大,下膠筒的接觸壓力比較小。顯然,隨著封隔器坐封載荷的增加,上膠筒與套管的接觸壓力明顯高于下膠筒的,所以上膠筒在密封套管環(huán)空上下壓差中起到關(guān)鍵作用。
4.2 坐封過程中軸向載荷與膠筒變形的分布規(guī)律
根據(jù)ANSYS分析輸出結(jié)果可以得到在150℃時(shí)不同坐封力下膠筒的軸向壓縮位移和軸向載荷的關(guān)系。穩(wěn)定溫度在150℃時(shí),依次施加10t、12t、13t、15t、17t、18t坐封力時(shí)膠筒的軸向壓縮位移分別為23.095mm、24.319mm、24.564mm、24.905mm、25.086mm、25.186mm。
實(shí)驗(yàn)發(fā)現(xiàn),隨著軸向載荷的逐步增大,膠筒軸向壓縮位移也漸漸增大,當(dāng)軸向載荷小于12t時(shí),膠筒軸向壓縮位移隨軸向載荷增大的梯度比較大,當(dāng)軸向載荷大于12t時(shí),膠筒軸向壓縮位移增大梯度慢慢減小。因此,當(dāng)軸向載荷過大時(shí),膠筒軸向壓縮位移的變化不是很大。
從數(shù)值模擬可以看出,封隔器的兩個(gè)膠筒與套管的接觸壓力分布為:下膠筒接觸壓力較小,上膠筒的接觸壓力較大,隨著軸向載荷的增
()()大,軸向壓縮量也增大,同時(shí)也可以發(fā)現(xiàn),隨著坐封力的增大,膠筒與套管接觸長度逐漸增加,膠筒外表柱面部份徑向變形受限制,膠筒內(nèi)表面變形如外表一樣向外鼓,當(dāng)載荷增加時(shí)膠筒被壓扁并在最后壓實(shí)。并且隨著軸向載荷的逐步增大,膠筒軸向壓縮位移一開始是漸漸增大,增大的梯度也比較大。當(dāng)軸向載荷繼續(xù)增加時(shí),膠筒軸向壓縮位移增大梯度會(huì)慢慢減小。因此,當(dāng)軸向載荷過大時(shí),膠筒軸向壓縮位移的變化不是很大,為了達(dá)到良好的密封,應(yīng)該計(jì)算合適的坐封載荷。
[1]曾憲平.對(duì)壓縮式封隔器膠筒變形分析[J].石油鉆采工藝,1983,13(4):65-68.
[2]楊秀娟,楊恒林.液壓封隔器膠筒座封過程數(shù)值分析[J].石油大學(xué)學(xué)報(bào),2003,05-0084-04.
[3]劉永輝.封隔器膠筒性能有限元分析及結(jié)構(gòu)改進(jìn)[D]西南石油大學(xué),2006
[4]楊曉翔.非線性橡膠材料的有限單元法[M].北京:石油工業(yè)出版社,1999-70-82.
[5] Ding Feng, Yongxin Yuan, Guanjun Xu. Finite Element Analysis Of The Packer Rubbers On Sealing Process[J].ISBN,2010