本发明属于应急救援辅助决策领域，涉及到自动化技术，特别是涉及一种化工事故多目标两阶段应急救援物资调度方法。
背景技术：
由于化工生产、储运等各环节多涉及危险化学品，以及化工园区多处于经济发达、人口密集的敏感地区，造成化工园区处于高风险状态。近年来，多起重大化工园区危险事故的发生，给人民的生命和财产带来了巨大的损害。当面对化工园区危险事故时，如何安全的规划救援车辆行驶路径以及合理的应急救援物资调度对救援工作的及时开展具有重要意义。
许多传统的应急救援物资调度方法在救援车辆的路径规划中只考虑路径长度因素，没有考虑到化工园区事故以及初始事故所引起的多米诺效应对园区以及周边道路产生的风险影响，因而容易导致规划的救援车辆路径具有较大风险。此外，因实际问题中，应急救援工作是系统而又复杂的科学问题，传统应急救援物资调度方法多以单个目标为救援调度目标，如救援时间最短、参与救援点最少，难以全面地考虑实际问题，可能造成救援工作受到影响。因此，应急救援物资调度非常困难。
技术实现要素：
本发明针对现有技术的不足，提出了一种化工事故多目标两阶段应急救援物资调度方法。
本发明的目标是针对应急救援物资调度中的一些难题，在满足所有给定的约束条件的背景下，分析化工园区事故后果以及多米诺效应，确定救援车辆路线，分配救援物资至可支配救援点，将救援物资从多个可支配救援点运输到受灾点，使救援开展时间、参与救援点数量和未满足救援物资数量最小化。
本发明的技术方案是通过充分考虑化工初始事故以及初始事故所引发的多米诺效应，化工园区道路信息、人口分布，提出一种基于道路风险的救援车辆路径规划模型，构建多目标两阶段应急救援物资调度模型。第一阶段根据路径特性，设计一种改进变异算子的遗传算法求解问题，第二阶段使用多目标遗传算法nsga-ii对模型进行求解，最终确立了化工事故多目标两阶段应急救援物资调度方法。
本发明的有益效果：本发明结合应急救援物资调度在传统的车辆路径规划方法的基础上，第一阶段考虑化工园区事故以及多米诺效应所产生的道路风险因素，构建了基于道路风险的救援车辆路径规划模型；第二阶段针对化工救援工作的复杂性，构建多目标应急救援物资调度模型，最终设计了一个化工事故多目标两阶段应急救援物资调度方法，本发明具有开放性、灵活性以及计算复杂度低等特点。
附图说明
图1为多米诺效应流程图；
图2为遗传算法流程图；
图3为交叉算子示意图；
图4为变异算子示意图；
图5为多目标遗传算法nsga-ii流程图。
具体实施方式
以下结合附图对本发明作进一步说明
本发明包括以下步骤：
步骤1:获取基础数据，其中包括化工园区事故信息、化工园区路段信息、周围可支配救援点物资储备信息、人口分布信息。
步骤2:分析化工园区以及周边道路风险，本发明将道路风险定义为化工园区初始事故以及多米诺效应下所产生个人风险的叠加。通过对初始事故类型以及造成伤害的主要因素进行分析，评价主要损伤因素对人体伤害以及周围设备所在造成的影响。其次分析初始事故所可能引发多米诺效应，从确定初始事故出发，针对热辐射，冲击波超压，抛射碎片等物理效应计算邻近目标设备的损坏几率，确定多米诺效应事故发生概率以及后果，确定失效单元，判断是否有下级单元。多米诺效应见图1，其中设备的损坏几率计算公式如下：
其中ym是不同损伤因素下的设备损伤概率单位，其计算公式如下所示：
ln1(ttf)＝-1.128ln(i)-2.667×10-5v 9.87
ln2(ttf)＝-0.97ln(i) 8.835v0.032
式中，tff为设备失效所需时间，s；v为二次设备容积，m3；ps为静态超压峰值，kpa；i为热辐射的强度，kw/m2。
通过上述两个部分可以计算初始事故所产生的后果，以及多米诺效应下事故产生的后果。事故后果带来的致死概率v公式如下所示:
其中yh为概率计算单位，其计算公式如下所示：
式中，i为热辐射强度，kw/m2；ps为静态超压峰值，kpa；c为毒气浓度，ppm；te为人体暴露时间，min；d为剂量。
则区域内任一点个人风险值计算表达式如下所示:
其中fa是第a个设备发生事故的频率，va是第a个设备发生事故造成的人员死亡概率,pab为设备a所引发的事故对设备b的多米诺效应拓展概率。
步骤3：构建化工事故多目标两阶段应急救援物资调度模型。
第一阶段：构建基于道路风险的救援车辆路径规划模型
本发明构建设化工园区路径网络节点无向图g＝(v,e),其中v为结点集,e为路径集，e中任意路段的长度与风险数值分别设为len、risk。因为len、risk之间量级差别较大，所以对道路长度和风险权重进行了标准归一化处理，使得两个权重对路径的选择具有相同的影响。路径规划模型如下式所示。
minz＝min(α1∑len α2∑risk)
其中α1和α2为权重系数，分别表示距离和安全两个不同的优化指标的重要性。该模型所计算出的最优路径为第二阶段多目标应急救援物资调度模型救援车辆运行时间参数提供根据。
第二阶段：构建多目标应急救援物资调度模型
本发明定义参数与变量分别为：j：表示事故灾害点的第j种应急资源需求，j＝1,2,...,m；dj：表示事故灾害点对第j种应急资源的需求量；vj：表示第j种应急资源的消耗速率；ai：表示第i个备选应急出救点，i＝1,2,...,n；ti：表示第i个备选应急出救点到事故灾害点的应急资源调配行程时间；rij：表示应急出救点i的第j种应急资源保有量；
为便于构建多目标应急资源调度模型，作基本假设：各应急救援点之间相互独立，不存在相互资源调度与补给；运输时间与车辆类型无关；各应急物资消耗速率为常量；救援过程中不考虑救援车辆装货与卸货时间；各救援点具有充足车辆用于救援物资调度；
本发明将优化目标函数定义为一个多目标的三维向量，minimize(z)＝[z1,z2,z3]，其中，z1是用于最小救援开展时间的目标函数；z2是用于最小化未满足救援物资数量的目标函数；z3是用于最小化救援参与点数量的目标函数。
目标函数z1：
上式表示为救援时效性最好
目标函数z2：
上式表示为救援安全性最好
目标函数z3：
上式表示为代表救援经济性最好
模型约束为：
上式表示为各应急救援点的资源派遣量不大于其资源保有量，且为非负整数；
上式中，ε为充分小数，此式确保当应急出救点ai派遣的第j种应急资源量
上式表示为当应急出救点ai参与应急资源调度时yi取1，否则取0
步骤4:模型求解
第一阶段模型求解：
(1)基于道路风险的车辆路径规划模型，应用了改进变异算子的的遗传算法进行求解，算法流程见图2所示，具体求解步骤如下：
①参数设置
设置种群数目n，最大迭代次数t,变异算子为0.05，交叉算子0.7；
②编码
一条路径为一个染色体，染色体编码没有采用传统遗传算法中的数值编码方式(如二进制和实数编码)而是采用该路径的结点序列,这样能更直接、具体地表达路径的实际意义,同时有利于适应度函数的选取与计算。结点号为染色体中的基因。由于路径不应出现环路,因此染色体中不允许有重复的基因码。染色体的第一个基因为路径起点,最后一个基因为路径终点,染色体的长度是变化的,但不可能超过最大长度n(n为无向图的结点数)。
③初始化种群
随机生成以救援点为起点，受灾点为终点的n条可行路径。
④适应度评价
根据染色体对应的路径风险数值以及长度数值对路径进行评价。
⑤选择
按照轮盘赌规则选择n个个体。
⑥交叉
交叉两个个体部分基因可以发现更好的个体。传统遗传算法中是随机进行一点或多点交叉的操作,但这样很容易产生断路或环路。因此本文在除首、尾结点之外的第一个重复结点进行交叉操作且只进行一点交叉操作。交叉算子示意图见图3。
⑦变异
本发明中对传统遗传算法变异进行改进，传统遗传算法在解决路径问题时，变异操作是随机选择路径中两个结点，并随机生成两结点之间可行路径对原路径进行替换。现对变异算子做出改进，随机选择路径某个结点，使用迪杰斯特拉算法生成路径起点到当前结点、当前结点到路径终点的最短路径，对原路径进行替换。变异算子示意图见图4。
⑧算法迭代
重复④到⑦，直至满足算法最大迭代数。最终，生成满足条件的路径。
第二阶段模型求解：
本发明应用了多目标遗传算法nsga-ii进行求解，算法流程图见图5，求解的具体步骤如下:
①初始化种群
设置初始种群数目大小为n，最大迭代次数t,变异算子为0.05，交叉算子为0.7，初始化种群，随机生成满足限制条件的n个个体；
②帕累托排序
在完成初始种群构建后，采用帕累托排序法划分种群层次。帕累托排序法采用帕累托支配关系比较种群中个体之间的优劣。假设种群为p，帕累托排序法需要计算种群内每一个个体p∈p的两个参数np和sp，np是种群中支配个体p的个体数，sp是种群中个体p支配的个体集合。当遍历种群内所有个体后，所有np＝0的个体将被划分到种群的第一层p1，对于p1内的个体l∈p1，其所支配的个体集合为sl，遍历sl中的个体m∈sl，执行nm＝nm-1，所有nm＝0的个体将被划分到种群的第二层p2，以此类推，直到整个种群被分层。序号为1的分层p1为非支配层，p1内所有个体均为当前种群的帕累托最优解。如图2所示，赋秩1的解优先于赋秩2的解，即，赋值秩为1的解在任何优化目标上均优于其余解。因此，分配给第1级的解是非支配的。
③进化算子
通过二进制锦标算法在当前种群内选择个体进行交叉变异操作，产生n个子代个体，分层序号越小的个体，被选中进行交叉变异操作的概率越大。其中交叉操作是通过随机选择矩阵的一行和一列，并随机交换由行和列组成的交叉区域之一，变异操作是通过随机选择矩阵的一个元素，设该点所对应的某种救援物资存储量为storage，并让其在[0,storage]间进行随机的取值。对父代种群和子代种群执行精英策略得到n个个体作为下一代的父代种群；
④算法迭代
重复②到③，直至满足算法最大迭代数。最终，基于变邻域搜索的混合多目标进化算法给出一组满足不同偏好的帕累托最优解。
1.一种化工事故多目标两阶段应急救援物资调度方法，其特征在于该方法包括以下步骤：
步骤1：获取基础数据，包括化工园区事故信息、化工园区路段信息、周围可支配救援点物资储备信息和人口分布信息；
步骤2：分析化工园区以及周边道路风险，将道路风险定义为化工园区初始事故以及多米诺效应下所产生个人风险的叠加；通过对初始事故类型以及造成伤害的主要因素进行分析，评价主要损伤因素对人体伤害以及周围设备所在造成的影响；
其次分析初始事故所可能引发多米诺效应，从确定初始事故出发，针对热辐射，冲击波超压，抛射碎片计算邻近目标设备的损坏几率，确定多米诺效应事故发生概率以及后果，确定失效单元，判断是否有下级单元；
计算初始事故所产生的后果，以及多米诺效应下事故产生的后果；事故后果带来的致死概率v；
则区域内任一点个人风险值计算表达式如下所示：
其中fa是第a个设备发生事故的频率，va是第a个设备发生事故造成的人员死亡概率，pab为设备a所引发的事故对设备b的多米诺效应拓展概率；
步骤3：构建化工事故多目标两阶段应急救援物资调度模型；
第一阶段：构建基于道路风险的救援车辆路径规划模型
构建设化工园区路径网络节点无向图g＝(v，e)，其中v为结点集，e为路径集，e中任意路段的长度与风险数值分别设为len、risk，路径规划模型如下式所示；
minz＝min(α1∑len α2∑risk)
其中α1和α2为权重系数，分别表示距离和安全两个不同的优化指标的重要性；
该模型所计算出的最优路径为第二阶段多目标应急救援物资调度模型救援车辆运行时间参数提供根据；
第二阶段：构建多目标应急救援物资调度模型
定义参数与变量分别为：j：表示事故灾害点的第j种应急资源需求，j＝1，2，...，m；dj：表示事故灾害点对第j种应急资源的需求量；vj：表示第j种应急资源的消耗速率；ai：表示第i个备选应急出救点，i＝1，2，...，n；ti：表示第i个备选应急出救点到事故灾害点的应急资源调配行程时间；rij：表示应急出救点i的第j种应急资源保有量；
设各应急救援点之间相互独立，不存在相互资源调度与补给；运输时间与车辆类型无关；各应急物资消耗速率为常量；救援过程中不考虑救援车辆装货与卸货时间；各救援点具有充足车辆用于救援物资调度；
将优化目标函数定义为一个多目标的三维向量，minimize(z)＝[z1，z2，z3]，其中，z1是用于最小救援开展时间的目标函数；z2是用于最小化未满足救援物资数量的目标函数；z3是用于最小化救援参与点数量的目标函数；
目标函数z1：
上式表示为救援时效性最好
目标函数z2：
上式表示为救援安全性最好
目标函数z3：
上式表示为代表救援经济性最好
模型约束为：
上式表示为各应急救援点的资源派遣量不大于其资源保有量，且为非负整数；
上式中，ε为充分小数，此式确保当应急出救点ai派遣的第j种应急资源量
上式表示为当应急出救点ai参与应急资源调度时yi取1，否则取0；
步骤4：模型求解
第一阶段模型求解：
(1)基于道路风险的车辆路径规划模型，应用改进变异算子的遗传算法进行求解，具体求解步骤如下：
①参数设置
设置种群数目，最大迭代次数，变异算子和交叉算子；
②编码
一条路径为一个染色体，染色体编码采用路径的结点序列，结点号为染色体中的基因；染色体的第一个基因为路径起点，最后一个基因为路径终点；
③初始化种群
随机生成以救援点为起点，受灾点为终点的n条可行路径；
④适应度评价
根据染色体对应的路径风险数值以及长度数值对路径进行评价；
⑤选择
按照轮盘赌规则选择n个个体；
⑥交叉
在除首、尾结点之外的第一个重复结点进行交叉操作且只进行一点交叉操作；
⑦变异
随机选择路径某个结点，使用迪杰斯特拉算法生成路径起点到当前结点、当前结点到路径终点的最短路径，对原路径进行替换；
⑧迭代
重复④到⑦，直至满足最大迭代数；最终，生成满足条件的路径；
第二阶段模型求解：
应用多目标遗传算法nsga-ii进行求解，求解的具体步骤如下：
①初始化种群
设置初始种群数目，最大迭代次数，变异算子，交叉算子，初始化种群，随机生成满足限制条件的n个个体；
②帕累托排序
采用帕累托排序法划分种群层次；帕累托排序法采用帕累托支配关系比较种群中个体之间的优劣；
③进化
通过二进制锦标算法在当前种群内选择个体进行交叉变异操作，产生n个子代个体；
④迭代
重复②到③，直至满足最大迭代数；最终给出一组满足不同偏好的帕累托最优解。
2.根据权利要求1所述的一种化工事故多目标两阶段应急救援物资调度方法，其特征在于：
所述目标设备的损坏几率计算公式如下：
其中ym是不同损伤因素下的设备损伤概率单位，其计算公式如下所示：
ln1(ttf)＝-1.128ln(i)-2.667×10-5v 9.87
ln2(ttf)＝-0.97ln(i) 8.835v0.032
式中，tff为设备失效所需时间；v为二次设备容积；ps为静态超压峰值；i为热辐射的强度。
3.根据权利要求1所述的一种化工事故多目标两阶段应急救援物资调度方法，其特征在于：
所述致死概率v公式如下所示：
其中yh为概率计算单位，其计算公式如下所示：
式中，i为热辐射强度；ps为静态超压峰值；c为毒气浓度；te为人体暴露时间；d为剂量。
4.根据权利要求1所述的一种化工事故多目标两阶段应急救援物资调度方法，其特征在于：
在帕累托排序法中，假设种群为p，计算种群内每一个个体p∈p的两个参数np和sp，np是种群中支配个体p的个体数，sp是种群中个体p支配的个体集合；当遍历种群内所有个体后，所有np＝0的个体将被划分到种群的第一层p1，对于p1内的个体l∈p1，其所支配的个体集合为sl，遍历sl中的个体m∈sl，执行nm＝nm-1，所有nm＝0的个体将被划分到种群的第二层p2，以此类推，直到整个种群被分层；序号为1的分层p1为非支配层，p1内所有个体均为当前种群的帕累托最优解；赋秩1的解优先于赋秩2的解，即，赋值秩为1的解在任何优化目标上均优于其余解，因此，分配给第1级的解是非支配的。
5.根据权利要求1所述的一种化工事故多目标两阶段应急救援物资调度方法，其特征在于：
所述的变异操作是通过随机选择矩阵的一个元素，设该点所对应的某种救援物资存储量为storage，并让其在[0，storage]间进行随机的取值；对父代种群和子代种群执行精英策略得到n个个体作为下一代的父代种群。
技术总结