CN115206135A - 一种不确定爬升率的航空器指令高度规划方法 - Google Patents
一种不确定爬升率的航空器指令高度规划方法 Download PDFInfo
- Publication number
- CN115206135A CN115206135A CN202210720445.1A CN202210720445A CN115206135A CN 115206135 A CN115206135 A CN 115206135A CN 202210720445 A CN202210720445 A CN 202210720445A CN 115206135 A CN115206135 A CN 115206135A
- Authority
- CN
- China
- Prior art keywords
- sector
- aircraft
- height
- cfl
- conflict
- 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.)
- Granted
Links
- 230000009194 climbing Effects 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000033001 locomotion Effects 0.000 claims description 15
- SLVOVFVZZFUEAS-UHFFFAOYSA-N 2-[2-[2-[bis(carboxymethyl)amino]ethoxy]ethyl-(carboxymethyl)amino]acetic acid Chemical compound OC(=O)CN(CC(O)=O)CCOCCN(CC(O)=O)CC(O)=O SLVOVFVZZFUEAS-UHFFFAOYSA-N 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 241000287196 Asthenes Species 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims 1
- 230000001105 regulatory effect Effects 0.000 claims 1
- 238000007726 management method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0047—Navigation or guidance aids for a single aircraft
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/003—Flight plan management
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提供了一种不确定爬升率的航空器指令高度规划方法,包括:步骤1,获取雷达航迹历史数据;步骤2,得到不同机型在不同方向不同时段的爬升率范围;步骤3,得到进扇区的指令高度集合和出扇区的指令高度集合;步骤4:计算航空器爬升率和下降率;步骤5:得到航空器可能的高度飞行剖面集合;步骤6:根据冲突探测识别结果对进扇区的指令高度集合和出扇区的指令高度集合进行分类;步骤7:进行冲突解脱。本发明为面向管制员在管制过程中提供了一种不确定爬升率的航空器指令高度规划方法,以辅助空中交通管理人员实现航空有效高效的运行,同时较少管制负荷,为下一代基于时间运行的空管自动化系统提供了技术支撑。
Description
技术领域
本发明涉及空中交通管理技术领域,尤其涉及不确定爬升率情况下航空器指令高度规划方法。
背景技术
在空中交通管理过程中,往往涉及到为了调节空域扇区内航空器之间的冲突,需要进行高度的调整,一般管制员会为航空器发出指令高度,让航空器在在指令高度保持。实际的操作过程中,对于如何选择最优化的指令高度,目前没有较好的方法,主要是通过管制员的经验来确定指令高度的值;这种方法存在以下问题:1)无法保证选择的指令高度是最优的;2)需要管制员自己计算并确定指令高度值,增加管制员的管制负荷;3)对于不确定爬升率性的情况,管制员很难通过大脑实现指令高度的计算和选择;因此基于上述挑战,需要提出不确定爬升率情况下航空器指令高度规划方法。
发明内容
发明目的:本发明所要解决的技术问题是针对当前或者下一代空管需求,从减少管制员管制负荷以及管制效率的角度出发,提供不确定爬升率情况下航空器指令高度规划方法,辅助管制员进行指令高度的选择,包括如下步骤:
步骤1,获取飞行计划信息中的航空器的关键航路点、经过的扇区编号等信息,进出扇区的高度限制、速度限制等信息;获取指定扇区的航空器的雷达航迹历史数据;
步骤2,通过统计分析的方法,得到不同机型在不同方向、不同时段的历史爬升率的均值和方差,以均值和方差得到不同机型在不同方向不同时段的爬升率范围;
步骤3,根据扇区空域的限制以及扇区移交协议中移交点的限制信息,选取能够用于指令高度设置的值,包括进扇区的指令高度集合和出扇区的指令高度集合;
步骤4:计算航空器爬升率和下降率;
步骤5:基于航空器爬升率和下降率,根据航空器当前高度、进扇区的指令高度集合和出扇区的指令高度集合,得到航空器可能的高度飞行剖面集合Trijk;
步骤6:计算扇区内各航空器潜在的运动区域,进行冲突探测识别;根据冲突探测识别结果,基于冲突发生的概率对进扇区的指令高度集合和出扇区的指令高度集合进行分类;
步骤7:针对存在潜在冲突风险指令高度CFL的值,基于调整高度、偏航和调整速度进行冲突解脱;
步骤8,重复步骤1~步骤7,直至管制扇区内各航空器具有无冲突的指令高度值,以辅助空中交通管理人员实现航空有效高效的运行,同时较少管制负荷。
步骤2包括:
步骤2-1:对获取的一段时间内的航空器雷达航迹历史数据按照航空器类型、经过的扇区、不同时间段以及不同的飞行方向(东单西双)进行爬升率数据分类,一段时间一般以一个月为周期;
步骤2-2:设样本数量为N,计算不同机型在不同方向、不同时段的历史爬升率的均值mean_ROCijk和方差st_ROCijk:
式中,ROCijk表示机型为i,方向为j,时段为k爬升率的值;
步骤2-3:得到不同机型在不同方向不同时段的爬升率范围,表示为:
[mean_ROCijk-st_ROCijk,mean_ROCijk+st_ROCijk]。
步骤3包括:
步骤3-1:管制员输入管制扇区空域的范围限制,包括扇区的编号、长、宽、高
步骤3-2:根据步骤3-1的输入,基于集合交集计算方法,计算管制扇区空间区域与管制扇区移交协议的限制区域的重叠区域,计算得到的交集记为QCFL;
步骤3-3:根据高度层的设置规则(东单西双,300米一个高度层),在交集QCFL中选取进扇区的指令高度集合CFL_set_in和出扇区的指令高度集合CFL_set_out,表示如下:
CFL_set_in={in_CFLijk|in_CFLijk满足扇区空域限制以及进扇区移交协议限制}
CFL_set_out={out_CFLijk|out_CFLijk满足扇区空域限制以及出扇区移交协议限制}
其中,in_CFLijk表示机型为i,方向为j,时段为k,且在交集QCFL中根据高度层的设置规则选取的进扇区的指令高度,所有的指令高度in_CFLijk组成进扇区的指令高度集合CFL_set_in;out_CFLijk表示机型为i,方向为j,时段为k,且在交集QCFL中根据高度层的设置规则选取的出扇区的指令高度,所有的指令高度out_CFLijk组成进扇区的指令高度集合CFL_set_out。
步骤4包括:为推算飞行高度剖面,通过如下公式计算航空器爬升率或下降率:
其中,ROCD表示航空器爬升率或下降率,T表示大气温度,单位为K;ΔT表示温度差,单位为K;VTAS表示真空速,由地速和风速的差可以获得;m表示航空器质量;g0=9.80665m/s2表示重力加速度。f{M}表示能量分享因子,表示在沿着固定的速度剖面爬升过程中,分配给爬升的能量与分配给加速的能量的比值。
步骤4中,按照BADA用户手册中提供的不同爬升段的能量分享因子的计算公式,分别计算能量分享因子,具体包括:
步骤4-1,根据BADA用户手册,计算加速爬升段的能量分享因子;
步骤4-2,根据BADA用户手册,计算等表速爬升段的能量分享因子;
步骤4-3,根据BADA用户手册,计算等马赫爬升段的能量分享因子。
步骤5包括:根据爬升率或下降率计算高度增量,计算公式如下:
Δh=ROCD·Δt
式中,Δt是系统计算采样时间周期;
根据航空器当前高度、进扇区的指令高度集合和出扇区的指令高度集合,根据步骤4爬升率计算公式,基于高度增量,计算得到航空器可能的高度飞行剖面集合Trijk。
步骤6包括:
步骤6-1:根据航空器当前的运行状态,具体包括通过雷达获取的当前的位置、速度信息,推算未来一段时间(15分钟-20分钟)航空可能的位置分布,位置分布的区域即为航空器的潜在的运动区域;
步骤6-2:计算航空器潜在运动区域与其他航空器潜在运动区域、特殊空域、飞行受限区的交集Q,计算公式为:
Q=STP1∩STP2…∩STPm
其中,STPi表示第i个航空器潜在运动区域、特殊空域、飞行限制空域,i=1,2,…m;
步骤6-3:计算潜在冲突概率P:
其中p(xu,yu,hu,tu)表示航空器位于坐标(xu,yu,hu,tu)的概率,u∈Q表示潜在冲突区域内的任意点,xu,yu,hu,tu分别表示航空器在任意点u的横轴的坐标值、纵轴的坐标值、竖轴的坐标值,以及过u的时间值,n表示相交区域内离散点个数;
步骤6-4:根据设定的阈值,如果潜在冲突概率大于阈值,则判定存在冲突,阈值一般取冲突的概率值,如取值0.5。
步骤6还包括:根据步骤6-4得到的结果,对进扇区的指令高度集合CFL_set_in和出扇区的指令高度集合CFL_set_out进行分类,具体包括如下三类:
当潜在冲突概率为零时,表示无冲突,对应的指令高度值是安全的,在系统设计中可用绿色表示;
当潜在冲突概率大于0但小于规定的阈值,存在潜在冲突风险,表示指令高度潜在安全风险;在系统设计中可用黄色表示;
潜在冲突概率大于规定的阈值,表示高冲突风险,指令高度值不安全,不可用,在系统设计中可用红色表示。
步骤7包括:
步骤7-1:基于速度控制实现冲突解脱,设当前速度为V,调速方式为然后计算水平和垂直方向上航空器之间的距离,如果满足飞行安全规则规定的最小安全间隔,则实现冲突解脱;如果不满足最小间隔,则执行步骤7-2和步骤7-3;
步骤7-2:针对水平方向有冲突且不能通过调速解决的情况,通过偏航的方式实现冲突解脱,允许的最大偏航角度为30度、45度和60度;允许的转弯时间T1计算公式如下:
其中DLmax表示最大偏航距离,trangle表示转弯角度,一般取值30度、45度或者60度;
步骤7-3:针对垂直方向有冲突且不能通过调速来解决的,通过调整高度实现冲突解脱;
步骤7-4:针对存在冲突的航迹实现冲突解脱后,执行步骤6-1~步骤6-4,进一步进行冲突探测;如果存在冲突,重复步骤7-1~步骤7-4,直至无冲突。
步骤7-3包括:如果需要到达的指令高度高于航空器当前的位置,则通过上升m个高度层,H的值由QCFL的值以及高度层划设规则计算得到,需要满足H∈QCFL且是300的整数倍;如果需要到达的指令高度低于航空器当前的位置,则通过下降m个高度层,m的获取方式同上;
相对于现有技术,本发明具有如下技术效果:首先,本发明为管制员选择指令高度提供了一种可实现的方法,基于该方法形成的辅助管制员的系统,能够减少管制员的管制负荷;其次,本发明为下一代空管系统提供了技术支持。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1是本发明实施部分提供的不确定性爬升率情况下航空器指令高度规划方法的工作流程示意图;
图2是本发明实施部分提供的不确定爬升率性情况下指令高度值的范围示意图;
图3是本发明实施部分提供的通过调整航向进行冲突解脱的示意图。
具体实施方式
如图1所示,本发明公开了一种不确定爬升率的航空器指令高度规划方法,属于空中交通管理领域。该方法首先获取飞行计划信息中的航空器的关键航路点、经过的扇区编号等信息,进出扇区的高度限制、速度限制等信息作为输入;通过统计分析的方法扇区的航空器的雷达航迹历史数据进行分析;得到不同机型在不同方向、不同时段的历史爬升率的均值和方差;然后,根据扇区空域的限制以及扇区移交协议中移交点的限制信息,选取可用于指令高度设置的值;引入航空器爬升率(或下降率)公式,计算航空器爬升下降率,根据航空器当前高度、进/出扇区的目标指令高度集,得到航空器可能的高度飞行剖面集合,计算扇区内各航空器潜在的运动区域,基于航空器的潜在运动区域之间,以及与特殊空域、飞行受限区的交集进行冲突探测识别;最后,根据冲突探测结果,基于冲突发生的概率对冲突进行分类,并基于调整高度、偏航以及调整速度进行冲突解脱;重复上述步骤,直至管制扇区内各航空器具有最优化的指令高度值。
本发明主要针对不确定爬升率情况下指令高度的规划选取提供一种技术,基于该方法形成的辅助管制员的系统,能够减少管制员的管制负荷;其次,本发明为下一代空管系统提供了技术支持,具体步骤如下:
步骤1,将获取的飞行计划信息中的航空器的关键航路点、经过的扇区编号等信息,进出扇区的高度限制、速度限制等信息;扇区的航空器的雷达航迹历史数据信息作为输入;
步骤2,通过统计分析的方法针对步骤1获取的扇区雷达航迹历史数据按照航空器类型、经过的扇区、不同时间段以及不同的飞行方向(东单西双)进行提取分类,根据如下公式计算爬升下降率的均值和方差:
式中,ROCijk表示机型为i,方向为j,时段为k爬升率的值;
根据计算得到均值与方差,将爬升率的不确定性表述为一个区间,如下所示:
[mean_ROCijk-st_ROCijk,mean_ROCijk+st_ROCijk]
步骤3,基于步骤1中获取的扇区空域的限制以及扇区移交协议中移交点的限制信息,计算管制扇区移交协议的限制信息(高度范围以及水平范围)与制扇区空域的范围限制的交集,表示为Q;根据高度层的设置规则(东单西双,300米一个高度层),在Q中选取可用的进/出扇区指令高度,表示如下:
CFL_set_in={in_CFLijk|in_CFLijk满足扇区空域限制以及进扇区移交协议限制}
CFL_set_out={out_CFLijk|out_CFLijk满足扇区空域限制以及出扇区移交协议限制}
步骤4:根据如下所示的航空器爬升率(或下降率)公式计算航空器的爬升率:
式中,T表示大气温度,单位为K;ΔT表示温度差,单位为K;VTAS表示真空速,由地速和风速的差可以获得;m表示航空器质量;g0=9.80665m/s2表示重力加速度。f(M)表示能量分享因子,其值按照BADA用户手册中提供的不同爬升段的能量分享因子的计算公式,分别计算能量分享因子;
步骤5:根据步骤4计算得到的爬升率,可以计算出高度增量,计算公式如下:
Δh=ROCD·Δt
根据航空器当前高度、进/出扇区的目标指令高度集CFL_set_in或者CFL_set_out,可以得到航空器可能的高度飞行剖面集合,计算结果如图2所示,图中蓝色线条表示不确定爬升率情况下可能的高度剖面集合,红色表示爬升率为1500ft/min的高度剖面(由于说明书附图只能是灰度图,看不出颜色,特此说明);
步骤6:根据步骤5,计算航空器的潜在的运动区域以及航空器潜在运动区域的与其他航空器潜在运动区域、特殊空域、飞行受限区的交集,计算公式为:
Q∩=STP1∩STP2…∩STPm
其中,STPi,i=1,2,…m表示潜在运动区域;
计算潜在冲突概率,计算公式为:
其中p(xu,yu,hu,tu)表示航空器位于坐标(xu,yu,hu,tu)的概率,u∈Q∩表示潜在冲突区域内的任意点,xu,yu,hu,tu分别表示航空器在任意点u的横轴的坐标值、纵轴的坐标值、竖轴的坐标值,以及过u的时间值,n表示相交区域内离散点个数;根据设定的阈值,如果潜在冲突概率大于阈值,则判定存在冲突,阈值一般取冲突的概率值,如取值0.5;
根据冲突探测结果,基于冲突发生的概率对CFL_set_in和CFL_set_out进行分类,具体包括三类:1)当发生冲突的概率为零时,表示无冲突,对应的指令高度值是安全的,在系统设计中可用绿色表示;2)当发生冲突的概率总冲突概率大于0但小于规定的阈值,存在潜在冲突风险,表示指令高度潜在安全风险;在系统设计中可用黄色表示3)总冲突概率大于规定的阈值,表示高冲突风险,指令高度值不安全,不可用,在系统设计中可用红色表示;
步骤7:针对步骤6中存在潜在冲突的指令高度剖面,进行冲突解脱,具体地,首先基于速度控制实现冲突解脱,设当前速度为V,调速方式为然后计算水平和垂直方向上航空器之间的距离,如果满足最小间隔,则实现冲突解脱;如果不满足最小间隔,则针对水平方向有冲突且不能通过调速解决的情况,通过偏航的方式实现冲突解脱,如图3所示,图中P1,P2,P3表示可偏航航路点,红色点表示目标航路点(由于附图只能是灰度图,看不出颜色,特此说明),允许的最大偏航角度为30度、45度以及60度;允许的转弯时间计算公式如下
式子中DLmax最大偏航距离,trangle转弯角度,一般取值30度、45度以及60度;
针对垂直方向有冲突且不能通过调速来解决的,通过调整高度实现冲突解脱;具体方法是:如果需要到达的指令高度高于航空器当前的位置,则通过上升一m个高度层,H的值由QCFL的值以及高度层划设规则计算得到,需要满足H∈QCFL且是300的整数倍;如果需要到达的指令高度低于航空器当前的位置,则通过下降m个高度层,m的获取方式同上;
针对存在冲突的航迹实现冲突解脱后,需进一步进行冲突探测;如果存在冲突,重复上述步骤,直至无冲突。
步骤8,重复步骤1~步骤7,直至管制扇区内各航空器具有最优化的指令高度值。
具体实现中,本发明还提供一种计算机存储介质,其中,该计算机存储介质可存储有程序,该程序执行时可包括本发明提供的一种不确定爬升率的航空器指令高度规划方法的各实施例中的部分或全部步骤。所述的存储介质可为磁碟、光盘、只读存储记忆体(read-only memory,ROM)或随机存储记忆体(random access memory,RAM)等。
本领域的技术人员可以清楚地了解到本发明实施例中的技术可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明实施例中的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
本发明提供了一种不确定爬升率的航空器指令高度规划方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (10)
1.一种不确定爬升率的航空器指令高度规划方法,其特征在于,包括如下步骤:
步骤1,获取飞行计划信息中的航空器的关键航路点、经过的扇区编号、进出扇区的高度限制、速度限制;获取指定扇区的航空器的雷达航迹历史数据;
步骤2,通过统计分析的方法,得到不同机型在不同方向、不同时段的历史爬升率的均值和方差,以均值和方差得到不同机型在不同方向不同时段的爬升率范围;
步骤3,根据扇区空域的限制以及扇区移交协议中移交点的限制信息,选取能够用于指令高度设置的值,包括进扇区的指令高度集合和出扇区的指令高度集合;
步骤4:计算航空器爬升率和下降率;
步骤5:基于航空器爬升率和下降率,根据航空器当前高度、进扇区的指令高度集合和出扇区的指令高度集合,得到航空器可能的高度飞行剖面集合Trijk;
步骤6:计算扇区内各航空器潜在的运动区域,进行冲突探测识别;根据冲突探测识别结果,基于冲突发生的概率对进扇区的指令高度集合和出扇区的指令高度集合进行分类;
步骤7:针对存在潜在冲突风险指令高度的值,基于调整高度、偏航和调整速度进行冲突解脱;
步骤8,重复步骤1~步骤7,直至管制扇区内各航空器具有无冲突的指令高度值。
3.根据权利要求2所述的方法,其特征在于,步骤3包括:
步骤3-1:输入管制扇区空域的范围限制,包括扇区的编号、长、宽、高
步骤3-2:根据步骤3-1的输入,基于集合交集计算方法,计算管制扇区空间区域与管制扇区移交协议的限制区域的重叠区域,计算得到的交集记为QCFL;
步骤3-3:根据高度层的设置规则,在交集QCFL中选取进扇区的指令高度集合CFL_set_in和出扇区的指令高度集合CFL_set_out,表示如下:CFL_set_in={in_CFLijk|in_CFLijk满足扇区空域限制以及进扇区移交协议限制}CFL_set_out={out_CFLijk|out_CFLijk满足扇区空域限制以及出扇区移交协议限制}
其中,in_CFLijk表示机型为i,方向为j,时段为k,且在交集QCFL中根据高度层的设置规则选取的进扇区的指令高度,所有的指令高度in_CFLijk组成进扇区的指令高度集合CFL_set_in;out_CFLijk表示机型为i,方向为j,时段为k,且在交集QCFL中根据高度层的设置规则选取的出扇区的指令高度,所有的指令高度out_CFLijk组成进扇区的指令高度集合CFL_set_out。
5.根据权利要求4所述的方法,其特征在于,步骤4中,按照BADA用户手册中提供的不同爬升段的能量分享因子的计算公式,分别计算能量分享因子,具体包括:
步骤4-1,根据BADA用户手册,计算加速爬升段的能量分享因子;
步骤4-2,根据BADA用户手册,计算等表速爬升段的能量分享因子;
步骤4-3,根据BADA用户手册,计算等马赫爬升段的能量分享因子。
6.根据权利要求5所述的方法,其特征在于,步骤5包括:根据爬升率或下降率计算高度增量,计算公式如下:
Δh=ROCD·Δt
式中,Δt是系统计算采样时间周期;
根据航空器当前高度、进扇区的指令高度集合和出扇区的指令高度集合,根据步骤4爬升率计算公式,基于高度增量,计算得到航空器可能的高度飞行剖面集合Trijk。
7.根据权利要求6所述的方法,其特征在于,步骤6包括:
步骤6-1:根据航空器当前的运行状态,具体包括通过雷达获取的当前的位置、速度信息,推算未来一段时间航空可能的位置分布,位置分布的区域即为航空器的潜在的运动区域;
步骤6-2:计算航空器潜在运动区域与其他航空器潜在运动区域、特殊空域、飞行受限区的交集Q,计算公式为:
Q=STP1∩STP2…∩STPm
其中,STPi表示第i个航空器潜在运动区域、特殊空域、飞行限制空域,i=1,2,…m;
步骤6-3:计算潜在冲突概率P:
其中p(xu,yu,hu,tu)表示航空器位于坐标(xu,yu,hu,tu)的概率,u∈Q表示潜在冲突区域内的任意点,xu,yu,hu,tu分别表示航空器在任意点u的横轴的坐标值、纵轴的坐标值、竖轴的坐标值,以及过u的时间值,n表示相交区域内离散点个数;
步骤6-4:根据设定的阈值,如果潜在冲突概率大于阈值,则判定存在冲突。
8.根据权利要求7所述的方法,其特征在于,步骤6还包括:根据步骤6-4得到的结果,对进扇区的指令高度集合CFL_set_in和出扇区的指令高度集合CFL_set_out进行分类,具体包括如下三类:
当潜在冲突概率为零时,表示无冲突,对应的指令高度值是安全的;
当潜在冲突概率大于0但小于规定的阈值,存在潜在冲突风险,表示指令高度潜在安全风险;
潜在冲突概率大于规定的阈值,表示高冲突风险,指令高度值不安全,不可用。
9.根据权利要求8所述的方法,其特征在于,步骤7包括:
步骤7-1:基于速度控制实现冲突解脱,设当前速度为V,调速方式为然后计算水平和垂直方向上航空器之间的距离,如果满足飞行规则规定的最小安全间隔,则实现冲突解脱;如果不满足最小间隔,则执行步骤7-2和步骤7-3;
步骤7-2:针对水平方向有冲突且不能通过调速解决的情况,通过偏航的方式实现冲突解脱,允许的最大偏航角度为30度、45度和60度;允许的转弯时间T1计算公式如下:
其中DLmax表示最大偏航距离,trangle表示转弯角度;
步骤7-3:针对垂直方向有冲突且不能通过调速来解决的,通过调整高度实现冲突解脱;
步骤7-4:针对存在冲突的航迹实现冲突解脱后,需执行步骤6-1~步骤6-4,进一步进行冲突探测;如果存在冲突,重复步骤7-1~步骤7-4,直至无冲突。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210720445.1A CN115206135B (zh) | 2022-06-16 | 2022-06-16 | 一种不确定爬升率的航空器指令高度规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210720445.1A CN115206135B (zh) | 2022-06-16 | 2022-06-16 | 一种不确定爬升率的航空器指令高度规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115206135A true CN115206135A (zh) | 2022-10-18 |
CN115206135B CN115206135B (zh) | 2024-02-13 |
Family
ID=83578723
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210720445.1A Active CN115206135B (zh) | 2022-06-16 | 2022-06-16 | 一种不确定爬升率的航空器指令高度规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115206135B (zh) |
Citations (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007004252A (ja) * | 2005-06-21 | 2007-01-11 | Oki Electric Ind Co Ltd | 航空路管制用管制卓における順序・間隔付けヒューマンインタフェース装置 |
US20090005960A1 (en) * | 2005-12-23 | 2009-01-01 | Alison Laura Udal Roberts | Air Traffic Control |
US20090012660A1 (en) * | 2005-12-23 | 2009-01-08 | Nats (En Route) Public Limited Company | Air Traffic Control |
CN101692315A (zh) * | 2009-09-25 | 2010-04-07 | 民航总局空管局技术中心 | 基于实时雷达数据的高精度4d飞机航迹分析方法 |
WO2013049875A1 (en) * | 2011-10-04 | 2013-04-11 | Thales Australia Limited | Portable device to control simulated aircraft in air traffic control training system |
US20130221164A1 (en) * | 2012-02-23 | 2013-08-29 | Ge Aviation Systems Llc | Method for flying an aircraft along a flight path |
CN105825720A (zh) * | 2016-04-11 | 2016-08-03 | 成都民航空管科技发展有限公司 | 一种有效减少短期飞行冲突告警虚警的优化方法 |
CN105956790A (zh) * | 2016-05-24 | 2016-09-21 | 南京航空航天大学 | 低空飞行态势安全性评估指标及其评估方法 |
CN108877298A (zh) * | 2018-07-25 | 2018-11-23 | 哈尔滨工业大学 | 一种基于四维航迹运行的航班计划安全性判断预警方法 |
CN110473431A (zh) * | 2019-09-17 | 2019-11-19 | 中国民用航空飞行学院 | 一种进近管制空域高度限制下的冲突预测方法 |
CN110930770A (zh) * | 2019-11-06 | 2020-03-27 | 南京莱斯信息技术股份有限公司 | 一种基于管制意图和飞机性能模型的四维航迹预测方法 |
CN111583724A (zh) * | 2020-05-08 | 2020-08-25 | 中国电子科技集团公司第二十八研究所 | 一种面向四维航迹运行的预战术阶段间隔管理方法 |
CN111613096A (zh) * | 2020-06-04 | 2020-09-01 | 成都民航空管科技发展有限公司 | 一种基于atc系统的cfl指令预先告警方法和系统 |
CN111640331A (zh) * | 2020-06-09 | 2020-09-08 | 成都民航空管科技发展有限公司 | 一种基于多高度层的cfl保护层计算方法及系统 |
CN112037583A (zh) * | 2020-09-02 | 2020-12-04 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN112330133A (zh) * | 2020-10-30 | 2021-02-05 | 南京航空航天大学 | 一种空中交通管制员培训效果智能评价方法及装置 |
CN112927562A (zh) * | 2021-01-28 | 2021-06-08 | 中国电子科技集团公司第二十八研究所 | 一种基于指定高度限制的航空器飞行高度剖面计算方法 |
CN113269991A (zh) * | 2021-04-22 | 2021-08-17 | 南京莱斯信息技术股份有限公司 | 基于实时航迹与管制意图的空中交通中期冲突探测方法 |
CN113470438A (zh) * | 2021-06-30 | 2021-10-01 | 中国电子科技集团公司第二十八研究所 | 一种基于逻辑时序推演仿真的无冲突航班轨迹生成方法 |
CN114115304A (zh) * | 2021-10-26 | 2022-03-01 | 南京航空航天大学 | 一种航空器四维爬升航迹规划方法及系统 |
CN114283624A (zh) * | 2021-12-21 | 2022-04-05 | 中国人民解放军93209部队 | 一种基于空地航迹信息共享的飞机位置预测方法 |
-
2022
- 2022-06-16 CN CN202210720445.1A patent/CN115206135B/zh active Active
Patent Citations (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007004252A (ja) * | 2005-06-21 | 2007-01-11 | Oki Electric Ind Co Ltd | 航空路管制用管制卓における順序・間隔付けヒューマンインタフェース装置 |
US20090005960A1 (en) * | 2005-12-23 | 2009-01-01 | Alison Laura Udal Roberts | Air Traffic Control |
US20090012660A1 (en) * | 2005-12-23 | 2009-01-08 | Nats (En Route) Public Limited Company | Air Traffic Control |
CN101692315A (zh) * | 2009-09-25 | 2010-04-07 | 民航总局空管局技术中心 | 基于实时雷达数据的高精度4d飞机航迹分析方法 |
WO2013049875A1 (en) * | 2011-10-04 | 2013-04-11 | Thales Australia Limited | Portable device to control simulated aircraft in air traffic control training system |
US20130221164A1 (en) * | 2012-02-23 | 2013-08-29 | Ge Aviation Systems Llc | Method for flying an aircraft along a flight path |
CN105825720A (zh) * | 2016-04-11 | 2016-08-03 | 成都民航空管科技发展有限公司 | 一种有效减少短期飞行冲突告警虚警的优化方法 |
CN105956790A (zh) * | 2016-05-24 | 2016-09-21 | 南京航空航天大学 | 低空飞行态势安全性评估指标及其评估方法 |
CN108877298A (zh) * | 2018-07-25 | 2018-11-23 | 哈尔滨工业大学 | 一种基于四维航迹运行的航班计划安全性判断预警方法 |
CN110473431A (zh) * | 2019-09-17 | 2019-11-19 | 中国民用航空飞行学院 | 一种进近管制空域高度限制下的冲突预测方法 |
CN110930770A (zh) * | 2019-11-06 | 2020-03-27 | 南京莱斯信息技术股份有限公司 | 一种基于管制意图和飞机性能模型的四维航迹预测方法 |
CN111583724A (zh) * | 2020-05-08 | 2020-08-25 | 中国电子科技集团公司第二十八研究所 | 一种面向四维航迹运行的预战术阶段间隔管理方法 |
CN111613096A (zh) * | 2020-06-04 | 2020-09-01 | 成都民航空管科技发展有限公司 | 一种基于atc系统的cfl指令预先告警方法和系统 |
CN111640331A (zh) * | 2020-06-09 | 2020-09-08 | 成都民航空管科技发展有限公司 | 一种基于多高度层的cfl保护层计算方法及系统 |
CN112037583A (zh) * | 2020-09-02 | 2020-12-04 | 中国航空无线电电子研究所 | 一种航空器高度层变更程序优化引导方法 |
CN112330133A (zh) * | 2020-10-30 | 2021-02-05 | 南京航空航天大学 | 一种空中交通管制员培训效果智能评价方法及装置 |
CN112927562A (zh) * | 2021-01-28 | 2021-06-08 | 中国电子科技集团公司第二十八研究所 | 一种基于指定高度限制的航空器飞行高度剖面计算方法 |
CN113269991A (zh) * | 2021-04-22 | 2021-08-17 | 南京莱斯信息技术股份有限公司 | 基于实时航迹与管制意图的空中交通中期冲突探测方法 |
CN113470438A (zh) * | 2021-06-30 | 2021-10-01 | 中国电子科技集团公司第二十八研究所 | 一种基于逻辑时序推演仿真的无冲突航班轨迹生成方法 |
CN114115304A (zh) * | 2021-10-26 | 2022-03-01 | 南京航空航天大学 | 一种航空器四维爬升航迹规划方法及系统 |
CN114283624A (zh) * | 2021-12-21 | 2022-04-05 | 中国人民解放军93209部队 | 一种基于空地航迹信息共享的飞机位置预测方法 |
Non-Patent Citations (9)
Title |
---|
PASINI, A等: "Direct route optimization for air traffic management improvement", 2018 IEEE/AIAA 37TH DIGITAL AVIONICS SYSTEMS CONFERENCE (DASC), pages 1 - 9 * |
YANG ZHANG等: "A Methodology for 4D Trajectory Optimization of Arrival Aircraft in Trajectory Based Operation", 2018 13TH WORLD CONGRESS ON INTELLIGENT CONTROL AND AUTOMATION (WCICA), pages 1284 - 1289 * |
刘岩等: "终端区空地数字化管制技术研究与系统设计", 信息化研究, vol. 43, no. 4, pages 13 - 18 * |
周建;王莉莉;AHMED RAHMANI;刘昕;: "分布式多agent系统在飞行冲突解脱中的应用", 广西师范大学学报(自然科学版), no. 03, pages 16 - 22 * |
张军峰;葛腾腾;陈强;王菲;: "离场航空器四维航迹预测及不确定性分析", 西南交通大学学报, no. 04, pages 800 - 806 * |
揭东等: "无人机冲突探测及解脱策略关键技术研究", 武汉理工大学学报(交通科学与工程版), vol. 42, no. 5, pages 776 - 782 * |
潘卫军;况金宏;王文博;李丘;: "程序管制下考虑高度层穿越的航路容量评估模型", 科学技术与工程, no. 01, 8 January 2015 (2015-01-08), pages 320 - 324 * |
王凯;隋东;胡京;: "基于规则库的航路飞行冲突探测与解脱", 哈尔滨商业大学学报(自然科学版), no. 01, pages 117 - 124 * |
陈子伟: "基于连续下降运行的航空器冲突探测研究", 中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑), pages 031 - 168 * |
Also Published As
Publication number | Publication date |
---|---|
CN115206135B (zh) | 2024-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110727287B (zh) | 用于确定爬升剖面的方法和系统以及预测爬升剖面的方法 | |
CN110930770B (zh) | 一种基于管制意图和飞机性能模型的四维航迹预测方法 | |
CN110243359B (zh) | 基于低空风预测模型的安全航迹规划方法 | |
JP6033549B2 (ja) | 航空機の飛行経路を予測するためのデータの供給 | |
CN110059863B (zh) | 一种基于所需到达时间的航空器四维航迹优化方法 | |
US9196165B2 (en) | Assessing feasibility of an aircraft trajectory | |
CN111583724A (zh) | 一种面向四维航迹运行的预战术阶段间隔管理方法 | |
CN114078339B (zh) | 基于飞行趋势推理的时空概率分布网格飞行冲突检测方法 | |
CN110727286B (zh) | 用于确定下降剖面的方法和系统 | |
CN112927562B (zh) | 一种基于指定高度限制的航空器飞行高度剖面计算方法 | |
CN112506227B (zh) | 一种用于民机全发失效迫降的辅助驾驶系统和方法 | |
JP2020077387A (ja) | 垂直飛行経路の最適化 | |
CN113593308A (zh) | 一种民航飞机智能化进场方法 | |
CN112669651A (zh) | 一种基于飞行动态情报中eet值修正预测过点时间的方法 | |
CN114969978B (zh) | 一种基于改进灰狼优化算法的飞机海上平台回收调度方法 | |
Çeçen et al. | Conflict-free en-route operations with horizontal resolution manoeuvers using a heuristic algorithm | |
Cai et al. | Mixed integer nonlinear programming for aircraft conflict avoidance by applying velocity and altitude changes | |
Robinso et al. | A concurrent sequencing and deconfliction algorithm for terminal area air traffic control | |
CN114637327A (zh) | 基于深度策略性梯度强化学习的在线轨迹生成制导方法 | |
Rosenow et al. | Multiple Aircraft in a multi-criteria Trajectory Optimization | |
CN115206135A (zh) | 一种不确定爬升率的航空器指令高度规划方法 | |
Vormer et al. | Optimization of flexible approach trajectories using a genetic algorithm | |
CN114721423B (zh) | 考虑避碰约束的多无人机协同到达预置目标的分配方法 | |
CN112258896A (zh) | 基于航迹的无人机融合空域运行方法 | |
EP3974937B1 (en) | Method and system for determining a descent profile |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |