徐昊 石慧 王遠(yuǎn)洋

摘要：油氣管道在復(fù)雜惡劣的運(yùn)行條件下會(huì)受到各種腐蝕因素與管道自身退化的影響，大幅降低管道運(yùn)行的安全性和可靠性。為提高管道壽命預(yù)測(cè)的準(zhǔn)確度，建立了一種融合腐蝕速率模型和數(shù)據(jù)驅(qū)動(dòng)的壽命預(yù)測(cè)方法。首先根據(jù)管道的歷史退化數(shù)據(jù)建立維納退化模型，其次利用腐蝕速率模型對(duì)漂移系數(shù)進(jìn)行描述，從而估計(jì)模型中的未知參數(shù)，最終獲得管道剩余壽命的概率密度函數(shù)及預(yù)測(cè)值。以某腐蝕油氣管道的退化數(shù)據(jù)為例，驗(yàn)證該模型的準(zhǔn)確性。研究結(jié)果表明：結(jié)合維納過程和數(shù)據(jù)驅(qū)動(dòng)模型，將影響管道腐蝕的外部環(huán)境因素和管道自身退化數(shù)據(jù)充分利用，能夠有效預(yù)測(cè)管道的壽命，具有較好的可行性與適用性。研究結(jié)果可為油氣管道維修策略制定提供參考。

關(guān)鍵詞：油氣管道；剩余壽命；維納過程；腐蝕速率；數(shù)據(jù)驅(qū)動(dòng)

0 引 言

管道是油氣的主要運(yùn)輸方式，在國(guó)民經(jīng)濟(jì)中也占相當(dāng)大的比重，管道的安全運(yùn)行也與人們的生活息息相關(guān)［1-2］。我國(guó)地域遼闊，地質(zhì)復(fù)雜，腐蝕環(huán)境差異顯著，管道極易受到外部環(huán)境的腐蝕，這樣就大幅降低了其安全性，縮短了使用壽命［3-4］。為有效預(yù)防油氣管道發(fā)生災(zāi)難性的事故，對(duì)管道進(jìn)行剩余壽命預(yù)測(cè)有著重要的意義。

長(zhǎng)輸管道腐蝕速率預(yù)測(cè)模型主要有經(jīng)驗(yàn)型預(yù)測(cè)模型、半經(jīng)驗(yàn)型預(yù)測(cè)模型和機(jī)理型預(yù)測(cè)模型［5-6］。曹乃寧［7］構(gòu)建的腐蝕速率模型首次考慮了CO2分壓的影響因素。劉熠［8］通過對(duì)油氣管道腐蝕因素的分類總結(jié)，建立了新的氣油兩相流腐蝕速率模型。近些年，對(duì)腐蝕速率的研究主要是建立新的腐蝕速率模型或者通過數(shù)據(jù)統(tǒng)計(jì)的方法建立數(shù)學(xué)模型，從而對(duì)管道進(jìn)行腐蝕速率預(yù)測(cè)［7］。腐蝕速率模型通?；诓煌母g機(jī)理和腐蝕因素建立，因?yàn)楣艿赖倪\(yùn)行狀態(tài)以及腐蝕因素復(fù)雜多樣，所以通過腐蝕速率模型預(yù)測(cè)腐蝕速率需要選擇相對(duì)應(yīng)的管道進(jìn)行管道的壽命預(yù)測(cè)，該方法具有一定的局限性。通過腐蝕速率模型計(jì)算出腐蝕速率數(shù)值，從而對(duì)管道剩余壽命進(jìn)行預(yù)測(cè)，雖然可以較快地預(yù)測(cè)管道剩余使用壽命，進(jìn)而采取相關(guān)的防護(hù)措施，但這種方法僅僅考慮了管道腐蝕的環(huán)境因素，并沒有考慮管道本身的退化量。

目前，國(guó)內(nèi)外學(xué)者通過數(shù)據(jù)驅(qū)動(dòng)的方法對(duì)管道進(jìn)行剩余壽命的預(yù)測(cè)已有大量研究，M.GASPERIN等［9］對(duì)設(shè)備進(jìn)行壽命分布的研究時(shí)充分考慮了線性退化模型；ZHOU M.等［10］總結(jié)了不同的模型對(duì)管道可靠性的影響，通過隨機(jī)過程對(duì)管道腐蝕的增長(zhǎng)進(jìn)行描述；張新生等［11］運(yùn)用逆高斯過程建立腐蝕退化模型；SUN B.等［12］結(jié)合了隨機(jī)過程和干擾理論，進(jìn)而對(duì)管道的腐蝕退化過程進(jìn)行建模；陳亮等［13］對(duì)管道進(jìn)行剩余壽命預(yù)測(cè)時(shí)充分考慮了測(cè)量誤差的影響；張鵬等［14］在考慮腐蝕油氣管道腐蝕穿孔、局部爆破和整體斷裂的條件下，建立了隨機(jī)變量相關(guān)的失效概率模型；FU G.Y.等［15］考慮多種失效模式及其隨機(jī)相關(guān)性，分析了埋地管道的可靠性；V.ARYAI等［16］在基于融合伽馬過程和Copula函數(shù)的條件下，估算了管道的失效概率。但是以上方法主要采用管道自身的歷史退化數(shù)據(jù)來(lái)進(jìn)行分析，缺乏一定的經(jīng)驗(yàn)指導(dǎo)，無(wú)法充分有效利用管道腐蝕模型指導(dǎo)壽命預(yù)測(cè)，不能對(duì)管道腐蝕的退化機(jī)理和規(guī)律進(jìn)行深入分析。

筆者將腐蝕速率物理模型與數(shù)據(jù)驅(qū)動(dòng)有機(jī)融合，應(yīng)用維納過程建立油氣管道退化模型，通過腐蝕速率模型估計(jì)模型中的未知參數(shù)，不僅可以充分考慮到腐蝕管道自身退化過程的不確定性，而且腐蝕速率模型可以考慮到溫度、CO2、H2S等外部環(huán)境影響因素，從而更加準(zhǔn)確地預(yù)測(cè)管道壽命。

1 模型描述

1.1 維納過程

維納過程可以描述設(shè)備的退化過程，具有良好的計(jì)算和分析能力［17-18］。一方面，如果設(shè)備退化過程隨時(shí)間線性變化，則可以描述線性漂移維納過程。另一方面，線性漂移維納過程模型的剩余壽命分布服從逆高斯過程。本節(jié)使用線性漂移維納過程模型來(lái)對(duì)油氣管道進(jìn)行失效模型的建立。

2 退化模型構(gòu)建

2.1 模型構(gòu)建

綜合上述對(duì)維納過程和腐蝕速率修正模型的分析，可通過維納過程對(duì)油氣管道退化數(shù)據(jù)進(jìn)行建模，并通過腐蝕速率模型描述模型中的漂移系數(shù)。其基本流程如圖2所示。

預(yù)測(cè)方法具體步驟如下：

（1）通過對(duì)檢測(cè)到的管道退化數(shù)據(jù)進(jìn)行正態(tài)分布檢驗(yàn)，判斷管道的腐蝕深度退化過程是否符合維納過程。若服從，則可對(duì)管道的腐蝕深度建立基于維納過程的性能退化模型。

（2）利用腐蝕速率修正模型得到不同檢測(cè)時(shí)刻的腐蝕速率，從而估計(jì)模型中的漂移系數(shù)θ；并通過極大似然估計(jì)進(jìn)一步獲取擴(kuò)散系數(shù)σ的估計(jì)值。

（3）將步驟（2）所得到的參數(shù)代入剩余壽命概率密度函數(shù)式中。根據(jù)概率密度函數(shù)定義，其最大值處所對(duì)應(yīng)的值即為管道所對(duì)應(yīng)的剩余壽命。

2.2 剩余壽命預(yù)測(cè)建模

3.3 壽命預(yù)測(cè)與腐蝕速率模型對(duì)比分析

通過腐蝕速率模型與極大似然估計(jì)法得到漂移系數(shù)與擴(kuò)散系數(shù)，并將不同檢測(cè)時(shí)刻的參數(shù)估計(jì)值代入該模型，結(jié)合式（6）和式（10）可得到管道剩余壽命的概率分布。概率密度函數(shù)的峰值作為不同檢測(cè)時(shí)刻的預(yù)測(cè)結(jié)果。圖5中，管道剩余壽命的實(shí)際值由紅色虛線表示，所求管道剩余壽命預(yù)測(cè)值由紅色實(shí)線表示，預(yù)測(cè)結(jié)果與實(shí)際值之間的誤差由兩條線之間的距離表示。

由圖5可知，隨著管道運(yùn)行年限的不斷延長(zhǎng)，管道剩余壽命概率密度函數(shù)的峰值也逐漸遞增，說明管道的剩余壽命不斷縮短，符合管道運(yùn)行的狀況。

將不同時(shí)刻下所檢測(cè)到的概率密度函數(shù)最大值所對(duì)應(yīng)的時(shí)間作為管道的預(yù)測(cè)值，可得到管道剩余壽命預(yù)測(cè)的結(jié)果，如圖6所示。

由圖6可知，隨著檢測(cè)時(shí)刻的不斷延長(zhǎng)，預(yù)測(cè)值與實(shí)際值之間的誤差逐漸縮小，預(yù)測(cè)精度也逐漸提高。

為驗(yàn)證本文方法的預(yù)測(cè)效果，與腐蝕速率模型進(jìn)行了對(duì)比分析?；诟g速率模型預(yù)測(cè)管道壽命的方法需要管道在不同的腐蝕環(huán)境下進(jìn)行數(shù)據(jù)采集，從而建立相對(duì)應(yīng)的數(shù)學(xué)表達(dá)式來(lái)預(yù)測(cè)管道壽命。該方法缺乏對(duì)管道自身退化過程的充分考慮，會(huì)導(dǎo)致模型預(yù)測(cè)誤差偏大?；诟g速率預(yù)測(cè)管道剩余壽命的公式：

通過本文提出的腐蝕速率模型對(duì)目標(biāo)管道進(jìn)行剩余壽命預(yù)測(cè)，并與實(shí)際值進(jìn)行對(duì)比，結(jié)果如圖7所示。

由圖7可知，本文提出的模型在不同檢測(cè)時(shí)刻的相對(duì)誤差均小于腐蝕速率模型，證明了前者在充分考慮管道腐蝕外部環(huán)境的影響下，提高了維納退化模型的預(yù)測(cè)精度，但該模型仍存在一定的誤差，所以在模型預(yù)測(cè)準(zhǔn)確性方面仍然需要進(jìn)一步地開展研究。

4 結(jié) 論

本文首先根據(jù)管道自身退化數(shù)據(jù)構(gòu)造退化模型，在考慮管道腐蝕外部環(huán)境因素的情況下，引入腐蝕速率模型估計(jì)維納過程模型中的未知參數(shù)，進(jìn)而進(jìn)行管道的壽命預(yù)測(cè)。最后以某油氣管道為例，將本文方法的預(yù)測(cè)結(jié)果與腐蝕速率模型進(jìn)行對(duì)比，驗(yàn)證了本文提出方法的準(zhǔn)確性。得到如下結(jié)論：

（1）通過分析油氣管道的退化機(jī)理和腐蝕影響因素，根據(jù)管道腐蝕退化特點(diǎn)建立了基于維納過程的腐蝕管道退化模型。

（2）結(jié)合維納過程和腐蝕速率模型，將影響管道腐蝕的外部環(huán)境因素建立起腐蝕速率模型和管道退化數(shù)據(jù)充分利用，有效提高了模型預(yù)測(cè)精度。

（3）與單一的腐蝕速率模型預(yù)測(cè)管道壽命結(jié)果對(duì)比，預(yù)測(cè)值在不同的檢測(cè)時(shí)刻的誤差值均較后者低，證明本文方法具有一定的準(zhǔn)確性，可為油氣管道維修策略制定提供參考。

參考文獻(xiàn)：

［1］辛佳興，陳金忠，李曉龍，等．油氣管道內(nèi)檢測(cè)技術(shù)研究前沿進(jìn)展［J］．石油機(jī)械，2022，50（5）：119-126．

XIN J X，CHEN J Z，LI X L，et al.Advances in oil and gas pipeline internal inspection technology［J］.China Petroleum Machinery，2022，50（5）： 119-126.

［2］李琳，王玉婷．油氣管道土壤腐蝕分級(jí)評(píng)價(jià)研究［J］．安全與環(huán)境學(xué)報(bào)，2019，19（5）：1497-1503．

LI L，WANG Y T.Probe into the grade-degree evaluation of the soil corrosion in the petro-gas transmission pipeline system［J］.Journal of Safety and Environment，2019，19（5）： 1497-1503.

［3］王小完，顧建榮，駱正山．基于可靠性的油氣管道系統(tǒng)維修策略研究［J］．消防科學(xué)與技術(shù)，2015，34（5）：655-659．

WANG X W，GU J R，LUO Z S.A research on maintenance strategy of oil and gas pipeline system based on reliability［J］.Fire Science and Technology，2015，34（5）： 655-659.

［4］趙一昭，柳勇，白凡，等．管道含外表面裂紋時(shí)的疲勞壽命預(yù)測(cè)研究［J］．石油機(jī)械，2020，48（8）：129-135．

ZHAO Y Z，LIU Y，BAI F，et al.Fatigue life prediction of pipe with cracks on external surface［J］.China Petroleum Machinery，2020，48（8）： 129-135.

［5］冷亞梅．含CO2多相流內(nèi)腐蝕研究［D］．東營(yíng)：中國(guó)石油大學(xué)（華東），2013．

LENG Y M.Study on CO2 internal corrosion of multiphase flow［D］.Dongying： China University of Petroleum （East China），2013.

［6］朱景龍，孫成，王佳，等．CO2腐蝕及控制研究進(jìn)展［J］．腐蝕科學(xué)與防護(hù)技術(shù)，2007，19（5）：350-353．

ZHU J L，SUN C，WANG J，et al.Advances on research of CO2 induced corrosion and its protection［J］.Corrosion Science and Protection Technology，2007，19（5）： 350-353.

［7］曹乃寧．油氣輸送腐蝕管道剩余壽命預(yù)測(cè)研究［D］．西安：西安建筑科技大學(xué)，2016．

CAO N N.Residual life prediction technology for oil & gas corroded pipeline［D］.Xian： Xian University of Architecture and Technology，2016.

［8］劉熠．油氣兩相流管道內(nèi)腐蝕研究［D］．成都：西南石油大學(xué)，2007．

LIU Y.Study on internal corrosion of oil gas two phase flow pipeline［D］.Chengdu： Southwest Petroleum University，2007.

［9］GAPERIN M，JURII，BOKOSKI P，et al.Model-based prognostics of gear health using stochastic dynamical models［J］.Mechanical Systems and Signal Processing，2011，25（2）： 537-548.

［10］ZHOU W，XIANG W，HONG H P.Sensitivity of system reliability of corroding pipelines to modeling of stochastic growth of corrosion defects［J］.Reliability Engineering & System Safety，2017，167： 428-438.

［11］張新生，趙夢(mèng)旭，王小完．尾段殘差修正GM（1，1）模型在管道腐蝕預(yù)測(cè)中的應(yīng)用［J］．中國(guó)安全科學(xué)學(xué)報(bào)，2017，27（1）：65-70．

ZHANG X S，ZHAO M X，WANG X W.Application of modified empennage residual error GM （1，1） model in prediction of pipeline corrosion［J］.China Safety Science Journal，2017，27（1）： 65-70.

［12］SUN B，LIAO B P，LI M M，et al.A time-variant reliability model for copper bending pipe under seawater-active corrosion based on the stochastic degradation process［J］.Materials，2018，11（4）： 507.

［13］陳亮，胡昌華．Gamma過程退化模型估計(jì)中測(cè)量誤差影響的仿真研究［J］．系統(tǒng)仿真技術(shù)，2010，6（1）：1-5．

CHEN L，HU C H.Simulation study on the effect of measurement error on estimation of Gamma process degradation model［J］.System Simulation Technology，2010，6（1）： 1-5.

［14］張鵬，彭?xiàng)睿紤]隨機(jī)變量相關(guān)性的腐蝕管道失效概率［J］．石油學(xué)報(bào)，2016，37（10）：1293-1301．

ZHANG P，PENG Y.Failure probability of corroded pipeline considering random variables correlation［J］.Acta Petrolei Sinica，2016，37（10）： 1293-1301.

［15］FU G Y，YANG W，LI C Q，et al.Reliability analysis of corrosion affected underground steel pipes considering multiple failure modes and their stochastic correlations［J］.Tunnelling and Underground Space Technology，2019，87： 56-63.

［16］ARYAI V，BAJI H，MAHMOODIAN M，et al.Time-dependent finite element reliability assessment of cast-iron water pipes subjected to spatio-temporal correlated corrosion process［J］.Reliability Engineering & System Safety，2020，197： 106802.

［17］WHITMORE G A.Estimating degradation by a wiener diffusion process subject to measurement error［J］.Lifetime Data Analysis，1995，1（3）： 307-319.

［18］WANG X.Wiener processes with random effects for degradation data［J］.Journal of Multivariate Analysis，2010，101（2）： 340-351.

［19］ROTA G C.The theory of stochastic processes［J］.Advances in Mathematics，1976，22（1）： 129.

［20］張啟龍，張曉誠(chéng)，韓耀圖，等．渤海多因素條件下的管材腐蝕規(guī)律試驗(yàn)研究［J］．石油機(jī)械，2021，49（7）：130-137．

ZHANG Q L，ZHANG X C，HAN Y T，et al.Experimental research on tubular goods corrosion under multi-factor conditions in Bohai oilfield［J］.China Petroleum Machinery，2021，49（7）： 130-137.

［21］程雅雯，趙會(huì)軍，彭浩平，等．H2S/CO2共存環(huán)境下管道內(nèi)腐蝕研究進(jìn)展［J］．常州大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，27（4）：63-68．

CHENG Y W，ZHAO H J，PENG H P，et al.Research progress of the pipeline under the coexistence of H2S/CO2［J］.Journal of Changzhou University （Natural Science Edition），2015，27（4）： 63-68.

［22］高國(guó)娟．基于陜北油田管道H2S/CO2腐蝕預(yù)測(cè)模型的應(yīng)用研究［D］．西安：西安石油大學(xué)，2020．

GAO G J.Application research of H2S/CO2 corrosion prediction model based on pipeline in northern Shan Xi oilfield［D］.Xian： Xian Shiyou University，2020.

［23］EBENUWA A U，TEE K F.Reliability estimation of buried steel pipes subjected to seismic effect［J］.Transportation Geotechnics，2019，20： 100242.

［24］胡昌華，施權(quán)，司小勝，等．?dāng)?shù)據(jù)驅(qū)動(dòng)的壽命預(yù)測(cè)和健康管理技術(shù)研究進(jìn)展［J］．信息與控制，2017，46（1）：72-82．

HU C H，SHI Q，SI X S，et al.Data-driven life prediction and health management： state of the art［J］.Information and Control，2017，46（1）： 72-82.

［25］任子強(qiáng)，司小勝，胡昌華，等．融合多源數(shù)據(jù)的非線性退化建模與剩余壽命預(yù)測(cè)［J］．中國(guó)測(cè)試，2020，46（2）：1-8．

REN Z Q，SI X S，HU C H，et al.Multi-source data fusion for nonlinear degradation modeling and remaining useful life prediction［J］.China Measurement & Test，2020，46（2）： 1-8.

［26］張新生，李亞云，王小完．基于逆高斯過程的腐蝕油氣管道維修策略［J］．石油學(xué)報(bào)，2017，38（3）：356-362．

ZHANG X S，LI Y Y，WANG X W.Maintenance strategy of corroded oil-gas pipeline based on inverse Gaussian process［J］.Acta Petrolei Sinica，2017，38（3）： 356-362.