姜明軍，魯欣衛(wèi)，朱魯方，李海壽，張偉娜，朱萬(wàn)林，馬國(guó)印，李晨安

(1.海洋石油工程股份有限公司，天津 300451； 2.中海石油(中國(guó))有限公司天津分公司，天津 300451；3.中法渤海地質(zhì)服務(wù)有限公司，天津 300451； 4.合肥通用機(jī)械研究院有限公司，合肥 230031)

0 引言

隨著我國(guó)海洋石油裝備的不斷發(fā)展，水下生產(chǎn)系統(tǒng)加采油船逐漸成為我國(guó)海洋油氣資源開(kāi)采的主要模式[1]。海洋油氣資源的開(kāi)采過(guò)程可簡(jiǎn)述如下：水下生產(chǎn)系統(tǒng)采集原油并通過(guò)海管外輸至采油船，在船上經(jīng)油氣分離、除沙、脫水、除硫、化學(xué)藥劑回收等一系列工藝處理流程后儲(chǔ)存在船體的儲(chǔ)油倉(cāng)內(nèi)。原油產(chǎn)品屬于易燃易爆物質(zhì)。當(dāng)儲(chǔ)油倉(cāng)內(nèi)游離的碳?xì)浠衔锶芏?、助燃氧氣溶度處于油氣爆炸范圍且引燃溫度達(dá)到燃爆下限時(shí)，會(huì)發(fā)生儲(chǔ)油倉(cāng)的燃爆現(xiàn)象[2]。為避免發(fā)生燃爆事故，必須為采油船設(shè)置惰氣發(fā)生器系統(tǒng)(inert gas generator system,IGGS)。

IGGS作為海洋油氣資源開(kāi)發(fā)的重要設(shè)備，在原油存儲(chǔ)與外輸環(huán)節(jié)均發(fā)揮重要作用。該系統(tǒng)一旦發(fā)生故障將會(huì)導(dǎo)致油田停產(chǎn)，為安全生產(chǎn)埋下巨大隱患。另一方面，IGGS的故障診斷往往過(guò)分依賴(lài)機(jī)修師的工作經(jīng)驗(yàn)[3]，且故障排查過(guò)程耗時(shí)久、效果差。綜上所述，快速、準(zhǔn)確的故障源定位，對(duì)于協(xié)助維修人員排除故障具有十分重要的意義。

為解決上述問(wèn)題，本文提出了一種利用貝葉斯網(wǎng)絡(luò)進(jìn)行IGGS故障診斷的方法：首先，簡(jiǎn)要介紹了IGGS的組成結(jié)構(gòu)、運(yùn)行機(jī)制；其次，闡述了有關(guān)貝葉斯網(wǎng)絡(luò)的基本知識(shí)并引入Netica軟件；最后，結(jié)合南海某采油船項(xiàng)目，給出IGGS故障診斷的案例分析。

1 IGGS的原理及組成

IGGS通過(guò)燃燒天然氣或燃油與空氣的混合物產(chǎn)生惰性氣體。空氣中主要包含了79%的氮?dú)夂?1%氧氣，燃燒過(guò)程消耗掉大部分氧氣，產(chǎn)生二氧化碳和水蒸氣，以及少量的一氧化碳和氫氣[4]。反應(yīng)原理如式(1)所示：

(CH)n+O2→CO2+H2O+CO+H2

(1)

IGGS產(chǎn)生的廢氣經(jīng)洗滌、除塵、降溫等工序處理后供采油船使用。此時(shí)，產(chǎn)生氣體的主要成分是氮?dú)饧岸趸?。雖然各廠家生產(chǎn)的IGGS外形差異很大，但其工作原理基本相同，功能、結(jié)構(gòu)十分類(lèi)似[5]。IGGS主要包括惰氣發(fā)生器、單元控制面板(unit control panel,UCP)控制柜、甲板水封撬及P/V破斷器。IGGS結(jié)構(gòu)如圖1所示。

圖1 IGGS結(jié)構(gòu)圖

IGGS中，最左側(cè)是燃?xì)夤┙o入口和燃油供給入口，為燃燒腔室提供燃?xì)饧叭加?。氣閥單元與燃油泵決定燃?xì)馀c燃油供給管路的通斷。風(fēng)機(jī)為燃燒腔室提供新鮮空氣。燃燒反應(yīng)在本地控制盤(pán)的調(diào)節(jié)下發(fā)生在惰氣發(fā)生器的燃燒腔室內(nèi)。燃燒模式分為燃油、燃?xì)夂突旌先N。

①燃油模式。燃油進(jìn)入燃燒腔室前被噴嘴噴散使其霧化，水平方向的霧狀燃油與切向新鮮空氣按照一定比例充分混合，并由引燃器點(diǎn)火引燃。

②燃?xì)饽Ｊ健Ｈ細(xì)獠恍枰F化，直接進(jìn)入燃燒腔室與新鮮空氣混合后燃燒。

③混合模式。燃油與燃?xì)馔瑫r(shí)作為燃料進(jìn)入燃燒腔室燃燒。

惰氣發(fā)生器的出口管線上裝有氧分析儀，測(cè)定排出氣的氧含量并將數(shù)據(jù)回傳至UCP控制柜。UCP控制柜與IGGS啟動(dòng)器、本地控制盤(pán)相互通信共同控制IGGS的啟停過(guò)程，使其安全、穩(wěn)定地運(yùn)行。排氣扇將燃燒腔室內(nèi)的惰氣排出。合格的惰氣通過(guò)B口經(jīng)甲板水封撬、P/V破斷器進(jìn)入船體原油倉(cāng)，起到封倉(cāng)的作用。當(dāng)惰氣的含氧量較高時(shí)，不合格的惰氣將直接排放到空氣中。Netica軟件僅支持英文字符輸入。為便于說(shuō)明，下文給出主要部件中英文名稱(chēng)對(duì)照表。

IGGS各部分中英文名稱(chēng)對(duì)照表如表1所示。

表1 IGGS各部分中英文名稱(chēng)對(duì)照表

2 貝葉斯網(wǎng)絡(luò)及Netica軟件

貝葉斯網(wǎng)絡(luò)(Bayesian network,BN)是概率圖模型的一種，結(jié)構(gòu)為一種有向無(wú)環(huán)圖[6]。BN由節(jié)點(diǎn)和節(jié)點(diǎn)之間的有向邊組成。節(jié)點(diǎn)代表了隨機(jī)變量其集合，記為I；有向邊代表了節(jié)點(diǎn)之間的條件概率，其集合記為E。則BN可簡(jiǎn)記為G=(I,E)。假設(shè)集合I中共有N個(gè)隨機(jī)變量Xi∈I(i=1,2,…N)，根據(jù)概率乘法的鏈?zhǔn)椒▌t，N個(gè)隨機(jī)變量的聯(lián)合概率分布如式(2)所示：

p(X1,X2,…,Xn)=p1(X1)p2(X2|X1)p3(X3|X1,X2)…

(2)

式中：pi(Xi)為節(jié)點(diǎn)Xi的概率。

根據(jù)有向無(wú)環(huán)圖的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)，網(wǎng)絡(luò)中任意節(jié)點(diǎn)的概率只與其父節(jié)點(diǎn)的概率有關(guān)。其聯(lián)合概率分布可簡(jiǎn)化為：

(3)

式中：F(Xi)為節(jié)點(diǎn)Xi的所有父節(jié)點(diǎn)；pi[Xi|F(Xi)]為父節(jié)點(diǎn)與節(jié)點(diǎn)Xi的條件概率。

式(3)表明：在BN中，聯(lián)合概率分布可表示為所有節(jié)點(diǎn)與其父節(jié)點(diǎn)條件概率的乘積。BN可以很方便地構(gòu)建事物間的不確定性、簡(jiǎn)化網(wǎng)絡(luò)推理過(guò)程，已經(jīng)被成功應(yīng)用于流程工業(yè)生產(chǎn)過(guò)程[7]、汽輪機(jī)[8]以及導(dǎo)彈發(fā)射過(guò)程的故障診斷[9]。

式(3)使概率計(jì)算程序化，可以很方便地利用計(jì)算機(jī)計(jì)算每一個(gè)結(jié)點(diǎn)的概率。Netica是Norsys軟件公司的一款基于BN的統(tǒng)計(jì)推斷軟件，已被各大公司、政府機(jī)構(gòu)及醫(yī)院用于決策分析、因果推斷，是各領(lǐng)域應(yīng)用較多的貝葉斯分析軟件[10-11]。Netica軟件圖形化的操作界面對(duì)用戶(hù)十分友好，工程人員稍加培訓(xùn)即可熟練掌握。軟件中可以很方便地定義節(jié)點(diǎn)、設(shè)置節(jié)點(diǎn)之間的條件概率、一鍵完成BN的推理過(guò)程，并以數(shù)字加圖形的方式顯示結(jié)果。針對(duì)先驗(yàn)概率不完整的情況，軟件支持導(dǎo)入測(cè)試數(shù)據(jù)訓(xùn)練網(wǎng)絡(luò)參數(shù)。使用軟件創(chuàng)建BN的過(guò)程主要包括：建立網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)、確定條件概率表(conditional probability tables,CPT)、使用網(wǎng)絡(luò)進(jìn)行推理。行業(yè)專(zhuān)家基于對(duì)研究對(duì)象工作機(jī)制、操作流程的理解，將其拆分為若干部分并確定各部分間的因果關(guān)系，在此基礎(chǔ)上建立BN的拓?fù)浣Y(jié)構(gòu)。條件概率表的確定或基于專(zhuān)家的先驗(yàn)知識(shí)，或基于測(cè)試數(shù)據(jù)。以下舉例說(shuō)明如何使用Netica創(chuàng)建一個(gè)簡(jiǎn)單的BN并進(jìn)行推理。

2.1 創(chuàng)建網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)

在IGGS中，當(dāng)有燃油供給且燃油泵正常工作時(shí)，惰氣發(fā)生器的燃燒腔室才會(huì)有燃油進(jìn)入。據(jù)此可創(chuàng)建BN的拓?fù)浣Y(jié)構(gòu)。燃油輸入系統(tǒng)的BN如圖2所示。

圖2 燃油輸入系統(tǒng)的BN示意圖

2.2 確定條件概率表

依據(jù)專(zhuān)家經(jīng)驗(yàn)及歷史數(shù)據(jù)，確定節(jié)點(diǎn)Fuel Oil Supply (S)的條件概率為p(S=Yes)=0.7、p(S=No)=0.3，表示燃油有供給的概率為70%、無(wú)供給的概率為30%。節(jié)點(diǎn)Fuel Oil Pump (P)的條件概率為p(P=Work)=0.8、p(P=Fault)=0.2，表示燃油泵正常工作的概率為80%、故障的概率為20%。節(jié)點(diǎn)Fuel Oil Input (I)的CPT如表2所示。

表2 節(jié)點(diǎn)Fuel Oil Input的CPT

表2中，第一行數(shù)據(jù)表示在燃油正常供給且燃油泵正常工作的條件下，燃燒腔室內(nèi)有燃油的概率為95%，無(wú)燃油的概率為5%。Netica軟件界面僅顯示節(jié)點(diǎn)的邊緣概率，如圖2所示，節(jié)點(diǎn)Fuel Oil Supply (S)與Fuel Oil Pump (P)沒(méi)有父節(jié)點(diǎn)其邊緣概率等于條件概率。節(jié)點(diǎn)Fuel Oil Input (I)的邊緣概率依據(jù)式(3)計(jì)算：

p(I=Y)=p(I=Y,S=Y,P=W)+p(I=Y,S=Y,P=F)+p(I=Y,S=N,P=W)+p(I=Y,S=N,P=F)=p(I=Y|S=Y,P=W)P(S=Y,P=W)+p(I=Y|S=Y,P=F)P(S=Y,P=F)+p(I=Y|S=N,P=W)P(S=N,P=W)+p(I=Y|S=N,P=F)P(S=N,P=F)

(4)

式中：Y、N、W、F分別為Yes、No、Work、Fault的首字母縮寫(xiě)。

圖2所示的BN為Head-to-Head結(jié)構(gòu)，在Fuel Oil Input (I)未知的情況下，隨機(jī)節(jié)點(diǎn)Fuel Oil Supply (S)與Fuel Oil Pump (P)相互獨(dú)立。式(4)可化簡(jiǎn)為：

p(I=Y)=p(I=Y|S=Y,P=W)P(S=Y)P(P=W)+p(I=Y|S=Y,P=F)P(S=Y)P(P=F)+p(I=Y|S=N,P=W)P(S=N)P(P=W)+p(I=Y|S=N,P=F)P(S=N)P(P=F)=0.95×0.7×0.8+0.02×0.7×0.2+0.02×0.3×0.8+0×0.3×0.2=0.539 6

(5)

根據(jù)概率的公理化定義，可得：

p(I=N)=1-p(I=Y)=1-0.539 6=0.460 4

(6)

2.3 使用網(wǎng)絡(luò)進(jìn)行推理

在Netica軟件界面中，單擊Fuel Oil Input (I)節(jié)點(diǎn)的Yes狀態(tài)，軟件自動(dòng)更改其值為100%，所有節(jié)點(diǎn)的邊緣概率值均發(fā)生變化。節(jié)點(diǎn)Fuel Oil Input置為Yes=100%如圖3所示。

圖3 節(jié)點(diǎn)Fuel Oil Input 置為Yes=100%

由圖3可知，若燃燒腔室有燃油輸入，則燃油有供給的概率是99.1%，燃油泵正常工作的概率是99.5%。

同理，單擊Fuel Oil Input (I)節(jié)點(diǎn)的No狀態(tài)，軟件自動(dòng)更改其值為100%。節(jié)點(diǎn)Fuel Oil Input置為No=100%如圖4所示。

圖4 節(jié)點(diǎn)Fuel Oil Input置為No=100%

由圖4可知，若燃燒腔室內(nèi)無(wú)燃油輸入，則燃油無(wú)供給的概率是64.1%、燃油泵故障的概率是42.8%。燃油無(wú)供給的概率高于燃油泵故障的概率，據(jù)此設(shè)備維修人員應(yīng)首先檢查燃油供給狀態(tài)，再查看燃油箱內(nèi)是否有燃油。若確定燃油有供給，則維修人員再確認(rèn)燃油泵是否故障。

3 案例分析

根據(jù)第一節(jié)介紹的IGGS的組成框圖創(chuàng)建BN。為便于描述，該案例僅對(duì)惰氣發(fā)生器撬塊TS-X-3803部分進(jìn)行建模，甲板水封撬(TS-T-3810)及P/V破斷器(TS-T-3820)暫不予考慮。

在Netica軟件中創(chuàng)建的惰氣發(fā)生器的BN如圖5所示。

圖5 惰氣發(fā)生器的BN示意圖

此時(shí)，惰氣發(fā)生器處于混合模式。節(jié)點(diǎn)Combustion Chamber的CPT如表3所示。

表3 節(jié)點(diǎn)Combustion Chamber的CPT

在Netica軟件界面中，單擊Inert Gas Output節(jié)點(diǎn)的No狀態(tài)，軟件自動(dòng)更改其值為100%，表示無(wú)惰氣輸出，惰氣發(fā)生器處于故障狀態(tài)。觀察其父節(jié)點(diǎn)，Combustion Chamber的故障概率變?yōu)?9.3%，而Extraction Fan的故障概率變?yōu)?2.6%。若能夠確定Combustion Chamber故障，則繼續(xù)單擊該節(jié)點(diǎn)的Fault狀態(tài)，其父節(jié)點(diǎn)的概率發(fā)生變化，Local Control Panel故障的概率變?yōu)?0.6%，燃燒腔室內(nèi)無(wú)燃?xì)夤?yīng)的概率為35.9%，無(wú)燃油供應(yīng)的概率為34.8%，無(wú)新鮮空氣供應(yīng)的概率為49.9%。

以故障概率值為基準(zhǔn)，從大到小依次確認(rèn)燃燒腔室內(nèi)新鮮空氣輸入、燃?xì)廨斎搿⑷加洼斎胧欠裾?。若確認(rèn)無(wú)新鮮空氣供應(yīng)，則單擊Fresh Air Input節(jié)點(diǎn)的No狀態(tài)。融合證據(jù)后的BN如圖6所示。圖6中，在惰氣發(fā)生器各部分的故障概率中Blower Fan節(jié)點(diǎn)的故障概率最大，為90.5%，設(shè)備維修人員需確認(rèn)風(fēng)機(jī)是否正常工作。按照上述故障診斷步驟，從結(jié)果開(kāi)始依次加入已經(jīng)確認(rèn)的各種證據(jù)，逐步倒推，直到排查出設(shè)備的故障部位。

圖6 融合證據(jù)后的BN示意圖

4 結(jié)論

針對(duì)設(shè)置有IGGS的采油船(如浮式生產(chǎn)儲(chǔ)卸油裝置、半潛式采油船)，可利用BN對(duì)IGGS進(jìn)行故障診斷。BN是概率模型與圖模型的有機(jī)結(jié)合體，便于人們以類(lèi)似于故障樹(shù)的模式排查IGGS的故障。在IGGS的維修過(guò)程中，相關(guān)維修人員可借助Netica軟件創(chuàng)建BN并逐步加入證據(jù)進(jìn)行推理，從而快速確定故障源頭、縮短維修時(shí)間，以節(jié)約海洋石油開(kāi)采成本。