基于分布式水文模型的水功能区设计流量研究.docVIP

基于分布式水文模型的水功能区设计流量 胡鹏周祖昊贾仰文周娜 中国水利水电科学研究院 流域水循环模拟与调控国家重点实验室，北京 100038 摘要：提出了基于分布式水文模型的水功能区设计流量计算方法。将此方法应用于松花江流域，对流域内300个水功能区的设计流量进行了计算，并取用水、下垫面和调水工程对流域在2020和2030年情景下水功能区设计流量变化情况进行了计算和分析。关键词：水功能区设计流量分布式水文模型松花江流域1 研究背景 为了促进水资源的合理开发和有效保护，我国于2011年颁布了《全国重要江河湖泊水功能区划》，将全国主要河流和湖泊划分为保护区、缓冲区、开发利用区和保留区共4类水功能一级区2888个，其中1133个开发利用区又被划分为饮用水源区、工业用水区、农业用水区、渔业用水区、景观娱乐用水区、过渡区和排污控制区共七类水功能二级区27个。2011年，中央一号文件《中共中央国务院关于加快水利改革发展的决定》中明确提出确立水功能区限制纳污红线，从严核定水域纳污容量，严格控制入河湖排污总量。水功能区纳污容量的确定主要是通过在设定水文条件下，某种污染物满足水功能区水质目标要求所能容纳的最大数量来计算[1-2]，设计流量的大小对于河流水功能区纳污容量有着直接的影响。目前，我国各大流域对于水功能区设计流量的具体核算方法有所差别，但主要均是根据长系列实测数据，选取一定保证率的枯水期流量作为设计流量。采用历史实测数据核算水功能区设计流量保障了结果的可靠性，但同时也存在着以下几方面的问题需要对每个水功能区单独进行计算，导致在实际工作和科学研究中工作量较大；完全通过数据推求，不利于变化环境下水功能区设计流量演变规律的研究；水文实测资料不足或无资料的地区处理。 目前国内外对于水域纳污容量的研究主要集中在给定水文条件下通过水质模型计算水域纳污能力上[3-7]，而对于给定水文条件即设计流量研究较少。徐小燕等[8]以珊溪水库为例探讨了不同取值方式对于设计流量大小的影响，提出了适合南方地区水功能区设计流量的计算方法。在人类活动的水文影响方面，有大量关于人类活动对于河川径流影响规律的研究[9-12]，但少有关于人类活动对于枯水期径流量以及水功能区设计流量的研究，Shalini Oogathoo[13]应用MIKE SHE模型对加拿大Canagagigue Creek 流域进行了模拟，并设定情景进行水资源演变分析，在城市化的情景下洪峰流量减少，枯水期流量增加。同样在分布式水文模型研究方面，绝大多数研究也只限于对于全年径流过程的模拟分析，少有专门针对枯水期径流模拟的研究。针对上述问题，本文提出了一种利用分布式水文模型计算水功能区设计流量的方法，并利用模型对研究区未来情景下枯水期径流量和水功能区设计流量变化情况进行了分析，提出了相应的应对措施。 2 研究区概况和资料收集 松花江流域位于中国东北地区北部，介于北纬41°42′51°38′、东经119°52′132°31′之间，行政区划包含黑龙江省全部、吉林省大部、辽宁省和内蒙古自治区一部分，总面积5万km2。2009年流域总人口近6000万，GDP约1万亿元，是中国重工业城市的集中地、重要农牧业生产基地。流域多年平均降雨量537mm，其中汛期69月的降水量占全年的60%～80%，冬季122月的降水量仅为全年的5%左右。流域地处中高纬度地区，冬季寒冷漫长，年负温期5个月左右，季节性冻土与多年冻土广泛分布，河流一般11月进入冰封期，年4月解冻。 流域水系主要包括北部的嫩江、南部的第二松花江以及二者汇合后的松花江干流三大干流水系。松花江流域共有水功能一级区207个，其中73个开发利用区又分为166个二级区，不同类型水功能区共有300个。松花江流域水功能区一级区及主要水文站点分布如图1所示。对研究区进行历史以及未来情景下分布式水文模拟所需资料主要包括地表高程、土地利用、土壤质地、气象、水文、社会经济、用水、水利工程等，相关数据来源及简要说明如表1所示。 图1 松花江流域水功能一级区及主要水文站点示意 表1 松花江流域历史水文过程模拟数据列表 数据类别 数据来源 格式与说明 地表高程 美国地质调查局 分辨率30 土地利用 中国科学院 1990年和2000年，分辨率30m 土壤质地 全国第二次土壤普查 比例尺1：100000 水文地质 《松辽流域水资源综合规划》、《中国水文地质分布图》 给水度、渗透系数、岩性分区、含水层厚度等 气象 中国气象局 35个气象站1956～2000逐日气温、日照、湿度、风速 降水 松辽委水文局、黑龙江、吉林省水文局、中国气象局 251个降水站1956～2000逐日降水观测资料 径流 《松辽流域水资源综合规划》 36个水文站1956～2000逐月实测径流和还原径流资料 大型灌区