[发明专利]停采状态下采空区自然发火模拟方法及应用

说明书

本发明提供了一种停采状态下采空区自然发火模拟方法及应用，其中，停采状态下采空区自然发火模拟方法，包括以下步骤：1)建立停采状态下采空区自然发火多场耦合数学模型，并确定其边界条件；2)编制停采状态下采空区自然发火多场耦合数学模型求解程序；3)数值模拟获取各物理场的动态演化结果；4)根据各物理场的动态演化结果，确定高氧浓度区以及高温区域，判断采空区是否存在自然发火危险。本发明所述的停采状态下采空区自然发火模拟方法，可以提高停采期间采空区自然发火预报的准确性和可靠性，对于采空区火灾防治具有积极意义。

技术领域

本发明属于煤矿采空区自然发火预测及防治技术领域，尤其是涉及一种停采状态下采空区自然发火模拟方法及应用。

背景技术

煤矿火灾一直是威胁井下安全生产的重大矿山灾害之一，而采空区煤炭自燃则是造成大部分矿井火灾的主要原因。由于采空区过大，并且地质条件极为复杂，氧浓度、温度等自燃指标不能直接测量，难以实施实时的监测，从而难以准确预测采空区煤自燃的具体位置。此外，在停采状态下，采空区也极易发生自燃现象，而目前关于停采状态下的预测采空区自然发火的有效方法少之又少。为了有效地预测停采状态下采空区自然发火变化规律，从而降低矿井危害，提出一种新的停采状态下的采空区自然发火模拟方法是具有重要的现实意义的。

发明内容

有鉴于此，本发明旨在提出一种停采状态下采空区自然发火模拟方法，以提高停采期间采空区自然发火预报的准确性和可靠性。

为达到上述目的，本发明的技术方案是这样实现的：

停采状态下采空区自然发火模拟方法，包括以下步骤：

1)建立停采状态下采空区自然发火的流场、氧浓度场、固体温度场和气体温度场多场耦合数学模型，并确定其边界条件；

2)编制停采状态下采空区自然发火多场耦合数学模型求解程序；

3)数值模拟获取各物理场的动态演化结果；

4)根据各物理场的动态演化结果，确定高氧浓度区以及高温区域，判断采空区是否存在自然发火危险。

优选的，步骤1)中，流场模型为式(1)：

边界条件为：

式中：V为流场中任一控制体的体积；S为流场中任一控制体的面积；ε₀为采空区内浮煤孔隙率，单位为％；ρ_g为空气密度，单位为kg/m³；τ为时间，单位为s；v_g空气渗流速度，单位为m/s；为单位法线向量；为空气的渗流速度的矢量形式，单位为m/s；K为渗透系数，单位为m/s；g为重力加速度，单位为m/s²；P为工作面总压差，单位为Pa；ρ₀为标况下的空气密度，单位为kg/m³；T₀为标况下的温度，单位为℃；T_g为气体的温度，单位为℃；P(x,y,z)是沿工作面的气压函数；Г₁为靠近工作面的边界，属于对流换热边界；Г₂、Г₃为上下两行煤柱，Г₄为采空区深部边界，Г₅为顶板边界，Г₆为底板边界。

优选的，步骤1)中，氧浓度场的模型为式(2)：

边界条件为：