一种大气污染物的污染源定位方法、装置、介质及设备与流程

1.本发明涉及污染源定位技术领域，尤其涉及一种大气污染物的污染源定位方法、装置、介质及设备。


背景技术：

2.大气污染物的监测是大气污染防治的基础，在环境保护工作中起着重要的作用，是不可缺少的环节之一。但传统的空气污染物监测方法无法快速准确地定位污染源的位置，例如，利用无人机搭载微电子传感器在工业区进行飞行监测，将监测数据通过无线电传输给主机，以定位污染源的位置，但是当一个区域内有多个污染源时，定位效率会大大下降。


技术实现要素：

3.本发明实施例的目的是提供一种大气污染物的污染源定位方法、装置、介质及设备，能够提高大气污染源的定位效率，同时提高定位的精确度。
4.为实现上述目的，本发明实施例提供了一种大气污染物的污染源定位方法，包括：
5.根据无人机上的气体传感器检测到的预设范围内的污染源的污染程度，构建污染源搜索模型；
6.利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，以根据所述最优搜索路径找到污染源的位置。
7.作为上述方案的改进，所述污染源搜索模型具体为：
[0008][0009]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a∈{lv′
,a
，lv′
,na
，l
v,a
，l
v,na
}为无人机搜索污染源的决策集合，lv′
,a
表示当污染源的污染程度为第二污染程度时，进行污染源的定位搜索，l
v,na
表示当污染源的污染程度为第一污染程度时，不进行污染源的定位搜索，l
v,a
表示当污染源的污染程度为第一污染程度时，进行污染源的定位搜索，lv′
,na
表示当污染源的污染程度为第二污染程度时，不进行污染源的定位搜索。
[0010]
作为上述方案的改进，当满足预设的第一参数条件时，a＝lv′
,a
∪l
v,na
；当满足预设的第二参数条件时，a＝l
v,a
∪lv′
,na
；
[0011]
其中，所述第一参数条件具体为：
[0012]
c《[p
0va
+(1-p0)u+f](pv′-pv)+(1-pv′
)f
[0013]
[0014][0015]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a＝lv′
,a
∪l
v,na
表示当污染源的污染程度为第二污染程度时，进行污染源的定位搜索，且当污染源的污染程度为第一污染程度时，不进行污染源的定位搜索；
[0016]
所述第二参数条件具体为：
[0017]
c》(1-pv′
)f+(pv′-pv)[p
0va
+(1-p0)u+f]
[0018][0019][0020]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a＝l
v,a
∪lv′
,na
表示当污染源的污染程度为第一污染程度时，进行污染源的定位搜索，且当污染源的污染程度为第二污染程度时，不进行污染源的定位搜索。
[0021]
作为上述方案的改进，所述气体传感器成功定位污染源的概率包括以下情况中的至少一种：
[0022]
当污染源的污染程度为第一污染程度时，气体传感器成功定位污染源的概率为p0(f)＝[1-exp(l-f)]pv，式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv为污染源的污染程度为第一污染程度的概率；
[0023]
当污染源的污染程度为第二污染程度时，气体传感器成功定位污染源的概率为p0(f)＝[1-exp(l-f)]pv′
，式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv′
为污染源的污染程度为第二污染程度的概率。
[0024]
作为上述方案的改进，所述气体传感器检测污染物的成本，具体为：
[0025]
c(h)＝kh
[0026]
式中，k表示气体传感器对单个污染源进行定位的成本，h表示污染源分布的分散程度。
[0027]
作为上述方案的改进，所述利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，包括：
[0028]
步骤a，获取大气污染物的浓度γ＝m
tiβ
并初始化参数；其中，t为污染物检测次数，i为大气污染物检出限，β为大气污染物，β∈[0,1]，m为污染物传感器检测形式，m
t
＝m
t-1
+(b/m
t-1
×
ngen)，m0取值为0.01，b取值为0.025，ngen表示最大迭代次数；
[0029]
步骤b，在预设范围内随机初始化污染源位置，将所述污染源搜索模型作为目标函
数，求解所述目标函数，计算当前污染源位置的最优值；
[0030]
步骤c，判断是否满足预设的终止条件，若是，输出最优值；若否，进入步骤d；
[0031]
步骤d，根据c1＝2exp(-(4t/ngen)2)计算收敛因子c1；其中，exp为指导函数，t为污染物检测次数，ngen表示最大迭代次数，d＝2-2t/ngen；
[0032]
步骤e，根据γ＝m
tiβ
计算预设范围内的每个污染源的大气污染物的浓度γi，根据a＝2dr-d计算收敛因子a；其中r∈[0,1]，d＝2-2t/ngen；
[0033]
步骤f，判断收敛因子c1是否大于1，若是，根据控制无人机进行全局搜索并获取最优置；若否，根据控制无人机进行局部搜索并获取最优置；其中，搜索并获取最优置；其中，搜索并获取最优置；其中，表示气体传感器在第i个污染源位置进行第t次污染物检测，r∈[0,1]，g
*
表示全局最优解，表示气体传感器检测第j维进行第t次污染物检测，表示气体传感器检测第k维进行第t次污染物检测，γi表示第i个污染源的污染物浓度，表示第j维气体传感器的位置，c3∈[0,1]，lbj、ubj为变量范围第j维的下界和上界；
[0034]
步骤g，执行t＝t+1，并返回步骤c。
[0035]
为实现上述目的，本发明实施例还提供了一种大气污染物的污染源定位装置，包括：
[0036]
污染源搜索模型构建模块，用于根据无人机上的气体传感器检测到的预设范围内的污染源的污染程度，构建污染源搜索模型；
[0037]
污染源定位模块，用于利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，以根据所述最优搜索路径找到污染源的位置。
[0038]
作为上述方案的改进，所述利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，包括：
[0039]
步骤a，获取大气污染物的浓度γ＝m
tiβ
并初始化参数；其中，t为污染物检测次数，i为大气污染物检出限，β为大气污染物，β∈[0,1]，m为污染物传感器检测形式，m
t
＝m
t-1
+(b/m
t-1
×
ngen)，m0取值为0.01，b取值为0.025，ngen表示最大迭代次数；
[0040]
步骤b，在预设范围内随机初始化污染源位置，将所述污染源搜索模型作为目标函数，求解所述目标函数，计算当前污染源位置的最优值；
[0041]
步骤c，判断是否满足预设的终止条件，若是，输出最优值；若否，进入步骤d；
[0042]
步骤d，根据c1＝2exp(-(4t/ngen)2)计算收敛因子c1；其中，exp为指导函数，t为污染物检测次数，ngen表示最大迭代次数，d＝2-2t/ngen；
[0043]
步骤e，根据γ＝m
tiβ
计算预设范围内的每个污染源的大气污染物的浓度γi，根据a＝2dr-d计算收敛因子a；其中r∈[0,1]，d＝2-2t/ngen；
[0044]
步骤f，判断收敛因子c1是否大于1，若是，根据控制
无人机进行全局搜索并获取最优置；若否，根据控制无人机进行局部搜索并获取最优置；其中，搜索并获取最优置；其中，搜索并获取最优置；其中，表示气体传感器在第i个污染源位置进行第t次污染物检测，r∈[0,1]，g
*
表示全局最优解，表示气体传感器检测第j维进行第t次污染物检测，表示气体传感器检测第k维进行第t次污染物检测，γi表示第i个污染源的污染物浓度，表示第j维气体传感器的位置，c3∈[0,1]，lbj、ubj为变量范围第j维的下界和上界；
[0045]
步骤g，执行t＝t+1，并返回步骤c。
[0046]
为实现上述目的，本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序；其中，所述计算机程序在运行时控制所述计算机可读存储介质所在的设备执行如上述所述的大气污染物的污染源定位方法。
[0047]
为实现上述目的，本发明实施例还提供了一种设备，包括处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序，所述处理器在执行所述计算机程序时实现如上述所述的大气污染物的污染源定位方法。
[0048]
与现有技术相比，本发明实施例提供一种大气污染物的污染源定位方法、装置、介质及设备，通过根据无人机上的气体传感器检测到的预设范围内的污染源的污染程度，构建污染源搜索模型；利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，以根据所述最优搜索路径找到污染源的位置，能够提高大气污染源的定位效率，同时提高定位的精确度。
附图说明
[0049]
图1是本发明实施例提供的一种大气污染物的污染源定位方法的流程图；
[0050]
图2是本发明实施例提供的一种大气污染物的污染源定位方法的又一流程图；
[0051]
图3是本发明实施例提供的一种大气污染物的污染源定位装置的结构框图；
[0052]
图4是本发明实施例提供的一种设备的结构框图；
具体实施方式
[0053]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0054]
参见图1-2，是本发明实施例提供的一种大气污染物的污染源定位方法的流程图，所述大气污染物的污染源定位方法，包括：
[0055]
s1、根据无人机上的气体传感器检测到的预设范围内的污染源的污染程度，构建污染源搜索模型；
[0056]
s2、利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，以根据所述最优搜索路径找到污染源的位置。
[0057]
具体地，所述污染源搜索模型具体为：
[0058][0059]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a∈{lv′
,a
，lv′
,na
，l
v,a
，l
v,na
}为无人机搜索污染源的决策集合，lv′
,a
表示当污染源的污染程度为第二污染程度时，进行污染源的定位搜索，l
v,na
表示当污染源的污染程度为第一污染程度时，不进行污染源的定位搜索，l
v,a
表示当污染源的污染程度为第一污染程度时，进行污染源的定位搜索，lv′
,na
表示当污染源的污染程度为第二污染程度时，不进行污染源的定位搜索。
[0060]
可以理解的是，预设范围具体为工业区，第一污染程度具体为重度污染、第二污染程度具体为轻度污染；在本发明实施例中，工业区内污染源的污染程度可以分为三个等级：无污染、轻度污染、重度污染，当一个工业区内存在多个不同污染程度的污染源时，无人机上的气体传感器在检测到多个污染源后会根据其污染程度、传播半径和距离进行优先级的判断，经过多次决策后计算出最优路径，其可简化为：搜索到多个污染源时，先对第一个污染源进行以下决策，即：“轻度污染”进行定位{lv′
,a
}、“轻度污染”不进行定位{lv′
,na
}、“重度污染”进行定位{l
v,a
}、“重度污染”不进行定位{l
v,na
}；然后再对下一污染源进行决策，最终决定出搜索和定位的路径。
[0061]
当气体传感器检测到大气污染物后未经过决策直接进行定位，期望收益为：
[0062]
ua＝pv[p
0va
+(1-p0)]+(1-pv′
)f
[0063]
当气体传感器检测到大气污染物后，先进行决策后再定位，其期望收益为：
[0064]
ua＝pv[p
0va
+(1-p0)u]+(1-pv′
)f-c
[0065]
当一个工业区内出现多个污染源且污染程度皆为轻度污染以上，气体传感器检测到某一污染物后决定“重度污染进行定位、轻度污染不进行定位”的期望收益为：
[0066]
ua＝[(p
0va
+(1-p0)u+f)p
v-f]pv′
+va(1-pv′
)pv′-c
[0067]
当气体传感器检测到某一污染物后决定“轻度污染进行定位、重度污染不进行定位”的期望收益为：
[0068]
ua＝[(p
0va
+(1-p0)u+f)pv′-f]+vapvpv′
(1-pv)-c
[0069]
可以理解的是，收益指的是根据污染源的情况所计算出的搜索路径，期望收益ua则指的是最优路径。在污染源搜索之后会根据计算出来的位置情况设定最优搜索路线。
[0070]
具体地，当满足预设的第一参数条件时，a＝lv′
,a
∪l
v,na
；当满足预设的第二参数条件时，a＝l
v,a
∪lv′
,na
；
[0071]
其中，所述第一参数条件具体为：
[0072]
c《[p
0va
+(1-p0)u+f](pv′-pv)+(1-pv′
)f
[0073]
[0074][0075]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a＝lv′
,a
∪l
v,na
表示当污染源的污染程度为第二污染程度时，进行污染源的定位搜索，且当污染源的污染程度为第一污染程度时，不进行污染源的定位搜索；
[0076]
所述第二参数条件具体为：
[0077]
c》(1-pv′
)f+(pv′-pv)[p
0va
+(1-p0)u+f]
[0078][0079][0080]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a＝l
v,a
∪lv′
,na
表示当污染源的污染程度为第一污染程度时，进行污染源的定位搜索，且当污染源的污染程度为第二污染程度时，不进行污染源的定位搜索。
[0081]
具体地，所述气体传感器成功定位污染源的概率包括以下情况中的至少一种：
[0082]
当污染源的污染程度为第一污染程度时，气体传感器成功定位污染源的概率为p0(f)＝[1-exp(l-f)]pv，式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv为污染源的污染程度为第一污染程度的概率；
[0083]
当污染源的污染程度为第二污染程度时，气体传感器成功定位污染源的概率为p0(f)＝[1-exp(l-f)]pv′
，式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv′
为污染源的污染程度为第二污染程度的概率。
[0084]
可以理解的是，由于污染源排放的污染物在经过空气传播后其浓度会有所下降，为了便于进行数值分析，当污染程度为重度污染时，无人机上的气体传感器成功定位污染源的概率：
[0085]
p0(f)＝[1-exp(l-f)]pv[0086]
式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv为污染源的污染程度为第一污染程度的概率；
[0087]
当污染程度为轻度污染时，气体传感器成功定位污染源的概率：
[0088]
p0(f)＝[1-exp(l-f)]pv′
[0089]
式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv′
为污染源的污染程度为第二污染程度的概率。
[0090]
可以理解的是，当气体传感器与污染源的距离越远，污染物经过传播后其浓度下
降得越多，成功定位污染源的概率越低。
[0091]
具体地，所述气体传感器检测污染物的成本，具体为：
[0092]
c(h)＝kh
[0093]
式中，k表示气体传感器对单个污染源进行定位的成本，h表示污染源分布的分散程度。
[0094]
可以理解的是，一般来说，当污染源分布的分散程度越高，无人机气体传感器检测与定位的成本则越高。
[0095]
具体地，所述利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，包括：
[0096]
步骤a，获取大气污染物的浓度γ＝m
tiβ
并初始化参数；其中，t为污染物检测次数，i为大气污染物检出限，β为大气污染物，β∈[0,1]，m为污染物传感器检测形式，m
t
＝m
t-1
+(b/m
t-1
×
ngen)，m0取值为0.01，b取值为0.025，ngen表示最大迭代次数；
[0097]
步骤b，在预设范围内随机初始化污染源位置，将所述污染源搜索模型作为目标函数，求解所述目标函数，计算当前污染源位置的最优值；
[0098]
步骤c，判断是否满足预设的终止条件，若是，输出最优值；若否，进入步骤d；
[0099]
步骤d，根据c1＝2exp(-(4t/ngen)2)计算收敛因子c1；其中，exp为指导函数，t为污染物检测次数，ngen表示最大迭代次数，d＝2-2t/ngen；
[0100]
步骤e，根据γ＝m
tiβ
计算预设范围内的每个污染源的大气污染物的浓度γi，根据a＝2dr-d计算收敛因子a；其中r∈[0,1]，d＝2-2t/ngen；
[0101]
步骤f，判断收敛因子c1是否大于1，若是，根据控制无人机进行全局搜索并获取最优置；若否，根据控制无人机进行局部搜索并获取最优置；其中，搜索并获取最优置；其中，搜索并获取最优置；其中，表示气体传感器在第i个污染源位置进行第t次污染物检测，r∈[0,1]，g
*
表示全局最优解，表示气体传感器检测第j维进行第t次污染物检测，表示气体传感器检测第k维进行第t次污染物检测，γi表示第i个污染源的污染物浓度，表示第j维气体传感器的位置，c3∈[0,1]，lbj、ubj为变量范围第j维的下界和上界；
[0102]
步骤g，执行t＝t+1，并返回步骤c。
[0103]
可以理解的，在蝴蝶算法中，设在某一个工业区内，有n个随机污染源，每个污染源都会排放大气污染物，大气污染源浓度的计算公式：
[0104]
γ＝m
tiβ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)式(1)中：
[0105]
排放的大气污染源的浓度为γ；
[0106]
污染物检测的次数为t；
[0107]
大气污染物检出限为i，其值由气体传感器的性能确定；
[0108]
大气污染物为β，β∈[0,1]；
[0109]
污染物传感器检测形式为m，计算公式如下：
[0110]mt
＝m
t-1
+(b/m
t-1
×
ngen)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0111]
式(2)中：m0取值为0.01；b取值为0.025；ngen表示最大迭代次数。
[0112]
当气体传感器检测到周围有浓度更高的大气污染物时，便会控制无人机朝向其移动，该阶段称为全局搜索，其计算公式为：
[0113][0114]
当气体传感器无法检测出周围的污染物时便进行随机移动，这一阶段称为局部搜索，其计算公式为：
[0115][0116]
式(3)、(4)中，表示气体传感器在第i个污染源位置进行第t次污染物检测；r∈[0,1]；g
*
表示全局最优解；γi表示第i个污染源的污染物浓度。
[0117]
为了进一步提高蝴蝶算法的探索能力和收敛的精度，分别将收敛因子a和c1引入蝴蝶算法的全局搜索和局部搜索公式中，则改进后的全局搜索计算公式为：
[0118][0119]
在污染物检测的前期a值的变化范围较大以增强全局搜索，在污染物检测的后期a值变化范围较小以增强局部搜索，在变化范围内，a值的计算公式如下：
[0120]
d＝2-2
×
t/ngen
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0121]
a＝2
×d×
r-d
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0122]
同理，将非线性收敛因子c1引入局部更新处，当c1》1时，算法进行全局搜索，否则进行局部搜索，c1会随着污染物检测的次数而增加，则改进后的局部搜索计算公式为：
[0123][0124]
式(8)中，c1的计算公式如下：
[0125]
c1＝2
×
exp(-(4
×
t/ngen)2)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0126]
式(8)(9)中，exp为指导函数，t为污染物检测次数，ngen表示最大迭代次数。
[0127]
经过双收敛因子的改进，发现局部搜索的能力有较大的提升，但是全局搜索的能力提升却不明显，这说明全局搜索和局部搜索并未达到有效平衡。借鉴于人工鱼群中算法中领导者策略以平衡全局搜索和局部搜索，领导者以食物为目标不断探索前进，因此考虑在局部搜索的过程中融合其领导者策略来平衡算法的全局搜索和局部搜索，领导者策略的计算公式如下：
[0128][0129]
式(10)中，表示第j维领导者(即气体传感器)的位置；c3∈[0,1]；lbj、ubj为变量范围第j维的下界和上界；
[0130]
因此在局部搜索过程中，蝴蝶算法的局部算子与人工鱼群中领导者机制同时进行
搜索，分别按式(8)和式(10)计算x1和x2即污染源搜索的位置结果，再按式(11)计算更新后的污染源位置：
[0131][0132]
本发明实施例所提供的大气污染物的污染源定位方法，通过根据无人机上的气体传感器检测到的预设范围内的污染源的污染程度，构建污染源搜索模型；利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，以根据所述最优搜索路径找到污染源的位置，能够提高大气污染源的定位效率，同时提高定位的精确度。由此可见，本发明实施例能够更好地平衡全局搜索和局部搜索，减少了气体传感器的能耗，加快了大气污染源的定位速度，提高了大气污染物的污染源定位的精确度。
[0133]
参见图3，图3是本发明实施例提供的一种大气污染物的污染源定位装置10的结构框图，所述大气污染物的污染源定位装置，包括：
[0134]
污染源搜索模型构建模块11，用于根据无人机上的气体传感器检测到的预设范围内的污染源的污染程度，构建污染源搜索模型；
[0135]
污染源定位模块12，用于利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，以根据所述最优搜索路径找到污染源的位置。
[0136]
优选地，所述污染源搜索模型具体为：
[0137][0138]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a∈{lv′
,a
，lv′
,na
，l
v,a
，l
v,na
}为无人机搜索污染源的决策集合，lv′
,a
表示当污染源的污染程度为第二污染程度时，进行污染源的定位搜索，l
v,na
表示当污染源的污染程度为第一污染程度时，不进行污染源的定位搜索，l
v,a
表示当污染源的污染程度为第一污染程度时，进行污染源的定位搜索，lv′
,na
表示当污染源的污染程度为第二污染程度时，不进行污染源的定位搜索。
[0139]
优选地，当满足预设的第一参数条件时，a＝lv′
,a
∪l
v,na
；当满足预设的第二参数条件时，a＝l
v,a
∪lv′
,na
；
[0140]
其中，所述第一参数条件具体为：
[0141]
c《[p
0va
+(1-p0)u+f](pv′-pv)+(1-pv′
)f
[0142][0143][0144]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a＝lv′
,a
∪l
v,na
表示当污染源的污染程度为
第二污染程度时，进行污染源的定位搜索，且当污染源的污染程度为第一污染程度时，不进行污染源的定位搜索；
[0145]
所述第二参数条件具体为：
[0146]
c》(1-pv′
)f+(pv′-pv)[p
0va
+(1-p0)u+f]
[0147][0148][0149]
式中，p0为气体传感器成功定位污染源的概率，va为预设范围内空气质量指数的下降值，u为无人机进行搜索污染源决策的次数，f为污染物的浓度随传播距离的增加而降低的值，pv为污染源的污染程度为第一污染程度的概率，pv′
为污染源的污染程度为第二污染程度的概率，c为气体传感器检测污染物的成本，a＝l
v,a
∪lv′
,na
表示当污染源的污染程度为第一污染程度时，进行污染源的定位搜索，且当污染源的污染程度为第二污染程度时，不进行污染源的定位搜索。
[0150]
优选地，所述气体传感器成功定位污染源的概率包括以下情况中的至少一种：
[0151]
当污染源的污染程度为第一污染程度时，气体传感器成功定位污染源的概率为p0(f)＝[1-exp(l-f)]pv，式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv为污染源的污染程度为第一污染程度的概率；
[0152]
当污染源的污染程度为第二污染程度时，气体传感器成功定位污染源的概率为p0(f)＝[1-exp(l-f)]pv′
，式中，f表示污染物的浓度随传播距离的增加而降低的值，l为无人机与所需定位的污染源的直线距离，pv′
为污染源的污染程度为第二污染程度的概率。
[0153]
优选地，所述气体传感器检测污染物的成本，具体为：
[0154]
c(h)＝kh
[0155]
式中，k表示气体传感器对单个污染源进行定位的成本，h表示污染源分布的分散程度。
[0156]
优选地，所述利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，包括：
[0157]
步骤a，获取大气污染物的浓度γ＝m
tiβ
并初始化参数；其中，t为污染物检测次数，i为大气污染物检出限，β为大气污染物，β∈[0,1]，m为污染物传感器检测形式，m
t
＝m
t-1
+(b/m
t-1
×
ngen)，m0取值为0.01，b取值为0.025，ngen表示最大迭代次数；
[0158]
步骤b，在预设范围内随机初始化污染源位置，将所述污染源搜索模型作为目标函数，求解所述目标函数，计算当前污染源位置的最优值；
[0159]
步骤c，判断是否满足预设的终止条件，若是，输出最优值；若否，进入步骤d；
[0160]
步骤d，根据c1＝2exp(-(4t/ngen)2)计算收敛因子c1；其中，exp为指导函数，t为污染物检测次数，ngen表示最大迭代次数，d＝2-2t/ngen；
[0161]
步骤e，根据γ＝m
tiβ
计算预设范围内的每个污染源的大气污染物的浓度γi，根据a＝2dr-d计算收敛因子a；其中r∈[0,1]，d＝2-2t/ngen；
[0162]
步骤f，判断收敛因子c1是否大于1，若是，根据控制
无人机进行全局搜索并获取最优置；若否，根据控制无人机进行局部搜索并获取最优置；其中，搜索并获取最优置；其中，搜索并获取最优置；其中，表示气体传感器在第i个污染源位置进行第t次污染物检测，r∈[0,1]，g
*
表示全局最优解，表示气体传感器检测第j维进行第t次污染物检测，表示气体传感器检测第k维进行第t次污染物检测，γi表示第i个污染源的污染物浓度，表示第j维气体传感器的位置，c3∈[0,1]，lbj、ubj为变量范围第j维的下界和上界；
[0163]
步骤g，执行t＝t+1，并返回步骤c。
[0164]
值得说明的是，本发明实施例所述的大气污染物的污染源定位装置10中各个模块的工作过程可参考上述实施例所述的大气污染物的污染源定位方法的工作过程，在此不再赘述。
[0165]
本发明实施例所提供的一种大气污染物的污染源定位装置10，通过根据无人机上的气体传感器检测到的预设范围内的污染源的污染程度，构建污染源搜索模型；利用改进的蝴蝶算法求解所述污染源搜索模型，得到污染源的最优搜索路径，以根据所述最优搜索路径找到污染源的位置，能够提高大气污染源的定位效率，同时提高定位的精确度。由此可见，本发明实施例能够更好地平衡全局搜索和局部搜索，减少了气体传感器的能耗，加快了大气污染源的定位速度，提高了大气污染物的污染源定位的精确度。
[0166]
本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序；其中，所述计算机程序在运行时控制所述计算机可读存储介质所在的设备执行如上述任一实施例所述的大气污染物的污染源定位方法。
[0167]
参见图4，图4是本发明实施例提供的一种设备20的结构框图，所述设备20包括：处理器21、存储器22以及存储在所述存储器22中并可在所述处理器21上运行的计算机程序。所述处理器21执行所述计算机程序时实现上述大气污染物的污染源定位方法实施例中的步骤。或者，所述处理器21执行所述计算机程序时实现上述各装置实施例中各模块/单元的功能。
[0168]
示例性的，所述计算机程序可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器22中，并由所述处理器21执行，以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序在所述设备20中的执行过程。
[0169]
所述设备20可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述设备20可包括，但不仅限于，处理器21、存储器22。本领域技术人员可以理解，所述示意图仅仅是设备20的示例，并不构成对设备20的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述设备20还可以包括输入输出设备、网络接入设备、总线等。
[0170]
所称处理器21可以是中央处理单元(central processing unit，cpu)，还可以是
其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述处理器21是所述设备20的控制中心，利用各种接口和线路连接整个设备20的各个部分。
[0171]
所述存储器22可用于存储所述计算机程序和/或模块，所述处理器21通过运行或执行存储在所述存储器22内的计算机程序和/或模块，以及调用存储在存储器22内的数据，实现所述设备20的各种功能。所述存储器22可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器22可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(smart media card,smc)，安全数字(secure digital,sd)卡，闪存卡(flash card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
[0172]
其中，所述设备20集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器21执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、u盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、电载波信号、电信信号以及软件分发介质等。
[0173]
需说明的是，以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外，本发明提供的装置实施例附图中，模块之间的连接关系表示它们之间具有通信连接，具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
[0174]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。