第八讲SWAT的非点源污染模拟.ppt

* * * * * * * * * * （4）模拟结果 水闸模拟 蚌埠闸 11座水闸模拟结果 （4）模拟结果 王家坝 23个水文站模拟结果 （4）模拟结果 王家坝 23个水文站模拟结果 （5）小结 通过对淮河流域的径流模拟表明，在率定期水库断面的平均相关系数和平均效率系数为0.72和0.47，水闸为0.83和0.67，水文断面为0.88和0.73。而在检验期内水库的平均相关系数和效率系数为0.68和0.40，水闸为0.81和0.51，水文断面为 0.80和0.56。水文站的模拟效果最好，其次为水闸，最差的为水库。但模拟结果均在可接受的范围内。 氨氮模拟相关系数在0.45以上的断面共19个，占53%，平均相关系数为0.46；CODMn模拟相关系数在0.45以上的断面共22个，占61%，平均相关系数为0.52。 温榆河流域水量水质联合调度 北京市水资源严重短缺，已成为经济社会发展的制约因素； 温榆河是唯一一条发源于北京市且常年有水的河流，流域生态水环境状况直接影响该区域经济社会的协调发展和功能定位。 河道水体污染严重，水生态系统退化明显；流域内所有河道水体水质全部为劣Ⅴ类水体 土壤类型 土地利用 排污口 气象信息 资料收集 子流域：39个 流域出口：北关闸 闸坝：10座 北关闸、上庄闸、沙河水库、尚信橡胶坝、郑各庄橡胶坝、十三陵水库、鲁疃闸、桃峪口水库、苇沟闸、辛堡闸 94个主要排污口 55座污水处理厂和再生水厂 参数 名称 定义 文件 取值范围 水文 CN2 第二种土壤湿润条件下的CN值 .mgt （43，96） ESCO 土壤蒸发修正系数 .hru （0.005，0.845） EPCO 植物蒸腾修正系数 .hru （0.005，0.845） SOL_AW 土壤有效蓄水能力 .sol （0.10，0.25） SOL_K 土壤的水力传导度 .sol （0.05，1.68） SOL_Z 土壤层厚度 .sol （96，1040） 水质 BSETLR COD在水库中的沉降速度 .lwq （-100，100） NSETLR 氮在水库的沉降速度 .lwq （-500，5.5） RK1 20°C 时CBOD 氧化速率 .swq （0.0，3.4） RK3 20°C 时CBOD的沉降速率 .swq （-0.10，0.36） BC1 20°C 时氨氮的氧化速率 .swq （0.05，0.10） BC3 20°C 时有机氮转化为氨氮的速率 .swq （0.20，0.40） CtoB 污水的可生化性系数 .swq （1.2，2.3） ERORGN 土壤中有机氮的富集率 .hru （1.8，3.5） ECBOD 土壤中有机污染物的富集率 .hru （0.1，1.5） 参数选择与率定 水量 氨氮 COD 水量水质模拟：沙河闸断面 水量 氨氮 COD 水量水质模拟：北关闸断面 断面 水量模拟 氨氮 CODMn 相对误差 相关系数 效率系数 相对误差 相对误差 沙河闸 8% 0.57 0.40 53% 30% 鲁疃闸 - - - 33% 56% 辛堡闸 - - - 34% 70% 苇沟闸 - - - 35% 44% 北关闸 -5% 0.86 0.84 15% 69% 断面 多年平均 不同保证率来水量 2004-2005年年平均径流模拟 P 50% P 75% P 95% 沙河闸 7830 6290 4530 3600 7348 鲁疃闸 13160 10870 8260 6880 13244 北关闸 52150 48320 41230 35490 62913 水量模拟：沙河闸和北关闸的年径流总量都接近于断面多年平均径流量；北关闸模拟相关系数和效率系数都在0.80以上。 水质模拟：沙河闸、鲁疃闸、辛堡闸、苇沟闸和北关闸的断面氨氮模拟平均相对误差为34%，COD的平均相对误差为54%。 水量水质模拟结果 基于耦合闸坝系统的温榆河流域SWAT模型，通过沙河闸、鲁疃闸、新堡闸、苇沟闸、北关闸等多闸坝的调控，进行水污染修复与整治，使温榆河干流水体水质浓度尽量达到所在水功能区划的浓度标准，从而提高流域内水资源的可利用量，实现流域水资源的高效利用。 闸坝优化调度模型与多情景分析 目标函数： 当流域污染严重时，下游水体的等标污染指数最小 库容约束 闸坝下泄流量约束 水量平衡约束 约束条件： 模型求解采用遗传算法 该指数数值小于等于1，说明该断面水体没有受到污染；数值大于1，说明该断面水体已受到污染,数值越大，超标越严重。 当流域污染较轻时，总达标水资源量最大 等标污染指数 max 氨氮 达标氨氮浓度 COD 达标COD浓度 ， 闸坝调度模型 采用Fortran语言编译QCmode模型，然后将其嵌入SWAT源程序中。二者之间以文本读入读出形