🔬 物理学

Dimensional effects on breakthrough behavior and pore-scale residence time distributions

本研究利用格子玻尔兹曼模拟证明，虽然二维孔隙尺度模型由于面外连通性的减少，会系统性地低估三维渗透率——这一差异通过欧拉示性数能得到最佳预测——但维度对突破曲线和停留时间等传输指标的影响，会根据介质的具体孔隙几何结构和孔隙率而产生显著差异。

原作者： Shahram Asgari, Christof Meile

发布于 2026-10-16

📖 1 分钟阅读☕ 轻松阅读

原作者： Shahram Asgari, Christof Meile

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

水与溶解物质在土壤、岩石和工程材料那隐秘而复杂的微观世界中移动，其方式往往是肉眼无法观测到的。这种被称为溶质运移的过程，对于从清理受污染的地下水到深层地下碳封存等诸多领域都至关重要。为了理解这些流体如何在沙粒或过滤器纤维之间的微小蜿蜒通道中穿行，科学家们构建了模拟单个孔隙尺度下流动的计算机模型。几十年来，研究人员一直面临一个实际的困境：虽然三维模型能提供最真实的复杂网络图景，但它们对计算能力的要求极高。二维切片——本质上是材料的平面横截面——运行起来要廉价且快速得多，因此许多人将其作为三维现实的替代品。长期以来的疑问在于，这种捷径是否以牺牲过多的准确性为代价，从而使其变得不再实用。

佐治亚大学的一个研究小组致力于测试这些平面模型究竟遗漏了多少信息。他们创建了三种不同类型多孔材料的数字版本：一种由缠绕的纤维组成，一种类似于堆叠的球体（如一堆大理石），第三种则模仿自然地质结构的随机、块状形态。针对每种材料，他们都构建了一个完整的三维块体以及一系列从其中切割出的二维切片。利用一种能够追踪流体颗粒和溶解物质运动的高效模拟方法，他们对比了水和示踪剂在平面切片与完整块体中的移动情况。他们特别观察了流体通过材料的难易程度、物质从一侧到达另一侧所需的时间，以及物质在途中扩散的程度。

结果证实，平面模型始终低估了流体通过材料的容易程度。在三维块体中，流体可以向各个方向绕过障碍物，寻找那些在平面切片中根本不存在的捷径。当研究人员观察二维切片时，发现由于流体无法通过“向上”或“向下”脱离平面的方式绕过障碍物，被迫走上了更长、更曲折的路径。这种效应在孔隙率较低的材料中最为显著，因为这类材料可供流动的空间本就狭窄。在这些情况下，对于纤维状材料，平面模型预测的渗透率值几乎比真实的的三维现实低了十倍；而对于类地质结构，则低了约四倍。研究人员发现，这种差异不仅在于存在多少空隙空间，更在于一种根本性的连接缺失。通过将世界扁平化，这些模型切断了让流体高效编织穿梭于网络中的关键的“面外”连接。

为了预测何时可以使用平面模型作为安全的替代方案，该团队转向了一种通过计算孔隙网络内连接数量的数学方法。他们使用了一种衡量指标，该指标本质上是在统计结构中的独立口袋、环路和封闭空间的数量。他们发现，当平面切片与完整块体的连接计数非常不同时，流体预测就会出现严重偏差。这种拓扑学度量可以作为一个快速的诊断工具：如果平面切片捕捉到了与完整块体相同的环路和连接模式，那么它很可能提供准确的流动图景。如果计数差异显著，则说明平面模型丢失了关键路径，此时必须进行完整的三维模拟。

除了研究水的流动速度外，该研究还考察了物质通过材料所需的时间以及它的扩散程度。这些运移指标的行为在很大程度上取决于材料的类型。在纤维介质中，无论材料的渗透率如何，平面模型与完整模型之间的差异都保持稳定，这可能是因为流动的控制因素是少数几条主通道，而这些通道在两个维度上的形态相似。在类地质结构中，随着材料渗透率的增加，两种模型之间的差距逐渐缩小，这表明在更松散、更开放的材料中，平面切片能够较好地捕捉流动路径。然而，在堆叠球体模型中，平面切片有时预测的旅行时间更快，有时则更慢，没有表现出清晰的规律，这表明单个切片中球体的特定排列可能会产生误导。

研究人员还测试了被追踪物质的化学性质是否会改变这些结果。他们对比了一种仅随水移动而不发生反应的保守示踪剂，以及一种会随移动与材料发生相互作用的反应性示踪剂。令人惊讶的是，示踪剂的类型对整体对比结果几乎没有影响。无论物质是惰性的还是反应性的，使用平面模型引入的误差主要受控于孔隙的形状和连通性，而非流体的化学性质。这表明，二维模型的结构局限性是一个基本的几何问题，而非化学问题。

最终，这项研究为必须在速度与精度之间做出抉择的科学家和工程师提供了明确的指南。虽然二维模型对于快速估算是有用的，但由于它们无法考虑到三维空间中可用的额外路径，因此会系统性地忽略流动的便捷性。研究表明，通过检查孔隙是如何连接的这一简单指标，人们可以在运行模拟之前判断一个平面切片是一个忠实的代理，还是一个具有误导性的简化模型。这一洞察力使研究人员能够更高效地利用计算资源，确保在关键时刻，自然界复杂的、三维性的本质得到充分尊重。

技术摘要：维度效应对突破行为及孔隙尺度停留时间分布的影响

问题陈述
在研究多孔介质中的溶质运移时，二维（2D）孔隙尺度模型常被用作三维（3D）模拟的计算高效代理。然而，由于连接性、曲折度和混合机制的差异，2D 和 3D 多孔结构表现出截然不同的流动和运移行为。虽然 3D 模拟能提供最真实的孔隙网络表征，但其计算强度极高。因此，研究人员经常依赖从微米级 CT 图像或薄切片中提取的 2D 切片。本研究旨在解决的核心问题是：这些 2D 类比模型在预测 3D 渗透率、突破曲线和流体年龄分布方面的准确程度如何，以及是否存在特定的拓扑或几何指标，可以决定何时 2D 模型可以作为全 3D 模拟的合适替代方案。

方法论
作者利用格子玻尔兹曼（Lattice Boltzmann, LB）方法，在离散体素网格上求解纳维-斯托克斯方程（Navier–Stokes）和对流-扩散-反应（ADR）方程。研究重点针对三种不同类型的合成多孔介质：

纤维状（Fibrous）： 代表生物组织和工程材料，通过放置圆柱形纤维生成。
高斯型（Gaussian）： 类似于土壤聚集体和砂岩等地质材料，通过滤波随机场生成。
重叠球体（Overlapping Spheres）： 代表颗粒材料和填充床，通过随机放置球体构建。

对于每种介质类型，生成了三种不同孔隙率（0.4、0.6 和 0.8）的 3D 计算域（100 × 100 × 100 体素）。2D 域作为 xz 平面切片从这些 3D 结构中提取，并确保流体贯通。

研究对比了以下指标：

渗透率： 分别通过达西定律计算 2D 切片及其 3D 母体域的渗透率。
运移指标： 分析了通量加权突破曲线，以导出第一阶时间矩（ $M_1$ ，平均到达时间）和变异系数（ $CV$ ，突破扩散度）。
示踪剂类型： 在对流主导条件下（Péclet 数分别为 4.1、12.5 和 37.5），模拟了保守型（非反应型）和反应型（代表具有一级反应速率的流体年龄）示踪剂。
拓扑分析： 计算了基于持续同调（Betti 数 $\beta_0, \beta_1, \beta_2$ ）得到的孔隙归一化欧拉特征数（ $\tilde{\chi}$ ），将其作为量化孔隙网络连接性的纯几何指标。

关键结果

渗透率的系统性低估：
2D 模拟在所有介质类型和孔隙率下均系统性地低估了真实的 3D 渗透率。这种偏差在低孔隙率时最为显著。例如，在 $\phi=0.4$ 的纤维状介质中，2D 渗透率几乎比 3D 值低 10 倍。这种偏差归因于 2D 切片中平面外连接性的缺失，这增加了有效曲折度并减少了贡献于活性流动的孔隙空间比例。Kozeny–Carman 关系证实，仅靠孔隙率差异无法解释这种差异；拓扑上的断裂才是主要驱动因素。
欧拉特征数（ $\tilde{\chi}$ ）的作用：
作为一种衡量连接性的拓扑度量，孔隙归一化欧拉特征数在很大程度上追踪了 2D 与 3D 渗透率之间的差异。负值的 $\tilde{\chi}$ （表示由环路主导的良好连接网络）与较高的渗透率和较小的 2D–3D 差异相关；而正值的 $\tilde{\chi}$ （破碎的网络）则与较低的渗透率和较大的差异相关。作者提出将 $\tilde{\chi}$ 作为一种低成本、仅基于几何形状的代理指标，用于在运行流动模拟前对 2D 切片进行筛选。
特定介质的运移趋势：
维度效应对运移指标（ $M_1$ 和 $CV$ ）的影响因介质类型而异：

纤维状介质： 平均旅行时间和突破扩散度的相对差异在整个渗透率范围内几乎保持不变。这是由于流动组织成主导通道，而这些通道在 2D 中得以保留。
高斯型介质： 2D–3D 差异随着渗透率的增加而减小。在低渗透率下，狭窄的孔喉限制了 2D 流动，但在 3D 中可以通过平面外路径绕行；随着渗透率上升，这些约束得到缓解，从而缩小了维度差距。
重叠球体介质： 平均到达时间的差异较大且具有双向性（部分 2D 切片高估，部分低估了 3D 结果），且与渗透率之间没有明确趋势。这反映了大型互连孔隙与单个 2D 切片特定孔隙率之间复杂的相互作用。

示踪剂反应性的影响：
保守型和反应型示踪剂表现出相似的维度偏差。虽然反应型示踪剂在平均旅行时间（ $M_1$ ）上的相对差异比保守型示踪剂系统性地大 20–40%，但整体模式是由孔隙几何结构而非示踪剂反应性决定的。突破扩散度（ $CV$ ）表现出更多变异性，其中高斯型介质显示出明显的偏差，即保守型示踪剂比反应型示踪剂表现出更大的 2D–3D 差异。

意义与主张
本文得出结论：孔隙几何结构和连接性是孔隙尺度建模中维度偏差的主要控制因素，而示踪剂反应性的作用次之。研究表明，尽管 2D 模型计算效率高，但它们会在渗透率和运移预测中引入系统误差，且这些误差无法仅通过匹配孔隙率来修正。

作者主张，孔隙归一化欧拉特征数（ $\tilde{\chi}$ ）提供了一种实用的预模拟筛选工具。通过将候选 2D 切片的 $\tilde{\chi}$ 与 3D 母体进行比较，研究人员可以识别哪些 2D 切片保留了 3D 系统的拓扑连接性。如果 $\tilde{\chi}$ 的值在符号和量级上一致，则该 2D 切片可以作为有效的代理；否则，则需要进行全 3D 模拟。这一方法表明，在构建代表 3D 系统的 2D 微流控器件（如微流控设备）时，通过拓扑描述符匹配孔隙尺度的连接性，与匹配孔隙率或表面积同样至关重要。

类似论文