一种基于改进D-S证据理论的水质判断方法与流程

本发明涉及环境监测领域，特别是指一种基于改进d-s证据理论的水质判断方法。

背景技术：

水是生命之源，是人们赖以生存的重要条件，随着社会工业化的步伐不断加快，国民经济迅速增长，社会发展与人口规模扩大造成的水资源污染成为了社会发展中不可忽视的问题。因此，加强水污染防治，建立水质工程对民生的改善，国家战略安全有很大的重要意义。而水质监测预警作为水污染防治中重要的一环，是环境监管中的重点和难点，只有正确分析监测点的水质状况，才能为上层行动决策做出最精准的估计，从而更好更快地控制污染的扩散。

多传感器数据融合就是通过对多个传感器的监测数据的支配和使用，利用多个传感器协同工作的优势，把传感器的监测数据在时间或空间上的互补数据和冗余数据按照融合算法融合出可以判定待测点环境状况的数据，从而得出对该监测点环境状况的一致性解释，提高整个水质监测预警的有效性与稳定性，从而可以更好地掌握水环境的情况，为水环境质量的评价估测提供有效的依据和手段。

目前，传统的多传感器数据融合算法大部分是由专家系统通过大量知识与经验计算出基本概率分配函数，再通过普通的融合算法对其融合。这种多传感器数据融合算法存在基本概率分配函数获取困难、计算量大、融合精确性不够、系统不适合大范围推广等问题。

技术实现要素：

本发明的主要目的在于针对传统的多传感器数据融合算法的基本概率分配函数获取困难、计算量大、融合精确性不够，提出一种基于改进d-s证据理论的水质判断方法。

本发明采用如下技术方案：

一种基于改进d-s证据理论的水质判断方法，其特征在于，包括如下步骤：

1)通过同类型水质传感器获取证据数据s＝{s1,s2,…,sm}，si为水质参数监测数据，1≤i＜m，m为水质参数监测数据的样本数量；

2)计算水质参数监测数据si基于五个不同水质类别的基本概率分配函数；

3)将m个水质参数监测数据所对应的基本概率分配函数进行融合；

4)将融合结果按照判决原则确定水质类别。

优选的，将五个水质类别设定为辨识框架θ＝{x1,x2,x3,x4,x5}，每个水质类别xj对应着一个限值sij，限值集s＝{si1,si2,si3,si4,si5}，步骤2)中，基本概率分配函数计算具体包括如下：

2.1)求出水质参数监测数据si与水质类别限值sij之间的海明距离:

hij＝|si-sij|；其中1≤j＜5

2.2)求出水质参数监测数据si对应的海明距离的平均值：

2.3)求出水质参数监测数据si与水质类别限值sij的相似度，利用正态曲线可得：

其中：a是大于0的系数，

2.4)求出水质参数监测数据si对应的水质类别xj的基本概率分配函数：

优选的，当水质参数监测数据满足以下条件

sij＜si＜si(j+1),1≤j＜5

则设水质参数监测数据样本所在区间的海明距离差值为：

h＝|hi(j+1)-hij|

新的海明距离为：

h'ik＝hik+h，k＝1,2,…,5,k≠j,k＝j+1

根据新的海明距离重复步骤2.2)-2.4)，得到基本概率分配函数。

优选的，步骤3)中，具体为：将同类型水质传感器获取的两个水质参数监测数据a和b的基本概率分配函数进行融合,a、b∈s，

融合后的基本概率分配函数为

其中f(x)是冲突证据概率分配函数，f(x)＝kq(x)，k为冲突因子,q(x)为证据对水质类别x的支持程度，将同类型水质传感器所获取证据数据的基本概率分配函数进行两两融合得到融合结果。

优选的，步骤4)中，所述判决的原则是：

规则一：该水质类别的基本概率分配函数是最大的；

规则二：该水质类别应该比其他的水质类别的基本概率分配函数要大于某一阙值；

规则三：不确定性的基本概率分配函数必须要小于某一阙值，；

规则四：该水质类别的基本概率分配函数要大于不确定性的基本概率分配函数。

由上述对本发明的描述可知，与现有技术相比，本发明具有如下有益效果：

本发明采用基于有序命题类问题的典型样本方法，简化了基本概率分配函数的获取，在大量减少基本概率分配函数获取的计算量情况下，也使得获取的数据更具有实时性和有效性。所有的基本概率分配函数进行融合，融合结果符合正常的思维逻辑，不像d-s证据理论存在“一票否决”的错误，也不像yager公式存在大量的不确定性，因此是合理有效的数据，可以按照判决原则来判定监测点的水质状况。

附图说明

图1是本发明方法在水质监测中的实现过程原理图。

以下结合附图和具体实施例对本发明作进一步详述。

具体实施方式

以下通过具体实施方式对本发明作进一步的描述。

一种基于改进d-s证据理论的水质判断方法，包括如下步骤：

1)获取证据数据

证据数据可以是溶解氧、高锰酸盐、氨氮等水质参数，以溶解氧水质参数为例，一个溶解氧传感器作为一个证据源，m为该证据源中证据的数量也即水质参数监测数据的样本数量，水质参数数据集s＝{s1,s2,…,sm}，对这水质参数数据si进行融合，1≤i＜m。

融合分为时间上的融合与空间上的融合。

设某一传感器t天内的水质参数监测数据stime＝{s1,s2,…,st},对这t个数据进行融合，可以获得t天内该监测点的水质状况。

设某一水域s个监测点传感器的水质参数监测数据sspace＝{s1,s2,…,ss}，对这s个数据进行融合，可以获得该水域在某时刻水质状况。

将时间融合的数据与空间融合的数据再次融合，便可得到监测水域的整体水质状况。

2)计算水质参数监测数据基于不同水质类别的基本概率分配函数

将五个水质类别设定为辨识框架θ＝{x1,x2,x3,x4,x5}，每个水质类别xj对应着一个限值，限值集s＝{si1,si2,si3,si4,si5}。则第i个水质参数监测数据在五个水质类别xj上的bpa的获取步骤如下：

2.1)求出水质参数监测数据si与水质类别限值sij之间的海明距离：

hij＝|si-sij|(1)

2.2)求出水质参数监测数据si对应的海明距离的平均值：

2.3)求出水质参数监测数据si与水质类别限值sij的相似度，利用正态曲线可得：

a是大于0的系数，为了保证正态曲线的斜度，取a＝5；b为：

2.4)求出水质参数监测数据si对应的水质类别xj的基本概率分配函数为：

另外，当水质参数监测数据满足以下条件时，

sij＜si＜si(j+1),1≤j＜5(6)

则设水质参数监测数据样本所在区间的海明距离差值为

h＝|hi(j+1)-hij|(7)

那么新的海明距离为hi'k＝hik+h，k＝1,2,…,5,k≠j,k＝j+1(8)

根据新的海明距离重复步骤2.2)-2.4)即重复公式(2)(3)(4)(5)的计算，即可获得基本概率分配函数。

采用典型方法来计算基本概率分配函数的目的是弥补传统基本概率分配函数需要依赖专家系统才能获取的缺点，使得获得的基本概率分配函数具有客观性与实时性。

3)将水质参数监测数据所对应的基本概率分配函数进行融合

将同类型水质传感器获取的两个水质参数监测数据a和b的基本概率分配函数进行融合，a、b∈s，

融合后的基本概率分配函数为

其中f(x)是冲突证据概率分配函数，f(x)＝kq(x)，k为冲突因子，q(x)为证据对水质类别x的支持程度，

将同类型水质传感器所获取证据数据的基本概率分配函数进行两两融合得到融合结果，五种不同水质类别下的融合结果可以按照判决原则来判定监测点的水质状况。

4)将融合结果按照判决原则确定水质类别

对水质类别进行判定，判决的原则是：

规则一：该水质类别的基本概率分配函数是最大的；

规则二：该水质类别应该比其他的水质类别的基本概率分配函数要大于某一阙值；

规则三：不确定性的基本概率分配函数m(θ)必须要小于某一阙值；

规则四：该水质类别的基本概率分配函数要大于不确定性的基本概率分配函数m(θ)；

将以上的基本规则用公式来表示，如下所示：

m(x1)＝max{m(x1),m(x2),m(x3),m(x4),m(x5)}(12)

m(x2)＝max{m(x2),m(x3),m(x4),m(x5)}(13)

m(x1)-m(x2)≥ε1(14)

m(θ)＜ε2(15)

m(x1)＞m(θ)(16)

那么x1为水质类别的判决结果，其中ε1,ε2为某一阙值。

不确定性的基本概率分配函数，是用1减去五个水质类别的确定性基本概率分配函数之和，得到的差值即为不确定性基本概率分配函数，之所以叫不确定性基本概率分配函数，也是相对前述步骤中所介绍的五个水质类别的确定性基本概率分配函数而言的，大多数正常情况下这五个水质类别的基本概率分配函数相加之和为1，此时不确定性基本概率分配函数就为0。如出现有不确定性的基本概率分配函数，则需要排除掉，不能用来判断水质的类别。

本发明采用对各类水质概率进行加权分配的融合方法，消除了传统d-s证据理论在水质参数数据融合中存在的逻辑错误问题，融合精度更高，说服力更强，从而可以对监测点的水质状况做出迅速、正确的判断。

应用举例

1)获取证据数据

通过溶解氧传感器，获取某一监测点的溶解氧水质参数监测数据s＝{4.02,4.35,4.88}，对这3个水质参数监测数据进行融合。

2)计算基本概率分配函数

水质类别限值s＝{7.5，6.0，5.0,3.0,2.0}

按照公式(1)～(8)计算基本概率分配函数

m1(x1)＝0.00000，m1(x2)＝0.00754，m1(x3)＝0.51915，m1(x4)＝0.46577，m1(x5)＝0.00754；

m2(x1)＝0.00000，m2(x2)＝0.00517，m2(x3)＝0.79520，m2(x4)＝0.19951，m2(x5)＝0.00012；

m3(x1)＝0.00037，m3(x2)＝0.00426，m3(x3)＝0.90278，m3(x4)＝0.09257，m3(x5)＝0.00000。

3)对基本概率分配函数进行融合

按照公式(9)～(11)对基本概率分配函数进行融合，

m(x1)＝0.00006，m(x2)＝0.00117，m(x3)＝0.92037，m(x4)＝0.07811，m(x5)＝0.00029。

4)针对融合结果按照判决原则判定水质类别

针对融合结果，按照公式(12)～(16)对监测点的水质类别进行判定。

《地表水环境质量标准》将地表水按功能高低依次划分为五类(i类、ii类、iii类、iv类、v类)，并给出了地表水环境质量标准五类水质基本项目标准限值表。由最终融合结果可以看出，m(x3)值最大，该监测点的水质属于iii类水质。故需要对该水质采取iii类水质相对应的污染处理措施。

上述仅为本发明的具体实施方式，但本发明的设计构思并不局限于此，凡利用此构思对本发明进行非实质性的改动，均应属于侵犯本发明保护范围的行为。