SWATMOD模型 耦合原理

3.几个典型耦合模型分析

3.1 “四水”转换模型 

“四水”转化是从水资源计算与评价的角度提出的，涉及大气水、地表水、土壤水与地下水相互转化关系问题。 “四水”转化模型把计算单元或流域（区域）看成一个系统，把降水、蒸发、地表水、土壤水、地下水及农业生态和灌溉等视为相互关联、相互制约的子系统。这个模型的理论基础是流域水循环和平衡原理，它利用物理方程和逻辑判断，计算出各部分水量变化过程。 

“四水”转化模型中的各部分水量转化层次明确而有序，但是该模型主要针对水资源统一评价问题建立，地表、地下水模型之间的相互作用和转换非常简化，因此应用时有很大的局限性。

3.2 SWATMOD模型

SWATMOD模型耦合了美国农业部农业研究局开发的半分布式水文学模型SWAT和美国地质调查局开发的MODFLOW模型。它不仅能模拟分析灌溉、施肥和耕作措施对农业生产和水资源利用方面的影响，而且能反映复杂含水层系统中地下水的动态[4]。 

SWATMOD模型中覆被层、土壤层、池塘、水库的水分运动由SWAT模拟，如图1虚线以上的部分；地下含水层的水分运动、河流与地下水的相互作用由MODFLOW模拟，如图1虚线以下部分；穿过虚线的箭头表示两者之间的交换。Perkins为了解决SWAT和MODFLOW之间时间、空间尺度不一致以及耦合的其他需要，修改了SWAT源代码。他将各地下水计算时段内SWAT的日模拟结果进行累计，借助GIS软件将每个HRU与其所对应的地下水模型网格联系起来，进而可以将SWAT模拟的潜在蒸散发、地下水补给量、河流水位转换为MOD-FLOW需要的尺度。借助GIS软件将河流以及它所经过的HRU和地下含水层网格联系起来，采用MODFLOW中的河流子程序包计算出地下水对河流的补给或排泄量，加入SWAT的河道入流中，用蒸发蒸腾子程序包计算潜水蒸发传给HRU。最近，Kim也提出了另一种将SWAT与MODFLOW相耦合的方法。该耦合模型以日为时间尺度，将MOD-FLOW和SWAT都分为输入部分和计算部分两大块，然后用MODFLOW取代SWAT中的地下水模块，由HRU在各时段调用。

SWATMOD模型将从陆地水文学角度建立的概念性水文模型和从水文地质学角度建立的地下水动力学模型相结合，能够更充分的利用水文气象和水文地质资料，而且耦合模型可以在两类模型中取长补短。例如，SWAT本身含有对地下水的描述，但并不能较准确反映河流与含水层之间的相互关系以及地下水抽水井的分布。用MODFLOW取代SWAT地下水模块，就能很好地解决这些问题。SWAT也能为MODFLOW提供更加准确的蒸散发、入渗补给量等的空间分布信息。

3.3 MODBRNCH模型

MODBRANCH模型是Swain E.D.在1994年提出的，它耦合了一维明渠不稳定地表水模型BRANCH和三维地下水流模型MODFLOW。MODBRANCH模型在空间上将地表水系细分为子支流，将子支流再分为子子支流，假定渠系为内联的，且平面上子支流的尺寸不能大于一个地下水单元的尺寸，垂向上河流不能跨越多个地下水模拟层。模型采用BRANCH，通过求解含侧向交换量的一维明渠运动波方程和连续性方程来模拟地表水流的运动；采用MODFLOW，通过求解根据连续性方程和达西定律建立的地下水运动方程对地下水流进行模拟。在计算地表水和地下水的交换量（基流量）时，根据MODFLOW中一个地下水时间步长内河网结点的地下水头始末值，内插得到BRANCH模型各相应时间段的地下水头，再根据达西定律计算该时段地表水和地下水的交换量。地表水和地下水的交换量被概化为通过低渗透性的河床底积层的垂向水流运动。

MODBRANCH模型将一维明渠不稳定地表水模型BRANCH和MODFLOW相耦合，对河流-含水层相互作用的模拟功能比MODFLOW更完善。此外，MODBRANCH中地表水和地下水模拟的时间尺度可以不一致，地下水的时间尺度可以是地表水时间尺度的好几倍，这样可适当减少需要收集的数据量[1]。

3.4 MIKE-SHE模型

MIKE-SHE系列模型是在SHE模型的基础上发展起来的模型，是一个综合性的、确定性的且具有物理意义的分布式水文模型。它能够模拟陆地水循环中几乎所有主要的水文过程，包括水流运动、水质和泥沙输移，可以用于解决与地表水和地下水相关的资源和环境问题。MIKE-SHE模型包括水流运动、水质、土壤侵蚀、灌溉模块。其中水流运动模块MIKE-SHE WM（如图2）是模型的核心，坡面/河道流通过求解连续性方程和扩散波方程来近似模拟；土壤水采用应用一维Richards方程进行模拟；地下水运动采用三维地下水的偏微分方程来描述，源漏项包括与非饱和层的水量交换、与河道的水量交换、地下水抽取/补给、潜水蒸发、排水管道的流量。模型应用时，流域在平面上被划分为许多矩形网格，通过离散化的坡面流、河道流、土壤渗流和地下水运动偏微分方程来建立相邻网格之间的时空关系，以有限差分法进行求解。各子模块可采用不同的时间尺度，但需要根据连接过程的尺度进行转换匹配。

MIKE-SHE模型采取较严格的水动力学瞬变偏微分方程描述水文过程，具有很好的物理基础，但需要大量精确的参数和数据支撑，因此建立和率定模型非常耗时，而且需要很深的专业知识，使模型在推广中面临着很多困难。此外，因为对土壤水流采取了一维垂向流的简化，它不太适用于求解坡面坡度过大的水流模拟。

结语

地表水和地下水耦合模型是分析地表水与地下水相互作用应用较广的方法。国内外已有的耦合模型大多是在对现实特定条件概化基础上建立起来的，有些是把一些成熟的地表水和地下水模型连接起来，如SWAT和MODFLOW的耦合，但由于两类模型在模拟方法、参数定义以及计算单元的划分等关键问题上存在巨大的差异，并且在理论上也存在一定的缺陷，在实际应用中有待于进一步发展和完善。系列软件MIKE-SHE和MODHMS模型功能较为完善，但它们需要大量的参数和数据支撑，建立和率定模型非常耗时，理论上的严格同操作上的困难之间的差距成为耦合模型应用的一道障碍。因此，综合水循环各组成部分的复杂系统，研究更为复杂的实际应用问题，适应更为广泛的需求是地表水和地下水耦合模型发展的必然趋势。此外，耦合模型所面临界面时空尺度转换、非饱和带土壤水分运动机理、参数的不确定性、河流与地下水相互转化机理等问题都需要进一步深入的研究。

参考文献：

[1] 胡立堂, 王忠静, 赵建世. 地表水和地下水相互作用及集成模型研究. 水利学报, 2007, 38(1): 54~59.

[2] 张德同, 付艳红. 白城平原区“四水”转化产流模型的探讨. 东北水利水电, 2004, 7(22): 6~9.

[3] Perkins S P, Sophocleous M A. Development of a comprehensive watershed model applied to study stream yield under drought conditions. Ground Water, 1999, 37(3): 418~426.

[4] Ross Mark, Geurink Jeffrey, Said, Ahmed. Evapotranspiration Conceptualization in the HSPF-Modflow Integrated Models. Journal of the American Water Resources Association, 2005, 41(5): 1013~1025.

[5] Steven L. Markstrom, Richard G. Niswonger. GSFLOW — Coupled Ground-Water and Surface-Water Flow Model Based on the Integration of the Precipitation-Runoff Modeling System (PRMS) and the Modular Ground-Water Flow Model (MODFLOW-2005). Chapter 1 of Section D, Ground-Water/Surface-Water Book 6, Modeling Techniques. Available at http://pubs.usgs.gov/tm/tm6d1/pdf/tm6d1.pdf.

[6] 贾仰文, 王浩, 倪广恒等. 分布式流域水文模型原理与实践. 北京: 中国水利水电出版社, 2005, 10~13.

[7] 王蕾, 倪广恒, 胡和平. 沁河流域地表水与地下水转换的模拟. 清华大学学报(自然科学版), 2006, 46(12): 1978~1981.

[8] 李兰, 钟名军. 基于 GIS 的 LL-Ⅱ 分布式降雨径流模型的结构. 水电能源科学, 2003, 21(4): 35~38.

[9] Yu Z, F W. Schwartz. Application of and Integrated Basin-Scale Hydrologic Model to Simulate Surface-Water and Ground-Water Interactions. Journal of the American Water Resources Association, 1998, 34(2): 409~425.

[10] Hsin-chi Jerry Lin, Jing-Ru Cheng. Using the Parallel WASH123D Code to Simulate Overland-Subsurface Interactions. World Water and Environmental Resources Congress, 2004, Salt Lake City, Utah, USA.

[11] (Hydro GeoLogic Software), MODHMS. Available at http://www.modhms.com/software/modhms.html.

[12] (Hydro GeoLogic Software), IHSim. Available at http://www.modhms.com/software/IHSim.html.

[13] Swain E D. Implementation and use of direct-flow connections in a coupled groundwater and surface-water model. Ground Water, 1994, 32(1): 139~144.

[14] Jobson H E, Harbaugh A W. Modifications to the diffusion analogy surface-water flow model (DAFLOW) for coupling to the modular finite-difference ground-water flow model (MODFLOW). U.S. Geological Survey Open-File Report, 1999.

[15] (WASY Software), IFMMike. Available at http://www.wasy.de/english/produkte/feflow/ifmmike11.html.

[16] Labolle EM, Ahmed AA, Fogg GE. Review of the integrated groundwater and surface-water model (IGSM). Ground Water, 2003, 41(2): 238~246.