[发明专利]基于栅格化处理的平原河网区河道面源污染负荷确定方法

权利要求书

1.一种基于栅格化处理的平原河网区河道面源污染负荷确定方法，其特征在于，包括以下步骤：

(1)绘制河网多边形结构图

根据计算区域的水系图，对河网进行概化处理，绘制由河道包围所形成的网状多边形结构图；

(2)计算域栅格化及空间运算

借助工具，对计算区域的土地利用图层和栅格图层进行空间叠合分析，生成每个网格单元不同土地利用类型的范围，统计相应土地利用的面积为河网多边形第i个网格第j种土地利用的面积，km²；

(3)计算汇流权重因子

按公式(1)进行计算：

式中：P_i^k为第i个网格到第k条河道上的权重因子，％；d_i^k为第i个网格到第k条河道的最小距离，km；A_k为第k条河道的断面面积，m²；m为包围该网格的河道数量；

(4)确定面源污染流向

统计每个网格到第k条河道的最大权重因子，按公式(2)计算：

M_i＝max{P_i^k,k＝1…m} (2)

式中：M_i为某个河网多边形的第i个网格汇流入第k条河道的最大权重因子，％；

如果第i个网格到第k条河道的权重因子最大，那么该网格的面源污染负荷全部汇入第k条河道；

(5)面源污染汇流区划分

依次计算河网多边形内每个网格的面源污染汇入周边河道的最大权重因子，所有汇入同一条河道的网格即构成该河道的汇流区；

(6)统计汇流区土地利用类型面积

按公式(3)统计河网多边形不同土地利用类型的面积；

式中：A_j分别为河网多边形第j种土地利用类型的面积，km²；n为某一河道汇流区的栅格数量；依次对所有河网多边形重复上述计算过程，得到计算区域全部河道的汇流区范围及对应的各种土地利用类型面积；

(7)将河道汇流区的各类土地利用面积输入面源污染模型，计算与每条河道对应的面源污染量WL_k；

(8)根据面源污染量WL_k，结合工业废水和生活污水的点源调查资料，统计每条河道各种污染物的污染负荷总量WL_s，将该数据与政府公布的水污染物限制排污总量WL_L进行比较：

如果WL_s≤WL_L，说明河道的实际污染负荷小于或等于限制排污总量，无需对该河道制定相应减排措施；

如果WL_s>WL_L，说明该河道的实际污染负荷已经超过限制排污总量，需要制定相应的污染治理和减排措施。

2.根据权利要求1所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法，其特征在于：所述步骤(2)中的工具为GIS软件平台的“Intersect”工具。

3.根据权利要求1所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法，其特征在于：所述步骤(7)中面源污染类型包含城市径流污染、旱地径流污染、稻田径流污染。