2.1 過程爆炸事故模式及致因分析 對蒸氣云爆炸和沸騰液體擴展蒸氣爆炸兩種典型的重大爆炸災害事故進行案例統計與分析,根據近40年來國內外已經正式發布的不同的工廠、工藝過程中所發生的重大氣云爆炸事故數據,根據各事故致因的不同,經統計分析,歸納提煉出的UVCE、BLEVE分為以下事故模式,見表2.1,2.2。
表2.1 蒸氣云爆炸事故模式統計分析
Table 2.1 Statistic and analysis of accident modes for UVCE
圖2.1表明,UVCE發生的場所中,煉油廠占比例最大,氣體加工、乙烯、石油化工、LPG罐區發生比率依次減少。
而導致UVCE發生的事故原因也是有多種多樣的,主要包括閥門泄漏、法蘭失效泄漏、管線失效(損壞、破裂、腐蝕)、儲罐失效(破裂、裂縫、腐蝕、超壓、沖擊作用)、閥門開啟、滿裝外溢等因素導致危險物質泄漏,形成氣云被引爆,以上主要致因所導致的UVCE事故模式所占比率,見圖2.2。
圖2.2表明,由于管線損壞導致的UVCE事故所占比率最大,而由于閥門開啟泄漏、滿裝外溢、儲罐失效(裂縫、腐蝕、破裂、超壓)等因素導致的UVCE事故比率依次減少。
從表2.2,可以看出,導致BELVE的主要原因外來火焰包圍、熱輻射沖擊(火災、太陽)、外來撞擊、貯罐破裂、管線泄漏、槽罐火災(罐車出軌等)、貯罐腐蝕泄漏、管線破裂(低溫、撞擊)、閥門泄漏(凍結、失效,無法關閉),引發BLEVE的重要事故致因所占比率,見圖2.3。
從圖2.3可以看出,在導致BLEVE的事故致因中,火焰包圍導致BLEVE的事故模式比率最大,而管線破裂或失效、外來熱腐蝕沖擊、閥門故障或失效泄漏、貯罐腐蝕泄漏、槽罐火災、外來撞擊等因素引發BLEVE的事故模式比率依次減少。以上統計分析是根據近40年來國內外已經正式公布的重大事故案例數據來進行的,沒有正式公布或事故案例數據不完整的不在統計之列。
2.2 氣云爆炸事故機理
氣云爆炸是化工裝置中發生的重大災難性事故,包括蒸氣云爆炸和沸騰液體擴展蒸氣爆炸兩種事故模式,氣云爆炸的共同特點是參與的爆炸物質量一般都較大,導致的危害和損失巨大。
2.2.1 蒸氣云爆炸事故機理 蒸氣云爆炸(UVCE)是由于氣體或易于揮發的液體燃料的大量快速泄漏,與周圍空氣混合形成覆蓋很大的范圍的“預混云”,在某一有限制空間遇點火而導致的爆炸。
UVCE發生有一定的條件,包括一定量的液化石油氣泄漏并與周圍空氣預混、延遲點火、局限化的空間[1](周圍環境如樹木、房屋、設備及其它建筑物等形成具有一定限制的空間)等。蒸氣云爆炸事故在液化氣等物質的運輸過程中并不常見,因為在運輸過程中即使發生液化氣泄漏,也不容易積聚起來。而在液化氣儲備場所,UVCE事故相對易于發生。導致UVCE發生的事故原因也是有多種的,主要包括上文統計分析出的閥門泄漏、法蘭失效泄漏、管線失效(損壞、破裂、腐蝕)、儲罐失效(破裂、裂縫、腐蝕、超壓、沖擊作用)、閥門開啟、滿裝外溢等因素導致危險物質泄漏,形成氣云被引爆。
UVCE的發生大多數是由于儲存液化氣等物質的設備罐體在機械、化學或熱作用下發生破壞而導致大量液化氣泄漏所引起的,此外工作人員在裝運取樣等日常業務中是否正確操作,也是導致罐內液化氣泄漏的一個重要因素。罐體破裂是導致UVCE發生的直接原因,因此研究罐體破裂的原因是研究UVCE事故機理的重點。液化氣容器在受到機械作用(如撞擊、打擊)、化學作用(如腐蝕)或熱作用(如火焰環境、熱沖擊)時,由于所收作用程度的不同,容器會發生以下幾種破裂模式。容器破裂導致蒸氣云爆炸機理,見圖2.4。
(1) 容器罐體突然炸裂,產生巨大的沖擊力,炸裂的碎片以很大的動量向四周拋射出去,比如容器受到物體猛烈的撞擊或震蕩,或遭受劇烈的熱沖擊等作用。
(2) 容器局部破裂,導致液化氣以氣液兩相的混合物從裂口噴出,比如容器在遭受到非均勻熱沖擊、小物體撞擊或化學腐蝕等作用,容器局部器壁材料軟化、失效或腐蝕成孔而導致容器局部破裂。
(3) 安全閥動作及失效,導致容器內壓力過高而爆炸。
蒸氣云爆炸特點。UVCE具有以下特點:一般由火災發展成爆燃,而不是爆轟;蒸氣云的形成是加壓存儲的可燃液體和液化氣大量泄漏的結果,存儲溫度一般大大高于它們的常壓沸點;參與蒸氣云爆炸的可燃氣體或蒸氣的量一般在5×103kg以上;參與蒸氣云爆炸的燃料最常見的是低分子碳氫化合物,偶爾也有其它物質,如氯乙烯、氫氣與異丙醇等;爆源初始尺寸與特征長度相當,并且蒸氣云爆炸的能量釋放速率比較小,是一種面源爆炸。
蒸氣云爆炸傷害形式。UVCE發生后,云霧區內的爆炸波作用、云霧區外的沖擊波作用、高溫燃燒作用和熱輻射作用,以及缺氧造成的窒息作用的主要因素,是造成對周圍人員、建筑物、儲罐等設備的傷害、破壞作用。蒸氣云爆炸的破壞效應表現為:形成相當大的火球,在大氣中形成爆轟波,其強度取決于氣云的燃燒速度;碎片效應通常可以忽略[2]。
蒸氣云爆炸氣云燃燒傳播及超壓形成機理。氣云點燃后的燃燒模式最可能是爆燃,而不是爆轟,爆燃是沿著波的前峰在壓力和密度上都減小的膨脹波,屬于亞音速的。當可燃氣云和空氣的預混物在低能量點火下就會發生爆燃。
2.2.2 沸騰液體擴展蒸氣爆炸事故機理
裝有液化氣的容器當處于火焰環境下、受到撞擊或機械失效等狀態時,容器突然破裂,壓力平衡破壞,LPG急劇氣化,大量的氣化的LPG釋放出來,并隨即被火焰點燃就會導致沸騰液體擴展蒸氣爆炸。根據實驗研究表明,BELVE可以分為熱BLEVE和冷BLEVE[3],根據BLEVE過熱極限理論,熱BLEVE是指BLEVE的發生是由于LPG在爆炸前高于LPG大氣過熱極限溫度,而冷BLEVE是指BLEVE的發生是由于LPG在爆炸前低于LPG大氣過熱極限溫度。冷BLEVE一般是由于強度比較差的容器由于機械或熱的作用而引起的災難性失效所導致的,而熱BLEVE一般是由于過熱液體在容器局部失效時發生噴射釋放而引發的過熱爆炸,它們的發生機理、條件及導致的后果分析與對比,見表2.3。
表2.3 熱BLEVE和冷BLEVE的比較
Table 2.3 Comparisons of Hot and Cold BLEVE
導致TLOC和BLEVE的能量來自于氣相區和液相區兩部分[4],對于快速的TLOC和BLEVE其能量主要來自于氣相區,而液相區的能量此時對BLEVE的發生并不重要,而慢速的TLOC和BLEVE,其能量主要來自于過熱液體的劇烈沸騰氣化,此時氣相區的能量只影響容器初始失效而對TLOC和BLEVE的發生影響不大,當氣相區能量不足以導致TLOC和BLEVE的發生時,液相區的能量此時就決定著TLOC和BLEVE的能否發生。
可見BLEVE的發生是由許多因素決定的[5],其中容器的壁厚、材料強度、外來熱輻射通量的大小及作用時間、泄壓閥的開啟面積和設定壓力以及熱保護措施等這些因素對BLEVE的發生起著決定性的作用。
BLEVE發生時,經常導致爆炸火球的產生,而影響爆炸火球形狀的影響因素有多種,包括容器的裝料系數、釋放的危險物質量、風速及風向、容器失效模式等。
2.3 氣云爆炸災害模擬計算模型
2.3.1 已有UVCE模擬計算模型的對比分析
根據氣云爆炸沖擊波預測模型的特點和復雜性,可以將它們分成數值模型、物理模型和相關模型三種。數值模型大多數是基于CFD(Computational Fluid Dynamics)方法,例如典型VCE數值模型有FLACS、REAGS,數值模型也許是用于模擬蒸氣云爆炸的理想方法,但由于其需要高性能的計算機以及模擬計算時間長等缺點,使得此方法在UVCE模擬方面的應用和推廣受到限制;而物理模型屬于簡化模型,例如典型VCE物理模型有SCOPE模型、CLICHÉ模型,物理模型用一個簡化的方法來描述VCE的物理過程,能夠預測大范圍內的爆炸超壓,但由于其簡化了VCE的爆炸過程,因而此類模型的模擬與預測精度受到影響。相關模型也就是縮放比率模型,依靠實驗結果而建立起來的,典型的UVCE相關模型包括TNT模型、TNO模型、ME模型、CAM模型[10];由于相關模型易于應用于災害風險評價方面,所以本文將以此類模型為研究重點。
TNT當量法和TNO(Multi-Energy)模型法[12]是蒸氣云爆炸(UVCE)模擬方法中的兩個典型模型。TNT當量法是把氣云爆炸的破壞作用轉化成TNT爆炸的破壞作用,從而把蒸氣云的量轉化成TNT當量。TNT當量法簡單易行,但有其明顯缺陷:
(1)TNT爆炸時爆源體積可忽略,而蒸氣云較大不能忽略,且隨著爆炸的進行,爆源體積在增大。
(2)TNT爆炸時能量是瞬間釋放的,而蒸氣云爆炸過程中能量的釋放速率是有限的。
(3)TNT爆炸過程形成的沖擊波強度大,但衰減速度快,而蒸氣云爆炸多屬爆燃過程,正壓作用時間較短,負壓作用時間較長。因而TNT當量法只適用于很強的蒸氣云爆炸且用以模擬爆炸遠場時偏差較小,模擬爆炸近場時高估蒸氣云爆炸產生的超壓。
(4)TNT當量法的當量系數難以確定,可變性大(0.02%-15.9%)。
TNT當量法關鍵模型:
WTNT=aWQ/QTNT (2-1)
z = R/(WTNT)1/3 (2-2)
Pi = (3.9/z1.85) + (0.5/z) (2-3)
WTNT,kg;a為LPG蒸氣云當量系數(統計平均值為0.04,占統計的60%);W為蒸氣云中LPG質量,kg;Q為LPG的燃燒熱,J/kg; QTNT為TNT的爆炸熱,J/kg; z為R處的爆炸特征長度;Pi為R處的爆炸超壓峰值。
由式(2-1)計算出對氣云爆炸有貢獻的爆炸物的當量,由式(2-2)可以確定目標R處的爆炸特征距離,再由式(2-3)計算出目標R處的爆炸超壓峰值。
TNO模型基本觀點是約束條件是增強氣云爆炸威力的關鍵因素,只有受約束的那部分氣云才對爆炸強度有作用,而不受約束的那部分蒸氣云幾乎對爆炸強度沒有貢獻;TNO模型以半球形氣云為模型,假設中心點火,火焰以恒定的速度傳播,從而以數值方法計算不同燃燒速度下的氣云爆炸強度,獲得一組爆炸強度曲線。TNO模型在理論上比較合理,但在具體應用中仍存在一些缺點:(1)如何確定受限區域的尺寸是個難以解決的問題。(2)忽略了處于開敞空間的氣云對爆炸強度的貢獻,而在實際中,當氣云形成時間較長,反應激烈的氣云,在開敞空間形成了均勻混合物,則對蒸氣云爆炸強度作用很大。(3)如果將整個氣云分成幾個爆炸源,它們的爆炸強度的疊加時的問題。(4)爆炸強度圖中的10個級別如何選取。圖中10條曲線代表不同的爆源強度(火焰速度不同),曲線10代表爆轟情況,曲線1-10依次增大;曲線6-7代表爆源強度居中的情況,用于一般氣云爆炸時,模擬遠場超壓較合適。曲線,見圖2.7。
TNO關鍵模型:
△ps、=△ps/p0 (2-4)
R、=R/(Ec/p0)1/3 (2-5)
Ec=3.5V (2-6)
r0、=R0/(Ec/p0)1/3 (2-7)
根據受限氣云的體積V就可以由式(2-6)計算出對爆炸有貢獻的燃燒能量,由式(2-5)計算出目標R處的R、,然后就可以從圖2.7查得△ps、,再由式(2-4)計算出目標R處的爆炸超壓值。
ME模型與TNO模型有很多的相似,ME模型認為由湍流火焰產生的沖擊波和爆炸過程中存在約束作用,在有障礙物的情況下,將產生巨大的沖擊波,而在沒有障礙物的條件下沖擊波將不會產生。ME模型能夠估算出UVCE的最大靜壓、最大動壓、正相超壓持續時間和超壓時間圖形。
CAM模型和ME模型有很多共同點,起初此模型建立時,爆炸的源強度是由決策樹得到的,而爆炸沖擊波的衰減是由一個簡單公式得到的,后來此模型進行了修改,考慮到了脈沖持續時間和形狀,現在此模型能夠對有障礙物的條件下的VCE進行較細致的預測。
為了檢驗上述幾個典型相關模型模擬結果的有效性,利用兩個實際案例來評價這些模型。由上述相關模型得到的爆炸超壓隨離爆源中心距離的空間分布與Flixborough事故(1974)數據的比較,如圖2.8,和La Mede 精煉廠事故 (1992) [14, 15] 的比較,如圖2.9。
從兩個對比圖中可見,上述的幾個相關模型的預測結果和事故案例觀察數據基本是相似的。然而,這些模型中的TNO模型的結果是較差的,它的結果一直比觀察到的結果要小;而TNT模型在爆炸近場高估了爆炸超壓值,而在爆炸遠場又低估了爆炸超壓值,ME模型和CAM模型的結果和觀察值吻合性較好。
2.3.2 已有BLEVE模擬計算模型的對比分析
火球模型包括火球最大直徑、持續時、抬升高度、火球表面熱輻射通量、目標接受熱劑量模型五個子模型。火球模型根據火球尺寸分為近地面火球模型和抬升火球模型,傳統的火球模型都忽略了由于LPG初始噴射使得火球中心有一定的抬升高度,即為近地面火球模型。對于慢速BLEVE來說由于有一定時間的初始噴射,火球有一個產生、燃燒擴大抬升和消失的過程,且火球有一定的抬升高度。近地面火球模型和抬升火球模型主要區別在于目標接受的火球熱輻射劑量模型不同,而在火球最大直徑和持續時間一樣。
根據火球熱輻射通量模型不同分為圓柱火焰模型(火球表面熱輻射通量假為常數)和點源模型(火球表面熱輻射通量依賴于火球中的LPG質量和持續時間)。其中典型的點源模型包括ILO、H. R. Greenberg & J . J. Cramer等模型,這些模型的模擬結果誤差較大,其中ILO模型沒有考慮容器壓力的影響,而H. R. Greenberg和J. J. Cramer模型沒有考慮容器形狀和空氣吸收的影響。但是,由于實驗數據的缺乏,這些模型的有效性有待進一步檢驗。而且輸入參數可變性大,沒有具體的用戶使用指南,妨礙了它們的普遍實際應用。
由于BLEVE發生后難以確定容器內殘余的LPG,因而在建立火球模型時,根據最大危險性原則,假設容器內全部LPG都消耗在火球中。火球最大直徑和持續時間與火球中的LPG質量成比例關系:
Dmax=aWb (2-8)
t=cWd (2-9)
a、b、c、d,不同的模型取值不同,典型的火球直徑及持續時間模型,見表2-4。
表2.4典型的火球直徑及持續時間模型
Table 2.4 Typical models for fireball diameter and duration
2.3.3 UVCE災害模擬評價模型
根據上述對蒸氣云爆炸模型的對比分析,本文將采用修正TNT模型來對UVCE進行災害模擬與預測。蒸氣云爆炸主要因沖擊波造成傷害,因而按超壓-沖量準則確定人員傷亡區域及財產損失區域。沖擊波超壓破壞準則[25],見表2.5。
表2.5 沖擊波超壓破壞、傷害準則
爆炸總能量:E=1.8aWQ (2-10)
E為LPG的爆炸總能量,J;式中1.8為地面爆炸系數
爆炸傷害半徑R: R=C(NE)1/3 (2-11)
C為爆炸實驗常數,取值;0.03~0.4;N為有限空間內爆炸發生系數,取10%
爆炸沖擊波正相最大超壓ΔP[25]:
ln(ΔP/P0)=-0.9216-1.5058ln(R、)+0.167 ln2(R、)-0.0320 ln3(R、) (2-12)
R、=D/(E/P0)1/3 (2-13)
R、無量綱距離;D為目標到蒸氣云中心距離,m;P0為大氣壓
死亡半徑指人在沖擊波作用下頭部撞擊致死半徑,由下式確定:
R1=1.98WP0.447 (2-14)
WP為LPG蒸氣云的丙烷當量(kg)
重傷半徑指人在沖擊波作用下耳鼓膜50%破裂半徑,由下式確定:
R2=9.187WP1/3 (2-15)
輕傷半徑指人在沖擊波作用下耳鼓膜1%破裂半徑,由下式確定:
R3=17.87WP1/3 (2-16)
財產損失半徑指在沖擊波作用下
上一篇:成都市學校安全工作管理實現網絡化 下一篇:沙塵暴天氣對輸電安全影響初步查明