甲烷爆炸對建筑物內(nèi)外壓力場分布的影響

李玉星，尹淵博，王雅真，劉翠偉

（1 中國石油大學(xué)（華東）山東省油氣儲運安全省級重點實驗室，山東青島266580； 2 中國石化青島安全工程研究院，山東青島266580； 3 中國石油大學(xué)（華東）儲運與建筑工程學(xué)院，山東青島266580）

引 言

近年來，管道泄漏爆燃事故頻繁發(fā)生，造成大量人員傷亡和財產(chǎn)損失。當(dāng)泄漏的甲烷氣體濃度處于爆炸極限內(nèi)，便存在爆炸的危險，因此需要首先明確甲烷氣體的爆炸極限。Vanderstraeten 等[1-5]研究發(fā)現(xiàn)，甲烷氣體在20℃，100 kPa 下的爆炸下限為4.6%±0.3%，上限15.8%±0.4%，但這一范圍在不同溫度壓力條件下會發(fā)生變化。Cui 等[6]測量了甲烷-空氣混合物在低初始溫度(113 K 以下)和高壓力下的爆炸特性，燃燒持續(xù)時間和火焰發(fā)展持續(xù)時間均隨初始壓力的增大而增大，隨初始溫度的升高而減小。Huang 等[7]研究表明甲烷在高壓下的燃燒下限略有降低，超過3 MPa 后燃燒上限顯著提高。Zheng 等[8]研究了存在濃度梯度(既非均勻混合)的甲烷-空氣混合物的燃燒，其結(jié)果表明，組分濃度梯度會使燃燒速度降低。

甲烷空氣預(yù)混氣體一旦發(fā)生爆炸，會產(chǎn)生極其嚴重的危害，為了能更好的減小危害程度，必須對其可能產(chǎn)生的后果進行分析。目前對于氣體爆炸超壓后果的計算通常使用TNT 當(dāng)量模型[9]和TNO 多能模型[10]，前者是將蒸氣云爆炸的能量當(dāng)量轉(zhuǎn)換為TNT炸藥的質(zhì)量，但TNT是理想爆源，導(dǎo)致該方法誤差較大。后者是將不同燃燒速率下的無量綱爆炸超壓與無量綱距離繪制成曲線圖查圖使用，使用時受人為的主觀影響很強，結(jié)果可信度并不高。Dobashi 等[11]對此進行了改進，得到了適用于氣體爆炸壓力場的計算公式。為更好的揭示氣體爆炸規(guī)律，部分學(xué)者對爆炸壓力場及超壓沖擊波的傳播做了深入研究：Bao等[12]通過設(shè)置不同體積和濃度的預(yù)混甲烷/空氣氣云爆炸實驗研究了開敞空間內(nèi)甲烷爆炸的相關(guān)特性。Gostintsev 等[13]基于自相似火焰?zhèn)鞑ダ碚?，提出了具有自發(fā)火焰不穩(wěn)定性的火焰?zhèn)鞑グ虢?jīng)驗公式。Kim[14]提出無約束氣體爆炸具有自相似傳播特征，并在研究中對Gostintsev 等[13]提出的公式做出了修正。Strehlow 等[15]利用火焰?zhèn)鞑サ暮愣ㄋ俣群图铀俣葋矸治鲇汕蛐位鹧娈a(chǎn)生的爆炸波。Li 等[16]分析了城市井蓋甲烷氣體爆炸超壓變化規(guī)律。Huang 等[17]實驗研究發(fā)現(xiàn)最大爆炸超壓與最大燃燒速度呈對數(shù)函數(shù)關(guān)系。

上述學(xué)者的研究工作大多是面向無約束條件下的爆炸工況，但實際生產(chǎn)過程中面對的往往是存在障礙物的情況。當(dāng)存在障礙物時，由于在傳播過程中氣體受到阻塞，火焰在通過障礙物時其鋒面扭曲褶皺形成湍流旋渦，部分學(xué)者針對障礙物存在條件下的爆炸進行了研究。Ciccarelli 等[18]研究發(fā)現(xiàn)，障礙物存在情況下沖擊波會限制最終的燃燒速度。Harrison等[19]通過實驗證明，爆炸火焰的傳播速度和超壓隨障礙物間距的增大出現(xiàn)先增加后減小的趨勢。Wang等[20]和Hall等[21]得出相似結(jié)論。不同形狀的障礙物(如圓柱形、三角形、立方體、菱形和平板形等)，都可以加速火焰的傳播[22]，Ibrahim 等[23]分析了不同形狀的障礙物對爆炸結(jié)果的影響，其實驗結(jié)果表明，板形障礙物所產(chǎn)生的超壓最高，而圓柱形障礙物對超壓的影響最低。Xu 等[24]通過實驗同樣得出類似的結(jié)論。

目前對于無約束條件下甲烷空氣混合氣體的燃燒過程以及爆炸超壓的研究已成體系，對于約束條件下，氣體預(yù)混爆炸在障礙物周圍產(chǎn)生壓力場的研究也取得了一定的進展，但目前尚未有對于氣體爆炸在障礙物內(nèi)部產(chǎn)生壓力場的研究，因此本文在經(jīng)實驗數(shù)據(jù)修正的無約束預(yù)混爆炸模型的基礎(chǔ)上，建立三維約束預(yù)混氣體爆炸模型，分析了障礙物高度和寬度對障礙物外部壓力場分布的影響，與之前相關(guān)文獻作比較驗證該模型與求解方法的正確性，在此基礎(chǔ)上研究了氣體預(yù)混爆炸在障礙物內(nèi)部產(chǎn)生壓力場的分布規(guī)律，給出了實際爆炸發(fā)生時相對安全的位置。

1 理論分析

1.1 傷害判定準(zhǔn)則

目前較為常用的爆炸傷害判斷標(biāo)準(zhǔn)主要有超壓準(zhǔn)則、沖量準(zhǔn)則以及二者結(jié)合使用的方法。超壓準(zhǔn)則認為，只有當(dāng)爆炸沖擊波超壓大于或等于某一臨界值時，才會對目標(biāo)造成一定的傷害。否則，爆炸沖擊波不會對目標(biāo)造成傷害。沖量準(zhǔn)則認為，只有當(dāng)作用于目標(biāo)的爆炸波沖量達到某一臨界值時才會引起傷害[25]。實踐表明，不同的爆炸源有不同的爆炸波形，同樣沖量值產(chǎn)生的破壞作用也可能會顯著不同[26]。文章選用超壓準(zhǔn)則作為衡量標(biāo)準(zhǔn)展開，表1 和表2 展示了超壓波對人體的傷害和對建筑物的破壞。

表1 沖擊波超壓對人體的傷害作用Table 1 Damage effect of shock wave overpressure on human body

表2 沖擊波超壓對建筑物的破壞作用Table 2 Damage effect of shock wave overpressure on buildings

1.2 方程及模型選擇

模型所利用到的基本空置方程有組分守恒方程、連續(xù)性方程、動量守恒方程以及能量守恒方程；經(jīng)過對不同湍流模型，最終發(fā)現(xiàn)使用SAS 湍流模型得出的結(jié)果更符合實際，SAS 模型能夠自動調(diào)整解析更小尺度的渦；燃燒模型分別用到了層流有限速率模型和渦-耗散-概念(EDC)模型。無約束爆炸過程中，氣流的運動不受阻礙，其湍流程度較低，層流有限速率燃燒模型既可以滿足計算精度又可以節(jié)約計算資源，但這一模型無法描述約束爆炸狀況。渦-耗散-概念(EDC)模型在湍流的流動結(jié)構(gòu)中納入詳細的化學(xué)反應(yīng)機理。它假設(shè)化學(xué)反應(yīng)均發(fā)生在小渦當(dāng)中，稱為精細渦。反應(yīng)時間由小渦生存時間和化學(xué)反應(yīng)本身所需要的時間來共同控制。其可以更好地描述湍流程度較高地燃燒過程。

在Fluent 中，假設(shè)小渦中的化學(xué)反應(yīng)發(fā)生在相同條件下，初始條件取為單元中當(dāng)前的組分和溫度。組分守恒方程中的化學(xué)反應(yīng)源項由式(1)得出。

式中，Ri為化學(xué)反應(yīng)源項；ε*為小渦的尺度；ρ 為密度，kg/m3;Y*為經(jīng)過一段時間的反應(yīng)后，小尺度組分的質(zhì)量分數(shù)；τ*為反應(yīng)經(jīng)過的時間尺度；Y 為同樣的物質(zhì)在反應(yīng)之前的濃度;上角標(biāo)*表示小尺度；下角標(biāo)i表示第i種物質(zhì)。

2 仿真模型建立與實驗驗證

首先建立無約束超壓沖擊波傳播模型，并使用實驗數(shù)據(jù)進行修正，得到一種可行的超壓傳播仿真模型，然后通過分析有無約束對爆炸超壓沖擊波傳播過程的影響，對無約束爆炸模型進行修正并應(yīng)用于約束條件下的爆炸過程，分析障礙物周圍及障礙物內(nèi)部在預(yù)混氣體爆炸時壓力場的分布。

2.1 無約束爆炸仿真計算

2.1.1 無約束爆炸仿真模型建立 足夠的網(wǎng)格密度是保證Fluent 計算精度的關(guān)鍵，但是受計算機計算能力和所需計算時間的限制，無法對網(wǎng)格進行無限加密。因此需要對不同網(wǎng)格密度進行計算分析，找到合理的網(wǎng)格密度，以確保計算效率[27]。本研究共設(shè)置三個不同網(wǎng)格密度，分別編為1#、2#、3#，具體的網(wǎng)格信息如表3所示，三種網(wǎng)格的質(zhì)量均在0.6以上，能夠較好的滿足Fluent的計算要求。

對邊長為2 m 的近地面甲烷/空氣預(yù)混爆炸進行數(shù)值模擬，使用中心點火方式，檢測距點火點水平位置為2 m處的超壓，得到的結(jié)果如圖1所示。

表3 網(wǎng)格信息Table3 Grid information

圖1 網(wǎng)格敏感性分析結(jié)果Fig.1 Sensitivity analysis of grid

由圖1 可見，2#與3#的計算結(jié)果，對于正超壓峰值和負超壓峰值相對誤差均小于5%，因此2#的網(wǎng)格密度比較適宜，能夠在保證計算結(jié)果可靠性的基礎(chǔ)上節(jié)約計算資源和計算時間，故本次研究采用2#的網(wǎng)格密度。

使用甲烷單步反應(yīng)方程式，其具體形式如下：

該化學(xué)模型的指前因子取值為1.35×1020，活化能為1.225×108J/kmol[28]。

建立一邊長為2 m 的甲烷/空氣預(yù)混立方體爆炸模型，預(yù)混氣體濃度為化學(xué)計量濃度9.5%(體積)，中心點火。計算監(jiān)測點為地面上距爆心水平距離為2 m 的位置?？紤]到無約束爆炸的對稱性以及節(jié)約計算資源，首先建立無約束爆炸二維仿真模型，并對中心點火區(qū)域另行加密，網(wǎng)格如圖2所示。

圖2 無約束爆炸模型網(wǎng)格示意圖Fig.2 Schematic diagram of unconstrained explosion model grid

考慮計算精度及計算機資源占用，時間步長取1×10-5s。在無約束爆炸中，由于沒有障礙物引起火焰鋒面的褶皺，故而燃燒反應(yīng)的湍流程度較低，因此層流燃燒模型更適合無約束爆炸工況[29-30]。邊界條件使用符合實際的壓力出口，采用PISO 算法，二階迎風(fēng)格式。

為選擇適合的湍流模型，分別對k-ε湍流模型、k-ε 湍流模型以及SAS 湍流模型進行模擬分析，使用不同湍流模型進行計算。表4總結(jié)了不同湍流模型的峰值超壓及峰值超壓到達時間。對比后續(xù)2.2節(jié)中實驗數(shù)據(jù)，可知SAS 湍流模型的計算結(jié)果更為接近實驗數(shù)據(jù)，因此選擇SAS 湍流模型來進行后續(xù)數(shù)值模擬。

表4 不同湍流模型峰值超壓與到達時間對比Table 4 Comparison of peak overpressure and arrival time of different turbulence models

采用Patch(定義某一局部條件)高溫區(qū)的方式進行點火，模型點火溫度選擇為能夠點著的最低溫度1500 K，使用式(2)計算以1500 K 高溫點火的點火能量，當(dāng)點火區(qū)域半徑為0.004 m 時，點火能量為100 mJ。

式中，Q 是點火能量，mJ；c 是比熱容，kJ/(kg·mol)；m是被點火區(qū)域氣體質(zhì)量，kg；T是溫度，℃。

保持點火溫度1500 K 不變，改變點火體積來改變點火能量，數(shù)值計算研究了球體點火區(qū)域半徑分別為0.004、0.006、0.008 和0.01 m 時點火能量對爆炸超壓的影響，結(jié)果列于表5中。由結(jié)果可知，在弱點火能情況下，點火能量的少許改變對爆炸超壓的影響并不大，因此后續(xù)模擬采用0.004 m 的點火球體半徑。

表5 點火能量對爆炸超壓的影響Table 5 Effect of ignition energy on explosion overpressure

2.1.2 無約束爆炸峰值超壓模擬計算 針對無約束條件下，不同體積的化學(xué)計量濃度下甲烷/空氣預(yù)混氣云，采用上述方法進行模擬。共設(shè)置3 個不同體積的立方形預(yù)混氣體爆炸工況，具體參數(shù)如表6所示。在爆炸中心線上，距爆炸中心點1、2、2.5、3和4 m處設(shè)置五個壓力監(jiān)測點。

表6 不同體積立方形預(yù)混氣體爆炸工況條件Table 6 Explosion conditions of different volume cubic premixed gases

2.2 無約束爆炸峰值超壓實驗驗證

由于開展甲烷氣體爆炸實驗成本高，難度大，因此本文使用Bao 等[12]的實驗數(shù)據(jù)進行驗證。Bao等[12]建立體積為V=1,8,27 m3的立方體形甲烷/空氣預(yù)混氣體進行了中心點火的無約束爆炸實驗，預(yù)混氣體滿足甲烷燃燒的化學(xué)計量濃度[甲烷濃度c=9.5%(體積)]，點火能量為100 mJ。實驗使用聚乙烯薄膜對預(yù)混氣體進行固定，甲烷的預(yù)混及穩(wěn)定均在薄膜內(nèi)進行。實驗驗證聚乙烯薄膜會在超壓產(chǎn)生瞬間破裂，對爆炸超壓的影響微乎其微，可以忽略不計。壓力監(jiān)測點設(shè)置在地面上與點火點水平距離分別為d=2,2.5,3 m 的位置。具體的實驗裝置如圖3 所示。圖4 為距離點火點水平距離分別為d=2,2.5,3 m 的地面處的超壓隨時間的變化。

選取模擬過程中，預(yù)混氣體邊長為2 m 的工況與實驗數(shù)據(jù)進行對比，模擬工況及監(jiān)測點位置與實驗相同，其模型示意圖如圖5 所示，圖6 繪制了預(yù)混氣體邊長為2 m 的工況在三個監(jiān)測點處的超壓隨時間變化圖。

將圖6 與圖4 中對應(yīng)實驗數(shù)據(jù)進行比較可見，F(xiàn)luent 數(shù)值計算結(jié)果一直處于波動狀態(tài)，但總體趨勢與實驗結(jié)果相同。對于爆炸正負峰值超壓的到達時間，F(xiàn)luent 模擬結(jié)果要略小于實驗數(shù)據(jù)，但在實際生產(chǎn)中，零點幾秒的時間可以忽略，不會產(chǎn)生影響。對于峰值數(shù)值差異較大的現(xiàn)象，提出如下假設(shè)：數(shù)值模擬中預(yù)混氣云內(nèi)所有的甲烷均發(fā)生了完全燃燒，而實際實驗工況中并非所有可燃物均能夠發(fā)生完全燃燒，部分可燃物發(fā)生不完全燃燒反應(yīng)，反應(yīng)放出的熱量小于模擬計算值；同時由于爆炸風(fēng)的吹散，使得壓力波附近的可燃氣云被壓縮，甲烷濃度升高，而外圍還未反應(yīng)的氣云向周圍擴撒，部分氣云與空氣再次發(fā)生混合稀釋，甲烷濃度不再是化學(xué)計量濃度，甚至超出爆炸極限范圍，故而實際爆炸工況產(chǎn)生的超壓強度要低于數(shù)值模擬結(jié)果。下面將證明此種假設(shè)的正確性。

圖7 繪制了預(yù)混氣云邊長2 m 工況中0.15、0.2和0.22 s時對應(yīng)的甲烷體積濃度云圖。

圖3 實驗裝置圖Fig.3 Diagram of experimental equipment

圖4 無約束甲烷/空氣預(yù)混爆炸在不同混合物體積、不同檢測距離的實驗數(shù)據(jù)結(jié)果Fig.4 Experimental results for unconstrained methane/air premixed explosions in different mixture volumes and different detection distances

圖5 無約束爆炸模型示意圖Fig.5 Schematic diagram of unconstrained explosion

由圖7 可見，甲烷的反應(yīng)一直在預(yù)混范圍內(nèi)進行，未曾發(fā)生氣云的吹散。研究同時監(jiān)測了距離爆心水平距離為1.1、1.3、1.5 m 處三點的甲烷濃度，在反應(yīng)過程中其甲烷濃度一直保持為零。以上分析印證了Fluent計算模擬值較大原因的兩條假設(shè)。

圖6 邊長2 m無約束預(yù)混立方氣體爆炸Fluent數(shù)值模擬結(jié)果Fig.6 Fluent numerical simulation results of 2 m unconstrained premixed cubic gas explosion

為了直觀地分析Fluent模擬計算值與實驗值關(guān)系，表7列出了原始模擬計算值、實驗值以及兩者的比值關(guān)系。正負峰值超壓的Fluent計算值與實驗值的比值均在2.3～3.5之間，大部分比值都接近3，所以這里選擇比值關(guān)系3 作為修正系數(shù)，即將Fluent 的原始模擬計算值乘以系數(shù)1/3 作為修正結(jié)果進行后續(xù)分析，且稱修正后體積為當(dāng)量體積。

圖7 2 m工況對應(yīng)甲烷體積濃度分布圖Fig.7 Methane volume concentration distribution of case 2 m

表7 Fluent模型正負超壓原始模擬計算值與實驗值對比Table 7 The positive and negative overpressure original Fluent simulation results and experimental data

2.3 約束爆炸仿真模型建立

無約束條件下發(fā)生的預(yù)混爆炸，由于沒有障礙物引起火焰鋒面的褶皺，故而燃燒反應(yīng)的湍流程度較低，而當(dāng)氣體爆炸影響區(qū)域內(nèi)存在障礙物時，由于在傳播過程中氣體受到阻塞，火焰在通過障礙物時其鋒面扭曲褶皺形成湍流旋渦，層流燃燒模型將不再適用。渦流耗散模型(EDC)能夠在湍流反應(yīng)流動中考慮詳細的化學(xué)反應(yīng)機理，計算精度高，可以用于描述由于障礙物存在產(chǎn)生的湍流旋渦。除此之外，障礙物存在不會改變預(yù)混爆炸的其他過程，因此，在無約束工況甲烷預(yù)混爆炸數(shù)值模型的基礎(chǔ)上，將層流燃燒模型換為EDC 模型，即可用來完成甲烷預(yù)混氣體約束爆炸模擬分析。

在上述研究的基礎(chǔ)上，使用選定的計算方法，建立約束條件下甲烷/空氣預(yù)混爆炸的三維仿真模型：立方預(yù)混氣體體積8 m3，處于化學(xué)計量濃度，距地面高度1 m；監(jiān)測區(qū)域為20 m×20 m×15 m 的長方體區(qū)域；障礙物長寬高為1 m×2 m×4 m，距爆心2 m，如圖8 所示。監(jiān)測點設(shè)置在障礙物的三個位置處：與爆心高度相同的正對爆炸面及背面，位于地面上的正對爆炸面及背面，以及頂部中心線處，具體分布如圖9 所示。網(wǎng)格尺寸與無約束工況相同，選為空氣區(qū)域0.2 m×0.2 m×0.2 m，可燃物區(qū)域0.02 m×0.02 m×0.02 m，點火區(qū)域加密0.005 m×0.005 m×0.005 m，在障礙物壁面附近同樣進行了網(wǎng)格加密，其他設(shè)置不變，建立用于計算約束爆炸的仿真模型。

圖8 障礙物的布置(a)與網(wǎng)格劃分(b)Fig.8 Schematic diagram of obstacle arrangement(a)and meshing(b)

圖9 監(jiān)測點布置示意圖Fig.9 Schematic diagram of monitoring point

3 障礙物存在工況下甲烷預(yù)混爆炸數(shù)值分析

主要從障礙物的高度、寬度兩個方面考慮，以數(shù)值模擬的方法來研究氣云爆炸在障礙物周圍壓力場的分布，并進一步探究了障礙物內(nèi)部的壓力場。

3.1 數(shù)值模擬工況設(shè)置

按照上述模型建立方法，通過改變障礙物幾何特征，得到表8 所示幾種仿真工況，分別研究了障礙物的高度和寬度對障礙物周圍壓力場分布的影響。

3.2 障礙物高度的影響結(jié)果分析

共設(shè)置5 個研究工況A-1～A-5，研究障礙物高度對爆炸后果的影響。障礙物的寬度均為2 m，厚度為1 m，障礙物正面到爆心的水平距離為2 m，障礙物的高度h 從2 m 增加到6 m，應(yīng)用2.3 節(jié)建立的模型進行計算，并對結(jié)果進行1/3 修正。因A-3、B-2 條件設(shè)置相同，選A-3 作為參比工況。圖10 展示了A-3(h=4 m)工況不同時刻通過爆心的XZ 平面壓力場分布圖。

由圖10 可見，當(dāng)預(yù)混氣云附近存在障礙物時，壓力場的分布不再是對稱的。正超壓場明顯偏向于障礙物側(cè)，負超壓場則偏向于無障礙物側(cè)，持續(xù)時間更長，但障礙物側(cè)的負超壓絕對值更大。障礙物的存在對正超壓的傳播起到了一定的阻擋作用，障礙物后的正超壓值相對較低，但后部會受到負壓的強烈沖擊，這在無障礙物側(cè)距離相等的位置是不存在的。爆炸產(chǎn)生的壓力場一直處在正負超壓的振蕩之中，短時間內(nèi)并不會回到零超壓。

表8 障礙物研究工況條件設(shè)置表Table 8 Obstacle research conditions setting table

圖10 A-3工況不同時刻通過爆心的XZ平面壓力場分布圖Fig.10 Pressure distribution of XZ plane passing through the ignition point at different times for case A-3

對于A 組障礙物工況，分別在障礙物正對爆炸面和背面的半側(cè)寬度設(shè)置監(jiān)測點，高度兩米處(與爆炸中心高度相同)4 個和地面處4 個。圖11 繪制了障礙物正對爆炸面上的超壓分布，實線為與爆炸中心同高度處的超壓分布，虛線為地面處的超壓分布。

由圖11 可知，從障礙物的中軸線向兩側(cè)方向正負超壓絕對值有所下降。對比實線簇與虛線簇，總體上地面處的超壓要大于與爆心水平高度處，這是由于地面與障礙物圍成的夾角阻礙了氣體的傳波，極大地增強了氣體的湍流強度，壓力波在拐角處的多次折射也增大了總體的壓力強度，造成地面處的超壓升高。這表明，在實際遇到爆炸事故時，不要盲目地從建筑物上層向下層疏散，障礙物下層受到正負超壓的反復(fù)沖擊，極易對人體造成損害[31]。

如圖12 所示，在障礙物背對爆炸源的一面，隨著障礙物高度的增加，障礙物后側(cè)的正峰值超壓減小，但負峰值超壓不斷增大。這說明障礙物的背面有一定的抵御爆炸正超壓的能力。同時，當(dāng)障礙物高度升高到一定值，也即高寬比較大時(h=6 m 工況)，障礙物造成的湍流對超壓的增強作于強于其抵御作用，障礙物后側(cè)超壓將會急劇增大。對比前后兩側(cè)的壓差可以得到同樣的特性，隨著障礙物高度的增加，兩側(cè)壓差逐漸增大，但當(dāng)障礙物高寬比較大時，壓差又出現(xiàn)減小的趨勢。這是由于正對爆炸面在地面與障礙物形成的拐角處由湍流增大的超壓更加容易繞過障礙物傳播至其后側(cè)空間。與障礙物正面的超壓場相同，背面同樣呈現(xiàn)地面處的正負超壓絕對值均大于爆心水平面的超壓值。綜合來看，障礙物背對爆炸面的中上部才是最安全的位置，在實際爆炸工況中不應(yīng)盲目向建筑物下部撤離。

圖11 A組障礙物正對爆炸面超壓分布曲線Fig.11 Overpressure distribution on the facing surface of obstacles of case group A

圖12 A組障礙物背對爆炸面超壓分布曲線Fig.12 Overpressure distribution on the back surface of obstacles of case group A

3.3 障礙物寬度的影響結(jié)果分析

設(shè)置5 個研究工況編號B-1～B-5，研究障礙物寬度對爆炸后果的影響。障礙物的高度均為4 m，厚度為1 m，障礙物正面到爆心的水平距離為2 m，障礙物的寬度b 從1 m 增加到5 m，使用與3.2 節(jié)相同的計算方法。為了更加直觀地分析障礙物寬度對爆炸超壓的影響，將各個工況最大正負超壓對應(yīng)時刻過爆心的XZ 平面壓力場繪制于圖13 和圖14中。

由圖可知，隨著障礙物寬度的增加，爆炸的影響范圍不斷增大。當(dāng)障礙物寬度增大時，地面與障礙物的拐角處氣體聚集并發(fā)生湍流增強作用，正超壓在拐角處急劇增大。障礙物對爆炸正超壓具有較好的抵御作用。負超壓的絕對值在障礙物寬度增加后同樣逐步升高，負壓場出現(xiàn)不對稱性，障礙物一側(cè)的負壓范圍增大。并且障礙物寬度增大后，障礙物后側(cè)出現(xiàn)大范圍的高強負壓場，說明障礙物對負壓的抵御能力較弱。

對B 組障礙物寬度工況，分別在障礙物正對爆炸面和背面的半側(cè)寬度上設(shè)置監(jiān)測點，監(jiān)測點高度分別設(shè)在2 m 高度處(與爆炸中心高度相同)和地面處。工況B-1、B-2 同一高度設(shè)置4 個監(jiān)測點，工況B-3、B-4、B-5 在同一高度設(shè)置5 個監(jiān)測點。圖15繪制了障礙物正對爆炸源面的超壓分布，實線為與爆炸中心同高度處的超壓分布，虛線為地面處的超壓分布。

圖13 B組工況爆炸最大正超壓對應(yīng)時刻壓力場Fig.13 Pressure fields at maximum positive overpressure corresponding time for case group B

圖14 B組工況爆炸最大負超壓對應(yīng)時刻壓力場Fig.14 Pressure fields at maximum negative overpressure corresponding time for case group B

圖15 B組障礙物正對爆炸面超壓分布曲線Fig.15 Overpressure distribution on the facing surface of obstacles of case group B

對比實線簇與虛線簇，總體上在障礙物正對爆炸的一側(cè)，地面附近無論正負超壓，其絕對值總是要大于上方與爆心同高度的位置，與高度工況相同，說明在地面與障礙物形成的拐角處可燃氣發(fā)生積聚，在增加湍流程度的同時增大了爆炸超壓。對于負壓，隨著障礙物寬度的增加，負超壓的絕對值呈現(xiàn)先上升后下降的趨勢，同時上下不同位置的壓差逐漸減小，說明隨著障礙物寬度的增大，對負壓具有削弱作用。在障礙物寬度方向上，正負超壓絕對值呈現(xiàn)出由中軸線向兩側(cè)遞減的趨勢。

如圖16 所示，在障礙物背對爆炸源的一面，隨著障礙物寬度的增加，障礙物后側(cè)的正峰值超壓不斷減小，但負峰值超壓總體上呈現(xiàn)小范圍的波動。

對比同一工況障礙物前后兩面的正超壓值可以發(fā)現(xiàn)壓差不斷增大，這說明障礙物寬度的增加，使得其抵御爆炸正超壓的能力不斷提高，能夠有效減輕正超壓對障礙物后方的破壞。但對于負超壓來說，兩側(cè)壓差并沒有明顯的變化，障礙物寬度的增加對其抵御負超壓的能力并無顯著效果，障礙物后方主要承受負壓破壞。與障礙物正面的超壓場相同，背面同樣呈現(xiàn)地面處的正負超壓絕對值均大于爆心水平面的超壓值，負壓強度要高于正壓。綜合來看，障礙物寬度是其抵御爆炸正超壓能力的重要因素。

3.4 障礙物內(nèi)部壓力場分布

圖16 B組障礙物背對爆炸源面超壓分布曲線Fig.16 Overpressure distribution on the back surface of obstacles of case group B

圖17 障礙物內(nèi)部監(jiān)測點布置Fig.17 The monitor points arrangement in obstacle

3.2節(jié)、3.3節(jié)中均假設(shè)障礙物是全封閉的，能夠?qū)?nèi)部起到有效的保護，但實際工況中的建筑物并非完全封閉。為了研究障礙物內(nèi)部的壓力場，設(shè)立工況C，障礙物尺寸與標(biāo)準(zhǔn)工況中相同(寬2 m，厚1 m，高4 m)，距離爆炸中心的距離為2 m，殼體厚度0.05 m。在障礙物下部開口0.6 m×1 m，模擬建筑物大門。障礙物內(nèi)部監(jiān)測點設(shè)置在中心YZ平面上，具體布置如圖17所示。

經(jīng)過Fluent 數(shù)值模擬計算，將封閉障礙物工況(A-3)和開門障礙物工況C 在爆炸過程中的最大正負超壓值列于表9中。當(dāng)障礙物存在一個下部開口時，對整個爆炸的最大正超壓影響并不大，但會明顯降低爆炸負超壓的最大絕對值。

在研究障礙物內(nèi)部壓力場之前，首先對障礙物下部開口后對其外部超壓場的影響進行分析。圖18 繪制了兩個工況障礙物正對爆炸面半側(cè)寬度上由中心向一側(cè)超壓的變化，其中實線代表z=2 m 高度上(與爆心同高)的超壓，虛線代表地面上(z=0 m)的超壓。

表9 工況A-3與工況C最大正負超壓對比Table 9 Maximum overpressure comparison between case A-3 and case C

由圖可見，當(dāng)障礙物下部開門后，爆心水平面上的正超壓略有降低，但并不明顯。地面處的正超壓相較于封閉式障礙物大大降低，這是因為原來積聚在此處的氣體能夠向障礙物內(nèi)部流動，大大降低了阻塞率。障礙物下部開口后，正側(cè)的負超壓變化并不明顯。

圖19 為兩個工況障礙物背對爆炸面半側(cè)寬度上由中心向一側(cè)超壓的變化，同樣實線代表z=2 m高度上(與爆心同高)的超壓，虛線代表地面上(z=0 m)的超壓。障礙物下部開門后，對于爆心水平面上的正負超壓均無明顯變化，但對于地面處，障礙物后側(cè)的正負超壓絕對值均顯著增大。

在此基礎(chǔ)上，對障礙物內(nèi)部的壓力場進行分析，圖20 為與爆炸中心同高度2 m 處，障礙物內(nèi)部半側(cè)寬度上由中軸線向一側(cè)的超壓變化，圖21則為障礙物中軸線上，由地面向上正負超壓的變化。

由障礙物中心至兩側(cè)正負超壓有少許下降，但并不明顯，總體上正負超壓的絕對值大致與障礙物背面相同，說明障礙物對內(nèi)部有較好的保護作用。隨著高度的增加，正負超壓的絕對值均明顯下降，正超壓的下降速度大于負超壓。

圖18 工況A-3和工況C障礙物正對爆炸側(cè)正負超壓分布Fig.18 Overpressure distribution on the facing surface of obstacles for case A-3 and case C

圖19 工況A-3和工況C障礙物背對爆炸側(cè)正負超壓分布Fig.19 Overpressure distribution on the back surface of obstacles for case A-3 and case C

圖20 障礙物內(nèi)部z=2 m半側(cè)寬度上正負超壓的變化Fig.20 Overpressure change diagram inside the obstacle on z=2 m height

總體來看，在建筑物不失效的前提下，即建筑物不坍塌能夠?qū)?nèi)部進行有效保護的情況下，越靠近障礙物的上層與兩側(cè)越安全。在緊急爆炸事故中，盲目的向下逃生，反而容易受到正負超壓的強烈沖擊。

圖21 障礙物內(nèi)部中軸線上正負超壓隨高度的變化Fig.21 Overpressure changes with height on the inner axis of the obstacle

4 結(jié) 論

文章針對天然氣泄漏甲烷預(yù)混爆炸的無約束工況，建立了峰值超壓模型，并利用實驗加以修正，分析對比有無障礙物時氣體預(yù)混爆炸的區(qū)別，將所建立模型修正后應(yīng)用于障礙物存在時的爆炸情形，分析了氣體爆炸時障礙物周圍及內(nèi)部的壓力場分布，主要得到以下結(jié)論。

(1)經(jīng)過對Fluent不同模型的分析以及與實驗數(shù)據(jù)的對比，SAS湍流模型、層流燃燒模型最為適合甲烷無約束預(yù)混爆炸，由于湍流程度不同，EDC 燃燒模型更為適用于約束爆炸情形。

(2)在弱點火能情況下，點火能量的少許改變對爆炸超壓的影響并不大。處于爆炸外圍的氣體會在沖擊波的作用下進一步與空氣發(fā)生混合，其濃度可以降低到爆炸下限之外，而壓力波附近的可燃氣云被壓縮，甲烷濃度升高，這部分甲烷將不完全燃燒甚至不燃燒。按照文章給出的仿真模型計算，模擬值與實驗值峰值之比近似為3。

(3)障礙物的存在對正超壓的傳播起到了一定的阻擋作用，障礙物后的正超壓值相對較低，但會受到負壓的強烈沖擊。隨著障礙物寬度的增加，其抵御正超壓的能力不斷提高。在障礙物豎直方向上，貼近地面處受地面與障礙物圍成的夾角的影響，湍流強度及超壓峰值會極大增強。

(4)綜合來看，背對爆炸面的中上部是最安全的位置，實際發(fā)生爆炸時，不要盲目自上而下撤離。在建筑物不失效的前提下，建筑物內(nèi)部越靠近障礙物的上層與兩側(cè)超壓值越小，更為安全。

符 號 說 明

c——比熱容，kJ/(kg·mol)

m——被點火區(qū)域氣體質(zhì)量，kg

Q——點火能量，mJ

T——溫度，℃

ρ——密度，kg/m3

下角標(biāo)

i——參與反應(yīng)的第i種物質(zhì)，i=1,2,3,…