一种采空区自然发火的预测方法
【专利摘要】本发明涉及一种采空区自然发火的预测方法，其步骤：设置静坐标下采空区自然发火模型；利用经坐标下采空区自然发火模型及采空区动态移动特点，引入移动坐标系建立移动坐标下采空区自然发火模型；在移动坐标下采空区自然发货模型基础上获得无因次的通用采空区自燃模型；根据无因次的通用采空区自燃模型获取自燃预测准则：耗氧准则和放热准则；求出耗氧准则和放热准则的取值范围，并利用无因次的通用采空区自燃模型对两准则取值范围内的不同取值组合进行自燃情况模拟，将模拟得到的采空区最高温度绘制成曲线得到采空区自燃判别图，完成对采空区自然发火的预测。本发明可广泛用于矿井生产中采空区自然发火的预测工作中。
【专利说明】
一种采空区自然发火的预测方法
技术领域
[0001] 本发明涉及一种煤矿采空区自燃预测方法，特别是关于一种煤矿采空区自燃预测 领域中的采空区自然发火的预测方法。
【背景技术】
[0002] 目前，采空区自然发火一直是我国煤矿的主要灾害之一，随着社会的发展和科技 的进步，采空区煤炭自燃预测的方法也越来越多。目前，采空区自然发火预测方法主要有： 测温法、气体分析法、数值模拟法等。测温法比较直观，受外界因素影响小，但是对预测预报 的要求很高，特别是传感器的选择，这就必然要增加经济成本，而且不排除出现误报的可 能。气体分析法，由于不同矿井，不同煤质也会有所不同，这给指标气体的选择和量化带来 一定的困难。数值模拟法，具有成本较低的优点，是进行采空区自燃预测预报的常用方法之 一。现有方法主要存在以下两方面问题:一是这些模型均是在静坐标下建立的，在考虑工作 面推进的情况时，虽然建立了数学模型，但由于采煤工作面移动，数学模型必然是移动边界 条件下的非稳态数学模型，因而过于复杂，以至不能对模型求解，或者提出大量假设对模型 进行化简，使得解算结果与实际情况相差太大。二是在没有对模型中各个物理量及其之间 的关系进行分析的前提下直接利用数值方法求解模型，固然该方法可用以解算采空区的温 度场、氧浓度场等相关数值，但它只能针对某个具体的采场空间，不能分析各个影响因素之 间的相互关系以及它们对自燃过程的影响，因此无法得到科学的判别自然发火指标。
[0003] 从以上列举的几种方法可以看出，每种方法都有各自的优缺点。因此，完善这些方 法，使之成为准确、实用的预测方法迫在眉睫。
【发明内容】
[0004] 针对上述问题，本发明的目的是提供一种更准确、实用的采空区自然发火的预测 方法。
[0005] 为实现上述目的，本发明采取以下技术方案:一种采空区自然发火的预测方法，其 特征在于，该预测方法包括以下步骤:1)设置静坐标下采空区自然发火模型;2)利用经坐标 下采空区自然发火模型及采空区动态移动特点，引入移动坐标系建立移动坐标下采空区自 然发火模型;3)在移动坐标下采空区自然发货模型基础上获得无因次的通用采空区自燃模 型;4)根据无因次的通用采空区自燃模型获取自燃预测准则:耗氧准则和放热准则;5)求出 耗氧准则和放热准则的取值范围，并利用无因次的通用采空区自燃模型对两准则取值范围 内的不同取值组合进行自燃情况模拟，将模拟得到的采空区最高温度绘制成曲线得到采空 区自燃判别图，完成对采空区自然发火的预测。
[0006] 优选地，所述步骤1)中，流场模型及边界条件：
[0007]
[0008]
[0009] 式中，Kx为采空区内浮煤的渗透系数，单位为m/s;p为空气静压和速压之和，单位 为Pa; Pg为采空区内空气的密度，单位为kg/m3; Po为采空区工作面空气的密度，单位为kg/ m3;g为重力加速度，单位为m/s2;a为煤层的倾角;η为采空区内的空隙率;τ为时间，单位为 s; vx、vy分别为沿X方向和y方向的渗流速度，单位为m/s; Γ^、Γ2、Γ3、Γ4为采空区边界，其中 为沿工作面切顶线侧，Γ2为采空区上部边界，Γ3为采空区下部边界，Γ4为采空区深部 边界;P (X，y)为沿r i边界上的风压函数，单位为Pa。
[0010] 优选地，所述步骤1)中，氧浓度场模型及边界条件：
[0011]
[0012]
[0013] 式中，C为氧气浓度，单位为mol/m3; d( 〇2)为氧气在空气中的扩散系数，单位为m2/ s;u(t)为松散煤体内的氧气消耗速，单位为mol/(s · m3) ;c(x，y)S Γ :边界下部的氧气浓度 函数，将Γ :边界分成上部Γ ls、下部Γ lx两部分。
[0014] 优选地，所述步骤1)中，温度场模型及边界条件：
[0015]
[0016]
[0017]式中，PeCFa-rOPsCs+nPgCwP^PhPj别为松散煤体、固体和气体的密度，单位 为kg/m3;Ce、Cs、匕分别为松散煤体、固体、气体的热容，单位为kJ/(kg · K);Ae = (I-n)As+n Ag入、\义分别为松散煤体、固体、气体的导热系数，单位为kj/(m· s ·Κ);?为温度，K;q (t)为单位时间、单位体积内固体颗粒的放热量，单位为kJ/(mol · 部的温度函数，单位为κ; r5~r1Q为假想的温度场边界。
[0018]优选地，所述步骤3)中，流场模型的无因次判别:设回采工作面长度L为定型尺寸， 工作面两端的压力差Δρο,其中回风侧的压力为PQ，工作面的推进速度VQ，进风流中的空气 密度PQ为基准密度，采空区的渗透系数为K x = Kof (V )，其中Ko的量纲与速度的量纲相同，故 取Ko为基准速度；
[0019] 在对移动坐标下的采空区自燃流场模型引入以下无因次准数:无因次坐标准数：
;无因次压力准数：
：无因次密度准数
无因次风速 准数
6因次压力常数准数：
[0020] 优选地，所述步骤3)中，氧浓度场模型的无因次判别:设采空区进风流中的氧气浓 度CQ为基准氧浓度，采空区进风流温度为t〇，煤升温氧化的临界温度为tc，△ tc = tc 一 t〇 ;在
对采空区自燃氧浓度场模型定义如下无因次准数:无因次氧气浓度准数= 无因次 温度准数
无因次氧气扩散系数准数： ￡因次耗氧速度系
~ a ~U c 数a准数
；无因次耗氧速度系数b准娄
[0021 ]优选地，所述步骤3)中，温度场模型的无因次判别:设松散煤体的导温系数
(Pe3Q3)Aag = W(PgCg);在对采空区自燃温度场模型再定义如下的无因次准数:采空区浮 煤的傅立叶数准数：
：;热容比准数
无因次放热强度系数e准数：
;无因次放热强度系数f准数：
；无因次工作面推进速度 准数：
[0022] 优选地，所述步骤4)中，耗氧准则Ro为：
[0023]
[0024]式中，MaS无因次耗氧速度系数a准数;L为回采工作面长度；Atc = tc - tQ，tQ为采 空区进风流温度，t。为煤升温氧化的临界温度;kb为与煤的比表面积有关的系数;kh为与煤 的厚度有关的系数;a为无因次耗氧速度系数;Co为进风风流中的氧气浓度;Ko为基准速度。 [0025]优选地，所述步骤4)中，放热准则Rt为：
[0026]
[0027]式中，E为无因次推进速度;Fos为实煤的傅立叶数;Ne为无因次放热强度系数;L为 回采工作面长度;kb为与煤的比表面积有关的系数;kh为与煤的厚度有关的系数;e为无因次 放热强度系数;VQ为工作面的推进速度；A = - to，to为采空区进风流温度，为煤升温 氧化的临界温度;Ce3为松散煤体的热容，单位为kj/(kg · K); Pe3为松散煤体的密度，单位为 kg/m3。
[0028] 优选地，所述步骤5)中，当采空区的最高无因次温度达到或超过1时，采空区就有 自然发火危险；当采空区的最高无因次温度小于1时，采空区不具有自然发火危险。
[0029] 本发明由于采取以上技术方案，其具有以下优点：本发明根据采空区动态移动的 特点，通过引入移动坐标系建立移动坐标下的采空区自然发火模型，简化了自然发火模型 求解，并对移动坐标下的采空区自然发火模型进行无因次分析。依据对模型的分析选取回 采工作面长度作为定型尺寸，选取基准速度、基准氧气浓度和基准温度等参考值，由这些参 考值定义了无因次压力、无因次氧扩散系数、傅立叶准数等无因次准则。以这些无因次准则 为基础建立了移动坐标下无因次的采空区自然发火模型。通过分析影响采空区自然发火的 无因次参数，解算不同无因次组合下的采空区自然发火的情况，得到了采空区的最高温度， 并绘制了采空区自然发火的判别曲线图。对于一个具体的采空区，可以通过查图方法来确 定它的自然发火可能性，该方法较为准确、实用。本发明可广泛用于矿井生产中采空区自然 发火的预测工作中。
【附图说明】
[0030]图1是本发明的采空区自然发火边界示意图；
[0031 ]图2是本发明的采空区自燃判别示意图。
【具体实施方式】
[0032]下面结合附图和实施例对本发明进行详细的描述。
[0033] 本发明提供一种采空区自然发火的预测方法，该方法为查图法，其包括以下步骤：
[0034] 1)设置静坐标下采空区自然发火模型：
[0035] 1.1)流场模型及边界条件：
[0036]
[0037]
[0038] 式中，Kx为采空区内浮煤的渗透系数，单位为m/s;p为空气静压和速压之和(不包 括位压），单位为Pa;p g*采空区内空气的密度，单位为kg/m3;Po为采空区工作面空气的密 度，单位为kg/m 3;g为重力加速度，单位为m/s2;α为煤层的倾角;η为采空区内的空隙率;τ为 时间，单位为s ;νχ、Vy分别为沿X方向和y方向的渗流速度，单位为m/s; Γ ^ Γ 2、Γ 3、Γ 4为采 空区边界，其中为沿工作面切顶线侧，Γ2为采空区上部边界，Γ 3为采空区下部边界，Γ 4 为采空区深部边界;P(x，y)为沿Γ:边界上的风压函数，单位为Pa。
[0039] 1.2)氧浓度场模型及边界条件：
[0040]
⑶
[0041]
(4)
[0042] 式中，c为氧气浓度，单位为mol/m3; d( 〇2)为氧气在空气中的扩散系数，单位为m2/ s; u (t)为松散煤体内的氧气消耗速，单位为mo I /(s · m3 )，它是与温度和氧气浓度有关的函 数，通过煤的氧化升温实验测得，并用线性回归成2段直线(分临界温度前和临界温度后）， 用 U(t)=Ckbkh(at+b)/CQ方程来表示，k b为与煤的比表面积有关的系数，kh为与煤的厚度有 关的系数，Co为进风风流中的氧气浓度，a、b分别为回归系数;c(X，y)为Γ :边界下部的氧气 浓度函数，由于采空区1^边界上一部分向采空区内漏风，一部分是从采空区向外漏风，所 以将Γ i边界分成上部Γ ls、下部Γ lx两部分，通过流场解算结果来判定上部Γ ls、下部Γ lx两 部分，通过流场解算结果来判定上部、下部的分界。
[0043] 1.3)温度场模型及边界条件：
[0044]
[0045]
[0046] 式中，PeCeia-rOPsCs+nPgC^PhPhPg*别为松散煤体、固体和气体的密度，单位 为kg/m3;C e、Cs、匕分别为松散煤体、固体、气体的热容，单位为kJ/(kg · K);Ae = (I-n)As+n Ag入、心入分别为松散煤体、固体、气体的导热系数，单位为kj/(m · s · K) ;t为温度，单位 为K;q(t)为单位时间、单位体积内固体颗粒的放热量，单位为kJ/(mol · s)，也是与温度和 氧气浓度有关的函数，通过煤的氧化升温实验测得，并用线性回归成2段直线（分临界温度 前和临界温度后），用q(t) = Ckbkh(et+f)/CQ方程来表示，e，f分别为回归系数。t(x，y)S「r 边界下部的温度函数，κ; r5~r1Q为假想的温度场边界，考虑到热量的传递不但在采空区 范围内进行，还要与周围围岩换热，所以假定了热量传递的范围（如图1所示）。
[0047] 2)利用经坐标下采空区自然发火模型及采空区动态移动特点，引入移动坐标系建 立移动坐标下采空区自然发火模型。
[0048] 引入动态坐标系，/)，动坐标随工作面推进而移动，动坐标与静坐标(x，y)的 关系为：
[0049]
[0050] 式中，V为工作面平均推进速度，m/s。
[0051 ]由式(7)可得：
[0052]
[0053] 式中，供为通用变量，可以代表压力、温度、氧气浓度。
[0054]将式(8)代入式（1)、式(3)和式(5)中并整理得移动坐标下的采空区自然发火数学 模型为：
[0055] 2.1)流场模型：
[0056]
[0057]
[0058]
[0059]
[0060]
[0061]在以上流场方程式(9)和氧浓度场方程式（10)中，由于回采工作面的推进速度远 小于采空区内的空气流动速度，故在两个方程中舍去了涉及到工作面推进速度的项。
[0062]相应的3个边界条件变为：
[0063] 2.1)流场：
[0064]
[0065]
[0066]
[0067]
[0068]
[0069] 由此可知，由流场模型式(9)和式（12)、氧浓度场模型式（10)和式（13)以及温度场 模型式(11)和式(14)共同组成了完整的移动坐标下采空区自然发火模型。
[0070] 3)在移动坐标下采空区自然发货模型基础上获得无因次的通用采空区自燃模型： 依据对移动坐标下采空区自然发火模型，选取回采工作面长度作为定型尺寸，选取基准速 度、基准氧气浓度和基准温度等参考值，由这些参考值定义无因次压力、无因次氧扩散系数 以及傅立叶准数等无因次准则;再以这些无因次准则为基础建立无因次的通用采空区自燃 模型。其具体过程如下：
[0071] 3.1)流场模型的无因次判别
[0072]设回采工作面长度L为定型尺寸，工作面两端的压力差Δρο,其中回风侧的压力为 PO，工作面的推进速度VO，进风流中的空气密度PO为基准密度，采空区的渗透系数为Kx = Kof (V )，其中Ko的量纲与速度的量纲相同，故取Ko为基准速度。
[0073] 对移动坐标下的采空区自燃流场模型式(9)引入以下无因次准数：
[0074] 无因次坐标准数：
[0075] 无因次压力准数
[0076] 无因次密度准数
[0077] 无因次风速准数 [0078]无因次压力常数
[0079] 移动坐标下的流场模型式(9)和边界条件式(12)用以上无因次量表示，得到：
[0080]
[0081]
[0082] 3.2)氧浓度场模型的无因次判别
[0083]设采空区进风流中的氧气浓度Co为基准氧浓度，采空区进风流温度为to,煤升温氧 化的临界温度为tc，Atc； = t。一to。在对采空区自燃氧浓度场模型式（10)定义如下无因次准 数：
[0084] 无因次氧气浓度
[0085] 无因次温度准数
[0086] 无因次氧气扩散
[0087] 无因次耗氧速度
[0088] 无因次耗氧速度
[0089]将氧浓度场方程式(10)和边界条件式(13)分别用以上无因次量表示，得到：
[0090]
[0091]
[0092] 3.3)温度场模型的无因次判别
[0093]设松散煤体的导温系数ae = Ae/(peCe)，令ag = Ae/(pgCg)。在对采空区自燃温度场模 型式(11)再定义如下的无因次准数：
[0094] 采空区浮煤的傅立叶数准数
[0095] 热容比准数：
[0096] 无因次放热强度系数e准数：
[0097] 无因次放热强度系数f准数：
[0098] 无因次工作面推进速度准数
[0099]将温度场方程式(11)和边界条件式(14)分别用以上无因次量表示，得到：
[0100]
[0101]
[0102] 通过以上分析，则无因次的流场模型式（15)和式（16)、无因次的氧浓度场模型式 (17)和式（18)以及无因次温度场模型式（19)和式(20)共同组成了完整的无因次的通用采 空区自燃模型。
[0103] 4)根据无因次的通用采空区自燃模型获取自燃预测准则，即耗氧准则和放热准 则。
[0104] 4.1)流场参数分析
[0105] 采空区流场模型主要影响到采空区内供风量大小，流场解算采空区压力的目的是 为计算采空区内风速大小，以便在解算采空区氧浓度场和温度场时使用。从式（15)可以看 出，影响采空区内风速大小的无因次参数有孔隙率分布函数f (XL)、无因次密度Ω、压力常 数Ψ。在这三个参数中孔隙率分布变化不大，无因次密度主要是受采空区温度变化影响，属 于内部因素。工作面的压力常数Ψ与工作面两端的压差和工作面的长度的比值有关，这个 比值的变化范围不大，故流场模型内的无因次量对采空区自然发火的影响不大。
[0106] 4.2)氧浓度场参数分析
[0107]从前面式（17)以看出影响氧气浓度的无因次参数主要有无因次氧扩散系数D、无 因次耗氧速度系数准数1和无因次耗氧速度系数准数Mb。在三个参数中无因次氧扩散系数D 不变，而无因次耗氧速度系数Mb约等于0,所以影响采空区氧气浓度的无因次参数主要是无 因次耗氧速度系数M a，称之为耗氧准则Ro，则
[0108]
(21)
[0109] 4.3)温度场参数分析
[0110]从前面式（19)可以看出影响采空区内温度的无因次参数主要有无因次推进速度 E、实煤的傅立叶数F〇s、无因次放热强度系数Ne、无因次放热强度系数Nf和热容比α。在四个 参数中无因次放热强度系数N f约等于0,热容比α变化不大，故影响采空区内温度的无因次 参数主要有无因次推进速度E、实煤的傅立叶数Fos和无因次放热强度系数Ne。将E、F〇s和Ne 重新组合成新的无因次参数，称为放热准则Rt，则
[0111]
(22)
[0112] 通过以上分析，把影响采空区自然发火的实际众多物理参数归结为两个主要无因 次准则，即耗氧准则R0和放热准则Rt，这两个准则统称为采空区自燃准则。对于一个实际回 采工作面的采空区，由与采空区自燃准则相关的实际物理量计算出采空区自燃准则，再由 这些自燃准则值的大小判断出采空区的自然发火情况。
[0113] 5)求出耗氧准则和放热准则的取值范围，并利用无因次的通用采空区自燃模型对 两准则取值范围内的不同取值组合进行自燃情况模拟，将模拟得到的采空区最高温度绘制 成曲线得到采空区自燃判别图。
[0114] 对于采空区自然发火的判断可以采用现有采空区自然发火解算系统来得到采空 区内的温度场分布，从而判断自然发火情况。但这样应用有个不足，就是对每个应用实例都 需要进行解算，不但费时，而且不方便。故本发明对各种情况下的采空区自然发火进行解 算，将解算的结果绘制成图。这样对于一个具体的采空区，可以通过查图方法来确定它的自 然发火可能性。
[0115] 分析影响采空区自然发火各物理量实际的取值范围，得到采空区自燃准则可能的 取值范围。在自燃准则的取值范围内按照它们对自然发火的影响不同分别取四个值。组合 两个自燃准则的不同取值分别进行采空区自然发火的模拟解算，得到不同准则组合下采空 区的最高无因次温度，并绘制成图形。
[0116] 经过分析采空区耗氧准则和放热准则的取值如表1所示。
[0117] 表1无因次准则取值
[0119] 将表中的无因次数值组合分别进行采空区自然发火模拟解算得到它们的不同组 合下采空区最高无因次温度，当采空区的最高无因次温度达到或超过1时，采空区就有自然 发火危险；当采空区的最高无因次温度小于1时，采空区不具有自然发火危险。将组合到的 最高无因次温度绘制成图（如图2所示），图中的横坐标为采空区自燃准则的耗氧准则Ro，图 的纵坐标为解算所得的采空区最高无因次温度，用0 max表示。图中的四条曲线分别为四个放 热准则对应不同Ro值所解算得采空区最高无因次温度的连线。由于傅立叶数的间隔大致为 等比的，所以纵坐标取为对数刻度。
[0120] 上述各实施例仅用于说明本发明，各部件的结构、尺寸、设置位置及形状都是可以 有所变化的，在本发明技术方案的基础上，凡根据本发明原理对个别部件进行的改进和等 同变换，均不应排除在本发明的保护范围之外。
【主权项】
1. 一种采空区自然发火的预测方法，其特征在于，该预测方法包括以下步骤： 1) 设置静坐标下采空区自然发火模型； 2) 利用经坐标下采空区自然发火模型及采空区动态移动特点，引入移动坐标系建立移 动坐标下采空区自然发火模型； 3) 在移动坐标下采空区自然发货模型基础上获得无因次的通用采空区自燃模型； 4) 根据无因次的通用采空区自燃模型获取自燃预测准则:耗氧准则和放热准则； 5) 求出耗氧准则和放热准则的取值范围，并利用无因次的通用采空区自燃模型对两准 则取值范围内的不同取值组合进行自燃情况模拟，将模拟得到的采空区最高温度绘制成曲 线得到采空区自燃判别图，完成对采空区自然发火的预测。2. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤1)中， 流场模型及边界条件：式中，Kx为采空区内浮煤的渗透系数，单位为m/s; p为空气静压和速压之和，单位为Pa; Pg为采空区内空气的密度，单位为kg/m3; pq为采空区工作面空气的密度，单位为kg/m3; g为 重力加速度，单位为m/s2;a为煤层的倾角;n为采空区内的空隙率;t为时间，单位为s ; Vx、Vy 分别为沿X方向和y方向的渗流速度，单位为m/s; h、r2、r3、r4为采空区边界，其中rd 沿工作面切顶线侧，r2为采空区上部边界，r3为采空区下部边界，r4为采空区深部边界;p (X，y)为沿边界上的风压函数，单位为Pa。3. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤1)中， 氧浓度场模型及边界条件：式中，c为氧气浓度，单位为mol/m3; d(〇2)为氧气在空气中的扩散系数，单位为m2/s ; u (t)为松散煤体内的氧气消耗速，单位为mol/(s ? m3) ;c(x，y)为r :边界下部的氧气浓度函 数，将r i边界分成上部r ls、下部r lx两部分。4. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤1)中， 温度场模型及边界条件：式中，P' = (1 -n) psCs+npgCg; pe、ps、Pg分别为松散煤体、固体和气体的密度，单位为kg/ m3; Ce、Cs、(^分别为松散煤体、固体、气体的热容，单位为kj/(kg ? K); = (1 -n)As+nAg入、 入s人分别为松散煤体、固体、气体的导热系数，单位为kj/(m? s ?!〇;〖为温度，K;q(t)为单 位时间、单位体积内固体颗粒的放热量，单位为kJ/(mol ? s);t(x，y)为「:边界下部的温度 函数，单位为K; r 5~r 1Q为假想的温度场边界。5. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤3)中， 流场模型的无因次判别：设回采工作面长度L为定型尺寸，工作面两端的压力差A PQ，其中 回风侧的压力为P〇,工作面的推进速度进风流中的空气密度P〇为基准密度，采空区的渗 透系数为K^Kof^)，其中Ko的量纲与速度的量纲相同，故取Ko为基准速度； 在对移动坐标下的采空区自燃流场模型引入以下无因次准数：无因次坐标准数：6. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤3)中， 氧浓度场模型的无因次判别：设采空区进风流中的氧气浓度co为基准氧浓度，采空区进风 流温度为to,煤升温氧化的临界温度为t c，A tc = tc - to;在对采空区自燃氧浓度场模型定义 如下无因次准数： 无因次氧气浓度准数；无因次温度准数： ；无因次氧气扩 散系数准数：；无因次耗氧速度系数a准数：无因次耗氧速度系数 b准数7. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤3)中， 温度场模型的无因次判别:设松散煤体的导温系数a e = Ae/(peCe)，$ag = Ae/(pgCg);在对采 空区自燃温度场模型再定义如下的无因次准数： 采空区浮煤的傅立叶数准数：；无因次放热强度 系数e准数；无因次放热强度系数f准数：无因次工作 面推进速度准数8. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤4)中， 耗氧准则Ro为：式中，无因次耗氧速度系数a准数;L为回采工作面长度；Atc = tc-to，t()为采空区 进风流温度，t。为煤升温氧化的临界温度;kb为与煤的比表面积有关的系数;kh为与煤的厚 度有关的系数;a为无因次耗氧速度系数;co为进风风流中的氧气浓度;Ko为基准速度。9. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤4)中， 放热准则Rt为：式中，E为无因次推进速度;F〇s为实煤的傅立叶数;Ne为无因次放热强度系数;L为回采 工作面长度;kb为与煤的比表面积有关的系数;kh为与煤的厚度有关的系数;e为无因次放热 强度系数;vo为工作面的推进速度；A tc = tc -to，to为采空区进风流温度，为煤升温氧化 的临界温度;Ce为松散煤体的热容，单位为kj/(kg ? K) ;pe为松散煤体的密度，单位为kg/m3。10. 如权利要求1所述的一种采空区自然发火的预测方法，其特征在于:所述步骤5)中， 当采空区的最高无因次温度达到或超过1时，采空区就有自然发火危险；当采空区的最高无 因次温度小于1时，采空区不具有自然发火危险。
【文档编号】G06F19/00GK106055880SQ201610351586
【公开日】2016年10月26日
【申请日】2016年5月25日
【发明人】朱建芳, 秦跃平, 郭文杰, 王晔, 段嘉敏, 马辉
【申请人】华北科技学院