基于Monte Carlo方法的SCV加熱s盤管可靠性評價

鄧 嬌 侯 磊 李麗敏 文亞成

1. 中國石油大學(xué)(北京)機械與儲運工程學(xué)院, 北京 102249;2. 中國市政工程中南設(shè)計研究總院, 湖北 武漢 430010;3. 大慶油田采油六廠, 黑龍江 大慶 163114

0 前言

浸沒燃燒式氣化器(Submerged Combustion Vaporizer,SCV)具有結(jié)構(gòu)緊湊、節(jié)省空間、換熱效率高、排煙溫度低、節(jié)省能源等優(yōu)點,是LNG接收站氣化外輸系統(tǒng)的重要設(shè)備,在冬季作為接收站的氣化設(shè)備,在其他季節(jié)作為備用氣化器,其可靠性直接影響氣化外輸系統(tǒng)的安全運行。SCV的組成復(fù)雜,調(diào)研結(jié)果表明,其他部件均可在調(diào)控之后平穩(wěn)工作,而加熱盤管所處的水浴環(huán)境pH值為4～6,會發(fā)生腐蝕現(xiàn)象,影響SCV的可靠性。目前國內(nèi)外對于SCV的結(jié)構(gòu)優(yōu)化設(shè)計和操作參數(shù)設(shè)置進(jìn)行了研究,未涉及其可靠性評價方面[1]。本文構(gòu)建SCV可靠性評價模型,建立加熱盤管的極限狀態(tài)方程,計算SCV可靠度,研究加熱盤管腐蝕隨運行時間和運行壓力的變化規(guī)律,對站場SCV的運行和維修具有指導(dǎo)意義。

1 SCV結(jié)構(gòu)和工作原理

圖1 SCV工作原理

2 可靠性評價模型的建立

為評價SCV加熱盤管產(chǎn)生腐蝕缺陷時能否在規(guī)定的壓力下繼續(xù)運行,利用應(yīng)力-強度干涉理論和管道腐蝕剩余強度評估標(biāo)準(zhǔn)建立可靠性評價模型。SCV加熱盤管的管材是316 L不銹鋼,適合于低合金鋼剩余強度評價的方法主要有五種[4-11]：美國機械工程師協(xié)會頒布的ASME B 31.G-1984適用于管道只有單一的軸向腐蝕缺陷；ASME B 31.G-1991在上一版的基礎(chǔ)上進(jìn)行了改進(jìn),對流變應(yīng)力、鼓脹系數(shù)和腐蝕缺陷面積三個影響因素做了修正,計算結(jié)果更加精確；ASME B 31.G-2009將評價方法分成四個等級,根據(jù)不同的腐蝕情況和收集的數(shù)據(jù)情況采用不同的評價方法,使得計算結(jié)果與實際情況更加符合；美國運輸部提出的RSTRENG標(biāo)準(zhǔn)需要非常詳細(xì)的腐蝕缺陷參數(shù),腐蝕參數(shù)主要包括腐蝕區(qū)域的形狀、長度和寬度；SY/T 6151《鋼制管道管體腐蝕損傷評價方法》將管道的腐蝕程度進(jìn)行層次劃分,分類進(jìn)行評定,這種方法也需要詳細(xì)且精確的管道腐蝕數(shù)據(jù)。由于SCV加熱盤管只承受內(nèi)壓,現(xiàn)場測得的腐蝕數(shù)據(jù)為腐蝕長度和腐蝕深度,從載荷類型、最佳適用材料和腐蝕數(shù)據(jù)的獲取情況三方面考慮,采用ASME B 31.G-2009的1級評價方法進(jìn)行SCV加熱盤管腐蝕剩余強度評價。

2.1 建立極限狀態(tài)方程

ASME B 31.G-2009中要求:當(dāng)失效壓力大于或等于操作壓力乘以1個合適的安全系數(shù)時,缺陷可以接受[9]。為了評估加熱盤管的可靠性,建立加熱盤管極限狀態(tài)函數(shù)[12]:

z=pf-ps

(1)

式中:z為加熱盤管極限狀態(tài)函數(shù)(z>0,加熱盤管可靠;z=0,加熱盤管處于極限狀態(tài);z<0,加熱盤管失效);pf為傳熱管失效時的壓力,MPa;ps為加熱盤管的最大安全運行壓力,MPa。

(2)

ps=S×p0

(3)

式中:SF為加熱盤管發(fā)生失效時的環(huán)向應(yīng)力,MPa;D為管道外徑,mm;t為管道壁厚,mm;S為安全系數(shù),取S=1.25;p0為加熱盤管運行壓力,MPa。

加熱盤管的失效應(yīng)力按下面公式計算:

(4)

式中:d為腐蝕深度,mm;M為鼓脹系數(shù);Sflow為流變應(yīng)力,MPa。

(5)

式中:L為腐蝕長度,mm。

(6)

流變應(yīng)力表示為：

Sflow=SMYS+69

(7)

式中:SMYS為最小屈服強度,MPa。

對于準(zhǔn)穩(wěn)態(tài)腐蝕,腐蝕的長度和深度隨時間擴大,T時刻腐蝕長度和腐蝕深度分別為：

L=L0+Va(T-T0)

(8)

d=d0+Vr(T-T0)

(9)

式中:T為現(xiàn)腐蝕檢測時間，a；T0為腐蝕檢測初始時刻,a;L0為腐蝕檢測初始時刻的腐蝕長度,mm;d0為腐蝕檢測初始時刻的腐蝕深度,mm;Va為長度腐蝕速率,mm/a;Vr為深度腐蝕速率,mm/a。

將式(2)～(6)帶入式(1),可得：

(10)

2.2 Monte Carlo方法仿真過程

Monte Carlo方法通常用于計算有隨機變量函數(shù)的函數(shù)值,通過對獨立隨機變量進(jìn)行抽樣實驗或隨機模擬計算,可用于隨機函數(shù)服從任意分布的函數(shù)問題。其收斂速度與隨機變量的維數(shù)無關(guān),只要所取的循環(huán)次數(shù)足夠多,就可以保證概率計算的準(zhǔn)確性和全面性[13-16]。加熱盤管的初始腐蝕尺寸、工況載荷、機械強度、軸向及徑向腐蝕增長速率等參數(shù)具有隨機性,本文采取Monte Carlo方法解決加熱盤管可靠性評價過程中的隨機性問題,定量分析加熱盤管的可靠性,仿真步驟如下:

1)確定模型中所有參數(shù)的分布規(guī)律,對腐蝕缺陷參數(shù)和傳熱管基本參數(shù)進(jìn)行隨機抽樣,抽樣次數(shù)N為106。

2)將抽樣所得數(shù)據(jù)代入式(10)中,計算z的值，得到傳熱管失效壓力與最大運行壓力的差值Δ p=z。

3)若z≤0,則記為1；若z>0,則記為0；統(tǒng)計失效次數(shù)n。

4)計算加熱盤管失效概率PF=n/N。

2.3 模型參數(shù)的確定

加熱盤管腐蝕剩余強度的極限狀態(tài)方程中包含的參數(shù)有管道外徑D、管道壁厚t、管道運行壓力p0、管道屈服強度SMYS、腐蝕深度d0、腐蝕長度L0、深度腐蝕速率Vr、長度腐蝕速率Va。為簡化計算,暫不考慮變量間的相關(guān)性,假定各參數(shù)相互獨立。

這些參數(shù)往往不是定值且服從一定的分布規(guī)律,通過調(diào)研文獻(xiàn)可知[17],工況載荷、盤管尺寸和腐蝕區(qū)域尺寸的分布規(guī)律均可用正態(tài)分布描述。選取所有變量為隨機變量且均服從正態(tài)分布。

3 計算實例

3.1 SCV加熱盤管剩余壽命計算

以某LNG接收站SCV加熱盤管為例,盤管采用316 L不銹鋼[18-20],通過現(xiàn)場調(diào)研,加熱盤管參數(shù)見表1。

表1加熱盤管參數(shù)

參數(shù)均值標(biāo)準(zhǔn)差分布形式外徑/mm3175001正態(tài)分布壁厚/mm2410001正態(tài)分布內(nèi)壓/MPa251正態(tài)分布屈服強度/MPa1702正態(tài)分布初始腐蝕深度/mm0015正態(tài)分布初始腐蝕長度/mm001正態(tài)分布深度腐蝕速率/(mm·a-1)0043001正態(tài)分布長度腐蝕速率/(mm·a-1)0801正態(tài)分布

對表1中的參數(shù)進(jìn)行抽樣,計算加熱盤管的失效概率。圖2為運行10、12、15、18 a的失效概率分布曲線。失效概率隨著使用年限的延長而增大,使用年限增長,腐蝕長度和腐蝕深度增加,加熱盤管承壓面積減少,失效概率增大。圖3為加熱盤管累積失效概率隨時間的變化曲線,累積失效概率隨運行時間的增加逐漸增加,運行前期累積失效概率增加緩慢,運行至某一時間,累積失效概率開始迅速增加。加熱盤管運行前10 a,累積失效概率較低,可靠度高,運行10 a時的累積失效概率為0.002 99；在運行10～12 a時,累積失效概率開始緩慢上升,可靠度逐漸降低,建議縮短檢修周期；從運行12 a開始,累積失效概率開始快速上升,可靠度迅速下降,應(yīng)密切檢查加熱盤管的運行情況,及時檢修；運行18 a時,其累積失效概率為0.411。

a)投入使用10 a的失效概率曲線

b)投入使用12 a的失效概率曲線

c)投入使用15 a的失效概率曲線

d)投入使用18 a的失效概率曲線圖2 加熱盤管失效概率分布曲線

圖3 加熱盤管累積失效概率隨時間的變化曲線

3.2 運行壓力對可靠性的影響

在加熱盤管的可靠性影響因素中,穩(wěn)定的腐蝕環(huán)境使得腐蝕參數(shù)變化較小,而運行壓力隨著工況的變化而變化,因此主要分析運行壓力對加熱盤管可靠度的影響。

隨著壓力的變化,加熱盤管失效概率發(fā)生變化。以盤管運行10 a為例,選取壓力分別為20、23、25 MPa,失效概率分布曲線見圖4。隨著壓力增大,失效概率曲線向坐標(biāo)橫軸負(fù)方向移動,失效概率增加。不同運行壓力下加熱盤管累積失效概率隨時間的變化見圖5,運行壓力為20～22 MPa時,累積失效概率隨壓力的變化較小,20 a內(nèi)累積失效概率均較低；運行壓力為23～26 MPa時,累積失效概率隨壓力的增加變化較大,第20 a不同運行壓力下累積失效概率均較高。運行20 a,運行壓力從22 MPa增加到26 MPa,累積失效概率從0.1上升至0.8。同時,當(dāng)運行時間較短時,壓力的變化對累積失效概率的影響較小,但當(dāng)運行時間較長時,累積失效概率隨壓力的變化較大。運行壓力從20 MPa增加到26 MPa,運行8 a,累積失效概率均接近0;運行20 a,累積失效概率大約從0上升至0.8。

圖4 不同運行壓力下加熱盤管失效概率曲線

圖5 不同運行壓力下加熱盤管累積失效概率隨時間的變化曲線

4 結(jié)論

1)依據(jù)應(yīng)力-強度干涉理論建立SCV加熱盤管的極限狀態(tài)方程,按照ASME B31.G-2009標(biāo)準(zhǔn)的流程計算加熱盤管的失效應(yīng)力。采用Monte Carlo方法對加熱盤管極限狀態(tài)函數(shù)中的參數(shù)進(jìn)行隨機抽樣,計算加熱盤管可靠度和預(yù)測加熱盤管剩余使用壽命。

2)對某LNG接收站SCV加熱盤管進(jìn)行腐蝕剩余壽命預(yù)測,隨著使用年限增加,加熱盤管失效概率增加,且增加速率逐漸加快。在盤管運行10 a時,其累積失效概率為0.002 99,可靠度為0.997;運行18 a時,其累積失效概率為0.411,可靠度為0.589。若以0.8為可靠度下限,則盤管使用壽命為14 a。

3)隨著盤管運行壓力的上升,其可靠度下降。當(dāng)運行年限較短時,運行壓力對盤管可靠度影響較小,當(dāng)運行年限較長時,盤管可靠度隨運行壓力的增加迅速上升。運行20 a,壓力20 ～22 MPa,加熱盤管可靠度高于0.9,壓力23 ～26 MPa,加熱盤管可靠度從0.8下降至0.2。

[1] 于國杰.LNG沉浸式燃燒型氣化器數(shù)值模擬[D].大連:大連理工大學(xué),2009:21-25.

Yu Guojie. Numerical Simulation of LNG Submerged-Combustion Vaporizer [D]. Dalian: Dalian University of Technology, 2009: 21-25.

[2] 劉世俊,郭 超.浸沒燃燒法在LNG氣化工藝中的應(yīng)用[J].城市燃?xì)?2015,12(490):22-27.

Liu Shijun, Guo Chao. Application of Immersion Combustion Method in LNG Gasification Process [J]. Urban Gas, 2015, 12 (490): 22-27.

[3] 程云東.浸沒燃燒式氣化器的設(shè)計缺陷及改進(jìn)方法[J].油氣儲運,2014,33(3):292-294.

Cheng Yundong. Design Defects and Modification of Submerged Combustion Vaporizer [J]. Oil & Gas Storage and Transportation,2014, 33 (3): 292-294.

[4] 俞樹榮,李建華,李淑欣,等.埋地管道腐蝕剩余壽命預(yù)測概率模型[J].中國安全科學(xué)學(xué)報,2008,18(6):11-15.

Yu Shurong, Li Jianhua, Li Shuxin, et al. Probability Model for the Prediction of Corrosion Remaining Life of Underground Pipelines [J]. China Safety Science Journal, 2008, 18 (6): 11-15.

[5] 帥 健.腐蝕管線的剩余壽命預(yù)測[J].中國石油大學(xué)學(xué)報:自然科學(xué)版,2003,27(4):91-93.

Shuai Jian. Prediction Method for Remaining Life of Corroded Pipelines [J]. Journal of the University of Petroleum, China: Edition of Natural Science, 2003, 27 (4): 91-93.

[6] 黃亞暉,周麗麗.B31.G腐蝕管道剩余強度評價比較分析[J].天然氣與石油,2012,30(4):62-66.

Huang Yahui, Zhou Lili. Evaluation on Residual Strength of Corroded B31. G Pipelines [J]. Natural Gas and Oil, 2012, 30 (4): 62-66.

[7] 何東升,郭 簡,張 鵬.腐蝕管道剩余強度評價方法及其應(yīng)用[J].石油學(xué)報,2007,28(6):125-128.

He Dongsheng, Guo Jian, Zhang Peng. Assessment Method for Remaining Strength of Corroded Pipeline and Its Application [J]. Acta Petrolei Sinica, 2007, 28 (6): 125-128.

[8] Kiefner J F, Vieth P H. A Modified Criterion for Evaluating the Remaining Strength of Corroded Pipe [M]. Columbus: Battelle Memorial Institute, 1989: 63-74.

[9] American Society of Mechanical Engineers. Manual of Det-ermining the Remaining Strength of Corroded Pipelines: ASME B31 G-1984 [S]. New York: ASME, 1984: 4-7.

[10] American Society of Mechanical Engineers. Manual of Det-ermining the Remaining Strength of Corroded Pipelines: ASME B31 G-1991[S]. New York: ASME, 1991: 4-7.

[11] American Society of Mechanical Engineers. Manual for Det-ermining the Remaining Strength of Corroded Pipelines: ASME B31 G-2009[S]. New York: ASME, 2009: 8-15.

[12] Petroleum and Natural Gas Industries. Pipeline Transportation Systems-Reliability-Based Limit State Methods: ISO 16708[S]. Geneva: ISO, 2006: 325-327.

[13] 趙 磊.重要性抽樣Monte-Carlo法在概率斷裂力學(xué)分析中的應(yīng)用研究[D].杭州:浙江工業(yè)大學(xué),2013: 25-27.

Zhao Lei. Research on the Application of the Monto-Carlo Simulation Based on Importance Sampling Method in the Analysis of Probabilistic Fracture Mechanics [D]. Hangzhou: Zhejiang University of Technology, 2013: 25-27.

[14] 張仙風(fēng),呂志鵬.基于MATLAB的蒙特卡羅方法在可靠性設(shè)計中的應(yīng)用[J].裝備制造技術(shù),2006,(4):82-83.

Zhang Xianfeng, Lv Zhipeng. The Application of MATLAB-based Monte Carlo Method in Reliability Design [J]. Equipment Manufacturing Technology, 2006, (4): 82-83.

[15] 孫春梅,李 琴,黃志強,等.基于Monte Carlo方法的腐蝕管道可靠性分析[J].油氣儲運,2015,34(8):811-816.

Sun Chunmei, Li Qin, Huang Zhiqiang, et al. Reliability Analysis of Corroded Pipelines Based on Monte Carlo Method [J]. Oil & Gas Storage and Transportation, 2015, 34 (8): 811-816.

[16] 張 崎,李興斯.結(jié)構(gòu)可靠性分析的模擬重要抽樣方法[J].工程力學(xué),2007,24(1):33-36.

Zhang Qi, Li Xingsi. Importance Sampling Approach in Structural Reliability Analysis Based on Kriging Simulation [J]. Engineering Mechanics, 2007, 24 (1): 33-36.

[17] 燕秀發(fā),謝禹鈞,戴 耀.壓力管道可靠性安全評定[J].機械設(shè)計與造,2004,2(1):6-8.

Yan Xiufa, Xie Yujun, Dai Yao. Pressure Pipings Reliability Safety Assessment [J]. Machinery Design and Manufacture, 2004, 2 (1): 6-8.

[18] 曾初升.316 L不銹鋼腐蝕性能電化學(xué)研究[D].昆明:昆明理工大學(xué),2006:65-85.

Zeng Chusheng. Electrochemical Study on Corrosion Behavior of 316 L Stainless Steel [D]. Kunming: Kunming University of Science and Technology, 2006: 65-85.

[19] 李 陽.316 L不銹鋼焊接接頭的微觀組織及點蝕性能[D].長春:吉林大學(xué),2015:21-28.

Li Yang. The Microstructure and Pitting Corrosion Property of the 316 L Stainless Steel Welded Joint [D]. Changchun: Jilin University, 2015: 21-28.

[20] 張艷麗.X 65碳鋼和316 L不銹鋼在模擬油田采出水中的腐蝕行為研究[D].青島:中國海洋大學(xué),2014.

Zhang Yanli. Corrosion Behavior of X 65 Carbon Steel and 316 L Stainless Steel in Simulated Oilfield Product Water [D]. Qingdao: Ocean University of China, 2014.