🔬 physics

A Streamline-Based Local Thermal Non-Equilibrium Framework for Heat Transport in Subsurface Porous Media

本研究引入了一种基于流线的地下热传输局部热非平衡（LTNE）框架，该框架将流体平流与多维传导及相间交换解耦，并通过基准测试和非均匀模拟证明，即使在宏观热响应表现出接近平衡的状态时，流体与岩石之间仍可能存在显著的局部温差。

原作者： Mohammad Nezam Uddin, Yin Feng, Boyun Guo

发布于 2026-07-07

📖 1 分钟阅读☕ 轻松阅读

原作者： Mohammad Nezam Uddin, Yin Feng, Boyun Guo

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

以下是对该研究论文的解释，通过使用日常类比将复杂概念拆解为简单易懂的内容。

大局观： “热岩石，冷水”问题

想象你正试图通过泵入冷水来冷却一个巨大的、滚烫的砖窑（地下岩石）。这就是地热能的工作原理：我们泵入冷水，它吸收岩石的热量，然后我们将热水泵出以产生电能。

长期以来，工程师们一直使用一个简单的经验法则来预测这一过程：“瞬间拥抱”理论（The "Instant Hug" Theory）。
他们假设一旦冷水接触到热岩石，两者的温度会立即变得相同。这就像拥抱一样；该理论认为你们会瞬间共享完全相同的体温。这让数学计算变得非常简单，但现实中，热量从岩石转移到水中是需要时间的。

新发现： “慢速握手”

这篇论文引入了一种更符合现实的观察问题的方法，称为局部热非平衡（Local Thermal Non-Equilibrium, LTNE）。

他们没有使用“瞬间拥抱”，而是提出了**“慢速握手”**。

现实情况： 当冷水流经热岩石时，水会保持一段时间的低温，而岩石也会保持一段时间的高温。它们不会立即达到温度一致。
旧方法的问题： 如果你假设它们瞬间达到平衡，你可能会误以为流出的地下水降温的速度比实际情况更快。这可能会导致对地热电站寿命或产能的错误判断。

他们是如何解决的： “河流地图”法

作者创建了一个新的计算机模型来追踪这种“慢速握手”。为了实现这一点，他们使用了一个巧妙的技巧，叫做流线模拟（Streamline Simulation）。

想象地下岩石是一个巨大的、混乱的迷宫。

旧方法（网格法）： 想象你在通过观察一个巨大的棋盘来追踪水的路径。你检查每一个方格以确定水的方位。这种方法很精确，但速度极慢，且对计算资源消耗巨大。
新方法（流线法）： 与其使用棋盘，不如想象画出一些看不见的河流（流线），展示水的精确流动方向。计算机只需像火车在轨道上行驶一样，沿着这些“河流”进行追踪。

神奇的技巧：
作者将问题拆分为两个同时发生的环节：

河流之旅： 他们计算冷水如何沿着“河流”（流线）移动。这既快速又简单。
握手过程： 他们计算水在经过岩石时“偷走”了多少热量。这个过程是在“棋盘”（背景网格）上进行的，以考虑到岩石的热量储存。

通过将“移动的水”与“静止的岩石”分离，他们可以运行得更快，同时保持极高的准确性。

他们的发现（结果）

1. “速度测试”（验证）
他们使用已知的数学问题测试了新模型。

结果： 他们的“河流地图”方法极其精确。它几乎完美地匹配了完美的数学解（准确度达 99.9%），甚至比旧的、更慢的方法还要精准。
“平衡”检查： 他们证明了如果热交换速度极快（就像一个非常紧密的拥抱），他们的模型会自动转化为旧的简单模型。这证明了他们的数学逻辑是正确的，因为旧的数学已被包含在其中作为一个特例。

2. “现实世界”测试（非均匀储层）
他们模拟了一个复杂的地下系统，其中包含“快车道”（高渗透通道）和“慢车道”（低渗透区域）。

令人惊讶的发现： 尽管在使用旧的“瞬间拥抱”模型和新的“慢速握手”模型时，流出的地下水看起来大致相同，但系统的内部情况却大不相同。
温度差距： 在新模型中，水和岩石在很长一段时间内都保持着显著的温差。
- 在简单模型中，温差很小（不到 0.5 度）。
- 在新模型中，温差巨大（高达 7.5 度）。
输出结果： 因为水在新的模型中保持较长时间的低温（因为它还没完成足够的“拥抱”），在 400 天后，新模型预测的流出地下水温度比旧模型预测的要热 13 度。

核心结论

论文得出结论：虽然“瞬间拥抱”模型（LTE）在预测井口出水温度的总体趋势方面尚可，但它无法洞察局部正在发生的情况。

类比：
想象一个拥挤的派对。

旧模型假设所有人瞬间对房间的温度达成共识。
新模型意识到，处于角落里的人（岩石）可能仍在流汗，而正冲进门口的人（水）却觉得冷得发抖。即使房间的平均温度看起来正常，角落里的人经历的现实却截然不同。

为什么这很重要：
如果你正在设计一座地热电站，知道岩石仍然比水热得多（即使流出的水看起来温度正常），这说明地下仍有大量的热能尚未被利用。新模型帮助工程师看到了这种“隐藏的能量”。

总结声明

方法： 他们构建了一个快速、准确的计算机模型，能够分别追踪水和岩石的温度。
准确性： 它完美匹配已知数学解，且比旧的基于网格的方法更精确。
关键发现： 在复杂的地下系统中，即使流出的水温变化不大，水和岩石之间仍可能存在巨大的温差（高达 7.5°C）。
局限性： 本文侧重于多孔介质（如沙子或岩石）中的热传输物理学，并未讨论除地热和石油开采建模需求之外的具体临床用途或未来应用。

技术摘要：一种基于流线法的地下多孔介质热传输局部热非平衡框架

问题陈述
地下多孔介质中的热传输对于地热能开采、热采油（EOR）以及含水层热能储存等应用至关重要。传统的油藏规模模拟通常依赖于局部热平衡（LTE）假设，该假设认为代表性体积单元内的流体相和岩石相共享单一温度。虽然 LTE 在计算上非常高效，但它假设相间热交换是迅速完成的。在流体停留时间相对于相间热交换较短的场景中（例如靠近注水井、沿优先流路径、裂缝地层或存在巨大温差时），该假设可能会失效。在这些条件下，局部热非平衡（LTNE）效应——即流体和岩石维持不同的温度——可能会显著改变热前缘的传播过程和生产温度的预测。然而，现有的 LTNE 模型通常局限于理想几何形状的解析解，或者计算量巨大的全域网格数值模型，后者并未显式利用对流主导系统中流向路径的结构特征。

方法论
作者开发了一种基于流线的 LTNE 框架，旨在捕捉独立的流体和岩石温度，同时保持流线模拟的高效性。该方法集成了以下组件：

共同的数值基础： LTE 和 LTNE 模型均使用相同的底层压力场、流线几何结构、映射算子和算子分裂工作流。这确保了两种表述之间可以直接进行“苹果对苹果”的对比。
算子分裂： 传输过程被分解为三个不同的步骤：
- 流线对流： 在正则化的飞行时间（TOF）坐标下，沿规则流线求解占主导地位的流向传输。在 LTNE 模型中，只有移动流体温度（ $T_f$ ）经历对流，并以全间隙速度（ $\beta = 1$ ）运动，而岩石温度（ $T_r$ ）则是静止的。
- 纵向传导： 在正则化的 TOF 网格上求解沿流线方向的热传导。
- 基于网格的更新： 多维传导和局部流-岩热交换在背景欧拉网格上进行处理。
解析交换更新： LTNE 实现的一个关键特征是在每个网格单元的交换子步骤中，使用精确且保守的解析解来更新耦合的流体和岩石温度。这避免了与高交换系数相关的数值刚性问题，并确保了能量的局部守恒。
热传递系数（ $h_a$ ）： 该框架允许将体积流-岩热传递系数指定为均匀参数或空间可变场。对于非均匀情况，通过基于局部渗透率对 $h_a$ 进行参数化，将高渗透区域与较小的基质块及更强的交换联系起来。

核心贡献

首个流线 LTNE 框架： 据作者所知，这是首个用于地下热传输的基于流线的 LTNE 框架，将流线法的效率扩展到了双温度系统。
物理机制分解： 该模型在物理上区分了移动流体对流（在 1D TOF 空间中求解）和静止岩石储热/交换（在背景网格上求解），在保留物理区别的同时保留了计算优势。
精确交换求解器： 实现了一种解析的、非迭代的交换更新，确保了即使在高交换系数下，局部交换子问题的无条件稳定性。

结果
该框架经过了解析基准测试验证，并应用于非均匀油藏模型（SPE10 第 36 层）：

1D 验证：
- LTE： 流线解与 Ogata–Banks 解析解匹配， $R^2$ 为 0.9928，优于同一网格上的有限差分法（FDM）（ $R^2 = 0.9669$ ）。
- LTNE： 流体和岩石温度解与 Schumann 解析解匹配， $R^2$ 为 0.9999，显著优于 FDM（ $R^2 \approx 0.9975$ ）。
- 收敛性： 随着体积交换系数（ $h_a$ ）从 10 增加到 100 W m⁻³ °C⁻¹，LTNE 热容量加权平均温度向 LTE 解收敛，其 $R^2$ 从 0.9442 增加到 0.9996。
非均匀油藏案例（500 天模拟）：
- 热突破： 在四分之一五点法配置下，突破发生在约 200 天时。
- 产出温度： 在 400 天左右，LTNE 模型预测的产出流体温度比 LTE 模型高出约 13 °C。这种差异源于 LTE 强制瞬时平衡，而 LTNE 则允许岩石保留热量并缓慢将其传递给流体。
- 空间变异性： 与均匀 $h_a$ 情况相比，引入基于渗透率的空间可变 $h_a$ 对宏观产出温度曲线的影响微乎其微。然而，它显著影响了局部的相间失衡。在 500 天时，最大流-岩温差从 0.43 °C（均匀情况）增加到 7.57 °C（可变情况），增幅达 18 倍。

意义与主张
本文主张，尽管宏观热响应（如产出温度曲线）可能表现出接近平衡的状态，但在非均匀系统中，特别是具有空间可变交换特性的系统中，显著的局部热非平衡现象仍然存在。研究表明，均匀的交换系数虽然可以准确重现宏观生产数据，但无法代表局部温度差的量级和空间分布。所提出的框架提供了一个计算高效的工具，用于判定何时 LTE 是足够的，以及何时必须保留双温度效应以实现准确的地下热传输模拟。