探索多孔介质中渐近菲克(asymptotically-Fickian)化学传输现象

《Journal of Contaminant Hydrology》:Exploring asymptotically-Fickian chemical transport in porous media

【字体: 时间:2026年02月13日 来源:Journal of Contaminant Hydrology 4.4

编辑推荐:

  化学传输在多孔介质中的非Fickian特性及CTRW-PT模型分析,提出基于连续时间随机游走(CTRW)的谱观点方法,可描述渐近Fickian行为并有效模拟多孔介质中化学物质传输的复杂现象,通过数值模拟和实验验证了参数空间扩展的必要性,并定义了可用于评估异常传输程度的简便参数。

  
马坦·科恩|布莱恩·伯科维茨
以色列魏茨曼科学研究所地球与行星科学系,雷霍沃特7610001

摘要

化学物质通过地质多孔介质的传输行为多种多样,从费克定律遵循到高度异常行为都有,其中包括一种渐近费克行为,即只有在长时间和/或大距离下传输才符合费克定律。渐近费克传输行为在文献中仍研究不足,在许多情况下,传输行为要么被简单地视为费克定律遵循,要么被视为异常行为。本文提出了一种基于光谱的化学传输视角,能够更好地描述渐近费克传输的细微差别。采用连续时间随机游走(CTRW-PT)模型来描述多孔介质中的传输过程,该模型考虑了代表广泛传输行为的参数空间,所有这些行为都属于渐近费克范畴。随后使用相同的CTRW-PT模型分析了柱状实验的实验结果。这两组结果——纯模拟和实验建模——共同证明了允许使用宽参数空间进行CTRW-PT建模的有效性和实用性。此外,还评估了几个易于计算的参数,作为判断传输异常程度趋势的潜在指标,特别是在演变系统的背景下。

引言

化学物质通过地质多孔介质的传输受多种因素影响,包括介质的结构和组成、流体速度以及所涉及化学物质的性质(无论是溶解态还是分散态,以及其反应能力、形成复合物的能力或吸附到介质上的能力)。具体来说,对于高度均匀的介质,传输通常是费克定律遵循的,即可以通过平流-扩散方程(ADE)在足够大的空间和时间尺度上进行量化。对于更不均匀的系统,传输行为会偏离费克定律(Berkowitz等人,2006年),即使在非常轻微不均匀的系统中也能观察到这种效应(Zhang和Lv,2007年)。非费克传输通常表现为不对称的、长尾的化学羽流,而费克传输下的化学羽流通常是“更紧密”的,并以对称的高斯方式发展(即,从初始注入点开始的浓度分布沿任何测量轴都是高斯的)。使用基于ADE的模型对不均匀系统进行建模需要对该系统进行映射,即在整个系统体积内分配水力传导率值,然后显式求解速度场。然而,精确的几何结构通常是未知的且难以调查的,因此更倾向于使用能够产生不均匀性效应(例如长尾化学羽流)的模型,同时将整个域视为一个单一的有效连续体(Berkowitz等人,2006年;Naseri等人,2020年)。
在研究污染物通过地质多孔介质的传输时,一个重要的问题是污染物在一个地点(例如井、含水层区域)停留的时间长度,直到浓度降低到安全或可处理的水平(Dedewanou等人,2015年)。这个问题的答案在很大程度上取决于长尾化学羽流的出现,特别是在考虑砷或重金属等有毒污染物时,极低的浓度就足以使地下水无法饮用(Nouri等人,2008年;Podgorski和Berg,2020年)。因此,有必要开发和使用能够描述实际地质系统中可能出现的传输行为谱系的模型。
在结构不均匀性的最大尺度相对于感兴趣的尺度足够小的情况下,最初的非费克传输可能在长时间和/或大距离上趋于费克行为(Dentz等人,2004年;Marseguerra和Zoia,2007年;Xiong等人,2006年);这种行为在这里被称为渐近费克。直观上,可以将其视为观察不均匀域的一个子部分,并逐渐“放大”,使得各个结构元素变得越来越小,直到整个域在宏观上看起来是均匀的(图1)。通过这样的系统迁移的化学羽流在足够长的时间和/或距离上会表现出越来越像费克定律的传输行为,就像在均匀域中迁移一样。尽管这种效应在概念上很容易理解,并且可以在数值模拟中展示,但在实验中却很难观察到(Berkowitz等人,2000年;Cortis等人,2004年;Dentz等人,2004年)。这主要是因为展示传输趋于费克行为可能需要从注入点开始跨越非常宽的距离进行测量,而所有这些测量都在结构不均匀性有限的受控系统中进行。总体而言,关于渐近费克传输的已发表文献很少(Dentz等人,2018年)。
在这项工作中,非费克异常这两个术语可以互换使用,用来描述任何不能被ADE很好地描述的溶质传输行为,包括在有限时间内甚至“短”时间内。术语传输异常程度用于描述传输行为与费克定律的差异程度。通常,系统越不均匀,传输异常程度越大。
虽然费克传输和异常传输之间的根本区别在于羽流形状的演变方式,但实际适用于现实世界系统的测量通常是系统内某一点(或几个点)处的浓度-时间曲线(BTCs)。在一些实验工作中,已经报道了对化学羽流的直接成像(Marica等人,2011年;Trabucchi等人,2023年)。通常,BTCs反映了羽流形状,因此不对称的长尾羽流会产生不对称的重尾BTCs。然而,对于相对较低的佩克莱特数,羽流形状有足够的时间在穿过测量位置时演变,这种情况下完美的高斯羽流也可能产生不对称的BTCs(Ogata和Banks,1961年)。
如果系统本身随时间演变,传输的性质也会发生变化。例如,沉淀和溶解(P&D)反应可以影响多孔介质的结构,从而导致速度场的变化。Cohen等人(2025年)测量了在60厘米长的沙填充柱中石膏沉淀程度增加时惰性示踪剂的BTCs,发现系统内沉淀程度的增加导致了传输异常程度的增加,表现为越来越不对称的重尾BTCs。这些BTC行为的趋势在实验和模拟中都有观察到(Deng等人,2025年;Saeibehrouzi等人,2025年;Trabucchi等人,2023年)。
如图1所示,将传输行为从高度异常到费克定律的平滑过渡(即没有单一时刻出现突变)概念化意味着可以将可能的传输行为范围视为一个连续谱,而不是费克/异常的二分法。这种基于光谱的方法不仅在考虑逐渐延长的时间或距离时有用,而且对于研究演变系统也很有用。特别是,如果系统的初始状态是高度均匀的(例如沙填充柱),并且它经历了逐渐的变化(例如由于沉淀),那么通过该系统的传输可以想象为表现出类似于图1中“放大”效应的效果,但方向相反:传输最初可能接近费克定律,然后逐渐变得更加异常,最终不可避免地进入渐近费克区域。
那么就出现了如何适当建模这种演变系统的问题,其中整体传输异常程度可能最初(或始终)相对较低,并且传输预计(至少对于演变系统的某些状态)是渐近费克的。使用基于连续时间随机游走(CTRW)方法的模型可以实现对广泛传输行为的建模,同时将整个域视为一个单一的有效介质。这类模型在概念上不同于传统的随机游走模型:虽然这两种类型的模型都使用模拟粒子来表示溶解或分散的化学物质包,并根据运动方程在系统内移动(即过渡;步进),但CTRW模型中的运动方程同时控制了时间上的过渡和空间上的过渡——也就是说,CTRW模型描述了粒子何时移动以及在哪里移动(Berkowitz等人,2006年;Bijeljic和Blunt,2006年;Rhodes等人,2008年)。某些CTRW模型公式能够描述渐近费克传输,如下所述。
在第2.1节中,简要介绍了CTRW模型,不仅描述了本工作中使用的模型,还介绍了其发展的原始模型和概念。这样做既是为了说明为什么某些模型中的CTRW幂律指数(在第2.1.1节中介绍)被限制在某些值范围内,也是为了向对CTRW建模感兴趣的读者提供可能有用的通用信息。使用CTRW进行的传输模拟在第3.1节中进行并展示,并与基于ADE的费克传输进行了比较。研究了模型变量值的影响,并讨论了使用宽参数范围对于幂律指数的有效性。然后,在第3.2节中,基于上述模拟的洞察,对Cohen等人(2025年)的一些结果进行了CTRW-PT分析。最后,在第3.3节中,定义并回顾了基于BTCs的易于计算的传输参数,以确定演变系统中传输异常程度趋势的有用度量。Cohen等人(2025年)表明,其中两个参数——洗脱时间窗口和最后一个峰的斜率——是有用的。

章节摘录

CTRW基础和连续体方法

CTRW模型可以采用连续体方法应用,描述浓度场根据适当的微分方程如何发展,或者采用粒子跟踪(PT)方法,描述单个粒子如何通过系统移动(Berkowitz等人,2006年;Dentz等人,2004年)。历史上,首先使用PT概念开发了一个由积分-微分方程控制的连续体方法模型,后来又开发了一个纯PT模型。
早期的CTRW模型使用了随机

突破曲线

图2显示了第2.1.3节中描述的所有数值模拟的突破曲线,包括线性和半对数(以y为坐标)尺度;包含半对数尺度是为了强调BTC的尾部行为。由于本工作的重点是BTC的不对称性,如果BTC可以通过ADE类型的模拟复制,并且具有相同的初始条件(即初始粒子位置)和任何vD值,则认为BTC是费克定律遵循的。

结论

历史上,通过多孔介质的溶质传输大多被认为是费克定律遵循或异常的,而渐近费克传输(即在足够长的时间和距离上趋于费克定律的行为)相对较少被研究。基于CTRW的模型用于描述化学物质的传输,采用PDF来表示粒子过渡时间,这些PDF由指数β控制的幂律形式表示。最早的CTRW模型关注渐近行为,对于这些行为

CRediT作者贡献声明

马坦·科恩:撰写——原始草稿,可视化,验证,方法论,调查,形式分析,概念化。布莱恩·伯科维茨:撰写——审阅与编辑,验证,监督,资源提供,方法论,调查,资金获取,概念化。

利益冲突声明

作者声明以下财务利益/个人关系可能被视为潜在的利益冲突。
布莱恩·伯科维茨报告称得到了以色列科学基金会的财务支持。如果有其他作者,他们声明没有已知的财务利益或个人关系可能影响本文报告的工作。

致谢

这项研究得到了以色列科学基金会(资助编号1008/20)和De Botton海洋科学中心的支持。B.B.担任Sam Zuckerberg水文学教授职位。
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普

知名企业招聘

热点排行

    今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号