天然氣凈化廠過程風(fēng)險監(jiān)測預(yù)警模型構(gòu)建方法研究

閻紅巧，曹瑜，吳順成，王順義，樊志強(qiáng)

(中國石油集團(tuán)安全環(huán)保技術(shù)研究院有限公司，北京 102200)

0 引言

天然氣作為高效的清潔能源，應(yīng)用市場廣闊。剛開采出的天然氣通常含有硫化氫、有機(jī)硫、二氧化碳、水及其他雜質(zhì)，需對其進(jìn)行凈化處理以滿足集輸要求，供用戶安全使用。保證天然氣凈化過程的安全平穩(wěn)運(yùn)行，是后續(xù)天然氣經(jīng)營的基礎(chǔ)。目前，天然氣凈化廠生產(chǎn)運(yùn)行風(fēng)險管控主要依靠DCS 和SCADA 等控制系統(tǒng)閾值報警，安全聯(lián)鎖系統(tǒng)(SIS)和緊急聯(lián)鎖系統(tǒng)(ESD) 聯(lián)鎖和人工定期巡檢的方式來防范生產(chǎn)運(yùn)行過程中的風(fēng)險[1]。但這些方法存在以下問題：(1)控制系統(tǒng)報警泛濫，報警有效性較低，導(dǎo)致操作員無法及時準(zhǔn)確響應(yīng)而發(fā)生事故；(2) 消警處置和定期巡檢受現(xiàn)場人員的專業(yè)程度、主觀思維和經(jīng)驗影響大，易因漏檢、誤判、誤操作、處置不及時而引發(fā)事故；(3)SIS系統(tǒng)、ESD 系統(tǒng)啟動會引起非計劃停機(jī)，導(dǎo)致生產(chǎn)力損失，生產(chǎn)運(yùn)行成本增加。

生產(chǎn)運(yùn)行風(fēng)險預(yù)警是通過深層挖掘生產(chǎn)過程參數(shù)與事故先兆之間的關(guān)聯(lián)關(guān)系，構(gòu)建表征風(fēng)險的預(yù)警模型，實現(xiàn)在裝置發(fā)生事故之前甚至是發(fā)生非正常工況之前預(yù)警，及早告知操作人員及管理人員，為風(fēng)險處置爭取更多時間，從而極大降低安全事故的發(fā)生和連帶損壞的風(fēng)險，提升系統(tǒng)的可靠性[2]。中國石化安全工程研究院多次開展了石化裝置的生產(chǎn)過程風(fēng)險的動態(tài)監(jiān)測預(yù)警技術(shù)的研究。包括基于數(shù)據(jù)驅(qū)動的異常監(jiān)測預(yù)警方法和基于專家系統(tǒng)的異常監(jiān)測預(yù)警方法[3-6]。中石油也開展了特定場景的預(yù)警技術(shù)研究，包括采用故障分類模型開展注水泵異常預(yù)警、時間序列預(yù)測模型進(jìn)行溢罐與剩余應(yīng)急時間預(yù)警、二元分類算法進(jìn)行脫硫塔發(fā)泡預(yù)警[7-8]。但這些研究主要聚焦某一特定場景進(jìn)行建模研究，缺乏整體系統(tǒng)的建模思路研究，或僅從工藝與設(shè)備的角度進(jìn)行簡單的劃分后分別選擇不同的建模方法，沒有闡明每種方法的適用條件和原因。本文結(jié)合天然氣凈化廠工藝和風(fēng)險特點，進(jìn)行系統(tǒng)的闡述和分析，形成完整的建模方法體系，為天然氣凈化廠全流程工藝風(fēng)險預(yù)警建模方法的選擇提供指導(dǎo)，也可為類似的石化裝置建模提供思路和參考。

1 天然氣凈化廠工藝概況

1.1 天然氣凈化廠工藝流程

天然氣凈化廠主體工藝裝置包括過濾分離裝置、脫硫裝置、脫水裝置、硫磺回收裝置、尾氣處理裝置[9]，總體工藝過程如圖1 所示。原料天然氣首先通過過濾分離裝置脫除原料天然氣中的游離水和大部分機(jī)械雜質(zhì)，再由脫硫過程脫除原料氣中的硫化物和部分二氧化碳，常用化學(xué)吸收劑甲基二乙醇胺(MDEA) 脫除，然后通過脫水裝置中的三甘醇溶液脫除氣體中的水分，減少管道和設(shè)備腐蝕，防止水合物的形成。凈化過程產(chǎn)生的酸性氣體由硫磺回收裝置進(jìn)行處理和回收。最后還需尾氣處理裝置來凈化硫磺回收裝置排出的尾氣，保證污染物含量達(dá)到規(guī)定排放要求。

圖1 常規(guī)天然氣凈化廠工藝簡圖

1.2 天然氣凈化廠工藝特征分析

結(jié)合天然氣凈化廠工藝，分析得到天然氣凈化廠生產(chǎn)運(yùn)行過程主要有以下特點[10-11]：

(1)原料氣具有易燃、易爆、有毒、腐蝕性。在生產(chǎn)運(yùn)行過程中，若發(fā)生設(shè)備異常或操作失誤等問題，輕者導(dǎo)致產(chǎn)品不合格、裝置停機(jī)，嚴(yán)重時可能發(fā)生火災(zāi)爆炸、人員中毒、傷亡等事故。

(2)凈化各個工藝過程均涉及高溫、高壓過程。若監(jiān)測及防范不到位，極易發(fā)生火災(zāi)、爆炸事故。

(3)凈化工藝過程復(fù)雜，涉及的設(shè)備、工藝流程、工藝參數(shù)、原料和產(chǎn)品多樣，整體耦合性強(qiáng)，一旦某處發(fā)生異常，容易導(dǎo)致連鎖效應(yīng)，引發(fā)重大安全風(fēng)險。

因此，有必要探索裝置異常發(fā)生前的特征和規(guī)律，建立系統(tǒng)性的預(yù)警建模思路，以便開展凈化廠全工藝流程的異常預(yù)警，實現(xiàn)動態(tài)的風(fēng)險感知與預(yù)測，降低事故發(fā)生可能性，提升天然氣凈化廠生產(chǎn)運(yùn)行過程風(fēng)險管控水平。

2 模型構(gòu)建方法分類

依據(jù)對建模場景的機(jī)理知識和數(shù)據(jù)的需求程度，常用的建模手段分別為機(jī)理建模、專家經(jīng)驗建模、數(shù)據(jù)驅(qū)動建模和混合建模方法，下面簡述這些建模方法的特點。

機(jī)理建模[12-13]是利用能量、質(zhì)量、動量守恒等各種物理化學(xué)原理來描述所研究的工藝場景的規(guī)律，模型的參數(shù)均具備明確的物理含義，最終形成參數(shù)的代數(shù)方程、微分方程、狀態(tài)方程等，從而明確參數(shù)之間的函數(shù)關(guān)系，過程如圖2 所示。通過這些數(shù)學(xué)表達(dá)式，可以精準(zhǔn)把控工藝流程的物理化學(xué)特性與變化規(guī)律，并且能夠預(yù)測各個參數(shù)之間的變化關(guān)系。機(jī)理建模是最傳統(tǒng)的一種建模方法，在理論研究或者工程應(yīng)用方面均有廣泛應(yīng)用。

圖2 機(jī)理模型

專家知識模型是充分獲取長期從事生產(chǎn)運(yùn)行實踐的專家積累的經(jīng)驗知識，將其轉(zhuǎn)化成數(shù)學(xué)規(guī)則，建立規(guī)則模型，可大幅降低建模難度和模型計算量，過程如圖3 所示。專家知識模型適用于能獲得有大量生產(chǎn)經(jīng)驗和工藝知識且在不斷實踐過程中被檢驗經(jīng)驗可靠的專家知識[14]。

圖3 專家知識模型

隨著工業(yè)信息化的發(fā)展，工廠的DCS/SCADA 等系統(tǒng)不斷完善，產(chǎn)生了大量的生產(chǎn)運(yùn)行數(shù)據(jù)。利用海量生產(chǎn)數(shù)據(jù)，融合統(tǒng)計學(xué)和機(jī)器學(xué)習(xí)等手段對數(shù)據(jù)進(jìn)行利用，挖掘數(shù)據(jù)之間的關(guān)系，可建立數(shù)學(xué)驅(qū)動模型[15]，過程如圖4 所示。數(shù)據(jù)驅(qū)動建模無需掌握所研究的工藝場景的內(nèi)部機(jī)理知識或參數(shù)間精確的數(shù)學(xué)表達(dá)式，重點是獲得研究對象的輸入輸出數(shù)據(jù)，進(jìn)行擬合，該方法建模周期短、求解模型難度小、非線性擬合能力強(qiáng)。因此在解決機(jī)理、經(jīng)驗知識不能解決的高階、非線性問題方面有較大的優(yōu)勢。數(shù)據(jù)驅(qū)動建模方法包括傳統(tǒng)的統(tǒng)計學(xué)方法和機(jī)器學(xué)習(xí)方法。

圖4 數(shù)據(jù)驅(qū)動模型

通過分析上述建模方法的特征，總結(jié)以上3 種建模方法的優(yōu)缺點，如表1 所示。

表1 常規(guī)建模方法對比

由表1 可知，三種方法也有各自的局限性，因此混合建模方法被提出。混合建模是將機(jī)理、專家經(jīng)驗和數(shù)據(jù)驅(qū)動方法結(jié)合，利用先驗知識建立機(jī)理或?qū)＜医?jīng)驗?zāi)Ｐ停儆脭?shù)據(jù)驅(qū)動建模方法表達(dá)未知或難以描述的部分，集三者所長，共同完成模型構(gòu)建，從而降低建模難度和過擬合風(fēng)險，提高求解速度，并賦予參數(shù)物理意義，提高可解釋性和外推泛化能力[16]。

3 天然氣凈化廠生產(chǎn)運(yùn)行風(fēng)險預(yù)警建模

根據(jù)天然氣凈化廠工藝概況，得到天然氣凈化廠生產(chǎn)運(yùn)行預(yù)警建模面臨的主要問題為：涉及多個工藝流程，物料數(shù)量、反應(yīng)復(fù)雜程度差異大；現(xiàn)場運(yùn)行專家經(jīng)驗水平參差不齊；通過DSC/SCADA 系統(tǒng)能獲得大量的正常生產(chǎn)運(yùn)行數(shù)據(jù)，但缺乏異常工況樣本。

從工程實際出發(fā)，預(yù)警建模方法應(yīng)按圖5 所示進(jìn)行選擇。根據(jù)本文第2 節(jié)可知，數(shù)據(jù)驅(qū)動建模方法作為黑箱模型，是快速簡單的模型構(gòu)建方法，只要獲取海量的數(shù)據(jù)樣本，就可對任何復(fù)雜過程建模。但由于化工過程監(jiān)管嚴(yán)格，工藝系統(tǒng)多數(shù)都在正常狀況下運(yùn)行，因此異常樣本量不足，導(dǎo)致數(shù)據(jù)驅(qū)動模型的準(zhǔn)確性和外推泛化能力不足，因此，需要在生產(chǎn)運(yùn)行過程中不斷的積累樣本，最終實現(xiàn)采用數(shù)據(jù)驅(qū)動建模方法構(gòu)建風(fēng)險預(yù)警模型這一終極目標(biāo)。目前，受制于樣本量，預(yù)警建模仍以過渡階段和主流方法為主，即采用機(jī)理建模、專家知識建模或者混合建模方法。

圖5 預(yù)警建模方法

當(dāng)然，機(jī)理建模、專家知識建模和混合建模方法的選擇也需結(jié)合風(fēng)險場景。本文結(jié)合西南油田某天然氣凈化廠的實際，開展凈化廠生產(chǎn)預(yù)警方法的深入探討和系統(tǒng)性分析。從保證模型構(gòu)建的準(zhǔn)確性與適用性的角度，將凈化廠預(yù)警場景進(jìn)行歸類，分別進(jìn)行預(yù)警建模方法的闡述并給出詳細(xì)示例，形成適用于凈化廠預(yù)警建模的方法體系，為快速開展凈化廠預(yù)警模型構(gòu)建提供指導(dǎo)，提高建模效率，最終實現(xiàn)風(fēng)險預(yù)警技術(shù)在天然氣凈化廠大幅推廣應(yīng)用，達(dá)到整體提升凈化廠的安全管理水平的目標(biāo)。

3.1 機(jī)理模型應(yīng)用場景—以儲罐溢罐為例

儲罐作為常用的存儲設(shè)備，主要承擔(dān)存儲介質(zhì)的作用，在生產(chǎn)運(yùn)行過程中，僅存在介質(zhì)的流入、流出，無其他物理化學(xué)反應(yīng)。儲罐常見的風(fēng)險為溢罐，因為進(jìn)出口流量不匹配，導(dǎo)致液位持續(xù)上升，而此時液位監(jiān)測儀表出現(xiàn)故障未及時報警處理，從而發(fā)生溢罐。將液位作為特征參數(shù)來表征儲罐溢罐，基于儲油罐D(zhuǎn)CS 系統(tǒng)實時監(jiān)測數(shù)據(jù)，利用流體力學(xué)理論和體積平衡原理，結(jié)合進(jìn)出口流量和液位變化速率，構(gòu)建機(jī)理模型，實時分析判斷儲罐是否發(fā)生溢罐事件。

預(yù)警模型如下：

式(1)～式(2)中：H為實時監(jiān)測液位值；H’為估算理論液位值；Hy為理論值與實際液位值的允許偏差；Q1為DCS 系統(tǒng)中進(jìn)罐實時監(jiān)測體積流量值；Q2為出罐實時監(jiān)測體積流量值；t為分析判斷周期；S為罐體的截面積。

設(shè)儲油罐進(jìn)口流量誤差為q1，出口流量誤差為q2，液位計誤差為h，理論液位值與實際監(jiān)測液位的差值應(yīng)將儀表誤差考慮在內(nèi)，因此Hy應(yīng)大于遠(yuǎn)傳儀表誤差帶來的液位變化。

預(yù)警的目的是提前發(fā)現(xiàn)風(fēng)險，及時處置。因此，儲罐溢罐的分析判斷周期t應(yīng)小于現(xiàn)場操作人員的巡檢周期t0。

此外，由于儀表即液位計和流量計本身存在誤差，分析判斷周期也應(yīng)考慮儀表誤差。

根據(jù)式(1)～式(4)，得到：

通過構(gòu)建上述的機(jī)理模型，當(dāng)實際測量值與預(yù)估值的差值超過了儀表誤差時，可進(jìn)行預(yù)警，在儲罐發(fā)生泄漏的初期及時發(fā)現(xiàn)以將事件的風(fēng)險降到最低。

表2 為模型驗證分析結(jié)果，結(jié)合現(xiàn)場液位計偏差和流量計偏差情況將閾值Hy設(shè)為5 mm，偏差小于5 mm，則認(rèn)為儲罐未泄漏；若偏差值大于5 mm，則認(rèn)為存在泄漏風(fēng)險。從表2 中可知，此時儲罐處于正常運(yùn)行范圍內(nèi)，無溢罐風(fēng)險。

表2 驗證分析結(jié)果表

表3 脫硫塔監(jiān)測數(shù)據(jù)

對諸如儲罐溢罐這類場景，主要特點有：(1)涉及物理化學(xué)反應(yīng)少甚至沒有理化反應(yīng)；(2) 僅有一種狀態(tài)的工藝介質(zhì)；(3)一個特征參數(shù)就可以表征其風(fēng)險。對于此類場景，則可直接采用機(jī)理建模，只需引入簡單的機(jī)理方程，如質(zhì)量守恒、能量守恒方程等，就可求得特征參數(shù)。在此類場景下，構(gòu)建的機(jī)理模型不僅計算量小、準(zhǔn)確度高、可解釋性強(qiáng)，還方便模型復(fù)用，使預(yù)警模型可以大范圍推廣應(yīng)用。

3.2 專家經(jīng)驗?zāi)Ｐ蛻?yīng)用場景—以脫硫塔發(fā)泡為例

脫硫系統(tǒng)是天然氣凈化廠的核心工藝，其中脫硫塔為最核心設(shè)備之一。在生產(chǎn)運(yùn)行過程中，脫硫塔異常發(fā)泡現(xiàn)象常發(fā)。根據(jù)脫硫塔的工藝可知，脫硫塔反應(yīng)過程不僅包括含硫天然氣、甲基二乙醇胺(MDEA)溶液兩種主要的混合介質(zhì)，還可能包含未除干凈的固體雜質(zhì)。脫硫塔設(shè)備內(nèi)部結(jié)構(gòu)復(fù)雜，包括多層塔板、溢流管、除沫器和溢流堰等。脫硫塔異常發(fā)泡表現(xiàn)為脫硫塔中的漿液發(fā)生溢流。受制于脫硫塔復(fù)雜的設(shè)備結(jié)構(gòu)和化學(xué)反應(yīng)，難以采用機(jī)理模型對發(fā)泡過程進(jìn)行表征。但因為脫硫塔發(fā)泡現(xiàn)象比較常見，使現(xiàn)場的操作人員積累了較成熟的經(jīng)驗知識。當(dāng)脫硫塔將要發(fā)生異常發(fā)泡時，脫硫塔差壓會偏離穩(wěn)定值并逐漸增大，而脫硫塔液位會偏離穩(wěn)定狀態(tài)并逐漸下降。因此，考慮根據(jù)獲取的專家經(jīng)驗構(gòu)建專家經(jīng)驗?zāi)Ｐ汀?/p>

選用脫硫塔差壓和液位這兩個特征參數(shù)來表征脫硫塔發(fā)泡。本例采用模糊評價方法將評價指標(biāo)數(shù)據(jù)化，并依據(jù)不同指標(biāo)對評價對象的影響程度來分配權(quán)重，實現(xiàn)脫硫塔發(fā)泡過程的量化分析，構(gòu)建脫硫塔異常發(fā)泡預(yù)警模型[17]。

根據(jù)特征參數(shù)對發(fā)泡的影響程度，確定偏差程度和權(quán)重。采用升半矩形分布和降半矩形分布的隸屬度函數(shù)，將工況參數(shù)進(jìn)行歸一化處理。脫硫塔“發(fā)泡”時差壓參數(shù)會偏離穩(wěn)定值表現(xiàn)為增加趨勢，設(shè)為偏大型因素，采用升半矩形分布，液位被稱為偏小型因素，采用降半矩形分布。

計算偏大型因素的偏差程度的公式為：

式中：x為該參數(shù)的實測值；[xmin,xmax]為該參數(shù)的正常范圍；xH為報警高限值。

計算其偏小型因素的偏差程度的公式為：

式中：x為該參數(shù)的實測值；[xmin,xmax]為該參數(shù)的正常范圍；xL為報警低限值；xH為報警高限值。生產(chǎn)現(xiàn)場《化工裝置操作工藝卡》中明確規(guī)定了xH和xL的設(shè)定值。

液位權(quán)重計算公式如下：

差壓權(quán)重計算公式：

特征參數(shù)權(quán)重之和為1。

式(9)～式(11)中：vA為液位參數(shù)的變化速率；vB為差壓參數(shù)的變化速率；εA為液位參數(shù)的權(quán)重；εB為差壓參數(shù)的權(quán)重。

依據(jù)特征參數(shù)偏差程度計算，得到特征參數(shù)的偏差程度矩陣μi：

依據(jù)特征參數(shù)在設(shè)備異常狀態(tài)下的變化速率，得到特征參數(shù)的權(quán)重向量εi：

根據(jù)特征參數(shù)的偏差程度矩陣μi以及特征參數(shù)間的權(quán)值向量εi，計算表征脫硫塔異常“發(fā)泡”狀態(tài)預(yù)警模型為：

式中：i為某一種異常狀態(tài)；c1、c2、c3分別對應(yīng)優(yōu)、中、差3 個等級；βi為表征某一種異常狀態(tài)的量化指標(biāo)，被稱為異常狀態(tài)置信度。

將預(yù)警模型應(yīng)用于某天然氣凈化廠，收集異常發(fā)泡時差壓和液位的監(jiān)測數(shù)據(jù)，計算得到液位參數(shù)權(quán)重εA為0.27，差壓參數(shù)權(quán)重εB為0.73。根據(jù)先驗知識，設(shè)備異常置信度β＞0.1 時即可確定脫硫塔發(fā)生異常“發(fā)泡”。

對比DSC 系統(tǒng)提醒和預(yù)警提醒的時間如表4 所示，發(fā)現(xiàn)預(yù)警時間比DCS 系統(tǒng)報警時間提前8 min，可達(dá)到過程安全管理關(guān)口前移的目的。

表4 各策略異常提醒區(qū)別

因此，對類似脫硫塔異常發(fā)泡的預(yù)警場景，主要特點有：(1) 涉及物理化學(xué)反應(yīng)；(2) 包含多種狀態(tài)的工藝介質(zhì)；(3) 風(fēng)險的設(shè)備載體結(jié)構(gòu)復(fù)雜；(4) 發(fā)生頻率高，風(fēng)險發(fā)生后，容易提煉出工藝參數(shù)的變化特征，可依賴專家經(jīng)驗；(5) 需要至少兩個及以上的特征參數(shù)來表征其風(fēng)險。對此類場景，優(yōu)先考慮獲取專家經(jīng)驗，將專家經(jīng)驗轉(zhuǎn)化成特征參數(shù)及參數(shù)間的關(guān)系規(guī)則，構(gòu)建經(jīng)驗?zāi)Ｐ停缺ＷC了模型的可靠性，又減少了建模的工作量。

3.3 混合建模應(yīng)用場景—以設(shè)備腐蝕預(yù)警為例

天然氣凈化過程中的容器、設(shè)備易因腐蝕導(dǎo)致器壁減薄，嚴(yán)重時會穿孔或失效，引發(fā)火災(zāi)爆炸風(fēng)險。但腐蝕機(jī)理、原因復(fù)雜多樣，受設(shè)備材質(zhì)、生產(chǎn)運(yùn)行環(huán)境、設(shè)備存儲或輸送的介質(zhì)類型及工藝參數(shù)變化的影響，難以構(gòu)建機(jī)理模型。另一方面，設(shè)備腐蝕是一個長期的過程，且在腐蝕初期，對設(shè)備的壓力、溫度影響不大，只能通過外觀觀測和設(shè)備壁厚監(jiān)測來識別，即使是長期從事生產(chǎn)運(yùn)行的專家也難以對腐蝕特征進(jìn)行可靠的預(yù)測，難以獲取可靠的專家經(jīng)驗。因此，考慮構(gòu)建混合模型來開展設(shè)備腐蝕預(yù)警研究。

在設(shè)備生產(chǎn)運(yùn)行過程中，會定期對設(shè)備的壁厚進(jìn)行檢測，積累了一定的壁厚檢測數(shù)據(jù)。可通過這些數(shù)據(jù)，構(gòu)建壁厚與時間變化的數(shù)據(jù)驅(qū)動模型，分析設(shè)備壁厚值隨時間的變化規(guī)律。然后結(jié)合設(shè)備承壓能力計算公式，推測設(shè)備承壓能力的變化情況，當(dāng)實際操作壓力大于承壓能力時，及時發(fā)出預(yù)警信息。通過構(gòu)建壁厚預(yù)測和承壓能力計算的混合模型，實現(xiàn)設(shè)備的腐蝕預(yù)警。

本文獲取了某凈化廠原料過濾器某點處的一組實際的壁厚測量值，如表5 所示，以此為例進(jìn)行腐蝕預(yù)警研究。本方法是基于壁厚測量數(shù)據(jù)，未對腐蝕的各種影響因素展開分析。因此，該方法僅適用于均勻腐蝕的壽命預(yù)測，通常在容器各部位腐蝕狀態(tài)差別不大的情況下選用。預(yù)測實際上就是用過去和現(xiàn)在的數(shù)據(jù)作為樣本集，構(gòu)建預(yù)測模型，對變量的未來行為狀態(tài)做出科學(xué)定量的估計。

表5 原料過濾器壁厚檢測數(shù)據(jù)

本文選用回歸方法來構(gòu)建設(shè)備壁厚預(yù)測模型，分別使用一元線性回歸和指數(shù)回歸方法進(jìn)行建模[18]。

一元線性回歸法的基本方程為：

式中：a、b分別為回歸系數(shù)；x為自變量時間；y為因變量預(yù)測值。

指數(shù)回歸方法的基本方程為：

使用上述方法擬合模擬設(shè)備的厚度變化趨勢并與預(yù)測結(jié)果進(jìn)行對比驗證，選擇預(yù)測誤差最小的方法預(yù)測設(shè)備的厚度。然后再結(jié)合設(shè)備承壓能力計算公式，進(jìn)行設(shè)備的腐蝕預(yù)警。設(shè)備承壓能力計算式如下：

式中：pc為計算壓力(MPa)；δ為壁厚(mm)；φ為焊接接頭系數(shù)；[σ]t為設(shè)計溫度下圓筒或球殼的計算壓力(MPa)；Di為筒體內(nèi)徑(mm)；P為該設(shè)備的操作壓力。

當(dāng)時pc＜P，進(jìn)行預(yù)警，則給出降壓運(yùn)行或增加設(shè)備強(qiáng)度提醒。

取表5 中前6 組數(shù)據(jù)為計算數(shù)據(jù)，第7 組數(shù)據(jù)為驗證數(shù)據(jù)。

一元線性回歸法得到的預(yù)測方程為：

指數(shù)回歸法得到的厚度預(yù)測方程為：

根據(jù)直線擬合厚度預(yù)測方程，可以得到2014 年4 月25 日的容器厚度應(yīng)為14.195 mm；根據(jù)指數(shù)擬合得到的2014 年4 月25 日的預(yù)測值14.143 mm。顯然，直線擬合的誤差較低，應(yīng)該采用直線擬合來預(yù)測設(shè)備的厚度。根據(jù)設(shè)備壁厚預(yù)測模型，再結(jié)合承壓能力機(jī)理模型，實時計算設(shè)備的承壓能力，并預(yù)測設(shè)備的失效時間，實現(xiàn)設(shè)備腐蝕風(fēng)險預(yù)警。

因此，對類似設(shè)備腐蝕這類場景，主要特點有：(1) 涉及復(fù)雜的物理化學(xué)反應(yīng)；(2) 受環(huán)境及工藝介質(zhì)類型和參數(shù)變化影響大，影響風(fēng)險發(fā)生的因素多，難以通過長期的運(yùn)行操作經(jīng)驗獲取風(fēng)險運(yùn)行規(guī)律；(3) 需要至少兩個及以上的特征參數(shù)來表征風(fēng)險。對于此類場景，可考慮可建立混合模型。例如使用機(jī)理或經(jīng)驗完成基礎(chǔ)模型搭建，數(shù)據(jù)驅(qū)動完成輔助參數(shù)計算、優(yōu)化等工作，綜合機(jī)理/ 專家經(jīng)驗和數(shù)據(jù)驅(qū)動的優(yōu)勢，克服單個建模方法的不足，建立準(zhǔn)確率高、預(yù)警效果佳的模型。

4 結(jié)語

開展生產(chǎn)運(yùn)行風(fēng)險預(yù)警研究是提高企業(yè)安全管理水平的重要手段，天然氣凈化廠也不例外。本文結(jié)合天然氣凈化廠工藝，形成覆蓋凈化工藝全流程的工藝風(fēng)險預(yù)警建模方法體系。未來可參考本文推薦的建模思路，開展具體風(fēng)險場景的預(yù)警模型構(gòu)建與應(yīng)用，以降低生產(chǎn)運(yùn)行風(fēng)險，提高凈化廠本質(zhì)安全水平。此外，還可將本文的思路應(yīng)用于其他石化工藝，開展生產(chǎn)運(yùn)行預(yù)警模型構(gòu)建。當(dāng)然，在涉及數(shù)據(jù)驅(qū)動建模時，仍舊存在著樣本數(shù)據(jù)不足，模型的準(zhǔn)確率具有不確定性的問題。因此，應(yīng)通過持續(xù)的預(yù)警建模研究不斷積累樣本，并結(jié)合機(jī)器學(xué)習(xí)技術(shù)的不斷發(fā)展，探索優(yōu)質(zhì)高效的數(shù)據(jù)驅(qū)動建模方法。