WO2018055825A1 - 性能推定装置及び性能推定方法 - Google Patents

性能推定装置及び性能推定方法 Download PDF

Info

Publication number
WO2018055825A1
WO2018055825A1 PCT/JP2017/015591 JP2017015591W WO2018055825A1 WO 2018055825 A1 WO2018055825 A1 WO 2018055825A1 JP 2017015591 W JP2017015591 W JP 2017015591W WO 2018055825 A1 WO2018055825 A1 WO 2018055825A1
Authority
WO
WIPO (PCT)
Prior art keywords
gain
value
operating point
estimation
sensor
Prior art date
Application number
PCT/JP2017/015591
Other languages
English (en)
French (fr)
Inventor
中村　恵子
大紀 垣内
麻衣 江本
祐司 小熊
Original Assignee
株式会社Ihi
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 株式会社Ihi filed Critical 株式会社Ihi
Publication of WO2018055825A1 publication Critical patent/WO2018055825A1/ja

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D25/00—Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
    • F02C7/00—Features, components parts, details or accessories, not provided for in, or of interest apart form groups F02C1/00 - F02C6/00; Air intakes for jet-propulsion plants
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
    • F02C9/00—Controlling gas-turbine plants; Controlling fuel supply in air- breathing jet-propulsion plants
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00—Testing or monitoring of control systems or parts thereof
    • G05B23/02—Electric testing or monitoring

Definitions

  • FIG. 1 is a diagram illustrating a schematic configuration of a performance estimation system including a performance estimation device according to an embodiment.
  • a performance estimation system 1 shown in FIG. 1 is a system that estimates the state of a target device based on sensor measurement values measured by a plurality of sensors provided in the target device that is an estimation target device.
  • a gas turbine engine will be described as an example of the target device.
  • (k) is added to the code representing the variable, but the description is not limited to a specific time, and is established at an arbitrary time.
  • the unit time of the performance estimation system 1 is “1”, and the larger the value in parentheses, the longer the time has elapsed.
  • the internal state quantity x e is an amount of the state inside the gas turbine engine 10 including a state quantity that cannot be measured by the sensor.
  • the internal state quantity x e is a ne- dimensional vector and includes the rotational speed of each element, the internal energy of each element, the air flow rate of each element, the enthalpy of each element, and the like.
  • the performance parameter value q is an engine performance parameter value that cannot be measured by the sensor.
  • the performance parameter value q is a p-dimensional vector and includes the fan flow rate, fan efficiency change, the flow rate of each compressor, the efficiency change of each compressor, and the like.
  • FIG. 4 is a diagram showing an example of the relationship between the operating point and the gain.
  • the gas turbine engine 10 is provided in an aircraft, and description is made using altitude and speed as engine operating conditions.
  • the gain storage unit 32 stores gains K1 to K4 and a map in which the ranges R1 to R4 are associated with the operating points included in each range.
  • the gain determining unit 33 gradually switches the set gain K (k) from the gain K1 to the gain K2.
  • the gain determination unit 33 changes the setting gain K (k) continuously (stepwise) from the gain K1 to the gain K2 by multiplying the switching time Ts.
  • the switching time Ts is set so that the amount of change of the setting gain K (k) per unit time does not diverge the estimation in the performance estimation device 30. That is, the gain determination unit 33 performs gain interpolation when switching the set gain K (k), and otherwise sets a fixed value as the set gain K (k).
  • the gain determination unit 33 outputs the determined set gain K (k) to the correction unit 34. Details of the setting gain K (k) determination process will be described later.
  • the estimation unit 35 uses the nonlinear simulation model of the gas turbine engine 10 to calculate the state estimated value x h (k) and the sensor estimated value y h (k) based on the control input value u (k) and the correction value. calculate.
  • the estimation unit 35 outputs the calculated sensor estimation value y h (k) to the error calculation unit 31, and outputs the calculated state estimation value x h (k) to the outside of the performance estimation device 30 (control device 2).
  • the estimated value updating unit 52 calculates the sensor estimated value y h (k) based on the prior estimated value x h ⁇ (k) output from the time updating unit 51. Specifically, the estimated value update unit 52 calculates the sensor estimated value y h (k) by using the nonlinear function g and Equation (2). The estimated value updating unit 52 outputs the calculated sensor estimated value y h (k) to the error calculating unit 31.
  • the estimation unit 35 calculates the state estimated value x h (k) based on the control input value u (k) and the correction value using the nonlinear simulation model of the gas turbine engine 10 (step S08). Specifically, the estimated value updating unit 53 of the estimating unit 35 adds the pre-estimated value x h ⁇ (k) and the correction value as shown in the equation (3), and the addition result is used as the state estimated value. Calculate as x h (k). Then, the estimated value updating unit 53 outputs the calculated state estimated value x h (k) to the outside and the delay unit 54 (step S09). As described above, a series of processes of the performance estimation method at time k ends.

Abstract

本開示の一側面に係る性能推定装置は、対象装置を制御するための制御入力値及び対象装置に設けられたセンサによって計測されたセンサ計測値に基づいて、センサによる計測が不可能な対象装置の状態量の推定値である状態推定値を出力する性能推定装置であって、対象装置の非線形シミュレーションモデルを用いて、状態推定値、及びセンサ推定値を算出する推定部と、センサ推定値とセンサ計測値との誤差を算出する誤差算出部と、第１作動点に対応する第１ゲイン、及び第２作動点に対応する第２ゲインを記憶するゲイン記憶部と、第１ゲイン及び第２ゲインに基づいて、対象装置が動作している作動点である設定作動点における設定ゲインを決定するゲイン決定部と、誤差と設定ゲインとに基づいて、補正値を生成する補正部と、を備え、推定部は、制御入力値及び補正値に基づいて、状態推定値及びセンサ推定値を算出し、状態推定値を出力する。

Description

性能推定装置及び性能推定方法
　本開示は、性能推定装置及び性能推定方法に関する。
　ガスタービンエンジン等では、センサによる計測が不可能な内部状態量及び性能パラメータが存在する。このようなセンサによる計測が不可能な内部状態量及び性能パラメータの値を推定する手法として、カルマンフィルタを用いた手法がある。例えば、特許文献１には、エンジンの非線形シミュレーションモデルを含むカルマンフィルタを用いて、センサ計測値から観測不能変数推定値を演算するガスタービンエンジンの性能推定システムが記載されている。
特開２００９－６８３５９号公報
　特許文献１に記載の性能推定システムでは、カルマンフィルタゲインが一定とされている。このカルマンフィルタゲインは、ガスタービンエンジンのある作動点において求められる。非線形シミュレーションモデルを用いた場合、非線形性が強い対象装置では、１つのゲインで全作動点に対応しようとすると以下のような問題がある。つまり、カルマンフィルタでは、推定の安定性と速応性とがトレードオフの関係にあるので、安定性を重視してカルマンフィルタゲインを設定すると応答が遅くなり、速応性を重視してカルマンフィルタゲインを設定すると、安定性が低下する。
　本開示は、センサによる計測が不可能な状態量の推定における安定性と速応性とを両立可能な性能推定装置及び性能推定方法を説明する。
　本開示の一側面に係る性能推定装置は、推定対象の装置である対象装置を制御するための制御入力値及び対象装置に設けられたセンサによって計測されたセンサ計測値に基づいて、センサによる計測が不可能な対象装置の状態量の推定値である状態推定値を出力する装置である。この性能推定装置は、対象装置の非線形シミュレーションモデルを用いて、状態推定値、及びセンサ計測値の推定値であるセンサ推定値を算出する推定部と、センサ推定値とセンサ計測値との誤差を算出する誤差算出部と、対象装置の１つの作動点である第１作動点に対応する第１ゲイン、及び対象装置の別の作動点である第２作動点に対応する第２ゲインを記憶するゲイン記憶部と、第１ゲイン及び第２ゲインに基づいて、対象装置が動作している作動点である設定作動点における設定ゲインを決定するゲイン決定部と、誤差と設定ゲインとに基づいて、補正値を生成する補正部と、を備える。推定部は、制御入力値及び補正値に基づいて、状態推定値及びセンサ推定値を算出し、状態推定値を出力する。
　本開示によれば、センサによる計測が不可能な状態量の推定における安定性と速応性とを両立できる。
図１は、一実施形態に係る性能推定装置を含む性能推定システムの概略構成を示す図である。 図２は、図１のガスタービンエンジンの各要素を模式的に示す図である。 図３は、図１の性能推定装置のハードウェア構成図である。 図４は、作動点とゲインとの関係の一例を示す図である。 図５は、図１の推定部の構成例を示す図である。 図６の（ａ）～（ｄ）は、設定作動点及び設定ゲインの時間変化の一例を示す図である。 図７は、図１の性能推定装置が行う性能推定方法の一連の処理を示すフローチャートである。
［１］実施形態の概要
　本開示の一側面に係る性能推定装置は、推定対象の装置である対象装置を制御するための制御入力値及び対象装置に設けられたセンサによって計測されたセンサ計測値に基づいて、センサによる計測が不可能な対象装置の状態量の推定値である状態推定値を出力する装置である。この性能推定装置は、対象装置の非線形シミュレーションモデルを用いて、状態推定値、及びセンサ計測値の推定値であるセンサ推定値を算出する推定部と、センサ推定値とセンサ計測値との誤差を算出する誤差算出部と、対象装置の１つの作動点である第１作動点に対応する第１ゲイン、及び対象装置の別の作動点である第２作動点に対応する第２ゲインを記憶するゲイン記憶部と、第１ゲイン及び第２ゲインに基づいて、対象装置が動作している作動点である設定作動点における設定ゲインを決定するゲイン決定部と、誤差と設定ゲインとに基づいて、補正値を生成する補正部と、を備える。推定部は、制御入力値及び補正値に基づいて、状態推定値及びセンサ推定値を算出し、状態推定値を出力する。
　この性能推定装置では、対象装置の非線形シミュレーションモデルを用いて、状態推定値及びセンサ推定値が算出される。この状態推定値及びセンサ推定値は、センサ推定値及びセンサ計測値の誤差と設定作動点における設定ゲインとに基づいて生成された補正値、並びに制御入力値に基づいて、算出される。設定ゲインは、第１ゲイン及び第２ゲインに基づいて決定される。このように、２つの作動点におけるゲインに基づいて、設定ゲインが決定されるので、全作動点の範囲を１つの固定ゲインで対応する場合と比較して、１つのゲインで対応する作動点の範囲が狭くなる。このため、全作動点の範囲を１つの固定ゲインで対応する場合と比較して、各ゲインが対応する範囲において、推定の安定性を向上させつつ、速応性を向上させるように第１ゲイン及び第２ゲインを設定することが可能となる。これらの第１ゲイン及び第２ゲインに基づいて、設定ゲインが決定されることから、全作動点の範囲を１つの固定ゲインで対応する場合と比較して、センサによる計測が不可能な状態量の推定における安定性と速応性とを両立することが可能となる。
　ゲイン決定部は、第１作動点を含む作動点の範囲である第１範囲に設定作動点が含まれる場合には第１ゲインを設定ゲインとしてもよく、第２作動点を含む作動点の範囲である第２範囲に設定作動点が含まれる場合には第２ゲインを設定ゲインとしてもよい。この場合、設定作動点が含まれる範囲に応じて、ゲイン記憶部に記憶されている第１ゲイン又は第２ゲインが設定ゲインとして用いられる。このため、ゲイン記憶部からゲインを取得して設定ゲインとすればよいので、設定ゲインの計算負荷を軽減することが可能となる。
　第１範囲と第２範囲とは互いに隣接していてもよい。ゲイン決定部は、設定作動点が第１範囲から第２範囲に移動した場合には、設定ゲインを第１ゲインから第２ゲインに徐々に切り替えてもよい。この場合、設定ゲインを徐々に切り替えることにより、設定ゲインの切り替えによる補正値の変化を抑えることができ、性能推定装置における推定が発散する可能性を低減することが可能となる。
　本開示の別の側面に係る性能推定方法は、推定対象の装置である対象装置に設けられたセンサによる計測が不可能な状態量の推定値である状態推定値を出力する方法である。この性能推定方法は、対象装置が動作している作動点である設定作動点を受信するステップと、対象装置を制御するための制御入力値を受信するステップと、対象装置に設けられたセンサによって計測されたセンサ計測値を受信するステップと、対象装置の非線形シミュレーションモデルを用いて、制御入力値に基づいて、センサ計測値の推定値であるセンサ推定値を算出するステップと、センサ推定値とセンサ計測値との誤差を算出するステップと、対象装置の１つの作動点である第１作動点に対応する第１ゲイン、及び対象装置の別の作動点である第２作動点に対応する第２ゲインに基づいて、設定作動点における設定ゲインを決定するステップと、誤差と設定ゲインとに基づいて、補正値を生成するステップと、非線形シミュレーションモデルを用いて、制御入力値及び補正値に基づいて、状態推定値を算出するステップと、状態推定値を出力するステップと、を含む。
　この性能推定方法では、対象装置の非線形シミュレーションモデルを用いて、制御入力値に基づいて、センサ推定値が算出される。対象装置の非線形シミュレーションモデルを用いて、センサ推定値及びセンサ計測値の誤差と設定作動点における設定ゲインとに基づいて生成された補正値、並びに制御入力値に基づいて、状態推定値が算出される。設定ゲインは、第１ゲイン及び第２ゲインに基づいて決定される。このように、２つの作動点におけるゲインに基づいて、設定ゲインが決定されるので、全作動点の範囲を１つの固定ゲインで対応する場合と比較して、１つのゲインで対応する作動点の範囲が狭くなる。このため、全作動点の範囲を１つの固定ゲインで対応する場合と比較して、各ゲインが対応する範囲において、推定の安定性を向上させつつ、速応性を向上させるように第１ゲイン及び第２ゲインを設定することが可能となる。これらの第１ゲイン及び第２ゲインに基づいて、設定ゲインが決定されることから、全作動点の範囲を１つの固定ゲインで対応する場合と比較して、センサによる計測が不可能な状態量の推定における安定性と速応性とを両立することが可能となる。
［２］実施形態の例示
　以下、本開示の実施形態について、図面を参照しながら説明する。なお、図面の説明において同一要素には同一符号を付し、重複する説明を省略する。
　図１は、一実施形態に係る性能推定装置を含む性能推定システムの概略構成を示す図である。図１に示される性能推定システム１は、推定対象の装置である対象装置に設けられる複数のセンサによって計測されたセンサ計測値に基づいて、対象装置の状態を推定するシステムである。以下、対象装置の一例として、ガスタービンエンジンを用いて説明する。なお、時間変化する変数の時刻ｋにおける値を表す場合に、当該変数を表す符号に（ｋ）を付して説明を行う場合があるが、特定の時刻に限られず、任意の時刻において成り立つ。また、性能推定システム１の単位時間を「１」としており、括弧内の値が大きいほど時間が経過していることを示している。
　性能推定システム１は、制御装置２と、ガスタービンエンジン１０と、を備えている。ガスタービンエンジン１０は、例えば、車両、船舶、及び航空機等のエンジンとして用いられる。
　図２は、ガスタービンエンジン１０の各要素を模式的に示す図である。図２に示されるように、ガスタービンエンジン１０は、ファン（ＦＡＮ）１１と、低圧圧縮機（ＬＰＣ）１２と、高圧圧縮機（ＨＰＣ）１３と、燃焼器（ＣＯＭＢ）１４と、高圧タービン（ＨＰＴ）１５と、低圧タービン（ＬＰＴ）１６と、ロータ１７と、ロータ１８と、を備えている。
　ファン１１は、ガスタービンエンジン１０の外部から空気を吸い込んで、吸い込んだ空気の一部を低圧圧縮機１２に供給する。低圧圧縮機１２は、ファン１１から供給された空気を圧縮することによって圧縮空気を生成し、圧縮空気を高圧圧縮機１３に供給する。高圧圧縮機１３は、低圧圧縮機１２から供給された圧縮空気をさらに圧縮して圧力をさらに高くすることによって高圧空気を生成し、高圧空気を燃焼器１４に供給する。燃焼器１４は、高圧圧縮機１３から供給された高圧空気に燃料を混ぜ合わせて燃焼する。燃焼器１４は、燃焼して得られた高温高圧の燃焼ガスを高圧タービン１５に供給する。
　高圧タービン１５は、燃焼器１４から供給された高温高圧の燃焼ガスによって回転し、ロータ１８を回転させる。低圧タービン１６は、高圧タービン１５を通過した燃焼ガスによって回転し、ロータ１７を回転させる。ロータ１７は、ファン１１、低圧圧縮機１２及び低圧タービン１６を一体的に回転可能に接続している。ロータ１８は、高圧圧縮機１３及び高圧タービン１５を一体的に回転可能に接続している。
　図１に戻って、ガスタービンエンジン１０には、複数のセンサ（不図示）が設けられている。ガスタービンエンジン１０は、後述するエンジン制御器２０によって制御される。ガスタービンエンジン１０は、ある時刻ｋにおいて、エンジン制御器２０から出力される制御入力値ｕ（ｋ）を入力し、ガスタービンエンジン１０に設けられているセンサによって検知されるセンサ計測値ｙ（ｋ）を出力する。制御入力値ｕ（ｋ）は、ガスタービンエンジン１０を制御するための値であり、エンジン動作条件に基づいて出力される。ガスタービンエンジン１０が航空機用のエンジンである場合、エンジン動作条件は、飛行条件ともいい、高度、機速及びパワーレバー角度（Power　Lever　Angle：ＰＬＡ）等である。各エンジン動作条件の値の組み合わせを作動点という。エンジン動作条件（作動点）は、飛行時において機体から送信される。制御入力値ｕ（ｋ）は、燃料流量及び抽気バルブの開度等を含む。センサ計測値ｙ（ｋ）は、ｒ次元であり、ガスタービンエンジン１０のロータ１７，１８の回転数、各部における温度及び圧力を含む。
　制御装置２は、エンジン制御器２０と、性能推定装置３０と、を備えている。エンジン制御器２０は、ガスタービンエンジン１０を制御する装置である。エンジン制御器２０は、エンジン動作条件（作動点）に基づいて、制御入力値ｕ（ｋ）を決定する。エンジン制御器２０は、制御入力値ｕ（ｋ）をガスタービンエンジン１０及び性能推定装置３０に出力する。
　性能推定装置３０は、ガスタービンエンジン１０の非線形シミュレーションモデルを用いたカルマンフィルタによって、ガスタービンエンジン１０の状態量ｘ（ｋ）を推定する。非線形シミュレーションモデルは、ガスタービンエンジン１０ごとに予め定められた非線形のシミュレーションモデルであり、非線形関数ｆ及び非線形関数ｇを用いて、後述の式（１）～式（２）で表される。非線形関数ｆは、対象装置（ガスタービンエンジン１０）ごとに予め定められており、状態量ｘ（ｋ）及び制御入力値ｕ（ｋ）の関数である。非線形関数ｇは、対象装置ごとに予め定められており、状態量ｘ（ｋ）の関数である。
　具体的には、性能推定装置３０は、制御入力値ｕ（ｋ）及びセンサ計測値ｙ（ｋ）に基づいて、ガスタービンエンジン１０の状態推定値ｘｈ（ｋ）を出力する。状態推定値ｘｈ（ｋ）は、ガスタービンエンジン１０の状態量ｘ（ｋ）の推定値である。状態量ｘ（ｋ）は、ガスタービンエンジン１０の状態量のうち、センサによる計測が不可能な状態量を含む状態量であり、内部状態量ｘｅ（ｋ）と性能パラメータ値ｑ（ｋ）とを合わせたｎ（ｎ＝ｎｅ＋ｐ）次元のベクトルである。状態推定値ｘｈ（ｋ）は、内部状態量ｘｅ（ｋ）の推定値と性能パラメータ値ｑ（ｋ）の推定値とを合わせたｎ（ｎ＝ｎｅ＋ｐ）次元のベクトルである。
　内部状態量ｘｅは、センサ計測が不可能な状態量を含むガスタービンエンジン１０内部の状態の量である。内部状態量ｘｅは、ｎｅ次元のベクトルであり、各要素の回転数、各要素の内部エネルギー、各要素の空気流量、各要素のエンタルピー等を含む。性能パラメータ値ｑは、センサ計測が不可能なエンジン性能パラメータの値である。性能パラメータ値ｑは、ｐ次元のベクトルであり、ファンの流量、ファンの効率変化、各圧縮機の流量、各圧縮機の効率変化等を含む。
　図３は、性能推定装置３０のハードウェア構成図である。図３に示されるように、性能推定装置３０は、物理的には、１又は複数のプロセッサ３０１、主記憶装置であるＲＡＭ（Random　Access　Memory）及びＲＯＭ（Read　Only　Memory)等の記憶装置３０２、ハードディスク装置等の補助記憶装置３０３、キーボード等の入力装置３０４、ディスプレイ等の出力装置３０５、並びに、データ送受信デバイスである通信装置３０６等のハードウェアを備えるコンピュータとして構成され得る。性能推定装置３０の図１に示される各機能は、記憶装置３０２等のハードウェアに１又は複数の所定のコンピュータプログラムを読み込ませることにより、１又は複数のプロセッサ３０１の制御のもとで各ハードウェアを動作させるとともに、記憶装置３０２及び補助記憶装置３０３におけるデータの読み出し及び書き込みを行うことで実現される。なお、上記の説明は性能推定装置３０のハードウェア構成として説明したが、制御装置２がプロセッサ３０１、記憶装置３０２、補助記憶装置３０３、入力装置３０４、出力装置３０５、及び通信装置３０６を備えるコンピュータとして構成されてもよい。
　図１に戻って、性能推定装置３０の詳細について説明する。性能推定装置３０は、誤差算出部３１と、ゲイン記憶部３２と、ゲイン決定部３３と、補正部３４と、推定部３５と、を備えている。
　誤差算出部３１は、ガスタービンエンジン１０から出力されたセンサ計測値ｙ（ｋ）と推定部３５によって算出されたセンサ推定値ｙｈ（ｋ）との誤差ｅ（ｋ）を算出する。この例では、誤差算出部３１は、センサ計測値ｙ（ｋ）からセンサ推定値ｙｈ（ｋ）を減算し、その減算結果を誤差ｅ（ｋ）として補正部３４に出力する。なお、センサ推定値ｙｈ（ｋ）は、センサ計測値ｙ（ｋ）の推定値である。
　ゲイン記憶部３２は、複数のゲインを記憶している。複数のゲインのそれぞれは、ガスタービンエンジン１０の特徴的な作動点を設計点として設定される。つまり、各ゲインは、対応する作動点において所定の計算によって算出される。各ゲインには、対応する作動点を含む所定の範囲が設定される。ゲイン記憶部３２は、作動点と範囲とを対応付けたマップ（テーブル）を記憶している。
　図４は、作動点とゲインとの関係の一例を示す図である。図４に示される例では、ガスタービンエンジン１０が航空機に設けられており、エンジン動作条件として、高度及び機速を用いて説明している。この例では、ゲイン記憶部３２は、ゲインＫ１～Ｋ４と、範囲Ｒ１～Ｒ４と各範囲に含まれる作動点とを対応付けたマップと、を記憶している。
　図４に示されるように、エンジン動作条件として作動点が取り得る作動範囲Ｒが設定されており、この作動範囲Ｒが４つの範囲Ｒ１～Ｒ４に分割されている。範囲Ｒ１（第１範囲）は、地上付近の作動点の範囲である。範囲Ｒ２（第２範囲）は、離陸（着陸）付近の作動点の範囲である。範囲Ｒ３は、巡航付近の作動点の範囲である。範囲Ｒ４は、高空低速付近の作動点の範囲である。範囲Ｒ１と範囲Ｒ２とは互いに隣接しており、範囲Ｒ２と範囲Ｒ３とは互いに隣接しており、範囲Ｒ２と範囲Ｒ４とは互いに隣接しており、範囲Ｒ３と範囲Ｒ４とは互いに隣接している。
　ゲインＫ１（第１ゲイン）は、範囲Ｒ１に対応するゲインである。ゲインＫ１は、範囲Ｒ１に含まれる静止状態の作動点Ｐ１（第１作動点）を設計点として設定される。ゲインＫ２（第２ゲイン）は、範囲Ｒ２に対応するゲインである。ゲインＫ２は、範囲Ｒ２に含まれる作動点Ｐ２（第２作動点）を設計点として設定される。作動点Ｐ２は、機体が離陸して巡航点に向かって上昇している状態又は機体が巡航点から着陸するために下降している状態の作動点である。ゲインＫ３は、範囲Ｒ３に対応するゲインである。ゲインＫ３は、範囲Ｒ３に含まれる巡航状態の作動点Ｐ３を設計点として設定される。ゲインＫ４は、範囲Ｒ４に対応するゲインである。ゲインＫ４は、範囲Ｒ４に含まれる高空低速の作動点Ｐ４を設計点として設定される。つまり、ゲインＫ１～Ｋ４は、それぞれ別の作動点を設計点として設計され、それぞれ別の範囲に割り当てられる。
　ゲイン決定部３３は、ゲイン記憶部３２に記憶されている複数のゲインに基づいて、設定作動点における設定ゲインＫ（ｋ）を決定する。設定作動点は、ガスタービンエンジン１０が動作している作動点、つまり現在の作動点である。ゲイン決定部３３は、ゲイン記憶部３２に記憶されているマップを参照し、設定作動点がいずれの範囲に含まれるかを判定し、設定作動点が含まれる範囲に基づいて設定ゲインＫ(ｋ)を決定する。
　図４の例では、ゲイン決定部３３は、例えば、作動点Ｐ１を含む範囲Ｒ１に設定作動点が含まれる場合には、ゲインＫ１を設定ゲインＫ(ｋ)とする。ゲイン決定部３３は、例えば、作動点Ｐ２を含む範囲Ｒ２に設定作動点が含まれる場合には、ゲインＫ２を設定ゲインＫ(ｋ)とする。ゲイン決定部３３は、例えば、作動点Ｐ３を含む範囲Ｒ３に設定作動点が含まれる場合には、ゲインＫ３を設定ゲインＫ(ｋ)とする。ゲイン決定部３３は、例えば、作動点Ｐ４を含む範囲Ｒ４に設定作動点が含まれる場合には、ゲインＫ４を設定ゲインＫ(ｋ)とする。
　ゲイン決定部３３は、ガスタービンエンジン１０のエンジン動作条件が変更されることによって、設定作動点がある範囲から異なる範囲に移動した場合（隣り合う範囲間の境界を越えた場合）には、設定ゲインＫ(ｋ)を徐々に切り替える。ゲイン決定部３３は、例えば、１時刻前に設定作動点が含まれていた範囲を記憶しておき、時刻ｋにおいて設定作動点が含まれる範囲と、時刻ｋ－１において設定作動点が含まれていた範囲とが異なる場合に、設定作動点が２つの範囲間を移動したと判定する。
　ゲイン決定部３３は、例えば、設定作動点が範囲Ｒ１から範囲Ｒ２に移動した場合には、設定ゲインＫ(ｋ)をゲインＫ１からゲインＫ２に徐々に切り替える。このとき、ゲイン決定部３３は、設定ゲインＫ(ｋ)をゲインＫ１からゲインＫ２まで切替時間Ｔｓを掛けて連続的（段階的）に変更する。設定ゲインＫ(ｋ)の切替処理において、単位時間ごとの設定ゲインＫ(ｋ)の変化量が、性能推定装置３０における推定を発散させない程度となるように、切替時間Ｔｓは設定される。つまり、ゲイン決定部３３は、設定ゲインＫ(ｋ)を切り替える際にゲイン補間を行い、それ以外では固定値を設定ゲインＫ(ｋ)とする。ゲイン決定部３３は、決定した設定ゲインＫ(ｋ)を補正部３４に出力する。設定ゲインＫ(ｋ)の決定処理の詳細は後述する。
　補正部３４は、誤差算出部３１によって算出された誤差ｅ（ｋ）とゲイン決定部３３によって決定された設定ゲインＫ（ｋ）とに基づいて、推定部３５に事前推定値ｘｈ －（ｋ）を補正させるための補正値を生成する。事前推定値ｘｈ －（ｋ）は、状態推定値ｘｈ（ｋ）の推定誤差を補正する前の推定値である。具体的には、補正部３４は、誤差ｅ（ｋ）が小さくなるように、補正値を生成する。補正部３４は、誤差ｅ（ｋ）と設定ゲインＫ（ｋ）とを乗算し、その乗算結果を補正値として推定部３５に出力する。
　推定部３５は、ガスタービンエンジン１０の非線形シミュレーションモデルを用いて、制御入力値ｕ（ｋ）及び補正値に基づいて、状態推定値ｘｈ（ｋ）、及びセンサ推定値ｙｈ（ｋ）を算出する。推定部３５は、算出したセンサ推定値ｙｈ（ｋ）を誤差算出部３１に出力し、算出した状態推定値ｘｈ（ｋ）を性能推定装置３０（制御装置２）の外部に出力する。
　図５は、推定部３５の構成例を示す図である。図５に示されるように、推定部３５は、時間更新部５１と、推定値更新部５２と、推定値更新部５３と、遅延部５４と、を備えている。
　時間更新部５１は、制御入力値ｕ（ｋ）と、遅延部５４から出力された１時刻前の状態推定値ｘｈ（ｋ－１）と、に基づいて、事前推定値ｘｈ －（ｋ）を算出する。具体的には、時間更新部５１は、非線形関数ｆを用い、式（１）によって事前推定値ｘｈ －（ｋ）を算出する。時間更新部５１は、算出した事前推定値ｘｈ －（ｋ）を推定値更新部５２及び推定値更新部５３に出力する。
Figure JPOXMLDOC01-appb-M000001
　推定値更新部５２は、時間更新部５１から出力された事前推定値ｘｈ －（ｋ）に基づいて、センサ推定値ｙｈ（ｋ）を算出する。具体的には、推定値更新部５２は、非線形関数ｇを用い、式（２）によってセンサ推定値ｙｈ（ｋ）を算出する。推定値更新部５２は、算出したセンサ推定値ｙｈ（ｋ）を誤差算出部３１に出力する。
Figure JPOXMLDOC01-appb-M000002
　推定値更新部５３は、時間更新部５１から出力された事前推定値ｘｈ －（ｋ）と、補正部３４から出力された補正値とに基づいて、状態推定値ｘｈ（ｋ）を算出する。具体的には、推定値更新部５３は加算器であり、式（３）に示されるように、事前推定値ｘｈ －（ｋ）と補正値とを加算し、その加算結果を状態推定値ｘｈ（ｋ）として算出する。推定値更新部５３は、算出した状態推定値ｘｈ（ｋ）を遅延部５４及び外部に出力する。
Figure JPOXMLDOC01-appb-M000003
　遅延部５４は、推定値更新部５３から出力された状態推定値ｘｈ（ｋ）を単位時間だけ遅延させ、次の時刻における推定に対し、状態推定値ｘｈ（ｋ－１）として時間更新部５１に出力する。
　ここで、図４及び図６を用いて、ゲイン決定処理を詳細に説明する。図６の（ａ）～（ｄ）は、設定作動点及び設定ゲインの時間変化の一例を示す図であり、図６の（ａ）は高度の時間変化の一例を示す図、図６の（ｂ）は機速の時間変化の一例を示す図、図６の（ｃ）はＰＬＡの時間変化の一例を示す図、図６の（ｄ）は設定ゲインの時間変化の一例を示す図である。各図の横軸は時間を示し、縦軸は作動点（高度、機速、及びＰＬＡ）及び設定ゲインの値を示している。図４の運用パターンＦに示されるように、航空機が離陸して巡航点まで上昇する場合を例にして説明する。
　まず、航空機の停止状態では、高度、機速及びＰＬＡはいずれも０である（作動点Ｐ１）。このため、設定作動点は範囲Ｒ１に含まれるので、ゲイン決定部３３は、ゲインＫ１を設定ゲインＫ(ｋ)として決定する。時刻０において航空機が動作を開始すると、機速が次第に大きくなり、航空機が離陸すると、ＰＬＡが増加するとともに高度が高くなる。そして、航空機が上昇を続け、高度、機速及びＰＬＡが増加し続ける。そして、時刻ｋ１において設定作動点が範囲Ｒ１と範囲Ｒ２との境界（ゲイン切替点Ｐｓ１）に達する。したがって、時刻０～ｋ１では、設定作動点は範囲Ｒ１内に留まっているので、ゲイン決定部３３は、設定ゲインＫ(ｋ)をゲインＫ１で維持する。
　時刻ｋ１において設定作動点がゲイン切替点Ｐｓ１に達すると、時刻ｋ１～ｋ２では、ゲイン決定部３３は、設定ゲインＫ(ｋ)をゲインＫ１からゲインＫ２まで切替時間Ｔｓを掛けて連続的（段階的）に変更する。時刻ｋ２は、時刻ｋ１に切替時間Ｔｓを加えた時刻である。ゲイン決定部３３は、例えば、式（４）を用いて、設定ゲインＫ(ｋ)を算出する。性能推定装置３０は、式（４）を用いて単位時間ごとに設定ゲインＫ(ｋ)を更新するので、単位時間では、設定ゲインＫ(ｋ)の変化量は（Ｋ２－Ｋ１）／Ｔｓとなる。設定ゲインＫ(ｋ)の変化量によって性能推定装置３０における推定が発散しないように、切替時間Ｔｓは設定されている。
Figure JPOXMLDOC01-appb-M000004
　そして、ＰＬＡが調整されながら、航空機がさらに上昇を続け、高度及び機速は増加し続ける。そして、時刻ｋ３において設定作動点が範囲Ｒ２と範囲Ｒ３との境界（ゲイン切替点Ｐｓ２）に達する。したがって、時刻ｋ１～ｋ３では、設定作動点は範囲Ｒ２内に留まっており、時刻ｋ２において設定ゲインＫ(ｋ)をゲインＫ１からゲインＫ２に切り替える処理が完了するので、時刻ｋ２～ｋ３では、ゲイン決定部３３は、設定ゲインＫ(ｋ)をゲインＫ２で維持する。
　時刻ｋ３において設定作動点がゲイン切替点Ｐｓ２に達すると、時刻ｋ３～ｋ４では、ゲイン決定部３３は、設定ゲインＫ(ｋ)をゲインＫ２からゲインＫ３まで切替時間Ｔｓを掛けて連続的（段階的）に変更する。時刻ｋ４は、時刻ｋ３に切替時間Ｔｓを加えた時刻である。ゲイン決定部３３は、例えば、式（５）を用いて、設定ゲインＫ(ｋ)を算出する。性能推定装置３０は、式（５）を用いて単位時間ごとに設定ゲインＫ(ｋ)を更新するので、単位時間では、設定ゲインＫ(ｋ)の変化量は（Ｋ３－Ｋ２）／Ｔｓとなる。設定ゲインＫ(ｋ)の変化量によって性能推定装置３０における推定が発散しないように、切替時間Ｔｓは設定されている。
Figure JPOXMLDOC01-appb-M000005
　そして、ＰＬＡが調整されながら、航空機がさらに上昇を続け、高度及び機速は増加し続ける。そして、設定作動点が巡航点（作動点Ｐ３）に達する。巡航点では、高度、機速及びＰＬＡは略一定である。したがって、時刻ｋ３以降では、設定作動点は範囲Ｒ３内に留まっており、時刻ｋ４において設定ゲインＫ(ｋ)をゲインＫ２からゲインＫ３に切り替える処理が完了するので、時刻ｋ４以降では、ゲイン決定部３３は、設定ゲインＫ(ｋ)をゲインＫ３で維持する。
　このように、ゲイン決定部３３は、設定作動点が、作動範囲Ｒにおいてゲインに対して設定されている範囲を超えた場合には、設定ゲインＫ(ｋ)を設定されているゲインから新たなゲインまで切替時間Ｔｓで徐々に変化させる。それ以外では、ゲイン決定部３３は、設定作動点が属する範囲に対応するゲインを設定ゲインＫ(ｋ)として決定する。また、ゲイン決定部３３は、設定ゲインＫ(ｋ)の切り替えを開始した時刻、つまり、設定作動点が範囲の境界を越えた時刻（時刻ｋ１，ｋ３）を記憶しておき、設定ゲインＫ(ｋ)の切り替えを開始した時刻から切替時間Ｔｓを経過していない場合には、切替処理を継続することとしてもよい。なお、航空機が巡航点から下降して地上に着陸する場合も同様に、設定ゲインＫ(ｋ)が決定される。
　次に、図７を参照して、性能推定装置３０が行う性能推定方法の一連の処理を説明する。図７は、性能推定装置３０が行う性能推定方法の一連の処理を示すフローチャートである。図７に示される処理は、単位時間ごとに繰り返し行われる。ここでは、ある時刻ｋにおける処理を説明する。この性能推定方法は、状態推定値ｘｈ（ｋ）を出力するための方法である。
　まず、機体から設定作動点が出力され、制御装置２がその設定作動点を受信する。受信した設定作動点は、エンジン制御器２０及び性能推定装置３０それぞれに入力される。性能推定装置３０では、ゲイン決定部３３が設定作動点を受信する（ステップＳ０１）。また、エンジン制御器２０は、設定作動点を受信すると、設定作動点に基づいて制御入力値ｕ（ｋ）を決定する。そして、エンジン制御器２０は、制御入力値ｕ（ｋ）をガスタービンエンジン１０及び性能推定装置３０に出力する。
　そして、性能推定装置３０では、推定部３５が、制御入力値ｕ（ｋ）を受信する（ステップＳ０２）。また、ガスタービンエンジン１０は、エンジン制御器２０から出力された制御入力値ｕ（ｋ）を入力し、ガスタービンエンジン１０に設けられているセンサによって計測されたセンサ計測値ｙ（ｋ）を性能推定装置３０に出力する。そして、性能推定装置３０の誤差算出部３１は、センサ計測値ｙ（ｋ）を受信する（ステップＳ０３）。
　続いて、推定部３５は、ガスタービンエンジン１０の非線形シミュレーションモデルを用いて、制御入力値ｕ（ｋ）に基づいて、センサ推定値ｙｈ（ｋ）を算出する（ステップＳ０４）。ステップＳ０４の処理を具体的に説明する。まず、推定部３５の時間更新部５１は、制御入力値ｕ（ｋ）と、遅延部５４から出力された状態推定値ｘｈ（ｋ－１）と、を用いて、式（１）によって事前推定値ｘｈ －（ｋ）を算出する。そして、時間更新部５１は、算出した事前推定値ｘｈ －（ｋ）を推定値更新部５２及び推定値更新部５３に出力する。そして、推定値更新部５２は、時間更新部５１から出力された事前推定値ｘｈ －（ｋ）を用いて、式（２）によってセンサ推定値ｙｈ（ｋ）を算出する。そして、推定値更新部５２は、算出したセンサ推定値ｙｈ（ｋ）を誤差算出部３１に出力する。
　続いて、誤差算出部３１は、センサ計測値ｙ（ｋ）とセンサ推定値ｙｈ（ｋ）との誤差ｅ（ｋ）を算出する（ステップＳ０５）。具体的には、誤差算出部３１は、センサ計測値ｙ（ｋ）からセンサ推定値ｙｈ（ｋ）を減算し、その減算結果を誤差ｅ（ｋ）として補正部３４に出力する。
　続いて、ゲイン決定部３３は、設定作動点に基づいて、設定ゲインＫ(ｋ)を決定する（ステップＳ０６）。ステップＳ０６の処理を具体的に説明する。まず、ゲイン決定部３３は、ゲイン記憶部３２に記憶されているマップを参照し、設定作動点がいずれの範囲（図４の例では、範囲Ｒ１～Ｒ４）に含まれるかを判定する。このとき、ゲイン決定部３３は、時刻ｋ－１から時刻ｋの間に、設定作動点が２つの範囲間を移動したか否か（隣り合う２つの範囲の境界を越えたか否か）を判定する。ゲイン決定部３３は、例えば、１時刻前に設定作動点が含まれていた範囲を記憶しておき、時刻ｋにおいて設定作動点が含まれる範囲と、時刻ｋ－１において設定作動点が含まれていた範囲とが異なる場合に、設定作動点が境界を越えて２つの範囲間を移動したと判定する。
　そして、ゲイン決定部３３は、例えば、設定作動点が範囲Ｒ１から範囲Ｒ２に移動した場合には、式（４）によって、設定ゲインＫ(ｋ)をゲインＫ１からゲインＫ２に切替時間Ｔｓを掛けて徐々に切り替える。一方、ゲイン決定部３３は、設定作動点が２つの範囲間を移動していないと判定した場合、前回設定ゲインＫ(ｋ)の切り替えを開始した時刻から切替時間Ｔｓが経過しているか否かを判定する。ゲイン決定部３３は、切替時間Ｔｓが経過していないと判定した場合には、設定ゲインＫ(ｋ)の切替処理を継続する。一方、ゲイン決定部３３は、切替時間Ｔｓが経過していると判定した場合には、設定作動点が含まれる範囲に対応付けられたゲインを設定ゲインＫ(ｋ)として決定する（つまり、設定ゲインＫ(ｋ)を変更しない）。そして、ゲイン決定部３３は、決定した設定ゲインＫ(ｋ)を補正部３４に出力する。
　続いて、補正部３４は、誤差算出部３１によって算出された誤差ｅ（ｋ）とゲイン決定部３３によって決定された設定ゲインＫ(ｋ)とに基づいて、補正値を生成する（ステップＳ０７）。具体的には、補正部３４は、誤差ｅ（ｋ）と設定ゲインＫ(ｋ)とを乗算し、その乗算結果Ｋ(ｋ)×ｅ（ｋ）を補正値として推定部３５に出力する。
　続いて、推定部３５は、ガスタービンエンジン１０の非線形シミュレーションモデルを用いて、制御入力値ｕ（ｋ）及び補正値に基づいて、状態推定値ｘｈ（ｋ）を算出する（ステップＳ０８）。具体的には、推定部３５の推定値更新部５３が、式（３）に示されるように、事前推定値ｘｈ －（ｋ）と補正値とを加算し、その加算結果を状態推定値ｘｈ（ｋ）として算出する。そして、推定値更新部５３は、算出した状態推定値ｘｈ（ｋ）を外部及び遅延部５４に出力する（ステップＳ０９）。以上のようにして、時刻ｋにおける性能推定方法の一連の処理が終了する。
　以上説明したように、性能推定装置３０では、ガスタービンエンジン１０の非線形シミュレーションモデルを用いて、状態推定値ｘｈ（ｋ）及びセンサ推定値ｙｈ（ｋ）が算出される。この状態推定値ｘｈ（ｋ）及びセンサ推定値ｙｈ（ｋ）は、センサ推定値ｙｈ（ｋ）及びセンサ計測値ｙ（ｋ）の誤差ｅ（ｋ）と設定作動点における設定ゲインＫ(ｋ)とに基づいて生成された補正値、並びに制御入力値ｕ（ｋ）に基づいて、算出される。設定ゲインＫ(ｋ)は、予め設定された複数のゲイン（図４の例では、ゲインＫ１～Ｋ４）に基づいて決定される。このように、複数の作動点におけるゲインに基づいて、設定ゲインＫ(ｋ)が決定されるので、全作動点の作動範囲Ｒを１つの固定ゲインで対応する場合と比較して、１つのゲインで対応する作動点の範囲（図４の例では、範囲Ｒ１～Ｒ４）が狭くなる。このため、全作動点の作動範囲Ｒを１つの固定ゲインで対応する場合と比較して、各範囲において、推定の安定性を向上させつつ、速応性を向上させるように各ゲイン（ゲインＫ１～Ｋ４）を設定することが可能となる。これらの複数のゲインに基づいて、設定ゲインＫ(ｋ)が決定されることから、全作動点の作動範囲Ｒを１つの固定ゲインで対応する場合と比較して、センサによる計測が不可能な状態量の推定における安定性と速応性とを両立することが可能となる。
　また、設定作動点が含まれる範囲（範囲Ｒ１～Ｒ４）に応じて、ゲイン記憶部３２に記憶されているゲインＫ１～Ｋ４が設定ゲインＫ(ｋ)として用いられる。このため、ゲイン記憶部３２からゲインＫ１～Ｋ４のいずれかを取得して設定ゲインＫ(ｋ)とすればよいので、設定ゲインＫ(ｋ)の計算負荷を軽減することが可能となる。
　設定ゲインＫ(ｋ)を急に切り替えると、補正値が大きく変化することがある。この場合、センサ推定値ｙｈ（ｋ）が大きく変化するので、誤差ｅ（ｋ）が大きく変化し、これにより補正値がさらに大きく変化する。このような動作が繰り返され、性能推定装置３０の推定が発散するおそれがある。これに対して、設定ゲインＫ(ｋ)を徐々に切り替えることにより、設定ゲインＫ(ｋ)の切り替えによる補正値の変化を抑えることができ、性能推定装置３０における推定が発散する可能性を低減することが可能となる。なお、航空機用のガスタービンエンジン等では、運用パターンがある程度決まっており、設定ゲインＫ(ｋ)の切り替え頻度は高くない。図４に示される例では、１回のフライトで、設定ゲインＫ(ｋ)は４回切り替えられるだけである。また、設定作動点に対応するゲインを逐次計算する計算負荷よりも、設定ゲインＫ(ｋ)の切替処理の計算負荷の方が低い。このため、設定ゲインＫ(ｋ)を徐々に切り替えたとしても、設定ゲインＫ(ｋ)の計算負荷を軽減することが可能となる。
　以上、本開示の実施形態について説明したが、本発明は上記実施形態に限定されない。例えば、上記実施形態では、対象装置の一例としてガスタービンエンジンが用いられているがこれに限られない。対象装置は、制御入力値ｕ（ｋ）に対して非線形な動作を行う装置であって、制御入力値ｕ（ｋ）を入力とし、センサ計測値ｙ（ｋ）及び状態量ｘ（ｋ）を出力とする非線形シミュレーションモデルで表現し得る装置であればよい。
　また、上記実施形態では、性能推定装置３０は、カルマンフィルタを用いて、ガスタービンエンジン１０（対象装置）の状態量ｘ（ｋ）を推定しているが、これに限られない。状態量ｘ（ｋ）を推定するための手法は、対象装置の非線形シミュレーションモデルとゲインとを用いる状態推定手法であればよい。
　また、切替時間Ｔｓは、互いに隣接する領域に対応する２つのゲインの差分の大きさに応じて設定されてもよい。例えば、上述の式（４）の切替時間Ｔｓと式（５）の切替時間Ｔｓとは異なる値に設定されてもよい。
　また、作動範囲Ｒはエンジン動作条件の取り得る範囲に応じて適宜設定され得る。所望の作動点が、各ゲインの設計点として採用され得る。採用された作動点に応じて作動範囲Ｒが分割され得る。つまり、作動範囲Ｒは４つの範囲に限られず、２以上の任意の数の範囲に分割され得る。例えば、作動範囲Ｒが２つの範囲Ｒ１及び範囲Ｒ２に分割される場合、ゲイン記憶部３２は、作動点Ｐ１に対応するゲインＫ１、及び作動点Ｐ２に対応するゲインＫ２を記憶する。ゲイン決定部３３は、ゲインＫ１及びゲインＫ２に基づいて、設定ゲインＫ(ｋ)を決定する。
　また、図７のステップＳ０６において、設定ゲインＫ（ｋ）の切り替え処理を中断しないように、設定ゲインＫ(ｋ)を決定してもよい。例えば、ゲイン決定部３３は、まず前回設定ゲインＫ(ｋ)の切り替えを開始した時刻から切替時間Ｔｓが経過しているか否かを判定し、切替時間Ｔｓが経過していないと判定した場合には、設定ゲインＫ(ｋ)の切替処理を継続してもよい。一方、ゲイン決定部３３は、切替時間Ｔｓが経過していると判定した場合には、ゲイン記憶部３２に記憶されているマップを参照し、設定作動点がいずれの範囲に含まれるかを判定するとともに、前回設定ゲインＫ(ｋ)の切り替えを開始した時刻から時刻ｋの間に、設定作動点が２つの範囲間を移動したか否かを判定してもよい。そして、ゲイン決定部３３は、設定作動点が２つの範囲間を移動した場合には、設定ゲインＫ(ｋ)を移動元の範囲に対応付けられたゲインから移動先の範囲に対応付けられたゲインに切替時間Ｔｓを掛けて徐々に切り替える一方、設定作動点が２つの範囲間を移動していないと判定した場合には、設定作動点が含まれる範囲に対応付けられたゲインを設定ゲインＫ(ｋ)として決定してもよい。
　本開示の性能推定装置及び性能推定方法によれば、センサによる計測が不可能な状態量の推定における安定性と速応性とを両立できる。
１０　ガスタービンエンジン
３０　性能推定装置
３１　誤差算出部
３２　ゲイン記憶部
３３　ゲイン決定部
３４　補正部
３５　推定部
５１　時間更新部
５２　推定値更新部
５３　推定値更新部
５４　遅延部
ｅ　誤差
Ｋ１　ゲイン
Ｋ２　ゲイン
Ｋ３　ゲイン
Ｋ４　ゲイン
Ｋ　設定ゲイン
Ｐ１　作動点
Ｐ２　作動点
Ｐ３　作動点
Ｐ４　作動点
Ｒ　作動範囲
Ｒ１　範囲
Ｒ２　範囲
Ｒ３　範囲
Ｒ４　範囲
Ｔｓ　切替時間
ｕ　制御入力値
ｘ　状態量
ｘｈ　状態推定値
ｘｈ －　事前推定値
ｙ　センサ計測値
ｙｈ　センサ推定値

Claims (4)

  1. 　推定対象の装置である対象装置を制御するための制御入力値及び前記対象装置に設けられたセンサによって計測されたセンサ計測値に基づいて、センサによる計測が不可能な前記対象装置の状態量の推定値である状態推定値を出力する性能推定装置であって、
    　前記対象装置の非線形シミュレーションモデルを用いて、前記状態推定値、及び前記センサ計測値の推定値であるセンサ推定値を算出する推定部と、
    　前記センサ推定値と前記センサ計測値との誤差を算出する誤差算出部と、
    　前記対象装置の１つの作動点である第１作動点に対応する第１ゲイン、及び前記対象装置の別の作動点である第２作動点に対応する第２ゲインを記憶するゲイン記憶部と、
    　前記第１ゲイン及び前記第２ゲインに基づいて、前記対象装置が動作している作動点である設定作動点における設定ゲインを決定するゲイン決定部と、
    　前記誤差と前記設定ゲインとに基づいて、補正値を生成する補正部と、
    を備え、
    　前記推定部は、前記制御入力値及び前記補正値に基づいて、前記状態推定値及び前記センサ推定値を算出し、前記状態推定値を出力する、性能推定装置。
  2. 　前記ゲイン決定部は、前記第１作動点を含む前記作動点の範囲である第１範囲に前記設定作動点が含まれる場合には前記第１ゲインを前記設定ゲインとし、前記第２作動点を含む前記作動点の範囲である第２範囲に前記設定作動点が含まれる場合には前記第２ゲインを前記設定ゲインとする、請求項１に記載の性能推定装置。
  3. 　前記第１範囲と前記第２範囲とは互いに隣接しており、
    　前記ゲイン決定部は、前記設定作動点が前記第１範囲から前記第２範囲に移動した場合には、前記設定ゲインを前記第１ゲインから前記第２ゲインに徐々に切り替える、請求項２に記載の性能推定装置。
  4. 　推定対象の装置である対象装置に設けられたセンサによる計測が不可能な状態量の推定値である状態推定値を出力する性能推定方法であって、
    　前記対象装置が動作している作動点である設定作動点を受信するステップと、
    　前記対象装置を制御するための制御入力値を受信するステップと、
    　前記対象装置に設けられたセンサによって計測されたセンサ計測値を受信するステップと、
    　前記対象装置の非線形シミュレーションモデルを用いて、前記制御入力値に基づいて、前記センサ計測値の推定値であるセンサ推定値を算出するステップと、
    　前記センサ推定値と前記センサ計測値との誤差を算出するステップと、
    　前記対象装置の１つの作動点である第１作動点に対応する第１ゲイン、及び前記対象装置の別の作動点である第２作動点に対応する第２ゲインに基づいて、前記設定作動点における設定ゲインを決定するステップと、
    　前記誤差と前記設定ゲインとに基づいて、補正値を生成するステップと、
    　前記非線形シミュレーションモデルを用いて、前記制御入力値及び前記補正値に基づいて、前記状態推定値を算出するステップと、
    　前記状態推定値を出力するステップと、
    を含む、性能推定方法。
PCT/JP2017/015591 2016-09-26 2017-04-18 性能推定装置及び性能推定方法 WO2018055825A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016-187017 2016-09-26
JP2016187017A JP2020003831A (ja) 2016-09-26 2016-09-26 性能推定装置及び性能推定方法

Publications (1)

Publication Number Publication Date
WO2018055825A1 true WO2018055825A1 (ja) 2018-03-29

Family

ID=61690290

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/015591 WO2018055825A1 (ja) 2016-09-26 2017-04-18 性能推定装置及び性能推定方法

Country Status (2)

Country Link
JP (1) JP2020003831A (ja)
WO (1) WO2018055825A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020142984A1 (zh) * 2019-01-10 2020-07-16 大连理工大学 一种基于误差区间观测器的航空发动机主动容错控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005248946A (ja) * 2004-03-02 2005-09-15 General Electric Co <Ge> ガスタービンエンジンのためのモデルベース制御システム及び方法
JP2009068359A (ja) * 2007-09-11 2009-04-02 Japan Aerospace Exploration Agency ガスタービンエンジンの性能推定システム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005248946A (ja) * 2004-03-02 2005-09-15 General Electric Co <Ge> ガスタービンエンジンのためのモデルベース制御システム及び方法
JP2009068359A (ja) * 2007-09-11 2009-04-02 Japan Aerospace Exploration Agency ガスタービンエンジンの性能推定システム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MASAFUMI HAMAGUCHI ET AL.: "Damping Control of Sloshing in a Liquid Container on a Cart with an Active Vibration Reducer of 6-DOF Parallel Linkage:-The Case of a Straight Path on Horizontal Plane", TRANSACTIONS OF THE SOCIETY OF INSTRUMENT AND CONTROL ENGINEERS, vol. 51, no. 2, 15 February 2015 (2015-02-15), pages 101 - 109 *
ZHANG, MENGTING ET AL.: "A Digital Dual-Mode LQR Optimal Controller for Boost Converter over a Wide Operation Range", PROCEEDINGS OF THE 35TH CHINESE CONTROL CONFERENCE, 27 July 2016 (2016-07-27), pages 5890 - 5895, XP032952331 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020142984A1 (zh) * 2019-01-10 2020-07-16 大连理工大学 一种基于误差区间观测器的航空发动机主动容错控制方法
US11635734B2 (en) 2019-01-10 2023-04-25 Dalian University Of Technology Interval error observer-based aircraft engine active fault tolerant control method

Also Published As

Publication number Publication date
JP2020003831A (ja) 2020-01-09

Similar Documents

Publication Publication Date Title
US8068997B2 (en) Continuous performance analysis system and method
US9878692B2 (en) Model-based optimal control for stall margin limit protection in an aircraft engine
US8849542B2 (en) Real time linearization of a component-level gas turbine engine model for model-based control
JP5583455B2 (ja) タービン部品にサージ保護を与えるためのシステム及び方法
JP2009008078A (ja) マルチ缶型燃焼器に対して燃焼ダイナミックス調整アルゴリズムを用いるためのシステム及び方法
KR102290578B1 (ko) 가스 터빈 제어 장치, 가스 터빈 플랜트 및 가스 터빈 제어 방법
EP2770390B1 (en) System and method for continuous performance analysis of systems that exhibit variable performance characteristics at different operating conditions
CN110647052B (zh) 一种变循环发动机模式切换自适应身份证模型构建方法
JP2009068359A (ja) ガスタービンエンジンの性能推定システム
JP6236979B2 (ja) ガスタービンエンジン最適制御装置
EP2993309A1 (en) Engine-induced aircraft cabin resonance reduction system and method
WO2018055825A1 (ja) 性能推定装置及び性能推定方法
US7111464B2 (en) Acceleration control in multi spool gas turbine engine
EP2886835A1 (en) Gas turbine firing temperature control system and method
EP2647811A1 (en) Gas turbine control device and power generation system
US11015531B2 (en) Engine-induced aircraft cabin resonance reduction system and method
JP6483571B2 (ja) 性能推定装置及び性能推定方法
US20120158374A1 (en) Method for real-time modeling of an n-dimensional surface
JP2018142200A (ja) 状態推定方法、及び状態推定装置
CN113075882A (zh) 基于非平衡态线性化的涡扇发动机多变量鲁棒变增益控制方法
RU2665011C2 (ru) Способ и система регулирования рабочей температуры устройства создания давления наддува транспортного средства и транспортное средство, содержащее систему для регулирования температуры
JP6816646B2 (ja) ゲイン決定装置、ゲイン決定方法、及び性能推定装置
JP5888947B2 (ja) 弁制御装置、ガスタービン、及び弁制御方法
US20230086029A1 (en) System and method for non-model based control utilizing turbine exit mach number surrogate
CN113741195A (zh) 一种航空发动机非线性控制方法及系统

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: 17852606

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: 17852606

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP