基于量纲分析理论的煤尘爆炸能量预测模型*

李雨成，刘天奇，周西华

(1.辽宁工程技术大学安全科学与工程学院,辽宁 阜新 123000;2.辽宁工程技术大学矿山热动力灾害与防治教育部重点实验室,辽宁 阜新 123000)

基于量纲分析理论的煤尘爆炸能量预测模型*

李雨成1,2，刘天奇1,2，周西华1,2

(1.辽宁工程技术大学安全科学与工程学院,辽宁 阜新 123000;2.辽宁工程技术大学矿山热动力灾害与防治教育部重点实验室,辽宁 阜新 123000)

为预测煤尘爆炸能量，基于量纲分析理论建立煤尘爆炸能量预测模型。选取爆炸能量E、空气密度ρ和大气压强p的量纲为导出量纲。根据量纲分析Π定理得出含有待定参数的具有普适性的能量预测模型。通过小型煤尘爆炸性实验设计，测定10次爆炸最长火焰长度平均值l0、10次最长火焰长度出现时间平均值t0与该小型煤尘爆炸中释放能量E0，确定模型中参数为0.467。对模型变量t、E、l的函数关系进行合理性检验。通过实测的15组不同时刻的火焰长度进行模型变量t、l幂指关系检验。检验结果表明：量纲选取完备，预测模型科学合理。

煤尘爆炸;能量预测;量纲分析;预测模型

煤矿井下作业中，一旦掘进面等处产生的浮尘遇到明火，极易发生爆炸事故[1-2]。煤尘爆炸在传播过程中，冲击波和火焰面释放的能量给巷道结构、生产设备及作业人员带来极大损害[3-4]。煤尘爆炸属于多次连环的有可燃悬浮颗粒参与的化学爆炸。对煤尘爆炸能量的研究主要集中于压力波和火焰波两方面。李庆钊等[5]利用20 L球形装置研究了煤尘浓度、环境初压等因素与爆炸能量的关系。程磊[6]研究了受限分岔管道、双向分岔管道、截面突变管道内煤尘爆炸冲击波的变化规律。张延松[7]、李润之[8]和王新等[9]研究了大型实验巷道中瓦斯诱导煤尘爆炸产生的压力波、火焰波的传播过程。李雨成等[10-11]研究了不同煤质的煤尘爆炸火焰传播特性，并建立了火焰长度预测模型。由于煤尘爆炸能量的数量级较大，不便于测试，因此，本文中基于量纲分析理论建立煤尘爆炸能量预测模型，以期为高数量级的爆炸能量预测提供理论基础。

1 基于量纲分析理论的煤尘爆炸能量预测模型

1.1 基本量纲与导出量纲的选取

量纲是自然界中物理量的单位，分为基本量纲和导出量纲[12-13]。本文中预测煤尘爆炸能量选取的基本量纲为煤尘爆炸火焰长度l、煤尘质量m、时间t的量纲，分别记为：

(1)

在基本量纲的基础上可形成导出量纲，本文中采用的导出量纲有煤尘爆炸能量E、空气密度ρ、大气压强p的量纲，分别记为:

(2)

根据量纲分析理论，可以找到与煤尘爆炸能量相关的变量及其之间的量纲关系。其中，量纲分析齐次原理是量纲分析理论中的重要原理，Π定理是利用量纲分析理论建模的基础，虽然根据Π定理往往不能求出模型的具体表达式，但并不影响建模过程及本文煤尘爆炸能量的预测结果。

1.2 基于量纲分析齐次原理与Π定理的煤尘爆炸能量预测模型

将煤尘爆炸中释放的能量记为E，在本文建模中，考虑的基本量纲有煤尘爆炸火焰传播距离(即火焰长度)l、质量m、时间t的量纲，还考虑了动力学范畴内与之相关的物理变量空气密度ρ与大气压强p的量纲。通过后文模型分析可验证所选取变量是准确充足有效的，其具体量纲表达式已由式(2)给出。则根据量纲分析Π定理，不妨设l与t、E、ρ、p间存在如下关系式l=(t,E,ρ,p)，将其写成隐函数形式为f(l,t,E,ρ,p)=0。根据量纲分析齐次原理中“导出量纲函数可写成基本量纲幂函数乘积形式”的结论，可把l、t、E、ρ、p表示成基本量纲L、M、T幂函数乘积形式:

(3)

根据量纲分析Π定理，由式(3)得到量纲矩阵为:

(4)

再根据Π定理，得到2个相互独立且量纲一量π1和π2，使得F(π1,π2)=0与f(l,t,E,ρ,p)=0等价。其中，π1和π2表达式分别为:

(5)

由于π1和π2均为量纲一量且彼此相互独立，因此两者必然可表达成某一函数ψ形式，即:

(6)

由此得到煤尘爆炸火焰长度l关于t、E、ρ、p的函数关系。在函数ψ中，由于煤尘爆炸火焰存在时间很短，而煤尘爆炸所释放能量很大，因此可做如下处理：

(7)

显然对某一固定函数形式ψ，ψ(0)必为一常数，不妨令ψ(0)=。由此将l推导成关于t、E、ρ的函数关系(空气密度ρ=1.25 kg/m3)如下：

(8)

式中:煤尘爆炸火焰持续时间t、煤尘爆炸火焰长度l均可测，只要再实测一组爆炸释放能量E就可估计出常数，进而利用式(8)即可针对不同时刻t的火焰长度l，估计煤尘爆炸不同时刻释放出的能量E，这对爆炸威力预测及爆炸后果评估具有重要意义。为确定常数，设计实验如下。

2 模型待定参数确定的小型实验设计与模型检验

2.1 小型煤尘爆炸实验设计

图1 煤尘爆炸性鉴定装置Fig.1 Identification equipment of coal dust explosibility

对于模型(8)中存在的3个未知参数l、t、E和一个待定参数，对于不同时刻t会对应存在煤尘爆炸火焰长度l和爆炸释放能量E，但对于时间t存在一个临界点，即为煤尘爆炸能量释放最大(即能量全部释放)的时刻t0，当t>t0时，火焰长度l和释放能量E均开始减小。为确定待定参数，可设计小型实验思路如下。

以t0时刻为基准，通过煤尘爆炸性鉴定实验，实验装置如图1所示，可以测得t0时刻对应的最长火焰长度l0。为便于看清数据，利用高速摄像机捕捉火焰图像，通过玻璃管上标尺刻度读取火焰长度。10次爆炸实验捕捉到最长火焰长度平均值为l0=0.65 m，对应10次实验的最长火焰长度出现的时间的平均值为t0=0.015 s。再根据所选煤样为无烟煤，其煤尘热值Q=2.9×107J/kg，煤样质量为1 g，进而得到该小型煤尘爆炸中释放能量为E0=2.9×104J。

2.2 模型待定参数计算

将2.1节中实验计算得到的l0、t0、E0代入式(8)，可计算得到待定参数为：

(9)

由此得到煤尘爆炸能量释放预测模型为:

(10)

采用该能量预测模型，可以在已知一组爆炸t时刻火焰长度l情况下，预测爆炸各不同时刻能量释放值。但对于该公式存在一定限制条件，即0t0时，火焰长度l开始逐渐衰减，所释放能量E也开始衰减。

2.3 模型检验

2.3.1 模型变量t、E、l合理性检验

对于由量纲分析原理和小型煤尘爆炸性鉴定实验设计推导出的煤尘爆炸能量预测模型(10)，当t=0时，l=0，显然与实际情况相符，即在爆炸起始时刻，煤尘爆炸火焰长度为0，煤尘爆炸释放能量为0。在此基础上，由式(10)可以看出，在t∈[0,t0]范围内，随着爆炸后时间的增大，火焰长度、爆炸释放能量也不断增大，即火焰长度、爆炸释放能量分别与爆炸时间呈幂指递增函数关系，也验证了模型的合理性。

2.3.2 模型变量t、l幂指关系检验

图2 拟合结果Fig.2 Fitted result

再进一步验证模型中幂指关系的合理性，首先暂不考虑变量E的作用，仅从l和t角度可知l正比于t的2/5次幂。为验证该关系，利用SA8独立型高速摄像机(每秒可最多拍摄30 000帧)捕捉的不同时刻火焰长度(为10次测试的平均值)(见表1)，进行l=atb(0

表1 不同时刻的火焰长度

2.4 应用前景分析

在预测模型(10)中，只要已知某时刻t的煤尘爆炸火焰传播距离l，就可以利用该模型预测t时刻爆炸释放的能量。而在煤矿现场，首先可通过技术手段确定爆源位置，在易爆区段安装可精确测量物体运动时间的时间传感器以及火焰传感器，以获取火焰传播最远距离和传播时间数据，以此进行爆炸能量估算。同时，即便测试条件受限，对于煤尘爆炸火焰传播规律的研究(包括火焰传播距离、火焰传播速度与时间的关系)仍然很多，文献[6-7]等均给出火焰传播距离与时间的关系表达式。因此，只要通过爆炸痕迹鉴定判断出火焰传播的最远距离，便可推导出爆炸经历的时间，再代入能量预测模型(10)同样可得出结果。预测模型(10)参数变量较少，量纲选取完备合理，同时使计算过程极大简化，省去以往复杂的偏微分方程的运算过程，降低了在模型求解时对数学理论掌握的高要求。此外，通过量纲分析方法和小型实验室煤尘爆炸性鉴定实验，便可预测计算出数量级相对较高的爆炸能量，大大简化运算量，同时，预测结果对爆炸事故处理有一定参考意义与应用价值。

3 结 论

(1)根据量纲齐次原理与量纲分析Π定理构建了煤尘爆炸能量模型。选取煤尘爆炸火焰传播距离l、质量m、时间t的量纲为基本量纲，选取爆炸能量E、空气密度ρ和大气压强p的量纲为导出量纲。由Π定理得出含有未知参数l、t、E和待定参数的具有普适性和创新性的能量预测模型。

(3)通过分析模型中变量t、E与变量l的函数关系验证了模型的合理性。利用实验实测的15组不同时刻t的火焰长度l检验模型中变量t、l的幂指关系，拟合方程结果为l正比于t的0.400 2次幂，而模型中l正比于t的2/5次幂，进一步验证了量纲选取的完备性和模型的科学性。从应用前景角度分析，该预测模型计算简单，变量参数少，对煤矿爆炸事故危害评估具有借鉴意义和应用价值。

[1] 金龙哲.矿井粉尘防治理论[M].北京:科学出版社,2010:45-56.

[2] 李雨成.矿井粉尘防治理论及技术[M].北京:煤炭工业出版社,2015:23-38.

[3] 蔡周全,罗振敏,程方明.瓦斯煤尘爆炸传播特性的实验研究[J].煤炭学报,2009,34(7):938-941. Cai Zhouquan, Luo Zhenmin, Cheng Fangming. Experimental study on propagation characteristics of gas and coal dust explosion[J]. Journal of Coal, 2009,34(7):938-941.

[4] 来诚锋,段滋华,张永发.煤粉末的爆炸机理[J].爆炸与冲击,2010,30(3):325-328. Lai Chengfeng, Duan Zihua, Zhang Yongfa. Explosion mechanism of carbon powder[J]. Explosion and Shock Waves, 2010,30(3):325-328.

[5] 李庆钊,翟成,吴海进.基于20L球形爆炸装置的煤尘爆炸特性研究[J].煤炭学报,2011,36(S1):119-124. Li Qingzhao, Zhai Cheng, Wu Haijin. Investigation on dust explosion characteristics using 20 L explosion sphere vessels[J]. Journal of Coal, 2011,36(S1):119-124.

[6] 程磊.受限空间煤尘爆炸冲击波传播衰减规律研究[D].焦作:河南理工大学,2011:79-91.

[7] 张延松.瓦斯爆炸诱导沉积煤尘爆炸研究[M].北京:中国矿业大学出版社,2011:19-31.

[8] 李润之.瓦斯爆炸诱导沉积煤尘爆炸的数值模拟[J].爆炸与冲击,2010,30(5):529-534. Li Runzhi. Numerical simulation of coal dust explosion induced by gas explosion[J]. Explosion and Shock Waves, 2010,30(5):529-534.

[91] 王新,李润之,张延松.瓦斯爆炸引起沉积煤尘爆炸传播实验研究[J].中国安全科学学报,2009,19(4):73-77. Wang Xin, Li Runzhi, Zhang Yansong. Experimental study on the propagation laws of gas explosion inducing coal dust explosion[J]. China Safety Science Journal, 2009,19(4):73-77.

[10] 李雨成,刘天奇,陈善乐.煤质指标对煤尘爆炸火焰长度影响作用的主成分分析[J].中国安全生产科学技术,2015,11(3):40-46. Li Yucheng, Liu Tianqi, Chen Shanle. Principal component analysis of impact of coal quality index on flame length in coal dust explosion[J]. Journal of Safety Science and Technology, 2015,11(3):40-46.

[11] Li Yucheng, Liu Tianqi, Chen Shanle. The prediction model research of flame length of coal dust explosion based on coal quality indexes’ PCRA[J]. International Journal of Earth Sciences and Engineering, 2014,7(5):1777-1787.

[12] Bridgman P W. Dimensional analysis[M]. New York: Classics Press, 2013:24-31.

[13] 谈庆明.量纲分析[M].北京:中国科学技术大学出版社,2005:48-59.

(责任编辑 张凌云)

An energy prediction model for coal dust explosion based on dimensional analysis

Li Yucheng1,2, Liu Tianqi1,2, Zhou Xihua1,2

(1.CollegeofSafetyScienceandEngineering,LiaoningTechnicalUniversity,Fuxin123000,Liaoning,China;2.KeyLaboratoryofMineThermo-motiveDisasterandPrevention,MinistryofEducation,LiaoningTechnicalUniversity,Fuxin123000,Liaoning,China)

To assess the coal dust explosive power, we used the dimensional analysis mathematical method to model the prediction of the energyreleased from coal dust explosion. The basic dimensions selected are those of the distance, the mass and the time, of coal dust explosion’s flame propagation. The selected output dimensions are those of explosive energy, the air density and the barometric pressure. According to the dimensional analysis Π theorem, we established the energy prediction model with unknown parameters and undetermined parameters, a model that is universally applicable. By carrying out a small-scale coal dust explosion experimental design, we determined the model parameters, and obtained the average maximum flame length of ten times of coal dust explosion, the average maximum flame time, and the released energy. The undetermined parameter was calculated as 0.467. It follows that the prediction model for coal dust explosive energy and its qualification were obtained. Based on our rationality analysis of this model, the experimental measured flame length of 15 different times of coal dust explosion was used to test the power exponent relationship. The test result verifies the completeness of the selected dimensions and the scientific rationality of the prediction model. Using fewer variable parameters, this model simplifies calculations and will provide important reference for explosion hazard assessment.

coal dust explosion; energy prediction; dimensional analysis; predicting model

10.11883/1001-1455(2017)03-0566-05

2015-09-22；

2015-12-28

国家自然科学基金项目(51204089,51274115);辽宁省教育厅科研项目(L2014131)

李雨成(1978— ),男,博士,副教授,硕士生导师,lyclntu@163.com。

O389;X936 国标学科代码： 13035

A