供水格局变化下海河平原区浅层地下水动态响应剖析.docVIP

供水格局变化下海河平原区浅层地下水动态响应剖析 　　【摘 要】南水北调中线工程通水后，海河平原区因水源置换与地下水压采，供水格局发生转变。基于水资源转化动态模拟模型MODCYCLE，在对2001～2010年现状地下水动态平衡模拟分析的基础上，设置不同供水方案情景，量化模拟未来浅层地下水的动态响应。结果表明：供水格局变化后，随着降水入渗量和地表灌溉渗漏量增加，地下水总补给量有所增加；随着人工开采量的减少，地下水总排泄量减少；地下水补排关系改善但仍呈现负均衡。研究可为今后建立海河平原区地下水合理开采模式提供依据，促进区域地下水可持续管理。 　　【关键词】海河平原区；MODCYCLE模型；浅层地下水；动态响应；供水格局 　　0 引言 　　海河平原区是我国经济发展的重要区域，地下水一直是主要供水水源，且供水比重也呈稳定增长趋势，近年来更高达66%。自80年代以来，在需水量迅速增加和降水衰减的共同作用下，海河平原区已经成为南水北调受水区地下水超采最为严重的区域[1]。长期无序过量的开采地下水资源，导致海河平原区地下水储量大量消耗，区域地下水水位持续下降，并引发严重的地面沉降、海水倒灌、水质污染等环境地质问题[2]。为确保未来海河平原区地下水的可持续利用，保障区域稳定健康发展，多年来学术界一直将当地地下水评价与研究作为关注热点。 　　韩瑞光研究建立了海河平原区浅层地下水概念模型，并提出今后模型建设建议[3]。费宇红等通过研究海河平原区地下水储量消耗过程，指出该区域地下水可开采利用的潜力已经十分有限，从长远看南水北调是解决缺水的理想途径[4]。何杉采用水量平衡的方法，研究分析了南水北调实施后，地下水开采量的减少与入渗补给量的增加，将促使海河平原浅层地下水局部得到恢复[5]。杜思思等联合运用MODFLOW与水资源配置模型ROWAS，模拟了有无南水北调两种对比情景下海河平原区地下水的演变[6]。 　　以上研究通过数据分析与模型模拟等方法对海河平原区的地下水资源作出了评价，但作为模拟情景水文条件的水文系列较短，考虑的情景方案较少。为从更完整的角度验证工程达效对海河平原区地下水循环恢复所起的作用，本文基于分布式水文模型MODCYCLE，结合多个典型的供水格局情景进行海河平原区地下水的详细模拟与动态响应分析。 　　1 海河平原区MODCYCLE模型的构建与验证 　　MODCYCLE模型是基于“自然――社会”二元特性开发的分布式水循环模拟模型[7]，充分考虑到对自然水循环过程与人工水循环过程的双重体现[8]，可用于人类活动干扰明显的海河平原区水循环系统的模拟量化。为保证水循环模拟的完整性，本文通过MODCYCLE构建海河流域水资源转化动态模拟模型，研究和辨析现状2001～2010年海河平原区浅层地下水动态平衡；选取5个代表性水资源配置方案，模拟预测不同水文系列条件（1956～2000年平水系列、1980～2005年近期枯水系列）和南水北调工程实施情况（南水北调中线工程一期达效、二期达效和加大中线一期引水20%）下海河平原区浅层地下水动态响应。 　　1.1 模型数据输入 　　按DEM将海河流域划分为2028个子流域，其中平原区子流域1165个。地下水数值模拟以4km为间距划分网格单元，有效单元格8383个。模拟气象数据采用收集的46个气象站点实测数据展布。地下水水位根据550个浅层地下水位观测井和210个深层地下水位观测井的观测数据插值计算。水文地质参数根据海河流域水文地质调查数据展布。 　　1.2 模型率定与验证 　　模型以2001～2005年为率定期，2006～2010年为验证期。考虑到海河流域水循环特性，选取地下水位、地下水蓄变量为验证指标。 　　1.2.1 地下水位检验 　　图1所示为2010年末（验证期末）的实测与模拟浅层地下水位等值线对比，从整体上看，模拟与实测地下水位等值线具有可比性，山前及中部地下水开采密集区的地下水位等值线变化幅度大。 　　1.2.2 浅层地下水蓄变量检验 　　2001～2010年海河流域浅层地下水蓄变过程统计值（根据2001～2010年《海河流域水资源公报》分析整理）与模拟值对比如图3。从蓄变模拟结果看，蓄变过程在变化趋势上一致。经计算得，浅层地下水蓄变量模拟与统计值之间相关系数为0.96，相关程度较高。 　　从总体上看，对于海河流域这种大空间尺度和长时期的水循环模拟研究，目前的率定验证结果基本满足要求。 　　2 地下水平衡现状与模拟情景设置 　　2.1 2001～2010年现状浅层地下水动态平衡 　　模拟现状年时段海河平原区浅层地下水年均补给总量约193.66亿m3。其中降水入渗量占总补给量的67.0%，为最主要的补给来源；灌溉渗漏补给量占8.7%。浅层地下水年均排泄总量22