本研究报告旨在向中文研究界介绍由Lichun Wang、Seonkyoo Yoon、Lizhi Zheng、Tiejun Wang、Xi Chen及Peter K. Kang合作完成，于2023年11月10日在线发表在《Journal of Hydrology》第627卷上的一项原创性研究成果。该研究题为“Flux exchange between fracture and matrix dictates late-time tracer tailing”，聚焦于裂隙岩体中溶质运移的核心挑战，即如何理解和量化溶质运移后期出现的幂律拖尾现象，并首次在三维空间中系统探究了裂隙与基质之间通量交换对拖尾行为的决定性影响。
研究的学术背景 该研究隶属于水文地质学与多孔介质传质领域。在地下水污染修复、核废物地质处置、二氧化碳地质封存等关键工程应用中，准确预测溶质在裂隙岩体中的运移行为至关重要。然而，裂隙介质的显著非均质性使得运移预测极具挑战。大量现场、实验与数值研究表明，裂隙岩体中的溶质运移常偏离经典的Fickian扩散定律，表现出非Fickian行为，其中幂律拖尾是一个典型特征，其溶质浓度衰减遵循 ( C \propto t^{-n} ) 的形式，n为拖尾指数。前人已提出多种机制解释幂律拖尾，如裂隙壁粗糙度导致的流速空间变异、裂隙网络中速度分布的慢衰减、以及基质多速率扩散等。然而，大量现场观测发现拖尾指数 n 常常小于1.5，这表明部分溶质可能永远无法从系统中排出。传统观点认为，对于溶质从裂隙表面垂直扩散进入无限大静止基质的情况，n的理论值为1.5。因此，n < 1.5现象的成因一直是开放性的科学问题。尽管近期研究提出了如速度分布慢衰减、大尺度径向扩散等机制，但裂隙与基质之间不仅存在扩散，还可能存在由压力差驱动的对流（即通量交换）。这种通量交换在多大程度上影响乃至主导拖尾过程，此前尚不明确，尤其是在真实三维粗糙裂隙的背景下。因此，本研究旨在填补这一知识空白，验证“裂隙与基质之间的对流交换通量可导致 n < 1.5 的严重拖尾”的核心假设。
详细的工作流程 本研究采用了基于高分辨率数值模拟（Direct Numerical Simulation, DNS）与随机升尺度建模相结合的系统性方法，工作流程主要包括以下几个步骤：
构建高分辨率三维裂隙-基质模型系统：研究以一块天然圣安娜凝灰岩裂隙为研究对象，该裂隙尺寸为80.3 mm × 94.8 mm，先前已通过高分辨率X射线计算机断层扫描（CT）获取其精确的三维形貌数据，空间分辨率极高（平面内0.038 mm²，法向0.05 mm）。基于此扫描数据，研究团队构建了两个模拟域：a) 纯单裂隙模型；b) 裂隙-基质耦合模型。后者通过在扫描得到的粗糙裂隙几何体外围包裹一个约100 mm高的多孔介质（基质）块体而构成，旨在研究两者的相互作用。裂隙的平均开度为2.14 mm，标准偏差为0.86 mm，粗糙度为0.40。
直接数值模拟裂隙-基质系统中的流体流动：采用商业软件COMSOL Multiphysics进行三维稳态流体流动模拟。对于裂隙域，直接求解Navier-Stokes方程；对于多孔基质域，求解Brinkman方程，该方程能够在自由流（裂隙）与Darcy流（基质）的界面处准确捕捉粘性剪切耗散。在裂隙的进出口施加固定的压力梯度以驱动流体从左向右流动，所有其他边界设为无滑移壁面。裂隙与基质界面处的压力和速度完全耦合，即两者相等。为研究通量交换强度的影响，系统性地改变基质渗透率 ( K_m )，范围从 ( 10^{-16} \, m^2 ) 到 ( 10^{-9} \, m^2 )，涵盖了七个数量级，共模拟12个不同渗透率情景，可能对应自然界中从双孔隙到双渗透率的各类裂隙-基质系统。模拟采用了自由四面体网格单元，针对裂隙区域采用精细网格（单元尺寸 ( 8\times10^{-5} \, m ) 至 ( 5\times10^{-4} \, m )），基质区域采用较粗网格（约 ( 1\times10^{-3} \, m )），最终裂隙-基质模型的网格单元数超过3.3亿个，确保了数值结果的收敛性与可靠性。
直接数值模拟溶质运移过程：在已求解的稳态流场基础上，进行溶质运移模拟。在裂隙域，求解平流-扩散方程；在基质域，求解平流-弥散/扩散方程。溶质的分子扩散系数设为固定值。对于仅考虑基质扩散（无对流交换）的情景，将基质中的流速设为零，此时基质内的方程简化为扩散方程。边界条件设置为：在裂隙入口处施加狄利克雷边界条件（浓度C=1），出口为开放边界（无扩散通量），所有其他边界均为无通量边界。每个溶质运移模拟耗时约1-2个实际工作日。
拖尾行为分析与指数 n 的量化：从数值模拟结果中提取出口处的通量加权浓度随时间的变化，即穿透曲线。为了更好地分析幂律拖尾行为，研究进一步计算了停留时间分布，即浓度对时间的导数在半对数坐标下的图形。对所有模拟结果，在一个特定的时间窗口内（约3至200个孔隙体积），拟合幂律函数 ( C \propto t^{-n} )，从而得到量化拖尾程度的指数 n。
基于连续时间随机游走（CTRW）的随机升尺度：为检验是否能用升尺度模型捕捉复杂的拖尾行为，研究应用了CTRW这一随机建模框架。具体方法是：从DNS结果中提取裂隙上游某位置（x=10 mm）的穿透曲线，将其视为溶质粒子移动一个固定步长所需过渡时间的累积分布函数。然后，通过逆变换采样算法从此分布中随机抽取一系列过渡时间，用于模拟大量溶质粒子（本研究为100万个）的随机运动，从而预测下游位置（x=80 mm）的停留时间分布，并与DNS结果进行对比。此方法新颖之处在于首次将其应用于考虑裂隙-基质通量交换的粗糙壁裂隙运移问题。
边界条件敏感性分析：为探究模型边界条件设置（如基质是否也具有平行于裂隙的纵向压力梯度）对结论的影响，研究额外进行了一组DNS模拟。在该模拟中，除了裂隙的进出口，也对基质的进出口施加了特定的压力和浓度边界条件。通过对比两种边界条件下（仅裂隙有压降 vs. 裂隙和基质均有压降）的穿透曲线和拖尾指数，评估边界条件对研究主要发现的影响。
主要结果 1. 流体流动与通量交换特征：在固定的整体压力梯度下（雷诺数Re=1.83，处于线性流态），单裂隙内由于粗糙度影响形成了明显的优势流道。在裂隙-基质系统中，压力梯度驱动流体在裂隙入口附近的区域从裂隙流入基质，流经一段距离后再返回裂隙，形成了裂隙与基质间的通量交换。随着基质渗透率 ( K_m ) 的增加，这种交换通量线性增加。一个有趣但常被忽视的现象是，裂隙出口的总流量（出流）也随着 ( K_m ) 的增加而单调线性增加，这验证了可渗透基质会增强裂隙-基质界面处的速度，从而提高整体流量的理论。
基质扩散与对流交换对拖尾的独立与联合影响：
CTRW升尺度模型的成功应用：对比研究表明，基于上游过渡时间分布的CTRW模型能够非常准确地预测不同 ( K_m ) 情景下（包括无基质、仅扩散、以及高渗透率导致强交换的情况）在下游观测到的溶质拖尾行为。这证明了CTRW模型具有捕捉裂隙-基质耦合系统中非Fickian运移的升尺度能力，也暗示该系统的拖尾行为主要由扩散-弥散过程主导，而连续跳跃间的速度相关性在低Peclet数条件下可以忽略。
边界条件影响的有限性：边界条件敏感性分析表明，即使在基质进出口也施加压力梯度的更复杂情景下，所得的穿透曲线和停留时间分布与原始模型（仅裂隙有压降）在整体上基本一致。虽然晚期拖尾部分存在细微差别，但由此引起的拖尾指数 n 的相对变化（约1.3%）远小于由改变基质渗透率所引起的变化。这强有力地支持了本研究的主要结论：裂隙与基质间的通量交换，而非特定的边界条件设置，主导了晚期拖尾行为。
研究的结论与意义 本研究的主要结论是：裂隙与岩体基质之间的通量交换（包括扩散和对流）是控制溶质运移晚期幂律拖尾行为的关键机制。特别是，当基质渗透率足够大（超过约 ( 10^{-12} \, m^2 ) ）时，对流交换通量成为主导，能够导致拖尾指数 n 显著下降至1.5以下，这为野外普遍观测到的严重拖尾现象（n < 1.5）提供了一种新的、机理性的解释。同时，研究首次成功地将连续时间随机游走模型应用于此类复杂系统，证明了其对拖尾行为进行有效升尺度预测的能力。
科学价值与应用价值：本研究深化了对裂隙岩体复杂运移过程的理解，强调了在预测非Fickian运移时，必须充分考虑基质渗透性及由此产生的裂隙-基质对流交换的重要性。这对于涉及污染物迁移、地下水修复、能源资源开采（如地热、页岩气）等领域的高精度数值模型构建具有重要指导意义。研究建立的 n 与基质渗透率之间的定量关系，为快速评估拖尾严重程度提供了一种潜在的简化工具。此外，将高精度DNS与CTRW升尺度相结合的方法论，为处理其他多尺度、多物理场耦合的渗流问题提供了范例。
研究的亮点 1. 研究视角新颖：首次在三维真实粗糙裂隙的框架下，系统且定量地揭示了裂隙-基质对流交换通量对溶质运移晚期幂律拖尾的决定性影响，并提出了解释 n < 1.5 现象的新机制。 2. 方法先进全面：结合了超高分辨率三维直接数值模拟（DNS）与随机理论升尺度（CTRW）两种方法，从机理揭示到模型预测，形成了完整的研究链条。DNS的规模（超过3.3亿网格单元）和精细度（同时求解裂隙Navier-Stokes流和基质Brinkman流）在同类研究中处于前沿。 3. 参数范围广泛：对基质渗透率进行了跨越七个数量级的系统性参数研究，涵盖了自然界可能存在的多种岩石类型和裂隙-基质相互作用模式，结论具有普适性。 4. 结论明确有力：明确指出了控制拖尾的渗透率阈值（~( 10^{-12} \, m^2 ) ），并建立了 n 与 ( \log_{10}(K_m) ) 的定量关系，同时通过边界条件敏感性分析增强了结论的稳健性。 5. 开创性应用：首次将CTRW升尺度模型应用于考虑通量交换的粗糙壁裂隙运移问题，并验证了其有效性，拓宽了该理论的应用范围。