基金项目:中国地震局震灾应急救援司专题课题“地震应急公共服务平台研究”资助.
(Earthquake Administration of Yunnan Province,Kunming 650224,Yunnan,China)
the Analytic Hierarchy Process; earthquake-site communications; emergency response mode
备注
基金项目:中国地震局震灾应急救援司专题课题“地震应急公共服务平台研究”资助.
对云南地震现场应急通信技术系统的现状进行梳理分析,以AHP层次分析法为基础,从系统的设备构成、功能特点、制约因素三个方面,阐述了各层面间的相互联系与系统自身存在的制约因素,确定了各级指标的权重,从而得出这些因素对实现总目标的影响或重要程度,并对决策模型的计算结果进行了详细的分析,最后应用该方法对比分析近几年云南典型地震现场应急通信技术保障实例。
In this paper,we analyzed the System of Yunnan Earthquake-site Emergency Communications which is in operation in the Earthquake Administration of Yunnan Province. On the basis of the Analytic Hierarchy Process(AHP),and aiming at the system's equipment,functions and features,and restricting factors,we studied relations between each hierarchy of the system and the restricting factors of the system itself. Then we decided the weights of indicators at all levels. Thus,according to these weights,we concluded how the restricting factors influence the achievement of the overall goal,and how important they were to the overall goal. Further,we analyzed the computed results from the decision-making model that we proposed in this paper. We found our results practicable compared with the real cases of the dispatch and deployment of the emergency communications on the sites of the earthquake events in Yunnan Province in recent years.
引言
自2011年3月24日发生缅甸孟帕亚7.2级地震以来,云南地区已进入新一轮强震活跃期(皇甫岗等,2014)。据统计,2011—2015年,云南共发生5.0级以上破坏性地震21次,平均每年4次,地震共造成731人死亡,5 315人受伤,直接经济总损失4 129 728万元①②③④⑤。云南地处欧亚大陆与印度洋板块碰撞带的东沿,地壳运动活跃,区域内横断山脉纵横交错,地质切割形态分布较广,高原山地环境山高坡陡。因此,云南地震活动频度高、人员伤亡多、灾害损失重、救援难度大。这些都给云南地震现场通信技术系统保障带来了严峻的考验。
近年来,我国发生了震级较大的汶川8.0级、芦山7.0级、鲁甸6.5级破坏性地震,从历次地震现场的震害可知,在中等以上破坏性地震的震中区,有线通信网络基本破坏,在应急期内难以立即恢复(张毓丰,邓宪民,2005),地震现场通信保障工作在地震应急处置中显得尤为重要。随着政府对防震减灾工作越来越重视,社会公众对地震应急处置的关注和要求越来越高,地震应急通信保障任务也就越来越重。为了不断提高地震应急通信保障能力,使地震现场应急通信技术系统的指挥调度和应用更加科学,一些学者对地震应急指挥技术系统硬件和软件(帅向华等,2009a,b)、地震应急指挥协同平台(姜立新,帅向华,2012)、地震救援现场应急通信系统的设计(洪惠群等,2012)进行了深入的研究,并对国家地震应急指挥技术系统建设中的关键技术(杨天青等,2010)进行了探索和应用分析,这些研究成果为云南地震现场通信技术系统的应用提供了很好的参考。
《国家地震应急预案》和《云南省地震应急预案》对地震应急响应级别的启动作了规定,《云南省地震系统应急工作方案》对现场工作队通信组的地震现场通信和网络保障、现场灾情获取和信息传输等工作内容作了说明。赵恒等(2007,2012)对破坏性地震的应急通信需求进行分析,认为不同震级所造成的灾区通信破坏程度不同,因而不同应急响应级别对应的地震现场通信技术系统保障模式也会有差异。云南地震现场通信技术系统由多套子系统组成,设备较多、功能复杂,在应急响应指挥调度决策的过程中,难免出现功能重复、装备搭配不合理、地震现场保障模式和灾区通信需求不匹配等问题。为了使该系统在地震应急中的指挥调度和应用更加科学,最大化发挥应急通信保障能力,笔者利用AHP层次分析法,根据不同的应急响应级别,对该系统的功能和制约因素进行分析研究,探索适应云南特殊地理环境的地震现场通信技术系统应急响应模式。
1 云南地震现场通信技术系统概况
经过国家十五项目、地震社会服务工程、云南省预防和处置地震灾害能力建设十项重大措施和十项重点工程等项目的规划建设,云南省地震局建成了由现场应急指挥技术系统、现场工作队通信指挥系统、前方移动指挥平台、无人机地震灾情获取系统、现场指挥所共5套系统(平台)组成的云南地震现场应急通信技术系统(表1)。
在历次地震应急响应中,云南省地震局地震现场应急通信技术系统都进行了不同程度的保障。据统计,在2011—2015年发生的21次5.0级以上破坏性地震中,云南省地震局启动Ⅳ级响应3次,Ⅲ级响应11次,Ⅱ级响应4次,Ⅰ级响应1次。其中,现场应急指挥技术系统进行了6次保障,现场工作队通信指挥系统进行了14次保障,现场指挥所进行了11次保障。无人机地震灾情获取系统建成时间较晚(2015年建成),在笔者的统计时间段内没有得到较好的应用,但在2014年“8·03”鲁甸地震和2014年“10·07”景谷地震应急响应中地震系统的其他单位采用了无人机对灾区进行了航拍,灾情获取效果较好。
从表1中可以看出,各分系统整体的装备和功能各有侧重,但也有相同和相似的部分,有的可独立使用,有的需要相互配合使用,有的装备可满足多种功能,有的功能需要多种装备配合才能实现,有些系统和装备又有各自的制约因素。
2 AHP方法的基本原理
AHP决策分析法由美国运筹学家Saaty于20世纪70年代提出,是一种定性与定量相结合的层次决策分析方法。它将复杂问题分解为多个组成因素,并将这些因素按支配关系进一步分解,按目标层、准则层、指标层排列起来,形成一个多目标、多层次的模型,形成有序的递阶层次结构。通过两两比较的方式确定层次中诸因素的相对重要性,然后综合评估主体的判断确定诸因素相对重要性的总顺序。
2.1 构造判断矩阵构造判断矩阵是AHP方法的关键步骤,表示针对上一层次中的某元素而言,评定该层次中各有关元素的相对重要性,其形式如表2所示(徐建华,2002)。
若以矩阵来表示,即
B=[b11 b12 … b1n
b21 b22 … b2n
bn1 bn2 … bnn](1)
式(1)中,B称为判断矩阵。其中,bij表示对于Ak而言,元素Bi对Bj的相对重要性的判断值,bij一般取1、3、5、7、9等5个等级指标,意义为:1表示Bi与Bj同等重要,3表示Bi比Bj重要一点,5表示Bi比Bj重要得多,7表示Bi比Bj重要很多,9表示Bi比Bj极端重要。2、4、6、8为相邻判断的中值。显然bii=1,bij=1/bji,(i,j=1,2,…,n)。
2.2 模型计算第一步,计算层次单排序。对矩阵B计算满足
BW=λmaxW(2)
的特征根和特征向量,并将特征向量归一化。在式(2)中,λmax为B的最大特征根,W为对应于λmax的正规化特征向量,W的归一化特征向量就是对应元素的单排序权重值。
检验判断矩阵的一致性,计算一致性指标
CI=(λmax-n)/(n-1)(3)
式中:当CI=0,即λmax=n时,判断矩阵具有完全一致性,反之,CI越大,判断矩阵的一致性越差。
计算一致性比例:
CR=(CI)/(RI)<0.10(4)
其中,RI为平均随机一致性指标(表3)。一般当CR<0.10时,判断矩阵具有满意的一致性,否则就要调整判断矩阵。
第二步,计算层次总排序。若上一层次所有元素A1,A2,…,Am的总排序权重值为a1,a2,…,am; 与Aj对应的本层次元素B1,B2,…,Bn的层次单排序结果为[bj1,bj2,…,bjn]T(若Bi与Aj无联系时,bji=0); 那么,B层次的总排序为表4。
对层次总排序进行一致性检验,分别计算下列指标:
CI=∑mj=1ajCIj(5)
RI=∑mj=1ajRIj(6)
CR=(CI)/(RI)<0.10(7)
同样的,当层次总排序的随机一致性比例CR<0.10时,层次总排序具有满意的一致性,否则就要调整该层次的判断矩阵。
3 云南地震现场通信技术系统AHP模型
3.1 建立层次结构根据云南地震现场通信技术系统的结构,按照总目标、系统层、功能层、制约因素4个层次列出决策因子(表5)。根据表5所述各个决策因素之间的相互关系(表1),建立如图1、2所示层次结构模型。
3.2 构造判断矩阵笔者选取多次参加云南地震现场的相关专家和专门从事地震现场通信的技术人员,发放了15份调查问卷,其中专家问卷4份,一般技术人员问卷11份,对系统本身、系统功能和制约因素进行问卷调查,对应功能层次结构图(图1)和制约因素层次结构图(图2),构造了基于云南地震现场通信技术系统的3组共33个判断矩阵。
层 Z6 箱体较重,运输和搬运困难 Z7 数据流量较大,网络要求高 Z8 便携系统电力持续时间有限 Z9 空域申请困难 Z10 技术要求较高 Z11 工作场地环境易受限制
3.3 模型计算(1)计算5个系统X1,X2,…,X5的相对权重,系统层数为1,因此计算出的相对权重既是层次单排序,也是层次总排序(表6),他们表示各系统对实现总目标的重要程度。
(2)计算每个功能G1,G2,…,G22对每个系统的相对权重(层次单排序),并用系统X1,X2,…,X5的权重对功能的相应权重加权后相加,计算功能层的组合权重(表7),它们表示各功能对实现总目标的重要程度。
(3)计算每个制约因素Z1,Z2,…,Z11对每个系统的相对权重(层次单排序),并用系统X1,X2,…,X5的权重对制约因素的相应权重加权后相加,计算制约因素层的组合权重(表8),它们表示各制约因素对实现总目标的制约程度。
利用Matlab程序(余胜威,2014)对该AHP模型进行计算,结果如表6~8所示。
4 计算结果与实例分析
5 结语
地震应急指挥决策中,现场通信技术系统的指挥调度是一个复杂的管理问题,笔者基于AHP方法对地震现场通信技术系统的应急响应模式进行分析,把系统的各种构成要素分模块分层次区别开来,赋予一定的权重后排序,在系统指挥调度时可避免响应过度和响应不足,使得地震应急响应更加合理。地震现场受灾情况复,不同地震现场对通信系统功能的保障需求也不同,在充分估计地震灾情的前提下可根据该模型计算的功能权重排序考虑功能配置,使地震现场通信功能的保障更加全面。地震应急响应时要考虑系统本身状况、灾区自然环境、受灾情况等因素的协调性,尽量避免通信系统受到制约不能充分发挥保障功能的情形。本文为复杂的问题提供了一种分层次的权重决策分析方法,计算结果为多样本计算结果,具有普遍性,但在地震应急响应中,一些特殊地震要考虑不同的现场通信保障需求,决策者在普遍性和特殊性之间还需做一些调整。另外,该模型的结果主要适用于云南地区,对于特大型城市、平原地区、经济发达地区的地震应急通信响应模式还需考虑不同的影响因素。
4.1 计算结果分析4.1.1 Ⅰ级响应结果分析当启动Ⅰ级响应时,系统的权重排序为X5>X2>X4>X3>X1(表6),此时的地震现场应急通信技术系统保障可优先考虑无人机地震灾情获取系统和现场工作队通信指挥系统,其次考虑现场指挥所系统和前方移动指挥平台。
从系统功能层次看,功能层次权重位于前5位的是G5、G12、G22、G2、G1(表7),此时可优先考虑卫星语音通话、地震灾情数据采集传输、无人机地震灾情图像拼接和灾害识别、数据/网络、音视频会议方面的功能保障。功能层次权重位于6~10位的是G3、G20/G21、G9、G13、G6(表7),其次考虑地震现场办公、无人机地面目标探查跟踪和定位、无人机地面信息航拍采集、地震现场照明、地震科普宣传、临时现场指挥所方面的功能保障。
从制约因素分析,Ⅰ级应急响应对应的震级较高,灾情严重,地震灾区滑坡分布普遍,道路中断的概率很大,通信瘫痪的可能性也很大。根据AHP决策分析的结果,权重排序位于前5位的制约因素是Z3、Z2、Z5、Z11、Z9(表8),此时应充分考虑大型车辆难于进入道路较差的灾区、车体较重不平衡、互联网瘫痪不能为通信系统提供互联网通信、工作场地环境限制等因素对现场通信技术系统相关设备的制约。
4.1.2 Ⅱ级响应结果分析当启动Ⅱ级响应时,系统的权重排序为X2>X5>X3>X4>X1(表6),此时的地震现场应急通信技术系统保障可优先考虑现场工作队通信指挥系统和无人机地震灾情获取系统,其次考虑前方移动指挥平台和现场指挥所系统。
从系统功能层次看,功能层次权重位于前5位的是G1、G22、G2、G5、G12(表7),此时可优先考虑音视频会议、无人机地震灾情图像拼接和灾害识别、数据/网络、卫星语音通话、地震灾情数据采集传输方面的功能保障。功能层次权重位于6~10位的是G20/G21、G3、G7、G8、G6(表7),其次考虑无人机地面目标探查跟踪和定位、无人机地面信息航拍采集、地震现场办公、地震现场通信保障、现场指挥协同调度、临时现场指挥所方面的功能保障。
从制约因素分析,Ⅱ级应急响应对应的震级高、灾情重,地震灾区有一定滑坡分布,道路可能中断,极震区通信可能瘫痪。根据AHP决策分析的结果,权重排序位于前5位的制约因素是Z3、Z2、Z5、Z9、Z4(表8),此时应充分考虑大型车
辆难于进入道路较差的灾区、车体较重不平衡、互联网瘫痪不能为通信系统提供互联网通信、无人机空域申请困难、卫星天线直径小等因素对现场通信技术系统相关设备的制约。
4.1.3 Ⅲ级响应结果分析当启动Ⅲ级响应时,系统的权重排序为X1>X4>X2>X3/X5(表6),此时的地震现场应急通信技术系统保障可优先考虑现场应急指挥技术系统、现场指挥所和现场工作队通信指挥系统。
从系统功能层次看,功能层次权重位于前5位的是G1、G5、G2、G7、G12(表7),此时可优先考虑音视频会议、卫星语音通话、数据/网络、地震现场通信保障、地震灾情数据采集传输方面的功能保障。
从制约因素分析,Ⅲ级应急响应对应的震级低,灾情不太重,地震灾区滑坡可能性也较小。此时主要考虑系统设备自身的状况和山区道路对通信车通行的影响。
4.2 实例分析笔者选取了2011年以后具有代表性的云南3次不同应急响应级别的地震应急通信技术系统保障实例进行对比,分别是2014年“8·03”鲁甸6.5级地震 云南省地震局.2014.2014年8月3日云南鲁甸6.5级地震灾害直接经济损失评估报告(云南灾区).、2014年“10·07”景谷6.6级地震 云南省地震局.2015.2014年10月7日云南景谷6.6级地震灾害直接经济损失评估报告.、2015年“3·01”沧源5.5级地震 云南省地震局.2015.2015年3月1日云南沧源5.5级地震灾害直接经济损失评估报告(中国灾区).,对比情
况见表9。从对比结果看,Ⅰ级和Ⅱ级响应的吻合性较好,云南地震现场应急通信技术系统的装备属于重型装备,设计时基于大震应急考虑而设计,因而在发生大地震时,系统的适用性会比较明显。Ⅲ级响应的吻合性稍差,原因是Ⅲ级响应一般针对较小破坏性地震而启动,灾区受灾情况相对较轻,通过灾区当地均能保障办公功能(G3),构造判断矩阵时考虑该功能权重小一些,沧源地震应急中出现该功能的保障属于现场指挥所系统(X4)本身具备的功能,因而模型计算结果更具有普偏性。Ⅲ级响应下,通常会考虑设备老化(Z1)、车体较重不平衡(Z2)、车体较大难于进入道路较差的灾区(Z3)等一些系统设备自身存在的制约因素,考虑地震造成灾区地质灾害的影响制约因素的权重会小一些。而对比实例中出现了动中通天线直径小功率低(Z4)、数据流量较大网络要求高(Z7)应属特例,主要原因是沧源位于中缅交界少数民族聚居区,卫星信号值弱 美国波音卫星系统公司.2003.亚洲4号卫星ku波段FSS频段中国波束EIRP覆盖图.,互联网基础建设滞后。由表中的对比情况和上述分析可得,模型计算结果与实例对比具有一定的吻合性。
- 洪惠群,林金发,吴楠楠,等.2012.地震救援现场应急通信系统的设计与集成[J].山西地震,(3):28-30.
- 皇甫岗,苏有锦,张建国,等.2013.新世纪滇缅交界大震活动研究[M].昆明:云南科技出版社.
- 姜立新,帅向华,聂高众,等.2012.地震应急指挥协同技术平台设计研究[J].震灾防御技术,7(3):294-302.
- 帅向华,姜立新,刘钦,等.2009a.地震应急指挥技术系统设计与实现[J].测绘通报,(7):38-41.
- 帅向华,姜立新,王栋梁.2009b.国家地震应急指挥软件系统研究[J].自然灾害学报,18(3):99-104.
- 徐建华.2002.现代地理学中的数学方法(第二版)[M].北京:高等教育出版社.
- 杨天青,帅向华.2010.国家地震应急指挥技术系统建设中的关键技术及应用[J].震灾防御技术,5(2):208-214.
- 余胜威.2014.MATLAB优化算法案例分析与应用[M].北京:清华大学出版社.
- 张毓丰,邓宪民.2005.地震现场应急救援的通信系统设计[J].灾害学,20(4):111-115.
- 赵恒,白仙富,张方浩,等.2012.破坏性地震的应急通信需求与应用初探[J].地震研究,35(1):139-144.
- 赵恒,周兰江,曹刻,等.2007.云南省地震局地震应急通信系统在2007年宁洱6.4级地震中的应用与发展[J].地震研究,30(4):395-400.