李 濤，丁愛芹，陳 強，李益良，裴曉含，張 濤

中國石油勘探開發(fā)研究院采油采氣裝備研究所 （北京 100083）

膨脹管技術(shù)是一種由特殊材料制成、具有良好塑性的金屬鋼管，下入井內(nèi)后進行膨脹，在冷作硬化效應(yīng)下，管材力學(xué)性能得到提高[1-3]，在管柱設(shè)計過程中，通常使用傳統(tǒng)的安全系數(shù)法，即將各因素的最大值（如外擠力）和名義值（如內(nèi)外徑）等固定值代入模型之中，得出一個確定的安全系數(shù)，但是受到塑性變形、殘余應(yīng)力等多重不確定因素影響，使用安全系數(shù)法設(shè)計膨脹管管柱存在較大的主觀隨意性。

通過分析不同膨脹率條件下膨脹管脹后抗外擠強度以及殘余應(yīng)力的變化規(guī)律，同時考慮到膨脹管膨脹過程中屈服強度、外徑、壁厚、橢圓度、殘余應(yīng)力等參數(shù)變化的隨機性，利用可靠性理論建立膨脹管外部載荷與失效概率的安全評價方法。

1 脹后膨脹管力學(xué)性能

通常認(rèn)為膨脹管膨脹屬于冷作硬化過程，脹后膨脹管屈服強度與抗拉強度提高，延展率下降，殘余應(yīng)力變大。膨脹率是描述膨脹管脹后尺寸的重要參數(shù)，對膨脹管力學(xué)性能以及可靠性有重要影響，利用有限元分析114mm×7mm規(guī)格碳鋼材料膨脹管在不同膨脹率條件下力學(xué)性能的變化規(guī)律。圖1為膨脹管塑性變形過程以及脹后等效應(yīng)力云圖。膨脹過程中在膨脹的部位即膨脹錐所在的位置出現(xiàn)等效應(yīng)力最大的紅色區(qū)，膨脹后應(yīng)力得到了釋放。

膨脹管膨脹過程中，在環(huán)向、軸向及徑向都將發(fā)生塑性變形，這種塑性變形往往是不均勻的，不均勻的變形在管體內(nèi)產(chǎn)生附加應(yīng)力，膨脹后會殘留在套管內(nèi)部形成殘余應(yīng)力[2-4]。在有限元模型上選擇一條路徑即z軸向，將模擬計算的結(jié)果數(shù)據(jù)映射到路徑上并繪出相關(guān)路徑圖，根據(jù)所得數(shù)據(jù)描繪膨脹后殘余應(yīng)力的分布情況。從圖2可以看出，不同膨脹率的膨脹管在膨脹過程中殘余應(yīng)力隨軸向的分布情況大體相同，殘余應(yīng)力隨膨脹率增加而變大；膨脹過程中殘余應(yīng)力隨軸向的分布情況相差很大,通常膨脹錐所在的位置殘余應(yīng)力值較大，膨脹后平均殘余應(yīng)力逐漸減小，應(yīng)力得到釋放。

圖2 不同膨脹率膨脹過程及脹后殘余應(yīng)力變化示意圖

膨脹管的失效方式一般有2種[4]，局部擠毀和整體屈服。利用有限元模擬膨脹管局部擠毀和整體屈服狀態(tài)，并獲取膨脹管屈服強度數(shù)據(jù)。從圖3中可以看到，套管擠毀后截面通常變?yōu)闄E圓形，局部或完全被擠扁。

圖3 膨脹管失效過程有限元分析

脹后膨脹管屈服強度受到2個因素的影響，一方面塑性變形過程中，膨脹管晶格扭曲、畸變，晶粒產(chǎn)生剪切、滑移，晶粒被拉長[5]，表面硬度增加；另一方面由于膨脹管膨脹過程中遵循體積不變原理[6]，隨著膨脹率的增加，膨脹管壁厚減少，強度降低。

相關(guān)試驗及有限元模擬分析表明[7-9]，隨著膨脹率的增加，脹后膨脹管強度逐漸增加，達到峰值后基本保持不變。變化規(guī)律如圖4所示。

圖4 膨脹管屈服失效變化規(guī)律

2 膨脹管可靠性模型

膨脹管的塑性變形會導(dǎo)致外徑、壁厚不均勻度、橢圓度以及殘余應(yīng)力等參數(shù)發(fā)生變化，同時由于膨脹率不同，這些參數(shù)的變化差異較大。實驗與模型計算表明，膨脹管抗外擠強度服從正態(tài)分布[8-10]。

2.1 膨脹管抗外擠強度的概率分布

針對膨脹管導(dǎo)致的脹后膨脹管強度變化，利用Klever-Tomano抗外擠強度模型[10-11]作為膨脹管抗外擠失效概率分布函數(shù)，其表達式為：

式中：PcR為脹后膨脹管抗外擠強度，MPa；Pe為理想圓管的極限彈性擠毀壓力，MPa；Pv為理想圓管的極限屈服擠毀壓力，MPa；E為彈性模量，MPa；γ為泊松比；D為膨脹管外徑，mm；t為膨脹管壁厚，mm；σy為屈服強度，MPa；ov表示膨脹管橢圓度，mm；ec表示膨脹管壁厚不均勻度，mm；rs為膨脹管殘余應(yīng)力，MPa；hn為膨脹管應(yīng)力應(yīng)變曲線形狀系數(shù)，取值0.009。

2.2 膨脹管抗內(nèi)壓的概率分布

膨脹管抗內(nèi)壓強度有3種表現(xiàn)形式[12-14]，屈服強度、塑性破裂強度和裂紋擴展破裂強度。如不考慮腐蝕、氫脆等因素影響，膨脹管在內(nèi)壓作用下的主要失效形式是塑性破裂。研究表明：Klever-Stewart[10,12]提出的塑性變形破壞內(nèi)壓強度模型精度較高，可以作為膨脹管抗外壓概率分布函數(shù)，表達式為：

式中：PiR為膨脹管抗內(nèi)壓強度，MPa；fu為膨脹管拉伸屈服強度，MPa；kα為內(nèi)壓屈服系數(shù),取值2.0;aN位膨脹管缺陷深度，mm，通常為5%膨脹管壁厚；n為材料應(yīng)力應(yīng)變強度硬化因子，無因次；D為膨脹管外徑，mm；tmin為膨脹管最小壁厚，mm。

3 膨脹管安全可靠性評價

根據(jù)前面建立的模型計算114mm×7mm規(guī)格膨脹管脹后抗外擠和抗內(nèi)壓強度的概率分布，可以得到膨脹管脹后抗外擠和抗內(nèi)壓強度服從正態(tài)分布規(guī)律，如圖5所示。

圖5 膨脹管脹后強度概率分布

同時利用模型計算不同載荷下膨脹管抗外擠與抗內(nèi)壓的失效概率，如表1所示。結(jié)果發(fā)現(xiàn)，膨脹管的失效概率與外部載荷并非正相關(guān)，不同膨脹率的膨脹管在同一載荷條件下失效概率有較大差異，用可靠性計算的失效概率評價指標(biāo)可以為膨脹管強度安全系數(shù)的科學(xué)設(shè)計提供依據(jù)。

4 結(jié)論

1）膨脹管脹后強度受到膨脹率與壁厚變化量雙重影響，通常隨膨脹率增加變大，達到峰值后基本保持不變。不同膨脹率的實體管在膨脹過程中殘余應(yīng)力隨軸向的分布情況大體相同，只是數(shù)值大小不同。

2）采用隨機抽樣法模擬膨脹管強度的隨機分布規(guī)律，發(fā)現(xiàn)受到橢圓度、壁厚不均勻度以及殘余應(yīng)力等因素影響，膨脹管抗外擠與抗內(nèi)壓強度服從正態(tài)分布規(guī)律。

表1 不同膨脹率的膨脹管失效概率與外部載荷數(shù)據(jù)統(tǒng)計

3）模型計算結(jié)果表明不同膨脹率條件的膨脹管失效概率不同，論文研究為膨脹管強度安全系數(shù)的科學(xué)設(shè)計提供依據(jù)。

[1]American Petroleum Institute.API Bulletin5C3-1994.Bulletin on formulas and calculations for casing,tubing,drill pipe and line pipe properties[S].Dallas：API Production Department,1994.

[2]張建兵，施太和，練章華.鉆井膨脹管膨脹過程中不均勻變形的試驗研究[J].石油機械，2004，32（7）：1-5.

[3]俞茂宏，何麗南，宋凌宇.雙剪應(yīng)力強度理論及推廣[J].中國科學(xué)（A），1985，28（12）：1113-1120.

[4]Lewis D B.Load and resistance factor design for oil country tubular goods[R].OTC7936,1995.

[5]Cunha J C,Demirdal B,Gui P.Use of quantitative risk analysisfor uncertainty quantification on drilling operations-review and lessons learned[R].SPE94980,2005.

[6]王國榮,劉清友,何霞.套管可靠性壽命預(yù)測[J].天然氣工業(yè),2002,22(5)：53-55.

[7]閆相禎,高進偉,楊秀娟.用可靠性理論解析API套管強度的計算公式[J].石油學(xué)報,2007,28(1)：122-126.

[8]何水清,王善.結(jié)構(gòu)可靠性分析與設(shè)計[M].北京：國防工業(yè)出版社,1993：14-20.

[9]王建東,史交齊,林凱,等.套管抗擠強度的可靠性評價方法[J].石油工業(yè)技術(shù)監(jiān)督,2005,21(11)：36-38.

[10]International Organization for Standardization.ISO/TR10400-2007Petroleum and natural gas industries-equations and calculations for the properties of casing,tubing,drill pipe and line pipe used as casing or tubing[S].Geneva：ISO Copyright Office,2007.

[11]Klever F J,Tamano T.A new OCTG strength equation for collapse under combined loads[R].SPE90904,2006.

[12]Klever F J.Formulas for rupture,necking,and wrinkling for OCTG under combined loads[R].SPE102585,2006.

[13]Adams A J,Parfitt S H L,Reeves T B,et al.Casing system riskanalysis using structural reliability[R].SPE25693,1993.

[14]Liang Q J.Application of quantitative risk analysis to pore pressure and fracture gradient prediction[R].SPE77354,2002.