王麗苑，杜曉亮，王興峰，吳世洋，李衛(wèi)東，王鑫羽

（西北礦冶研究院環(huán)境資源研究所，甘肅 白銀730900）

1 引言

環(huán)境風(fēng)險是由自發(fā)的自然原因和人類活動引起的，通過環(huán)境介質(zhì)傳播的，能對人類社會及自然環(huán)境產(chǎn)生破壞、損害及至毀滅性作用等不幸事件發(fā)生的概率及其后果［1，2］。環(huán)境風(fēng)險評價是分析和預(yù)測建設(shè)項目存在的潛在危險、有害因素，項目建設(shè)和運行期間可能發(fā)生的突發(fā)性事件和事故，引起有毒有害和易燃易爆等物質(zhì)泄漏，所造成的人身安全與環(huán)境影響及損害程度，提出合理可行的防范、應(yīng)急與減緩措施，以使建設(shè)項目事故率、損失和環(huán)境影響達到可接受的水平［2］。為了達到環(huán)境風(fēng)險評價的目的，環(huán)境風(fēng)險預(yù)案的計算與后果預(yù)測就成為環(huán)境風(fēng)險評價中的重中之重。

2 硫酸儲罐發(fā)生泄漏的可能性分析

硫酸貯罐發(fā)生泄漏的事故樹分析見圖1。硫酸儲罐泄漏事故樹定量分析結(jié)果見表1。

事故樹結(jié)構(gòu)函數(shù)為：T＝（q1＋q2＋q3）q4＋q5＋q6q7q8＝q1q4＋q2q4＋q3q4＋q5＋q6q7q8；

事故樹有5個最小割集：｛｛q1q4｝、｛q2q4｝、｛q3q4｝、｛q5｝、｛q6q7q8｝｝；

表1 硫酸儲罐泄漏事故樹定量計算結(jié)果

事故概率的估算：

圖1 硫酸儲罐泄漏事故樹

通過對以上事故樹的計算與分析得知，硫酸貯罐發(fā)生泄漏時的危險重要系數(shù)順序為：q（7）＝q（8）＞q（6）＝q（4）＝q（5）＞q（1）＞q（3）＝q（2）。

事件名稱是：自動控制系統(tǒng)失靈＝未及時關(guān)閉閥門＞生產(chǎn)報警系統(tǒng)事故＝毒氣檢測事故＝外觀檢測不合格＞閥門故障＞管道破裂＝貯罐破裂。

造成硫酸貯罐發(fā)生泄漏的因素很多，其中基本事件q（7）與q（8）的危險度最大，是造成硫酸貯罐發(fā)生泄漏的直接原因；其次q（6）、q（4）、q（5）是結(jié)構(gòu)危險度相等的3個原因，再下來是q（1）、q（3）、q（2）等3個原因。由此可以看出，要控制基本事件——硫酸貯罐泄漏的發(fā)生，必須從最危險的事件q（7）與q（8）入手。因此，要控制硫酸貯罐泄漏事故的發(fā)生，在設(shè)計階段就要按標(biāo)準(zhǔn)的要求對各個元件進行選型與設(shè)計，以確保達到安全的本質(zhì)；其次在生產(chǎn)過程中要經(jīng)常檢查，以確保各類元件處于正常的狀態(tài)，要經(jīng)常對職工進行安全教育，使職工處于一個高度的認(rèn)識狀態(tài)，以確保安全地運行。

由事故樹的定量分析可知，儲罐泄漏引發(fā)火災(zāi)爆炸事故的發(fā)生概率約為10－6。此類事故雖然發(fā)生概率較低，但危害后果較嚴(yán)重。

3 項目事故源強確定

本項目中硫酸泄漏的主要原因是儲運設(shè)施缺乏維護，造成罐體或管道開裂引起硫酸泄漏；裂口面積0.01m2，儲罐泄漏后，安全系統(tǒng)報警，操作人員在20min內(nèi)使儲罐泄漏得到制止，并在泄漏物料上方噴灑泡沫，覆蓋泄漏物料阻止泄漏液體的揮發(fā)，同時采取有效的收集措施，在20min內(nèi)將泄漏物料收集到備用儲罐［3］。

3.1 泄漏量計算

硫酸泄漏速率計算公式［4］可采用下式：

式中：QL為液體泄漏速度，kg/s；Cd為液體泄漏系數(shù)，此值常用0.6～0.64，取值為0.62；A為裂口面積，m2，取值0.2×0.01＝0.002m2；P為容器內(nèi)介質(zhì)壓力，178200Pa；P0為環(huán)境壓力，78520Pa；g為重力加速度，9.8m/s2；h為 裂 口 之 上 液 位 高 度，9m；ρ為 密 度，取1.84g/mL。

由計算可知，硫酸泄漏速度為0.752kg/s，20min泄漏量約為902.4kg。以濃硫酸密度為1.84t/m3，硫酸在地面形成的液體厚度0.005m推算，泄露的硫酸在地面形成的面積為735m2［5］。

3.2 泄漏液體蒸發(fā)量計算

泄漏液體的蒸發(fā)分為閃蒸蒸發(fā)、熱量蒸發(fā)和質(zhì)量蒸發(fā)［6］3種，由于濃硫酸常壓下的沸點為330℃，而項目儲罐儲存溫度和環(huán)境溫度均不高于40℃，當(dāng)液體泄漏時不發(fā)生閃蒸和熱量蒸發(fā)，因此不考慮閃蒸蒸發(fā)量和熱量蒸發(fā)量。

質(zhì)量蒸發(fā)速度Q3按下式計算［7］：

式中：Q3為質(zhì)量蒸發(fā)速度，kg/s；a，n為大氣穩(wěn)定度系數(shù)（液池蒸發(fā)模式參數(shù)見表2）；p為液體表面蒸氣壓，Pa（取值1300Pa）；R為氣體常數(shù)，J/mol·K（取值為8.31）；T0為環(huán)境溫度，K（按282.8K計算）；u為風(fēng)速，m/s（取值1/2.2/5m/s）；r為液池半徑，m（經(jīng)計算本項目罐區(qū)液池等效半徑為15.3m）；M為液體摩爾質(zhì)量，kg/mol（取值0.098kg/mol）。

表2 液池蒸發(fā)模式參數(shù)

表3 不同條件下的排放速率 g/s

根據(jù)公式（2）計算出項目不同風(fēng)速、不同穩(wěn)定度下儲罐泄漏時溶劑油的排放速率見表3。

4 硫酸泄漏事故的后果計算

4.1 后果計算預(yù)測模式

有毒有害氣體事故排放時間短，并且具有煙團排放的特點，故硫酸泄漏采用《建設(shè)項目環(huán)境風(fēng)險評價技術(shù)導(dǎo)則》（HJ/T169－2004）中的多煙團擴散模式進行預(yù)測，計算事故后釋放的氣體擴散。在事故后果評價中采用下列煙團公式［8］：

式中：C（x，y，0）為下風(fēng)向地面坐標(biāo)處的空氣中污染物濃度（mg/m3）；x0，y0，z0為煙團中心坐標(biāo)；Q為事故期間煙團的排放量；σX、σy、σz為X、Y、Z方向的擴散參數(shù)（m），常取σX＝σy。

4.2 硫酸霧連續(xù)泄漏擴散濃度分布及危害距離預(yù)測計算與評價

硫酸泄漏預(yù)測時間為20min，第20min時刻，下風(fēng)向硫酸霧擴散濃度分布，預(yù)測對環(huán)境的影響情況。

根據(jù)查閱相關(guān)資料［9～11］，風(fēng)險評價指標(biāo)體系見表4。

表4 風(fēng)險評價指標(biāo)系數(shù)一覽 mg/m3

4.3 預(yù)測結(jié)果

為了說明不同氣象條件下溶劑油儲罐對周圍空氣環(huán)境的影響情況，分別選取小風(fēng)（1.0m/s）、年均風(fēng)速（1.92m/s）及大風(fēng)（5m/s）3種風(fēng)速和B、D、E這3種穩(wěn)定度下，預(yù)測發(fā)煙硫酸泄漏后硫酸霧下風(fēng)向地面濃度，結(jié)果見表5及圖2～圖4。

表5 硫酸儲罐泄露時下風(fēng)向地面濃度結(jié)果一覽

圖2 風(fēng)速1m/s時，不同大氣穩(wěn)定度條件下硫酸泄漏時濃度變化曲線

圖3 風(fēng)速1.92m/s時，不同大氣穩(wěn)定度條件下硫酸泄漏時濃度變化曲線

圖4 風(fēng)速5m/s時，不同大氣穩(wěn)定度條件下硫酸泄漏時濃度變化曲線

由表5的預(yù)測結(jié)果可知，在硫酸儲罐泄漏事故發(fā)生20min時，靜風(fēng)條件，E穩(wěn)定度下，下風(fēng)向最大落地濃度110.12g/m3，出現(xiàn)距離1.6m處，居住區(qū)最高容許濃度范圍183.2m；年均風(fēng)速1.92m/s條件，E穩(wěn)定度下，下風(fēng)向最大落地濃度1870.47g/m3，出現(xiàn)距離為2.8m處，居住區(qū)最高容許濃度范圍317.9m；大風(fēng)條件，E穩(wěn)定度條件下，最大落地濃度703.15g/m3，出現(xiàn)在7.3m處，居住區(qū)最高容許濃度范圍317.9m。如圖5所示，挑選年平均風(fēng)速時的預(yù)測圖形進行分析。

如圖5所示在年均風(fēng)速條件下影響范圍較小，對附近800m范圍內(nèi)居民和建筑物基本無影響，大于800m后有輕微影響。

項目敏感點居民區(qū)距離廠區(qū)最近距離1000m左右，位于廠址上風(fēng)向，基本不會受到大的影響，硫酸泄漏只是短時的，通過大氣的擴散和稀釋，影響會逐漸消失。

圖5 1.92m/s時，E穩(wěn)定度下污染物擴散范圍

5 結(jié)論與建議

化工行業(yè)環(huán)境風(fēng)險評價一般應(yīng)包括風(fēng)險識別、源項分析、后果計算、風(fēng)險防范措施及應(yīng)急預(yù)案等主要內(nèi)容，其中事故發(fā)生頻率及事故源強的確定是風(fēng)險評價的難點和重點。根據(jù)本文的分析與計算，本項目結(jié)論如下。

（1）本項目事故源強為20min的泄漏量902.4kg；

（2）多煙團擴散模式適用于硫酸泄漏時的后果計算與預(yù)測；

（3）應(yīng)用環(huán)境風(fēng)險評價軟件計算年平均風(fēng)速下E穩(wěn)定度下硫酸儲罐泄漏影響范圍為1273m。

［1］曹希壽.區(qū)域環(huán)境風(fēng)險評價與管理初探［J］.中國環(huán)境科學(xué)，1994，14（6）：456～470.

［2］陸雍森.環(huán)境評價［M］.2版.上海：同濟大學(xué)出版社，1999：531～558.

［3］邵超峰.突發(fā)性大氣污染事件的環(huán)境分析按評估與管理［J］.環(huán)境科學(xué)與技術(shù)，2009，32（6）：200～204.

［4］Kamal M.AI－Subhi AI－Harbi.Application of the AHP in project Management［J］.Project Management，2001（19）：19～27.

［5］王 超.液體泄漏形成液池擴散面積的計算方法綜述［J］.安全與環(huán)境工程，2012，19（6）：125～128.

［6］D Colbournea，K O Suenb.Appraising the flammability hazardsof hydrocarbon refrigerants using quantitative risk assessmentmodel［J］.International Journal of Refrigeration，2004，33（1）：774～783.

［7］潘旭海，蔣軍成.事故泄漏源模型研究與分析［J］.南京工業(yè)大學(xué)學(xué)報：自然科學(xué)版，2002，24（1）：105～110.

［8］劉 蕾，彭 量.多煙團模式在環(huán)境風(fēng)險評價中的應(yīng)用［J］.氣象科技，2012，38（6）：684～688.

［9］Asian Development Bank.Environmentla Risk Assessment［R］.Office of the Environment，ADB，1990.

［10］付 鐵.石化項目大氣環(huán)境風(fēng)險評價研究［D］.蘭州：蘭州大學(xué)，2006.

［11］林 莉，孔慶波.淺析化工企業(yè)環(huán)境風(fēng)險評價［J］.環(huán)境科學(xué)與管理，2008，33（4）：78～82.