王建 王冬 張培理 梁建軍 劉沖 劉俊家

摘 ?????要：多分支結(jié)構(gòu)坑道在石油化工、天然氣工業(yè)中應(yīng)用十分廣泛。研究多分支結(jié)構(gòu)坑道內(nèi)的油氣爆炸過程，對(duì)避免或減小人員傷亡和財(cái)產(chǎn)損失、提高多分支結(jié)構(gòu)坑道的安全防護(hù)水平具有重要意義。針對(duì)三種多分支結(jié)構(gòu)坑道（一字分布型、交錯(cuò)分布型、相對(duì)分布型），采用 WALE 湍流模型和Zimout 預(yù)混燃燒模型對(duì)上述三種坑道內(nèi)的油氣爆炸過程進(jìn)行了模擬分析，重點(diǎn)研究了火焰?zhèn)鞑ニ俣取⒒鹧婷娣e、超壓峰值、超壓上升速率等爆炸過程的關(guān)鍵參數(shù)。數(shù)值模擬結(jié)果表明：總?cè)莘e不變時(shí)，一字分布型和交錯(cuò)分布型坑道內(nèi)的火焰速度、升壓速率均大于相對(duì)分布型坑道;支坑道數(shù)量不變時(shí)，一字分布型和交錯(cuò)分布型火焰?zhèn)鞑ヒ?guī)律相似，但交錯(cuò)型分布坑道內(nèi)較一字型分布作用更強(qiáng);數(shù)值分析結(jié)果表明，火焰?zhèn)鞑ニ俣取⒒鹧婷娣e與爆炸超壓特性關(guān)系密切，三者之間相互耦合、相互激勵(lì)，存在顯著的正反饋機(jī)制。

關(guān) ?鍵 ?詞：分布形式; 多分支結(jié)構(gòu);密閉空間;油氣;大渦模擬

中圖分類號(hào)：TE88 ???????文獻(xiàn)標(biāo)識(shí)碼： A ???????文章編號(hào)：1671-0460（2019）08-1811-05

Abstract： Multi-branched structure tunnels are widely used in petrochemical and natural gas industries. Studying the oil and gas explosion process in multi-branch tunnels is of great significance to avoid or reduce casualties and property losses， and to improve the safety level of multi-branch tunnels. In this paper， the WALE turbulence model and the Zimont premixed combustion model were used to simulate the oil and gas explosion process in three kinds of multi-branched structure tunnels. The key parameters of the explosion process were studied， such as flame propagation speed， flame area， overpressure peak value， overpressure rise rate and so on. The numerical simulation results showed that： when the total volume was constant， the flame propagation speed， overpressure peak and boost rate in the one-line distribution and staggered-type tunnels were larger than those of the relative distribution tunnels; when the number of branch tunnels was unchanged， the flame propagation laws of one-line distribution and staggered-type tunnels were similar， but the staggered distribution tunnel had stronger effect than the one-line distribution tunnel. The numerical analysis results showed that the relationship among flame propagation speed， flame area and explosion overpressure characteristics was close， they were mutually coupled and mutually motivated， and there was a significant positive feedback mechanism.

Key words： Distribution form; Multi-branch structure; Confined space; Oil and gas; Large eddy simulation

在石油化工、天然氣工業(yè)中，多分支結(jié)構(gòu)坑道（包括坑道、 巷道等）應(yīng)用較為廣泛，比如坑道工程、城市天然氣管網(wǎng)、煤礦開采巷道等。 油氣爆炸發(fā)生在具有多分支結(jié)構(gòu)的受限空間中時(shí)，多分支結(jié)構(gòu)會(huì)對(duì)燃燒爆炸過程中的熱輸運(yùn)和組分輸運(yùn)等過程產(chǎn)生較大影響，與直管或球形容器中的油氣爆炸不同，爆炸后爆炸波會(huì)往復(fù)疊加，壓力會(huì)明顯上升，在傳播過程中很可能發(fā)展成爆轟波，造成了重大人員傷亡和財(cái)產(chǎn)損失[1]。例如，2017年7月，貴州省青龍沙區(qū)天然氣隧道泄漏爆炸，造成8人死亡，35人受傷; 同年7月，天然氣在吉林省松原市寧江區(qū)繁華路發(fā)生泄漏爆炸。因此，如何控制多分支結(jié)構(gòu)坑道中油氣爆燃事故的發(fā)生，降低油氣爆燃的威脅，減少可能造成的人員傷亡和財(cái)產(chǎn)損失，是當(dāng)前亟待解決的主要課題。

相關(guān)領(lǐng)域的專家學(xué)者對(duì)涉及分支結(jié)構(gòu)的密閉空間可燃?xì)怏w爆炸傳播過程進(jìn)行了大量研究，取得了一定的成果。王漢良等[2]研究了L型管道，楊志等[3]研究了Z型管道，ZHAI、ZHU等[4，5]研究了U型管道，杜揚(yáng)、李國(guó)慶、蔣新生等[6-9]研究了T型管道，研究發(fā)現(xiàn)，上述涉分支結(jié)構(gòu)均對(duì)爆炸強(qiáng)度、壓力有強(qiáng)化作用，對(duì)火焰?zhèn)鞑ニ俣扔忻黠@加強(qiáng)，結(jié)果表明了分支空間結(jié)構(gòu)與可燃?xì)怏w爆炸特性密切相關(guān)。

上述涉分支結(jié)構(gòu)容器內(nèi)可燃?xì)怏w爆炸研究，雖然部分涉及了分支坑道，但是，涉及的分支坑道數(shù)量較少、形式較單一，研究的系統(tǒng)性和深度也不足以推廣至多分支結(jié)構(gòu)坑道內(nèi)的油氣爆炸環(huán)境。因此，在實(shí)驗(yàn)的基礎(chǔ)上，建立了多分支結(jié)構(gòu)封閉坑道油氣爆炸的數(shù)值模擬模型，深入研究了油氣爆炸燃燒特性。這項(xiàng)研究將有助于探索多分支結(jié)構(gòu)坑道中油氣的爆炸、燃燒和傳播過程規(guī)律，極大地豐富了油氣爆炸燃燒領(lǐng)域的理論體系，為相關(guān)領(lǐng)域提供了理論參考。

1 ?數(shù)值模擬模型

1.1 ?基本假設(shè)

在具有復(fù)雜分支結(jié)構(gòu)坑道內(nèi)的油氣爆炸燃燒過程是一個(gè)具有強(qiáng)湍流場(chǎng)特征和燃燒反應(yīng)劇烈的過程，包括熱傳遞和組分傳輸、流體流動(dòng)、化學(xué)反應(yīng)和許多其他影響因素。由于上述因素影響，難以建立數(shù)值模擬模型。為了保證數(shù)值模擬準(zhǔn)確可行，在理論模型的構(gòu)建中，本文做以下假設(shè)和簡(jiǎn)化[10]：

（1）封閉坑道內(nèi)各組分氣體均為理想氣體，滿足理想氣體狀態(tài)方程;

（2）油氣燃燒反應(yīng)涉及到復(fù)雜的多步基元反應(yīng)，由Lawrence國(guó)家實(shí)驗(yàn)室發(fā)現(xiàn)的C1～C4詳細(xì)反應(yīng)機(jī)理包括155種成分、689個(gè)基元反應(yīng)[11]。但是為了研究該過程中壓力和火焰?zhèn)鞑ミ^程中的變化，需要根據(jù)烴類化合物燃燒模型的研究將燃燒反應(yīng)簡(jiǎn)化為成兩步反應(yīng)：

（3）充斥密閉空間內(nèi)的油氣混合物處于常溫常壓狀態(tài)，呈均勻分布;

（4）將壁面邊界視為剛性、絕熱、光滑和密閉的。

1.2 ?湍流模型

目前，求解可燃?xì)怏w爆炸問題的數(shù)值方法主要有有直接模擬法（DNS）、雷諾平均法（RANS）和大渦模擬法（LES）[12-14]。在針對(duì)可燃?xì)怏w爆炸的數(shù)值模擬計(jì)算中，RANS模型和 LES 模型都被廣泛采用，為更精確爆燃過程和精細(xì)化流場(chǎng)特征，項(xiàng)目組前期對(duì)fluent軟件的四種常用湍流模型（WALE模型、Dynamic magorinsky-Lilly 模型、RNG k-ε模型和 Reynolds Stress 模型）進(jìn)行試運(yùn)算，并將這四種湍流模型的計(jì)算結(jié)果同相應(yīng)的實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，以此判定這四種湍流模型的精確度。結(jié)果表明，WALE大渦模型計(jì)算所得各特征參數(shù)與實(shí)驗(yàn)結(jié)果相比較，誤差比其余三種湍流模型計(jì)算結(jié)果最小[15，16]。因此下文將選用LES計(jì)算方法對(duì)爆燃過程的湍流信息和流場(chǎng)精細(xì)化特征進(jìn)行仿真研究。

預(yù)混氣體燃燒中大渦模擬的基本原理是模擬復(fù)雜的燃燒化學(xué)反應(yīng)過程和復(fù)雜爆炸壓力波系統(tǒng)的形成。大渦模型主要由連續(xù)性方程、動(dòng)量方程、能量方程、反應(yīng)進(jìn)程變量方程和理想氣體狀態(tài)方程組成：

2 ?幾何模型構(gòu)建和網(wǎng)格劃分

為全面反映不同分布形式的多分支結(jié)構(gòu)坑道油氣爆炸火焰-流場(chǎng)-超壓的耦合機(jī)理以及火焰特殊形態(tài)結(jié)構(gòu)特征形成機(jī)理，在本文的數(shù)值模擬部分，主要建立了3個(gè)幾何模型工況形式：分別為一字分布型（the lined layout）、相對(duì)分布型（the relative layout）、交錯(cuò)分布型（the staggered layout），如圖1-a、1-b、1-c所示。

在執(zhí)行模擬計(jì)算之前對(duì)計(jì)算的幾何模型進(jìn)行網(wǎng)格劃分。為保證求解精度，本文采用邊長(zhǎng)為4 mm 的六面體結(jié)構(gòu)網(wǎng)格對(duì)坑道模型進(jìn)行區(qū)域劃分，總網(wǎng)格數(shù)量為 125 萬個(gè)。由于爆炸反應(yīng)持續(xù)時(shí)間極短，坑道的壁面與火焰鋒面之間的熱交換極其有限，因而將坑道壁面邊界（Wall）都設(shè)置為絕熱無滑移。初始條件設(shè)為：計(jì)算區(qū)域壓力為大氣壓力，即超壓P0=0 Pa;初始反應(yīng)物選用體積濃度為1.5%的汽油蒸氣，層流火焰速度假設(shè)為定值0.34 m/s;點(diǎn)火區(qū)域的初始溫度T0設(shè)置為3 000 K，其他區(qū)域初始溫度T0=300 K;流場(chǎng)初速度以及反應(yīng)進(jìn)程變量均初始化為0。在左側(cè)坑道中心位置的封閉端，Patch設(shè)定半徑為2 mm的半球形（Sphere）點(diǎn)火區(qū)域（反應(yīng)進(jìn)度變量c=1），用以模擬點(diǎn)火功能[17]。

3 ?結(jié)果分析

3.1 ?火焰?zhèn)鞑ミ^程分析

圖2、圖3為三種分布形式的多分支坑道內(nèi)油氣爆炸火焰鋒面位置以及火焰?zhèn)鞑ニ俣鹊姆抡娼Y(jié)果的對(duì)比，如圖所示。從火焰鋒面位置曲線可見，三種分布形式的變化趨勢(shì)相對(duì)均勻和一致。隨著坑道內(nèi)反應(yīng)進(jìn)程的進(jìn)行，坑道內(nèi)火焰前鋒距離右側(cè)點(diǎn)火端的距離逐漸增大，增長(zhǎng)趨勢(shì)均呈現(xiàn)單調(diào)遞增規(guī)律，一字分布型和交錯(cuò)分布型曲線近似重合0～40 ms（即火焰從點(diǎn)火段傳播到第三個(gè)分支的時(shí)間區(qū)間）一字分布型斜率略大，60 ms（火焰從點(diǎn)火段傳播到第四個(gè)分支的時(shí)間區(qū)間）后兩者重合，兩者斜率皆大于相對(duì)分布型。

一字分布型和交錯(cuò)分布型曲線變化趨勢(shì)相似，火焰?zhèn)鞑ニ俣犬a(chǎn)生四個(gè)峰值，對(duì)應(yīng)四個(gè)分支結(jié)構(gòu)處的條件變化，其中一字分布型速度峰值高于交錯(cuò)分布型，且20～60 ms區(qū)間的鋒面位置早于交錯(cuò)分布型，根據(jù)仿真結(jié)果推測(cè)火焰徑向傳播至某一側(cè)管壁產(chǎn)生的渦旋會(huì)影響火焰的徑向傳播速度和輕微阻礙火焰的軸向傳播速度，因此使得同側(cè)分布分支比交錯(cuò)兩側(cè)分布分支的主坑道火焰先傳播到下一個(gè)分支位置，因此交錯(cuò)分布型主坑道火焰速度峰值大于一字分布型。相對(duì)分布型較前兩種分布形式主坑道火焰?zhèn)鞑ニ俣让黠@更小，說明總?cè)莘e不變時(shí)，說明總?cè)莘e不變時(shí)，一字分布型和交錯(cuò)分布型較相對(duì)分布型火焰爆炸傳播過程顯著加強(qiáng)。

3.2 ?爆炸超壓與火焰?zhèn)鞑ヱ詈详P(guān)系

結(jié)合圖4、5進(jìn)行討論分析，根據(jù)壓力上升速率的變化規(guī)律曲線來分析，火焰?zhèn)鞑ミ^程可劃分為從點(diǎn)火到發(fā)展的三個(gè)階段。

第一階段（0到20 ms）為點(diǎn)火初期，火焰鋒面還未接觸多分支結(jié)構(gòu)坑道，此時(shí)反應(yīng)初始，僅在點(diǎn)火頭周圍反應(yīng)，升壓速率較低，超壓變化幅度非常小。三種分布形式的變化規(guī)律呈現(xiàn)單調(diào)遞增規(guī)律，一字分布型和交錯(cuò)分布型升壓速率均大于相對(duì)分布型。第二階段（20～126 ms），一字分布型和交錯(cuò)分布型兩者曲線基本一致，將其視為等同一并分析。20 ms時(shí)，此時(shí)一字分布型和交錯(cuò)分布型的火焰鋒面前端恰好傳播至第一個(gè)分支，火焰形狀開始變得不規(guī)則。并且在20～60 ms 之間，隨著火焰?zhèn)鞑ソ?jīng)過分支數(shù)量的增多，超壓上升速率也逐漸增大。并且當(dāng)火焰?zhèn)鞑ソ?jīng)過四個(gè)分支之后，兩者的升壓速率于70 ms處達(dá)到最大峰值。相對(duì)分布型軸向火焰?zhèn)鞑ニ俣嚷谇皟烧?，故升壓速率加速過程慢于前兩者，于126 ms達(dá)到峰值。第三階段（126 ms之后）由于從該時(shí)刻后火焰只在主坑道內(nèi)傳播且多分支結(jié)構(gòu)坑道未燃?xì)怏w減少引起的火焰局部熄滅，使得火焰面積大幅減少，導(dǎo)致火焰?zhèn)鞑ニ俣妊杆贉p小，從而使得升壓速率驟然下降，相對(duì)分布型升壓速率小于于前兩者。