王海濤,孔明慧

(大連交通大學(xué) 土木與安全工程學(xué)院，大連 116028)

?

三次指數(shù)平滑法預(yù)測管道腐蝕速率的應(yīng)用

王海濤,孔明慧

(大連交通大學(xué) 土木與安全工程學(xué)院，大連 116028)

為準(zhǔn)確預(yù)測管道腐蝕速率并得到精度較高的腐蝕速率預(yù)測模型，采用已有統(tǒng)計數(shù)據(jù)作為研究對象。首先，采用三次指數(shù)平滑法，建立管道腐蝕速率預(yù)測模型，對腐蝕速率數(shù)據(jù)進(jìn)行擬合與預(yù)測分析；然后，得出預(yù)測模型中最合理的權(quán)重系數(shù)α；最后，通過與一次指數(shù)平滑法、二次指數(shù)平滑法進(jìn)行對比分析，得出了三次指數(shù)平滑法預(yù)測精度較高，預(yù)測值與實際值相符合的結(jié)論。結(jié)果表明，三次指數(shù)平滑法預(yù)測模型符合管道腐蝕速率的特點，對有效預(yù)測管道腐蝕具有一定的參考價值和指導(dǎo)意義。

管道腐蝕；腐蝕速率預(yù)測；預(yù)測模型；三次指數(shù)平滑法；權(quán)重系數(shù)；對比分析

腐蝕預(yù)測一直是腐蝕領(lǐng)域中重要的研究內(nèi)容。國內(nèi)外學(xué)者在腐蝕機理、腐蝕剩余強度評價等方面開展了大量的研究工作，并取得了一定的研究成果，但對于腐蝕分析和評價中有著重要作用的預(yù)測模型的研究卻相對較少。管道腐蝕速率的預(yù)測模型也僅限于de Waard預(yù)測模型[1]、B.Mishra預(yù)測模型[2]、BP神經(jīng)網(wǎng)絡(luò)模型[3]、灰色模型[4-6]等幾個。大部分方法普遍存在“長周期、大區(qū)域、低信度”的缺陷。如神經(jīng)網(wǎng)絡(luò)需要較為全面、能代表各種情況的樣本數(shù)據(jù)，否則訓(xùn)練出的神經(jīng)網(wǎng)絡(luò)外延性不強。而采用灰色模型的預(yù)測精度不理想，原因是模型要求數(shù)據(jù)序列必須呈指數(shù)規(guī)律，而且對隨機波動性較大的數(shù)據(jù)序列的擬合較差。

針對以上預(yù)測方法的不足和缺陷，本工作提出了一種基于三次指數(shù)平滑法的預(yù)測模型，該模型具有使用方便、操作簡單，采用“重近輕遠(yuǎn)”的策略，使其跟蹤時序的變化，能夠較好地顯現(xiàn)時序的變化趨勢，適用于數(shù)據(jù)波動較大，而且數(shù)據(jù)呈現(xiàn)非線性變化規(guī)律的情況。

1　基于三次指數(shù)平滑法的預(yù)測模型

指數(shù)平滑法是布朗(Robert Brown)所提出的一種對時間序列的態(tài)勢具有穩(wěn)定性或規(guī)則性的預(yù)測方法。他認(rèn)為時間序列是可以被合理順勢推延的，最近的過去態(tài)勢在某種程度上會持續(xù)到最近的未來，所以將較大的權(quán)數(shù)放在最近的資料中[7]。

指數(shù)平滑法具有使用方便、操作簡單等特點，近年來被廣泛應(yīng)用于諸多研究領(lǐng)域[8-11]。指數(shù)平滑法是在移動平均法基礎(chǔ)上發(fā)展起來的一種時間序列分析預(yù)測法，通過計算指數(shù)平滑值，配合一定的時間序列預(yù)測模型對現(xiàn)象的未來進(jìn)行預(yù)測。指數(shù)平滑法包括一次指數(shù)，二次指數(shù)，三次指數(shù)。簡單一次指數(shù)平滑法只適用于數(shù)據(jù)波動不大的模型，對有些在某一時間突然有較大變化的數(shù)據(jù),預(yù)測時會產(chǎn)生較大的誤差。而三次指數(shù)能夠跟蹤時序的非線性變化趨勢，而且隨著計算機技術(shù)的發(fā)展，該方法可以在Excel表格中實現(xiàn)，可操作性和實用性很強，被廣泛應(yīng)用于各個領(lǐng)域[12]。

1.1三次指數(shù)平滑法的預(yù)測模型

三次指數(shù)平滑法是在二次指數(shù)平滑值的基礎(chǔ)上又進(jìn)行了一次平滑，可以用它來估算二次多項式的參數(shù)值，三次指數(shù)平滑法的預(yù)測模型如下：

(1)

式中：Yt+T為預(yù)測值；T為預(yù)測的超前時間數(shù)，t=1,2,3,…n,(n為原始數(shù)據(jù)個數(shù)),at，bt，ct為預(yù)測系數(shù)。at，bt，ct的表達(dá)式如下：

(2)

(3)

(4)

(5)

(6)

(7)

式中：Xt為實際值，α為權(quán)重系數(shù)，在0～1.0取值。

三次指數(shù)平滑法與一次、二次指數(shù)平滑法相比要復(fù)雜很多，但修正預(yù)測值的目的是相同的，三次指數(shù)平滑能夠跟蹤時序的非線性變化趨勢[13]。

1.2初始值和權(quán)重系數(shù)的選取

當(dāng)原數(shù)列的項數(shù)較多時(大于15項)，可以選用第一期的觀察值或選用比第一期前一期的觀察值作為初始值。當(dāng)原數(shù)列的項數(shù)較少時(小于15項)，可以選取最初幾期(一般為前三期)的平均數(shù)作為初始值。

權(quán)重系數(shù)α是新舊數(shù)據(jù)在預(yù)測中所起不同作用的比例因子[14]。α越大，新數(shù)據(jù)所起的作用越大，模型的靈敏度就越高，適應(yīng)新水平較快，但容易過敏；反之，α過小，則比較保守，可能會落后新的發(fā)展趨勢。一般來說，α值應(yīng)按以下情況處理：

(1) 如果觀察值的長期趨勢變動接近穩(wěn)定的常數(shù)，應(yīng)取居中α值(一般取0.4～0.6)，使觀察值在指數(shù)平滑中具有大小接近的權(quán)數(shù)；

(2) 如果觀察值呈現(xiàn)明顯的變動時，則宜取較大的α值(一般取0.6～0.9)，從而使近期觀察值能迅速反映在未來的預(yù)測值中；

(3) 如果觀察值的長期趨勢變動較緩慢，則宜取較小的α值(一般取0.1～0.4)，使遠(yuǎn)期觀察值的特征也能反映在指數(shù)平滑值中。

2　實例分析

2.1數(shù)據(jù)來源

本試驗?zāi)M了在常溫常壓下管線鋼的腐蝕速率隨時間的變化規(guī)律[15]。利用失重法對腐蝕試件進(jìn)行分析，通過計算得出其腐蝕速率，共8個數(shù)據(jù)，作為原始樣本(見表1)。根據(jù)試驗得出的腐蝕速率數(shù)據(jù)，對其進(jìn)行擬合與預(yù)測分析。

表1　失重量與平均腐蝕速率隨時間的變化

2.2基于三次指數(shù)平滑法模型在管道腐蝕速率預(yù)測中的應(yīng)用

首先把表1中的數(shù)據(jù)轉(zhuǎn)換成曲線圖，見圖1?？梢钥闯?數(shù)據(jù)具有比較明顯的波動性，而且具有明顯的非線性變化趨勢。應(yīng)用三次指數(shù)平滑預(yù)測方法，對實際值與預(yù)測值進(jìn)行擬合。具體步驟如下:

(1) 初始值的選取

(2) 權(quán)重系數(shù)的選取

如果觀察值的長期趨勢變動較緩慢，則宜取較小的α值(一般取0.1～0.4)。

從圖1中可以看出，預(yù)測所需的原始數(shù)據(jù)即管線腐蝕試驗得出的腐蝕速率值漸漸趨于緩慢變化，故權(quán)重系數(shù)α取較小值，在0.1～0.4中選取。取α為0.2,0.3,0.32,0.33,0.35,0.4分別對腐蝕速率進(jìn)行預(yù)測，如圖2所示，最后選用預(yù)測值與實測值之間最小的誤差時的α作為權(quán)重系數(shù)。

在預(yù)測應(yīng)用時，α需根據(jù)觀察值的具體變化情況，來大致確定α取值范圍，然后取幾個α值進(jìn)行試算，然后用管線鋼腐蝕試驗所得數(shù)據(jù)的預(yù)測值與實際值進(jìn)行誤差分析，經(jīng)過誤差分析以及圖2所示，最終取α=0.33進(jìn)行預(yù)測,見圖2(d)。

(3) 預(yù)測結(jié)果

將初始值分別代入式(5)～(7)，求出所對應(yīng)得一次、二次、三次指數(shù)平滑修勻新數(shù)列值，然后再把新值代入式(2)～(4)，求出預(yù)測系數(shù)，最后代入三次指數(shù)平滑預(yù)測式(1)，以此類推，求得最終的預(yù)測值，見表2。

經(jīng)計算:

S15(1)=0.371 734 08，S15(2)=0.437 980 946，S15(3)=0.459 842 412;a15=0.261 101 813，b15=-0.081 736，c15=-0.005 384。

所以，Y15+T=0.261 101 813-0.081 736T-0.005 384T2，分別將T=1，T=2代入可得第16天、17天的管線腐蝕速率的預(yù)測值。

3　實例驗證

參考相關(guān)管道腐蝕統(tǒng)計資料，以我國四川某管道不同時間測得的平均腐蝕速率(表3)為初始預(yù)測數(shù)據(jù)建立預(yù)測模型。

根據(jù)上述計算過程，應(yīng)用布朗三次指數(shù)平滑法對該組數(shù)據(jù)進(jìn)行分析和預(yù)測，1999-2007年的預(yù)測值分別由該年份之前數(shù)據(jù)取T=1計算。預(yù)測結(jié)果趨勢圖如圖3所示。由圖3可見，針對管道腐蝕進(jìn)行的布朗三次指數(shù)平滑法預(yù)測精度較高，當(dāng)α=0.33時預(yù)測值與實際值基本符合，預(yù)測值曲線與實際值曲線擬合成功，且2005、2006 和2007 年的預(yù)測值與實際值的相對誤差分別為1.08%、1.24% 和1.86%，表明應(yīng)用這種方法進(jìn)行管道腐蝕速率預(yù)測是合理可行的?？梢酝ㄟ^這種方法進(jìn)行預(yù)測，預(yù)報未來短期內(nèi)的管道發(fā)展趨勢，為決策者提供決策依據(jù)。

表2　三次指數(shù)平滑法預(yù)測計算表

表3　不同時間測得的16Mn管道鋼的平均腐蝕速率

4　結(jié)論

(1) 提出了將三次指數(shù)平滑法預(yù)測模型應(yīng)用于管道腐蝕速率預(yù)測的方法。結(jié)合我國埋地管線腐蝕的實際情況和事故特性，探討了腐蝕速率預(yù)測的理論。

(2) 基于匯總統(tǒng)計的管道腐蝕速率數(shù)據(jù)，應(yīng)用布朗三次指數(shù)平滑法,對1999-2007年之前數(shù)據(jù)取T=1計算，預(yù)測、繪制出了管道腐蝕速率的發(fā)展趨勢，預(yù)測值與實際值基本符合,誤差較小。

(3) 布朗三次指數(shù)平滑法預(yù)測精度較高，而且用Excel很容易實現(xiàn)，易于理解，便于接受。

[1]WAARD C D,LOTZ U,MILLIAMS D E. Predictive Model for CO2corrosion engineering in wet natural gas pipelines[J]. Corrosion,1991,47(12):976-985.

[2]MISHRA B,HASSAN S A,OLSON D L,et al. Development of a predictive model for activation-controlled corrosion of steel in solutions containing carbon dioxide[J]. Corrosion,1997,53(11):852-859.

[3]章玉婷,楊劍鋒. 基于BP神經(jīng)網(wǎng)絡(luò)的管道腐蝕速率預(yù)測[J]. 腐蝕研究,2013,27(9):67-71.

[4]喻西崇,趙金洲,鄔亞玲,等. 利用灰色理論預(yù)測注水管道腐蝕速率的變化趨勢[J]. 腐蝕與防護,2003,24(2):51-54.

[5]李秀娟,梅甫定. 灰色神經(jīng)網(wǎng)絡(luò)模型在油氣管道腐蝕速度預(yù)測中的應(yīng)用[J]. 安全與環(huán)境工程,2006,13(4):77-80.

[6]鄧聚龍. 灰色系統(tǒng)理論及其應(yīng)用[M]. 武漢:華中理工大學(xué)出版社,1988.

[7]齊穎. 一次指數(shù)平滑模型預(yù)測法及實際應(yīng)用[J]. 北京機械工業(yè)學(xué)院學(xué)報,2004,19(4):50-53,59.

[8]顏柳,麻鳳海. 三次指數(shù)平滑法在城市地鐵變形預(yù)測中的應(yīng)用[J]. 交通科技與經(jīng)濟,2007(5):62-71.

[9]俞志峰,謝正文. 基于指數(shù)平滑技術(shù)的灰色沉降預(yù)測模型及應(yīng)用[J]. 中南公路工程,2007,32(3):120-123.

[10]顧曉葉,于殿寶,朱玉生. 事故預(yù)測數(shù)學(xué)模型的研究與實踐[J]. 中國安全科學(xué)學(xué)報,2009,19(12):23-27.

[11]袁昌明,崔曉君. 事故預(yù)測模型的建立與應(yīng)用[J]. 中國計量學(xué)院學(xué)報,2004,15(4):314-316.

[12]王洪德,曹英浩. 道路交通事故的三次指數(shù)平滑預(yù)測法[J]. 遼寧工程技術(shù)大學(xué)學(xué)報,2014,33(1):42-46.

[13]朱慶明,張浩. 三次指數(shù)平滑法在煤礦事故預(yù)測中的應(yīng)用研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù),2012,8(4):103-106.

[14]王長江. 指數(shù)平滑法中平滑系數(shù)的選擇研究[J]. 中北大學(xué)學(xué)報,2006,27(6):558-561.

[15]李彬. 海底集輸管線鋼的腐蝕行為研究[D]. 西安:西安工業(yè)大學(xué)，2011.

Application of Cubic Exponential Smoothing Method to Pipeline Corrosion Rate Forecasting

WANG Hai-tao, KONG Ming-hui

(School of Civil and Safety Engineering, Dalian Jiaotong University, Dalian 116028, China)

In order to predict the pipeline corrosion rate accurately and get high precision of corrosion rate forecasting model, the existing experimental data were considered as the research object. Firstly, the cubic exponential smoothing method was used to establish pipeline corrosion rate forecasting model to fit and forecast the data of the corrosion rate. Secondly, the most reasonable weight coefficient α in the prediction model was obtained. Finally, a comparison with the single exponential smoothing method and the double exponential smoothing indicates that the forecasting precision of the cubic exponential smoothing method was preferable, and the predicted values conformed to the true values. Results show that the cubic exponential smoothing forecasting model is consistent with the characteristics of pipeline corrosion rate, and the method can provide a certain reference and guidance for effective forecasting the pipeline corrosion.

pipeline corrosion; corrosion rate forecasting; forecasting model; cubic exponential smoothing method; weight coefficient; contrastive analysis

10.11973/fsyfh-201601002

2014-12-23

國家自然科學(xué)基金項目(51208073)； 遼寧省博士啟動基金項目(20121061)

王海濤(1982-),副教授,博士,從事巖土及地下工程研究工作,15998581208,wht@djtu.edu.cn

TE988

B

1005-748X(2016)01-0008-04