🔬 physics

Exploring Particle Clogging Mechanisms through Deterministic Pore Network Modeling

本研究开发并验证了一个确定性孔隙网络建模框架，用于模拟多孔介质中的颗粒堵塞动力学，揭示了颗粒尺寸和电导率比如何影响渗透率下降，同时识别了当前的局限性并为地下水修复和过滤领域的应用提出了未来的改进方案。

原作者： Qianjing Tang, Amir Raoof

发布于 2026-07-07

📖 1 分钟阅读☕ 轻松阅读

原作者： Qianjing Tang, Amir Raoof

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

想象一种多孔材料，比如海绵或一层沙子，不要把它看作一个实心的块状物，而要把它看作一座由微小隧道组成的巨大且复杂的城市。有些隧道是宽阔的大道，而有些则是狭窄的小巷。这就是科学家所说的“多孔介质”。

这篇论文关于建立这个城市的数字孪生（Digital Twin），以观察当微小的旅行者（颗粒）试图穿过它时会发生什么。有时，这些旅行者会被卡住，从而堵塞隧道。这被称为“堵塞（clogging）”。乌得勒支大学的 Qianjing Tang 和 Amir Raoof 开发了一个计算机模型，用以精确预测这些堵塞是如何以及何时发生的，以及它们如何改变系统中的水流。

以下是他们研究结果的拆解，使用了简单的类比：

1. 数字城市（模型）

与其尝试模拟一个巨大的真实海绵中每一滴水和每一个颗粒（这会让超级计算机运行很久），作者将海绵简化为了一个相互连接的管道网络。

管道： 这些代表了材料中的“喉咙”（狭窄的连接处）和“孔隙”（较宽的腔室）。
交通： 水流通过这些管道，带着微小的颗粒随之移动。
目标： 他们想要观察交通拥堵（堵塞）是如何形成的，以及这些拥堵如何减慢整个城市的交通流量（渗透率）。

2. “交通拥堵”规则（确定性堵塞）

在以往的模型中，科学家经常使用掷骰子（随机性）来猜测何时会发生交通拥堵。这篇论文引入了一套严格的规则手册（确定性准则）。

把它想象成一个收费站。模型会精确计算需要有多少辆车（颗粒）通过特定的狭窄隧道，该隧道才会完全被堵死。这不是一种猜测；这是一个基于数学计数的逻辑。一旦达到这个数值，该隧道就被视为“堵塞”。

3. “滤饼”效应

这是他们发现的一个巧妙之处：即使一个隧道被堵塞了，也不一定意味着所有的交通都停止了。

类比： 想象一条高速公路车道被一堆沙子堵住了。汽车无法从沙子中穿过，但如果还有空间，它们仍然可以绕着沙子行驶；或者沙子本身可能足够多孔，允许少量车辆缓慢挤过去。
发现： 模型显示，当一个隧道被堵塞时，颗粒会形成一个“滤饼（filter cake）”。水仍然可以流过这个滤饼，但颗粒会被拦截。这意味着隧道并没有完全废弃，它只是变成了一个缓慢的、经过过滤的车道。

4. 大小决定游戏

研究人员测试了当“旅行者”（颗粒）的大小不同时会发生什么。

小旅行者（1-2 微米）： 它们就像微小的蚂蚁。它们可以钻过许多狭窄的小巷。当它们被卡住时，它们会先堵塞细小的路径，但大的干道仍保持开放。交通会逐渐减慢，就像一个缓慢的渗漏。
大旅行者（4-5 微米）： 它们就像送货卡车。它们会在狭窄的小巷中非常快地被卡住。因为它们如此迅速地封锁了主要路线，整个城市网格可能会被切断成孤立的岛屿。交通停止得更快，且“堵塞”会停留在入口处（入口处），而不是深入到城市的内部。

核心洞察： 较大的颗粒堵塞系统的速度更快，但在某些情况下，与较小的颗粒相比，它们造成的总流量下降反而没那么剧烈，因为它们在入口处就迅速封锁了，导致系统的其余部分甚至没有机会变得混乱。这是在堵塞速度与减少流量程度之间的一场复杂舞蹈。

5. “交通流”的变化

随着堵塞的发生，水的流动并不只是均匀地减慢。它变得混乱了。

类比： 想象一条河流。如果你在河里放几块石头，水会在石头周围流得更快。如果你筑起一座大坝，水就必须寻找新的路径。
发现： 随着隧道被堵塞，水会在少数剩余的开放隧道中加速，而在被堵塞的隧道中减速。由于最快路径和最慢路径之间的差异（称为“变异系数”）变得巨大，流量变得非常不均匀且难以预测。

6. 模型目前还做不到的事情（局限性）

作者坦诚地说明了他们的数字城市目前缺失了什么：

“粘性”因素： 模型假设颗粒只有在物理上堵塞管道时才会卡住。但在现实中，颗粒可能会在堵塞之前就粘在管道壁上，就像架子上的灰尘一样。模型目前尚未考虑这种“预堵塞”的粘性。
“混合装”问题： 模型目前只能处理一次一种尺寸的颗粒（比如只处理红色的车）。现实生活中通常会有大卡车和小汽车的混合物（双峰混合物），而模型目前还无法模拟这种情况。

总结

这篇论文展示了一种新的、基于规则的计算机模型，它就像是微观海绵的交通模拟器。它表明：

堵塞从小处开始（在狭窄的隧道中）并向外扩散。
较大的颗粒堵塞更快，但它们产生的流量模式与较小的颗粒不同。
被堵塞的隧道并非完全失效； 它们充当了过滤器，让水流过但拦截颗粒。
流量变得混乱， 随着系统堵塞，水会在少数剩余的开放路径中快速奔涌。

作者希望这一工具能帮助工程师设计更好的过滤器、清理地下水，或通过预测系统何时以及如何被堵塞来改进石油开采。

技术摘要：通过确定性孔隙网络建模探索颗粒堵塞机制

问题陈述
颗粒在多孔介质中的堵塞是一个关键过程，它会改变渗透率和流体流动，影响从地下水修复、过滤系统设计到地热能作业等广泛的应用领域。虽然直接数值模拟 (DNS) 和格子玻尔兹曼方法 (LBM) 可以提供高分辨率的见解，但其计算受到限制，仅适用于较小的区域，可能无法代表代表性单元体积 (REV)。孔隙网络建模 (PNM) 为中尺度分析提供了一种计算高效的替代方案，但通常依赖于随机的、基于概率的标准来确定孔隙喉道的堵塞时机。这些随机方法为预测引入了不确定性。此外，现有模型往往难以准确捕捉一旦开始堵塞后，流体通量与颗粒通量之间的动态相互作用，特别是在部分堵塞区域内流体通量与颗粒通量分离的问题。

方法论
本研究利用 PoreFlow 框架（由 Raoof 等人于 2013 年开发的综合 PNM 工具）来模拟颗粒堵塞动力学。该方法建立在三个核心组件之上：

确定性堵塞准则： 模型并未依赖随机概率分布，而是采用了基于 Wyss 等人 (2006) 提出并由 Dersoir (2015) 完善的临界颗粒数 $N^*$ 的确定性准则。 $N^*$ 定义了堵塞特定喉道所需的颗粒数量。模型根据以下因素动态计算每个喉道的 $N^*$ ：
- 附着势： 通过 DLVO 理论进行估算，并受离子强度（例如 NaCl 浓度）的影响。
- 几何因素： 特别是喉道宽度 ( $W$ ) 与颗粒直径 ( $D$ ) 的比例，遵循幂律关系。
- 流速： 引入了阈值效应，即高于临界值 ( $Q_c$ ) 的流速会影响 $N^*$ 。
- 模型假设一旦通过某一喉道的累积颗粒数超过 $N^*$ ，该喉道即转变为“堵塞状态”。
流体-颗粒解耦（自过滤）： 在发生堵塞时，模型将流体流动与颗粒通量分离。虽然流体仍会通过被堵塞的喉道（尽管电导率降低），但颗粒通量被拦截，从而有效地将下游流体中的颗粒浓度降至零。这被建模为在喉道内生长出的“滤饼”或塞子。该塞子的长度 ( $L_{cl}$ ) 根据截留颗粒的体积进行计算，并且喉道的水力电导率使用串联电阻模型（结合开放部分和堵塞部分）进行更新。
模拟设置： 在 $1000 \times 1000 \times 2000$ $\mu$ m 的区域内生成了一个 3D 孔隙网络（包含 1656 个孔隙，3082 个喉道）。模拟在单相流条件下运行，采用恒定入口压力 (2.3 kPa) 和恒定的入口颗粒浓度。研究系统地改变了颗粒尺寸 (1–5 $\mu$ m) 和堵塞电导率比（即堵塞状态下的电导率与开放状态下的电导率之比）来进行敏感性分析。

主要贡献

确定性框架： 本研究建立了一个透明、模块化且确定性的喉道堵塞准则，在保持物理可解释性的同时，减少了与随机模型相关的不确定性。
动态电导率更新： 它实现了一种机制，使堵塞的喉道保持部分渗透性，从而允许模拟“滤饼”效应，即即使颗粒传输停止，流体流动仍然持续。
敏感性分析： 该工作对颗粒尺寸和堵塞喉道电导率比如何影响堵塞的演化过程以及最终的渗透率变化进行了系统性的调查。

结果

渗透率演化： 模型成功捕捉到了与实验观察一致的渗透率指数级下降。堵塞从较窄的喉道开始，并随着时间的推移向较大的喉道传播。
颗粒尺寸效应：
- 较大颗粒： 加速了堵塞过程，导致电导率迅速丧失并形成“完全堵塞”状态，此时网络连通性受到显著破坏。这些场景导致堵塞集中在入口附近（平均堵塞深度较浅），并随着流体重新分布导致速度异质性（变异系数，$CV$）下降。
- 较小颗粒： 导致“部分堵塞”状态，其中电导率下降较为缓慢。这些场景在较长时间内保持较高的网络连通性，并由于优先流路径的形成表现出持续增加的速度异质性（$CV$）。
电导率比的影响： 较高的电导率比（表明堵塞后的喉道仍允许显著流量）会延迟堵塞的进展。这种效应在具有快速堵塞潜力（如较大颗粒）的场景中更为明显，表明颗粒尺寸与堵塞塞子渗透性之间存在关键的相互作用。
空间分布： 随着堵塞的进行，平均堵塞深度减小（向入口移动），且堵塞深度的方差减小，表明堵塞事件向入口点集中。

意义与局限性
论文认为，该 PNM 框架为预测多孔介质中的堵塞行为提供了一个强大的工具，在地下水修复、过滤设计和地热能领域具有广泛的应用前景。确定性方法增强了对孔隙尺度现象及其对宏观尺度属性集体影响的可解释性。

然而，作者也谦虚地承认了特定的局限性：

当前模型未考虑在发生完全堵塞事件之前颗粒向喉道壁的附着，这可能低估了总体的颗粒截留量。
该模型目前仅限于单分散颗粒和单相流，缺乏处理多分散颗粒混合物或多相流条件（例如在非饱和带中）的能力。
研究建议，未来的改进应纳入颗粒附着动力学以及针对多分散系统的概率场景，以更好地反映现实世界的复杂性。