模型参数由与命运和传输建模相关的材料属性组成,这些建模与独立于要建模的组成物种的地质单元相关。 每个模拟所需的具体模型参数由延迟模型和反应模型定义。
通用模型参数:
IRCTOP
IRCTOP 是一个整数标志,指示如何输入反应变量:
· 按单元 (IRCTOP=2) [默认] - 所有反应变量均指定为逐个单元的 3-D 数组。
· 按层 (IRCTOP=0) - 所有反应变量均指定为一维数组,数组中的每个值对应于单个层。该选项主要是为了保留与MT3D之前版本的兼容性。
IRCTOP 不适用于 MODFLOW-6 模拟
IGETSC
IGETSC 是一个整数标志,指示非平衡吸附时是否应读取所有物质的非平衡吸附或固定相的初始浓度。
模拟 (ISOTHM = 4) 或双域传质 (ISOTHM = 5 或 6)。 该标志是通过回答以下问题来设置的:是否应该读取吸附相或固定相浓度?
· 否 - 则 IGETSC = 0,并且不会读取吸附相或固定相的初始浓度。 默认情况下,假设吸附相与溶解相平衡 (ISOTHM = 4),并且假设固定域的浓度为零 (ISOTHM = 5 或 6)。
· 是 - 则 IGETSC =1,并且将读取所有物质的吸附相或固定液相的初始浓度(如适用)。
IGETSC 不适用于 MODFLOW-6 模拟; 始终读取移动域(以及适用的固定域)中溶解溶质浓度的初始浓度,并且仅支持非动力学/平衡吸附模型。
延迟模型 - 模型参数
Bulk Density (RHOB)
土壤容重是土壤的干质量除以土壤的体积,用于计算延迟因子。 创建的任何新模型的默认土壤容重值为 1700 kg/m3。 如果在当前传输变量中未选择吸附方法,则模拟不需要体积密度值,并且左侧工具栏中的所有选项都将被禁用。 对于大多数传输引擎
(MT3DMS、RT3D、SEAWAT 和 SURFACT),为移动和非移动域的每个模型单元指定一个体积密度。 然而,对于 MODFLOW-6,可以根据适用情况为每个模型单元分别为移动域
(RHOB) 和每个固定域(RhoIm1、RhoIm2、RhoIm3)定义干堆积密度。 固定域的干堆积密度用于计算每个固定域中吸附成分衰减的质量平衡项。 有关更多信息,请参阅 Langevin 等人 (2022) 的 MODFLOW-6 地下水输送模型的技术文档。
RHOB 是土壤的干容重 [M/L3]。 创建的任何新模型的默认土壤容重值为 1700 kg/m3
RhobIm(或 RhoIm1)是第一固定域中土壤的干容重 [M/L3]
RhoIm2 是第二固定域中土壤的干容重 [M/L3]
RhoIm3 是第三固定域土壤的干容重 [M/L3]
Immobile Domain Porosity (PhiIm)
Immobile domain porosity is used to calculate the mass balance of constituents in the immobile domain.
PhiIm (or PhiIm1) is porosity of the (first) immobile domain as the volume of the (first) immobile domain pore space per total volume [L3/L3] → [-]
PhiIm2 is porosity of the second immobile domain as the volume of the second immobile domain pore space per total volume [L3/L3]→ [-]
PhiIm3 is porosity of the third immobile domain as the volume of the second immobile domain pore space per total volume [L3/L3] → [-]
PhiF is the fraction of the mobile domain relative to the total volume [L3/L3] → [-]
For more information, please see the property zones Theory section.
Reaction Model/Modules - Model Parameters
The reaction modules associated with RT3D may require model parameter values and settings to be selected at the Define Modeling Objectives workflow step, these include:
· Instantaneous aerobic degradation of BTEX [IREACT=1];
· Six-Species, First-Order, Rate-Limited, BTEX Degradation using Sequential Electron Acceptors [IREACT=3];
· Rate-Limited Sorption [IREACT=4];
· Double-Monod Model [IREACT=5];
· Sequential First-Order Decay [IREACT=6]; and
· Aerobic/Anaerobic PCE/TCE Dechlorination [IREACT=7]
where IREACT is a flag in the RT3D RCT package file whose value determines which reaction module is active in the transport model run.Note that the developers of RT3D did not implement a reaction model for IREACT=2 (Clement, 2003).
Please Note: yield values associated with several reaction modules in RT3D are defined on a mg/L basis. To be consistent the user must use mg/L units for all concentrations when using these reaction modules. See additional notes for each reaction module.