重大危險源分級標準

1、v1.0可編輯可修改重大危險源分級標準（征求意見稿）1適用范圍本規(guī)范規(guī)定了重大危險源評估分級的方法和程序。本規(guī)范為重大危險源評估分級技術規(guī)范， 適用于包括儲罐區(qū)、庫區(qū)、生產(chǎn)場所等重大危險源。2規(guī)范性引用文件下列文件中的條款，通過本規(guī)范的引用而成為本標準的條款。凡是標注日期的引用文件，其隨后所有的修改單（不包括勘誤的內容）或修訂版 均不適用于本規(guī)范，然而，鼓勵根據(jù)本規(guī)范達成協(xié)議的各方研究是否可使 用這些文件的最新版本。凡是不注日期的引用文件，其最新版本適用于本 規(guī)范。中華人民共和國安全生產(chǎn)法危險化學品安全管理條例安全生產(chǎn)許可證條例重大危險源辨識（GB18218安全評價通則關于規(guī)范重大危險源監(jiān)督與

2、管理工作的通知（安監(jiān)總協(xié)調字2005125號）3術語和定義卜列術語和定義適用于本規(guī)范重大危險源 major hazard installations重大危險源是指長期地或者臨時地生產(chǎn)、搬運、使用或者儲存危險物品，且 危險物品的數(shù)量等于或超過臨界量的單元（包括場所和設施）。4重大危險源分級判據(jù)重大危險源分級判據(jù)如表1所示。表1重大危險源分級判據(jù)危險源等級一級重大危險源二級重大危險源三級重大危險源四級重大危險源具體判別的依據(jù)如下:一級重大危險源：可能造成死亡二級重大危險源：可能造成死亡三級重大危險源：可能造成死亡四級重大危險源：可能造成死亡分級判據(jù)死亡人數(shù)可能造成30人（含30人）以上可能造成10

3、 29人可能造成39人可能造成1-2人30人（含30人）以上的重大危險源;10-29人的重大危險源；3-9人的重大危險源；1-2人的重大危險源。5重大危險源死亡人數(shù)及財產(chǎn)損失計算方法可能造成的死亡人數(shù)評價程序為：將重大危險源的周邊區(qū)域劃分成等間隔的網(wǎng)格區(qū)，用一笛卡爾坐標體系的 網(wǎng)格覆蓋城市的區(qū)域地圖（如圖1所示），網(wǎng)格間距大小取決于當?shù)厝丝诿芏龋?以 不影響計算結果為準。確定每一網(wǎng)格內的人員數(shù)量，通過火災（室內火災除外）、爆炸、毒物泄漏 擴散事故后果模型計算重大危險源事故在每一網(wǎng)格中心處產(chǎn)生的熱輻射、超壓或毒物濃度的數(shù)值，然后通過熱輻射、沖擊波超壓、中毒概率函數(shù)將其其轉化為造 成死亡的概率。將

4、每一網(wǎng)格中心的死亡率與人口數(shù)量相乘，即得到死亡的人數(shù)。將所有網(wǎng)格的死亡人數(shù)求和，即得到總的死亡人數(shù)。具體用下式表?。?1)nDi S Vi i 1式中，N為總的死亡人數(shù)，D為第i個網(wǎng)格的人口密度，S為網(wǎng)格面積，Vi 為第i個網(wǎng)格的個人死亡率，n為網(wǎng)格的數(shù)目。FeHRlzen industry: Petroleum refinery iFMflfe 1: 十皿PeSroleunS re-finiryiRfBanffukhadlM扁：Y,二 We| lendb-wDbmt& 后二 dChemical industry圖1死亡人數(shù)計算原理示意圖采用財產(chǎn)損失半徑的方法評估事故后果造成的損失，

5、并假定此半徑內沒有損失的財產(chǎn)與此半徑外損失的財產(chǎn)相互抵消， 或者說此半徑內的財產(chǎn)完全損失，此半徑外的財產(chǎn)完全無損失。財產(chǎn)損失半徑通過火災、爆炸事故后果模型確定。財產(chǎn)損失半徑按下式計算：19Ri，1/3K iWTNT217631751WTNT(2)式中，R為i區(qū)半徑，mi K為常量熱輻射對建筑物的影響直接取決于熱輻射強度的大小及作用時間的長短，以引燃木材的熱通量作為對建筑物破壞財產(chǎn)損失半徑，計算公式如下:q 6730t 4/5 25400(3)t W/Mc(4)式中，q為引燃木材的熱通量（W/m）, t為熱輻射作用時間，即火災持續(xù)時問（s）。6重大危險源評價分級程序重大危險源的評價分級程序如下圖

6、所示。如果一種危險物質具有多種事故形 態(tài)，按照后果最嚴重的事故形態(tài)考慮，即遵循“最大危險原則” 。各類重大危險 源具體事故情景選擇、后果計算及死亡概率計算過程參見附錄 A。圖2重大危險源評價分級程序附錄A:重大危險源事故后果模型儲罐區(qū)重大事故后果分析A.1.1儲罐區(qū)的主要事故后果類型A.1.1.1 池火災易燃液體如汽油、苯、甲醇、乙酸乙酯等，一旦從儲罐及管路中泄漏到地面 后，將向四周流淌、擴展，形成一定厚度的液池，若受到防火堤、隔堤的阻擋， 液體將在限定區(qū)域（相當于圍堰）內得以積聚，形成一定范圍的液池。這時，若 遇到火源，液池可能被點燃，發(fā)生地面池火災。A.1.1.2蒸氣云爆炸易燃易爆氣體如H

7、2、天然氣等，泄漏后隨著風向擴散，與周圍空氣混合成易 燃易爆混合物，在擴散擴過程中如遇到點火源，延遲點火，由于存在某些特殊原 因和條件，火焰加速傳播，產(chǎn)生爆炸沖擊波超壓，發(fā)生蒸氣云爆炸。易燃易爆的液化氣體如液化石油氣、 液化丙烷、液化丁烷等，其沸點遠小于 環(huán)境溫度，泄漏后將會由于自身的熱量、地面?zhèn)鳠?、太陽輻射、氣流運動等迅速 蒸發(fā)，在液池上面形成蒸氣云，與周圍空氣混合成易燃易爆混合物，并且隨著風向擴散，擴散擴過程中如遇到點火源，也會發(fā)生蒸氣云爆炸。A.1.1.3噴射火對于易燃易爆氣體如H2、天然氣，以及易燃易爆的液化氣體來說，泄漏后可 能因摩擦產(chǎn)生的靜電立即點火，產(chǎn)生噴射火。A.1.1.4沸騰

8、液體擴展蒸氣云爆炸易燃易爆的液化氣體容器在外部火焰的烘烤下可能發(fā)生突然破裂，壓力平衡被破壞，液體急劇氣化，并隨即被火焰點燃而發(fā)生爆炸，產(chǎn)生巨大的火球。這種 事故被稱為沸騰液體擴展為蒸氣云爆炸。A.1.1.5中毒事故毒性的液化氣體如液氯、液氨等，由于沸點小于環(huán)境溫度，泄漏后會因自身 熱量、地面?zhèn)鳠?、太陽輻射、氣流運動等迅速蒸發(fā)，生成有毒蒸氣云，密集在泄 漏源周圍，隨后由于環(huán)境溫度、地形、風力和湍流等因素影響產(chǎn)生漂移、擴散， 范圍變大，濃度減小。A.1.2儲罐區(qū)主要事故后果模型A.1.2.1池火災事故后果模型池火災火焰的幾何尺寸及輻射參數(shù)按如下步驟計算。計算池直徑根據(jù)泄漏的液體量和地面性質，按下式

9、可計算最大可能的池面積。S亞此訪（1）式中，S為液池面積（m）,W為泄漏液體的質量（kg）,為液體的密度（kg/m3） 鼻皿為最小油層厚度（m。最小物料層厚度與地面性質對應關系見表 1。表1不同性質地面物料層厚度表地面性質最小物料層厚度（成草地粗糙地面平整地面混凝土地面平靜的水面確定火焰高度計算池火焰高度的經(jīng)驗公式如下：h LD 42 mf/（ 0'gD）0.61（2）式中：L為火焰高度（m, D為池直徑（m）, m為燃燒速率（kg/m2s）, P 0為空氣密度（kg/m3）, g為引力常數(shù)。計算火焰表面熱通量假定能量由圓柱形火焰?zhèn)让婧晚敳肯蛑車鶆蜉椛?，用下式計算火焰表面的熱通?q

10、(o0.25 D2 HCmf f0.25 D 2 DL(3)式中，q。為火焰表面的熱通量（kw/m2）, A H為燃燒熱（kJ/kg ）,冗為圓周 率，f為熱輻射系數(shù)（可取為）,m為燃燒速率（kg/m2s）,其它符號同前目標接收到的熱通量的計算目標接收到的熱通量q（r）的計算公式為：q（r） q°（1 0.058lnr）V（4）式中，q（r）為目標接收到的熱通量（kw/m2）, q0為由式（3）計算的火焰表 面的熱通量（kw/m2）, r為目標到油區(qū)中心的水平距離（V為視角系數(shù)。視角系數(shù)的計算角系數(shù)V與目標到火焰垂直軸的距離與火焰半徑之比s,火焰高度與直徑之比h有關。V,(Vv2 V

11、H2)VHA BA (b 1/s), 1tan0.5/(b21)0.50.5B (a1/s), 1tan/(a2 1)0.5Vvtan 1 (h /(s21)0.5)/s h(JK)/s0.5K tan 1s 1/s1 0.5(5)(6)(7)(8)(9)(10)(11)22-a (h s 1)/(2s)(b (1 s2)/(2s)(13)其中A、R J、K、VH、VV是為了描述方便而引入的中間變量，冗為圓周率。A.1.2.2蒸氣云爆炸事故后果模型蒸氣云爆炸產(chǎn)生的沖擊波超壓是其主要危害。沖擊波超壓可通過傳統(tǒng)的TNT當量系數(shù)法進行計算，將事故爆炸產(chǎn)生的爆炸能量等同于一定當量的TNT也可根據(jù)爆炸能

12、量直接計算。（1） TNT當量法確定閃蒸系數(shù)在熱力學數(shù)據(jù)資料的基礎上，用下式估算燃料的閃蒸部分。Cp TF 1 exp L(14)式中，F(xiàn)為蒸發(fā)系數(shù)，C為燃料的平均比熱（kJ/kgK） , AT為環(huán)境壓力下容器內溫度與沸點的溫差（K） , L為汽化熱（kJ/kg） o計算云團中燃料的質量:Wf 2FW(15)式中，W為云團中燃料的質量（kg） , W為泄漏的燃料的質量（kg） , F為閃蒸系數(shù)。計算TNT當量WTNTeWf H f / H tnt（16）式中，WNt為燃料的TNT 當量 （kg） ,W為云團中燃料的質量（kg）, H為燃料的燃燒熱（MJ/kg） , Hnt為1Q1J1O2無因次

13、距離 mkg -1/3. Wtnt口 1 無因次超壓IX圖1 A Ps R曲線TNT的爆熱(MJ/kg) , a e為TNT當量系數(shù)，推薦 y將實際距離轉化為無因次距離:R R/Wtnt(17)式中，R為離爆炸點的實際距離(m), R為無因次距離(m)。在離爆炸點距離為R處，根據(jù)相應的R值，查圖1得到超壓，進而預測人員 受傷害和建筑受破壞的情況。(2)直接計算法在得到云團中燃料的質量的情況下，可按下式直接計算爆炸沖擊波超壓A p23ln( ps/pa)0.9126 1.5058(lnZ) 0.1675(ln Z) 0.0320(ln Z) (18)(WZ0 12)z R/(馬 13(19)Pa

14、E 1.8 WQC(20)式中，A ps為沖擊波正相最大超壓(Pa), Z為無量綱距離，Pa為環(huán)境壓力， R為目標到爆源的水平距離(項，E為爆源總能量(J), a為蒸氣云當量系數(shù)， 一般取，W為蒸氣云中對爆炸沖擊波有實際貢獻的燃料質量(Kg), Q為燃料的燃燒熱(J/Kg)。A.1.2.3噴射火事故后果模型加壓的可燃物泄漏時形成射流，如果在泄漏裂口處被點燃，則形成噴射火。假定火焰為圓錐形，并用從泄漏處到火焰長度4/5處的點源模型來表示?；鹧骈L度計算噴射火的火焰長度可用如下方程得到：0 444 L (Hem).161.66(21)式中，L為火焰長度（m, HC為燃燒熱（J/kg ）, m為質量流

15、速（kg/s ）熱輻射的通量計算距離火焰點源為X（m）處接收到的熱輻射通量可用下式表示:fH Cmq24 X2 1000(22)式中，q為距離X處接收的熱輻射的通量（KW/市，f為熱輻射率，為大 氣傳輸率。大氣傳輸率p按下式計算：1 0.0565 In X(23)A.1.2.4沸騰液體擴展為蒸氣云爆炸事故后果模型計算主要包括如下步驟?；鹎蛑睆紻 2.665W0.327(24)式中，D為火球直徑（m1 , W為火球中消耗的可燃物質量（Kg）。對單罐儲存, W取罐容量的50%對雙罐儲存，W取罐容量的70%對多罐儲存，W取罐容量的90%火球持續(xù)時間t 1.089W0.327（25）式中，t為火球持續(xù)

16、時間（s）, W同式（24）?；鹎蛱叨然鹎蛟谌紵龝r，將抬升到一定高度?；鹎蛑行木嚯x地面的高度H由下式估計：H D（26）火球表面熱輻射能量假設火球表面熱輻射能量是均勻擴散的?；鹎虮砻鏌彷椛淠芰縎EP由下式計_2SEP FsmHa /( D t)(27)式中，F(xiàn)s為火球表面輻射的能量比，H為火球的有效燃燒熱(J/Kg)。Fs與儲罐破裂瞬間儲存物料的飽和蒸氣壓力P (MPa有關：Fs 0.27P0.32(28)對于因外部火災引起的BLEVE*故，上式中的P值可取儲罐安全閥啟動壓力 Pv (MPa的倍，即：P 1.21FV(29)H由下式求得：Ha Hc Hv CpT(30)式中，H為燃燒熱(J

17、/kg), T為常沸點下的蒸發(fā)熱(J/kg), Cp為恒壓比熱 (J/), T為火球表面火焰溫度與環(huán)境溫度之差(K), 一般來說T=1700K視角系數(shù)視角系數(shù)F的計算公式如下：2F (D/2)/r)(31)式中，r為目標到火球中心的距離(m。令目標與儲罐白水平距離為 x(m),則：r (X2 H 2)0.5(32)大氣熱傳遞系數(shù)火球表面輻射的熱能在大氣中傳輸時，由于空氣的吸收及散射作用，一部分 能量損失掉了。假定能量損失比為a,則大氣熱傳遞系數(shù)a 1 o a和大氣中 的CO和H2O的含量、熱傳輸距離及輻射光譜的特性等因素有關。a可由以下的經(jīng)驗公式來求?。篴 2.02(pwr) 0.09(33)

18、式中，國為環(huán)境溫度下空氣中的水蒸氣壓(N/m), r為目標到火球表面的距離(mj)。Pw pW RH(34)0式中，Pw為環(huán)境溫度下的飽和水蒸氣壓(N/m), RH為相對濕度。r' r D/2(35)火球熱輻射強度分布函數(shù)在不考慮障礙物對火球熱輻射產(chǎn)生阻擋作用的條件下，距離儲罐X處的熱輻射強度q (W/m)可由下式計算：q SEP F a(36)A.1.2.5中毒事故后果模型(1)泄漏模型液體泄漏速率模型液體泄漏可根據(jù)流體力學中的柏努力方程計算泄漏量。 當裂口不規(guī)則時，可 采取等效尺寸代替；當泄漏過程中壓力變化時，則往往采用經(jīng)驗公式。柏努力方 程如下：c2(P P。)Q CdA , 2

19、gh'(37)式中，Q為液體泄漏速率(kg/s) , Cd為無量綱泄漏系數(shù)，是液體密度(kg/m 3) , A是泄漏孔面積(m2), P為罐壓(Pa) , R為大氣壓力(Pa) , g為引力常 數(shù)(9.8m/s 2) , h為液壓高度(m)。液體出口速度可按下式計算：Qu Cd A(38)式中，u為液體出口速度(m/s),其他符號如前。持續(xù)時間按下式計算：ts u0/(Cd g)(AT/A)(39)v1.0可編輯可修改式中，U0為初始流速（m/s）, Ar為罐內液面積（m2）。泄漏系數(shù)C的取值通?？蓮臉藴驶瘜W工程手冊中查到。對于管道破裂，Cd的典型取值為。表2為常用的液體泄漏系數(shù)數(shù)據(jù)。

20、表2液體泄漏系數(shù)Cd雷諾數(shù)Re裂 口形狀圓形（多邊形）三角形長方形>100<100這個方法沒有考慮泄漏速率對時間的依賴關系 （壓力隨時間而降低以及液壓 高度下降）。因此，計算出的泄漏速率是保守的最大可能泄漏速率氣體泄漏模型壓力氣體泄漏通常以射流的方式發(fā)生，泄漏的速度與其流動的狀態(tài)有關，其 特征可用臨界流（最大出口速度等于聲速）或亞臨界流來描述。Perry等人用如下的關系式作為臨界流的判斷準則：當式（40）成立時，氣體流動屬音速流動；當式（41）成立時，氣體流動屬亞音速流動25k(40)Po21P k 1k(41)Po2 此P k 1式中，Po為環(huán)境大氣壓力（Pa）, P為容器壓力（

21、Pa）, k為氣體的絕熱指數(shù),即定壓比熱G和定容比熱C之比。對于很多氣體，臨界比值（P/P。）c近似為2,也就是說儲壓近似等于大氣 壓力的兩倍，此時流體泄漏的出口速度近似等于聲速。 臨界流的質量泄漏速率可按下式計算:Q CdAPk 1Mk 2 k 1RT k 1氣體呈亞音速流動時，其泄漏量為：k 1Mk 2 k 1、Q YCdAPj （43） RT k 1式中，Q是氣體泄漏速率（kg/s）, Cd為氣體泄漏稀疏，A為裂口面積（m2）, M是氣體相對分子質量，R是普適氣體常數(shù)（J mol -1K1）, T是氣體的儲存溫度2P k 1P0（K）, Y為氣體膨脹因子，按式（44）計算k 1(44)P

22、okP上述考慮的為理想氣體的不可逆絕熱擴散過程。 止匕外，沒有考慮氣體泄漏速率隨時間的變化，因此使用初始儲存條件必然導致保守的結果。兩相流泄漏模型Cude在1975年建議了兩相流泄漏關系式。假設源容器和泄漏點之間的管道長度和管道直徑之比L/ D>12泄漏點壓力與泄漏點上流壓力之比Pc/P=。具體計算方法如下：第一步，按下式計算兩相流的質量分數(shù)：(45)(T Tc)CpH式中，M為蒸發(fā)的液體占液體總量的比例，Tc是對應于泄漏點壓力 R的平衡 溫度（K）, T是對應于泄漏點上流壓力 P的平衡溫度（K）, G是液體的定壓比熱 J / （kg K），T是液體的蒸發(fā)熱（J/kg）。第二步，按下式計

23、算兩相流的平均密度：(46)(47)1M V 1 MVV1式中，、v和1分別是兩相流、蒸氣和液體的密度（kg/m3）第三步，按下式計算兩相流的質量泄漏速率 Q （kg/S）:Q ACd 1 2 (P Pc)式中，Cd為泄漏系數(shù)，多數(shù)情況下，取 Cd=, A為裂口面積（m2）, P為兩相 混合物的壓力（Pa）, PC為臨界壓力（Pa）。如果L/D<12,先按前面介紹的方法計算純液體泄漏速率和兩相流泄漏速 率，再用內插法加以修正。兩相流實際泄漏速率的計算公式為：Q Qvi （Qi Qvi）（12 L/D）/10（48）式中，Q、Qvi和Qi分別為兩相流實際泄漏速率、按式（43）計算出的兩相

24、流泄漏速率和純液體泄漏速率（kg/S）。如果L/D02, 一般認為泄漏為純液體泄漏。（3）非重氣云擴散模型瞬間泄漏擴散模型2QC(x, y,z,t,He)(2 )3/2 -exp 喂 expexpx y z(z H)22 Zx (x Vt)22y2 2乙y(49)連續(xù)泄漏擴散模型C(x,y,z,t,He)2Q'exo2 V y zexo(zH)22 2exp(z H)22 Z(50)式中，C為氣云中危險物質濃度(kg/mi),He為泄漏源有效高度（源瞬間泄漏量（kg）, Q'為源連續(xù)泄漏速率（kg/s）, V為風速（m）Zs）, t為泄漏后的時間（s）,y和z分別為x、y和z方

25、向的擴散系數(shù)（m）對于連續(xù)泄漏，平均時間取i0min。其中（Tx, （Ty, （Tz與地面的有效粗糙度 有關。地面有效粗糙度長度如下表所示。表3地面有效粗糙度長度表地回類型Z0/m地回類型Z0/m草原、平坦開闊地<農(nóng)作物地區(qū)村落、分散的樹林1分散的高矮建筑物(城市)14密集的高矮建筑物(大城市)4有效粗糙度ZoW 0.1m地區(qū)的擴散參數(shù)按下表選取。表4 Z 0W 0.1m地區(qū)的擴散參數(shù)大氣穩(wěn)定度0- y/m0- z/mA/人-1/2(1+B/人-1/2(1+C/人-1/2(1+/人-1/2(1+D-1/2(1+-1/2(1+E(1+-1/2(1+-1/2F/A -1/2(1+(1+-1/

26、2有效粗糙度Zo> 0.1m的粗糙地形擴散系數(shù)為：yy0 fyzz0 fzfy(Zo) 1a°ZoIn xfz(x, Zo)(bo Colnx)(do eolnx) 1Z0f0 g0式中，(T y0、(Tz0按表4中的數(shù)值取值。其他系數(shù)按表 5取值表5不同大氣穩(wěn)定度下的系數(shù)值穩(wěn)定度ABCDEFv1.0可編輯可修改a0b0Codoeo fogo式（49）和式（5。）中泄漏源有效高度是指泄漏氣體形成的氣云基本上變成 水平狀時氣云中心的離地高度。在大多數(shù)問題中，泄漏源有效高度難以與泄漏源 實際高度相一致。事實上，它等于泄漏源實際高度加泄漏源抬升高度。泄漏源抬升高度可以用下面的公式近似

27、計算：_1,、H Vsd1.5 o.268Fa（Ts Ta）Ts d/V（51）或H 2.4Vsd/V（52）式中：AH是泄漏源抬升高度（m, Vs是氣云出口速度（用7 s）, d是出口直 徑（m）, V是環(huán)境風速（m/s） , pa是環(huán)境大氣壓力（Pa）, Ts是氣云出口溫度（K）, Ta是環(huán)境大氣溫度（K）。計算出泄漏源抬升高度以后，將泄漏源抬升高度與泄漏源實際高度相加就得 到了泄漏源有效高度。（3）重氣云擴散模型常用模型有盒子模型和平板模型兩類。盒子模型用來描述瞬間泄漏形成的重氣云團的運動，平板模型用來描述連續(xù)泄漏形成的重氣云羽的運動。這兩類模型的核心是因空氣進入而引起的氣云質量增加速率

28、方程。盒子模型盒子模型使用如下假設：I、重氣云團為正立的坍塌圓柱體，圓柱體初始高度等于初始半徑的一半II、在重氣云團內部，溫度、密度和危險氣體濃度等參數(shù)均勻分布。III、重氣云團中心的移動速度等于風速。重氣擴散的盒子模型示意圖如下圖所示。/W/Z/A7/» 空氣通過界面進入云困> 云地向外的徑向增動43圖2重氣云團盒子模型坍塌圓柱體的徑向蔓延速度由下式確定：1 / 2，一、Vf dr/dt g( a)/ ah(53)式中，Vf為圓柱體的徑向蔓延速度(m/s) , r為圓柱體半徑(m), h為圓柱體高度(m), t為泄漏后時間(s)等式兩邊同時乘以2r,上式變成下面的形式:(54

29、)2 1/2 a)/ a "dr2/dt 2 g ( a a 1/2由于假設重氣云團和環(huán)境之間沒有熱量交換，重氣云團的浮力將守恒，即:g ( a)/ a V g ( a)/ a V。(55)將上式代入式(53),積分后得到：991 /2rr。2 g ( ° a)/ a V。/ t(56)式中，r。為重氣云團的初始半徑(m), Vo為重氣云團的初始體積(m3),。為重氣云團的初始密度(kg/m 3) 0由于假設重氣云團是圓柱體，初始高度等于初始半徑的一半，因此重氣云團 初始半徑的計算公式為：r。 (2V。/ )1/3(57)隨著空氣的不斷進入，云團的高度和體積也將不斷變化。云

30、團體積隨時間的 變化速率由下式確定：dV/dt ( R2)VT (2 Rh)VE(58)式中，重氣云團體積VR2h , VT和VP分別為空氣從頂部和邊緣進入重氣云團的速率(m/s)。由于重氣云團內部危險氣體質量守恒，因此，在重氣云團擴散過程中，下式存立：_22C/Co Vo/V (h°r°)/(hr)(59)3式中，G和C分別為初始時刻和t時刻重氣女團內部危險物質濃度(kg/m)。任意時刻重氣云團的半徑按式(56)計算。如果知道任意時刻重氣云團高度的計算公式，利用上式就可計算任意時刻重氣云團內部危險物質濃度。但這里不準備采用先推導重氣云團高度的計算公式，然后計算重氣云團體積

31、和危險物質濃度 的方法。而是先采用量綱分析法求重氣云團的體積和濃度，然后利用上式反推重 氣云團的高度。無量綱量V / Vo與x / VO/3之間存在如下函數(shù)關系：V Vo(x/Vo1/3)1.5,x Vo1/3(60)式中，x為下風向距離(m)。它與時間、風速之間的關系為：x Vt(61)將上式代入式(52),得到：CCo(x/Vo1/3) 1.5,x Vo1/3(62)將圓柱形重氣云團的體積V r2h代入式(60),可推導出1/3、1.521 /3h Vo(x/Vo ) /( r ),x Vo(63)隨著空氣的不斷進入，重氣云團的密度將不斷減小，重氣坍塌引起的擴散將 逐步讓位于環(huán)境湍流引起的

32、擴散。目前，判斷重氣坍塌過程終止的準則主要有:I、e準則定義 (p a)/ a。e準則認為，如果e小于或等于某個臨界值(在之 間)，重氣坍塌引起的擴散將讓位于環(huán)境湍流引起的擴散。下面推導轉變發(fā)生的位置，令：E gV( p a)/ a(64)gV將式(60)代入式(64),得到：E g Vo(x/Vo1/3)1.5(65)從上式求出x,得到：x E2/3Vo 1/3(g ) 2/3(66)由于不考慮云團與環(huán)境之間的熱交換，云團浮力守恒，故E=Eo。代入上式得到轉變點對應的下風向距離為： 2/31/32/3Xf Eo Vo(g cr)(67)II、Ri準則對于瞬間泄漏，定義Richardson數(shù)R

33、i (g( p a)/ a)V1/3/V*2。Ri準則 認為，如果Ri小于或等于某個臨界值(在110之間)，重氣坍塌引起的擴散將 讓位于環(huán)境湍流引起的擴散。由于云團內部浮力守恒，因此Eo E gV( a)/ a(68)對上式進行包等變換，得到：( a)/ a Eo/(gV)(69)將式(69)代入Richardson數(shù)的定義式，得到：Ri Eo/(xVo1/3V*2)(70)從上式求出轉變點下風向距離x,得到：xf Eo /(RicrVo1/3V*2)(71)III、Vf準則定義重氣云團徑向蔓延速度 Vf = dr/dt 。Vf準則認為，如果V小于或等于某 個臨界值，重氣坍塌引起的擴散將讓位于

34、環(huán)境湍流引起的擴散。不同研究人員提出的重氣云擴散階段終止時的臨界V值相差很大。例如，van Ulden認為，重氣云擴散階段終止的條件是 Vf = 2V*, Germeles和Drake認 為，重氣云擴散階段終止的條件是 Vf=V, Cox和Carpenter認為，重氣云擴散 階段終止的條件是Vf = d(ry/dt ,h為橫風向擴散系數(shù)(m), Eidsvik (1980)認 為，重氣云擴散階段終止的條件是M = *。這些準則覆蓋范圍很寬，從 9 =到V=V(假設摩+ S速度M=V715)。Germeles和Drake提出的準則太嚴，按照他們 提出的準則，即使存在重氣云擴散階段，這個階段持續(xù)的

35、時間也很短。不過，多 數(shù)研究人員認為，V的臨界值具有與M相同的數(shù)量別。平板模型平板模型使用了如下假設：I 、重氣云羽橫截面為矩形，橫風向半寬為 b(m),垂直方向高度為h(m)。在 泄漏源點，云羽半寬為高度的兩倍，即：bo= 2h。II 、重氣云羽橫截面內，濃度、溫度、密度等參數(shù)均勻分布。III 、重氣云羽的軸向蔓延速度等于風速。在重氣云羽的擴散過程中，橫截面半寬的變化由下式確定：Vdb/dx (gh( a)/ a)1/2(72)由于假設重氣云羽與環(huán)境之間無熱量交換， 重氣云羽的浮力通量在擴散過程 中守恒,即：2gVbh( a)/ a 2gVbo%( o a)/ a(73)將式(73)代入式(

36、72),積分后得到：b bo1 1.5gho( o a)/ a1/2X(Vbo) 12/3(74)由于重氣云羽初始半寬等于初始高度的兩倍，重氣云羽的初始體積通量為:2Vo' 2bohoV boV(75)從上式可以求出重氣云羽的初始半寬：bo 2ho (Vo/V)1/2(76)隨著空氣的進入，不僅重氣云羽的橫風向水平尺寸要增大，重氣云羽的高度也要增加。重氣云羽高度的變化與下風向距離間的關系由下式確定：dh (we/v)dx(77)式中，W為空氣卷吸系數(shù)，且假設空氣卷吸系數(shù)由下式確定： ' We 3.5V* /(11.67 Ri)(78)式中，Ri為當?shù)豏ichardson數(shù)。上式

37、表明，隨著Richardson 數(shù)的增加，空氣卷吸系數(shù)減小。Richardson數(shù)的定義為：Ri g( p a)/ ah/V*'2(79)式中，V*是垂直萬向的特征湍流速度(m/s),由下式確定：V*' 1.3(V*/V)(4/9)(db/dt)2 V21/2(80)式中，V*是摩擦速度。由于x Vt ,因此，db/dt Vdb/dx,結合式(74),得到： _ _，1/21-db/dt (2/3)(Vebo)1.5gho( o a)/ a Web。) (81)1 1.5gho( o a)/ a1/2X(Vebo) 1 1/3求解由式(78)式(81)組成的聯(lián)立方程組，可以求得

38、任意下風向距離重 氣云羽的高度。由于重氣云羽橫截面上危險物質通量守恒，因此2bhVC 2bohoVCo(82)上式兩邊同時除以2bhV,得到重氣云羽中危險物質濃度的計算公式：C bohoCo /(bh)(83)式中，C表示重氣云羽內危險物質濃度(kg/m3),下標0指初始條件無論是重氣云團擴散，還是重氣云羽擴散，一旦滿足前面討論過的轉變條件，重力驅動擴散將轉變?yōu)榄h(huán)境湍流驅動擴散。為了將轉變前后兩個不同的擴散過程有機銜接起來，需要進行虛源計算。所謂虛源，是指位于轉變點上游某處的虛擬泄漏源。虛源計算的目的是確定虛源與轉變點之間的距離。進行虛源計算時應遵循下面的原則：在相同的泄漏和擴散條件(相同源強

39、、相同地形、相同氣象條件等)下，利用重氣云擴散模型對實源泄漏進行擴散分析得到的轉變點所在位置危險物質濃 度等于利用高斯模型對虛源泄漏進行擴散分析得到的轉變點所在位置危險物質 濃度。虛源計算時假設轉變點的下風向擴散系數(shù)x、橫風向擴散系數(shù)y和垂直方向擴散系數(shù)z分別由下面三個公式計算：xy (瞬間泄漏才需要)(84)yb/21/2(85)zh/21/2(86)因此，如果知道擴散系數(shù)與下風向距離的關系，就可以計算出虛源與轉變點 之間的距離。例如，如果y 0.1X,那么，虛源與轉變點之間的距離 Xv 10b/21/2。庫區(qū)重大事故后果分析A.2.1庫區(qū)的主要事故類型根據(jù)儲存場所的不同以及儲存危險品特性的

40、不同，庫區(qū)主要的事故后果類型如下：(1)倉庫中存儲TNT等爆炸性物品時，容易發(fā)生凝聚相含能材料爆炸，產(chǎn) 生非常嚴重的后果。由于爆炸性物品不得露天堆放，只能存儲在倉庫中，因此通 常發(fā)生的是受限空間的爆炸。(2)易燃、易爆的氣體(包括液化氣體)如液化石油氣鋼瓶等在倉庫中存 儲時，發(fā)生泄漏并在擴散過程中遇到點火源，則很容易發(fā)生蒸氣云爆炸事故。(3)有毒氣體(包括液化氣體)如液氯、液氨鋼瓶在倉庫中存儲時，發(fā)生 泄漏并擴散很容易引起中毒事故。(4)易燃液體如苯、甲醇等瓶裝、桶裝的化工原料在倉庫中存儲時，泄漏 后很容易引發(fā)室內池火災。(5)易燃固體、自燃物品、遇濕易燃物品等瓶裝、桶裝、袋裝的物品在倉 庫中

41、存儲時，容易發(fā)生室內固體火災事故。(6)易燃液體的桶裝的化工原料在堆場中存儲時，泄漏后很容易引發(fā)大面 積的池火災。(7)易燃固體物品在堆場中存儲時，容易發(fā)生固體火災事故。A.2.2庫區(qū)的主要事故后果模型A.2.2.1凝聚相含能材料爆炸凝聚相含能材料爆炸能產(chǎn)生多種破壞效應， 如熱輻射、一次破片作用、有毒 氣體產(chǎn)物的致命效應，但破壞力最強，破壞區(qū)域最大的是沖擊波的破壞效應， 因 此，爆炸模型主要考慮沖擊波的傷害作用。凝聚相含能材料的爆炸沖擊波超壓A p可按下式計算：(87)(88)(89)(90)1 0.1567Z 3( p 5)Ps0.137Z 3 0.119Z 2 0.269Z 1 0.019

42、(1 p 10)PsP/P0z R/( )13 P。E 1.8WQC式中，Z為無量綱距離，A p為目標處的超壓值(Pa), p0為環(huán)境壓力，R為 目標到爆源的水平距離(m, E為爆源總能量(J), W為含能材料的質量(Kg), Q為爆炸物的爆熱(J/Kg)。A.2.2.2蒸氣云爆炸見A.1.3節(jié)蒸氣云爆炸事故后果模型。A.2.2.3毒物的泄漏擴散見A.1.3節(jié)中毒事故后果模型。A.2.2.4池火災見A.1.3節(jié)池火災事故后果模型。A.2.2.5固體火災固體火災的熱輻射參數(shù)按點源模型估計。此模型認為火焰射出的能量為燃燒 的一部分，并且輻射強度與目標至火源中心距離的平方成反比。q(r)fMcHc4

43、 r2(91)式中，q（r）為目標接收到的輻射強度（W/rri）, f為輻射系數(shù)，可取f= , M為燃燒速率（kg/s ）, H為燃燒熱（J/kg ）, r為目標至火源中心的水平距離（rr）0生產(chǎn)場所重大事故后果分析A.3.1生產(chǎn)場所主要事故類型生產(chǎn)場所的事故類型非常復雜，因反應介質、工藝設備與機器、操作條件的 不同而不同，常見的危害較大的主要包括以下幾類：（1）爆炸物理爆炸化工容器及設備由于設計、制造、腐蝕或低溫、材料缺陷、交變載荷的作用， 使得器壁的平均應力超過材料的屈服點或強度極限，導致脆性疲勞、疲勞破裂和 應力腐蝕破裂發(fā)生物理爆炸，也可因安全泄放裝置失靈、液化氣體充裝過量、嚴 重受熱膨

44、脹、違章超負荷運行等發(fā)生物理爆炸。常見的如水夾汽包、化工容器、 液化器氣瓶等的爆炸。化學爆炸化工設備和機器內的物質發(fā)生極迅速、劇烈的化學反應而產(chǎn)生高溫高壓可引 發(fā)瞬間的爆炸現(xiàn)象。一般可分為簡單分解爆炸、復雜分解爆炸和爆炸性混合物爆 炸。在化工、石油化工生產(chǎn)中發(fā)生的化學爆炸絕大部分是爆炸性混合物爆炸。例如由于負壓操作、系統(tǒng)用氣、水封不嚴或失效，空氣用入到裝置中，形成爆炸性 混合物，發(fā)生化學爆炸；再如硝化反應過程中，由于溫度控制不良，很容易引起 爆炸。蒸氣云爆炸化工設備和機器由于密封裝置失效、 設備管道腐蝕、磨損或疲勞破裂、斷裂 以及安裝檢修不良、操作失誤等原因，可燃性氣體從化工裝置、設備、管道內

45、泄 漏或噴射，擴散到周圍環(huán)境中，達到爆炸極限，若遇到明火或高溫很就會發(fā)生蒸 氣云爆炸。(2)燃燒當化工設備和機器由于密封裝置失效、 設備管道腐蝕、磨損或疲勞破裂、斷 裂以及安裝檢修不良、操作失誤等原因導致物料泄漏時， 對于易燃液體而言，泄 漏后形成一定范圍的液池，若遇到火源，液池可能被點燃，發(fā)生 地面池火災；對 于易燃氣體而言，泄漏后立即遇到明火或高溫，或因高速摩擦產(chǎn)生靜電而產(chǎn)生 噴 射火，也可延時點火產(chǎn)生 閃火。含有易燃易爆液化氣體的容器在外部火焰的烘烤下可能沸騰液體擴展為蒸氣爆炸，產(chǎn)生巨大的火球。(3)中毒在化工、石化生產(chǎn)中，由于設備密封不嚴、嚴重腐蝕穿孔、疲勞破裂、磨損、 超壓引起的設備

46、與管道突然斷裂、開錯閥門、閥門密封不嚴、水封失效等原因， 很容易造成毒性氣體的泄漏，向周圍環(huán)境擴散，造成人員的中毒事故。主要化工設備和機器的事故類型及損壞尺寸見下表。表6主要化工設備和機器的事故類型及損壞尺寸序號設備類型事故類型損壞尺寸1塔（吸收塔、蒸福塔、萃取塔、干燥塔）物理爆炸、化學爆炸、嚴重泄漏內部爆炸全部泄漏；孔蓋泄漏20%T徑；噴嘴斷裂100%徑；管路破裂20M 100%管徑。2儲槽、壓縮空氣儲罐、緩沖罐物理爆炸、化學爆炸、嚴重泄漏物理爆炸全部泄漏； 噴嘴斷裂100%r徑；管路破裂20% 100%f徑。3中間儲罐物理爆炸、嚴重泄漏物理爆炸全部泄漏；接頭泄漏20% 100%管徑。4換熱

47、器、冷凝器、冷卻器、再沸器嚴重泄漏并可導致爆炸內部爆炸全部泄漏；管路破裂20M 100%管徑。5反應釜、合成塔、流化床、物理爆炸、化學爆炸、嚴重泄漏內部爆炸全部泄漏；孔蓋泄漏20%T徑；噴嘴斷裂100%f徑；管路破裂20M 100%管徑。6管道、法蘭、接頭撓性連接器、過濾器嚴重泄漏法蘭泄漏20%f徑；管道泄漏20% 100%管徑；接頭損壞20M 100%T徑。7閥門、泄放閥嚴重泄漏殼體泄漏20% 100%f徑；蓋子、桿損壞20%世漏。8壓縮機、泵嚴重泄漏殼體泄漏20M 100%!徑；密封蓋損壞20%世漏。9風機負壓空氣進入化學爆炸，泄漏內部爆炸全部泄漏；殼體泄漏20% 100%管徑。10余熱鍋

48、爐、加熱爐爐管泄漏、爐體爆炸內部爆炸全部泄漏；管路破裂20% 100%管徑。11火炬嚴重泄漏接頭泄漏20% 100%f徑。A.3.2生產(chǎn)場所主要事故后果模型A.3.2.1物理爆炸物理爆炸就是物質狀態(tài)參數(shù)（溫度、壓力、體積）迅速發(fā)生變化，在瞬間放 出大量能量并對外做功的現(xiàn)象。物理爆炸的特點是在爆炸現(xiàn)象發(fā)生過程中，導致 爆炸發(fā)生的介質的化學性質不發(fā)生變化，發(fā)生變化的僅是介質的狀態(tài)參數(shù)。物理爆炸如壓力容器破裂時，氣體膨脹所釋放的能量（即爆破能量）不僅與 氣體壓力和容器的體積有關，而且與介質在容器內的物性相態(tài)有關。當壓力容器中介質為壓縮氣體，即以氣態(tài)形式存在而發(fā)生物理爆炸時， 其釋 放的爆破能量為：k

49、 10.1013-Eg -pV1 （-） 103（92）k 1p式中，Eg為氣體的爆破能量（kJ）, p為容器內氣體的絕對壓力（MPa, V 為容器的容積（m3）, k為氣體的絕熱指數(shù)，即氣體的定壓比熱與定容比熱之比。當介質全部為液體時，鑒于通常用液體加壓時所做的功作為常溫液體壓力容器爆炸時釋放的能量，爆破能量計算模型如下：Ei ”VA（93）2式中，Ei為常溫液體壓力容器爆炸時釋放的能量（kJ）, p為液體的絕對壓 力（Pa）, V為容器的體積（mi）, t為液體在壓力p和溫度T下的壓縮系數(shù)（Pa I。而液化氣體一般在容器內以氣液兩態(tài)存在。當容器破裂發(fā)生爆炸時，除了氣體的急劇膨脹做功外，還有

50、過熱液體激烈的蒸發(fā)過程。在大多數(shù)情況下，這類容 器內的飽和液體占有容器介質重量的絕大部分，它的爆破能量比飽和氣體大得 多，一般計算時不考慮氣體膨脹做的功。過熱狀態(tài)下液體在容器破裂時釋放出爆破能量可按下式計算：E (Hi H2) (S1 S2)TiW(94)式中，E為過熱狀態(tài)液體的爆破能量(kJ), H為爆炸前液化液體的始(kJ/kg ), 代為在大氣壓力下飽和液體的始(kJ/kg ) , Si為爆炸前飽和液體的嫡 (kJ/(kg ? C), &為在大氣壓力下飽和液體的嫡(kJ/(kg ? C), Ti為介質在 大氣壓力下的沸點(C) , W為飽和液體的質量(kg)。壓力容器爆破時，爆破

51、能量向外釋放以沖擊波能量、碎片能量和容器殘余變 形能量三種形式表現(xiàn)出來。根據(jù)介紹，后二者所消耗的能量只占總爆破能量的 3 %15%也就是說大部分能量產(chǎn)生空氣沖擊波。A.3.2.2化學爆炸對于簡單分解爆炸、復雜分解爆炸產(chǎn)生的爆炸可用凝聚相含能材料爆炸模型計算見A.2.3節(jié);對于爆炸性混合物氣體爆炸，可用蒸氣云爆炸事故后果模型， 見節(jié)。A.3.2.3池火災見A.1.3節(jié)池火災事故后果模型。A.3.2.4沸騰液體擴展蒸氣云爆炸見A.1.3節(jié)沸騰液體擴展蒸氣云爆炸事故后果模型。A.3.2.5蒸氣云爆炸見A.1.3節(jié)蒸氣云爆炸事故后果模型。A.3.2.6噴射火見A.1.3節(jié)噴射火事故后果模型。A.3.2.7閃火模型閃火是可燃性氣體或蒸氣泄漏到空氣中，與之混合后被點燃而發(fā)生的一種非 爆炸性的燃燒過程。閃火的主要危害來自熱輻射和火焰直接接觸?？扇嘉镌茍F的大小決定了可能造成直接火焰接觸危害的面積，而云團的大小則部分取決于擴散 和泄漏條件。在閃火模型中，假定閃火是一個火焰以恒定速度傳播的過程。閃火模型的幾何關系如圖3所小。圖3閃火模型火焰高度火焰