基于空间马尔可夫模型的流域非点源优先管理区识别
Identification of Non-Point Source Priority Management Areas in Watershed Based on Spatial Markov Model
DOI: 10.12677/aep.2025.153040, PDF,    科研立项经费支持
作者: 秦俏华, 王小胜, 李慧芩：河北工程大学数理科学与工程学院，河北 邯郸；侯 越, 王 帅：河北省大清河河务中心，河北 保定；王明哲, 贾茂平：河北省保定市水文勘测研究中心，河北 保定；庞树江*：河北工程大学数理科学与工程学院，河北 邯郸；中国科学院东北地理与农业生态研究所，吉林 长春
关键词: 非点源污染；优先管理区；SWAT模型；空间马尔可夫算法；最佳管理措施；Non-Point Source Pollution； Priority Management Areas； SWAT Model； Spatial Markov Algorithm； Best Management Practices
摘要: 非点源污染是导致地表和地下水环境恶化的重要原因之一。优先管理区(Priority Management Areas, PMAs)识别对非点源污染精准高效管控至关重要。然而，在已有PMAs识别研究中很少考虑河流系统中的传输损失，且缺乏对识别结果准确性和合理性的验证。基于此，本文通过耦合SWAT模型与空间Markov算法，构建了一个考虑河流滞留效应和削减负荷的多级PMAs识别框架，应用于海河水系的王快水库流域，并验证PMAs识别结果在不同管理措施下的准确性。主要结果如下：(1) 在校准和验证期SWAT模型性能良好，纳什效率系数NSE > 0.60、R2 > 0.79；(2) 河道氮磷滞留系数与径流量呈非线性关系，存在显著的径流阈值。当径流量介于1.14~8.23 m3/s，总氮(TN)和总磷(TP)的滞留系数趋于稳定。(3) TN的PMAs均为#29和#30号子流域，负荷贡献为74.35% (枯水年)、58.45% (平水年)和76.46% (丰水年)，面积占流域总面积的14.22%。(4) TP的PMAs为#33 (枯水年)、#28和#33 (平水年)、#29和#30 (丰水年)号子流域，负荷贡献分别达到53.19%、78.83%和69.69%，面积占1.20%、2.95%和14.22%。(5) 仿真结果显示，PMAs削减效率约为非PMAs的4.25倍(TN)和5.37倍(TP)，证实了PMAs识别结果的合理性。该研究结果可为相关流域的非点源污染控制提供有效方法支撑。
Abstract: Non-point source (NPS) pollution is one of the important contributors to the deterioration of the surface and groundwater environment. The identification of priority management areas (PMAs) is crucial for effective management of NPS pollution. However, transport losses within river ecosystems are seldom considered when identifying PMAs, and there is a lack of validation for the accuracy of the results. Therefore, a framework for multi-level PMAs identification was conducted, considering the effect of river retention and retention load by integrating the SWAT model with the Markov algorithm. This framework is applied to the Wangkuai Reservoir Watershed in the Haihe River system, and the accuracy of the PMAs identification results under different best management measures (BMPs) is verified. The main results are illustrated as follows: (1) The SWAT model showed strong performance with high statistical indicators (NSE > 0.60, R2 > 0.79) for the calibration and validation periods. (2) A non-linear relationship between the retention coefficient and runoff with a significant runoff threshold was detected. For both total nitrogen (TN) and total phosphorus (TP), the retention coefficients held steady when runoff surpassed 1.14~8.23 m3/s. (3) The PMAs for TN are sub-watersheds #29 and #30, with load contributions of 74.35% (dry year), 58.45% (normal year), and 76.46% (wet year), and the area accounts for 14.22% of the total basin area. (4) The PMAs for TP are sub-watershed #33 (dry year), sub-watershed #28 and sub-watershed #33 (normal year), and sub-watershed #29 and #30 (wet year), with load contributions reaching 53.19%, 78.83%, and 69.69% respectively, and the areas account for 1.20%, 2.95%, and 14.22% respectively. (5) The reduction efficiency of PMAs is approximately 4.25 times (TN) and 5.37 times (TP) that of non-PMAs, confirming the accuracy of the PMAs identification results. The results will provide a reference for NPS pollution control in similar watersheds.
文章引用：秦俏华, 王小胜, 侯越, 李慧芩, 王帅, 王明哲, 贾茂平, 庞树江. 基于空间马尔可夫模型的流域非点源优先管理区识别[J]. 环境保护前沿, 2025, 15(3): 321-336. https://doi.org/10.12677/aep.2025.153040

参考文献

[1] Tong, X., Lai, X. and Liang, Q. (2023) An Improved Non-Point Source Pollution Model for Catchment-Scale Hydrological Processes and Phosphorus Loads. Journal of Hydrology, 621, Article ID: 129588. [Google Scholar] [CrossRef] 
[2] 李洪庆, 陈明慧, 程飞飞, 等. 基于“源-汇”理论的深圳河湾流域非点源污染风险评价[J]. 农业工程学报, 2024, 40(19): 217-224.
[3] Carpenter, S.R., Caraco, N.F., Correll, D.L., Howarth, R.W., Sharpley, A.N. and Smith, V.H. (1998) Nonpoint Pollution of Surface Waters with Phosphorus and Nitrogen. Ecological Applications, 8, 559-568. [Google Scholar] [CrossRef] 
[4] 任衍淦, 吴明艳, 张亚群. 基于SWAT模型祖厉河流域非点源污染优先管理区识别研究[J]. 中国农村水利水电, 2024(8): 33-40.
[5] Wu, T., Zhu, L., Shen, S., Zhu, A., Shi, M. and Qin, C. (2023) Identification of Watershed Priority Management Areas Based on Landscape Positions: An Implementation Using SWAT+. Journal of Hydrology, 619, Article ID: 129281. [Google Scholar] [CrossRef] 
[6] Chen, L., Zhong, Y., Wei, G., Cai, Y. and Shen, Z. (2014) Development of an Integrated Modeling Approach for Identifying Multilevel Non-Point-Source Priority Management Areas at the Watershed Scale. Water Resources Research, 50, 4095-4109. [Google Scholar] [CrossRef] 
[7] Guo, Y., Wang, X., Melching, C. and Nan, Z. (2022) Identification Method and Application of Critical Load Contribution Areas Based on River Retention Effect. Journal of Environmental Management, 305, Article ID: 114314. [Google Scholar] [CrossRef] [PubMed]
[8] 桑灿, 李兆华, 冯学高, 等. 白莲河流域总磷关键源区识别与最佳管理措施(BMPs)研究[J]. 长江流域资源与环境, 2024, 33(2): 374-386.
[9] Ahsan, A., Das, S.K., Khan, M.H.R.B., Ng, A.W.M., Al-Ansari, N., Ahmed, S., et al. (2023) Modeling the Impacts of Best Management Practices (BMPs) on Pollution Reduction in the Yarra River Catchment, Australia. Applied Water Science, 13, Article No. 98. [Google Scholar] [CrossRef] 
[10] 郑博福, 刘海燕, 吴汉卿, 等. 中国农田磷流失风险评价及其关键驱动因素[J]. 农业工程学报, 2024, 40(2): 332-343.
[11] Tan, M.L., Gassman, P.W., Yang, X. and Haywood, J. (2020) A Review of SWAT Applications, Performance and Future Needs for Simulation of Hydro-Climatic Extremes. Advances in Water Resources, 143, Article ID: 103662. [Google Scholar] [CrossRef] 
[12] Arnold, J.G., Srinivasan, R., Muttiah, R.S. and Williams, J.R. (1998) Large Area Hydrologic Modeling and Assessment Part I: Model Development. JAWRA Journal of the American Water Resources Association, 34, 73-89. [Google Scholar] [CrossRef] 
[13] 刘伟, 李逢港, 李沁芮, 等. 基于SWAT模型的釜溪河流域非点源污染负荷模拟研究[J]. 环境污染与防治, 2023, 45(8): 1176-1182.
[14] Moriasi, D.N., Arnold, J.G., Van Liew, M.W., Bingner, R.L., Harmel, R.D. and Veith, T.L. (2007) Model Evaluation Guidelines for Systematic Quantification of Accuracy in Watershed Simulations. Transactions of the ASABE, 50, 885-900. [Google Scholar] [CrossRef] 
[15] Moriasi, D.N., Gitau, M.W., Pai, N. and Daggupati, P. (2015) Hydrologic and Water Quality Models: Performance Measures and Evaluation Criteria. Transactions of the ASABE, 58, 1763-1785. [Google Scholar] [CrossRef] 
[16] Ma, Q., Pang, Y. and Mu, R. (2021) Water Environmental Capacity Calculation Based on Control of Contamination Zone for Water Environment Functional Zones in Jiangsu Section of Yangtze River, China. Water, 13, Article 587. [Google Scholar] [CrossRef] 
[17] Grimvall, A. and Stålnacke, P. (1996) Statistical Methods for Source Apportionment of Riverine Loads of Pollutants. Environmetrics, 7, 201-213. [Google Scholar] [CrossRef] 
[18] 曲玲, 于志强, 郝启勇, 等. 基于SWAT模型的尼山水库面源污染最佳管理措施研究[J]. 水资源与水工程学报, 2022, 33(4): 32-39, 49.
[19] 王敏, 张雨桐, 李奇宸, 等. 安徽省兆河流域非点源污染模拟及最佳管理措施[J]. 水土保持通报, 2021, 41(5): 127-136, 153.
[20] Liu, R., Xu, F., Zhang, P., Yu, W. and Men, C. (2016) Identifying Non-Point Source Critical Source Areas Based on Multi-Factors at a Basin Scale with Swat. Journal of Hydrology, 533, 379-388. [Google Scholar] [CrossRef] 
[21] Risal, A. and Parajuli, P.B. (2022) Evaluation of the Impact of Best Management Practices on Streamflow, Sediment and Nutrient Yield at Field and Watershed Scales. Water Resources Management, 36, 1093-1105. [Google Scholar] [CrossRef]