WO2020129315A1 - プログラム、情報処理装置及び情報処理方法 - Google Patents

プログラム、情報処理装置及び情報処理方法 Download PDF

Info

Publication number
WO2020129315A1
WO2020129315A1 PCT/JP2019/033863 JP2019033863W WO2020129315A1 WO 2020129315 A1 WO2020129315 A1 WO 2020129315A1 JP 2019033863 W JP2019033863 W JP 2019033863W WO 2020129315 A1 WO2020129315 A1 WO 2020129315A1
Authority
WO
WIPO (PCT)
Prior art keywords
arrival time
target point
time
probability distribution
value
Prior art date
Application number
PCT/JP2019/033863
Other languages
English (en)
French (fr)
Inventor
賀宣 松野
Original Assignee
国立研究開発法人宇宙航空研究開発機構
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 国立研究開発法人宇宙航空研究開発機構 filed Critical 国立研究開発法人宇宙航空研究開発機構
Priority to US17/416,265 priority Critical patent/US20220051575A1/en
Publication of WO2020129315A1 publication Critical patent/WO2020129315A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
    • G08G5/0017—Arrangements for implementing traffic-related aircraft activities, e.g. arrangements for generating, displaying, acquiring or managing traffic information
    • G08G5/0021—Arrangements for implementing traffic-related aircraft activities, e.g. arrangements for generating, displaying, acquiring or managing traffic information located in the aircraft
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64C—AEROPLANES; HELICOPTERS
    • B64C13/00—Control systems or transmitting systems for actuating flying-control surfaces, lift-increasing flaps, air brakes, or spoilers
    • B64C13/02—Initiating means
    • B64C13/16—Initiating means actuated automatically, e.g. responsive to gust detectors
    • B64C13/18—Initiating means actuated automatically, e.g. responsive to gust detectors using automatic pilot
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64D—EQUIPMENT FOR FITTING IN OR TO AIRCRAFT; FLIGHT SUITS; PARACHUTES; ARRANGEMENT OR MOUNTING OF POWER PLANTS OR PROPULSION TRANSMISSIONS IN AIRCRAFT
    • B64D45/00—Aircraft indicators or protectors not otherwise provided for
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F1/00—Ground or aircraft-carrier-deck installations
    • B64F1/36—Other airport installations
    • 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
    • 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
    • G08G5/0052—Navigation or guidance aids for a single aircraft for cruising
    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
    • G08G5/0073—Surveillance aids
    • G08G5/0091—Surveillance aids for monitoring atmospheric conditions

Definitions

  • the present invention takes into consideration uncertainties during flight (such as weather prediction errors and airframe characteristic errors), and optimizes control inputs (speed, altitude, route, etc.) that satisfy the required arrival time at any target point on the flight route.
  • the present invention relates to a program, an information processing apparatus, and an information processing method for determining a corner, a cost index, a rise/fall rate, etc.).
  • an object of the present invention is to reduce the number of speed adjustments and the like to suppress the deterioration of fuel efficiency and comfort while suppressing the influence of uncertainty on the arrival time, an information processing device, and To provide an information processing method.
  • a program quantifies an uncertain predicted value that affects the arrival time of a moving body at a target point by using an actual measurement value to determine the uncertain predicted value.
  • the probability distribution of the estimated arrival time to the target point when the control input to the moving body that affects the arrival time at the target point is calculated, and the current control input is used, Calculated using the probability distribution of the uncertain predicted value, in the probability distribution of the estimated arrival time to the target point, the estimated arrival time for the arrival request time to the target point when using the current control input Is determined to be outside a first threshold value, and when the first statistic is outside the first threshold, the second statistic of the estimated arrival time with respect to the arrival request time is
  • the computer is caused to execute the step of determining the control input so as to be within the second threshold value.
  • the present invention by robustly managing the time with respect to the uncertain prediction value that affects the arrival time of the moving body at the target point, that is, the uncertainty with respect to the arrival time, the number of speed adjustments, etc. is reduced and the fuel consumption is reduced.
  • the influence of uncertainty on the arrival time can be suppressed while suppressing the deterioration of comfort.
  • the control input in the step of determining the control input, when the first statistic is outside the first threshold, the second statistic is within the second threshold. And the control input may be determined so as to optimize a predetermined cost.
  • the first and second statistics may be at least one of a probability, a moment, and a representative value in a probability distribution of the estimated arrival time with respect to the arrival request time. ..
  • the moving body is an aircraft
  • the uncertain prediction value is a weather prediction value and/or aircraft characteristics
  • the actual measurement value is relative to the weather prediction value.
  • the moving body is an aircraft
  • the control input to the aircraft is at least one of speed, altitude, route angle, cost index, and ascent/descent rate. Good.
  • An information processing apparatus quantifies an uncertain predicted value that affects the arrival time of a moving body at a target point by using a measured value to calculate a probability distribution of the uncertain predicted value.
  • An uncertainty estimation unit that is a control input to the moving body that affects the arrival time at the target point, and the probability distribution of the estimated arrival time to the target point when the current control input is used.
  • An estimated time-of-arrival estimation unit that calculates using the probability distribution of the uncertain predicted value, and a probability distribution of estimated time of arrival to the target point, up to the target point when the current control input is used.
  • An information processing method is to quantify an uncertain predicted value that affects the arrival time of a mobile body at a target point by using a measured value to calculate a probability distribution of the uncertain predicted value.
  • the probability distribution of the estimated arrival time to the target point which is a control input to the moving body that affects the arrival time to the target point, is calculated using the current control input.
  • the control input is determined so that
  • the present invention it is possible to reduce the number of speed adjustments and the like to suppress the deterioration of fuel efficiency and comfort while suppressing the influence of uncertainty on the arrival time.
  • FIG. 6 is a diagram for explaining the content (continuous type) of processing by the arrival time control unit 13.
  • FIG. 6 is a diagram for explaining the content (discrete type) of processing by the arrival time control unit 13.
  • 6 is a flowchart showing the operation of the information processing device 1.
  • 9 is a graph for comparing the number of speed changes of the information processing apparatus 1 according to the present embodiment with the related art.
  • 7 is a graph (No. 1) for comparing the time management accuracy of the information processing apparatus 1 according to the present embodiment with the related art.
  • 6 is a graph (No. 2) for comparing the time management accuracy of the information processing apparatus 1 according to the present embodiment with the related art.
  • FIG. 1 is a block diagram showing the configuration of an information processing apparatus according to an embodiment of the present invention.
  • the information processing apparatus 1 shown in FIG. 1 is used for, for example, an existing on-board device or a peripheral device obtained during flight through an existing on-board device or a data link by weather forecast correction and error generation technology (see Non-Patent Document 1). In-flight uncertainty is corrected and residual error is quantified using flight data.
  • the information processing apparatus 1 utilizes the correction information and the error information generated during the flight as described above, considers the influence of the quantified uncertainty, and ensures that the arrival time falls within the target time management accuracy. , Determine the optimum control input (speed, altitude, route angle, cost index, climb/descent rate, etc.). That is, the information processing device 1 takes into consideration uncertainty during flight (such as weather prediction error and airframe characteristic error), and the optimum control input (speed, speed, Altitude, route angle, cost index, climb/descent rate, etc.).
  • the onboard equipment is used to calculate the optimum control input, and the determined optimum control input is directly transmitted to the control system of the aircraft to automatically guide and control the aircraft, or, for example, the onboard electronic flight
  • the information is transmitted to the pilot through the display device of the bag (EFB) to assist the pilot in making a flight decision.
  • EFB display device of the bag
  • the information processing device 1 is typically composed of an application program of an aircraft on-board device or an EFB application program, and includes an uncertainty estimation unit 11 and an estimated arrival time estimation unit 12. , And an arrival time control unit 13.
  • the uncertainty estimation unit 11 quantifies the uncertain predicted value in flight, which affects the arrival time of the aircraft at the target point, using the actual measurement value to calculate the probability distribution of the uncertain predicted value.
  • the uncertain predicted value that affects the arrival time is typically weather forecast, more specifically, for example, predicted wind speed, wind direction, and temperature
  • the actual measurement value is, for example, meteorological observation. This is a value, more specifically, the wind speed, wind direction, and temperature actually observed by the aircraft.
  • the uncertain predicted value includes not only the weather predicted value but also aircraft characteristics. More specifically, the airframe characteristics are, for example, fuel flow rate, engine parameters, thrust, and aerodynamic force.
  • the actual measurement value is a machine body measurement value for the machine body characteristic. In the following description, only the weather forecast will be described as an example of the predicted value, but a combination with the airframe characteristic or only the airframe characteristic may be used as the uncertain predicted value.
  • the control input to the aircraft is typically the flight speed of the aircraft. Control inputs to the aircraft include the flight speed of the aircraft, altitude, route angle, cost index, and climb/descent rate. Further, the control input to the aircraft may be two or more of these. Furthermore, it may be time-varying or time-invariant.
  • Fig. 2 shows an example of an uncertain predicted value, that is, a probability distribution of uncertainties.
  • the probability distribution of the uncertainties is generated by quantifying the weather forecast error, which is the difference between the weather forecast value and the actual observed value. For example, using the meteorological observation value (x_obs), the probability (p(x_1), p(x_2),..., Of the weather forecast (x_1, x_2,..., x_N) of the N patterns input in advance is used. p(x_N)) is calculated.
  • the continuous type probability distribution is illustrated in FIG. 2, a discrete type probability distribution may be generated as shown in FIG.
  • the estimated arrival time estimation unit 12 uses the probability distribution of the uncertainties shown in FIG. 2 or FIG. 3 to estimate the estimated arrival time to the target point when the flight route is planned at the current flight speed. Calculate the probability distribution.
  • Figure 4 shows an example of the probability distribution of estimated arrival times.
  • the estimated time of arrival (ETA) when flying at the current speed (v) is calculated for each of the above N patterns of weather forecast (x_1, x_2,..., x_N).
  • ETA_1, ETA_2,..., ETA_N and their probabilities (p′(ETA_1), p′(ETA_2),..., P′(ETA_N)) are calculated.
  • the estimated arrival time when flying at the speed v under the weather forecast x_1 is ETA_1
  • its probability p′(ETA_1) is the probability p(x_1) that the weather forecast is x_1.
  • the continuous type probability distribution is illustrated in FIG. 4, a discrete type probability distribution may be generated as shown in FIG.
  • the arrival time control unit 13 determines whether the statistical amount of the estimated arrival time with respect to the arrival request time at the current flight speed v is outside a predetermined threshold value based on the probability distribution of the estimated arrival time shown in FIG. 4 or 5.
  • the flight speed is determined so that the statistic of the estimated arrival time with respect to the arrival request time is within the predetermined threshold.
  • the probability Pr that satisfies the arrival request time (RTA) within the allowable range (-tol1 to +tol2). (The area of the diagonal line in the figure) is used, and the probability Pr is used for the threshold determination.
  • the arrival time control unit 13 maximizes the probability Pr, that is, maximizes the area of the diagonal line in the figure, as shown in the lower diagram of FIG. To determine the flight speed.
  • FIG. 6 exemplifies a continuous type probability distribution, it can be similarly implemented even if it is a discrete type probability distribution as shown in FIG. 7.
  • the arrival time control unit 13 sets the flight speed to v, the estimated arrival time to ETA, the arrival request time to RTA, the allowable range to -tol to +tol2, the probability Pr, and the probability distribution of the estimated arrival time to the target point to p.
  • '(ETA) if the probability Pr is less than or equal to a predetermined threshold value, the set speed v opt that maximizes the probability Pr is maximized by the following formula (1) when arriving at the arrival request time, If it is desired to arrive by the arrival request time, it is maximized by the following equation (2), and if it is desired to arrive after the arrival request time, it is maximized by the following equation (3).
  • FIG. 8 is a flowchart showing the operation of the information processing device 1 configured as above.
  • the uncertainty estimation unit 11 uses the meteorological observation value and the aircraft measurement value during the flight to obtain the probability distribution (see FIG. 2 or FIG. 3) of the uncertainty element during the flight (the weather forecast, the aircraft characteristic, etc.). It is generated (step 801).
  • the estimated arrival time estimation unit 12 uses the probability distribution of the uncertainties (see FIG. 2 or 3) to obtain the current control input (speed, altitude, route angle, cost index, ascent/descent rate, etc.).
  • the probability distribution (see FIG. 4 or FIG. 5) of the estimated arrival time at the target point when the flight is performed in step S802 is generated (step 802). To generate.
  • the arrival time control unit 13 determines the threshold value of the estimated arrival time of the current control input (step 803), and if it is outside the threshold value, optimizes the control input that satisfies the arrival request time within the allowable range. (Step 804).
  • the flight speed is maximized in consideration of the influence of the error so that the arrival time falls within the time management accuracy. Can be determined.
  • the arrival time control unit according to the present invention is different from the first statistic when the first statistic of the estimated arrival time with respect to the arrival request time is outside the first threshold which is a predetermined threshold, or The control input is determined such that the second statistic of the estimated arrival time with respect to the arrival request time, which is the first statistic, is different from the first threshold value or is within the second threshold value which is the first threshold value. It is a thing.
  • the arrival time control unit 13 determines the set speed so as to maximize the probability Pr.
  • the arrival time control unit 13 determines that the probability Pr is equal to or less than a predetermined threshold value.
  • the control input may be determined so as to optimize a predetermined cost under the condition that the probability Pr is within a predetermined threshold.
  • the optimization of the predetermined cost includes, for example, minimizing the fuel consumption amount, but it may be not only the fuel consumption amount but also another cost factor, or a combination thereof.
  • the probability Pr that satisfies the arrival request time (RTA) within the allowable range (-tol1 to +tol2) is used as the statistic.
  • Value, variance, skewness, etc.), representative value (median value, percentile value, mode value, etc.), or the like, or a combination thereof may be used.
  • the arrival time control unit uses the following equation (4) when the expected value E[ETA] of the estimated arrival time (ETA) is less than or equal to a predetermined threshold value. Then, the control input is determined so as to minimize the square of the difference between the expected value E[ETA] and the arrival request time RTA. ..
  • the arrival time control unit optimizes the predetermined cost under the condition that the expected value E[ETA] is within the predetermined threshold when the expected value E[ETA] of the estimated arrival time (ETA) is less than or equal to the predetermined threshold.
  • the control input is determined so that For example, consider the fuel consumption up to the target point as the predetermined cost. If you want to arrive at the arrival request time: Minimize fuel consumption Constraint: RTA-tol1 ⁇ E[ETA] ⁇ RTA + tol2 If you want to arrive by the arrival request time: Minimize fuel consumption Constraint: E[ETA] ⁇ RTA + tol2 If you want to arrive after the arrival request time: Minimize fuel consumption Constraint: RTA-tol1 ⁇ E [ETA]
  • the present invention can be applied not only to implementation on an aircraft but also to a ground air traffic control system.
  • the feasible arrival time can be determined and transmitted to each aircraft.
  • the present invention can be applied not only to aircraft but also to ships, trains, automobiles, etc., as long as the operation manages the arrival time.
  • uncertainties such as tidal currents in the case of ships, and uncertainties such as traffic jams in the case of automobiles.
  • the optimum speed can be determined by applying the present invention in consideration of such uncertainties in the operation.
  • Information processing device 11 Uncertainty estimation unit 12: Estimated arrival time estimation unit 13: Arrival time control unit

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Mechanical Engineering (AREA)
  • Traffic Control Systems (AREA)

Abstract

【課題】速度調整回数等を低減して燃費と快適性の悪化を抑えつつ到着時刻に対する不確定性の影響を抑えること。 【解決手段】情報処理装置１は、移動体の目標地点への到着時刻に影響する不確定な予測値を、実測値を用いて定量化して前記不確定な予測値の確率分布を算出する不確定性推定部１１と、前記目標地点への到着時刻に影響する移動体への制御入力であって、現在の前記制御入力を用いた場合の前記目標地点までの到着予定時刻の確率分布を、前記不確定な予測値の確率分布を用いて算出する到着予定時刻推定部１２と、前記目標地点までの到着予定時刻の確率分布において、前記現在の制御入力を用いた場合の前記目標地点までの到着要求時刻に対する到着予定時刻の第１の統計量が第１の閾値外かどうかを判定し、前記第１の統計量が前記第１の閾値外の場合に、前記到着要求時刻に対する到着予定時刻の第２の統計量が第２の閾値内となるように前記制御入力を決定する到着時刻制御部１３とを具備する。

Description

プログラム、情報処理装置及び情報処理方法
　本発明は、例えば飛行中の不確定性（気象予測誤差や機体特性誤差等）を考慮し、飛行経路上の任意の目標地点への到着要求時刻を満たす最適な制御入力（速度、高度、経路角、コストインデックス、上昇・降下率等）を決定するプログラム、情報処理装置及び情報処理方法に関する。
　航空交通需要の増大に伴い高密度運航を行うためには、時間管理を導入した運航、つまり航空機が目標地点を要求時刻に通過する必要があり、到着予定時刻を予測し到着要求時刻を満たすように飛行速度を制御する到着時刻制御技術が要求される（特許文献１参照）。
　また、飛行中に生じる気象予測誤差等の不確実性の影響から到着予定時刻を正確に予測することは困難なため、到着時刻の不確実性を定量化し、不確実性を使用した到着時刻制御技術も提案されている（特許文献２参照）。
米国特許第５１２１３２５号公報 米国特許第８１５０５８８号公報
「Ｎｏｗｃａｓｔｉｎｇ　ａｌｇｏｒｉｓｍ　ｆｏｒ　ｗｉｎｄ　ｆｉｅｌｄｓ　ｕｓｉｎｇ　ｅｎｓｅｍｂｌｅ　ｆｏｒｅｃａｓｔｉｎｇ　ａｎｄ　ａｉｒｃｒａｆｔ　ｆｌｉｇｈｔ　ｄａｔａ」　Ｒｙｏｔａ　Ｋｉｋｕｃｈｉ　ｅｔ　ａｌ．２０１８　Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ　Ａｐｌｉｃａｔｉｏｎｓ
　特許文献１に記載された技術では、飛行中の不確実性の影響から速度調整が高頻度化し加減速を繰り返すため、燃費と快適性が悪化する。また、特許文献２に記載された技術でも、事前に定めた最低使用可能速度と最高使用可能速度との切り替えのみを想定したバンバン制御であり、燃費と快適性が悪化する。そのため、特許文献１及び２に記載された技術では、実用性に乏しかった。
　以上のような事情に鑑み、本発明の目的は、速度調整回数等を低減して燃費と快適性の悪化を抑えつつ到着時刻に対する不確定性の影響を抑えることができるプログラム、情報処理装置及び情報処理方法を提供することにある。
　上記目的を達成するため、本発明の一形態に係るプログラムは、移動体の目標地点への到着時刻に影響する不確定な予測値を、実測値を用いて定量化して前記不確定な予測値の確率分布を算出し、前記目標地点への到着時刻に影響する移動体への制御入力であって、現在の前記制御入力を用いた場合の前記目標地点までの到着予定時刻の確率分布を、前記不確定な予測値の確率分布を用いて算出し、前記目標地点までの到着予定時刻の確率分布において、前記現在の制御入力を用いた場合の前記目標地点までの到着要求時刻に対する到着予定時刻の第１の統計量が第１の閾値外かどうかを判定し、前記第１の統計量が前記第１の閾値外の場合に、前記到着要求時刻に対する到着予定時刻の第２の統計量が第２の閾値内となるように前記制御入力を決定するステップをコンピュータに実行させる。
　本発明では、移動体の目標地点への到着時刻に影響する不確定な予測値、つまり到着時刻に対する不確定性に対してロバストに時間管理を行うことで、速度調整回数等を低減して燃費と快適性の悪化を抑えつつ到着時刻に対する不確定性の影響を抑えることができる。
　本発明の一形態に係るプログラムでは、前記制御入力を決定するステップは、前記第１の統計量が前記第１の閾値外の場合に、前記第２の統計量が前記第２の閾値内となりかつ所定のコストを最適化するように前記制御入力を決定してもよい。
　本発明の一形態に係るプログラムでは、前記第１及び第２の統計量は、前記到着要求時刻に対する前記到着予定時刻の確率分布における確率、モーメント及び代表値のうち少なくとも１つであってもよい。
　本発明の一形態に係るプログラムでは、前記移動体は、航空機であり、前記不確定な予測値は、気象予測値及び／又は機体特性であり、前記実測値は、前記気象予測値に対しては気象観測値、前記機体特性に対しては機体計測値であってもよい。
　本発明の一形態に係るプログラムでは、前記移動体は、航空機であり、前記航空機への制御入力は、速度、高度、経路角、コストインデックス及び上昇・降下率のうち少なくとも一つであってもよい。
　本発明の一形態に係る情報処理装置は、移動体の目標地点への到着時刻に影響する不確定な予測値を、実測値を用いて定量化して前記不確定な予測値の確率分布を算出する不確定性推定部と、前記目標地点への到着時刻に影響する移動体への制御入力であって、現在の前記制御入力を用いた場合の前記目標地点までの到着予定時刻の確率分布を、前記不確定な予測値の確率分布を用いて算出する到着予定時刻推定部と、前記目標地点までの到着予定時刻の確率分布において、前記現在の制御入力を用いた場合の前記目標地点までの到着要求時刻に対する到着予定時刻の第１の統計量が第１の閾値外かどうかを判定し、前記第１の統計量が前記第１の閾値外の場合に、前記到着要求時刻に対する到着予定時刻の第２の統計量が第２の閾値内となるように前記制御入力を決定する到着時刻制御部とを具備する。
　本発明の一形態に係る情報処理方法は、移動体の目標地点への到着時刻に影響する不確定な予測値を、実測値を用いて定量化して前記不確定な予測値の確率分布を算出し、前記目標地点への到着時刻に影響する移動体への制御入力であって、現在の前記制御入力を用いた場合の前記目標地点までの到着予定時刻の確率分布を、前記不確定な予測値の確率分布を用いて算出し、前記目標地点までの到着予定時刻の確率分布において、前記現在の制御入力を用いた場合の前記目標地点までの到着要求時刻に対する到着予定時刻の第１の統計量が第１の閾値外かどうかを判定し、前記第１の統計量が前記第１の閾値外の場合に、前記到着要求時刻に対する到着予定時刻の第２の統計量が第２の閾値内となるように前記制御入力を決定する。
　本発明によれば、速度調整回数等を低減して燃費と快適性の悪化を抑えつつ到着時刻に対する不確定性の影響を抑えることができる。
本発明の一実施形態に係る情報処理装置１の構成を示すブロック図である。 不確定要素の確率分布（連続型の確率分布）の一例を示すグラフである。 不確定要素の確率分布（離散型の確率分布）の一例を示すグラフである。 到着予定時刻の確率分布（連続型の確率分布）の一例を示すグラフである。 到着予定時刻の確率分布（離散型の確率分布）の一例を示すグラフである。 到着時刻制御部１３による処理の内容（連続型）を説明するための図である。 到着時刻制御部１３による処理の内容（離散型）を説明するための図である。 情報処理装置１の動作を示すフローチャートである。 従来技術と本実施形態に係る情報処理装置１の速度変更回数を比較するためのグラフである。 従来技術と本実施形態に係る情報処理装置１の時間管理精度を比較するためのグラフ（その１）である。 従来技術と本実施形態に係る情報処理装置１の時間管理精度を比較するためのグラフ（その２）である。
　以下、図面を参照しながら、本発明の実施形態を説明する。
　図１は、本発明の一実施形態に係る情報処理装置の構成を示すブロック図である。　
　図１に示す情報処理装置１は、例えば気象予測補正及び誤差生成技術（非特許文献１参照）等により、既存の機上装置もしくはデータリンクを介して飛行中に得られる自機や周辺機の飛行データ等を用いて、飛行中の不確実性を補正すると共に、残留誤差を定量化する。
　情報処理装置１は、このように飛行中に生成される補正情報及び誤差情報を活用し、定量化された不確実性の影響を考え、到着時刻が目標とする時間管理精度内に収まるように、最適な制御入力（速度、高度、経路角、コストインデックス、上昇・降下率等）を決定する。すなわち、情報処理装置１は、飛行中の不確定性（気象予測誤差や機体特性誤差等）を考慮し、飛行経路上の任意の目標地点への到着要求時刻を満たす最適な制御入力（速度、高度、経路角、コストインデックス、上昇・降下率等）を決定する。
　最適な制御入力の算出には機上機器を使用し、決定された最適な制御入力は、航空機の制御システムに直接伝達され航空機を自動で誘導制御するか、もしくは、例えば搭載したエレクトロニック・フライト・バッグ（ＥＦＢ）の表示機器を介してパイロットに伝達し、パイロットの運航判断を支援する。
　つまり、情報処理装置１は、典型的には、航空機の機上機器のアプリケーションプログラムからなるもの、或いはＥＦＢのアプリケーションプログラムからなるもので、不確定性推定部１１と、到着予定時刻推定部１２と、到着時刻制御部１３とを有する。
　不確定性推定部１１は、航空機の目標地点への到着時刻に影響する飛行中の不確定な予測値を、実測値を用いて定量化して不確定な予測値の確率分布を算出する。
　ここで、到着時刻に影響する不確定な予測値とは、典型的には、気象予測であり、より詳細には例えば予測される風速や風向、気温であり、実測値とは、例えば気象観測値であり、より詳細には当該航空機で実際に観測された風速や風向、気温である。前記の不確定な予測値は、気象予測値の他、機体特性等もある。機体特性とは、より詳細には、例えば燃料流量やエンジンパラメータ、推力、空気力である。前記の実測値は、前記の機体特性に対する機体計測値である。以下の説明では、予測値として気象予測だけを例にして説明するが、機体特性等との組み合わせ、或いは機体特性等のみを不確定な予測値としてもよい。また、航空機への制御入力とは、典型的には、航空機の飛行速度である。航空機への制御入力は、航空機の飛行速度の他、高度や経路角、コストインデックス、上昇・降下率等がある。また、航空機への制御入力は、これらのうち二つ以上であってもよい。さらに、時変でも時不変であってもよい。
　図２に不確定な予測値、つまり不確定要素の確率分布の一例を示す。不確定要素の確率分布は、気象予測値と実測値である気象観測値との差である気象予測誤差を定量化して生成される。例えば、気象観測値(ｘ＿ｏｂｓ)を用いて、事前に入力されたＮパターンの気象予測（ｘ＿１，ｘ＿２，・・・，ｘ＿Ｎ）の確率（ｐ（ｘ＿１），ｐ（ｘ＿２），・・・，ｐ（ｘ＿Ｎ））を算出する。なお、図２には連続型の確率分布を例示しているが、図３に示すように離散型の確率分布を生成するようにしてもよい。
　到着予定時刻推定部１２は、図２又は図３に示した不確定要素の確率分布を用いて、現在の飛行速度で計画している飛行経路を通った場合の目標地点までの到着予定時刻の確率分布を算出する。
　図４に到着予定時刻の確率分布の一例を示す。例えば、現在の速度（ｖ）で飛行した場合の到着予定時刻（ＥＴＡ）を、上記のＮパターンの気象予測（ｘ＿１，ｘ＿２，・・・，ｘ＿Ｎ）各々の場合で算出し（この算出結果を（ＥＴＡ＿１，ＥＴＡ＿２，・・・，ＥＴＡ＿Ｎ）とする）、それらの確率（ｐ'（ＥＴＡ＿１），ｐ'（ＥＴＡ＿２），・・・，ｐ'（ＥＴＡ＿Ｎ））を算出する。ここで、気象予測ｘ＿１のもと速度ｖで飛行した場合の到着予定時刻はＥＴＡ＿１であり、その確率ｐ'（ＥＴＡ＿１）は気象予測がｘ＿１となる確率ｐ（ｘ＿１）となる。なお、図４には連続型の確率分布を例示しているが、図５に示すように離散型の確率分布を生成するようにしてもよい。
　到着時刻制御部１３は、図４又は図５に示した到着予定時刻の確率分布に基づいて、現在の飛行速度ｖでの到着要求時刻に対する到着予定時刻の統計量が所定の閾値外かどうかを判定し、統計量が所定の閾値外の場合に、到着要求時刻に対する到着予定時刻の統計量が所定の閾値内となるように飛行速度を決定する。例えば、図６の上図に示すように、到着要求時刻（ＲＴＡ）に対する到着予定時刻（ＥＴＡ）の統計量として、到着要求時刻（ＲＴＡ）の許容範囲内（－ｔｏｌ１～＋ｔｏｌ２）を満たす確率Ｐｒ（図中の斜線の面積）を用い、その確率Ｐｒを閾値判定に用いる。
　到着時刻制御部１３は、その確率Ｐｒが所定の閾値以下の場合には、図６の下図に示すように、確率Ｐｒが最大化するように、すなわち図中の斜線の面積が最大化するように飛行速度を決定する。なお、図６には連続型の確率分布を例示しているが、図７に示すように離散型の確率分布であっても同様に実施が可能である。
　到着時刻制御部１３は、例えば、飛行速度をｖ、到着予定時刻をＥＴＡ、到着要求時刻をＲＴＡ、許容範囲を－ｔｏｌ１～＋ｔｏｌ２、確率をＰｒ、目標地点までの到着予定時刻の確率分布をｐ'（ＥＴＡ）としたとき、確率Ｐｒが所定の閾値以下の場合に、確率Ｐｒを最大化する設定速度ｖｏｐｔを、到着要求時刻に到着したい場合には以下の式（１）で最大化、到着要求時刻までに到着したい場合には以下の式（２）で最大化及び到着要求時刻以降に到着したい場合には以下の式（３）で最大化して決定する。
Figure JPOXMLDOC01-appb-I000001
 　　　　　　　　　　　　　　　　　　　　　　　　　　　　式（１）
Figure JPOXMLDOC01-appb-I000002
　　　　　　　　　　　　　　　　　　　　　　　　　　　　式（２）
Figure JPOXMLDOC01-appb-I000003
　　　　　　　　　　　　　　　　　　　　　　　　　　　　　式（３）
　図８はこのように構成された情報処理装置１の動作を示すフローチャートである。　
　まず、不確定性推定部１１が、飛行中の気象観測値や機体計測値を用いて、飛行中の不確定要素（気象予測や機体特性等）の確率分布（図２又は図３参照）を生成する（ステップ８０１）。　
　次に、到着予定時刻推定部１２が、不確定要素の確率分布（図２又は図３参照）を用いて、現在の制御入力（速度、高度、経路角、コストインデックス、上昇・降下率等）で飛行した場合の目標地点への到着予定時刻の確率分布（図４又は図５参照）を生成する（ステップ８０２）。
を生成する。　
　次に、到着時刻制御部１３が、現在の制御入力での到着予定時刻の閾値を判定し（ステップ８０３）、閾値外の場合には、到着要求時刻を許容範囲内で満たす制御入力を最適化する（ステップ８０４）。
　図９に示すように、これまで提案された技術（特許文献１、２等参照）では、速度調整回数が高頻度化（図９のａ）するため、現実の航空機の運航では使用され難かったが、本実施形態に係る情報処理装置１は、不確定性に対してロバストに時間管理を行うことができるので、速度調整回数を低減し（図９のｂ）快適性や燃費が改善され、実用性が高くなる。また、本実施形態に係る情報処理装置１は、到着時間誤差においても、図１０に示すように、これまで提案された技術（図１０のａ）と比べて向上する（図１０のｂ）。ちなみに、図１１に示すように、これまで提案された技術（図１１のａ）では、現在速度で飛行を続けると、到着時刻は大きくずれ、時間管理精度内（例えば目標時刻精度±５秒）に収まる確率が低下するが、本実施形態に係る情報処理装置１（図１１のｂ）では、誤差の影響を考慮し、到着時刻が時間管理精度内に収まる確率を最大化するように飛行速度を決定することができる。
　本発明は、上記の実施形態には限定されず、本発明の技術思想の範囲内で様々な変形や応用による実施が可能であり、その実施の範囲も本発明の技術的範囲に属するものである。　
　例えば、本発明に係る到着時刻制御部は、到着要求時刻に対する到着予定時刻の第１の統計量が所定の閾値である第１の閾値外の場合に、第１の統計量とは異なる、或いは第１の統計量である到着要求時刻に対する到着予定時刻の第２の統計量が第１の閾値とは異なる、或いは第１の閾値である第２の閾値内となるように制御入力を決定するものである。
　従って、上記の実施形態では、到着時刻制御部１３は、確率Ｐｒを最大化するように設定速度を決定していたが、例えば到着時刻制御部１３は、その確率Ｐｒが所定の閾値以下の場合には、確率Ｐｒを所定の閾値内にするという条件下で所定のコストを最適化するように制御入力を決定してもよい。所定のコストを最適化とは、例えば燃料消費量を最小化することが挙げられるが、燃料消費量だけでなく他のコスト要因であってもよく、或いはこれらの組み合わせであってもよい。
　また、上記の実施形態では、統計量として、到着要求時刻（ＲＴＡ）の許容範囲内（－ｔｏｌ１～＋ｔｏｌ２）を満たす確率Ｐｒを用いていたが、統計量としては、確率Ｐｒ以外にモーメント（期待値や分散、歪度等）、代表値（中央値、パーセンタイル値、最頻値等）等を用いてもよく、或いはこれらの組み合わせを用いてもよい。例えば、統計量として期待値を用いた場合に、到着時刻制御部は、到着予定時刻（ETA）の期待値E[ETA]が所定の閾値以下の場合に、下記の式（４）に示すように期待値E[ETA]と到着要求時刻RTAとの差の二乗を最小化するように制御入力を決定する。　
　(E[ETA]－RTA)２　　　　　　　　　　　　　　　　　　　式（４）
　到着時刻制御部は、到着予定時刻（ETA）の期待値E[ETA]が所定の閾値以下の場合に、期待値E[ETA]を所定の閾値内にするという条件下で所定のコストを最適化するように制御入力を決定する。例えば、所定のコストとして、目標地点までの燃料消費量を考える。　
　到着要求時刻に到着したい場合：燃料消費量を最小化
　　　　　　　　　　制約条件：RTA－tol1≦E[ETA]≦RTA＋tol2
到着要求時刻までに到着したい場合：燃料消費量を最小化
　　　　　　　　　　制約条件：E[ETA]≦RTA＋tol2
到着要求時刻以降に到着したい場合：燃料消費量を最小化
　　　　　　　　　　制約条件：RTA－tol1≦E[ETA]
　本発明は、航空機上での実施に限らず、地上の航空管制システムへも適用できる。この場合に、地上から機上に対して到着時刻を指定する際、地上で得られる飛行中の不確実性情報及び各機体の使用可能速度範囲をもとに、各機体が調整可能な到着時刻幅を算出することで、実現可能な到着時刻を決定し各機上へ伝達することができる。このように管制情報処理システムを統合することにより、複数の航空機が整然と飛行すると共に、位置や時間の予見を向上し、航空機の間隔を短縮することが可能となり、混雑空域において効率的な降下方式や高密度運航を実現でき、また運航者の希望する自由度の高い飛行機道を実現でき、更に運航前からの協調的な軌道調整により遅延の最小化を図ることができる。
　また、本発明は、到着時刻が管理される運航であれば、航空機に限らず、船舶や電車、自動車等に適用可能である。例えば、船舶であれば潮流等の不確実性が存在し、自動車であれば渋滞等の不確実性がある。このような運航等における不確実性を考えた上で、本発明を適用し、最適速度を決定することができる。
　１　　　　：情報処理装置
１１　　　：不確定性推定部
１２　　　：到着予定時刻推定部
１３　　　：到着時刻制御部

Claims (7)

  1. 　移動体の目標地点への到着時刻に影響する不確定な予測値を、実測値を用いて定量化して前記不確定な予測値の確率分布を算出し、
    　前記目標地点への到着時刻に影響する移動体への制御入力であって、現在の前記制御入力を用いた場合の前記目標地点までの到着予定時刻の確率分布を、前記不確定な予測値の確率分布を用いて算出し、
    　前記目標地点までの到着予定時刻の確率分布において、前記現在の制御入力を用いた場合の前記目標地点までの到着要求時刻に対する到着予定時刻の第１の統計量が第１の閾値外かどうかを判定し、
    　前記第１の統計量が前記第１の閾値外の場合に、前記到着要求時刻に対する到着予定時刻の第２の統計量が第２の閾値内となるように前記制御入力を決定する
    　ステップをコンピュータに実行させるプログラム。
  2. 　前記制御入力を決定するステップは、前記第１の統計量が前記第１の閾値外の場合に、前記第２の統計量が前記第２の閾値内となりかつ所定のコストを最適化するように前記制御入力を決定する
    　請求項１に記載のプログラム。
  3. 　前記第１及び第２の統計量は、前記到着要求時刻に対する前記到着予定時刻の確率分布における確率、モーメント及び代表値のうち少なくとも１つである
    　請求項１又は２に記載のプログラム。
  4. 　前記移動体は、航空機であり、
    　前記不確定な予測値は、気象予測値及び／又は機体特性であり、前記実測値は、前記気象予測値に対しては気象観測値、前記機体特性に対しては機体計測値である
    　請求項１～３のうちいずれか一項に記載のプログラム。
  5. 　前記移動体は、航空機であり、
    　前記航空機への制御入力は、速度、高度、経路角、コストインデックス及び上昇・降下率のうち少なくとも１つである
    　請求項１～４のうちいずれか一項に記載のプログラム。
  6. 　移動体の目標地点への到着時刻に影響する不確定な予測値を、実測値を用いて定量化して前記不確定な予測値の確率分布を算出する不確定性推定部と、
    　前記目標地点への到着時刻に影響する移動体への制御入力であって、現在の前記制御入力を用いた場合の前記目標地点までの到着予定時刻の確率分布を、前記不確定な予測値の確率分布を用いて算出する到着予定時刻推定部と、
    　前記目標地点までの到着予定時刻の確率分布において、前記現在の制御入力を用いた場合の前記目標地点までの到着要求時刻に対する到着予定時刻の第１の統計量が第１の閾値外かどうかを判定し、前記第１の統計量が前記第１の閾値外の場合に、前記到着要求時刻に対する到着予定時刻の第２の統計量が第２の閾値内となるように前記制御入力を決定する到着時刻制御部と
    　を具備する情報処理装置。
  7. 　移動体の目標地点への到着時刻に影響する不確定な予測値を、実測値を用いて定量化して前記不確定な予測値の確率分布を算出し、
    　前記目標地点への到着時刻に影響する移動体への制御入力であって、現在の前記制御入力を用いた場合の前記目標地点までの到着予定時刻の確率分布を、前記不確定な予測値の確率分布を用いて算出し、
    　前記目標地点までの到着予定時刻の確率分布において、前記現在の制御入力を用いた場合の前記目標地点までの到着要求時刻に対する到着予定時刻の第１の統計量が第１の閾値外かどうかを判定し、
    　前記第１の統計量が前記第１の閾値外の場合に、前記到着要求時刻に対する到着予定時刻の第２の統計量が第２の閾値内となるように前記制御入力を決定する
    　情報処理方法。
PCT/JP2019/033863 2018-12-20 2019-08-29 プログラム、情報処理装置及び情報処理方法 WO2020129315A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US17/416,265 US20220051575A1 (en) 2018-12-20 2019-08-29 Program, information processing apparatus, and information processing method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018-237900 2018-12-20
JP2018237900A JP7228229B2 (ja) 2018-12-20 2018-12-20 プログラム、情報処理装置及び情報処理方法

Publications (1)

Publication Number Publication Date
WO2020129315A1 true WO2020129315A1 (ja) 2020-06-25

Family

ID=71100849

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/033863 WO2020129315A1 (ja) 2018-12-20 2019-08-29 プログラム、情報処理装置及び情報処理方法

Country Status (3)

Country Link
US (1) US20220051575A1 (ja)
JP (1) JP7228229B2 (ja)
WO (1) WO2020129315A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220051573A1 (en) * 2020-08-14 2022-02-17 Honeywell International Inc. Method and system for initiating and managing required time of arrival constraints for an aircraft for all phases of flight
JP2024042218A (ja) * 2022-09-15 2024-03-28 株式会社日立製作所 飛行体位置監視システム、および、飛行体位置監視方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012510108A (ja) * 2008-11-25 2012-04-26 ゼネラル・エレクトリック・カンパニイ 到着時刻不確定性を使用する到着時刻制御の方法およびシステム
JP2012174270A (ja) * 2011-02-22 2012-09-10 General Electric Co <Ge> 航空交通を管理する方法およびシステム
JP2013177120A (ja) * 2012-02-27 2013-09-09 Ge Aviation Systems Llc 飛行計画を飛行中に調節するための方法
US8862287B1 (en) * 2010-05-17 2014-10-14 The Boeing Company Four dimensional trajectory based operation flight plans
JP2017503272A (ja) * 2013-12-31 2017-01-26 ザ・ボーイング・カンパニーＴｈｅ Ｂｏｅｉｎｇ Ｃｏｍｐａｎｙ 航空機トラジェクトリの画定および予測のためのシステムおよび方法
US20170249849A1 (en) * 2016-02-29 2017-08-31 The Boeing Company Method and electronic device for optimizing a flight trajectory of an aircraft subject to time of arrival control constraints

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8781651B2 (en) * 2012-09-21 2014-07-15 Georgia Tech Research Corporation Systems and methods providing a fuel-efficient RTA implementation with uncertain winds
US9193442B1 (en) * 2014-05-21 2015-11-24 Rockwell Collins, Inc. Predictable and required time of arrival compliant optimized profile descents with four dimensional flight management system and related method
US10718885B2 (en) * 2016-08-12 2020-07-21 The Boeing Company System and methods for aviation winds uncertainty prediction

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012510108A (ja) * 2008-11-25 2012-04-26 ゼネラル・エレクトリック・カンパニイ 到着時刻不確定性を使用する到着時刻制御の方法およびシステム
US8862287B1 (en) * 2010-05-17 2014-10-14 The Boeing Company Four dimensional trajectory based operation flight plans
JP2012174270A (ja) * 2011-02-22 2012-09-10 General Electric Co <Ge> 航空交通を管理する方法およびシステム
JP2013177120A (ja) * 2012-02-27 2013-09-09 Ge Aviation Systems Llc 飛行計画を飛行中に調節するための方法
JP2017503272A (ja) * 2013-12-31 2017-01-26 ザ・ボーイング・カンパニーＴｈｅ Ｂｏｅｉｎｇ Ｃｏｍｐａｎｙ 航空機トラジェクトリの画定および予測のためのシステムおよび方法
US20170249849A1 (en) * 2016-02-29 2017-08-31 The Boeing Company Method and electronic device for optimizing a flight trajectory of an aircraft subject to time of arrival control constraints

Also Published As

Publication number Publication date
JP2020101879A (ja) 2020-07-02
JP7228229B2 (ja) 2023-02-24
US20220051575A1 (en) 2022-02-17

Similar Documents

Publication Publication Date Title
US8788189B2 (en) Aircraft control system and method for reaching a waypoint at a required time of arrival
JP5289581B2 (ja) 到着時刻不確定性を使用する到着時刻制御の方法およびシステム
US8165734B2 (en) Method and device for guiding an aircraft
JP6034130B2 (ja) 航空機パラメータを推定するための方法およびシステム
EP3232292B1 (en) Aircraft systems and methods with multiple sap speed profiles
US9330574B1 (en) In-flight generation of RTA-compliant optimal profile descent paths
US7283895B2 (en) Method and device for calculating a target speed for an aircraft
US20140277853A1 (en) System and method for determining aircraft operational parameters and enhancing aircraft operation
US20120004837A1 (en) method and system of controlling air traffic
US20100131125A1 (en) Method for assisting in the management of the flight of an aircraft in order to keep to a time constraint
CN104077928B (zh) 适用于航路外无规律飞行的通航过点位置信息更新方法
WO2020129315A1 (ja) プログラム、情報処理装置及び情報処理方法
US9507053B2 (en) Using aircraft trajectory data to infer atmospheric conditions
US8386174B2 (en) Method and device for calculating a time window
Nikoleris et al. Effect of trajectory prediction and stochastic runway occupancy times on aircraft delays
US20120253650A1 (en) Method And Device For Checking That An Aircraft Fulfils Successive Time Constraints
Dalmau et al. Fast sensitivity-based optimal trajectory updates for descent operations subject to time constraints
WO2019192615A1 (zh) 基于大数据的空中交通流量管理方法及系统
JP6456442B2 (ja) 移動計画管理システム、プロファイル管理システム、プロファイル管理プログラム及び移動計画管理プログラム
Huang et al. Probabilistic estimation of state-dependent hybrid mode transitions for aircraft arrival time prediction
US11486710B2 (en) System and method for determining a three-dimensional travel path for a vehicle
US20240169110A1 (en) Calculating, storing, and retrieving high-accuracy aircraft performance data
Kasai Car-following model with multiple predicting and controlling modules based on assumptions of anticipation behavior

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19897836

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19897836

Country of ref document: EP

Kind code of ref document: A1