4. 流域空间单元划分与数据关联框架

流域土地利用变化与污染负荷的关联分析，第一步不是建模，而是把“空间”和“数据”对齐。流域是一个连续的地理系统，而土地利用、气象、土壤、排污等数据往往来自不同来源、不同分辨率、不同坐标体系。如果不先建立统一的空间单元体系，后续的负荷估算和关联分析就会因为“空间错位”而失去意义。本节围绕三个层次展开：子流域划分、水文响应单元（HRU）构建，以及土地利用与污染负荷的空间对应关系，最后落到“土地利用—负荷”关联分析数据库的实操建库。

一、子流域划分：把连续流域切成可计算单元

子流域划分的核心依据是地形与水系。常用方法是基于数字高程模型（DEM）进行流向分析、汇流累积量计算和河网提取，再在河网节点或汇水口处设置出水口，将流域切分为若干子流域。划分时需关注三个参数：集水面积阈值、最小子流域面积和出水口位置。集水面积阈值决定河网密度，阈值越小河网越密、子流域越多；最小子流域面积用于合并过碎单元，避免产生大量无数据支撑的小单元。

以一个约 1200 km² 的示例流域为例，DEM 分辨率为 30 m，若集水面积阈值取 5 km²，可提取出约 40～60 个子流域；若阈值提高到 20 km²，子流域数量会明显减少。这里的关键判断是：子流域数量应与可用监测数据、土地利用数据的分辨率相匹配。如果只有 3 个水质监测断面，却划分出 80 个子流域，多数子流域将缺乏验证数据，关联分析的可信度反而下降。因此，子流域划分不是越细越好，而是“计算精度”与“数据支撑”之间的平衡。

二、HRU 构建：在子流域内部叠加下垫面异质性

子流域内部往往同时存在耕地、林地、建设用地等不同下垫面，其产流和污染物输出特征差异很大。水文响应单元（HRU）就是在子流域内部，按“土地利用类型 + 土壤类型 + 坡度等级”的组合划分出的均质单元。同一 HRU 内，认为其产流机制和污染负荷输出系数近似一致。

构建 HRU 的实操步骤是：先对土地利用栅格重分类（如耕地、林地、草地、建设用地、水域），再对土壤类型按水文分组（A/B/C/D 四类）重分类，最后对坡度分级（如 0～5%、5～15%、>15%）。三者叠加后，统计每个子流域内各组合的面积占比，并设定面积阈值（如小于子流域面积 1% 的组合并入邻近主导类型），避免 HRU 数量爆炸。示例中，一个 30 km² 的子流域经叠加后可得到 8～15 个有效 HRU，其中“耕地—C 类土壤—5～15% 坡度”这一组合可能占 35%，是后续污染负荷估算的重点对象。

三、土地利用与污染负荷的空间对应关系

建立空间对应关系的本质，是让每一条负荷数据都能找到它“来自哪块地”。常用做法有两种：一是“HRU 汇总法”，即先计算每个 HRU 的污染负荷，再按子流域汇总；二是“监测断面反推法”，即利用断面监测浓度和流量估算子流域负荷，再按 HRU 面积占比分摊。前者适合无监测资料或资料稀少的流域，后者适合有连续监测的流域。实际项目中常将两者结合：用输出系数法给出 HRU 负荷初值，再用监测断面数据校正子流域总量。

以总氮（TN）为例，假设示例流域采用输出系数：耕地 12 kg/(ha·a)、林地 2.5 kg/(ha·a)、建设用地 9 kg/(ha·a)、草地 4 kg/(ha·a)。某子流域面积 30 km²，其中耕地 35%、林地 22%、建设用地 18%、草地 12%，其余为水域和未利用地。则子流域 TN 负荷约为：耕地 1050 ha×12=12600 kg，林地 660 ha×2.5=1650 kg，建设用地 540 ha×9=4860 kg，草地 360 ha×4=1440 kg，合计约 20.55 t/a。这一数值再与断面监测估算值对比，若偏差超过 30%，需检查输出系数取值或 HRU 面积统计是否准确。

四、“土地利用—负荷”关联分析数据库的实操建库

数据库建议采用“空间主键 + 属性表”的结构。核心表包括：子流域表（SUB_ID、面积、出水口坐标）、HRU 表（HRU_ID、SUB_ID、土地利用类型、土壤组、坡度级、面积、占比）、土地利用变化表（HRU_ID、年份、类型、面积、转移来源）、负荷表（HRU_ID、指标、输出系数、负荷量、监测校正值）。四张表通过 HRU_ID 和 SUB_ID 关联，形成可查询、可追溯的分析底库。

建库时需特别注意三点：第一，统一坐标系统和投影，所有面积统计必须在投影坐标系下进行，避免经纬度直接算面积带来的偏差；第二，保留“年份”字段，因为土地利用变化分析需要多期数据对比，同一 HRU 在不同年份的类型和面积可能不同；第三，记录数据来源和版本，便于后续复核。完成建库后，即可用 SQL 或 GIS 属性查询快速回答“某子流域建设用地增加 10% 对应 TN 负荷增加多少”这类关联问题，为后续情景分析打下基础。