考虑多因素影响的无人机碰撞风险评估_童亮.pdf

第 44 卷第 4 期兵 器 装 备 工 程 学 报2023 年 4月收稿日期:2022 06 08;修回日期:2022 07 24基金项目:陕西省自然科学基础研究计划(2022JM-412)作者简介:童亮(1995)，男，硕士研究生，助理工程师，E-mail:1044149827 qq com。通信作者:甘旭升(1971)，男，博士，副教授，E-mail:ganxusheng123163 com。doi:1011809/bqzbgcxb202304040考虑多因素影响的无人机碰撞风险评估童亮1，甘旭升1，张宏宏1，金剑勇2(1 空军工程大学 空管领航学院，西安710051;2 中国人民解放军 95140 部队，广东惠州516000)摘要:采用蒙特卡洛法评估多因素影响下的无人机与民航客机的碰撞风险，从而对进入融合空域的无人机冲突解脱概率进行计算。首先建立了基于速度或者航向调整方法的无人机避撞数学模型，随后对融合空域内影响无人机与民航客机定位误差的关键因素构成与分布进行了分析，并利用蒙特卡洛法对典型冲突场景下的碰撞风险进行仿真。结果表明:航向解脱与速度解脱 2 种策略在一定程度上均能避免无人机与民航客机之间危险接近，在同等误差条件下，速度解脱策略造成的严重事故往往更低;民航客机定位误差对避障效果影响较大;得到了给定误差条件下，不同冲突场景允许的最小安全间隔标准，能够为未来融合空域内无人机与民航客机、电子围栏安全间隔的设置提供理论依据。关键词:安全评估;速度障碍法;冲突解脱;蒙特卡洛本文引用格式:童亮，甘旭升，张宏宏，等 考虑多因素影响的无人机碰撞风险评估 J 兵器装备工程学报，2023，44(4):282 289Citation format:TONG Liang，GAN Xusheng，ZHANG Honghong，et al isk assessment of UAV collision consideringmultiple factors J Journal of Ordnance Equipment Engineering，2023，44(4):282 289中图分类号:V271 4;TB24;TP391 9文献标识码:A文章编号:2096 2304(2023)04 0282 08isk assessment of UAV collision considering multiple factorsTONG Liang1，GAN Xusheng1，ZHANG Honghong1，JIN Jianyong2(1 Air Traffic Control and Navigation College，Air Force Engineering University，Xi an 710051，China;2 Unit 95140 of the Chinese People s Liberation Army，Huizhou 516000，China)Abstract:This paper uses the Monte Carlo method to evaluate collision risks between a UAV and a civilairliner under the influence of multiple factors so as to evaluate and calculate the conflict resolutionprobability of UAV entering the fused airspace Firstly，a mathematical model of UAV collision avoidancebased on speed or course adjustment methods is established Then，the composition and distribution of keyfactors affecting the positioning errors of the UAV and the civil airliner in the fused airspace are analyzed，and the simulation of the collision risks in typical conflict scenarios are conducted using the Monte Carlomethod The results show that:Both course release and speed release strategies can avoid dangerousapproaches between a UAV and a civil airliner to a certain extent，and the serious accidents caused by thelatter strategy are less likely to occur under the same error conditions The positioning error of the civilairliner has a greater impact on obstacle avoidance The minimum safety interval criteria allowed fordifferent conflict scenarios under a given error condition are obtained，which can be used as the theoreticalbasis for the future setting of safety intervals in fused airspace between a UAV and a civil airliner as well asbetween electronic fencesKey words:safety assessment;velocity obstacle method;conflict resolution;Monte Carlo0引言随着航空工业技术的快速发展，无人机凭借其机动能力强、成本低、风险小等特点，在军事、农业、交通等领域发挥越来越重要的作用。由于缺乏可靠的冲突解脱能力，无人机目前通常在隔离空域内运行1。融合空域是指无人机与有其他有人机同时运行的空域。为了保证无人机在融合空域执行任务的可靠性，降低对空域带来的不安全影响，需要对无人机碰撞风险概率进行评估2。国内外对航空器飞行碰撞风险评估进行了一系列研究，最早用于碰撞风险评估的模型是 eich 模型3 4，分析了由于航空器位置误差造成的碰撞风险。EVENT 模型用来评估纵向间隔与侧向间隔的碰撞风险5。石磊等基于三维布朗运动模型，计算出航空器碰撞概率6。王松涛等人综合考虑高空风等环境因素对碰撞概率的影响7。杨新湦以碰撞风险为量化指标，基于尾涡强度对无人机与民机的纵向安全间隔进行评估8。张宏宏为解决无人机在复杂低空环境内运行安全性较低的现状，提出了基于飞行环境的风险评估方法，为低空无人机航路规划方法提供参考9。谷志鸣利用蒙特卡洛对冲突解脱算法可靠性进行评估10。张兆宁提出基于考虑所需 CNS(communication，navigation and surveillance)性能、防撞系统和高空风等影响因素作用下的碰撞风险评估11。目前，国内外对民航客机之间的碰撞风险已有一定的研究，但是对融合空域内无人机避撞风险的研究相对较少。传统研究中，有人航空器之间的碰撞风险集中于考虑人为因素的影响，学者们从不同的角度、维度对人为因素的进行了定性、定量的阐述。无人机没有机载的飞行驾驶员对飞机进行直接操纵，缺乏有人机与有人机之间“看见即躲避”的避撞方式，使得融合空域内无人机与有人机在避撞过程中主要依赖于无线信号/指令的传输，无线信号/指令的传输过程中易受到多种环境因素的干扰，导致信号/指令出现中断、丢包等现象，进而影响无人机避撞行为的执行，从而引发相撞事故。无人机进入融合空域与有人机共同运行是一种趋势，尤其是无人机与民航班机在融合空域中的共同运行。因此，需要对无人机在融合空域内的碰撞风险评估进行深入研究。本文针对无人机进入融合空域的现实需求，引入适用于无人机的速度障碍法模型。对融合空域内影响无人机与民航客机定位误差的关键因素构成与分布进行了分析。利用蒙特卡洛方法对经典冲突场景下的碰撞风险进行仿真，利用具体的理论数据得出无人机与有人机之间碰撞风险评估结果，为无人机进入融合空域提供一种理论参考。1无人机避撞模型针对无人机在融合空域内避障问题，国内外学者将智能优化算法12 13、势场法14 15、几何法16 17 等方法应用在冲突解脱领，均取得不错效果。Jrmy Omer18 将基于时间离散的航迹优化算法和基于空间离散化的航迹优化算法以及单次调整的航迹优化算法进行了比较。从综合表现来说单次调整优化算法能够解决更为复杂的情况，并且算法得到的冲突消解机动产生的消耗较少，计算时间适中。因此，本文选取基于速度障碍法的集中式单次调整优化方法作为评估对象。速度障碍法定义了一个障碍区域，当无人机与入侵障碍物的相对速度落在障碍区域之外，便可完成冲突解脱19。如图 1 所示，A 与 O 点分别为无人机与动态障碍物的位置，速度分别为 v1、v2，相对速度为 v=v1 v2，无人机与障碍物之间的安全距离设为 d1，无人机与障碍物之间的连线 AO 与v1、v2的夹角分别为 1、2。无人机与障碍物之间的动态关系可转化成以相对速度 v运行的无人机与静止障碍物之间的关系。图 1速度障碍模型示意图Fig 1 Diagram of speed barrier model无人机可以通过改变航向与速度，使得相对速度 v落在速度障碍区之外，实现避障。1 1航向调整如图 2 所示，无人机仅通过调整航向进行避障，无人机航向调整前后与坐标轴的夹角分别为 1与 1，相对速度调整前后分别为 v与 v，无人机速度调整前后分别为 v1与 v1，则有，根据文献 20，冲突解脱过程中，无人机航向调整量1为:1=1 1=arcsin(v2sin(2)/v1)1(1)图 2无人机航向调整前后相对速度矢量示意图Fig 2 Schematic diagram of relative velocity vector ofUAV before and after course adjustment1 2速度调整如图 3 所示，无人机仅通过调整速度大小进行避障，根382童亮，等:考虑多因素影响的无人机碰撞风险评估据文献 20，冲突解脱过程中，无人机速度调整量 v1为:v1=v1 v1=v2sin(2+)/sin(1)v1(2)图 3无人机速度调整前后相对速度矢量示意图Fig 3 Schematic diagram of relative velocity vector ofUAV before and after speed adjustment基于速度障碍法的集中式单次调整优化方法，无人机实现对障碍物规避并及时恢复到原航迹20，效果较好，但假设过于理想，忽略了无人机运行不确定性因素对避障的影响，因此需要综合考虑影响因素，对碰撞风险进行评估。2无人机运行影响因素分析无人机在融合空域内运行时，受到多随机因素影响，定位误差随着影响因素变化，造成一定程度的偏航，甚至导致无人机与空域内动态障碍物的危险接近、空中相撞，威胁无人机安全运行。本文选取影响无人机偏航的关键因素进行研究，假设造成无人机误差的各个环境要素之间相互独立，对多因素影响下的无人机碰撞风险进行评估。2 1导航性能对无人机定位影响目前，无人机定位方式主要是采用 GPS 卫星导航，其误差来源主要有卫星误差、传播误差与接收误差。卫星误差主要包括星历误差、星钟误差以及相对论效应影响等，传播误差主要包括电离层延迟误差、对流层折射误差与多路径误差，接收误差主要包括接收机时钟误差、噪声误差与接收机天线相位中心误差21。GPS 一般采用三星定位原理，除 3个用于定位的卫星外，还需要一个原子钟，由第四个卫星代替，假设各个卫星对无人机定位影响相互独立，互不干涉。本文以 P 码编码方式为例，图 4 给出了各种误差来源造成的定位误差以及系统总误差21。根据定位原理，无人机 GPS 定位误差 x在运行速度方向服从均值 x为 0，方差为 x的高斯正态分布，标准差 x=4 16 12=32 2 m=0 032 2 km，定位误差 x密度函数表示为:f(x)=12xexp2x22()x(3)图 4各类误差对定位误差的影响Fig 4 The effect of various errors on positioning errors2 2气象因素对无人机定位影响空域内运行中的无人机极易受到不确定气象因素的影响，造成定位误差，假设无人机飞行所受顺风及逆风频率相同，无人机在高空风影响下的定位误差 w符合正态分布22，密度函数表示为:f(w)=12wexp(x w)222()w(4)式中，w与 w分别表示单独考虑气象因素(高空风)时，无人机定位误差的均值与方差。2 3无人机状态改变造成的位置误差融合空域内的无人机通过改变航向或速度，实现冲突解脱，而在实际运行时，无人机航向与速度的改变并非瞬间完成，存在一定的延迟与位置误差。无人机系统的时间延迟不仅与机身操纵系统到响应时间、数据链信息传输延迟等因素有关，还与无人机地面站信息处理时间、链路上下传输延迟等因素有关23。2 3 1改变航向时的定位误差根据速度障碍法避障原理，无人机在 AB 航线上运行，到B 点进行瞬间转向，再向 BC 方向继续运行。而实际运行的无人机受性能约束，为实现等效转向效果，需要在 A 点以 P为圆心，为半径进行转向，运行到 C 点之后恢复直线航行，为理论转向角度，如图 5 所示。图 5无人机改变航向示意图Fig 5 Diagram of UAV changing course482兵 器 装 备 工 程 学 报http:/bzxb cqut edu cn/由于无人机转弯时受到指令延迟、转弯速度等因素影响，无人机期望转弯角度与实际转角之间产生偏差，其实际的转弯角度 服从正态分布，概率密度表示为:f()=12exp()222()(5)2 3 2改变速度时的定位误差根据速度障碍法避障原理，无人机可改变自身速度大小实现避撞。实际无人机速度变化受自身性能约束，不可能瞬间完成变速，如图6 所示，无人机运行速度由 v1变为 v1时，飞行航程 S 为图中阴影面积，可表示为:S=12(v1+v1)tf(6)式中，tf为无人机变速结束时间，受自身性能约束，可表示为:tf=v1 v1amax=v1amax(7)式中，amax为无人机最大加速度。图 6无人机改变速度示意图Fig 6 Diagram of UAV changing speed无人机按照最大加速度 amax指令改变速度时，受到不确定性因素干扰影响，实际最大加速度 amax服从正态分布，表示为:f(a)=12aexp(a amax)222()a(8)式中，a是加速度的方差，取 0 01。2 4民航客机定位误差无人机在融合空域内运行时动态障碍物类型众多，这里选取民航客机作为动态障碍物，无人机通过执行避撞策略实现冲突解脱。根据文献 11，运行中的航空器由于定位误差，会偏离理论航行路线，造成碰撞风险。民航飞机导航误差通常由导航精度引起的，密度函数服从期望值为 0 的双指数分布，可表示为:y=12bexpy a()b(9)式中:a 为位置参数，b 为尺度参数且 b 0，二者为偏航误差概率密度中的参数。根据所需导航性能(required navigational performance，NP)定义，NPn 指航空器在 n 海里之内，以 95%的概率偏离航线的距离11。则有:nn12aexpy()ady=0 95(10)3仿真分析为验证速度障碍法避障的有效性，本节运用环境进行仿真，其中无人机与民航飞机的初始位置、航向与运行速度给定，应用蒙特卡洛算法，分析运行中的不确定性因素对避障算法的影响。根据无人机运行特点，无人机安全圆半径设为15 km，机载传感器探测半径设为 100 km。相对于巡航阶段，民机进离场与起降阶段时飞行高度较低，也是与无人机发生冲突，产生危险接近最多的阶段。本文构建的冲突场景为无人机与进离场与起降阶段时民航客机航线交叉相遇下的冲突解脱场景，结合我国航空器飞行间隔标准，根据航空器运行最小间隔距离，对碰撞风险进行分类评估24，如表 1 所示，当两机最小间隔距离小于 1 76 km时，定为严重事故，当两机最小间隔距离处于 1 76 3 km时，定为危险接近，当两机最小间隔距离处于 3 10 km 时，定为飞行冲突，当两机最小间隔距离处于 10 15 km 时，定为存在风险，当两机最小间隔距离大于 15 km 时，定为安全运行。表 1事故风险等级划分Table 1 Classification of accident risk levels运行最小间隔/km碰撞风险等级 0，1 76严重事故(1 76，3危险接近(3，10飞行冲突(10，15存在风险(15，)安全运行3 1冲突场景仿真选取大型无人机与民航航机作为碰撞对象，融合空域内不确定因素参数具体数据如表 2 所示。无人机与民航客机初始信息如表 3 所示。表 2所需参数值Table 2 equired parameter value物理量变量分布参数取值无人机 GPS定位误差x正态分布x=0，x=0 032 2高空风影响下的定位误差w正态分布w=0，w=0 01实际转弯角度正态分布=，=0 01加速度a正态分布a=amax，a=0 01偏航误差概率密度y双指数分布a=0 2582童亮，等:考虑多因素影响的无人机碰撞风险评估场景 1、场景 2 分别为无人机任务航线与民航客机航线夹角分别为锐角、钝角的冲突场景，根据速度障碍法避障策略，无人机通过单独调整航向，改变无人机与民航客机之间的相对速度，实现融合空域内的避撞。场景 1 中无人机与民航客机之间的理论航迹图、双方理论间隔与实际间隔分别如图 7 所示。场景 3、场景 4 分别为无人机任务航线与民航客机航线夹角分别为锐角、钝角的冲突场景，根据速度障碍法避障策略，解脱策略采用速度解脱，无人机通过单独调整速度大小，改变无人机与民航客机之间的相对速度，实现融合空域内的避撞。场景 3 中无人机与民航客机之间的理论航迹图、双方理论间隔与实际间隔分别如图 8 所示。场景 2、场景 4 避撞场景、间隔对比仿真图与场景 1、场景 3 类似，限于篇幅不再罗列。表 3无人机与民航客机初始信息Table 3 UAV and civil aircraft initial information场景解脱策略冲突类型(航线夹角)冲突方起点/km航向/()速度/(kmh1)仿真步长/h12航向解脱锐角(45)钝角(45)无人机(0，200)0300民航客机(60，0)75600无人机(50，100)15300民航客机(350，20)1506000 0010 00134速度解脱锐角(45)钝角(45)无人机(0，200)0300民航客机(60，0)75600无人机(50，100)15300民航客机(350，20)1506000 0010 001图 7场景 1 仿真示意图Fig 7 Scenario 1 simulation diagram图 8场景 3 仿真示意图Fig 8 Scenario 3 simulation diagram682兵 器 装 备 工 程 学 报http:/bzxb cqut edu cn/采用蒙特卡洛仿真模拟，对场景 1 场景 4 分别进行 10万次重复实验，统计无人机与民航客机之间的风险等级分别处于严重事故、危险接近、飞行冲突、存在风险、安全运行等 5个等级下的发生次数，具体数据如表 4 所示，分布如图 9所示。表 4不同场景下事故率统计Table 4 Accident rate statistics under different scenarios风险等级航向解脱场景 1场景2速度解脱场景3场景4严重事故0002 480006 040000 240000 38危险接近0007 120004 240000 540000 84飞行冲突0144 140227 230006 450004 81存在风险0351 050346 370339 710455 75安全运行0495 210416 120653 240538 22图 9不同场景下事故率分布Fig 9 Accident rate distribution under different scenarios分析场景 1 场景 4 中各风险等级发生概率，可以发现航向解脱与速度解脱 2 种策略在一定程度上均能避免无人机与民航客机之间的严重事故，但离国际民航组织提出的有人机总体安全目标水平(target level of safety，TLS)次事故/飞行小时还有一定差距;同时在同等误差条件下，速度解脱策略造成的严重事故往往更低，因此在同等冲突场景下，应优先选择速度解脱策略进行避撞。3 2民航客机定位误差对安全目标水平的影响为更精确描述民航客机定位误差对安全目标水平的影响程度，针对场景 1，保持其他误差参数不变，对不同 NP 性能环境下的避撞效果进行仿真，如图 10 所示，无人机与民航客机的碰撞严重事故率随着避撞间隔的增大而降低，同时，民航客机导航性能越精确，达到同等安全目标时所需的避撞间隔越小。3 3无人机转向误差对安全目标水平的影响将无人机转向误差方差值 分别设为 0 05、0 01 与0.15，对在 3 种转向误差下的无人机安全目标水平进行实验，得到如图 11 所示结果，仿真结果表明:在合理误差下，无人机转向误差对安全目标水平影响程度较小。图 10民航客机定位精度对安全目标水平的影响Fig 10 Influence of positioning accuracy on safetytarget level of civil aircraft图 11无人机转向误差对安全目标水平的影响Fig 11 Influence of UAV steering error on safety target level3 4高空风影响下的定位误差对安全目标水平的影响选取高空风影响下的定位误差为研究对象，分别研究物理量在、2、3 下的误差，具体分布如图 12 所示。可知，同一避撞间隔下，定位误差取值增大，安全目标水平的取值相应增大，安全性越低;随着避撞间隔的增大，安全目标水平取值减小，安全性提高。图 12不同高空风误差下的安全目标水平示意图Fig 12 Safety target level diagram under differentupper air wind errors782童亮，等:考虑多因素影响的无人机碰撞风险评估3 5无人机碰撞间隔设置目前，无人机与民航客机之间的安全间隔距离暂无明确标准，本文将安全间隔距离分别设为 10 20 km，以场景 1 为研究对象，在同等误差条件下进行1 000 万次蒙特卡洛仿真，各安全间隔下的严重事故发生概率如图 13 所示。可以看出，无人机与民航客机发生严重事故概率随安全间隔的增大而减小，当安全间隔变为 25 km 左右时，严重事故发生概率大约与有人机总体安全目标水平(次事故/飞行小时)相等，为无人机融入有人机空域的安全间隔设置提供了评估方式。图 13严重事故概率变化趋势Fig 13 Trend of serious accident probability4结论当前国内外对无人机在融合空域内冲突解脱技术的研究较多，实验案例往往仅验证解脱算法的有效性，很少有从安全性的角度对解脱算法进行考虑，并且目前对融合空域内无人机运行安全评估体系还很缺乏。因此，利用蒙特卡洛技术对无人机与民航客机的冲突解脱算法进行大规模仿真，得到如下结论:1)航向解脱与速度解脱 2 种策略在一定程度上均能避免无人机与民航客机之间危险接近，避障间隔设置越小，离国际民航组织提出的有人机总体安全目标水平差距越大。在同等误差条件下，速度解脱策略造成的严重事故往往更低，因此在同等冲突场景下，应优先选择速度解脱策略进行避撞。2)无人机与民航客机的碰撞严重事故率随着避撞间隔的增大而降低，同时，民航客机导航性能越精确，达到同等安全目标时所需的避撞间隔越小。在合理误差下，无人机转向误差对安全目标水平影响程度较小。3)无人机与民航客机发生严重事故概率随安全间隔的增大而减小，当安全间隔设为 25 km 左右时，严重事故发生概率大约与有人机总体安全目标水平相等，为无人机与障碍物、电子围栏安全间隔的设置提供了一种理论计算方式。参考文献:1 张泽京，张曙光，柳旭，等 无人机系统安全目标水平预估方法 J 航空动力学报，2018，33(4):1017 1024ZHANG Zejing，ZHANG Shuguang，LIU Xu，et al Safetytarget level prediction method for unmanned aerial vehiclesystems J Journal of Aerospace Power，2018，33(4):10171024 2 邓力 无人机与民航客机碰撞概率研究 J 南京理工大学学报，2019，43(1):122 128DENG Li Study on collision probability between UAV andcivil aircraftJ Journal of Nanjing University of Scienceand Technology，2019，43(1):122 128 3 EICH P G Analysis of long-range air traffic systems:Sep-aration standards J Journal of the Institute of Navigation，1966，19(1):88 98 4 韩松臣，裴成功，隋东，等 平行区域导航航路安全性分析 J 航空学报，2006，27(6):1023 1027HAN Songchen，PEI Chenggong，SUI Dong，et al Safety A-nalysis of parallel regional navigation routeJ Journal ofAviation，2006，27(6):1023 1027 5PETE B Longitudinal collision risk for ATC track sys-tems:a hazardous event modelJ Journal of Navigation，2006，59(1):55 70 6 石磊，吴仁彪，黄晓晓 基于总体冲突概率和三维布朗运动的冲突探测算法 J 电子与信息学报，2015，37(2):360 366SHI Lei，WU enbiao，HUANG Xioaxiao Conflict detectionalgorithm based on total conflict probability and three-di-mensional Brownian motionJ Journal of ElectronicsInformation Technology，2015，37(2):360 366 7 王松涛，吕宗平，张兆宁 基于环境要素的自由飞行条件下航空器碰撞风险模型J 中国科技论文，2016，11(19):2183 2186WANG Songtao，LV Zongping，ZHANG Zhaoning Collisionrisk model of aircraft under free flight conditions based onenvironmental factors J Science and Technology Papers ofChina，2016，11(19):2183 2186 8 杨新湦，任治 基于尾涡强度的无人机与民机纵向安全间隔评估研究J 中国安全生产科学技术，2018，14(10):64 69YANG Xinnie，EN Zhi Longitudinalsafety interval evalua-tion of unmannedaerial vehicle and civil aircraft based onwake vortex intensity J China Safety Science and Tech-nology，2018，14(10):64 69 9 张宏宏，甘旭升，李双峰，等 复杂低空环境下考虑区域风险评估的无人机航路规划 J 仪器仪表学报，2021，42(1):257 266ZHANG Honghong，GAN Xusheng，LI Shuangfeng，et alouteplanning of UAV considering regional risk assessment882兵 器 装 备 工 程 学 报http:/bzxb cqut edu cn/in complex low-altitude environment J Chinese Journal ofScientific Instrument，2021，42(1):257 266 10 梁智韬，贾高伟，王建峰 无人机集群任务分配中群智能算 法 性 能 对 比J 火 力 与 指 挥 控 制，2022，47(6):28 37LIANG Zhitao，JIA Gaowei，WANG Jianfeng The compari-son of swarm intelligence algorithm aboutthe mission assign-ment of uav swarm J Fire Control Command Control，2022，47(6):28 37 11 张兆宁，余晟 自由飞行下基于航空器定位误差修正的碰撞风险研究 J 安全与环境学报，2017，17(5):16331636ZHANG Zhaoning，YU Sheng Collision risk study based onaircraft positioning error correction in free flight J Journalof Safety and Environment，2017，17(5):1633 1636 12 LI P Y Multi-objective path planning of unmanned aerialvehicle based on genetic algorithm J Agricultural Equip-ment Vehicle Engineering，2019，57(1):68 70，86 13 LI F L，CHEN S，FAN X J，et al Path planning based onfirefly algorithm in dynamic unknown environment J Au-tomation Instrumentation，2019，34(6):53 58 14 管祥民，吕人力 基于混合人工势场与蚁群算法的多飞行器冲突解脱方法 J 武汉理工大学学报(交通科学与工程版)，2020，44(1):28 33GUAN Xiangming，LV enli Multi-aircraft conflict resolu-tion method based on hybrid artificial potential field and antcolony algorithm J Journal of Wuhan University of Tech-nology(Transportation Science and Engineering)，2020，44(1):28 33 15 YANG X，YANG W，ZHANG H，et al A new method for ro-bot path planning based artificial potential field C/2016IEEE 11th Conference on Industrial Electronics and Appli-cations(ICIEA)IEEE，2016:1299 16 DUAND N Constant speed optimal reciprocal collision a-voidanceJ Transportation esearch Part C:EmergingTechnologies，2018，96(11):366 379 17 杨秀霞，张毅，周硙硙，等 基于空间障碍球冠的 UAV 保角映射避碰决策 J 华中科技大学学报(自然科学版)，2019，47(2):127 132YANG Xiuxia，ZHANG Yi，ZHOU Weiwei，et al Collisionavoidance decision of UAV conformal mapping based onspace obstacle ball crown J Journal of HuazhongUniversi-ty of Science and Technology(Natural Science Edition)，2019，47(2):127 132 18 JMYOME Comparison of mixed-integer linear modelsfor fuel-optimal air conflict resolution with recoveryJ IEEE Transactionson Intelligent Transportation Systems，2015(16):3126 3137 19 DUAND N，BANIE N Does ATM need centralized co-ordination?Autonomous conflict resolution analysis in aconstrained speed environmentJ Air Traffic ControlQuarterly，2015，23(4):1 9 20 张宏宏，甘旭升，李昂，等 基于速度障碍法的无人机避障与航迹恢复策略 J 系统工程与电子技术，2020，42(8):1759 1767ZHANG Honghong，GAN Xusheng，LI Ang，et al Obstacleavoidance and track recovery strategy of uav based on speedobstacle method J Systems Engineering and Electronics，2020，42(8):1759 1767 21 游广芝 GPS 导航定位中的误差分析与修正 D 哈尔滨:哈尔滨工业大学，2007YOU Guangzhi Error analysis and correction in GPS navi-gation D Harbin:Harbin Institute of Technology，2007 22 吕宗平，王松涛，张兆宁 基于 CEAM 和 IDAC 综合分析法的自由飞行下航空器碰撞风险评估 J 安全与环境工程，2017，24(1):111 114LV Zongping，WANG Songtao，ZHANG Zhaoning Collisionrisk assessment of aircraft under free flight based onCEAM and IDAC comprehens