随着社会的高速发展，我国水环境质量面临的压力持续升高，水环境污染负荷也严重超标，在此背景下，我国近年来先后制定了《水污染防治行动计划》、《“十四五”重点流域水环境综合治理规划》等决策部署，其中多条细则指出目前国内的水环境的管理控制应该逐渐由浓度、目标总量控制转变为容量总量控制。黄石市的湖泊众多，水资源丰富，其中保安湖为《湖北省湖泊保护名录》中划定的省控重点湖泊，然而近年来保安湖水质状况并不稳定，不能满足生态文明建设需求，因此对保安湖开展基于总量控制的水环境容量研究至关重要。

  本研究以保安湖为研究对象，通过对保安湖流域水文气象、自然资源、地形地貌和社会经济等信息进行调查分析，建立保安湖流域基础数据库，并以此基础开展基于输出系数法的保安湖污染源解析工作。随后根据保安湖历年常规监测断面水质监测数据，基于WASP模型构建保安湖水质模型。使用模型试错法推算保安湖纳污染能力、核算日最大污染负荷总量，在等比例法的基础上结合当地实际情况制定出TMDL管理计划，最终为保安湖提出针对性污染控制对策。

  主要研究结果如下：

（1）保安湖的总氮、总磷污染现状均未能达到Ⅲ类水体的水质管理目标，近年水质呈现恶化趋势，仅达到Ⅳ类标准。总氮、总磷污染均呈现湖周高，湖心低的分布情况，且丰水期水质明显劣于枯水期，枯水期整体离散程度更稳定。

（2）基于输出系数法得出的保安湖污染源解析结果显示，总氮入湖量为 927.15 t/a、总磷入湖量为 138.21 t/a。其中总氮入湖污染负荷主要来源于城镇生活，总磷入湖污染负荷主要来源于水产养殖。

（3）基于WASP模型建立了保安湖水质模型，与实测值进行对比，总氮模拟值平均相对误差为15.03%，总磷模拟值平均相对误差为15.21%，可以认为基于WASP模型建立的保安湖水质模型具有良好的水质模拟效果。

（4）保安湖总氮TMDL为1797kg/d,总磷TMDL为107kg/d。将TMDL进行二级分配，总氮点源、面源分配量分别为363.08 kg/d和1253.88 kg/d；总磷分别为10.31 kg/d和95.97 kg/d。二级分配中保安镇为重点分配区域。结合保安湖流域现状，从技术和管理角度提出相应削减措施以保障TMDL计划的实施。

关键词：WASP模型; 保安湖；TMDL；水环境容量；污染负荷分配

[1] BORAH D K, YAGOW G, SALEH A, et al. Sediment and nutrient modeling for TMDL development and implementation [J]. Transactions of the Asabe, 2006, 49(4): 967-86.

[2] BORAH D K, AHMADISHARAF E, PADMANABHAN G, et al. Watershed Models for Development and Implementation of Total Maximum Daily Loads [J]. Journal of Hydrologic Engineering, 2019, 24(1).

[3] 王健, 向峰, 邱飞，等. 水质预测模型研究进展 [J]. 环境科学导刊, 2018, 37(04): 63-7.

[4] 谢伟. 美国TMDL制度发展及启示 [J]. 社会科学家, 2017, No.247(11): 100-6.

[5] 吴玲玲. 基于水质目标管理的颍河TMDL计划研究 [D]; 南京大学, 2018.

[6] 郑军, 李乐, 郑静. 欧美水生态环境管理历程和现状研究 [J]. 中国生态文明, 2022, No.52(04): 70-5.

[7] 沙健, 李青, 石春力，等. 基于美国TMDL计划下的湖库流域模型化管理概述; proceedings of the 中国环境科学学会2009年学术年会, 中国湖北武汉, F, 2009 [C].

[8] 侯智. 基于WASP模型的大冶湖TMDL管理研究 [D]; 武汉理工大学, 2018.

[9] 孙文章, 曹升乐, 徐光杰. 应用WASP对东昌湖水质进行模拟研究 [J]. 山东大学学报(工学版), 2008, No.168(02): 83-5+100.

[10] 崔慧萍, 彭森, 吴卿，等. 绿色校园雨水利用及WASP模型模拟评价 [J]. 中国给水排水, 2018, 34(21): 127-32.

[11] 王亚宁, 李一平, 程月，等. 大型浅水湖泊水质模型边界负荷敏感性分析 [J]. 环境科学, 2021, 42(06): 2778-86.

[12] 朱瑶, 梁志伟, 李伟，等. 流域水环境污染模型及其应用研究综述 [J]. 应用生态学报, 2013, 24(10): 3012-8.

[13] 陶子夜. 负荷历时曲线法在通顺河武汉段流域TMDL计划中的应用 [D]; 武汉理工大学, 2017.

[14] 赖格英, 于革. 太湖流域营养物质输移的模拟评估研究 [J]. 河海大学学报(自然科学版), 2007, No.163(02): 140-4.

[15] 范翻平. 基于Delft3D模型的鄱阳湖水动力模拟研究 [D]; 江西师范大学, 2010.

[16] 刘雪梅. 基于水动力-水质-水生态综合模型的查干湖多水源调控 [D]; 中国科学院大学(中国科学院东北地理与农业生态研究所), 2021.

[17] RUEDA-BAYONA J G, HORRILLO-CARABALLO J, CHAPARRO T R. Modelling of surface river plume using set-up and input data files of Delft-3D model [J]. Data in brief, 2020, 31: 105899.

[18] PENG S, FU G Y Z, ZHAO X H. Integration of USEPA WASP model in a GIS platform [J]. Journal of Zhejiang University-Science A, 2010, 11(12): 1015-24.

[19] 李军, 张明华. WASP富营养化水质模型在GIS中的完全集成研究 [J]. 长江流域资源与环境, 2012, 21(S1): 93-8.

[20] 程晨健, 李天文, 陈靖，等. 利用WASP模型和GIS可视化集成的水质监测与模拟——以渭河为例 [J]. 地下水, 2011, 33(02): 52-5.

[21] ZUO X Q, WU X Q, YU J. Extraction of Watershed Characteristics of Red River Based on SRTM DEM; proceedings of the 2nd International Conference on Computer and Automation Engineering (ICCAE), Singapore, SINGAPORE, F Feb 26-28, 2010 [C]. 2010.

[22] XU R, HUANG X X, LUO L, et al. A new grid-associated algorithm in the distributed hydrological model simulations [J]. Science China-Technological Sciences, 2010, 53(1): 235-41.

[23] 胡晞. 基于WASP模型的湘江湘潭段水质目标管理研究 [D]; 湘潭大学, 2013.

[24] 陈美丹, 姚琪, 徐爱兰. WASP水质模型及其研究进展 [J]. 水利科技与经济, 2006, (07): 420-2+6.

[25] BOYACIOGLU H, ALPASLAN M N. Total maximum daily load (TMDL) based sustainable basin growth and management strategy [J]. Environmental Monitoring and Assessment, 2008, 146(1-3): 411-21.

[26] 蔡明, 李怀恩, 庄咏涛，等. 改进的输出系数法在流域非点源污染负荷估算中的应用 [J]. 水利学报, 2004, (07): 40-5.

[27] 李俊晓, 李朝奎, 殷智慧. 基于ArcGIS的克里金插值方法及其应用 [J]. 测绘通报, 2013, No.438(09): 87-90+7.

[28] 刘佳明. 湖北梁子湖近45年来蓄洪容积变化遥感监测分析 [J]. 中国防汛抗旱, 2019, 29(04): 29-35.

[29] 王小菊, 王仁忠, 陈表娟，等. 浅析水环境监测的质量保证和质量控制 [J]. 能源环境保护, 2019, 33(04): 49-51+5.

[30] 孙崇亮, 王卷乐. 基于DEM的水系自动提取与分级研究进展 [J]. 地理科学进展, 2008, (01): 118-24.

[31] 周茹. 遥感影像在地理教学中的应用研究 [D]; 华中师范大学, 2021.

[32] 邱斌, 李萍萍, 钟晨宇，等. 海河流域农村非点源污染现状及空间特征分析 [J]. 中国环境科学, 2012, 32(03): 564-70.

[33] 廖安平, 彭舒, 武昊，等. 30m全球地表覆盖遥感制图生产体系与实践 [J]. 测绘通报, 2015, No.463(10): 4-8.

[34] 李卉, 苏保林. 平原河网地区农业非点源污染负荷估算方法综述 [J]. 北京师范大学学报(自然科学版), 2009, 45(Z1): 662-6.

[35] 朱梅. 海河流域农业非点源污染负荷估算与评价研究 [D]; 中国农业科学院, 2011.

[36] 夏菁, 张翔, 朱志龙，等. TMDL计划在长湖水污染总量控制中的应用 [J]. 环境科学与技术, 2015, 38(07): 176-81.

[37] 李静. 太湖流域农业污染压力评估 [D]; 山东师范大学, 2012.

[38] 孙辰. 汉江襄阳段水环境容量及总量控制研究 [D]; 华中科技大学, 2013.

[39] 连国奇. 百花湖富营养化调查及沉积物磷形态研究 [D]; 贵州师范大学, 2009.

[40] 彭森. 基于WASP模型的不确定性水质模型研究 [D]; 天津大学, 2010.

[41] 郭卫鹏. 基于WASP模型的黄家湖TMDL管理研究 [D]; 华中科技大学, 2016.

[42] 金梦. 基于WASP水质模型的水环境容量计算 [D]; 辽宁大学, 2011.

[43] 张永祥, 王磊, 姚伟涛，等. WASP模型参数率定与敏感性分析 [J]. 水资源与水工程学报, 2009, 20(05): 28-30.

[44] 崔广柏, 陈星, 向龙，等. 平原河网区水系连通改善水环境效果评估 [J]. 水利学报, 2017, 48(12): 1429-37.

[45] 嵇灵烨. 基于环境容量的总量控制方法比较研究 [D]; 浙江大学, 2018.

[46] 杨龙, 王晓燕, 孟庆义. 美国TMDL计划的研究现状及其发展趋势 [J]. 环境科学与技术, 2008, (09): 72-6.

[47] DILKS D W, FREEDMAN P L. Improved consideration of the margin of safety in total maximum daily load development [J]. Journal of Environmental Engineering-Asce, 2004, 130(6): 690-4.

[48] 王双玲. 基于流域水生态承载力的污染物总量控制技术研究 [D]; 武汉大学, 2014.

[49] 杜云彬. 洪泽湖生态系统健康评价与污染负荷总量控制研究 [D]; 重庆交通大学, 2021.

[50] LIU Q Q, LIU D, ZHANG Y, et al. Study on spatial distribution of pollutants and total amount reduction in the Duliujian River (Tianjin, China) [J]. Water Supply, 2020, 20(4): 1380-95.

[51] 杨水化, 彭正洪, 焦洪赞，等. 城市富营养化湖泊的外源污染负荷与贡献解析——以武汉市后官湖为例 [J]. 湖泊科学, 2020, 32(04): 941-51.

[52] TIAN Y, HU H Y, ZHANG J. Solution to water resource scarcity: water reclamation and reuse [J]. Environmental Science and Pollution Research, 2017, 24(6): 5095-7.

[53] 杨帆, 张明俊, 胡晓晖. 城镇化污水治理中对人工湿地的运用 [J]. 资源节约与环保, 2023, No.254(01): 63-6.

[54] 杨娜, 王趁义, 徐园园，等. 黑臭小微水体治理技术的研究现状与发展趋势 [J]. 工业水处理, 2021, 41(05): 15-21.

[55] 沈峥, 刘洪波, 张亚雷. 中国“厕所革命”的现状、问题及其对策思考 [J]. 中国环境管理, 2018, 10(02): 45-8.

[56] 罗其友, 刘洋, 唐华俊，等. 新时期我国农业结构调整战略研究 [J]. 中国工程科学, 2018, 20(05): 31-8.

[57] 朱蒋洁, 曾艳, 陈敬安，等. 我国农业面源污染治理技术研究进展 [J]. 四川环境, 2014, 33(03): 153-61.

[58] 高庆标, 黄力华. 临沂市农村环境污染现状及综合整治 [J]. 农业环境与发展, 2009, 26(03): 46-51.

[59] 汪劲. 中国生态补偿制度建设历程及展望 [J]. 环境保护, 2014, 42(05): 18-22.

[60] Brian C.P. Mercury minimization measures to meet total maximum daily load requirements [J]. Federal Facilities Environment Journal, 2003.