CN111505205A - 一种气体泄漏源强反算改进型搜索算法 - Google Patents
一种气体泄漏源强反算改进型搜索算法 Download PDFInfo
- Publication number
- CN111505205A CN111505205A CN202010063258.1A CN202010063258A CN111505205A CN 111505205 A CN111505205 A CN 111505205A CN 202010063258 A CN202010063258 A CN 202010063258A CN 111505205 A CN111505205 A CN 111505205A
- Authority
- CN
- China
- Prior art keywords
- leakage
- algorithm
- source
- calculation
- optimal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 43
- 238000010845 search algorithm Methods 0.000 title claims abstract description 15
- 241000167857 Bourreria Species 0.000 title claims abstract 5
- 239000002245 particle Substances 0.000 claims abstract description 49
- 238000009792 diffusion process Methods 0.000 claims abstract description 27
- 238000005259 measurement Methods 0.000 claims abstract description 25
- 238000012614 Monte-Carlo sampling Methods 0.000 claims abstract description 22
- 230000002068 genetic effect Effects 0.000 claims abstract description 18
- 238000005070 sampling Methods 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000002708 enhancing effect Effects 0.000 claims description 5
- 108090000623 proteins and genes Proteins 0.000 claims description 5
- 239000000779 smoke Substances 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000012544 monitoring process Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 244000257022 tick clover Species 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008786 sensory perception of smell Effects 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 241000208125 Nicotiana Species 0.000 description 1
- 235000002637 Nicotiana tabacum Nutrition 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000383 hazardous chemical Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0027—General constructional details of gas analysers, e.g. portable test equipment concerning the detector
- G01N33/0036—General constructional details of gas analysers, e.g. portable test equipment concerning the detector specially adapted to detect a particular component
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M3/00—Investigating fluid-tightness of structures
- G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0027—General constructional details of gas analysers, e.g. portable test equipment concerning the detector
- G01N33/0031—General constructional details of gas analysers, e.g. portable test equipment concerning the detector comprising two or more sensors, e.g. a sensor array
- G01N33/0034—General constructional details of gas analysers, e.g. portable test equipment concerning the detector comprising two or more sensors, e.g. a sensor array comprising neural networks or related mathematical techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display
- G01N33/0068—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display using a computer specifically programmed
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Medicinal Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Food Science & Technology (AREA)
- Mathematical Physics (AREA)
- Combustion & Propulsion (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biomedical Technology (AREA)
- Probability & Statistics with Applications (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种气体泄漏源强反算改进型搜索算法,包括以下步骤:使用传感器网络获取泄漏现场的浓度测量数据,先验扩散模型计算泄漏现场的浓度计算数据,通过建立测量数据和先验扩散模型间的函数关系,将反算问题转化为最优化求解问题,利用改进型遗传算法完成泄漏源参数预估值;将改进型遗传算法的预估值结果作为马尔科夫链蒙特卡洛抽样初始抽样点,利用马尔科夫链蒙特卡洛抽样算法得到泄漏源位置的概率分布范围;以马尔科夫链蒙特卡洛抽样算法提供的泄漏源位置的分布范围,使用移动机器人进场主动嗅探进行处理,利用改进型保守收敛粒子群算法用于完成泄漏源精确定位。先利用基于传感器网络的定位技术确定疑似的泄漏源区域,而后在疑似区域内有针对性的释放移动机器人进行进一步搜索,取长补短联合处理危害气体泄漏事故中源强反算问题。
Description
技术领域
本发明涉及危险化学品气体泄漏检测技术领域,尤其是一种气体泄漏源强 反算改进型搜索算法。
背景技术
危险化学品气体泄漏源位置和强度的确定是危险化学品气体泄漏事故应 急救援的基石,开展此领域的研究将具有其重要性和紧迫性。危险化学品气体 泄漏往往是突发性的,导致泄漏时间、泄漏地点和泄漏环境均具有不确定性。 在突发情况下,要求在信息缺失的有限条件下,以最短的时间确定未知泄漏源, 以进一步确定应急疏散区域和安全距离,为应急决策提供依据。危险化学品气 体泄漏发生后,利用监测设备对泄漏现场及周围进行气体浓度的监测,进而利 用所监测数据实现泄漏源位置的定位。泄漏源定位形式主要分为两种:传感器 网络监测计算定位和移动机器人协同搜索定位。
传感器网络监测计算定位的思路是融合先验的气体扩散模型和实时采集 的气体浓度数据对气源参数迭代求解的分布式估计问题。该方法的实现方式是 首先通过先验模型依据候选参数集得到模拟浓度,再通过贝叶斯理论将模拟浓 度和传感器网络测量浓度建立匹配函数,通过不断优化匹配函数结果,从而得 到气源参数的后验分布。
移动机器人协同搜索定位的思路是利用移动机器人安装相应的气体浓度 传感器形成的具有嗅觉感知的个体,故而应有于源强反算问题的移动机器人搜 索又被成为主动嗅觉问题,机器人利用风向信息和气体浓度信息,采用启发性 搜索策略实现定位搜索。基于传感器网络的泄漏源定位受限于传感器网络的数 据有限性和网络布置的固定性,基于移动机器人协同的泄漏源定位则可以充分 发挥移动机器人的运动灵活性以实现高细致更有针对性的局部区域搜索,从而 可以使得移动机器人的场景适应性更强,可以应对更为复杂动态的环境。
对于基于无线传感器网络的气体源监测定位技术存在如下的研究缺点:
(1)属于统计概率分布类别的数据处理算法需要基于事故信息以及待求 参数的先验信息,虽然可以有效地评估反演结果的不确定性,但是对于后验分 布的抽样过程极为耗时,且初始抽样点的选取直接关系到抽样结果的准确与 否;
(2)属于优化最优解类别的反算算法对先验扩散模型存在较强的依赖, 且由于源强反算问题的高度非线性,只能采用直接优化算法进行反算,但是以 遗传算法为主的直接优化算法需要极大的种群规模才能得到较为理想的结果, 而这将直接导致计算成本指数增加;
(3)传感器网络规模与搜索目标准确性之间的平衡问题:显而易见,要 获得精度更高的疑似源参数,必然需要增加传感器网络的规模,这是典型的成 本和成果的平衡。
对于基于移动机器人的气体源监测定位技术则存在如下技术缺点:
(1)湍流造成的烟羽动态分布特性。小尺度漩涡将气味/气体烟羽撕裂成 许多细丝,造成局部微小区域浓度分布数值振荡性,容易给移动机器人带来局 部极值陷阱;由于实际气体扩散场的破碎现象明显,故而移动机器人面对规模 较大的搜索空间的适应性较低;
(2)当前的指导多机器人***的算法效率普遍较低,高效的启发策略缺 失,并且少有研究涉及到多源场景下的并行定位,这极大地限制移动机器人协 同搜索的场景适应性。
发明内容
本发明的目的是提供一种气体泄漏源强反算改进型搜索算法,先利用基于 传感器网络的定位技术确定疑似的泄漏源区域,而后在疑似区域内有针对性的 释放移动机器人进行进一步搜索,取长补短联合处理危害气体泄漏事故中源强 反算问题。
为实现上述目的,本发明采用下述技术方案:
一种气体泄漏源强反算改进型搜索算法,包括以下步骤:
使用传感器网络获取泄漏现场的浓度测量数据,先验扩散模型计算泄漏现 场的浓度计算数据,通过建立测量数据和先验扩散模型间的函数关系,将反算 问题转化为最优化求解问题,利用改进型遗传算法完成泄漏源参数预估值;
将改进型遗传算法的预估值结果作为马尔科夫链蒙特卡洛抽样初始抽样 点,利用马尔科夫链蒙特卡洛抽样算法得到泄漏源位置的概率分布范围;
以马尔科夫链蒙特卡洛抽样算法提供的泄漏源位置的分布范围,使用移动 机器人进场主动嗅探进行处理,利用改进型保守收敛粒子群算法用于完成泄漏 源精确定位。
进一步地,所述使用传感器网络获取泄漏现场的浓度测量数据,先验扩散 模型计算泄漏现场的浓度计算数据,通过建立测量数据和先验扩散模型间的函 数关系,将反算问题转化为最优化求解问题,具体包括:
进一步地,所述利用改进型遗传算法完成泄漏源参数预估值,具体包括:
收集每步迭代在选择步骤中被丢弃的泄漏源单位样本个体,将被丢弃的泄 漏源单位样本个体单独存入淘汰者基因库用于扩大交叉的范围;
收集每步迭代最优的泄漏源单位样本个体,引入“跟随和搜索”策略,加强 种群向最优点方向演变概率,同时当最优点更新陷入停滞时,以当前最优点为 起点进行主动局部搜索操作,直到新产生的泄漏源粒子数目符合要求。
进一步地,所述利用马尔科夫链蒙特卡洛抽样算法得到泄漏源位置的概率 分布范围,具体包括:
使用自适应Metropolis算法作为马尔科夫链蒙特卡洛抽样的抽样算法,采 用“对数平衡”数据处理和“部分匹配,协同搜索”策略,直到满足抽样终止 条件。
进一步地,所述利用改进型保守收敛粒子群算法用于完成泄漏源精确定 位,具体包括:
初期按照标准粒子群算法运行,直到粒子群历史最优值更新停滞代数超过 设定限值;
从同代中适应度较小的子群中征调部分粒子与历史最优值粒子组成“最优 集团”,用于加强最优值点附近的局部搜索力度;
远离最优点重置部分适应度较小的粒子的位置和速度,以保证算法尽可能 地跳出局部极值;
动态更新最优集团,直到最优集团成功收缩定位或迭代代数达到上限。
发明内容中提供的效果仅仅是实施例的效果,而不是发明所有的全部效 果,上述技术方案中的一个技术方案具有如下优点或有益效果:
本发明按照“源强参数初步预估-源位置范围确定-泄漏源精准定位”的三 段式研究思路,提出更高效源强反算问题的搜索算法。初阶段使用传感器网络 获取泄漏现场的浓度测量数据,先验扩散模型计算泄漏现场的浓度计算数据, 通过建立测量数据和先验扩散模型间的函数关系,将反算问题转化为最优化求 解问题,而后利用优化算法求解最优源强参数。遗传算法用于源强反算问题存 在早熟收敛和计算成本高等缺点,本申请采用改进型遗传算法,可以较低的计 算成本下取得较优的参数解。通过引入“淘汰者基因库”扩充算法前期的搜索 空间,引入“跟随和搜索”策略强化收敛域内的局部搜索力度。
其次,以马尔科夫链蒙特卡洛抽样完成源位置范围的确定。改进型遗传算 法的结果作为马尔科夫链蒙特卡洛抽样初始抽样点,以解决马尔科夫链蒙特卡 洛抽样初始点依赖性问题,提出“部分匹配,协同搜索”的运算策略,来提升 其误差适应性。
基于移动机器人进场主动嗅探的方案进行处理,传统的粒子群算法存在搜 索速度较慢,定位结果偏离大的缺点。因此提出一种改进型保守收敛粒子群算 法来实现泄漏源精确定位。通过引入“征调和重置”策略使得算法具有更强的
附图说明
图1是本发明气体泄漏源强反算改进型搜索算法流程图;
图2是本发明改进型遗传算法MGA的新交叉流程图;
图3是图2中步骤S2的一个实施例;
图4是图2中步骤S3的一个实施例。
具体实施方式
为能清楚说明本方案的技术特点,下面通过具体实施方式,并结合其附图, 对本发明进行详细阐述。下文的公开提供了许多不同的实施例或例子用来实现 本发明的不同结构。为了简化本发明的公开,下文中对特定例子的部件和设置 进行描述。此外,本发明可以在不同例子中重复参考数字和/或字母。这种重 复是为了简化和清楚的目的,其本身不指示所讨论各种实施例和/或设置之间 的关系。应当注意,在附图中所图示的部件不一定按比例绘制。本发明省略了 对公知组件和处理技术及工艺的描述以避免不必要地限制本发明。
如图1所示,一种气体泄漏源强反算改进型搜索算法,包括以下步骤:
S1、使用传感器网络获取泄漏现场的浓度测量数据,先验扩散模型计算泄 漏现场的浓度计算数据,通过建立测量数据和先验扩散模型间的函数关系,将 反算问题转化为最优化求解问题,利用改进型遗传算法完成泄漏源参数预估 值;
S2、将改进型遗传算法的预估值结果作为马尔科夫链蒙特卡洛抽样初始抽 样点,利用马尔科夫链蒙特卡洛抽样算法得到泄漏源位置的概率分布范围;
S3、以马尔科夫链蒙特卡洛抽样算法提供的泄漏源位置的分布范围,使用 移动机器人进场主动嗅探进行处理,利用改进型保守收敛粒子群算法用于完成 泄漏源精确定位。
作为本申请的一个实施例,步骤S1中,使用传感器网络获取泄漏现场的 浓度测量数据,先验扩散模型计算泄漏现场的浓度计算数据,通过建立测量数 据和先验扩散模型间的函数关系,将反算问题转化为最优化求解问题,具体包 括:
式中,
Q:瞬时泄漏的物料质量,单位为kg;
u:风速;
t:泄漏时间;
如图2所示,作为本申请的一个实施例,改进型遗传算法MGA的新交叉流 程如下:
(1)从“选择”阶段保留下来的子群中随机选取两个不同个体,分别记 为P1,P2,即选择保留下来的泄漏物单位样本。
(2)产生随机数r1~U(0,1),记交叉率为α,如果r1<α,则P1,P2依旧按照GA 传统的交叉流程进行交叉操作;如果r1≥α,则转入步骤(3);
(3)记***当前的最优点更新停滞代数为S,设置的最优点停滞代数上 限为SetMax,如果S<SetMax,则进入步骤(4),否则转入步骤(5);
(4)从“选择”淘汰的淘汰者基因库中随机选取两个个体,分别记为 Ps1,Ps2,即从选择淘汰的泄漏物单位样本,用于分别匹配P1,P2进行交叉操作, 其中P1,P2作为母体进行交叉,记母系继承率为β;
(5)产生随机数r1~U(0,1),记追随率为γ,如果r1<γ,则进入步骤(6), 否则转入步骤(7);
(6)追随过程:记当前最优粒子为Pbest,P1,P2按照如下操作进行跟随,其 中f为随机数且f~N(0,1):
Pnew_i=Pbest+f×(Pbest-Pi)
(7)主动搜索:记当前最优粒子为Pbest,P1,P2按照如下操作进行主动搜索:
Pnew_i=β×Pi+Pbest×(1-β)
(8)重复步骤(1)-(7),直到新产生的泄漏物粒子数目得到要求。
相比于传统GA的交叉流程设计,MGA的主要改进为:收集每步迭代在 选择步骤中被丢弃的泄漏源单位样本个体,将被丢弃的泄漏源单位样本个体单 独存入淘汰者基因库用于扩大交叉的范围;
收集每步迭代最优的泄漏源单位样本个体,引入“跟随和搜索”策略,加强 种群向最优点方向演变概率,同时当最优点更新陷入停滞时,以当前最优点为 起点进行主动局部搜索操作,直到新产生的泄漏物粒子数目符合要求。
如图3所示,作为本申请的一个实施例,步骤S2中,利用马尔科夫链蒙 特卡洛抽样算法得到泄漏源位置的概率分布范围,具体包括:
S21、初始化马氏链抽样起点为θ0=(x0,y0,Q0),将第t步得到目标参数点记为 θt=(xt,yt,Qt),其中Xt,Yt,为测量位置的坐标。
MCMC方法将模型误差、测量误差以及参数的先验分布等通过概率密度函数 的方式表达,通过抽样结果的分布更合理地描述源强反算问题的不确定性。当 目标对象的搜索空间较复杂时,MCMC对初始点和提议分布的选择很敏感,不合 理的初始点会导致马氏链在收敛域外进行大量的无效抽样。因此采用MGA的预 估结果为MCMC的抽样初始起点。
S24、重复S22、S23直到抽样终止条件满足,如抽样代数达到要求。
如图4所示,作为本申请的一个实施例,步骤S3中,以马尔科夫链蒙特 卡洛抽样算法提供的泄漏源位置的分布范围,使用移动机器人进场主动嗅探进 行处理,利用改进型保守收敛粒子群算法用于完成泄漏源精确定位,具体包括:
S31、初始化粒子群的参数,分别为泄漏物样本的粒子总数、坐标点和粒 子初始速度;
S32、初期按标准粒子群算法运行,直到粒子群历史最优值更新停滞代数 超过设定限制;
式中,
Vi:粒子的飞翔速度;
Pi:粒子个体经历过的最好位置;
Pg:整个群体所有粒子经过的最好位置;
W:惯性因子,其值为非负;其值较大,全局寻优能力强,局部寻优能力弱; 其值较大,全局寻优能力弱,局部寻优能力强;因此,w越大,粒子飞翔的速 度和位置更新的幅度也就越大,因此偏离原先寻优轨道的程度也就越大;
ρ:范围参数;
C1,C2:加速常数,C1是根据个体自身的经验进行判断的常数,C2是根据群体 的经验;二者是学习因子,通常C1=C2=2;
r1,r2:为[0,1]范围内变换的随机数;
S35、根据已更新的速度更新粒子群的新位置,然后计算粒子群的新适应 值,更新粒子的自身历史最优,而后更新粒子群的全局历史最优,如果普通子 群中出现了新的全局最优值,则将当前已有的“最优集团”子群直接全部设置 为新的最优粒子的辅助粒子群,此为“最优集团”动态更新操作;
S36、如果普通子群中没有出现新的全局最优值,则根据“最优集团”搜 索新的最优值是否成功,更新“最优集团”搜索结果记录,根据搜索结果记录 是否满足连续成功或是连续失败次数更新局部搜索范围参数ρ;
S37、重复步骤S34、S35、S36直到“最优集团”成功收缩定位或是迭代 代数达到上限。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明 保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上, 本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发 明的保护范围以内。
Claims (6)
1.一种气体泄漏源强反算改进型搜索算法,其特征是,包括以下步骤:
使用传感器网络获取泄漏现场的浓度测量数据,利用先验扩散模型计算泄漏现场的浓度计算数据,通过建立测量数据和先验扩散模型间的函数关系,将反算问题转化为最优化求解问题,应用改进型遗传算法完成泄漏源参数预估值计算;
将改进型遗传算法的预估值结果作为马尔科夫链蒙特卡洛抽样算法的初始抽样点,利用马尔科夫链蒙特卡洛抽样算法得到泄漏源位置的概率分布范围;
以马尔科夫链蒙特卡洛抽样算法提供的泄漏源位置的分布范围,使用移动机器人在泄漏源分布范围内寻找泄漏源,利用改进型保守收敛粒子群算法用于完成泄漏源精确定位。
5.如权利要求4所述的气体泄漏源强反算改进型搜索算法,其特征是,所述利用马尔科夫链蒙特卡洛抽样算法得到泄漏源位置的概率分布范围,具体包括:
使用自适应Metropolis算法作为马尔科夫链蒙特卡洛抽样的抽样算法,采用“对数平衡”对候选点进行数据处理,直到满足抽样终止条件。
6.如权利要求5所述的气体泄漏源强反算改进型搜索算法,其特征是,所述利用改进型保守收敛粒子群算法用于完成泄漏源精确定位,具体包括:
初期按照标准粒子群算法运行,直到粒子群历史最优值更新停滞代数超过设定限值;
从同代中适应度较小的子群中征调部分粒子与历史最优值粒子组成“最优集团”,用于加强最优值点附近的局部搜索力度;
远离最优点重置部分适应度较小的粒子的位置和速度,以保证算法尽可能地跳出局部极值;
动态更新最优集团,直到最优集团成功收缩定位或迭代代数达到上限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010063258.1A CN111505205A (zh) | 2020-01-20 | 2020-01-20 | 一种气体泄漏源强反算改进型搜索算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010063258.1A CN111505205A (zh) | 2020-01-20 | 2020-01-20 | 一种气体泄漏源强反算改进型搜索算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111505205A true CN111505205A (zh) | 2020-08-07 |
Family
ID=71868919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010063258.1A Pending CN111505205A (zh) | 2020-01-20 | 2020-01-20 | 一种气体泄漏源强反算改进型搜索算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111505205A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113326655A (zh) * | 2021-05-25 | 2021-08-31 | 广西电网有限责任公司电力科学研究院 | 辐射型配电网可靠性与经济性综合评估方法及装置 |
CN113586968A (zh) * | 2021-06-01 | 2021-11-02 | 北京市燃气集团有限责任公司 | 一种天然气泄漏源定位方法及装置 |
CN114154383A (zh) * | 2021-12-01 | 2022-03-08 | 中国人民解放军国防科技大学 | 基于认知搜索策略的多机器人源搜索方案生成方法及*** |
WO2022153921A1 (ja) * | 2021-01-13 | 2022-07-21 | Arithmer株式会社 | モデル生成装置、モデル生成方法、拡散源位置出力装置、拡散範囲出力装置及び拡散源位置出力方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104834017A (zh) * | 2015-04-30 | 2015-08-12 | 华南理工大学 | 一种有毒有害气体泄漏扩散事故源定位方法 |
CN106779185A (zh) * | 2016-11-30 | 2017-05-31 | 江苏省邮电规划设计院有限责任公司 | 一种基于危化品扩散模型的交通应急指挥方法 |
CN107402283A (zh) * | 2016-05-20 | 2017-11-28 | 北京化工大学 | 一种借助移动机器人实现的气体泄漏溯源方法 |
CN108920818A (zh) * | 2018-06-29 | 2018-11-30 | 中国安全生产科学研究院 | 危险化学品气体泄漏等浓度线质心轨迹溯源算法 |
-
2020
- 2020-01-20 CN CN202010063258.1A patent/CN111505205A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104834017A (zh) * | 2015-04-30 | 2015-08-12 | 华南理工大学 | 一种有毒有害气体泄漏扩散事故源定位方法 |
CN107402283A (zh) * | 2016-05-20 | 2017-11-28 | 北京化工大学 | 一种借助移动机器人实现的气体泄漏溯源方法 |
CN106779185A (zh) * | 2016-11-30 | 2017-05-31 | 江苏省邮电规划设计院有限责任公司 | 一种基于危化品扩散模型的交通应急指挥方法 |
CN108920818A (zh) * | 2018-06-29 | 2018-11-30 | 中国安全生产科学研究院 | 危险化学品气体泄漏等浓度线质心轨迹溯源算法 |
Non-Patent Citations (1)
Title |
---|
张儒: "危险化学品事故泄漏源强反算算法研究", 《中国优秀硕博士学位论文全文数据库(硕士) 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022153921A1 (ja) * | 2021-01-13 | 2022-07-21 | Arithmer株式会社 | モデル生成装置、モデル生成方法、拡散源位置出力装置、拡散範囲出力装置及び拡散源位置出力方法 |
CN113326655A (zh) * | 2021-05-25 | 2021-08-31 | 广西电网有限责任公司电力科学研究院 | 辐射型配电网可靠性与经济性综合评估方法及装置 |
CN113586968A (zh) * | 2021-06-01 | 2021-11-02 | 北京市燃气集团有限责任公司 | 一种天然气泄漏源定位方法及装置 |
CN114154383A (zh) * | 2021-12-01 | 2022-03-08 | 中国人民解放军国防科技大学 | 基于认知搜索策略的多机器人源搜索方案生成方法及*** |
CN114154383B (zh) * | 2021-12-01 | 2024-07-12 | 中国人民解放军国防科技大学 | 基于认知搜索策略的多机器人源搜索方案生成方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111505205A (zh) | 一种气体泄漏源强反算改进型搜索算法 | |
WO2018218832A1 (zh) | 风电场的预测风速校正方法及装置 | |
CN110147610B (zh) | 一种用于河流突发水污染事件的溯源方法 | |
CN110232471B (zh) | 一种降水传感网节点布局优化方法及装置 | |
CN110796299A (zh) | 一种雷电预测方法 | |
US20110313958A1 (en) | System and method for empirical ensemble-based virtual sensing of particulates | |
CN111080070A (zh) | 一种基于空间误差的城市土地利用模拟元胞自动机方法 | |
CN111797572A (zh) | 一种城市事故灾害演化模拟及风险预测预警方法 | |
CN112270112A (zh) | 有毒有害气体泄露溯源新方法 | |
Lopez et al. | Effect of direction on wind speed estimation in complex terrain using neural networks | |
CN114662344B (zh) | 基于连续在线观测数据的大气污染溯源预测方法及*** | |
CN111428942A (zh) | 基于变网格技术提取微地形因子的线路覆冰厚度预测方法 | |
CN114611801A (zh) | 一种基于改进的鲸鱼优化算法的旅行商问题求解方法 | |
Park et al. | Source term estimation using deep reinforcement learning with Gaussian mixture model feature extraction for mobile sensors | |
CN114219345A (zh) | 一种基于数据挖掘的二次空气质量预报优化方法 | |
KR101313822B1 (ko) | 신경망 분석 기법을 이용한 풍력 밀도 예측 방법 | |
CN108594241B (zh) | 一种auv声隐身态势评估方法 | |
Emeksiz | Multi-gen genetic programming based improved innovative model for extrapolation of wind data at high altitudes, case study: Turkey | |
Fenech et al. | Selecting a global climate model for understanding future projections of climate change | |
Materukhin et al. | Simulation of Spatio-Temporal Data Streams from Geosensors Located On mobile Objects | |
CN110210072B (zh) | 基于近似模型及差分进化算法求解高维优化问题的方法 | |
CN114441397B (zh) | 基于高斯烟羽和质心定位的pm2.5污染源定位方法 | |
CN115034159A (zh) | 一种海上风电场的功率预测方法、装置、存储介质及*** | |
CN112434857B (zh) | 一种基于ga-pso优化的无人机机场网络选址布置方法 | |
CN113570120A (zh) | 基于改进粒子群算法的风电功率预测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200807 |
|
RJ01 | Rejection of invention patent application after publication |