地中海沿海含水层中与盐水入侵过程相关的地下水污染程度的多变量评估

《Journal of Contaminant Hydrology》:Multivariate assessment of groundwater contamination levels associated with saline intrusion processes in Mediterranean coastal aquifers

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

编辑推荐:

  连续时间随机游走粒子追踪模型用于描述地质多孔介质中渐近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。然而,对于相对较低的佩克莱特数(Péclet number),羽流形状有足够的时间在穿过测量位置时演变,这种情况下完美的高斯羽流也可能产生不对称的BTCs(Ogata和Banks,1961年)。
如果系统本身随时间演变,通过系统的传输性质也会发生变化。例如,沉淀和溶解(P&D)反应可以影响多孔介质的结构,从而导致速度场的变化。Cohen等人(2025年)测量了在一个60厘米长的沙填充柱中石膏沉淀程度增加时惰性示踪剂的BTCs,发现系统内沉淀程度的增加导致了异常传输程度的增加,表现为越来越不对称的重尾BTCs。这些BTC行为的趋势在实验和模拟中都有观察到,无论是以溶解为主导的系统还是以侵蚀为主导的系统(Deng等人,2025年;Saeibehrouzi等人,2025年;Trabucchi等人,2023年)。
如图1所示,将传输行为从高度异常到费克定律的平滑过渡(即没有单一时刻出现突变)概念化意味着可以将可能的传输行为范围视为一个连续谱,而不是费克定律/异常的二分法。这种基于谱的方法不仅在考虑逐渐延长的时间或距离时有用,而且也适用于研究动态系统。特别是,如果系统的初始状态是高度均匀的(例如沙填充柱),并且它经历了逐渐的变化(例如由于沉淀),则可以想象通过该系统的传输会表现出类似于图1中“放大”效应的效果,但方向相反:传输最初可能接近费克定律,然后逐渐变得更为异常,最终不可避免地进入渐近费克定律范围。
那么就出现了如何适当地对这样的动态系统进行建模的问题,其中整体异常传输程度可能最初(或始终)相对较低,且传输预计(至少在系统的某些状态下)是渐近费克定律的。使用基于连续时间随机游走(CTRW)方法的模型可以实现对广泛传输行为的建模,同时将整个域视为一个单一的有效介质。这类模型在概念上不同于传统的随机游走模型:虽然这两种类型的模型都使用模拟粒子来表示溶解或分散的化学物质包,并根据运动方程在系统内移动(即进行步进),但CTRW模型中的运动方程不仅控制了粒子的位置移动,还控制了时间移动(Berkowitz等人,2006年;Bijeljic和Blunt,2006年;Rhodes等人,2008年)。某些CTRW模型公式能够描述渐近费克定律传输,如下所述。
在第2.1节中,简要介绍了CTRW模型,不仅描述了本工作中使用的模型,还介绍了其发展的原始模型和概念基础。这样做既是为了说明为什么某些模型中的CTRW幂律指数(在第2.1.1节中介绍)被限制在特定范围内,也是为对CTRW建模感兴趣的读者提供可能有用的通用信息。第3.1节进行了使用CTRW的传输模拟,并将其与基于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号