羅振敏，李 逵，毛文龍，張 江，程方明，王 濤，傅 文

0 引言

以甲烷為主要成分的天然氣作為一種清潔能源日益廣泛地應用在工業(yè)和生活中，一直以來，由含甲烷的可燃氣體引發(fā)的爆炸事故時有發(fā)生，造成嚴重的財產(chǎn)損失和人員傷亡。乙烯作為可燃氣體，具有爆炸極限小、壓力峰值高和危險度大等特性，若與甲烷混合會增加混合氣體的活性，使得甲烷爆炸極限范圍明顯增大，加劇甲烷爆炸，而造成更大的危害［1］。因此，開展C2H4/CH4混合氣體爆炸的研究對安全防爆具有一定指導意義。

目前，國內(nèi)外學者主要對初始條件、點火能量、障礙物存在、惰性及部分可燃性氣體等影響甲烷爆炸的因素進行研究［2－7］，在爆炸極限方面，陳曉坤等采用近球型氣體爆炸反應罐測試了在空氣中C2H6對甲烷爆炸極限參數(shù)的影響［8］;周利華對礦井火區(qū)可燃性混合氣體爆炸三角形判別法及其爆炸危險性進行了研究［9］;Van den Schoor和Verplaetsen研究了在溫度和壓力升高時，空氣中C2H6和 C3H8爆炸極限的變化規(guī)律［10－11］;Bolk 等研究了流動狀態(tài)下 C2H6—Air—N2混合氣體的爆炸極限［12］。而涉及到乙烯對甲烷爆炸特性的研究較少，僅部分學者開展了考察，研究表明，少量的乙烯氣體可使空氣中甲烷爆炸的界限加寬，危險度增大［1］;C2H4與CO之間存在著協(xié)同氧化誘導效應以及C2H4/CH4混合氣體的爆炸極限臨界氧濃度降低，使甲烷爆炸危險性增大［13－14］，這些研究也只是在宏觀上進行了定性分析，并未從微觀角度深入研究。隨著鏈式反應理論的提出和化學動力學的發(fā)展，發(fā)現(xiàn)可燃氣體爆炸并不是單一化學反應，而是一種復雜的熱爆炸和鏈式反應，自由基的生成和消失直接影響反應的進行［15］。因此，如果從微觀分子角度來定量地分析C2H4/CH4混合爆炸反應的機理，那么將會深入了解C2H4影響甲烷爆炸的規(guī)律。為此，文中基于氣體爆炸極限測試實驗，研究乙烯對甲烷爆炸極限的影響，同時應用Gaussian 09程序，采用密度泛函理論的B3LYP/6－31G方法優(yōu)化分子構型，并對相關基元反應作定量分析，從熱力學和動力學角度，研究C2H4對甲烷鏈式爆炸的微觀作用機理，得出C2H4/CH4混合體系步簡化反應機理，從微觀角度分析乙烯對甲烷爆炸的影響規(guī)律。

1 實驗與計算方法

1.1 實驗系統(tǒng)

實驗系統(tǒng)采用標準的可燃氣體爆炸極限測試裝置，如圖1所示，主要由爆炸反應管、點火器、循環(huán)泵、真空壓力計、泄壓口等組成。爆炸反應管長為1 400±50 mm，管內(nèi)徑60±5 mm，底部裝置通徑為25 mm的泄壓閥;配氣、循環(huán)攪拌、點火等步驟均在計算機上自動控制完成，配氣精確度為0.1%，采用電壓為300 VA電壓互感器為點火電源，電極間距設置3～4mm，放電時間0．1～1 s可調(diào);實驗在常壓下進行，初始溫度18～26℃，環(huán)境濕度56～64%RH，分別測試加入不同體積分數(shù)(0%，0.4%，0．8%，1．2%，1．6%，2．0%)的 C2H4氣體后CH4在空氣中的爆炸極限，循環(huán)攪拌時間為5 min，點火放電時間為500 ms．

1.2 計算方法

運用Gaussian 09量子化學程序包，采用密度泛函理論的 B3LYP/6－31G方法［16］，優(yōu)化反應物(Re)、中間體(IM)、過渡態(tài)(TS)和產(chǎn)物(P)的幾何構型，計算出相關基元反應的熱力學數(shù)據(jù)如活化能、反應焓變、吉布斯自由能變，并用反應內(nèi)稟坐標法驗證反應路徑［17］，來考察C2H4對甲烷鏈式爆炸的微觀機理。

圖1 可燃氣體爆炸極限測試裝置Fig．1 Test device for explosion limits of combustible gases

2 實驗結果及計算分析

2.1 C 2 H 4對空氣中CH 4爆炸極限的影響

實驗測得 CH4/air中加入不同體積分數(shù)的C2H4氣體時CH4的爆炸極限，且計算甲烷的爆炸危險度及爆炸極限范圍如圖2和圖3所示，來研究C2H4對甲烷爆的影響規(guī)律。其中，可燃氣體的爆炸危險度F是爆炸上下限的平均值和下限值的差再除以平均值［18］，其表達式 F=(槡UL －L)/，U為爆炸上限，%;L為爆炸下限，%，F(xiàn)值范圍為0～1，F(xiàn)值越接近于1則表明CH4爆炸危險性越大。

由圖2可知，隨著C2H4濃度的增加，爆炸極限范圍增大，爆炸下限下降比較明顯。圖3表明C2H4濃度增加，空氣中CH4爆炸危險度F值逐漸增大，2%的 C2H4使得 CH4爆炸危險度增大約59.1%．表明，C2H4在一定程度上增大了CH4在空氣中的爆炸危險性。

2.2 C2 H 4/CH 4鏈式爆炸反應過程計算及分析

對于甲烷燃爆反應機理已有了很多相關報道，若某物質的存在能生成主要中間產(chǎn)物或優(yōu)先于某些鏈反應發(fā)生，從而生成更多中間產(chǎn)物，那么就提高了反應體系活化中心濃度，促進了甲烷爆炸［19－23］?；诖?，現(xiàn)借鑒敏感性分析甲烷燃爆反應的簡化機理［24］，結合C2H4/CH4混合體系中可能發(fā)生的反應，選取其中11個基元反應，從反應發(fā)生的難易程度和概率出發(fā)，進行反應機理的簡化，來分析C2H4對甲烷爆炸影響的微觀作用機理。

圖2 C2 H4對CH4爆炸極限的影響Fig．2 Influence of ethylene on explosion limits of methane

圖3 C2H4對CH4爆炸危險度的影響Fig．3 Influence of ethylene on explosive hazard degree of methane

2．2．1 反應歷程分析

圖4為各反應中間體、過渡態(tài)及生成物構型，紅色虛線表示原子間空間距離，并未實際成鍵。由圖4可以看出，IM1的碳碳雙鍵直接裂解到生成物P1，并不存在過渡態(tài)TS1，IM2中的一個氧原子與·OH達到成鍵距離后形成·HO2;IM3中一個氫原子發(fā)生斷裂，與O2結合生成P3;IM4中碳碳雙鍵斷裂后與O2的2個氧原子分別結合形成P4，IM5碳氧雙鍵發(fā)生斷裂后一個氧原子與羥基結合形成P5，IM6中碳氧鍵發(fā)生斷裂后氧原子被氫原子取代，形成P6，IM7中的一個碳氫鍵發(fā)生斷裂后，2個游離的氫原子結合，形成P7;IM8中游離的氧原子攻擊碳氫鍵，并與斷裂后的氫原子結合形成P8，IM9中同樣是游離的氧原子攻擊碳氫鍵，并與斷裂后的氫原子結合形成P9，IM10中碳氫鍵發(fā)生斷裂后的氫原子與游離的羥基結合形成P10，IM11中游離羥基攻擊碳氫鍵，碳氫單鍵發(fā)生斷裂后氫原子與羥基結合生成P11．通過圖4所示的計算得知，各中間體和生成物均存在。

圖4 B3LYP/6－31G水平下各反應中間體、過渡態(tài)、及生成物的構型(鍵長單位/10－10m，鍵角單位/(°))Fig．4 Geometries of intermediate，transition states and products at the B3LYP/6 －31G level(Bond lengths/10 －10 m，bond angles/(°))

2．2．2 熱力學與動力學分析

11個基元反應的反應物、過渡態(tài)和產(chǎn)物以及正逆反應活化能、反應焓變和吉布斯自由能見表1～4．其中，Ea為正反應活化能;Eb為逆反應活化能;ΔHθ為反應焓變;ΔGθ為反應吉布斯自由能變。

在初始的CH4/C2H4混合體系中，首先可能發(fā)生的反應只有反應1，2，3，從反應熱力學角度出發(fā)，由表1可知，3個反應的ΔGθ均大于0，因此常溫常壓下均不能自發(fā)進行。對比反應(1)C2H4→·C2H3+·H和反應(3)C2H4+O2→·C2H3+·HO2的熱力學數(shù)據(jù)可看出，2個反應的ΔHθ都小于0屬于放熱反應，反應(3)的正反應活化能Ea為141．331 kJ/mol，逆反應活化能 Eb為 328．608 kJ/mol，此反應Ea＜Eb，反應有利于正向進行。反應(1)ΔE 為 28．959 kJ/mol，不存在過渡態(tài)，雖然反應(3)有利于反應正向進行，但對比化學反應平衡，對于從反應物開始反應時的次序幾乎無影響，而反應(3)的反應活化能小于反應(1)的反應活化能，同樣條件下，在此為了方便對比分析，簡化描述為率先發(fā)生是反應(1)。反應(2)CH4+O2→·CH3+HO2的焓變ΔHθ大于0為吸熱反應，正反應活化能 Ea為 289．645 kJ/mol，逆反應活化能 Eb為36．494 kJ/mol，此反應 Ea＞ Eb，反應有利于逆向進行。對比反應(2)和反應(3)的熱力學數(shù)據(jù)可看出，反應(3)有利于反應正向進行，也就是說更有利于C2H4與O2結合;反應(3)為放熱反應，同時反應(3)的正反應活化能Ea小于反應(2)，反應(3)發(fā)生所需跨越的能壘較小，因此，同樣條件下，反應(3)比反應(2)優(yōu)先發(fā)生。故此處可以說反應(1)首先引發(fā)了C2H4/CH4混合體系爆炸鏈鎖反應，為后續(xù)爆炸鏈反應提供熱量以及關鍵自由基·H 和·C2H3．

表1 反應(1)～(4)的熱力學數(shù)據(jù)Tab． Thermodynamic data of reaction(1)～(4)

由于反應1比反應3的優(yōu)先發(fā)生，可以使體系中的C2H4大部分反應生成·C2H3，因此在后續(xù)的鏈式反應中，體系中與O2的反應主要體現(xiàn)為反應2，4．由表1可知，鏈反應(4)·C2H3+O2→·CH2O+HCO正反應活化能Ea為22．133 kJ/mol遠小于反應 2的 289．645 kJ/mol，逆反應活化能 Eb為199．853 kJ/mol，焓變 ΔHθ為 －177．720 kJ/mol，為放熱反應，ΔGθ＜0，此反應在常溫常壓下能自發(fā)進行，同時，此反應 Ea＜Eb，反應有利于正向進行。比較反應(2)的熱力學數(shù)據(jù)可以發(fā)現(xiàn)，反應(4)比反應(2)更容易發(fā)生，也就是說在反應(1)首先引發(fā)初始混合體系鏈鎖反應后，·C2H3可以優(yōu)先或代替CH4去與O2發(fā)生反應，同時為體系提供熱量和鏈分支反應(5)的反應物CH2O以及爆炸過程關鍵自由基·HO2的生成，促進了CH4爆炸的鏈引發(fā)和鏈分支反應。因此，C2H4的加入使得CH4在空氣中更易發(fā)生爆炸反應，爆炸危險性將增大，這也解釋了C2H4的加入使得空氣中CH4爆炸下限出現(xiàn)明顯下降和爆炸危險度F值的增大。

由表2可知，反應(5)CH2O+O2→·HO2+·HCO正反應活化能Ea為193．184 kJ/mol遠小于反應(2)的 289．645 kJ/mol，逆反應活化能 Eb為50.515 kJ/mol，ΔHθ為 141．646 kJ/mol，是吸熱反應，ΔGθ＞0，此反應在常溫常壓下不能自發(fā)進行，同時，此反應Ea＞Eb，反應有利于逆向進行。同理，反應(6)·H+CH2O→·O+CH3為放熱反應，在常溫常壓下能自發(fā)進行，而反應(7)CH4+·H→·CH3+H2是吸熱反應且不能自發(fā)進行，且反應(6)的Ea小于反應(7)的Ea，因此反應(6)可以優(yōu)先發(fā)生，可以說在反應4優(yōu)先反應生成更多的CH2O后，CH2O仍可優(yōu)先與O2發(fā)生反應，增加了CH4爆炸過程關鍵自由基·O，這也就是說C2H4的加入間接增加了CH4燃爆鏈分支反應的反應物，進而促進了整個鏈式反應的進行，宏觀上可以表現(xiàn)為促進了CH4的爆炸。

表2 反應(5)～(7)的熱力學數(shù)據(jù)Tab．2 Thermodynamic data of reaction(5) ～(7)

表3 反應(8)～(9)的熱力學數(shù)據(jù)Tab．3 Thermodynamic data of reaction(8) ～(9)

表4 反應(10)～(11)的熱力學數(shù)據(jù)Tab．4 Thermodynamic data of reaction(10) ～(11)

由表3可知，反應(8)CH4+·→O·CH3+·OH為吸熱反應，在常溫常壓下不能自發(fā)進行，反應(9)C2H4+·→O·C2H3+·O屬于吸熱反應且能自發(fā)進行，有利于正向反應，同時反應(9)的Ea遠小于反應(8)的 Ea，因此，率先發(fā)生的是反應(9)。故此處可以說C2H4優(yōu)先CH4與·O反應，增加了CH4爆炸過程關鍵自由基·C2H3和·OH，促進了甲烷爆炸。

由表4可知，反應(10)CH4+·H→·CH3+H2O和反應(11)C2H4+·OH→·C2H3+H2O均為放熱反應且在常溫常壓下都能自發(fā)進行，但反應(10)的正反應活化能小于反應(11)的正反應活化能，因此，率先發(fā)生是反應(10)。雖然CH4優(yōu)先C2H4與·OH反應，但C2H4的存在增加了關鍵自由基·C2H3，促進了爆炸鏈分支反應物CH2O的生成。

為進一步驗證所研究反應歷程的準確性，利用內(nèi)稟反應坐標法計算反應通道。受限于篇幅，僅以反應(1)C2H4→·C2H3+·H和反應(2)CH4+O2→·CH3+·HO2為例展示反應路徑追蹤圖，如圖5，圖6所示。從各反應IRC反應路徑分析，皆經(jīng)歷了從初始低能量點構型到生成低能量穩(wěn)定構型產(chǎn)物的過程，直接驗證了文中11個反應的正確性，反應的能量變化情況與實際反應相符。

圖5 C2H4 →·C2H3+·H IRC反應路徑追蹤Fig．5 C2 H4 →·C2 H3+·H IRC reaction path tracking

圖6 CH4+O2 →·CH3+·HO2 IRC反應路徑追蹤Fig．6 CH4+O2 →·CH3+·HO2 IRCreaction path tracking

2．2．3 C2H4/CH4混合體系鏈式爆炸反應簡化機理

在C2H4/CH4混合體系爆炸反應過程中，一開始整個反應系統(tǒng)中充滿著C2H4，CH4和空氣。在點火的一瞬間，混合系統(tǒng)中首先可能發(fā)生如下反應

反應(1)的正反應活化能最小，反應發(fā)生所需跨越的能壘較小，不存在過渡態(tài)，所需吸收的熱量極小，所以在同樣的條件下，反應(1)可以優(yōu)先發(fā)生，或者說反應(1)發(fā)生的概率遠大于反應(2)和(3)，在此為了方便對比分析，簡化描述為在C2H4/CH4混合體系中率先發(fā)生是反應(1)。當反應(1)、(3)、(2)發(fā)生后，混合體系中除了 C2H4，CH4和 O2外，還有反應(1)產(chǎn)生的新的自由基·H，·C2H3，CH3，HO2．接著可能發(fā)生如下反應

同理上述分析，對比各基元反應熱力學和反應動力學參數(shù)，從反應發(fā)生的難易程度和概率出發(fā)，分析得出C2H4/CH4混合體系鏈式爆炸引發(fā)階段反應簡化機理如下

甲烷燃爆的基元反應是一個復雜的動態(tài)過程，并不能簡單認為每個基元反應在體系中存在狹義的順序。上述基元反應順序并不能說明這些反應是一個接一個的發(fā)生，而是就體系中存在的自由基可能發(fā)生的反應和反應的難易程度及概率進行了歸納。

3 結論

1)隨著C2H4體積分數(shù)的不斷增大，空氣中CH4爆炸上限和下限均下降，且爆炸下限下降比較明顯，而對上限的影響幅度不大，CH4爆炸極限范圍逐漸變寬，爆炸危險度F值增大;

2)C2H4→·C2H3+·H首先引發(fā)了C2H4/CH4混合體系爆炸鏈鎖反應并釋放能量，為后續(xù)鏈式反應提供所需的自由基，降低了混合體系爆炸發(fā)生的困難程度，解釋了C2H4加入后CH4爆炸下限下降的原因;

3)C2H4→·C2H3+·H產(chǎn)物中·C2H3優(yōu)先CH4與O2發(fā)生反應，增加了CH4爆炸鏈分支反應所需的中間產(chǎn)物CH2O，促進了甲烷爆炸，同時C2H4的存在增加了CH4爆炸過程中關鍵自由基·O，·C2H3，·OH，·H，·OH 以及·HCO，提高了反應中中間關鍵產(chǎn)物的活化中心濃度，在一定程度上促進了甲烷的爆炸;

4)通過分析C2H4對甲烷/空氣爆炸的微觀作用機理，得出了C2H4/CH4混合體系鏈式爆炸的16步簡化反應機理。

［1］ 賈寶山，胡如霞，皮子坤，等．采空區(qū)遺煤自燃產(chǎn)生的C2H4促進瓦斯爆炸特性［J］．遼寧工程技術大學學報:自然科學版，2015，34(6):677 －682．JIA Bao-shan，HU Ru-xia，PI Zi-kun，et al．Promoting gas explosion characteristics of C2H4produced by spontaneous combustion of residual coal in gob［J］．Journal of Liaoning Technical University:Natural Science Edition，2015，34(6):677 －682．

［2］ 左 哲，姚志強，高進東，等．受限空間內(nèi)天然氣爆炸反應過程模擬分析［J］．天然氣工業(yè)，2015，35(6):131－137．ZUO Zhe，YAO Zhi-qiang，GAO Jin-dong，et al．Simulation analysis of natural gas explosion within confined space［J］．Natural Gas Industry，2015，35(6):131－137．

［3］ 劉 丹，司榮軍，李潤之．環(huán)境濕度對瓦斯爆炸特性的影響［J］．高壓物理學報，2015，29(4):307 －312．LIU Dan，SI Rong-Jun，LI Run-zhi．Ambient humidity influence on explosion characteristics of methane-air mixture［J］．Chinese Journal of High Pressure Physics，2015，29(4):307 －312．

［4］ 李潤之，司榮軍．點火能量對瓦斯爆炸壓力影響的實驗研究［J］．礦業(yè)安全與環(huán)保，2010，37(2):14 －16．LI Run-zhi，SI Rong-Jun．Experimental study on impact of ignition energy on gas explosion pressure［J］．Mining Safety and Environmental Protection，2010，37(2):14－16．

［5］ 張應安，劉振翼，王 峰，等．含CO2天然氣燃燒爆炸特性實驗研究［J］．天然氣工業(yè)，2009，29(6):110－112．ZHANG Ying-an，LIU Zhen-yi，WANG Feng，et al．Experimental studies on fire and explosion characteristics of natural gas with carbon dioxide［J］．Natural Gas Industry，2009，29(6):110 －112．

［6］ 任韶然，李海奎，李磊兵，等．惰性及特種可燃氣體對甲烷爆炸特性的影響實驗及分析［J］．天然氣工業(yè)，2013，33(10):110 －115．REN Shao-ran，LI Hai-kui，LI Lei-bing，et al．An experimental study of effects of inert and special flammable gases on methane’s explosion characteristics［J］．Natural Gas Industry，2013，33(10):110 －115．

［7］ 程方明，鄧 軍．一氧化碳影響二氧化碳惰化甲烷爆炸的實驗研究［J］．西安科技大學學報，2016，36(3):315－319．CHENG Fang-ming，DENG Jun．Experimental study on the influence of carbon monoxide on carbon dioxide-inerted methane explosion［J］．Journal of Xi’an University of Science and Technology，2016，36(3):315 －319．

［8］ 陳曉坤，張 娥．可燃性氣體對CH4爆炸特性影響的實驗研究［J］．煤礦安全，2009，40(10):7 －10．CHEN Xiao-kun，ZHANG E．Experimental study on the methane mixed combustive gas explosive characteristic［J］．Safety in Coal Mines，2009，40(10):7 － 10．

［9］ 周利華．礦井火區(qū)可燃性混合氣體爆炸三角形判斷法及其爆炸危險性分析［J］．中國安全科學學報，2001，11(2):47 －51．ZHOU Li-hua．Explosion triangle of mixed inflammable gases and explosion risk analysis in the mine fire district［J］．China Safety Science Journal，2001，11(2):47－51．

［10］ Van Den Schoor F，Verplastsen F．The upper explosion limit of lower alkanes and alkenes in air at elevated pressures and temperatures［J］．Journal of Hazardous Materials，2006，128(1):1 －9．

［11］ Norman F，Van Den Schoor F，Verplastsen F．Auto-ignition and upper explosion limit of rich propane-air mixtures at elevated pressures［J］．Journal of Hazardous Materials，2006，137(2)，666 －671．

［12］ Bolk J W，Siccama N B，Westerterp K R．Flammability limits in flowing ethane-air-nitrogen mixtures:An experimental study［J］．Chemical Engineering Science，1996，51(10):2231－2239．

［13］鄧 軍，吳曉春，程 超．CH4，CO，C2H4多元可燃氣體爆炸實驗研究［J］．煤炭現(xiàn)代化，2007(5):63－65．DENG Jun，WU Xiao-chun，CHENG Chao．Experimental study on explosion of multicomponent combustibility gases containing CH4，CO and C2H4［J］．Coal Mine Modernization，2007(5):63 －65．

［14］鄧 軍，程 超，吳曉春．煤礦可燃性氣體爆炸氧濃度的實驗研究［J］．煤礦安全，2007，38(6):5－7．DENG Jun，CHENG Chao，WU Xiao-chun．Experimental study on explosive oxygen concentration of coal mine combustive gases［J］．Safety in Coal Mines，2007，38(6):5－7．

［15］高 娜，張延松，胡毅亭．受限空間瓦斯爆炸鏈式反應動力學分析［J］．中國安全科學學報，2014，24(1):60－65．GAO Na，ZHANG Yan-song，HU Yi-ting．Dynamics analysis of gas explosion chain reaction in restricted space［J］．China Safety Science Journal，2014，24(1):60－65．

［16］ Karunakaran Remya，Cherumuttathu H Suresh．Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small non-covalent dimers?A benchmark study using gaussian 09［J］．Journal of Computational Chemistry，2013，34(15):1341－1353．

［17］ Bakken V，Helgaker T．The efficient optimization of molecular geometries using redundant internal coordinates［J］．Journal of Chemical Physics，2002，117(20):9160－9174．

［18］近藤重雄，韓 美．新的爆炸危險度F值［J］．低溫與特氣，1995(2):57－60．JINTENG Zhong-xiong，HAN Mei．A new explosive hazard degree F［J］．Low Temperature and Specialty Gases，1995(2):57 －60．

［19］羅振敏，鄧 軍，郭曉波．基于Gaussian的瓦斯爆炸微觀反應機理［J］．遼寧工程技術大學學報，2008，27(3):325－328．LUO Zhen-min，DENG Jun，GUO Xiao-bo．Microcosmic mechanism of gas explosion based on Gaussion［J］．Journal of Liaoning Technical University，2008，27(3):325－328．

［20］羅振敏．鏈式反應理論在礦井瓦斯抑爆中的應用［J］．煤礦安全，2009(2):67 －69．LUO Zhen-min．Application of chain reaction theory in mine gas explosion suppression［J］．Safety in Coal Mines，2009(2):67 － 69．

［21］羅振敏，張 江，王 濤．CO/CH4鏈式爆炸反應機理研究［J］．兵工學報，2017，38(S1):49 －59．LUO Zhen-min，ZHANG Jiang，WANG Tao．The study of CO/CH4explosion chain reaction mechanism［J］．Acta Armamentarii，2017，38(S1):49 －59．

［22］羅振敏，張 江，任軍瑩，等．NH4H2PO4粉體熱分解產(chǎn)物在瓦斯爆炸中的作用［J］．煤炭學報，2017，42(6):1489－1495．LUO Zhen-min，ZHANG Jiang，REN Jun-ying，et al．Role of thermal decomposition products of NH4H2PO4powder in gas explosion［J］．Journal of China Coal Society，2017，42(6):1489 －1495．

［23］羅振敏，康 凱，任軍瑩．NH3對甲烷鏈式爆炸的微觀作用機理［J］．煤炭學報，2016，41(4):876 －883．LUO Zhen-min，KANG Kai，REN Jun-ying．Microscopic mechanism of NH3on chain of methane explosion［J］．Journal of China Coal Society，2016，41(4):876 －883．

［24］侯金麗，金 平，蔡國飆．基于敏感性分析的氧/甲烷燃燒反應簡化機理［J］．航空動力學報，2012，27(7):1549－1554．HOU Jin-li，JIN Ping，CAI Guo-biao．Reduced mechanism for oxygen/methane combustion based on sensitivity analysis［J］．Journal of Aerospace Power，2012，27(7):1549－1554．