鄭運(yùn)虎， 姜　峰， 杜超飛

(蘭州理工大學(xué) 石油化工學(xué)院 ， 甘肅 蘭州 730050)

?

基于灰色關(guān)聯(lián)理論的海洋立管分段風(fēng)險(xiǎn)評(píng)估

鄭運(yùn)虎， 姜峰， 杜超飛

(蘭州理工大學(xué) 石油化工學(xué)院 ， 甘肅 蘭州 730050)

針對(duì)整根立管評(píng)定的可行性差的問(wèn)題，提出了將整根立管進(jìn)行分段計(jì)算、綜合評(píng)定的方法。通過(guò)對(duì)海洋立管常見(jiàn)風(fēng)險(xiǎn)因素進(jìn)行收集整理，建立風(fēng)險(xiǎn)評(píng)價(jià)模型，采用灰色關(guān)聯(lián)度理論對(duì)海洋立管進(jìn)行分段風(fēng)險(xiǎn)評(píng)估。根據(jù)各主要風(fēng)險(xiǎn)因素對(duì)發(fā)生風(fēng)險(xiǎn)概率的關(guān)聯(lián)度，計(jì)算出各段立管在各主要風(fēng)險(xiǎn)因素下的失效權(quán)重，將所得權(quán)重結(jié)合立管各區(qū)段在多種主要風(fēng)險(xiǎn)因素下的失效關(guān)聯(lián)系數(shù)，進(jìn)一步計(jì)算出立管各區(qū)段的失效概率。這種評(píng)估方法為海洋立管的風(fēng)險(xiǎn)評(píng)估提供了一種新的思路，對(duì)加強(qiáng)海洋立管各區(qū)段管理與維護(hù)提供一定依據(jù)。

灰色理論；關(guān)聯(lián)度；海洋立管；風(fēng)險(xiǎn)；權(quán)重

隨著海上油田的大力開(kāi)發(fā)，各種石油開(kāi)采和集輸設(shè)備迅速投入生產(chǎn)和使用。海上油氣開(kāi)發(fā)的重點(diǎn)和難點(diǎn)是確保油氣集輸設(shè)備的安全可靠性，而海洋立管作為油氣輸送的關(guān)鍵設(shè)備，在整個(gè)石油開(kāi)采和集輸設(shè)備中占很大比例，其在役的安全可靠性是保證石油順利開(kāi)采和集輸?shù)年P(guān)鍵。自從灰色理論在1982被國(guó)外學(xué)者提出后，在國(guó)內(nèi)得到較快的應(yīng)用和推廣，但是我國(guó)對(duì)海上油田的開(kāi)發(fā)起步較晚，對(duì)海上油氣開(kāi)發(fā)設(shè)備可靠性相關(guān)評(píng)價(jià)方法也極為匱乏。由于海洋立管的每段工作海域都有它所對(duì)應(yīng)的特定工作環(huán)境，從整體上進(jìn)行宏觀評(píng)價(jià)可行性差。因此，只有結(jié)合我國(guó)海洋環(huán)境具體環(huán)境和設(shè)備在役狀況，有針對(duì)性的進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià)和理論研究，將定性分析轉(zhuǎn)變?yōu)槎糠治觯拍軓恼嬲饬x上實(shí)現(xiàn)海洋管道設(shè)備的可靠性分析。

海洋立管在油氣輸送過(guò)程中受環(huán)境因素影響較大，不確定風(fēng)險(xiǎn)因素較多，對(duì)海洋立管工作穩(wěn)定性造成較大影響。對(duì)于少數(shù)據(jù)、小樣本、信息不完全和經(jīng)驗(yàn)缺乏的不確定性問(wèn)題，通常采用灰色理論來(lái)處理[1-5]?；疑P(guān)聯(lián)理論中的灰色評(píng)價(jià)法利用灰色關(guān)聯(lián)作為直接測(cè)度，充分利用白化信息，減少評(píng)估誤差，使得評(píng)價(jià)結(jié)果更加科學(xué)合理，與實(shí)際情況更符合。本文通過(guò)運(yùn)用灰色關(guān)聯(lián)理論，結(jié)合海洋立管各段的具體工況和自身特點(diǎn)將立管進(jìn)行分段計(jì)算，得到各段立管的定量分析結(jié)果，充分反映出管道整體可靠性概況[6-10]。

1　關(guān)聯(lián)度理論與權(quán)重

(1) 選取參考序列x0={x0(1),x0(2),…,x0(n)}，把與x0進(jìn)行比較的序列xj叫做比較序列xj={xj(1),xj(2),…,xj(n)},j=1,2,…,m。通常將序列中具有代表性的數(shù)據(jù)作為參考數(shù)列，也稱(chēng)為最優(yōu)集或母序列。

(3) 求差序列：

(4) 求兩極差：

(5) 求關(guān)聯(lián)系數(shù)：

式中，ξ為分辨系數(shù)，一般取ξ=0.5。為方便下一步計(jì)算可以將關(guān)聯(lián)系數(shù)以矩陣R的形式表示[11]。

2　風(fēng)險(xiǎn)因素模型與灰色關(guān)聯(lián)分析

影響立管穩(wěn)定工作的風(fēng)險(xiǎn)因素有很多，結(jié)合歷史記錄與實(shí)際經(jīng)驗(yàn)找出事故類(lèi)型，收集整理風(fēng)險(xiǎn)因素，將各風(fēng)險(xiǎn)因素按照W. K. Muhlbauer[6]提供的方法，將立管風(fēng)險(xiǎn)因素劃分為：腐蝕因素[7-9]、載荷因素、第三方破壞、溫度和磨損。每類(lèi)風(fēng)險(xiǎn)因素都包含著相應(yīng)的影響因素。海洋立管的風(fēng)險(xiǎn)因素模型如圖1所示。

(2) 根據(jù)灰色關(guān)聯(lián)理論具體分析步驟，建立海洋立管灰色關(guān)聯(lián)分析流程圖，如圖2所示。

運(yùn)用灰色關(guān)聯(lián)分析法計(jì)算，計(jì)算各種主要事故因素的關(guān)聯(lián)度，通過(guò)比較各個(gè)因素關(guān)聯(lián)度，對(duì)各因素的風(fēng)險(xiǎn)概率進(jìn)行大小排序。依據(jù)所選用的參考序列(最優(yōu)集)的性質(zhì)，得出的各因素關(guān)聯(lián)度越大表示該因素與最優(yōu)集的性質(zhì)越接近，若以最小風(fēng)險(xiǎn)概率作為參考序列，因素關(guān)聯(lián)度越大則表示該因素發(fā)生風(fēng)險(xiǎn)的概率越小[10-13]，與最小概率的性質(zhì)接近。反之，表示風(fēng)險(xiǎn)發(fā)生概率越大。

圖1　海洋立管風(fēng)險(xiǎn)因素模型

Fig.1Marine riser risk factor model

圖2　灰色關(guān)聯(lián)分析流程

Fig.2Grey relational analysis flowchart

3　實(shí)例應(yīng)用

某深海立管根據(jù)立管作業(yè)環(huán)境分為海上大氣區(qū)、浪濺區(qū)、潮差區(qū)、全浸區(qū)和海底泥層區(qū)。根據(jù)海洋立管各區(qū)段失效記錄數(shù)據(jù)和實(shí)際經(jīng)驗(yàn)，給定管道在不同區(qū)段不同風(fēng)險(xiǎn)因素下的失效概率如表1所示。用灰色關(guān)聯(lián)理論計(jì)算各個(gè)風(fēng)險(xiǎn)因素的關(guān)聯(lián)度以及各區(qū)段在這些風(fēng)險(xiǎn)影響下的失效概率。

3.1風(fēng)險(xiǎn)因素關(guān)聯(lián)分析與權(quán)重計(jì)算

(1) 建立參考序列。選取各風(fēng)險(xiǎn)因素最小概率數(shù)據(jù)做為參考序列，則x0=(0.000 479,0.000 023,0.000 019,0.000 021,0.000 020)

(2) 根據(jù)均值像公式將比較序列均值化處理,得：

表1　海洋立管失效概率數(shù)據(jù)

(3) 根據(jù)求差序列公式，求得差序列：

Δ1=(3.069 07,1.278 48,1.297 26,0.773 17,0.720 18);

Δ2=(4.169 57,1.377 17,1.359 06,2.578 06,0.144 73);

Δ3=(2.505 17,1.692 77,1.536 06,0.210 56,0.065 66);

Δ4=(3.194 07,0.568 37,0.640 76,1.027 87,1.956 96);

Δ5=(2.172 97,1.397 87,1.410 66,0.382 77, 0.018 44)。

(4) 由兩極差公式，求得兩極差：

M=4.169 57；m=0.018 44。

(5) 求關(guān)聯(lián)系數(shù)，根據(jù)關(guān)聯(lián)系數(shù)求解公式得：

ξ1=(0.408 08,0.625 35,0.621 88,0.735 92,0.749 82)；

ξ2=(0.336 28,0.607 52,0.610 71,0.451 06,0.943 35)；

ξ3=(0.458 22,0.516 77,0.530 86,0.716 29,0.978 04)；

ξ4=(0.398 42,0.792 72,0.871 67,0.975 70,0.520 37)；

ξ5=(0.493 97,0.603 90,0.601 70,0.852 35,1)。

(6) 運(yùn)用關(guān)聯(lián)度公式計(jì)算關(guān)聯(lián)度和第一層各風(fēng)險(xiǎn)因素權(quán)重：

r1=0.628 212；r2=0.589 786；r3=0.640 037；

r4=0.711 778；r5=0.710 388。

通過(guò)比較得：

r4>r5>r3>r1>r2。

計(jì)算風(fēng)險(xiǎn)權(quán)重：

根據(jù)權(quán)重計(jì)算公式通過(guò)歸一化處理得權(quán)重：

w1=0.207 7；w2=0.221 3；w3=0.203 9；

w4=0.183 4；w5=0.183 7。

最優(yōu)比較數(shù)列是根據(jù)最小風(fēng)險(xiǎn)概率數(shù)據(jù)確定的，數(shù)據(jù)越大表示與所比較數(shù)列的關(guān)聯(lián)度越高，發(fā)生風(fēng)險(xiǎn)的可能性就越小，反之發(fā)生風(fēng)險(xiǎn)的可能性就越大。因此，根據(jù)所得結(jié)果可判斷發(fā)生風(fēng)險(xiǎn)的可能性依次為：載荷>缺陷>第三方破壞>磨損>溫度。根據(jù)實(shí)際應(yīng)用中管道失效情況，可知所得結(jié)果與實(shí)際基本相符。依此算法同樣可計(jì)算其他層次的風(fēng)險(xiǎn)因素的關(guān)聯(lián)度。

3.2海洋立管區(qū)段風(fēng)險(xiǎn)關(guān)聯(lián)度計(jì)算

(1) 建立參考序列。為確定海洋立管各區(qū)段風(fēng)險(xiǎn)大小情況，選取風(fēng)險(xiǎn)因素中的風(fēng)險(xiǎn)概率最大值作為參考值，根據(jù)表1中的數(shù)據(jù)，確定參考序列x0=(0.000 791,0.001 082,0.000 148,0.000 058,0.000 275)

(2) 均值化處理比較序列:

(3) 求差序列:

Δ1=(1.179 95, 2.136 49,0.301 09,0.007 09,0.651 29);

Δ2=(0.532 30, 0.569 17,0.099 05, 0.064 54,0.005 27);

Δ3=(0.566 07,0.582 79,0.067 17, 0.060 08,0.012 56);

Δ4=(0.199 99,0.823 53,0.224 92,0.027 98,0.367 72);

Δ5=(2.389 02,2.187 87,0.153 07,0.369 17,0.414 33)。

(4) 求兩極差:

M=2.389 02;m=0.005 27。

(5) 求關(guān)聯(lián)系數(shù):

ε1=(0.505 29,0.360 19,0.802 21,0.998 48,0.650 01);

ε2=(0.694 79,0.680 27,0.927 50,0.952 92,1);

ε3=(0.681 47,0.675 06,0.950 94,0.956 31,0.993 96);

ε4=(0.860 36,0.594 53,0.845 25,0.981 42,0.767 99);

ε5=(0.334 80,0.354 71,0.890 32,0.707 28,0.645 74);

(6) 求關(guān)聯(lián)度:

根據(jù)關(guān)聯(lián)系數(shù)的系數(shù)矩陣R

wj=(0.207 7,0.221 3,0.203 9,0.183 4,0.183 7)

為方便對(duì)比，按照均值化求關(guān)聯(lián)度得到rj，rj=(0.663 236,0.851 096,0.851 548,0.809 91,0.618 57)。

4　結(jié)論

海洋立管的每段工作區(qū)域都有它對(duì)應(yīng)的特定海況，缺乏直接可行的宏觀評(píng)價(jià)方法。本文針對(duì)整根立管評(píng)定可行性差的問(wèn)題，提出了一套將整根立管進(jìn)行分段計(jì)算，綜合評(píng)定的方法。通過(guò)深入了解具體環(huán)節(jié)，建立合適的評(píng)價(jià)內(nèi)容和方法，得到最終的評(píng)價(jià)結(jié)果。在計(jì)算過(guò)程中，為了準(zhǔn)確的掌握各個(gè)風(fēng)險(xiǎn)因素對(duì)海洋立管失效影響的大小，通過(guò)構(gòu)建海洋立管風(fēng)險(xiǎn)因素模型，利用灰色關(guān)聯(lián)理論，確定了各個(gè)風(fēng)險(xiǎn)的關(guān)聯(lián)度，進(jìn)而算得各風(fēng)險(xiǎn)因素對(duì)海洋立管失效的概率權(quán)重。將所得各風(fēng)險(xiǎn)概率權(quán)重與海洋立管不同區(qū)段風(fēng)險(xiǎn)因素的關(guān)聯(lián)系數(shù)結(jié)合，得到海洋立管在不同區(qū)段的失效概率。在計(jì)算各個(gè)風(fēng)險(xiǎn)因素對(duì)海洋立管失效概率權(quán)重過(guò)程，摒棄了傳統(tǒng)的專(zhuān)家打分或者主觀臆斷法，而是將實(shí)際應(yīng)用中各風(fēng)險(xiǎn)概率數(shù)據(jù)與灰色關(guān)聯(lián)理論相結(jié)合。這種計(jì)算權(quán)重方法，大大減小了主觀判斷帶的誤差，提高了計(jì)算精度。計(jì)算結(jié)果表明，風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果與海洋立管實(shí)際具體失效情況基本符合，該方法具有較高的應(yīng)用價(jià)值，可作為海洋立管風(fēng)險(xiǎn)評(píng)估的一種手段，為海洋立管的可靠性評(píng)價(jià)以及立區(qū)段維護(hù)理提供可靠的依據(jù)。

[1]李建華,元東民,俞樹(shù)榮,等.基于灰色理論的長(zhǎng)輸管道風(fēng)險(xiǎn)分析[J].石油機(jī)械,2006,34(11):24-26.

Li Jianhua,Yuan Dongmin,Yu Shurong,et al. Risk analysis of oil and gas pipelines based on grey theory[J]. China Petroleum Machinery, 2006,34(11): 24-26.

[2]何學(xué)秋.安全工程學(xué)[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,2000.

[3]劉思峰.灰色系統(tǒng)理論與應(yīng)用[M].北京:科學(xué)出版社,2007.

[4]鄧聚龍.灰色系統(tǒng)基本方法[M].武漢:華中科技大學(xué)出版社,2005.

[5]鄧聚龍.灰色系統(tǒng)理論教程[M].武漢:華中理工大學(xué)出版社,1992.

[6]Muhlbauer W K. Pipeline Risk Management Manual.2nded[M].Houston:Gulf Publishing Company,1996.

[7]程浩力,劉德俊.模型預(yù)測(cè)海底管道腐蝕速率軟件的開(kāi)發(fā)[J].遼寧石油化工大學(xué)學(xué)報(bào),2011, 31(3):58-61.

Cheng Haoli,Liu Dejun. Application of GM (1, 1) on the prediction software of corrosive rate of submarine pipelines[J]. Journal of Liaoning Shihua University, 2011,31(3):58-61.

[8]孫春一,石彬.基于分形插值方法的油氣管道腐蝕灰色預(yù)測(cè)[J].遼寧石油化工大學(xué)學(xué)報(bào),2009,29(3):52-54.

Sun Chunyi,Shi Bin. Grey predication of corrosion of the oil and gas pipes based on fractal interpolation method[J]. Journal of Liaoning University of Petroleum & Chemical Technology, 2009,29(3):52-54.

[9]代真,牛曉光,沈士明.流體力學(xué)因素對(duì)空冷器沖蝕影響的定量分析[J].石油化工高等學(xué)校學(xué)報(bào),2007,20(4):85-88.

Dai Zhen,Niu Xiaoguang,Shen Shiming. Quantitative analysis for effect of hydrodynamic factors on erosion-corrosion of air cooler[J]. Journal of Petrochemical Universities,2007,20(4):85-88.

[10]劉展,王智平,俞樹(shù)榮,等.基于灰色理論的壓 力管道風(fēng)險(xiǎn)分析技術(shù)研究[J].石油化工設(shè)備,2008,37(2):1-3.

Liu Zhan,Wang Zhiping,Yu Shurong,et al.Study of pressure pipeline risk analysis technology based on gray theory—theory and key technology of pressure pipeline risk management(4)[J]. Petro-Chemical Equipment, 2008,37(2):1-3.

[11]張圣柱,吳宗之,張健,等.油氣管道選線和風(fēng)險(xiǎn)評(píng)價(jià)相關(guān)法規(guī)與方法[J].油氣儲(chǔ)運(yùn),2012, 31(9):663-669.

Zhang Shengzhu, Wu Zongzhi, Zhang Jian,et al.Relevant regulations and methods for route selection of oil and gas pipeline and risk assessment[J].Oil & Gas Storage and Transportation,2012,31(9):663-669.

[12]李小波.油其長(zhǎng)輸管道風(fēng)險(xiǎn)管理的數(shù)據(jù)管理[J].石油礦場(chǎng)機(jī)械,2013,42(5):27-30.

Li Xiaobo. Data management in oil-and-gas pipeline risk supervision[J].Oil and Field Equipment, 2013,42(5):27-30.

[13]譚成仟,宋子齊,吳向紅.儲(chǔ)層油氣產(chǎn)能的灰色理論預(yù)測(cè)方法[J].系統(tǒng)理論與實(shí)踐,2010, 10(10):101-106.

Tan Chengqian,Song Ziqi,Wu Xianghong. The predicting method of reservoir production capacity by grey system theory[J]. Systems Engineering-Theory & Practice,2010, 10(10):101-106.

[14]王一飛,潘志遠(yuǎn),黃小平,等.深海立管渦激振動(dòng)疲勞損傷影響因素分析[J].船舶力學(xué),2006,10(5):76-83.

Wang Yifei,Pan Zhiyuan,Huang Xiaoping,et al. Effects of various factors on the VIV-induced fatigue damage in deep sea risers[J]. Journal of Ship Mechanics, 2006,10(5):76-83.

(編輯王亞新)

The Risk Assessment of the Segmented Marine Riser Based on Gray Relevance Theory

Zheng Yunhu， Jiang Feng， Du Chaofei

(CollegeofPetrochemicalEngineering,LanzhouUniversityofTechnology,LanzhouGansu730050,China)

Each work area of the marine riser has its corresponding specific sea conditions, and the feasibility of macro-evaluate is poor. For the feasibility of the whole problem of poor riser assessment, a method was proposed to segment the entire riser calculations, and assessment comprehensively. The marine riser common risk factors were collected to establish risk assessment model, and the gray correlation theory was used to assessment segmented riskof marine riser. According to the major risk factors associated with the risk of occurrence probability, for each segment, the invalidation problem and unreasonable weightswere calculated, the probability of failure of each section of the riser was further calculated by combining the weight with failure correlation coefficient of each section in a variety of major risk factors . This method of risk assessment of the marine riser provides a new way to strengthen the management and maintenance for marine riser of each section.

Gray theory; Correlation; Marine riser; Risk; Weights

1006-396X(2016)03-0069-05

2013-12-06

2014-03-12

甘肅省教育廳資助項(xiàng)目(1205ZTC067)；國(guó)家質(zhì)檢公益項(xiàng)目(201210026，201310152)。

鄭運(yùn)虎(1986-),男,碩士，從事承壓設(shè)備安全可靠性研究；E-mail:921238461@qq.com。

TQ050.2

Adoi:10.3969/j.issn.1006-396X.2016.03.014

投稿網(wǎng)址：http://journal.lnpu.edu.cn