TracerLPM（版本1）：用于解释环境示踪剂数据中地下水年龄分布的Excel®工作簿

___Y1

已于 2024-02-23 14:40:33 修改

阅读量1.2k

点赞数 12

文章标签：数据库

于 2024-02-23 11:59:32 首次发布

本文链接：https://blog.csdn.net/weixin_44259522/article/details/136249149

版权

By Bryant C. Jurgens, J.K. Böhlke, and Sandra M. Eberts

Abstract

TracerLPM是一个交互式的Excel®（2007或更高版本）工作簿程序，用于通过使用集总参数模型（LPMs）评估环境示踪剂数据中的地下水年龄分布。集总参数模型是基于简化的含水层几何和流动配置的运输数学模型，考虑了含水层、井孔或排水区内的流体动力学弥散或混合效应。

工作簿中包括五个主要的LPMs：活塞流模型（PFM），指数混合模型（EMM），指数活塞流模型（EPM），部分指数模型（PEM）和弥散模型（DM）。可以通过以不同的组合结合主要的LPMs来创建二元混合模型（BMM）。通过在参数中包括未饱和带中的旅行时间，可以进一步扩展模型。

TracerLPM还允许用户输入从其他方法确定的年龄分布，例如数值地下水流模型的粒子跟踪结果或不包括在此程序中的其他LPMs。同时解释年轻地下水（指示20世纪40年代后的补给的人为大气气体和同位素物质）和更古老地下水（碳-14和氦-4）的示踪剂，从而可以对具有广泛年龄范围的样本的地下水年龄分布进行约束。

TracerLPM被组织成一个全面的解释性方法，包括水文地质概念化、数据和模型的视觉检查以及最佳参数估计。可以通过比较两种方式评估地下水年龄分布：（1）同时分析多个示踪剂，用于检查其在建模浓度方面的一致性（示踪剂-示踪剂应用），或（2）评估示踪剂时间序列数据以确定其与建模趋势的一致性（示踪剂-时间应用）。对于仅在某一时间点测量单个示踪剂的样本，也可以获得地下水年龄估计；然而，通常需要先了解适当的LPM，因为平均年龄通常不是唯一确定的。

LPM的输出浓度取决于模型参数和样本日期。所有LPMs都有一个平均年龄参数。EPM、PEM和DM还有一个描述样本中年龄混合程度的额外参数。BMMs有一个用于混合中第一个组分的分数的参数。一个LPM以及其参数值一起提供了包含在样本中的每个补给年龄的年龄分布或水分的分数贡献的描述。对于PFM，年龄分布是一个在一个明确年龄上的单位脉冲。对于其他LPMs，年龄分布可以更广泛，跨足十年、世纪、千年或更长时间。对于具有混合地下水年龄的样本，示踪剂数据的报告解释包括LPM名称、平均年龄以及任何其他独立模型参数的值。

TracerLPM还可以用于模拟井、泉、河流或其他地下水排放接受体对引入的非点源污染物（如硝酸盐）的响应。这是通过将LPM或用户定义的年龄分布与水位的污染物负荷信息相结合来实现的。有关历史污染物负荷的信息可用于帮助评估模型匹配实际条件的能力并了解观察到的污染物趋势，而有关未来污染物负荷情景的信息则可用于预测潜在的污染物趋势。

Introduction

环境示踪剂数据的收集已经成为世界各地地下水质量调查的例行工作。示踪剂提供关于地下水年龄、补给速率、补给源、地下水团的混合以及地下水-地表水相互作用的信息。从示踪剂测量中获得的信息提供了对地下水流系统的更完整的图像。

通常情况下，环境示踪剂测量的目的是为了确定地下水的年龄。在某些情况下，通过将样品中测量的示踪剂浓度与充入含水层的水中的示踪剂输入历史相关联，假定示踪剂从补给区到排水区（例如井）运移时不受流体动力学弥散或混合的影响，来估算地下水年龄。因此，含有示踪剂的水被假定具有一个明确的年龄。这个年龄通常被称为活塞流年龄或表观年龄。

然而，从井中抽取的水往往包含许多不同年龄和补给历史的水包裹，因为井通常横跨几英尺或更多的含水层。此外，异质含水层中的弥散和混合可能导致即使在屏幕间隔较短的井中也出现年龄分布较广的情况（Weissmann等人，2002）。因此，在许多水文地质环境中，假设活塞流是不太可能成立的。这些过程也适用于泉和河流的地下水排放调查。

自20世纪50年代以来，已经开发并应用了几种数学模型或集总参数模型（LPMs）来解释环境示踪剂（Vogel，1967；Eriksson，1971；Maloszewski和Zuber，1982；Amin和Campana，1996）。这些模型基于简化的含水层几何和流动配置，考虑了含水层或井孔内的弥散和混合效应。在实践中，模型是基于其对正在研究的含水层和井（泉）的概念相关性而应用的。尽管这些模型已经被证明能够解释不同水文地质环境中观察到的示踪剂浓度，但使用集总参数模型解释环境示踪剂并不是例行的。

最近（2011年），美国地质调查局的国家水质评估（NAWQA）计划研究了对比水文地质环境中向公共供水井输送自然和人为污染物（TANC）的问题（Eberts等人，2005）。首次完成的四项研究分别位于加利福尼亚州莫德斯托附近的东圣华金谷的中央山谷含水层；佛罗里达州坦帕附近的弗洛里达含水层；康涅狄格州伍德伯里附近的庞佩劳格河流域；以及内布拉斯加州林肯附近的高原含水层。TANC研究在每个研究区域为地下水年龄测定收集了多个示踪剂，并使用集总参数模型（LPMs）解释了示踪剂浓度。这些结果与应用于每个研究区域的详细局部尺度地下水流模型的粒子跟踪方法的结果进行了比较（Eberts等人，2012）。该研究的主要发现包括：（1）表观（活塞流）年龄经常是对来自公共供水井和一些短屏幕监测井的水的平均年龄的误导性描述；（2）基于年龄分布的建模污染物响应可能与基于表观年龄或平均年龄的响应大不相同；（3）当基于相似的概念模型并根据相似的示踪剂数据进行校准时，LPMs可能为混合年龄水的井提供与粒子跟踪结果类似的年龄分布。因此，通过用测量的示踪剂浓度评估LPMs相对于开发三维地下水流和运移模型而言，这是一种相对简便和经济的方法，可以获得关于井中地下水样本的年龄分布的有价值信息。

目前可用于根据环境示踪剂数据估算地下水年龄的程序包括FLOWPC（Maloszewski和Zuber，1996），BOXMODEL（Kinzelbach等人，2002；Zoellmann等人，2002），LUMPED和LUMPEDUS（Ozyurt和Bayari，2002；2005），以及TRACERMODEL（Böhlke，2006）。这些程序主要设计用于年轻地下水的示踪剂，如氚（3H），氯氟烃（CFCs）和六氟化

硫（SF6）。其中一些程序缺少一个或多个LPM。BOXMODEL包含弥散模型，但不包括指数活塞流模型或二元混合模型。TRACERMODEL的初始版本包含指数活塞流模型和二元混合模型，但没有弥散模型。此外，其中一些程序缺少用于协助校准在样品中测量的多个示踪剂的最小化算法。

本报告中介绍的程序，TracerLPM，基于工作簿程序TRACERMODEL（Böhlke，2006），该程序已在美国地质调查局进行的多个研究项目中使用过（Böhlke和Denver，1995；Focazio等人，1998；Katz等人，1999；Katz等人，2001；Plummer等人，2001；Böhlke，2002：Böhlke和Krantz，2003；Lindsey等人，2003；Landon等人，2008；Jurgens等人，2008；Katz等人，2009；Brown等人，2009）。TRACERMODEL工作簿强调对多个示踪剂进行同时建模，以评估样本的平均年龄和年龄分布，并具有“示踪剂-示踪剂”图的图形分析功能。TRACERMODEL的另一个工作簿版本是为比较示踪剂时间序列数据和模型而开发的。在新的工作簿程序TracerLPM中，对TRACERMODEL进行了许多改进和改进，包括增加了额外的二元模型、老年地下水示踪剂（碳-14、氦-4）、最小化方案和示踪剂浓度的预测。

本报告的目的是记录Excel®工作簿程序TracerLPM，该程序旨在基于地下水中的环境示踪剂评估地下水年龄分布和平均年龄，适用于来自井和泉的样品。本报告分为四个部分。第一部分描述了模型及其在特定水文地质环境中的适用性。第二部分描述了当前包含在此工作簿中的环境示踪剂。第三部分概述了工作簿并详细描述了各个工作表。第四部分包含说明程序的使用和功能的示例。

Lumped Parameter Model Calculations

TracerLPM工作簿包含五个可用于确定样本的年龄分布和平均年龄的LPMs：活塞流模型（PFM）、指数混合模型（EMM）、指数活塞流模型（EPM）、部分指数模型（PEM）和分散模型（DM）。这些模型中的每一个都可以与另一个模型结合，以创建二元混合模型（BMM）。这导致了25种可能的BMM组合。这些LPMs对应于地下水从含水层的入口位置（补给区）流向含水层的出口位置（井或泉），并在数学上表示为过渡时间分布函数或出口年龄分布函数 [g(t)]（Maloszewski和Zuber，1982）。水样被设想为由许多“小包裹”组成，它们沿不同的流动路径到达样本点；每个小包裹代表相对离散的地下水年龄和示踪剂浓度。对于稳态地下水系统，可以使用出口年龄分布函数和示踪剂的衰减函数（例如放射性衰变或生物降解）从含水层入口位置的示踪剂输入历史来计算含水层出口位置的模拟示踪剂浓度：

$C_{out}=\int_{-\infty }^{t}C_{in}(t')e^{-\lambda (t-t')}g(t-t')dt' \qquad (1)$

其中 $C_{out}$ 是出口示踪剂浓度， $C_{in}$ 是入口示踪剂在时间 $t'$ 的浓度， $t$ 是样本日期， $t'$ 是水包裹进入系统的日期， $\lambda$ 是衰减常数，即单位时间的分数损失， $t-t'$ 是水包裹的年龄。

模型集总参数方法假定示踪剂被注入并在流体通量中检测到，并且示踪剂在行为上是保守的（除了可能的放射性衰变或降解），并随水流动。因此，由示踪剂浓度推断出的平均年龄等于从系统或筛选间隔中排放的水的平均年龄。对于一些示踪剂或水文地质情境，示踪剂行进时间与水行进时间之间的偏差可能导致对样本中水的年龄分布和平均年龄的不可靠估计。因此，重要的是要考虑扩散和地球化学交换等过程如何影响那些不属于水分子的特定示踪剂，并导致基于示踪剂的年龄与水龄之间存在较大差异。

样本的平均年龄（τs）是从描述样本中示踪剂浓度的出口年龄分布函数中导出的：
$\tau _s=\int_{-\infty }^{t}(t-t')g(t-t')dt \qquad(2)$
其可以通过数值逼近：
$\tau_s =\sum_{i=1}^{\infty }t_iX_i(\Delta t)\qquad (3)$
其中， $t_i$ 是水包裹的年龄 $( t-t' )$ ， $X_i$ 是由给定年龄增量对应的水包裹在样本中所占比例， $\Delta t$ 是时间步长（年龄增量）。

通常情况下，样本的平均年龄 $\tau_s$ 不等于整个含水层中地下水的平均年龄 $\tau_{aq}$ ，但在一些情况下，比如在EMM或EPM中，样本的平均年龄也对应于整个含水层中地下水的平均年龄。

对于每个LPM年龄分布函数 $g(t - t')$ ，方程1的解析闭合形式的解都是通过使用名为“TracerLPMfunctions”的Excel插件（为32位和64位版本的Excel编译）导出的。该XLL插件包含了LPM的工作表函数，使用C++语言在Microsoft® Visual Studio®和Microsoft Office Excel® 2010软件开发工具包中进行编程（有关更多详细信息，请参见附录B中的兼容性部分）。

活塞流模型（PFM）

活塞流模型（PFM）假设示踪物质从入口位置（补给区）流向出口位置（例如井或泉），不经历水动力弥散或混合。PFM可适用于水文地质环境中弥散较低、平均线性流速较高或从补给到排泄的流动路径较短的情况。在不受限制的含水层中，从浅层、短筛选的监测井或在有小型补给区的受限含水层中测得的示踪物质可以近似地遵循活塞流行为（见图1）。活塞流模型的出口年龄分布函数如下：

$PFM_{g(t-t')}=\delta(t-t'-\tau_s)\qquad(4)$

其中
$\delta$ 是狄拉克（Dirac）δ函数。

通过以下公式计算PFM：

$C_{out}(t)=\begin{cases}C_{in}(t-\tau_s)e^{-\lambda-\tau_s}, & \text{for } t=\tau_s\\ 0, & \text{for } t \neq \tau_s\end{cases}\qquad(5)$

其中 $\tau_s$ 是延迟时间， $C_{out}(t)$ 是出口处的浓度， $C_{in}(t)$ 是入口处的浓度， $\lambda$ 是衰减常数。

图1.理想化水文地质含水层结构的示意图，其中活塞流模型可能适用。顶部图示一个接受面源补给的不受限制的含水层，其中短筛选的井位于含水层的浅部。第二个图（改编自Maloszewski和Zuber，1982）显示一个有小型补给区的受限制含水层，以及一个位于补给区下游的井。假设弥散或扩散过程对示踪物浓度梯度几乎没有影响。右指向的箭头表示采样点。底部图显示了示踪物从补给区到井的活塞流输送所期望的出口年龄频率分布，平均年龄为25年，以1年递增。

指数混合模型（EMM）

指数混合模型（EMM）适用于均匀、不受限制、厚度恒定的含水层，该含水层接收均匀的补给（见图2）。

这种情况导致地下水年龄的垂直分层，从水位表为零开始，逐渐呈对数增长，直到在含水层底部趋近于无穷大（Vogel，1967；Appelo和Postma，1996）。该模型适用于完全贯穿井或向泉或河流排放的含水层。EMM模型适用于流线沿着纵向和横向均不发生弥散的情况，混合发生在井筒或泉眼内而不是在含水层内。EMM还描述了完全混合水库的年龄分布（Eriksson，1971）。EMM出口年龄分布函数如下：

$EMM_{g(t-t')}=\frac{1}{\tau_s}e^{(\frac{t-t'}{\tau_s})} \qquad(6)$

通过使用对于每个年龄增量（Δt）的卷积积分的以下闭式解，从采样日期减去年龄增量并向后迈进时间，直到输出浓度变化不超过10^{–6}为止，可以计算EMM：

$C_{out}(t) = \sum_{t'=-\infty}^{t-\Delta t}\frac{1}{\tau_s}\frac{1}{\frac{1}{\tau_s}+\lambda}\left[e^{-\left(\frac{1}{\tau_s}+\lambda\right)(t-t'-\Delta t)} - e^{-\left(\frac{1}{\tau_s}+\lambda\right)(t-t')}\right]\qquad(7)$
其中 $\tau_s$ 是延迟时间， $C_{out}(t)$ 是出口处的浓度， $\lambda$ 是衰减常数。

图2. 在指数模型可能适用的理想化含水层配置的示意图。顶部图示一个不受限制的含水层，接收均匀的补给，样品取自整个含水层厚度的筛选良好的井或泉。假设混合发生在井或泉中而不是在含水层中。右指向的箭头表示采样点。底部图显示了从井或泉中示踪物的指数混合所期望的出口年龄频率分布，平均年龄为25年，以1年递增。

Exponential Piston-Flow Model (EPM)

The exponential piston-flow model (EPM) can be used to describe an aquifer that has two segments of flow in series: a segment of exponential flow followed by a segment of piston-flow. This model can be used to describe discharge from an aquifer of constant thickness with an upgradient unconfined portion receiving areally distributed recharge (the exponential part) connected to a downgradient confined portion or an unconfined portion receiving little to no recharge (piston-flow part; fig. 3). The EPM also can be used to describe piston-flow transport within the unsaturated zone followed by exponential mixing; however, piston-flow transport through the unsaturated zone has been implemented separately for all models in TracerLPM. In addition, some tracers, like tritiogenic helium-3 for example, cannot be modeled correctly using an EPM to describe unsaturated zone (UZ) transport. Therefore, the EPM, as implemented in this program, is intended for situations where exponential flow precedes piston-flow within the saturated zone. The EPM age distribution function is as follows:

总平均年龄、平均年龄和非饱和带行程时间

在这个程序中，样本的总平均年龄被定义为饱和带地下水的平均年龄加上通过非饱和带（UZ）的行程时间。正如这暗示的那样，当通过非饱和带的行程时间为零或小于一年时，总平均年龄大致等于饱和带行程时间的平均年龄。对于BMMs，总平均年龄由以下公式给出：

$\tau_{total} = f_1\tau_1 + (1-f_1)\tau_2 + UZ_{tt}\qquad(22)$

其中：

$\tau_{total}$ 是混合物的总平均年龄，

$f_1$ 是第一组分的分数，

$\tau_1$ 是第一组分的平均年龄，

$\tau_2$ 是第二组分的平均年龄，以及

$UZ_{tt}$ 是非饱和带行程时间。

对示踪物的非饱和带输运建模时，假定了活塞流行为。因此，5年的非饱和带行程时间导致3H在降水形式进入地下后5年进入饱和带。由于3H是一种放射性示踪物，进入含水层的浓度可能明显低于降水中的初始3H浓度。此外，由于3Hetrit通常会在水包穿越水位并与大气交换隔离之前失去到大气中，非饱和带对3H的输运的长行程时间导致含水层中3Hetrit浓度较低的峰值。

重要的是要注意，在非饱和带中弥散可能会影响3H和其他示踪物的输运。

尽管TracerLPM中没有明确将非饱和带中的弥散编程为选项，但可以通过将弥散模型应用于原始示踪物输入历史并将模型结果输入到TracerInput工作表中，创建一个考虑弥散的示踪物输入函数。

TracerLPM允许用户为每个示踪物独立地将非饱和带行程时间设置为TracerInput工作表上指定的常数值。这使用户能够为在地下水中具有年龄不一致性的示踪物定义不同的非饱和带行程时间，这不能通过LPM模型参数的变化来解释。例如，3H和其他可溶性示踪物可以作为含水物种随渗透水一起穿越非饱和带，而相对不溶性的大气气体示踪物，如CFCs，可以在非饱和带空气中更快地到达水位。因此，为了正确地模拟地下水中的3H和CFC浓度，必须将水性示踪物与大气气体示踪物的非饱和带行程时间设置为不同的值。

TracerLPM还允许用户将一个或多个示踪物的非饱和带行程时间设置为在其他几个工作表中定义的非饱和带行程时间参数。这使用户在从输出和图形工作表测试方案时无需在TracerInput工作表上为每个示踪物单独更改非饱和带行程时间。每个示踪物的非饱和带行程时间的处理可以设置为TracerInput工作表上指定的常数值，或者设置为由其他工作表上的非饱和带行程时间参数控制。