SWAT模型在我国山区小流域水质预测中的应用——以琼江流域为例.pdf

第42 卷第5期2023年10 月水环境SWAT模型在我国山区小流域水质预测中的应用-以琼江流域为例一四川环境SICHUAN ENVIRONMENTD0I:10.14034/ki.schj.2023.05.020Vol.42,No.5October 2023胡鑫，刘兵，杨永安”，王远铭4，李克锋4（1.四川省环境应急与事故调查中心，成都6 10 0 0 0；2.遂宁市生态环境安全应急中心，四川遂宁6 2 9 0 0 0；3.西华师范大学环境科学与工程学院，四川南充6 37 0 0 2；4.四川大学水力学与山区河流开发保护国家重点实验室，成都6 10 0 6 5）摘要：小流域水环境受人类活动的影响显著，为解析我国典型山区小流域的水质受人类活动的影响规律，选取琼江流域遂宁段为研究对象，构建了针对琼江流域的SWAT模型，在模型验证的基础上并预测了2 0 30 年社会经济发展条件与相应治污条件下琼江流域遂宁段的水质时空变化特征。研究结果表明，2 0 30 年琼江流域大安断面逐月流量变化范围为2.8240.44m/s，年均流量为16.2 0 m/s，其中8 月流量最大；各支流中蟠龙河的流量最大，年均流量达到1.13m/s；2 0 30 年琼江干流NH,-N与TP浓度较高，在枯水季节尤其1 2 月水质较差，难以满足I类水质标准；干流控制安居城区以上的高滩断面除1月的水质能达到类水质要求，流经安居城区后，干流萝卜园断面在2 月份TP未能达到类水质要求，至安居城区以下大安断面，水质进一步变差，13月份难以满足类水质要求。成果可为我国山区小流域水质的研究提供借鉴，同时为琼江流域水环境保护和政府决策，提供科学依据和技术支持。关键词：小流域；SWAT模型；水质预测；时空变化中图分类号：X522Application of SWAT Model in Water Quality Prediction of Small Watershed in theMountainous Area of China:A Case Study of Qiongjiang River BasinHU Xin,LIU Bing,YANG Yong-an,WANG Yuan-ming,LI Ke-feng(1.Center of Environmental Emergency&Accident Investigation of Sichuan Province,Chengdu 610000,China,2.Emergency Center for Ecological Environmental Safety of Suining,Suining,Sichuan 629000,China;3.School of Environmental Science&Engineering,China West Normal University,Nanchong,Sichuan 637002;4.State Key Laboratory of Hydraulics&Mountain River Engineering,Sichuan University,Chengdu 610065,China)Abstract:The water environment of small river basin is always significantly affected by human activities.In order to analyzethe law of water quality affected by human activities in typical small mountain basins in China,this paper selected Suining sectionof Qiongjiang River Basin as the research object to establish a calibrated SWAT model.Based on the model verification,thetemporal and spatial changes of water quality in the Suining section of Qingjiang River basin under the social and economicdevelopment conditions and corresponding pollution control conditions in the year 2030 were predicted in this study.The resultsshowed that the monthly water flow at the Da-an section ranges from 2.84 40.44 m/s in the year 2030,and the annual waterflow is 16.20 m/s,with the largest flow in August.Among the tributaries,Panlong River has the largest annual discharge,withan average value of 1.13 m/s.The concentrations of NH,-N and TP in the mainstream are relatively high in 2030,and the waterquality is poor in the dry season,especially January and February,which is difficult to achieve the water quality standard of classII.The Gaotan section controls the water quality upstream of the urban area,and it has a relatively good water quality except inJanuary.In February,the water quality can not reach the standard of class IIl after flowing to the Luoboyuan section.Waterquality at the Da-an section further deteriorates,and it is dfficult to meet the class II water quality requirements from January toMarch.The results of this study can provide reference for the research on the water quality in small mountain river basins inChina,and provide scientific basis and technique support for water environment protection and government decision-making of theQiongjiang River.Keywords:Small river basin;SWAT model;water quality prediction;temporal and spatial variation.收稿日期：2 0 2 2-0 6-10基金项目：四川省省级科技计划项目（2 0 18 SZYZF0001）。作者简介：胡鑫（19 8 5）男，四川乐至人，毕业于西南大学环境科学专业，工程师，主要从事环境风险防控与预警研究工作。通讯作者：杨永安，。文献标识码：A文章编号：10 0 1-36 44(2 0 2 3)0 5-0 143-10144前言小流域水环境受人类活动的影响显著。长江上游山区小流域经济发展水平低，流域内面源污染占据主导地位且时空分布不均；此外，我国长江上游地区经济发展迅速，未来区域用地结构、产业结构调整变化大，当前水污染管控尚处于较低水平，增加了山区小流域水质不达标的风险，为长江上游山区绿色小流域的构建带来了挑战。当前，针对山区小流域水环境的研究集中于流域山洪灾害、流域水文模拟及流域水土保持等，对山区小流域水质的分析及预测较为缺乏 2-3。琼江是我国西南地区典型的流程短、流量小、面源广、沿程工业源少的小流域，具有很强的代表性。近年来琼江流域沿岸城镇及农村社会经济快速发展的同时，也给琼江带来了巨大的生态环境压力 4。琼江流域遂宁市内的大安断面，近年来水质不能稳定达到类水质标准，被列为四川省挂牌整治的小流域达标考核断面 5。目前关于琼江流域的水质研究集中于关键水质指标的监测与评估(6-7 1，或是流域点、面源污染负荷的统计与复核，针对琼江流域特征的水质预测模型至今没有建立，流域水质随着未来社会经济发展的时空变化规律尚未被揭示，已成为琼江流域水环境保护的重要制约因素之一。分布式水文模型（Soil and Water AssessmentTool，简称SWAT）是美国农业部农业研究中心研发的一个基于GIS平台的分布式流域水文模型，采用日为时间步长进行连续模拟计算，能够模拟和预测流域内气候、土地利用方式、土壤类型和农业管四川环境理措施等对径流量、泥沙和化学污染物质输出的影响 8-9)。由于SWAT模型具有适用性好、模型性能稳健、运算效率高、连续模拟时间长等优点，在模拟流域内营养物质输移和时空演变 10 ，湿地对污染物的降解能力，农业非点源污染空间分布 12 14，流域内氮、磷等污染物流失 15-16 ，农药等污染物的运移与转换等方面表现出了良好的模拟效果 17 。本文选取琼江流域遂宁段为研究对象，通过构建针对琼江流域的SWAT模型，预测了社会经济发展条件下琼江流域遂宁段的水质时空变化特征。本文的工作拓展了我国山区小流域的水环境研究，研究结果可为琼江流域水环境保护和政府决策，提供科学依据和技术支持。1材料与方法1.1研究区域琼江是江右岸一级支流，源于四川省资阳市乐至县三星乡金马寺，经该县的棺材函后流人遂宁市安居区，于大安乡陈家坝下半边山出安居境，进人重庆市潼南区，至铜梁县安居镇汇人江。琼江全长2 35km，全流域面积456 0 km。市境内流程117km，流域面积10 0 8.4km，天然落差40 m，河道平均比降0.34%o，琼江流域遂宁段水系分布如图1所示。琼江主要流经的安居区下辖2 个街道和21个乡镇。2 0 17 年末，全区总人口为7 8.36 万，其中农业人口7 1.0 9 万，占9 0.7 3%，城镇人口19.2万人，城镇化率为2 9.9 1%。42卷四川省車庆市大安断面S0FIOUrN图例资阳诚城点源入510微城出口20干米控制点105-400数字代表划分的子流域，星号代表流域点源图1琼江流域遂宁段水系图Fig.1Water system map of Suining section of Qiongjiang River Basin5期1.2研究方法基于流域数字高程地图、土地利用类型、土壤类型和土地坡度，本文首先构建了针对琼江流域的SWAT模型。根据实测气象数据与区域统计年鉴，分别完成SWAT模型的气象数据库构建与管理措施方案输入，将琼江实测流量数据、水质数据及统计的污染源数据作为SWAT模型边界条件的输人。琼江实测流量数据与水质数据也将用于SWAT模型的参数率定与验证。在完成模型的构建、参数的率定与验证后，采用建立的SWAT模型，根据琼江流域已有社会经济发展规划对琼江流域2 0 30 年河道水文、水质过程进行模拟预测，解析2 0 30 年社会经济发展条件下琼江流域遂宁段的水质时空变化特征。2结果与讨论2.1模型的构建与验证2.1.1模型建立SWAT模型输人的空间数据包括数字高程地图胡鑫等：SWAT模型在我国山区小流域水质预测中的应用一一以琼江流域为例145（D EM）、土地利用类型和土壤类型。DEM采用ASTERGDEM提取的30 m分辨率的数字高程数据图，土地利用数据采用遂宁市安居区精度为30 m的2 0 19 年土地利用类型栅格文件。将收集到的土地利用信息转化为SWAT模型能够识别的类型后，进行土地利用重分类（图2 b），琼江流域遂宁段内，耕地占比最多，为59.43%；林地其次，占比为2 0.37%；水田第三，占比为15.52%；水域占比为1.6 6%；其他建设用地占比1.40%；城镇用地占比1.18%；裸土地占比1.0 0%；农村居民点占比0.40%。土壤数据采用精度为1:10 0 万的HWSD-V1.2世界土壤数据栅格文件，建立琼江流域遂宁段土壤类型数据库。研究区域土壤主要包括石灰冲积土和人为土，面积超过10%的土壤类型如下图2 c所示。琼江流域遂宁段为典型的丘陵地貌，坡度类型复杂，坡度为0 10、10 2 0、和2 0 以上的土地分别占2 8.16%、38.8 6%和32.9 8%，其中坡度为10 2 0 的土地占比最高（图2 d）。大英养(b)乐图例重庆市(c)图创20+*(d)图例图2 SWAT模型构建琼江流域遂宁段（a）数字高程、（b）土地利用、（c）土壤类型与（d）土地坡度Fig.2 The(a)DEM,(b)land use,(c)soil category and(d)land slopeof Suining section of Qiongjiang River based on the SWAT model气象数据采用2 0 112 0 18 年的实测值（图3）用于模型参数率定，并结合气候预测再分析数据模型（Climate Forecast System Reanalysis，C FSR）中的未来气候数据构建分布式SWAT模型的气象数据库，包括降雨、最高和最低气温、太阳辐射、相对湿度和风速。依据遂宁市安居区2 0 17 2 0 18 年统计年鉴，琼江流域遂宁段内主要粮食作物为春小麦、大麦、玉米、水稻等，主要经济作物为油料和棉花等，将上述作物所需肥料用量按照耕作时间输人模型，完成管理措施方案。146四川环境日110日1042卷日10400.00350.00300.00(uu)鲁250.00200.00150.00100.0050.000.004540353025(.）120151050-5降雨量(mm)一最低气温()一最高气温()日19 10日1$8 10350.0%300.0250.0200.0150.0100.050.00.0加载DEM，划分子流域，设置流域入口、子流域出口和流域总出口。考虑流域内水文地质、行政区域和重点监测断面等情况，将大安断面设置为子流域出口，其余子流域出口由模型自动生成节点，大安断面在流域末端和四川省边界，设置为流域总出口。共划分30 个子流域，见图1。琼江遂宁段干流控制的子流域从上游到下游依次为：17、9、5、2、3、8、11、15、2 0、2 2、2 5、2 6、2 7、28、30；琼江遂宁段支流盐井河控制的子流域为18，莲花河控制的子流域为6，中心场河控制的子流域为10，双龙河控制的子流域为1，白马河控制6.005.004.00(s/叫)单区3.002.001.000.00太阳辐射(W/)相对湿度(%)图3琼江流域2 0 11 2 0 18 年气象数据实测值Fig.3 Meteorological data measured in Qiongjiang River basin from 2011 to 2018风速(m/s)的子流域为4，观音河控制的子流域为7，石洞河控制的子流域为1、13和16，胡家河控制的子流域为14，蟠龙河控制的子流域为2 4，会龙河和玉丰河控制的子流域为19、2 1和2 3。研究收集了大安水文站点的流量和水质数据（NH,-N、T P）的逐月值，污染源数据用于边界条件的输人，主要包括点源污染（城镇生活污水、工业源及规模化畜禽养殖）和面源污染（农村生活污水、散养畜禽和农业种植）。根据琼江流域污染负荷现场调查与计算，各子流域点源污染统计见表1。非点源污染负荷统计见表2。5期表1琼江流域遂宁段子流域点源污染统计Tab.1 The point source pollution at different sub-areasof Suining section of Qiongjiang River子流域名称及编号琼江干流2 0胡家河14蟠龙河2 4石洞河16琼江干流15琼江干流2 7玉丰河2 1琼江干流9盐井河18琼江干流2中心场河10会龙河2 1琼江干流30琼江干流2 9莲花河6观音河7石洞河12.1.2模型参数率定与验证利用收集到的水文和水质多年实测数据，通过SWAT-CUP率定SWAT模型参数。将多年实测数据划分为模型需要的不同时期，其中，流量参数的预热期（即初始计算时期）为2 0 0 8 2 0 11年，率定期为2 0 142 0 16 年，验证期为2 0 17 年，水质参Tab.3 The result of parameter calibration in Suining section of the Qiongjiang River basin参数选代方式CH_K2.IteVCN2.mgtR湿润条件下的初始SCS径流曲线值-12.38 8ALPHA_BF.gWVREVAPMN.gWVEPCO.hruVGW_SPYLD.gWBIOMX.mgtHRU_SLP.hruSOL_BD().sol SOL_AWC().solCH_N2.rteSOL_Z().sol SURLAG.bsnGW_REVAP.gWGWHT.gWCANMX.hruANION_EXCL.sol胡鑫等：SWAT模型在我国山区小流域水质预测中的应用以琼江流域为例NH;-N198.0397.63025.93912.63717.62441.37518.54636.85610.99125.6468.69522.54825.5945.2417.3556.8026.023主河道河床有效水力传导率基流系数浅层地下水再蒸发水深阈值植物吸收补偿因子浅层含水层产流量比V生物混合效应R平均坡度R土壤饱和容重R土壤层有效水容重V主河道河床的曼宁系数R土壤表层到底层的深度V地表径流滞后时间V地下水再蒸发系数V地下水初始高度V最大冠层截留量V阴离子孔隙率147表2 王琼江流域遂宁段子流域非点源污染统计Tab.2 Statistics of non-point source pollution in(kg/d)Suining section of Qiongjiang River BasinTP子流域编号43.107盐井河182.010石洞河16.245会龙河192.421中心场河104.137琼江干流2 014.870琼江干流23.450琼江干流306.861蟠龙河2 42.883石洞河163.460观音河71.892胡家河143.362玉丰河2 110.707琼江干流91.000莲花河62.720琼江干流2 91.458琼江干流2 71.662数的预热期为2 0 0 8 2 0 11年，率定期为2 0 17 年，验证期为2 0 18 年。根据SWAT输入/输出手册和文献调研及局部敏感性分析选择参数，利用SWATC U P提供的SUFI-2算法率定参数并分析灵敏度，其中地表径流、基流、地下水、坡度、蒸散发、土壤含水等相关参数较敏感，各参数率定成果见表3。表3琼江流域遂宁段主要参数率定成果物理意义值P值35.7230.0000.0003.0670.002-2.9830.0032.7150.0072.4120.0172.2050.029-1.7710.078-1.4470.1501.2380.2171.1230.263-1.0890.2781.0730.2850.8300.4080.5180.6050.4750.6350.3570.721(kg/ha)素磷素27.4663.52352.9087.954286.92580.144179.60814.963231.55567.846231.65527.606536.03660.292219.69263.033131.51611.969197.10461.414307.61844.415239.24130.569171.45668.106160.51944.839151.83032.497104.96342.516初始范围最终值-0.01 500351.996735-0.5 0.80.540517010.2953810500436.227570010.0109610 0.40.049728010.755239-0.1 0.10.049782-0.1 0.10.000298-0.1 0.1-0.053216-0.01 0.30.108397-0.1 0.1-0.0635000.05 2416.0279370.02 0.20.1258810250.7663870 1000.0179470.01 10.292567148续表3参数OV_N.hruSOL_K().solGW_DELAY.gWGWQMN.gWSOL_CRK.sol四川环境选代方式物理意义V地表径流曼宁系数R饱和导水率V地下水迟滞时间V浅层含水层产生基流的阈值深度V土壤裂缝体积势42卷值P值0.2970.7670.2600.7950.1290.8980.0760.9400.0550.956初始范围0.01 10.2 0.202000 500001最终值0.0222990.1337570.0497283157.5903320.625978采用纳什效率系数（NS）和决定系数（R）两个指标评价实测数据和最佳模拟结果的拟合度，认为NS0.5且R大于0.5时模型精度可以接受。NS和R?公式如下：_ E1-1(Qm.i-Q.)2NS=1Z1(Qm.;-Qm)2Z1-（Q m;-Q m)（Q s,;-Q.)?R=Z=l(Qmi-Qm)?Zi=1(Qsi-Q.)式中：Qm为实测数据，m/s；Q m,是第i个实测数据，m/s；Q s 是第i个模拟数据，m/s；Q m是实测数据平均值，m/s；n 是实测数据量。1201-100806040200大安断面率定期和验证期NS分别为0.7 3和0.62，决定系数R分别为0.7 2 和0.6 7，根据模拟结果评价标准，率定效果良好，因此率定所得流量参数满足本流域研究需要（图4）。考虑到大安断面作为琼江流域遂宁段控制断面，同时也是琼江(1)水质在遂宁境内的重要考核断面，模型对大安断面的径流量与水质进行了预测，结果表明2 0 14年1(2)月2 0 17 年12 月的径流量预测值能实测值较好地与实测值吻合，模型2 0 18 年水质预测的结果表明，除2 月外，其余各月份TP与NH,-N值都能较好地与实测值吻合（图5）。月平均降雨量一观测值一模拟值E+9102015年1月051015202530354045501$L10时间图4大安断面流量率定及验证效果图Fig.4The calibration and validation of water flow at Da-an section2.0(T/Au)dL1.51.00.5011/261/152.5(T/u)N-HN2.01.5。观测值模拟值Od3/64/256/148/39/2211/11 12/31时间图5大安断面水质率定及验证效果图Fig.5 The calibration and validation of water flow at Da-an section一模拟值观测值1.0O0.5011/2651/153/64/256/148/39/2211/1112/31时间O9O5期2.2琼江流域遂宁段水质预测2.2.1预测因子及工况采用本文建立的分布式水文模型SWAT对琼江流域河道水文、水质过程进行模拟预测。根据遂宁市针对琼江流域已有社会经济发展规划 18 ，预测水平年选择2 0 30 年，预测工况考虑基于模拟的2030年流量条件下的相关治污规划水平、以及三个控制节点（高滩、萝卜园、大安）断面在已有治污规划水平下的水质年内变化过程。至2 0 30 年，安居区区域常住人口将达到7 7.2 0 万人，城镇化水平为56.9%，中心城区人均综合用水量指标为450升（人天），其他镇区为40 0 升（人天）。2030年安居城区将分片建设三个污水处理厂，即玉丰河污水处理厂（设计规模1m/d）、乌木厅污水处理厂（设计规模8 万m/d）与龙眼井污水处理厂（设计规模2.8 万m/d），其余乡镇各设一座胡鑫等：SWAT模型在我国山区小流域水质预测中的应用一一以琼江流域为例149污水处理厂。中心城区污水处理厂出水水质达到城镇污水处理厂污染物排放标准（CB189182002）的一级A标准，各乡镇根据受纳水体环境功能要求执行一级B标准。在已有规划治污水平下，预测年各乡镇污水厂污水设计处理能力将大于污水排放量。根据对琼江流域的水质分析，预测因子选择 NH,-N、T P。2.2.2河道流量根据流域产汇流过程，干流范围从上游至下游流量呈现增加趋势，大安断面处流量达到最大值，年均流量为16.2 0 m/s。各支流中蟠龙河的流量最大，年均流量达到1.13m/s。从时间角度来看，大安断面逐月流量变化范围为2.8 2 40.44m/s，其中8 月流量最大。河道流量在2 0 30 年每个季度代表月沿程变化如图6 所示。2121控制点Q_FEB_cms0.069590-0.2566000.256601-0.4905000.490501-1.1880001.188001-1.7390001.739001-4.017000控制点Q_MAY_cmS300.336000-0.5824000.582401-1.4020001.402001-3.50100002.5510千米(2月)303.501001-8.45400002.558.454001-14.40000010千米(5月)2121控制点Q_AUG_cms1.022000-2.5560002.556001-6.1390006.139001-16.97000016.970001-26.47000026.470001-40.440000Fig.6 The changes of water flow along Qiongjiang River in typical months of forecast year(2030)2.2.3琼江流域遂宁段主要干支流沿程水质模拟干流NH,-N浓度较高，其中城区以上干流河控制点Q_NOV_cms0.216300-0.4770000.477001-1.4930001.493001-5.14900002.5510干米(8月)图6 预测水平年（2 0 30 年）典型月份河道流量沿程变化5.149001-10.02000010.020001-15.630000段在12 月超过类水质标准，城区以下干流河段在13月超过类水质标准。各支流中，中心02.55(11 月)10千米150场河和半边河NH,-N浓度在1月超过类水质标准；双龙河和观音河在1月和6 月超过类水质标准；盐井河和莲花河NH3-N浓度在12 月和12月超过类水质标准；胡家沟河NH,-N浓度在1月、2 月、6 月和12 月超过类水质标准；石洞河和蟠龙河NH,-N浓度在13月、6 月、10 月、12四川环境月超过类水质标准；玉丰河NH,-N浓度仅在6月达到类水质标准。从时间角度来看，研究区域内4 5月、7 9 月水质相对较好，仅玉丰河NH3-N浓度超过类水质标准，其余月份均有多个河段不能达标。已有治污规划水平下各子流域每个季度代表月NH,-N沿程变化如图7 所示。42卷2121控制点NH3_FEB_mg/L0.166250-0.2500000.250001-0.5000000.500001-0.7500000.750001-1.000001.000001-2.745783控制点NH3_MAY_mg/L0.166250-0.250000.250001-0.5000000.50001-0.75000002.5510千米(2月)0.750001-1.0000001.000001-1.0699250.2.55(5月)10千米21控制点NH3_AUG_mg/L0.166250-0.2500000.250001-0.5000000.500001-0.7500000.750001-1.00001.000001-1.021517Fig.7 The changes of NH,-N along Qiongjiang River at typical months of forecast year(2030)从空间角度来看，干流TP浓度较高，其中城区以上干流河段在12 月超过类水质标准，城区范围干流河段在1月超过类水质标准，城区以下干流河段在1 3月、1112 月超过类水质标准。各支流中，中心场河和半边河TP浓度在1月超过类水质标准；双龙河和观音河在12 月和6月超过类水质标准；盐井河TP浓度在12 月和12 月超过类水质标准；莲花河和胡家沟河TP浓度在12 月、6 月和12 月超过类水质标准；控制点NH3_NOV_mg/L0.166250-0.2500000.250001-0.5000000.50001-0.75000002.5510千米(8月)图7 预测水平年（2 0 30 年）已有治污规划下典型月份流域NH,-N水质沿程变化290.750001-1.000002.551.000001-1.379849(11月)石洞河TP浓度在13月、6 月和12 月超过类水质标准；玉丰河TP浓度在1 4月、10 12 月超过类水质标准；蟠龙河TP浓度在13月、56月、9 12 月超过类水质标准。从时间角度来看，研究区域内干支流TP在7 8 月能够达到类水质标准，其余月份均有部分河段不能达标。已有治污规划水平下每个季度代表月各子流域TP沿程变化如图8 所示。10千米5期胡鑫等：SWAT模型在我国山区小流域水质预测中的应用一一以琼江流域为例151控制点TP_FEB_mg/L0.030901-0.0500000.050001-0.100000.100001-0.1500000.1500010.2000000.200001-0.463455控制点TP_MAY_mg/L300.030901-0.0500000.050001-0.100000.10001-0.15000002.5.510千米(2月)300.150001-0.2000000.2.550.200001-0.40000010干米(5月)21控制点TP_AUG_mg/L0.030901-0.0500000.050001-0.10000.100001-0.1500000.150001-0.2000000.200001-0.250000Fig.8 The changes of TP along Qiongjiang River at typical months of forecast year(2030)2.2.4重点断面全年水质模拟预测选择城区以上单元控制断面高滩断面、城区范围控制断面萝卜园断面和城区以下控制断面大安断面作为水质年内变化过程的重点分析断面。各典型断面NH,-N在已有治污规划水平下的年内变化过程如图9 所示。已有治污规划水平下，各典型断面NH3-N均不能稳定达标。高滩断面NH,-N浓度变化范围为0.2 53 1.2 8 4mg/L，除1月超过类水质标准外，其余月份均能达到类水质要求；萝卜园断面NH,-N浓度变化趋势与高滩断面一致，变化范围为0.2 9 9 1.0 12 mg/L，1月浓度超过类水质标准；大安断面较高滩、萝卜园断面水质稍差，NH,-N浓度在0.40 8 1.8 6 3mg/L之间波动，其中13月超过类水质标准，1月浓度最大，其余月份均能达到类水质标准。各典型断面TP在已有治污规划水平下的年内变化过程如图10 所示。已有治污规划水平下，各典型断面TP均不能稳定达标。高滩断面TP浓度控制点TP_NOV_mg/L0.030901-0.0500000.050001-0.100000.100001-0.150002.5510干米(8月)图：预测水平年（2 0 30 年）已有治污规划下典型月份子流域TP水质沿程变化类水质线1.00.80.60.40.2012 34 56789101112时间（月）图9预测水平年（2 0 30 年）已有治污规划下重点断面NH,-N年内变化Fig.9The changes of NH,-N at key sections under theexiting plan of pollution control in the forecast year(2030)变化范围为0.0 53 0.30 0 mg/L，除1月超过类水质标准外，其余月份均能达到类水质要求；萝卜园断面TP浓度变化趋势与高滩断面一致，变化范围为0.0 6 7 0.2 34mg/L，其中1 2 月超过类水质标准，1月浓度最大，其余月份均能达到类水质标准；大安断面较高滩、萝卜园断面水质稍0.150001-0.2000000.200001-0.4568312.01.8(T/u)N-HN1.61.41.22.55(11月)一高滩断面日一萝卜园断面。大安断面10干米152差，TP浓度在0.117 0.47 7 mg/L之间波动，其中13月超过类水质标准，1月浓度最大，其余月份均能达到类水质标准。0.6r0.5F(T/au)dL0.40.30.20.1F012234567891101112时间（月）图10 预测水平年（2 0 30 年）已有治污规划下重点断面TP年内变化Fig.10The changes of TP at key sections under theexiting plan of pollution control in the forecast year(2030)3结 论（1）基于琼江流域遂宁段构建的SWAT模型受地表径流、基流、地下水、坡度、蒸散发与土壤含水等参数的影响较为敏感；模型在率定期和验证期纳什效率系数分别为0.7 3和0.6 2，决定系数分别为0.7 2 和0.6 7，率定效果良好；模型对大安断面2 0 14年1月 2 0 17 年12 月的径流量预测值能实测值较好地与实测值吻合，除2 月外，2 0 18 年其余各月份TP与NH,-N值都能较好地与实测值吻合。（2）2 0 30 年琼江流域遂宁段干流从上游至下游流量呈现增加趋势，大安断面年均流量为16.2 0m/s，各支流中蟠龙河的流量最大，年均流量达到1.13m/s。从时间角度来看，大安断面逐月流量变化范围为2.8 2 40.44m/s，其中8 月流量最大。（3）2 0 30 年琼江干流NH,-N与TP浓度较高，在枯水季节尤其12 月水质较差，难以满足类水质标准，各支流在枯水季节水质也有不同程度超标。干流控制安居城区以上的高滩断面除1月的其余各月份水质能达到类水质要求，流经安居城区以后，干流萝卜园断面在2 月份TP也未能达到类水质要求，至出遂宁的大安断面，水质进一步变差，13月份难以满足类水质要求。（4）本文构建了基于琼江流域遂宁段构建的SWAT模型，解析了社会经济发展条件下琼江流域遂宁段的水质时空变化特征，为我国山区小流域水质的研究提供了思路，研究结果可为琼江流域水环境保护和政府决策，提供科学依据和技术支持。四川环境参考文献：【1刘凌雪，敖天其，胡正，等。琼江流域（安居段)水质及面源污染综合评价J.水土保持研究，2 0 19，2 6（6）：37 2-376.一一高滩断面一日一萝卜园断面一大安断面类水质线42卷2 杨玲玲，覃光华，张利平，等降雨移