６．継続上昇運航(CCO)に関する研究

干渉点

2016年 6月 9日 電子航法研究所研究発表会
横浜国立大学 上野誠也
航空交通管理領域 虎谷大地
1

本日の発表
１．研究の背景と課題
課題１ 他機との時間的・空間的間隔の確保
２-１ 時空間座標系 ２-２ 単機の軌道最適化

課題２ 干渉点での通過順序の最適化
３-１ 評価関数曲線

３-２ 軌道と順序の最適化

課題３ 実運用システムへの展開
４-１ 継続的最適化 ４-２ 計算量の低減
５ 研究の成果
2

1. 研究の背景と課題 １/５
航空交通管理(ATM)の変革• 現行のATM
Airspace B• 将来のATM
Airspace A

Airspace C
Integrated
ATM system

Airport A
Airport B

空域ベースのATM運用

軌道ベースのATM運用! 4次元（位置と時間）軌道運用が必要とされる
3

1. 研究の背景と課題 ２/５◆継続上昇運航（CCO）：離陸からの上昇区間において、水平飛行区間を用いずに、上昇する運航方式消費燃料の削減
運航時間の短縮

高度✓環境影響の軽減✓管制作業の低減
高度

消費燃料

飛行時間[ft]

[kg]

[s]

CCO 31,000

5878

STP

76102

CCO

STP

水平距離10,000

Roachら AIAA-2010-9379より引用

4

1. 研究の背景と課題 ３/５
関空のおける上昇機と降下機の干渉

降下機１
関空
降下機２

上昇機１

上昇機と降下機の干渉点がある。

5

1. 研究の背景と課題 ４/５
高度15000

上昇機１は高度８０００ｆｔで巡航し、降下機２を通した。
降下機１

Hp [ft]

10000
降下機２

5000

上昇機１
高度8000ftで巡航

0-40-20x [NM]0206

１. 研究の背景と課題 ５/５
課題１ 他機との時間的・空間的間隔の確保⇒単機の軌道最適化

課題２ 干渉点での通過順序の最適化⇒複数機の軌道と順序の最適化
課題３ 実運用システムへの展開⇒連続的な軌道生成⇒計算処理時間の短縮
7

課題１ 他機との時間的・空間的間隔の確保

時空間座標系を導入し、時間を空間座標と同じ状態変数として扱うことで、時間的間隔と空間的間隔のバランスを考慮した最適軌道が生成できる。単機の軌道最適化でその有効性を示す。
8

２-1 時空間座標系 １/３
航空機の垂直面内運動方程式
V V   (Thr  D  W sin  ) / m  d     ( L  W cos  ) /( mV ) 
V cos dt x Hp  V sin  9

２-1 時空間座標系 ２/３
 V   (Thr  D  W sin  ) / m  d     ( L  W cos  ) /( mV ) V cos dt x
 
Hp  V sin   V   (T  D  W sin  ) /(mV )  
2    ( L  W cos ) /(mV ) d 
x cos 
ds sin Hp 
 t  

1 /V
 

ds  Vdt10

２-1 時空間座標系 ３/３◆時空間座標系の利点
位置

一般の座標系

時空間座標系

A機

位置
A機

B機

時間
異なる時間の比較が困難である

時間

位置時間

B機

時間的間隔と空間的間隔を同等に扱うことができる

２-２ 単機の軌道最適化 １/4◆問題の設定降下機： CDOによる降下（優先度が高い）
上昇機： CCOによる上昇→最適回避

干渉点

12

２-２ 単機の軌道最適化 ２/４◆回避の条件
空間的間隔： RCR=5[nm] (約9260[m])時間的間隔： ΔtCR=44.0[s] (=9260[m]/210[m/s])

従来の座標系ではどちらかを選んで最適化計算を行った。
13

２-２ 単機の軌道最適化 ３/４◆時空間座標系における空間的間隔と時間的間隔楕円体で表現

x
空間的間隔による回避

時間的間隔による回避yt

14

２-２ 単機の軌道最適化 ４/４
・楕円体上の一点を通過する経路を求める
・出発時間も最適化する計算を行う

消費燃料消費燃料の比較（出発時間の調整）

出発遅延時間

出発時間を変更する時間的回避を行うと消費燃料が減る
15

課題２ 干渉点での通過順序の最適化

各機体の干渉点を通過する時間に対する評価関数値を示した評価関数曲線を導入し、時間間隔を確保した順序を混合
線形計画法で解いた。複数機の軌道最適化でその有効性を示す。
16

３-１ 評価関数曲線 １/３
軌道順序同時最適化問題（CCO機とCDO機の干渉）
des1, des2: 降下機1, 2
cl1, cl2 : 上昇機1, 2

干渉点• 降下機と上昇機がCDOとCCOで飛行．速度の変更は可能．
• CDOとCCOなので，干渉点における回避を時間間隔で回避．
• 干渉点における通過時間と終端における燃料消費量で評価．
17

３-１ 評価関数曲線 ２/３
評価関数曲線の導入
通過時間に対する燃料消費量（評価関数）を示した曲線を「評価関数曲線」と呼ぶ。

干渉点

評価関数曲線は最適な通過時間（tfOpt）で最小となる下に凸の曲線になる。
◆求めるためには最適制御問題を多数回解く必要がある。
計算量が多いので実時間システムには不利となる。
⇒高速化アルゴリズムで解決
18

３-１ 評価関数曲線 ３/３◆評価関数曲線の複数機へ適用
干渉点の時間間隔を保ちながら最適軌道を採用する

干渉点• 混合整数計画法を用いて、軌道と順序を同時に最適化することができる．

19

３-２ 軌道と順序の最適化 １/５
計算条件
des1, des2: 降下機1, 2
cl1, cl2 : 上昇機1, 2干渉点

• CCO機とCDO機が2機ずつ飛行• 干渉点にて最小時間（60秒）間隔を確保• 燃料消費量で評価

ケース１ 4機が同一機種（B772）
ケース２ CCO機の1機が異なる機種（B744）
ケース３ 出発時刻の遅延を許容
20

３-２ 軌道と順序の最適化 ２/５
ケース1（全てB772）

最適順序 des1 – cl1 – des2 – cl2

60 [s]
評価関数曲線

21

３-２ 軌道と順序の最適化 ３/５
計算条件
des1, des2: 降下機1, 2
cl1, cl2 : 上昇機1, 2

干渉点• CCO機とCDO機が2機ずつ飛行• 交差点にて最小時間（60秒）間隔を確保• 燃料消費量で評価

ケース１ 4機が同一機種（B772）
ケース２ CCO機の1機が異なる機種（B744）
ケース３ 出発時刻の遅延を許容
22

３-２ 軌道と順序の最適化 ４/５
ケース2（cl2がB744，軌道順序同時最適化）

最適順序
des1 – cl1 – cl2 – des2

評価関数曲線• 通過時間で燃料消費量が急変するB744を優先した順序．
23

３-２ 軌道と順序の最適化 ５/５
ケース3（比較： cl2がB744，順番指定）
燃料消費量 [lbs]
ケース2-1

31266.5

ケース2-231271.4- 4.9 [lbs]指定順序
des1 – cl1 – des2 – cl2

評価関数曲線• 軌道と順序を同時に最適化した方が燃料消費量が少なかった．
24

課題３ 実運用システムへの展開

実運用システムでは常に新しい離陸機と着陸機が領域に加わり、継続的な最適化が短時間で行う必要がある。継続的な最適化にはレシーディングホライズンを、計算時間の短縮にはホモトピー法を導入した。
25

４-１ 継続的最適化 １/２◆ レシーディングホライズン
後退ホライズン

x目標値
実行最適化
t

実行ホライズン• 後退ホライズンの区間の最適化を行い、得られた最適入力の実行ホライズンの部分区間のみを用いる。
26

４-１ 継続的最適化 ２/２
後退ホライズン

xx

目標値

目標値

実行最適化

tt

実行ホライズン••問題を分割することで解くべき最適化問題の計算時間を短縮
繰り返し最適化による無限時間の最適化
HpAirportt27

４-２ 計算量の低減 １/３
最適化計算の問題点••

評価関数曲線を得るためには複数点の評価関数値が必要．
1点の最適解を得るためには、膨大な繰り返し計算が必要。
分割点数
J

1

7263

4 5

評価関数曲線の計算- 各点は軌道最適化計算により得られる．
- 左図では7回の最適化計算が必要．

→ ホモトピー法の特徴を利用- 非線形連立方程式の数値解法- 近似を用いないので、計算時間が一定- 初期仮定解の近傍解を得るには有力△：初期仮定
○：最適解

28

４-２ 計算量の低減 ２/３連続変化（ホモトピー）法- ある最適解を初期仮定解に用いると，その近傍の解を少ない計算回数で求めることができる．
- 解は滑らかに変化する．
- 右図で，一番下の点のみ通常通り
求め，その点を初期解として他の点を求めると数値積分は16回．

J

積分回数1111111029

４-２ 計算量の低減 ３/３
計算時間比較 （軌道順序同時最適化 v.s. 順序総当たり）
総当たり計算計算時間[s]64ホモトピー法

２１機数
ホモトピー法 [s]総当たり計算 [s]

２４３機数
21.37
2.19

５3
1.85
15.00

42.49
106.5154.85
804.50

一般的PCを使用
30

５. 研究の成果 １/２
課題１ 他機との時間的・空間的間隔の確保⇒時空間座標系の導入で両者のバランスを考慮した軌道生成を可能とした。

課題２ 干渉点での通過順序の最適化⇒評価関数曲線の導入で混合線形計画問題として解くことができるようになった。
課題３ 実運用システムへの展開⇒レシーディングホライズンの導入で連続的な軌道生成を可能とし、ホモトピー法の利用で処理時間の短縮ができた。
31

５. 研究の成果 ２/２
公表論文[1] 虎谷大地、上野誠也、樋口丈浩、“空港周辺空域における上昇機のための最適な回避軌道に関する研究”、航空宇宙技術、Vol.14、pp.105-112、2015．
[2] D. Toratani, S. Ueno, T. Higuchi, “SimultaneousOptimization Method for Trajectory and Sequence
for Receding Horizon Guidance in Terminal Area”,
SICE Journal of CMSI, Vol. 8, No. 2, pp.144-153,
2015.
[3] 虎谷大地，上野誠也，樋口丈浩，“継続上昇・降下運航のための干渉回避手法に関する研究”，航空宇宙技術，（掲載可）．
その他：国際会議にて４件の口頭発表
32