🔬 physics

The Statistical physics of unsaturated soil water: kinetic theory and non commutative pore water dynamics

本文提出了一个多尺度统计力学框架，将非饱和土壤水的动力学统一为一个单一的动力学连续体方程，证明了诸如 Richards 方程等经典模型以及优先流和滞后现象，都是由控制非平衡孔隙占据率的孔隙解析达姆科勒数（Damkohler number）所导致的特定极限或几何结果。

原作者： Riccardo Rigon

发布于 2026-07-13

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Riccardo Rigon

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

不要把土壤仅仅想象成一团泥土，而要把它想象成一个巨大的、三维的迷宫，由数以十亿计、粗细各异的微小中空管状结构组成。有些管子像人的头发一样细；有些则像饮用吸管一样宽。在这个迷宫中，水正试图寻找出路。

一个多世纪以来，科学家们一直试图用一个著名的规则——**理查兹方程（Richards equation）*来描述水的运动方式。你可以把这个方程看作一本“慢动作”规则书。它假设如果你把水倒在土壤上，每一滴水都会瞬间找到最舒适、最节省能量的位置。它假设水总是处于完美的、平衡的状态，并且唯一重要的因素仅仅是土壤中水的总量*。

重大发现：土壤拥有记忆（以及一个速度限制）

里卡多·里贡（Riccardo Rigon）的这篇论文认为，旧的规则书是不完整的。这就像仅通过观察人群的平均密度来描述一场混乱的“开场舞”（mosh pit）一样——你忽略了人们正在奔跑、互相碰撞以及被困在特定位置的事实。

该论文引入了一个更详细的新视角。它不再仅仅询问“土壤有多湿？”，而是询问：“哪些特定的管子是满的，哪些是空的？”

作者提出了一种新的“动力学”理论。他们将土壤中的水视为一个位于拥有许多房间的建筑中的人群。

微观尺度（走廊）： 水根据物理定律（哈根-泊肃叶定律/Hagen–Poiseuille rates）在管子之间移动。这就像水流经管道一样。
介观尺度（房间）： 当你观察一整个房间（一小块土壤）时，水正在不断地交换位置，试图稳定到最舒适的排列状态。
宏观尺度（建筑）： 当你缩放到整个田野时，你会看到水从一个房间流向另一个房间。

“达姆科勒数”（Damköhler Number）：土壤的速度计

这篇论文中最令人兴奋的部分是一个被称为**达姆科勒数（Da）**的神奇数字。你可以把它看作土壤的“速度计”。它比较了两个速度：

水想要重新排列自己以寻找完美位置的速度（弛豫速度）。
雨水击中地面的速度（强迫速度）。

如果雨下得很慢（Da 很小）： 水有足够的时间重新排列自己。它能找到完美的位点，就像旧的理查兹方程所预测的那样。此时土壤表现得像一种平静、缓慢流动的流体。
如果是一场突如其来的暴雨（Da 很大）： 水被灌入的速度超过了它重新排列的速度。它没有时间找到那些微小且舒适的管子。相反，它会卡在大而快的管子里，产生“优先流”（preferential flow）（就像水通过消防水带喷射而出一样）。

“优先流”之谜解开了吗？

长期以来，科学家们认为“理查兹流”（缓慢、稳定）和“优先流”（快速、混乱）是两种完全不同的现象，需要两套不同的规则书。

这篇论文指出：不。它们是同一个方程，只是根据降雨速度的不同而表现出不同的行为。

理查兹流是当土壤有时间喘息时发生的情况。
优先流是当土壤被压倒时发生的情况。

并不是土壤改变了主意，而是事件的速度改变了结果。论文表明，如果你知道降雨的速度和管子的尺寸，你就可以准确预测土壤何时会从“平静”转向“混乱”。

“滞后”环：几何幽灵

你可能已经注意到，土壤在变干的过程中和在变湿的过程中保持水分的方式是不同的。这被称为滞后现象（hysteresis）。

旧观点： 科学家认为这是因为土壤具有“双稳态”管子——就像取决于历史记录而卡在“开”或“关闭”位置的灯开关。
新观点： 论文认为这并非关于卡住的开关，而是关于几何形状与时间。

想象你在森林里行走。如果你走得很慢，你可以跨过每一根树枝。如果你奔跑，你可能会被树枝绊倒并走上一条不同的路径。如果你绕圈行走（先湿后干），你可能不会回到最初的起点，因为“奔跑”的路径与“行走”的路径不同。

作者称之为几何相位（geometric phase，或称全纯性/holonomy）。这就像一个“几何幽灵”。土壤记住的是它达到当前湿度所经过的路径，而不只是湿度本身。论文指出，这种“记忆”的大小（滞后环）随降雨强度的平方增长。如果你将降雨速度增加一倍，记忆效应会扩大四倍。这是一个具体的、可测试的预测，而旧的“卡住开关”模型无法解释这一点。

什么是被证明的，什么是被建议的，以及什么是被排除的

被排除的： 论文明确反对认为“优先流”和“理查兹流”是两种完全无关的现象。它们是一个连续的交叉过程。它也排除了滞后现象是由单个孔隙进入“二元状态”导致的观点；相反，它是整个网络几何形状和强迫速度的结果。
被证明的（在其数学推导中）： 作者在数学上证明了他们的新方程符合物理定律（如质量守恒），并且系统始终向着更低能量状态移动（H-定理）。他们展示了如果将一切都放慢速度，新方程在数学上可以转化为旧的理查兹方程。
被建议/模拟的： 滞后环面积随降雨强度平方增长（ $H \propto I^2$ ）这一具体预测是基于其理论的一个猜想。他们建议这可以在实验室中进行测试。他们还建议，“田间持水量”（土壤停止排水的点）实际上是一个“渗透阈值”（水网络破碎的点），这是一个合理的想法，得到了其模型的支持，但仍需更多的现实验证。
仅通过模拟展示： 关于“非对易性”（即先湿后干与先干后湿得到不同结果这一事实）是滞后现象根源的观点，是通过他们的数学推导以及对小型管状网络进行的计算机模拟来展示的。

底线

这篇论文不仅提供了一个新公式，它还提供了一种看待世界的新方式。它表明土壤不是一个静态的海绵，而是一个动态的、具有时间敏感性的迷宫。水是缓慢均匀地流动，还是在混乱的流中冲刷而过，完全取决于降雨的速度与土壤内部重新排列的速度之间的关系。

作者们确信，他们的新方程是开启“缓慢稳定流”与“快速混乱流”的“万能钥匙”，将两者统一成一个单一、优美且具有严密数学逻辑的故事。他们邀请科学家们去测试他们的预测：如果你测量不同降雨速度下的滞后环，它应该呈二次方增长。如果确实如此，那么“几何幽灵”就是真实的。如果不是，那么理论需要重写。

技术摘要：非饱和土壤水的统计物理学

问题陈述
非饱和多孔介质的动力学本质上是不可逆的，其特征是依赖于强迫作用（润湿与干燥）历史的非平衡构型。经典的模型（主要是 Richards 方程）假设存在局部热力学平衡，并使用单一标量状态变量（含水量 $\theta$ 或基质势 $\psi$ ）。这一框架无法解释以下现象：

滞后现象（Hysteresis）： 相同的宏观状态对应于不同的内部构型，这取决于润湿或干燥的路径。
优势流（Preferential Flow）： 在剧烈强迫作用下，存在缓慢的毛细管控制再分配与通过大孔隙的快速旁路流并存的现象，这违反了 Richards 方程的单域假设。
速率依赖性（Rate Dependence）： 速率无关模型无法捕捉滞后现象和流动机制对强迫速度的依赖性。

现有的双域或双渗透率模型通过引入额外的状态变量和交换项来处理这些问题，但这些模型属于现象学描述，并未从统一的物理描述中推导而来。

方法论
本文通过三个尺度开发了一种统计力学理论，将孔隙占据分数 $g(r, x, t)$ （在位置 $x$ 、时间 $t$ 时，半径为 $r$ 的孔隙被水填充的比例）视为基本状态变量。

微观尺度（动力学速率）：
- 孔隙间的液体转移由 Hagen–Poiseuille 速率 $\kappa(r)$ 和驱动势 $\Phi(r, r')$ 控制，后者定义为不同孔径类别的化学势之差（结合了毛细作用和重力项）。
- 推导出了状态相关的速率常数 $\kappa(r, g)$ ，该常数随着填充分数 $g$ 的增加，实现了从 Lucas–Washburn 动力学（快速初始入侵）到 Hagen–Poiseuille 流（稳态）的过渡。
中观尺度（主方程）：
- 对代表性单元体积（REV）进行平均，得到关于 $g$ 的主方程。
- 该方程包含代表孔隙类别间水交换的增益（ $G$ ）项和损失（ $L$ ）项。这些项利用了“分子混沌假设”（Stosszahlansatz），假设供体和受体占据状态在统计上是独立的，从而导致双线性形式 $g(r')[1-g(r)]$ 。
- 系统被证明是一个 Onsager 梯度流，在具有状态相关迁移率的 Gibbs 自由能泛函 $F[g]$ 下进行下降。这确保了在等温、无强迫系统中的 H-定理（ $dF/dt \leq 0$ ）。
连续介质尺度（场方程）：
- 将 REV 收缩为一个点，使 $g$ 提升为连续场。
- 生成的根本方程是一个平衡律：
  $\partial_t g + \nabla \cdot \mathbf{F} = C[g] - E - T$
  其中 $C[g]$ 是 REV 内部的再分配算子， $\nabla \cdot \mathbf{F}$ 是输运散度（REV 间通量）， $E$ 和 $T$ 是汇项（蒸发、根系吸收）。
- 通量 $\mathbf{F}$ 在孔径半径上是非局域的，通过连通性矩阵 $C(r, r')$ 将不同的孔径类别耦合在一起。

核心贡献与结果

统一流动机制： 本理论引入了一个单一的无量纲控制参数——孔隙解析达姆科勒数（Damköhler number） $Da(r, x) = \tau_{redis}(r) / \tau_{forcing}(x)$ ，它比较了内部弛豫时间与强迫时间尺度。
- $Da \ll 1$ （准静态）： 系统追踪局部平衡。通过对动力学方程进行 Chapman–Enskog 约化，可以恢复 Richards 方程，其中基质势 $\psi$ 和水力传导度 $K(\theta)$ 作为极限出现。
- $Da \gtrsim 1$ （动力学）： 再分配无法跟上强迫作用的速度。水占据的是几何上可达的路径（大孔隙），而非热力学最优路径，从而自然地重现了优势流，无需引入人为的 dual-domain 假设。
滞后的几何起源： 滞后现象并非归因于单个孔隙的双稳态（如墨水瓶效应），而是归因于“强迫束”（forcing bundle）的全纯性（holonomy）。
- 润湿算子（ $L_W$ ）和干燥算子（ $L_D$ ）由于非选择性填充（降水）与受拓扑结构约束的选择性排空（排水/蒸发）之间的不对称性，并不对易（ $[L_W, L_D] \neq 0$ ）。
- 这种非对易性在封闭的强迫循环后产生了一个几何相位（残余状态）。
- 预测： 对于慢速循环，滞后环面积 $H$ 与强迫强度呈二次方关系（ $H \propto I^2$ ），这是速率无关模型所不具备的特征。
源自微观结构的本构关系：
- 水力传导度 $K(\theta)$ 被推导为一种输运系数（迁移率核的谱极限），而非拟合参数。
- 理论恢复了经典极限：毛细管束模型（对角极限）、临界路径模型（渗流附近的谱极限）以及动态毛细压力模型（一阶矩近似）。
- 田间持水量被识别为水网络破碎时的渗透阈值（ $\theta_c$ ），此时宏观传导度消失。
不可约状态变量： 在非准静态极限下，标量含水量 $\theta$ 不足以表征系统。完整的分布 $g(r, x, t)$ 才是不可约的状态变量。

意义与主张
本文声称提供了一个建设性的、基于热力学的框架，将此前截然不同的现象学（Richards 流、优势流、滞后现象）统一到一个单一的动力学方程中。

理论严密性： 该理论被呈现为多孔介质版的 Boltzmann 方程，其中 $g(r)$ 对应于速度分布 $f(v)$ ，其 Chapman–Enskog 约化过程镜像了 Navier–Stokes 方程的推导过程。
预测能力： 所有输入（孔径分布 $f(r)$ 、连通性 $C(r, r')$ 、扭曲度 $\bar{\tau}$ ）均为可通过微型 CT 图像测量的几何属性，无需针对宏观持水曲线或传导度数据进行校准。
可证伪性： 该理论提出了具体的、可测试的预测，最显著的是滞后环面积随强迫强度变化的 $H \propto I^2$ 标度律，以及基于土壤结构和流体性质的优势流起始阈值强度的存在。
适用范围： 该框架具有普适性，可应用于其他具有受限且依赖历史的转变过程的系统（如塑性变形、磁滞现象），并且可以通过优化驱动势 $\Phi$ 来扩展到吸附、渗透或热力学项。

文章得出结论：平衡流与非平衡流之间的分歧不在于选择不同的模型，而在于认识到在单一统一的动力学理论中，由达姆科勒数控制的连续交叉过程。

技术摘要：非饱和土壤水的统计物理学

类似论文