(57)【特許請求の範囲】
【請求項１】
非定常状態の運転状態を含む監視対象物の異常診断方法であって、
前記監視対象物のシミュレーションモデルを作成するモデル作成ステップと、
前記監視対象物の運転状態における内部状態量を計測して実測値を抽出する計測ステップと、
前記監視対象物の運転状態と同一の制御入力値を前記シミュレーションモデルにインプットして前記監視対象物の内部状態量の予測値を算出する予測ステップと、
前記実測値と前記予測値との差分からマハラノビス距離を算出するマハラノビス距離算出ステップと、
前記マハラノビス距離に基づいて前記監視対象物の運転状態が異常であるか否か診断する異常診断ステップと、
を備えることを特徴とする異常診断方法。
【請求項２】
前記マハラノビス距離算出ステップは、前記差分とその積分値とを成分とするエラーベクトルを算出するステップを含む、ことを特徴とする請求項１に記載の異常診断方法。
【請求項３】
前記予測ステップは、時系列的に一つ前の実測値に基づいて前記予測値を算出する、ことを特徴とする請求項２に記載の異常診断方法。
【請求項４】
非定常状態の運転状態を含む監視対象物の異常診断システムであって、
前記監視対象物を模擬したシミュレーションモデルと、
前記監視対象物の運転状態における内部状態量を計測する計測手段と、
前記シミュレーションモデルにより求められた予測値と前記計測手段から抽出された実測値との差分からマハラノビス距離を算出するとともに該マハラノビス距離に基づいて前記監視対象物の運転状態が異常であるか否か診断する診断装置と、
少なくとも前記監視対象物及び前記シミュレーションモデルに同一の制御入力値を送信する制御装置と、
を備えることを特徴とする異常診断システム。
【請求項５】
前記診断装置は、前記差分とその積分値とを成分とするエラーベクトルに基づいて前記マハラノビス距離を算出する、ことを特徴とする請求項４に記載の異常診断システム。
【請求項６】
前記シミュレーションモデルは、時系列的に一つ前の実測値に基づいて前記予測値を算出する、ことを特徴とする請求項５に記載の異常診断システム。
【請求項７】
前記監視対象物は、再使用型宇宙機用エンジンである、ことを特徴とする請求項４に記載の異常診断システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、異常診断方法及び異常診断システムに関し、特に、動的変化を有する非定常状態における異常の発見に適した異常診断方法及び異常診断システムに関する。
【背景技術】
【０００２】
ガスタービン発電プラント、原子力発電プラント、火力発電プラント等の各種プラントやジェットエンジン等の内燃機関の分野では、安定した運転や出力を保持するために、運転状態（試験を含む）を監視して異常診断することが行われている。
【０００３】
例えば、特許文献１には、監視対象から所定の監視対象データを取得し、そのマハラノビス距離を算出して監視対象の異常を検知する監視手段と、異常の予兆を示す異常信号と関連する監視対象データである関連信号とを抽出して所定の入力信号を生成するデータ処理手段と、入力信号に基づいて監視対象の故障診断を行う故障診断手段と、を備え、監視対象における異常の予兆の監視から故障診断までの一連の処理を自動化できる監視システムが開示されている。
【０００４】
また、特許文献２には、過去に入力した各変数の値を含む蓄積データを記憶する蓄積データ記憶部と、各変数について新たに入力した値と蓄積データに含まれる所定期間の最大値及び最小値を抽出して中間の値を中央値として決定する決定部と、各変数について新たに入力した値と中央値との差分値を求める第１算出部と、求められた各変数の差分値と所定の単位空間のデータを利用してマハラノビス距離を求める第２算出部と、マハラノビス距離が予め定められる閾値の範囲内であるかを判定し異常を診断する判定部と、を備え、プラントから新たに入力する複数の変数の値を所定の単位空間と比較して前記プラントの異常を診断する異常診断装置が開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１１−０９０３８２号公報
【特許文献２】特開２０１４−０３５２８２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、プラントや内燃機関等の監視対象物は、一般に、安定した運転状態を示す定常状態と、定常状態に至るまでの過渡的な不安定な運転状態を示す非定常状態と、を有している。非定常状態は、同一の監視対象物であっても、その時の環境条件や運転条件等によって異なるものであり、同一の動的変化を示すことはほとんどない。
【０００７】
特許文献１に記載された監視システムでは、監視対象データのマハラノビス距離を算出することにより異常の予兆を示す異常信号及び関連信号から故障診断の入力信号を生成している。ここで、監視対象データのマハラノビス距離を算出した後、異常又はその予兆を示しているか否か判断するには、予め基準データを準備しておく必要がある。このとき、定常状態の場合には、運転状態や出力状態が安定していることから、基準データを準備することが可能である。しかしながら、動的変化を伴う非定常状態の場合には、監視対象データのみからは基準データを作成することができず、異常診断を行うことができない。
【０００８】
また、特許文献２に記載された異常診断装置においても、過去の蓄積データを用いてマハラノビス距離を算出していることから、特許文献１と同様に、定常状態の場合には過去のデータと比較して異常診断できるものの、非定常状態の場合には異常診断することができない。
【０００９】
本発明は、上述した問題点に鑑みて創案されたものであり、監視対象物の定常状態のみならず非定常状態についても異常診断することができる、異常診断方法及び異常診断システムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明によれば、非定常状態の運転状態を含む監視対象物の異常診断方法であって、前記監視対象物のシミュレーションモデルを作成するモデル作成ステップと、前記監視対象物の運転状態における内部状態量を計測して実測値を抽出する計測ステップと、前記監視対象物の運転状態と同一の制御入力値を前記シミュレーションモデルにインプットして前記監視対象物の内部状態量の予測値を算出する予測ステップと、前記実測値と前記予測値との差分からマハラノビス距離を算出するマハラノビス距離算出ステップと、前記マハラノビス距離に基づいて前記監視対象物の運転状態が異常であるか否か診断する異常診断ステップと、を備えることを特徴とする異常診断方法が提供される。
【００１１】
前記マハラノビス距離算出ステップは、前記差分とその積分値とを成分とするエラーベクトルを算出するステップを含んでいてもよい。さらに、前記予測ステップは、時系列的に一つ前の実測値に基づいて前記予測値を算出するようにしてもよい。
【００１２】
また、本発明によれば、非定常状態の運転状態を含む監視対象物の異常診断システムであって、前記監視対象物を模擬したシミュレーションモデルと、前記監視対象物の運転状態における内部状態量を計測する計測手段と、前記シミュレーションモデルにより求められた予測値と前記計測手段から抽出された実測値との差分からマハラノビス距離を算出するとともに該マハラノビス距離に基づいて前記監視対象物の運転状態が異常であるか否か診断する診断装置と、少なくとも前記監視対象物及び前記シミュレーションモデルに同一の制御入力値を送信する制御装置と、を備えることを特徴とする異常診断システムが提供される。
【００１３】
前記診断装置は、前記差分とその積分値とを成分とするエラーベクトルに基づいて前記マハラノビス距離を算出するようにしてもよい。さらに、前記シミュレーションモデルは、時系列的に一つ前の実測値に基づいて前記予測値を算出するようにしてもよい。また、前記監視対象物は、例えば、再使用型宇宙機用エンジンである。
【発明の効果】
【００１４】
本発明に係る異常診断方法及び異常診断システムによれば、監視対象物の内部状態を模擬するシミュレーションモデルを作成し、監視対象物の実測値とシミュレーションモデルの予測値との差分を用いて監視対象物の異常の有無を診断するようにしたことから、シミュレーションモデルにより異常診断時の環境条件や運転条件に適応した予測値を算出することができるとともに、差分を取ることにより監視対象物の実測値を正常値に対する変動値に置換することができる。したがって、監視対象物の運転状態が非定常状態の場合であっても、その動的変化に追従して対応することができ、監視対象物の定常状態のみならず非定常状態についても異常診断することができる。また、異常診断にマハラノビス距離を用いることにより、異常診断の簡素化及び高速化を図ることができる。
【図面の簡単な説明】
【００１５】
【図１】本発明に係る異常診断システムを示す概略全体構成図である。
【図２】本発明に係る異常診断方法を示すフロー図である。
【図３】マハラノビス距離算出ステップの説明図であり、（ａ）はエラーベクトル、（ｂ）は予測値の算出方法の一例、を示している。
【図４】異常診断ステップの説明図であり、（ａ）はマハラノビス距離の概念図、（ｂ）は異常診断の概念図、を示している。
【図５】本発明を再使用型宇宙機用エンジンに適用した場合の有効性を検証するための説明図であり、（ａ）は制御入力値、（ｂ）は異常な数値を意図的に含む模擬データ、（ｃ）はマハラノビス距離による異常診断結果、を示している。
【発明を実施するための形態】
【００１６】
以下、本発明の実施形態について図１〜図５を用いて説明する。ここで、図１は、本発明に係る異常診断システムを示す概略全体構成図である。図２は、本発明に係る異常診断方法を示すフロー図である。図３は、マハラノビス距離算出ステップの説明図であり、（ａ）はエラーベクトル、（ｂ）は予測値の算出方法の一例、を示している。図４は、異常診断ステップの説明図であり、（ａ）はマハラノビス距離の概念図、（ｂ）は異常診断の概念図、を示している。
【００１７】
本発明の一実施形態に係る異常診断システム１は、図１に示したように、非定常状態の運転状態を含む監視対象物２の異常診断システムであって、監視対象物２を模擬したシミュレーションモデル３と、監視対象物２の運転状態における所定の内部状態量を計測する計測手段４と、シミュレーションモデル３により求められた予測値ｘと計測手段４から抽出された実測値ｘ^（実際の表記はｘの上に^。以下、同じ。）との差分（ｘ^−ｘ）からマハラノビス距離ＭＤを算出するとともにマハラノビス距離ＭＤに基づいて監視対象物２の運転状態が異常であるか否か診断する診断装置５と、監視対象物２及びシミュレーションモデル３に同一の制御入力値ｕを送信する制御装置６と、を備えている。
【００１８】
監視対象物２は、例えば、再使用型宇宙機用エンジンであるが、これに限定されるものではなく、ジェットエンジン等の他の内燃機関や、ガスタービン発電プラント、原子力発電プラント、火力発電プラント、化学プラント等の各種プラント等であってもよい。特に、監視対象物２は、運転状態において、安定した運転状態を示す定常状態と、定常状態に至るまでの過渡的な不安定な運転状態を示す非定常状態と、を含むことが好ましい。
【００１９】
シミュレーションモデル３は、監視対象物２の内部状態量を推定可能な内部推定モデルであり、例えば、数値シミュレーション技術を応用して作成される。シミュレーションモデル作成時には、リアルタイムでの処理を考慮し、漸化式表現(ARMA)で記述してもよい。例えば、監視対象物２が再使用型宇宙機用エンジンの場合には、内部状態量として、例えば、燃焼圧Ｐｃ、再生冷却出口温度Ｔｊｍｆ、燃料ポンプ回転数Ｎｆ、酸化剤ポンプ回転数Ｎｏ、燃料ポンプ出口圧力Ｐｄｆ、酸化剤ポンプ出口圧力Ｐｄｏ等が選択される。したがって、これらの内部状態量を算出可能なシミュレーションモデルが作成される。シミュレーションモデル３は、監視対象物２の全体を模擬した一つのシミュレーションモデルであってもよいし、内部状態量を個別に算出可能な複数のシミュレーションモデルにより構築されていてもよい。
【００２０】
計測手段４は、監視対象物２に設置され、例えば、燃焼圧Ｐｃ、再生冷却出口温度Ｔｊｍｆ、燃料ポンプ回転数Ｎｆ、酸化剤ポンプ回転数Ｎｏ、燃料ポンプ出口圧力Ｐｄｆ、酸化剤ポンプ出口圧力Ｐｄｏ等の内部状態量を計測するセンサである。例えば、計測手段４は、圧力計、温度計、ロータリエンコーダ等であるが、これらに限定されるものではなく、監視対象物２によって、計測する内部状態量によって、適宜選択される。
【００２１】
制御装置６は、監視対象物２を運転するために必要な制御入力値ｕを送信する装置である。この監視対象物２の運転は、実運用であってもよいし、監視対象物２の試験であってもよい。また、制御装置６は、監視対象物２の運転に必要な制御入力値ｕをシミュレーションモデル３にも送信する。シミュレーションモデル３は、この制御入力値ｕに基づいて内部状態量を計算し、各内部状態量について予測値ｘを算出する。なお、制御入力値ｕにより運転された監視対象物２の出力値ｙを計測して外部に抽出するようにしてもよい。
【００２２】
診断装置５は、計測手段４により計測された実測値ｘ^のデータと、シミュレーションモデル３により算出された予測値ｘのデータと、を受信し、これらのデータを用いて監視対象物２の異常診断を行う装置である。かかる診断装置５では、例えば、図２に記載したフロー図に基づいて処理される。なお、診断結果及び診断データは、診断装置５からモニタ出力、紙出力、データ出力等によって外部出力するようにしてもよい。
【００２３】
ここで、図２に記載したフロー図は、監視対象物２のシミュレーションモデル３を作成するモデル作成ステップＳｔｅｐ１と、監視対象物２の運転を開始する運転開始ステップＳｔｅｐ２と、監視対象物２の運転状態における内部状態量を計測して実測値ｘ^を抽出する計測ステップＳｔｅｐ３と、監視対象物２の運転状態と同一の制御入力値ｕをシミュレーションモデル３にインプットして監視対象物２の内部状態量の予測値ｘを算出する予測ステップＳｔｅｐ４と、実測値ｘ^と予測値ｘとの差分（ｘ^−ｘ）からマハラノビス距離ＭＤを算出するマハラノビス距離算出ステップＳｔｅｐ５と、マハラノビス距離ＭＤに基づいて監視対象物２の運転状態が異常であるか否か診断する異常診断ステップＳｔｅｐ６と、を備えている。
【００２４】
上述した診断装置５では、マハラノビス距離算出ステップＳｔｅｐ５及び異常診断ステップＳｔｅｐ６の処理を行う。本実施形態における異常診断方法では、マハラノビス距離を用いた多変量解析に基づいて、得られたデータ（実測値ｘ^）が異常であるか否かを診断している。このマハラノビス距離を用いることにより、複数の変数の相関関係を一度に処理することができ、個々の変数について個別に異常であるか否かを診断する必要がなく、異常診断の簡素化及び高速化を図ることができる。
【００２５】
また、マハラノビス距離算出ステップＳｔｅｐ５は、図２に示したように、実測値ｘ^と予測値ｘとの差分（ｘ^−ｘ）を算出する差分算出ステップＳｔｅｐ５１と、差分（ｘ^−ｘ）と誤差の積分値Σεとを成分とするエラーベクトルεを算出するエラーベクトル算出ステップＳｔｅｐ５２と、エラーベクトルεに基づいてマハラノビス距離ＭＤを計算するマハラノビス距離計算ステップＳｔｅｐ５３と、を含んでいてもよい。
【００２６】
ここで、エラーベクトルεは、例えば、図３（ａ）に示したように表記することができる。エラーベクトルεの一成分を構成する積分値Σεは、時々刻々と変化するエラーベクトルを連続的に計算すれば、いわゆる積分値として算出できるものであり、一定の時間（スパン）毎にエラーベクトルεを計算する場合には、積分値Σεは差分（ｘ^−ｘ）の総和として算出することができる。このように、誤差（差分）の積分値Σεを用いることにより、同一方向への累積誤差評価感度が脆弱となることを防ぐことができる。
【００２７】
例えば、内部状態量として、燃焼圧Ｐｃ、再生冷却出口温度Ｔｊｍｆ、燃料ポンプ回転数Ｎｆ、酸化剤ポンプ回転数Ｎｏ、燃料ポンプ出口圧力Ｐｄｆ、酸化剤ポンプ出口圧力Ｐｄｏを選択した場合には、エラーベクトルεは、図３（ａ）に示したように、（ΔＰｃ，ΔＴｊｍｆ，ΔＮｆ，ΔＮｏ，ΔＰｄｆ，ΔＰｄｏ，ΣΔＰｃ，ΣΔＴｊｍｆ，ΣΔＮｆ，ΣΔＮｏ，ΣΔＰｄｆ，ΣΔＰｄｏ）の行列として表記することもできる。この場合、エラーベクトルεは、１２個の変数を含んでいることから、これらの変数により形成されるベクトル空間Ｒ１２に含まれる。
【００２８】
また、予測ステップＳｔｅｐ４は、監視対象物２の運転と同一の制御入力値ｕをシミュレーションモデル３にインプットする入力ステップＳｔｅｐ４１と、制御入力値ｕに基づいて内部状態量の予測値ｘを算出する予測値算出ステップＳｔｅｐ４２と、を含んでいる。予測値算出ステップＳｔｅｐ４２（予測ステップＳｔｅｐ４）は、図３（ｂ）に示したように、時系列的に一つ前の実測値ｘｎ−１^に基づいて予測値ｘｎを算出するようにしてもよい。すなわち、予測値ｘｎは実測値ｘｎ−１^に基づいて算出され、予測値ｘｎ＋１は実測値ｘｎ^に基づいて算出される。かかる処理により、誤差の累積を抑制することができ、予測値ｘｎの精度を向上させることができ、異常診断の精度を向上させることができる。
【００２９】
マハラノビス距離計算ステップＳｔｅｐ５３において、エラーベクトルεからマハラノビス距離ＭＤを計算するには、まず、数１を用いてエラーベクトルεを規格化し、物理量単位に依存しない状態に変換する。エラーベクトルεの規格化には、運転期間中の全平均値ベクトルε￣（実際の表記はεの上に￣。以下、同じ。）及び偏差σεを用いる。
【００３０】
【数１】
0006610987.tif 000002
【００３１】
次に、数２を用いて、マハラノビス距離ＭＤを計算する。ここで、εＴは、エラーベクトルεの転置行列を意味し、ｄｉｍ（ε）は、エラーベクトルεの次元を意味している。また、共分散行列は、例えば、正常であると診断された過去の蓄積データから導出することができる。
【００３２】
【数２】
0006610987.tif 000003
【００３３】
マハラノビス距離ＭＤを算出し、等距離の点を結ぶことにより、例えば、図４（ａ）に示したような内部状態量の相関関係を求めることができ、図示した略楕円領域の中心から離れるに従って誤差が大きく、この領域から逸脱した場合に異常であると診断することができる。ここで、図４（ａ）に示した相関関係は、直感的理解を促すために、内部状態量Ｄ１，Ｄ２の２変数のみの相関関係を示している。この相関関係によれば、略楕円領域の長径方向に対しては誤差の許容量が大きく、略楕円領域の短径方向に対しては誤差の許容量が小さいことが理解できる。なお、図示しないが、上述したように、１２個の変数を用いた場合には、１２次元の相関関係を求めることとなる。
【００３４】
異常診断ステップＳｔｅｐ６では、例えば、図４（ｂ）に示したように、時々刻々と変化する誤差（差分値）に対して、診断時毎にマハラノビス距離ＭＤを計算し、その都度、誤差（差分値）がマハラノビス距離ＭＤの範囲内であるかを判断する。例えば、時間ｔ１におけるマハラノビス距離ＭＤ１、時間ｔ２におけるマハラノビス距離ＭＤ２、時間ｔ３におけるマハラノビス距離ＭＤ３、時間ｔ４におけるマハラノビス距離ＭＤ４、時間ｔ５におけるマハラノビス距離ＭＤ５は、その時の環境条件や運転条件等によって時々刻々と変化するものである。なお、図４（ｂ）に示した図は、本実施形態に係る異常診断方法の直感的理解を促すために図示したものである。
【００３５】
上述した本実施形態に係る異常診断方法及び異常診断システム１によれば、監視対象物２の内部状態を模擬するシミュレーションモデル３を作成し、監視対象物２の実測値ｘ^とシミュレーションモデル３の予測値ｘとの差分（ｘ^−ｘ）を用いて監視対象物２の異常の有無を診断するようにしたことから、シミュレーションモデル３により異常診断時の環境条件や運転条件に適応した予測値ｘを算出することができるとともに、差分を取ることにより監視対象物２の実測値ｘ^を正常値に対する変動値に置換することができる。したがって、監視対象物２の運転状態が非定常状態の場合であっても、その動的変化に追従して対応することができ、監視対象物２の定常状態のみならず非定常状態についても異常診断することができる。
【００３６】
ここで、図５は、本発明を再使用型宇宙機用エンジンに適用した場合の有効性を検証するための説明図であり、（ａ）は制御入力値、（ｂ）は異常な数値を意図的に含む模擬データ、（ｃ）はマハラノビス距離による異常診断結果、を示している。図５（ａ）及び（ｂ）において、推力の数値を実線、燃料の数値を一点鎖線、酸化剤の数値を点線、燃焼圧の数値を二点鎖線、で表示している。なお、図５（ａ）において、推力が上に凸となっている部分（略台形部分）は非定常状態を模擬したものである。
【００３７】
図５（ａ）に示した推力を得るために、図示したように、燃料及び酸化剤の分量を制御するものとする。いま、マハラノビス距離ＭＤによる異常診断の有効性を検証するために、図５（ｂ）に示したように、正常な実測値に対してオフセット値（図中のαの部分）を与えて、異常な数値を意図的に含む模擬データを作成した。そして、この模擬データとシミュレーションモデル３により求められる予測値を用いて、上述したマハラノビス距離算出ステップＳｔｅｐ５の処理を行ったところ、図５（ｃ）に示した結果が得られた。
【００３８】
図５（ｃ）において、実線はマハラノビス距離ＭＤの時間的変化を示しており、●（黒丸）は異常であると診断された点を示している。この検証結果によれば、意図的に異常な数値を付与したオフセット部分に対応した部分のマハラノビス距離ＭＤについて異常であると診断されていることが理解できる。したがって、本実施形態に係る異常診断方法及び異常診断システム１は、非定常状態を含む運転状態に対して異常診断可能な応答性を有していることが認められる。
【００３９】
本発明は上述した実施形態に限定されず、本発明の趣旨を逸脱しない範囲で種々変更が可能であることは勿論である。
【符号の説明】
【００４０】
１  異常診断システム
２  監視対象物
３  シミュレーションモデル
４  計測手段
５  診断装置
６  制御装置
Ｓｔｅｐ１  モデル作成ステップ
Ｓｔｅｐ２  運転開始ステップ
Ｓｔｅｐ３  計測ステップ
Ｓｔｅｐ４  予測ステップ
Ｓｔｅｐ４１  入力ステップ
Ｓｔｅｐ４２  予測値算出ステップ
Ｓｔｅｐ５  マハラノビス距離算出ステップ
Ｓｔｅｐ５１  差分算出ステップ
Ｓｔｅｐ５２  エラーベクトル算出ステップ
Ｓｔｅｐ５３  マハラノビス距離計算ステップ
Ｓｔｅｐ６  異常診断ステップ
【図１】
0006610987.tif 000004
【図２】
0006610987.tif 000005
【図３】
0006610987.tif 000006
【図４】
0006610987.tif 000007
【図５】
0006610987.tif 000008