🔬 凝聚态物理

A Conservative Darcy-Brinkman-Stokes Formulation for PHREEQC-Coupled Reactive Transport with Cell-Based Surface Complexation

本文提出了一种在 GeoChemFoam 中实现的质量守恒 Darcy-Brinkman-Stokes 公式化方法，并将其与 PHREEQC 相耦合，从而能够在解析及非解析的多孔介质结构中，实现具有基于单元的表面络合作用的精确反应输运模拟。

原作者： Hannah P. Menke, Florian Doster, Julien Maes

发布于 2026-10-15

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Hannah P. Menke, Florian Doster, Julien Maes

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

在深层地下，水的流动并非像汽车在高速公路上那样在空旷的隧道中穿行。相反，它是在由孔隙组成的岩石迷宫中移动，这意味着岩石充满了微小的、相互连通的孔洞和裂缝。在许多这类地下环境中，岩石不仅仅是一个被动的容器；其表面具有化学活性。随着水的移动，溶解的矿物质和污染物会粘附在这些表面上，这一过程减缓了它们的迁移速度并改变了它们的浓度。这种相互作用对于理解污染物如何在地下水层中移动、我们如何从盐水中提取锂等有价值的矿物质，以及如何将放射性废物安全地储存数千年，都至关重要。科学家面临的挑战在于，岩石的内部结构极其复杂。有些部分拥有水流湍急的大型开放通道，而另一些部分则是由具有微观孔隙的细密、海绵状材料组成的，这些孔隙蕴含着巨大的内部表面积。传统的计算机模型往往难以同时捕捉到大型通道中的快速流动以及微小、不可见的孔隙内缓慢且具有粘性的化学反应。

一个研究团队开发了一种模拟这些复杂相互作用的新方法，使他们能够在单一的统一网格上同时模拟水的流动和颗粒的化学粘附。他们的方法将岩石视为一种连续的材料，而不是将其视为一系列分离的大孔洞和微型海绵；在他们的模型中，水的含量和可用于化学反应的表面积可以在模型的每一个微小单元之间发生变化。在这种新方法中，计算机将岩石图像划分为许多小单元（cells）。对于每个单元，模型都会计算实际流经其中的水量，以及可供化学物质粘附的表面积，即使该表面积在原始图像中过于微小而无法被看见。通过追踪化学物质的总质量——包括溶解在水中的部分和粘附在岩石上的部分——研究人员确保在计算过程中不会因失误导致物质的丢失或凭空产生。这种质量守恒对于准确预测污染物需要多久才能迁移，或者可以回收多少有价值的矿物质，都至关重要。这是因为质量守恒对于做出关于污染物迁移时间或可回收矿物量的准确预测至关重要。

研究人员使用一系列严格的基准测试来检验他们的新方法是否有效。他们首先从简单的、可以通过手工计算答案的一维场景开始，确认他们的计算机模型在化学物质移动速度以及由于粘附导致的减速效应方面，与预期结果相匹配。随后，他们转向了使用真实岩石和合成材料的三维图像。在其中一项测试中，他们模拟了一个具有光滑球形形状的岩石，以观察其方法能否准确测量表面积。他们发现，即使在图像分辨率相对较低的情况下，他们的方法也能以小于百分之一的误差计算出表面积。在另一项测试中，他们将这种新方法与一种更详细、计算成本更高的模拟方法进行了对比，后者明确绘制了每一个微小的孔隙。他们的新方法速度更快，且能够重现详细模拟的结果，在捕获锂的总量方面差异小于百分之一，这证明了该方法无需解析每一个微观细节，也能处理现实世界岩石的复杂性。

随后，该团队将他们的方法应用于三个不同的现实场景，以观察其实际表现。首先，他们模拟了一个裂隙岩石系统，其中铅离子在涂有氧化铁的裂隙中移动。模拟展示了铅是如何由于粘附在涂层上而减速的，从而清晰地描绘出污染物如何在裂隙网络中扩散。其次，他们模拟了从含有大量镁的咸盐水中捕获锰的过程。由于锰和镁会竞争岩石表面的相同粘附位点，模型必须计算出谁能胜出。结果显示了锰是如何逐渐置换掉镁的，其在岩石表面积累的模式取决于水的流动。最后，他们模拟了锶在石灰岩层中的运动，这一场景与深层地质处置库的安全性相关。在这种情况下，锶必须与水中已存在的钙进行竞争。模拟成功追踪了锶在向前移动时如何将钙从岩石表面挤压下来，形成了一个明显的化学交换波。

这项工作的核心发现是：无需解析每一个微观孔隙，即可对复杂的多尺度岩层进行反应输运建模。通过将反应表面积视为单元本身的一种属性，而非仅仅是可见的岩石边界，研究人员可以将那些对化学反应至关重要的、不可见的内部表面纳入考虑。这使得模拟污染物和有价值矿物质在地下环境中的移动变得更加准确，即使岩石包含大型开放空间和微小海绵状区域的混合体。该方法具有保守性，这意味着它严格遵循质量平衡定律，确保整个模拟过程中物质总量保持不变。这种可靠性使其成为预测地下系统长期行为的强大工具，无论是清理受污染的含水层，还是设计更安全的放射性废物储存库。研究人员证明，他们的方法适用于广泛的流速和岩石类型，为理解深层地球隐藏的化学作用提供了一个稳健的框架。

技术摘要：一种用于 PHREEQC 耦合反应输运的保守型 Darcy-Brinkman-Stokes 公式化方法

问题陈述
多尺度多孔介质（如碳酸盐岩）中的反应输运提出了显著的模型挑战。这些系统通常包含开放的大孔隙以及微孔基质区域。虽然基于图像的模拟可以解析大孔隙，但微孔相的内部反应表面积往往低于图像分辨率（亚体素尺度）。传统的 Darcy 尺度模型依赖于本构上尺度假设，无法捕捉空间异质性；而完全解析的孔隙尺度模型（例如显式边界 Darcy-Brinkman-Stokes 或 DBS）在主导反应面积位于粗糙计算单元内部而非网格边界时会失效。现有的 PHREEQC 耦合工作流通常将化学作为通量条件施加在解析的界面上，这导致未解析的内部表面积被忽略，或者被表示为破坏质量守恒和单元尺度库存物理定义的权宜性的体积汇项。

方法论
作者提出了一种在 GeoChemFoam（一个基于 OpenFOAM 的工具箱）中实现的质量守恒 DBS 公式化方法，并通过 PhreeqcRM 接口与 PHREEQC 相耦合。其核心创新在于一种基于单元的表面络合方法，其中反应表面是单元体积的属性，而非网格边界的属性。

控制方程： 流动受 Darcy-Brinkman-Stokes 方程支配，其中单元中心的孔隙率场 ( $\epsilon$ ) 在固体与开放孔隙之间平滑变化。对物理移动水组分 ( $\theta$ ) 进行输运求解，该组分控制着水相存储与扩散。
表面积定义： 反应表面积 ( $a_{surf}$ $a_{s u r f}$ ) 被定义为一个单元中心的体密度 ( $m^2 m^{-3}_{bulk}$ $m^{2} m_{b u l k}^{- 3}$ )，由两个相加项组成：
1. 界面项 ( $a_{interface}^{surf}$ )：通过平滑界面算法（如 occupancyGaussian）从解析的图像界面（孔隙相与固体相之间）重建。
2. 内禀项 ( $a_{intrinsic}^{struct}$ )：一种给定的亚体素贡献，由基于相位特定表面密度 ( $s_k$ ) 并结合几何相位组分 ( $f_k$ ) 加权得出。这使得能够表示图像分割无法解析的内部微孔表面积。
保守耦合： 局部质量库存定义为 $M = \theta C + a_{surf} \Gamma$ $M = θ C + a_{s u r f} Γ$ ，其中 $C$ $C$ 是水相浓度， $\Gamma$ $Γ$ 是表面浓度。求解器使用一阶 Lie 分裂方案：
1. 输运步骤：使用有限体积 DBS 输运法推进水相组分的质量守恒过程。
2. 化学步骤：将更新后的水相状态、移动水组分及基于单元的表面库存映射到 PHREEQC。PHREEQC 计算局部热力学平衡（表面络合），并返回更新后的水相和吸附态组分，确保每个单元内水相与表面相之间的严格质量守恒。
实现： 反应表面积和位点密度在模拟过程中被视为固定场。移动水组分低于阈值 ( $\theta_{mask}$ ) 的单元不参与化学计算。

主要贡献

解耦流动与化学体积：该公式化方法区分了数值 DBS 孔隙率（控制流动阻力）与物理移动水组分（控制化学溶液体积），尽管在所呈现的饱和案例中两者被设为相等。
基于单元的未解析表面：该方法成功地将未解析的内部表面积集成作为一种内禀体积密度，从而实现了在无需显式解析内部孔隙网络的情况下，对微孔基质区域进行反应输运建模。
保守质量平衡：该方法确保了溶解态和表面结合态组分进入相同的局部质量平衡，避免了源项的启发式重新分配。

结果与验证
该公式化方法通过了一套全面的测试套件进行了验证：

1D 滞后基准测试：耦合的平流-扩散-吸附响应与解析解及有限域参考值相匹配，证明了正确的滞后因子和空间与时间的一阶收敛性。
封闭系统批次平衡：质量平衡测试确认，PHREEQC 内部水相与表面相之间的质量转移保持了总库存，残差低于 $10^{-9}$ 。
表面积重建：使用平滑界面算法针对解析几何体（球体和圆柱柱体）进行了测试。与体拟合的 snappyHexMesh 参考值相比，误差被限制在亚百分比水平（例如，在足够分辨率下，2D 柱体的误差 $<0.2\%$ ，3D 球体的误差 $<1\%$ ）。
解析-DBS 两尺度基准测试：通过锂微流控芯片模型，将完全解析的孔隙尺度模拟与上尺度 DBS 公式化方法进行对比，结果显示在 Péclet 数 ( $Pe$ ) 为 0.1 至 10 的范围内，该方法对水相和吸附态库存的重建准确度均在 0.6% 以内。虽然总摄取量非常准确，但吸附态流向剖面存在高达 4.75% 的形状差异，这归因于孔隙尺度浓度变化的平滑处理。
应用案例：该模型被应用于三个复杂的基于图像的场景：
- 2D 断裂网络：氢氧化铁 (HFO) 涂层对 Pb $^{2+}$ 的滞后作用，捕捉了强位点占据的空间异质性。
- 3D Kelvin Cell：竞争性的 Mn $^{2+}$ /Mg $^{2+}$ 交换，展示了模型处理微孔晶格中多组分竞争的能力。
- 3D Estaillades 石灰岩：模拟放射性核素迁移的 Sr $^{2+}$ /Ca $^{2+}$ 交换，展示了 Sr $^{2+}$ 在水相和吸附相中置换 Ca $^{2+}$ 的层析交换模式。

意义与声明
本文声称该框架支持在解析及未解析孔隙结构中的保守表面络合建模。通过将反应表面视为单元属性而非网格边界，该方法克服了传统孔隙尺度模型在多尺度介质中的分辨率限制。这对于由选择性和表面可及性控制性能的应用（如污染物拦截、地质处置库中的放射性核素输运以及关键矿物回收，如直接锂提取）具有重要意义。作者指出，内禀面积表示的准确性依赖于外部表征（如 BET 吸附）来约束亚体速表面密度，且目前的实现假设为局部平衡和固定表面积，不包括动力学反应及矿物溶解/沉淀。

技术摘要：一种用于 PHREEQC 耦合反应输运的保守型 Darcy-Brinkman-Stokes 公式化方法

类似论文