巷道围岩与风流不稳定换热准数计算方法

Calculation method of unsteady heat-transfer dimensionless number for surrounding rock and airflow in roadway

  • 摘要: 围岩放热是深部高温矿井最主要的热源,其关键参数不稳定换热准数Kuτ直接决定风温预测精度。现有获取Kuτ的方法(解析级数、查图、查表、分段回归、逐巷数值模拟)或计算繁琐,或插值误差大,或无法批量实施,难以满足工程快速、准确、统一计算的需求。为解决上述问题,对巷道围岩与风流不稳定换热准数计算方法展开研究。研究以圆形巷道为对象,建立极坐标下无量纲瞬态导热微分方程,将Kuτ表征为毕渥数Bi与傅里叶数Fo的函数;采用有限体积法离散,节点半径与时间步长按等比级数布置,兼顾剧烈非稳态阶段的精度与后期计算效率,基于Visual Studio自主开发求解程序,在不同Bi数、0.01≤Fo1000的范围内计算高精度离散解;通过变量转换发现(BiKuτ)−1Fo−0.5呈高度线性关系,据此构建统一的二次多项式回归公式。计算结果与经典理论解相吻合,验证了数值模型与程序的正确性;回归公式形式简洁、无分段;回归公式平均相对误差1.8%,最大相对误差 4.2%,优于既有分段回归;算例显示,巷道通风5 d、1 a、10 a后,不稳定换热准数分别降至初始值的29.0%、7.7%、5.0%,与现实规律一致。该研究给出了 BiFo 全域连续变化条件下Kuτ的统一显式表达式,突破了分段拟合的局限,形式简洁,可直接嵌入风温预测软件,实现上百条巷道围岩散热秒级精确计算,为深部矿井热害快速评估及降温设计提供了可靠、高效的工具;该方法可扩展至非圆形断面及非均质围岩,后续将耦合风流温湿过程,实现全矿井热环境动态预测。

     

    Abstract: Heat release from surrounding rock is the dominant heat source in deep high-temperature mines, and its key parameter—the unsteady heat-transfer dimensionless number Kuτ—directly determines the accuracy of air-temperature prediction. Existing methods for obtaining Kuτ (analytical series, charts, tables, piecewise regressions, or case-by-case numerical simulations) are either computationally cumbersome, suffer from large interpolation errors, or cannot be implemented in batches, thus failing to meet engineering demands for rapid, accurate and unified calculations. To address the above issues, a study on the calculation method of the unsteady heat-transfer dimensionless number between roadway surrounding rock and airflow has been carried out. Taking a circular roadway as the research object, a dimensionless transient heat-conduction equation in polar coordinates is established and expresses Kuτ as a function of the Biot number Bi and the Fourier number Fo. The finite-volume method is employed for discretisation; the nodal radii and time steps are arranged in geometric progressions to guarantee accuracy during the violent unsteady stage while maintaining computational efficiency at later times. An in-house solver developed in Visual Studio is used to obtain high-precision discrete solutions over the ranges of different Bi and 0.01 ≤ Fo ≤ 1000. Through a variable transformation, a highly linear relationship between (BiKuτ)−1 and Fo−0.5 is discovered, on the basis of which a unified quadratic-polynomial regression formula is constructed. The numerical results agree well with classical analytical solutions, validating the model and the code. The regression formula is concise and free of piecewise definitions; its average relative error is 1.8% and its maximum relative error is 4.2%, outperforming existing piecewise regressions. A case study shows that the unsteady heat-transfer dimensionless number drops to 29.0%, 7.7% and 5.0% of its initial value after 5 days, 1 year and 10 years of ventilation, respectively, consistent with field observations. A unified explicit expression of Kuτ valid over the entire continuous ranges of Bi and Fo is presented, eliminating the limitations of piece-wise fitting. Owing to its compact form, the formula can be directly embedded in mine air-temperature prediction software to achieve second-level accurate calculation of heat dissipation from surrounding rock for hundreds of roadways, offering a reliable and efficient tool for rapid thermal-hazard assessment and cooling design in deep mines. The method is readily extendable to non-circular cross-sections and heterogeneous rock, and future work will couple air temperature–humidity processes to enable dynamic prediction of the thermal environment of entire mines.

     

/

返回文章
返回