第12章. 超越基本建模概念

我迄今为止尚未见到任何问题，无论多么复杂，只要你以正确的方式看待，它就不会变得更加复杂。

—— 保尔·安德森

哦，那条路通向疯狂；让我远离那条路。

——《李尔王》（第三幕，第四场）

12.1 介绍

地下水流模型，如第1-10章所述，可以解决许多实际问题；然而，地下水模拟人员最终将遇到涉及超出简单地下水流的过程的问题。在第3.1节中，我们简要提到了不饱和流、可变密度流和多相流。模拟这三个过程中的任何一个都需要比基本地下水流方程（Eqn (3.12)）更复杂的流动方程。模拟溶质和热量传输以及对地下水进行严格的模拟表面水路需要将地下水流模型与传输模型或降雨径流模型链接或耦合，分别。此外，复杂的地下过程，如管道和断裂的流动以及含水层压缩，通常不包括在基本的地下水流模拟中，尽管一些流动代码具有针对这些过程的选项或附加包。此外，还有用于模拟具有随机非均质性含水层特性的地下水流和传输的代码，以及专为在管理问题中优化决策而设计的代码。

本书的第一版（Anderson和Woessner，1992）包括一个最终章节，简要总结了模拟复杂地下水问题的理论和方法。然而，现在有许多模拟地下水模型中的复杂过程的方法。例如，在地下水流代码中有不饱和流和将地下水和河流流量模型耦合的选项。现代编程语言，如Python，非常适合处理地下水问题（例如，Bakker，2014），可以扩展建模工具以包含其他科学的前沿和经过良好测试的方法。还有许多其他教材、报告和期刊论文，可帮助从业者学习如何模拟地下水是重要但不是唯一感兴趣的过程的复杂问题。鉴于自第一版以来信息的爆炸性增长，关于高级主题的单一全面章节是不可能的。相反，本章旨在为那些希望超越第1-10章中涵盖的地下水流建模的基本概念的人提供起点。在本章中，我们：（1）将读者引回我们的书中简要讨论高级过程模拟的章节；（2）提供其他一些思考和一些参考资料，以帮助读者开始模拟地下水系统中的复杂过程。有关本书末尾问题的Python版本等其他参考资料和资源，请访问http://appliedgwmodeling.elsevier.com。

首先，我们涵盖了复杂的地下流动过程：通过管道和断裂的流动；含水层压缩；不饱和流动；和可变密度流。然后，我们讨论地下水中热量和溶质的传输，这需要将地下水流模型与一个或多个传输模型链接或耦合。接下来，我们讨论地表水过程；地下水代码内的简单模拟地表水可以通过选项和包来处理，但严格的地表水建模需要将地下水代码与地表水代码耦合。最后，我们讨论随机地下水建模和决策支持和优化框架中的地下水建模。

在进行复杂过程的建模之前，模型制定者首先应考虑是否有必要模拟地下水流之外的其他过程。模拟超出地下水流的附加过程是否合理？适当的模型（第9.2节）更简单、成本更低吗？模拟复杂过程需要模型制定者更多的知识和技能，以及额外的参数和这些参数的适当值。在没有信息来约束新参数的情况下，添加额外的参数会在模型校准期间（图10.2）增加成本。否则，一个最初不恰当和不确定的地下水问题（第9.2节）变成了一个更加复杂的问题，但仍然是基本上不适当和不确定的问题。此外，模拟复杂过程的复杂性增加了计算负担；解决方案的潜在数值不稳定性；要分析的输出量；以及结果的不确定性。美国国家研究委员会（NRC，1990，第14页）于1990年得出的结论仍然相关：地下水模型的复杂性在变化。用于分析特定场地的模型的复杂性应由分析的问题类型确定。虽然更复杂的模型增加了可以描述的情况范围，但增加的复杂性需要更多的输入数据，需要模型制定者更高水平和范围的技能，并且如果输入数据不可用或质量不足以指定模型的参数，则可能引入输出的更大不确定性。

然而，即使是由观测值较差约束的复杂模型在价值上也可能具有价值，因为只有探索复杂过程对预测及其相关不确定性的影响才能有价值。使用第9章和第10章的工具，例如参数可辨识性（图10.10）或在Box 9.6中描述的更先进的工具，模型制定者可以识别

并专注于对预测最重要的因素，即主导预测不确定性的因素。

在整本书中，我们使用有限差分（FD）代码MODFLOW和有限元（FE）代码FEFLOW来说明地下水流建模概念。

幸运的是，MODFLOW和FEFLOW都可以通过在代码中打开高级选项、添加模块或将其链接或耦合到其他代码来处理许多复杂的过程。在下面的讨论中，我们提到了这些代码中的特定高级选项，但建议读者咨询MODFLOW和FEFLOW的开发人员维护的网站，以获取功能和相关信息的最新列表。

虽然我们用MODFLOW和FEFLOW进行说明，但本书中包括本章在内的概念对于任何地下水建模代码都是相关的。MODFLOW和FEFLOW的摘要描述，以及由Bear和Cheng（2010，第583-591页）提供的其他一些代码，包括流行的Transport Of Unsaturated Groundwater and Heat（TOUGH）代码的摘要描述，Barnett等人（2012年；表4-1，第42-44页）。 TOUGH是一个多过程控制体积FD代码，模拟“多孔和断裂介质中水汽、蒸汽、非冷凝气体和热量的耦合传输”（劳伦斯伯克利国家实验室：http://esd.lbl.gov/research/projects/ tough/）。 TOUGH旨在模拟“渗透性介质中耦合热、水文、地球化学和机械过程” ，适用于“核废物处置、环境修复问题、地热能生产、油气储层以及气体水合物沉积的地质储层，地质碳封存、非饱和带水文。” TOUGH目前对于表示其他潜在重要因素（如地表水过程）的能力有限；用于模拟化学反应的选项也仅限于与其主要应用有关的模拟地热储层和高放射性核废物地质存储库周围流体的化学反应。有关TOUGH的更多信息，请参见Finsterle等人（2014）的评论论文以及TOUGH网站。

12.2 复杂的地下水流过程

12.2.1 裂缝和管道中的流动

当裂缝均匀分布且相互连接时，可以将流动系统适当地模拟为等效多孔介质（EPM）（Box 2.2），使用本书第1-10章中描述的方法。有效的含水层参数（对于水力导数、储存和孔隙度）被定义为使模型能够再现现场观察到的流动模式和速率。EPM方法假定裂缝多孔材料可以被视为连续体，并且可以定义具有有效含水层参数的代表性基元体积（REV；第3.2节）。

EPM方法的应用困难在于确定定义等效水力特性所需的REV的适当大小（例如，Muldoon和Bradbury，2005年）。EPM方法可能足以表示区域流系统