(57)【特許請求の範囲】
【請求項１】
発電プラントから状態量を計測した計測信号に基づいてプラントの運転状態を診断し、診断結果を画像表示装置に表示する発電プラントの診断装置において、
発電プラントの診断装置に発電プラントの状態量を計測した計測信号を用いて診断に用いるモデルを構築するモデル構築手段と、前記モデルで診断する運転条件と計測信号の正規化方法を定義するモデル定義手段と、前記モデル構築手段で構築したモデルを用いて発電プラントの運転状態を診断する診断手段を備え、
前記モデル定義手段に、発電プラントの運転条件を判定する運転条件判定部と、運転条件判定部で判定した運転条件毎に計測信号の正規化条件を決定する正規化条件決定部を備え、
前記診断手段では、運転条件に合わせて診断モデルを切り替えて診断することを特徴とする発電プラントの診断装置。
【請求項２】
請求項１に記載した発電プラントの診断装置において、
前記運転条件判定部は、発電プラントの運転条件を一定負荷運転中、起動運転中、停止運転中、負荷変化運転中のいずれかに分類する演算装置と、運転条件が一定負荷運転中の場合は負荷帯を抽出する演算装置と、運転条件が起動運転中の場合は起動モードの種類と負荷変化率を抽出する演算装置と、停止運転中の時は停止モードの種類と停止動作開始出力を抽出する演算装置と、負荷変化運転中の時は負荷変化率、負荷変化開始時の出力、負荷変化終了時の出力を抽出する演算装置を備えることを特徴とする発電プラントの診断装置。
【請求項３】
請求項１に記載した発電プラントの診断装置において、
前記正規化条件決定部では、前記運転条件判定部で判定した運転条件毎にモデル構築用データを分割する演算装置と、運転条件毎のデータ変化幅の情報を抽出する演算装置と、前記データ変化幅の情報に基づいて正規化範囲を決定する演算装置を備えることを特徴とする発電プラントの診断装置。
【請求項４】
請求項３に記載した発電プラントの診断装置において、
負荷変化中における計測信号の正規化範囲は、出力指令値の関数とする演算装置を備えることを特徴とする発電プラントの診断装置。
【請求項５】
請求項３に記載した発電プラントの診断装置において、
前記正規化条件決定部では、負荷変化中の計測信号について無駄時間を推定する演算装置と、推定した無駄時間分計測信号をシフトさせる演算装置を動作させた後、
請求項３に記載した演算装置を動作させることを特徴とする発電プラントの診断装置。
【請求項６】
請求項３〜請求項５に記載した発電プラントの診断装置において、
前記正規化条件決定部では、計測信号の最大値、最小値付近での正規化範囲を広くする演算装置を備えたことを特徴とする発電プラントの診断装置。
【請求項７】
請求項１に記載した発電プラントの診断装置において、
計測信号の経時変化を示すトレンドグラフと、前記運転条件判定手段で判定した運転条件が切り替わる時刻と、前記正規化条件決定手段で決定した正規化範囲を重ねて表示させ、
前記運転条件が切り替わる時刻と、前記正規化範囲を任意に変更するための画面表示装置を備えたことを特徴とする発電プラントの診断装置。
【請求項８】
発電プラントから状態量を計測した計測信号に基づいてプラントの運転状態を診断し、診断結果を画像表示装置に表示する発電プラントの診断方法において、
発電プラントの運転条件を判定するステップと、運転条件判定部で判定した運転条件毎に計測信号の正規化条件を決定するステップと、発電プラントの診断装置に発電プラントの状態量を計測した計測信号を用いて診断に用いるモデルを構築するステップと、前記モデルで診断する運転条件と計測信号の正規化方法を定義するステップと、運転条件に合わせて診断モデルを切り替えて診断するステップを備えたことを特徴とする発電プラントの診断方法。
【請求項９】
請求項８に記載した発電プラントの診断方法において、
前記発電プラントの運転条件を判定するステップは、発電プラントの運転条件を一定負荷運転中、起動運転中、停止運転中、負荷変化運転中のいずれかに分類するステップと、運転条件が一定負荷運転中の場合は負荷帯を抽出するステップと、運転条件が起動運転中の場合は起動モードの種類と負荷変化率を抽出するステップと、停止運転中の時は停止モードの種類と停止動作開始出力を抽出するステップと、負荷変化運転中の時は負荷変化率、負荷変化開始時の出力、負荷変化終了時の出力を抽出するステップを備えることを特徴とする発電プラントの診断方法。
【請求項１０】
請求項８に記載した発電プラントの診断方法において、
前記正規化条件を決定するステップでは、前記運転条件判定部で判定した運転条件毎にモデル構築用データを分割するステップと、運転条件毎のデータ変化幅の情報を抽出するステップと、前記データ変化幅の情報に基づいて正規化範囲を決定するステップを備えることを特徴とする発電プラントの診断方法。
【請求項１１】
請求項１０に記載した発電プラントの診断方法において、
負荷変化中における計測信号の正規化範囲は、出力指令値の関数とするステップを備えることを特徴とする発電プラントの診断方法。
【請求項１２】
請求項１０に記載した発電プラントの診断方法において、
前記正規化条件を決定するステップでは、負荷変化中の計測信号について無駄時間を推定するステップと、推定した無駄時間分計測信号をシフトさせるステップを実行した後、請求項１０に記載したステップを実行することを特徴とする発電プラントの診断方法。
【請求項１３】
請求項１０〜請求項１２に記載した発電プラントの診断方法において、
前記正規化条件を決定するステップでは、計測信号の最大値、最小値付近での正規化範囲を広くするステップを備えたことを特徴とする発電プラントの診断方法。
【請求項１４】
請求項８に記載した発電プラントの診断方法において、
計測信号の経時変化を示すトレンドグラフと、前記運転条件判定手段で判定した運転条件が切り替わる時刻と、前記正規化条件決定手段で決定した正規化範囲を重ねて表示させ、
前記運転条件が切り替わる時刻と、前記正規化範囲を任意に変更できることを特徴とする発電プラントの診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、発電プラントの診断装置、及び発電プラントの診断方法に関する。
【背景技術】
【０００２】
発電プラントの診断装置は、プラントに異常な過渡事象や事故等が生じた際に、プラントからの計測信号を基にその異常や事故の発生を検知する。
【０００３】
公知例のプラントの診断装置として特開２００５−１６５３７５号公報には、適応共鳴理論（Adaptive Resonance Theory：ＡＲＴ）を用いたプラントの診断装置が開示されている。ここで、ＡＲＴとは、多次元のデータをその類似度に応じてカテゴリーに分類する技術である。
【０００４】
この特開２００５−１６５３７５号公報に記載されたプラントの診断装置の技術においては、まず、プラントの運転データを記録した過去の計測信号から、プラントの状態が正常と考えられる期間の計測信号をモデル構築用データとして抽出する。そして、ＡＲＴを用いて、モデル構築用データを複数のカテゴリー（正常カテゴリー）に分類して診断に用いる正常モデルを作成する。次に、プラントの現在の計測信号をＡＲＴでカテゴリーに分類する。この現在の計測信号が正常モデルと一致しない時、すなわち正常カテゴリーに分類できない時は、新しいカテゴリー（新規カテゴリー）が生成される。つまり、新規カテゴリーの発生は、計測信号の傾向が変化し、プラントの状態が変化したことを意味する。そこで、異常の発生を新規カテゴリーの発生で判断することとし、新規カテゴリーの発生割合が閾値を越えた場合に異常と診断する技術である。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００５−１６５３７５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
発電プラントは、起動、停止、一定負荷、負荷変化など、様々な条件で運転する。運転条件が異なると、計測信号の変化する範囲も異なる。
【０００７】
また、診断に用いる計測信号の前処理として、正規化処理がある。正規化処理では、正規化の下限値が０、正規化の上限値が１になるように、計測信号を処理する。正規化の下限値、上限値は、事前に設定する必要がある。先行技術の診断装置では、運転条件に関わらず、同じ正規化条件で計測信号を処理する。そのため、正規化範囲は全ての運転条件における計測信号の変化範囲を包含するように、広く設定する必要がある。
【０００８】
計測信号の変化範囲に比べて正規化範囲が広いと、計測信号を正規化した後の値の変化が小さくなる。そのため、異常発生時の計測信号の傾向変化をとらえることができず、異常発生時にも新規カテゴリーが発生しない場合があった。これは、失報の原因となる。
【０００９】
本発明の目的は、運転条件に合わせて正規化範囲を適切に決定することで失報を抑制し、診断精度を向上させた診断装置を提供することにある。
【課題を解決するための手段】
【００１０】
発電プラントから状態量を計測した計測信号に基づいてプラントの運転状態を診断し、診断結果を画像表示装置に表示する発電プラントの診断装置において、発電プラントの診断装置に発電プラントの状態量を計測した計測信号を用いて診断に用いるモデルを構築するモデル構築手段と、前記モデルで診断する運転条件と計測信号の正規化方法を定義するモデル定義手段と、前記モデル構築手段で構築したモデルを用いて発電プラントの運転状態を診断する診断手段を備え、前記モデル定義手段に、発電プラントの運転条件を判定する運転条件判定部と、運転条件判定部で判定した運転条件毎にデータの正規化条件を決定する正規化条件決定部を備え、前記診断手段では、運転条件に合わせて診断モデルを切り替えて診断することを特徴とする。
【発明の効果】
【００１１】
本発明の発電プラントの診断装置を用いることにより、異常発生時に異常を検知できない失報を低減し、診断精度を向上できる。また、正規化範囲を自動的に決定することができ、診断装置の調整期間を短縮できる。
【図面の簡単な説明】
【００１２】
【図１】本発明の一実施例である発電プラントの診断装置の構成を示す制御ブロック図。
【図２】図１に示した発電プラントの診断装置の基本動作を示すフローチャート図、及び動作のタイミングを示した説明図。
【図３】図１に示したプラントの診断装置におけるモデル構築手段、診断手段において、データを分類する機能の実装例を示す説明図。
【図４】図１に示した発電プラントの診断装置におけるモデル構築手段で計測信号を分類する実例を示す説明図。
【図５】図４（ａ）に示した発電プラントにおける起動モードとプロセス値の関係の説明図。
【図６】図１に示した発電プラントの診断装置における運転条件判定部５００の動作の説明図。
【図７】図１に示した発電プラントの診断装置における正規化条件決定部６００の第１の実施例の説明図。
【図８】図１に示した発電プラントの診断装置における正規化条件決定部６００の第２の実施例の説明図。
【図９】図１に示した発電プラントの診断装置における正規化条件決定部６００の第３の実施例の説明図。
【図１０】図１に示した発電プラントの診断装置におけるモデル構築モードと診断モードの動作フローチャートの説明図。
【図１１】図１に示した発電プラントの診断装置におけるデータベースに保存されているデータ態様の説明図。
【図１２】図１に示した発電プラントの診断装置の適用効果の説明図。
【図１３】図１に示した発電プラントの診断装置における画像表示装置に表示される画面の説明図。
【発明を実施するための形態】
【００１３】
次に、本発明の実施例である発電プラントの診断装置について図面を参照して以下に説明する。
【００１４】
図１は、本発明の一実施例である発電プラントの診断装置を説明するブロック図である。図１に示した発電プラントの診断装置においては、プラント１００の状態を診断装置２００により診断する。
【００１５】
診断装置２００は、診断装置２００を構成する演算装置として、モデル定義手段４００、モデル構築手段７００、及び診断手段８００を備えている。この診断装置２００は、データベースとして、計測信号データベース３１０、モデル定義データベース３２０、及び診断モデルデータベース３３０を備えている。なお、図１においては、データベースをＤＢと略記している。
【００１６】
計測信号データベース３１０、モデル定義データベース３２０、及び診断モデルデータベース３３０のデータベースには、電子化された情報が記録されており、通常、電子ファイル（電子データ）と呼ばれる。
【００１７】
モデル構築手段７００は、プラント１００の運転状態を計測した計測信号から、プラント１００の過去の運転状態の計測信号を蓄積した蓄積データを基に、プラントの正常状態を学習した診断モデルを作成する。
【００１８】
モデル定義手段４００は、診断モデルで診断する運転条件と、モデル構築時における計測信号の正規化方法を定義する。
【００１９】
診断手段８００は、モデル構築手段７００で作成した診断モデルの値のデータと、計測したプラント１００の計測信号のデータとを比較する。計測信号が、正常状態を学習した診断モデルの値と一致すれば、プラントの状態を正常と判定し、一致しなければ異常と判定する。
【００２０】
また、診断装置２００は、外部とのインターフェイスとして外部入力インターフェイス２１０及び外部出力インターフェイス２２０を備えている。
【００２１】
そして、外部入力インターフェイス２１０を介してプラント１００の運転状態である各種状態量を計測した計測信号１と、運転管理室９００に備えられているキーボード９２０及びマウス９３０で構成される外部入力装置９１０の操作で作成する外部入力信号２が診断装置２００に取り込まれる。また、外部出力インターフェイス２２０を介して、画像表示情報１１を運転管理室９００に備えられている画像表示装置９４０に出力する。
【００２２】
なお、本実施例のプラントの診断装置においては、モデル定義手段４００、モデル構築手段７００、診断手段８００、計測信号データベース３１０、モデル定義データベース３２０、診断モデルデータベース３３０が診断装置２００の内部に備えられているが、これらの一部を診断装置２００の外部に配置し、データのみをこれらの装置の間で通信するようにしてもよい。
【００２３】
また、本実施例のプラントの診断装置において診断対象とするプラント１００は１基である場合を示しているが、１台の診断装置２００で複数基のプラント１００を診断することも可能である。
【００２４】
次に、本実施例のプラントの診断装置に備えられた診断装置２００の動作を説明する。
【００２５】
図１に示した本実施例のプラントの診断装置において、プラント１００の各種状態量を計測した計測信号１は外部入力インターフェイス２１０を介して取り込まれる。計測信号３は、診断装置２００に設置された計測信号データベース３１０に保存される。
【００２６】
モデル定義手段４００には、運転条件判定部５００、および正規化条件決定部６００がそれぞれ備えられている。モデル定義手段４００は、計測信号４の入力に対し、モデル定義情報５をモデル定義データベース３２０に出力する。
【００２７】
発電プラントは出力を一定にして運転する一定負荷運転、プラントを起動する起動運転、プラントを停止する停止運転、出力を変化させる負荷変化運転がある。運転条件判定部５００では、計測信号データベース３１０に蓄積されているプラント１００の蓄積データを用いて、これらのデータを一定負荷運転中、起動運転中、停止運転中、負荷変化運転中の運転モードに分割する。また、それぞれの運転モードの特徴量を抽出する。この機能の詳細は、図６を用いて後述する。
【００２８】
また、モデル構築手段７００では、モデルを構築するための前処理として、計測信号を正規化処理する。正規化条件決定部６００では、各運転モードに対して、適切な正規化条件を決定する。この機能の詳細は、図７〜図９を用いて後述する。前述のモデル定義情報５は、運転条件データと正規化条件データで構成される。
【００２９】
モデル構築手段７００では、計測信号データベース３１０に蓄積されているプラント１００の計測信号４と、モデル定義データベース３２０に保存されているモデル定義情報６を用いて、診断に用いるモデルを構築する。モデル構築手段７００で作成したモデル情報８は、診断モデルデータベース３３０に保存する。
【００３０】
モデル構築手段７００を実装する技術として、適応共鳴理論、ベクトル量子化などのクラスタリング技術がある。また、診断に用いるモデルは上記のクラスタリング方法に限定されるわけではなく、物理式を用いたモデル、ニューラルネットワークなどの統計モデルを用いることもできる。
【００３１】
前記診断装置２００に設置された診断手段８００では、計測信号４の入力に対して、モデル定義データベース３２０のモデル定義情報７と、診断モデルデータベース３３０のモデル情報９を参照することによってプラント１００の運転状態を診断し、その診断結果１０を出力する。
【００３２】
診断手段８００が診断したプラント１００の現在の運転状態に対する診断結果１０は、外部出力インターフェイス２２０を介して、画像表示情報１１として運転管理室９００に設置した画像表示装置９４０に送信され、表示される。これにより、運転管理室９００にいるオペレータに、プラント１００の運転状態に対する診断結果を通知する。
【００３３】
このようにして、本実施例のプラントの診断装置２００では、プラントの状態が変化したことを、オペレータに通知する。
【００３４】
また、診断装置２００に設置された計測信号データベース３１０、モデル定義データベース３２０、診断モデルデータベース３３０に保存されている診断装置情報５０は、任意に運転管理室９００の画像表示装置９４０に表示できるようになっている。また、これらの情報は、キーボード９２０とマウス９３０で構成する外部入力装置９１０を操作して生成する外部入力信号２で、修正できる。
【００３５】
次に、本実施例のプラントの診断装置の動作について説明する。以下に診断装置２００の動作フローチャートを、図１に示したプラントの診断装置の基本動作を示すフローチャート図である図２（ａ）を用いて説明する。
【００３６】
図２（ａ）のフローチャート図に示すように、診断装置２００の基本動作はステップ２０１、２０２、２０３を組み合わせて実行する。
【００３７】
まず、ステップ２０１では診断装置２００の動作モードがモデル構築モードか診断モードかを判定する。そして、モデル構築モードの場合はステップ２０２に進み、診断モードの時はステップ２０３に進む。
【００３８】
ステップ２０２を動作させると、モデル定義手段４００、モデル構築手段７００が動作する。その結果、モデル定義情報５とモデル情報８が生成され、作成した情報はそれぞれモデル定義データベース３２０、診断モデルデータベース３３０に保存される。モデル構築モードの動作の詳細は、図１０（ａ）を用いて後述する。
【００３９】
また、ステップ２０３を動作させると、診断手段８００にてプラント１００の運転状態を診断し、診断結果１０の含まれる画像表示情報１１を画像表示装置９４０送信することにより、プラント１００の運転状態を画像表示装置９４０に表示する。診断モードの動作の詳細は、図１０（ｂ）を用いて後述する。
【００４０】
診断装置２００のモデル構築モードと診断モードを動作させるタイミングはオペレータによって任意に指定することができる。以下、図２（ｂ）〜（ｄ）を用いて、モデル構築モードと診断モードを動作させるタイミングの各種の実施例をそれぞれ説明する。
【００４１】
図２（ｂ）に示した実施例においては、計測信号のサンプリング周期毎にモデル構築モードと診断モードの両方を動作させて診断する。
【００４２】
計測信号を取得する毎に診断モデルを更新することで、常に最新のモデルを用いた診断が可能である。
【００４３】
しかし、モデル構築に用いるデータ量が多い時は、モデル構築に時間を要するため、サンプリング周期内に計算が終了しない可能性がある。
【００４４】
このような場合は、図２（ｃ）に示した実施例のように、所定の設定期間毎に正常状態モデル構築モードを動作させ、サンプリング周期毎に診断モードのみを動作させて診断することもできる。図２（ｂ）及び図２（ｃ）に示した実施例の方法では、サンプリング周期毎に診断モードが実行され、オンラインでプラントの状態を診断できる。
【００４５】
また、図２（ｄ）に示した実施例のように、オペレータがモデル構築、診断を実施する外部入力信号２を診断装置２００に入力することで、任意のタイミングでモデル構築モードと診断モードを動作できる。つまり、色々な条件を変えてプラント１００の運転状態を診断することが可能となる。
【００４６】
次に、図３、図４を用いて本実施例のプラントの診断装置を構成する診断装置２００を構成するモデル構築手段７００、及び診断手段８００に備えた、プラント１００の計測信号４を分類する機能を説明する。
【００４７】
本実施例のプラントの診断装置では、データ分類機能に適応共鳴理論（Adaptive Resonance Theory：ＡＲＴ）を適用した場合について述べる。尚、データ分類機能として、ベクトル量子化等、他のクラスタリング手法を用いることもできる。
【００４８】
図３（ａ）に示すように、データ分類機能はデータ前処理装置７１０とＡＲＴモジュール７２０で構成する。データ前処理装置７１０は、運転データをＡＲＴモジュール７２０の入力データに変換する。
【００４９】
以下に、前記データ前処理装置７１０及びＡＲＴモジュール７２０によるそれらの手順（工程）について説明する。
【００５０】
まず、データ前処理装置７１０において、モデル定義データベース３２０に保存されている正規化条件の情報を用いて、計測項目毎にデータを正規化する。計測信号を正規化したデータＮxi（ｎ）及び正規化したデータの補数ＣＮxi（ｎ）（＝１−Ｎxi（ｎ））を含むデータを入力データＩｉ（ｎ）とする。この入力データＩｉ（ｎ）が、ＡＲＴモジュール７２０に入力される。
【００５１】
ＡＲＴモジュール７２０においては、入力データであるプラント１００の計測信号４を複数のカテゴリーに分類する。
【００５２】
ＡＲＴモジュール７２０は、Ｆ０レイヤー７２１、Ｆ１レイヤー７２２、Ｆ２レイヤー７２３、メモリ７２４及び選択サブシステム７２５を備え、これらは相互に結合している。Ｆ１レイヤー７２２及びＦ２レイヤー７２３は、重み係数を介して結合している。重み係数は、入力データが分類されるカテゴリーのプロトタイプ（原型）を表している。ここで、プロトタイプとは、カテゴリーの代表値を表すものである。
【００５３】
次に、ＡＲＴモジュール７２０のアルゴリズムについて説明する。
【００５４】
ＡＲＴモジュール７２０に入力データが入力された場合のアルゴリズムの概要は、下記の処理１〜処理５のようになる。
【００５５】
処理１：Ｆ０レイヤー７２１により入力ベクトルを正規化し、ノイズを除去する。
【００５６】
処理２：Ｆ１レイヤー７２２に入力された入力データと重み係数との比較により、ふさわしいカテゴリーの候補を選択する。
【００５７】
処理３：選択サブシステム７２５で選択したカテゴリーの妥当性がパラメータρとの比により評価される。妥当と判断されれば、入力データはそのカテゴリーに分類され、処理４に進む。一方、妥当と判断されなければ、そのカテゴリーはリセットされ、他のカテゴリーからふさわしいカテゴリーの候補を選択する（処理２を繰り返す）。パラメータρの値を大きくするとカテゴリーの分類が細かくなり、ρの値を小さくすると分類が粗くなる。このパラメータρをビジランス（ｖｉｇｉｌａｎｃｅ）パラメータと呼ぶ。
【００５８】
処理４：処理２において全ての既存のカテゴリーがリセットされると、入力データが新規カテゴリーに属すると判断され、新規カテゴリーのプロトタイプを表す新しい重み係数を生成する。
【００５９】
処理５：入力データがカテゴリーＪに分類されると、カテゴリーＪに対応する重み係数ＷＪ（ｎｅｗ）は、過去の重み係数ＷＪ（ｏｌｄ）及び入力データｐ（又は入力データから派生したデータ）を用いて下記式（１）により更新される。
【００６０】
【数１】
0005337853.tif 000002
【００６１】
ここで、Ｋｗは、学習率パラメータ（０＜Ｋｗ＜１）であり、入力ベクトルを新しい重み係数に反映させる度合いを決定する値である。
【００６２】
尚、式（１）及び後述する式（２）乃至式（１２）の各演算式は前記ＡＲＴモジュール７２０に組み込まれている。
【００６３】
ＡＲＴモジュール７２０のデータ分類アルゴリズムの特徴は、上記の処理４にある。
【００６４】
処理４においては、学習した時のパターンと異なる入力データが入力された場合、記録されているパターンを変更せずに新しいパターンを記録することができる。このため、過去に学習したパターンを記録しながら、新たなパターンを記録することが可能となる。
【００６５】
このように、入力データとして予め与えた運転データを与えると、ＡＲＴモジュール７２０は与えられたパターンを学習する。したがって、学習済みのＡＲＴモジュール７２０に新たな入力データが入力されると、上記アルゴリズムにより、過去におけるどのパターンに近いかを判定することができる。また、過去に経験したことのないパターンであれば、新規カテゴリーに分類される。
【００６６】
図３（ｂ）は、Ｆ０レイヤー７２１の構成を示すブロック図である。Ｆ０レイヤー７２１では、入力データＩiを各時刻で再度正規化し、Ｆ１レイヤー７２１、及び選択サブシステム７２５に入力する正規化入力ベクトルｕi0を作成する。
【００６７】
始めに、入力データＩiから、式（２）に従ってｗi0を計算する。ここで、ａは定数である。
【００６８】
【数２】
0005337853.tif 000003
【００６９】
次に、ｗi0を正規化したｘi0を、式（３）を用いて計算する。ここで、‖・‖はノルムを表す記号である。
【００７０】
【数３】
0005337853.tif 000004
【００７１】
そして、式（４）を用いて、ｘi0からノイズを除去したｖi0を計算する。ただし、θはノイズを除去するための定数である。式（４）の計算により、微小な値は０となるため、入力データのノイズが除去される。
【００７２】
【数４】
0005337853.tif 000005
【００７３】
最後に、式（５）を用いて正規化入力ベクトルｕi0を求める。ｕi0はＦ１レイヤーの入力となる。
【００７４】
【数５】
0005337853.tif 000006
【００７５】
図３（ｃ）は、Ｆ１レイヤー７２２の構成を示すブロック図である。Ｆ１レイヤー７２２では、式（５）で求めたｕi0を短期記憶として保持し、Ｆ２レイヤー７２２に入力するｐiを計算する。Ｆ２レイヤーの計算式をまとめて式（６）〜式（１２）に示す。ただし、ａ、ｂは定数、f（・）は式（４）で示した関数、ＴjはＦ２レイヤー７２２で計算する適合度である。
【００７６】
【数６】
0005337853.tif 000007
【００７７】
【数７】
0005337853.tif 000008
【００７８】
【数８】
0005337853.tif 000009
【００７９】
【数９】
0005337853.tif 000010
【００８０】
【数１０】
0005337853.tif 000011
【００８１】
【数１１】
0005337853.tif 000012
【００８２】
ただし、
【００８３】
【数１２】
0005337853.tif 000013
【００８４】
次に、図４を用いて本実施例のプラントの診断装置を構成する診断装置２００を構成するモデル構築手段７００に備えた、プラント１００の計測信号４でモデルを構築する機能を説明する。
【００８５】
まず、図４（ａ）を用いてプラント１００の実施例を説明し、計測信号４に含まれる情報を述べる。次に、図４（ｂ）及び図４（ｃ）を用いて、計測信号４をカテゴリーに分類する様子を述べる。
【００８６】
図４（ａ）は、プラント１００の実施例である火力発電プラントを示すブロック図である。
【００８７】
図４（ａ）において、火力発電プラント１００は、ガスタービン発電機１１０、制御装置１２０及びデータ送信装置１３０を含む。ガスタービン発電機１１０は、発電機１１１、圧縮機１１２、燃焼器１１３及びタービン１１４を含む。
【００８８】
発電に際しては、圧縮機１１２にて吸い込んだ空気を圧縮して圧縮空気とし、この圧縮空気を燃焼器１１３に送り、燃料と混合して燃焼する。燃焼により発生した高圧ガスを用いてタービン１１４を回転させ、発電機１１１により発電を行う。
【００８９】
制御装置１２０においては、電力需要に応じてガスタービン発電機１１０の出力を制御する。また、制御装置１２０は、ガスタービン発電機１１０に設置されたセンサ（図示せず）で計測した運転データ１０２を入力データとしている。運転データ１０２は、吸気温度、燃料投入量、タービン排ガス温度、タービン回転数、発電機発電量、タービン軸振動などの状態量であり、サンプリング周期毎に計測している。また、大気温度などの気象情報も計測している。
【００９０】
制御装置１２０においては、これらの運転データ１０２を用いて、ガスタービン発電機１１０を制御するための制御信号１０１を算出する。また、制御装置１２０では、運転データ１０２の値が予め設定した範囲を逸脱した時に警報を発生させる処理を実施している。警報信号は、運転データ１０２が予め設定した範囲を逸脱した時に「１」、範囲内の時は「０」のデジタル信号として処理する。警報信号が「１」の時は、音や画面表示などで、警報の内容をオペレータに通知する。
【００９１】
信号データ送信装置１３０は、制御装置１２０で計測した運転データ１０２、及び制御装置１２０で算出した制御信号１０１、および警報信号を含む計測信号１を診断装置２００に送信する。
【００９２】
図４（ｂ）は、プラント１００から取得した計測信号１を、カテゴリーに分類した結果を説明する図である。横軸は、時間、縦軸は計測信号、カテゴリー番号である。図４（ｃ）は、プラント１００の計測信号１を、カテゴリーに分類した分類結果の一例を示す図である。
【００９３】
図４（ｃ）は、一例として、計測信号のうちの２項目を表示したものであり、２次元のグラフで表記した。また、縦軸及び横軸は、それぞれの項目の計測信号を規格化して示した。
【００９４】
計測信号は、図３（ａ）のＡＲＴモジュール７２０によって複数のカテゴリー１０００（図４（ｃ）に示す円）に分割される。１つの円が、１つのカテゴリーに相当する。
【００９５】
本実施例では、計測信号は４つのカテゴリーに分類されている。カテゴリー番号１は、項目Ａの値が大きく、項目Ｂの値が小さいグループ、カテゴリー番号２は、項目Ａ、項目Ｂの値が共に小さいグループ、カテゴリー番号３は項目Ａの値が小さく、項目Ｂの値が大きいグループ、カテゴリー番号４は項目Ａ、項目Ｂの値が共に大きいグループである。
【００９６】
図４（ｂ）に示すように、診断開始前の正常期間のデータは、カテゴリー１〜３に分類された。診断開始後の前半のデータはカテゴリー２に分類されており、正常期間と同じカテゴリーである。この場合、データの傾向が正常期間と同じであることから、正常であると診断する。一方、診断開始後の後半のデータはカテゴリー４に分類されており、正常期間とは異なるカテゴリーである。データの傾向が異なることから、プラントの状態が変化し、異常が発生した可能性がある。この場合、本発明の診断装置２００では、プラントのオペレータに異常発生の可能性があることを画像表示装置９４０に表示し、オペレータに通知する。
【００９７】
尚、本実施例においては、２項目の計測信号をカテゴリーに分類する例を述べたが、３項目以上の計測信号について多次元の座標を用いてカテゴリーに分類することもできる。
【００９８】
図５は、図４（ａ）に示した発電プラントにおける起動モードとプロセス値の関係を説明する図であり、出力指令値とタービン機器温度の経時変化を示している。
【００９９】
代表的な起動モードとして、ホットスタート、コールドスタートがある。タービンや圧縮機が熱い状態で再起動する場合をホットスタートと呼ぶ。また、比較的長期間停止し、タービン、圧縮機が冷えた状態で再起動する場合とコールドスタートと呼ぶ。図５（ｂ）に示すように、起動モードが異なると、起動開始時のタービン機器温度、負荷変化中の温度の変化範囲も異なる。
【０１００】
従来の診断装置では、運転条件に関わらず１つの診断モデルを構築していたため、正規化の範囲はデータを包含するように決定する必要がある。すなわち、例えば、正規化範囲の上限値を１０１０、下限値を１０１１としていた。
【０１０１】
例えば、ホットスタートの低負荷時では、プロセス値が変化する範囲１０２０に比べて、正規化範囲１０２２が広くなる。
【０１０２】
データの変化範囲に比べて正規化範囲が広いと、正規化後の値変化が小さくなる。そのため、異常発生時のデータの変化をとらえることができず、異常発生時にも新規カテゴリーが発生しない場合があった。これは、失報の原因となる。
【０１０３】
本発明に備えられているモデル定義手段４００を用いることで、運転条件に合わせて正規化範囲を適切に決定する。本機能により、失報を抑制して診断精度を向上できる。以下、その具体的な方法について説明する。
【０１０４】
図６は、モデル定義手段４００の構成要素である運転条件判定部５００の動作を説明するフローチャート図である。図６に示すように、本アルゴリズムはステップ５１０、５２０、５３０、５４０、５５０を組み合わせて実行する。
【０１０５】
まず、ステップ５１０では、計測信号４に蓄積されているデータを、一定負荷運転中、起動運転中、停止運転中、負荷変化運転中の期間毎に分割する。
【０１０６】
出力が変化していない時、すなわち、出力の計測信号の変化率が小さい場合は、一定負荷運転中とする。また、発電プラントの計測信号に含まれている、起動中、停止中、負荷変化中を区別する信号を基に起動運転中、停止運転中、負荷変化運転中の期間に分割する。
【０１０７】
一定負荷運転中の場合はステップ５２０に、起動運転中の場合はステップ５３０に、停止運転中の場合はステップ５４０に、負荷変化運転中の場合はステップ５５０にそれぞれ進む。
【０１０８】
ステップ５２０では、負荷帯に関する情報を抽出する。例えば、定格出力の０〜５０％を低出力、５０％〜１００％を高出力など、出力に応じて負荷帯を分類する。ステップ５３０では、起動モードの種類、起動中の負荷変化率に関する情報を抽出する。起動モードとしては、ホットスタートモード、コールドスタートモードなどがある。ステップ５４０では、停止モード、停止動作開始時出力に関する情報を抽出する。停止モードとしては、通常運用の中でプラントを停止するモード、異常発生時に緊急遮断するモードなどがある。ステップ５５０では、負荷変化率、負荷変化開始時出力、負荷変化終了時出力に関する情報を抽出する。
【０１０９】
尚、本実施例ではステップ５２０、５３０、５４０、５５０で、一定負荷運転中、起動運転中、停止運転中、負荷変化運転中の状態を区別するため、先に述べた情報を抽出しているが、この情報量を増やすことも可能である。例えば、タービン回転数、昇速率、プラントに供給する燃料の種類、大気温度など、プラント１００の計測信号を処理して得られる情報であれば、ステップ５２０、５３０、５４０、５５０で抽出する情報に追加してもよい。
【０１１０】
図６（ｂ）に示すように、同じ計測項目を用いた診断モデル７１０は、図６（ａ）で抽出した情報に基づいてサブ診断モデル７２０に分類したモデルの集合で構築する。サブ診断モデル毎に異なる正規化条件を定義して診断する。
【０１１１】
以下、図７〜図９を用いて、正規化条件を決定する正規化条件決定部６００の実施例を説明する。正規化条件決定部６００で決定した正規化条件の情報を用いて、データ前処理装置７１０で計測信号を正規化する。計測信号ｘｉのデータ項目数がＮ個でｎ番目の計測信号をｘ（ｎ）とする。正規化したデータＮxi（ｎ）は、下記式（１３）で表される。ただし、Ｎmin（ｎ）は正規化の下限値、Ｎmax（ｎ）は正規化の上限値である、
【０１１２】
【数１３】
0005337853.tif 000014
【０１１３】
正規化条件決定部６００では、（１３）式で示したＮmin（ｎ）、Ｎmax（ｎ）を決定する。
【０１１４】
図７は、正規化条件決定部６００の第１の実施例の説明図である。
【０１１５】
図７（ａ）は、正規化条件決定部６００の第１の実施例のフローチャート図である。図７（ａ）に示すように、本アルゴリズムはステップ６１１、６１２、６１３を組み合わせて実行する。
【０１１６】
まず、ステップ６１１では、運転条件毎にモデル構築用データを分割する。ステップ６１２では、運転条件毎に計測信号変化幅の情報を抽出する。
【０１１７】
ステップ６１３では、正規化範囲の上限値Ｎmax１（ｎ）、正規化範囲の下限値Ｎmin１（ｎ）を決定する。
【０１１８】
一定負荷運転中の場合は、（１４）（１５）式を用いる。ここで、Ｄmax１（ｎ）は計測信号の最大値、Ｄmin１（ｎ）は最小値であり、α、βは定数である。
【０１１９】
【数１４】
0005337853.tif 000015
【０１２０】
【数１５】
0005337853.tif 000016
【０１２１】
負荷変化中、起動中、停止中の場合は、出力指令値に合わせて、式（１６）（１７）のように正規化範囲を決定する。
【０１２２】
【数１６】
0005337853.tif 000017
【０１２３】
【数１７】
0005337853.tif 000018
【０１２４】
ここで、ＭＷは出力指令値であり、ａ、ｃは、負荷変化中の計測信号を１次近似した時の傾きである。また、ｂ、ｄは、式（１８）（１９）で計算する値である。
【０１２５】
【数１８】
0005337853.tif 000019
【０１２６】
【数１９】
0005337853.tif 000020
【０１２７】
ここで、ｆは負荷変化中の計測信号を１次近似した時の切片であり、Ｄmax２は計測信号を１次近似した直線と計測信号との偏差の最大値、Ｄmin２は偏差の最小値である。
【０１２８】
図７（ｂ）は、低出力から高出力に負荷変化した時の、出力と計測信号の経時変化を説明する図である。時刻Ｔ０ａまでは低出力、時刻Ｔ０ａから時刻Ｔ０ｂまでは負荷変化、時刻Ｔ０ｂ以降は高出力である。
【０１２９】
図７（ａ）に示したフローチャートを動作させることにより、低出力時の正規化範囲の上限値は１０３０、下限値は１０４０、負荷変化中の正規化範囲の上限値は１０５０、下限値は１０６０、高出力時の正規化範囲の上限値は１０７０、下限値は１０８０に決定される。このように、本発明の正規化範囲は、運転条件に合わせて変化する。
【０１３０】
図８は、正規化条件決定部６００の第２の実施例の説明図である。図８（ａ）に示すように、本アルゴリズムはステップ６３１、６３２、６３３を組み合わせて実行する。
【０１３１】
ステップ６３１では、負荷変化モードのデータを抽出する。ステップ６３２では、データ項目毎に、無駄時間を推定する。プロセス信号を入出力データとしてインパルス応答によりむだ時間を推定する。具体的な計算方法として参考文献「制御のための上級システム同定」（東京電機大学出版局）に記載されている方法が挙げられる。ステップ６３３では、ステップ６３２で推定した無駄時間分データをシフトした後、データ項目毎に正規化範囲を決定する。
【０１３２】
この様子を、図８（ｂ）を用いて説明する。出力指令値が増加し始める時刻Ｔ１ａより、計測信号が変化し始める時刻が遅くなる。この時間が無駄時間である。データ項目Ａの無駄時間は（Ｔ２ａ−Ｔ１ａ）、データ項目Ｂの無駄時間は（Ｔ３ａ−Ｔ１ａ）である。
【０１３３】
図８（ａ）のフローチャートを用いて計測信号を無駄時間分ずらした後、図７のフローチャートを用いて正規化範囲を決定する。これにより、出力指令値に合わせて計測信号の正規化範囲を決定できる。
【０１３４】
図９は、図１に示した発電プラントの診断装置における正規化条件決定部６００の第３の実施例の説明図である。
【０１３５】
計測信号の最大値、最小値付近での正規化範囲を広くして、最大値、最小値近傍での計測信号の変化を検知しやすくする。
【０１３６】
尚、正規化条件決定部６００の正規化方法は上述した内容に限定されるものではなく、運転条件毎に正規化条件を決定するものであれば良い。たとえば、プラントの設計情報や計測器の仕様からデータ値が変化する範囲を推定し、その範囲を正規化範囲としてもよい。
【０１３７】
図１０は、図１に示した発電プラントの診断装置２００の動作モードを説明するフローチャート図である。図１０（ａ）は図２におけるモデル構築モードの動作フローチャートであり、図１０（ｂ）は診断モードの動作フローチャートである。
【０１３８】
図１０（ａ）に示すように、モデル構築モードは、ステップ１２００、１２１０、１２２０、１２３０を組み合わせて実行する。ステップ１２００では、計測信号データベース３１０からモデル構築に用いる期間のデータを抽出する。この期間は、プラント１００のオペレータにより任意に設定できる。次に、ステップ１２１０では運転条件判定部５００を動作させる。図６（ａ）に示したフローチャートが動作し、モデル構築に用いる期間のデータを、一定負荷運転中、起動運転中、停止運転中、負荷変化運転中の運転モード毎に分割し、さらにステップ５２０、５３０、５４０、５５０を動作させて各モードのデータを特徴ごとに分割する。次に、ステップ１２２０では正規化条件決定部６００を動作させる。図７（ａ）、図８（ａ）で説明したフローチャートを動作させ、ステップ１２１０で分割したグループ毎に正規化条件を決定する。ステップ１２１０、１２２０を動作させて得られたモデル定義情報５は、モデル定義データベース３２０に保存する。最後に、ステップ１２３０ではモデル構築手段７００を動作させる。ステップ１２１０で分割した各グループは、サブ診断モデル７２０（図６（ｂ）参照）として定義され、サブ診断モデル毎に定義されている正規化条件を用いて計測信号を処理し、図３で述べたＡＲＴを用いて診断モデルを構築する。
【０１３９】
図１０（ｂ）に示すように、診断モードは、ステップ１３００、１３１０、１３２０を組み合わせて実行する。ステップ１３００では、計測信号データベース３１０から診断する期間のデータを抽出する。次に、ステップ１３１０では、診断期間のデータを運転条件判定部５００で処理する。モデル構築モードで構築した複数のサブ診断モデルから、運転条件が一致するサブ診断モデルを抽出する。最後に、ステップ１３２０では診断手段８００を動作させる。診断手段８００では、モデル定義データベース３２０からステップ１３３０で抽出したサブ診断モデル情報を抽出する。サブ診断モデルのカテゴリー番号と、診断期間のデータをＡＲＴで分類して得られたカテゴリー番号を比較する。モデル構築モードを動作させた時と同じカテゴリーの時は正常と診断し、異なるカテゴリー番号が発生していた時は異常と診断する。診断結果１０を外部出力インターフェイス２２０に出力する。
【０１４０】
図１１は、本発明のデータベースに保存されているデータの態様を説明する図である。
【０１４１】
図１１（ａ）に示すように、計測信号データベース３１０には、プラント１００に対して計測した運転データである計測信号１（図では、データ項目Ａ、Ｂ、Ｃを記載）の値が、サンプリング周期（縦軸の時刻）毎に保存される。
【０１４２】
表示画面３１１において縦横に移動可能なスクロールボックス３１２及び３１３を用いることにより、広範囲のデータをスクロール表示することができる。
【０１４３】
モデル定義データベース３２０には、図１１（ｂ）に示すように、運転条件と正規化範囲の情報が対応付けられて保存される。
【０１４４】
診断モデルデータベースには、図１１（ｃ）に示すように、カテゴリー番号と重み係数の関係が保存されている。ここで、重み係数とは、カテゴリーの中心座標のことである。
【０１４５】
図１２は、図１に示した発電プラントの診断装置の適用効果を説明する図である。
【０１４６】
時刻Ｔ４で負荷変化を開始し、時刻Ｔ５で負荷変化を終了した。負荷変化終了後、時刻Ｔ６で異常が発生した。
【０１４７】
図１２は、出力指令値と、タービン温度の関係をカテゴリー番号と、温度の計測信号の関係も合わせて図示してある。異常発生に伴い温度が上昇したが、カテゴリー４の範囲内であるため正常状態と同じカテゴリーに分類された。従来方式では、正規化範囲が広く、異常発生に伴う温度上昇をとらえることができず、異常を検知できなかった。
【０１４８】
本発明では、運転条件に合わせて診断モデルが切り替わる。本実施例では、時刻Ｔ４までは運転条件１診断モデル（一定負荷運転中、低出力モデル）、時刻Ｔ４〜Ｔ５の間は運転条件２診断モデル（負荷変化運転中）、時刻Ｔ５以降は運転条件３診断モデル（一定負荷運転中、高出力モデル）で診断するように切り替わる。
【０１４９】
運転条件３診断モデルにおいて、正常状態はカテゴリー番号１〜５に分類されている。異常発生前は正常状態と同じカテゴリーに分類されているが、異常発生後のデータはカテゴリー番号６の新規カテゴリーに分類され、異常を検知できた。すなわち、従来方式よりも状態を細かく分類できるので、異常を検知できない失報を抑制できた。
【０１５０】
図１３は、図１に示した発電プラントの診断装置２００における画像表示装置９４０に表示される画面を説明する図である。
【０１５１】
カーソル９５１をマウス９３０で操作して、正規化範囲の上限値（９５４、９５６、９５８）、正規化範囲の下限値（９５５、９５７、９５９）、運転条件の境界時刻（９５２、９５３）を任意に調整できる。調整した結果をモデルに反映するには、図１３の画面上で実行ボタン９６０をクリックする。この操作により、図１０（ｂ）に示したモデル定義データベースの運転条件と正規化範囲の情報が変更され、調整結果をモデル構築、及び診断動作に反映できる。
【０１５２】
なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明をわかりやすく説明するために詳細に記載したものであり、必ずしも説明したすべての構成を備えるものに限定されるものではない。
【０１５３】
また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を集積回路で設計する等によりハードウェアで実現してもよい。また上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行するソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル、計測信号、算出情報等の情報は、メモリやハードディスク等の記憶装置、またはＩＣカード、ＳＤカード、ＤＶＤ等の記憶媒体に置くことができる。よって、各処理、各構成は処理ユニットやプログラムモジュールとして実現可能である。
【０１５４】
また、情報線は説明上必要と考えられるものを示しており、製品上必ずしもすべての制御線や情報線を示しているとは限らない。実際には、ほとんど全ての構成が相互に接続されていると考えてもよい。
【０１５５】
本実施例によれば、発電プラントの異常を高精度で検知する発電プラントの診断装置及びプラントの診断方法が得られる。
【産業上の利用可能性】
【０１５６】
本発明は、プラントの診断装置及びプラントの診断方法として各種プラント等に広く適用可能である。
【符号の説明】
【０１５７】
１、３、４  計測信号
２  外部入力信号
５、６、７  モデル定義情報
８、９  モデル情報
１０  診断結果
１１  画面表示情報
５０  診断装置情報
１００  プラント
２００  診断装置
２１０  外部入力インターフェイス
２２０  外部出力インターフェイス
３１０  計測信号データベース
３２０  モデル定義データベース
３３０  診断モデルデータベース
４００  モデル定義手段
５００  運転条件判定部
６００  正規化条件決定部
７００  モデル構築手段
８００  診断手段
９００  運転管理室
９１０  外部入力装置
９２０  キーボード
９３０  マウス
９４０  画像表示装置
【図１】
0005337853.tif 000021
【図２】
0005337853.tif 000022
【図３】
0005337853.tif 000023
【図４】
0005337853.tif 000024
【図５】
0005337853.tif 000025
【図６】
0005337853.tif 000026
【図７】
0005337853.tif 000027
【図８】
0005337853.tif 000028
【図９】
0005337853.tif 000029
【図１０】
0005337853.tif 000030
【図１１】
0005337853.tif 000031
【図１２】
0005337853.tif 000032
【図１３】
0005337853.tif 000033