基于EFDC模型的重庆库区水环境容量研究

三峡生态环境监测

EcologyandEnvironmentalMonitoringofThreeGorges2019年6月Jun.2019□研究论文DOI：10.19478/j.cnki.2096-2347.2019.02.02

基于EFDC模型的重庆库区水环境容量研究

郑瑶*，胡学斌，何

（重庆大学

强，何劼韫

三峡库区生态环境教育部重点实验室，重庆400045）

摘要：基于水环境数学模型EFDC（environmentfluiddynamicscode），建立了重庆库区长江干流的水动力水质模型，通过

2017年监测水文、水质资料对模型进行了率定和验证，模拟结果与监测资料吻合度较好，表明水动力水质模型计算结果较为合理可靠。根据该模型量化各个排污口对水质控制断面的响应系数，以地表水Ⅱ类水标准为规划目标，用线性规划法求解水环境容量。计算结果表明，重庆库区化学需氧量（COD）、氨氮（NH3-N）、总磷（TP）的水环境容量分别为776t/a，11038t/a和922t/a。重庆库区干流水动力水质模型的建立与水环境容量的计算，为三峡库区的水环境保护工作提供了一个重要的管理和分析平台，对于三峡库区的污染物总量控制具有重要的意义。关键词：EFDC；重庆库区；水环境容量；长江干流

引用格式：郑瑶，胡学斌，何强，等.基于EFDC模型的重庆库区水环境容量研究[J].三峡生态环境监测，2019，4（2）：9-16.Citationformat:ZHENGY,HUXB,HEQ，etal.WaterenvironmentcapacityinThreeGorgesReservior(Chongqingsection)basedonEFDCmodel[J].EcologyandEnvironmentalMonitoringofThreeGorges，2019，4（2）：9-16.中图分类号：X321

文献标识码：A

文章编号：2096-2347（2019）02-0009-08

WaterEnvironmentCapacityinThreeGorgesReservior(Chongqing

Section)BasedonEFDCModel

(KeyLaboratoryofEco-environmentsofThreeGorgesReservoirRegion(MinistryofEducation),ChongqingUniversity,Chongqing400045,China)

ZHENGYao*,HUXuebin,HEQiang,HEJieyun

Abstract:BasedonthewaterenvironmentmathematicalmodelEFDC(environmentalfluiddynamicscode),the

hydrodynamicandwaterqualitymodeloftheYangtzeRiverinChongqingsectionwasestablished.Themodelwascalibratedandverifiedbyusingthehydrologyandwaterqualitydatamonitoredin2017.Thesimulationre⁃sultsareingoodagreementwiththemeasureddata,indicatingthatthemodelcalculationsarereasonableandre⁃liable.Accordingtothemodel,theresponsecoefficientwasquantifiedforeachsewageoutlettothewaterqualitycontrolsection,andthelinearenvironmentmethodwasusedtodeterminethecapacityofwaterenvironmentbytakingClassIIofthenationalstandardsforsurfacewaterqualityastheplanningtarget.Thecalculationresultsshowthatthewaterenvironmentalcapacityforchemicaloxygendemand,ammonianitrogenandtotalphosphorusinChongqingsectionare776t/a,11038t/aand922t/a,respectively.TheestablishmentofthehydrodynamicwaterqualitymodelfortheYangtzeRiverinChongqingsectionandthesimulatedcapacityofwaterenvironmentprovideanimportantplatformformanagementandanalysisofwaterenvironmentalprotection,andisofgreatsig⁃nificanceforthetotalmasscontrolofpollutantsintheThreeGorgesReservoirarea.

Keywords:EFDC;Chongqingreservoirarea;waterenvironmentcapacity;theYangtzeRiver

收稿日期：2019-02-26

基金项目：国家水体污染控制与治理科技重大专项（2014ZX07323001）。作者简介：郑瑶（1995—），女，硕士研究生，主要从事流域水环境治理方案研究。E-mail:20161702013@cqu.edu.cn

10三峡生态环境监测http://sxhc.cbpt.cnki.net/第4卷了《三峡后续工作规划》，正式进入后三峡时期，三峡库区主要任务由三峡工程建设和移民转变为库区生态环境保护和库区移民安稳致富。由于三峡工程关闸蓄水，水位上升流速减缓，三峡库区前端更接近人工湖泊的性质而不是河道的性质，2011年5月18日，常务会议讨论通过

型，根据库区不同蓄水位水库调度运行方式下水文特点设计模拟边界条件，科学分析识别重庆库区不同控制单元河段水环境容量。

1

1.1

研究区域与EFDC模型概述

研究区域概况

自净能力减弱，水环境容量减少，三峡库区内的污染输入对水质的影响程度增加，从而导致了支流部分区域富营养化和水华等亟待解决的环境问题，库区水环境较脆弱[1-2]。

治行动计划》2015年，党、发布了，明确提出强化环境治理目标管理，《水污染防

而水环境容量计算是环境目标管理的核心内容，它在建立污染物排放与水质响应关系的基础上，计算水体可以容纳的最大污染物量，从真正意义上将污染物控制与水质目标紧密结合起来[3]。近年来已有学者对三峡库区进行了水环境容量研究分析。刘昭伟等[4]通过模拟三峡水库万州段水质状况，提出一种岸边水环境容量的计算方法，计算出蓄水前后岸边水环境容量，分析了三峡工程对水库岸边水环境容量的影响。重庆交通大学陈培帅[5]运用MIKE21软件建立重庆主城区两江水质模型，以排污口污染混合区总长度为控制指标，计算出重庆主城两江水环境容量，并结合重庆主城各区经济结构，得到与主城各区经济发展相适应的水环境容量。丁世敏等[6]采用允许排污量法，计算了重庆市涪陵区“十三五”期间的水环境容量，为三峡库区城镇的水污染控制和水环境管理提供更有效的支持。王晓青等[7-8]采用垂线平均的二维水流、水质数学模型，对小江在三峡库区蓄水前后化学需氧量（COD）、总氮（TN）和总磷（TP）水环境容量变化进行了研究，发现蓄水后TP水环境容量减小幅度最大，故实施污染物总量控制要特别注意磷污染物的削减。

从以上研究可以看出，关于三峡库区小流域水环境容量研究常见，而大流域的水环境容量鲜有模型方法的应用研究。为更清晰地了解三峡库区水质状况及纳污能力，并制定合理的流域水环境治理整体解决方案，本研究以对三峡库区有重要影响的重庆段干流部分为研究对象，基于水环境数学模型EFDC（environmentfluiddynamicscode）软件建立了重庆库区干流的水动力和水质模

三峡库区重庆段东起巫山县、西至江津区、

南起武隆区、北至开州区，地理范围为北纬28°31′个重庆区县，长约～31°44′、东经560km105°49′，面积约～110°12′4.62×10，包括4km2，约22占整个三峡库区面积的85.6%，由此可以看出，它占据着整个三峡库区至关重要的生态位置[9-10]。通过对历年三峡水库水体水质状况统计分析，自蓄水以来，重庆库区干流水质总体较好，且近年来有向好的趋势；沿程各断面水质变化趋势不完全一致，但从上游到下游，水质评价结果优于Ⅲ类的比例总体上呈现上升的趋势[11]。1.2

EFDC模型概述

水环境研究直接面临的问题是空间尺度大、污染物成分多及环境过程复杂，如果只是通过实验分析或现场监测具有较大的难度，但通过建立数学模型模拟水环境的数值，借助数学工具和方法可以很好地解决这个困难[12]。EFDC模型是一个开源的地表水建模系统，有机地将水动力、泥沙、污染物和水质等模块用单一的源代码完全集成。该模型由美国弗吉尼亚海洋科学研究所（VIMS）和威廉玛丽学院的海洋科学学院于1988年共同开发，后由美国国家环境保护局资助，改用Fortran95进行再次开发优化。在这之后又先后经历了EFDC_1D，EFDC_Hydro进，直到现在的，EFDC_GVCEFDC_Explore，EFDC_DSI（EE），并且逐渐发等版本的改展为各研究所、高校、环保部门广泛应用的水环境模拟工具[13]。

有毒物质及沉积物等模块组成。其中，水动力模

EFDC模型主要由水动力、水质、标量输运、

块是EFDC模型建立的基础，其余模块都是建立在水动力模块的基础上。在水动力模块中主要涉及了淡水流、大气作用、水深、表面高程、底摩擦力、流速、湍流混合、盐度、水温9大部分，在本研究中，涉及了淡水流、水深、底摩擦力、流速、湍流混合以及水温部分。作为研究主要部分的水质模块，包括了碳的各种组成，以及氮、磷、溶解

第4卷第2期郑瑶，等：基于EFDC模型的重庆库区水环境容量研究

11氧、悬浮藻类等21个水质变量，在本研究中主要涉

及了氨氮（NH3-N）、TP以及COD这3个水质变量。

算的网格为10611个。

由于三峡库区河底地形复杂且水深，所以水下地形的获取有较大难度。江辽等[14]在研究感潮河段数值模拟时，利用模型对概化后河道的水动力状态如河长、河宽以及河深进行数值试验，并与实际观测值进行了比较。得出参数中河宽对潮位和流量的影响最大，其次是河长。因此，本研究在建立模型时为了优化计算方法，提高计算效率，在考虑河底地形时忽略水下地形起伏并将河底视为平底。重庆库区干流纵断面如图1所示。

2

2.12.1.1

重庆库区干流水环境模型

模型构建计算网格生成

本研究选取三峡库区重庆段干流为模拟的计

算区域，采用由荷兰DelftHydraulic公司开发的

Delft3D的RGFGRID模块作为贴体网格的生成工网格文件（*.grd）导入EFDC软件中，总共参与计

具，调整优化后再将生成的重庆库区河流的二维

ZYX

8.991

WaterColumn

[Time138.000]

45

SecchiDepth(m)

VerticalExaggeration:10

图1

Fig.1

河底概化后的纵断面

Thelongitudinalsectionoftheriverbottom

2.1.2

工构筑物以及开边界。在本研究中，将上游来水、重庆库区沿线流域面积超过100km2的河流、沿线的日废水流量超过1.0×105t的污染源这三部分均按照水流入流输入的方式设定；出流边界设定为开力水质模型的输入条件。

三峡库区重庆段干流共62个入流口，包括进水口、40处支流入流、21处点源污水入流。进水口是江津上游的区域来水，其流量与水质数据来自重庆市水利局公布的今日八时水情以及环保部公布的国控断面水质数据；支流流量与污染负荷是根据水利监测站提供的监测数据或文献资料计算分析得出[15]；排污口流量与污染物排放量来自重庆市生态环境监测中心公布的污染源监督性监边界，充分利用流域监测数据作为重庆库区水动

EFDC模型的边界条件有三类：水流入流、水

边界条件设定

测信息；下游出水口采用巫山水文站的水位数据。将这些数据以时间序列格式输入EFDC中，模拟时长为1年，设定的边界时间序列以2017年1月1日为起点，模拟时间步长为1s，CPU运算时间为67h。2.1.3

EFDC具有很好的通用性，数值计算能力强，相关参数的取值

大多数情况下，EFDC模型中的许多参数不需要修率。参照相关文献资料，三峡库区的河道糙率取

改[16]。水动力模型中所涉及的关键参数为河道糙值范围为0.032~0.037[17-18]，经多次调试和测算后，确定河道糙率均考虑为0.035。水质模型各项初始参数值主要根据相关文献以及EFDC的模拟实例确定[19]。模型选用第1~183天（即2017年1—6月）的模拟数据与丰收坝、和尚山、李渡取水口3个国控断面的监测数据（环保部公示的监测数据）进行参数率定，结果见表1。

12三峡生态环境监测http://sxhc.cbpt.cnki.net/第4卷表1

EFDC内部参数名KROKTRKRPKLPKDPKRNKLNKDNKHNitNKHNitDOTNitKNit2rNitMKHCODKCD

Tab.1

水质模块主要参数

文献参考范围

1~5.321~1.20.001~0.10.001~0.10.01~0.10.001~0.010.001~0.10.001~0.10.2~10.01~110~300.0045~0.0060.01~0.80.5~20.01~0.6

率定结果1.61.030.0050.0750.10.0050.020.030.20.220

0.00450.151.490.01

单位——d-1d-1d-1d-1d-1d-1

gN·m-3gO2·m-3

℃—gN·m-3gO2·m-3

d-1

Themainparametersofwaterqualitymodule

代表意义复氧速率

复氧速率的温度效应系数稳定颗粒态有机磷最小水解率不稳定颗粒态有机磷最小水解率

溶解有机磷最小水解率稳定颗粒态有机氮最小水解率不稳定颗粒态有机氮最小水解率

溶解有机氮最小水解率NH3-N对硝化作用半饱和常数DO对硝化作用半饱和常数

硝化参考温度

硝化反应理想温度效应常数

最大硝化速率

COD降解需氧半饱和常数

COD降解速率

2.2模型验证

2.2.1水动力模型验证

在许多河流、河口、海湾及河漫滩模型中，水位是进行水动力模型验证的首选指标[20]。本研

177.50175.00172.50170.00167.50

QingxichangElevation1:1Correlation

究选取涪陵区清溪场、万州区万县断面水位监测值，与模型数据进行了对比，图2为各断面实测水位与模拟水位相关性对比结果，表明实测水位与模拟水位之间的相关系数均达0.9以上。

177.50175.00172.50170.00167.50

WanxianElevation1:1Correlation

水位计算数据162.50160.00157.50155.00152.50150.00147.50

NumberofData/ModelPairs:3Statistic:R-Squared:0.9932

水位计算数据165.00

165.00162.50160.00157.50155.00152.50150.00147.50

NumberofData/ModelPairs:3Statistic:R-Squared:0.9932145.00

145.00147.50150.00152.50155.00157.50160.00162.50165.00167.50170.00172.50175.00177.50145.00

145.00147.50150.00152.50155.00157.50160.00162.50165.00167.50170.00172.50175.00177.50水位实测数据（a）涪陵区清溪场

水位实测数据（b）万州区万县

图2

Fig.2

水位计算数据和实测数据相关性对比结果

Comparisonandcorrelationbetweenthesimulated(verticalaxis)andobserved(horizontalaxis)waterlevels

2.2.2

水质模型验证

模型的验证采用第183~365天丰收坝、和尚

相对误差分别为21.83%，18.81%，14.96%。从总体上看，三峡水库水质模型模拟结果与监测值误差不超过30%，在合理误差范围内，证明该模型较好地体现了污染物的时空分布规律，基本反映了重庆库区干流主要水质变化过程，可用于实际的水质预测分析和水环境容量的计算[21]。

山、李渡取水口3个国控断面的监测数据对模型进TP浓度模拟值与实测值对比验证的结果，结果显

行验证。图3展示了和尚山监测站点COD、NH3-N、示，COD、NH3-N、TP浓度模拟值与实测值的平均

第4卷第2期郑瑶，等：基于EFDC模型的重庆库区水环境容量研究

133.23.0浓度模拟值浓度实测值

1）-L2.8·g（m2.6/度2.4浓DO2.2C2.01.81.61.4

180195210225240255时间270/d

285300315330345360

0.200.18浓度模拟值浓度实测值

1）0.16-L·g0.14（m/度0.12浓N-0.103HN0.080.060.04

180195210225240255时间270/d

285300315330345360

0.140.130.12浓度模拟值浓度实测值

1）-L·0.11g（m0.10/度浓0.09PT0.080.070.060.050.04

0.03

180195210225240255时间270/d

285300315330345360

图3和尚山监测站点COD，NH3-N，TP浓度模拟值与实

测值对比验证结果

Fig.3

centrationsComparisonofCOD,betweenNHthesimulatedandmeasuredcon⁃

3-N,TPatHeshangshanstation

2.3模型结果分析

重庆库区入流量较大，多年平均入流量为

8速受季节、地理位置以及较大支流汇入影响较大。670m3/s，从模型模拟结果发现，库区内水流流

在冬季枯水期，库区外来输入的流量较小且运行水位较高、储水量大，使得库区流速整体较缓；夏季与冬季条件相反，输入的水量较大且库区范围在夏季属于降水量较大且较为集中的区域，故而夏季库区整体流速较快。在成坝后，水流速度总体上是随监测站点位置向三峡大坝方向推移逐渐变缓，江津至重庆主城段常年流速为0.5～1.5m/s，涪陵至万州段流速为0.08～1m/s，云阳至巫山段十分接近湖泊型水库，流速为0.03～0.2m/s，水流相对较缓。重庆库区沿程都有支流汇入长江干流。比如，在重庆主城段水流速度突然增加，主要原因是在朝天门区域有嘉陵江汇入，导致流量迅速增加；在涪陵段因为乌江的汇入也有类似的情况。嘉陵江以及乌江等主要支流的共同作用，在一定程度上增大了库区内重庆主城至丰都段的水流速度。

通过对污染物浓度场的模拟发现，在不同运行情况及地理位置下，水质状况以及污染带状况也有所不同。江津至重庆主城段区域与云阳至巫山段区域比较，污染带的长度逐渐变长，面积逐渐增大，影响的范围也逐渐增加。145浓度进入江水后削减的幅度较小，形成的污染带m水位运行条件相比，COD、NH175m水位与3-N及TP的面积更大。云阳至巫山段区域在175m水位时的支

流入河口处的污染物浓度更高，形成的污染带范围更广，有多个污染带叠加形成一个大面积的污染带的情况。

3

水环境容量计算

3.1

计算方法

随着水环境数学模型的广泛应用及计算机技

术的不断进步，水环境容量逐渐形成了解析公式法、模型试错法、系统最优化法、概率稀释模型法、未确知数学法等计算方法。其中系统最优化法包括线性规划、非线性规划、动态规划及随机规划等[22]。水环境容量计算中所采用的主要是解析公式法[23]与线性规划法[24-25]。解析公式法计算简便，适用于水文要素、水力条件变化不大的水域；基于动态水环境模型的线性规划法，自动化程度高、精度高、对边界条件及设计条件的适应能力强，适应水域多。因此，本研究基于EFDC模型选用线性规划法计算重庆库区水环境容量，以

14三峡生态环境监测http://sxhc.cbpt.cnki.net/第4卷期为重庆库区水质保持和水环境治理方案提供科学依据。

线性规划法基本思路为：①基于河流水环境模型，建立排污口污染物排放量和控制断面水质标准浓度之间的动态响应关系；②以污染物最大允许排放量为目标函数，控制断面达标为约束条件，利用线性规划最大最优解的方式求解水环境容量[26]。

目标函数：maxL=∑xiQi(i=1,2,⋯,n)（1）

ni=1

n

质响应系数，量纲为1；xi为排污口允许排放浓度，单位为mg/L；L为所有排污口的污染物负荷排放总量，单位为g/s；Cj0为水质控制断面的污染背景浓度，单位为mg/L；Cj为水质控制断面的环境标准值，单位为mg/L；Qi为排污口设计排污流量，单位为m3/s；Cmax为第i个排污口的浓度上限值，单位为mg/L。3.2计算条件

3.2.1水质保护目标

2008年批准实施的《三峡库区及其上游水污染防治规划（修订本）》，要求三峡库区及

式中：i为排污口编号；n为排污口数目；j为水质控制断面编号；m为水质控制断面数目；aij为第i

ì（2）

ï∑aijxi+Cj0≤Cj(j=1,2,⋯,m)约束条件：íi=1ïî0其上游主要控制断面水质整体上基本达到国家地表水环境质量Ⅱ类标准。根据《全国重要江河湖泊水功能区划（2011—2030年）》、长江流域水资源质量公报及重庆市地表水环境质量市控断面位置，确定长江干流重庆库区段水功能区划、相关控制断面现状水质与水质目标，见表2。

个排污口的单位负荷量对第j个水质控制断面的水

表2

Tab.2

序号12345671011121314151617

水功能区名称

长江上游珍稀特有鱼类自然保护区

大渡口饮用、工业用水区渝中饮用、景观用水区江北工业用水区巴南、长寿保留区长寿工业用水区涪陵李渡饮用、工业用水区涪陵饮用、景观用水区丰都饮用水源区忠县工业、景观用水区忠县、万州保留区万州饮用、工业用水区云阳饮用、工业用水区云阳、奉节保留区奉节工业、景观娱乐用水区

巫山饮用水源区巫山保留区

重庆库区水功能区划及水质目标

起点朱坨丰收坝和尚山寸滩鱼嘴扇陀黄草峡李渡取水点清溪场镇江镇苏家新田镇晒网坝苦草沱龙洞镇白帝城巫山长江大桥

终点丰收坝和尚山寸滩鱼嘴扇陀黄草峡李渡取水点清溪场镇江镇苏家新田镇晒网坝苦草沱龙洞镇白帝城巫山长江大桥

培石

现状水质

ⅡⅡⅡⅡⅡⅡⅢⅡⅡⅢⅡⅡⅡⅡⅡⅡⅡ

水质目标ⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡⅡ

WaterfunctionzoningandwaterqualitytargetforChongqingReservoirareaoftheYangtzeRiver

3.2.2

设计流量

作为水环境容量计算，设计流量的确定主要

各断面平均流速最小，本研究选择枯水期90%保证率的最枯月平均流量作为最不利设计流量，上

28]

游最不利设计流量取3000m3/s[5，。

考虑最不利情况[27]。结合三峡库区不同蓄水位水库调度运行方式对水库水文条件的影响，以及枯水期上游来水量较小和库区运行较高水位使流域

3.2.3

水质响应系数求解

水质响应系数的求解，是以线性规划法求解

第4卷第2期郑瑶，等：基于EFDC模型的重庆库区水环境容量研究

15水环境容量的关键。根据定义，在每个排污口设定1个单位负荷量，而其余排污口无负荷量排出，然后用水质模型计算出在这种情况下的浓度分布，找出对应水质控制断面的浓度平均值，即为水质响应系数aij。改变排污口，重复以上步骤，即可求

表3

Tab.3水功能区

长江上游珍稀特有鱼类自然保护区

大渡口饮用、工业用水区渝中饮用、景观用水区江北工业用水区巴南、长寿保留区长寿工业用水区涪陵李渡饮用、工业用水区涪陵饮用、景观用水区丰都饮用水源区忠县工业、景观用水区忠县、万州保留区万州龙宝饮用、工业用水区云阳饮用、工业用水区云阳、奉节保留区奉节工业、景观娱乐用水区

巫山饮用水源区巫山保留区合计

出每个排污口的水质响应系数[29]。3.3

计算结果与分析

通过线性最优方法求解水环境容量，采用

LINGO程序计算线性规划方程，得到重庆库区各区划COD、NH3-N、TP的水环境容量，见表3。

重庆库区水环境容量计算结果（t/a）

COD

NH3-N

TP155.2043.5276.48107.4481.9915.48-56.9262.3440.59-50.8696.7614.9375.69141.4247.9050.4620.50922.91

ResultsofwaterenvironmentcapacityofChongqingReservoirarea(t/a)

6391.462118.5912665.0611100.673616.5160.481859.365757.033901.79-1017.291935.231752.617450.388199.362739.8213.21819.94776.28

902.333.221292.982234.33452.07504.584.84599.18478.48-101.73483.811.221343.721024.92326.71406.81102.4911038.26

总量COD为776.28t/a，NH3-N为11038.26t/a，TP为922.91t/a，其中涪陵李渡饮用、工业用水区TP背景浓度已高于水质管理目标，水环境容量为

枯水期设计流量下重庆库区干流水环境容量提供决策依据。

4结论

（1）为了研究重庆库区水环境主要污染物及

负值，需削减57t/a，忠县工业、景观用水区分别需削减1017t/a，102t/a，51t/a，才能满足控制断面达标要求。

COD、NH3-N、TP背景浓度均高于水质管理目标，

剩余水环境容量，本研究通过EFDC软件构建重庆库区长江干流的水动力水质模型，研究该河流的水动力场特性及水质变化趋势，模型验证结果较好，可以用于模拟分析重庆库区的水环境状况。EFDC水动力水质模型计算出各个排污口对水质控制断面的响应系数，以地表水Ⅱ类水标准为规划排放量，为重庆库区流域污染物合理削减量、水环境改善及水质保障打下理论基础，可为库区水环境规划管理提供决策依据。

目标，运用线性规划法求解各个控制单元的允许

（2）将重庆库区分为17个控制单元，基于

在运用线性优化法计算水环境容量过程中，时常出现某些排污口被“优化掉”的现象，即某些排污口的允许排放量为零，即对于水质响应系数大的排污口的允许排放量为零，未考虑各个排污口自身的经济、环境、技术发展状况，忽略了公平问题、效率问题，并不能作为最终控制单元-排污口总量分配方案，但可以结合库区经济、社会、技术等多种因素，为排污总量分配

16参考文献

三峡生态环境监测http://sxhc.cbpt.cnki.net/第4卷[16]李一平，邱利，唐春燕，等.湖泊水动力模型外部输入条

410-416.

件不确定性和敏感性分析[J].中国环境科学，2014，34（2）：

[1]文传浩，秦方鹏，王钰莹，等.从库区管理到流域治理：三

峡库区水环境管理的战略转变[J].西部论坛，2017，27（2）：

[17]黄庆超.三峡库区典型支流的水文水质特征研究[D].武

[2]58-62.

andFUB,WUB,LUY,etal.ThreeGorgesproject：[3]Geographychallenges单保庆，，2010for，34（the6environment[J].ProgressinPhysicalEfforts

方案制定方法及其案例研究王超，李叙勇，）等：741-7.

.基于水质目标管理的河流治理

[J].环境科学学报，2015，35（8）：

[4]2314-2323.

刘昭伟，的计算与分析陈永灿，[J].付健，水力发电学报，等.三峡水库万州段岸边环境容量

2006，25（4）：51-56.[5]交通大学，陈培帅.重庆主城区两江水环境容量研究2013.

[D].重庆：重庆[6]丁世敏，量控制研究熊竹，—封享华，—以重庆市涪陵区为例等.三峡库区城镇水环境容量及总

[J].三峡生态环境

监测，2016，1（1）：60-.

[7]王晓青，[J].环境科学研究，郭劲松.三峡水库蓄水后小江水环境容量的变化[8]2012，25（1）：36-42.

耦合模型研究王晓青.三峡库区澎溪河[D].重庆：重庆大学，（小江）富营养化及水动力水质

2012.

[9]区生态安全情景模拟及指标阈值确定张梦婕，官冬杰，苏维词.基于系统动力学的重庆三峡库

[J].生态学报，

[10]2015[D].费红，35重庆：.（三峡库区14）：4880-40.

重庆师范大学，（重庆段）2013.

农村面源污染空间分异研究[11]究卓海华，[J].长江流域资源与环境，吴云丽，刘旻璇，等2017.三峡水库水质变化趋势研

，26（6）：925-936.[12]重庆：张文时重庆大学，.基于EFDC2014.

模型的山地河流水动力水质模拟[D].[13]环境影响评价，张以飞，王玉琳，2015汪靓，37（.EFDC3）：70-72模型概述与应用分析，92.

[J].

[14]江辽，试验[J].尹小玲，水电能源科学，朱磊.潮汐河口概化感潮河段的数值模拟

2015，33（8）：9-12.

[15]水质王聪（.水系基于）GIS影响研究和RS[D].的三峡库区土地利用变化对支流

重庆：重庆大学，2017.汉：华中农业大学，2016.

[18]国环境科学出版社，孟伟.流域水污染物总量控制技术与示范2008.

[M].北京：中

[19]modelWUG,：CaseXUZ.studyPredictionintheofalgalbloomingusingEFDC

Modelling，2011，222（6）：1245-1252.

DaoxiangLake[J].Ecological

[20]刘夏明，拟中的应用及进展李俊清，豆小敏，[J].环境科学与技术，等.EFDC模型在河口水环境模

2011，34（S1）：

[21]136-140研究卫志宏，，360.

[J].生态科学，杨振祥，吕兴菊，2013，32等（3.）洱海动态水环境容量模拟

：282-2.

[22]回顾与展望董飞，，[J].彭文启，水科学进展，等.地表水水环境容量计算方法

2014，25（3）：451-463.[23]forZHAOF，LIC，CHENL，etal.Anintegrated[24]–reservoiraccountingcombinationforwatersystem[J].Waterenvironmentalcapacity，2018，ofthemethod

riverofHANChinanitrogenH，LI：ModelingandK,WANGphosphorusX,etpollutionsal.Environmental10（4）：483.

inJiaozhoucapacity

Bay,tin，2011，63（5-12and）：262-266.

assessing[J].MarinePollutionBulle⁃[25]ronmentalRANJBARMH,HADJIZADEHviacapacityofphosphorusZinN.GorganEstimationBayof，Iranenvi⁃

，[26]tal董一博Monitoringa3Decological-hydrodynamic.浑太河水环境容量计算研究andAssessment，2016model[J].，188[D].（11邯郸：）Environmen⁃：9.程大学，2015.

河北工[27]何羽究[D]..三峡库区主要污染物总量控制方案及安全边际研

重庆：重庆大学，2012.

[28]张卫东，控制技术研究邓春光，[M].王孟北京：.三峡库区流域水环境容量与总量

中国环境出版社，2013.[29]总量的控制分配卫志宏，唐雄飞，[J].杨振祥，湖泊科学，等.洱海主要污染物允许排放

2013，25（5）：665-673.