Experimental estimation of the combustion regime in the oil-gas explosion process
-
摘要: 首先分析讨论了油气爆炸过程中火焰燃烧模式的估计方法,然后在激波管内进行了低、中、高3次不同初始油气浓度条件下的油气爆炸实验,通过实验数据分别计算出了低、中、高初始油气浓度条件下油气爆炸在初期、中期和后期的丹姆克尔数和湍流雷诺数,最后依靠丹姆克尔数-湍流雷诺数图对低、中、高初始油气浓度条件下油气爆炸初期、中期和后期的火焰燃烧模式进行了定量估计。结果表明:低、中、高初始油气浓度条件下激波管油气爆炸过程初期、中期和后期的火焰燃烧模式均为漩涡内小火焰模式。Abstract: In this artical, firstly, estimation method of flame combustion regime of the oil-gas explosion was discussed and three oil-gas explosion experiments under the conditions of low, middle and high initial gas vapor concentration were carried out, and then the Damköhler number and the turbulent Reynolds number for the early, interim and late stage of the oil-gas explosion at low, middle and high initial gas vapor concentration conditions were calculated according to the experimental data. Finally, through the Damköhler number vs. Reynolds number diagram, the combustion regimes for each stage of the oil-gas explosion at low, middle and high initial gas vapor concentration conditions were quantitative estimated. Results show that the combustions at early, interim and late stage of the gas-air explosion under the conditions of low, middle and high initial gas vapor concentration in the tube have the same regime of flameletes-in-eddies. The conclusions of this paper can provide some useful reference for the further study of combustion regime and the numerical analysis model selection of the gas-oil explosion.
-
油气(包括天然气、瓦斯和原油、成品油挥发物等)属于典型的易燃易爆气体,稍有不慎就可能引发火灾爆炸安全事故,造成重大人员伤亡和巨额经济损失。近年来,随着国民经济与社会生产的快速发展,受限空间油气火灾爆炸事故时有发生,如2007年我国西南某油库发生的油库爆炸事故;2013年11月,山东青岛原油输油管线发生爆炸事故;2014年7月,台湾高雄发生燃气爆炸事故,以上事故均带来严重的人员伤亡和财产损失。为了避免此类事故的发生,减少人员伤亡与财产损失,尚需加强对油气爆炸机理的研究,从而提出有效的油气爆炸防治措施。
受限空间油气爆炸火焰属于典型的湍流预混火焰,湍流火焰的燃烧模式与油气爆炸过程的化学反应过程、湍流流动以及传热传质过程密切相关。一般来说,不同的油气爆炸火焰燃烧模式会导致不同的油气爆炸特性,如爆炸超压、超压上升速率、火焰速度等,究其内在的原因之一就是湍流会使火焰前锋面发生褶皱和扭曲[1],进而对油气爆炸过程的传热传质[2-3],甚至化学反应历程[4-5]产生影响。因此,研究受限空间油气爆炸火焰的燃烧模式对于进一步认识油气爆炸过程的机理和规律、构建油气爆炸过程数学模型具有十分重要的意义。
1. 油气爆炸火焰燃烧模式的估计方法
对湍流火焰燃烧模式的研究一直是可燃气体燃烧和爆炸领域的热点和难点[6]。F.A.Williams[7]和J.Abraham等[8]根据湍流柯尔莫格洛夫微尺度lk、湍流积分尺度l0和层流火焰厚度δL等参数构成的判据,将湍流预混火焰分为3种模式:褶皱层流火焰模式、漩涡内小火焰模式和分布反应模式。具体方法是先计算湍流预混火焰的丹姆克尔数Da和湍流雷诺数Rel0,然后根据反映湍流预混火焰3种模式的湍流雷诺数-丹姆克尔数图(如图 1所示),来确定湍流火焰的燃烧模式。
图 1中,粗实线(lk/δL=1)上方是满足威廉斯-克里莫夫判据的褶皱层流火焰存在的条件,粗实线(l0/δL=1)下方是满足丹姆克尔判据的分布反应模式存在的条件,两条粗实线之间是满足漩涡内小火焰模式的分布条件。因此只要计算出油气爆炸过程的丹姆克尔数Da和湍流雷诺数Rel0,根据图 1就可以判断油气爆炸火焰燃烧的模式。
丹姆克尔数Da在燃烧中是一个很重要的量纲一参数,定义为:
Da=τflowτchem=l0/v∗rmsδL/SL=(l0δL)(SLv∗rms) (1) 式中:τflow为流动特征时间,τchem为化学特征时间,vrms*为湍流脉动速度,SL为层流火焰速度。l0/SL为火焰几何尺度比,SL/vrms*为火焰相对湍流强度。
湍流雷诺数Rel0的定义为:
Rel0=ρv∗rmsl0μ (2) 式中:ρ为燃烧火焰的密度,μ可燃气体的为动力粘度。
从上面的分析可以看出,要想判断油气爆炸火焰的模式,必须得到油气爆炸过程的湍流脉动速度vrms*,湍流积分尺度l0,层流火焰厚度δL和层流火焰速度SL等参数,显然在实验室要想直接获得这些参数是非常困难的。
本文拟通过激波管中的油气爆炸实验,测量油气爆炸过程中已燃气体和未燃气体的温度、平均流速、爆炸压力等数据来估计油气爆炸过得的丹姆克尔数Da和湍流雷诺数Rel0,进而估计油气爆炸过程的火焰模式。
2. 实验装置和实验方法
本文中的实验装置主要包括激波管、真空循环泵、油气雾化装置、数据采集系统、点火系统、高速摄影仪、电脑等。数据采集系统,主要包括压力采集系统、浓度采集仪、温度采集设备和流场速度设备。实验装置布置示意图如图 2所示。
激波管尺寸为200 mm×200 mm×6 100 mm,其中包括一段长300 mm的玻璃观察段。油气雾化装置的结构和工作原理见参考文献[9-10]。
压力采集系统主要由压力传感器、采集卡和电脑组成。沿激波管共布置了3个压力传感器,实验中以这3个压力传感器采集到数据的平均值作为激波管内的压力值。浓度采集系统主要由GXH-1050型红外分析仪和NHA-502型汽车尾气分析仪组成。实验中油气浓度由GXH-1050型红外分析仪测试,而其他气体组分,如O2、CO、CO2等,均由NHA-502型汽车尾气分析仪测量。
实验中未燃气体和已燃气体的温度由热电偶、采集卡和电脑组成的温度采集设备测量。根据实验需要,沿激波管共设置了3组热电偶,分别布置在激波管50、305和560 cm处。
油气爆炸过程中湍流的脉动速度的测量较为困难,本文中拟先用6162型高温热线风速仪来记录油气爆炸过程中的瞬时气体流速,并根据整个油气爆炸持续的时间得到平均流速,然后根据任一时刻的瞬时速度和平均流速得到湍流脉动速度。该热线风速仪探头设置在激波管中间位置,可以每隔5~10 ms测一次气体的瞬时流速,并储存在计算机内。
不同的初始条件会导致不同的实验结果。为了使实验结果具有代表性,选取低(1.20%)、中(1.65%)、高(2.10%)3种初始油气体积分数进行实验,经计算,其油气当量比ϕ(按异辛烷计算)分别为0.72、1.00和1.28。其他初始条件设置如下:初始温度为常温,压力为常压,点火能设为5 J。
本文先以ϕ=1.00(即油气初始体积分数为1.65%)时的实验结果为例,给出火焰燃烧模式估计的详细步骤与结果;然后再列出初始油气体积分数分别为1.20%和2.10%时的实验数据,并分别估计其燃烧模式;最后根据3次燃烧模式计算的结果,分析和讨论激波管内油气爆炸火焰的燃烧模式。
3. 实验结果与火焰燃烧模式计算
3.1 实验结果与分析
图 3为初始油气体积分数为1.65%时(ϕ=1),流场速度v和爆炸超压p随时间变化的曲线。由图 3可以看出整个油气爆炸持续了约288 ms,在此过程中激波管内部的气体速度在-60~100 m/s之间不断振荡波动,其平均流速v≈21.20 m/s。激波管内气体流速任一时刻的脉动量vrms*可以由瞬时流速v(t)和平均流速v计算:
v∗rms=|v(t)−ˉv| (3) 图 4为初始油气体积分数为1.65%时(ϕ=1),3组热电偶采集到的流场温度T随时间的变化曲线。从图 4可以看出,未燃气体的温度Tu约为300 K,1#、2#和3#热电偶采集到的已燃气体的最大温度Tmax分别约为1 127.39、1 203.62和1 324.05 K。已燃气体的最大温度正是火焰经过热电偶时的燃烧气体的温度,火焰离开热电偶后,热电偶处的气体温度由于剧烈的热损失开始下降,之后温度趋于平稳,此时的温度即为燃烧后的气体温度为Tb。可以认为热电偶采集到的最大温度Tmax即为油气爆炸已燃气体的瞬时温度。由于1#、2#和3#热电偶分别设置在激波管的前端、中部和尾部,因此其各自采集到的最大温度代表了油气爆炸初期、充分发展和后期的燃烧气体的温度。
本文中拟分别对油气爆炸初期(50 ms)、充分发展(150 ms)和后期(250 ms)的火焰燃烧模式进行估计。从图 3和图 4可以得到ϕ=1,t=50、150和250 ms时激波管内气体流速、压力和已燃气体温度等数据,如表 1所示。
表 1 各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1)Table 1. Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1)t/ms v/(m·s-1) vrms* p/Pa Tmax/K 50 55.82 34.62 98 530 1 127.39 150 72.76 51.56 340 800 1 203.62 250 29.14 7.94 486 190 1 324.05 3.2 油气爆炸火焰燃烧模式计算
湍流积分尺度代表了湍流涡的平均尺度,根据杨辉等[11]和C.Sak等[11]的研究,管道或坑道内湍流的平均积分尺度约为管道或坑道横截面积的1/20~1/5。为了便于估算,设积分尺度为激波管横截面尺寸的1/10,即l0=0.02 m。根据式(1),先计算50 ms时的丹姆克尔数。
其中流动特征时间为τflow=l0/vrms*=0.000 6 s。
化学特征时间τchem=δL/SL,根据文献[13]:
SL=SL, ref(TuTu, ref)γ(ppref)β (5a) SL, ref=BM+B2(ϕ−ϕM)2 (5b) 式中:Tu, ref=298 K;pref为标准大气压;BM,B2和ϕM是由燃料类型确定的常数,由表 2给出;温度指数γ和压力指数β是当量比ϕ的函数,其经验公式如下[13]:
{γ=2.18−0.8(ϕ−1) β=−0.16+0.22(ϕ−1) (6) 表 2 参数BM,B2和ϕM取值Table 2. Value of BM, B2 and ϕM燃料 ϕM BM/(cm·s-1) B2/(cm·s-1) 甲醇 1.11 36.92 -140.51 丙烷 1.08 34.22 -138.65 异辛烷 1.13 26.32 -84.72 RMFD-303 1.13 27.58 -78.34 根据表 2,选异辛烷的BM,B2和ϕM参考值来估算SL。根据式5(b),SL, ref=24.89 cm/s;油气体积分数为1.65%时当量比ϕ=1,则根据式(6),γ和β为2.18和-0.16。因此根据表 1,SL=0.253 7 m/s。
层流火焰厚度的计算采用D.B.Spalding的计算式[14]:
{{\rm{ \mathsf{ δ} }}_{\rm{L}}} = \frac{{2\alpha }}{{{S_{\rm{L}}}}} (7) 式中:α为热扩散率,用油气爆炸50 ms时平均温度0.5(Tb+Tu)=713.7 K下空气的物性[15]来近似,并进行压力修正,得α=α0(pref/p)=1.04×10-4 m2/s,α0为空气的热扩散率。因此,δL=8.20×10-4 m。由此可得,τchem=3.23×10-3 s。根据式(1),Da=0.19。
油气混合物的动力粘度μ仍以空气在Tb下的动力粘度μb来近似,根据文献[15],μb=4.56×10-5 N·s/m2,则湍流雷诺数
R{e_{{l_0}}} = \frac{{\rho v_{{\rm{rms}}}^*{l_0}}}{\mu } = \frac{{\left( {p/R{T_{\rm{b}}}} \right)v_{{\rm{rms}}}^*{l_0}}}{{{\mu _b}}} = 4\;603 (8) 式中:R为理想气体常数。
同理,计算出油气爆炸150和250 ms时的丹姆克尔数Da和湍流雷诺数Rel0,如表 3所示。
表 3 各时刻的Da和Rel0计算值(ϕ=1)Table 3. Calculated values of Da & Rel0 at different times (ϕ=1)t/ms Da Rel0 50 0.19 4 603 150 0.26 31 366 250 21 4 012 表 4~5分别是油气当量比为0.72和1.28时激波管内50、150和250 ms时刻气体流速、压力和已燃气体温度的实验数据。根据上文的方法和步骤,当油气当量比分别为0.72和1.28时50、150和250 ms时刻的丹姆克尔数Da和湍流雷诺数Rel0计算结果如表 6所示。
表 4 各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=0.72)Table 4. Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=0.72)t/ms v/(m·s-1) vrms* p/Pa Tmax/K 50 65.13 45.45 115 120 1 158.42 150 -5.74 25.42 322 410 1 201.77 250 11.04 8.64 357 380 1 238.51 表 5 各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1.28)Table 5. Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1.28)t/ms v/(m·s-1) vrms* p/Pa Tmax/K 50 29.69 7.10 11 510 1 251.63 150 86.91 64.32 227 780 1 296.15 250 12.64 9.95 504 310 1 357.36 表 6 各时刻的Da和Rel0计算值Table 6. Calculated values of Da & Rel0 at different timest/ms Da Rel0 ϕ=0.72 ϕ=1.28 ϕ=0.72 ϕ=1.28 50 0.03 1.31 6 767 93 150 0.10 0.15 9 993 15 838 250 0.31 1.76 3 589 4 970 将表 3和表 6中的数据,绘入丹姆克尔数-湍流雷诺数图中,如图 5所示。由图 5可以看出,油气当量比分别为0.72, 1.00和1.28时,油气爆炸50、150和250 ms时的火焰燃烧模式均落入了漩涡内小火焰模式区域。也就是说,低中高初始油气浓度下,油气爆炸初期、充分发展和后期的火焰燃烧模式均为漩涡小火焰模式。这一结论与D.R.Ballal等[16-17]对预混丙烷火焰燃烧模式的实验研究结论一致,他们通过实验发现大部分受限流动中的预混丙烷火焰燃烧模式与漩涡内小火焰模式相关。
图 6是当初始油气体积分数为1.65%时油气爆炸在148、152和156 ms时刻火焰的高速摄影照片。由图 6可以清楚地看出油气爆炸火焰前锋产生了明显的褶皱,并且随着油气爆炸的发展,火焰前锋后生成了许多小火焰点。这些小火焰点很可能与火焰前锋褶皱的不断加强破碎过程有关。
在漩涡内小火焰模式下,油气爆炸火焰的燃烧支持漩涡破碎模型的理论[18],即认为燃烧速度取决于未燃气体破碎成更小微团的速度,由于不断地破碎,使得未燃混合物与已燃热烟气之间有足够的界面进行反应[4],这表明,并不是化学反应的速率决定着燃烧速度,而是湍流混合速度控制着燃烧过程。这一点在对受限空间油气爆炸火焰燃烧的数值分析模型进行选择时至关重要,该结果为受限空间油气爆炸火焰燃烧的数值分析模型的选择提供了依据。
3. 结论
本文首先根据实验先得到了低、中、高3种初始油气浓度条件下油气爆炸过程的气体流速、压力以及火焰温度。通过对实验数据的分析和处理,分别确定了低、中、高3种初始油气浓度条件下油气爆炸过程初期(50 ms)、中期(150 ms)和后期(250 ms)时的气体脉动速度、压力、已燃气体温度和未燃气体温度,最后根据丹姆克尔数和湍流雷诺数对低、中、高3种初始油气浓度条件下油气爆炸过程中的初期、中期和后期火焰燃烧模型分别进行了定量估计。结果表明,低、中、高3种初始油气浓度条件下,激波管油气爆炸过程初期、中期和后期的火焰燃烧模式均为漩涡内小火焰模式;在漩涡内小火焰模式下,油气爆炸过程中火焰的燃烧支持漩涡破碎模型理论,即燃烧速度取决于未燃气体破碎成更小微团的速度,这也表明并不是化学反应的速率决定着燃烧速度,而是湍流混合速度控制着燃烧过程。本文的研究方法和结论可为进一步研究油气爆炸火焰燃烧模式及其数值分析模型的选择提供的参考。
-
表 1 各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1)
Table 1. Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1)
t/ms v/(m·s-1) vrms* p/Pa Tmax/K 50 55.82 34.62 98 530 1 127.39 150 72.76 51.56 340 800 1 203.62 250 29.14 7.94 486 190 1 324.05 表 2 参数BM,B2和ϕM取值
Table 2. Value of BM, B2 and ϕM
燃料 ϕM BM/(cm·s-1) B2/(cm·s-1) 甲醇 1.11 36.92 -140.51 丙烷 1.08 34.22 -138.65 异辛烷 1.13 26.32 -84.72 RMFD-303 1.13 27.58 -78.34 表 3 各时刻的Da和Rel0计算值(ϕ=1)
Table 3. Calculated values of Da & Rel0 at different times (ϕ=1)
t/ms Da Rel0 50 0.19 4 603 150 0.26 31 366 250 21 4 012 表 4 各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=0.72)
Table 4. Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=0.72)
t/ms v/(m·s-1) vrms* p/Pa Tmax/K 50 65.13 45.45 115 120 1 158.42 150 -5.74 25.42 322 410 1 201.77 250 11.04 8.64 357 380 1 238.51 表 5 各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1.28)
Table 5. Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1.28)
t/ms v/(m·s-1) vrms* p/Pa Tmax/K 50 29.69 7.10 11 510 1 251.63 150 86.91 64.32 227 780 1 296.15 250 12.64 9.95 504 310 1 357.36 表 6 各时刻的Da和Rel0计算值
Table 6. Calculated values of Da & Rel0 at different times
t/ms Da Rel0 ϕ=0.72 ϕ=1.28 ϕ=0.72 ϕ=1.28 50 0.03 1.31 6 767 93 150 0.10 0.15 9 993 15 838 250 0.31 1.76 3 589 4 970 -
[1] Steinberg A M, Driscoll J F. Straining and wrinkling processes during turbulence-premixed flame interaction measured using temporally-resolved diagnostics[J]. Combustion and Flame, 2009, 156(12):2285-2306. doi: 10.1016/j.combustflame.2009.06.024 [2] Shin D H, Lieuwen T. Flame wrinkle destruction processes in harmonically forced, laminar premixed flames[J]. Combustion and Flame, 2012, 159(11):3312-3322. doi: 10.1016/j.combustflame.2012.06.015 [3] Yi Y, Geng L, Jing G. Experimental study on the fractal characteristic of methane explosion flame[J]. Safety Science, 2012, 50(4):679-683. doi: 10.1016/j.ssci.2011.08.052 [4] Won S H, Windom B, Jiang B, et al. The role of low temperature fuel chemistry on turbulent flame propagation[J]. Combustion and Flame, 2014, 161(2):475-483. doi: 10.1016/j.combustflame.2013.08.027 [5] Zhang M, Wang J, Wu J, et al. Flame front structure of turbulent premixed flames of syngas oxyfuel mixtures[J]. International Journal of Hydrogen Energy, 2014, 39(10):5176-5185. doi: 10.1016/j.ijhydene.2014.01.038 [6] Mukaiyama K, Shibayama S, Kuwana K. Fractal structures of hydrodynamically unstable and diffusive-thermally unstable flames[J]. Combustion and Flame, 2013, 160(11):2471-2475. doi: 10.1016/j.combustflame.2013.05.017 [7] Williams F A. Asymptotic methods in turbulent combustions[J]. AIAA Journal, 1986, 24:867-875. doi: 10.2514/3.9361 [8] Abraham J, Williams F A, Bracco F V. A discussion of turbulent flame structure in premixed charges[R]. Warrendale, Pennsylvania, USA: Society of automotive engineers, 1985. [9] Zhang P, Du Y, Zhou Y, et al. Explosions of gasoline-air mixture in the tunnels containing branch configuration[J]. Journal of Loss Prevention in the Process Industries, 2013, 26(6):1279-1284. doi: 10.1016/j.jlp.2013.07.003 [10] Yang D, Li Z P, Hong O Y. Effects of humidity, temperature and slow oxidation reactions on the occurrence of gasoline-air explosions[J]. Journal of Fire Protection Engineering, 2013, 23(3):226-238. doi: 10.1177/1042391513486464 [11] 杨辉, 崔鑫, 郑昕等.管道中湍流强度及湍流积分尺度随时间的变化研究[J].安全与环境工程, 2013, 20(4):102-104. doi: 10.3969/j.issn.1671-1556.2013.04.023Yang H, Cui X, Zheng X, et al.Investigation of variation of turbulence intensity and turbulent integral scale with time in tube[J]. Safety and Environmental Engineering, 2013, 20(4):102-104. doi: 10.3969/j.issn.1671-1556.2013.04.023 [12] Sak C, Liu R, Ting D S-K, et al. The role of turbulence length scale and turbulence intensity on forced convection from a heated horizontal circular cylinder[J]. Experimental Thermal and Fluid Science, 2007, 31(4):279-289. doi: 10.1016/j.expthermflusci.2006.04.007 [13] Metghalchi M, Keck J C. Burning velocities of mixtures of air with methanol, isooctane, and indolene at high pressure and temperature[J]. Combustion and Flame, 1982, 48:191-210. doi: 10.1016/0010-2180(82)90127-4 [14] Spadling D B. Combustion and Mass Transfer[M]. New York: Pergamon, 1979:59. [15] Stephen R T. An introduction to combustion:concepts and application[M]. New York: The McGraw-Hill Companies, 2000:appendix C. [16] Ballal D R, Lefebvre A H. The structure and propagation of turbulent flames[J]. Proceedings of the Royoal Society of London, 1979, 344(1637):217-234. http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_1312.1890 [17] Ballal D R, Lefebvre A H. The structure of a premixed turbulent flame[J]. Proceedings of the Royoal Society of London, 1979, 367(1730):253-280. http://www.jstor.org/stable/79908 [18] Poinsot T, Candel S, Trouvé A. Applications of direct numerical simulation to premixed turbulent combustion[J]. Progress in Energy and Combustion Science, 1995, 21(6):531-576. doi: 10.1016/0360-1285(95)00011-9 期刊类型引用(6)
1. 张培理,肖俊,王建,李玉玺. 分支管道长度和数量对油气爆炸超压特性的影响. 安全与环境学报. 2022(05): 2412-2419 . 百度学术
2. 蒋新生,谢威,赵亚东,李静野,李进,余彬彬. 不同长径比的狭长管道油气爆炸实验. 油气储运. 2020(08): 879-884 . 百度学术
3. 蒋新生,谢威,魏树旺,徐建楠,周毅. 有无泄压条件下的管道油气爆轰实验研究. 振动与冲击. 2018(16): 261-266 . 百度学术
4. 张培理,齐圣,王世茂,程顺国. 93号汽油蒸气绝热燃烧温度的计算分析. 后勤工程学院学报. 2017(02): 33-39 . 百度学术
5. 韦世豪,杜扬,王世茂,李蒙. 非均匀容积式受限空间油气爆炸超压与火焰特征. 后勤工程学院学报. 2017(04): 30-36 . 百度学术
6. 韦世豪,杜扬,王世茂,李蒙. 不同形状受限空间内油气爆燃特性的实验研究. 中国安全生产科学技术. 2017(05): 41-47 . 百度学术
其他类型引用(2)
-