1. 引言
近年来，随着我国经济社会持续高速发展和人口大幅增加，生产、生活及生态用水矛盾突出，生产、生活用水挤占河道生态环境用水现象较为普遍。许多河道断面的生态需水得不到满足，水环境承载力下降，加之经济社会发展带来的污染物大幅增加，河流污染日趋严重。水资源与水环境承载力不足问题凸显，这已成为制约流域经济社会可持续发展和水生态文明建设的瓶颈 [1] 。
河流水生态环境问题诊因复杂，不仅涉及岸上、水下、上下游，还涉及水利、环保、住建等部门，水环境的治理具有流域性、系统性及跨部门性等特点。传统的水环境治理存在局部化、零散化、线性化的问题，缺乏系统治理思路和跨部门的统筹，缺乏流域尺度关联性的思考，缺乏水质水量一体化的考虑，措施单一，难以做到精准治理。治理结果往往存在局部性、短效性的问题，导致屡治屡败，不能达到真正改善流域水环境的目的 [2] 。因此，唯有从流域的视角看待入河污染物的产生、输移及其对水环境的影响，按照“节水优先、空间均衡、系统治理、两手发力”新时期治水思路，坚持“共抓大保护、不搞大开发”，严守以水量、水质、水生态要素为核心的流域生态保护红线，提出跨部门、跨专业的流域水环境治理模式，才是解决流域水生态环境问题的根本策略 [3] 。
本文以成都市沱江流域水环境综合治理为例，提出“大流域统筹、小流域分区、差异化施策”的流域水环境综合治理思路和分区治理方案，并采用MIKE11一维河网水环境模型，以考核断面水质达标为约束，利用模型的反馈机制，对规划方案进行优化，提出成都市沱江流域水环境综合治理分区方案和分区治理标准，用以指导成都市沱江流域的水环境综合治理项目的实施。
2. 研究区基本情况
成都市沱江流域跨越成都市东北部和东南部，范围涉及都江堰、彭州、郫都等11个区(市)县，沱江水系分支众多(见图1)，沱江流域成都段主要支流有青白江、湔江、毗河、北河、资水河(阳化河)、绛溪河等河流。沱江属长江一级支流，北河、青白江、毗河在金堂县境内赵镇汇流后始称沱江。成都市沱江流域自金堂县官仓镇满天星起，至简阳市新市镇界牌村，河段长146 km，流域面积6458 km2，涉及沱江干流、34条沱江一级支流、176条二级支流、219条三级支流和181条四级支流。
Figure 1. Water system map of Tuojiang river basin of Chengdu city
图1. 成都市沱江流域水系
3. 流域综合治理规划思路
规划思路：以现状调查为基础，以问题为导向，以研究范围内河流国考、省考及市考断面水质达标为刚性要求。坚持“大流域统筹、小流域分区、差异化施策”的流域治理思路和“水环境、水生态、水资源、水安全”全要素覆盖的系统治理体系，实现对入河污染物“源头减排、过程控制、末端治理”的全过程控制，并通过建立入河污染物—水质目标响应关系模型，以水环境容量为约束，科学提出分区、分期治理措施及治理标准，真正做到“精准施治，一河一策”。
技术路线(见图2)：1) 开展基础数据收集，并在此基础上系统分析现状水问题；2) 进行控制小流域单元划分，把复杂的流域水环境问题分解到各小流域单元，结合各小流域单元水污染成因、污染源组成及国家、地方相关政策规划，按照“本地产污，本地削减”原则，制定污染物分区减排方案及治理标准(初始方案)；3) 按照该方案及治理标准构建水量水质数学模型，对河流水质进行动态模拟；4) 进行目标可达性分析，坚持“以量(排放量)定标(治理标准)”，不可达则反馈回方案层，重新制定方案、优化措施；达标则形成流域水环境分区治理方案及治理标准(优化方案)，依靠河湖长制实施执行减排方案，实现精准治理。
Figure 2. Technology roadmap of the planning strategy on water environment comprehensive treatment for the river basin
图2. 流域水环境综合治理规划方案制定技术路线图
4. 流域一维水环境模型
4.1. 治理目标
根据成都市人民政府与四川省环境保护厅签订的《成都市水污染防治目标责任书(2016年4月)》的断面水质考核目标，并结合《成都市水务“十三五”规划》、《沱江流域水污染防治规划2017~2020》、《成都市环境保护“十三五”规划》、《四川省沱江“一河一策”管理保护方案》等相关规划及文件要求，确定驾虹、致和电站、三邑大桥、拦河堰、西河天平、毗河二桥、三皇庙、宏缘、爱民桥、新市、临江寺、红日大桥断面2020、2035年的水质目标均为地表水III类(宏缘断面总磷除外，为0.22 mg/L)，具体分布见图1。
4.2. 模型选取
为建立流域污染物排放与水质考核断面之间的输入-响应关系，需选用合适的水动力和水质模型进行模拟预测。沱江流域面积大、河网水系复杂，本研究采用目前世界上领先且经过实际工程验证最多的、被水环境研究人员广泛认同的DHI公司MIKE11模型 [4] [5] [6] 。
4.3. 模型搭建及率定
4.3.1. 河网概化
在建立模型时，需要按照控制断面及河网拓扑关系进行概化。河网概化的原则是能基本反映天然河网的水力特性，建立模型时将流量较小的河道进行合并、概化，河道纵比降通过控制断面的高程进行控制，并根据模型需要进行适当平顺处理。概化时将主要的输水河道纳入计算范围，将次要的河道和水体根据等效原理，归并为单一河道和节点，使概化前后河道的输水能力相等、调蓄能力不变。
经概化，成都市沱江流域计算区域的主要河流为沱江干流、支流水系，包括沱江干流、毗河、青白江、西江河、清溪河及绛溪河等42条河流。基于对流域内现状工业、污水处理厂、城镇生活、农村生活、畜禽养殖以及农田面积分布等信息的调查，对研究区域排污口进行概化，共概化排污口42个，位置分布及河网概化见图1。
4.3.2. 模型输入条件
为全面分析枯水条件下流域主要河道水质状况及变化规律，本次计算选取枯水年(90%)进行分析，典型年选择是通过成都市60多年降雨资料进行频率分析确定。根据边界水质现状条件，都江堰、石堤堰分水质边界取II类，跨界河流水质边界均按照地表功能区要求按III类选取，跨界断面水质超标不予考虑。
4.3.3. 模型参数率定
采用2016年逐月水文实测资料及通过降雨、蒸发资料采用产汇流模型计算得到的边界条件对水动力模型进行率定，采用试错法调试各河道的糙率。采用2016年逐月水文实测资料进行模型主要断面的率定结果分析，三皇庙断面流量平均相对误差小于9.56%，所建模型模拟水文变化过程的效果较好。
在水动力模型的基础上，采用2016年成都市环境监测站进行的逐月水质监测结果对水质模型进行率定，得到COD、NH3-N、TP降解系数分别为0.2、0.15、0.1 d−1。去除监测数据中个别极值数据，主要断面COD、NH3-N、TP模型计算值与实测值平均相对误差分别为16.99%、28.88%、28.17%，主要污染物计算值与实测值吻合性较好，模型基本可靠，可用于预测。
5. 规划方案模拟分析
5.1. 规划方案
在成都市沱江流域水环境综合治理规划中，针对各流域片所存在的问题，制定了以下规划方案：1) 沱江上游青白江片区，以乡镇污水与农村面源治理为主；2) 沱江中游毗河及西江河流域片，以城市点面源治理、农村面源治理、河道清淤、生态补水为主；3) 沱江下游以城市控源截污、农村面源污染治理并重。
5.2. 模拟情景
本论文根据水文年型(枯水年)、水平年(2020年、2035年)以及规划方案(初始方案、优化方案)，设定四种模拟情景，分别是：初始方案下的2020年情景和2035年情景；优化方案下的2020年和2035年情景。
5.3. 水质模拟分析
1) 初始方案
在成都市“十三五”沱江流域监测断面中，将入境断面与源头水断面以外的断面纳入本次水质模拟，共选取了12个监测断面，即驾虹、致和电站、三邑大桥、拦河堰、西河天平、毗河二桥、三皇庙、宏缘、爱民桥、新市、临江寺、红日大桥断面。各监测断面不同水平年的污染物浓度模拟值见图3。
模拟结果表明，青白江流域驾虹、致和电站断面、沱江流域爱民桥断面在不同水平年，均能达到地表水III类标准；西江河流域西河天平、毗河流域拦河堰、毗河二桥断面在2020年、2035年均不能达标，超标因子为总磷、氨氮；青白江流域三邑大桥、阳化河流域红日大桥、沱江流域三皇庙、宏缘、新市、临江寺断面在2020年不能达标，超标因子主要为总磷，在2035年均能达标。
Figure 3. Simulation value of pollutant concentration at each inspection section for the initial scheme (left: 2020, right: 2035)
图3. 初始方案中各考核断面2020、2035水平年的污染物浓度模拟值(左：2020年，右：2035年)
2) 优化方案
初始方案中，青白江流域三邑大桥、阳化河流域红日大桥、沱江流域三皇庙、宏缘、新市、临江寺等6个断面在2020年不达标，西江河和毗河的3个考核断面在2020年和2035年均未达标，考虑到毗河和西江河水质情况不容乐观，因此，规划提出对初始方案进行如下改进：
a) 在2020年，毗河枯水期(11月~4月)从石堤堰补充7 m3/s的水量，西江河全年从西江河源头补充3 m3/s的水量，补水水质均为II类；西江河、毗河流域所有日处理水量大于1万m3的污水处理厂尾水经人工湿地等处理后达到地表水III类标准；毗河、西江河流域全域禁止畜禽养殖，毗河、西江河流域农业面源污染削减分别为15%和30%；沱江流域散养畜禽养殖削减30%。
b) 在2035年，毗河枯水期(11月~4月)从石堤堰补充7 m3/s的水量，西江河枯水期(11月~4月)从西江河源头补充3 m3/s的水量，补水水质均为II类；西江河、毗河流域所有日处理水量大于1万m3的污水处理厂尾水经人工湿地等处理后达到地表水III类标准；毗河、西江河流域全域禁止畜禽养殖；毗河、西江河流域农业面源污染削减50%；沱江流域散养畜禽养殖削减30%。
实施上述优化措施后，原未达标的9个断面均能达到地表水III类标准。2020年、2035年原未达标断面污染物浓度年均模拟值见图4。
Figure 4. Simulation value of pollutant concentration at each inspection section for the optimization scheme (left: 2020, right: 2035)
图4. 优化方案中各考核断面2020、2035水平年的污染物浓度模拟值(左：2020年，右：2035年)
6. 结语
本文通过对沱江流域水环境综合治理规划研究，得到几点认识：1) 针对大流域尺度，系统而完整地提出了“大流域统筹、小流域分区、差异化施策、全过程控制、全要素覆盖”的流域水环境治理思路。2) 采用MIKE11建立流域一维水环境模型，利用模型对入河污染物与控制断面水质间的响应与反馈机制，可对规划方案进行精细模拟分析。3) 本研究思路应用于沱江流域水环境综合治理工作实践表明，通过方案优化后，提出的增加毗河流域生态环境补水、提高尾水排放标准等综合治理措施，可使各考核断面水质达标，为流域内各区县河长开展水环境治理工作提供重要抓手。4) 本文提出的流域水环境综合治理思路，统筹了环保、水利及住建等部门的要求，体现了水环境治理的流域性和系统性，对类似工作具有积极的借鉴意义。
基金项目
国家重点研发计划资助项目(2016YFC0402005)。