3. 污染负荷估算方法体系

流域土地利用变化对污染负荷的影响，最终要落到“量”的判断上。估算方法选得对不对，直接决定负荷结果能不能支撑管理决策。本节围绕三类主流方法——输出系数模型、排污系数法、水文水质机理模型（SWAT、HSPF等）——讲清各自的适用条件与选择依据，并用一组示例数据完成一次负荷估算的初步实战。

三类方法的核心逻辑与适用边界

输出系数模型把每一类土地利用视为一个“负荷发生单元”，用单位面积年输出系数乘以面积再求和，形式为 L = Σ(Eᵢ·Aᵢ) + 点源。它假设负荷与土地利用面积近似成正比，对降雨、地形、输移过程的刻画很弱。优点是数据需求低（只需土地利用面积和文献或本地实测系数），适合资料匮乏流域的年度尺度估算和情景对比；缺点是难以反映年内变化和迁移路径，系数外推风险大。

排污系数法从污染源端出发，按人口、畜禽存栏、农田施肥量等活动强度乘以相应排污系数，逐源累加。它擅长核算点源与分散畜禽、农村生活等源项的贡献结构，适合污染源普查、总量分配和减排方案设计；但对源到受纳水体的输移衰减考虑不足，通常需要配合入河系数修正。

水文水质机理模型（SWAT、HSPF）以日或小时步长模拟降水—产流—汇流—污染物迁移转化全过程，能输出逐日负荷、区分地表径流与地下路径贡献、评估管理措施（如退耕、缓冲带）的响应。代价是数据门槛高：需要DEM、土壤、气象、土地利用、农业管理与水质监测数据，且需率定验证。它适合有监测基础、需要过程机制和措施情景分析的流域。

方法选择的三条判断线

实际工作中可按三条线快速定位：数据线——只有土地利用图和文献系数，选输出系数模型；有污染源普查或统计数据，选排污系数法；有DEM、土壤、气象与逐日水质，才具备机理模型条件。时间线——只需年度总量，前两类足够；需逐日或汛期非点源响应，必须上机理模型。问题线——问“哪类土地利用贡献大”，输出系数模型配合情景面积调整即可；问“减排措施能削减多少”，机理模型的过程机制不可替代。三者并非互斥，常见做法是用排污系数法核算点源、输出系数模型估算非点源、再用机理模型校核总量。

初步实战：某子流域总氮负荷估算

设某子流域面积 120 km²，土地利用为耕地 45 km²、林地 50 km²、建设用地 15 km²、草地 10 km²。采用输出系数模型估算非点源总氮，系数取自同类流域文献值（示例假设）：耕地 12 kg/(hm²·a)、林地 2.5、建设用地 9、草地 4。注意单位换算：1 km² = 100 hm²。

非点源 TN 负荷（示例计算）
耕地：45 × 100 × 12   = 54 000 kg/a
林地：50 × 100 × 2.5  = 12 500 kg/a
建设用地：15 × 100 × 9 = 13 500 kg/a
草地：10 × 100 × 4    =  4 000 kg/a
非点源合计 = 84 000 kg/a = 84.0 t/a

点源（排污系数法，示例）：
生活污水：1.2 万人 × 1.8 g/(人·d) × 365 ≈ 7.9 t/a
畜禽养殖：折算后约 5.6 t/a
点源合计 ≈ 13.5 t/a

入河系数取 0.65（示例假设）：
入河负荷 ≈ (84.0 + 13.5) × 0.65 ≈ 63.4 t/a

结果解读要抓住两点：一是贡献结构，本例耕地贡献约 64%，是首要控制对象，这与“耕地占比高、系数大”的设定一致；二是敏感性，输出系数的不确定性通常在 ±30% 以上，应做系数上下限的区间估算，而不是只报一个点值。若该流域有出口断面逐月监测，可用实测通量反算入河系数，替代假设值，显著提升可信度。

综上，方法体系的选择本质是“数据条件—时间尺度—管理问题”三者的匹配。输出系数模型求快、排污系数法求源、机理模型求过程，实战中常组合使用：以系数法搭骨架、以机理模型校核关键断面，再用情景面积变化推演土地利用调整的负荷响应。