Calpuff模型对轻质气体扩散模拟的改进方法和装置与流程

本公开涉及本发明涉及轻质气体泄漏扩散的动态模拟技术领域，具体而言，涉及一种calpuff模型对轻质气体扩散模拟的改进方法和装置。

背景技术：

如今石油、化工企业在生产、储存、运输和使用过程中存在着各种危险物质。这些危险物质其中就包含轻质气体，由于排放过后比空气密度小，在浮力的作用下，相比普通的气体会有明显的抬升；通过来自边界层的湍流反射作用也会对轻质气体的输送和扩散作用产生一定影响。详细过程如下：当轻气发生连续泄漏时，它所形成的羽流由于与空气之间存在密度差，在浮力的作用下羽流会有明显上升的现象；而且当泄漏轻气的初始温度高于空气时，这一羽流上升作用将更为明显。在主导风向的作用下，上升的羽流向下风向方向弯曲，形成一个典型的弯曲羽流。此时需要对羽流的有效泄漏源高进行计算；来自于地表面和边界层的瑞流反射作用显然也对轻气的输送和扩散作用产生一定的影响。因此地表状况的不同与所需要计算的轻气浓度是息息相关的，如高密度城市区域和旷野地区的扩散，最后影响的区域范围显然是不同的。

calpuff模型是一个用来模拟不稳定状态的多层、多物种污染的高斯型烟团扩散模式，它适用于模拟时空都在变化的气象条件下污染物的迁移、转化和清除。它考虑了复杂地形的影响，海岸的交界影响，建筑物的下洗影响，干湿沉降以及简单的化学转化，可以计算出在预设点的浓度和沉降量。但是，针对目前的calpuff模型具有一定的限制性。具体的，针对轻质气体排放后的抬升：calpuff模型在计算气态污染物排放高度的时候，会考虑污染物由于温度小于空气温度造成一定程度上的抬升，但是对于轻质气体，不仅温度其密度与空气存在密度差，在浮力作用下会有明显上升，calpuff模型无法模拟这种特性；针对轻质气体的扩散阶段：calpuff模型在计算完气体的抬升之后，将污染物气团密度看作与空气密度相当，污染物气团随大气迁移扩散，但是对于轻质气体，更易受到来自于地表面和边界层的湍流反射作用，导致轻质气体经过边界层的反射后有一定程度的下沉速度，calpuff模型无法模拟这种特性。

技术实现要素：

为了解决现有技术中的技术问题，本公开针对不同特征的轻质气体，对calpuff模型的参数进行调整，使得calpuff能够实现模拟轻质气体迁移扩散。即本公开实施例提供了一种calpuff模型对轻质气体扩散模拟的改进方法和装置，获取事故源信息与环境信息，并对事故源信息与环境信息进行收集；对收集过程中的排放高度参数与沉降速度参数进行优化操作；将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。经过排放高度和沉降速度两个参数优化后的calpuff模型，能够适用于模拟轻质气体的迁移扩散，从而能够更方便快速的对轻质气体的泄漏事故进行动态模拟；提出的参数优化的公式里面的参数均比较容易获取，从而能够在泄漏事故中以较短的时间内搜集到必须的数据，以实现泄漏事故的动态模拟。

第一方面，本公开实施例提供了一种calpuff模型对轻质气体扩散模拟的改进方法，包括以下步骤：获取事故源信息与环境信息，并对所述事故源信息与所述环境信息进行收集；对收集过程中的排放高度参数与沉降速度参数进行优化操作；将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。

在其中一个实施例中，所述事故源信息包括：轻质气体的排放高度、轻质气体的排放温度、轻质气体的排放速度以及轻质气体的排放量。

在其中一个实施例中，所述环境信息包括：环境风速、莫氏长度、摩擦速度、季节/下垫面。

在其中一个实施例中，所述对收集与预处理过程中的排放高度参数与沉降速度参数进行优化操作包括：针对轻质气体计算浮力通量以及针对轻质气体计算上升的抬升高度。

在其中一个实施例中，所述对收集与预处理过程中的排放高度参数与沉降速度参数进行优化操作还包括：

根据轻质气体的泄漏源参数，通过第一预设公式计算上升的抬升高度，其中，所述第一预设公式为：

其中，x为下风向距离(m)，u为环境风速(m/s)，轻质气体额外的抬升高度随下风向距离增大而增大，δh为由于轻质气体造成的上升的抬升高度；

通过第二预设公式计算calpuff模型中的排放高度的输入值，其中，所述第二预设公式为：hi＝hs+δh；

其中，hs为实际排放高度(m)，hi为calpuff模型中输入的排放高度(m)；

根据第三预设公式计算浮力通量，其中，所述第三预设公式为：其中，vs为排放出口处烟气排放速率(m/s)，g为重力加速度(m/s²)，ts为烟气出口温度(k)，ta为环境大气温度(k)。

在其中一个实施例中，所述对收集与预处理过程中的排放高度参数与沉降速度参数进行优化操作还包括：

根据所述环境信息，基于第四预设公式计算气体的沉降速度，计算的沉降速度vd作为calpuff模型中的气体干沉降速度的输入值，所述第四预设公式为：

其中，ra——湍流层空气动力学阻力，rb——准层流副层阻力；ra是由湍流运动引起的污染物从大气输送至地面附近时受到的阻力，根据近地层由污染物的质量输送和热量输送类似的假定获取第五预设公式，其中，所述第五预设公式为：

其中，zr是计算vd时选择的参考高度(m)，ua是摩擦速度(m/s)，k是vonkarman常数，取0.4，z0是下垫面的粗糙度长，由局地地形特征经验地获得(m)，ψc为与质量通量有关的稳定度修正函数，

其中，1/l是monin-obukhov长度的倒数(1/m)；

rb是污染物向地表沉积时经过近地面片流层所受的阻力，其和通过与地表相接的类片流层的传输相关联；

其中，ua是摩擦速度(m/s)，k是vonkarman常数，取0.4，pr为prandtl数，取值0.72，sc为schmidt数，sc＝μ0/di，μ0为空气运动粘度，取值0.15cm²s^-1，di为与组分相关的第1种污染物的分子扩散度。

在其中一个实施例中，还包括：通过判断多个季节多个下垫面类型获取受体表面阻力。

第二方面，本公开实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述的方法的步骤。

第三方面，本公开实施例提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述的方法的步骤。

第四方面，本公开实施例提供了一种calpuff模型对轻质气体扩散模拟的改进装置，所述装置包括：获取与收集模块，用于获取事故源信息与环境信息，并对所述事故源信息与所述环境信息进行收集；优化模块，用于对收集过程中的排放高度参数与沉降速度参数进行优化操作；结果获取模块，用于将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。

本发明提供的一种calpuff模型对轻质气体扩散模拟的改进方法和装置，获取事故源信息与环境信息，并对事故源信息与环境信息进行收集；对收集过程中的排放高度参数与沉降速度参数进行优化操作；将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。该方法经过排放高度和沉降速度两个参数优化后的calpuff模型，能够适用于模拟轻质气体的迁移扩散，从而能够更方便快速的对轻质气体的泄漏事故进行动态模拟；提出的参数优化的公式里面的参数均比较容易获取，从而能够在泄漏事故中以较短的时间内搜集到必须的数据，以实现泄漏事故的动态模拟。

附图说明

为了更清楚地说明本公开实施例的技术方案，下面对实施例描述中所需要使用的附图作简单地介绍：

图1为本发明一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法的步骤流程示意图；

图2为本发明另一实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法的步骤流程示意图；

图3为本发明一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法中的下风向中心线上最大时均浓度变化示意图；

图4为本发明一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法中的下风向地面时均浓度场示意图；

图5为本发明一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进装置的结构示意图。

具体实施方式

下面结合附图和实施例对本申请进行进一步的详细介绍。

在下述介绍中，术语“第一”、“第二”仅为用于描述的目的，而不能理解为指示或暗示相对重要性。下述介绍提供了本公开的多个实施例，不同实施例之间可以替换或者合并组合，因此本申请也可认为包含所记载的相同和/或不同实施例的所有可能组合。因而，如果一个实施例包含特征a、b、c，另一个实施例包含特征b、d，那么本申请也应视为包括含有a、b、c、d的一个或多个所有其他可能的组合的实施例，尽管该实施例可能并未在以下内容中有明确的文字记载。

为了使本发明的目的、技术方案及优点更加清楚明白，以下通过实施例，并结合附图，对本发明一种calpuff模型对轻质气体扩散模拟的改进方法和装置的具体实施方式进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

如图1所示，为一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法的流程示意图，具体包括以下步骤：

步骤102，获取事故源信息与环境信息，并对事故源信息与环境信息进行收集。具体的，事故源信息包括：轻质气体的排放高度、轻质气体的排放温度、轻质气体的排放速度以及轻质气体的排放量。由此，提高了针对轻质气体的事故源信息获取的多样性与选择性。此外，环境信息包括：环境风速、莫氏长度、摩擦速度、季节/下垫面。由此，提高了环境信息的多选择性。

步骤104，对收集过程中的排放高度参数与沉降速度参数进行优化操作。需要说明的是，排放高度参数优化中包括计算浮力通量以及计算上升的抬升高度；沉降速度参数优化包括计算湍流运动阻力、计算近地面阻力以及计算受体表面阻力。具体的，对收集与预处理过程中的排放高度参数与沉降速度参数进行优化操作包括：针对轻质气体计算浮力通量以及针对轻质气体计算上升的抬升高度。

进一步地，对收集与预处理过程中的排放高度参数与沉降速度参数进行优化操作还包括：

根据轻质气体的泄漏源参数，通过第一预设公式计算上升的抬升高度，其中，第一预设公式为：

其中，x为下风向距离(m)，u为环境风速(m/s)，轻质气体额外的抬升高度随下风向距离增大而增大，δh为由于轻质气体造成的上升的抬升高度；

通过第二预设公式计算calpuff模型中的排放高度的输入值，其中，第二预设公式为：hi＝hs+δh；

其中，hs为实际排放高度(m)，hi为calpuff模型中输入的排放高度(m)；

根据第三预设公式计算浮力通量，其中，所述第三预设公式为：其中，vs为排放出口处烟气排放速率(m/s)，g为重力加速度(m/s²)，ts为烟气出口温度(k)，ta为环境大气温度(k)。

更进一步地，对收集与预处理过程中的排放高度参数与沉降速度参数进行优化操作还包括：

根据环境信息，基于第四预设公式计算气体的沉降速度，计算的沉降速度vd作为calpuff模型中的气体干沉降速度的输入值，第四预设公式为：

其中，ra——湍流层空气动力学阻力，rb——准层流副层阻力；ra是由湍流运动引起的污染物从大气输送至地面附近时受到的阻力，根据近地层由污染物的质量输送和热量输送类似的假定获取第五预设公式，其中，所述第五预设公式为：

其中，zr是计算vd时选择的参考高度(m)，ua是摩擦速度(m/s)，k是vonkarman常数，取0.4，z0是下垫面的粗糙度长，由局地地形特征经验地获得(m)，ψc为与质量通量有关的稳定度修正函数，

其中，1/l是monin-obukhov长度的倒数(1/m)；

rb是污染物向地表沉积时经过近地面片流层所受的阻力，其和通过与地表相接的类片流层的传输相关联；

其中，ua是摩擦速度(m/s)，k是vonkarman常数，取0.4，pr为prandtl数，取值0.72，sc为schmidt数，sc＝μ0/di，μ0为空气运动粘度，取值0.15cm²s^-1，di为与组分相关的第1种污染物的分子扩散度。

此外，还需要说明的是，本公开所涉及的calpuff模型对轻质气体扩散模拟的改进方法还包括：通过判断多个季节多个下垫面类型获取受体表面阻力。

步骤106，将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。

本发明提供的一种calpuff模型对轻质气体扩散模拟的改进方法，获取事故源信息与环境信息，并对事故源信息与环境信息进行收集；对收集过程中的排放高度参数与沉降速度参数进行优化操作；将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。该方法经过排放高度和沉降速度两个参数优化后的calpuff模型，能够适用于模拟轻质气体的迁移扩散，从而能够更方便快速的对轻质气体的泄漏事故进行动态模拟；提出的参数优化的公式里面的参数均比较容易获取，从而能够在泄漏事故中以较短的时间内搜集到必须的数据，以实现泄漏事故的动态模拟。

为了更精准、灵活地理解与运用本公开所提出的一种calpuff模型对轻质气体扩散模拟的改进方法，进行以下示例。需要说明的是，本公开所保护的范围不限于以下示例。

需要说明的是，通常气体的schmidt和prandtl数与分子量相关，具体参考表1，即表1为几种气体的schmidt和prandtl数及相对于水分子扩散系数比率。

此外，还需要说明的是，rc是受体表面阻力，包含其他多种阻力，通常通过判断不同季节不同下垫面类型获得。如表2所示意。表2各种下垫面沉积时所受到的表面阻力。

具体的，结合图2-4所示，即图2为本发明另一实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法的步骤流程示意图；图3为本发明一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法中的下风向中心线上最大时均浓度变化示意图；图4为本发明一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进方法中的下风向地面时均浓度场示意图。针对文献中有毒氨气的试验，有毒氨气在白天因意外事故发生连续泄漏，源强为20kg/s，泄漏初始速率为1.5m/s，实际泄漏高度为3m，定义氨气与大气环境温差因子为0.1，事故发生时风速为4m/s，大气稳定度为c，地面粗糙度为0.5m。

首先，关于排放高度计算，根据现场收集到的参数，首先计算浮力通量fb：

然后计算由于轻质气体导致的排放高度增加量δh。

δh＝1.6×1.47^1/3×4^-1×(3×10)^2/3＝4.4；则经过参数优化的排放高度hi＝hs+δh＝3+4.4＝7.4。

进一步地，关于沉降速度计算，根据现场收集到的参数，首先计算由湍流运动引起的污染物从大气输送至地面附近时受到的阻力ra：

然后计算污染物向地表沉积时经过近地面片流层所受的阻力rb。即然后根据现场季节及下垫面类型判断受体表面阻力rc；rc＝0；最后计算轻质气体干沉降速度vd；

综上所述，calpuff计算结果及与实测结果对比，将上述计算的排放高度和沉降速度作为calpuff模型参数输入，最终得到本次泄漏事故液氯的浓度场。根据国家卫生部门相关规定，当氨气浓度达到30mg/m³时，人体出现轻微反应；达到350mg/m³时，人体出现中度反应；达到550mg/m³时，人体出现强烈反应；达到1500mg/m³时，接触一段时间后可能导致死亡；当浓度大于3500mg/m³时，则可能吸入后立即死亡。根据这五个毒性阈值对模拟计算得到的氨气浓度影响区域进行划分，距泄漏源从近到远依次可以分为：立即致死区、死亡区、强烈反应区、中度反应区、以及轻微反应区。如表3所示，即表3模拟计算浓度与文献试验浓度对比。

从以上结果可以看出，除轻微反应区外，改进前calpuff估计的不同阈值范围的长度均大于文献试验值，而改进后该长度基本与文献试验值吻合，说明了进过改进的calpuff模型能够很好的模拟轻质气体。

综上可知，本发明公开了一种基于calpuff模型模拟轻质气体的改进方法，提出对calpuff模型的输入参数排放高度和干沉降速度进行参数优化，最后利用calpuff模型动态模拟轻质气体的迁移扩散。试验结果表明：该算法能很好的模拟出轻质气体迁移扩散特性。本公开通过一种算法，针对不同特征的轻质气体，对calpuff模型的参数进行调整，使得calpuff能够实现模拟轻质气体迁移扩散。对calpuff模型中排放高度和沉降速度两个参数进行优化，以实现calpuff模型对于轻质气体的模拟。

基于同一发明构思，还提供了一种calpuff模型对轻质气体扩散模拟的改进装置。由于此装置解决问题的原理与前述一种calpuff模型对轻质气体扩散模拟的改进方法相似，因此，该装置的实施可以按照前述方法的具体步骤实现，重复之处不再赘述。

如图5所示，为一个实施例中的一种calpuff模型对轻质气体扩散模拟的改进装置的结构示意图。该calpuff模型对轻质气体扩散模拟的改进装置10包括：获取与收集模块200、优化模块400和结果获取模块600。

其中，获取与收集模块200用于获取事故源信息与环境信息，并对事故源信息与环境信息进行收集；优化模块400用于对收集过程中的排放高度参数与沉降速度参数进行优化操作；结果获取模块600用于将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。

本发明提供的一种calpuff模型对轻质气体扩散模拟的改进装置，首先通过获取与收集模块获取事故源信息与环境信息，并对事故源信息与环境信息进行收集；再通过优化模块对收集过程中的排放高度参数与沉降速度参数进行优化操作；最终通过结果获取模块将优化后的参数输入至calpuff模型中，获取轻质气体浓度场结果。该装置经过排放高度和沉降速度两个参数优化后的calpuff模型，能够适用于模拟轻质气体的迁移扩散，从而能够更方便快速的对轻质气体的泄漏事故进行动态模拟；提出的参数优化的公式里面的参数均比较容易获取，从而能够在泄漏事故中以较短的时间内搜集到必须的数据，以实现泄漏事故的动态模拟。

本发明实施例还提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该程序被图1和图2中处理器执行。本发明实施例还提供了一种包含指令的计算机程序产品。当该计算机程序产品在计算机上运行时，使得计算机执行上述图1和图2的方法。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read-onlymemory，rom)或随机存储记忆体(randomaccessmemory，ram)等。以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。为了示例和描述的目的已经给出了以上描述。此描述不意图将本公开的实施例限制到在此公开的形式。尽管以上已经讨论了多个示例方面和实施例，但是本领域技术人员将认识到其某些变型、修改、改变、添加和子组合。