Evaluation of the Application Efficiency of the Earthquake Emergency Collaborative Technology Platform
-
摘要: 根据地震灾害应急协同技术平台的保障条件、功能架构、性能指标、成果产出和实震检验等方面的应用情况和特点,采用德尔菲法从影响地震灾害应急协同技术平台状态和应用效能的诸多因素中,梳理出3项准则层指标、13项子准则层指标和52项具体指标,构建应用效能评估指标体系,并运用层次分析法确定各层级指标相对于上一层指标的重要程度排序及权重。运用模糊综合评价法取优、良、差3个等级,以云南省地震灾害应急协同技术平台为例开展应用效能评估。评估结果表明,该平台应用效能处于中等偏上的水平。实例应用结果表明,可利用德尔菲法、层次分析法和模糊综合评价法相结合的方式进行地震灾害应急协同技术平台应用效能评估。Abstract: According to the application and characteristics of the guarantee conditions, functional architecture, performance indicators, achievement output and actual earthquake inspection of the earthquake disaster emergency collaborative technology platform, the Delphi method was used to sort out 3 criterion-level indicators, 13 sub-criterion-level indicators and 52 specific indicators from the many factors that affect the status and application efficiency of the earthquake emergency collaborative technology platform. And an application efficacy evaluation index system is constructed. The analytic hierarchy process is used to determine the importance ranking and weight of each level index relative to the previous level index. The fuzzy comprehensive evaluation method is used to get three grades of “excellent”, “good”, and “poor”, and the application efficiency evaluation is carried out by taking the earthquake emergency collaborative technology platform in Yunnan Province as an example. The evaluation results show that the application performance of the platform is at an upper-middle level. The application results of the example show that the method proposed in this paper combines the Delphi method, the analytic hierarchy process and the fuzzy comprehensive evaluation method to evaluate the application efficiency of the earthquake disaster emergency coordination technology platform. The platform or system to carry out application performance evaluation provides universal experience and has certain reference significance.
-
Key words:
- Earthquake disaster /
- Collaborative service /
- Efficacy evaluation /
- Emergency
-
引言
道路是生命线系统的重要组成部分,道路系统不仅对国民经济发展至关重要,也是灾区震后得以及时救援和恢复重建的保障(刘金龙等,2013)。至2019年,中国公路总里程居世界第一,共484.7万km,其中高速公路达14.3万km,汶川地震灾区范围内公路总里程达62671km(刘爱文等,2008),受损公路总里程达31412km,近一半的公路受损,仅道路破坏带来的直接经济损失高达612亿。对地震引起的道路破坏机理进行分析,了解道路震害影响因素,从而针对特定的地质条件修复和加固道路,对提高道路抗震水平和震后恢复具有重要的现实意义(王伟等,2014;顾全等,2017;李帅等,2017)。
周德培等(2010)结合工程震害实例,根据震害现象分析各类边坡和相应支挡结构的震害机制;陈乐生(2012a)通过对汶川地震公路震害的调查统计,得到了地基条件、地基类型、所在位置、道路、断裂带等因素与路基破坏的关系;胡衡(2018)总结了路基、支挡结构和边坡震害,并给出了发生道路震害的主要因素。目前,关于道路震害分析的研究大都仅给出导致道路震害的主要因素。
本文选取道路构件(挡土墙、边坡和路基路面)典型破坏现象,分析和总结相应的震害特点,并将道路构件震害按破坏形式进一步将挡土墙分为墙身破坏和倾斜破坏两类;根据边坡坡度和岩石类型,将其分为崩塌型滑坡和塌陷滑移型滑坡两类;路基发生永久变形的三种情况,包括路基差异沉降、路堤边坡失稳、断层破裂引起的差异位移。总结每类破坏形式的常见工程和自然条件,讨论地震动导致道路构件破坏的原因,可加深对公路系统震损特征的了解,有利于提高道路抗震能力和震后恢复能力。
1. 挡土墙破坏机理分析
挡土墙指用于支承路基填土或山坡土体、防止填土或土体变形失稳的构造物。地震直接作用引起的挡土墙破坏包括垮塌、墙身剪断、整体倾斜、倾覆、开裂变形,挡土墙震害现象如图 1-7(中华人民共和国交通运输部等,2009;陈乐生,2012a),地震的间接破坏由次生灾害(如滑坡、土体崩塌等)导致挡土墙被掩埋或砸坏,所以可认为间接破坏的挡土墙是由边坡破坏引起的。汶川地震中挡土墙震害现象包括垮塌、墙身剪断、整体倾斜、倾覆、开裂变形、挡土墙被掩埋或砸坏等。
汶川地震中,烈度为Ⅵ-Ⅺ度的调查地区存在395处挡土墙破坏(陈乐生,2012a),直接破坏的震害数存在364处,占总震害数的92.2%,间接破坏的震害数仅占总震害数的7.8%。在直接破坏的挡土墙中,垮塌、墙身剪断、开裂变形破坏存在291处;73处为挡土墙倾斜破坏,所以按挡土墙墙身结构是否破坏分为挡土墙因荷载过大导致承载力不足而产生的墙身结构破坏和挡土墙整体倾斜或倾覆。
1.1 挡土墙墙身结构破坏
挡土墙垮塌、中部剪断、开裂和鼓胀裂缝均使墙身结构发生了破坏。不同砌筑方式挡土墙破坏数量调查结果显示,汶川地震调查的405处挡土墙震害中,发生破坏的浆砌挡土墙占80%,发生破坏的混凝土挡土墙占17%,发生破坏的干砌挡土墙占2%,发生破坏的其他类型挡土墙占1%。干砌挡土墙在汶川地震中普遍发生了破坏,仅部分浆砌挡土墙基本完好,可知干砌和浆砌挡土墙是挡土墙墙身破坏的主要类型。根据各地震烈度下挡土墙破坏数量统计结果(图 8),可知随着地震烈度的增大,挡土墙破坏数量呈增多的趋势,且破坏主要集中于高烈度区,Ⅷ度区挡土墙破坏数量少于Ⅶ度区,Ⅷ度区挡土墙数量少,且多为混凝土挡土墙,这说明混凝土挡土墙抗震性能较浆砌和石砌挡土墙好。
对浆砌和石砌挡土墙抗震性能进行分析,由于挡土墙墙身结构破坏是由挡土墙结构构件承受的应力过大引起的,如果重力式混凝土挡土墙结构设计恰当,对于静态荷载承载力计算,一般取荷载系数为1.7,因此对于一般强度的地震动,挡土墙不会发生墙身结构破坏,除非挡土墙内部存在缺陷(如截面不足、施工缝不良和缺乏加强内部支撑)。以浆砌、石砌施工方法修筑的挡土墙,在石块间的砂浆位置易形成软弱区,整体性较差,是墙身破坏的主要原因,山区挡土墙施工质量较差也是墙身易破坏的原因。
综上可知,在地震的反复作用下,挡土墙在结构缺陷部位形成裂缝,对于浆砌挡土墙,存在多处软弱部位,易发生墙面大面积开裂、鼓胀和隆起。随着地震动强度的增大,混凝土挡土墙裂缝处的抗剪切力不足抵抗裂缝上部的土压力,直接使挡土墙被剪断,而浆砌、干砌挡土墙因砂浆材料的全面屈服整体垮塌。
1.2 挡土墙滑移、倾斜或倾覆
挡土墙发生倾斜破坏时墙身整体结构较好,破坏特征为顶部位移过大,整体失稳。调查结果显示(陈乐生,2012a),挡土墙倾斜震害发生在岩石地基上的仅占22.4%,发生在土质地基上的高达58.4%,在119处路肩墙破坏结果中,挡土墙未在岩质地基上发生过度倾斜或倾覆破坏。由此可知,相比岩石地基,在土质地基上挡土墙更易发生破坏。考虑土质地基承载力较低,地震发生后挡土墙墙趾位置处的土压力增大,土质地基易发生变形,因此导致土质地基上的挡土墙更易发生滑移和倾斜破坏。
挡土墙试验证明(杨长卫等,2015),随着地震强度的增大,挡土墙背后土压力分布不均匀,土压力强度基本与加速度成正比,从整体土压力分布规律来看,加速度越大,土压力分布越接近三角形。这可能是因为挡土墙上部约束少,更易产生位移,而挡土墙下部被土体约束,位移小,上部产生位移后,下部将受到更大的土压力。对于高烈度区的地震动峰值加速度,当挡土墙正面墙趾处土压力不足时,挡土墙在墙背主动土压力作用下发生滑移;当挡土墙正面墙趾处土压力足够时,导致墙体内、墙趾和墙踵处弯矩显著增加,发生倾斜破坏。墙身高度对挡土墙破坏的影响见表 1,由长度占比可知,随着高度的增加,挡土墙越易发生破坏(图 4),这说明地震作用在挡土墙背后填土上时,挡土墙越高,地震惯性力的放大效应越明显,转动力矩越大。对于结构完好的挡土墙,随着地震强度的增大,挡土墙墙身越高,其上部位移越易积累,因而导致挡土墙发生倾斜甚至倾覆破坏。
表 1 G213都映高速挡土墙破坏情况Table 1. The damage situation of retaining wall on highway G213 from Dujiangyan to Yingxiu挡土墙墙身高度/m 总长度/m 破坏长度/m 破坏长度占比/% <6 3784.4 1131.5 0.299 6-8 3217.9 1030.0 0.320 8-10 1259.4 537.6 0.427 >10 488.2 229.0 0.469 2. 边坡破坏机理分析
汶川地震中,由边坡破坏造成多处交通堵塞,地震诱发滑坡和路堤失稳已引起公路系统大规模瘫痪,这种破坏不仅造成严重的直接经济损失(赵纪生等,2008),且阻塞道路,冲击和掩埋车辆,滑坡滚落的岩石将路基路面和挡土墙砸坏,对震后抢修造成了极大困难。汶川地震造成的典型边坡破坏现象如图 9-13所示(徐锡伟,2009)。
对发生破坏的边坡岩土类型进行调查,结果见表 2,由表 2可知岩质更易发生边坡震害。边坡坡度与震害数量的关系如图 14所示,由图 14可知,随着坡度的增加,边坡损坏数量显著增加,在坡度>45°的陡坡上,边坡震害较集中。边坡震害现象表明边坡地形结构与边坡破坏规模具有很大关系。综上分析,边坡岩土类型、边坡坡度、地形结构等因素对汶川地震边坡震害的影响较大,根据不同震害现象,按边坡震害形成模式分为崩塌型滑坡和塌陷滑移型滑坡。
表 2 边坡岩土类型与震害数量Table 2. Types of slope rock and quantity of earthquake damage边坡岩土类型 震害数量/个 占比/% 岩质 251 47.0 土质 154 28.8 上层土质下层岩质 129 24.2 2.1 崩塌型滑坡
崩塌型滑坡指岩石土体在自重或地震等其他作用下从陡峭山坡急速下落,有时伴随着破碎岩体沿边坡大规模下落现象。由图 9-12所示边坡情况可知,此类边坡破坏规模普遍较大,属浅层岩石破坏,发生在边坡较高的位置,且发生在坡度>45°的陡坡上。较坚硬的岩石易形成较大的陡坡,这些陡坡岩石属于硬性岩石。地震发生后,边坡上的岩石沿边坡下落至地面,形成堆积阻塞,掩埋道路(图 10、图 12)。玄武岩和花岗岩等发育良好的岩石,部分岩体甚至重心外伸,且风化程度严重、岩石间胶结强度弱,在地震作用下,因风化或雨水侵蚀存在薄弱部位的岩体进一步开裂,从陡坡滑落或坠落,较大的岩体形成落石(图 11)。
调查发现,震中附近约150m高的水坝记录到峰值加速度高达1.5g-2.0g,而相同地点水坝底地震动峰值加速度较小,高差悬殊、斜坡山脊等地形对地震动的放大效应更明显。此外,考虑震中区较长时间的强震动作用,分析认为强震动作用时间长也是震中区边坡发生严重破坏的原因。综上可知,地震对边坡具有以下影响:①边坡地带对地震动有放大效应,尤其是陡峭的斜坡;②强震作用时间越长,坡体破裂越严重。边坡崩塌与上盘效应也有密切关系,北川映秀段公路位于汶川地震断层上盘,以花岗岩等硬性岩石为主,岩石受风化影响严重,边坡岩体发生了大规模松动和崩塌。
综上所述,边坡发生大规模崩塌滑坡现象的原因是发育良好的硬性岩石在严重的风化影响下,岩层层理间的胶结强度降低,产生微小裂缝,尤其是竖向层理的边坡岩体裂缝较明显。地形高差悬殊、坡度较陡的边坡对地震动强度有放大效应,遭遇强震时,地震波在已开裂的岩层中不断发生折射和反射,加大了破裂面上的拉应力,且岩层间产生的位移随之累积。这会导致结构面岩层强度进一步削弱,最终超过岩层破裂强度。上述破坏过程不仅发生在破裂界面上,岩石内部也受上述因素影响,导致岩体从震裂到松弛,抗剪能力较弱的岩层首先发生断裂,破裂的岩石从陡坡上被挤出,下落过程较迅速。
2.2 塌陷滑移型滑坡
塌陷滑移型滑坡指斜坡上岩石土体在重力或地震等其他作用下,沿着滑动面移动形成的滑坡。土体塌陷(陈乐生,2012b)是发生在冲积或海洋洪积平原天然沉积层内的破坏,常出现在溪流、渠道沿岸等地势较陡峭的边坡上,有时也出现在有浅层或上层滞水面的中等陡峭边坡上。图 13所示的岩土塌陷滑移,滑动岩体基本保持为整体,调查发现,发生此类破坏的边坡主要集中于坡度为30°-45°相对陡峭的边坡上,因坡度较崩塌型滑坡破坏的边坡小,移动速度相对较慢。通过滑动岩体基本保持了其相对位置,可知滑动面在地震发生前即为相对软弱的部位。
在调查的179个塌陷边坡中,塌陷破坏主要集中于Ⅸ-Ⅺ度区,约占破坏总数的80%,可知塌陷滑移型滑坡受地震强度的影响较大。通过岩石性质调查,可知地震诱发的土体塌陷多发生于土质边坡上,这些岩石多为弱胶结沉积岩,且岩体为软弱和软硬互层。部分边坡破坏发生在河道或沼泽沉积物等软土地基上,这些地基压实度不足,且高路堤人工填土较松散,所以压实性及地基条件对地震边坡破坏的影响较大。
由以上震害特点,对塌陷滑移型滑坡进行理论计算分析,为简化计算将边坡转化为二维平面,将地震动作为拟静力处理,水平和竖直两个方向的地震动影响相互独立,认为水平地震动系数kH用0.1、0.25和0.4分别表示小震(7度)、中震(8度)和大震(9度),相应的竖向地震动系数kV用0.05、0.125和0.2表示。如图 15所示,α表示滑动体滑动面与边坡坡面之间的角度;β表示水平向与边坡坡面的角度;FN表示对滑动体的支持力;FS表示平行于滑动面的力。
对滑动体受力分析,建立沿滑动面平行和垂直方向建立平衡方程。
$$\left\{\begin{array}{l}F_{\mathrm{N}}=\left[\left(1-k_{\mathrm{V}}\right) \cos \alpha-k_{\mathrm{H}} \sin \alpha\right] \times W \\ F_{\mathrm{S}}=\left[\left(1-k_{\mathrm{V}}\right) \sin \alpha+k_{\mathrm{H}} \sin \alpha\right] \times W\end{array}\right. $$ (1) 式中,为方便计算将滑动体简化为三角形,W为滑动体自重,用W=γ(tan α-H2/tan β)表示;γ为滑动体重度。
沿滑动面向上的力由土体的粘聚力和摩擦力提供,表示为:
$$F_{\mathrm{S}}=\left(c l+F_{\mathrm{N}} \tan \varphi\right) $$ (2) FS为土体表面具有的抗力;cl为破坏面粘聚力总和;FN为土体摩擦力;φ为滑动面上的土体内摩擦角。
当滑动体处于临界状态(即有滑动趋势的岩体在粘聚力以及静摩擦作用下恰好保持静止),FS用公式(3)表示为:
$$F_{\mathrm{S}}=\left(c l+F_{\mathrm{N}} \tan \varphi\right)=\left[\left(1-0.5 k_{\mathrm{H}}\right) \sin \alpha+k_{\mathrm{H}} \sin \alpha\right] \times \gamma\left(H^{2} / \tan \alpha-H^{2} / \tan \beta\right) $$ (3) 通常将竖向地震动系数用水平地震动系数乘以相应的系数表示,即kV=ξkH,最终可以
得出:
$$\begin{aligned} k_{\mathrm{H}}=&\left[\frac{2 c}{\gamma H} \sin \beta+(\cos \alpha \tan \varphi-\sin \alpha) \times \sin (\beta-\alpha)\right]/ \\ &\{[(\cos \alpha+\sin \alpha \tan \varphi)-\xi(\sin \alpha-\cos \alpha \tan \varphi)] \times \sin (\beta-\alpha)\} \end{aligned} $$ (4) 式(4)为边坡将要发生滑移时的临界条件,改变边坡坡度、岩土的性质与地震动强度可能会打破平衡条件,当右侧值大于左侧值时,边坡会发生破坏。
综上分析,塌陷滑移型滑坡发生破坏时,软弱或软硬互层岩土边坡在重力作用及雨水等因素的影响下,产生相对薄弱层,有滑动趋势的岩体在接触面抗剪切力及静摩擦力作用下保持静止;当地震动强度超过某临界值时,薄弱层抗剪切力和摩擦力不足以承受滑动体重力分量及加速度,滑动体沿薄弱层被拉裂,薄弱层产生贯穿的裂缝,滑动体整体沿裂缝向下滑移。滑坡规模与边坡坡度有关,在相对陡峭的边坡上,随着坡度的增加,滑坡规模越来越大。
3. 路基路面破坏机理分析
随着大地震的发生,道路路面往往发生破坏。地层变形或地层破坏是导致道路破坏的直接原因(李杰,2012),由于强烈的地震动,地基局部隆起或面层断裂时有发生,但相比其他形式的破坏,此类震害现象较少。汶川地震造成的路基路面常见震害现象如图 16-20所示(中华人民共和国交通运输部等,2009)。
地震导致的面层破坏不仅是铺砌表面的路面破坏,多由路基永久性地层变形造成。稳定且不易受地震诱发永久变形的路基支撑上的面层,在地震作用下较少破坏。在地基或路堤特别脆弱的位置,路基发生严重破坏,无承载能力的面层也会随之破坏。基于汶川地震出现的路基路面破坏现象,根据路基产生永久变形的原因,将路基破坏分为以下情况:
(1)压实度不足导致路基差异沉降
图 16所示的路基差异沉降中,裂缝延伸至边坡一侧,路基类型为半填半挖型。在汶川地震调查的区段中,Ⅵ-Ⅺ度区内路基路面震害共579处,Ⅸ-Ⅺ度区震害较集中,震害受损程度较大,共498处,占总震害数量的86%,地震强度对路基的影响较大。路基破坏情况如图 21、图 22所示,由图可知半填半挖是发生破坏的主要路基类型,土质是发生破坏的主要路基地质条件。
考虑半填半挖路基人工填筑部分及土质地基压实度不足且承载力较低,易发生塑性变形。由此认为,路基震害除与地震烈度有关外,下卧层存在易受地震诱发地层位移影响的压实不良的路堤材料或压实度存在明显差异的路基也易产生差异沉降破坏。综上分析,在地震动的反复作用下,路基根据压实度不同发生不同的塑性形变,压实度不足的路基沉降较大,导致上层路面随之破坏,如地震发生后路桥结合处因压实度明显不同造成错台现象和半填半挖路基局部沉降错台现象。
(2)路堤边坡失稳引起路基破坏
图 17、图 18所示为路堤边坡失稳引起的面层竖向位移震害现象,由图 18可知,路基在临路堤边坡一侧大范围塌陷,失去了支撑作用的路面不足以保证车辆安全通行。汶川地震中,仅受地震作用而破坏的路基较少,道路严重破坏多由浅层路基边坡塌陷所致。调查显示土质路基和半填半挖形式路基震损情况最多,且多为下边坡一侧临空路基下滑,导致路基发生严重毁坏,沿河公路出现此类破坏较多。分析认为半填半挖路基存在天然软弱结合面,沿河路堤边坡在河水的冲刷作用下易造成底部空虚,这是路堤边坡失稳的主要原因。
综上可知,此类破坏类似边坡塌陷滑移破坏,不合理填筑的过陡路堤及土质路基在雨水等外因影响下底部产生较小的破裂层或天然存在薄弱层,在强烈的地震作用下,破裂层或薄弱层处裂缝进一步开展,抗剪强度不足时,路堤边坡发生整体塌滑,进而失稳,造成路基竖向剪切破坏,有时会造成路基大范围整体滑移。小范围的路堤边坡塌陷往往仅使面层发生位移和错动,大规模路堤边坡塌陷或深层塌陷,均会造成大规模路基竖向断裂,有时路基和路面沉降量高达几米。
(3)断层破裂引起差异位移
考虑断层破裂附近的道路破坏严重,且常出现较大裂缝(图 19),有时甚至出现高差几米的陡坎(图 20)。断层走向与公路走向的夹角统计结果显示,随着平行至垂直的变化,路基破坏数量逐渐减少,当与断层走向平行时,破坏集中,且出现多处隆起、错台破坏。根据断层作用方式的不同,将断层破裂对附近路基的影响分为以下方面:①对于走滑断层,断层的相对位移直接导致断层延伸至地表处,产生较大裂缝甚至高差悬殊的陡坎。此外,根据错台发生的不同位置可知主断层破裂会导致附近次生断层的出现。利用有限元模型对断层导致的公路破坏情况进行模拟(邓龙胜等,2009),结果表明断层导致路基路面的破坏模式包括张拉、剪切、弯压和复合破坏,主要以张拉破坏为主。②断层不仅对断裂带产生影响,走滑断层上下盘效应的影响也不能忽视。断层上盘效应对上断层PGA有明显的增大作用(王栋,2010;范优铭等,2017),距断层裂缝越近,对道路上盘效应的影响越大。北川至映秀镇公路位于断层上盘,与下盘距断层相同距离的道路相比,震害较严重。由此可知,对于斜向走滑断层,主断层面两侧的相对移动导致道路路基被地震裂缝拉断,与此同时,主断层断裂导致次生断层破裂,随着断层的相对移动,路基随之产生位移,最终出现路基隆起、陡坎等破坏。对于水平方向的断层,上部路基路面整体随着地震产生位移,但不产生相对位移,因此水平方向的断层造成路堤破裂破坏的情况少。断层上盘效应使位于断层上盘的路基承受更大的地震动峰值加速度,相当于增大了地震动强度,从而使压实度不足的路基更易发生不均匀沉降,使土质路基边坡易发生失稳塌陷和整体滑移破坏。
4. 结语
本文针对汶川地震震损道路进行机理分析,通过总结道路构件(挡土墙、边坡和路基路面)破坏类型,对破坏形式进行合理分类,给出每类破坏形式常见的自然地质条件和工程因素,得出以下结论:
(1)对挡土墙墙身破坏和倾斜倾覆破坏进行分析,对于浆砌和石砌挡土墙结构,忽略施工质量和材料自身等因素,易在石块间的砂浆位置产生软弱区,这是挡土墙开裂、墙身剪断甚至整体垮塌的主要原因。混凝土挡土墙一般不会发生结构自身破坏,除非挡土墙内部存在缺陷,包括截面不足、施工缝不良和缺乏内部支撑。挡土墙倾斜破坏常发生在土质地基中,因地震动作用于挡土墙背后回填土上,导致主动土压力增大,墙趾在基础中被约束,且土质地基易变形,挡土墙越高墙趾处的转动力矩越大,挡土墙上部易产生位移,导致挡土墙倾斜,整体失稳。
(2)根据不同边坡岩石类型和坡度产生的破坏现象,分为以下破坏模式:①崩塌型滑坡,发生大规模崩塌型滑坡现象的边坡具有硬质岩石、地形高差悬殊、坡度较陡、风化严重、处于断层上盘位置的特点。在地震作用下,抗剪能力较弱的岩层发生断裂,破裂的岩石从陡坡上被挤出。②塌陷滑移型滑坡,常见的人工填土中,发生大部分破坏的为填土松散,且压实度较差,或被填在河道、沼泽沉积物等软土地基上,在地震作用下,坡度为30°-45°的陡坡上中软岩石边坡在较深层沿软弱面断开,滑动体内部保持相对位置,整体下滑。
(3)路基产生永久变形的情况包括:①路基差异沉降,土体在地震作用下,压实度不同的位置发生差异塑性变形;②路堤边坡失稳,类似边坡塌陷滑移破坏形式,在强烈的地震作用下,路基破裂层或薄弱层处裂缝进一步开展,路堤边坡发生整体塌滑,进而失稳,造成路基竖向剪切破坏;③断层造成路基永久变形,斜向走滑断层两侧有相对移动,导致道路路堤被地震裂缝拉断,且断层上盘效应明显增大地震动强度。对于道路路基路面,虽可通过分析进行抗震设计,可在震前有减小或消除大多数面层破坏的改造或加固方法,但震后修复价格往往令人无法接受,且震前改造或加固仅对面层结构具有可行性,所以建议采用震后面层破坏快速修理的被动策略,而不是在预防路基破坏上加大投入的主动策略。
-
表 1 标度定义
Table 1. The scale definition
标度aij 含义 1 表示元素Ai与Aj相比一样重要 3 表示元素Ai与Aj相比略为重要 5 表示元素Ai与Aj相比明显重要 7 表示元素Ai与Aj相比非常重要 9 表示元素Ai与Aj相比绝对重要 2,4,6,8 为以上两判断之间的中间状态对应的标度值 表 2 平均随机一致性指标
Table 2. Mean values of the stochastic consistency indicator
指标 数值 n 1 2 3 4 5 6 7 8 9 RI 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 表 3 A-B判断矩阵及其权重
Table 3. A-B judgment matrix and its weights
影响因素 影响因素 权重W 一致性检验结果 B1 B2 B3 B1 1 1/2 1/3 0.163 4 CR=0.007 9<0.1,
通过一致性检验B2 2 1 1/2 0.297 0 B3 3 2 1 0.539 6 表 4 B1-C判断矩阵及其权重
Table 4. B1-C judgment matrix and its weights
子准则层 子准则层 权重W 一致性检验结果 C11 C12 C13 C11 1 1/2 1/2 0.2 CR=0<0.1,
通过一致性检验C12 2 1 1 0.4 C13 2 1 1 0.4 表 5 B2-C判断矩阵及其权重
Table 5. B2-C judgment matrix and its weights
子准则层 子准则层 权重W 一致性检验结果 C21 C22 C23 C24 C25 C21 1 1/3 1/2 1 1/3 0.097 6 CR=0.021 2<0.1,
通过一致性检验C22 3 1 1 3 1 0.269 9 C23 2 1 1 2 1/3 0.184 2 C24 1 1/3 1/2 1 1/4 0.092 1 C25 3 1 3 4 1 0.356 2 表 6 B3-C判断矩阵及其权重
Table 6. B3-C judgment matrix and its weights
子准则层 子准则层 权重W 一致性检验结果 C31 C32 C33 C34 C35 C31 1 1/2 1/4 1/3 1/5 0.062 7 CR=0.014 5<0.1,
通过一致性检验C32 2 1 1/3 1/2 1/4 0.099 4 C33 4 3 1 2 1/2 0.268 5 C34 3 2 1/2 1 1/2 0.177 1 C35 5 4 2 2 1 0.392 4 表 7 C11-D判断矩阵及其权重
Table 7. C11-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D111 D112 D113 D111 1 3 2 0.547 8 CR=0.013 2<0.1,通过一致性检验 D112 1/3 1 1/2 0.150 7 D113 1/2 2 1 0.301 5 表 8 C12-D判断矩阵及其权重
Table 8. C12-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D121 D122 D123 D124 D125 D121 1 1/2 4 2 3 0.265 0 CR=0.037 8<0.1,
通过一致性检验D122 2 1 5 3 4 0.420 0 D123 1/4 1/5 1 1/2 1/3 0.061 9 D124 1/2 1/3 2 1 1/2 0.112 7 D125 1/3 1/4 3 2 1 0.140 4 表 9 C13-D判断矩阵及其权重
Table 9. C13-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验 D131 D132 D133 D134 D131 1 3 3 4 0.512 0 CR=0.030 1<0.1,
通过一致性检验D132 1/3 1 2 2 0.224 6 D133 1/3 1/2 1 2 0.158 8 D134 1/4 1/2 1/2 1 0.104 5 表 10 C21-D判断矩阵及其权重
Table 10. C21-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D211 D212 D213 D214 D215 D211 1 1/3 1/4 1/2 1/3 0.073 5 CR=0.007 4<0.1,
通过一致性检验D212 3 1 1/2 2 1 0.215 4 D213 4 2 1 3 2 0.375 0 D214 2 1/2 1/3 1 1/2 0.120 8 D215 3 1 1/2 2 1 0.215 4 表 11 C22-D判断矩阵及其权重
Table 11. C22-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D221 D222 D223 D224 D225 D221 1 1/2 1 1/2 2 0.157 8 CR=0.003<0.1,
通过一致性检验D222 2 1 2 1 3 0.297 9 D223 1 1/2 1 1/2 2 0.157 8 D224 2 1 2 1 3 0.297 9 D225 1/2 1/3 1/2 1/3 1 0.088 5 表 12 C23-D判断矩阵及其权重
Table 12. C23-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D231 D232 D233 D234 D231 1 2 2 3 0.423 1 CR=0.038<0.1,
通过一致性检验D232 1/2 1 1 2 0.227 4 D233 1/2 1 1 2 0.227 4 D234 1/3 1/2 1/2 1 0.122 2 表 13 C24-D判断矩阵及其权重
Table 13. C24-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D241 D242 D243 D241 1 2 1 0.400 0 CR=0<0.1,
通过一致性检验D242 1/2 1 1/2 0.200 0 D243 1 2 1 0.400 0 表 14 C25-D判断矩阵及其权重
Table 14. C25-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D251 D252 D253 D254 D255 D251 1 3 2 2 3 0.3682 CR=0.003<0.1,
通过一致性检验D252 1/3 1 1/2 1/2 1 0.1094 D253 1/2 2 1 1 2 0.2065 D254 1/2 2 1 1 2 0.2065 D255 1/3 1 1/2 1/2 1 0.1094 表 15 C31-D判断矩阵及其权重
Table 15. C31-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D311 D312 D313 D311 1 1/2 1/3 0.168 0 CR=0.013 2<0.1,
通过一致性检验D312 2 1 1/2 0.277 3 D313 3 2 1 0.554 7 表 16 C32-D判断矩阵及其权重
Table 16. C32-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D321 D322 D323 D324 D321 1 3 3 4 0.515 0 CR=0.007 6<0.1,
通过一致性检验D322 1/3 1 1 2 0.190 0 D323 1/3 1 1 2 0.190 0 D324 1/4 1/2 1/2 1 0.105 1 表 17 C33-D判断矩阵及其权重
Table 17. C33-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D331 D332 D333 D334 D331 1 1/3 1 1/4 0.109 4 CR=0.007 6<0.1,
通过一致性检验D332 3 1 3 1/2 0.296 7 D333 1 1/3 1 1/4 0.109 4 D334 4 2 4 1 0.484 5 表 18 C34-D判断矩阵及其权重
Table 18. C34-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D341 D342 D343 D341 1 2 3 0.524 0 CR=0.013 2<0.1,
通过一致性检验D342 1/2 1 2 0.317 4 D343 1/3 1/2 1 0.158 7 表 19 C35-D判断矩阵及其权重
Table 19. C35-D judgment matrix and its weights
指标层 指标层 权重W 一致性检验结果 D351 D352 D353 D354 D351 1 1/2 3 2 0.277 6 CR=0.011 5<0.1,
通过一致性检验D352 2 1 4 3 0.466 8 D353 1/3 1/4 1 1/2 0.095 3 D354 1/2 1/3 2 1 0.160 3 表 20 地震灾害应急协同技术平台应用效能评估指标在各隶属度等级的概率
Table 20. The probability for different membership degrees of application efficiency evaluation index system of the collaborative technology platform for earthquake emergency
指标层 隶属度 优 良 差 标准 概率 标准 概率 标准 概率 应急指挥场所
具备情况D111有专用场所 1 与其他功能共用场所 0 无场所 0 应急指挥场所建筑物设防标准D112 高于本地最新设防要求 1 达到本地最新设防要求 0 低于本地最新设防要求 0 应急指挥场所可
容纳人数D11350人及以上 1 30~50人 0 30人以下 0 网络信道D121 具备互联网、地震行业网、电子政务网3种信道联通 1 具备互联网、地震行业网、电子政务网中的2种信道联通 0 仅具备互联网联通 0 运算设备D122 5~10台专业图形工作站进行评估运算,有数据库服务器和应用服务器 0.71 2~4台专业图形工作站进行评估运算,有1种服务器 0.29 没有或仅有1台专业图形工作站进行评估运算,无服务器 0 打印条件D123 可打印各类尺寸图幅(A0、A1、A3、A4) 1 可打印常用尺寸图幅(A1、A3、A4) 0 仅能打印A3、A4图幅 0 显示条件D124 高精度大屏显示 0 中低分辨率大屏显示,有辅屏显示 1 仅有电脑显示 0 控制设备D125 有中控系统,可远程控制指挥场所所有音视频信号,且可预制模板,一键切换 1 无中控系统,仅在指挥场所内控制音视频信号 0 多源信号控制困难 0 数据库建设D131 已建成且运维正常 1 有,但存在运维问题 0 没有 0 丰富性D132 基础地理信息、地震行业信息、社会经济信息等9大类42子类数据齐全 0 基础地理信息、地震行业信息、社会经济信息等9大类42子类数据基本齐全 1 基础地理信息、地震行业信息、社会经济信息等9大类42子类数据有较多缺失 0 颗粒度D133 可按州(市)、县、乡(镇)、村公里格网统计 0 可按州(市)、县、乡(镇)公里格网统计 1 可按州(市)、县公里格网统计 0 数据时效性D134 实时更新 0 1~3年更新 1 ≥5年更新 0 友好性D211 软件系统易理解、易操作、人工干预 0.29 软件系统较易理解、操作不复杂,不可人工干预 0.71 软件系统难理解,不易操作 0 时效性D212 软件系统在2~5 min内完成启动、计算、传输、显示 0.57 软件系统在5~10 min内完成启动、计算、传输、显示 0.43 软件系统完成启动、计算、传输、显示时间超过10 min 0 稳定性D213 软件系统全年24小时运行,各功能完全正常运行 0.71 软件系统全年24小时运行,各功能基本正常运行 0.29 软件系统部分功能经常不能正常运行 0 可扩展性D214 软件系统可拓展性强,可根据用户需求对系统功能、产出结果进行灵活开发调整 0.29 软件系统具有一定可拓展性,可根据用户需求对部分系统功能进行开发调整 0.71 软件系统部分功能、产品不可调整 0 安全性D215 系统设置了防火墙、访问控制机制、身份认证机制和数据加密机制等安全机制,采用网络监控与入侵防范系统实现入侵检测、病毒和木马扫描、安全扫描、安全审计 1 系统设置了防火墙、访问控制机制、身份认证机制和数据加密机制等安全机制 0 系统仅有身份认证等较少的安全机制 0 自动响应功能D221 自动获取地震三要素,并完成启动、计算、传输、显示 0.57 自动获取地震三要素,并完成启动、计算,需要人工传输和显示功能 0.29 不可自动获取地震三要素,需人工完成启动、计算、传输、显示功能 0.14 快速评估功能D222 在2~5 min完成1次7级地震评估,且产出评估结果完备 1.00 在5~30 min完成1次7级地震,且评估结果完备 0 完成1次7级地震时间超过30 min,评估结果完备 0 人工修正功能D223 平台收到修正信息后,可通过手机远程自动对评估结果进行修正 0.29 平台收到修正信息后可人工对评估结果进行修正 0.43 不能实现平台收到修正信息后对评估结果的修正 0.28 信息发布功能D224 预评估结果自动采用多种手段进行信息发布 0.43 预评估结果自动采用单一手段进行信息发布 0.29 预评估结果仅可通过快速人工手动进行信息发布服务 0.28 产品展示功能D225 实现对基础数据与评估结果数据的可定制化可视化展示,可适应不同显示场景,界面美观 0.43 实现对基础数据与评估结果数据的固定场景可视化展示,且可视化界面美观 0.29 未实现数据远程展示功能 0.29 地震台网数据
自动获取分析修正
功能D231实现速报台网和预警台网产出余震信息、仪器烈度等数据自动获取,并自动分析修正评估影响场 0 实现速报台网和预警台网产出余震信息、仪器烈度等数据自动获取,但须人工分析修正评估影响场 0 未实现速报台网和预警台网产出余震信息和仪器烈度等数据自动获取分析功能 1 震后趋势分析数
据的自动获取分析
修正功能D232自动获取震源机制解、破裂过程、余震时空分布等震后趋势分析数据,并自动根据获取的震后趋势分析数据进行快速动态修正 0 自动获取震源机制解、破裂过程、余震时空分布等震后趋势分析数据,但须人工根据获取的震后趋势分析数据进行快速动态修正 0 未实现震后趋势分析数据的自动获取分析修正功能 1 多源灾情上报
分析功能D233实现移动终端、PC端、APP等多种方式的灾情上报,并对上报的灾情信息进行自动分析,自动灾情研判生成地震影响范围 0.29 实现移动终端、PC端、APP等多种方式的灾情上报,须人工对上报的灾情信息进行统计分析,手动灾情研判生成地震影响范围 0.57 仅实现单一手段的灾情上报功能,对上报信息未进行任何处理 0.14 动态大数据获取
与自动分析修正
功能D234实现互联网(微博等)或物联网(通信基站等)大数据自动获取信息,进行灾情研判与自动分析修正预评估结果 0.29 实现少量或单一动态网大数据获取,自动灾情研判,须手动分析修正预评估结果 0.71 未实现网络大数据获取与自动分析修正功能 0 数据管理功能D241 实现对平台数据库中数据的统计、查询、管理、审核等操作,并对统计、查询结果进行可视化显示与数据导出 1 实现对平台数据库中数据的统计、查询、管理、审核等操作,不可对统计、查询结果进行可视化显示与数据导出 0 无数据管理操作功能 0 接口管理功能D242 实现对平台全功能接口的代码调试、接口设置、接口授权、接口统计等功能 0.43 实现对平台部分功能接口的代码调试、接口设置、接口授权、接口统计等功能 0.29 无法对接口进行管理 0.29 用户管理功能D243 实现对平台用户分级、分层管理,可对软件系统功能权限进行设置管理 0.71 实现对平台用户分级管理,无权限控制管理 0.29 无用户管理功能,仅可对用户信息进行查询 0 系统并发人数D251 ≥50人同时登录访问 0.29 20~50人同时登录访问 0.43 ≤20人同时登录访问 0.28 视频系统连通数D252 实现4个及以上视频系统的互联互通 0.29 实现2、3个系统视频的互联互通 0.71 仅与行业内进行互联互通 0 资源实时共享D253 实现数据、计算资源、计算模型、产出结果模板等资源的实时共享 0.57 实现数据、计算资源、计算模型、产出结果模板等部分资源的实时共享 0.43 不可进行资源实时共享 0 专家协同功能D254 实现地震行业专家文字、语音、视频等多方式的在线协同应急功能 0.29 实现地震行业专家较单一方式
的在线协同应急功能0.71 未实现行业专家系统应急功能 0 部门协同功能D255 实现抗震救灾指挥部成员单位文字、语音、视频等多方式的在线协同应急功能 0.29 实现抗震救灾指挥部成员单位较单一方式的在线协同应急功能 0.71 未实现部门系统应急功能 0 采用技术手段D311 应用系统(APP)、微信、短信推送,全部技术手段 0.29 2种技术手段 0.43 单一技术手段 0.29 服务产品模式D312 电子网页版、纸质版、短信版全部服务模式 0.43 2种服务模式 0.57 单一服务模式 0 服务对象范围D313 行业用户、主要指挥部成员单位、信息发布部门,全部服务对象 0.29 地震行业和主要指挥部成员单位 0.71 单一服务对象 0 产品推送时间D321 产品推送到服务对象的时间,产出后自动推送,10 min内 0.43 人工推送,30 min内 0.57 人工推送,超过30 min 0 产品持续更新D322 72 h内按照应急响应规程和实际的应急救援需求全时程动态更新 0.57 48 h内按照应急响应规程全时程动态更新 0.43 在规定时间内仅能静态产出,无法动态更新 0 专家协同时效D323 1 h内专家快速查询、修改、反馈结果 0.86 2 h内专家查询、修改、反馈
结果0.14 4 h内专家查询、修改、反馈结果 0 部门协同时效D324 2 h内迅速联通多部门视频会议,进行协同 0.71 4 h内逐步联通多部门视频会议,进行协同 0.29 8 h内联通多部门视频会议,进行协同或无法连通 0 专题图件质量D331 图件主题突出、要素完整、配色规范、图幅合适 0.71 图件主题基本突出、要素基本完整、配色基本规范、图幅基本合适 0.29 图件主题不突出、要素不完整 0 评估报告质量D332 评估报告结果准确,评估要素完整 0.29 评估报告结果基本准确,评估要素基本完整 0.71 评估报告结果偏差大,评估要素不完整 0 对策建议产出
质量D333对策建议具有针对性和可行性 0.57 对策建议具有一定的针对性和可行性 0.43 对策建议的针对性和可行性较差 0 灾情研判的
可靠性D334灾情研判结果可靠,可提供决策建议 0.29 灾情研判结果基本可靠 0.71 灾情研判结果可靠性较差 0 地震行业用户D341 满意度90%以上 0.71 满意度60%~90% 0.29 满意度60%以下 0 抗震救灾指挥部成员单位用户D342 满意度85%以上 0.57 满意度60%~85% 0.43 满意度60%以下 0 信息发布部门D343 满意度80%以上 0.57 满意度60%~80% 0.43 满意度60%以下 0 烈度评估结果D351 烈度等级、范围、方向基本一致 0.57 烈度等级、范围、方向有偏差,但在误差范围内 0.43 烈度等级、范围、方向偏差较大 0 人员伤亡评估
结果D352评估结果基本对应,起到了重要的决策支撑作用 0.57 评估结果有误差,与应急响应较匹配 0.29 评估结果存在数量级上的误差 0.14 经济损失评估
结果D353评估结果基本对应,起到了重要的决策支撑作用 0.57 评估结果有误差,与应急响应较匹配 0.43 评估结果存在数量级上的误差 0 应急响应等级
建议D354建议等级与实际相符 0.71 建议等级与实际基本相符 0.29 建议等级与实际相差较大 0 表 21 子准则层评价结果
Table 21. Evaluation results at the secondary rule level
子准则层 评价结果 优 良 差 C11 1.00 0 0 C12 0.77 0.23 0 C13 0.51 0.49 0 C21 0.66 0.34 0 C22 0.60 0.22 0.18 C23 0.10 0.22 0.68 C24 0.77 0.17 0.06 C25 0.51 0.28 0.21 C31 0.32 0.63 0.05 C32 0.57 0.43 0 C33 0.36 0.64 0 C34 0.65 0.35 0 C35 0.59 0.34 0.07 表 22 准则层评价结果
Table 22. Evaluation results at the criterion level
准则层 评价结果 优 良 差 B1 0.71 0.29 0 B2 0.50 0.25 0.25 B3 0.52 0.45 0.03 -
曹彦波, 李兆隆, 李永强等, 2015. 云南地震应急快速评估模型本地化集成研究. 地震研究, 38(1): 148—154 doi: 10.3969/j.issn.1000-0666.2015.01.020Cao Y. B. , Li Z. L. , Li Y. Q. , et al. , 2015. Research on the localization integration of Yunnan earthquake emergency rapid assessment model. Journal of Seismological Research, 38(1): 148—154. (in Chinese) doi: 10.3969/j.issn.1000-0666.2015.01.020 邓树荣, 余庆坤, 张方浩等, 2021. 基于ADC方法的云南地震现场指挥部视频会议集成系统应用效能评估. 地震研究, 44(2): 292—300Deng S. R. , Yu Q. K. , Zhang F. H. , et al. , 2021. Effectiveness evaluation of the integrated video conference system in earthquake field headquarters based on ADC method. Journal of Seismological Research, 44(2): 292—300. (in Chinese) 邓玉莲, 2021. 城市排水管网状态和运行效能评估方法的研究与应用. 北京: 北京建筑大学.Deng Y. L. , 2021. Research and application of evaluation method for state and operation efficiency of urban drainage network. Beijing: Beijing University of Civil Engineering and Architecture. (in Chinese) 姜立新, 帅向华, 聂高众等, 2012. 地震应急指挥协同技术平台设计研究. 震灾防御技术, 7(3): 294—302Jiang L. X. , Shuai X. H. , Nie G. Z. , et al. , 2012. Study on the design of earthquake emergency command collaboration technology platform. Technology for Earthquake Disaster Prevention, 7(3): 294—302. (in Chinese) 李罡风, 赵建和, 朱生水等, 2010. 安徽地震前兆台网监测效能综合分析. 地震研究, 33(3): 280—286 doi: 10.3969/j.issn.1000-0666.2010.03.007Li G. F. , Zhao J. H. , Zhu S. S. , et al. , 2010. Synthetic study of the observational efficiency of Anhui earthquake precursory network. Journal of Seismological Research, 33(3): 280—286. (in Chinese) doi: 10.3969/j.issn.1000-0666.2010.03.007 李敏, 郑川, 李永强, 2018. 云南省省、市、县级地震应急信息协同技术系统的设计与实现. 地震研究, 41(1): 139—147 doi: 10.3969/j.issn.1000-0666.2018.01.018Li M. , Zheng C. , Li Y. Q. , 2018. Design and implement of research on earthquake emergency information collaborative technology system of provincal-municipal-county in Yunan province. Journal of Seismological Research, 41(1): 139—147. (in Chinese) doi: 10.3969/j.issn.1000-0666.2018.01.018 李兆隆, 吕佳丽, 郑川等, 2019. 云环境下的云南地震应急指挥技术软件系统设计与实现. 地震研究, 42(2): 196—203Li Z. L. , Lv J. L. , Zheng C. , et al. , 2019. Design and implementation of Yunnan earthquake emergency command system in cloud environment. Journal of Seismological Research, 42(2): 196—203. (in Chinese) 曲凯, 陈西选, 徐珞等, 2017. 可定制的效能评估系统研究与实现. 软件工程, 20(2): 19—21, 11 doi: 10.3969/j.issn.1008-0775.2017.02.006Qu K. , Chen X. X. , Xu L. , et al. , 2017. Research and implementation of the customizable effectiveness evaluation system. Software Engineering, 20(2): 19—21, 11. (in Chinese) doi: 10.3969/j.issn.1008-0775.2017.02.006 徐国祥, 2008. 统计预测和决策. 3版. 上海: 上海财经大学出版社.Xu G. X. , 2008. Statistical forecasting and decision-making. 3 rd ed. Shanghai: Shanghai University of Finance & Economics Press. (in Chinese) 徐建华, 2002. 现代地理学中的数学方法. 2版. 北京: 高等教育出版社.Xu J. H. , 2002. Mathematical methods in contemporary geography. 2 nd ed. Beijing: Higher Education Press. (in Chinese) 张方浩, 李永强, 余庆坤等, 2015. 省级地震应急信息公共服务平台设计探讨. 震灾防御技术, 10(3): 657—663Zhang F. H. , Li Y. Q. , Yu Q. K. , et al. , 2015. Design of the public service platform for earthquake-emergency information at the provincial level. Technology for Earthquake Disaster Prevention, 10(3): 657—663. (in Chinese) 张方浩, 李永强, 曹彦波等, 2019. 基于云技术的云南地震现场应急指挥技术系统优化. 地震研究, 42(2): 257—264Zhang F. H. , Li Y. Q. , Cao Y. B. , et al. , 2019. Optimization of the Yunnan post-earthquake field emergency command technical system based on cloud technology. Journal of Seismological Research, 42(2): 257—264. (in Chinese) 张杰, 唐宏, 苏凯等, 2009. 效能评估方法研究. 北京: 国防工业出版社. 周森, 2012. 基于层次分析与概率计算综合法的广州市白云区岩溶塌陷风险分析. 广州: 华南理工大学.Zhou S. , 2012. Risk analysis of the karst collapse in Guangzhou Baiyun district based on analytical hierarchy process combined with probability calculation. Guangzhou: South China University of Technology. (in Chinese) Saaty T. L., Kearns K. P., 1985. The analytic hierarchy process. In: Saaty T. L., Kearns K. P., eds., Analytical Planning. Amsterdam: Elsevier. -