1 概述

據(jù)不完全統(tǒng)計，2019年我國共發(fā)生燃氣事故722起，其中發(fā)生在餐飲店的事故逾100起。臨街餐飲店通常將就餐區(qū)布置在臨街一側(cè)，廚房位于內(nèi)側(cè)，許多是沒有對外門窗的內(nèi)廚房，存在通風(fēng)不良隱患。對這些餐飲店進行燃氣泄漏后果分析對控制燃氣事故的發(fā)生有重要意義。

范堅強突然打電話給一杭，約他去辦公室談《真相》的創(chuàng)作進度。一杭不太喜歡談?wù)撜趯懽鞯男≌f，他認為創(chuàng)作是很私密的事情，別人的意見容易造成不必要的干擾，何況，還沒成型的東西，說出去也是對自己的不負責。當然，這一次他準備去，就算范堅強不找他，他也要尋個機會去一趟一風(fēng)公司。他需要確認，核桃臉記事本上提到的范總是不是范堅強。

目前已有許多學(xué)者對燃氣泄漏或者爆炸過程進行了研究。Huh等人

分別采用CFD方法和簡化的室內(nèi)氣體流量模型對建筑內(nèi)氣體泄漏擴散進行模擬計算，對比驗證和評價室內(nèi)釋放模型的性能。Litto等人

模擬建筑物內(nèi)貧甲烷泄漏的氣流規(guī)律，研究門、通風(fēng)口和泄漏位置的影響，結(jié)果表明建筑主門保持關(guān)閉時，大部分無組織甲烷通過通風(fēng)口排出。陳琪等人

以民用廚房為研究對象，模擬在門窗關(guān)閉的條件下微量連續(xù)點源泄漏的甲烷氣體擴散運動特性，得出室內(nèi)危險區(qū)域范圍。Sun等人

在鋼制燃燒室中進行甲烷-空氣混合物的爆炸試驗，研究了初始壓力、濃度、初始湍流對火焰發(fā)展的影響。Ji等人

對建筑內(nèi)天然氣泄漏爆炸事故進行調(diào)查研究，將不同條件下的爆炸模擬結(jié)果與實際情況對比，得出事故發(fā)生前建筑內(nèi)可能存在的天然氣填充量。Pedersen等人

采用CFD方法對兩室封閉空間內(nèi)天然氣爆炸進行模擬研究，結(jié)果表明超壓曲線有2個壓力峰值。第1個壓力峰值與防爆通風(fēng)板泄壓有關(guān)，第2個壓力峰值隨排氣口結(jié)構(gòu)、填充程度而變化。Sezer等人

設(shè)計各種尺寸和幾何形狀的通風(fēng)口，研究其對通風(fēng)室爆炸后果產(chǎn)生的影響。Li等人

提出了一種用FLACS模擬封閉空間內(nèi)湍流火焰的CFD模擬方法，并用FLUENT軟件計算了低湍流爆炸波在自由空氣中的傳播，將其應(yīng)用于預(yù)測由外部爆炸引起的爆炸室內(nèi)的峰值壓力，計算結(jié)果與實驗結(jié)果吻合較好。李紅培

建立開放式廚房燃氣泄漏擴散的模型，利用FLUENT模擬不同影響因素下開放式廚房內(nèi)燃氣擴散，結(jié)果表明開放式廚房開窗有風(fēng)情況和加隔墻情況會降低事故風(fēng)險發(fā)生的概率。Baalisampang等人

建立一種分析模型來評估火災(zāi)、爆炸產(chǎn)生的超壓和熱輻射對人體形成的危害后果。目前幾乎沒有學(xué)者對臨街餐飲店內(nèi)天然氣泄漏爆炸事故進行模擬研究，筆者對此進行研究，并對比分析產(chǎn)生的后果。

本文采用計算流體力學(xué)方法，建立兩種長寬比的餐飲店物理模型，對門、風(fēng)機開閉狀態(tài)下的天然氣泄漏擴散進行模擬，計算得到天然氣質(zhì)量分數(shù)分布。然后在FLACS軟件中設(shè)置點火點進行爆炸模擬，對比研究不同泄漏時間進行點火的超壓分布情況，研究餐飲店門開閉狀態(tài)對超壓分布的影響，分析爆炸過程中餐飲店內(nèi)超壓分布情況、氣流運動情況。針對最危險的模擬工況，分析對人體傷害以及對建筑物的破壞效應(yīng)。

2 臨街餐飲店模型與模擬參數(shù)設(shè)置

本文坐標單位均為m。

2.1 物理模型

① 泄漏擴散物理模型

利用FLUENT前處理模塊ICEM建立典型臨街餐飲店泄漏擴散物理模型，見圖1，模型平面圖見圖2(圖中尺寸單位為mm)，模型尺寸見表1，模型內(nèi)物體

坐標見表2。餐飲店門和窗口位于餐廳

軸中間位置，餐廳門寬2 m，高2 m。窗口(是個孔洞)位于餐廳門正上方，寬1 m，高0.3 m。廚房門距離

平面2.2 m，寬0.8 m，高3 m。模型主要由就餐區(qū)和廚房組成，廚房門始終保持開啟，厚度為0.1 m。矩形送、排風(fēng)管道寬0.35 m，中心距地面2.625 m，送、排風(fēng)口為圓形向下，直徑為0.3 m。燃氣管道內(nèi)直徑為32 mm，貼附地面設(shè)置，泄漏點位于水平燃氣管道端頭的一段豎直管段上。

圖中

——餐飲店總長度，m

④ 格雷馬斯.行動元、角色和形象.選自張寅德.敘事學(xué)研究.王國卿譯.中國社會科學(xué)出版社，1989年:第119-125頁.

——就餐區(qū)長度，m

——就餐區(qū)和廚房的隔墻厚度，m

——廚房長度，m

為了探討農(nóng)村地區(qū)收入不平等究竟是由城市群與非城市群間(組間)，還是其組內(nèi)收入差距擴大所引起的，本文采用 Shorrocks[24]的子群體分解法來進行分析。由于使用基尼系數(shù)進行子群體分解時，存在“剩余項”，不適合用來進行子群體分解[25]。因此，本文將中國農(nóng)村劃分為城市群地區(qū)與非城市群地區(qū)，用泰爾指數(shù)來度量不平等，結(jié)果如表1所示。通過泰爾指數(shù)分解，可知：

——餐飲店寬度，m

——餐飲店高度，m

《治譜》主要依據(jù)從政流程與行政事務(wù)分門別類進行編排，其條目如：初選門；到任門；堂事門；詞訟門；錢糧門；人命門；賊盜門；獄囚門；待人門；(雜事門、補遺)；(憲約、慎刑說)[26](P85)。其中《賊盜門·截殺法》曰：

進行爆炸模擬時，為便于了解沖擊波對餐飲店外部空間的影響程度，在泄漏擴散物理模型的基礎(chǔ)上擴大餐飲店外部空間模擬區(qū)域，得到爆炸物理模型(采用FLACS軟件的前處理器CASD建立)，見圖3。模型1擴大后稱為模型3，模型2擴大后稱為模型4。模型3長27 m，寬24 m，高度為8 m；模型4長36 m，寬24 m，高度為8 m。當餐飲店門關(guān)閉時，將餐飲店門設(shè)置成POPOUT類型泄壓板，其開啟壓力為0.01 MPa。

第五，因為地震后原有家園被摧毀，一部分羌族村寨被迫遷移。原汶川縣龍溪鄉(xiāng)直臺村和垮坡村夕格組原有的生存空間不再適宜生存，于是被迫遷往了位于邛崍市的南寶農(nóng)場。由原本的羌族聚居區(qū)遷往漢人區(qū)，不僅是生存空間的遷移，更是文化生存空間的改變。

2.2 網(wǎng)格劃分

① 泄漏擴散物理模型

采用ICEM劃分非結(jié)構(gòu)網(wǎng)格，泄漏點附近通常具有較大的速度梯度、濃度梯度和壓力梯度，因此對其上方區(qū)域進行局部網(wǎng)格加密。模型1網(wǎng)格數(shù)量為1 141 234個，模型2網(wǎng)格數(shù)量為1 699 798個。

② 爆炸物理模型

采用CASD劃分網(wǎng)格，網(wǎng)格類型為三維笛卡爾網(wǎng)格，是邊長為0.10 m的六面體網(wǎng)格。

2.3 參數(shù)設(shè)置

① 泄漏擴散

調(diào)控進度 創(chuàng)造過程經(jīng)歷準備期、醞釀期、明朗（頓悟）期和驗證期．其中，醞釀期對獲得創(chuàng)造性想法至關(guān)重要，大量觀念的組合和選擇就發(fā)生在這個時期．大量研究支持醞釀期的意識思維休息假說及遺忘固著假說．因此，教師一方面應(yīng)提供富有挑戰(zhàn)性的任務(wù)讓學(xué)生體驗醞釀期，另一方面應(yīng)調(diào)控教學(xué)的進度，保證學(xué)生有足夠的醞釀時間，如在課堂話語和活動中進行短暫的休息，讓學(xué)生在一段時間內(nèi)同時處理多個任務(wù)或完成長時作業(yè)等．

采用FLUENT軟件進行泄漏擴散計算時，湍流模型采用

-

模型，求解模型采用PISO模型。模型中餐桌、灶臺、吸油煙機、冰箱等家具家電以及墻面、風(fēng)管、燃氣管道均設(shè)置為壁面邊界條件，窗口邊界條件為Outflow。餐廳門關(guān)閉時，餐廳門設(shè)置為壁面邊界條件；餐廳門開啟時，餐廳門設(shè)置為Outflow。風(fēng)機關(guān)閉時，送風(fēng)口與排風(fēng)口為壁面邊界條件；風(fēng)機開啟時，送風(fēng)口與排風(fēng)口邊界條件為Velocity，送、排風(fēng)速均為8.96 m/s，湍流強度為3.47%。天然氣假設(shè)為純甲烷。初始絕對壓力為101.325 kPa，初始溫度為25 ℃。

模型1和模型2中均設(shè)置了1個天然氣體積分數(shù)監(jiān)測點，當監(jiān)測到天然氣體積分數(shù)達到1%時開啟送風(fēng)機和排風(fēng)機。模型1監(jiān)測點坐標為(3.0，7.0，2.7)，模型2監(jiān)測點坐標為(3.0，14.0，2.7)。

泄漏點坐標為(3.1，7.9，0.1)，面積為50 mm

，泄漏方向沿

方向。泄漏點邊界條件為velocity-inlet，泄漏速度為54.18 m/s，采用小孔泄漏模型流量計算方法

計算得泄漏流量為9.8 m

/h。

要想更好地對公司進行管理，就需要從內(nèi)到外地進行人力資源管理工作。首先調(diào)整人力資源管理的策略，不僅要從客戶方面考慮，還要結(jié)合市場以及其他部門，人力資源管理應(yīng)該最大限度地為企業(yè)創(chuàng)造利益和價值。此外，人力資源管理者需要具備提高企業(yè)吸引度，增加企業(yè)投資者的能力，站在他人或受益者的角度進行考慮的能力。

② 爆炸物理模型

由圖8可知，當餐飲店內(nèi)氣云燃燒形成的壓力未達到餐飲店門破壞壓力(10 kPa)時，主要在餐飲店內(nèi)部形成超壓。隨著燃燒產(chǎn)物逐漸增加，餐飲店內(nèi)壓力積聚到10 kPa時，餐飲店門開啟，超壓傳播到外部，超壓峰值達18 kPa以上。模型1泄漏1 h點火后200 ms、574 ms時

=0.25 m、

=0.75 m平面上的速度分布見圖9。

采用FLACS軟件進行爆炸計算時，根據(jù)泄漏擴散模擬得到的結(jié)果設(shè)置餐飲店內(nèi)的天然氣體積分數(shù)，模型采用歐拉(EULER)邊界條件，點火點設(shè)于灶臺附近，模型3點火點坐標(3.2，7.9，0.2)，模型4點火點坐標(3.2，14.9，0.2)。設(shè)置初始絕對壓力為101 325 Pa，初始溫度為20 ℃。輸出參數(shù)為爆炸超壓和氣流速度。

③ 壓力測點設(shè)置

定義直線1的方程為：

對于模型1、3，沿直線1，在

=-15、-10、-6、-3、0、2、4、6、8處布置壓力測點。

對于模型2、4，沿直線1，在

=-15、-10、-6、-3、0、2、4、6、8、11、13、15處布置壓力測點。

3 泄漏擴散模擬結(jié)果分析

3.1 餐飲店門開閉對泄漏擴散的影響

當討論到由“數(shù)據(jù)分析”做主導(dǎo)的企業(yè)面臨的最大挑戰(zhàn)時，與會者們一致認為投資和管理的焦點過于看重科技，卻忽略了人和分析流程。

模型1風(fēng)機啟動時刻為167 s，模型2風(fēng)機啟動時刻為138 s。由圖7可知，風(fēng)機啟動前，天然氣體積分數(shù)大于1%的氣云體積隨時間呈增加趨勢；風(fēng)機啟動后，體積分數(shù)迅速下降，各模型天然氣體積分數(shù)大于5%的氣云體積都在風(fēng)機啟動后的20 s內(nèi)降低至0?？芍L(fēng)能在短時間內(nèi)有效控制餐飲店內(nèi)泄漏天然氣的擴散。

對比圖4、5可知，在餐飲店門開啟和關(guān)閉條件下，天然氣泄漏后的擴散趨勢總體相同。但由于餐飲店門開啟時，可通過餐飲店門直接和外面大氣進行質(zhì)量交換，因此天然氣質(zhì)量分數(shù)比餐飲店門關(guān)閉時明顯降低。圖5中天然氣等質(zhì)量分數(shù)線明顯不水平，餐飲店門附近存在較大的質(zhì)量分數(shù)梯度；而餐飲店門關(guān)閉時，等質(zhì)量分數(shù)線基本水平，餐飲店門附近不存在大的質(zhì)量分數(shù)梯度。對比圖4a、5a可知，泄漏10 min時，餐飲店門開啟時就餐區(qū)的天然氣質(zhì)量分數(shù)大于餐飲店門關(guān)閉時。對比圖4c、5c可知，泄漏2 h，餐飲店門開啟時，模型1內(nèi)天然氣質(zhì)量分數(shù)總體低于0.07，而餐飲店門關(guān)閉時，模型1內(nèi)質(zhì)量分數(shù)均在0.15以上。對比圖4f、5f可知，泄漏2 h，餐飲店門開啟時，模型2內(nèi)質(zhì)量分數(shù)總體在0.04以下；餐飲店門關(guān)閉時，模型2內(nèi)質(zhì)量分數(shù)大于0.05?？傮w來說，由于餐飲店門開啟狀態(tài)下天然氣擴散速度更快，泄漏剛開始時，天然氣質(zhì)量分數(shù)更高。一段時間后，天然氣通過餐飲店門不斷擴散至外部，天然氣質(zhì)量分數(shù)比餐飲店門關(guān)閉時低。

農(nóng)村空巢老人經(jīng)濟供養(yǎng)水平較低，貧困程度相對較高。由于年齡限制了身體各方面的機能，年齡越高，身體機能越弱，勞動能力隨之下降。但是，空巢老人家庭條件差，老人的收入來源不固定，大多數(shù)老人選擇進行農(nóng)作物耕種來獲取收入，因此老人投入農(nóng)耕的時間和精力便增多，勞動強度大，一旦遇到天災(zāi)人禍，農(nóng)作物必須及時收割或直接放棄來年繼續(xù)耕作，這對老人來說無疑是巨大的挑戰(zhàn)。

餐飲店門開閉時，不同天然氣體積分數(shù)范圍對應(yīng)的氣云體積隨時間的變化曲線見圖6。

原始數(shù)據(jù)中除面狀河流外、同時存在線狀河流數(shù)據(jù)。其中線狀河流與面狀河流共有以下三種拓撲關(guān)系：相交、相接、相離。圖 8（a）、（b）、（c）即為原始河流數(shù)據(jù)中線面三種狀態(tài)的可視化。圖 8（d）、（e）、（f）為運用本文算法、并對應(yīng)三種拓撲關(guān)系的中軸線提取結(jié)果。對于原本由線狀河流柵格化得到的單線條河流都能夠很好的保留。當線狀河流與面狀河流存在重疊或相連時，面狀河流提取中軸線與單線條河流在重疊或相連處仍可以保持完整連通的狀態(tài)，不會對提取結(jié)果造成不良影響，實現(xiàn)了全自動化的中軸提取。

定義 1[10] 對簡單圖G,設(shè)映射E(G)→{1,2,…,k}滿足:對?uv,uE(G)(v≠w),有f(uv)≠f(uw),其中N，則稱f為G的k-正常邊染色,記為G的k-PEC,稱χ′(G)=min{k|G有k-PEC}是圖G的邊色數(shù)。

圖6中，無論餐飲店門開閉，兩種模型中天然氣體積分數(shù)大于1%的氣云體積均隨時間增加至餐飲店體積后保持恒定。對比圖6a、6c，餐飲店體積越小，空間內(nèi)天然氣體積分數(shù)達到1%所需的時間越短。由圖6a、6c可知，在餐飲店門關(guān)閉條件下，天然氣體積分數(shù)在爆炸極限范圍內(nèi)的體積隨時間增加至餐飲店體積后保持一段時間，然后逐漸減小至0。由圖6b、6d可知，對于模型1，餐飲店門開啟，泄漏剛開始時天然氣體積分數(shù)大于1%的氣云體積稍高于體積分數(shù)大于5%的氣云體積，而1 h后兩者幾乎相同，泄漏1～4 h餐飲店內(nèi)天然氣體積分數(shù)都保持在爆炸極限范圍內(nèi)。對于模型2，餐飲店門開啟，泄漏4 h內(nèi)天然氣體積分數(shù)大于1%的氣云體積均大于體積分數(shù)大于5%的氣云體積，表明泄漏4 h內(nèi)餐飲店內(nèi)僅部分氣云的天然氣體積分數(shù)達到爆炸極限，并且在爆炸極限范圍內(nèi)的氣云體積呈先增加后降低的趨勢。

3.2 送、排風(fēng)機啟停對泄漏擴散的影響

餐飲店門保持關(guān)閉狀態(tài)，保持其他條件不變，當體積分數(shù)監(jiān)測點監(jiān)測到天然氣體積分數(shù)達到1%時，開啟送、排風(fēng)機對餐廳進行通風(fēng)。不同天然氣體積分數(shù)對應(yīng)的氣云體積隨時間的變化曲線見圖7。

由圖4可知，天然氣泄漏后，以泄漏點為中心，受浮力作用向四周擴散。天然氣由廚房擴散至就餐區(qū)的過程中，先貼附屋頂從廚房門擴散到就餐空間，在浮力和動量的作用下，沿屋頂和墻壁擴散至整個就餐區(qū)，泄漏2 h后兩個模型內(nèi)天然氣質(zhì)量分數(shù)均在0.07以上?？傮w來看，廚房內(nèi)下部空間的天然氣質(zhì)量分數(shù)平均值高于上部空間，而就餐區(qū)下部空間的天然氣質(zhì)量分數(shù)明顯低于上部空間。兩種模型天然氣擴散規(guī)律基本一致，但同時刻模型1天然氣質(zhì)量分數(shù)均大于模型2，泄漏2 h后模型1內(nèi)天然氣質(zhì)量分數(shù)幾乎在0.15以上，模型2內(nèi)天然氣質(zhì)量分數(shù)在0.15以上的區(qū)域較小?？芍惋嫷牦w積(指餐飲店內(nèi)部空間的體積)對天然氣擴散的影響較大，體積越小，泄漏相同時間后天然氣質(zhì)量分數(shù)越高。

4 爆炸模擬結(jié)果分析

4.1 餐飲店爆炸過程分析

僅以風(fēng)機關(guān)閉、餐廳門關(guān)閉狀態(tài)下，模型3內(nèi)泄漏流量為9.8 m

/h，持續(xù)泄漏1 h被點燃爆炸為例。點火574 ms后餐飲店內(nèi)達到超壓峰值，點火后各時刻(200 ms、450 ms、574 ms)

=0.75 m平面上的超壓分布見圖8。

② 爆炸

由圖9可知，點火后燃燒反應(yīng)釋放的熱量以及燃燒產(chǎn)物膨脹使附近氣體向四周運動，其速度沿反應(yīng)前沿向外圍依次遞減，氣流經(jīng)過廚房門時節(jié)流，使得此處速度增加。隨著燃燒范圍不斷擴大，氣流速度也不斷增加，到點火發(fā)生574 ms時，大量燃燒產(chǎn)物通過餐飲店門節(jié)流運動至外部大氣空間，此時最大氣流速度在420 m/s以上。

4.2 餐飲店門開閉對超壓的影響

不通風(fēng)時，分別在兩種模型泄漏0.5 h、1 h、1.5 h、2 h點火，對餐飲店門開啟和關(guān)閉兩種狀態(tài)進行模擬。由圖6a可知，模型1泄漏2 h后餐飲店內(nèi)天然氣體積分數(shù)均大于爆炸上限，因此模型3泄漏2 h點火時，模擬空間的初始天然氣體積分數(shù)按爆炸上限(15%)進行設(shè)置，模型內(nèi)各測點超壓峰值見圖10。

由于長柄雙花木不同發(fā)育階段的植物經(jīng)濟策略不同，因此對其進行保護和保育時應(yīng)針對其不同發(fā)育階段采取相應(yīng)的保護措施。長柄雙花木幼苗需要較高的葉含水量，有研究表明長柄雙花木幼苗葉片水含量和光合速率隨遮蔭程度的增加而上升[46]，其幼苗適宜生長在較濕潤且遮蔭的環(huán)境中。因此在遷地保護時應(yīng)注意對幼苗進行適當?shù)恼谑a，以促進幼苗的生長；而在幼樹和成樹階段，葉厚度和葉面積逐漸增大，植株需要較強的光照進行光合作用。因此對于就地保護和遷地保護種群的幼樹和成樹，應(yīng)對其周圍樹木進行適當修剪與清理，以保證其對光照的需求。

以泄漏開始時刻為0時刻，時間用

表示。不啟動送風(fēng)機和排風(fēng)機，餐飲店門關(guān)閉、開啟時，模型1和2內(nèi)各時刻天然氣質(zhì)量分數(shù)分布云圖分別見圖4、5。

點火后，溫度迅速升高引起氣體膨脹，當受到墻體、家具等障礙物的阻礙時，火焰?zhèn)鞑ニ俣燃涌烨覊毫Ψe聚。同一種模型，泄漏時間越長，點火形成的超壓越大。爆炸發(fā)生后，兩種模型餐飲店內(nèi)最靠近點火位置的測點處的超壓峰值最大，各測點超壓峰值隨著與點火點距離增加而逐漸減小。對比圖10a、10c可知，在餐飲店門關(guān)閉的條件下，泄漏0.5 h、1 h進行點火，此時模型3內(nèi)形成的可燃氣云體積更大，模型3內(nèi)的超壓峰值比模型4大；隨著泄漏時間的延長，在泄漏1.5 h、2 h進行點火，模型3內(nèi)各測點超壓峰值均小于模型4。對比圖10a、10b可知，餐飲店門關(guān)閉時，氣云燃燒形成的超壓以及燃燒產(chǎn)物未能有效泄放，在餐飲店內(nèi)部積聚直到餐飲店門被破壞；餐飲店門開啟時，壓力波迅速運動至外部空間，因此相同測點處餐飲店門關(guān)閉形成的超壓峰值比餐飲店門開啟時大?？傮w來說，泄漏剛開始的一段時間內(nèi)，體積較小的餐飲店內(nèi)進行點火爆炸產(chǎn)生的超壓更大；天然氣泄漏時間較長時，體積較大的餐飲店點火爆炸產(chǎn)生的超壓更大。

4.3 爆炸后果分析

① 超壓對人體的影響

超壓對人體產(chǎn)生傷害，本文采用概率函數(shù)法

評估超壓對人體各部位產(chǎn)生的傷害概率。由前文可知，餐飲店門關(guān)閉的條件下，模型4泄漏2 h后點火爆炸內(nèi)部產(chǎn)生的超壓峰值最大，對這種情況下爆炸過程餐飲店內(nèi)部和外部各測點處超壓對人體各部位產(chǎn)生的傷害概率進行計算。

計算結(jié)果顯示，餐飲店內(nèi)部超壓導(dǎo)致的人體傷害概率大于餐飲店外部，餐飲店內(nèi)部各測點處的超壓對人體頭部的傷害概率接近100%，而餐飲店外部超壓基本不對人體產(chǎn)生傷害。人體位移傷害概率隨距點火點距離增加而逐漸減小。在餐飲店內(nèi)外，對肺部造成損傷的可能性幾乎為0。

② 超壓對建筑物的影響

采用沖擊波對建筑物(4層以下)影響的概率函數(shù)法

，計算超壓對建筑物的損壞概率。仍針對模型4泄漏2 h進行點火爆炸的情況，對各測點超壓對建筑物的損壞概率進行計算。結(jié)果表明，餐飲店內(nèi)部各測點處發(fā)生玻璃窗損壞、主要結(jié)構(gòu)破壞以及建筑倒塌概率都在99.9%以上；餐飲店外部，計算域內(nèi)的玻璃窗損壞概率在99.9%以上，主要結(jié)構(gòu)損壞概率在32%以下，建筑倒塌概率在5%以下。

③ 其他影響

爆炸過程中，除了沖擊波引起的超壓變化外，高速爆炸風(fēng)產(chǎn)生的沖擊動壓具有與超壓同樣的破壞效應(yīng)。模型3泄漏2 h點火爆炸，最大氣流速度達552 m/s，模型4泄漏2 h點火爆炸，最大氣流速度達799 m/s。高速爆炸風(fēng)產(chǎn)生的沖擊動壓導(dǎo)致餐飲店內(nèi)設(shè)備發(fā)生位移，也會使人體受到拋擲和磕碰致傷。

此外，爆炸過程中急劇燃燒產(chǎn)生高溫。模型3泄漏2 h點火產(chǎn)生的高溫逾2 200 K，模型4泄漏2 h點火產(chǎn)生的高溫逾2 300 K，會使人體受到高溫傷害、燒傷或致死。

5 結(jié)論

① 典型臨街餐飲店內(nèi)發(fā)生天然氣泄漏時，餐飲店體積對天然氣擴散產(chǎn)生影響，在泄漏流量相同的情況下，餐飲店體積越小，天然氣體積分數(shù)越高。

② 相同泄漏源條件下發(fā)生天然氣泄漏，相比餐飲店門關(guān)閉，餐飲店門開啟時天然氣擴散速度更快，導(dǎo)致泄漏初始時段天然氣體積分數(shù)更高。泄漏一段時間后，天然氣通過餐飲店門不斷擴散至餐飲店外部，餐飲店內(nèi)天然氣體積分數(shù)比門關(guān)閉時低。

③ 一定條件下，風(fēng)機開啟后20 s內(nèi)可顯著消除餐飲店內(nèi)的可燃氣云，因此合理設(shè)置送排風(fēng)系統(tǒng)可以有效控制餐飲店內(nèi)天然氣泄漏事故后果。

④ 在餐飲店門關(guān)閉條件下，模型3泄漏1 h后點火，爆炸過程中壓力波破壞餐飲店門，將燃燒產(chǎn)物泄放至外部空間，點火后574 ms產(chǎn)生超壓峰值18 kPa以上，同時燃燒形成的熱量以及燃燒產(chǎn)物膨脹使得計算區(qū)域內(nèi)產(chǎn)生高速氣流，爆炸計算域內(nèi)最大氣流速度達420 m/s。

⑤ 對于同一模型，無論餐飲店門開閉，餐飲店內(nèi)天然氣體積分數(shù)達到爆炸上限以前，點火爆炸產(chǎn)生的超壓峰值隨泄漏時間延長而增大。對于不同體積餐飲店，泄漏剛開始的一段時間內(nèi)，體積較小的餐飲店內(nèi)點火爆炸產(chǎn)生的超壓更大；天然氣泄漏時間較長時，體積較大的餐飲店點火爆炸產(chǎn)生的超壓更大。

Craig原本在金融業(yè)工作，2014年7月他辭職了。離開了原來的崗位后，他橫跨半個地球，搬到了澳大利亞生活。他告訴我們，他對攝影的熱情就是從那時培養(yǎng)起來的，并且從此一發(fā)不可收拾。起初，他主要從事風(fēng)光攝影和旅行攝影的工作。而且他一直熱衷拍攝各地的風(fēng)光和多樣的人文。如今，他的攝影愛好已經(jīng)足以謀生，現(xiàn)在他在Bristol從事肖像攝影和商業(yè)攝影的工作。

⑥ 人體位移傷害概率隨距點火點距離的增加而逐漸減小。高速爆炸風(fēng)產(chǎn)生的動壓會使人體受到拋擲和磕碰致傷。爆炸過程中急劇燃燒產(chǎn)生的高溫會使人體受到高溫傷害、燒傷或致死。

[ 1 ] HUH C，CHOI S，LEE J M. Concentration model for gas releases in buildings and the mitigation effect[J]. Journal of Loss Prevention in the Process Industries，2020，65:104135-1-9.

[ 2 ] LITTO R，HAYES R E，LIU B. Capturing fugitive methane emissions from natural gas compressor buildings[J]. Journal of Environmental Management，2007(3): 347-361.

[ 3 ] 陳琪，陳彩霞，吳亭亭，等. 室內(nèi)天然氣管道微量泄漏氣體擴散特性的CFD模擬[J]. 安全與環(huán)境學(xué)報，2018(6): 2224-2229.

[ 4 ] SUN S，QIU Y，XING H，et al. Effects of concentration and initial turbulence on the vented explosion characteristics of methane-air mixtures[J]. Fuel，2020(1):1-9.

[ 5 ] JI T，QIAN X，YUAN M，et al. Case study of a natural gas explosion in Beijing，China[J]. Journal of Loss Prevention in the Process Industries，2017，49: 401-410.

[ 6 ] PEDERSEN H H，TOMLIN G，MIDDHA P，et al. Modelling large-scale vented gas explosions in a twin-compartment enclosure[J]. Journal of Loss Prevention in the Process Industries，2013(6):1604-1615.

[ 7 ] SEZER H，KRONZ F，AKKERMAN V，et al. Methane-induced explosions in vented enclosures[J]. Journal of Loss Prevention in the Process Industries，2017，48: 199-206.

[ 8 ] LI J，HAO H. Internal and external pressure prediction of vented gas explosion in large rooms by using analytical and CFD methods[J]. Journal of Loss Prevention in the Process Industries，2017，49: 367-381.

[ 9 ] 李紅培. 開放式廚房燃氣泄漏爆炸模擬研究(碩士學(xué)位論文)[D]. 濟南：山東建筑大學(xué)，2019:10-60.

[10] BAALISAMPANG T，ABBASSI R，GARANIYA V，et al. Modelling an integrated impact of fire，explosion and combustion products during transitional events caused by an accidental release of LNG[J]. Process Safety and Environmental Protection，2019，128: 259-272.

[11] 黃小美，彭世尼，徐海東，等. 燃氣管道泄漏流量的計算[J]. 煤氣與熱力，2008(3):B11-B16.

[12] ASSAEL M J，KAKOSIMOS K E. Fires，explosions，and toxic gas dispersions:effects calculation and risk analysis[M]. Boca Raton(US): A Wiley Company，2010：23-25.

[13] TNO. Methods for the determination of possible damage to people and objects resulting from releases of hazardous materials[M]. Voorburg：Director-general of Labour，1992:37-44.