本发明公开了一种基于随机Petri网的不正常航班衍生事件链式效应预测方法，步骤为：对不正常航班事件中的不确定性因素进行分析；建立不确定性因素的形式化表达方法；构建不确定性因素扰动下不正常航班事件的结构模型；建立基于航班执行过程解析的不正常航班波及模型；建立不正常航班衍生事件链的层次结构模型；建立不正常航班衍生事件链式效应解析模型，建立不正常航班衍生事件链式演化系统；构建不正常航班衍生事件演化的随机事件链Petri网模型；基于马尔科夫链标识转移概率预测不正常航班衍生事件链式效应。本发明能够提升不正常航班的预测预警精度及其引发衍生事件应急处置和航班恢复的合理性和可靠性，具有广泛的推广应用前景。
1.一种基于随机Petri网的不正常航班衍生事件链式效应预测方法，其特征在于：包括以下步骤：步骤a，对不正常航班事件中的不确定性因素进行分析；具体步骤为：根据不确定性因素的影响范围，将不确定性因素分为飞机自身不确定性因素、机场及航空公司内部不确定性因素、航线网络不确定性因素；
步骤b，在步骤a的不确定性因素分析的基础上，对不确定性因素扰动下不正常航班事件进行共性分析与事件描述，建立不确定性因素的形式化表达方法；具体步骤为：将不正常航班事件的共性特征描述为不正常航班事件的作用过程导致航班正常运行秩序的破坏，采用危险供体、危险受体和映射方式对不正常航班事件进行定义：危险供体在映射方式的作用下，最终导致危险受体遭到破坏的过程；将不正常航班事件E抽象为危险供体D、危险受体R和映射方式M三个要素，进而将不正常航班事件E＝表示为：D＝(d1,d2,…,di)，R＝(r1,r2,…,rj)，M＝(m1,m2,…,mk)，其中，D表示危险供体di的集合，其中，i＝1,2,…，即各种不确定性因素；R表示危险受体rj的集合，其中，j＝1,2,…，即受到不确定性因素影响的航班、旅客；M表示映射方式mk的集合，其中，k＝1,2,…，即延误、取消、返航、碰撞、冲突；用相关物理量表达危险供体与危险受体之间的映射关系，其中，同一个映射方式能够用不同的参数进行描述，参数数值反映了不正常航班的映射程度，由导致航班不正常的不确定性因素和航班状态决定，不同的参数之间通过相互推导，能够确定不正常航班事件中危险供体和危险受体之间的映射规律；
步骤c，分析不正常航班事件的结构特征，构建不确定性因素扰动下不正常航班事件的结构模型；具体步骤为：首先，将不正常航班事件中危险供体属性集Φ、危险受体属性集Ψ和映射集Ξ表示为：
其中，Φi代表不正常航班事件危险供体的属性，Ψj代表不正常航班事件危险受体的属性，mk代表不正常航班事件的映射方式；Φ包含不正常航班事件危险供体的全部属性，Ψ包含不正常航班事件危险受体的全部属性，Ξ包含不正常航班事件的全部映射方式；
其次，对不正常航班事件E的危险供体集合D和危险受体集合R建立属性函数 和ψ
其中， 表示不正常航班事件E中危险供体D的属性，ψ(R)表示不正常航班事件E中危险受体R的属性；
危险供体D和危险受体R的属性决定了其是否能产生映射，进而为不正常航班事件的映射方式m构建属性映射判定函数γm(P,Q):其中，P为不正常航班事件中危险供体的属性集，Q为不正常航班事件中危险受体的属性集；γm＝1时说明包含属性集P的不确定性因素能够和包含属性集Q的航班产生映射m；γm＝0时说明不确定性因素和航班不能产生映射m；
相应的，对任意不正常航班事件的映射集M＝{mk|k＝1,2,…}构建属性映射判定函数γM(P,Q):其中，γM(P,Q)表示，若具有P属性集的危险供体D能够在映射集M中m的映射下与具有Q属性集的危险受体R发生相互作用，则称危险供体D与危险受体R产生映射M；
最终将不确定性因素扰动下的不正常航班事件结构模型描述为：
步骤d，考虑航班之间的关联性，将不确定性因素融入不正常航班波及过程中，根据不确定性因素的触发时机对航班执行过程进行细分，建立基于航班执行过程解析的不正常航班波及模型；具体步骤为：根据各种不确定因素对不正常航班带来的连锁反应，考虑不确定性因素的触发时机和不正常航班在时间轴上的波及规律，将航班执行过程细分为滑出阶段、爬升阶段、巡航阶段、下降阶段、滑入阶段和服务保障阶段，计算实际关舱门、起飞、着陆、开舱门、过站服务保障时间与计划时间的差值，精确定位航班延误的发生阶段，所述航班延误的发生阶段包括关舱门延误DTb、起飞延误DTd、滑出延误DTto、爬升延误DTcl、巡航延误DTcr、下降延误DTla、开舱门延误DTo、滑入延误DTti、过站服务保障时间TS，根据不正常航班演化过程，推导出如下不正常航班波及模型：N次不正常航班波及的数学形式：
由此可推导出整条航班链中最后一个航班的关舱门延误时间
步骤e，分析不正常航班执行过程中可能引发的衍生事件，根据衍生事件的演化模式，分析不正常航班衍生事件的链式波及效应，建立不正常航班衍生事件链的层次结构模型，描述不正常航班衍生事件之间的逻辑关系；具体步骤为：将不正常航班事件的波及描述为产生了不正常航班衍生事件，根据航班不正常的发生阶段，分析不正常航班衍生事件的表现形式及特征，抽象出不正常航班衍生事件的演化方式，即直链式演化、发散式演化和集中式演化和循环式演化；对不正常航班衍生事件进行层次结构划分，根据不正常航班事件E＝，其衍生事件为ES，将不正常航班可能引发的全部衍生事件集表示为：上式中，当i＝1时，说明先前事件只能引发一个衍生事件ES，先前事件E与衍生事件ES形式一种最基本的事件链单元，描述如下：L(E→ES)
当不正常航班衍生事件形成一种多层直链式结构时，则将上一层事件Ei与其衍生事件Ei+1的链式关系描述为：当不正常航班事件同时诱发多个衍生事件时，形成一种发散式演化结构，则将这种情况下的航班延误衍生事件链式关系结构描述为：在不正常航班衍生事件演化过程中，当多种不同的事件同时诱发了相同的衍生事件时，将源事件集表示为：OE＝{E|Ei→ES,i＝1,2,…,n}，将多个源事件导致相同衍生事件发生的链式关系结构描述为：L(Ei→ES),i＝2,…,n
步骤f，挖掘衍生事件的各种作用形式及作用过程，探究衍生事件发展演化机理，建立不正常航班衍生事件链式效应解析模型，推演不正常航班引发的各种衍生事件，建立不正常航班衍生事件链式演化系统；具体步骤为：将可能引发的各种不正常航班事件E，记作E1,E2,…,En，在不确定性因素的作用下，任意两个不正常航班事件或衍生事件Ei和Ej之间关系表示为：由上式建立能够描述存在于不正常航班事件中的衍生关系的布尔型矩阵：
根据布尔运算法则，将不正常航班事件的n级衍生关系矩阵表示为 其
中， a’ij＝(ai1∩a1j)∪(ai2∩a2j)∪…∪(aik∩akj),k＝1,2,…,n；a’ij代表不正常航班事件之间的二级衍生关系，进而根据逻辑和运算，得到最终的不正常航班事件衍生关系矩阵U：其中，衍生关系矩阵U中的元素uij代表航班延误事件Ei能否经过很多次的衍生后引发另一事件Ej；
针对不正常航班事件发生后可能引发的n个事件E1,E2,…,En，构建事件状态转移矩阵E＝[e1,e2,…,en]，矩阵中元素ei含义如下：不正常航班事件发生后，必然会相继引发各种衍生事件，则事件状态转移矩阵E＝[e1,e2,…,en]将随时间而变化；t时刻的事件状态转移矩阵由该时刻的不正常航班事件状态决定，结合可能发生的n个不正常航班事件的衍生关系矩阵A，推演出不正常航班事件中可能引发的各级衍生事件：其中，矩阵 中的每个元素 代表由t时刻的不正常航班事件推演出的k级衍生事件，其中，i＝1,2,…,n；同理得到任何时刻不正常航班事件的衍生趋势；
根据上述不正常航班衍生趋势的推演，梳理不正常航班及其波及导致的下游航班或下游机场航班不正常的一系列衍生事件链，形成一种多级、多层次的不正常航班链式关系，分析衍生事件之间相互耦合作用，构建不正常航班衍生事件链式演化系统；
步骤g，引入随机Petri网模型，利用随机Petri网离散事件动态表达和并行计算方法，构建不正常航班衍生事件演化的随机事件链Petri网模型；
步骤h，根据随机Petri网与马尔科夫链的同构关系，将不正常航班衍生事件演化的随机事件链Petri网模型转化为马尔科夫链，基于马尔科夫链标识转移概率预测不正常航班衍生事件链式效应。
2.根据权利要求1所述的基于随机Petri网的不正常航班衍生事件链式效应预测方法，其特征在于：所述步骤g的具体步骤为：根据不正常航班衍生事件链式演化过程，定义连续时间随机事件链Petri网SECPN＝(P,T,F,K,W,M,λ)，将库所P＝{P1,P2,…,Pn}表示为不正常航班衍生事件演化过程中的危险供体和危险受体；变迁T＝{t1,t2,…,tm}表示为危险供体和危险受体之间的映射方式，也即触发条件；有向弧集 I表示变迁输入弧的集合， O表示变迁输出弧的集合， F中允许有禁止弧，禁止弧仅存在与从库
所到变迁的弧；库所容量函数K＝1；弧函数W＝1；SECPN网标识M:P→N,为一向量，表示各库所中的托肯数量，M0表示系统初始状态，即初始标识；变迁平均实施速率集合λ＝{λ1,λ2,…,λm}，λt表示变迁tt∈T的平均实施速率；建立相应的不正常航班衍生事件演化的随机事件链Petri网模型。
3.根据权利要求2所述的基于随机Petri网的不正常航班衍生事件链式效应预测方法，其特征在于：所述不正常航班衍生事件演化的随机事件链Petri网模型由42个库所和44个变迁组成，各库所与变迁的具体含义表示为：库所的有限集P:
P1：飞机自身不确定性因素；P2：机场和航空公司内部不确定性因素；P3：航线网络中不确定性因素；P4：不正常航班事件；P5：谣言；P6：核心参与者；P7：情绪感染者；P8：旁观者；P9：旅客群体聚集；P10：旅客权利诉求；P11：旅客群体性暴力事件；P12：肢体冲突；P13：围堵登机门；P14：占机、罢机；P15：登机口封闭；P16：参与者被控制；P17：旅客情绪缓和；P18：航班堆积；
P19：机务人员工作负荷增加；P20：飞机机械故障；P21：飞机故障排除；P22：空中流量加剧；P23：管制员工作负荷增加；P24：空中流量管制；P25：空中流量管制解除；P26：机场设施设备故障；
P27：设施设备故障排除；P28：飞机损坏；P29：跑道冲突；P30：滑行道冲突；P31：停机位冲突；P32：航空器相撞；P33：异物入侵跑道道面；P34：异物入侵滑行道道面；P35：异物入侵停机位道面；
P36：跑道关闭；P37：滑行道关闭；P38：停机位关闭；P39：跑道正常使用；P40：滑行道正常使用；
t1：各种不确定性因素触发；t2：信息发布不及时；t3：信息传播；t4：旅客情绪相互感染；
t5：谣言扩散；t6：谣言肆意传播；t7：相关部门干预不力；t8：处置不及时或处置效果不佳；t9：旅客群体性事件恶化；t10：有效的情绪疏导和警察威慑；t11：疏散旅客，实时发布航班状态；
t12：不正常航班未及时处理；t13：航班保障资源不足；t14：空域容量不足；t15：机场设施设备供应不足；t16：航空器未能及时避让；t17：飞行区运行冲突；t18：航班保障能力下降；t19：飞机检修；t20：管制员指挥不正常航班运行；t21：未能及时疏散控制运行航班；t22：及时疏散空中运行航班；t23：设施设备检修；t24：设施设备故障未及时发现，并用于航班保障；t25：航空器撞击机场设施设备；t26：航空器之间相撞；t27：航空器搬移并检修；t28：在跑道上产生碰撞碎片；t29：在滑行道上产生碰撞碎片；t30：在停机位产生碰撞碎片；t31：跑道上碎片清除；t32：跑道上碎片未及时清除；t33：滑行道上碎片清除；t34：滑行道上碎片未及时清除；t35：停机位上碎片清除；t36：停机位上碎片未及时清除；t37：跑道冲突未有效处置；t38：航空器搬移，跑道冲突解除；t39：滑行道冲突未有效处置；t40：航空器搬移，滑行道冲突解除；t41：停机位冲突未有效处置；t42：航空器搬移，停机位冲突解除；t43：指挥航班正常运行；t44：善后处理。
4.根据权利要求1所述的基于随机Petri网的不正常航班衍生事件链式效应预测方法，其特征在于：所述步骤h的具体步骤为：通过分析基于随机事件链Petri网的不正常航班衍生事件演化模型，得出该模型的可实施与实施，即变迁激发规则：令SECPN＝(P,T,F,K,W,M,λ)为不正常航班衍生事件演化的随机事件链Petri网系统，M是(P,T,F)上的标识，t∈T；
(2)若t在M有发生权，则t可发生并将M改为M'，记作M[t>M',M'由下式给出：其中，*t＝{p|(p∈P∪T)∧((p,t)∈F)}称为t的前置集，t*＝{p|(p∈P∪T)∧((t,p)∈F)}称为t的后置集；
根据上述随机Petri网的变迁激发规则，对不正常航班衍生事件演化模型建立初始标识，得到不同变迁所得到的可达集；
在连续时间内，基于随机事件链Petri网的不正常航班衍生事件演化模型中的每一个变迁t从可实施到实施需要延时，对 t变成可实施的时刻到它实施时刻之间被看成一个连续随机变量xt，其中xt取正实数值，且服从与指数分布函数，即其中实参数λt>0是变迁t的平均实施速率，变量x≥0；
将不正常航班衍生事件演化的随机事件链Petri网模型同构于齐次马尔可夫链，其各库所转换的触发条件为时间随机变量，并服从于指数分布，变迁t1,t2,…,t44的平均实施速率分别为λ1,λ2,…,λ44；根据不正常航班衍生事件演化的随机事件链Petri网模型的状态标识M1,M2,…,Mn，得到其同构的马尔科夫链；马尔科夫链中有向弧表示不正常航班衍生事件演化过程的随机事件链Petri网模型从一个状态到另一个状态的转换；
通过对现实生活中不正常航班及其衍生事件的观测，统计出t1,t2,…,t44的变迁概率，即得到不正常航班及其衍生事件的马尔可夫链中各标识状态发生变迁的概率λ1,λ2,…,λ44，根据各状态标识Mi的发生序列ta,tb,tc,…，进行归一化处理后得到ta,tb,tc,…的各个变迁权重ωa,ωb,ωc,…及其从可实施到实施所需要的近似时间τa,τb,τc,…；
由于模型中τ服从于指数分布，Fi(τ)表示变迁ti发生所经历时间τ的概率分布，其概率分布函数为：若变迁只有一个可达集，由于变迁ti的时间为τ，则状态标识Mu变迁到Mv的概率为：
若状态标识有两个或两个以上变迁，在时间τ内，变迁ti发生，tj,…,tk不发生，则状态标识Mu变迁到Mv的概率为：对于一个给定的发生序列L＝{MuMi,…,MjMv}，其中u≠i≠j≠v,u,i,j,v∈L,系统由状态u到状态v的发生概率为：
其中，|L|＝n，n为变迁个数，h为系统由状态u到状态v所有可能的变迁个数；
由此，对不正常航班事件可能引发的衍生事件进行预测和概率计算。
基于随机Petri网的不正常航班衍生事件链式效应预测方法
技术领域
[0001] 本发明涉及信息系统技术领域，具体涉及一种基于随机Petri网的不正常航班衍生事件链式效应预测方法。
背景技术
[0002] 据不完全统计，近年来伴随民航运输量的逐年攀升，航班正常率水平却表现出下降趋势，在2011年，由于民航总局颁布了一连串的航班延误整治办法，使得航班正常率有所回升，但也只有77.9％，而随着航班量的进一步增长，2014年全国平均航班正常率已降落到68％。不正常航班治理一直是我国乃至世界民航发展中的热点问题。从民航运行实际来看，单架航空器依次执行多个航段而导致众航班之间存在天然的关联关系，加之航班计划编排十分紧凑，空域和机场容量日趋饱和，延误一旦产生，将以航班链作为传播媒介，向后续航班传播，航班延误的波及范围和延误程度在航班链中迅速扩大，此时极易发生大面积航班延误事件。此外，现有不正常航班恢复效率低下，特别是大面积航班延误事件极易引发一系列包括旅客群体性事件、跑道入侵、航空器或地面车辆相撞、航班保障不及时和机场保障资源紧张等在内的衍生事件，形成了一种多级、多层次的链式关系，即衍生事件的链式效应。
链式效应的产生加剧了不正常航班应急处置的不确定性，使航班恢复的难度成倍增加。若能对不正常航班引发的衍生事件进行有效预测，则可提前部署衍生事件的应急处置措施，尽可能减小衍生事件的影响程度和范围，提高不正常航班恢复效率。
[0003] 频繁发生的不正常航班事件逐渐上升为社会各界关心的主要问题。国内外学者的关注重点主要集中在不正常航班及其延误波及效应、不正常航班恢复方法等研究领域，而针对不正常航班引发衍生事件的链式效应则研究较少。然而，当不正常航班引发衍生事件时，通常会发现这些事件之间相互联结，相互影响，会发生一系列连锁反应，形成复杂的链式效应。目前关于不正常航班引发衍生事件的研究集中于不正常航班引发旅客群体性事件(简单的两层结构事件链)的因素分析和应对方法上。如刘小平基于突变论分析不正常航班导致旅客群体性事件的致因原理，认为航班不正常事件的外部因素和民航运输系统的内部环境共同作用并达到阈值时发生突变，引发旅客群体性事件；赵斌应用社会燃烧理论分析不正常航班引发的群体性事件的原因，认为民航法制建设滞后、服务不到位、信息不畅通是群体性事件的燃烧材料，谣言传播形成事件的助燃条件，延误后的群体冲突是不正常航班引发群体性事件的导火索；杨磊着重研究了不正常航班引发群体性事件的应急处置措施，提出了包括加强单位协调、提高航空公司人员责任意识、统一各航空公司补偿标准等应对方法。而国内外关于其他类型突发事件引发衍生事件链式效应的研究已成为研究热点，主要有：自然灾害触发技术灾害研究(Na-tech)、事故灾难的多米诺效应研究和灾害链与事件链理论研究。
发明内容
[0004] 发明目的：为了克服现有不正常航班预测方法的不足，本发明的目的是提供一种基于随机Petri网的不正常航班衍生事件链式效应预测方法，以增加不正常航班预测的全面性和准确性，提高不正常航班恢复的全局性。
[0005] 技术方案：为实现上述目的，本发明采用的技术方案是：
[0006] 一种基于随机Petri网的不正常航班衍生事件链式效应预测方法，包括以下步骤：
[0007] 步骤a，对不正常航班事件中的不确定性因素进行分析；
[0008] 步骤b，在步骤a的不确定性因素分析的基础上，对不确定性因素扰动下不正常航班事件进行共性分析与事件描述，建立不确定性因素的形式化表达方法；
[0009] 步骤c，分析不正常航班事件的结构特征，构建不确定性因素扰动下不正常航班事件的结构模型；
[0010] 步骤d，考虑航班之间的关联性，将不确定性因素融入不正常航班波及过程中，根据不确定性因素的触发时机对航班执行过程进行细分，建立基于航班执行过程解析的不正常航班波及模型；
[0011] 步骤e，分析不正常航班执行过程中可能引发的衍生事件，根据衍生事件的演化模式，分析不正常航班衍生事件的链式波及效应，建立不正常航班衍生事件链的层次结构模型，描述不正常航班衍生事件之间的逻辑关系；
[0012] 步骤f，挖掘衍生事件的各种作用形式及作用过程，探究衍生事件发展演化机理，建立不正常航班衍生事件链式效应解析模型，推演不正常航班引发的各种衍生事件，建立不正常航班衍生事件链式演化系统；
[0013] 步骤g，引入随机Petri网模型，利用随机Petri网离散事件动态表达和并行计算方法，构建不正常航班衍生事件演化的随机事件链Petri网模型；
[0014] 步骤h，根据随机Petri网与马尔科夫链的同构关系，将不正常航班衍生事件演化的随机事件链Petri网模型转化为马尔科夫链，基于马尔科夫链标识转移概率预测不正常航班衍生事件链式效应。
[0015] 进一步的，所述步骤a的具体步骤为：根据不确定性因素的影响范围，将不确定性因素分为飞机自身不确定性因素、机场及航空公司内部不确定性因素、航线网络不确定性因素。
[0016] 进一步的，所述步骤b的具体步骤为：将不正常航班事件的共性特征描述为不正常航班事件的作用过程导致航班正常运行秩序的破坏，采用危险供体、危险受体和映射方式对不正常航班事件进行定义：危险供体在映射方式的作用下，最终导致危险受体遭到破坏的过程；将不正常航班事件E抽象为危险供体D、危险受体R和映射方式M三个要素，进而将不正常航班事件E＝表示为：D＝(d1,d2,…,di)，R＝(r1,r2,…,rj)，M＝(m1,m2,…,mk)，其中，D表示危险供体di(i＝1,2,…)的集合，即各种不确定性因素；R表示危险受体rj(j＝1,2,…)的集合，即受到不确定性因素影响的航班、旅客；M表示映射方式mk(k＝1,2,…)的集合，即延误、取消、返航、碰撞、冲突；用相关物理量表达危险供体与危险受体之间的映射关系，其中，同一个映射方式能够用不同的参数进行描述，参数数值反映了不正常航班的映射程度，由导致航班不正常的不确定性因素和航班状态决定，不同的参数之间通过相互推导，能够确定不正常航班事件中危险供体和危险受体之间的映射规律。
[0017] 进一步的，所述步骤c的具体步骤为：
[0018] 首先，将不正常航班事件中危险供体属性集Φ、危险受体属性集Ψ和映射集Ξ表示为：
[0019] Φ＝{Φi|i＝1,2,…}
[0020] Ψ＝{Ψj|j＝1,2,…}
[0021] Ξ＝{mk|k＝1,2,…}
[0022] 其中，Φi代表不正常航班事件危险供体的属性，Ψj代表不正常航班事件危险受体的属性，mk代表不正常航班事件的映射方式；Φ包含不正常航班事件危险供体的全部属性，Ψ包含不正常航班事件危险受体的全部属性，Ξ包含不正常航班事件的全部映射方式；
[0023] 其次，对不正常航班事件E的危险供体集合D和危险受体集合R建立属性函数和ψ(X)，则不正常航班事件E＝有：
[0024]
[0025]
[0026]
[0027] 其中， 表示不正常航班事件E中危险供体D的属性，ψ(R)表示不正常航班事件E中危险受体R的属性；
[0028] 危险供体D和危险受体R的属性决定了其是否能产生映射，进而为不正常航班事件的映射方式m构建属性映射判定函数γm(P,Q):
[0029]
[0030]
[0031]
[0032] 其中，P为不正常航班事件中危险供体的属性集，Q为不正常航班事件中危险受体的属性集；γm＝1时说明包含属性集P的不确定性因素能够和包含属性集Q的航班产生映射m；γm＝1时说明不确定性因素和航班不能产生映射m；
[0033] 相应的，对任意不正常航班事件的映射集M＝{mk|k＝1,2,…}构建属性映射判定函数γM(P,Q):
[0034]
[0035]
[0036]
[0037] 其中，γM(P,Q)表示，若具有P属性集的危险供体D能够在映射集M中m的映射下与具有Q属性集的危险受体R发生相互作用，则称危险供体D与危险受体R产生映射M；
[0038] 最终将不确定性因素扰动下的不正常航班事件结构模型描述为：
[0039] E＝
[0040]
[0041] 进一步的，步骤d的具体步骤为：根据各种不确定因素对不正常航班带来的连锁反应，考虑不确定性因素的触发时机和不正常航班在时间轴上的波及规律，将航班执行过程细分为滑出阶段、爬升阶段、巡航阶段、下降阶段、滑入阶段和服务保障阶段，计算实际关舱门、起飞、着陆、开舱门、过站服务保障时间与计划时间的差值，精确定位航班延误的发生阶段，所述航班延误的发生阶段包括关舱门延误DTb、起飞延误DTd、滑出延误DTto、爬升延误DTcl、巡航延误DTcr、下降延误DTla、开舱门延误DTo、滑入延误DTti、过站服务保障时间TS，根据不正常航班演化过程，推导出如下不正常航班波及模型：
[0042] N次不正常航班波及的数学形式：
[0043]
[0044] 由此可推导出整条航班链中最后一个航班的关舱门延误时间
[0045]
[0046] 进一步的，步骤e的具体步骤为：将不正常航班事件的波及描述为产生了不正常航班衍生事件，根据航班不正常的发生阶段，分析不正常航班衍生事件的表现形式及特征，抽象出不正常航班衍生事件的演化方式，即直链式演化、发散式演化和集中式演化和循环式演化；对不正常航班衍生事件进行层次结构划分，根据不正常航班事件E＝，其衍生事件为ES，将不正常航班可能引发的全部衍生事件集表示为：
[0047]
[0048] 上式中，当i＝1时，说明先前事件只能引发一个衍生事件ES，先前事件E与衍生事件ES形式一种最基本的事件链单元，描述如下：
[0049] L(E→ES)
[0050] 当不正常航班衍生事件形成一种多层直链式结构时，则将上一层事件Ei与其衍生事件Ei+1的链式关系描述为：
[0051]
[0052] 当不正常航班事件同时诱发多个衍生事件时，形成一种发散式演化结构，则将这种情况下的航班延误衍生事件链式关系结构描述为：
[0053]
[0054]
[0055] 在不正常航班衍生事件演化过程中，当多种不同的事件同时诱发了相同的衍生事件时，将源事件集表示为：OE＝{E|Ei→ES,i＝1,2,…,n}，将多个源事件导致相同衍生事件发生的链式关系结构描述为：
[0056] L(Ei→ES),i＝2,…,n
[0057]
[0058] 进一步的，所述步骤f的具体步骤为：将可能引发的各种不正常航班事件E，记作E1,E2,…,En，在不确定性因素的作用下，任意两个不正常航班事件或衍生事件Ei和Ej之间关系表示为：
[0059]
[0060] 由上式建立能够描述存在于不正常航班事件中的衍生关系的布尔型矩阵：
[0061]
[0062] 根据布尔运算法则，将不正常航班事件的n级衍生关系矩阵表示为其中，AοA＝A2＝(a'ij)，a'ij＝(ai1∩a1j)∪(ai2∩a2j)∪…∪(aik∩akj),k＝1,2,…,n；a'ij代表不正常航班事件之间的二级衍生关系，进而根据逻辑和运算，得到最终的不正常航班事件衍生关系矩阵U：
[0063]
[0064]
[0065] 其中，衍生关系矩阵U中的元素uij代表航班延误事件Ei能否经过很多次的衍生后引发另一事件Ej；
[0066] 针对不正常航班事件发生后可能引发的n个事件E1,E2,…,En，构建事件状态转移矩阵E＝[e1,e2,…,en]，矩阵中元素ei含义如下：
[0067]
[0068] 不正常航班事件发生后，必然会相继引发各种衍生事件，则事件状态转移矩阵E＝[e1,e2,…,en]将随时间而变化；t时刻的事件状态转移矩阵由该时刻的不正常航班事件状态决定，结合可能发生的n个不正常航班事件的衍生关系矩阵A，推演出不正常航班事件中可能引发的各级衍生事件：
[0069]
[0070] 其中，矩阵 中的每个元素 代表由t时刻的不正常航班事件推演出的k级衍生事件；同理得到任何时刻不正常航班事件的衍生趋势；
[0071] 根据上述不正常航班衍生趋势的推演，梳理不正常航班及其波及导致的下游航班或下游机场航班不正常的一系列衍生事件链，形成一种多级、多层次的不正常航班链式关系，分析衍生事件之间相互耦合作用，构建不正常航班衍生事件链式演化系统。
[0072] 进一步的，所述步骤g的具体步骤为：根据不正常航班衍生事件链式演化过程，定义连续时间随机事件链Petri网SECPN＝(P,T,F,K,W,M,λ)，将库所P＝{P1,P2,…,Pn}表示为不正常航班衍生事件演化过程中的危险供体和危险受体；变迁T＝{t1,t2,…,tm}表示为危险供体和危险受体之间的映射方式，也即触发条件；有向弧集 I表示变迁输入弧的集合， O表示变迁输出弧的集合， F中允许有禁止弧，禁止弧仅存在与
从库所到变迁的弧；库所容量函数K＝1；弧函数W＝1；SECPN网标识M:P→N,为一向量，表示各库所中的托肯数量，M0表示系统初始状态，即初始标识；变迁平均实施速率集合λ＝{λ1,λ2,…,λm}，λt表示变迁tt∈T的平均实施速率；建立相应的不正常航班衍生事件演化的随机事件链Petri网模型。
[0073] 其中，所述不正常航班衍生事件演化的随机事件链Petri网模型由42个库所和44个变迁组成，各库所与变迁的具体含义表示为：
[0074] 库所的有限集P:
[0075] P1：飞机自身不确定性因素；P2：机场和航空公司内部不确定性因素；P3：航线网络中不确定性因素；P4：不正常航班事件；P5：谣言；P6：核心参与者；P7：情绪感染者；P8：旁观者；P9：旅客群体聚集；P10：旅客权利诉求；P11：旅客群体性暴力事件；P12：肢体冲突；P13：围堵登机门；P14：占机、罢机；P15：登机口封闭；P16：参与者被控制；P17：旅客情绪缓和；P18：航班堆积；P19：机务人员工作负荷增加；P20：飞机机械故障；P21：飞机故障排除；P22：空中流量加剧；
P23：管制员工作负荷增加；P24：空中流量管制；P25：空中流量管制解除；P26：机场设施设备故障；P27：设施设备故障排除；P28：飞机损坏；P29：跑道冲突；P30：滑行道冲突；P31：停机位冲突；
P32：航空器相撞；P33：异物入侵跑道道面；P34：异物入侵滑行道道面；P35：异物入侵停机位道面；P36：跑道关闭；P37：滑行道关闭；P38：停机位关闭；P39：跑道正常使用；P40：滑行道正常使用；P41：停机位正常使用；P42：事件平息；
[0076] 变迁的有限集T:
[0077] t1：各种不确定性因素触发；t2：信息发布不及时；t3：信息传播；t4：旅客情绪相互感染；t5：谣言扩散；t6：谣言肆意传播；t7：相关部门干预不力；t8：处置不及时或处置效果不佳；t9：旅客群体性事件恶化；t10：有效的情绪疏导和警察威慑；t11：疏散旅客，实时发布航班状态；t12：不正常航班未及时处理；t13：航班保障资源不足；t14：空域容量不足；t15：机场设施设备供应不足；t16：航空器未能及时避让；t17：飞行区运行冲突；t18：航班保障能力下降；t19：飞机检修；t20：管制员指挥不正常航班运行；t21：未能及时疏散控制运行航班；t22：及时疏散空中运行航班；t23：设施设备检修；t24：设施设备故障未及时发现，并用于航班保障；t25：航空器撞击机场设施设备；t26：航空器之间相撞；t27：航空器搬移并检修；t28：在跑道上产生碰撞碎片；t29：在滑行道上产生碰撞碎片；t30：在停机位产生碰撞碎片；t31：跑道上碎片清除；
t32：跑道上碎片未及时清除；t33：滑行道上碎片清除；t34：滑行道上碎片未及时清除；t35：停机位上碎片清除；t36：停机位上碎片未及时清除；t37：跑道冲突未有效处置；t38：航空器搬移，跑道冲突解除；t39：滑行道冲突未有效处置；t40：航空器搬移，滑行道冲突解除；t41：停机位冲突未有效处置；t42：航空器搬移，停机位冲突解除；t43：指挥航班正常运行；t44：善后处理。
[0078] 进一步的，所述步骤h的具体步骤为：通过分析基于随机事件链Petri网的不正常航班衍生事件演化模型，得出该模型的可实施与实施：
[0079] 令SECPN＝(P,T,F,K,W,M,λ)为不正常航班衍生事件演化的随机事件链Petri网系统，M是(P,T,F)上的标识，t∈T；
[0080] (1)t在M有发生权的条件是：
[0081]
[0082] t在M有发生权的事实记为M[t>。
[0083] (2)若t在M有发生权，则t可发生并将M改为M'，记作M[t>M',M'由下式给出：
[0084]
[0085] 其中，*t＝{p|(p∈P∪T)∧((p,t)∈F)}称为t的前置集，t*＝{p|(p∈P∪T)∧((t,p)∈F)}称为t的后置集；
[0086] 根据上述随机Petri网的变迁激发规则，对不正常航班衍生事件演化模型建立初始标识，得到不同变迁所得到的可达集；
[0087] 在连续时间内，基于随机事件链Petri网的不正常航班衍生事件演化模型中的每一个变迁t从可实施到实施需要延时，对 t变成可实施的时刻到它实施时刻之间被看成一个连续随机变量xt，其中xt取正实数值，且服从与指数分布函数，即
[0088]
[0089] 其中实参数λt>0是变迁t的平均实施速率，变量x≥0；
[0090] 将不正常航班衍生事件演化的随机事件链Petri网模型同构于齐次马尔可夫链，其各库所转换的触发条件为时间随机变量，并服从于指数分布，变迁t1,t2,…,t44的平均实施速率分别为λ1,λ2,…,λ44；根据不正常航班衍生事件演化的随机事件链Petri网模型的状态标识M1,M2,…,Mn，得到其同构的马尔科夫链；马尔科夫链中有向弧表示不正常航班衍生事件演化过程的随机事件链Petri网模型从一个状态到另一个状态的转换；
[0091] 通过对现实生活中不正常航班及其衍生事件的观测，统计出t1,t2,…,t44的变迁概率，即得到不正常航班及其衍生事件的马尔可夫链中各标识状态发生变迁的概率λ1,λ2,…,λ44，根据各状态标识Mi的发生序列ta,tb,tc,…，进行归一化处理后得到ta,tb,tc,…的各个变迁权重ωa,ωb,ωc,…及其从可实施到实施所需要的近似时间τa,τb,τc,…；
[0092] 由于模型中τ服从于指数分布，Fi(τ)表示变迁ti发生所经历时间τ的概率分布，其概率分布函数为：
[0093]
[0094] 若变迁只有一个可达集，由于变迁ti的时间为τ，则状态标识Mu变迁到Mv的概率为：
[0095]
[0096] 若状态标识有两个或两个以上变迁，在时间τ内，变迁ti发生，tj,…,tk不发生，则状态标识Mu变迁到Mv的概率为：
[0097]
[0098] 对于一个给定的发生序列L＝{MuMi,…,MjMv}，其中u≠i≠j≠v,u,i,j,v∈L,[0099] 系统由状态u到状态v的发生概率为：
[0100]
[0101] 其中，|L|＝n，n为变迁个数，h为系统由状态u到状态v所有可能的变迁个数；
[0102] 由此，对不正常航班事件可能引发的衍生事件进行预测和概率计算。
[0103] 有益效果：
[0104] 与现有的技术相比，本发明的有益效果是：基于不正常航班事件三要素理论建立不确定性因素扰动下不正常航班事件的结构模型，更加精确地描述和表征各类不确定性因素对不正常航班事件的影响作用；基于航班执行过程解析的不正常航班波及模型，准确刻画了不正常航班在不同时间段和不同节点上对相关航班的波及影响，提高了不正常航班波及模型的准确性和延误时间预测精度；从事件链角度建立不正常航班衍生事件层次结构模型和链式效应解析模型，反映了不正常航班事件间的相互作用机理，实现了不正常航班及其衍生事件之间逻辑关系的描述和可能引发的各种衍生事件的推演；利用随机Petri网和马尔科夫链的建模分析方法，建立了不正常航班衍生事件演化的随机事件链Petri网模型和马尔科夫链，充分考虑了多种衍生事件对不正常航班事件的反馈作用，更加全面地预估不正常航班衍生事件的发展态势和发生概率，提高了不正常航班事件后果的预测精度。
附图说明
[0105] 图1为实施方式中基于随机Petri网的不正常航班衍生事件链式效应预测方法的流程图；
[0106] 图2为实施方式中不正常航班事件的三要素关系图；
[0107] 图3为实施方式中不正常航班波及过程示意图；
[0108] 图4为实施方式中不正常航班衍生事件链式演化系统图；
[0109] 图5为实施方式中不正常航班衍生事件演化的随机事件链Petri网模型。
具体实施方式
[0110] 下面结合附图对本发明做更进一步的解释。
[0111] 如图1所示，本发明的一种基于随机Petri网的不正常航班衍生事件链式效应预测方法，包括以下步骤：
[0112] 第一步：不正常航班事件中不确定性因素分析；
[0113] 根据不确定性因素的影响范围，将不确定性因素分为飞机自身不确定性因素、机场及航空公司内部不确定性因素、航线网络不确定性因素。
[0114] 第二步：不确定性因素扰动下不正常航班事件共性分析与事件描述；
[0115] 将不正常航班事件的共性特征描述为不正常航班事件的作用过程导致航班正常运行秩序的破坏，采用危险供体、危险受体和映射方式对不正常航班事件进行定义：危险供体在映射方式的作用下，最终导致危险受体遭到破坏的过程。将不正常航班事件E抽象为危险供体D、危险受体R和映射方式M三个要素，进而将不正常航班事件E＝表示为：D＝(d1,d2,…,di)，R＝(r1,r2,…,rj)，M＝(m1,m2,…,mk)，如图2所示。其中，D表示危险供体di(i＝1,2,…)的集合，即各种不确定性因素；R表示危险受体rj(j＝1,2,…)的集合，即受到不确定性因素影响的航班、旅客等；M表示映射方式mk(k＝1,2,…)的集合，即延误、取消、返航、碰撞、冲突等。用相关物理量表达危险供体与危险受体之间的映射关系，如表1所示。
[0116] 表1不正常航班事件的映射方式表
[0117]
[0118] 其中，同一个映射方式可用不同的参数进行描述，参数数值反映了不正常航班的映射程度，由导致航班不正常的不确定性因素和航班状态决定，不同的参数之间通过相互推导，可确定不正常航班事件中危险供体和危险受体之间的映射规律。
[0119] 第三步：分析不正常航班事件的结构特征，构建不确定性因素扰动下不正常航班事件的结构模型；
[0120] 首先，将不正常航班事件中危险供体属性集Φ、危险受体属性集Ψ和映射集Ξ表示为：
[0121]
[0122] 其中，Φi代表不正常航班事件危险供体的属性，Ψj代表不正常航班事件危险受体的属性，mk代表不正常航班事件的映射方式。Φ包含不正常航班事件危险供体的全部属性，Ψ包含不正常航班事件危险受体的全部属性，Ξ包含不正常航班事件的全部映射方式。
[0123] 其次，对不正常航班事件E的危险供体集合D和危险受体集合R建立属性函数和ψ(X)，则不正常航班事件E＝有：
[0124]
[0125] 其中， 表示不正常航班事件E中危险供体D的属性，ψ(R)表示不正常航班事件E中危险受体R的属性。
[0126] 危险供体D和危险受体R的属性决定了其是否能产生映射，进而为不正常航班事件的映射方式m构建属性映射判定函数γm(P,Q):
[0127]
[0128] 其中，P为不正常航班事件中危险供体的属性集，其中危险供体为不确定性因素，Q为不正常航班事件中危险受体的属性集，其中危险受体为航班。γm＝1时说明包含属性集P的不确定性因素能够和包含属性集Q的航班产生映射m；γm＝1时说明不确定性因素和航班不能产生映射m。
[0129] 相应的，对任意不正常航班事件的映射集M＝{mk|k＝1,2,…}构建属性映射判定函数γM(P,Q):
[0130]
[0131] 其中，γM(P,Q)表示，若具有P属性集的危险供体D能够在映射集M中m的映射下与具有Q属性集的危险受体R发生相互作用，则称危险供体D与危险受体R产生映射M。
[0132] 最终将不确定性因素扰动下的不正常航班事件结构模型描述为：
[0133]
[0134] 第四步：根据不确定性因素的触发时机对航班执行过程进行细分，建立基于航班执行过程解析的不正常航班波及模型；
[0135] 根据各种不确定因素对不正常航班带来的连锁反应，考虑不确定性因素的触发时机和不正常航班在时间轴上的波及规律，将航班执行过程细分为滑出阶段、爬升阶段、巡航阶段、下降阶段、滑入阶段和服务保障阶段，如图3所示，计算实际关舱门、起飞、着陆、开舱门、过站服务保障时间与计划时间的差值，精确定位航班延误的发生阶段，航班延误的发生阶段包括关舱门延误DTb、起飞延误DTd、滑出延误DTto、爬升延误DTcl、巡航延误DTcr、下降延la o ti误DT 、开舱门延误DT 、滑入延误DT 、过站服务保障时间TS，根据不正常航班演化过程，推导出如下不正常航班波及模型：
[0136] N次不正常航班波及的数学形式：
[0137]
[0138] 由此可推导出整条航班链中最后一个航班的关舱门延误时间 其中，根据国内航班运行规定，一架飞机一天内执行的航班任务数量不得超过10个，因此N最大取值为10：
[0139]
[0140] 第五步：分析不正常航班执行过程中可能引发的衍生事件及链式波及效应，建立不正常航班衍生事件链的层次结构模型，描述不正常航班衍生事件之间的逻辑关系；
[0141] 将不正常航班事件的波及描述为产生了不正常航班衍生事件，根据航班不正常的发生阶段，分析不正常航班衍生事件的表现形式及特征，抽象出不正常航班衍生事件的演化方式，即直链式演化、发散式演化和集中式演化和循环式演化。对不正常航班衍生事件进S行层次结构划分，根据不正常航班事件E＝，其衍生事件为E ，将不正常航班可能引发的全部衍生事件集表示为：
[0142]
[0143] 上式中，当i＝1时，说明先前事件只能引发一个衍生事件ES，先前事件E与衍生事件ES形式一种最基本的事件链单元，描述如下：
[0144] L(E→ES) (9)
[0145] 当不正常航班衍生事件形成一种多层直链式结构时，则将上一层事件Ei与其衍生事件Ei+1的链式关系描述为：
[0146]
[0147] 当不正常航班事件同时诱发多个衍生事件时，形成一种发散式演化结构，则将这种情况下的航班延误衍生事件链式关系结构描述为：
[0148]
[0149]
[0150] 在不正常航班衍生事件演化过程中，当多种不同的事件同时诱发了相同的衍生事件时，将源事件集表示为：OE＝{E|Ei→ES,i＝1,2,…,n}，将多个源事件导致相同衍生事件发生的链式关系结构描述为：
[0151] L(Ei→ES),i＝2,…,n (13)
[0152]
[0153] 第六步：挖掘衍生事件的各种作用形式及作用过程，探究衍生事件发展演化机理，建立不正常航班衍生事件链式效应解析模型，推演不正常航班引发的各种衍生事件，建立不正常航班衍生事件链式演化系统；
[0154] 将可能引发的各种不正常航班事件E，记作E1,E2,…,En，在不确定性因素的作用下，任意两个不正常航班事件或衍生事件Ei和Ej之间关系表示为：
[0155]
[0156] 由上式建立能够描述存在于不正常航班事件中的衍生关系的布尔型矩阵：
[0157]
[0158] 根据布尔运算法则，将不正常航班事件的n级衍生关系矩阵表示为其中，AοA＝A2＝(a'ij)，a'ij＝(ai1∩a1j)∪(ai2∩a2j)∪…∪(aik∩akj),k＝1,2,…,n。a'ij代表不正常航班事件之间的二级衍生关系，进而根据逻辑和运算，得到最终的不正常航班事件衍生关系矩阵U：
[0159]
[0160]
[0161] 其中，衍生关系矩阵U中的元素uij代表航班延误事件Ei能否经过很多次的衍生后引发另一事件Ej。
[0162] 针对不正常航班事件发生后可能引发的n个事件E1,E2,…,En，构建事件状态转移矩阵E＝[e1,e2,…,en]，矩阵中元素ei含义如下：
[0163]
[0164] 不正常航班事件发生后，必然会相继引发各种衍生事件，则事件状态转移矩阵E＝[e1,e2,…,en]将随时间而变化。t时刻的事件状态转移矩阵由该时刻的不正常航班事件状态决定，结合可能发生的n个不正常航班事件的衍生关系矩阵A，可以推演出不正常航班事件中可能引发的各级衍生事件：
[0165]
[0166] 其中，矩阵 中的每个元素 代表由t时刻的不正常航班事件推演出的k级衍生事件。同理可得任何时刻不正常航班事件的衍生趋势。
[0167] 梳理不正常航班及其波及导致的下游航班或下游机场航班不正常等一系列衍生事件链，形成一种多级、多层次的不正常航班链式关系，分析衍生事件之间相互耦合作用，构建不正常航班衍生事件链式演化系统，如图4所示。
[0168] 第七步：引入随机Petri网模型，利用随机Petri网离散事件动态表达和并行计算方法，构建不正常航班衍生事件演化的随机事件链Petri网模型；
[0169] 根据不正常航班衍生事件链式演化过程，定义连续时间随机事件链Petri网SECPN＝(P,T,F,K,W,M,λ)，将库所P＝{P1,P2,…,Pn}表示为不正常航班衍生事件演化过程中的危险供体和危险受体；变迁T＝{t1,t2,…,tm}表示为危险供体和危险受体之间的映射方式，也即触发条件；有向弧集 I表示变迁输入弧的集合， O表示变迁输出弧的集合， F中允许有禁止弧，禁止弧仅存在与从库所到变迁的弧；库所容量函数K＝
1；弧函数W＝1；SECPN网标识M:P→N,为一向量，表示各库所中的托肯数量，M0表示系统初始状态，即初始标识；变迁平均实施速率集合λ＝{λ1,λ2,…,λm}，λt表示变迁tt∈T的平均实施速率。建立相应的不正常航班衍生事件演化的随机事件链Petri网模型，如图5。图5所示模型由42个库所和44个变迁组成，各库所与变迁的具体含义表示为：
[0170] 库所的有限集P:
[0171] P1：飞机自身不确定性因素；P2：机场和航空公司内部不确定性因素；P3：航线网络中不确定性因素；P4：不正常航班事件；P5：谣言；P6：核心参与者；P7：情绪感染者；P8：旁观者；P9：旅客群体聚集；P10：旅客权利诉求；P11：旅客群体性暴力事件；P12：肢体冲突；P13：围堵登机门；P14：占机、罢机；P15：登机口封闭；P16：参与者被控制；P17：旅客情绪缓和；P18：航班堆积；P19：机务人员工作负荷增加；P20：飞机机械故障；P21：飞机故障排除；P22：空中流量加剧；
P23：管制员工作负荷增加；P24：空中流量管制；P25：空中流量管制解除；P26：机场设施设备故障；P27：设施设备故障排除；P28：飞机损坏；P29：跑道冲突；P30：滑行道冲突；P31：停机位冲突；
P32：航空器相撞；P33：异物入侵跑道道面；P34：异物入侵滑行道道面；P35：异物入侵停机位道面；P36：跑道关闭；P37：滑行道关闭；P38：停机位关闭；P39：跑道正常使用；P40：滑行道正常使用；P41：停机位正常使用；P42：事件平息。
[0172] 变迁的有限集T:
[0173] t1：各种不确定性因素触发；t2：信息发布不及时；t3：信息传播；t4：旅客情绪相互感染；t5：谣言扩散；t6：谣言肆意传播；t7：相关部门干预不力；t8：处置不及时或处置效果不佳；t9：旅客群体性事件恶化；t10：有效的情绪疏导和警察威慑；t11：疏散旅客，实时发布航班状态；t12：不正常航班未及时处理；t13：航班保障资源不足；t14：空域容量不足；t15：机场设施设备供应不足；t16：航空器未能及时避让；t17：飞行区运行冲突；t18：航班保障能力下降；t19：飞机检修；t20：管制员指挥不正常航班运行；t21：未能及时疏散控制运行航班；t22：及时疏散空中运行航班；t23：设施设备检修；t24：设施设备故障未及时发现，并用于航班保障；t25：航空器撞击机场设施设备；t26：航空器之间相撞；t27：航空器搬移并检修；t28：在跑道上产生碰撞碎片；t29：在滑行道上产生碰撞碎片；t30：在停机位产生碰撞碎片；t31：跑道上碎片清除；
t32：跑道上碎片未及时清除；t33：滑行道上碎片清除；t34：滑行道上碎片未及时清除；t35：停机位上碎片清除；t36：停机位上碎片未及时清除；t37：跑道冲突未有效处置；t38：航空器搬移，跑道冲突解除；t39：滑行道冲突未有效处置；t40：航空器搬移，滑行道冲突解除；t41：停机位冲突未有效处置；t42：航空器搬移，停机位冲突解除；t43：指挥航班正常运行；t44：善后处理。
[0174] 第八步：将不正常航班衍生事件演化的随机事件链Petri网模型转化为马尔科夫链，基于马尔科夫链标识转移概率预测不正常航班衍生事件链式效应。
[0175] 通过分析基于随机事件链Petri网(SECPN)的不正常航班衍生事件演化模型，得出该模型的可实施与实施，即变迁激发规则：
[0176] 令SECPN＝(P,T,F,K,W,M,λ)为不正常航班衍生事件演化的随机事件链Petri网(SECPN)系统，M是(P,T,F)上的标识，t∈T。
[0177] (1)t在M有发生权的条件是：
[0178]
[0179] t在M有发生权的事实记为M[t>。
[0180] (2)若t在M有发生权，则t可发生并将M改为M'，记作M[t>M',M'由下式给出：
[0181]
[0182] 其中，*t＝{p|(p∈P∪T)∧((p,t)∈F)}称为t的前置集，t*＝{p|(p∈P∪T)∧((t,p)∈F)}称为t的后置集。
[0183] 根据上述随机事件链Petri网的变迁激发规则，对不正常航班衍生事件演化模型建立初始标识，得到不同变迁所得到的可达集。
[0184] 在连续时间内，基于随机事件链Petri网的不正常航班衍生事件演化模型中的每一个变迁t从可实施到实施需要延时，对 t变成可实施的时刻到它实施时刻之间被看成一个连续随机变量xt，其中xt取正实数值，且服从与指数分布函数，即
[0185]
[0186] 其中实参数λt>0是变迁t的平均实施速率，变量x≥0。
[0187] 将不正常航班衍生事件演化的随机事件链Petri网模型同构于齐次马尔可夫链，其各库所转换的触发条件为时间随机变量，并服从于指数分布，变迁t1,t2,…,t44的平均实施速率分别为λ1,λ2,…,λ44。根据不正常航班衍生事件演化的随机事件链Petri网模型的状态标识M1,M2,…,Mn，可得到其同构的马尔科夫链。马尔科夫链中有向弧表示不正常航班衍生事件演化过程的随机事件链Petri网模型从一个状态到另一个状态的转换。
[0188] 通过对现实生活中不正常航班及其衍生事件的观测，统计出t1,t2,…,t44的变迁概率，即得到不正常航班及其衍生事件的马尔可夫链中各标识状态发生变迁的概率λ1,λ2,…,λ44，根据各状态标识Mi的发生序列ta,tb,tc,…，进行归一化处理后得到ta,tb,tc,…的各个变迁权重ωa,ωb,ωc,…及其从可实施到实施所需要的近似时间τa,τb,τc,…。
[0189] 由于模型中τ服从于指数分布，Fi(τ)表示变迁ti发生所经历时间τ的概率分布，其概率分布函数为：
[0190]
[0191] 若变迁只有一个可达集，由于变迁ti的时间为τ，则状态标识Mu变迁到Mv的概率为：
[0192]
[0193] 若状态标识有两个或两个以上变迁，在时间τ内，变迁ti发生，tj,…,tk不发生，则状态标识Mu变迁到Mv的概率为：
[0194]
[0195] 对于一个给定的发生序列L＝{MuMi,…,MjMv}，其中u≠i≠j≠v,u,i,j,v∈L,[0196] 系统由状态u到状态v的发生概率为：
[0197]
[0198] 其中，|L|＝n，n为变迁个数，h为系统由状态u到状态v所有可能的变迁个数。
[0199] 由此，可对不正常航班事件可能引发的衍生事件进行预测和概率计算。
[0200] 本发明将不确定性因素引入到不正常航班事件中，利用事件链式效应思想，构建不确定性因素扰动下的不正常航班事件结构模型，定量描述不确定因素的作用机制；提出基于航班执行过程解析的不正常航班波及模型，描述不正常航班波及过程及演化规律，刻画了不正常航班在不同时间段和不同节点上对相关航班的波及影响；而后根据不正常航班衍生事件的链式波及效应，建立不正常航班衍生事件链的层次结构模型和链式效应解析模型，描述事件之间逻辑关系并推演不正常航班可能引发的衍生事件；在构建不正常航班衍生事件链式演化系统的基础上，建立不正常航班衍生事件演化的随机事件链Petri网模型，利用随机Petri网与马尔科夫链的同构关系，预测不正常航班可能引发的衍生事件及发生概率。借助事件链式效应模型梳理出不正常航班及其波及导致的下游航班或下游机场航班不正常等一系列衍生事件链，能够使单一航班延误对其他航班的影响得以凸显，在此基础上结合随机事件链Petri网模型预测不正常航班衍生事件的演化趋势和演化概率，增加不正常航班预测的全面性和准确性，提高不正常航班恢复的全局性，得到一种预测精度、预测范围均优于以往算法的不正常航班预测方法。
[0201] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
