汽油蒸氣與異辛烷定容爆炸壓力的當量比等效

張培理　宋哲　張德洋　劉沖　王建

摘????? 要：在理論計算異辛烷定容絕熱爆炸壓力的基礎上，通過標準20 L球形爆炸裝置內的汽油蒸氣爆炸實驗，獲得了不同汽油蒸氣體積濃度下的爆炸壓力。通過分析異辛烷定容爆炸壓力理論計算值與汽油蒸氣定容爆炸壓力實驗值之間的誤差，提出了汽油蒸氣與異辛烷定容爆炸壓力的當量比等效關系式。該關系式不僅可以方便地用異辛烷來預測汽油蒸氣爆炸壓力，還可以修正汽油蒸氣爆炸壓力數值計算模型，使模型計算值與實驗值之間的誤差顯著降低。

關? 鍵? 詞：汽油蒸氣;爆炸壓力;當量比;等效關系式

中圖分類號：TQ038.4 ??????文獻標識碼： A?????? 文章編號： 1671-0460（2020）06-1013-04

Equivalent Ratio Equivalence for the Explosion Pressure of

Gasoline-Air and Isooctane-Air Mixture

ZHANG Pei-li¹， SONG Zhe²，ZHANG De-yang³， LIU Chong¹， WANG Jian¹

（1. Dept. of Fuel， Army Logistics University， Chongqing 401331， China;

2. Chongqing Fire Safety Research & Service Co.， Ltd.， Chongqing 401120， China;

3. 92274 Troops， Qingdao Shandong 266400， China）

Abstract： On the basis of calculating isooctane adiabatic explosion pressures under different equivalent ratios， the explosion pressures of the gasoline-air mixture under different gasoline vapor concentrations were obtained by carrying explosion experiments in the standard 20 L spherical explosive device. By analyzing the error between the theoretical value of isooctane explosion pressure and the experiment pressure of the gasoline-air mixture， the equivalence relation formula of the equivalent ratio for the gasoline vapor and isooctane explosion pressure was put forward. This relation formula can be used not only to predict gasoline vapor explosion pressure， but also to amend the gasoline vapor explosion pressure numerical calculation model， so that the error between the calculated value and the experimental value can be significantly reduced.

Key words： Gasoline vapor; Explosion pressure; Equivalent ratio; Equivalence relation formula

汽油作為重要的動力燃料之一，廣泛應用于各類動力設備。同時，汽油又是典型的混合物，成分構成復雜，氣相色譜-質譜分析表明其至少有94種成分，質量濃度由高到低依次為烷烴、烯烴、苯系物、環烷烴等^[1]。由于汽油成份的復雜性，因此其揮發成的油蒸氣成分也十分復雜。

目前相關領域在對汽油蒸氣的燃燒放熱、汽油蒸氣氧化動力學模型等開展研究時，為了簡化模型和便于分析，一般以汽油蒸氣替代物來開展相關實驗和計算，最常見的汽油蒸氣替代物有異辛烷（iso-octane）和PRF（Primary Reference Fuels）兩種。PRF是異辛烷和正庚烷的二元混合物，主要用于研究內燃機點火過程及辛烷值的確定等^[2-3]。異辛烷常常被國內外學者用來研究汽油蒸氣氧化過程的放熱^[4-5]、化學動力學模型^[6-13]以及火焰速度等特性^{[12， 14-15]}，因此在計算汽油蒸氣爆炸超壓、火焰溫度時，常用異辛烷來代替油蒸汽成份，符合相關領域的處理習慣。然而這樣的簡化處理，常常會造成理論計算值與實際數據之間存在一定誤差。為了消除這類誤差，必須深入研究異辛烷與油蒸汽之間的等效關系。

本文以異辛烷和汽油蒸氣為研究對象，一方面將在理論分析的基礎上，計算不同當量比下異辛烷的定容爆炸壓力。另一方面，通過20 L球形標準容器內的汽油蒸氣爆炸實驗，獲得不同汽油蒸氣體積濃度下汽油蒸氣定容爆炸壓力數據。通過分析異辛烷定容爆炸壓力計算值與實際汽油蒸氣爆炸實驗壓力數據之間的關系，試圖找出汽油蒸氣與異辛烷定容爆炸壓力的當量比等效關系式。從而為準確計算汽油蒸氣爆炸壓力以及汽油蒸氣爆炸過程數值分析建模等提供重要參考。

1 ?異辛烷定容爆炸壓力的理論計算

根據熱力學第一定律，當燃料-空氣混合物在定容條件下進行絕熱燃燒時，有如下的熱力學方程：

式中：V —容器體積，m³;

P_init和P_f —分別表示初始和終了狀態的壓力，Pa;

H_reac和H_prod —分別表示反應物和生成物的絕對焓，kJ/kmol。

其數值可由式（2）計算。

式中： ?—組分i標準狀態下的生成焓;

—組分i的定容比熱。

應用理想氣體狀態方程PV=nRT，消去P_fV項，即

式中： T_init—初始溫度，取值298 K;

T_ad，v—定容絕熱燃燒溫度，K;

n_reac和n_prod—分別表示反應物和生成物的摩爾數總和，其值與實際當量比或燃料體積濃度有關，可根據實際化學反應方程式求得。

根據式（3）及初始狀態和終態時的氣體狀態方程就可以計算出P_f。

上述計算過程可借助計算機來完成。表1為異辛烷（C₈H₁₈）當量比F 取值0.3～2.0時的計算結果。圖1為根據計算結果得到的T_ad，V和P_f隨F 的變化曲線。

由圖1可以看出，計算得到的異辛烷定容絕熱燃燒溫度T_ad，p和終態壓力P_f均隨當量比F 的增加呈現先上升后下降的變化規律。當F=1.1時定容絕熱燃燒溫度T_ad，p達到峰值（2 666.66 K），當F=1.2時終態壓力P_f達到峰值（989 207 Pa）。終態壓力P_f峰值對應的當量比F既不是化學計量濃度下的F=1.0，也不是最大定容絕熱燃燒溫度對應的F=1.1。造成這一現象的原因主要與燃燒熱、產物熱容、產物總摩爾數等因素有關，具體原因參見前期研究文獻^[16]，在此不再贅述。

2? 汽油蒸氣定容爆炸壓力實驗測定

2.1? 裝置與方法

本文采用標準20 L球形爆炸裝置（如圖2所示）測定汽油蒸氣爆炸壓力。相對于其他爆炸裝置來說，20 L球形爆炸裝置具有體積小、操作靈活方便、實驗消耗低等優勢，已成為實驗室研究可燃氣體（液體蒸氣）、粉塵爆炸的主流裝置。

實驗時，將左右兩側接口改成帶閥門雙通，分別用于連接配氣系統和濃度測試系統，壓力傳感器安裝在抽真空口。另外，為了模擬絕熱條件，實驗前在球體外側包上保溫絕熱材料以減少熱量的損失，改進后的實物如圖3所示。

實驗方法步驟如下：

1）將配氣系統、濃度測試系統與爆炸裝置連接好;

2）用清潔空氣對爆炸腔體進行通風吹掃，安裝好壓力傳感器，蓋上密封蓋;

3）開啟汽油蒸氣霧化循環系統，充入汽油蒸氣并進行循環攪拌;

4）循環2 min左右后停泵靜置1 min，利用濃度測試系統分析汽油蒸氣及氧濃度，若汽油蒸氣濃度未達到設定值再進行配氣動作直至符合要求;

5）關閉爆炸裝置所有閥門，連接點火器，調好壓力采集系統，而后點火;

6）記錄實驗數據;

7）打開密封蓋，排出廢氣，用清潔空氣進行吹掃，進行下一組實驗。

2.2? 結果與分析

基于上文實驗裝置與實驗步驟，控制初始油蒸汽體積分數步長在0.2%之內，汽油蒸氣初始體積分數在0.9%～2.4%之間，共完成了10組實驗。其典型的壓力曲線如圖4所示。

由圖4可以看出，汽油蒸氣在被點燃后，壓力迅速上升至峰值P_max，之后壓力緩慢降低。這主要是由于點火后，火焰呈球形由球形容器中心向容器內壁面迅速發展，當火焰燃燒至內壁面時，汽油蒸氣與空氣混合物燃燒過程結束，此時容器內溫度最大，壓力達到峰值。之后，溫度逐漸降低，壓力隨之緩慢降低。因此，壓力曲線上的壓力峰值P_max即是理論分析中的壓力終值P_f。

表2為10組實驗得到的壓力峰值P_max以及相同初始體積濃度下異辛烷定容爆炸壓力終值P_f。圖5為根據表2得到的對比圖。

由圖5可以看出，①汽油蒸氣爆炸壓力實驗值和異辛烷爆炸壓力理論值均隨初始體積分數（當量比）呈先上升后下降的趨勢。該趨勢可用開口向下的二次曲線描述。②相同當量比下，理論計算值均高于實驗值。這是因為，盡管實驗過程中對20 L球形容器外壁面采用了絕熱處理，但與理論計算時的理想絕熱條件相比仍有較大的熱損失。另一方面，汽油蒸氣的實際燃燒化學反應過程與異辛烷燃燒反應過程在放熱量、放熱速率等方面必然存在差異，這導致了實驗壓力與理論計算值之間的差異。③當Φ≤1時，實驗值與理論值之間的差值較小且較為一致，約為0.126 MPa;當Φ>1時，實驗值與理論值之間的差值隨Φ的增大而增大。

為了便于對定容絕熱條件下汽油蒸氣與異辛烷爆炸壓力進行等效分析，將實驗值與理論值的差值按對應的當量比繪圖并擬合，如圖6所示。

結合圖5和圖6，針對汽油蒸氣爆炸壓力實驗值與異辛烷爆炸壓力理論計算值之間的差異，現提出以下等效關系式：

式（5）即為汽油蒸氣定容爆炸壓力的異辛烷當量比等效關系式。該等效關系式的意義主要體現在以下幾個方面。

1）只要確定了實驗時的汽油蒸氣體積濃度，計算出對應的當量比后，便可以方便地預測汽油蒸氣與空氣預混后的定容爆炸壓力。

2）汽油蒸氣爆炸領域數值分析計算建模時，常常以異辛烷來代替汽油蒸氣的成份，然而由模型計算得到的壓力值往往比實際實驗值高。式（5）可以修正汽油蒸氣爆炸壓力計算模型，使模型計算值與實際實驗值之間的誤差顯著降低。

3? 結論

1）碳氫燃料與空氣預混后，其當量比可表示為燃料體積濃度、碳氫燃料分子式中碳原子和氫原子個數的函數。

2）當Φ≤1時，汽油蒸氣定容爆炸壓力實驗值與理論值之間的差值較小且較為一致;當Φ>1時，其實驗值與理論值之間的差值隨Φ的增大而增大。

3）汽油蒸氣定容爆炸壓力的異辛烷當量比等效關系式可由式（5）描述。只要確定了實驗時的汽油蒸氣體積濃度，便可以方便地預測汽油蒸氣與空氣預混后的定容爆炸壓力。

致謝：本文得到后勤工程學院青年科研基金資助（YQ16-420802）。

參考文獻：

[1]鄒麗，李傳真.火焰原子吸收光譜法測定汽油中錳含量[J].當代化工，2016，45（5）：1050-1051.

[2]SILEGHEM L， ALEKSEEV V A， VANCOILLIE J， et al. Laminar burning velocity of gasoline and the gasoline surrogate components iso-octane， n-heptane and toluene[J].Fuel， 2013， 112（3）： 355-365.

[3]CURRAN H J， PITZ W J， WESTBROOK C K， et al. Oxidation of automotive primary reference fuels at elevated pressures[C]. Twenty-Seventh Symposium （International） on Combustion， The Combustion Institute， 1998： 27（1）： 379–387.

[4]HALSTEAD M P， KIRSCH L J， QUINN C P. The autoignition of hydrocarbon fuels at high temperatures and pressures-fitting of a mathematical model[J].Combustion and Flame， 1977， 30（77）： 45-60.

[5]BOHM H，LACAS F. On extinction limits and polycyclic aromatic hydrocarbon formation in strained counter flow diffusion flames from 1 to 6 bar[J].Proceedings of the Combustion Institute， 2000， 28（2）： 2627–2634.

[6]LIU Y， JIA M，XIE M，et al. Improvement on a skeletal chemical kinetic model of iso-octane for internal combustion engine by using a practical methodology[J].Fuel， 2013， 103（1）： 884-891.

[7]TANAKA S， AYALA F， JAMES C. K. A reduced chemical kinetic model for HCCI combustion of primary reference fuels in a rapid compression machine[J].Combustion and Flame， 2003， 133（4）： 467–481.

[8]TSURUSHIMA T. A new skeletal PRF kinetic model for HCCI combustion[J].Proceedings of the Combustion Institute， 2009， 32（2）： 2835-2841.

[9]LU T.， CHUNG K. L. Linear time reduction of large kinetic mechanisms with directed relation graph：n-Heptane and iso-octane[J].Combustion and Flame， 2006， 144（1-2） ：24–36.

[10]梁振龍. 汽油燃料替代混合物化學動力學模型研究[D]. 重慶： 重慶大學， 2015.

[11]賈明， 解茂昭. 適用于HCCI發動機的異辛烷氧化的化學動力學模型（I）—現有模型的比較分析[J]. 內燃機學報， 2006 3 （24）： 227-234.

[12]COX R A，COLE J A，Chemical aspects of the autoignition of hydrocarbon-air mixture[J].Combustion and flame， 1985， 60 （2）： 109-123.

[13]CALLAHAN C V， HELD T J， DRYER F L， et al. Experimental data and kinetic modeling of primary reference fuel mixture[J].Proceedings of the Combustion Institute， 1996， 26（1）：739-746.

[14]JERZEMBECK S， PETERSA N， PEPIOT-DESJARDINS P，et al. Laminar burning velocities at high pressure for primary reference fuels and gasoline： experimental and numerical investigation[J].Combustion and Flame， 2009， 156（2）： 292-301.

[15]KELLEY A P， LIU W， XIN Y X，et al. Laminar flame speeds non-premixed stagnation ignition and reduced mechanisms in the oxidation of iso-octane[J].Proceedings of the Combustion Institute， 2011， 33（1）： 501-508.

[16]張培理， 齊圣， 王世茂， 程順國. 93號汽油蒸氣絕熱燃燒溫度的計算分析[J]. 后勤工程學院學報， 2017， 33 （2）： 33-39.