SWAT模型与MODFLOW模型的耦合计算及应用.pdf

《SWAT模型与MODFLOW模型的耦合计算及应用.pdf》由会员分享，可在线阅读，更多相关《SWAT模型与MODFLOW模型的耦合计算及应用.pdf（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第26卷第3期2010年5月水资源保护WATER RESOURCES PROTECTIONVol.26 No.3May 2010 基金项目:河海大学理科基金项目作者简介:张雪刚(1971),男,甘肃会宁人,讲师,主要从事水文水资源方面的研究。E2mail:zhangxg SWAT模型与MODFLOW模型的耦合计算及应用张雪刚1,毛媛媛2,董家瑞1,李致家1(1.河海大学水文水资源学院,江苏 南京 210098;2.江苏省水利工程规划办公室,江苏 南京 210029)摘要:为了准确模拟区域地下水动态变化,将地表水SWAT模型与地下水MODFLOW模型进行耦合计算。根据SWAT分布式模型的特点,以

2、其水文响应单元(HRU)和MODFLOW模型的有限差分网格(cell)作为基本交换单元,将SWAT模型计算的地下水补给量和潜水蒸发量引入MODFLOW模型的地下水补给(RCH)模块和潜水蒸发(EVT)模块中,并将其应用于徐州市张集地区的地下水模拟计算。结果表明,SWAT模型与MODFLOW模型的耦合计算能准确模拟和预测该地区的地下水水情及其地表水和地下水之间的相互作用。关键词:地表水;地下水;SWAT模型;MODFLOW模型;模型耦合;徐州市中图分类号:P339 文献标识码:A 文章编号:10046933(2010)03004904A coupled model simulation and

3、application of SWAT2MODFLOWZHANG Xue2gang1,MAO Yuan2yuan2,DONG Jia2rui1,LI Zhi2jia1(1.College of Hydrology and Water Resources,Hohai University,Nanjing210098,China;2.Jiangsu Water ResourcesPlanning Office,Nanjing210029,China)Abstract:In order to simulate change in regional groundwater,the quasi2dist

4、ributed watershed model,SWAT wascoupled with the fully2distributed groundwater model MODFLOW,in simulation.The hydrologic response units(HRU)in the SWAT model,its basic elements,were exchanged with cells in the MODFLOW model.By using the HRUcellconversion interface,the distributed groundwater rechar

5、ge rate and the groundwater evapotranspiration calculated bySWATwere introduced into the RCH and EVT modules of MODFLOW.The application of groundwater simulation in theZhangji area of Xuzhou demonstrated that the coupled SWAT2MODFLOW was capable of predicting the groundwatersituation,and of being us

6、ed for investigating the interaction between surface water and groundwater.Key words:surface water;groundwater;SWAT model;MODFLOW model;model coupled;Xuzhou City 地表水和地下水之间存在密切的联系和相互作用,地表水和地下水之间的相互作用规律一直是水文学和水文地质学研究的热点和难点。目前,国内外有关地表水和地下水模拟的模型很多,如以模拟地表径流为主的TOPMODEL1、SWAT模型23、新安江 模 型 等4,以 模 拟 饱 和 地 下 水

7、 流 为 主 的MODFLOW5、FEFLOW6等,在不同区域的水资源管理和规划中得到了广泛应用。在地表水和地下水的联合模拟方面也有许多相关研究,开发的模型有SWATMOD78、IGSM9、MODHMS10、MIKE2SHE11模型等。这些耦合与集成模型功能完善,但是由于结构复杂,同时需要大量参数和数据支撑,模型建立和率定非常耗时,在应用上受到一定限制。因此根据实际情况,研究简单、实用的地表水与地下水模型耦合计算的方法,准确、有效地模拟地表水与地下水系统之间的相互作用,对于水资源的准确评价,科学开发利用、有效管理以及生态环境的保护具有重要意义。张集地区位于江苏省徐州市,区内地下水资源丰富,建有

8、日供水量10万t的地下水水源地,是向徐州市提供饮用水的水源地之一。由于地区内地质条件特点,地表水渗漏作用强烈,对地下水存在很强94的补给作用。准确模拟和预测地区地下水的动态变化,对保障徐州市供水至关重要,也是实现区域水资源有效保护和可持续利用的关键。在地下水模型的模拟计算中,地下水的补给量和潜水蒸发量直接影响地下水动态模拟的结果,但是通常很难对其进行准确估算和赋值。笔者根据张集地区地下水模拟的要求,在地区SWAT模型模拟和MODFLOW模型建立的基础上,根据SWAT分布式模型的特点,通过SWAT模型水文响应单元(HRU)与MODFLOW差分网格(cell)的对应,构建SWAT2MODFLOW的

9、交互界面。应用ARCGIS软件,将SWAT模型地下水补给量与潜水蒸发量的模拟计算结果及空间分布引入MODFLOW相应的地下水补给子模块RCH和潜水蒸发子模块EVT中,实现模型的耦合计算。对张集地区地下水水位进行模拟,并与实测结果进行对比分析。1 模型方法1.1SWAT模型SWAT(soil and water assessment tool)模型是美国农业部(USDA)农业研究所(ARS)开发的适用于大流域尺度的分布式非点源污染计算模型,用于模拟预测在具有多种土壤类型、土地利用和管理条件的大尺度流域里,土地管理措施对水、沙和化学物质的长期影响。SWAT是基于物理过程的模拟模型,模型主要输入要素

10、包括气象、水文、土壤特性、植被覆盖、化学物质和土地利用情况等。SWAT模型计算效率较高,可以模拟和研究长期的变化过程。在国外水资源管理、水土流失等方面得到了成功应用,在国内的研究和应用主要侧重于水文模拟方面。SW AT模型是具有很强物理机制的半分布式流域水文模型,在SW AT模型中,流域被划分为多个子流域,在每个子流域内,根据不同的土地利用和土壤类型,再将每个子流域划分为1个或多个下垫面相对单一的水文响应单元(HRU),作为模型的基本计算单元。每个HRU的计算结果在子流域范围内进行累加,并演算到支流,最后通过河道汇流计算得到流域出口处的结果。笔者应用A VSW AT 2000进行计算。1.2M

11、ODFLOWMODFLOW(the modular finite2difference groundwaterflow model)是由美国地质调查局(USGS)开发的、用来模拟地下水流动和地下水中污染物迁移等特性的计算机程序,是基于达西定律和地下水质量平衡的具有物理意义的三维地下水模拟模型,使用有限差分方法,采用模块化的程序结构。Visual MODFLOW12是由加拿大waterloo水文地质公司在MODFLOW软件的基础上,综合已有的MODFLOW、MODPATH、MT3D、RT3D和WinPEST等地下水模型软件开发的,可进行三维水流模拟和溶质运移模拟。Visual MODFLOW以其

12、简单实用的求解方法、广泛的适用范围及强大的可视化功能等特点,成为目前国际上应用最广泛的三维地下水流和溶质运移模拟评价的标准可视化专业软件系统。笔者应用Visual MODFLOW410进行地下水模拟计算。1.3SWAT2MODFLOW耦合方法在SWAT模型中,水文响应单元(HRU)是模型最基本的计算单元,反映了子流域内土地利用和土壤覆盖情况的综合影响,是分布式模拟计算的基础。通过模拟,可以得到每个HRU内的地表径流、地下径流、壤中流、蒸发、渗漏量、地下水的补给量等计算结果,将HRU计算结果在子流域和整个流域进行累加,就得到计算模型要素在整个流域的空间分布。MODFLOW的基本计算单元是有限差分

13、网格(cell)。根据SW AT中HRU和MODFLOW中有限差分网格的特点,构建HRU2cell的交互界面13(图1)。图1中上半部分显示了SW AT模型1个子流域中HRU的空间分布,下半部分显示了MODFLOW中有限差分网格(cell)的划分情况,网格中的数字分别对应SW AT模型中HRU的编号。通过这种对应,可以将SW AT模型计算的模型要素的空间分布赋值给MODFLOW,从而实现SW AT与MODFLOW的耦合计算。图1SWAT模型与MODFLOW的HRU2cell耦合示意图2SWAT2MODFLOW耦合计算在张集地区的应用2.1 研究区域概况张集地区位于江苏省徐州市东南,距徐州市中心

14、25 km,总面积360 km2。地区地下水储量丰富。地区总体地形西北高、东南低,黄泛冲击平原是区内主要地貌形态,在水源地北部及南部边界地区有碳酸岩盐组成的低山丘陵区,山峰海拔高60237 m,其余地区为平原区,高程2838 m。张集地区数字高程模型DEM及其地理位置见图2。张集地区的地下水主要有孔隙水含水层和岩溶05图2 张集地区DEM及其在徐州市的地理位置示意图水含水层两种。孔隙含水层广泛分布在平原区,组成以砂土、粉土、黏性土为主,厚度一般不足40 m。水源地的岩溶含水层以震旦系和寒武系的碳酸盐岩为主,岩性多变,厚度大、范围广,是水源地地下水的主要开采层区。2.2SWAT模型率定及模拟计算

15、张集地区面积较小,且缺少实测地表水径流站,因此选择与张集地区相邻的具有相似条件的栏杆集地区,采用栏杆集水文站19801987年8年的月实测流量对模型参数进行率定。栏杆集地区位于张集水源地西部,流域面积600km2,与张集地区同属淮河流域,地区内的气候条件、土地利用及土壤覆盖情况与张集地区基本相同。模型率定的参数及地表径流量的拟合情况见表1和图3。表1SWAT模型率定的主要参数值名称含义数值名称含义数值CN2径流曲线系数-7GW_REVAP地下水蒸发系数011ESCO土壤蒸发补偿系数018GWQMN回归流产生所需地下水深阀值100mmALPHA_BF基流因子017REVAPMN浅层含水层中再蒸发

16、或下渗的地下水深度阀值2mmSOL_AWC土壤有效含水率01035引用Nash Sutcliffe模型效率系数Ens和线性回归系数R2评估模型的模拟效果。模型率定结果的确定性系数为0185,相关系数为0186。将率定的SWAT模型应用于张集地区。根据地区土地利用和土壤类型,共划分12个子流域。由于张集区域面积较小,且土地利用以农田为主,土地开发利用方式及土壤覆盖情况比较单一,因此HRU划分采用优势土被占有法,每个子流域对应1个HRU,共分为12个HRU(图3)。选用地区2005年的实测降雨和蒸发资料,进行SWAT月模拟计算,得到区域内每个HRU的地下水补给量和潜水蒸发量值。2.3SWAT与MO

17、DFLOW的耦合计算根据张集地区的水文地质条件,建立区域的MODFLOW地下水模型。模型分为3层,第1层为非承压的孔隙含水层,第2层和第3层为承压岩溶含水层。模型网格划分为每层46行 60列,每个单元格(cell)的大小为500m500m。对每个含水层,设置其渗透系数、地下水贮水率、给水度等参数和模型的边界条件14。应用ARCGIS软件,将SW AT模型计算的地下水补给量和潜水蒸发量结果引入MODFLOW模型边界条件中的地下水补给子模块RCH和潜水蒸发子模块EVT中,进行地下水模拟计算。图4为模拟得到的区域内35口地下水观测井2005年4月5日的实测地下水位与模型计算地下水位的散点图,模拟结果

18、的确定性系数为0178,模拟结果与实测值的相关系数R2为0181。图5显示了位于第7个HRU内的孔隙含水层一观测井实测与计算地下水水位过程线,计算得到的地下水水位随地下水补给量的变化趋势与实测地下水位的变化基本吻合。图3SWAT模型率定的栏杆集月径流量与实测流量对比图4 张集地区2005年4月5日观测井计算水位与实测水位散点拟合图3 结 论本文将地表水SWAT模型与地下水MODFLOW模型进行耦合并应用于徐州张集地区的地下水水位模拟计算。根据张集地区地下水模拟的要求,在地区SWAT模型模拟和MODFLOW模型建立的基础上,根据SWAT分布式模型的特点,将其水文响应单元(HRU)和MODFLOW

19、的有限差分网格(cell)相对15图5 张集地区孔隙含水层观测井2005年地下水位过程线应,应用ARCGIS软件,将SWAT模型计算的具有空间分布特征的模型要素地下水补给量和潜水蒸发量引入MODFLOW的地下水补给模块(RCH)和蒸发模块(EVT)中,实现模型的耦合计算。对区域2005年地下水模拟计算结果表明,计算得到的区域地下水位与实测结果基本吻合。SWAT模型与MODFLOW的耦合计算对于更准确地模拟和预测地下水水情,进一步研究地表水与地下水之间的相互作用提供了科学方法,对于区域水资源的科学规划和管理,实现区域水资源的可持续利用具有实际意义。参考文献:1 BEVEN K.T opmodel

20、:a critiqueJ.Hydrological Processes,1997,11:10691085.2张银辉.SWAT模型及其研究进展J.地理科学进展,2005(5):121130.3王中根,刘昌明.SWAT模型的原理、结构及应用研究J.地理科学进展,2003(1):2124.4赵人俊.流域水文模拟:新安江模型与陕北模型M.北京:中国水利水电出版社,1984.5 MCDONALD M G,HARBAUGH A W.A modular three2dimensional finite differnce ground2water flow model R.Washington:U S Ge

21、ological Survey Open2File Report,1988.6 DIERRSCH HJ G,FEFLOW users manualM.Berlin:Institutefor Water Resources Planning and System Research Ltd,1998.7 JOBSON H E,HARBAUGH A W.Modifications to the difusionanalogy surface2water flow model(DAFLOW)for coupling tothemodularfinitediferenceground2waterflow

22、model(MODFLOW)R.Washington:U S Geological Survey OpenFile Report,1999.8 SOPHOE LEOUS M A,K OE LIKERJ K,G OVINDARAJU R S,eta1.Integratednumericalmodelingforbasin2widewatermanagement:the case of the Rattlesnake Creek basin insouth central KansasJ.Journal of Hydrology,1999,214(14):179196.9 LABOLLE E M,

23、AHMED A A,FOGG G E.Review of theintegrated groundwater and surface2water model(IGSM)J.Groundwater,2003,41(2):238246.10 PANDAY S,HUY AK ORN P S.A fully coupled physically2based spatially2distributed model for evaluating surface flowJ.Advances in Water Resources,2004,27(4):361382.11贾仰文,王浩,倪广恒,等.分布式流域水

24、文模型原理与实践M.北京:中国水利水电出版社,2005.12 Anon.Visual MODFLOW V.21812 users manual forprofessional application in three2dimensional groundwater flowand contaminant transport modeling M.Ontario:WaterlooHydro Geologic Inc,2000.13 NAM W K,MOON C,Y OO S W,et al.Development andapplication of the integrated SWAT2MOD

25、FLOW model J.Journal of Hydrology,2008,(356):116.14江苏省水文地质工程地质勘察院.江苏省徐州市张集水源地供水水文地质详查报告R.南京:江苏省水文地质工程地质勘察院,2001.(收稿日期:20091022 编辑:高渭文)(上接第44页)9WANGYuan2yuan,CHENFei2zhou.Decompositionandphosphorus release from fourdifferentsize fractions ofmicrocytis spp.take from Lake Taihu,China J.Journal ofEnviro

26、nmental Science,2007,7:892899.10谢平,陈宜瑜.加强淡水生态系统中生物多样性的研究与保护C/中国环境与发展国际合作委员会.保护中国的生物多样性.北京:中国环境科学出版社,1996:208217.11丁笑生,段红英.微囊藻毒素及其生殖毒性的研究进展J.水产科学,2007,26(1):6162.12左金龙,崔福义,刘智晓,等.饮用水中蓝藻毒素污染研究进展J.环境污染治理技术与设备,2006,7(3):813.(收稿日期:20081102 编辑:高渭文)勘 误2010年第2期第76页图1(a)、(b)横坐标标值有误,正确标值见下图。特此更正。图1 进样温度对待测物响应值的影响25