一种燃气管线泄漏突发事件安全疏散方法及系统

1.本发明涉及突发事件安全疏散领域，更具体涉及一种燃气管线泄漏突发事件安全疏散方法及系统。


背景技术：

2.近年来，燃气泄漏突发事件频繁发生，灾害放大效应也越加明显。燃气突发事件具有连锁性、复杂性和放大性的特点，容易造成严重的人员伤亡后果。因此，当燃气泄漏事件发生时，需要迅速做出决策，及时、正确地进行人员疏散，才有可能将伤害降至最低。因此，有必要研究一种燃气管线泄漏突发事件的安全疏散方法。
3.1.中国专利申请号cn201810348059.8，公开了一种危化企业应急疏散线路优选方法及装置，该方法是通过构建真实泄漏场景集合进行泄漏模拟，实时记录各泄漏场景中各预设监测点气体浓度，确定可能泄漏扩散区域；判断不同泄漏场景中各待疏散点的气体浓度是否达到预设报警阈值，将可触发报警的泄漏场景作为该待疏散点的疏散场景集合，确定疏散路线，并从中选择到达所有可能泄漏扩散区域累加边界上任意路线出口点的疏散路线为备选疏散路线集合；计算备选疏散路线各疏散场景人员中毒风险得到累积个人中毒风险，遍历备选疏散路线集合，以最小累积个人中毒风险为优化函数得到最优疏散路线。该方法的缺点是没有考虑有毒气体浓度不同伤害程度对于人员疏散路径选择的影响，这可能会增加人员受伤害的风险。
4.2.中国专利申请号cn201610320272.9，公开了毒气泄漏事故下的人员疏散安全评估方法、装置及设备，该方法是采用泄漏扩散模型进行模拟，得出有毒气体时间序列的空间浓度分布状况，即浓度场。模拟沿设定疏散路径的人员疏散过程，获取疏散个体在不同时间的空间位置，计算设定疏散路径的疏散个体的毒性负荷，根据毒性负荷与有毒气体的种类计算疏散个体的致死概率，进而得到疏散个体的死亡率，最终判断在设定疏散路径的人员疏散安全程度。该方法通过评估在疏散路径上人员受有毒有害气体的影响，进而确定疏散路径的安全性，但考虑的因素比较片面，路径选择方法比较单一，很多场景下可能并不适用。
5.现有关于突发事件安全疏散的方法在研究思路上都比较一致，通常都是先确定疏散范围，再确定疏散范围内的人员数量、道路分布与庇护场所位置，最后研究最佳疏散路径。但是这些方法首先在考虑突发事件人员疏散范围时，对于事件后果类型的分析都不够充分；其次在确定最佳疏散路径时，选取的影响因素也多是道路的拥挤度、人员疏散速度等，而没有考虑事件后果不同伤害程度对于疏散路径选择的影响；最后在应对突发事件选择庇护场所方面，没有考虑到实际上庇护场所距离事发地的距离比较远，要经历多个疏散路段才能完成疏散，这就有可能会造成虽然选择的疏散路径比较短，但由于疏散人群在该路径中处于疏散范围内的时间较长，而大大增加了人员受伤害的风险。综上所述现有技术关于突发事件安全疏散的方法考虑不全面，难以实现安全高效的人员疏散。


技术实现要素：

6.本发明所要解决的技术问题在于现有技术关于突发事件安全疏散的方法存在考虑不全面，难以实现安全高效的人员疏散的问题。
7.本发明通过以下技术手段实现解决上述技术问题的：一种燃气管线泄漏突发事件安全疏散方法，所述方法包括：
8.步骤一：根据燃气泄漏事故后果模型获取燃气泄漏突发事件的各类后果影响范围以及伤害分区，根据各类后果影响范围以及伤害分区获取燃气泄漏突发事件的疏散范围；
9.步骤二：获取疏散出发点、疏散出口节点、警戒点、疏散路网节点以及庇护场所；
10.步骤三：获取疏散路网内各路段的疏散时间，根据各路段的疏散时间、通行道路所处的事故伤害分区及周边危化企业数量等信息获取疏散效果；
11.步骤四：获取各疏散出发点经过中间节点到对应的疏散出口的所有路径，根据各路径的疏散效果确定最佳疏散效果路径；
12.步骤五：获取合适的庇护场所。
13.本发明综合考虑燃气突发事件的各类后果影响范围及伤害分区，将伤害分区叠加到疏散路径上，对各个疏散路径上的风险与时间进行综合性评估，进而得到最佳疏散效果的路径及疏散出口，并给出适合的庇护场所位置，从而实现高效、安全的人员疏散。
14.优选的，所述步骤二包括：
15.步骤201：根据地理信息系统识别出疏散范围内的各类建筑物位置信息，可由城市人口密度热力图或相关政府部门获取各类建筑物的人口数量，若无法获取，则根据建筑物的属性预估该影响范围内的人口数量；
16.步骤202：以各类建筑物的出口作为疏散范围内人群疏散出发点a
o
且
17.步骤203：识别出疏散范围内的所有道路交通控制点即交叉路口作为疏散路径的中间节点a
m
且识别出疏散范围外所有最近的道路交叉口作为疏散出口节点a
e
且中间节点a
m
、疏散出口节点a
e
即为燃气管线泄漏突发事件疏散需要设置警戒的警戒点；
18.步骤204：由疏散的出发点、中间节点、疏散出口节点共同构成了燃气泄漏突发事件的疏散路网节点a
ω
且a
ω
＝{a
o
,a
m
,a
e
}，则用表示a
ω
内任意两个相邻节点之间的道路距离，其中a，p＝1,2,3
‥
，b+c+d；
19.步骤205：识别出事发地疏散范围外所有的庇护场所a
d
且对应的容量为p
d
且
20.优选的，所述步骤三包括：
21.通过公式获取人群速度和密度关系，其中，a表示高度自由行走速度，ρ
tran
表示人群运动速度减小时的人群密度，ρ
crit
表示人群临界密度，ρ
max
表示人群最大密度，ρ表示人群密度变量；
22.根据人群速度和密度关系式利用公式获取各个路段的疏散时间；
23.利用公式获取疏散人群经过某相邻路段的疏散效果，其中，α为第一权重取0.7，β为第二权重取值0.3；r1为所在伤害分区的评分值；r2为道路周边危化企业数量的对应评分值；γ为疏散人群通过上一交通控制点的时间，可连接道路交通系统实时获取，若无法实时获取，一般取30s。
24.优选的，所述步骤四包括：利用0
‑
1规划找出某一疏散出发点经过中间节点到对应某一疏散出口的所有路径，从中确定疏散效果最佳的路径。
25.优选的，所述步骤四还包括：
26.安全疏散路径是沿着方向进行，假设开始节点编号为1，到某个出口编号为9，路过节点为1，2,，3
‥
q；
27.对于任意路线x
kj
，k，j＝1,2
‥
q，如果经过线路x
kj
，则x
kj
＝1，否则x
kj
＝0，以任一节点为平衡节点，进入节点为+1，离开节点为
‑
1，则
28.对于起点节点1，所有节点都是离开此节点，所以
‑
x
14
‑
x
13
‑
x
12
＝
‑
1。
29.对于结束节点9，所有都是进入该节点，所以x
49
+x
89
＝1，故约束条件为：
[0030][0031]
满足约束条件的情况下，构建最佳疏散效果的0
‑
1规划模型：
[0032][0033]
根据最佳疏散效果的0
‑
1规划模型，获取某一疏散出发点对应所有的疏散出口a
e
的最佳疏散效果对应的路径，则必然存在最佳疏散效果对应的疏散出口使得出发点到疏散出口的疏散效果最佳，最终得到每一个疏散出发点所对应的最佳疏散出口
[0034]
优选的，所述步骤五包括：筛选距离疏散范围道路距离最近的庇护场所
对比其人口容量与待疏散人群数量的大小，若该庇护场所可以满足需求，则以该庇护场所作为本次安全疏散的终点；若该庇护场所的人口容量不能满足疏散需求，则继续筛选距离疏散范围道路距离较近的庇护场所进行比较，以此类推，直至最后选择的庇护场所人口容量满足疏散需求。
[0035]
本发明还提供一种燃气管线泄漏突发事件安全疏散系统，所述系统包括：
[0036]
影响范围获取模块，用于根据燃气泄漏事故后果模型获取燃气泄漏突发事件的各类后果影响范围以及伤害分区，根据各类后果影响范围以及伤害分区获取燃气泄漏突发事件的疏散范围；
[0037]
疏散处理模块，用于获取疏散出发点、疏散出口节点、警戒点、疏散路网节点以及庇护场所；
[0038]
疏散效果获取模块，用于获取疏散路网内各路段的疏散时间，根据各路段的疏散时间、通行道路所处的事故伤害分区及周边危化企业数量等信息获取疏散效果；
[0039]
最佳疏散效果模块，用于获取各疏散出发点经过中间节点到对应的疏散出口的所有路径，根据各路径的疏散效果确定最佳疏散效果路径；
[0040]
庇护场所获取模块，用于获取合适的庇护场所。
[0041]
优选的，所述疏散处理模块还用于：
[0042]
步骤201：根据地理信息系统识别出疏散范围内的各类建筑物位置信息，可由城市人口密度热力图或相关政府部门获取各类建筑物的人口数量，若无法获取，则根据建筑物的属性预估该影响范围内的人口数量；
[0043]
步骤202：以各类建筑物的出口作为疏散范围内人群疏散出发点a
o
且
[0044]
步骤203：识别出疏散范围内的所有道路交通控制点即交叉路口作为疏散路径的中间节点a
m
且识别出疏散范围外所有最近的道路交叉口作为疏散出口节点a
e
且中间节点a
m
、疏散出口节点a
e
即为燃气管线泄漏突发事件疏散需要设置警戒的警戒点；
[0045]
步骤204：由疏散的出发点、中间节点、疏散出口节点共同构成了燃气泄漏突发事件的疏散路网节点a
ω
且a
ω
＝{a
o
,a
m
,a
e
}，则用表示a
ω
内任意两个相邻节点之间的道路距离，其中a，p＝1,2,3
‥
，b+c+d；
[0046]
步骤205：识别出事发地疏散范围外所有的庇护场所a
d
且对应的容量为p
d
且
[0047]
优选的，所述疏散效果获取模块还用于：
[0048]
通过公式获取人群速度和密度关系式，其中，a表示高度自由行走速度，ρ
tran
表示人群运动速度减小时的人群密度，ρ
crit
表示人群临界密度，ρ
max
表示人群最大密度，ρ表示人群密度变量；
[0049]
根据人群速度和密度关系式利用公式获取各个路段的疏散时间；
[0050]
利用公式获取疏散人群经过某相邻路段的疏散效果，其中，α为第一权重取0.7，β为第二权重取值0.3；r1为所在伤害分区的评分值；r2为道路周边危化企业数量的对应评分值；γ为疏散人群通过上一交通控制点的时间，可连接道路交通系统实时获取，若无法实时获取，一般取30s。
[0051]
优选的，所述最佳疏散效果模块还用于：利用0
‑
1规划找出某一疏散出发点经过中间节点到对应某一疏散出口的所有路径，从中确定疏散效果最佳的路径。
[0052]
优选的，所述最佳疏散效果模块还用于：
[0053]
根据最佳疏散效果的0
‑
1规划模型，获取某一疏散出发点对应所有的疏散出口a
e
的最佳疏散效果对应的路径，则必然存在最佳疏散效果对应的疏散出口使得出发点到疏散出口的疏散效果最佳，最终得到每一个疏散出发点所对应的最佳疏散出口
[0054]
优选的，所述庇护场所获取模块还用于：筛选距离疏散范围最近的庇护场所且e＝1,2
…
,f对比其人口容量与待疏散人群数量的大小，若该庇护场所可以满足需求，则以该庇护场所作为本次安全疏散的终点；若该庇护场所的人口容量不能满足疏散需求，则继续筛选距离疏散范围道路距离较近的庇护场所且r＝1,2
…
,f进行比较，以此类推，直至最后选择的庇护场所人口容量满足疏散需求。
[0055]
本发明的优点在于：本发明综合考虑燃气突发事件的各类后果影响范围及伤害分区，将伤害分区叠加到疏散路径上，对各个疏散路径上的风险与时间进行综合性评估，进而得到最佳疏散效果的路径及疏散出口，并给出适合的庇护场所位置，从而实现高效、安全的人员疏散。
附图说明
[0056]
图1为本发明实施例所公开的一种燃气管线泄漏突发事件安全疏散方法的流程图；
[0057]
图2为本发明实施例所公开的一种燃气管线泄漏突发事件安全疏散方法中管道泄漏模型示意图；
[0058]
图3为本发明实施例所公开的一种燃气管线泄漏突发事件安全疏散方法中疏散路径的网络图。
具体实施方式
[0059]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出
创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0060]
实施例1
[0061]
如图1所示，一种燃气管线泄漏突发事件安全疏散方法，所述方法包括：
[0062]
步骤s1：根据燃气泄漏事故后果模型获取燃气泄漏突发事件的各类后果影响范围以及伤害分区，根据各类后果影响范围以及伤害分区获取燃气泄漏突发事件的疏散范围；具体过程为：
[0063]
(1)燃气泄漏事故模型分析
[0064]
1)根据燃气事故位置、泄漏类型描述，获取泄漏燃气管线压力、管径等基本参数，基于燃气突发事件后果危险最大化原则，根据管道泄漏模型计算燃气管道泄漏量。
[0065]
方法1：燃气管线泄漏时，管道内和泄漏孔内的气体分别存在两种流态，即临界流和亚临界流，最为普遍的情形为管道内亚临界流，孔口临界流或亚临界流。孔口亚临界流的情形一般是指随着泄漏孔径的增大，气体由临界流转为亚临界流，此外还包括管道断面破裂的情况。根据上述情形并结合经典的“孔”和“管道”模型开发了一种新模型，用于计算中低压下分配系统中的气体释放量。该模型兼顾了小孔泄漏和管道断裂泄漏两种情形。模型在使用时考虑到最大泄漏量的情况，即泄漏孔直径为管道直径时进行计算。
[0066]
图2为管道泄漏模型示意图，以下详细介绍具体计算过程：
[0067]
第一步：计算ma2(ma2>0)
[0068][0069][0070]
式中，d
d
为泄漏孔直径，m；当泄漏口直径缺省时，可取直径为d，管道直径，m；a
or
为泄漏孔面积，m2；a
c
为管道截面积，m2，k为气体绝热指数，默认取1.28；ma2为管道内泄漏口处的气流马赫数。
[0071]
第二步：计算中间参量y1，y2，
[0072][0073][0074]
式中，ma1为管道初始位置气流的马赫数；
[0075]
第三步：计算流体雷诺数：
[0076]
[0077]
式中，
[0078]
式中：ρ为流体密度，kg/m3；v为流体流速，m/s，对于城市天然气管道中压按12m/s，低压管道按6m/s计算；d为管道直径，m；μ为天然气动力粘度系数，取1.01
×
105pa
·
s；p1为管道初始点压力，pa；m为天然气摩尔质量，kg/mol，取0.1748kg/mol；r为通用气体常数，8.314j/(mol
·
k)；t1为初始点温度，k。
[0079]
第四步：计算管道范宁摩擦系数
[0080]
若re<100000，
[0081]
f＝0.079re
‑
0.25
ꢀꢀ
(6)
[0082]
若re≧100000，
[0083]
f＝0.0232re
‑
0.1507
ꢀꢀꢀꢀ
(7)
[0084]
式中，f为管道范宁摩擦系数；
[0085]
第五步：求解ma1(ma1>0)。初值定为0.0001，
[0086][0087]
式中，l为初始点到泄漏点的管道长度，m。
[0088]
第六步：根据泄漏孔面积a
or
计算质量流量q1为：
[0089][0090]
式中，q1为泄漏孔处气体的质量流量，kg/s；p1初始位置管道压力，pa；t1为初始位置气体温度，k。
[0091]
第七步，p2，t2计算
[0092][0093][0094]
式中，p2，t2分别为泄漏孔正下方管道内气体的压强和温度，pa，k。
[0095]
方法2：燃气管道泄漏速率的计算方法如下式：
[0096][0097]
式中：
[0098]
p
atm
——标准大气压，pa；
乡村条件f0.04x(1+0.0001x)
‑
1/2
0.016x(1+0.0003x)
‑1城镇条件a
‑
b0.32x(1+0.0004x)
‑
1/2
0.24x(1+0.001x)
‑
1/2
城镇条件c0.22x(1+0.0004x)
‑
1/2
0.20x城镇条件d0.16x(1+0.0004x)
‑
1/2
0.14x(1+0.003x)
‑
1/2
城镇e
‑
f0.11x(1+0.0004x)
‑
1/2
0.08x(1+0.0015x)
‑
1/2
[0114]
为简化计算，令he＝y＝z＝0，则
[0115][0116]
根据风险最大化原则，燃气泄漏扩散系数一般按照城镇条件f计算。
[0117]
3)取燃气扩散体积浓度值(可由质量浓度转换)不小于3％vol(甲烷爆炸下限)的范围，为人员疏散区域s1。为了防止发生火灾或爆炸事故，应该尽快通知危险区域内人员停止可能产生点火源的行为，并组织危险区域内居民疏散，基于此危险区边界浓度选择燃气爆炸下限浓度3％，通过将体积浓度转换为质量浓度c＝c
v
·
m/22.4，带入式16、式17即可得到疏散范围s1。
[0118]
(2)燃气泄漏事故后果模型分析
[0119]
1)喷射火模型s2[0120]
喷射火是指天然气管道发生泄漏，带压气体会从破裂处高速的喷出，遇有点火源形成的喷射火灾。喷射火主要通过热辐射作用周边的环境和人群，强烈的热辐射可能烧毁设备甚至造成人员伤亡等。
[0121]
假定火焰为圆锥形，并用从泄漏处到火焰长度4/5处的点源模型来表示。
[0122]
喷射火火焰长度可用如下方程得到
[0123][0124]
化简可得l
i
＝17q
0.444
ꢀꢀꢀꢀꢀꢀ
(19)
[0125]
点热源的热辐射通量按下式计算
[0126][0127]
式中：q为入射热辐射剂量，kw/m2；η为热辐射率，对于喷射火，取0.2或0.35，此处取0.35；τ
a
为大气传输率，τ可简化取1；q为气体释放率，kg/s；h
c
为燃烧热，对于甲烷，取55600kj/kg；r1为火源与被辐射点的水平距离，m。
[0128]
被辐射点至泄漏点距离可认为：
[0129][0130]
对于喷射火灾，根据表2不同事故伤害阈值，将热通量值带入上式即可得到不同伤害分区对应的影响范围；
[0131]
2)蒸气云爆炸模型s3[0132]
蒸气云爆炸是指燃气大量快速泄漏，与周围空气混合形成覆盖范围很大的“预混云”，在某一有限空间遇点火源而导致蒸气云爆炸。蒸气云爆炸产生的超压冲击波会对周边
的人员和造成一定的伤害。
[0133]
蒸气云的tnt当量：
[0134][0135]
式中m
f
为蒸气云中燃料的总质量，kg；m
f
＝q
×
t；a
tnt
为蒸气云的tnt当量系数，一般a
tnt
＝10％；q
tnt
为tnt的爆热，取q
tnt
＝4520kj/kg。对于质量为m
tnt
的tnt标准爆源，在地面发生爆炸时，爆炸场冲击波超压描述如下：
[0136][0137]
式中：r—爆炸场某点至爆源的距离，m；
△
p—爆炸波的入射超压，pa。
[0138]
根据表2不同事故伤害阈值，将超压值带入上式即可得到不同伤害分区对应的影响范围。
[0139]
3)火球模型s4[0140]
燃气泄漏爆炸释放巨大的能量加热周围的空气形成向外猛烈膨胀的火球。爆炸火球的伤害形式主要为皮肤伤害和视网膜烧伤，严重时能致人死亡。
[0141]
燃气在泄漏后如果在尚未与空气充分混合即被点燃，属于扩散性的燃烧，即形成球形或半球形的火体，称为火球。火球能够产生强烈的热辐射，对管道周围人员和财产构成伤害。火球燃烧可以认为是瞬态燃烧，因此采用热剂量准则。对处于火球内部的目标，显然将被热辐射毁伤，而对于火球以外的目标，其毁伤程度可以根据目标所接受的热毁伤阈值来判定，同样根据给定的阈值即可求出毁伤半径。
[0142]
距离火球中心某距离处的目标受到的入射热辐射剂量，具体为：
[0143][0144]
式中：q—入射热辐射剂量，kj/m2；tc为热传导系数，取1.0；r2为目标距离火球中心的距离，m；q
t
为火球的总辐射量，kj。可以用参与燃烧的泄漏天然气m1(质量为90％的泄漏量)燃烧释放的所有热量表示，即q
t
＝0.9
×
q
×
t
×
hc，t泄漏时间按照30s计算。
[0145]
根据表2不同事故伤害阈值，将热辐射值带入上式即可得到不同伤害分区对应的影响范围。
[0146]
4)地下空间爆炸模型s5[0147]
当燃气管线因地质灾害导致管线弯曲变形，发生微小泄漏后，泄漏气体存在通过土壤扩散至相邻地下空间聚集爆炸的可能。考虑到爆炸后果，地下空间主要指连通管线(污水管线、雨水管线、暗渠及同功能管线、电力管线)。结合相关研究，地下空间爆炸致灾模式以超压为主。
[0148]
当燃气管线因地质灾害导致管线弯曲变形，发生微小泄漏后，泄漏气体存在通过土壤扩散至相邻地下空间聚集爆炸的可能。考虑到爆炸后果，这里地下空间主要指连通管线(污水管线、雨水管线、暗渠及同功能管线、电力管线)线。根据相关实验，爆炸破片影响半径小于超压伤害半径，在此不再考虑破片影响。地下空间爆炸致灾模式主要为超压。
[0149]
以待评估燃气管段12.5m范围连通管线为计算对象。
[0150]
首先需要考虑燃气扩散至相邻空间概率，概率由下式计算
[0151][0152]
式中p
di
为燃气管线距离周边第i个连通管线的扩散概率，d
i
为燃气管线距第i个地下空间的最短距离。根据相关实验结果，超压对人员伤害半径为
[0153]
其中s
6i
为周边第i个地下空间(管线)截面积。
[0154]
对于连通管线截面积s
6i
：
[0155]
1)对于雨污水管线，其截面积可表示为：
[0156][0157]
式中r
y
为雨污水管线的直径。
[0158]
2)对于暗渠等长方形管道，其截面积可表示为：
[0159]
s
6i
＝ab
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(27)
[0160]
式中a，b分别为暗渠的宽，高。
[0161]
因此，利用燃气泄漏喷射火模型、蒸气云爆炸模型、火球模型、地下空间爆炸模型对人员的伤害程度来划分人员疏散范围。根据表2不同事故伤害阈值得出燃气爆炸对应的不同伤害后果范围s
2i
，s
3i
，s
4i
，s
5i
，(i＝1，2，3)，其中i＝1，2，3，分别对应死亡区、重伤区、轻伤区的范围，以轻伤区的边界作为人员疏散范围划分标准。则燃气泄漏突发事件的伤害分区可表示为：
[0162]
死亡区范围：s
h1
＝s
21
∪s
31
∪s
41
∪s
51
；
[0163]
重伤区范围：s
h2
＝s
22
∪s
32
∪s
42
∪s
52
；
[0164]
轻伤区范围：s
h3
＝s
23
∪s
33
∪s
43
∪s
53
；
[0165]
由此可得燃气泄漏突发事件的疏散范围为：
[0166]
s＝s1∪s
h3
＝s1∪s
23
∪s
33
∪s
43
∪s
53
[0167]
表2不同事故伤害阈值
[0168]
伤害分区热通量/(kw
·
m
‑2)超压δp/(mpa)热剂量/(kj
·
m
‑2)死亡区>37.5>0.075>592重伤区12.5～37.50.025～0.075392～592轻伤区1.6～12.50.01～0.025172～392安全区<1.6<0.01<172
[0169]
步骤s2：获取疏散出发点、疏散出口节点、警戒点、疏散路网节点以及庇护场所；具体过程为：
[0170]
步骤201：根据地理信息系统识别出疏散范围内的各类建筑物位置信息，可由城市人口密度热力图或相关政府部门获取各类建筑物的人口数量，若无法获取，则根据建筑物的属性预估该影响范围内的人口数量，具体参阅文献《王海波,蔡铭,钟舒琦.一种基于建筑
物属性的城市人口分布方法[b]或者专利文献cn201611086311.x,2016,11》；
[0171]
步骤202：以各类建筑物的出口作为疏散范围内人群疏散出发点a
o
且
[0172]
步骤203：识别出疏散范围内的所有道路交通控制点即交叉路口作为疏散路径的中间节点a
m
且识别出疏散范围外所有最近的道路交叉口作为疏散出口节点a
e
且中间节点a
m
、疏散出口节点a
e
即为燃气管线泄漏突发事件疏散需要设置警戒的警戒点；
[0173]
步骤204：由疏散的出发点、中间节点、疏散出口节点共同构成了燃气泄漏突发事件的疏散路网节点a
ω
且a
ω
＝{a
o
,a
m
,a
e
}，则用表示a
ω
内任意两个相邻节点之间的道路距离，其中a，p＝1,2,3
‥
，b+c+d；
[0174]
步骤205：识别出事发地疏散范围外所有的庇护场所a
d
且对应的容量为p
d
且
[0175]
步骤s3：获取疏散路网内各路段的疏散时间，根据各路段的疏散时间、通行道路所处的事故伤害分区及周边危化企业数量等信息获取疏散效果；具体过程为：
[0176]
在进行燃气管线泄漏突发事件安全疏散时，将疏散范围内的人群从a
o
疏散到a
e
所历经的道路a
o
a
e
，即为该疏散事件可供选择的疏散路径。为了有效评估每条路径的疏散效果，首先需要计算人群通过a
o
→
a
e
路径上各个路段的通行时间，再基于通行道路所处的事故伤害分区及周边危化企业数量等信息，计算单个路段的疏散效果，最后综合得到整条路径的疏散效果。
[0177]
人群密度反映的是一个空间内人员的数量和稠密程度，通常用单位面积上的人数表示。一般情况下，最大人群密度是7～8人/m2，这取决与人体的生理尺寸。人群的移动速度与人群密度有很大关系，当人群密度很小的时候，人群可以以最大的速度移动，随着密度的增加，人群的移动速度会减低，当人群密度继续增大，人群的移动速度达到临界最小值，甚至由于踩踏事故的发生，人群的移动速度甚至会降低为零。目前为止，还没有一个统一的关系式可以表达人群速度与人群密度之间的变化关系，根据目前的研究成果，为了确保疏散的安全有序性，人群速度和密度关系主要采用以下公式：其中，a表示高度自由行走速度，ρ
tran
表示人群运动速度减小时的人群密度，ρ
crit
表示人群临界密度，ρ
max
表示人群最大密度，ρ表示人群密度变量，a＝1.4m/s，ρ
tran
＝0.8人/m2，ρ
crit
＝a/2b＝2.8人/m2，ρ
max
＝a/b＝5.6人/m2。
[0178]
根据人群速度和密度关系式利用公式获取各个路段的疏散时间；
[0179]
单个路段的风险主要考虑两方面：路段所在的伤害分区和周边周边危化企业数量等信息，如下表所示：
[0180]
表3单个路段风险评估值
[0181][0182]
利用公式获取疏散人群经过某相邻路段的疏散效果，其中，α为第一权重取0.7，β为第二权重取值0.3；r1为所在伤害分区的评分值；r2为道路周边危化企业数量的对应评分值；γ为疏散人群通过上一交通控制点的时间，可连接道路交通系统实时获取，若无法实时获取，一般取30s。，值越小，则该路径的疏散效果越好。
[0183]
步骤s4：获取各疏散出发点经过中间节点到对应的疏散出口的所有路径，根据各路径的疏散效果确定最佳疏散效果路径；利用0
‑
1规划找出某一疏散出发点经过中间节点到对应某一疏散出口的所有路径，从中确定疏散效果最佳的路径。以下主要采用0
‑
1规划的方法，具体参阅文献《李映红.线性0
‑
1规划模型的排序解法[j].西南交通大学学报,2001.》以及《王裕婷.基于0
‑
1规划的高层建筑物疏散路径优化[j].设计与案例,2020,7.》。
[0184]
如图3所示，安全疏散路径是沿着方向进行，假设开始节点编号为1，到某个出口编号为9，路过节点为1，2,3
‥
q；
[0185]
对于任意路线x
kj
，k，j＝1,2
‥
q，如果经过线路x
kj
，则x
kj
＝1，否则x
kj
＝0，以任一节点为平衡节点，进入节点为+1，离开节点为
‑
1，则
[0186]
如图3所示，对于起点节点1，所有节点都是离开此节点，所以
‑
x
14
‑
x
13
‑
x
12
＝
‑
1。
[0187]
对于结束节点9，所有都是进入该节点，所以x
49
+x
89
＝1，故约束条件为：
[0188][0189]
满足约束条件的情况下，构建最佳疏散效果的0
‑
1规划模型：
[0190][0191]
根据最佳疏散效果的0
‑
1规划模型，获取某一疏散出发点对应所有的疏散出口a
e
的最佳疏散效果对应的路径，则必然存在最
佳疏散效果对应的疏散出口使得出发点到疏散出口的疏散效果最佳，最终得到每一个疏散出发点所对应的最佳疏散出口
[0192]
步骤s5：获取合适的庇护场所，具体过程为：筛选距离疏散范围道路距离最近的庇护场所对比其人口容量与待疏散人群数量的大小，若该庇护场所可以满足需求，则以该庇护场所作为本次安全疏散的终点；若该庇护场所的人口容量不能满足疏散需求，则继续筛选距离疏散范围道路距离较近的庇护场所进行比较，以此类推，直至最后选择的庇护场所人口容量满足疏散需求。
[0193]
通过以上技术方案，本发明提供的一种燃气管线泄漏突发事件安全疏散方法，在燃气泄漏事件发生进行人员安全疏散时，不仅要考虑路径最短、时间最少等问题，更要考虑人员在疏散道路上所承受的事故伤害影响，防止人员踏入风险程度更高的道路，从而造成更大的人员伤亡。燃气泄漏突发事件发生后，首要任务就是要把影响范围内的人群疏散出去，若在确定庇护场所时，只单一考虑疏散点到庇护场所的道路距离情况，很有可能造成疏散人群在影响范围内滞留时间相对过长，也就使得人群遭受伤害的可能性随之升高，因此确定疏散范围外最近的道路控制点作为中转节点就显得很有必要。基于此，本发明综合考虑燃气突发事件的各类后果影响范围及伤害分区，将伤害分区叠加到疏散路径上，对各个疏散路径上的风险与时间进行综合性评估，进而得到最佳疏散效果的路径及中转点，并给出适合的庇护场所，从而实现高效、安全的人员疏散。
[0194]
实施例2
[0195]
本发明还提供一种燃气管线泄漏突发事件安全疏散系统，所述系统包括：
[0196]
影响范围获取模块，用于根据燃气泄漏事故后果模型获取燃气泄漏突发事件的各类后果影响范围以及伤害分区，根据各类后果影响范围以及伤害分区获取燃气泄漏突发事件的疏散范围；
[0197]
疏散处理模块，用于获取疏散出发点、疏散出口节点、警戒点、疏散路网节点以及庇护场所；
[0198]
疏散效果获取模块，用于获取疏散范围内各路段的疏散时间，根据各路段的疏散时间、通行道路所处的事故伤害分区及周边危化企业数量等信息获取疏散效果；
[0199]
最佳疏散效果模块，用于获取各疏散出发点经过中间节点到对应的疏散出口的所有路径，根据各路径的疏散效果确定最佳疏散效果；
[0200]
庇护场所获取模块，用于获取合适的庇护场所。
[0201]
具体的，所述疏散处理模块还用于：
[0202]
步骤201：根据地理信息系统识别出疏散范围内的各类建筑物位置信息，可由城市人口密度热力图或相关政府部门获取各类建筑物的人口数量，若无法获取，则根据建筑物的属性预估该影响范围内的人口数量；
[0203]
步骤202：以各类建筑物的出口作为疏散范围内人群疏散出发点a
o
且
[0204]
步骤203：识别出疏散范围内的所有道路交通控制点即交叉路口作为疏散路径的
中间节点a
m
且识别出疏散范围外所有最近的道路交叉口作为疏散出口节点a
e
且中间节点a
m
、疏散出口节点a
e
即为燃气管线泄漏突发事件疏散需要设置警戒的警戒点；
[0205]
步骤204：由疏散的出发点、中间节点、疏散出口节点共同构成了燃气泄漏突发事件的疏散路网节点a
ω
且a
ω
＝{a
o
,a
m
,a
e
}，则用表示a
ω
内任意两个相邻节点之间的道路距离，其中a，p＝1,2,3
‥
，b+c+d；
[0206]
步骤205：识别出事发地疏散范围外所有的庇护场所a
d
且对应的容量为p
d
且
[0207]
具体的，所述疏散效果获取模块还用于：
[0208]
通过公式获取人群速度和密度关系式，其中，a表示高度自由行走速度，ρ
tran
表示人群运动速度减小时的人群密度，ρ
crit
表示人群临界密度，ρ
max
表示人群最大密度，ρ表示人群密度变量；
[0209]
根据人群速度和密度关系式利用公式获取各个路段的疏散时间；
[0210]
利用公式获取疏散人群经过某相邻路段的疏散效果，其中，α为第一权重取0.7，β为第二权重取值0.3；r1为所在伤害分区的评分值；r2为道路周边危化企业数量的对应评分值；γ为疏散人群通过上一交通控制点的时间，可连接道路交通系统实时获取，若无法实时获取，一般取30s。
[0211]
具体的，所述最佳疏散效果模块还用于：利用0
‑
1规划找出某一疏散出发点经过中间节点到对应某一疏散出口的所有路径，从中确定疏散效果最佳的路径。
[0212]
具体的，所述最佳疏散效果模块还用于：
[0213]
如图3所示，安全疏散路径是沿着方向进行，假设开始节点编号为1，到某个出口编号为9，路过节点为1，2,，3
‥
q；
[0214]
对于任意路线x
kj
，k，j＝1,2
‥
q，如果经过线路x
kj
，则x
kj
＝1，否则x
kj
＝0，以任一节点为平衡节点，进入节点为+1，离开节点为
‑
1，则
[0215]
对于起点节点1，所有节点都是离开此节点，所以
‑
x
14
‑
x
13
‑
x
12
＝
‑
1。
[0216]
对于结束节点9，所有都是进入该节点，所以x
49
+x
89
＝1，故约束条件为：
[0217][0218]
满足约束条件的情况下，构建最佳疏散效果的0
‑
1规划模型：
[0219][0220]
根据最佳疏散效果的0
‑
1规划模型，获取某一疏散出发点对应所有的疏散出口a
e
的最佳疏散效果对应的路径，则必然存在最佳疏散效果对应的疏散出口使得出发点到疏散出口的疏散效果最佳，由此最终可以得到每一个疏散出发点所对应的最佳疏散出口
[0221]
具体的，所述庇护场所获取模块还用于：筛选距离疏散范围最近的庇护场所且e＝1,2
…
,f对比其人口容量与待疏散人群数量的大小，若该庇护场所可以满足需求，则以该庇护场所作为本次安全疏散的终点；若该庇护场所的人口容量不能满足疏散需求，则继续筛选距离疏散范围道路距离较近的庇护场所且r＝1,2
…
,f进行比较，以此类推，直至最后选择的庇护场所人口容量满足疏散需求。
[0222]
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。