Vulnerability Analysis of En-route Network Based on Cascading Failure
-
摘要:
为探究级联失效所造成的航路网络脆弱性变化规律,基于不同攻击方式的级联失效过程对航路网络脆弱性进行分析. 首先,结合航路网络超容运行实际,提取非失效状态节点的过载运行、节点失效的概率性和节点所经负载的可分流性3个特点,构建航路网络级联失效模型和不同攻击方式的级联失效过程;然后,从航路点运行能力损失角度提出航路网络级脆弱性指标以及与级联失效模型相结合的脆弱性分析方法;最后,以民航华东区域为实例,分析航路网络脆弱性与所建模型参数的关联性,探究不同参数条件下航路网络脆弱性指标的变化规律,并设计了3种攻击实验. 研究表明:各航路点运行能力随着负载流量的超容范围增加而提升,航路网络脆弱性随之降低;航路网络对选择性攻击方式较为敏感,尤其是基于介数的攻击.
Abstract:In order to explore the variation of the vulnerability of the en-route network caused by cascading failure, the vulnerability of the en-route network was analyzed based on the cascading failure process of different attack methods. By combining with the actual overcapacity operation of the en-route network, the overload operation of the nodes in the non-failure state, the probability of node failure, and the divertability of the load passed by the nodes were extracted. The cascading failure model of the en-route network and the cascading failure process based on different attack methods were constructed, and the en-route network-level vulnerability index was proposed from the perspective of the loss of waypoint operation capacity, as well as the vulnerability analysis method combined with the cascading failure model. The correlation between the vulnerability of the en-route network and the proposed model parameters was analyzed by studying the operation data in Eastern China, and the change rule of the vulnerability index of the en-route network under different parameters was explored. Three kinds of attack experiments were designed, and the results show that the operation capacity of each waypoint increases with the overcapacity range of the load flow, and the vulnerability of the en-route network decreases; the en-route network is more sensitive to selective attack methods, especially those based on the betweenness.
-
Key words:
- air traffic /
- en-route network /
- complex networks /
- cascading failure /
- vulnerability
-
柴油机有热效率高、功率范围宽和可靠性高等优点,作为动力源广泛应用于各类交通工具. 然而,随着环境污染和化石资源的不断消耗,世界各国对柴油机的经济性和排放性提出了越来越高的要求[1]. 从喷孔喷出油束的喷油特性直接影响柴油和空气的混合特性,特别是喷孔的喷油稳定性会直接影响每次喷射的循环喷油量以及油束喷雾在缸内的分布,从而对柴油机缸内燃烧过程产生重要的影响 [2]. 因此,喷油器的喷油规律以及喷油的稳定性对柴油机进行缸内喷雾、燃烧性能优化具有重要的意义[3].
学者们针对喷油器喷孔喷油速率开展了大量研究工作,取得了丰富的研究成果[4-6]. Suh等[7]采用Bosch喷油规律测量方法,开展了电磁驱动式和压电驱动式喷油器喷油规律实验研究,实验发现,在同一工况条件下,压电式喷油器的循环喷油量高于电磁式喷油器. Payri等[8]基于Bosch法开展了柴油和汽油为燃料的喷油速率实验研究发现,在相同的喷射条件下,柴油的喷射速率高于汽油的,这是由于柴油和汽油密度不同造成的. Postrioti等[9]采用Zeuch和Bosch 2种方法进行喷油规律测量结果准确性的比较,结果表明,Bosch法测量的喷油开始时刻比Zeuch法测量的喷油始点延迟50~70 μs. Cavicchi等[10]采用Zeuch法结合激光多普勒测振仪开展了喷油速率与针阀升程的实验研究发现,在喷油器电流信号加载后第70 μs时才测量到针阀的位移,且在针阀抬起的瞬间喷油速率出现一个下降的趋势,这是因为在针阀抬起后,喷孔内的燃油出现了短暂的回流,在针阀开启8~10 μs后喷油速率才得以上升. Desantes等[11]基于动量法开展了生物质-柴油的混合油喷油规律实验,研究发现,在不同的工况下,喷油器喷孔的最大喷油速率随生物质掺混比的上升而增大. Luo等[12]基于动量法对多喷孔喷油器的喷油规律进行了同步测量,对于直径相同的喷孔,喷孔的水平夹角越小,喷油起始时间会越晚,且喷油量会越小. Faghani等[13]采用动量法对天然气直喷喷射器的稳定性进行研究,随着喷射脉宽的增加,相同工况下喷射的稳定性会上升,喷气规律曲线间的差异性会降低.
和传统的喷油规律测量方法相比,动量法是唯一可以测量多喷孔喷油器各孔喷油速率的实验方法,但较少有研究涉及到利用动量法对多孔喷油器在不同喷油压力下的稳定性. 因此,本文首先根据动量法开发了多孔喷油器各孔喷油规律测量装置;然后,基于研发的测量装置,在不同的负荷工况下,测量喷油器各喷孔喷油量的差异;在不同的喷油压力下,针对单一喷孔进行了喷孔喷油特性的研究,揭示了喷油压力对喷油量稳定性的影响机制.
1. 动量法测量原理和喷油稳定性评价
1.1 基于动量法的喷油规律测量原理
图1所示为动量法测量喷油规律原理示意. 图中:u为燃油的流速,˙m为喷孔内燃油瞬态质量, F为冲击力. 采用动量法时,在喷孔出口处一定距离上布置力传感器,力传感器表面与喷孔表面平行,此时假设油束完全撞击到传感器表面.
根据质量守恒定理,喷孔内燃油瞬态质量流量为
˙m=ρuA, (1) 式中:ρ为燃油密度,A为喷孔的流通面积.
通过传感器采集到瞬态冲击力,根据动量守恒定理可知,冲击力F为
F=˙mu. (2) 联立式(1)、(2),得到
˙m=√ρAF. (3) 1.2 喷油量稳定性评价指标
对于共轨多孔喷油器而言,喷孔与喷孔之间由于加工误差和流动的不一致性会导致喷油器各喷孔间燃油喷油量的差异. 此外,对于同一个喷孔,不同循环次数下,每次喷射的喷油量也存在不一致性现象. 因此,本文定义喷油量波动率σ来评价喷孔喷油量的稳定性,如式(4).
σ=|mi−mave|mave, (4) 式中:mi为第i个喷孔的喷油量,mave为总的喷油量的平均值.
2. 动量法测量系统设计
2.1 力传感器参数确定
据伯努利方程和动量定理,可计算出在最大压差下喷孔油束喷雾的理论最大冲击力为
Fmax=ρu2maxA=2(P1−P2)A, (5) 式中:P1为喷油压力,P2为喷孔出口背压,umax为油束在喷孔处的流速.
在本文的实验条件下,P1=120 MPa,P2=0.1 MPa,计算得到Fmax=17.03 N. 由于喷油持续时间为毫秒级,瞬态力传感器需要具有高的响应频率. 考虑到测量范围和响应频率要求,选择了大洋公司DYZ-100型力传感器作为本文实验用的瞬态力传感器,具体参数如表1所示.
表 1 力传感器主要性能指标Table 1. Main performance indicators of force sensor指标 数值 量程/N 0~20 输出精度/% 0.05 响应频率/kHz 10 工作温度/℃ −20~80 2.2 喷油规律计算程序
在进行喷油规律计算时,首先通过数据实时采集子程序进行油束力信号采集,为保证对信号的完整采集,采样率设置为100 KS/s,同时根据2 Hz的喷油周期设置0.5 s的采样时间,并采用循环结构进行信号的连续采集;完成信号采集后,通过滤波子程序对采集到的数据进行滤波,以减少干扰杂波,在进行滤波时,根据FFT(Fast Fourier Transform)频谱分析结果发现,试验中主要的波形集中在0~
1500 Hz频率段,因此,选择1500 Hz低通滤波以抑制高频噪声;最后,根据式(3)编程得到瞬态质量流量,在喷油持续期内进行积分处理,得到循环喷油量.3. 装置测量结果准确性验证
图2所示为本文建立的基于动量法的喷油器喷孔喷油规律测量装置. 该测量装置包括高压共轨喷油系统、高压共轨喷油系统控制台、共轨喷油器、实验柜、瞬态力传感器、油雾过滤器、数据实时采集系统、实验柜.
图3为法国EFS公司的EMI 2型单次喷射仪,在进行喷油规律测量时,首先将喷油器安装在单次喷射测量仪机械项部,调整背压关闭排油阀,然后向底部安装有可上下移动活塞的密封盛油腔内喷射燃油,从而导致控制活塞产生位移. 通过位移传感器将位移以电压信号输出,在对电压信号进行处理后,便可以得到相应位移下的循环喷油量以及相应的喷油规律曲线. EMI 2型单次喷射仪是EFS公司较为成熟的商业化测试仪器,测量误差小于满量程的± 0.1%,广泛应用在喷油规律测量的研究中 [14-15].
图4为所开发的实验装置在喷油压力100 MPa和喷油脉宽1.5 ms的条件下重复进行20次喷射的喷油规律. 由图4可知:每次喷射过程中喷油规律曲线的变化趋势均相同,并且不同喷射次数下的波谷、波峰出现的时间也一致;从循环喷油量来看,第9次喷射的循环喷油量的最大为131.23 mg,最小的循环喷油量出现在第16次为122.87 mg,循环喷油量的最大测量误差为6.8%,本文所开发的实验装置具有较好的复现性.
图5所示为喷油压力100 MPa和喷油脉宽1.5 ms,本文开发的测量装置和EFS公司的EMI 2型商业单次喷射仪在重复进行20次试验下测量结果均值的对比. 从图5看到:在喷油速率上升阶段,本文测量装置的测量结果和EMI 2测量结果的一致性较好,但在最大喷油速率上存在一定的差异,本文测量的最大值为80.96 g/s,而EIM2测量值为88.34 g/s,两者之间的偏差为8.4%,满足不超过10%的工程误差的要求[16];在喷油速率下降阶段,本文测量的喷油速率下降较EMI 2更缓,喷油结束时刻测量值晚于EMI 2测量值. Payri等[17]在采用动量法进行喷油规律测量时同样出现了上述喷油速率下降段延迟的现象,在采用动量法进行喷油规律测量时,喷孔和传感器之间存在一定的距离,在针阀关闭后会有一部分油束需要一定的时间才能到达传感器表面,造成了基于动量法测量出现的后端喷油规律变缓的现象,但这对喷油量几乎没有影响.
实际上,循环喷油量是保证柴油机所需动力性能的重要参数. 从循环喷油量的测量结果来看,本文测量值为126.38 mg,EMI2测量值为126.53 mg,两者的相对测量偏差仅为0.12%. 因此,从喷油规律曲线和循环喷油量可以看到,本文研发的基于动量法的喷孔喷油规律测量装置具有较高的测量准确性,可以采用该测量装置进行后续喷油特性测量研究.
4. 结果与讨论
4.1 各喷孔喷油差异性实验研究
图6为本文采用的Bosch多孔喷油器结构示意,在喷嘴上均匀分布4个直径为0.3 mm,水平夹角为20° 的喷孔,针阀位于喷油器的中心,通过电磁阀控制针阀的起落改变燃油的流通面积.
喷油压力60 MPa,喷油脉宽0.5 ms时各喷孔喷油速率实验结果对比如图7. 从图7观察到:喷孔间喷油速率差异较大,孔4喷油速率的峰值最小,孔2喷油速率的峰值最大,而孔1 (1.84 mg)和孔3(1.91 mg)喷油速率曲线几乎重合在一起;从循环喷油量角度分析发现,孔4喷油量为1.40 mg,是循环喷油量最小的,而孔2的循环喷油量为2.18 mg,为最大的,孔4和孔2循环喷油量的测量偏差达到55.71 %. 孔4和孔2喷油规律曲线显著差异的原因可能是针阀在抬升和落座过程中发生了向孔4方面的偏心,从而使得孔4实际的喷孔流通面积减小.
4.2 低喷油压力下喷油稳定性实验研究
通过3.1部分研究发现,喷油压力和喷油脉宽共同作用会影响喷孔喷油的稳定性. 图8 为低喷油压力下(60 MPa)喷油速率随喷油脉宽的变化. 从图8看到:在脉宽为0.5 ms时,由于针阀没有完全打开而出现小喷油量;随着喷油脉宽从0.5 ms增加到2.0 ms过程中,喷油速率最大值不断增大,喷油速率曲线出现波动,这是由于针阀运动导致了高压共轨喷油器内产生了明显的压力波动,而压力波动又反过来影响喷油速率曲线的波动.
为了更清楚获得不同脉宽下喷孔喷油量波动率,针对孔3,采用式(4)计算50次喷射实验下喷油量波动率,如图9. 从图9看到:随着喷油脉宽的增加,喷孔喷油量波动率从15.88%降低9.81%,降幅为37.97%. 这个现象的原因在于,对于60 MPa小喷油压力,本文共轨喷油器理论上可实现的最大喷油速率为15.41 g/s,在喷油脉宽<2 ms条件下,针阀并没有完全打开,此时,针阀抬起后便落座,未达到最大限位,使得针阀及针阀座之间的燃油流动成为主导,针阀运动不稳定引起了较大的喷油量波动率;随着针阀完全打开,针阀运动的不稳定性降低,喷油量波动率随着喷油脉宽的增加而减小.
4.3 高喷油压力下喷油稳定性实验研究
图10为高喷油压力下(120 MPa)喷油速率随喷油脉宽的变化. 从图10看到:当喷油脉宽为0.5 ms时,针阀并没有达到最大升程位置,而当喷油脉宽增加到1.0 ms时,喷油速率出现了稳定段,并且其最大喷油速率为22.61 g/s,和理论值22.23 g/s几乎一致,从而判断出喷油脉宽1.0 ms下针阀完全打开;随着喷油脉宽继续增加到2.0 ms,喷油速率中间平缓段变长,尽管存在一定的波动,主要是由于喷油诱发的压力波动导致.
与图9所示的低喷油压力下脉宽对喷油量波动率的影响相比,从图11看到:在更高的喷油压力下,不同脉宽的喷油量波动率均降低,且随着脉宽的增加,喷油量波动率降低变缓. 这是因为,随着喷油压力的上升,相同的脉宽下针阀更容易达到最大升程. 结合图10可以观察到:在120 MPa喷油压力下,脉宽大于1.0 ms时,喷油速率出现了稳定段,对应的喷油波动率也降低;由于压力波动的影响,喷油量的循环波动率也产生波动,因此才出现了在脉宽为2.0 ms时循环波动率增大的现象, 但仅增加了1%. 再次说明,针阀无法实现最大升程是高压共轨燃油喷射系统喷射不稳定性的主要原因.
5. 结 论
1) 基于动量法开发了共轨喷油器喷孔喷油规律测量装置,与商业单次喷射仪测量结果进行了对比,在最高喷油速率方面,两者测量偏差小于10%;在循环喷油量方面,两者测量偏差为0.12%,证明了本文研发的测量装置测量的准确性.
2) 在低喷油压力条件下,随着喷油脉宽从0.5 ms增加到2.0 ms,喷孔喷油量波动率随之降低,在上述低喷油压力和脉宽条件下,针阀无法完全打开,针阀—针阀座之间燃油的不稳定流动成为影响喷油量波动的主要因素.
3) 更高的喷油压力下,喷油脉宽相同时,针阀更容易达到最大升程,对喷孔喷油量波动起到决定作用的是喷孔参数. 因此,在0.5~2.0 ms喷油脉宽范围内喷孔喷油量波动率处于5%之内. 这说明,针阀无法实现最大升程是高压共轨系统燃油喷射不稳定性的主要原因.
-
表 1 航路网络节点度分布情况
Table 1. Distribution of node degree of en-route network
度数 节点数/个 比例/% 1 24 7.0 2 151 44.0 3 67 19.5 4 65 19.0 5 16 4.7 6 13 3.8 7 3 0.9 8 1 0.3 9 2 0.6 10 1 0.3 -
[1] 程光权,陆永中,张明星,等. 复杂网络节点重要度评估及网络脆弱性分析[J]. 国防科技大学学报,2017,39(1): 120-127.CHENG Guangquan, LU Yongzhong, ZHANG Mingxing, et al. Node importance evaluation and network vulnerability analysis on complex network[J]. Journal of National University of Defense Technology, 2017, 39(1): 120-127. [2] HU Y Q, ZHOU D, ZHANG R, et al. Percolation of interdependent networks with intersimilarity[J]. Physical Review E, 2013, 88(5): 052805.1-052805.7. [3] WANG J W, SUN E H, XU B, et al. Abnormal cascading failure spreading on complex networks[J]. Chaos, Solitons & Fractals, 2016, 91: 695-701. [4] KIM M, KIM J S. A model for cascading failures with the probability of failure described as a logistic function[J]. Scientific Reports, 2022, 12: 989-998. doi: 10.1038/s41598-021-04753-z [5] 徐伟举. 基于复杂网络的美国航空线路网络的抗毁性研究[D]. 成都:西南交通大学,2013. [6] 王兴隆,贺敏,刘明学. 空中交通CPS级联失效与缓解策略[J]. 北京航空航天大学学报,2021,47(12): 2426-2433.WANG Xinglong, HE Min, LIU Mingxue. Air traffic CPS cascading failure and mitigation strategy[J]. Journal of Beijing University of Aeronautics and Astronautics, 2021, 47(12): 2426-2433. [7] 汪凯强,王布宏,曾乐雅,等. 不同负载条件下的航空网络抗毁性研究[J]. 电光与控制,2022,29(7): 86-90.WANG Kaiqiang, WANG Buhong, ZENG Leya, et al. Survivability of airline network under different load models[J]. Electronics Optics & Control, 2022, 29(7): 86-90. [8] MIRZASOLEIMAN B, BABAEI M, JALILI M, et al. Cascaded failures in weighted networks[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 2011, 84: 046114.1-046114.8. [9] CUMELLES J, LORDAN O, SALLAN J M. Cascading failures in airport networks[J]. Journal of Air Transport Management, 2021, 92: 102026.1-102026.10. [10] 齐雁楠,高经东. 空域扇区网络级联失效抗毁性及优化策略[J]. 航空学报,2018,39(12): 349-357.QI Yannan, GAO Jingdong. Cascading failure invulnerability and optimization strategy of airspace sector network[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(12): 349-357. [11] 王兴隆,刘明学,高经东. 基于复杂网络的空域扇区网络建模与抗毁性研究[J]. 武汉理工大学学报,2018,40(11): 50-55.WANG Xinglong, LIU Mingxue, GAO Jingdong. Research on modeling analysis and invulnerability of sector network based on complex network[J]. Journal of Wuhan University of Technology, 2018, 40(11): 50-55. [12] 王兴隆,潘维煌,赵末. 华北航空相依网络统计特征与鲁棒性分析[J]. 科学技术与工程,2018,18(13): 180-185.WANG Xinglong, PAN Weihuang, ZHAO Mo. Analysis of network characteristics and robustness of aeronautical interdependent network in North China[J]. Science Technology and Engineering, 2018, 18(13): 180-185. [13] 王兴隆,潘维煌. 航空相依网络的弹性测度与分析[J]. 计算机仿真,2019,36(11): 40-43,58.WANG Xinglong, PAN Weihuang. The resilience measurement and analysis of aeronautical interdependent network[J]. Computer Simulation, 2019, 36(11): 40-43,58. [14] 王兴隆,潘维煌,赵末. 航空相依网络故障传播机理研究[J]. 武汉理工大学学报,2018,40(6): 41-46.WANG Xinglong, PAN Weihuang, ZHAO Mo. Research on fault propagation mechanism of aviation interdependent network[J]. Journal of Wuhan University of Technology, 2018, 40(6): 41-46. [15] ZHANG J, ZHANG Z H. A review of the research methods on vulnerability of transportation system[C]// CICTP 2019. Nanjing: American Society of Civil Engineers, 2019: 4294-4305. [16] GUZE S. Graph theory approach to the vulnerability of transportation networks[J]. Algorithms, 2019, 12(12): 270.1-270.12. [17] 王兴隆,潘维煌,赵末. 空中交通相依网络的脆弱性研究[J]. 航空学报,2018,39(12): 268-277.WANG Xinglong, PAN Weihuang, ZHAO Mo. Vulnerability of air traffic interdependent network[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(12): 268-277. [18] 周子丰. 基于复杂网络理论的民航陆基空管保障系统脆弱性分析[D]. 广汉:中国民用航空飞行学院,2018. [19] WILKINSON S M, DUNN S, MA S. The vulnerability of the European air traffic network to spatial hazards[J]. Natural Hazards, 2012, 60(3): 1027-1036. doi: 10.1007/s11069-011-9885-6 [20] 李航,胡小兵. 一种改进的民用航空网络空间脆弱性模型[J]. 交通运输系统工程与信息,2018,18(4): 202-208.LI Hang, HU Xiaobing. An improved spatial vulnerability model for civil aviation network[J]. Journal of Transportation Systems Engineering and Information Technology, 2018, 18(4): 202-208. [21] JANIĆ M. Modelling the resilience, friability and costs of an air transport network affected by a large-scale disruptive event[J]. Transportation Research Part A: Policy and Practice, 2015, 71: 1-16. doi: 10.1016/j.tra.2014.10.023 [22] VOLTES-DORTA A, RODRÍGUEZ-DÉNIZ H, SUAU-SANCHEZ P. Vulnerability of the European air transport network to major airport closures from the perspective of passenger delays: ranking the most critical airports[J]. Transportation Research Part A: Policy and Practice, 2017, 96: 119-145. doi: 10.1016/j.tra.2016.12.009 [23] LORDAN O, KLOPHAUS R. Measuring the vulnerability of global airline alliances to member exits[J]. Transportation Research Procedia, 2017, 25: 7-16. doi: 10.1016/j.trpro.2017.05.189 [24] Roy S, Xue M, Sridhar B. Vulnerability metrics for the airspace system[C]//Twelfth USA/Europe Air Traffic Management Research and Development Seminar (ATM2017). Seattle: Curran Associates, 2017: 1-26. [25] TAMIMI A, HAHN A, ROY S. Cyber threat impact analysis to air traffic flows through dynamic queue networks[J]. ACM Transactions on Cyber-Physical Systems, 2020, 4(3): 1-22. [26] 田文,方琴,周雪芳,宋津津. 航路网络关键节点识别方法研究[J]. 西南交通大学学报,2025,60(1): 233-242.TIAN Wen, ZHOU Xuefang, FANG Qin, et al. Research on identification method of key nodes in en-route network[J]. Journal of Southwest Jiaotong University,2025,60(1): 233-242. [27] 种鹏云. 基于复杂网络的危险品运输网络拓扑特性、级联失效机制及抗毁性研究[D]. 成都:西南交通大学,2015. [28] 欧阳敏,费奇,余明辉,等. 复杂网络的功效性与脆弱性研究综述[J]. 计算机科学,2008,35(6): 1-4.OUYANG Min, FEI Qi, YU Minghui, et al. Survey on efficiency and vulnerability of complex network[J]. Computer Science, 2008, 35(6): 1-4. -