CN111142555A - 一种基于碰撞风险的机场无人机管控区域划设方法 - Google Patents

一种基于碰撞风险的机场无人机管控区域划设方法 Download PDF

Info

Publication number
CN111142555A
CN111142555A CN201911306980.7A CN201911306980A CN111142555A CN 111142555 A CN111142555 A CN 111142555A CN 201911306980 A CN201911306980 A CN 201911306980A CN 111142555 A CN111142555 A CN 111142555A
Authority
CN
China
Prior art keywords
area
vertical surface
lsf
control area
unmanned aerial
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
Application number
CN201911306980.7A
Other languages
English (en)
Other versions
CN111142555B (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.)
Second Research Institute of CAAC
Original Assignee
Second Research Institute of CAAC
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 Second Research Institute of CAAC filed Critical Second Research Institute of CAAC
Priority to CN201911306980.7A priority Critical patent/CN111142555B/zh
Publication of CN111142555A publication Critical patent/CN111142555A/zh
Application granted granted Critical
Publication of CN111142555B publication Critical patent/CN111142555B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Traffic Control Systems (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提供了一种基于碰撞风险的机场无人机管控区域划设方法,包括如下步骤:S1根据碰撞风险概率e,划设管控区内核区;S2基于管控区内核区,划设与其对应的管控区缓冲区;S3划设无人机管控区,无人机管控区是管控区内核区和管控区缓冲区的并集。本发明提供的方法是国内首次以无人机为对象的机场周边飞行区域分类划设方法,利用该方法划设出了管控区内核区、管控区缓冲区、以及无人机的可飞区,为无人机相关管理规定提供技术支持。确保在不影响民航飞机运行安全与效率的前提下,实现对无人机的精细化运行管理。

Description

一种基于碰撞风险的机场无人机管控区域划设方法
技术领域
本发明涉及空中交通管制技术领域,尤其涉及无人机空中运行管理。
背景技术
无人机是由控制站管理(包括远程操纵或自主飞行)的航空器,也称远程驾 驶航空器。由于无人机具有机动性好、生存能力高、过载能力强等特点,在物流 运输、电力巡线、农林植保、抢险救灾、地理测绘等领域拥有广阔的应用前景。
然而,在近几年民用无人机市场规模不断扩大的同时,由于无人机易制造、易获取、使用者广泛,由此导致无人机运行中潜藏的风险也多种多样,其中较为典型的风险便是干扰民航机场航班正常起降。2014年7月,正准备降落在伦敦希思罗机场的民航客机险些与黑飞无人机相撞,事故调查,该事件被当局定性为“相撞风险严重”,等级“A”;2016年4月,一架载有132名乘客和5名机组人员的民航客机,在由日内瓦飞往伦敦的途中与一架遥控无人机相撞;而据统计,仅2015-2016两年时间内,国内民用无人机扰飞事件就高达27起。由于目前绝大部分无人机在超低空飞行区域运行,民航航班在超低空空域运行时对应航班起降阶段,即此时民航航班处于机场周边超低空飞行区域,因此,为防止无人机干扰民航航班正常起降,相关部门对机场周边超低空飞行区域的无人机运行实施了严格的管控。
最初,我国绝大部分地区根据民航局颁布的《民用机场运行安全管理规定》(CCAR-140) 中第一百六十六条所规定的机场周边障碍物限制范围,将机场周边包含机场障碍物限制面在内的距机场跑道中心线两侧各10公里,跑道端外20公里的区域,划设为无人机管控区域。
2017年5月,民航局根据国际民航组织《国际民用航空公约附件14-机场》中对障碍物限制面的相关规定,陆续公布我国运输机场障碍物限制面保护范围,该范围由如图1所示的12个坐标点依次连接而成,其中圆弧部分的半径为7070米,随后多地将该范围(障碍物限制面保护范围)上方空域作为民用轻型无人机管控空域。该工作提出了无人机管控区域容差缓冲区的概念,容差缓冲区是机场障碍物限制面向外延展一定距离对应的区域。诚然,在机场周边超低空飞行区域划设无人机管控区,对无人机在机场周边的运行进行限制,是保障民航航班安全运行的必要措施。但若机场周边针对无人机的管控区域划设过大,也会使得无人机可飞区域减少,从而对无人机产业发展造成一定的阻碍。目前的无人机管控区域划设缺乏科学性,主要体现如下:
首先,机场周边无人机管控区域划设参照的《民用机场运行安全管理规定》所规定的机场周边障碍物限制范围,是针对静止障碍物作出的限制,而目前参照的国际民航组织《国际民用航空公约附件14-机场》所划定的包括内水平面、锥形面等障碍物限制面及其保护范围 (含容差缓冲区),同样是针对静止障碍物。无人机作为移动物体,套用上述针对静止障碍物所划定的限制范围或保护范围对无人机飞行活动进行管控,显然是不尽合理的。
其次,上述两类方案是针对所有跑道的固定方案,欠缺对不同跑道运行特性的适应,理由如下:一,不同的跑道对应的进离场程序截然不同,用固定方案去适应不同的跑道未必科学合理,无人机管控空域的划设应该适应具体的飞行程序。例如就目前方案而言,机场障碍物限制面是跑道两侧对称的区域,然而航空器起降阶段未必是沿着跑道方向,甚至有可能是明显偏向于跑道一侧,无人机管控区域应该结合实际运行情况划设,未必对称;对称的无人机管控区域对于上述跑道来说,可能会导致管控区的范围过大,从而限制无人机飞行区域。二,就目前方案而言,所有跑道对应的容差缓冲区的圆弧半径均为7070米,该固定的划设缺少相关依据。事实上,根据《国际民用航空公约附件14-机场》规定,不同类型跑道对应的内水平面和锥形面大小不一,导致不同类型跑道的障碍物限制面不同。因此,不同类型跑道在划设不同障碍物限制面的同时却又划设相同的容差缓冲区圆弧半径,缺乏科学依据。
事实上,无人机作为移动物体,其管控区域的划设应客观、充分的考虑无人机与航空器的碰撞风险。碰撞风险评估是民航保障安全底线的最重要工作之一,风险本质上是随机发生的,不可能根本的杜绝风险发生的可能性。为此国际民航组织对空中交通碰撞风险提出了安全目标水平的概念(每对飞机发生碰撞频率为5×10-9次/飞行小时),不同的学者也尝试着从各个方面度量碰撞风险,旨在合理的判断空中交通系统是否满足安全目标水平的要求、确定航空器间安全间隔等。
碰撞风险评估的核心在于飞行误差的确立,通常涉及到纵向误差、侧向误差、竖向误差三个维度,现有方法侧重于宏观分析,对于某些微观场景精度不足。例如,航空器飞行轨迹目前多采取直线模型描述,在航空器转弯的时候该模型将不再适用,合理的航空器飞行轨迹曲线模型是进行碰撞风险评估的重大改进,特别是应用在航空器进离场阶段。此外,大量研究假设航空器的飞行误差在一段航段内满足固定的分布,方便碰撞风险评估的计算的同时对风险评估的精度有一定的损失。更为精确的航空器飞行误差分布应该与航空器飞行轨迹中对应的空间位置有关,在分析航空器爬升降落等微观场景时尤为重要。
综上,由于目前专门针对无人机的机场周边管控区域划设的法规还处于缺失状态,亟需一套聚焦更为精确的航空器飞行误差分布,基于碰撞风险的机场周边无人机管控区域划设方法,为相关法规的制定提供技术支持,在既保障民航航班安全起降的同时,尽可能缩小机场周边无人机管控区域范围,为机场周边运行的无人机提供更多可飞区域,从而进一步促进无人机产业的发展。
发明内容
本发明的目的是:提供一种机场周边飞行区域划设无人机管控区域的方法,该方法是以无人机为对象的,以每条跑道进离场程序特点为出发点,基于碰撞风险在机场周边飞行区域分类划设的方法,为无人机能够飞起来、实现无人机各项功能奠定基础,为实现对无人机的精细化运行管理提供技术支持。
为确保在不影响民航飞机运行安全与效率的前提下,
一种基于碰撞风险的机场无人机管控区域划设方法,包括如下步骤:
S1根据碰撞风险概率e,划设无人机管控区内核区,包括:
S101划设每条跑道管控区内核区
划设跑道离场程序管控区内核区;
划设跑道进场程序管控区内核区;
划设每条跑道管控区内核区,即为跑道离场程序管控区内核区和跑道进场程序管控区内核区的并集;
S102划设机场管控区内核区
机场管控区内核区是每条跑道管控区内核区的并集;
S2基于管控区内核区,划设无人机管控区缓冲区
S201划设每条跑道管控区缓冲区,其为每条跑道管控区内核区水平向外延伸不少于ρ3;ρ3为意图控制无人机到无人机被反制落地期间飞行的最大水平距离;
S202划设机场管控区缓冲区,即为每条跑道管控区缓冲区的并集减去管控区内核区;
S3划设无人机管控区
无人机管控区是管控区内核区和管控区缓冲区的并集;
S4划设无人机可飞区
划设所述无人机管控区的补集,即为无人机可飞区。
优选的,所述S101还包括步骤如下,
所述跑道离场程序管控区内核区nfzd为离场航段曲线
Figure BDA00023221664900000412
对应的后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面nplane(s2)围成的包含
Figure BDA0002322166490000042
的三维空域;
其中:离场航段曲线
Figure BDA0002322166490000043
的离场水平投影曲线L(s),s为离场水平投影曲线L(s)的弧长参数,s2<s1
nplane(s2)是离场水平投影曲线L(s)在点s2处的法平面,参数s2=infCS(s1),其中CS(s1)为满足约束条件:
Figure BDA0002322166490000044
的s组成的集合,inf为下确界函数;
ssf(s1)为与离场航段曲线
Figure BDA0002322166490000045
平行的曲面;
Figure BDA0002322166490000046
后侧竖直面bsf(s1)是曲线L(s)在起点处的法平面,bsf(s1)=nplane(0);
左侧竖直面lsf(s1)是位于曲线L(s)左侧,且距L(s)上每一点距离均相等的竖直面,即以曲线
Figure BDA0002322166490000047
为底边的竖直面,
Figure BDA0002322166490000048
右侧竖直面rsf(s1)为位于曲线L(s)右侧,且距L(s)上每一点距离均相等的竖直面,即以曲线
Figure BDA0002322166490000049
为底边的竖直面,
Figure BDA00023221664900000410
dz是离场航段曲线与底面ssf(s1)之间的距离;
dl是离场航段曲线与其对应的左侧竖直面之间的距离
dr是离场航段曲线与其对应的右侧竖直面rsf(s0)之间的距离;
Figure BDA00023221664900000411
为单位法向量;
区域Ω是高h以下的三维空域。
优选的,所述S101还包括步骤如下,
所述跑道进场程序管控区内核区nfza为进场航段曲线
Figure BDA0002322166490000051
对应的由后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面 nplane(s2)围成的包含进场航段曲线
Figure BDA0002322166490000052
的三维空域;
其中:进场航段曲线
Figure BDA0002322166490000053
的进场水平投影曲线L(s),s为进场水平投影曲线L(s)的弧长参数,s2<s1
nplane(s2)是离场水平投影曲线L(s)在点s2处的法平面,参数s2=infCS(s1),其中CS(s1)为满足约束条件:
Figure BDA0002322166490000054
的s组成的集合,inf为下确界函数;
ssf(s1)为与进场航段曲线
Figure BDA00023221664900000511
平行的曲面;
Figure BDA0002322166490000055
后侧竖直面bsf(s1)是曲线L(s)在起点处的法平面,bsf(s1)=nplane(0);
左侧竖直面lsf(s1)是位于曲线L(s)左侧,且距L(s)上每一点距离均相等的竖直面,即以曲线
Figure BDA0002322166490000056
为底边的竖直面,
Figure BDA0002322166490000057
右侧竖直面rsf(s1)为位于曲线L(s)右侧,且距L(s)上每一点距离均相等的竖直面,即以曲线
Figure BDA0002322166490000058
为底边的竖直面,
Figure BDA0002322166490000059
dz是进场航段曲线与底面ssf(s1)之间的距离;
dl是进场航段曲线与其对应的左侧竖直面之间的距离
dr是进场航段曲线与其对应的右侧竖直面rsf(s1)之间的距离;
Figure BDA00023221664900000510
为单位法向量;
区域Ω是高h以下的三维空域。
优选的,所述S201还包括步骤如下:
每条跑道离场程序管控区缓冲区界限包括:
内侧竖直面:每条跑道离场程序管控区内核区围成的三维空域的四个竖直面;
外后侧竖直面Bsf(s1):与后侧面竖直面bsf(s1)平行的竖直平面
Figure BDA0002322166490000061
外左侧竖直面Lsf(s1):位于左侧竖直面lsf(s1)左侧,且距lsf(s1)上每一点距离均相等的竖直面,即以曲线
Figure BDA0002322166490000062
为底边的竖直面
Lsf(s1)={(x,y,z)|(x,y,0)∈LL(s),z≥0}
外右侧竖直面Rsf(s1):位于右侧竖直面rsf(s1)右侧,且距rsf(s1)上每一点距离均相等的垂直面,即以曲线
Figure BDA0002322166490000063
为底边的竖直面
Rsf(s1)={(x,y,z)|(x,y,0)∈LR(s),z≥0}
外前侧竖直面Fsf(s1):离场航段水平投影曲线L(s)在点s3处的法平面nplane(s3);
Db是外后侧竖直面Bsf(s1)与后侧面竖直面bsf(s1)之间的距离;
DL是外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离;
DR是外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离。
优选的,所述S201还包括步骤如下:
每条跑道离场程序管控区缓冲区为满足下述约束条件的最小区域:
约束条件1:外后侧竖直面Bsf(s1),外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)分别与后侧竖直面bsf(s1)、左侧竖直面lsf(s1)、右侧竖直面rsf(s1)的距离不小于ρ3,即min{Db,DL,DR}≥ρ3
约束条件2:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s2)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ3
约束条件3:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s4)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ2
其中:
ρ2是无人机被反制达最大高度时对应的水平距离;ρ3为意图控制无人机到无人机被反制落地期间飞行达到的最大水平距离;
参数值s4是满足条件
Figure BDA0002322166490000071
的s的最大值;
h2是无人机被反制落地期间飞行的最大高度;
Figure BDA0002322166490000072
为h2的等高面。
优选的,所述S201还包括步骤如下:
每条跑道进场程序管控区缓冲区界限包括:
内侧竖直面:每条跑道进场程序管控区内核区围成的三维空间对应的竖直面;
外后侧竖直面Bsf(s1):与后侧面竖直面bsf(s1)平行的竖直平面
Figure BDA0002322166490000073
外左侧竖直面Lsf(s1):位于左侧竖直面lsf(s1)左侧,且距lsf(s1)上每一点距离均相等的竖直面,即以曲线
Figure BDA0002322166490000074
为底边的竖直面
Lsf(s1)={(x,y,z)|(x,y,0)∈LL(s),z≥0}
外右侧竖直面Rsf(s1):位于右侧竖直面rsf(s1)右侧,且距rsf(s1)上每一点距离均相等的垂直面,即以曲线
Figure BDA0002322166490000075
为底边的竖直面
Rsf(s1)={(x,y,z)|(x,y,0)∈LR(s),z≥0}
外前侧竖直面Fsf(s1):进场航段水平投影曲线L(s)在点s3处的法平面nplane(s3);
Db是外后侧竖直面Bsf(s1)与后侧竖直面bsf(s1)之间的距离;
dl是进场航段曲线与左侧竖直面lsf(s1)之间的距离;
dr是进场航段曲线与右侧竖直面rsf(s1)之间的距离;
DL是外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离;
DR是外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离;
Figure BDA0002322166490000076
为单位法向量。
优选的,所述S201还包括步骤如下:
每条跑道进场程序管控区缓冲区为满足下述约束条件的最小区域:
约束条件1:外后侧竖直面Bsf(s1),外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)分别与后侧竖直面bsf(s1)、左侧竖直面lsf(s1),右侧竖直面rsf(s1)的距离不小于ρ3,即min{Db,DL,DR}≥ρ3
约束条件2:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s2)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ3
约束条件3:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s4)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ2
其中:
ρ2是无人机被反制达最大高度时对应的水平距离;ρ3为意图控制无人机到无人机被反制落地期间飞行的最大水平距离;
参数值s4是满足条件
Figure BDA0002322166490000081
的s的最大值;
h2是无人机被反制落地期间飞行的最大高度;
Figure BDA0002322166490000082
为h2的等高面。
本发明的有益效果是:
(1)本发明提供的方法是国内首次以无人机为对象的机场周边飞行区域分类划设方法,利用该方法划设出了管控区内核区、管控区缓冲区、以及无人机的可飞区。克服了现有以静止障碍物为划设依据的缺陷;
(2)本发明提供的划设方法以民航航班进离场程序特点为依据,基于碰撞风险划设无人机管控区域,实现了无人机管控区的精确划设;
(3)本发明划设方法大大缩小了机场周边无人机管控区域范围,较现有规定的机场周边无人机管控区区域范围缩小了25%以上。
(4)本发明提供的方法为无人机能够飞起来、实现无人机各项功能奠定基础,为确保在不影响民航飞机运行安全与效率的前提下,实现对民用无人机管控区的精准划设,无人机的精细化运行管理提供技术支持。
(5)本发明是基于碰撞风险的管控区划设,是无人机空中交通运行管理的重要内容,保证了无人机不会干扰有人航空器正常运行,对无人机在机场周边超低空空域的运行进行限制,是保障民航航班安全运行的必要措施,也是无人机空管的一项关键技术,对无人机产业化至关重要。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1为中国民航机场保护范围示意图;
图2为碰撞风险区结构示意图;
图3为离场程序下无人机管控区结构示意图;
图4为无人机飞行特性示意图;
图5为无人机管控区缓冲区划设原理图;
图6为进场程序下无人机管控区结构示意图;
图7为某机场无人机管控区示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不限定本发明。
本发明所考虑的空域为机场周边某一范围内的空域
Figure BDA0002322166490000091
(例如沿机场各跑道中心线两侧各 10千米,跑道端外20千米的区域),重点关注该空域真高h以下的低空空域Ω。所涉及空域包含无人机管控区(RFZ)和无人机可飞区(FFZ)。其中无人机管控区将依据功能的不同分为管控区内核区(NFZ)以及管控区缓冲区(LFZ)。
无人机可飞区(FFZ):民用轻型无人机(同时满足空机重量不超过4千克,最大起飞重量不超过7千克,最大平飞速度不超过100千米/小时,具备符合空域管理要求的空域保持能力和可靠被监视能力的遥控航空器)无需申请,直接进行飞行活动的空域。
无人机管控区(RFZ):无人机需经过有关部门审批之后,方可按照预先提交的飞行计划进行飞行活动的空域,无人机管控区以及无人机可飞区构成了整个机场周边空域Ω,且两者没有交集,即:
FFZ∪RFZ=Ω
Figure BDA0002322166490000101
无人机管控区由管控区内核区和管控区缓冲区组成:
管控区内核区(NFZ):管控区内核区为管控区中间的内核区域,该区域通常情况下是没有无人机飞行的,是管控区划设的重要依据。
管控区内核区NFZ包括所有跑道runway离场程序管控区内核区nfzd与进近程序管控区内核区nfza。即对应每条跑道runway,跑道运行管控区内核区nfzrunway为:
nfzrunway=nfzd∪nfza
机场跑道运行管控区内核区NFZ3为:
Figure BDA0002322166490000102
管控区缓冲区(LFZ):一般划设在机场管控区内核区和可飞区之间。机场周边无人机管控区缓冲区主要目的是留出足够的时间和空间余度对可飞区内发生特情(无人机失控飞向机场管控区内核区)或企图恶意入侵机场管控区内核区的无人机进行干扰反制,或通知机场内航班对该类无人机进行避让,保证航班起降的安全,主要为管控区内核区向外部延伸的空域。管控区内核区以及管控区缓冲区构成无人机管控区的划分,即:
NFZULFZ=RFZ
Figure BDA0002322166490000103
本发明的机场无人机运行管理区域划设方法,包括:S1划设无人机管控区,S2划设无人机可飞区,所述无人机可飞区为无人机管控区的补集。
航段碰撞风险区
如图2所示,航段
Figure BDA0002322166490000104
对应的碰撞风险区rz(s1)为一个包含航段曲线
Figure BDA0002322166490000108
Figure BDA0002322166490000109
在内的三维区域,其边界由以下各个部分组成:
后侧竖直面bsf(s1):即水平投影曲线L(s)在起点处的法平面
bsf(s1)=nplane(0)
左侧竖直面lsf(s1):位于曲线L(s)左侧,且距L(s)上每一点距离均相等的竖直面,即以曲线
Figure BDA0002322166490000106
为底边的竖直面
Figure BDA0002322166490000107
右侧竖直面rsf(s1):位于曲线L(s)右侧,且距L上每一点距离均相等的竖直面。即以曲线
Figure BDA0002322166490000111
为底边的竖直面
Figure BDA0002322166490000112
前侧竖直面fsf(s1):包含曲线L(s)在点L(s1)处法平面
fsf(s1)=nplane(0)
一个底面ssf(s1):垂直于左、右侧面竖直面且与航段曲线
Figure BDA0002322166490000113
平行的曲面
Figure BDA0002322166490000114
碰撞风险区边界要求
假设航班实际运行航迹在点L(s)的侧向,垂向误差εy(s)的分布函数为Fy(s,εy),εz(s)的分布函数为Fz(s,εz),可接受的碰撞概率为e,那么航段
Figure BDA0002322166490000115
对应的碰撞风险区为满足下述条件的最小区域:
条件一:假设该航段运行航班飞机的最大侧宽为2λy,那么左、右侧竖直面参数dl、dr应满足:
Figure BDA0002322166490000116
Figure BDA0002322166490000117
Figure BDA0002322166490000118
P(x)是事件x发生的概率
条件二:假设该航段运行飞机的最大高度为2λz,那么底面参数应满足:
Figure BDA0002322166490000119
P(εz(s)<-dzz)=F(s,-dzz)≤e
P(x)是事件x发生的概率
本发明称满足以上两个条件的参数dl,dr,dz为碰撞风险概率e规避的,将碰撞风险概率e规避的参数dl,dr,dz的集合分别表示为RAl(e),RAr(e),RAz(e),那么碰撞风险区边界对应的参数分别为:
dl=infRAl(e)
dr=infRAr(e)
dz=infRAz(e)
其中:inf代表实数集合的下确界;
dl是航段曲线与左侧竖直面lsf(s1)之间的距离;
dz是航段曲线曲线与底面ssf(s1)之间的距离;
dr是航段曲线与右侧竖直面rsf(s1)之间的距离。
S1划设无人机管控区
运行程序一:跑道离场程序下管控区,如图3所示
空间结构
记dep为所考虑跑道对应的某离场程序,该离场程序对应的航段曲线
Figure BDA0002322166490000121
其水平投影曲线 L(s),那么离场程序管控区内核区nfz由离场程序某航段
Figure BDA0002322166490000122
的碰撞风险区rzd(s1) 确定,其中航段的起点
Figure BDA0002322166490000123
为跑道中线端点(同时认为是起飞离场程序航线起点)。本发明分别记航段曲线
Figure BDA0002322166490000124
的碰撞风险区对应的前侧竖直面fsf(s1),后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),那么该航段对应的管控区内核区nfz(s1)为由航段曲线
Figure BDA0002322166490000125
对应的后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面 nplane(s2)围成的包含航段曲线
Figure BDA0002322166490000126
的三维空域。管控区内核区的空间结构主要由参数值s1,s2确定,参数值s2的选取是内核区nfz(s1)的关键因素。
管控区内核区nfz(s1)约束条件
管控区内核区区域nfz(s1)主要是描述碰撞风险区rzd(s1)在区域Ω内的部分以及其向上向下的延展区域,参数s2应满足相交曲线nplane(s2)∩ssf(s1)不小于h,即:
Figure BDA0002322166490000127
记满足以上约束条件的所有满足条件的s组成的集合为CS(s1),那么s2应设置为:
s2=infCS(s1)
管控区缓冲区1fz划设模型
跑道离场程序管控区rfz是为了保障民航在空域融入无人机飞行时民航航班安全运行的空域,其以限飞内核区nfz为中心,以一定的缓冲距离向外延展得到缓冲区lfz。并保证在保有反制系统的空域结构内,可飞区ffz内的无人机在没有经过飞行申请的情况下无法接近管控区内核区nfz,以保证非合作类无人机与民航航班的运行隔离,杜绝威胁产生。
管控区缓冲区1fz组成
管控区缓冲区lfz(s1)界限由以下几部分组成:
一个顶面:区域Ω的上边界面
内侧竖直面:管控区内核区围成的三维空域的四个竖直面,即后侧竖直面
bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面nplane(s2);
外后侧竖直面Bsf(s1):与内后侧面bsf(s1)平行的竖直平面
Figure BDA0002322166490000131
外左侧竖直面Lsf(s1):位于内左侧竖直面lsf(s1)左侧,且距lsf(s1)上每一点距离均相等的竖直面。即以曲线
Figure BDA0002322166490000132
为底边的竖直面
Lsf(s1)={(x,y,z)|(x,y,0)∈LL(s),z≥0}
外右侧竖直面Rsf(s1):位于内右侧竖直面rsf(s1)右侧,且距rsf(s1)上每一点距离均相等的垂直面,即以曲线
Figure BDA0002322166490000133
为底边的竖直面
Rsf(s1)={(x,y,z)|(x,y,0)∈LR(s),z≥0}
外前侧竖直面Fsf(s1):水平投影曲线L(s)在点L(s3)处法平面nplane(s3)
管控区缓冲区空间几何结构主要由参数Db,DL,DR,s3确定,缓冲区lfz(s1)的划设关键在于给出上述参数,其中,
Db是外后侧竖直面Bsf(s1)与后侧竖直面bsf(s1)之间的距离;
DL是外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离;
DR是外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离;
参数值s3<s1
缓冲区lfz(s1)限制条件
参数Db,DL,Dr的选取主要依赖于轻微型无人机的极限运行性能,s3的选取同时依赖于无人机性能以及民航航班的运行数据。假设轻型无人机最大平飞速度为vh,最大爬升速度为vz,从探测到未经审批的无人机UA进入无人机限制区到干扰反制成功所用时间为Tr,那么,无人机于t=0时在可飞区域边界以最大运行性能冲入管控区,那么其可能在限制区内飞行最远的情况为:t∈[0,Tr]时,无人机仍然以最大性能(垂直速度以及水平速度为vz,vh)飞行;在Tr时刻,无人机失去动力系统,将做最大性能初速度的自由落体运动。
如图4所示,根据物理知识可知在t=Tr时,无人机上升高度Vz×Tr,水平飞行距离ρ1=Vh×Tr,此时无人机所在高度为h1=h+vz×Tr;在经过时间
Figure BDA0002322166490000141
后,无人机UA 达到制高点,此时垂直速度均转换为重力势能,故此时UA所在高度
Figure BDA0002322166490000142
UA在经过时间t3后着陆,此时UA的垂直速度为
Figure BDA0002322166490000143
Figure BDA0002322166490000144
综上所述,无人机UA在整个过程中上升的最大高度为:
Figure BDA0002322166490000145
此时UA飞行的水平距离ρ2=vht2;无人机UA在空中飞行总时间为:
Figure BDA0002322166490000146
故全过程中飞行的最大水平距离ρ3为:
Figure BDA0002322166490000147
记s4<s3为满足条件
Figure BDA0002322166490000148
的最大的s值,那么管控区缓冲区lfz(s1) 为同时满足下述条件的最小区域,如图5所示:
条件一:外后侧竖直面Bsf(s1)与后侧竖直面bsf(s1)之间的距离、外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离、外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离均不小于ρ3,即:
min{Db,DL,DR}≥ρ3
条件二:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s2)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ3条件三:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s4)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ2
参数s1的约束方程
离场程序管控区与某一条航段
Figure BDA0002322166490000151
的碰撞风险区成对应关系,也就是说管控区rfz 的空间几何属性依赖于参数s1的选取,本发明给出参数s1满足的约束。当已知航班飞行误差的分布与无人机极限飞行性能时,参数s2,s4均可由参数s1固定,即可以认为参数s2,s4均为参数s1的函数,可以表示为s2(s1),s4(s1),参数s3为参数s2,s4的函数,故而依然为参数s1的函数。在此时,参数s1应该由下述方程组给出:
s3(s1)=min{s2(s1)+ρ3,s4(s1)+ρ2}
s1=s3(s1)
上述方程的具体形式依赖于飞机飞行误差的分布规律以及无人机的极限飞行性能。
运行程序二:进场程序管控区划设,如图5所示
管控区内核区nfz划设模型
空间结构
记app为所考虑跑道对应的某进场程序,该进场程序对应的航段曲线
Figure BDA0002322166490000152
那么进场程序管控区内核区nfz由进场程序某航段
Figure BDA0002322166490000153
的碰撞风险区rzd(s1)确定,其中航段的起点
Figure BDA0002322166490000154
为跑道中线端点(同时认为是起飞进场程序航线起点)。本发明分别记航段曲线
Figure BDA0002322166490000155
的碰撞风险区对应的前侧竖直面fsf(s1),后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),那么该航段对应的管控区内核区nfz(s1)为由航段曲线
Figure BDA0002322166490000161
对应的后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面nplane(s2)围成的包含航段曲线
Figure BDA0002322166490000164
,s∈[0,s2]的三维空域。管控区内核区的空间结构主要由参数值s1,s2确定,参数值s2的选取是内核区 nfz(s1)的关键因素。
管控区内核区nfz(s1)约束条件
管控区内核区域nfz(s1)主要是描述碰撞风险区rzd(s1)在区域Ω内的部分以及其向上向下的延展区域,参数s2应满足相交曲线nplane(s2)∩ssf(s1)不小于h,即:
Figure BDA0002322166490000162
记满足以上约束条件的所有满足条件的s组成的集合为CS(s1),那么s2应设置为:
s2=infCS(s1)
管控区缓冲区lfz划设模型
跑道进场程序管控区rfz是为了保障民航在空域融入无人机飞行时民航航班安全运行的空域,其以内核区nfz为中心,以一定的缓冲距离向外延展得到缓冲区lfz。并保证在保有反制系统的空域结构内,可飞区ffz内的无人机在没有经过飞行申请的情况下无法接近管控区内核区 nfz,以保证非合作类无人机与民航航班的运行隔离,杜绝威胁产生。
管控区缓冲区组成
管控区缓冲区lfz(s1)界限由以下几部分组成:
一个顶面:区域Ω的上边界面;
内侧竖直面:管控区内核区围成的三维空域的四个竖直面,即后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面nplane(s2);
外后侧竖直面Bsf(s1):与内后侧面bsf(s1)平行的竖直平面
Figure BDA0002322166490000163
外左侧竖直面Lsf(s1):位于左侧竖直面lsf(s1)左侧,且距lsf(s1)上每一点距离均相等的竖直面。即以曲线
Figure BDA0002322166490000171
勾底边的竖直面
Lsf(s1)={(x,y,z)|(x,y,0)∈LL(s),z≥0}
外右侧竖直面Rsf(s1):位于曲线rsf(s1)右侧,且距rsf(s1)上每一点距离均相等的垂直面,即以曲线
Figure BDA0002322166490000172
为底边的竖直面
Rsf(s1)={(x,y,z)|(x,y,0)∈LR(s),z≥0}
外前侧竖直面Fsf(s1):曲线L(s)在点L(s3)处法平面nplane(s3)
管控区缓冲区空间几何结构主要由参数Db,DL,DR,s3确定,缓冲区lfz(s1)的划设关键在于给出上述参数。其中,
Db是外后侧竖直面Bsf(s1)与后侧竖直面bsf(s1)之间的距离;
DL是外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离;
DR是外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离;
参数值s3<s1
缓冲区lfz(s1)限制条件
参数Db,DL,DR选取主要依赖于轻微型无人机的极限运行性能,s3的选取同时依赖于无人机性能以及民航航班的运行数据。假设轻型无人机最大平飞速度为vh,最大爬升速度为vz,从探测到未经审批的无人机UA进入无人机限制区到干扰反制成功所用时间为Tr,那么,无人机于t=0时在可飞区域边界以最大运行性能冲入管控区,那么其可能在限制区内飞行最远的情况为:t∈[0,Tr]时,无人机仍然以最大性能(垂直速度以及水平速度为vz,vh)飞行;在Tr时刻,无人机失去动力系统,将做最大性能初速度的自由落体运动。
根据物理知识可知在t=Tr时,无人机上升高度Vz×Tr,水平飞行距离ρ1=Vh×Tr,此时无人机所在高度为h1=h+vz×Tr;在经过时间
Figure BDA0002322166490000173
后,无人机UA达到制高点,此时垂直速度均转换为重力势能,故此时UA所在高度
Figure BDA0002322166490000174
UA在经过时间t3后着陆,此时UA的垂直速度为
Figure BDA0002322166490000175
Figure BDA0002322166490000176
综上所述,无人机UA在整个过程中上升的最大高度为:
Figure BDA0002322166490000181
此时UA飞行的水平距离ρ2=vht2;无人机UA在空中飞行总时间为:
Figure BDA0002322166490000182
故全过程中飞行的最大水平距离ρ3为:
Figure BDA0002322166490000183
记s4<s3为满足条件
Figure BDA0002322166490000184
的最大的s值,那么管控区缓冲区lfz(s1) 为同时满足下述条件的最小区域:
条件一:外后侧竖直面Bsf(s1)与后侧竖直面bsf(s1)之间的距离、外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离、外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离均不小于ρ3,即:
min{Db,DL,DR)≥ρ3
条件二:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与 nplane(s2)位于左侧竖直面lsf(si),右侧竖直面rsf(s1)的部分子集的距离Df不小于ρ3
条件三:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与 nplane(s4)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离Df不小于ρ2
参数s1的约束方程
通过前述讨论,进场程序管控区与某一条航段
Figure BDA0002322166490000185
,s∈[0,s1]的碰撞风险区成对应关系,也就是说管控区rfz的空间几何属性依赖于参数s1的选取,本发明给出参数s1满足的约束。当已知航班飞行误差的分布与无人机极限飞行性能时,参数s2,s4均可由参数s1固定,即可以认为参数s2,s4均为参数s1的函数,可以表示为s2(s1),s4(s1),参数s3为参数s2,s4的函数,故而依然为参数s1的函数。在此时,参数s1应该由下述方程组给出:
s3(s1)=min{s2(s1)+ρ3,s4(s1)+ρ2}
s1=s3(s1)
某机场实例:
如图6所示结合某机场近一个月的实际飞行数据,其具体实施如下:
某机场特点如下:
所论述的机场为单跑道机场,该跑道对应有一个进场程序和一个离场程序,机场附近均为平地;
进场航段和离场航段对应的参数s1均为10000米,其航段曲线的表达式以及对应的误差分布均已知;
S1根据碰撞风险概率e,划设管控区内核区
管控区内核区主要依据航段运行性能,运行数据经风险评估后确立的。
1.1进场程序
不失一般性,航班在最后进近阶段航段
Figure BDA0002322166490000191
为直线,根据仪表进近图信息可知水平速度 vx=93.05米/秒,vy=0米/秒,vz=4.9米/秒,故可知航空器在进近阶段的下降梯度为
Figure BDA0002322166490000192
故航段
Figure BDA0002322166490000193
满足条件:
x(s)=s
y(s)=0
z(s)=ka×s
对于碰撞风险概率e=5×10-9,根据碰撞风险区边界条件可取dl=dr=283.887,dz= 129.615。
假设h=120米,此刻,管控区内核区nplane(s2)对应的参数s2满足:
z(s2)=h+dz=120+129.615=249.615
由于航段对应的梯度ka=0.0526,故
Figure BDA0002322166490000194
故管控区内核区nfza为由
Figure BDA0002322166490000195
对应的后侧竖直面bsf(10000),左侧竖直面lsf(10000),右侧竖直面rsf(10000),法平面nplane(4740.42)围成的三维空域。
1.2离场程序
部分跑道离场程序对应的航段为曲线(航班起飞的过程中有转弯的情况)。本例中对应的航段
Figure BDA0002322166490000201
为:
Figure BDA0002322166490000202
Figure BDA0002322166490000203
z(s)=0.05s
该航段描述的是起飞2千米后沿着半径为6千米的圆弧转弯角度
Figure BDA0002322166490000204
后继续直线飞行的情况。本例之中,对于碰撞风险概率e=5×10-9,根据碰撞风险区边界条件可取dl=dr=283.887,dz=133.5。
管控区内核区边界nplane(s2)对应的参数s2满足:
z(s2)=h+dz=120+133.5=253.5
由于航段对应的梯度kd=0.05,故
Figure BDA0002322166490000205
故而,管控区内核区nfzd为由航段曲线
Figure BDA0002322166490000206
后侧竖直面bsf(10000),左侧竖直面lsf(10000),右侧竖直面rsf(10000),法平面nplane(5070)围成的三维空域。
S2根据管控区内核区,划设与其对应的管控区缓冲区
管控区缓冲区将结合无人机运行性能以及航班运行数据共同确立。假设允许自由飞行的无人机(UA)的运行性能满足vz=2米/秒,vh=100千米/小时=27.778米/秒,反制系统的探测以及反制响应时Tr=0.5秒,取重力加速度g=9.8米/秒2那么根据物理知识可知在t=Tr时,无人机上升高度h0为:
h0=vz×Tr=2×0.5=1米
水平飞行距离ρ1为:
ρ1=vh×Tr=27.778×0.5=13.889米
此时无人机所在高度h1为:
h1=h0+h=1+120=121米
由于无人在垂直方向为匀加速直线运动,可知当无人机达到最高点是速度为零,此刻由匀加速直线运动公式:
vend=vinit+g×t
可知,无人机爬升到最高点所经历的时间为:
Figure BDA0002322166490000211
此时垂直速度均转换为重力势能,由能量守恒公式:
Figure BDA0002322166490000212
故此时UA所在高度h2为:
Figure BDA0002322166490000213
在无人机从至高点落到地面的过程中,再次利用能量守恒公式可知落地时垂直速度vz3为:
Figure BDA0002322166490000214
无人机从至高点下落至地面所经历的时间t3为:
Figure BDA0002322166490000215
综上所述,无人机UA在整个过程中上升的最大高度为:
Figure BDA0002322166490000221
UA在达到之高点时,UA飞行的水平距离ρ2为:
ρ2=vh(Tr+t2)=19.558米
无人机UA在空中飞行总时间为:
Figure BDA0002322166490000222
将Tr=0.5,vz=2,g=9.8,h=120代入上式可知t=5.678s故全过程中飞行的最大水平距离ρ3为:
Figure BDA0002322166490000223
2.1进场程序管控区缓冲区
进场阶段管控区内核区已经由1.1节讨论清楚;结合无人机空中运行的性能分析,可知由于航段为直线,由模型可知,
Figure BDA0002322166490000224
由管控区条件可知参数s3应为:
s3=min{s23,s32}=max{4894.13,4600.248}=4894.13米
综上所述,管控区缓冲区对应的参数s4=4580.69米,Db=DL=DR=157.71米,s3=4894.13米。
2.2离场程序管控区缓冲区
离场阶段管控区内核区已经由1.2节讨论清楚;结合无人机空中运行的性能分析,可知由于航段在s>s2的情况下为直线,由模型可知,
Figure BDA0002322166490000231
由管控区条件可知参数s3应为:
s4=min{s23,s32}=max{5227.71,5112.638}=5227.71米
综上所述,管控区缓冲区对应的参数s3=5094.08米,Db=DL=DR=157.71米,s3=5227.71米。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。尽管本发明已进行了一定程度的描述,明显地,在不脱离本发明的精神和范围的条件下,可进行各个条件的适当变化。可以理解,本发明不限于所述实施方案,而归于权利要求的范围,其包括所述每个因素的等同替换。

Claims (7)

1.一种基于碰撞风险的机场无人机管控区域划设方法,其特征在于,包括如下步骤:
S1根据碰撞风险概率e,划设无人机管控区内核区,包括:
S101划设每条跑道管控区内核区
划设跑道离场程序管控区内核区;
划设跑道进场程序管控区内核区;
划设每条跑道管控区内核区,即为跑道离场程序管控区内核区和跑道进场程序管控区内核区的并集;
S102划设机场管控区内核区
机场管控区内核区是每条跑道管控区内核区的并集;
S2基于S1划设的管控区内核区,划设无人机管控区缓冲区
S201划设每条跑道管控区缓冲区,其为每条跑道管控区内核区水平向外延伸不少于ρ3;ρ3为意图控制无人机到无人机被反制落地期间飞行的最大水平距离;
S202划设机场管控区缓冲区,即为每条跑道管控区缓冲区的并集减去管控区内核区;
S3划设无人机管控区
无人机管控区是管控区内核区和管控区缓冲区的并集。
2.如权利要求1所述划设方法,其特征在于,所述S101还包括步骤如下,
所述跑道离场程序管控区内核区nfzd为离场航段曲线
Figure FDA0002322166480000011
s∈[0,s1]对应的后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面nplane(s2)围成的包含
Figure FDA0002322166480000012
s∈[0,s2]的三维空域;
其中:离场航段曲线
Figure FDA0002322166480000013
的离场水平投影曲线L(s),s为离场水平投影曲线L(s)的弧长参数,参数值s2<s1
nplane(s2)是离场水平投影曲线L(s)在点s2处的法平面,参数s2=infCS(s1),其中CS(s1)为满足约束条件:
Figure FDA0002322166480000014
的s组成的集合,inf为下确界函数;
ssf(s1)为与离场航段曲线
Figure FDA0002322166480000015
平行的曲面:
Figure FDA0002322166480000016
后侧竖直面bsf(s1)是曲线L(s)在起点处的法平面,bsf(s1)=nplane(0);
左侧竖直面lsf(s1)是位于曲线L(s)左侧,且距L(s)上每一点距离均为dl的竖直面,即以曲线
Figure FDA0002322166480000021
为底边的竖直面,
Figure FDA0002322166480000022
右侧竖直面rsf(s1)为位于曲线L(s)右侧,且距L(s)上每一点距离均为dr的竖直面,即以曲线
Figure FDA0002322166480000023
为底边的竖直面,
Figure FDA0002322166480000024
dz是离场航段曲线与其对应的底面ssf(s1)之间的距离;
dl是离场航段曲线与其对应的左侧竖直面lsf(s1)之间的距离
dr是离场航段曲线与其对应的右侧竖直面rsf(s1)之间的距离;
Figure FDA0002322166480000025
为曲线L(s)的单位法向量;
Ω是真高h以下的三维空域。
3.如权利要求1或2所述划设方法,其特征在于,所述S101还包括步骤如下,
所述跑道进场程序管控区内核区nfza为进场航段曲线
Figure FDA0002322166480000026
对应的由后侧竖直面bsf(s1),左侧竖直面lsf(s1),右侧竖直面rsf(s1),法平面nplane(s2)围成的包含进场航段曲线
Figure FDA0002322166480000027
s∈[0,s2]的三维空域;
其中:进场航段曲线
Figure FDA0002322166480000028
的进场水平投影曲线L(s),s为进场水平投影曲线L(s)的弧长参数,参数值s2<s1
nplane(s2)是离场水平投影曲线L(s)在点s2处的法平面,参数s2=infCS(s1),其中CS(s1)为满足约束条件:
Figure FDA0002322166480000029
的s组成的集合,inf为下确界函数;
ssf(s1)为与进场航段曲线
Figure FDA00023221664800000210
平行的曲面;
Figure FDA00023221664800000211
后侧竖直面bsf(s1)是曲线L(s)在起点处的法平面,bsf(s1)=nplane(0);
左侧竖直面lsf(s1)是位于曲线L(s)左侧,且距L(s)上每一点距离均为dl竖直面,即以曲线
Figure FDA00023221664800000212
为底边的竖直面,
Figure FDA00023221664800000213
右侧竖直面rsf(s1)为位于曲线L(s)右侧,且距L(s)上每一点距离均为dr的竖直面,即以曲线
Figure FDA0002322166480000031
为底边的竖直面,
Figure FDA0002322166480000032
dz是进场航段曲线与底面ssf(s1)之间的距离;
dl是进场航段曲线与其对应的左侧竖直面lsf(s1)之间的距离
dr是进场航段曲线与其对应的右侧竖直面rsf(s1)之间的距离;
Figure FDA0002322166480000033
为水平投影曲线L(s)的单位法向量;
区域Ω是真高h以下的三维空域。
4.如权利要求2所述划设方法,其特征在于,所述S201还包括步骤如下:
每条跑道离场程序管控区缓冲区界限包括:离场航段曲线
Figure FDA0002322166480000034
s∈[0,s1]对应的
内侧竖直面:每条跑道离场程序管控区内核区围成的三维空域对应四个的竖直面;
外后侧竖直面Bsf(s1):与后侧面竖直面bsf(s1)平行的竖直平面
Figure FDA0002322166480000035
外左侧竖直面Lsf(s1):位于左侧竖直面lsf(s1)左侧,且距lsf(s1)上每一点距离均为DL的竖直面,即以曲线
Figure FDA0002322166480000036
为底边的竖直面
Lsf(s1)={(x,y,z)|(x,y,0)∈LL(s),z≥0}
外右侧竖直面Rsf(s1):位于右侧竖直面rsf(s1)右侧,且距rsf(s1)上每一点距离均相等的垂直面,即以曲线
Figure FDA0002322166480000037
为底边的竖直面
Rsf(s1)={(x,y,z)|(x,y,0)∈LR(s),z≥0}
外前侧竖直面Fsf(s1):离场航段水平投影曲线L(s)在点L(s3)处的法平面nplane(s3);Db是外后侧竖直面Bsf(s1)与后侧面竖直面bsf(s1)之间的距离;
DL是外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离;
DR是外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离。
5.如权利要求4所述划设方法,其特征在于,所述S201还包括步骤如下:每条跑道离场程序管控区缓冲区为满足下述约束条件的最小区域:
约束条件1:外后侧竖直面Bsf(s1)与后侧竖直面bsf(s1)的距离,外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)的距离,外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)的距离均不小于ρ3,即min{Db,DL,DR}≥ρ3
约束条件2:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s2)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ3
约束条件3:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s4)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ2
其中:
ρ2是无人机被反制达最大高度时对应的水平距离;ρ3为意图控制无人机到无人机被反制落地期间飞行的最大水平距离;
s4是满足条件
Figure FDA0002322166480000041
的s的最大值;
h2是无人机被反制落地期间飞行达到的最大高度;
Figure FDA0002322166480000042
为h2的等高面。
6.如权利要求3所述划设方法,其特征在于,所述S201还包括步骤如下:
每条跑道进场程序管控区缓冲区界限包括:进场航段曲线
Figure FDA0002322166480000043
s∈[0,s1]对应的
内侧竖直面:每条跑道进场程序管控区内核区围成的三维空间对应的四个竖直面;
外后侧竖直面Bsf(s1):与后侧面竖直面bsf(s1)平行的竖直平面
Figure FDA0002322166480000044
外左侧竖直面Lsf(s1):位于左侧竖直面lsf(s1)左侧,且距lsf(s1)上每一点距离均为DL的竖直面,即以曲线
Figure FDA0002322166480000045
为底边的竖直面
Lsf(s1)={(x,y,z)|(x,y,0)∈LL(s),z≥0}
外右侧竖直面Rsf(s1):位于右侧竖直面rsf(s1)右侧,且距rsf(s1)上每一点距离均为DR的竖直面,即以曲线
Figure FDA0002322166480000051
为底边的竖直面
Rsf(s1)={(x,y,z)|(x,y,0)∈LR(s),z≥0}
外前侧竖直面Fsf(s1):进场航段水平投影曲线L(s)在点s3处的法平面nplane(s3);
Db是外后侧竖直面Bsf(s1)与后侧面竖直面bsf(s1)之间的距离;
DL是外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)之间的距离;
DR是外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)之间的距离。
7.如权利要求6所述划设方法,其特征在于,所述S201还包括步骤如下:
每条跑道进场程序管控区缓冲区为满足下述约束条件的最小区域:
约束条件1:外后侧竖直面Bsf(s1)与后侧竖直面bsf(s1)的距离,外左侧竖直面Lsf(s1)与左侧竖直面lsf(s1)的距离,外右侧竖直面Rsf(s1)与右侧竖直面rsf(s1)的距离均不小于ρ3,即min{Db,DL,DR}≥ρ3
约束条件2:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s2)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ3
约束条件3:法平面nplane(s3)位于外左侧竖直面Lsf(s1),外右侧竖直面Rsf(s1)的部分子集与法平面nplane(s4)位于左侧竖直面lsf(s1),右侧竖直面rsf(s1)的部分子集的距离不小于ρ2
其中:
ρ2是无人机被反制达最大高度时对应的水平距离;ρ3为意图控制无人机到无人机被反制落地期间飞行的最大水平距离;
参数值s4是满足条件
Figure FDA0002322166480000052
的s的最大值;
h2是无人机被反制落地期间飞行的最大高度;
Figure FDA0002322166480000053
为h2的等高面。
CN201911306980.7A 2019-12-17 2019-12-17 一种基于碰撞风险的机场无人机管控区域划设方法 Active CN111142555B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911306980.7A CN111142555B (zh) 2019-12-17 2019-12-17 一种基于碰撞风险的机场无人机管控区域划设方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911306980.7A CN111142555B (zh) 2019-12-17 2019-12-17 一种基于碰撞风险的机场无人机管控区域划设方法

Publications (2)

Publication Number Publication Date
CN111142555A true CN111142555A (zh) 2020-05-12
CN111142555B CN111142555B (zh) 2020-12-22

Family

ID=70518716

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911306980.7A Active CN111142555B (zh) 2019-12-17 2019-12-17 一种基于碰撞风险的机场无人机管控区域划设方法

Country Status (1)

Country Link
CN (1) CN111142555B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112911225A (zh) * 2021-01-19 2021-06-04 深圳科盾量子信息科技有限公司 基于量子加密的视频监控方法
CN113359834A (zh) * 2021-06-22 2021-09-07 中国民用航空总局第二研究所 一种无人机运行监控方法、系统及监控平台
CN115630531A (zh) * 2022-12-19 2023-01-20 中南大学 无人机控制系统自动化安全评估方法
CN118506618A (zh) * 2024-07-17 2024-08-16 北京航空航天大学杭州创新研究院 一种低空适飞空域边界划设方法、装置、存储介质及设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105247593A (zh) * 2014-04-17 2016-01-13 深圳市大疆创新科技有限公司 飞行禁区的飞行控制
CN107689166A (zh) * 2017-10-16 2018-02-13 哈工大机器人集团重庆普耀信息产业发展有限公司 一种机场飞行区监控管理系统及管理方法
CN108733073A (zh) * 2018-05-21 2018-11-02 厦门安胜网络科技有限公司 一种区域内无人机管控系统、方法及可读介质
CN110942672A (zh) * 2018-09-21 2020-03-31 郑州信大捷安信息技术股份有限公司 用于禁飞区的无人机管控系统及方法
CN111065983A (zh) * 2018-10-18 2020-04-24 深圳市大疆创新科技有限公司 限飞区生成方法和设备、控制无人飞行器飞行的方法和设备
CN111145597A (zh) * 2019-12-17 2020-05-12 中国民用航空总局第二研究所 一种基于碰撞风险的飞行航段的无人机管控区缓冲区划设方法
CN111192481A (zh) * 2019-12-17 2020-05-22 中国民用航空总局第二研究所 一种基于碰撞风险的进离场程序无人机管控区边界确定方法
CN111191350A (zh) * 2019-12-17 2020-05-22 中国民用航空总局第二研究所 一种飞行航段碰撞风险区的划设方法及应用

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105247593A (zh) * 2014-04-17 2016-01-13 深圳市大疆创新科技有限公司 飞行禁区的飞行控制
CN107689166A (zh) * 2017-10-16 2018-02-13 哈工大机器人集团重庆普耀信息产业发展有限公司 一种机场飞行区监控管理系统及管理方法
CN108733073A (zh) * 2018-05-21 2018-11-02 厦门安胜网络科技有限公司 一种区域内无人机管控系统、方法及可读介质
CN110942672A (zh) * 2018-09-21 2020-03-31 郑州信大捷安信息技术股份有限公司 用于禁飞区的无人机管控系统及方法
CN111065983A (zh) * 2018-10-18 2020-04-24 深圳市大疆创新科技有限公司 限飞区生成方法和设备、控制无人飞行器飞行的方法和设备
CN111145597A (zh) * 2019-12-17 2020-05-12 中国民用航空总局第二研究所 一种基于碰撞风险的飞行航段的无人机管控区缓冲区划设方法
CN111192481A (zh) * 2019-12-17 2020-05-22 中国民用航空总局第二研究所 一种基于碰撞风险的进离场程序无人机管控区边界确定方法
CN111191350A (zh) * 2019-12-17 2020-05-22 中国民用航空总局第二研究所 一种飞行航段碰撞风险区的划设方法及应用

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JIANPING ZHANG 等: "A Sector Capacity Assessment Method Based on Airspace Utilization Efficiency", 《JOURNAL OF PHYSICS》 *
任家龙: "终端区无人机飞行限制区划设研究", 《CNKI》 *
闫少琨: "无人机运行安全风险评价", 《CNKI》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112911225A (zh) * 2021-01-19 2021-06-04 深圳科盾量子信息科技有限公司 基于量子加密的视频监控方法
CN113359834A (zh) * 2021-06-22 2021-09-07 中国民用航空总局第二研究所 一种无人机运行监控方法、系统及监控平台
CN113359834B (zh) * 2021-06-22 2022-06-10 中国民用航空总局第二研究所 一种无人机运行监控方法、系统及监控平台
CN115630531A (zh) * 2022-12-19 2023-01-20 中南大学 无人机控制系统自动化安全评估方法
CN118506618A (zh) * 2024-07-17 2024-08-16 北京航空航天大学杭州创新研究院 一种低空适飞空域边界划设方法、装置、存储介质及设备

Also Published As

Publication number Publication date
CN111142555B (zh) 2020-12-22

Similar Documents

Publication Publication Date Title
Watkins et al. Ten questions concerning the use of drones in urban environments
CN110243359B (zh) 基于低空风预测模型的安全航迹规划方法
CN111142555B (zh) 一种基于碰撞风险的机场无人机管控区域划设方法
US20200393852A1 (en) Three dimensional aircraft autonomous navigation under constraints
CN111192481B (zh) 一种基于碰撞风险的进离场程序无人机管控区边界确定方法
CN103116360B (zh) 一种无人机避障控制方法
US8275499B2 (en) Device for calculating a flight plan of an aircraft
CN202221566U (zh) 一种基于性能导航的飞行程序设计系统及验证平台
CN102867073A (zh) 一种基于性能导航的飞行程序设计系统及验证平台和验证方法
KR20170080354A (ko) 무인기 안전 비행을 위한 가상적 스카이웨이와 이를 적용한 관제시스템 및 무인기 항법장치와 서비스
Bulusu et al. Capacity estimation for low altitude airspace
CN104615849A (zh) 一种适用于通用航空的飞行计划评估系统及实现方法
CN111145597B (zh) 一种基于碰撞风险的飞行航段的无人机管控区缓冲区划设方法
CN112053593B (zh) 一种基于风险评估的机场无人机管控区域划设方法
CN110322733A (zh) 一种侧向跑道保护区到达起飞窗建立方法
CN113655812A (zh) 一种大型货运无人机空中发动机熄火的着陆方法
Bulusu Urban air mobility: Deconstructing the next revolution in urban transportation-feasibility, capacity and productivity
Zhang et al. Empirical study of airport geofencing for unmanned aircraft operation based on flight track distribution
Kulik et al. Intelligent transport systems in aerospace engineering
Lou et al. An RRT* based method for dynamic mission balancing for urban air mobility under uncertain operational conditions
KR20230078097A (ko) 도심항공교통 실증 관리를 위한 디지털 트윈 기술에 기반한 3d 가시화 방법
CN111191350B (zh) 一种飞行航段碰撞风险区的划设方法及应用
US20230312116A1 (en) Aerial vehicle and control method thereof, using hybrid distributed propulsion system
Shen et al. A dynamic airspace planning framework with ads-b tracks for manned and unmanned aircraft at low-altitude sharing airspace
US20230410666A1 (en) 3d space data generation method, device and computer program for flight guidance of aircraft

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