石油罐區(qū)應急操作仿真系統(tǒng)設計

許 祺 夏 濤

（北京化工大學信息科學與技術學院）

由于石油化工安全事故頻發(fā)，因此石化企業(yè)的安全生產越來越重要。 “十二五”規(guī)劃明確提出要求企業(yè)員工提高應急處置能力。 早在20世紀80年代，西方國家就用仿真培訓系統(tǒng)來訓練工人并取得了卓越的效果。

傳統(tǒng)化工仿真［1］側重于“工藝流程建?！?，具體更偏向于職業(yè)學校學生和新員工熟悉工藝流程和DCS操作，缺少關于應急救援的培訓。國內也有專注于石化生產安全的建模仿真平臺［2］，在建模過程中對設備生產安全屬性進行抽取，對事故擴散演變也進行了基于概率的建模。 然而缺乏災害機理模型，無災害救援描述，多不支持三維展示，培訓方面價值有限。 國內目前建立的基于3D或VR的應急救援演練系統(tǒng)［3］，幾乎都以腳本式操作為主，更大程度上是將應急預案內容在三維環(huán)境里進行展示，太流于表面化，對災害模擬、演練推演等涉及較少，沒有脫離游戲的本質。 而多數(shù)基于事故機理模型的災害仿真平臺［4］多用于事故分析，不適合應急演練，在訓練層面上表現(xiàn)得不夠好。 因此，研究并實現(xiàn)VR與災害仿真模型有機結合的應急演練平臺是國內化工安全領域迫切的需求。

筆者設計的應急操作仿真系統(tǒng)， 按照機理化、全息化、空間化和多樣化的原則，在化工仿真的基礎上加入災害建模，引入虛擬現(xiàn)實地理信息系統(tǒng) （Virtual Reality Geographic Information System，VRGIS），支持多客戶端、內外操協(xié)同演練，仿真場景多變，應急演練更真實有效。

1 事故仿真支持

為了更好地支持事故機理建模與更多要素，在傳統(tǒng)建模的基礎上，做了一些災害仿真與應急演練方面的研究，其主要內容包括以下3個部分。

1.1 災害場概念

災害場［5］其實就是能量場，多個相同類型的事故會疊加形成一個能量場，不同類型的事故會形成不同的能量場。 筆者設計的仿真系統(tǒng)支持3種災害場：溫度場、氣體場和爆炸超壓。 溫度場和氣體場對周邊模型的影響是持續(xù)的，而爆炸超壓對周邊模型的影響是一次性的。 火災事故會形成一個溫度場，爆炸事故會形成沖擊波，毒氣泄漏事故會形成氣團。 所有設備或多或少都會受到災害場的影響。

1.2 三維模型

近年來三維仿真技術逐漸成熟，國內大部分化工廠都有自己的三維模型，但大部分工廠的三維模型并沒有被很好地利用起來。 虛擬現(xiàn)實比動畫更加強調參演人員的互動性，所以它在應急演練中占了很重要的比重。 三維場景給仿真帶來的不止是視覺上的沖擊，在應急演練過程中，三維引擎中的碰撞檢測扮演了很重要的角色。

1.3 事故多米諾理論

事故多米諾理論 （Accident Domino Theory）是海因里希于1929年提出的一種解釋和分析事故成因的理論。 初始事故發(fā)生后進一步擴展影響到相鄰設備設施， 從而持續(xù)引發(fā)連鎖事故發(fā)生，造成更為嚴重的事故后果［6］。為了輔助研究，筆者所設計的仿真系統(tǒng)也支持事故多米諾模型。

2 隱患模型建模方法

化工設備是化工廠中實現(xiàn)化工生產所采用的必需工具，長期處于擱置易燃易爆、有毒有害化學物品的狀態(tài)。 設備缺陷是指設備存在威脅設備和人身安全、影響經濟運行、降低設備“健康”水平、影響設備性能或壽命等的隱患，或已經造成化工產品質量不合格等的異?，F(xiàn)象。

化工設備仿真建模是構建化工仿真系統(tǒng)的基石， 針對化工裝置中出現(xiàn)的每一個實體設備，都應該有一個該設備的仿真模型與之對應。 傳統(tǒng)仿真模型可復用性和擴展性都比較差，筆者設計的仿真模型因采用面向對象的思想，可復用性很高， 并在傳統(tǒng)建模的基礎上加入了多種要素，例如設備自身“健康”值、空間坐標及天氣風向等。模型的多樣性使得應急演練仿真更加接近真實場景，增加了仿真的要素，仿真效果更真實。

2.1 隱患模型構成

設備隱患產生的原因是多方面的，有設計缺陷、管理不當、維修不當及工藝參數(shù)變化等。 為了響應設備的各種隱患，需要在設備上綁定一個或多個災害模型，組成隱患設備，即“工藝模型+災害模型=隱患模型”。 工藝模型與災害模型單獨建模。 以閥門模型來舉例說明，閥門可能會產生靜電，也可能會發(fā)生泄漏，泄漏后形成的液池可能會經由靜電點火燃燒產生池火（圖1）。 工藝模型建模已經有比較成熟的機理模型，能夠為系統(tǒng)的實現(xiàn)提供科學的理論依據(jù)。

圖1 閥門隱患模型構成

2.2 災害模型建模

以池火災為例，可燃液體泄漏后流到地面形成液池， 遇到火源燃燒或達到燃點而形成池火。經查閱文獻，Heskestad G經過廣泛的實驗數(shù)據(jù)進行數(shù)學處理［7］，得到了一系列池火以及傷害的相關計算公式，包括液體泄漏速率、燃燒速度、火焰高度、熱輻射通量及目標入射熱輻射強度等。

液體泄漏速率公式由伯努利方程導出：

式中 A——泄漏裂口面積，m2；

Cd——泄漏系數(shù)；

g——重力加速度，m/s2；

h——裂口處液位高度，m；

p——容器內壓力，Pa；

p0——環(huán)境壓力，一般取大氣壓，Pa；

QL——液體泄漏速率，kg/s；

ρ——液體密度，kg/m3。

當液池可燃物的沸點高于周圍環(huán)境溫度時，液池表面單位面積燃燒速率為：

式中 Cp——液體的定壓比熱容，J/（kg·K）；

dm/dt——表面單位面積燃燒速率，kg/（m2·s）；

H——液體蒸發(fā)熱，J/kg；

HC——液體燃燒熱，J/kg；

T0——環(huán)境溫度，K；

Tb——液體沸點，K。

火焰高度計算公式為：

式中 Hf——火焰高度，m；

r——液池半徑，m；

ρ0——周圍空氣密度，一般取2.93kg/m3。

液池燃燒時放出的總熱輻射通量為：

式中 Q——總熱輻射通量，W；

η——效率因子，可取0.13～0.35。

假設全部輻射熱量由液池中心點的小球面輻射出來，則在距離池中心某一距離處的入射熱輻射強度為：

式中 I——熱輻射強度，W/m2；

tc——熱傳導系數(shù)， 在無相對理想的數(shù)據(jù)時，可取1；

X——目標點到液池中心的距離，m。

熱輻射是引起二次事故的重要原因，Cozzani V等提出了熱輻射導致設備損壞的概率模型［8］：

式中 P——設備損壞概率；

x——概率單位值的微元；

Y——概率單位變量。

2.3 應急模型建模

系統(tǒng)中還有一類應急模型， 它包含滅火器、消防車、防護服及人模型等。 應急模型的作用是讓仿真演練更為真實，操作性更強。 應急模塊也是安全機理模型，需通過查閱文獻、使用仿真軟件等方式讓模塊建模更加合理準確。

以泡沫滅火器［9］為例，儲罐火災泡沫滅火劑用量公式為：

式中 q——混合液供給強度，L/（min·m2）；

Qf——泡沫液需求量，L；

S——燃燒液面面積，m2；

t——泡沫連續(xù)供給時間，min；

α——泡沫液的混合比。

以人為例建立人受熱輻射影響的模型，根據(jù)式（5） 可計算出人受到的熱輻射強度。 Acciarri M等提出了人在未裝備防火服等防火裝備時暴露于火焰熱輻射下的死亡概率公式［10，11］，該公式較為復雜，在仿真中，為保證實時性，應事先將人員受熱輻射影響的后果計算出來。 根據(jù)該公式，再基于合理假設，可確定人員在熱輻射場中多久會死亡，假定人員總生命值為100點，表1為人員受熱輻射傷害情況。

表1 人員受熱輻射傷害

3 系統(tǒng)實現(xiàn)及開發(fā)

系統(tǒng)提供仿DCS頁面、3D頁面和頭戴式VR仿真，多樣且直觀的表現(xiàn)方式能夠滿足多工種員工訓練的需求。 系統(tǒng)具有強大的可擴展性和復用性，除了仿真教學與培訓外，還可用于事故演變、事故驗證等的研究。 同時系統(tǒng)還可提供強大的接口， 可以作為數(shù)據(jù)中心與第三方平臺進行通信，為其他化工生產的研究提供實時仿真數(shù)據(jù)。

3.1 系統(tǒng)功能模塊

系統(tǒng)功能模塊如圖2所示，包括工程師站、教師站、學員站和系統(tǒng)管理。

圖2 系統(tǒng)功能模塊框圖

工程師站是由專業(yè)的化工工程師，通過簡單的系統(tǒng)語言進行模型構建，并通過系統(tǒng)平臺進行裝置的組裝調試。 工程師站的主要功能圍繞建模展開，目前支持設備模型導入、圖形組模、工段建模及工廠建模等功能，并增加了數(shù)據(jù)查看、趨勢顯示及指定設備運行等輔助調試功能，方便工程師快速完成建模工作。

學員站用來提供給客戶進行裝置的培訓與考核， 培訓方式有DCS操作、 三維與DCS結合操作。 目前可以同時接入兩種不同的學員操作站：仿DCS操作站和3D操作站。 這兩種操作站分別對應內操與外操的訓練，數(shù)據(jù)互聯(lián)互通，在同一個工藝模型里工作， 實現(xiàn)了內操與外操的協(xié)同訓練。

3.2 仿真實例

3.2.1 案例對象工藝概況

2002年1月7日某縣一加油站發(fā)生爆燃事故，加油站內一臺加油機被燒得面目全非， 加油站屋頂被掀開，一側的圍墻也被炸倒，造成1人死亡，2人受傷。據(jù)調查，事故原因為油罐車司機駕駛油罐車為加油站送油，未按規(guī)程采用密封式輸油法，而直接將輸油管插入儲油罐中。 加油站一員工前去關閥門時，所穿的衣服產生靜電引起爆燃。

以此事故為例， 對加油站進行建模仿真，并說明系統(tǒng)如何對事故演變進行仿真，同時給出了進行應急演練操作的步驟。 表2和圖3為加油站的基本情況，站內存有5塊滅火氈，一臺25kg推車式泡沫滅火器。

表2 儲罐結構參數(shù)

3.2.2 異常與事故處理模擬

基于第2節(jié)提出的隱患設備建模方法對加油站進行建模。

仿真流程如下：油罐車進入加油站時，未做除靜電處理。 罐車?？吭?2#汽油罐V01右側，在卸油途中，操作工未采用密封式輸油法，致使油管意外脫落并發(fā)生油品泄漏。 此時因靜電引燃液池，形成池火。

圖3 加油站平面布置圖

已知汽油泄漏點位于西側靠墻區(qū)域，距東側的加油槍15m；西側因為有墻體阻擋，汽油只能往其他方向漫延。 風速為3m/s，風向為西，加油機整體位于火災的下風向。泄漏口直徑為8cm；汽油液位最高點距泄漏口垂直高度為2m； 汽油密度為800kg/m3； 汽油的單位面積燃燒速率為0.055kg/（m2·s）。 忽略罐內外空氣壓差，根據(jù)液體泄漏速率計算公式計算可得： 液位高度在1.0m時，泄漏速率為10.4kg/s；當液位高度在0.1m時，泄漏速率為3.2kg/s。

分兩種情況進行仿真：泄漏發(fā)生時，罐車汽油液位距泄漏口高度分別為1.0m和0.1m。

情況一：第270s發(fā)生泄漏，此時罐車汽油液位距泄漏口高度1.0m，剩余汽油超過10t，泄漏量較大，液池面積因持續(xù)泄漏擴大，第295s發(fā)生火災。因站內只有小型推車式滅火器，操作工判斷無法救援，迅速逃離并報警。 火勢蔓延，引燃并燒毀加油機，若操作工不及時逃離也會受到傷害。 圖4為泄漏量與池火等效半徑趨勢。

圖4 泄漏量與池火等效半徑趨勢

情況二：第110s發(fā)生泄漏，此時罐車汽油液位距泄漏口高度0.1m，剩余汽油質量約為1t，第145s時液池被靜電火花點著，發(fā)生火災。 在第880s時，操作工找到滅火器后開始對池火進行撲救并成功撲滅。 因罐車內已無汽油，不存在復燃可能性，此次應急救援結束。 圖5為對池火進行撲救的三維仿真畫面，圖6為此次滅火救援關鍵參數(shù)趨勢。

圖5 對池火進行撲救三維仿真畫面

圖6 滅火救援關鍵參數(shù)趨勢

3.2.3 結果

通過對加油站罐區(qū)卸油事故仿真驗證了所設計的仿真系統(tǒng)。 仿真系統(tǒng)以安全機理模型為基礎，在仿真系統(tǒng)的統(tǒng)一調度下，模型自行判斷引入災害模型， 所有模型都會對災害場作出反應，進而出現(xiàn)事故多米諾現(xiàn)象。 用戶可以在仿DCS畫面和三維場景中對流程中出現(xiàn)的事故災害進行協(xié)同處理救援。

4 結束語

基于隱患模型建模方法，搭建了罐區(qū)應急操作仿真平臺，將數(shù)據(jù)渲染到虛擬現(xiàn)實場景中。 目前實現(xiàn)的仿真場景，涵蓋了常見的幾類化工災害事故，實現(xiàn)了災害仿真與應急演練相結合，不僅可用于員工培訓，還可用于事故分析以及輔助制定應急預案等。 現(xiàn)在災害模型計算并不很精確，在以后的工作中，需要收集更多的真實事故數(shù)據(jù)來繼續(xù)矯正災害模型，災害對于災害場的疊加作用也要更精細。