2. 井筒传热基础与温度场控制方程

井筒温度场是油井生产动态分析、水合物防治、结蜡预测、热采效果评价以及分布式光纤测温解释的共同基础。要建立可用的温度剖面，必须先弄清热量在井筒中究竟以什么方式、沿什么方向传递，再把这种传递写成可求解的数学方程。本节从三种基本传热方式入手，逐步过渡到井筒径向与轴向传热机制，最后给出温度场控制方程的建立过程。

一、三种基本传热方式在井筒中的具体表现

热传导是依靠分子热运动、不发生宏观物质位移的热量传递，遵循傅里叶定律：热流密度与温度梯度成正比，方向与温度梯度相反，即 q = −λ·dT/dx，λ 为导热系数（W/(m·K)）。在井筒中，热传导主要发生在三个位置：油管壁、套管壁与水泥环的径向传热；地层岩石内部的径向与垂向传热；以及井筒内静止流体段（如关井后的油套环空液体）的传热。钢材导热系数约 45～50 W/(m·K)，水泥环约 0.8～1.5 W/(m·K)，而地层岩石通常在 1.5～3.5 W/(m·K) 之间，这种数量级差异决定了“薄薄的水泥环”往往是径向热阻的重要贡献者。

对流是流体宏观运动携带热量的传递方式，分为强迫对流与自然对流。牛顿冷却定律给出对流换热量 q = h·(T_w − T_f)，h 为对流换热系数（W/(m²·K)）。井筒内油流向上运动属于强迫对流，其换热强度与流速、流体物性、流态密切相关；油套环空中的液体在温差驱动下会产生自然对流，其等效换热系数通常低于强迫对流。对流是井筒内“把热量从地层带到井口”的主要载体，也是轴向传热的核心机制。

热辐射是以电磁波形式传递能量，服从斯特藩—玻尔兹曼定律 q = εσ(T₁⁴ − T₂⁴)。在常规油井中，油管与套管之间的环空若充满液体或存在管壁屏蔽，辐射换热量占比很小，工程上常并入环空等效热阻处理；但在热采井、注蒸汽井、环空为气体或真空隔热油管的场合，辐射不可忽略，隔热油管正是通过降低环空气体导热与辐射来大幅提高热阻。

二、井筒径向与轴向传热机制

井筒传热可以分解为两个方向。径向传热描述热量从地层深处经水泥环、套管、环空、油管进入井筒流体（或反向）的过程，其路径是“地层 → 水泥环 → 套管壁 → 环空 → 油管壁 → 油管内流体”，每一段都对应一个热阻，串联构成总热阻。工程上常用总传热系数 U 或单位长度热阻 R 来概括这一串联路径，R 越大，井筒与地层之间的热交换越弱，井口温度保持得越好。

轴向传热则描述热量随流体沿井深方向的输运，包括流体本身的对流携带（主导项）、油管壁与流体的对流交换、以及管壁与环空之间的导热。轴向传热的强弱直接决定井口温度：产量越高、流速越大，流体在井筒中停留时间越短，向地层散失的热量比例越小，井口温度越接近井底温度。反之，低产井或关井工况下，轴向对流减弱，径向散热占主导，井口温度会明显下降。

需要特别说明方向约定：本节统一取井深 z 向下为正、温度 T 随深度增加而升高，径向取由井筒中心指向地层为正。后续方程与图示均按此约定，若切换参照系须重新标注符号。

三、温度场控制方程的建立

建立控制方程的基本思路是：在井筒中取一个微元控制体（轴向长度 dz、径向厚度 dr），对控制体做能量守恒——流入热量减去流出热量，加上源项，等于控制体内能的变化。对井筒内流体，取轴向一维、径向与井壁换热的简化模型，可得：

ρ·c_p·(∂T/∂t + v·∂T/∂z) = ∂/∂z(λ_eff·∂T/∂z) + q_wall
q_wall = U·(T_formation − T_fluid) · (2π r / A)

其中 ρ 为流体密度，c_p 为定压比热，v 为流速（向上流动时取负号，与 z 向下为正的约定一致），λ_eff 为有效轴向导热系数（通常远小于对流项，可忽略），q_wall 为通过井壁的径向换热源项，U 为总传热系数，T_formation 为地层温度。稳态条件下 ∂T/∂t = 0，方程退化为常微分方程，可沿井深积分求解。

总传热系数 U 由串联热阻确定：1/U = 1/h₁ + r₁·ln(r₂/r₁)/λ_tube + 1/(h_ann + h_rad) + …，工程上常将其折算为“单位长度热阻 R（m·K/W）”，再乘以井段长度得到该段热阻。地层温度则通常按地温梯度线性近似：T_f(z) = T_surface + g_T·z，g_T 为地温梯度（约 0.025～0.035 ℃/m）。

四、一个可复算的教学示例

设某直井井深 2000 m，地温梯度 0.03 ℃/m，地表温度 15 ℃，则井底地层温度约 75 ℃。取油管内径 0.062 m、产量折算流速 0.5 m/s、流体比热 2200 J/(kg·K)、密度 850 kg/m³，单位长度总热阻取 0.12 m·K/W（示例假设值）。若忽略轴向导热，稳态方程可写为 v·dT/dz = (T_f − T)/ (ρ·c_p·R·A)，沿井深自下而上分段积分，即可得到井口温度。按上述参数粗算，井口温度约比井底低 8～12 ℃，产量越低、热阻越小，温降越大。这一结果属于示例性估算，实际应用需以实测物性与热阻标定为准。

掌握上述传热机制与方程结构后，后续的多相流计算才能把“温度”与“压力、含气率、物性”耦合起来，形成完整的井筒剖面求解框架。