基于定位周期的泄漏源定位方法、装置和计算机设备

摘要

本申请涉及一种基于定位周期的泄漏源定位方法、装置和计算机设备。包括：根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；监测当前周期的气体浓度信息，根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值；监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值；根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。上述方法能够根据目标函数模型计算值获得泄漏源强度和泄漏源位置，实现对危化气体泄漏源进行精确定位。

权利要求

1.一种基于定位周期的泄漏源定位方法，其特征在于，所述方法包括：

根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；

监测当前周期的气体浓度信息，根据所述当前周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述当前周期的目标函数模型计算值；

监测下一周期的气体浓度信息，根据所述下一周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述下一周期的目标函数模型计算值；

其中，所述监测当前周期的气体浓度信息和/或所述监测下一周期的气体浓度信息，是按照气体浓度升高的方向监测的；

根据所述当前周期的目标函数模型计算值和所述下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

2.根据权利要求1所述的方法，其特征在于，所述当前周期的气体浓度信息包括监测时间和与所述监测时间对应的气体浓度位置；

所述监测当前周期的气体浓度信息，根据所述当前周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述当前周期的目标函数模型计算值，包括：监测所述当前周期的气体浓度信息，根据所述当前周期的气体浓度信息和所述大气扩散模型，得到所述当前周期的大气扩散模型计算浓度序列；

根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到所述当前周期的目标函数模型计算值。

3.根据权利要求2所述的方法，其特征在于，所述监测所述当前周期的气体浓度信息，根据所述当前周期的气体浓度信息和所述大气扩散模型，得到所述当前周期的大气扩散模型计算浓度序列，包括：将获取到的与所述大气扩散模型相关的模型参数、所述监测时间和与所述监测时间对应的气体浓度位置输入所述大气扩散模型，得到所述当前周期的大气扩散模型计算浓度序列。

4.根据权利要求3所述的方法，其特征在于，所述将获取到的与所述大气扩散模型相关的模型参数、所述监测时间和与所述监测时间对应的气体浓度位置输入所述大气扩散模型，之前还包括：根据所述当前周期的气体在水平方向、竖直方向的扩散系数以及与所述监测时间对应的风速信息，得到与所述当前周期的大气扩散模型相关的模型参数。

5.根据权利要求2所述的方法，其特征在于，所述当前周期的气体浓度信息包括监测时间、与所述监测时间对应的气体浓度大小和与所述监测时间对应的气体浓度位置；

所述根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到所述当前周期的目标函数模型计算值，包括：根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列；其中，所述第一目标函数适应度序列同与所述监测时间对应的气体浓度大小和与所述监测时间对应的气体浓度位置存在第一对应关系；

对所述第一目标函数适应度序列进行排序及筛选，得到第二目标函数适应度序列；

对所述第二目标函数适应度序列执行预设次数遗传算法，得到中间种群适应度序列；

其中，所述中间种群适应度序列与所述第二目标函数适应度序列存在第二对应关系；

对所述中间种群适应度序列进行局部寻优，得到最优中间种群适应度；

根据所述第一对应关系、所述第二对应关系、所述最优中间种群适应度和所述目标函数模型，得到所述当前周期的目标函数模型计算值。

6.根据权利要求5所述的方法，其特征在于，所述根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列，包括：根据与所述监测时间对应的气体浓度大小和与所述监测时间对应的气体浓度位置，确定初始泄漏源强度集和初始泄漏源位置集；

根据所述初始泄漏源强度集和所述大气扩散模型计算浓度序列，得到初始大气扩散模型计算浓度集；

根据所述初始泄漏源强度集、所述初始泄漏源位置集、所述当前周期的气体浓度信息和所述初始大气扩散模型计算浓度集，得到所述第一目标函数适应度序列。

7.根据权利要求1所述的方法，其特征在于，所述根据所述当前周期的目标函数模型计算值和所述下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置，包括：将与所述当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与所述下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与所述下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置；

确定所述第一泄漏源强度和所述第二泄漏源强度之间的强度偏差值，确定所述第一泄漏源位置和所述第二泄漏源位置之间的距离偏差值；

若所述强度偏差值与所述距离偏差值的加权和小于或等于第一预设阈值，则将所述第二泄漏源强度确定为所述泄漏源强度，将所述第二泄漏源位置确定为所述泄漏源位置。

8.根据权利要求7所述的方法，其特征在于，所述方法还包括：

若所述强度偏差值与所述距离偏差值的加权和大于所述第一预设阈值，则将所述下一周期作为更新后的当前周期，并返回执行所述监测下一周期的气体浓度信息，根据所述下一周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述下一周期的目标函数模型计算值的步骤。

9.根据权利要求7所述的方法，其特征在于，所述将与所述当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与所述下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与所述下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置，包括：若所述当前周期的目标函数模型计算值小于或等于第二预设阈值，则将与所述当前周期的目标函数模型计算值对应的气体浓度大小确定为所述第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为所述第一泄漏源位置；

若所述下一周期的目标函数模型计算值小于或等于所述第二预设阈值，则将与所述下一周期的目标函数模型计算值对应的气体浓度大小确定为所述第二泄漏源强度，将与所述下一周期的目标函数模型计算值对应的气体浓度位置确定为所述第二泄漏源位置。

10.一种基于定位周期的泄漏源定位装置，其特征在于，所述装置包括：

构建模块，用于根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；

第一确定模块，用于监测当前周期的气体浓度信息，根据所述当前周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述当前周期的目标函数模型计算值；

第二确定模块，用于监测下一周期的气体浓度信息，根据所述下一周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述下一周期的目标函数模型计算值；

监测控制模块，用于按照气体浓度升高的方向监测所述当前周期的气体浓度信息和/或所述下一周期的气体浓度信息；

定位模块，用于根据所述当前周期的目标函数模型计算值和所述下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

11.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至9中任一项所述方法的步骤。

12.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至9中任一项所述的方法的步骤。

说明书

基于定位周期的泄漏源定位方法、装置和计算机设备

技术领域

[0001] 本申请涉及定位技术领域，特别涉及一种基于定位周期的泄漏源定位方法、装置和计算机设备。

背景技术

[0002] 危险化学品往往具有扩散性和毒害性，一旦发生泄漏，将造成严重的财产损失和人员伤亡。在危化品泄漏事故发生后，第一时间对危化品泄漏源进行定位，是事故应急处置、应急救援以及评估事故潜在危害的基础。危化品泄漏事故发生后，可以根据危化品泄漏事故中的气体泄漏情况来定位危化品泄漏源。

[0003] 传统的危化气体泄漏源定位方法是根据区域内危化气体的浓度梯度和扩散方向对危化气体泄漏源进行定位，常常会理想化泄漏事故现场环境。然而，危化品泄漏事故中往往伴随着爆炸、燃烧和复杂气流的发生，这些爆炸、燃烧和复杂气流的发生使得现场环境多数情况下并不理想化，非理想化的现场环境致使对危化气体泄漏源无法进行精确定位。

发明内容

[0004] 基于此，有必要针对上述技术问题，提供一种基于定位周期的泄漏源定位方法、装置和计算机设备。

[0005] 一种基于定位周期的泄漏源定位方法，所述方法包括：

[0006] 根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；

[0007] 监测当前周期的气体浓度信息，根据所述当前周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述当前周期的目标函数模型计算值；

[0008] 监测下一周期的气体浓度信息，根据所述下一周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述下一周期的目标函数模型计算值；

[0009] 根据所述当前周期的目标函数模型计算值和所述下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

[0010] 在其中一个实施例中，所述当前周期的气体浓度信息包括监测时间和与所述监测时间对应的气体浓度位置；

[0011] 所述监测当前周期的气体浓度信息，根据所述当前周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述当前周期的目标函数模型计算值，包括：

[0012] 监测所述当前周期的气体浓度信息，根据所述当前周期的气体浓度信息和所述大气扩散模型，得到所述当前周期的大气扩散模型计算浓度序列；

[0013] 根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到所述当前周期的目标函数模型计算值。

[0014] 在其中一个实施例中，所述监测所述当前周期的气体浓度信息，根据所述当前周期的气体浓度信息和所述大气扩散模型，得到所述当前周期的大气扩散模型计算浓度序列，包括：

[0015] 将获取到的与所述大气扩散模型相关的模型参数、所述监测时间和与所述监测时间对应的气体浓度位置输入所述大气扩散模型，得到所述当前周期的大气扩散模型计算浓度序列。

[0016] 在其中一个实施例中，所述将获取到的与所述大气扩散模型相关的模型参数、所述监测时间和与所述监测时间对应的气体浓度位置输入所述大气扩散模型，之前还包括：

[0017] 根据所述当前周期的气体在水平方向、竖直方向的扩散系数以及与所述监测时间对应的风速信息，得到与所述当前周期的大气扩散模型相关的模型参数。

[0018] 在其中一个实施例中，所述当前周期的气体浓度信息包括监测时间、与所述监测时间对应的气体浓度大小和与所述监测时间对应的气体浓度位置；

[0019] 所述根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到所述当前周期的目标函数模型计算值，包括：

[0020] 根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列；其中，所述第一目标函数适应度序列同与所述监测时间对应的气体浓度大小和与所述监测时间对应的气体浓度位置存在第一对应关系；

[0021] 对所述第一目标函数适应度序列进行排序及筛选，得到第二目标函数适应度序列；

[0022] 对所述第二目标函数适应度序列执行预设次数遗传算法，得到中间种群适应度序列；其中，所述中间种群适应度序列与所述第二目标函数适应度序列存在第二对应关系；

[0023] 对所述中间种群适应度序列进行局部寻优，得到最优中间种群适应度；

[0024] 根据所述第一对应关系、所述第二对应关系、所述最优中间种群适应度和所述目标函数模型，得到所述当前周期的目标函数模型计算值。

[0025] 在其中一个实施例中，所述根据所述目标函数模型、所述当前周期的气体浓度信息和所述当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列，包括：

[0026] 根据与所述监测时间对应的气体浓度大小和与所述监测时间对应的气体浓度位置，确定初始泄漏源强度集和初始泄漏源位置集；

[0027] 根据所述初始泄漏源强度集和所述大气扩散模型计算浓度序列，得到初始大气扩散模型计算浓度集；

[0028] 根据所述初始泄漏源强度集、所述初始泄漏源位置集、所述当前周期的气体浓度信息和所述初始大气扩散模型计算浓度集，得到所述第一目标函数适应度序列。

[0029] 在其中一个实施例中，所述根据所述当前周期的目标函数模型计算值和所述下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置，包括：

[0030] 将与所述当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与所述下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与所述下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置；

[0031] 确定所述第一泄漏源强度和所述第二泄漏源强度之间的强度偏差值，确定所述第一泄漏源位置和所述第二泄漏源位置之间的距离偏差值；

[0032] 若所述强度偏差值与所述距离偏差值的加权和小于或等于第一预设阈值，则将所述第二泄漏源强度确定为所述泄漏源强度，将所述第二泄漏源位置确定为所述泄漏源位置。

[0033] 在其中一个实施例中，所述方法还包括：

[0034] 若所述强度偏差值与所述距离偏差值的加权和大于所述第一预设阈值，则将所述下一周期作为更新后的当前周期，并返回执行所述监测下一周期的气体浓度信息，根据所述下一周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述下一周期的目标函数模型计算值的步骤。

[0035] 在其中一个实施例中，所述将与所述当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与所述下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与所述下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置，包括：

[0036] 若所述当前周期的目标函数模型计算值小于或等于第二预设阈值，则将与所述当前周期的目标函数模型计算值对应的气体浓度大小确定为所述第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为所述第一泄漏源位置；

[0037] 若所述下一周期的目标函数模型计算值小于或等于所述第二预设阈值，则将与所述下一周期的目标函数模型计算值对应的气体浓度大小确定为所述第二泄漏源强度，将与所述下一周期的目标函数模型计算值对应的气体浓度位置确定为所述第二泄漏源位置。

[0038] 一种基于定位周期的泄漏源定位装置，所述装置包括：

[0039] 构建模块，用于根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；

[0040] 第一确定模块，用于监测当前周期的气体浓度信息，根据所述当前周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述当前周期的目标函数模型计算值；

[0041] 第二确定模块，用于监测下一周期的气体浓度信息，根据所述下一周期的气体浓度信息、所述大气扩散模型和所述目标函数模型，得到所述下一周期的目标函数模型计算值；

[0042] 定位模块，用于根据所述当前周期的目标函数模型计算值和所述下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

[0043] 一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现上述任一项所述方法的步骤。

[0044] 一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述方法的步骤。

[0045] 上述基于定位周期的泄漏源定位方法、装置和计算机设备，通过根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型，监测当前周期的气体浓度信息，并根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值，监测下一周期的气体浓度信息，并根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值，进一步根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，实现对泄漏源强度和泄漏源位置的定位。上述方法能够根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值获得泄漏源强度和泄漏源位置，实现对危化气体泄漏源进行精确定位。

附图说明

[0046] 图1为一个实施例中基于定位周期的泄漏源定位方法的流程示意图；

[0047] 图2为一个实施例中步骤S200的细化步骤的流程示意图；

[0048] 图3为一个实施例中步骤S220的细化步骤的流程示意图；

[0049] 图4为一个实施例中步骤S221的细化步骤的流程示意图；

[0050] 图5为一个实施例中基于定位周期的泄漏源定位装置的结构框图；

[0051] 图6为一个实施例中计算机设备的内部结构图。

具体实施方式

[0052] 为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

[0053] 可以理解本申请中所使用的术语“第一”、“第二”等可在本文中用于描述各种序列关系及对应关系，但这些序列关系及对应关系不受这些术语限制。这些术语仅用于将一个序列关系及对应关系与另一个序列关系及对应关系区分开来。

[0054] 在一个实施例中，如图1所示，提供了一种基于定位周期的泄漏源定位方法，包括以下步骤：

[0055] 步骤S100，根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型。

[0056] 在本发明实施例中，大气扩散模型是一种计算模拟实际情况下的气体扩散迁移状况的数学模型。其中，大气扩散模型可以根据泄漏源所处的环境状态进行选择，例如，针对特殊气象条件和地形的扩散模型、封闭型扩散模型、熏烟型扩散模型、山区大气扩散模型和沿海大气扩散模式型，在此不做具体限定。

[0057] 其中，气体浓度信息指危化品泄漏源处释放出的气体在泄漏源周围空间中的浓度信息，此处的气体浓度信息是一种标识符号，不是实测的数据。

[0058] 步骤S200，监测当前周期的气体浓度信息，根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值。

[0059] 本发明实施例中，当监测到一定时长的监测数据之后才进行一次目标函数模型的计算，当前周期中的“周期”是指监测气体浓度信息的时间段，其中，并不对当前周期的时长做具体限定，可选地，当前周期的时长可以为1分钟到10分钟。当前周期的气体浓度信息包括监测到的监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置等信息。其中，监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置存在着一一对应关系。可选地，与监测时间对应的气体浓度大小是由离子迁移谱仪监测到的，与监测时间对应的气体浓度位置是由定位装置定位到的。离子迁移谱仪监测在其周围空间中不断地监测气体浓度大小，当离子迁移谱仪首次捕获到泄漏源气体的浓度信息时，则将首次捕获到泄漏源气体对应的监测时间视作当前周期的开始时间。其中，需要说明的是当前周期的气体浓度信息是离子迁移谱仪按照气体浓度升高的方向进行监测的。

[0060] 在本发明实施例中，目标函数模型计算值是指根据监测到的当前周期的气体浓度信息计算得到的目标函数模型的值。可选地，将前位置的气体浓度信息和大气扩散模型输入目标函数模型可以得到目标函数模型计算值。

[0061] 在本发明实施例中，假设周期为5分钟，离子迁移谱仪首次捕获到泄漏源气体浓度信息时的时间点为10:00，则当前周期为10:00-10:05，离子迁移谱仪按照气体浓度升高的方向进行监测，监测到10:05的时候，结束当前周期的监测，并得到当前周期的目标函数模型计算值。

[0062] 步骤S300，监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值。

[0063] 在本发明实施例中，下一周期中的“周期”与步骤S200中的“周期”相似，是指监测气体浓度信息的时间段，此处不做详细介绍。

[0064] 在本发明实施例中，假设下一周期为10:05-10:10，根据当前周期的目标函数模型计算值，会反推出一个与当前周期的气体浓度信息对应的泄漏源位置，离子迁移谱仪沿着反推出的与当前周期的气体浓度信息对应的泄漏源位置的方向移动并进行监测，监测到10:10的时候，结束所述下一周期的监测。

[0065] 其中，下一周期的气体浓度信息与步骤S200中的当前周期的气体浓度信息所监测的对象相同，下一周期的气体浓度信息也包括在下一周期中监测到的监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置。其中，监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置存在着一一对应关系。同样地，下一周期的目标函数模型计算值与步骤S200中当前周期的目标函数模型计算值的计算方式相同，在此不再赘述。

[0066] 步骤S400，根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

[0067] 具体地，可以根据当前周期的目标函数模型计算值和下一周期的目标函数计算值反推得到与当前周期的目标函数计算值对应和与下一周期的目标函数计算值对应的气体浓度大小和气体浓度位置。进一步，根据与当前周期的目标函数计算值对应的气体浓度大小和与下一周期的目标函数计算值对应的气体浓度大小之间的强度关系，以及与当前周期的目标函数计算值对应的气体浓度位置和与下一周期的目标函数计算值对应的气体浓度位置之间的距离关系，确定出泄漏源强度和泄漏源位置。

[0068] 上述实施例中，目标函数模型可表示为：

[0069]

[0070] 其中，L0表示待求的泄漏源位置，Q0表示待求的气体泄漏源强度，ti表示第i个时刻，Li表示时间为ti时监测点的位置Li＝(xi，yi)，CMol(Li，ti)表示在时间ti时、位置(xi，yi)处的大气扩散模型计算浓度，CIMS(Li，ti)表示在时间ti时、位置(xi，yi)处离子迁移谱仪监测到的气体浓度大小。可选地，采用全球定位系统(Global Position System，GPS)定位时(xi，yi)表示监测点在全球坐标系下的经度和纬度。

[0071] 上述基于定位周期的泄漏源定位方法，通过根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型，监测当前周期的气体浓度信息，并根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值，监测下一周期的气体浓度信息，并根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值，进一步根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，实现对泄漏源强度和泄漏源位置的定位。上述方法能够根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值获得泄漏源强度和泄漏源位置，实现对危化气体泄漏源进行精确定位。

[0072] 在其中一个实施例中，如图2所示，为步骤S200的细化步骤的流程示意图，其中，当前周期的气体浓度信息包括监测时间和与监测时间对应的气体浓度位置；步骤S200，监测当前周期的气体浓度信息，根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值，具体包括以下步骤：

[0073] 步骤S210，监测当前周期的气体浓度信息，根据当前周期的气体浓度信息和大气扩散模型，得到当前周期的大气扩散模型计算浓度序列。

[0074] 其中，大气扩散模型计算浓度序列是指根据已知的当前周期的气体浓度信息和大气扩散模型计算得到的关于大气扩散模型的序列。可选地，将当前周期的气体浓度信息输入大气扩散模型，可以得到大气扩散模型计算浓度序列。

[0075] 步骤S220，根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到当前周期的目标函数模型计算值。

[0076] 在本发明实施例中，将监测到的当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列输入目标函数模型，得到的是一个包含若干函数模型计算值的序列，对该序列进行排序和搜索，可以得到一个最优的函数模型计算值，将这个最优的函数模型计算值作为上述目标函数模型计算值。

[0077] 上述实施例中，计算得到目标函数模型计算值是定位泄漏源的关键步骤之一，获取目标函数模型计算值为定位泄漏源提供基础。

[0078] 在其中一个实施例中，步骤S210，监测当前周期的气体浓度信息，根据当前周期的气体浓度信息和大气扩散模型，得到当前周期的大气扩散模型计算浓度序列，包括：

[0079] 将获取到的与大气扩散模型相关的模型参数、监测时间和与监测时间对应的气体浓度位置输入大气扩散模型，得到当前周期的大气扩散模型计算浓度序列。

[0080] 可选地，大气扩散模型可以为高斯烟团模型。高斯烟团模型可表示为:

[0081]

[0082] 其中，C(x，y，t)表示时间为t时位置(x，y)处的大气模型计算浓度，v为风速，σx、σy、和σz分别为气体在x(水平轴)、y(水平轴)、z(竖直方向)方向的扩散系数。其中风速v包含风的速度和方向等信息，由超声波风速仪测得。

[0083] 上述实施例中，选取高斯烟团模型作为实施对象，将与高斯烟团模型相关的模型参数(例如，气体在水平方向、竖直方向的扩散系数以及与监测时间对应的风速信息等)、监测时间和与监测时间对应的气体浓度位置输入高斯烟团模型，得到高斯烟团模型计算浓度序列，此时的高斯烟团模型计算浓度序列中仅包含Q0一个未知参量，其中，将仅包含未知参量Q0的高斯烟团模型计算浓度序列作为大气扩散模型计算浓度序列，将大气扩散模型计算浓度序列作为确定目标函数中的泄漏源强度和泄漏源位置的基础。

[0084] 在其中一个实施例中，将获取到的与大气扩散模型相关的模型参数、监测时间和与监测时间对应的气体浓度位置输入大气扩散模型，之前还包括：

[0085] 根据当前周期的气体在水平方向、竖直方向的扩散系数以及与监测时间对应的风速信息，得到与当前周期的大气扩散模型相关的模型参数。其中，与监测时间对应的风速信息是指在当前监测时间下，与监测时间对应的气体浓度位置处的风的速度和方向。

[0086] 上述实施例中，大气扩散模型与当前气体在水平方向和竖直方向的扩散成都相关，并且与监测点当时的风向和风速相关，因此，在计算过程中，需要先获取与大气扩散模型相关的模型参数。

[0087] 在其中一个实施例中，如图3所示，为步骤S220的细化步骤的流程示意图，其中，当前周期的气体浓度信息包括监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置；步骤S220，根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到当前周期的目标函数模型计算值，具体包括以下步骤：

[0088] 步骤S221，根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列。

[0089] 其中，第一目标函数适应度序列同与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置存在第一对应关系。

[0090] 具体地，将监测到的当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列输入目标函数模型，可以得到第一目标函数适应度序列，因此第一目标函数适应度序列同当前周期的气体浓度信息(监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置)之间存在着一个一一对应的关系，可以把这个一一对应的关系称作第一对应关系。例如，当前周期的气体浓度信息序列可以表示为Aa1、Aa2……Aan，第一目标函数适应度序列可以表示为Tt1、Tt2……Ttn，其中，an等于tn，Aa1与Tt1对应，Aa2与Tt2对应……Aan与Ttn对应。

[0091] 步骤S222，对第一目标函数适应度序列进行排序及筛选，得到第二目标函数适应度序列。

[0092] 将第一目标函数适应度序列进行排序，剔除极值后，得到第二目标函数适应度序列，此时的第二目标函数适应度序列与第一目标函数适应度序列中的元素相对应。因此，第二目标函数适应度序列同当前周期的气体浓度信息(监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置)之间也存在着一个对应的关系。例如，得到的第二目标函数适应度序列可以表示为Ttt1、Ttt2……Tttn，其中，ttn小于等于tn，Ttt1、Ttt2……Tttn中的每一元素均与第一目标函数适应度序列Tt1、Tt2……Ttn中的唯一一个元素相对应，而第一目标函数适应度序列Tt1、Tt2……Ttn中的每个元素也与当前周期的气体浓度信息序列Aa1、Aa2……Aan中的唯一一个元素相对应，因此，Ttt1、Ttt2……Tttn中的每一元素均与当前周期的气体浓度信息序列Aa1、Aa2……Aan中的唯一一个元素相对应。

[0093] 步骤S223，对第二目标函数适应度序列执行预设次数遗传算法，得到中间种群适应度序列。

[0094] 其中，中间种群适应度序列与第二目标函数适应度序列存在第二对应关系。

[0095] 可选地，在本步骤中当执行了60代遗传算法迭代后，得到一个中间种群，在实际操作中对执行遗传算法的预设次数并不做具体限定。其中，得到的中间种群适应度序列与第二目标函数适应度序列中的元素相对应，也存在着一个对应的关系，而第二目标函数适应度序列与第一目标函数适应度序列中的元素相对应，第一目标函数适应度序列与当前周期的气体浓度信息(监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置)之间也存在着一个对应的关系。因此，中间种群适应度与当前周期的气体浓度信息(监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置)也存在着一个一一对应的关系。例如，得到的中间种群适应度序列可以表示为Mm1、Mm2……Mmn，其中mn小于等于ttn，ttn小于等于tn，Mm1、Mm2……Mmn中的每一元素均与第二目标函数适应度序列Ttt1、Ttt2……Tttn中的唯一一个元素相对应，Ttt1、Ttt2……Tttn中的每一元素均与第一目标函数适应度序列Tt1、Tt2……Ttn中的唯一一个元素相对应，也与当前周期的气体浓度信息序列Aa1、Aa2……Aan中的唯一一个元素相对应，因此，Mm1、Mm2……Mmn中的每一元素均与当前周期气体浓度信息序列Aa1、Aa2……Aan中的唯一一个元素相对应。

[0096] 步骤S224，对中间种群适应度序列进行局部寻优，得到最优中间种群适应度。

[0097] 在本发明实施例中，对中间种群适应度序列进行排序，选取中间种群适应度处于预设比例的个体执行模式搜索算法进行局部寻优，得到最优的中间种群适应度，由步骤S223中的解释可知，此时的最优中间种群适应度与当前周期气体浓度信息序列Aa1、Aa2……Aan中的唯一一个元素相对应，即对应唯一一组与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置。可选地，中间种群适应度处于预设比例的个体可以根据具体实施情况设定，可以为处于前20％的中间种群适应度。

[0098] 步骤S225，根据第一对应关系、第二对应关系、最优中间种群适应度和目标函数模型，得到目标函数模型计算值。

[0099] 其中，步骤S224得到最优中间种群适应度后，因为最优中间种群适应度包含于中间种群适应度序列中，且中间种群适应度序列与第二目标函数适应度序列存在第二对应关系，则根据第二对应关系可以从第二目标函数适应度序列中查找出与最优中间种群适应度对应的第二目标函数适应度中的元素(其中，为了便于后续描述，将第二目标函数适应度序列中与最优中间种群适应度对应的第二目标函数适应度中的元素统称为第二目标元素)，进一步地，因为第二目标函数适应度序列与第一目标函数适应度序列中的元素相对应，则可以从第一目标函数适应度序列中查找出与第一目标函数适应度中的元素(其中，为了便于后续描述，将第一目标函数适应度序列中与第一目标元素对应的元素统称为第一目标元素)，进一步地，因为第一目标函数适应度序列同与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置之间存在第一对应关系，则根据第一对应关系，可以确定与第一目标元素对应的当前周期气体浓度信息(监测时间、与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置)，将与第一目标元素对应的当前周期气体浓度信息输入目标函数模型，得到确定的目标函数模型计算值。

[0100] 上述实施例中，根据目标函数模型、气体的浓度大小和大气扩散模型计算浓度序列，得到第一目标函数适应度序列，并对第一目标函数适应度序列进行预处理后执行预设次数的遗传算法，得到中间种群适应度序列，进一步，运用局部寻优得到最优中间种群适应度。根据第一对应关系、第二对应关系、最优中间种群适应度和目标函数模型得到特定的与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置对应的目标函数模型计算值。

[0101] 在其中一个实施例中，如图4所示，为步骤S221的细化步骤的流程示意图，根据目标函数模型、当前周期的气体浓度信息和大气扩散模型计算浓度序列，得到第一目标函数适应度序列，具体包括以下步骤：

[0102] 步骤S2211，根据与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置，确定初始泄漏源强度集和初始泄漏源位置集。

[0103] 其中，对泄漏源位置L＝(x，y)和泄漏源强度Q0进行实数编码，并进行随机初始化，可以得到初始泄漏源强度集和初始泄漏源位置集。可选地，对泄漏源位置L＝(x，y)和泄漏源强度Q0进行实数编码和随机初始化这一步骤，是根据与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置进行实数编码和随机初始化的，此时与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置是指监测到当前周期的气体浓度信息后，与任意监测时间对应的气体浓度大小和与任意监测时间对应的气体浓度位置。

[0104] 步骤S2212，根据初始泄漏源强度集和大气扩散模型计算浓度序列，得到初始大气扩散模型计算浓度集。

[0105] 其中，大气扩散模型计算浓度序列仅包含未知参量泄漏源强度Q0，步骤S2211中对泄漏源强度Q0进行了初始化，将初始化后的Q0输入大气扩散模型计算浓度序列，能够得到不含未知参量的初始大气扩散模型计算浓度集。

[0106] 步骤S2213，根据初始泄漏源强度集、初始泄漏源位置集、当前周期的气体浓度信息和初始大气扩散模型计算浓度集，得到第一目标函数适应度序列。

[0107] 将泄漏源强度集、初始泄漏源位置集、当前周期的气体浓度信息和初始大气扩散模型计算浓度集输入目标函数模型式(1)中，得到第一目标函数适应度序列。

[0108] 上述实施例中，得到第一目标函数适应度序列为确定目标函数中的泄漏源强度和泄漏源位置提供基础。

[0109] 在其中一个实施例中，步骤S400，根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置，包括：

[0110] 将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置。

[0111] 在本发明实施例中，离子迁移谱仪以列维分布生成的路径进行对泄漏源气体进行监测，一旦监测到气体浓度信息，就将得到气体浓度信息的第一时刻作为当前周期的初始时刻，进一步，按照气体浓度升高的方向进行监测，直至时刻tn，其中n为所确定的定位周期的长度，n可以为1分钟到10分钟。由当前周期监测到的数据可以得到与当前周期的目标函数模型计算值，并将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置。

[0112] 可选地，下一周期的确定方式为以tn的监测位置为基准，向着第一泄漏源位置行进，直至时刻t2n，tn至t2n为下一周期的监测时段。

[0113] 可选地，确定第一泄漏源强度和第二泄漏源强度之间的强度偏差值，确定第一泄漏源位置和第二泄漏源位置之间的距离偏差值；若强度偏差值与距离偏差值的加权和小于或等于第一预设阈值，则将第二泄漏源强度确定为泄漏源强度，将第二泄漏源位置确定为泄漏源位置。

[0114] 其中，确定泄漏源强度和泄漏源位置的依据为，当强度偏差值与距离偏差值的加权和小于或等于一个预设的阈值时，说明两次反推出的泄漏源位置非常相近，表明泄漏源估算结果已足够稳定，则将第二泄漏源强度确定为泄漏源强度，将第二泄漏源位置确定为泄漏源位置。

[0115] 可选地，第一泄漏源强度与第二泄漏源强度之间的强度偏差值可以表示为其中， 表示第二泄漏源强度， 表示第一泄漏源强度，ΔQk表示第二泄漏源强度与第一泄漏源强度的差值的绝对值，将ΔQk称为强度偏差值。其中，第一泄漏源位置与第二泄漏源位置之间的距离偏差值可以表示为
其中，ΔLk表示距离偏差值， 表示表示第
二泄漏源位置， 表示第一泄漏源位置。

[0116] 可选地，强度偏差值与距离偏差值的加权和可以表示为ΔSk＝ξΔLk+(1-ξ)ΔQk，其中，ΔSk表示强度偏差值与距离偏差值的加权和，ξ为预设的权值参数，ξ的大小为0到1，(当ξ大于0.5时，表示泄漏源位置的权重大于泄漏源强度的权重；当ξ小于0.5时，表示泄漏源强度的权重大于泄漏源位置的权重；当ξ等于0.5时，表示泄漏源位置的权重与泄漏源强度的权重相当)，可以根据需求对ξ进行设置。当ΔSk低于某一误差限δ时，表明反推结果非常接近真实泄漏源的实际信息，则将第二泄漏源位置所对应的气体浓度大小确定为预估的泄漏源强度，将第二泄漏源位置确定为预估的泄漏源位置。

[0117] 可选地，若上述加权和大于第一预设阈值，则将所述下一周期作为更新后的当前周期，并返回执行监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值的步骤。

[0118] 在本发明实施例中，若加权和大于第一预设阈值时，说明两次反推出的泄漏源位置相距很远，表明泄漏源估算结果还不够稳定。此时，则将所述下一周期的气体浓度信息作为更新后的当前周期的气体浓度信息，将所述下一周期的目标函数模型计算值作为更新后的当前周期的目标函数模型计算值，下一周期确定的第二泄漏源强度和第二泄漏源位置作为更新后的第一泄漏源强度和第一泄漏源位置。上述信息更新之后，离子迁移谱仪以当前所处位置为基准，向着更新前的第一泄漏源位置和更新后的第一泄漏源位置(即更新前的第二泄漏源位置)的中点位置行进，重新监测气体浓度信息并执行确定目标函数模型计算值等步骤，直到得到泄漏源强度和泄漏源位置。

[0119] 其中，为了方便理解，对上述方案举例说明，例如，第一周期时段为0到tn，第二周期时段为tn到t2n，第三周期时段为t2n到t3n，第四周期时段为t3n到t4n……第k周期时段为tkn-1到tkn。以第一周期为当前周期，第二周期为下一周期，此时，离子迁移谱仪从0时刻对应的位置，沿着气体浓度升高的方向行进至tn时刻，用0到tn时段内监测到的气体浓度信息反推出一个泄漏源强度和泄漏源位置，可以称作第一周期泄漏源强度和第一周期泄漏源位置。接下来，离子迁移谱仪以tn时刻对应的位置为起点，以第一周期泄漏源位置为终点，行进至t2n时刻，用tn到t2n时段内监测到的气体浓度信息反推出一个泄漏源强度和泄漏源位置，可以称作第二周期泄漏源强度和第二周期泄漏源位置。此时，若第一周期与第二周期确定的强度偏差值与距离偏差值的加权和(其中，强度偏差值表示第二周期确定的泄漏源强度与第一周期确定的漏源强度差值的绝对值，距离偏差值表示第二周期确定的泄漏源位置与第一周期确定的漏源位置的距离)小于或等于第一预设阈值，则确定第二周期泄漏源强度为泄漏源强度，确定第二周期泄漏源位置为泄漏源位置。

[0120] 若第一周期与第二周期确定的强度偏差值与距离偏差值的加权和大于第一预设阈值，则返回执行监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值的步骤。具体地：以第二周期为更新后的当前周期，第三周期为更新后的下一周期，离子迁移谱仪以t2n时刻对应的位置为起点，以第一周期泄漏源位置和第二周期泄漏源位置的中点为终点，行进至t3n时刻，用t2n到t3n时段内监测到的气体浓度信息反推出一个泄漏源强度和泄漏源位置，可以称作第三周期泄漏源强度和第三周期泄漏源位置。

[0121] 此时，若第二周期与第三周期确定的强度偏差值与距离偏差值的加权和小于或等于第一预设阈值，则确定第三周期泄漏源强度为泄漏源强度，确定第三周期泄漏源位置为泄漏源位置。

[0122] 若第二周期与第三周期确定的强度偏差值与距离偏差值的加权和大于第一预设阈值，则返回执行监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值的步骤。具体地：以第三周期为更新后的当前周期，第四周期为更新后的下一周期，此时，离子迁移谱仪以t3n时刻对应的位置为起点，以第一周期泄漏源位置、第二周期泄漏源位置和第三周期泄漏源位置的中心点为终点，行进至t4n时刻，用t3n到t4n时段内监测到的气体浓度信息反推出一个泄漏源强度和泄漏源位置，可以称作第四周期泄漏源强度和第四周期泄漏源位置。

[0123] 此时，若第三周期与第四周期确定的强度偏差值与距离偏差值的加权和小于或等于第一预设阈值，则确定第四周期泄漏源强度为泄漏源强度，确定第四周期泄漏源位置为泄漏源位置。

[0124] 若第三周期与第四周期确定的强度偏差值与距离偏差值的加权和大于第一预设阈值，则返回执行监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值的步骤。具体地，按照上述方式依次更新，将以第四周期为更新后的当前周期，第五周期为更新后的下一周期进行数据监测和反推，在此不再赘述。

[0125] 上述实施例中，能够根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，获得泄漏源强度和泄漏源位置，实现对危化气体泄漏源进行精确定位，并且这种当得到一定时长的监测数据之后才进行一次泄漏源定位，可以大幅降低了运算量，并避免了因环境噪声引起的监测数据非平稳波动对泄漏源定位反推的不利影响。

[0126] 在其中一个实施例中，将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置，包括：

[0127] 若当前周期的目标函数模型计算值小于或等于第二预设阈值，则将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置。

[0128] 在本发明实施例中，当前周期的目标函数模型计算值小于或等于第二预设阈值为一个预设的收敛条件，表示大气扩散模型计算浓度CMol(Li，ti)与离子迁移谱仪监测浓度CIMS(Li，ti)差的平方和小于或等于设定的误差限，平方和小于或等于误差限表明反推结果非常接近真实泄漏源的实际信息，则将满足该收敛条件所对应的气体浓度大小确定为第一泄漏源强度，将满足该收敛条件所对应的气体浓度位置确定为第一泄漏源位置。

[0129] 可选地，若当前周期的目标函数模型计算值大于第二预设阈值，则返回执行监测当前周期气体浓度信息的步骤，直到当前周期的目标函数模型计算值小于或等于第二预设阈值，则将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置。

[0130] 具体地，若当前周期的目标函数模型计算值大于第二预设阈值，则更新当前周期的气体浓度信息，其中，更新当前周期的气体浓度信息是指重新采集或者更新用于计算的当前周期的气体浓度信息的采集时间和位置，并将新的当前周期的气体浓度信息输入预设大气扩散模型，得到更新后的大气扩散模型计算浓度序列，并将更新后的当前周期的气体浓度信息和更新后的大气扩散模型计算浓度序列输入目标函数模型，得到更新后的目标函数模型计算值。比较更新后的目标函数模型计算值与第二预设阈值的大小关系，直到更新后的目标函数模型计算值小于或等于第二预设阈值，则将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置。可选地，若下一周期的目标函数模型计算值小于或等于第二预设阈值，则将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置。在本发明实施例中，第二泄漏源强度、第二泄漏源位置的确定方式与第一泄漏源强度、第一泄漏源位置的方式相同，在此不再赘述。

[0131] 上述实施例中，能够根据当前周期的目标函数模型计算值和预设收敛条件获得第一泄漏源强度和第一泄漏源位置，并在当前周期的目标函数模型计算值不满足预设收敛条件的情况下不断修正定位方法中涉及模型的参数，获得第一泄漏源强度和第一泄漏源位置。

[0132] 在一个实施例中，如图5所示，提供了一种基于定位周期的泄漏源定位装置，包括：构建模块501、第一确定模块502、第二确定模块503和定位模块504，其中：

[0133] 构建模块501，用于根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；

[0134] 第一确定模块502，用于监测当前周期的气体浓度信息，根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值；

[0135] 第二确定模块503，用于监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值；

[0136] 定位模块504，用于根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

[0137] 在其中一个实施例中，第一确定模块502还用于执行监测当前周期的气体浓度信息，根据当前周期的气体浓度信息和大气扩散模型，得到当前周期的大气扩散模型计算浓度序列；根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到当前周期的目标函数模型计算值。

[0138] 在其中一个实施例中，第一确定模块502还用于执行将获取到的与大气扩散模型相关的模型参数、监测时间和与监测时间对应的气体浓度位置输入大气扩散模型，得到当前周期的大气扩散模型计算浓度序列。

[0139] 在其中一个实施例中，第一确定模块502还用于执行根据当前周期的气体在水平方向、竖直方向的扩散系数以及与监测时间对应的风速信息，得到与当前周期的大气扩散模型相关的模型参数。

[0140] 在其中一个实施例中，第一确定模块502还用于执行根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列；其中，第一目标函数适应度序列同与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置存在第一对应关系；对第一目标函数适应度序列进行排序及筛选，得到第二目标函数适应度序列；对第二目标函数适应度序列执行预设次数遗传算法，得到中间种群适应度序列；其中，中间种群适应度序列与第二目标函数适应度序列存在第二对应关系；对中间种群适应度序列进行局部寻优，得到最优中间种群适应度；根据第一对应关系、第二对应关系、最优中间种群适应度和目标函数模型，得到当前周期的目标函数模型计算值。

[0141] 在其中一个实施例中，第一确定模块502还用于执行根据与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置，确定初始泄漏源强度集和初始泄漏源位置集；根据初始泄漏源强度集和大气扩散模型计算浓度序列，得到初始大气扩散模型计算浓度集；根据初始泄漏源强度集、初始泄漏源位置集、当前周期的气体浓度信息和初始大气扩散模型计算浓度集，得到第一目标函数适应度序列。

[0142] 在其中一个实施例中，定位模块504还用于执行将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置；确定第一泄漏源强度和第二泄漏源强度之间的强度偏差值，确定第一泄漏源位置和第二泄漏源位置之间的距离偏差值；若强度偏差值与距离偏差值的加权和小于或等于第一预设阈值，则将第二泄漏源强度确定为泄漏源强度，将第二泄漏源位置确定为泄漏源位置。

[0143] 在其中一个实施例中，定位模块504还用于执行若强度偏差值与距离偏差值的加权和大于第一预设阈值，则将下一周期作为更新后的当前周期，并返回执行监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值的步骤。

[0144] 在其中一个实施例中，定位模块504还用于执行若当前周期的目标函数模型计算值小于或等于第二预设阈值，则将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置；若下一周期的目标函数模型计算值小于或等于第二预设阈值，则将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置。

[0145] 关于基于定位周期的泄漏源定位装置的具体限定可以参见上文中对于基于定位周期的泄漏源定位方法的限定，在此不再赘述。上述基于定位周期的泄漏源定位装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。

[0146] 在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图6所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储泄漏源定位数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于定位周期的泄漏源定位方法。

[0147] 本领域技术人员可以理解，图6中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

[0148] 在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器存储有计算机程序，处理器执行计算机程序时实现以下步骤：

[0149] 根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；

[0150] 监测当前周期的气体浓度信息，根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值；

[0151] 监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值；

[0152] 根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

[0153] 在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：监测当前周期的气体浓度信息，根据当前周期的气体浓度信息和大气扩散模型，得到当前周期的大气扩散模型计算浓度序列；根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到当前周期的目标函数模型计算值。

[0154] 在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：将获取到的与大气扩散模型相关的模型参数、监测时间和与监测时间对应的气体浓度位置输入大气扩散模型，得到当前周期的大气扩散模型计算浓度序列。

[0155] 在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：根据当前周期的气体在水平方向、竖直方向的扩散系数以及与监测时间对应的风速信息，得到与当前周期的大气扩散模型相关的模型参数。

[0156] 在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列；其中，第一目标函数适应度序列同与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置存在第一对应关系；对第一目标函数适应度序列进行排序及筛选，得到第二目标函数适应度序列；对第二目标函数适应度序列执行预设次数遗传算法，得到中间种群适应度序列；其中，中间种群适应度序列与第二目标函数适应度序列存在第二对应关系；对中间种群适应度序列进行局部寻优，得到最优中间种群适应度；根据第一对应关系、第二对应关系、最优中间种群适应度和目标函数模型，得到当前周期的目标函数模型计算值。

[0157] 在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：根据与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置，确定初始泄漏源强度集和初始泄漏源位置集；根据初始泄漏源强度集和大气扩散模型计算浓度序列，得到初始大气扩散模型计算浓度集；根据初始泄漏源强度集、初始泄漏源位置集、当前周期的气体浓度信息和初始大气扩散模型计算浓度集，得到第一目标函数适应度序列。

[0158] 在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置；确定第一泄漏源强度和第二泄漏源强度之间的强度偏差值，确定第一泄漏源位置和第二泄漏源位置之间的距离偏差值；若强度偏差值与距离偏差值的加权和小于或等于第一预设阈值，则将第二泄漏源强度确定为泄漏源强度，将第二泄漏源位置确定为泄漏源位置。

[0159] 在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：若强度偏差值与距离偏差值的加权和大于第一预设阈值，则将下一周期作为更新后的当前周期，并返回执行监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值的步骤。在其中一个实施例中，处理器执行计算机程序时还实现以下步骤：若当前周期的目标函数模型计算值小于或等于第二预设阈值，则将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置；若下一周期的目标函数模型计算值小于或等于第二预设阈值，则将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置。

[0160] 在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：

[0161] 根据预设的大气扩散模型构建与气体浓度信息关联的目标函数模型；

[0162] 监测当前周期的气体浓度信息，根据当前周期的气体浓度信息、大气扩散模型和目标函数模型，得到当前周期的目标函数模型计算值；

[0163] 监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值；

[0164] 根据当前周期的目标函数模型计算值和下一周期的目标函数模型计算值，确定泄漏源强度和泄漏源位置。

[0165] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：监测当前周期的气体浓度信息，根据当前周期的气体浓度信息和大气扩散模型，得到当前周期的大气扩散模型计算浓度序列；根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到当前周期的目标函数模型计算值。

[0166] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：将获取到的与大气扩散模型相关的模型参数、监测时间和与监测时间对应的气体浓度位置输入大气扩散模型，得到当前周期的大气扩散模型计算浓度序列。

[0167] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据当前周期的气体在水平方向、竖直方向的扩散系数以及与监测时间对应的风速信息，得到与当前周期的大气扩散模型相关的模型参数。

[0168] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据目标函数模型、当前周期的气体浓度信息和当前周期的大气扩散模型计算浓度序列，得到第一目标函数适应度序列；其中，第一目标函数适应度序列同与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置存在第一对应关系；对第一目标函数适应度序列进行排序及筛选，得到第二目标函数适应度序列；对第二目标函数适应度序列执行预设次数遗传算法，得到中间种群适应度序列；其中，中间种群适应度序列与第二目标函数适应度序列存在第二对应关系；对中间种群适应度序列进行局部寻优，得到最优中间种群适应度；根据第一对应关系、第二对应关系、最优中间种群适应度和目标函数模型，得到当前周期的目标函数模型计算值。

[0169] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据与监测时间对应的气体浓度大小和与监测时间对应的气体浓度位置，确定初始泄漏源强度集和初始泄漏源位置集；根据初始泄漏源强度集和大气扩散模型计算浓度序列，得到初始大气扩散模型计算浓度集；根据初始泄漏源强度集、初始泄漏源位置集、当前周期的气体浓度信息和初始大气扩散模型计算浓度集，得到第一目标函数适应度序列。

[0170] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置，以及将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置；确定第一泄漏源强度和第二泄漏源强度之间的强度偏差值，确定第一泄漏源位置和第二泄漏源位置之间的距离偏差值；若强度偏差值与距离偏差值的加权和小于或等于第一预设阈值，则将第二泄漏源强度确定为泄漏源强度，将第二泄漏源位置确定为泄漏源位置。

[0171] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：若强度偏差值与距离偏差值的加权和大于第一预设阈值，则将下一周期作为更新后的当前周期，并返回执行监测下一周期的气体浓度信息，根据下一周期的气体浓度信息、大气扩散模型和目标函数模型，得到下一周期的目标函数模型计算值的步骤。

[0172] 在其中一个实施例中，计算机程序被处理器执行时还实现以下步骤：若当前周期的目标函数模型计算值小于或等于第二预设阈值，则将与当前周期的目标函数模型计算值对应的气体浓度大小确定为第一泄漏源强度，将与当前周期的目标函数模型计算值对应的气体浓度位置确定为第一泄漏源位置；若下一周期的目标函数模型计算值小于或等于第二预设阈值，则将与下一周期的目标函数模型计算值对应的气体浓度大小确定为第二泄漏源强度，将与下一周期的目标函数模型计算值对应的气体浓度位置确定为第二泄漏源位置。

[0173] 为实现危化品泄漏源的定位，需要采集足够的监测数据，利用基于监测数据与危化品大气扩散模型构成的目标函数模型计算出泄漏源位置。在采集危化品监测信息方面，本申请利用离子迁移谱仪作为监测手段，能对极低浓度的危化品气体进行移动监测，利用离子迁移谱测得的当前周期的气体浓度信息作为目标函数模型的输入数据。在基于定位周期的泄漏源定位方法方面，以某一时刻在某一位置处离子迁移谱仪测得的危化品浓度与危化品大气扩散模型计算浓度的匹配度作为目标函数，建立基于离子迁移谱仪监测的泄漏源定位模型，利用基于遗传算法和模式搜索混合优化算法不断修正定位模型的参数，得到的最优泄漏源参数即为所估计的泄漏源强度和泄漏源位置。

[0174] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。

[0175] 以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

[0176] 以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。