成 誠,夏 天,龐奇志
(中國地質(zhì)大學(xué)(武漢)工程學(xué)院,湖北 武漢 430074)
環(huán)氧乙烷是一種常被用于制造殺菌劑的有毒的致癌物質(zhì),化學(xué)式為C2H4O,分子量為44.05,是一種有機(jī)化合物,因其特殊的性質(zhì)被廣泛應(yīng)用于許多行業(yè),如口罩的殺菌消毒。環(huán)氧乙烷具有非?;顫姷幕瘜W(xué)性質(zhì),易燃易爆炸,爆炸極限為3%~100%,并能與許多化合物發(fā)生化學(xué)反應(yīng)。環(huán)氧乙烷的熔點(diǎn)為-111℃,沸點(diǎn)為10.7℃,在低溫下呈現(xiàn)為無色透明的液態(tài),在常溫下呈現(xiàn)為無色帶有醚刺激性氣味的氣態(tài),易溶于水,密度為0.882 g/cm3。
目前,化工儲(chǔ)罐泄漏事故頻發(fā),且事故后果受環(huán)境因素的影響,急需對(duì)儲(chǔ)罐泄漏事故后果以及環(huán)境因素對(duì)儲(chǔ)罐泄漏事故后果的影響進(jìn)行研究。
對(duì)于儲(chǔ)罐泄漏事故造成的危害范圍的模擬研究,相艷景等[1]利用ALOHA軟件對(duì)環(huán)氧乙烷儲(chǔ)罐泄漏事故進(jìn)行了模擬,計(jì)算得出了事故的危害范圍;朱云峰等[2]利用ALOHA軟件模擬了環(huán)氧丙烷儲(chǔ)罐泄漏事故造成的危害范圍;田水承等[3]對(duì)氯乙烯儲(chǔ)罐泄漏事故后果的影響范圍進(jìn)行模擬,得到了可能事故場(chǎng)景下氯乙烯毒氣擴(kuò)散區(qū)域、閃火可燃區(qū)域和蒸汽云爆炸超壓沖擊波影響區(qū)域。
對(duì)于環(huán)境因素對(duì)?;沸孤┦鹿屎蠊挠绊懷芯?,孔大令[4]通過模擬不同環(huán)境條件下液氨泄漏擴(kuò)散的警戒范圍,分析了環(huán)境因素對(duì)液氨泄漏警戒范圍的影響;胡秀甲[5]利用ALOHA軟件模擬了內(nèi)外因?qū)σ郝刃孤U(kuò)散范圍的影響;鐘岸等[6]利用FLUENT軟件模擬了環(huán)氧乙烷隨泄漏速率、自然風(fēng)速和地面粗糙度的動(dòng)態(tài)擴(kuò)散規(guī)律;葛安然等[7]利用ALOHA軟件對(duì)甲醇儲(chǔ)罐泄漏后毒氣擴(kuò)散范圍的影響因素進(jìn)行了敏感性分析。
但是,目前現(xiàn)有研究中針對(duì)不同季節(jié)場(chǎng)景下環(huán)氧乙烷儲(chǔ)罐泄漏事故后果的模擬研究相對(duì)較少。鑒于此,本文分別對(duì)夏季和冬季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的中毒、池火災(zāi)和蒸氣云爆炸事故后果進(jìn)行了數(shù)值模擬研究,并對(duì)得出的事故危害范圍進(jìn)行了對(duì)比分析,以為環(huán)氧乙烷儲(chǔ)罐泄漏事故的預(yù)防與控制提供依據(jù)。
環(huán)氧乙烷泄漏后易發(fā)生氣化,人體暴露于該環(huán)境中會(huì)造成人員中毒傷亡,若點(diǎn)火延遲時(shí)會(huì)發(fā)生閃火或蒸氣云爆炸,未氣化的液態(tài)部分遇火源會(huì)形成池火災(zāi)[2]。
本文選擇某地某化工廠的環(huán)氧乙烷儲(chǔ)罐進(jìn)行數(shù)值模擬研究。環(huán)境因素會(huì)隨著季節(jié)的不同發(fā)生改變,不同季節(jié)的平均溫度不同,尤其是夏季與冬季的溫差很大,除了溫差外,平均風(fēng)速、風(fēng)向、相對(duì)濕度等環(huán)境因素也不相同。由于環(huán)氧乙烷儲(chǔ)罐發(fā)生泄漏后造成的事故后果會(huì)受到環(huán)境因素的影響,所以有必要對(duì)夏季和冬季時(shí)環(huán)氧乙烷儲(chǔ)罐泄漏事故分別進(jìn)行模擬,模擬結(jié)果可以為環(huán)氧乙烷儲(chǔ)罐不同季節(jié)的安全管理、事故預(yù)防和應(yīng)急救援工作提供參考,具有一定的實(shí)際應(yīng)用價(jià)值。
《危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲(chǔ)存設(shè)施外部安全防護(hù)距離確定方法》(GB/T 37243—2019)中規(guī)定10 cm孔徑泄漏為大孔泄漏,本文選取1~10 cm泄漏孔徑,利用ALOHA軟件對(duì)夏冬季場(chǎng)景下環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏后可能產(chǎn)生的事故后果進(jìn)行了仿真模擬。
該化工廠環(huán)氧乙烷儲(chǔ)罐區(qū)位于東經(jīng)114°04′、北緯30°38′,儲(chǔ)罐為臥式,直徑為2.75 m,長(zhǎng)為8.4 m,儲(chǔ)罐尺寸為50 m3,儲(chǔ)存溫度為5℃,儲(chǔ)存壓力為1.1 MPa,充裝系數(shù)為69%。該地區(qū)夏季平均溫度取最高月平均溫度為29℃,平均風(fēng)速取2.5 m/s,風(fēng)向?yàn)闁|南風(fēng),相對(duì)濕度為80%;冬季平均溫度取最低月平均溫度為4℃,平均風(fēng)速取2.6 m/s,風(fēng)向?yàn)槠憋L(fēng),相對(duì)濕度為75%。該儲(chǔ)罐區(qū)為水泥地面,設(shè)定儲(chǔ)罐底部距離地面為0.3 m處連接管破裂,并導(dǎo)致環(huán)氧乙烷連續(xù)性泄漏。表1為環(huán)氧乙烷儲(chǔ)罐區(qū)的環(huán)境因素。
表1 環(huán)氧乙烷儲(chǔ)罐區(qū)的環(huán)境因素
ALOHA軟件最初由美國環(huán)保署化學(xué)制品突發(fā)事件和預(yù)備辦公室、美國國家海洋和大氣管理局響應(yīng)和恢復(fù)辦公室共同開發(fā),2000年起我國研究者開始引進(jìn)[2]。ALOHA軟件采用的數(shù)學(xué)模型有高斯模型、DEGADIS重氣擴(kuò)散模型等。
由于環(huán)氧乙烷的分子量大于空氣的平均分子量,其相對(duì)空氣密度為1.52,比空氣重,因此環(huán)氧乙烷泄漏后擴(kuò)散會(huì)形成重氣云團(tuán),且能與空氣混合形成大型薄霧云并擴(kuò)散,故本次模擬時(shí)選用DEGADIS重氣擴(kuò)散模型。
DEGADIS重氣擴(kuò)散模型假設(shè)氣體擴(kuò)散濃度在側(cè)風(fēng)向上采用修正的高斯分布,而豎直方向上采用指數(shù)分布[2],其數(shù)學(xué)表達(dá)式如下:
c(x,y,z)=
(1)
式中:c(x,y,z)為任一點(diǎn)處氣體的質(zhì)量濃度(kg/m3);Cc(x)為在表面中心線處氣體的質(zhì)量濃度(kg/m3);b(x)為水平中間范圍的半寬(m);Sy(x)為水平方向上氣體擴(kuò)散的濃度比率;Sz為垂直方向上氣體擴(kuò)散的濃度比率;n為風(fēng)力常數(shù);z0為水平方向和垂直方向風(fēng)力圖基準(zhǔn)高度(m);ux為x方向的風(fēng)速(m/s);u0為z0高度上的實(shí)際風(fēng)速(m/s)。
在未遇到明火的情況下對(duì)環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散所造成的中毒事故的毒性影響區(qū)域范圍進(jìn)行了仿真模擬。并根據(jù)ALOHA軟件自動(dòng)引用的數(shù)據(jù)標(biāo)準(zhǔn),即美國緊急響應(yīng)計(jì)劃指南(ERPG)中規(guī)定的ERPG-2、 ERPG-3標(biāo)準(zhǔn)濃度值來確定中毒事故的毒性影響區(qū)域,其中ERPG-2為50 ppm(1 ppm=10-6),ERPG-3為500 ppm。當(dāng)毒氣濃度在50 ppm和500 ppm之間時(shí),會(huì)嚴(yán)重影響人體健康,會(huì)對(duì)人體某些器官造成傷害;當(dāng)毒氣濃度達(dá)到500 ppm(ERPG-3)以上時(shí),會(huì)對(duì)人體健康造成嚴(yán)重的影響并威脅生命[1]。本文利用ALOHA軟件對(duì)環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的中毒事故后果進(jìn)行了模擬計(jì)算,并確定了中毒事故的毒性影響范圍。
選取夏季環(huán)境因素參數(shù),通過ALOHA軟件模擬得到夏季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的中毒事故的毒性危害距離以及發(fā)生大孔泄漏(10 cm)時(shí)的毒性影響范圍,見表2和圖1。
表2 夏季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的中毒事故毒性危害距離
圖1 夏季環(huán)氧乙烷儲(chǔ)罐發(fā)生大孔泄漏(10 cm)時(shí)中毒事 故的毒性影響范圍Fig.1 Toxic influence range of large hole leakage(10 cm) of ethylene oxide storage tank in summer
由表2和圖1可以看出,夏季場(chǎng)景下隨著泄漏孔徑的增大,環(huán)氧乙烷儲(chǔ)罐發(fā)生泄漏擴(kuò)散所造成的中毒事故的毒性危害距離也隨之增大,且毒性危害距離隨泄漏孔徑增大的變化范圍在環(huán)氧乙烷濃度為50 ppm時(shí)比500 ppm時(shí)大。
選取冬季環(huán)境因素參數(shù),通過ALOHA軟件模擬得到冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的中毒事故的毒性危害距離以及發(fā)生大孔泄漏(10 cm)時(shí)的毒性影響范圍,見表3和圖2。
表3 冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的中毒事故毒性危害距離
圖2 冬季環(huán)氧乙烷儲(chǔ)罐發(fā)生大孔泄漏(10 cm)時(shí)中 毒事故的毒性影響范圍Fig.2 Toxic influence range of large hole leakage (10 cm) of ethylene oxide storage tank in winter
由表3和圖2可以看出,冬季場(chǎng)景下隨著泄漏孔徑的增大,環(huán)氧乙烷儲(chǔ)罐發(fā)生泄漏擴(kuò)散所造成的中毒事故的毒性危害距離也隨之增大,且毒性危害距離隨泄漏孔徑增大的變化范圍在環(huán)氧乙烷濃度為50 ppm時(shí)比500 ppm時(shí)大。
根據(jù)夏冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的中毒事故毒性危害距離的模擬結(jié)果(見表2和表3),得到夏季和冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的中毒事故的毒性危害距離折線圖,見圖3。
圖3 夏季和冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散 造成的中毒事故毒性危害距離折線圖Fig.3 Toxic hazard distance caused by leakage and diffusion of different apertures of ethylene oxide storage tank in different seasons
由圖3可見,夏季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散所造成的中毒事故的毒性危害距離均大于冬季環(huán)氧乙烷儲(chǔ)罐泄漏,且毒性危害距離的變化范圍在環(huán)氧乙烷濃度為50 ppm時(shí)比500 ppm時(shí)更明顯。
環(huán)氧乙烷儲(chǔ)罐泄漏后的液態(tài)成分遇點(diǎn)火源會(huì)形成池火災(zāi),泄漏孔徑大小不同,形成的池火災(zāi)的熱輻射強(qiáng)度也不相同[1]。由于池火災(zāi)的主要危害是熱輻射對(duì)人體的傷害,故在火災(zāi)事故后果分析中,以熱輻射強(qiáng)度來劃分危害區(qū)域。假設(shè)受到熱輻射影響的時(shí)間為1 min,當(dāng)熱輻射強(qiáng)度等級(jí)為3級(jí),即2 kW/m2時(shí),會(huì)導(dǎo)致人體疼痛;當(dāng)熱輻射強(qiáng)度等級(jí)為2級(jí),即5 kW/m2時(shí),會(huì)導(dǎo)致人體燒傷;當(dāng)熱輻射強(qiáng)度等級(jí)為1級(jí),即10 kW/m2時(shí),會(huì)導(dǎo)致人體死亡。本文采用這三個(gè)熱輻射強(qiáng)度來劃分夏季和冬季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的池火災(zāi)事故的熱輻射影響范圍。
選取夏季環(huán)境因素參數(shù),通過ALOHA軟件模擬得到夏季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的池火災(zāi)事故的熱輻射危害距離以及發(fā)生大孔泄漏(10 cm)時(shí)池火災(zāi)事故的熱輻射影響范圍,見表4和圖4。
表4 夏季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的池火災(zāi)事故熱輻射危害距離
圖4 夏季環(huán)氧乙烷儲(chǔ)罐發(fā)生大孔泄漏(10 cm)時(shí)池 火災(zāi)事故的熱輻射影響范圍Fig.4 Influence range of pool fire in case of large hole leakage(10 cm) of ethylene oxide storage tank in summer
由表4和圖4可以看出,夏季場(chǎng)景下隨著泄漏孔徑的增大,環(huán)氧乙烷儲(chǔ)罐發(fā)生泄漏擴(kuò)散所造成的池火災(zāi)事故的熱輻射危害距離也隨之增大,且池火災(zāi)事故的熱輻射危害距離在不同熱輻射強(qiáng)度(10 kW/m2、5 kW/m2、2 kW/m2)下隨泄漏孔徑增大的變化范圍也依次增大。
選取冬季環(huán)境因素參數(shù),通過ALOHA軟件模擬得到冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的池火災(zāi)事故的熱輻射危害距離以及發(fā)生大孔泄漏(10 cm)時(shí)的池火災(zāi)事故的熱輻射影響范圍,見表5和圖5。
表5 冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的池火災(zāi)事故熱輻射危害距離
圖5 冬季環(huán)氧乙烷儲(chǔ)罐發(fā)生大孔泄漏(10 cm)時(shí)池火 災(zāi)事故的熱輻射影響范圍Fig.5 Influence range of pool fire in case of large hole leakage(10 cm) of ethlene oxide storage tank in winter
由表5和圖5可以看出,冬季場(chǎng)景下隨著泄漏孔徑的增大,環(huán)氧乙烷儲(chǔ)罐發(fā)生泄漏擴(kuò)散所造成的池火災(zāi)事故的熱輻射危害距離也隨之增大,且熱輻射危害距離在不同熱輻射強(qiáng)度(10 kW/m2、5 kW/m2、2 kW/m2)下隨泄漏孔徑增大的變化范圍也依次增大。
根據(jù)夏、冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的池火災(zāi)事故熱輻射危害距離的模擬結(jié)果(見表4和表5),得到夏季和冬季場(chǎng)景下環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成池火災(zāi)事故的熱輻射危害距離折線圖,見圖6。
圖6 夏季和冬季場(chǎng)景下環(huán)氧乙烷儲(chǔ)罐不同孔徑泄 漏擴(kuò)散造成的池火災(zāi)事故熱輻射危害距離Fig.6 Heat radiation hazard distance of pool fire caused by leakage and diffusion of different apertures of ethylene oxide storage tank in different seasons
由圖6可見:冬季環(huán)氧乙烷儲(chǔ)罐泄漏所造成的池火災(zāi)事故的熱輻射危害距離略微大于夏季環(huán)氧乙烷儲(chǔ)罐泄漏。
池火災(zāi)是環(huán)氧乙烷泄漏后的液態(tài)成分遇點(diǎn)火源形成的,且環(huán)氧乙烷的沸點(diǎn)為10.7℃,平時(shí)以液態(tài)儲(chǔ)存在儲(chǔ)罐中,但當(dāng)發(fā)生環(huán)氧乙烷儲(chǔ)罐泄漏時(shí),冬季儲(chǔ)罐泄漏后汽化的環(huán)氧乙烷相比于夏季較少。
蒸氣云爆炸(VCE)是可燃性氣體混合物經(jīng)過一段延遲時(shí)間后,可燃蒸氣云被點(diǎn)火源點(diǎn)燃發(fā)生著火或爆炸,產(chǎn)生了危險(xiǎn)的沖擊波超壓。由于蒸氣云爆炸的主要危害來自沖擊波超壓,故利用ALOHA軟件對(duì)環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的蒸氣云爆炸事故后果進(jìn)行了數(shù)值模擬計(jì)算。當(dāng)沖擊波超壓強(qiáng)度為0.055 MPa時(shí),會(huì)造成建筑物毀壞;當(dāng)沖擊波超壓強(qiáng)度為0.024 MPa時(shí),會(huì)造成人員重傷;當(dāng)沖擊波超壓強(qiáng)度為0.007 MPa時(shí),會(huì)造成玻璃破碎[1]。
選取夏季環(huán)境因素參數(shù),通過ALOHA軟件模擬得到夏季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的蒸氣云爆炸事故產(chǎn)生的沖擊波超壓危害距離以及發(fā)生大孔泄漏(10 cm)時(shí)蒸氣云爆炸事故產(chǎn)生的沖擊波超壓的影響范圍,見表6和圖7。
表6 夏季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的蒸氣云爆炸事故沖擊波超壓危害距離
圖7 夏季環(huán)氧乙烷儲(chǔ)罐發(fā)生大孔泄漏(10 cm)時(shí) 蒸氣云爆炸事故產(chǎn)生的沖擊波超壓影響范圍Fig.7 Influence range of overpressure of shock wave of vapor cloud explosion in case of large hole leakage (10 cm) of ethylene oxide storage tank in summer
由表6和圖7可以看出,夏季場(chǎng)景下隨著泄漏孔徑的增大,環(huán)氧乙烷儲(chǔ)罐發(fā)生泄漏擴(kuò)散所造成的蒸氣云爆炸事故的沖擊波超壓危害距離也隨之增大,且沖擊波超壓危害距離在不同沖擊波超壓強(qiáng)度(0.055 MPa、0.024 MPa、0.007 MPa)下隨泄漏孔徑增大的變化范圍也依次增大。
選取冬季環(huán)境因素參數(shù),通過ALOHA模擬得到冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的蒸氣云爆炸事故產(chǎn)生的沖擊波超壓危害距離以及發(fā)生大孔泄漏(10 cm)時(shí)蒸氣云爆炸事故產(chǎn)生的沖擊波超壓的影響范圍,見表7和圖8。
表7 冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的蒸氣云爆炸事故沖擊波超壓危害距離
圖8 冬季環(huán)氧乙烷儲(chǔ)罐發(fā)生大孔泄漏(10 cm)時(shí) 蒸氣云爆炸事故產(chǎn)生的沖擊波超壓影響范圍Fig.8 Influence range of overpressure of shock wave of vapor cloud explosion in case of large hole leakage (10 cm) of ethylene oxide storage tank in winter
由表7和圖8可以看出,冬季場(chǎng)景下隨著泄漏孔徑的增大,環(huán)氧乙烷儲(chǔ)罐發(fā)生泄漏擴(kuò)散所造成的蒸氣云爆炸事故的沖擊波超壓危害距離也隨之增大,且沖擊波超壓危害距離在不同沖擊波超壓強(qiáng)度(0.055 MPa、0.024 MPa、0.007 MPa)下隨泄漏孔徑增大的變化范圍也依次增大。
根據(jù)夏、冬季環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的蒸氣云爆炸事故沖擊波超壓危害距離的模擬結(jié)果(見表6和表7),得到夏季和冬季場(chǎng)景下環(huán)氧乙烷儲(chǔ)罐不同孔徑泄漏擴(kuò)散造成的蒸氣云爆炸事故產(chǎn)生的沖擊波超壓危害距離折線圖,見圖9。
圖9 夏季和冬季場(chǎng)景下環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的 蒸氣云爆炸事故沖擊波超壓危害距離Fig.9 Hazard distance of vapor cloud explosion caused by leakage and diffusion of ehtylene oxide storagae tank in summer and winter
由圖9可見,夏季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的蒸氣云爆炸事故產(chǎn)生的沖擊波超壓危害距離均明顯大于冬季環(huán)氧乙烷儲(chǔ)罐泄漏,且沖擊波超壓危害距離的變化范圍在沖擊波超壓強(qiáng)度為0.007 MPa時(shí)比0.024 MPa和0.055 MPa時(shí)更明顯。
本文基于ALOHA軟件對(duì)夏季和冬季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的中毒、池火災(zāi)和蒸氣云爆炸事故后果進(jìn)行了數(shù)值模擬研究和對(duì)比分析,得出如下結(jié)論:
(1) 夏季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散所造成的中毒事故的毒性危害距離均大于冬季環(huán)氧乙烷儲(chǔ)罐泄漏,且毒性危害距離的變化范圍在環(huán)氧乙烷濃度為50 ppm時(shí)比500 ppm時(shí)更明顯。
(2) 夏季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散造成的蒸氣云爆炸事故所產(chǎn)生的沖擊波超壓危害距離均明顯大于冬季環(huán)氧乙烷儲(chǔ)罐泄漏,且沖擊波超壓危害距離的變化范圍在沖擊波超壓強(qiáng)度為0.007 MPa時(shí)比0.024 MPa和0.055 MPa更明顯。
(3) 模擬結(jié)果表明:不同于中毒事故和蒸氣云爆炸事故,冬季環(huán)氧乙烷儲(chǔ)罐泄漏擴(kuò)散所造成的池火災(zāi)事故的熱輻射危害距離略微大于夏季環(huán)氧乙烷儲(chǔ)罐泄漏。