一种基于软测量的河道排污口网格化溯源排查方法与流程

1.本发明涉及河道排污溯源技术领域，更具体的说是涉及一种基于软测量的河道排污口网格化溯源排查方法。


背景技术：

2.河道排污口排查整治是水环境治理的基础性工作。目前我国城市化地区虽然统计污水收集率平均达到90％以上，但是按照污染负荷收集率计算，平均只有60％，表明仍有大量污染物进入河道。入河排污口是污染物进入河道的最后一道“闸口”，只有全面摸清入河排污口的数量和污染物排放量，才能切实提高污染物收集能力，推动水环境质量从根本上好转。但河道排污口情况复杂，在水面以上排污口排查取得明显成效的现状下，当前面临的难点是水下排污口的排查。
3.基于此，提出基于河道网格化水位、水质监测的排污口溯源方法，来确定排污口调查的重点河段。其特点是水位和水质的监测不需要水下作业，水位的监测操作简便，且可在长时段内实施灵活性的动态监测。在水质监测模块，涉及到软测量模型的构建，以解决部分重要水质参数不能进行原位实时监测、不能获得高频监测数据的问题。软测量技术的核心就是建立数学模型。由于水质参数之间存在复杂的内在联系，天然水体中污染物的生化反应是一个强耦合的多输入、多输出的动态系统，具有高度非线性、不确定性等特点，传统的数学建模的方法受到了挑战。综上，基于排污口全面排查的需求，研发一种操作简单，能够快速进行排污口溯源的方法是本领域技术人员亟需解决的问题。


技术实现要素：

4.有鉴于此，本发明提供了一种基于软测量的河道排污口网格化溯源排查方法，克服了上述的缺陷。
5.为了实现上述目的，本发明提供如下技术方案：
6.一种基于软测量的河道排污口网格化溯源排查方法，具体步骤为：
7.将待测河流均匀划分为多个河段，并根据划分河段确定监测点；
8.获取各个监测点的水位监测数据和水质监测数据；
9.基于水位监测数据和水质监测数据，通过软测量方法获取各个监测点的污染物浓度数据；
10.通过各个监测点的污染物浓度数据，计算各个监测点的综合水质标识指数；
11.构建基于综合水质标识指数的污染评价体系，判断存在污染排放的河段；
12.对存在污染排放的河段采取监测点加密布设，直至达到预设的排污口定位精度要求。
13.其中，水质监测数据为常规可高频监测的水质数据，包括：固体悬浮物浓度ss、溶解氧do、氨氮nh3-n、浊度等数据。
14.可选的，软测量方法构建的步骤为：
15.确定输入输出变量；
16.根据输入输出变量确定子神经网络的个数；
17.根据输入输出变量以及子神经网络的个数，建立子神经网络模型；
18.对建立的子神经网络模型分别进行训练；
19.根据训练好的子神经网络模型对预测样本进行预测，将预测样本数据作为子神经网络模型的输入，获得监测点的污染物数据。
20.可选的，子神经网络模型的训练的具体步骤为：
21.基于训练数据对子神经网络分别进行训练，并在训练过程中利用快速efast对子神经网络的冗余隐神经元进行修剪；
22.直至各个子神经网络的误差均小于规定误差，停止训练，获得子神经网络模型。
23.可选的，判断存在污染排放的河段的具体步骤为：
24.构建基于综合水质标识指数的污染评价体系；
25.通过污染评价体系计算相邻监测点之间的综合水质标识指数变化情况；
26.若相邻监测点之间的综合水质标识指数发生变化，确定相邻监测点间的河段为存在污染排放的河段。
27.可选的，综合水质标识指数的获取公式为：
[0028][0029]
式中，ii为第i个监测点的综合水质指标；p
nh3’为第i个监测点的nh
3-n的单因子水质标识指数；p
tn’为第i个监测点的tn的单因子水质标识指数；p
tp’为第i个监测点的tp的单因子水质标识指数；p
cod’为第i个监测点的cod的单因子水质标识指数；p
do’为第i个监测点的do的单因子水质标识指数。
[0030]
可选的，单因子水质标识指数的获取公式为：
[0031]
p
wq
＝x1·
x2；
[0032]
式中，p
wq
为某一项水质指标的单因子水质指数；x1为待评价水质指标的水质类别；x2为监测数据在x1类水质变化区间中所处的位置，根据公式按四舍五入的原则计算确定。
[0033]
可选的，当水质为i～v类水时x1和x2的获取方式为：
[0034]
根据水质监测数据与国家标准的比较确定x1：x1＝1，表示该指标为i类水；x1＝2，表示该指标为ii类水；x1＝3，表示该指标为iii类水；x1＝4，表示该指标为iv类水；x1＝5，表示该指标为v类水。
[0035]
非溶解氧指标的x2获取公式为：
[0036][0037]
式中，cj为第j项水质指标的实测浓度，c
jk下
≤cj≤c
jk上
；c
jk下
为第j项水质指标第k类水区间浓度的下限值，k＝x1；c
jk上
为第j项水质指标第k类水区间浓度的上限值，k＝x1。
[0038]
溶解氧指标的x2获取公式为：
[0039][0040]
式中，c
do
为溶解氧的实测浓度；c
dok上
为第k类水中溶解氧浓度的区间边界上限值，k
＝x1；c
dok下
为第k类水中溶解氧浓度的区间边界下限值，k＝x1。
[0041]
可选的，当水质劣于v类时x1和x2的获取方式为：
[0042]
非溶解氧指标x1·
x2的获取公式为：
[0043][0044]
式中，cj为第j项水质指标实测浓度；c
j5上
为第j项指标v类水浓度上限值，x1·
x2根据计算结果按四舍五入的原则取小数点后一位确定。
[0045]
溶解氧指标x1·
x2的获取公式为：
[0046][0047]
式中，c
do
为溶解氧的实测浓度；c
do，5下
为溶解氧的v类水浓度下限值，c
do，5下
＝2.0mg/l；m为计算公式修正系数，取m＝4。
[0048]
可选的，相邻监测点之间的综合水质标识指数变化情况的计算公式为：
[0049]ii
》i
i-1
；
[0050]
式中，ii为第i个监测点的综合水质指标；i
i-1
为第i-1个监测点的综合水质指标。
[0051]
经由上述的技术方案可知，与现有技术相比，本发明公开了一种基于软测量的河道排污口网格化溯源排查方法，通过采用快速efast方法对神经网络的结构进行修剪，简化了神经网络的结构，加快了神经网络的反应速度的同时，也避免了拟合精度低的问题；本发明以综合水质标识指数的变化情况为依据，确定排污口的位置，降低了检测难度，提高了排污口辨识的速度与准确性。
附图说明
[0052]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
[0053]
图1为本发明的方法流程示意图；
[0054]
图2为软测量神经网络建模拓扑结构图；
[0055]
图3为单因子水质指数的结构示意图；
[0056]
图4(a)为水质类别为i～v类水时单因子水质标识指数为非溶解氧指标的x2的计算示意图；图4(b)为水质类别为i～v类水时单因子水质标识指数为溶解氧指标的x2的计算示意图；
[0057]
图5(a)为水质类别劣于v类水时单因子水质标识指数为非溶解氧指标的x2的计算示意图；图5(b)为水质类别劣于v类水时单因子水质标识指数为溶解氧指标的x2的计算示意图；
[0058]
图6为本发明实施例的河道排污口溯源的具体实例示意图。
具体实施方式
[0059]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0060]
本发明实施例公开了一种基于软测量的河道排污口网格化溯源排查方法，步骤如图1所述，具体为：
[0061]
步骤1、将8km的待测河流均匀划分为5个河段，根据河段确定6个监测点，开展网格化水位、水质监测；
[0062]
其中，确定监测点时以河段的划分位置为监测点位。
[0063]
水质监测数据包括固体悬浮物浓度ss、溶解氧do、氨氮nh
3-n、浊度等常规可高频监测的水质数据。
[0064]
步骤2、通过软测量获得各个监测点的污染物浓度数据，具体步骤为：
[0065]
步骤21、确定输入输出变量；
[0066]
污染物总氮tn、污染物总磷tp、化学需氧量cod等数据的监测较为复杂，尽管传感器技术的发展，这些重要的水质参数仍不能进行原位实时测量，需要建立基于常规高频水质监测参数(水位h、ph值、溶解氧do、浊度等)的软测量，来实现这些水质参数的高频“监测”。本实施例中将输入输出变量划分为三组，分别为：输出水质数据cod的输入辅助变量是水位h、固体悬浮物浓度ss、溶解氧do的监测数据；输出水质数据tn的输入辅助变量是水位h、do、nh
3-n、浊度的监测数据；输出水质数据tp的输入辅助变量是水位h、浊度的监测数据。
[0067]
步骤22、输入输出变量，分别建立三个子神经网络模型；
[0068]
根据输入输出变量分别建立子神经网络，本实施例以任一子神经网络为例：
[0069]
初始化神经网络：确定l-p-1的连接方式，对神经网络的权值进行随机赋值；即一个输入层有l个神经元，隐含层有p个神经元的单输出三层前向神经网络；
[0070]
隐含层的输出表达式为：
[0071][0072]
式中，zj(t)为隐含层第j个神经元的输出；为输入层权值；xi为神经网络的输入；i为输入层第i个神经元；
[0073]
隐含层神经元输出和神经网络输出的关系为：
[0074][0075]
式中，为输出层权值；
[0076]
定义误差函数为：
[0077][0078]
式中，t为训练样本的个数；yd(t)为第t为个训练样本的期望输出。
[0079]
步骤23、对建立的神经网络模型进行训练，直至子神经网络的误差小于规定误差；
[0080]
根据训练样本训练神经网络，并在训练过程中利用快速傅立叶幅值灵敏度检验扩展法(extend fourier amplitude sensitivity test，efast)对神经网络的冗余隐神经元进行修剪，以简化神经网络结构，增加神经网络的泛化能力和预测速度；
[0081]
在实施例中为计算各个隐含层神经元的灵敏度，对神经网络进行修剪，即删除低于5％灵敏度的隐含层神经元。
[0082]
当子神经网络的误差小于规定误差时，通过确定的各个网络的输入变量，输出变量，网络层数，网络隐含层个数，即得到三个子神经网络，结构如图2所示。
[0083]
步骤24、根据训练好的神经网络模型对预测样本进行预测，将预测样本数据(h、ss、do、nh
3-n、浊度的监测数据)作为网络的输入，利用构建的三个子神经网络输出获得监测点的cod、tn、tp水质浓度数据。
[0084]
步骤3、通过各个监测点的污染物浓度数据，计算各个监测点的综合水质标识指数ii变化情况，根据上下游监测点的综合水质标识指数ii变化情况判断划分的5个河段是否存在污染排入；
[0085]
具体为：根据监测点上下游的综合水质标识指数ii变化情况确定存在污染排放的河段。
[0086]
各监测点位的综合水质标识指数通过计算获得：
[0087][0088]
式中，ii为第i个监测点位的综合水质指标；p
nh3’为第i个监测点位的nh
3-n的单因子水质标识指数；p
tn’为第i个监测点位的tn的单因子水质标识指数；p
tp’为第i个监测点位的tp的单因子水质标识指数；p
cod’为第i个监测点位的cod的单因子水质标识指数；p
do’为第i个监测点位的do的单因子水质标识指数。
[0089]
可选的，单因子水质标识指数的获取公式为：
[0090]
p
wq
＝x1·
x2；
[0091]
式中，p
wq
为某一项水质指标的单因子水质指数；x1为待评价水质指标的水质类别；如图3所示，x2为监测数据在x1类水质变化区间中所处的位置，根据公式按四舍五入的原则计算确定。
[0092]
当水质为i～v类水时x1和x2的获取方式为：
[0093]
根据水质监测数据与国家标准的比较确定x1：x1＝1，表示该指标为i类水；x1＝2，表示该指标为ii类水；x1＝3，表示该指标为iii类水；x1＝4，表示该指标为iv类水；x1＝5，表示该指标为v类水。
[0094]
非溶解氧指标的x2获取公式为：
[0095][0096]
式中，cj为第j项水质指标的实测浓度，c
jk下
≤cj≤c
jk上
；c
jk下
为第j项水质指标第k类水区间浓度的下限值，k＝x1；c
jk上
为第j项水质指标第k类水区间浓度的上限值，k＝x1，获得的非溶解氧指标的x2如图4(a)所示。
[0097]
溶解氧指标的x2获取公式为：
[0098][0099]
式中，c
do
为溶解氧的实测浓度；c
dok上
为第k类水中溶解氧浓度高的区间边界值，k＝x1；c
dok下
为第k类水中溶解氧浓度低的区间边界值，k＝x1；获得的溶解氧指标的x2如图4(b)所示。
[0100]
当水质劣于v类时x1和x2的获取方式为：
[0101]
非溶解氧指标x1·
x2的获取公式为：
[0102][0103]
式中，cj为第j项水质指标实测浓度；c
j5上
为第j项指标v类水浓度上限值，x1·
x2根据计算结果按四舍五入的原则取小数点后一位确定；获得的非溶解氧指标x1·
x2如图5(a)所示。
[0104]
溶解氧指标x1·
x2的获取公式为：
[0105][0106]
式中，c
do
为溶解氧的实测浓度；c
do，5下
为溶解氧的v类水浓度下限值，c
do，5下
＝2.0mg/l；m为计算公式修正系数，取m＝4；获得的溶解氧指标x1·
x2如图5(b)所示。
[0107]
当上下游相邻监测点位的综合水质指标升高，即ii》i
i-1
；即可说明监测点位i-1与i间的河段存在污水排放；
[0108]
式中，ii为第i个监测点位的综合水质指标；i
i-1
为第i-1个监测点位的综合水质指标。
[0109]
步骤4、对存在污染物排放的河段，加密监测点位的布设，逐步缩小排查范围，直至达到预设的排污口定位精度要求。
[0110]
本实施例中所述的上下游均依据河流流向确定。
[0111]
如图6所示，本实施例中在网格化监测的5个河段中，河段3上下游两个点位间综合水质标识指数升高0.8，判断该河段存在污染排放，故单独对该1.6km河段进行监测点位的加密布设。采用对折法逐步缩小排污口排查范围。具体为：在河段3的中部布设水位、水质监测点位，计算综合水质标识指数，与上下游监测点位对比，判断污染排放位于河段前半段或是后半段。如此重复操作4次，即可将排污口定位精度控制为100m，达到预设的排污口定位精度要求，实现排污口位置的溯源。
[0112]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
[0113]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。