物理审计与不确定性感知的反应性硝酸盐输运替代建模

《Science of The Total Environment》：Physics audited and uncertainty aware surrogate modeling for reactive nitrate transport in groundwater

　　

地下水中硝酸盐（NO₃^-）污染常因电子供体限制、氧化还原竞争和长驻留时间而持续存在，这些因素制约了自然衰减。本研究提出一个物理审计与不确定性感知的替代建模框架，用于多孔介质中反应性硝酸盐输运。研究人员生成了四个一维基准的高保真度PHREEQC模拟，其地球化学复杂度递增，从线性异养反硝化到具有显式NO₃^-–Fe(III)竞争的双底物Monod动力学。使用CatBoost梯度提升算法直接在空间–时间浓度场上训练替代模型，通过贝叶斯优化调优超参数，并采用模拟级别数据划分以防止信息泄漏。替代模型接受水动力学和地球化学输入以及空间和时间坐标，预测整个域内的硝酸盐浓度场。采用子采样集成来表征认知变异性，并使用共形分位数回归提供校准后的预测区间，从而实现不确定性量化。应用事后物理审计来验证与主导的对流–弥散–反应平衡、溶质质量守恒和非负性的一致性。结果表明，替代模型准确再现了非线性反应锋面和供体–受体竞争，测试R²高达0.997，质量平衡误差通常低于0.5%。认知不确定性相对于偶然变异性较小，后者集中于反应过渡带。该框架是一个可靠且计算高效的工具，适用于硝酸盐污染含水层的情景分析和基于风险的地下水质量管理。

**论文解读文章**

**研究背景**
地下水硝酸盐（NO₃^-）污染对饮用水安全和水生生态系统完整性构成广泛威胁。尽管地表氮输入有所减少，但在许多流域中，地下水硝酸盐浓度并未相应降低，这反映了水文滞后效应，如长驻留时间、包气带和低渗透弱透水层中的储存以及从次级储层的缓慢释放。此外，硝酸盐衰减受到生物地球化学条件的强烈限制：在氧化至亚氧化含水层中，反硝化作用取决于合适电子供体的可利用性，而硝酸盐去除往往反映出含水层基质有限的电子供体容量，而非内在反应动力学。与其他电子受体（如Fe(III)）的竞争进一步分割了供体分配，可能限制硝酸盐还原。这些耦合的水文和生物地球化学控制强调了在预测地下水硝酸盐持久性时，必须显式表征供体限制和氧化还原竞争的必要性。

**研究动机**
传统的反应性输运模型（RTM：Reactive Transport Model）能够耦合对流–弥散与平衡和动力学反应，但因其计算成本高昂，难以进行大规模集合模拟、情景分析和不确定性传播。数据驱动的替代模型可以以更低成本近似RTM的输入–输出行为，但目前大多数应用于硝酸盐输运的替代模型采用化学上简化的公式（如一阶衰减），未探索多底物动力学或显式氧化还原竞争。本研究旨在开发一个结合物理审计和不确定性感知的替代建模框架，在高保真度RTM模拟基础上，高效模拟反应性硝酸盐输运，同时保留机理水地球化学保真度。

**研究内容与意义**
研究人员构建了四个一维反应输运基准，涵盖从线性异养反硝化到具有显式NO₃^-–Fe(III)竞争的双底物Monod动力学，使用PHREEQC通过拉丁超立方采样（LHS：Latin Hypercube Sampling）生成模拟。采用CatBoost梯度提升算法在完整的空间–时间浓度场上训练替代模型，并结合贝叶斯优化调优超参数。通过子采样集成和共形分位数回归分解预测不确定性为认知（epistemic）和偶然（aleatoric）成分。事后应用物理审计验证对流–弥散–反应平衡、质量守恒和非负性。结果表明，替代模型准确复制了非线性反应锋面和供体–受体竞争，测试R²高达0.997，质量平衡误差低于0.5%，认知不确定性较小。该框架为硝酸盐污染含水层的情景分析和基于风险的地下水质量管理提供了可靠且计算高效的工具。论文发表在《Science of The Total Environment》。

**关键技术方法**（不超过250字）
本研究采用的关键技术方法包括：(1) 使用PHREEQC（v3.8.6-17100）结合phreeqc.dat热力学数据库，通过拉丁超立方采样生成水动力学和地球化学参数空间填充设计，构建四个一维反应输运基准模拟；(2) 采用CatBoost梯度提升算法，在空间–时间浓度场上训练替代模型，并通过贝叶斯优化调整超参数；(3) 不确定性量化：通过基于模拟标识符的子采样（80%训练模拟，无放回）集成估计认知方差，同时使用共形分位数回归（0.1和0.9分位数）经留一模拟外出（LOSO：Leave-One-Simulation-Out）校准后得到标准化预测区间，区分认知与偶然变异性；(4) 事后物理审计：通过有限差分计算局部对流–弥散–反应（ADR）残差、全局质量平衡误差、非负性和边界残差指数（BRI：Boundary Residual Index）验证物理一致性。无样本队列来源，全部数据来自数值模拟。

**研究结果**

**3.1 Benchmark 1: 线性硝酸盐反硝化**
通过比较CatBoost、LightGBM和随机森林三种集成算法，CatBoost取得最低测试误差（RMSE=6.1×10^-5 mol kgw^-1，R²=0.95），优于其他算法。物理审计表明ADR残差和质量平衡误差均较小，无负浓度，替代模型准确再现了指数型硝酸盐衰减行为。

**3.2 Benchmark 2: 双线性氧化还原竞争**
在引入Fe(III)作为竞争电子受体后，CatBoost保持强预测能力，测试R²=0.9385，残差和误差仍较小。替代模型精确捕捉了由供体分配锋面后的硝酸盐衰减前沿，残差局限于反应过渡带。

**3.3 Benchmark 3: 双底物Monod反硝化**
采用非线性双底物Monod动力学时，CatBoost测试R²=0.9472，RMSE=6.9×10^-5 mol kgw^-1。物理审计确认ADR平衡满足，虽出现少量负预测（极小），但不影响全局质量平衡。替代模型再现了指数–S形衰减特征。

**3.4 Benchmark 4: 双Monod动力学与氧化还原竞争**
在最复杂的基准情景（双底物Monod+Fe(III)竞争）中，CatBoost表现最优，测试R²=0.9965，RMSE=1.7×10^-5 mol kgw^-1。物理审计显示残差、质量平衡误差和边界残差均处于低水平。精度提升归因于供体共享和氧化还原竞争使模拟数据流形更结构化，学习更容易。

**3.5 比较分析**
通过分析所有基准中最大硝酸盐浓度梯度的分布，发现锋面陡度范围类似，说明动力学对最大锋面陡度影响较小，锋面位置和推进是主要差异来源。替代模型残差主要集中于陡梯度区域，而非系统性偏差。

**3.6 不确定性量化**
针对Benchmark 4的量化表明：认知方差高度右偏且集中于零，而偶然方差分布更广且量级更大，表明总预测不确定性的主导来源是偶然变异性（源于子网格过程或固有异质性）。共形化80%预测区间（PI80）宽度沿反应衰减锋面升高，其他区域较小，经验覆盖率接近名义水平。

**3.7 综合、启示与局限**
综合四个基准，CatBoost替代模型在不同复杂度下均能高保真模拟一维反应性硝酸盐输运。物理审计验证了物理一致性。不确定性量化将认知与偶然成分分离，表明替代模型不确定性小于数据固有变异性。局限包括：所有基准均为一维均匀介质模拟，未考虑真实含水层中的异质性、优先流和多维混合；少量负预测可通过后处理修正；进一步扩展需考虑二维/三维域、时变边界条件、氧气抑制等。

**总结讨论与结论翻译**
讨论部分强调，该物理审计和不确定性感知的替代框架对地下水质量评估和管理具有重要意义。训练后的替代模型评估硝酸盐浓度场及相关不确定性的速度比底层PHREEQC模型快3–5个数量级，支持原本计算不可行的正向不确定性传播、情景筛选和敏感性分析，能快速评估不同负荷、供体可用性或修复情景下的硝酸盐持久性，并支持在不确定地下条件下的概率监管超标风险筛查。基准结构、物理审计和不确定性量化协议具有通用性，可转移至其他供体限制和多电子受体氧化还原系统，为开发物理可解释的风险决策替代模型提供模板。

**研究结论**（翻译自4. Conclusions）：
• 在完整空间–时间场上训练的CatBoost替代模型高保真地再现了硝酸盐浓度（Benchmark 4：测试R²=0.9965，RMSE≈1.7×10^-5 mol kgw^-1（约1 mg L^-1以NO₃^-计））。
• 误差局限于狭窄的反应过渡带；残差结构反映小的锋面位置和梯度失配，而非系统性输运或反应偏差。
• 事后物理审计验证了对流–弥散–反应平衡和全局质量守恒，域内残差和质量平衡误差<0.5%，边界残差指数<0.05。
• 不确定性量化协议区分了认知与偶然成分；校准后的80%预测区间（PI80）达到了名义覆盖率，总预测方差由集中于衰减锋面的偶然变异性主导。这些不确定性估计基于底层模拟，不捕捉模型结构误差或真实异质性的额外不确定性。
• 利用训练后的替代模型，大规模正向模拟和不确定性量化分析运行速度比直接RTM评估快多个数量级，使原本计算上不可行的参数探索和不确定性研究成为可能。
• 基准结构和工作流程可转移至其他供体限制或多电子受体氧化还原系统（如硫酸盐、锰或铁驱动的衰减），为开发将反应性输运过程与地下水质量和生态系统保护相联系的环境可解释替代模型提供了通用模板。