第二章：污染物传输（ACT）模块分析

2.1概述

ACT模块集成到MODFLOW-SURFACT/MODHMS中，用于对多达五种溶质的瞬态或稳态流场中的平流-分散输运方程进行数值求解。输运分析中使用的数值格式是严格质量守恒的。每种物种都允许线性或非线性（Freundlich）平衡吸附。允许在液相和土壤上进行一级生化降解。还实现了直链降解（或放射性衰变）和转化产物的生成。提供了仅允许在活动阶段进行运输的选项；气相和水相中的传输（其中一个对流动方程不起作用，其中扩散是传输机制）；以及在空气和水相中传输污染物，所述污染物从不动的消耗NAPL相溶解和挥发。与MODFLOW公式完全兼容的块中心有限差分空间离散化与最先进的TVD通量限制方案一起用于求解输运方程。本章讨论了运输模拟器的制定和实现。

2.2多物种运输问题的提法

2.2.1活动相输运的控制方程

控制污染物种类k在可变饱和多孔介质中的三维传输的偏微分方程可以用原始形式写成如下（Bear，1979）：

式中，Dij 为表观水力色散张量，ck 为活性相中 k 组分的溶质浓度，csk 为土壤吸附的 k 组分浓度，vi 为达西速度，φ 为有效孔隙度，Sα 为土壤饱和度。 活性流体相 α、λαk 和 λsk 分别是活性流体相 α 和土壤中组分 k 的一阶衰减系数，q 是每单位体积多孔介质通过源（或汇）的体积流量，c *k 是源（或汇）的溶质浓度，Γ 是组分 k 从活性相到非活性相的传质速率（仅对于活性相中的污染物传输为零，即，它显示在此处） ，以便稍后将方程推广到包括质量传递）。 分量 (j) 是发生变换产物时子产物 k 的父分量，方程 (1) 中的最后两项表示分量 k 的生成，这是所有父项 j = 1 衰减到 NPAR 的结果（其中 NPAR 是 k) 的亲本数量，分别在水中和土壤中。 xikj 是母体成分 j 转化为成分 k 的分数（这是直链衰减的统一，其中一个单位的母体产生一个单位的子体。）。 支链衰变详细信息将在第 2.2.6 节中提供。 活性相 α 被认为是求解流动方程的相，如第 1 卷第 2 章中所述。方程 (1) 表示由于吸附多孔介质中的平流和菲克色散而导致的衰变污染物的传输。 活跃期。
当活性流体相 α 为水 (w) 且 Sw=1（即，对于完全承压含水层或饱和含水层）时，可变饱和多孔介质中污染物传输的控制方程（方程 (1)）可简化为饱和传输方程 细胞）。

回想一下，BCF4 输入包中的索引 IREALSL 决定了流动模拟的活动阶段。 如果 IREALSL=0，则伪本构关系确定水活性系统中的地下水位。 如果 IREALSL= 1 或 2，则水在可变饱和系统中处于活性相，如果 IREALSL = 3 或 4，则空气 (a) 是可变饱和系统的活性相。 对于无约束条件（即，使用伪本构关系），地下水位上方的网格块的饱和值等于零（残差），并且在稳态条件下在干电池中发生传输，进入其中的质量通量等于质量 通量从底部流出，从而与这种情况下流动方程的 GSVE 考虑保持一致。 对于 0<Sw<1，流量和污染物迁移发生在无约束条件下，回想起来，当使用伪土壤关系确定地下水位时（即当 IREALSL=0 时），Sw 表示饱和网格块厚度。
活性相污染物浓度的存储期限（第一项）
等式（1）的右侧）针对瞬态传输的情况进行了扩展
不饱和流场，其中用于流动的比屈服（Sy）值不同于（小于）等式（1）的孔隙率（φ）值。 具体来说，污染物的储存发生在任何网格块的孔隙空间的两个区域内。 第一区域包括由Sy项占据的空间，第二区域包括剩余的孔隙空间(φ-Sy)。 例如，使用拟土函数的无侧限流解决方案将 Sy 视为可排水孔隙度，剩余孔隙空间被不可移动的水填充。 因此，扩展存储期限给出：

式中，Sachive是剩余孔隙空间中活性相的饱和度（φ-Sy）。这种膨胀是必要的，因为流体饱和度的变化是根据空间Sy计算的，这是由于净流入减去流出（用于传输方程中的达西速度）而发生的。剩余空间不包含在平流项中，并且可以任选地假设其被活跃相（Sactivity＝1）或被非活跃相（Sachive＝0）填充。利用基本传输输入读取标志IPHSFLG，以指示剩余空间是由活动（IPHSFLG=1）阶段还是由非活动（IPHSFLG=0）阶段填充。传输方程的原始形式的质量守恒性质，即使用平流项中的体积流体通量来移动污染物（而不是传输方程的简化形式所使用的速度），对这种物理情况是尖锐的。对于方程（1）的父项的衰变或增益，不进行类似的展开，因为当φ与Sy不同时，它们不构成质量守恒问题。