4. 围岩与风流热交换原理
在矿井高温热害计算中，围岩是最主要、也是最“隐蔽”的热源。风流穿过井筒、大巷和采区巷道时，并不直接接触原始岩体，而是通过巷道壁面与围岩发生换热。要算清风流温升，必须回答三个问题：围岩内部温度如何分布（温度场）、壁面与风流之间换热强弱如何量化（不稳定换热系数）、以及用什么模型把二者与风流温升联系起来。本节围绕这三条主线展开，并说明工程计算中常用的简化假设及其适用边界。
一、围岩温度场：从原始岩温到壁面温降
未受扰动的深部岩体处于原始岩温状态，其数值由地表恒温带温度与地温梯度决定，常用线性近似：tr0 = ts + g·H，其中 ts 为恒温带温度（约 10～15 ℃），g 为地温梯度（常见 0.02～0.035 ℃/m），H 为埋深。例如取 ts=12 ℃、g=0.03 ℃/m、H=800 m，则原始岩温约 36 ℃，这正是深井热害的根源。
巷道开掘后，壁面被风流冷却，围岩内部形成随时间向深部推进的温度降。以巷道中心为轴、径向坐标为 r，在“准稳态”近似下，围岩温度场可写为：
t(r,τ) = t_r0 - (t_r0 - t_b) · [ln(r/r_0)] / [ln(r_f(τ)/r_0)] (r_0 ≤ r ≤ r_f)式中 tb 为壁面温度，r0 为巷道等效半径，rf(τ) 为“冷却影响半径”，随通风时间 τ 增大而向外扩展，近似满足 rf ∝ √(aτ)，a 为围岩导温系数（约 1.0×10⁻⁶ m²/s 量级）。物理含义很直观：通风初期影响圈很薄，壁面附近温度梯度大、放热强烈；通风时间越长，影响圈越厚，温度梯度变缓，放热量逐渐衰减。这一“随时间衰减”的特性，正是不稳定换热系数要刻画的核心。
二、不稳定换热系数：把复杂温度场压缩成一个参数
若直接按牛顿冷却定律 q = α(tb − tf) 计算壁面热流，需要知道壁面温度 tb，而 tb 又取决于围岩内部导热，二者耦合。工程上更常用的做法是改用“原始岩温—风流温度”为驱动温差，引入不稳定换热系数 Kτ：
q = K_τ · (t_r0 - t_f) [W/m²]
K_τ = 1 / [ (1/α) + (r_0/λ) · ln(r_f(τ)/r_0) ]其中 α 为壁面与风流之间的对流换热系数（W/(m²·K)），λ 为围岩导热系数（W/(m·K)）。该式的结构说明：Kτ 由“对流热阻 1/α”与“围岩导热热阻”串联而成，通风时间越长，导热热阻越大，Kτ 越小。示例取值：α=8 W/(m²·K)、λ=2.5 W/(m·K)、r₀=2 m，通风 1 年时 r_f≈6 m，则导热热阻项约 (2/2.5)·ln3≈0.88，1/α=0.125，得 Kτ≈1.0 W/(m²·K)；若通风 5 年 r_f≈13 m，则 Kτ 降至约 0.55 W/(m²·K)。可见长期通风巷道围岩放热显著减弱，这一趋势对风温预测至关重要。
三、热交换计算模型：从壁面热流到风流温升
把不稳定换热系数用于一段长度为 dx 的巷道微元，可建立风流能量平衡。设风量为 Q（m³/s）、空气密度 ρ、比热容 cp，微元周长为 U，则：
ρ·Q·c_p·dt_f = K_τ·U·(t_r0 - t_f)·dx
积分得：t_f(x) = t_r0 - (t_r0 - t_f0)·exp[ -K_τ·U·x / (ρ·Q·c_p) ]该指数模型表明：风流温度沿巷道按指数规律趋近原始岩温，趋近速度由综合参数 KτU/(ρQcp) 决定。工程上常定义“热交换准数”或直接计算温升：
Δt_f = (t_r0 - t_f0)·[1 - exp(-K_τ·U·L/(ρ·Q·c_p))]示例：某大巷 L=500 m、U=12 m、Q=20 m³/s、ρ=1.2 kg/m³、cp=1005 J/(kg·K)、Kτ=1.0 W/(m²·K)、tr0=36 ℃、入口风温 26 ℃。计算指数项：KτUL=1.0×12×500=6000 W/K；ρQcp=1.2×20×1005≈24120 W/K；比值≈0.249，1−e−0.249≈0.220，故 Δtf≈10×0.220≈2.2 ℃，出口风温约 28.2 ℃。若风量提高到 30 m³/s，比值降至 0.166，温升约 1.5 ℃——这正是“加大风量降温”的定量依据。
四、简化假设及其适用边界
上述模型之所以可用，依赖一组明确的简化假设，使用时必须清楚其边界：
- 围岩均质各向同性、物性恒定：λ、a 不随温度与岩性变化；层状或含水岩层需分段取值。
- 原始岩温按地温梯度线性分布：忽略局部地热异常、岩浆岩体与断层带的影响。
- 巷道断面等效为圆形：用等效半径 r₀ = √(A/π) 或水力半径折算，非圆形断面误差一般在 10%～20%。
- 壁面温度沿周向均匀、沿轴向缓变：忽略局部淋水、支护结构与壁面粗糙度的差异。
- 风流沿巷道充分混合、物性不变：不计压缩功、水分蒸发潜热与风机温升，这些需另行叠加。
- 准稳态导热：忽略围岩热容量引起的瞬态项，适用于通风时间远大于热扩散时间尺度的巷道。
综上，围岩温度场描述“热从哪里来、如何分布”，不稳定换热系数把复杂导热过程压缩为一个随通风时间衰减的参数，指数型热交换模型则把壁面热流转化为风流沿程温升。三者串联，构成了矿井风流温度计算中最基础、也最常用的一环。掌握各简化假设的适用条件，才能在实际矿井中既算得快、又算得准。