美国地下水模拟软件MODFLOW的开发.docVIP

一、简 介 从上世纪70年代初开始，地下水模型的开发和应用一直是美国地质调查局（USGS）水文研究的一个重要组成部分。在过去的35年间，USGS已经开发出了一系列计算机模型来模拟饱和的和非饱和的地下水流、溶质运移和化学反应。其中应用最广泛的是MODFLOW模型，它利用有限差分法模拟三维（3D）地下水流（GWF）（Harbaugh，2005）。尽管最初的设想是将MODFLOW作为一种地下水流模型，但其模块化结构为结合其它模拟性能提供了一个稳健的结构，这无疑拓展了MODFLOW最初的设计理念。目前，与MODFLOW相关的模型系列能够模拟耦合的地下水和地表水系统、溶质运移、非饱和流、参数估计和地下水管理（GWM）（优化模型）。本文旨在总结当前USGS在MODFLOW开发方面的一些指令。在过去的两年间，几个与MODFLOW相关的计划的出版为本次总结提供了一个良好的契机。 USGS把模型开发作为自己职责的一部分，为美国水资源的有效管理提供了公正准确的水文资料和科学分析。比较有代表性的是，为解决水资源评价计划中不断涌现出来的一些特定问题，已经为MODFLOW开发出了一些新性能，如需要在非饱和带厚度很大的地区加强地下水与地表水相互关系的模拟或者沿海地区海水入侵的模拟。USGS的MODFLOW开发工作分布在几个机构组织中，通过地下水部（Office of Ground Water）进行整体协调、优先项目确定以及为代码开发提供支持。代码开发的这种分布式的形式，很大程度上是MODFLOW最初设计理念的一个副产物，它为模型新功能的增加提供了一种容易改进和完善的程序结构。我们认为，USGS为开发和测试计算模型提供了一个独特的环境，因为它是一个国家机构，且研究和应用程序之间也能得到最可能多的反馈，数据收集和资源评价计划也具有时间连续性。USGS的科学家们能够与调查局之外的同事们一起工作，在过去的几年间，对MODFLOW新功能的开发而言，这种合作一直是非常重要的。 二、最近的开发 MODFLOW的最新版本是MODFLOW-2005（Harbaugh，2005）。这一版本只包含了GWF和观测程序（Observation Processes），但是预计它的应用范围会不断拓展。与MODFLOW-2000（Harbaugh等，2000）相比，MODFLOW-2005的主要变化是使用了Fortran模块来声明可以为子程序共享的数据。尽管这些模块不会影响MODFLOW的一般用户，但是它们提供了一种在一次MODFLOW模拟中定义多重模型网格的手段，这是Mehl和Hill（2005）开发的新的局部网格加密（LGR）功能的基础。LGR允许用户在网格更大更粗糙的母模型中模拟更高精度的局部网格（称为子模型）中的水流。在短距离内的水力梯度发生明显变化的条件下，经常需要进行局部加密，以对含水层的非均质性进行更详细的刻画或者对污染物运移进行场地尺度的模拟。LGR利用两步迭代求解方法来耦合母模型和子模型，这样两者共同边界上的水头和流量是平衡的（Mehl和Hill，2005）。尽管目前的MODFLOW-2005版本没有直接支持形式参数估计和模型灵敏度分析，但是UCODE的新版本为MODFLOW-2005（和其它模型）提供了模型自动校正功能（UCODE_2005；Poeter等，2005）。可以应用JUPITE应用程序接口（JUPITER Application Programming Interface）使用其它一些程序进行灵敏度分析、模型校正、不确定性评价和优化（Banta等，2006）。 过去的两年间已经开发出几个新的MODFLOW软件包和程序，拓宽了MODFLOW的应用范围。由于篇幅限制，我们在此只介绍一些最新的进展，即与非饱和流、GWM、变密度流和溶质运移相关的那些开发工作。 最近开发的两个新软件包和一个新程序使MODFLOW的非饱和带水流模拟能力得到了显著提升。这些新功能为MODFLOW刻画地表与地下的水文联系（如土壤带的入渗，地下水补给和地表水/地下水的关系）提供了一些改进方法。Streamflow-Routing（SFR2）和Unsaturated-Zone（UZF1）这两个新软件包利用三维Richards流动方程的运动波近似，模拟非饱和带的一维（1D）垂向流动和储水。两个软件包都基于这样的假设，即非饱和流只受重力势能梯度的作用，而忽略负的势能梯度。SFR2软件包（Niswonger和Prudic，2005）模拟的是河流下面的非饱和流，其中河流与下部的含水层没有水力联系。UZF1软件包（Niswonger等，2006）模拟的是地表入渗条件下，非饱和带中水的流动、储存和蒸发。与SFR2和UZF1软件包相反，变饱和流（Variably Saturated Flow，VSF