CN104143170B - 低空救援空中交通调度指挥系统及其调度指挥方法 - Google Patents

低空救援空中交通调度指挥系统及其调度指挥方法 Download PDF

Info

Publication number
CN104143170B
CN104143170B CN201410361394.3A CN201410361394A CN104143170B CN 104143170 B CN104143170 B CN 104143170B CN 201410361394 A CN201410361394 A CN 201410361394A CN 104143170 B CN104143170 B CN 104143170B
Authority
CN
China
Prior art keywords
mrow
disaster
rescue
msub
information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410361394.3A
Other languages
English (en)
Other versions
CN104143170A (zh
Inventor
张明
周毅
王磊
王硕
李翔宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201410361394.3A priority Critical patent/CN104143170B/zh
Publication of CN104143170A publication Critical patent/CN104143170A/zh
Application granted granted Critical
Publication of CN104143170B publication Critical patent/CN104143170B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种低空救援空中交通调度指挥系统及其调度指挥方法,根据初始的航空器和物资分布信息,初始灾情信息,对于空中交通救援飞行计划和物资调度进行规划,获得航空器和物资的总调度方案确定救援总需求,在每个滚动的救援信息更新周期,通过信息熵对灾情信息融合实现短期灾情预测,采用模糊聚类算法建立灾区灾情的优先级划分,以空中交通信息和物资供需信息为约束,建立基于空中交通和物资信息的物资和航空器的调度模型,获取航空器和物资调度方案,本发明将空中交通信息和物资供需信息融入空中交通资源调度的决策中,提高航空救援资源调度的有效性和准确性。

Description

低空救援空中交通调度指挥系统及其调度指挥方法
技术领域
本发明涉及一种低空救援空中交通调度指挥系统及其调度指挥方法,属于民航航空救援领域。
背景技术
在抗灾救援及处置突发事件的各项措施中,航空救援具有快速、高效、受地理空间限制少等优势,是世界上许多国家普遍采用的最有效手段。在航空应急救援中,施救的航空器在真高1000米以下的低空运行环境由于受到地形环境、密集飞行等不确定因素的干扰,以及对于航空运力资源和物资的优化调度缺乏统筹考虑,往往存在安全风险大、施救效率低和实施不合理等问题。
目前,对于基于空中交通和物资信息的低空救援空中交通调度指挥的研究主要存在如下两个问题:
第一,在航空救援的航空器和物资的优化调度研究方面,目前的救援物资调度研究,一方面尽管有不少文献涉及物资调度模型,但并没有涉及救援信息变化和多物资下的航空应急救援调度,另一方面,目前较少将航空器调度和物资调度结合考虑,并且没有划分灾区灾情的优先级,以空中交通信息和物资供需信息为约束,建立基于空中交通和物资信息的物资和航空器的调度模型和方法。
第二,在基于空中交通和物资信息的低空救援空中交通调度指挥系统研究方面,目前没有空中交通信息和物资信息的动态变化下,同时实现航空器救援飞行计划优化以及航空救援物资和航空器优化调度方案的空中交通调度指挥系统,以至于存在航空器调度的利用率低,物资需求分析及投放的准确性低的问题。
具体参考文献如下:
1、朱新平,基于工作流技术的机场应急救援系统研究,南京航空航天大学硕士论文,2008年。
2、尧丰,韩松臣,朱新平.低空空域综合管理系统软件架构研究,武汉理工大学学报·信息与管理工程版,2010,32(2):257-260。
3、尧丰,基于RUP的低空空域管理系统关键技术初步研究,南京航空航天大学硕士论文,2010年。
发明内容
本发明所要解决的技术问题是:提供一种低空救援空中交通调度指挥系统及其调度指挥方法,首先根据初始灾情信息及物资信息,设计总调度方案,然后,根据实时获取的灾情及物资、航空器信息实时调整调度方案,解决了现有技术中航空器调度的利用率低,物资需求分析及投放的准确性低的问题。
本发明为解决上述技术问题,采用如下技术方案:
低空救援空中交通调度指挥方法,包括如下步骤:
步骤1、根据初始的航空器和物资分布信息,初始灾情信息,对于空中交通救援航空器飞行计划和物资调度进行规划,并确定航空器和物资调度信息的总调度方案,设定滚动的救援信息更新周期;
步骤2、首先,在每个滚动的救援信息更新周期内,实时采集各灾区的灾情信息,并采用信息熵的方法预测下一个救援信息更新周期的灾情信息,根据预测的灾情信息,预测物资的需求信息,然后,对灾区进行分组,确定各灾区对物资的需求信息,最后,根据当前的物资供应信息、航空器运力信息及各灾区对物资的需求信息,获取物资和航空器调度规划方案;
步骤3、将物资和航空器调度规划方案通过ADS-B地空数据链发布到各个航空器中;
步骤4、修正步骤1中的总调度方案中的航空器和物资调度信息,并将该航空器和物资调度信息通过网络发布。
所述步骤1包括:预先设定航空器的调度总架次Z1航空器参与救援总时间Z2其中,xfdcth为第f架航空器是否在第t天进行第h次由出救点d满载运往受灾点c的出程飞行,如果是xfdcth为1,否则为0,C为受灾点集合,F为航空器集合,D为出救点集合,k为总调度周期,s为航空器调度次数,ta表示每天实际作业时间,k、s均为正整数,航空器的调度总架次Z1、航空器参与救援总时间Z2满足如下约束条件:
其中rc为受灾点c调度周期内总的耐用品需求架次,vc为受灾点c快速消耗品每天的消耗架次,k'为耐用品的调度周期,k'为正整数;
其中,αc为受灾点c快速消耗品每天最大需求满足度;
其中,b为航空器每天允许的标准作业时间,ta为航空器每天实际作业时间,tL为航空器所需装载时间,tU为航空器所需卸载时间,tf为航空器所需加油时间,为q型航空器完成出程飞行所需时间,为q型航空器完成出程飞行所需时间;
其中,yfcdth为第f架航空器是否在第t天进行第h次由受灾点c到出救点d的入程飞行,如果是yfcdth为1,否则为0;
其中,efd为第f架航空器是否在制定调度计划前停靠在出救点d,如果是efd为1,否则为0,vfd,t-1为第f架航空器是否在第t-1天的结束调度飞行任务后停靠在出救点d,如果是vfd,t-1为1,否则为0;;
其中,vfdt为第f架航空器是否在第t天结束调度飞行任务后停靠在出救点d,如果是vfdt为1,否则为0;
所述步骤2包括如下步骤:
步骤2-1、在每个滚动的救援信息更新周期内,实时采集各灾区的灾情信息,并采用信息熵的方法预测下一个救援信息更新周期的灾情信息,预测模型为:
其中,表示在给定时间间隔t灾区i对物资l的时变需求量;al分别表示灾区i对于物资l每人次平均小时需求量以及对应的缓冲需求,其中l分别对应快速消耗品和耐用品;表示在给定的时间间隔t-ε内灾区i对物资l的时变需求量;L表示任意灾区连续两次物资分配预设的上限值;Rfmcg和Rdura分别表示快速消耗品和耐用品两类物资;z1-α表示选定统计值,其对应的时变物资需求短缺率被设定为α;δi(t)表示在给定的时间间隔t灾区i受困人数的估值;表示灾区i对物资l的时变需求量标准差;
步骤2-2、通过多准则模糊聚类方法对灾区进行分组;
步骤2-3、基于当前的物资供应信息、航空器运力信息及各灾区对物资的需求信息,建立物资调度和航空器飞行计划的规划模型,获取物资和航空器调度规划方案,其中,物资调度规划模型为:
航空器飞行计划的规划模型为:
其中,为各受灾点时变物资需求满足率,为物资配送的飞行总里程,M表示分配中心的数量,索引为m;Km表示分配中心m航空器架次,k为索引;ig为各个分组的受灾点数目,为决策变量,表示在时间间隔t由物资分配中心m运往受灾点ig物资l的数量,为受灾点ig物资l的需求量;nmk表示分配中心m处第k架航空器服务受灾点个数,i为索引;Rmj表示分配中心m区域内第j条路径;rmji表示受灾点ig∈g在分配中心m中的路径j中的顺序为i,令rmj0为配送中心,表示受灾点之间的距离,ig,ig'∈g,表示受灾点与分配中心之间的距离,m、k、l、i、i均为正整数。
所述步骤2-2中通过多准则模糊聚类方法对灾区进行分组的具体过程如下:
(1)初始数据处理,根据灾区受灾程度设定五种紧急属性,分别为:很高、高、中等、低、很低;
(2)模糊聚类分组,
首先,将五种紧急属性判别向量转化为二元数据,对于受灾点i,4×1的紧急属性判别向量为其中为在时间间隔t内受灾点i该地区的震级与烈度;为在时间间隔t内受灾点i受灾人数中伤亡人数的比例;为受灾点i当前时间t与上一次救援时间的差值;为在时间间隔t内受灾点i内的建筑物的损毁程度;
其次,生成I×I维的模糊关联矩阵其中I为受灾点个数,模糊关联矩阵中任意一个元素wpq(t)表示受灾点p和受灾点q之间的相互关系的紧密程度,p、q为正整数,
最后,通过聚类将具有相似紧急程度的受灾点划分到同一紧急属性组中,确定分组结果。
所述步骤2-3中物资调度规划模型、航空器飞行计划的规划模型均满足如下约束条件:
总的物资配送量其中,
每架航空器的最大载重量
每架航空器正常运输飞行状态下所能飞行的最大航程
其中,
每条路径的受灾点的组成Rnj={rnki|rnki∈{1,2,…,Lm},i=1,2,…,nmk}
航空器参与救援状况
对于任意受灾点,其任何时刻的配送量其中Ul表示物资l的单位体积,Qmk、fmk、cmk、vmk分别表示分配中心m处第k架航空器的最大载荷量、最大装载油量,平均耗油率和平均飞行速度;Lm为分配中心m的援助受灾点数目,N为所有受灾点总数目。
低空救援空中交通调度指挥系统,包括系统管理指挥终端、服务器、主控模块、ADS-B工作站、飞行状态数据库、航空器机载终端,以及救援物资信息源;飞行状态数据库通过ADS-B工作站接收航空器机载终端的空中交通信息,所述空中交通信息包括航空器性能数据和位置数据,所述服务器包括飞行计划优化服务器、救援调度管理服务器,其中,飞行计划优化服务器根据灾情信息、运力信息和物资信息,完成各个航空器初始路径集的优化,并将优化的路径集提供给救援调度管理服务器,救援调度管理服务器包括救援物资信息处理模块、航空器性能数据库、救援地理信息处理模块、应急预案规则库、救援调度模块和救援仿真显示模块;救援调度管理服务器根据空中交通信息及救援物资信息处理各模块的数据获取救援航空器和物资分布信息,并通过网络将救援航空器和物资分布信息传输至远程救援客户端显示,低空救援空中交通调度指挥系统管理指挥终端根据权限修改救援方案。
与现有技术相比,本发明具有如下有益效果:
1、根据初始的航空器和物资分布信息,初始灾情信息,对于空中交通救援飞行计划和物资调度进行规划,获得航空器和物资的总调度方案确定救援总需求,为进一步救援开展提供了航空运力和物资保障,也争取了宝贵的救援时间。
2、在每个滚动的救援信息更新周期内,通过信息熵对灾情信息融合实现短期灾情预测,采用模糊聚类算法建立灾区灾情的优先级划分,以空中交通信息和物资供需信息为约束,建立基于空中交通和物资信息的物资和航空器的调度模型,确定航空器和物资调度方案,相对于传统救援交通工具和物资调度,将空中交通信息和物资供需信息融入空中交通资源调度的决策中,提高航空救援资源调度的有效性和准确性。
附图说明
图1为本发明系统的结构示意图;
图2为本发明系统模块关系图;
图3为本发明方法流程图;
图4为采用模糊层次聚类方法依据灾情对灾区进行分组流程图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
本发明所述系统的结构如图1所示,包括系统管理指挥终端、服务器、主控模块、ADS-B工作站、飞行状态数据库、航空器机载终端,以及救援物资信息源;所述服务器包括飞行计划优化服务器、救援调度管理服务器。
航空器机载终端具备ADS-B OUT和IN功能的机载综合收发系统,救援信息物资信息源主要由救援与物资信息收发器完成对于救援物资和航空器的接收和发送。
飞行计划优化服务器根据灾情信息、运力信息和物资信息,完成对于各个航空器运行初始路径优化,并将优化的路径集提供给救援调度管理服务器。
飞行状态数据库通过ADS-B工作站接收航空器机载终端的空中交通信息,空中交通信息包括航空器性能数据和位置数据。
救援调度管理服务器提供救援物资信息处理模块、航空器性能数据库、救援地理信息处理模块、应急预案规则库、救援调度模块和救援仿真显示模块;救援调度管理服务器通过飞行状态数据库的空中交通信息和救援物资信息网络获取的物资信息,通过救援调度模块和应急预案规则库,结合初始路径集确定航空器和物资调度方案,通过网络传输到远程救援客户端(如电脑等固定远程终端和手机,笔记本等移动远程终端)显示救援航空器和物资分布信息,低空救援空中交通调度指挥系统管理指挥终端可根据权限修改救援方案。
具体系统的各个模块关系图,如图2所示,具体包括如下;
救援调度模块是系统的核心模块,接受救援信息和动态物资信息,同时通过航空器性能数据库和航空救援的应急预案库的检索,完成对于航空救援的调度计算,并通过输出端,完成物资和航空器调度信息的输出。
航空救援地理信息处理模块,通过救援网络,获取救援点与出救点的位置信息,灾情信息;通过ADS-B工作站,获取航空器的位置信息,结合GIS,实现对于航空救援运力和灾情的动态信息处理。
救援物资信息处理模块,通过救援通信网络,确定各个出救点的物资信息(如存储量,存储种类,救灾点物资需求量和需求种类信息等),为救援调度模块提供动态的调度物资供需信息。
航空器性能数据库,通过ADS-B,获取航空器的地速,燃油消耗,航向,高度和携带物资量和种类。应急预案规则库是根据航空救援特点和应急演练,确定的航空器救援基本规则。仿真显示模块,提供救援调度方案下航空器、出救点和救灾点物资和航空器的动态运行仿真,并提供救援演练下的仿真复现。
本发明所述低空救援空中交通调度指挥方法的流程如图3所示,具体包括如下:
(1)首先对于需要救援的灾区进行物资和航空器预先设定总体调度方案
通过对自然灾害条件下低空应急救援的分析,结合实际的航空器救援调度实施场景,建立了航空器调度的整数规划模型。该模型规定在一个调度周期内,使得出救航空器架次最大且总的调度时间最小,此外分别针对耐用品对应急时间要求高以及快速消费品持续消耗的不同特性,从多个出救点通过航空器向多个受灾点进行物资组合调度,满足调度周期内航空器出程飞行架次达到最大化,从而提高单位时间内的物资调运量,进而提高航空应急救援的效率;当出救架次达到饱和后,即无法出救更多架次的航空器时,使得参与救援的航空器总的调度时间最短,从而降低出救成本。
预先设定航空器的调度总架次Z1满足调度周期内航空器出程飞行架次达到最大化,从而提高单位时间内的物资调运量,进而提高航空应急救援的效率;航空器参与救援总时间Z2当出救架次达到饱和后,即无法出救更多架次的航空器时,使得参与救援的航空器总的调度时间最短,从而降低出救成本其中,xfdcth为第f架航空器是否在第t天进行第h次由出救点d满载运往受灾点c的出程飞行,如果是xfdcth为1,否则为0,C为受灾点集合,F为航空器集合,D为出救点集合,k为总调度周期,单位为天,s为航空器调度次数,ta表示每天实际作业时间,单位为分钟,k、s均为正整数,航空器的调度总架次Z1、航空器参与救援总时间Z2满足如下约束条件:
其中rc为受灾点c调度周期内总的耐用品需求架次,vc为受灾点c快速消耗品每天的消耗架次,k'为耐用品的调度周期,k'为正整数,保证两类物资中的耐用品在调度周期(k天)开始前k'天内完成所需资源调度,并且同时满足受灾点对快速消耗品的物资需求;
确保消耗性物资在耐用品调度完成后的周期内,耐用品的总的需求架次可以得到保障;
其中,αc为受灾点c快速消耗品每天最大需求满足度;确保调度周期内每天消耗品连续性需求得到满足,即避免出现消耗品耗尽;
其中,b为航空器每天允许的标准作业时间,ta为航空器每天实际作业时间,tL为航空器所需装载时间,tU为航空器所需卸载时间,tf为航空器所需加油时间,为q型航空器完成出程飞行所需时间,为q型航空器完成出程飞行所需时间;b、ta、tL、tU、tf的单位均为分钟;确保每天的实际调度时间在允许的调度作业时间限制之内;
其中,yfcdth为第f架航空器是否在第t天进行第h次由受灾点c到出救点d的入程飞行,如果是yfcdth为1,否则为0;确保总的出入程飞行架次相同,且航空器只有在上一次执行了入程飞行作业后,本次方可执行出程飞行任务;
其中,efd为第f架航空器是否在制定调度计划前停靠在出救点d,如果是efd为1,否则为0,vfd,t-1为第f架航空器是否在第t-1天的结束调度飞行任务后停靠在出救点d,如果是vfd,t-1为1,否则为0;确保航空器执行出程飞行的条件得到满足;
其中,vfdt为第f架航空器是否在第t天结束调度飞行任务后停靠在出救点d,如果是vfdt为1,否则为0;确保航空器执行调度任务后停靠在出救点的条件得到满足;
确保每架航空器在执行完当天的调度任务后只得停靠在一个出救点处;
所有变量均为二元变量,即变量数值仅可以取0和1两个数值。
预先设定滚动的救援信息更新周期,本发明每个滚动的救援信息更新周期设定为4小时,本发明所述的救援调度周期即为所述的滚动的救援信息更新周期。
通过预先设定上述总体调度,可实现对于初始灾情下,根据初始物资和航空器分布,实现对航空器和物资的总体调度,输出调度的计划和总体方案。
(2)救援调度周期内动态物资和航空器调度
在对上述航空救援总体调度的基础上,完成在救援信息变化下的救援调度周期内短期调度规划,首先需要根据多个信息源,对于物资需求信息的准确预测。
第一步,救援调度周期内动态物资需求预测,在预测过程中,通常需要多个信息源对搜集到的信息进行综合评判,以期对目标对象进行合理的预测,此时需要解决各个指标的相应权重赋值问题,从而真实反映各个信息源在评判和决策过程中所占有的地位和所起到的作用。而信息熵从信息的不确定性和概率测度的角度来表征信源的不定度,该方法可尽量消除各项指标权重计算的人为干扰,使得评价结果更为实际。对于每一个受灾点(用i表示),总共有ji种类型的信息源(比如记者,公共部门,私人救援队),且彼此之间具有相互独立性,伤亡人数观测值分为其中表示对于受灾点i,在时间间隔t内,对于给定类型ji的观测者的总人数。为了获取置信度的信息,我们进一步假设在时间间隔t内,受灾点i处,假定观测类型为ji的观测者所提供的观测值服从一个随时间变化的均值为方差为的高斯随机分布。随后,我们采用上述定义的高斯过程,总共将整个置信区间划分为M个信息置信水平以便分析不同渠道获得的信息的可靠性并对数据进行分级。
通过不同角度测量得到的伤亡数据,由于其获取途径不同,在统计预测中为了提高结果的准确度,须对不同两侧手段结果进行赋权,可信度不同加不同的权重。根据问题中多信息源特性,选择信息熵法。则下一步骤为熵值估算,首先,在时间间隔t内,受灾点i处,对于给定的信息源ji,我们可以采用上式计算给定置信水平为m的后验概率随后我们采用Shannon的信息熵,即熵值时间间隔t,信息源ji,有:并采用熵值越大,所获得的权重越小,确定合适的权重接下来,在时间间隔t内,对于受灾点i,累计伤亡总数Xi(t)即为各个组成部分的加权和,即到此为止,所有的实时伤亡数据均可以通过上述多信息源模糊模型求得,则相应的被困人数δi(t)为:最后,时变物资需求与该区的生还人数高度相关,则快速消费品和耐用品物资的需求预测模型为:
其中表示在给定时间间隔t灾区i对物资l的时变需求量;al分别表示灾区i对于物资l每人次平均小时需求量以及对应的缓冲需求,其中l分别对应快速消耗品和耐用品;表示在给定的时间间隔t-ε内灾区i对物资l的时变需求量;L表示任意灾区连续两次物资分配预设的上限值;Rfmcg和Rdura分别表示快速消耗品和耐用品两类物资;z1-α表示选定统计值,其对应的时变物资需求短缺率被设定为α;δi(t)表示在给定的时间间隔t灾区i受困人数的估值;表示灾区i对物资l的时变需求量标准差。
第二步,救援调度周期内基于灾情的灾区分组
为了便于物资分配,该阶段通过多准则模糊聚类方法对灾区进行分组,具体流程如图4所示,对于各受灾点供需不匹配下的动态物资配送具有重要意义,既可以更为有效的判定各受灾点需求程度,又能达到更为有效的救灾效果。具体步骤如下:
第一分步骤,初始数据处理,根据灾区受灾程度设定五种紧急属性,分别为:很高、高、中等、低、很低;
第二分步骤,模糊聚类分组,
首先,将五种紧急属性判别向量转化为二元数据,对于受灾点i,4×1的紧急属性判别向量为Vi(t):其中为在时间间隔t内受灾点i该地区的震级与烈度;为在时间间隔t内受灾点i受灾人数中伤亡人数的比例;为受灾点i当前时间t与上一次救援时间的差值;为在时间间隔t内受灾点i内的建筑物的损毁程度;
其次,生成I×I维的模糊关联矩阵其中I为受灾点个数,模糊关联矩阵中任意一个元素wpq(t)表示受灾点p和受灾点q之间的相互关系的紧密程度,p、q为正整数,
最后,通过聚类将具有相似紧急程度的受灾点划分到同一紧急属性组中,确定分组结果且使得各组之间的紧急属性差异较为明显。具体看附图4。
第三步,救援调度周期内航空器和物资调度
通过数据融合完成对各个受灾点的伤亡人数统计以及所需各种物资量的判定,通过模糊聚类方法对各个受灾点进行灾区分组,确定了各个受灾点的物资救灾响应优先级。本发明采用一个复合加权多目标模型来处理多个应急物资配送中心向多个受灾点的多类应急物资分配问题。即在各个调度周期各个受灾点优先级及各类物资需求量确定的情况下,通过合理安排航空器的出救架次和时间,使得各个受灾点的物资满足率最大化,且同时使得调度飞行航程最小化。
对于每一个受灾点分组g,在给定的时间间隔t内,使得各个受灾点时变物资需求满足率最大化,即其二是物资配送的飞行总里程最短,即其中M表示分配中心的数量,索引为m;Km表示分配中心m航空器架次,k为索引;ig为各个分组的受灾点数目,为决策变量,表示在时间间隔t由物资分配中心m运往受灾点ig物资l的数量。为受灾点ig物资l的需求量;nmk表示分配中心m处第k架航空器服务受灾点个数,i为索引;Rmj表示分配中心m区域内第j条路径;rmji表示受灾点ig∈g在分配中心m中的路径j中的顺序为i(不含配送中心),令rmj0为配送中心,(ig,ig'∈g)表示受灾点之间的距离,表示受灾点与分配中心之间的距离。同时,满足如下几个约束条件:
①确保在任意救援时段,对于任意受灾点,其总的配送量不得超过其实际时变需求量,即
②确保每架航空器不得超载运输飞行,即
③确保每架航空器的飞行距离不得超过其正常运输飞行状态下所能飞行的最大航程
④确保各个分配中心辖下的所有受灾点都获得救援
⑤确保所有受灾点都得到救助
⑥表示每条路径的受灾点的组成Rnj={rnki|rnki∈{1,2,…,Lm},i=1,2,…,nmk};
⑦限制每个受灾点只能由一架航空器进行救援
⑧表示当分配中心m处第k架航空器的救援受灾点数量大于等于1时,则说明该航空器参与了救援,取sign(nmk)=1,否则未参与救援,
⑨确保对于任意受灾点,其任何时刻的配送量的非负性,使得问题具有实际意义,即
其中Ul表示物资l的单位体积,Qmk、fmk、cmk、vmk分别表示分配中心m处第k架航空器的最大载荷量、最大装载油量,平均耗油率和平均飞行速度;Lm为分配中心m的援助受灾点数目,N为所有受灾点总数目。
以上仅为本发明的较佳实施例,但本发明的保护范围并不局限于此,任何熟悉本领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求所界定的保护范围为准。

Claims (5)

1.低空救援空中交通调度指挥方法,其特征在于:包括如下步骤:
步骤1、根据初始的航空器和物资分布信息,初始灾情信息,对于空中交通救援航空器飞行计划和物资调度进行规划,并确定航空器和物资调度信息的总调度方案,设定滚动的救援信息更新周期;
步骤2、首先,在每个滚动的救援信息更新周期内,实时采集各灾区的灾情信息,并采用信息熵的方法预测下一个救援信息更新周期的灾情信息,根据预测的灾情信息,预测物资的需求信息,然后,对灾区进行分组,确定各灾区对物资的需求信息,最后,根据当前的物资供应信息、航空器运力信息及各灾区对物资的需求信息,获取物资和航空器调度规划方案;
所述步骤2包括如下步骤:
步骤2-1、根据救援需求信息多信息源特性,采用信息熵法确定救援信息:首先,在时间间隔t内,受灾点i处,对于给定的信息源ji,计算给定置信水平为m的后验概率其中,表示伤亡人数观测值,表示对于受灾点i,在时间间隔t内,对于ji的观测者的总人数;其次,采用Shannon的信息熵,即熵值时间间隔t,信息源ji,有:并采用熵值越大,所获得的权重越小,确定合适的权重最后,在时间间隔t内,对于受灾点i,累计伤亡总数Xi(t)即为各个组成部分的加权和,即其中,表示观测类型为ji的观测者所提供的观测值服从一个随时间变化的高斯随机分布的均值;所有的实时伤亡数据均能够通过上述内容求得,则相应的被困人数δi(t)为:其中,Si表示受灾点i处的社会统计总人口;
步骤2-2、在每个滚动的救援信息更新周期内,实时采集各灾区的灾情信息,并采用信息熵的方法预测下一个救援信息更新周期的灾情信息,预测模型为:
<mrow> <msubsup> <mi>D</mi> <mi>i</mi> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>{</mo> <msup> <mi>a</mi> <mi>l</mi> </msup> <mo>&amp;times;</mo> <msub> <mi>&amp;delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mi>L</mi> <mo>+</mo> <msub> <mi>z</mi> <mrow> <mn>1</mn> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>&amp;times;</mo> <msubsup> <mi>STD</mi> <mi>i</mi> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <msqrt> <mi>L</mi> </msqrt> <mo>,</mo> <mn>0</mn> <mo>}</mo> </mrow> </mtd> <mtd> <mrow> <mi>l</mi> <mo>&amp;Element;</mo> <msub> <mi>R</mi> <mrow> <mi>f</mi> <mi>m</mi> <mi>c</mi> <mi>g</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>{</mo> <msup> <mi>a</mi> <mi>l</mi> </msup> <mo>&amp;times;</mo> <msub> <mi>&amp;delta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>b</mi> <mi>i</mi> <mi>l</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>&amp;epsiv;</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msubsup> <mi>A</mi> <mi>i</mi> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;epsiv;</mi> <mo>)</mo> </mrow> <mo>,</mo> <mn>0</mn> <mo>}</mo> </mrow> </mtd> <mtd> <mrow> <mi>l</mi> <mo>&amp;Element;</mo> <msub> <mi>R</mi> <mrow> <mi>d</mi> <mi>u</mi> <mi>r</mi> <mi>a</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,表示在给定时间间隔t灾区i对物资l的时变需求量;al分别表示灾区i对于物资l每人次平均小时需求量以及对应的缓冲需求,其中l分别对应快速消耗品和耐用品;表示在给定的时间间隔t-ε内灾区i对物资l的时变需求量;L表示任意灾区连续两次物资分配预设的上限值;Rfmcg和Rdura分别表示快速消耗品和耐用品两类物资;z1-α表示选定统计值,其对应的时变物资需求短缺率被设定为α;δi(t)表示在给定的时间间隔t灾区i受困人数的估值;表示灾区i对物资l的时变需求量标准差;
步骤2-3、通过多准则模糊聚类方法对灾区进行分组;
步骤2-4、基于当前的物资供应信息、航空器运力信息及各灾区对物资的需求信息,建立物资调度和航空器飞行计划的规划模型,获取物资和航空器调度规划方案,其中,物资调度规划模型为:
<mrow> <mi>max</mi> <mi> </mi> <msubsup> <mi>F</mi> <mi>g</mi> <mn>1</mn> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mo>&amp;Sigma;</mo> <mrow> <mo>&amp;ForAll;</mo> <mi>l</mi> </mrow> </msub> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </msubsup> <msub> <mo>&amp;Sigma;</mo> <mrow> <mo>&amp;ForAll;</mo> <msub> <mi>i</mi> <mi>g</mi> </msub> <mo>&amp;Element;</mo> <mi>g</mi> </mrow> </msub> <msubsup> <mi>X</mi> <mrow> <mi>m</mi> <mo>,</mo> <msub> <mi>i</mi> <mi>g</mi> </msub> </mrow> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mo>&amp;ForAll;</mo> <mi>l</mi> </mrow> </munder> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <munder> <mo>&amp;Sigma;</mo> <mrow> <mo>&amp;ForAll;</mo> <msub> <mi>i</mi> <mi>g</mi> </msub> <mo>&amp;Element;</mo> <mi>g</mi> </mrow> </munder> <msubsup> <mi>D</mi> <msub> <mi>i</mi> <mi>g</mi> </msub> <mi>l</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mo>&amp;ForAll;</mo> <mrow> <mo>(</mo> <mi>g</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
航空器飞行计划的规划模型为:
<mrow> <mi>min</mi> <mi> </mi> <msubsup> <mi>F</mi> <mi>g</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mo>{</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>K</mi> <mi>m</mi> </msub> </munderover> <mo>&amp;lsqb;</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mrow> <mi>m</mi> <mi>k</mi> </mrow> </msub> </munderover> <msub> <mi>d</mi> <mrow> <msub> <mi>r</mi> <mrow> <mi>m</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>m</mi> <mi>k</mi> <mi>i</mi> </mrow> </msub> </mrow> </msub> <mo>+</mo> <msub> <mi>d</mi> <mrow> <msub> <mi>r</mi> <mrow> <msub> <mi>mjn</mi> <mrow> <mi>m</mi> <mi>k</mi> </mrow> </msub> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>m</mi> <mi>j</mi> <mn>0</mn> </mrow> </msub> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mrow> <mi>m</mi> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>}</mo> <mo>,</mo> </mrow>
其中,g表示受灾点分组,为各受灾点时变物资需求满足率,为物资配送的飞行总里程,M表示分配中心的数量,索引为m;Km表示分配中心m航空器架次,k为索引;ig为各个分组的受灾点数目,为决策变量,表示在时间间隔t由物资分配中心m运往受灾点ig物资l的数量,为受灾点ig物资l的需求量;nmk表示分配中心m处第k架航空器服务受灾点个数,i为索引;rmj表示分配中心m区域内第j条路径;rmji表示受灾点ig∈g在分配中心m中的路径j中的顺序为i,令rmj0为配送中心,m、k、l、i均为正整数;
步骤3、将物资和航空器调度规划方案通过ADS-B地空数据链发布到各个航空器中;
步骤4、修正步骤1中的总调度方案中的航空器和物资调度信息,并将该航空器和物资调度信息通过网络发布。
2.根据权利要求1所述低空救援空中交通调度指挥方法,其特征在于:所述步骤1包括:预先设定航空器的调度总架次Z1航空器参与救援总时间Z2其中,xfdcth为第f架航空器是否在第t天进行第h次由出救点d满载运往受灾点c的出程飞行,如果是xfdcth为1,否则为0,C为受灾点集合,F为航空器集合,D为出救点集合,k为总调度周期,s为航空器调度次数,ta表示每天实际作业时间,k、s均为正整数,航空器的调度总架次Z1、航空器参与救援总时间Z2满足如下约束条件:
其中rc为受灾点c调度周期内总的耐用品需求架次,vc为受灾点c快速消耗品每天的消耗架次,k'为耐用品的调度周期,k'为正整数;
其中,αc为受灾点c快速消耗品每天最大需求满足度;
其中,b为航空器每天允许的标准作业时间,ta为航空器每天实际作业时间,tL为航空器所需装载时间,tU为航空器所需卸载时间,tf为航空器所需加油时间,为q型航空器完成出程飞行所需时间,为q型航空器完成出程飞行所需时间;
其中,yfcdth为第f架航空器是否在第t天进行第h次由受灾点c到出救点d的入程飞行,如果是yfcdth为1,否则为0;
其中,efd为第f架航空器是否在制定调度计划前停靠在出救点d,如果是efd为1,否则为0,vfd,t-1为第f架航空器是否在第t-1天的结束调度飞行任务后停靠在出救点d,如果是vfd,t-1为1,否则为0;
其中,vfdt为第f架航空器是否在第t天结束调度飞行任务后停靠在出救点d,如果是vfdt为1,否则为0;
3.根据权利要求1所述低空救援空中交通调度指挥方法,其特征在于:所述步骤2-3中通过多准则模糊聚类方法对灾区进行分组的具体过程如下:
(1)初始数据处理,根据灾区受灾程度设定五种紧急属性,分别为:很高、高、中等、低、很低;
(2)模糊聚类分组,
首先,将五种紧急属性判别向量转化为二元数据,对于受灾点i,4×1的紧急属性判别向量为Vi(t):其中为在时间间隔t内受灾点i该地区的震级与烈度;为在时间间隔t内受灾点i受灾人数中伤亡人数的比例;为受灾点i当前时间t与上一次救援时间的差值;为在时间间隔t内受灾点i内的建筑物的损毁程度;
其次,生成I×I维的模糊关联矩阵其中I为受灾点个数,模糊关联矩阵中任意一个元素wpq(t)表示受灾点p和受灾点q之间的相互关系的紧密程度,p、q为正整数,
最后,通过聚类将具有相似紧急程度的受灾点划分到同一紧急属性组中,确定分组结果。
4.根据权利要求1所述低空救援空中交通调度指挥方法,其特征在于:所述步骤2-4中物资调度规划模型、航空器飞行计划的规划模型均满足如下约束条件:
总的物资配送量其中,
每架航空器的最大载重量
每架航空器正常运输飞行状态下所能飞行的最大航程
其中,
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>K</mi> <mi>m</mi> </msub> </munderover> <msub> <mi>n</mi> <mrow> <mi>m</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>L</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mo>&amp;ForAll;</mo> <mi>m</mi> <mo>)</mo> </mrow> </mrow>
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>L</mi> <mi>m</mi> </msub> <mo>=</mo> <mi>N</mi> </mrow>
每条路径的受灾点的组成Rnj={rnki|rnki∈{1,2,…,Lm},i=1,2,…,nmk}
航空器参与救援状况
对于任意受灾点,其任何时刻的配送量其中Ul表示物资l的单位体积,Qmk、fmk、cmk、vmk分别表示分配中心m处第k架航空器的最大载荷量、最大装载油量、平均耗油率和平均飞行速度;Lm为分配中心m的援助受灾点数目,N为所有受灾点总数目。
5.低空救援空中交通调度指挥系统,其特征在于:包括系统管理指挥终端、服务器、主控模块、ADS-B工作站、飞行状态数据库、航空器机载终端,以及救援物资信息源;飞行状态数据库通过ADS-B工作站接收航空器机载终端的空中交通信息,所述空中交通信息包括航空器性能数据和位置数据,所述服务器包括飞行计划优化服务器、救援调度管理服务器,其中,飞行计划优化服务器根据灾情信息、运力信息和物资信息,完成各个航空器初始路径集的优化,并将优化的路径集提供给救援调度管理服务器,救援调度管理服务器包括救援物资信息处理模块、航空器性能数据库、救援地理信息处理模块、应急预案规则库、救援调度模块和救援仿真显示模块;救援调度管理服务器根据空中交通信息及救援物资信息处理各模块的数据获取救援航空器和物资分布信息,并通过网络将救援航空器和物资分布信息传输至远程救援客户端显示,低空救援空中交通调度指挥系统管理指挥终端根据权限修改救援方案;航空器机载终端具备ADS-B OUT和IN功能的机载综合收发系统,救援信息物资信息源由救援与物资信息收发器完成对于救援物资和航空器的接收和发送。
CN201410361394.3A 2014-07-25 2014-07-25 低空救援空中交通调度指挥系统及其调度指挥方法 Active CN104143170B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410361394.3A CN104143170B (zh) 2014-07-25 2014-07-25 低空救援空中交通调度指挥系统及其调度指挥方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410361394.3A CN104143170B (zh) 2014-07-25 2014-07-25 低空救援空中交通调度指挥系统及其调度指挥方法

Publications (2)

Publication Number Publication Date
CN104143170A CN104143170A (zh) 2014-11-12
CN104143170B true CN104143170B (zh) 2018-01-30

Family

ID=51852338

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410361394.3A Active CN104143170B (zh) 2014-07-25 2014-07-25 低空救援空中交通调度指挥系统及其调度指挥方法

Country Status (1)

Country Link
CN (1) CN104143170B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105513432B (zh) * 2016-01-14 2017-11-17 中国民用航空总局第二研究所 用于加油调度的航班信息处理方法及系统
CN109540166A (zh) * 2018-11-30 2019-03-29 电子科技大学 一种基于高斯过程的安全路径规划方法
CN109636282B (zh) * 2018-12-10 2020-04-10 中南大学 一种航空货邮配送方法
CN109948918B (zh) * 2019-03-08 2021-04-20 北京交通大学 属地应急众储物资的综合分配方法
CN110222919A (zh) * 2019-04-18 2019-09-10 天津科技大学 应急救援平台
CN110322169B (zh) * 2019-07-12 2020-11-10 河北冀联人力资源服务集团有限公司 一种任务发布的方法及装置
CN112187717B (zh) * 2020-08-27 2021-10-29 中国人民解放军军事科学院国防科技创新研究院 一种基于救援场景的数据链消息处理方法
CN113159397B (zh) * 2021-03-31 2022-10-11 上海城市地理信息系统发展有限公司 一种救灾物资预存管理方法、装置和电子设备
CN115392814A (zh) * 2022-11-01 2022-11-25 中科雨辰科技有限公司 一种配置任务的数据处理系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102549558A (zh) * 2009-08-12 2012-07-04 Dal全球服务公司有限责任公司 实时服务管理与调度系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050197848A1 (en) * 2004-03-08 2005-09-08 Chou Y. H. Airport customer support dispatch system and method for operation for the same

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102549558A (zh) * 2009-08-12 2012-07-04 Dal全球服务公司有限责任公司 实时服务管理与调度系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《基于GIS的空管应急指挥辅助显示系统研究》;陈文涛;《中国优秀硕士学位论文全文数据库信息科技辑》;20111215;论文摘要,第3页,第31页,第41页,第44页 *

Also Published As

Publication number Publication date
CN104143170A (zh) 2014-11-12

Similar Documents

Publication Publication Date Title
CN104143170B (zh) 低空救援空中交通调度指挥系统及其调度指挥方法
Shavarani et al. A congested capacitated multi-level fuzzy facility location problem: An efficient drone delivery system
Chowdhury et al. Drones for disaster response and relief operations: A continuous approximation model
Ploetner et al. Long-term application potential of urban air mobility complementing public transport: an upper Bavaria example
Macias et al. Optimal hub selection for rapid medical deliveries using unmanned aerial vehicles
CN114066062B (zh) 一种城市空中交通的物流需求预测方法及系统
CN103116825A (zh) 智慧城市管理系统
CN105894847A (zh) 一种云平台环境下无监督学习的实时公交动态调度系统及方法
CN110288118A (zh) 一种基于无人机的空中救援响应方法
CN109977066A (zh) 属地智慧应急投送一体机
Szeto et al. Link-based multi-class hazmat routing-scheduling problem: A multiple demon approach
van Steenbergen et al. A simulation framework for UAV-aided humanitarian logistics
Zhang et al. A multi-trip electric bus routing model considering equity during short-notice evacuations
Chen et al. Multi-objective optimisation model of emergency material allocation in emergency logistics: A view of utility, priority and economic principles
Lu et al. Optimizing distribution of droneports for emergency monitoring of flood disasters in China
CN110059893A (zh) 一种含规划机场的大型繁忙机场终端区起降量预测方法
Radzki et al. Periodic planning of UAVs' fleet mission with the uncertainty of travel parameters
Hansen et al. Multiple airport systems in the United States: current status and future prospects
Yang et al. A novel integration platform to reduce flight delays in the National Airspace System
Yang et al. Machine learning-enabled adaptive air traffic recommendation system for disaster evacuation
Arreeras et al. An improvement on shelter airport selection model during large-scale volcanic disasters: A case study of Hakoneyama Japan
Ma et al. A data-driven optimal method for massive passenger flow evacuation at airports under large-scale flight delays
Pak et al. Can Urban Air Mobility become reality? Opportunities, challenges and selected research results
Zhu Decision making under uncertainties for air traffic flow management
Maheshwari Enabling Ride-Sharing in on-Demand Air Service Operations Through Reinforcement Learning

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant