多情景下航空公司航线网络抗毁性研究_胡小兵.pdf
《多情景下航空公司航线网络抗毁性研究_胡小兵.pdf》由会员分享,可在线阅读,更多相关《多情景下航空公司航线网络抗毁性研究_胡小兵.pdf(7页珍藏版)》请在咨信网上搜索。
1、文章编号:1009-6094(2023)07-2223-07多情景下航空公司航线网络抗毁性研究*胡小兵1,魏媛2,李航1(1 中国民航大学安全科学与工程学院,天津 300300;2 中国民航大学电子信息与自动化学院,天津 300300)摘要:航空运输易受多种外部因素干扰,评估不同情景下航线网络的运行水平对有效减轻外部干扰具有重要意义。运用复杂网络方法,分别构建国内 3 家航空公司(A、Z、M)无向加权航线网络。基于网络拓扑结构,进一步考虑空间灾害影响,利用加权度中心性、加权介数中心性和机场灾害风险度识别重要节点,并设计随机、蓄意和自然灾害 3 种灾害情景,借助网络效率和最大连通子图相对大小评估
2、航线网络抗毁性。结果表明:3 家航空公司的航线网络对随机灾害具有较强的抗毁性,而对蓄意灾害表现出脆弱性;其中,基于加权介数中心性的蓄意灾害对航线网络的破坏性最大,且蓄意灾害下 Z 航线网络更容易瘫痪;暴雨对 3 家航空公司航线网络的破坏力大于沙尘暴大于随机灾害,且自然灾害下 M 航线网络更容易瘫痪。关键词:安全工程;航空运输;加权网络;随机灾害;蓄意灾害;自然灾害;抗毁性中图分类号:X949文献标志码:ADOI:10.13637/j issn 1009-6094.2022.0665*收稿日期:2022 04 13作者简介:胡小兵,教授,从事复杂系统安全和风险、智能计算、空中交通管理研究,hux
3、btg ;李航(通信作者),讲师,博士,从事复杂网络系统安全和风险研究,lih cauc edu cn。基金项目:天津市教委科研计划项目(2021KJ045)0引言航空运输是国家和地区经济快速发展的基础,但同时,航空运输易受多种不确定因素影响,这些突发事件极有可能导致航线停飞或机场关闭,从而影响航空运输网络效率;且网络面向突发事件的抗毁能力是评估航空运输系统安全的核心要素。因此,研究多情景下航空运输网络抗毁性具有重要意义。近年来,很多学者借助复杂网络理论对航空网络展开研究,当前研究主要集中在两方面:一是关于航空网络拓扑结构的研究1 2,大量研究发现我国航空网络具有小世界特征3,度分布服从幂律分
4、布4;二 是 有 关 航 空 网 络 特 性 的 研 究,如 抗 毁性5 6。所谓抗毁性,目前公认的是 1997 年 Ellison等7 提出的:指网络系统在遭受攻击、故障和意外事故时仍能够及时完成其关键任务的能力。2000 年,Albert 等8 最先开始复杂网络抗毁性研究,研究发现,随机攻击下,无标度网络比随机网络具有更强的抗毁性;在选择性攻击下,无标度网络要比随机网络更容易崩溃。在此基础上,不同学者对不同网络展开研究。2012 年,曾小舟9 对中国航空网络结构及抗毁性进行了实证研究,发现中国航空网络是小世界网络,且对随机攻击具有较强的抗毁性,而对蓄意攻击具有脆弱性,整个网络的可靠性由少数
5、关键机场决定。卓志强等10 分析了“一带一路”沿线航空网络特征并构建了随机攻击和蓄意攻击的仿真系统。柴星思11 研究中国航空网络发现,蓄意攻击下,移除约 21%的节点后,网络结构迅速瓦解,直至移除约 80%的节点后,网络整体运输效率完全丧失。但上述研究对象多为不同区域范围内的网络或单独某家航空公司网络12 13,较少将多家航空公司航线网络进行比较并评估其抗毁性。虽然,也有学者们对航空公司航线网络进行分析,如罗心雨等14 比较了随机攻击和蓄意攻击下国有三大航司航空网络的抗毁性。但这些研究主要基于随机和蓄意攻击两种情景,对现实灾害情景的考虑不够,如恶劣天气等。而在航线网络拓扑结构基础上,进一步考虑
6、灾害本身特点,从整体上分析系统的抗毁性水平对提升航空运输网络服务的稳定性和可靠性具有重要意义15。鉴于此,本文根据2020 年民航业发展统计公报16,对比各航空公司运输周转量和运输规模后,从国内主要航空公司中选取 3 个代表性航空公司A、Z、M,全面分析和比较多种灾害情景下,3 家航空公司航线网络的抗毁性水平,以期为航空公司网络优化提供一定的理论依据。1加权航线网络及抗毁性为系统研究航空公司航线网络抗毁性,本文将从 3 个方面入手进行分析:首先,构建无向加权航线网络模型;其次,确定航空运输可能遭受的灾害情景;最后,选取网络抗毁性测度指标,用于评估不同灾害情景下航线网络的抗毁性。1.1航线网络建
7、模考虑到航班量是衡量航线重要性的关键因素,为更准确地分析航空公司航线网络,构造加权航线网络模型。对于包含 N 个机场的民航运输系统,可3222第 23 卷第 7 期2023 年 7 月安全 与 环 境 学 报Journal of Safety and EnvironmentVol 23No 7Jul,2023建立其加权航线网络 G 为:G=(V,L,W)。其中,V代表网络的节点集,即网络中所有的机场;L 表示网络中所有航线组成的边集合。用 N N 邻接矩阵 A表示航线网络中各机场间航线连接情况,其中 A(i,j)=1 表示机场 i、j 之间有直通航线,否则为 0。W为航线间的航班量,即权值,可
8、用矩阵 W=(Wmn)表示,其中 m,nN。1.2航线网络拓扑结构特性分析1.2.1度与加权度节点度是指目标节点与其他节点存在的连边数量,能够反映节点在局部区域内的重要性,度越大表示节点越重要。这里,将节点 i 的加权度定义为与该节点相连的所有边的权值之和,记为 DF(i),计算方法如下。DF(i)=Nj=1Wij(1)平均加权度k是指网络中所有节点加权度的平均值,计算式为k=Ni=1DF(i)N(2)1.2.2聚类系数节点 i 的聚类系数 C(i)是指网络中节点 i 的 k个邻居节点间实际存在的航线数与总的最大可能航线数之比,计算方法如下。C(i)=2Eiki(ki 1)(3)式中Ei表示节
9、点 i 的 ki个相邻节点间实际存在的航线数。整个网络的聚类系数 C 是指网络中所有节点聚类系数的平均值,计算为C=1NC(i)(4)1.3灾害情景确定在建立航线网络的基础上,要研究航线网络抗毁性,需要先确定航空运输可能面临的灾害情景。航空运输易受多种因素干扰,随机攻击和蓄意攻击是研究复杂网络抗毁性时常用的两种基本方式,为更全面地评估航线网络的抗毁性,本文进一步引入对航空运输正常服务造成经常性影响的自然灾害情景,如灾害性天气。然后,要确定具体情景的节点移除顺序,具体方式如下。1)随机灾害情景17:随机的删除网络中的某些节点。随机灾害下节点被删除的顺序与其影响力无关,每个节点被删除的概率相同。2
10、)蓄意灾害情景17:指按照某种策略对节点进行重要性排序,优先删除对网络影响力较大的节点。本文利用加权度中心性和加权介数中心性对节点进行重要性排序,依次对节点进行攻击。以两城市间的航班量作为权重,基于传统节点重要度评价指标进行加权改进。a.加权度中心性。度中心性是刻画节点中心性最直接的度量指标,通常一个节点与网络中其他节点连接越多越重要。传统节点度中心性指通过计算网络中各节点的度,并由大到小排序,排名越靠前的节点越重要。基于航班量对传统节点度指标进行加权改进,见式(1),于是,加权度中心性指基于上述加权方法得到加权度排名并进行节点重要性分析。b.加权介数中心性。节点介数中心性用网络中经过某个节点
11、的最短路径的数目反映节点重要性。传统上,通常将节点 i的介数中心性定义为节点 i 通过网络中所有起点 终点(Origin-Destination,OD)对最短路径数与网络总路径数的比值。于是,基于航班量的加权介数中心性 BF(i)可表示为BF(i)=sitGFst(i)Fst(5)式中Fst表示节点 s 到节点 t 所有基于航班量加权的最短路径数,Fst(i)表示节点 s 到节点 t 中经过节点 i 的相应加权最短路径数。实际上是通过用各条最短路径上的航班量对不同最短路径进行加权,即运用加权最短路径区分不同最短路径的重要性。3)自然灾害情景:指根据机场灾害风险度对节点进行重要性排序,依次对节点
12、进行攻击。机场灾害风险度是指在特定时间和特定区域内,可能发生的特定自然灾害对机场造成预期损失的程度。依据灾害系统论,机场灾害风险度取决于致灾因子危险性、承载体脆弱性和设防水平等因素,这里,将机场 i 的灾害风险度(i)表示为(i)=HiVCi1Fi(6)式中H 表示致灾因子危险性,H=fs,f 为灾害发生频率,s 为灾害强度,这里将灾害强度 s 简单定为 1,直接用灾害发生频率反映灾害的危险性;VC表示机场脆弱性,这里用机场暴露的旅客量,也即机场吞吐量近似反映;F 反映某一机场的设防水平,这里用机场自身等级(l)与当地国内生产总值(Gross DomesticProduct,GDP,PGD)的
13、乘积反映,即 Fi=liPGDi,机场设防水平越高,受自然灾害影响越小,机场的灾害风4222Vol 23No 7安全 与 环 境 学 报第 23 卷第 7 期险度也越低,因此,与 F 呈反比的关系。1.4航线网络抗毁性测度确定航空运输面临的灾害情景后,需要提出适合的航线网络抗毁性测度指标,以分析系统在不同灾害情景下的抗毁性水平。网络效率和最大连通子图是常用的网络性能指标。本文通过分析网络遭受攻击前后网络性能指标的变化情况评估航线网络抗毁性。1)网络效率。网络效率指节点之间最短路径倒数和的平均值,用 E 表示,数学表达式为E=1N(N 1)i,jV;ij1dij(7)式中E 的取值范围为 0,1
14、;dij表示机场 i 到机场 j 的最短路径长度。网络连通性与网络效率呈正相关。令E=Ebe Eaf(8)式中Ebe表示网络遭受攻击前的网络效率;Eaf表示网络遭受攻击后的网络效率;E 表示网络遭受攻击前后网络效率 E 的下降幅度,E 越大,表明网络效率下降幅度越大,网络损毁程度越大,网络越容易瘫痪。图 1航线网络拓扑图Fig 1Airline network topology2)最大连通子图相对大小。最大连通子图指网络遭到攻击后分裂出的所有子图中节点数最多的子图。最大连通子图相对大小是指最大连通子图中节点数占原始网络节点数的比例,用 S 表示,数学表达式为S(G)=nN(9)式中N 为网络中
15、节点总数;n 为移除节点后最大连通子图所包含的节点数量。网络遭受攻击前后 S 的下降幅度可直观地反映网络被破坏的程度。2案例研究通过比较各航空公司航空运输周转量和运输规模,从国内 64 家运输航空公司中选取 3 个代表 A、Z、M 进行仿真试验。2.1数据收集及处理航线网络构建:运用世界航空运输资讯数据库(Cirium Dashboard)收集了 2020 年 8 月 3 家航空公司(A、Z、M)的航班数据,并基于此构建无向加权航线网络。需要说明的是,这里仅考虑直飞和大陆航线。图 1 为国内 3 家航空公司航线网络拓扑图,总体来看,3 家航空公司航线东西分布不均,呈现“东密西疏”的特点。利用
16、MATLAB 软件对 3 家航空公司航线网络进行指标计算,得到各航空公司航线网络基本特征值,见表 1。通过相关特征值对比发现,M 的通航城市和航线数最多,网络规模最大;A 的通航城市和航线数最少,聚类系数最大,表明其整体运行效率最高;M 的平均加权度高于 A 和 Z,为 427.61,表明在1 个月内,每个机场与其他机场平均有 427 班的航班量。机场灾害风险度确定:这里,主要考虑暴雨和沙尘暴两类灾害情景。本文首先通过中国气象数据网收集了 20112020 年的两类灾害发生频次数据,并表 1三家航空公司航线网络基本特征Table 1Basic characteristics of the ro
17、ute networkof the three major airlines航空公司节点数边数平均加权度聚类系数A 航120356309.650.393 5Z 航143536423.930.169 2M 航147614427.610.224 052222023 年 7 月胡小兵,等:多情景下航空公司航线网络抗毁性研究Jul,2023对空间分布的频次等级进行简单划分。总体来说,我国暴雨主要集中在南方和东部沿海地区,沙尘暴易发生在西北内陆地区。机场吞吐量以 2020 年各机场数据为准;机场等级根据我国民用机场等级划分方法18 进行划分,各机场所在地区 GDP 通过国家统计局19 收集了 2020
18、年各地 GDP。于是,运用前面所提机场灾害风险度公式,可分别计算得到各机场在暴雨和沙尘暴灾害下的风险度。2.2仿真试验根据 1.3 节确定的灾害情景,即随机灾害、节点加权度中心性优先的蓄意灾害、节点加权介数中心性优先的蓄意灾害、暴雨灾害和沙尘暴灾害,分别进行网络抗打击试验。具体设计以 10%节点打击为步长,每个步长进行 100 次打击试验,并计算网络效率和最大连通子图相对大小的变化平均值。具体仿真试验结果见图 2 和 3。图 2 为 3 家航空公司航线网络节点遭受随机灾害、蓄意灾害和自然灾害后,网络效率 E 的变化情况。总体来看,无论哪种灾害情景,3 家航空公司航线网络的网络效率均随节点失效比
19、例增加而减小,直至为 0。蓄意灾害、自然灾害和随机灾害下,失效节点数分别在20%、70%和80%左右时 E 趋于0,表明不同灾害情景对航线网络的破坏性不同,即蓄意灾害对网络的破坏力大于自然灾害大于随机灾害。这是因为由少数加权度或介数较大的机场构成的子网络是整体航线网络的重要组成部分,在航线网络中具有枢纽性地位,对航线网络的整体连通性起着重要作用;而机场的灾害风险度与其地理位置和灾害的空间分布有关,灾害风险度较大的机场并不一定在网络拓扑结构的中心连接位置,对网络整体连通性的影响相对较小;随机灾害情景下,节点失效具有随机性,并不是每次都攻击网络中最为重要的节点,因此对网络影响较小。进一步来看,当
20、3 家航空公司航线网络遭受不同类型的蓄意灾害时,10%的节点失效时网络效率E 急剧下降,且通过网络遭受灾害前后网络效率 E的下降幅度可以看出加权介数中心性值较大的节点失效对航线网络性能影响更大。当 3 家航空公司航线网络遭受自然灾害时,网络效率下降较缓慢,当节点失效数在 70%左右时,网络效率趋近于 0,网络基本瘫痪;相比于沙尘暴灾害,暴雨对航线网络的影响更大。原因是,我国沙尘暴主要集中在西北和华北地区,暴雨多发生在我国东南沿海地区,而各航空公司航线均具有“东密西疏”的特点,因此暴雨灾害对各航空公司航线网络的破坏性大于沙尘暴。就航空公司而言,蓄意灾害下,相比于 M,A 航和 Z 航的网络效率
21、E 下降更快;Z 航线网络节点失效数达 20%时 E 已趋于 0,而 A 和 M 航线网络的节点失效数达 30%时 E 才趋于 0,此时整个网络中几乎只有孤立节点,即网络瘫痪。这表明相比于 A 和M,蓄意灾害对 Z 航线网络的破坏性更大。暴雨和沙尘暴两种自然灾害下,3 家航空公司航线网络在第二次攻击前后网络效率均出现了差值大于 0.1 的情况,其中 Z 的网络效率差值最大,A 的最小。而随机灾害情景下,当节点失效数在 80%左右时,3 家航空公司航线网络的 E 趋近于 0,网络基本瘫痪;这表明随机灾害情景下,3 家航空公司航线网络的抗打击能力较强。图 2航空公司航线网络效率与节点删除比例的关系
22、Fig 2elationship between route network efficiencyand node deletion ratio of airlines图 3 给出了 3 家航空公司航线网络分别遭受随机灾害、蓄意灾害和自然灾害后,最大连通子图相对6222Vol 23No 7安全 与 环 境 学 报第 23 卷第 7 期大小 S 的变化情况。总体来看,3 家航空公司航线网络初始状态下均没有孤立节点,因此初始状态下最大连通子图相对大小 S 为 1,随着节点失效比例的增加,最大连通子图的规模不断减小。且通过对比网络瘫痪时的节点失效数可以看出无论哪个航空公司,蓄意灾害下最大连通子图相对
23、大小 S 的下降幅度均大于自然灾害大于随机灾害。因为少数度或介数较大的节点大多是网络中的枢纽节点,这些节点一旦失效,与其有连接的节点都会受到干扰甚至破坏,从而导致网络中孤立节点增多,网络规模急剧减小。图 3航空公司航线网络最大连通子图相对大小与节点删除比例的关系Fig 3elationship between relative size of themaximum connected subgraph and nodedeletion ratio in the airline network进一步来看,蓄意灾害下,3 家航空公司航线网络的最大连通子图相对大小 S 均快速减小,且节点加权介数中心
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 多情 航空公司 航线 网络 抗毁性 研究 小兵
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。