金雪蓮

（遼寧工業(yè)大學 理學院，遼寧 錦州 121001）

1 問題的提出

當發(fā)生有毒氣體突發(fā)性泄露事故時，有關部門需要快速對泄漏源進行定位和識別，并科學預測有毒氣體的蔓延及影響范圍。其中，分析有毒氣體在大氣中的擴散是泄露事故后果分析的重要內(nèi)容，其目的在于定量地描述泄漏事故對人員和環(huán)境造成傷害的程度，并預測危害后果。2010年，一輛裝載環(huán)氧乙烷的運輸車輛在某國道上側(cè)翻，造成有毒氣體外泄。雖經(jīng)消防人員緊急處置，在等待救援和處置過程中，仍有大約1 000 個單位質(zhì)量的環(huán)氧乙烷氣體擴散到周邊區(qū)域。為了分析泄露事故可能引發(fā)的后果，在以事發(fā)點為中心東西南北各3 km、高度150 m 范圍的空間區(qū)域中，監(jiān)測部門對毒物的濃度進行了抽樣測量(見表1)。由于公路北側(cè)是農(nóng)田而南側(cè)是綠化林帶，可能導致公路南北兩側(cè)的污染程度有一定差異，因此抽樣測量分為2 個部分進行，其中前3 列為取樣點的坐標，第4 列為取樣點處毒物的濃度。

表1 2 個地區(qū)毒氣濃度的抽樣結(jié)果

2 基本假設

假設1有毒氣體釋放瞬間是向各個方向等強度釋放的。

假設2氣體在整個擴散過程中不考慮溫度、風速、風向、熱傳遞、熱輻射的因素。

假設3氣體在整個擴散過程中遵守理想氣體狀態(tài)方程。

假設4氣體的擴散服從擴散定律。

假設5遵守質(zhì)量守恒定律。

3 模型建立

3.1 有毒氣體濃度的變化規(guī)律

建立空間直角坐標系，將毒氣泄漏瞬間記作t=0，泄漏點選為坐標原點，根據(jù)假設4，每秒通過法向面積為1 m2的流量為：

其中：grad 表示梯度；σ表示氣體的擴散系數(shù)；q表示毒氣的流量；負號表示氣體的擴散方向是從濃度高的地方向濃度低的方向擴散。

圖1 為模型的理想化圖形。

圖1 模型的理想化圖形

當時間間隔為Δt時，通過毒氣域Ω的流量為：

其中：Q釋放出有毒氣體的總量；S表示包圍毒氣域所形成的空間曲面；n包圍毒氣域所形成的空間曲面的法向量。

時間間隔為Δt時，毒氣域Ω內(nèi)氣體的增量為：

其中：C(x,y,z,t)表示t時刻氣體在任意一點處的濃度。

由假設5 可知：

由高斯公式可以得到：

其中：div 為散度。

由積分中值定理得：

根據(jù)假設1，在t=0 時刻氣體泄漏點處的氣體濃度為：

其中δ(x,y,z)表示點源函數(shù)，由(7)可以解得偏微分方程(6)的解為：

這樣就計算出氣體在任意時刻、任意點的濃度，從結(jié)果中可以看出有毒氣體濃度C的等值面是以泄漏點為球心的球面上，并且球的半徑越大，濃度值越小。

4 模型的求解

對上節(jié)建立的模型使用MATLAB 程序進行求解，對于數(shù)據(jù)1，由于該組數(shù)據(jù)是在某一時刻對不同數(shù)據(jù)點測量的濃度值，因此時間t可以取為常數(shù)1。毒氣總量為1 000 個單位，利用該組數(shù)據(jù)可對擴散系數(shù)進行擬合，得到區(qū)域1 的擴散系數(shù)為1.6。對于第二組數(shù)據(jù)，利用同樣的方法可得其擴散系數(shù)為0.75。將其分別帶入到模型當中得到：

農(nóng)田毒氣擴散模型

綠化帶毒氣擴散模型

下面分別對這2 種擴散模型進行MATLAB 仿真，并對仿真結(jié)果進行分析。

4.1 農(nóng)田毒氣濃度圖形分析

圖2 第1 h 農(nóng)田高度為0 m 即地表毒氣濃度擴散圖

圖3 第1 h 農(nóng)田高度為20 m 毒氣濃度擴散圖

圖4 第12 h 農(nóng)田高度為0 m 即地表毒氣濃度擴散圖

圖5 第12 h 農(nóng)田高度為20 m 毒氣濃度擴散圖

從圖2～圖5 分析來看：(1)在固定時間內(nèi)隨著高度的增加毒氣濃度降低，從整個區(qū)域看，周邊區(qū)域的毒氣濃度是呈現(xiàn)降低趨勢的。(2)12 h 后，地表毒氣濃度明顯降低，同一高度的濃度雖然有所增加，但是總體濃度已經(jīng)明顯降低小于0.28。(3)12 h后，毒氣范圍縮小。

4.2 綠化帶毒氣濃度圖形分析

從圖6～圖9 分析看：(1)在固定時間內(nèi)隨著高度的增加毒氣濃度降低，從整個區(qū)域看，周邊區(qū)域的毒氣濃度是呈現(xiàn)降低趨勢的。(2)12 h 后，地表毒氣濃度明顯降低，同一高度的濃度雖然有所增加，但是總體濃度已經(jīng)明顯降低小于0.9。(3)12 h后，毒氣范圍縮小。

經(jīng)比較分析發(fā)現(xiàn)，綠化帶區(qū)域經(jīng)過12 h 后與農(nóng)田區(qū)域進行比較，農(nóng)田區(qū)域的毒氣濃度降低得更快。

5 結(jié)束語

圖6 第1 h 綠化帶高度為0 m 即地表毒氣濃度擴散圖

圖7 第1 h 綠化帶高度為20 m 毒氣濃度擴散圖

圖8 第12 h 綠化帶高度為0 m 即地表毒氣濃度擴散圖

圖9 第12 h 綠化帶高度為20 m 毒氣濃度擴散圖

比較農(nóng)田和綠化帶2 種不同區(qū)域的毒氣擴散程度可以看出，盡管2 種植被均有凈化作用，但農(nóng)田域毒氣濃度降低得更快。盡管如此毒氣擴散的程度和植被的空間體積也存在明顯關系，與風速、風向等外界條件也有關系，所以本文還是有一定局限性，仍然需要進一步研究。