張 銳.

(中石化石油工程設(shè)計有限公司，山東東營 257000)

某氣藏已開發(fā)至中期，大都氣、水同產(chǎn)，且氣、液中含H2S、CO2等腐蝕介質(zhì)。它們?nèi)苡谒笮纬扇跛?，對井下管柱具有極強的腐蝕性，不僅影響了油氣井的正常作業(yè)，降低了油管的使用期限，而且給修井作業(yè)帶來不必要的麻煩。大量的現(xiàn)場經(jīng)驗表明，含H2S、CO2共存時的腐蝕行為遠比它們單獨作用時復(fù)雜的多。由于不同區(qū)塊CO2和H2S含量各不相同，井下腐蝕環(huán)境各異，即使在相同的區(qū)塊，井下油管的腐蝕狀況也會存在差別，因此，有效地開展預(yù)測油管腐蝕速度及其影響因素和相互作用關(guān)系的科學(xué)研究，對預(yù)測油管使用壽命，保證氣井的安全生產(chǎn)，降低修井成本，具有重要意義。

氣井井下油管腐蝕是一個復(fù)雜的體系，腐蝕速度也是由各種影響因素相互作用、共同影響的，它不但與腐蝕介質(zhì)濃度及分壓因素有關(guān)，而且與環(huán)境條件、油管材質(zhì)及產(chǎn)水量等有關(guān)。由于各影響因素具有隨機性和不確定性，因而，用簡單的數(shù)理模型來表達各影響因素與腐蝕速度之間對應(yīng)的關(guān)系是不切實際的。像這種具有模糊性和復(fù)雜性的腐蝕速度的預(yù)測問題，可以把井下管柱的腐蝕看成是灰色系統(tǒng)，通過使用滑動平均法對傳統(tǒng)的灰色模型進行優(yōu)化，結(jié)果表明具有很好的適用性。

1 滑動平均的灰色預(yù)測的原理及方法

(1)確定需要預(yù)測的灰色系統(tǒng)中，某項指標在連續(xù)若干個時間內(nèi)按順序排列的原始數(shù)列。X(0)=[X(0)(1),X(0)(2),X(0)(3),…,X(0)(n)]

(2)采用滑動平均法處理原始數(shù)列。

滑動平均法的計算公式為：

(1)

兩端點的計算公式為：

(2)

(3)

(3)累加生成1-AGO序列。

經(jīng)滑動平均法計算后，將所得的新序列作一次累加處理，得到X′(0)的1-AGO序列X(1)。

i=1,2,…,n

X(1)=[x(1)(1),x(1)(2),…,x(1)(n)]

=[x′(0)(1),x′(0)(1)+x′(0)(2),…,x′(0)(1)

+…+x′(0)(n)]

(4)

(4)將累加后的數(shù)列X(1)用微分方程(5)近似表示。

(5)

其中a為發(fā)展灰數(shù)，u為內(nèi)生控制灰數(shù)。

(5)構(gòu)造向量Yn和矩陣B。

Yn=[x(0)(2),x(0)(3),…,x(0)(n)]T

(6)

(7)

(8)

(7)確定模型。

(9)

(8)對計算所得的擬合值進行一次累減處理，得到1-IAGO序列，得到原始數(shù)據(jù)的還原值。

(10)

2 模型精度檢驗

灰色系統(tǒng)理論有三種模型精度檢驗方法，分別為：殘差檢驗、關(guān)聯(lián)度檢驗和后驗差檢驗。

(1)殘差檢驗。

殘差檢驗包括絕對誤差檢驗和相對誤差檢驗。

絕對誤差：

(11)

相對誤差：

(12)

(2)關(guān)聯(lián)度檢驗。

分別計算原始數(shù)據(jù)還原值與其真實值之間絕對誤差的絕對值，比較得出其最小值和最大值：

(13)

(14)

計算關(guān)聯(lián)系數(shù)w(i)：

(15)

k為分辨系數(shù)，一般取0.5。

關(guān)聯(lián)度R計算：

(16)

檢驗標準：k=0.5時，R>0.6有意義。

(3)后驗差檢驗。

原始數(shù)列的均值：

(17)

殘差均值：

(18)

原始數(shù)列的方差：

(19)

殘差的方差：

(20)

后驗差比值C：

(21)

小誤差概率P：

(22)

(4)精度等級。

按照上面兩個指標，可以將預(yù)測精度劃分為四個等級，等級劃分如表1所示:

表1 預(yù)測精度等級劃分Table 1 Classification of prediction accuracy

3 實例應(yīng)用

高溫、氣液混輸和天然氣中普遍含CO2、H2S是某氣田生產(chǎn)腐蝕環(huán)境的主要特點，由此產(chǎn)生的腐蝕問題較為嚴重。該氣田由于井下油管投運時間長，井下流體腐蝕性強，導(dǎo)致井下油套管腐蝕嚴重，影響氣井生產(chǎn)和修井作業(yè),部分油管腐蝕嚴重導(dǎo)致氣井修死。M氣井為中度腐蝕的一口井，采用光油管完井，即油套環(huán)空連通，對套管采用多臂井徑儀檢測發(fā)現(xiàn)在741米井段腐蝕變形。

為了真實可靠地評價井下的腐蝕狀態(tài)，解決井下掛片檢測實施的不便性、局限性和與真實腐蝕狀態(tài)的誤差性，采用了CFL-1動態(tài)腐蝕試驗環(huán)道裝置進行實驗研究，該裝置操作方便，真實可靠，大大降低了腐蝕檢測成本，減少作業(yè)，為安全生產(chǎn)提供指導(dǎo)。根據(jù)M井的溫度、壓力、流速等參數(shù)確定出M井最近七年的腐蝕速度，如表2所示。

(1)原始數(shù)列為:

X(0)=[0.262,0.247,0.217,0.183,0.157,0.142,0.131]利用滑動平均法對原始數(shù)據(jù)進行處理,得到的滑動平均序列為：X′(0)=[0.25825,0.24325,0.216,0.185,0.15975,0.143,0.13375]經(jīng)滑動平均處理進行一次累加后得到的新序列X′(0)：X(1)=[0.25825,0.5015,0.7175,0.9025,1.06225,1.20525,1.339]

表2 某氣田M井最近七年的腐蝕速度實測值Table 2 Measured corrosion rates of M well in a gas field in the last seven years

(2)參數(shù)估計并構(gòu)造預(yù)測模型。

用最小二乘法解得：

則對應(yīng)的時間響應(yīng)函數(shù)，即預(yù)測模型為：

X(1)(k+1)=2.194074-1.935824e-0.135k

(3)累減還原預(yù)測結(jié)果。

(4)計算相對誤差。

q(i)=[-0.0145,-0.0102,-0.0159,0.0193,0.0368,0.0035,-0.0522]K取0.5時，計算得到關(guān)聯(lián)度R=0.79，滿足R>0.6的檢驗標準。

(5)后驗差檢驗。

殘差均值：ξ(0)=0.02177

原始數(shù)列的方差：S12=0.002679

殘差的方差：S22=0.000285

后驗差比值C：C=0.326

計算得出0.6745S1=0.03491，由表3可以看出i=1,2，……，7時，|ξ(0)(i)-`ξ(0)|的值均小于0.6745S1，所以P值為1。

(6)模型精度檢驗。

以某氣田M井最近七年測得的腐蝕速度為研究對象，對七年的腐蝕速度進行擬合，擬合情況如表4所示：

表3 |ξ(0)(i)-ξ(0)|值的求解Table 3 The solution of |ξ(0)(i)-ξ(0)|

表3 |ξ(0)(i)-ξ(0)|值的求解Table 3 The solution of |ξ(0)(i)-ξ(0)|

i|ξ(0)(i)-ξ(0)|10.072 720.011 5730.005 8740.002 4750.015 0360.018 2770.030 43

表4 擬合值和實際測量值、殘差和相對殘差Table 4 Fitting values , actual measured values, residuals and relative residuals

根據(jù)預(yù)測精度等級劃分結(jié)果表3可知，此滑動平均的灰色預(yù)測模型的預(yù)測精度等級為好。且由表4可見，實際測量值和擬合值很接近，殘差和相對殘差都較小，相對誤差控制在5%以內(nèi)，具有較高的擬合精度。因此，可以用該模型對該氣田M井未來五年的腐蝕速度進行預(yù)測。預(yù)測結(jié)果如表5所示。

4 結(jié)論

(1)針對井下油管腐蝕問題，提出了滑動平均的灰色理論預(yù)測模型，為井下腐蝕監(jiān)測提供依據(jù)，是處理油管腐蝕速度預(yù)測問題的有效工具。

(2)根據(jù)M井最近七年的油管腐蝕速度實測值，建立油管腐蝕預(yù)測模型，經(jīng)驗證，該模型的預(yù)測精度等級為好(小誤差概率P=1>0.95，后驗差比值C=0.325<0.35)，預(yù)測精度控制在5%以內(nèi)。

表5 氣田M井未來五年腐蝕速度的預(yù)測值Table 5 Predicted corrosion rate of M well in gas field in the next five years

(3)建立的油管腐蝕速度預(yù)測模型隨時間呈指數(shù)函數(shù)變化，并預(yù)測得出M井未來五年的油管腐蝕速度分別為0.108 8 mm/s，0.094 5 mm/s，0.083 mm/s，0.072 5 mm/s，0.063 4 mm/s。

(4)推廣該模型用于油管腐蝕預(yù)測，可根據(jù)反饋的數(shù)據(jù)制定更為有效的防腐措施，或為油管的更換使用提供指導(dǎo)作用。