第八讲 SWAT非点源污染模拟

SWAT的非点源污染模拟主要内容1、非点源污染模拟原理2、模拟过程与关键步骤3、点源和非点源模拟应用实例一、非点源模拟原理模型概述流域中非点源污染的来源主要是农业化肥和农药的施用。氮、磷等营养物和农药的迁移和转化取决于化合物在土壤环境中经历的转换过程。SWAT模型可以模拟水文响应单元内氮和磷完整的营养物循环过程和任何一种农药的降解过程。模型概述SWAT模拟流域内泥沙、藻类、溶解氧、有机污染、多种不同形式的氮、磷以及农药等污染物质的运移与转换。物质循环模型可以分为非点源污染模块、河道水质模块以及湖泊水库水质模块。据统计，目前已在国际期刊上发表相关文献近200篇，中国知网收录文献292篇。流域污染物迁移转化过程示意图模型概述非点源污染模块河道水质模块湖泊、水库水质模块模型概述非点源污染模块农药模拟模块河道水质模块模型扩展应用作物生长模块河道演算模块模型非点源污染模拟原理物质运移内容：1、营养物，包括有机氮、铵氮、亚硝氮、硝氮、有机/无机磷；2、农药；3、重金属；4、病原菌；5、生化需氧量、溶解氧过程：与水循环过程类似，多数先在水文响应单元计算各种产出量（陆面过程），再输出到子流域的滞蓄水体以及河网系统进行循环过程氮素的土壤循环SWAT模型监测土壤中的五种不同氮库:无机氮库（NH4+和NO3-）有机氮库（新有机氮、活性有机氮和稳定有机氮）新有机氮（FreshorganicN）是指以农作物残余和微生物的形式存在的有机氮;活性有机氮和稳定有机氮是指土壤中以腐殖质形式存在两种不同类型的氮。土壤初始氮

腐殖质有机氮初始硝酸氮矿化作用、分解作用/固氮作用腐殖质矿化作物残余的分解和矿化Nminf,ly

=0.8·δntr,ly·orgNfrsh,ly硝化作用和氨挥发硝化挥发

新有机氮orgNfrsh,surf=0.0015·rsdsurf

反硝化作用如果γsw,ly≥0.95，则如果γsw,ly<0.95，则Ndenit,ly=0.0

Ndenit,ly

=NO3ly•(1-exp[-1.4γtmp,lyorgCly])降雨中的氮Nrain=0.01·RNO3·Rday

氮运移磷素的土壤循环化肥收割绿肥/垃圾/淤泥绿肥/垃圾/淤泥矿化土壤有机质吸附态无机磷固持解吸吸附矿物性土壤中磷主要以三种形式存在:腐殖质中的有机磷难溶的矿物磷土壤溶液中的植物可利用磷

土壤初始磷

稳定矿物磷活性矿物磷矿化作用、分解作用/固磷作用腐殖质矿化作物残余的分解和矿化无机磷的吸附作用溶液中迁移到活性矿物库活性矿物库中迁移到溶液腐殖质有机磷淋溶作用Pmina,ly=1.4•βmin•γtmp,lyγsw,ly•orgPact,lyPsol┃actly

=Psolutiontly-minPactly·Psol┃actly

=0.1·(Psolutiontly-minPactly·)磷运移农药的地表水体循环主河道中（仅模拟一种农药）：1、水中：降解、挥发、底泥沉降、出流；2、底泥中：降解、再悬浮、解吸附、掩埋；滞蓄水体中（仅水库，过程与主河道基本一致）农药的地表水体循环叶面施药挥发降解冲洗入渗淋溶挥发产流带走分解表层和亚表层施药淋洗作用：降解作用：pstf,wsh=frwsh·pstf

psts,ly,t=psts,ly,0·exp[-kp,sol·t]病原菌运移病原菌的来源：SWAT仅概念性地模拟两种菌群，一种为可长久生存的，一种为相对短命的。病原菌来自于绿肥(动物粪便），SWAT的肥料类型数据库中有各种绿肥的含菌量参数。在施肥过程中，病原菌也一同施入，并在植物叶面和表层10mm的土层中进行分配。另外可通过点源输入。病原菌的运移模拟：HRU的病原菌：1、植物叶面；2、土壤溶液；3、吸附在土壤颗粒中1、植物叶面的雨水冲刷；2、菌群的死亡-生长(3种介质）；3、土壤溶液中细菌的淋溶（淋溶的细菌认为死亡）；4、地表产流带走；5、在地表水体中（主河道、仅水库）随水流演进，并仅考虑死亡过程河道水质模型SWAT模型中的河流水质模型采用QUAL2E模型。该模型是以溶解氧为中心的多变量的综合性河流水质模型。它能按照使用者的要求，模拟河道中泥沙、叶绿素——藻类、有机氮、氨氮、有机磷、可溶性磷、亚硝酸盐、硝酸盐、矿物质磷、生化需氧量、溶解氧、农药、三种任选的重金属等不同水质指标的浓度变化过程。营养元素的地表水体循环主河道中：1、藻类的死亡->增加有机氮和有机磷；2、藻类的生长->氮磷固持（包括铵氮、硝氮、无机磷）；3、有机氮/有机磷的底泥沉降；3、有机氮的矿化及底泥吸附态铵氮的释放->铵氮；4、铵氮的硝化->亚硝氮->硝氮；5、有机磷的矿化及底泥吸附态无机磷的释放->无机磷滞蓄水体中（池塘、湿地、水库,不包括洼地）：1、随地表径流进入滞蓄水体，并出流进入主河道2、仅考虑营养元素的沉降，不模拟转化过程河道水质模拟市政和工厂排放非点源排放溶解物运移生物降解和转化颗粒运移稀释和扩散沉积作用沉积和再悬浮沉积和积累泥沙吸附作用河道水质模型1—复氧作用；2—河底生物(包括底泥)的耗氧；3—碳化合物BOD耗氧；4—光合作用产氧；5—氨氮氧化耗氧；6—亚硝酸氮氧化耗氧；7—碳化BOD的沉淀；8—浮游植物对硝酸氮的吸收；9—浮游植物对磷(磷酸盐磷)；10—浮游植物呼吸产生磷；11—浮游植物的死亡和沉淀；12—浮游植物呼吸产生氨氮；13—底泥释放氨氮；14—氨氮转化为亚硝酸氮；15—亚硝酸氮转化为硝酸氮；16—底泥释放磷。藻类模拟藻类（或叶绿素a）的生长和分解可以通过与生长率、呼吸速率、沉降率以及当前河流中藻体数量建立函数计算得出。藻类模拟呼吸或死亡速率：沉降速率：生长速率：光氮磷氮元素

在含氧水体中氮是逐级转化的，从有机氮转成氨氮再转为氮，亚硝酸态氮，最后成为硝酸氮。有机氮、氨氮、亚硝酸氮和硝酸氮通过沉降和底泥吸附从河道中除去。氮元素有机氮：氨氮：亚硝酸态氮：硝态氮：磷元素

磷循环和氮循环相似。藻体死亡后，藻体内的磷转化为有机磷。有机磷再经过矿化成为可被藻类吸收的溶解态磷。同时，有机磷也可能会通过沉降作用从河流中去除。有机磷：无机（可溶性）磷生化需氧量及溶解氧生化需氧量(CBOD)：指水体中有机物质分解所需要的氧气的总量。由点源输入到河网系统，在主河道和水库中循环。1、CBOD的底泥沉降；2、CBOD的氧化分解（耗氧）生化需氧量及溶解氧

溶解氧(DO)：水体中的溶解氧浓度和大气复氧、光合作用、动植物呼吸作用、底泥需氧量、生化需氧量、硝化作用、盐度、温度等相关。1、复氧过程。大气复氧、藻类的光合作用、水体扰动复氧(坝体)2、消耗过程。藻类的呼吸作用、CBOD耗氧、底泥耗氧、铵氮硝化、亚硝氮硝化重金属的运移

允许模拟3种SWAT的重金属来源于点源输入。重金属运移是SWAT的物质运移模拟中唯一一个不考虑HRU因素的物质运移，仅是简单地计算随河网系统水流演进过程和物质量平衡。蓄水体水质模拟SWAT湖泊、水库水质模型采用了一个简单的污染负荷平衡模型。同时在模拟较为复杂的湖泊时，SWAT采用了分布式WASP水质模型。蓄水体水质模拟

模拟水体中营养物质的转化过程时，SWAT假设在一个完全均匀的水体，当营养物质进入水体的同时，就会均匀分布于整个水体。这个假设忽略了水体中的分层现象以及由于水体表层浮游生物引起的表层和下层之间的差异的存在。水体初始水量：水体初始氮磷含量：沉积作用导致水体中营养物质的损失量等于区域内沉积物－水界面通量和面积的乘积。在稳态情形下，水库污染负荷总平衡方程为：二、建模过程与关键步骤——以AVSWAT2005模型自带数据库为例影响因素DEM精度、降水不均匀性、子流域划分输入参数：土壤物化性质、植被生理生态水资源开发利用：水利工程、工农业取用水其他因素：对流域的熟悉程度模拟：径流模拟：基础泥沙模拟：迁移机理污染物模拟：与径流和泥沙的关系主要步骤主要步骤数据准备矢量图：DEM、土地利用、土壤图、气象水文站点位置等文本：土地利用和土壤数据库链接表表格：气象水文站点（包含高程、经纬度信息）输入气象资料表：各站点日最高最低气温、日降水量、日蒸发蒸腾、日风速露点温度等其他：取用水、引水地点等校核资料：控制水文站日径流量、水库调节方式、实测污染物负荷数据准备点排放数据表－年度负荷点源或者入口排放数据可以用四种方法进行总结：稳定日负荷，平均年负荷，平均月负荷，日负荷。用稳定日负荷对排放数据进行总结，则要在点排放数据对话框中输入数据。对其他三种方法，界面要求预先生成包含点排放数据的文件。Name定义Name定义YEAR测量数据的年份NH3YR年日均氨负荷(kg/day)FLOYR年日均水负荷量(m3/day)NO2YR年日均亚硝酸盐负荷(kg/day)SEDYR年日均泥沙负荷(kg/day)CMTL1YR年日均保守金属＃1负荷(kg/day)ORGNYR年日均有机氮负荷(kg/day)CMTL2YR年日均保守金属＃2负荷(kg/day)ORGPYR年日均有机磷负荷(kg/day)CMTL3YR年日均保守金属＃3负荷(kg/day)NO3YR年日均硝酸盐负荷(kg/day)BACTPYR年日均持久细菌负荷(#bact/100ml)MINPYR年日均可溶性磷负荷(kg/day)BACTLPYR年日均非持久细菌负荷(#bact/100ml)主要步骤开始界面主要步骤流域描述主界面主要步骤流域描述主要步骤添加土地利用和土壤数据-ArcSWAT中有坡度划分主要步骤添加气象数据-划分子流域降水和气象数据应在相应目录下主要步骤写入数据主要步骤设定参数-运行模型主要步骤输出结果查看手动修改输入参数数据库土壤气象土地利用杀虫剂耕作城市点源污染排放入口排放水库子流域点源污染排放排放情景选择与输入稳定日负荷点源污染排放排放情景选择与输入年均日负荷月均日负荷日负荷流程相似入口排放修改稳定日负荷年均日负荷月均日负荷日负荷与点源污染输入方式类似水库输入数据修改水库输入数据修改泄流方式选择年均泄流速度日水库出流月均水库出流模拟-设定目标泄流量水库输入数据修改水库水质数据子流域数据修改子流域气象站点常规数据水文响应单元地下水河道汇流用水……子流域数据修改气象站点子流域数据修改常规输入数据修改水文相应单元子流域数据修改主河道汇流地下水用水子流域数据修改农业管理措施一般管理信息水文参数-CN值子流域数据修改农业管理措施特定操作必需以日或者植物生长阶段（热量单元）来安排种植开始、灌溉、施肥、杀虫剂子流域数据修改农业管理措施耕作收获和割除收获子流域数据修改农业管理措施放牧自动灌溉自动施肥子流域数据修改农业管理措施土地清理圈围和释放长期施肥选择数据适用区域子流域数据修改土壤化学参数河流水质子流域数据修改水库坑塘输出结果查看结果并作图输出结果子流域1在1977年的产水量保存工程回到主界面保存模拟步骤重点水量、泥沙模拟精度高点源污染输入正确：位置、量农业管理措施：施肥量、日期（生长期）灌溉耕作土壤化学参数准确三、应用实例新安江流域非点源污染模拟淮河流域闸坝调控下水量水质模拟温榆河流域水量水质联合调度68在中国干旱半干旱、湿润地区等10多个小、中、大流域尺度得到了应用。主要应用于径流模拟、非点源污染模拟、闸坝评估及调度、以及气候变化等方面。新安江流域非点源污染模拟69数据收集新安江流域非点源污染模拟数据类型尺度数据属性DEM网格:30m×30m高程,坡度,坡长土地利用网格:30m×30m土地利用分类土壤类型1:1,000,000土壤分类及其理化数据雨量站48个(2000-2010年)日降水量气象站3个(2000-2010年)日最高最低气温、风速、相对湿度、日照时数水文站6个(2000-2010年)月径流过程水质断面3个(2000-2010年)月泥沙、总磷、总氮等水质指标浓度社会经济13个县级行政区化肥使用量(折纯)、牲畜存栏数、人口等数据收集71参数选择及率定-水质①地表产流过程影响泥沙、N和P负荷的参数：SPCON：泥沙输移线性系数；SPEXP：泥沙输移指数系数；USLE_K：土壤侵蚀力因子ERORGN：氮在土壤中的富集率；ERORGP：磷在土壤中的富集率；NPERCO：氮的下渗系数，地表径流中氮的浓度与下渗水流中氮的浓度比；PPERCO：磷的下渗系数，地表径流中氮的浓度与下渗水流中氮的浓度比，；PHOSKD：磷的土壤分离系数，指表层10mm土壤中可溶性磷的浓度和地表径流中可溶性磷的浓度之比；②影响氮、磷在河道水质模拟参数：在水体中有机氮可以转化到氨氮、亚硝酸盐氮然后到硝酸盐氮。磷循环过程与氮循环有相似的转化关系，故考虑如下几个参数的率定：BC3：20˚C时有机氮转化为氨氮的速度常数；BC1：20˚C时氨氮生物氧化速度常数；BC2：20˚C时亚硝酸盐氮的生物氧化速度常数；BC4：20˚C时有机磷的矿化速度常数；RK1：20°C时CBOD氧化速率RK3：20°C时CBOD的沉降速率细菌分配系数Mgt和fert.dat文件的制作编号矿物质氮的比例矿物质磷的比例有机磷的比例有机氮的比例氨氮的比例稳定细菌比例弱稳定细菌比例名称MgtFert.dat73流量泥沙总氮总磷屯溪模拟结果74流量泥沙总氮总磷渔梁模拟结果75模拟指标时段评价指标三阳坑临溪新亭月潭渔梁屯溪水量率定期体积误差-0.020.00-0.09-0.090.02-0.02相关系数0.950.960.950.950.930.95效率系数0.910.920.900.890.870.91检验期体积误差0.010.010.090.080.070.07相关系数0.970.950.910.970.960.97效率系数0.920.890.810.930.880.92泥沙率定期体积误差-----0.09-0.06相关系数----0.720.85效率系数----0.500.60检验期体积误差-----0.050.09相关系数----0.900.89效率系数----0.750.76总氮率定期体积误差----0.03-0.06相关系数----0.670.50检验期体积误差----0.000.03相关系数----0.380.30总磷率定期体积误差-----0.060.02相关系数----0.850.74检验期体积误差----0.100.14相关系数----0.530.80新安江非点源污染估算模型模拟结果76降水-径流空间分布图77污染负荷空间分布78泥沙总氮总磷行政区县各年产污量淮河流域闸坝调控下水量水质模拟嵌入闸坝群运行的流域水循环分布式SWAT模型

淮河流域是我国人口最密集、闸坝最多和污染最为严重的流域。淮河流域目前已兴建约1.1万余座闸坝,占全国的一半以上。闸坝在防洪、灌溉供水等方面发挥巨大效益。但工程修建后对河流生态与环境造成的影响，已成为国内外关注的焦点。淮河流域淮河流域概况（1）SWAT模型与闸坝调度方式的耦合 SWAT模型中水库水平衡方程:SWAT模型??SWAT模型所提供的四种模拟水库：出流方式中实测日/月出流量需要输入模拟期间的所有水库日/月观测径流数据；目标泄水控制出流需要调查获得每个水库每个月的目标蓄水量以及最小/最大下泄流量这三种模式方式都需要非常详细的闸坝运行资料；无控制年平均放水率只需要主溢洪道和紧急溢洪道水位对应的蓄水量和水面面积，并没有包含闸坝的调度，将直接导致径流过程模拟的误差较大。?（1）SWAT模型与闸坝调度方式的耦合

在我国闸坝的设计时，闸坝的水位、库容和下泄流量存在对应关系。闸坝的运行调度时，闸坝的下泄流量主要依据闸上水体达到某一水位（或库容）。因此三者之间存在一定的相关关系：白龟山水库：白龟山水库的泄流量、坝上蓄水量、库容面积及水位之间的相关关系

（1）SWAT模型与闸坝调度方式的耦合 淮河流域重点水库和水闸泄流量、坝上蓄水量及水位之间的相关关系

相关性非常显著（2）SWAT中水质模块的改进

（1）在水库水质模块中，采用污染负荷平衡方程考虑BOD的沉淀、降解，添加模型中水库BOD模拟功能。（2）建立COD与BOD之间的关系。从COD和BOD氧化机理上分析，两者存在一定的线性关系。（3）修改模型中水库、河道和子流域输出程序，以浓度形式输出不同水质指标的变化过程。

研究选取蚌埠闸为流域出口，建立了闸坝运行环境下SWAT模型，共划分了129个子流域和468个水文响应单元。主要大型闸坝按其所在地理位置作为相应子流域的出口控制断面加入模型中。（3）模型的建立 淮河流域闸坝、河网、子流域耦合分布图子流域：129hru：468闸坝：31参数的选择和率定1963-2000年（率定期1963-1990年，检验期1991-2000年）40个闸控断面和水文断面水量模拟：水质模拟：1984-2000年期间不连续月浓度资料。模拟负荷为氨氮和CODMn。36个断面敏感性分析指标：类别指标敏感性I0.00≤|I|<0.05可以忽略II0.05≤|I|<0.20适中III0.20≤|I|<1.00敏感IV|I|≥1.00非常敏感参数的选择和率定参数河流水库径流NH3-NCODMn径流NH3-NCODMnⅠClassⅠClassⅠClassⅠClassⅠClassⅠClass水量CN0.32Ⅲ(p)0.44Ⅲ(p)1.07Ⅳ(p)0.18Ⅱ(p)0.22Ⅲ(p)0.61Ⅲ(p)ESCO0.34Ⅲ(p)0.19Ⅱ(p)0.32Ⅲ(p)0.25Ⅲ(p)0.82Ⅲ(p)0.10Ⅱ(p)EPCO-0.05Ⅱ(n)-0.04Ⅰ(n)-0.03Ⅰ(n)-0.03Ⅰ(n)-0.01Ⅰ(n)0.02Ⅰ(p)Sol_AW-0.45Ⅲ(n)-0.27Ⅲ(n)-1.42Ⅳ(n)-0.32Ⅲ(n)-0.15Ⅱ(n)-0.44Ⅲ(n)Sol_k0.10Ⅱ(p)-0.01Ⅰ(n)-0.40Ⅲ(n)0.06Ⅱ(p)1.92Ⅳ(p)-0.06Ⅱ(n)Sol_z-0.89Ⅲ(n)-0.22Ⅲ(n)-1.14Ⅳ(n)-0.44Ⅲ(n)-0.11Ⅱ(n)-0.54Ⅲ(n)水质NSETLR-NS-0.11Ⅱ(n)-NS-NS-1.13Ⅳ(n)-NSBSETLR-NS-NS-0.19Ⅱ(n)-NS-NS-0.64Ⅲ(n)Chla-NS0.00Ⅰ(p)0.00Ⅰ(p)-NS0.00Ⅰ(p)0.00Ⅰ(p)RK1-NS-NS-4.02Ⅳ(n)-NS-NS-2.73Ⅳ(n)RK3-NS-NS-1.61Ⅳ(n)-NS-NS-0.87Ⅲ(n)BC1-NS-0.88Ⅲ(n)-NS-NS-1.55Ⅳ(n)-NSBC3-NS0.28Ⅲ(p)-NS-NS0.20Ⅲ(p)-NSCtoB-NS-NS0.79Ⅲ(p)-NS-NS0.93Ⅲ(p)ERORGN-NS0.97Ⅲ(p)-NS-NS0.09Ⅰ(p)-NSErcbod-NS-NS0.76Ⅲ(p)-NS-NS0.41Ⅲ(p)（4）模拟结果

水库模拟白龟山水库12座水库模拟结果（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）水质BSETLRCOD在水库中的沉降速度.lwq（-100，100）NSETLR氮在水库的沉降速度.lwq（-500，5.5）RK120°C时CBOD氧化速率.swq（0.0，3.4）RK320°C时CBOD的沉降速率.swq（-0.10，0.36）BC120°C时氨氮的氧化速率.swq（0.05，0.10）BC320°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.570.4053%30%鲁疃闸---33%56%辛堡闸---34%70%苇沟闸---35%44%北关闸-5%0.860.8415%69%断面多年平均不同保证率来水量2004-2005年年平均径流模拟P=50%P=75%P=95%沙河闸78306290453036007348鲁疃闸13160108708260688013244北关闸5215048320412303549062913水量模拟：沙河闸和北关闸的年径流总量都接近于断面多年平均径流量；北关闸模拟相关系数和效率系数都在0.80以上。水质模拟：沙河闸、鲁疃闸、辛堡闸、苇沟闸和北关闸的断面氨氮模拟平均相对误差为34%，COD的平均相对误差为54%。水量水质模拟结果基于耦合闸坝系统的温榆河流域SWAT模型，通过沙河闸、鲁疃闸、新堡闸、苇沟闸、北关闸等多闸坝的调控，进行水污染修复与整治，使温榆河干流水体水质浓度尽量达到所在水功能区划的浓度标准，从而提高流域内水资源的可利用量，实现流域水资源的高效利用。闸坝优化调度模型与多情景分析目标函数：

当流域污染严重时，下游水体的等标污染指数最小库容约束闸坝下泄流量约束水量平衡约束约束条件：模型求解采用遗传算法该指数数值小于等于1，说明该断面水体没有受到污染；数值大于1，说明该断面水体已受到污染,数值越大，超标越严重。

当流域污染较轻时，总达标水资源量最