0 引言
1 低空空域优化模型
1.1 问题描述
1.2 基于流量驱动的空域分层优化模型
1.2.1 目标函数
1.2.2 决策变量
1.2.3 约束条件
2 FGA-ASO算法设计
2.1 求解方法选择
2.2 设计思路
2.2.1 编码与初始化策略
2.2.2 约束处理与变量修复机制
2.2.3 算子设计
2.3 FGA-ASO算法流程
3 实验设计
3.1 仿真环境建模与关键要素提取
3.1.1 高层建筑物约束与地理围栏生成
3.1.2 出行需求建模与起降场网络生成
3.2 仿真环境及参数配置
表1 仿真参数配置表 |
| 参数 类别 | 参数名称 | 参数值 | 设置依据与说明 |
|---|---|---|---|
| 空域 结构 参数 | 空域高度 范围 | 70~870 m | 覆盖我国G类空域 (真高300 m以下)并向上延伸,留出缓冲空间 |
| 分层 数量 | 8层 | 在管理复杂性和空域 利用率之间取得平衡 | |
| 层高 | 100 m | 参考典型无人机垂直性能,提供安全的垂直间隔 | |
| 飞行 器性 能参 数 | 巡航 速度 | 100 km/h (≈27.8 m/s) | 参考主流载物/载人 eVTOL机型(如EH216-S)的典型巡航速度 |
| 初始飞行 策略 | 直线飞行 (大圆航线) | 简化模型,以最直接路径连接起降点,用于基础 冲突评估 | |
| 冲突 安全 参数 | 水平安全 间隔 | 50 m | 综合考虑无人机导航精度、通信延迟与安全余量 |
| 垂直安全 间隔 | 25 m | 约为1/2水平间隔, 符合航空器间隔标准惯例 | |
| 建筑物安全容差 | 30 m | 为飞行器与建筑物之间 预留的额外安全缓冲区 | |
| 冲突探测前瞻时间 | 300 s | 提供足够的时间用于 预测冲突 | |
| 优化 算法 参数 | 种群规模 | 100 | 在计算效率与解空间探索能力之间取得平衡 |
| 迭代代数 | 120 | 通过前期实验验证, 足以使算法收敛至满意解 | |
| 基准交叉概率 | 0.9 | 采用较高的概率以促进 优良基因的交换 | |
| 变异概率 | 0.2 | 设置适当的扰动概率以 维持种群多样性, 避免早熟收敛 | |
| 交叉分布 指数 | 20 | 控制SBX交叉操作,数值越大,子代越接近父代 | |
| 变异分布 指数 | 20 | 控制多项式变异操作,数值越大,变异越局部化 |
3.3 评估指标
3.3.1 冲突数(Conflict Count, CC)
3.3.2 侵入数(Intrusion Count, IC)
3.3.3 侵入严重性(Intrusion Severity, IS)
3.4 空域结构优化算法对比实验
3.4.1 对比基线
3.4.2 实验流程
4 结果分析
4.1 算法收敛特性分析
4.2 解质量综合评价
表2 空域结构优化算法对比实验结果 |
| 优化算法 | 目标函数值 f | 总冲突数 | 飞行器间冲突 | 飞行器-建筑冲突数 | 平均侵入严重性 | 平均计算时间 /s |
|---|---|---|---|---|---|---|
| FGA-ASO(本文算法) | 435.31 | 386 | 385 | 1 | 0.493 1 | 88 605.83 |
| 标准遗传算法 | 543.53 | 493 | 490 | 3 | 0.505 3 | 96 922.77 |
| 精英保留的遗传算法(EGA) | 473.62 | 423 | 421 | 2 | 0.506 2 | 96 046.98 |
| 差分进化算法(DE) | 485.98 | 433 | 431 | 2 | 0.529 8 | 97 818.96 |
| 稳态遗传算法(steadyGA) | 525.74 | 475 | 473 | 2 | 0.507 4 | 95 823.53 |
| 种马遗传算法(studGA) | 488.64 | 437 | 436 | 1 | 0.516 4 | 93 715.30 |
| 代际差分遗传算法(GGAP_SGA) | 460.25 | 410 | 408 | 2 | 0.502 5 | 92 103.42 |
| 均匀分层基准 | 581.02 | 528 | 525 | 3 | 0.530 2 | — |
4.3 FGA-ASO改进机制效果分析
4.3.1 消融实验与对比设置
表3 各算法变体性能对比 |
| 算法变体 | 总冲突数 | 平均侵入 严重性 | 收敛 速度 |
|---|---|---|---|
| FGA-ASO-Full | 386.0±8.2 | 0.493±0.012 | 40±5 |
| FGA-ASO-noHIS | 398.5±10.1 | 0.502±0.015 | 55±7 |
| FGA-ASO-noFWC | 415.2±12.5 | 0.510±0.018 | 45±6 |
| FGA-ASO-noAMS | 394.3±9.8 | 0.513±0.014 | 42±5 |
4.3.2 流量加权交叉(FWC)的作用机制分析
表4 空域结构分层对比分析 |
| 方案 | 层 号 | 层边界 | 层宽 /° | 层内 航班 数 | 占总 流量比 例(%) | 层内平均 密度 (架次/°) |
|---|---|---|---|---|---|---|
| FGA-ASO (带FWC) | 1 | 0°~42° | 42 | 30 | 11.7 | 0.714 |
| 2 | 42°~116° | 74 | 57 | 22.3 | 0.770 | |
| 3 | 116°~185° | 69 | 33 | 12.9 | 0.478 | |
| 4 | 185°~223° | 38 | 35 | 13.7 | 0.921 | |
| 5 | 223°~242° | 19 | 26 | 10.2 | 1.368 | |
| 6 | 242°~295° | 53 | 44 | 17.2 | 0.830 | |
| 7 | 295°~320° | 25 | 11 | 4.3 | 0.440 | |
| 8 | 320°~360° | 40 | 20 | 7.8 | 0.500 | |
| 无FWC 变体 | 1 | 0°~55° | 55 | 39 | 15.2 | 0.709 |
| 2 | 55°~110° | 55 | 38 | 14.8 | 0.691 | |
| 3 | 110°~165° | 55 | 41 | 16.0 | 0.745 | |
| 4 | 165°~200° | 35 | 17 | 6.6 | 0.486 | |
| 5 | 200°~245° | 45 | 46 | 18.0 | 1.022 | |
| 6 | 245°~290° | 45 | 33 | 12.9 | 0.733 | |
| 7 | 290°~335° | 45 | 29 | 11.3 | 0.644 | |
| 8 | 335°~360° | 25 | 13 | 5.1 | 0.520 |
注:表中数据经四舍五入处理,分项之和存在微小误差。 |