一种水体污染类型的界定方法

摘要

本发明公开一种水体污染类型的界定方法，包括S1收集监测数据，形成数据库；S2根据数据库，确定枯水边界和丰水边界；S3根据枯水边界及区域降雨量，确定参考基流；S4将全年分组，建立分组通量数据库；S5在通量核定时间边界范围内，建立污染物通量核定模型；S6利用污染物通量核定模型，将污染物通量分为枯水时段污染物通量和丰水时段污染物通量；S7根据丰水时段污染物通量和丰水时段污染物通量，计算污染强度指数，根据污染强度指数及占比对水体污染类型进行界定。本发明通过对污染类型的定量划分，提高了污染类型界定的准确性，实现了区域水环境精准化治理，为水体污染控制管理提供有力的决策支持，对水环境保护具有重要意义。

权利要求

1.一种水体污染类型的界定方法，其特征在于，包括以下步骤：S1、收集监测数据，形成数据库，所述监测数据包括断面流量Q、污染物浓度C和区域降雨量R；

S2、根据步骤S1得到的断面流量Q和区域降雨量R，确定枯水边界Tk和丰水边界Tf；

S3、根据步骤S2确定的枯水边界Tk及区域降雨量R，确定参考基流Q0和步长n；

S4、根据步骤S3确定的步长n、将全年分组，根据步骤S1中的监测数据断面流量Q和污染物浓度C，建立分组通量数据库；

2

所述通量数据库中流量浓度分布为：Qn～N(μn,σn)；

2

所述通量数据库中污染物浓度分布为：Cn～N(μn,σn)；

S5、在通量核定时间边界T范围内，建立污染物通量核定模型；

所述污染物通量核定模型的公式如下：其中，P表示t时段内通过断面的污染物通量，单位为吨；

Q(t)为t时段中流量分布；C(t)为t时段中污染物浓度分布；

t为通量核定时间边界T内的时间，k为单位转换系数；

S6、利用步骤S5得到的污染物通量核定模型，根据流量Q与参考基流Q0大小关系，将污染物通量分为枯水时段污染物通量Pk和丰水时段污染物通量Pf；

S7、根据步骤S6得到的枯水时段污染物通量Pk和丰水时段污染物通量Pf，计算污染强度指数，根据污染强度指数 及概率对水体污染类型进行界定，所述界定采用的标准为：

所述污染强度指数 的计算公式为：其中，Pf为丰水时段污染物通量；

Pk为枯水时段污染物通量；

在(0,1)区间范围内，越趋近于1，面源污染越重；

当Pf‑Pk≤0时，

a为1年。

2.根据权利要求1所述的界定方法，其特征在于，步骤S1中所述监测数据可通过实地监测或自动监测站数据收集获取。

3.根据权利要求1所述的界定方法，其特征在于，步骤S1中所述污染物包括但不限于化学需氧量COD、氨氮NH3‑N、总氮TN、总磷TP。

4.根据权利要求1所述的界定方法，其特征在于，步骤S2中所述枯水边界Tk为流量Q和区域降水量R均出现最小值的时段；所述丰水边界Tf为流量Q和区域降水量R均出现最大值的时段。

5.根据权利要求1所述的界定方法，其特征在于，步骤S3中所述参考基流Q0为枯水边界2

Tk中，连续最大时段区域降雨量为零时的流量Q，参考基流分布记为：Q0～N(μ0,σ0)。

6.根据权利要求1所述的界定方法，其特征在于，步骤S3中所述步长n为在枯水边界Tk中，区域降雨量为零时的连续最大天数。

7.根据权利要求1所述的界定方法，其特征在于，步骤S4中全年分组的公式为：INT＝N/n      (3)；

其中，INT为组数；N为365天或366天，n为步长。

8.根据权利要求1所述的界定方法，其特征在于，步骤S5中通量核定时间边界T的计算方法为：将步骤S1得到的检测数据断面流量Q和降雨R，按照时间顺序依次排列，选取枯水边界Tk和丰水边界Tf所包含的时段，即为所述通量核定时间边界T。

9.根据权利要求1所述的界定方法，其特征在于，步骤S6中所述枯水时段污染物通量Pk为枯水边界Tk内，Q≤Q0时的污染物通量；所述丰水时段污染物通量Pf为丰水边界Tf内，Q＞Q0时的污染物通量。

10.根据权利要求1‑9任一项所述的界定方法在水体污染治理中的应用。

说明书

一种水体污染类型的界定方法

技术领域

[0001] 本发明属于水体污染防治技术领域，具体来说，涉及一种水体污染类型的界定方法。

背景技术

[0002] 在水污染防治领域，环境污染可根据其污染排放是否明确可分为“点源污染”与“面源污染”。点源污染是指有固定排放点的污染源，在数学模型中常用一点表示以简化计
算，如工业废水及城市生活污水，由排放口集中汇入江河湖泊。面源污染相对于点源来说则
没有固定污染排放点，污染物以广域的、分散的、微量的形式进入水体。

[0003] 目前，对于水体污染类型的界定，通常采用专家评估的经验法，或者是基于理论排放系数概算方法加以区分。然而，该类技术方法存在较大的不确定问题，导致点源污染和面
源污染类型界定不清，无法实现区域水环境精准化治理。

[0004] 因此，开发一种简便、高效的水体污染类型的界定方法，可以针对不同水体污染类型采用不同的水污染治理措施，为水体污染控制管理提供有力的决策支持，提高了污染防
治决策的成效性，对水环境保护具有重要意义。

发明内容

[0005] 为了解决以上技术问题，本发明基于大数据挖掘技术，提供了一种简便、高效的水体污染类型的界定方法。

[0006] 为了实现上述目的，本发明采用以下技术方案：

[0007] 一种水体污染类型的界定方法，包括以下步骤：

[0008] S1、收集监测数据，形成数据库，所述监测数据包括断面流量Q、污染物浓度C和区域降雨量R；

[0009] S2、根据步骤S1得到的断面流量Q和区域降雨量R，确定枯水边界Tk和丰水边界Tf；

[0010] S3、根据步骤S2确定的枯水边界Tk及区域降雨量R，确定参考基流Q0和步长n；

[0011] S4、根据步骤S3确定的步长n、将全年分组，根据步骤S1中的监测数据断面流量Q和污染物浓度C，建立分组通量数据库；

[0012] 所述通量数据库中流量浓度分布为：Qn～N(μn,σn2)；

[0013] 所述通量数据库中污染物浓度分布为：Cn～N(μn,σn2)；

[0014] S5、在通量核定时间边界T范围内，建立污染物通量核定模型；

[0015] 所述污染物通量核定模型的公式如下：

[0016]

[0017] 其中，P表示t时段内通过断面的污染物通量，单位为吨；

[0018] Q(t)为t时段中流量分布；C(t)为t时段中污染物浓度分布；

[0019] t为通量核定时间边界T内的时间，k为单位转换系数；

[0020] S6、利用步骤S5得到的污染物通量核定模型，根据流量Q与参考基流Q0大小关系，将污染物通量分为枯水时段污染物通量Pk和丰水时段污染物通量Pf；

[0021] S7、根据步骤S6得到的枯水时段污染物通量Pk和丰水时段污染物通量Pf，计算污染强度指数，根据污染强度指数 及概率两个指标对水体污染类型进行界定，所述界定采用
的标准为：

[0022]

[0023] 所述污染强度指数 的计算公式为：

[0024]

[0025] 其中，Pf为丰水时段污染物通量；

[0026] Pk为枯水时段污染物通量；

[0027] 在(0,1)区间范围内，越趋近于1，面源污染越重；

[0028] a为1年。

[0029] 优选地，步骤S1中所述监测数据可通过实地监测或自动监测站数据收集获取。

[0030] 优选地，步骤S1中所述污染物包括但不限于化学需氧量(COD)、氨氮(NH3‑N)、总氮(TN)、总磷(TP)。

[0031] 优选地，步骤S2中所述枯水边界Tk为流量Q和区域降水量R均出现最小值的时段；所述丰水边界Tf为流量Q和区域降水量R均出现最大值的时段。

[0032] 优选地，步骤S3中所述参考基流Q0为枯水边界Tk中，连续最大时段区域降雨量为零2
时的流量Q，参考基流Q0分布为：Q0～N(μ0,σ0)。

[0033] 优选地，步长n为在枯水边界Tk中，区域降雨量为零时的连续最大天数。

[0034] 优选地，步骤S4中全年分组的公式为：

[0035] INT＝N/n    (3)；

[0036] 其中，INT为组数；N为365天或366天，n为步长。

[0037] 优选地，步骤S5中通量核定时间边界T的计算方法为：将步骤S1得到的检测数据断面流量Q和降雨R，按照时间顺序依次排列，选取枯水边界Tk和丰水边界Tf所包含的时段，即
为通量核定时间边界T。

[0038] 优选地，步骤S6中所述枯水时段污染物通量Pk为枯水边界Tk内，Q≤Q0时的污染物通量；所述丰水时段污染物通量Pf为丰水边界Tf内，Q＞Q0时的污染物通量。

[0039] 本发明还提供了上述水体污染类型的界定方法在水体污染防治中的应用。

[0040] 所述应用为：根据界定的风险等级污染类型，采取不同的污染控制策略，具体如下：

[0041]

[0042] 其中，1为点源控制措施：针对断面上游区域，对具体的污染排放口进行污染治理。接入市政污水治理系统，增加末端污水处理设备，或者采取关停等方式；

[0043] 2为点面混合型控制措施：采用综合治理措施，既要对点源型污染进行治理，也要对面源污染进行防控；此外，混合型污染控制还要对点面源污染控制措施的选择、布局加以
优化，发挥工程效益的合力作用；

[0044] 3为面源控制措施；针对断面上游区域，加强面源污染的防控，采取生态沟渠拦截、湿地生态净化系统和/或水源涵养林修复措施。

[0045] 本发明的有益效果为：

[0046] (1)本发明的河流水体污染类型的界定方法，基于大数据挖掘技术，对水体污染类型进行了定量划分，精确度高，提高了水体污染类型界定的准确性。

[0047] (2)本发明根据污染强度指数和概率两个指数，将水体污染类型及风险等级进行了定量划分，可以针对不同水体污染类型采用不同的水污染治理措施，实现了区域水环境
精准化治理，为水体污染控制管理提供有力的决策支持，提高了污染防治决策的成效性，对
水环境保护具有重要意义。

具体实施方式

[0048] 以下结合特定的具体实例说明本发明的实施方式，在进一步描述本发明具体实施方式之前，应理解，本发明的保护范围不局限于下述特定的具体实施方案；还应当理解，本
发明实施例中使用的术语是为了描述特定的具体实施方案，而不是为了限制本发明的保护
范围。

[0049] 除非另外定义，本文中使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同意义。

[0050] 一种河流水体污染类型的界定方法，包括以下步骤：

[0051] S1、收集监测数据断面流量Q、污染物浓度C和区域降雨量R，形成数据库；所述监测数据可通过实地监测或自动监测站数据收集获取；所述污染物包括但不限于化学需氧量
(COD)、氨氮(NH3‑N)、总氮(TN)、总磷(TP)。

[0052] S2、枯水边界和丰水边界由断面流量和区域降雨量指标共同确定，根据步骤S1得到的断面流量Q和区域降雨量R，确定枯水边界Tk和丰水边界Tf；

[0053] 其中，枯水边界Tk为QMin∩RMin，流量和区域降雨量均出现最小值的时段；

[0054] 丰水边界Tf为QMax∩RMax，流量和区域降雨量均出现最大值的时段。

[0055] S3、基本径流是河道中能常年存在的那部分径流，主要是由地下水补给，基流大小与流域面积和自然地理条件有关，它是水资源评价中一个很重要的特征值；

[0056] 根据步骤S2确定的枯水边界Tk及区域降雨量R，确定参考基流Q0和步长n；

[0057] 在枯水边界Tk中，连续最大时段区域降雨量为零时的流量Q，确定为参考基流Q0,参2
考基流Q0分布记为：Q0～N(μ0,σ0)；其中，在枯水边界Tk中，区域降雨量为零时的连续最大天
数，记为步长n，单位为天。

[0058] S4、根据步骤S3确定的步长n、将全年分组，根据步骤S1中的监测数据断面流量Q和污染物浓度C，建立分组通量数据库；

[0059] 全年分组的公式为：

[0060] INT＝N/n    (3)；

[0061] 其中，INT为组数；N为365天或366天，n为步长，单位为天；

[0062] 所述通量数据库中流量浓度分布为：Qn～N(μn,σn2)；

[0063] 所述通量数据库中污染物浓度分布为：Cn～N(μn,σn2)；

[0064] S5、将步骤S1得到的检测数据断面流量Q和降雨R，按照时间顺序依次排列，选取枯水边界Tk和丰水边界Tf所包含的时段，记为通量核定时间边界T；

[0065] 在通量核定时间边界T范围内，建立污染物通量核定模型；

[0066] 所述污染物通量核定模型的公式如下：

[0067]

[0068] 其中，P表示t时段内通过断面的污染物通量，单位为吨；

[0069] Q(t)为t时段中流量分布；C(t)为t时段中污染物浓度分布；

[0070] t为通量核定时间边界T内的时间，k为单位转换系数；

[0071] S6、利用步骤S5得到的污染物通量核定模型，根据流量Q与参考基流Q0大小关系，将污染物通量分为枯水时段污染物通量Pk和丰水时段污染物通量Pf；

[0072] 枯水边界Tk内，Q≤Q0时的污染物通量，记为枯水时段污染物通量Pk；

[0073] 丰水边界Tf内，Q＞Q0时的污染物通量，记为丰水时段污染物通量Pf。

[0074] S7、根据步骤S6得到的枯水时段污染物通量Pk和丰水时段污染物通量Pf，计算污染强度指数，根据污染强度指数 及概率对水体污染类型进行界定，所述界定采用的标准为：

[0075]

[0076] 所述污染强度指数 的计算公式为：

[0077]

[0078] 其中，Pf为丰水时段污染物通量；

[0079] Pk为枯水时段污染物通量；

[0080] 在(0,1)区间范围内，越趋近于1，面源污染越重；

[0081] a为1年。

[0082] 本发明还提供了上述水体污染类型的界定方法在水体污染防治中的应用。

[0083] 所述应用为：根据界定的风险等级和污染类型，采取不同的污染控制策略，具体如下：

[0084]

[0085] 其中，1代表点源控制措施：针对断面上游区域，对具体的污染排放口进行污染治理。接入市政污水治理系统，增加末端污水处理设备，或者采取关停等方式；

[0086] 2代表点面混合型控制措施：采用综合治理措施，既要对点源型污染进行治理，也要对面源污染进行防控；此外，混合型污染控制还要对点面源污染控制措施的选择、布局加
以优化，发挥工程效益的合力作用；

[0087] 3代表面源控制措施；针对断面上游区域，加强面源污染的防控，采取生态沟渠拦截、湿地生态净化系统和/或水源涵养林修复措施。

[0088] 本发明根据污染强度指数和概率两个指数，将水体污染类型及风险等级进行了定量划分，可以针对不同水体污染类型采用不同的水污染治理措施，实现了区域水环境精准
化治理，为水体污染控制管理提供有力的决策支持，提高了污染防治决策的成效性，对水环
境保护具有重要意义。

[0089] 实施例1一种河流水体污染类型的界定方法

[0090] 一种河流水体污染类型的界定方法，包括以下步骤：

[0091] S1、监测断面流量(Q)、污染物浓度(C)以及区域降雨量(R)等数据形成数据库，其中，所述污染物包括但不限于化学需氧量(COD)、氨氮(NH3‑N)、总氮(TN)、总磷(TP)，如表1
所示。

[0092] 表1断面流量、污染物浓度及区域降雨量数据库

[0093]年 月 日 20‑20时降水量mm 流量数据m3/s 水质浓度
2020 1 1 0 1.00 2.00
2020 1 2 0 1.20 1.80
2020 1 3 0 0.80 1.70
2020 1 4 0 0.90 2.20
2020 1 5 0 0.70 2.10
2020 1 6 0 0.90 2.00
2020 1 7 5 0.83 1.60
2020 1 8 0 1.02 1.96
2020 1 9 0 0.91 1.97
2020 1 10 3 1.14 2.00
2020 1 11 2 0.82 2.23
2020 1 12 0 1.07 2.07
… … … … … …
… … … … … …
2020 12 26 0 0.87 1.77
2020 12 27 0 1.03 1.86
2020 12 28 0 1.30 1.81
2020 12 29 2 1.06 1.67
2020 12 30 0 0.96 1.86
2020 12 31 0 0.93 1.82

[0094] S2、确定枯水边界和基本径流：

[0095] 2020年1月1日至2020年1月6日，持续降雨量RMin＝0时最大时段，参考基流记作Q0～N(0.92,0.16)，n＝6。

[0096] S3、建立通量数据库：

[0097] 根据时间步长n＝6，将2020年1月1日至12月31日，共计分为61组。

[0098] 分组1：Q1～N(0.92,0.16)，C1～N(1.97,0.17)；

[0099] 分组2：Q2～N(0.97,0.12)，C2～N(1.97,0.19)；

[0100] ……

[0101] 分组61：Q61～N(1.02,0.14)，C61～N(1.80 0.08)。

[0102] S4、计算污染强度指数：

[0103] 根据分组通量数据库(分组流量与浓度分布特征)，分别抽样数据10万次，形成10万组计量数据组，形成污染物通量核定模型，根据流量Q与参考基流Q0大小关系，将污染物
通量分为枯水时段污染物通量Pk和丰水时段污染物通量Pf，并计算污染强度指数。根据计
算，其污染强度指数 为0.6，概率为75％，污染控制型为点面混合型，中高风险，应采取点
面混合型控制措施。

[0104] 以上是结合具体实施例对本发明进一步的描述，但这些实施例仅是范例性的，并不对本发明的范围构成任何限制。本领域技术人员应该理解的是，在不偏离本发明的精神
和范围下可以对本发明技术方案的细节和形式进行修改或替换，但这些修改和替换均落入
本发明的保护范围内。