(57)【特許請求の範囲】
【請求項１】
発電設備に設けたセンサによるセンサ情報と該設備の運転に関する運転情報及びメンテナンス情報に基づいて、該発電設備の機器の損傷を有限要素法の数値解析モデルにより評価する診断装置において、
前記センサ情報と前記運転情報とメンテナンス情報を前記数値解析モデルの境界条件に変換する境界条件算出手段と、算出された境界条件を与えられる数値解析モデルにより損傷の数値解析を実施する数値解析手段と、
かつ、実験計画法を用いて、前記センサ情報と前記運転情報とメンテナンス情報から算出された境界条件を変数として直交表に割付け、この直交表に従って前記数値解析モデルにより運転条件を模擬した熱解析及び構造解析を行い、その解析結果と実機の検査データを比較し、その差が最小となるように前記変数である前記境界条件を決定する逆問題解析手段、または、前記差を最小化する最適化手法を用いて、前記境界条件を決定する逆問題解析手段とを設け、
前記逆問題解析手段により求められた変数である前記境界条件により同定された数値解析モデルを前記数値解析手段に再設定した後、前記境界条件算出手段を用いてセンサ情報と運転情報とメンテナンス情報を境界条件に変換し、前記同定された損傷の数値解析を実施し、損傷結果を出力することを特徴とする発電設備の損傷診断装置。
【請求項２】
請求項１において、前記数値解析手段が予め決められた設備の起動停止回数毎の損傷を算出すると共に、算出した損傷を累計し、損傷の累計がある基準値を超えたときに損傷に関する情報を出力する損傷評価手段を設けることを特徴とする発電設備の損傷診断装置。
【請求項３】
発電設備に設けたセンサによるセンサ情報と該設備の運転に関する運転情報及びメンテナンス情報に基づいて、該発電設備の機器の損傷を有限要素法の数値解析モデルにより評価する診断装置において、
実験計画法により、前記センサ情報と前記運転情報とメンテナンス情報を前記数値解析モデルの境界条件の変数として直交表に割付ける直交表作成手段と、この直交表に従い、前記数値解析モデルにより運転条件を模擬した熱解析及び構造解析を行い、その解析結果の物理量を用いた損傷解析を行う数値解析手段と、前記数値解析手段で得られた損傷解析結果について各変数の効果を算出して有意な変数を、分散分析を用い、予め設定された有意水準を基準として、選択する分散分析手段と、前記有意な変数に基いて損傷に係わる評価項目の近似式を作成する近似式作成手段とを設け、該近似式に前記センサ情報と前記運転情報を入力して損傷診断を行うことを特徴とする発電設備の損傷診断装置。
【請求項４】
請求項１〜３のいずれかにおいて、前記センサ情報と前記運転情報を収集する計算機を設け、該計算機と前記診断装置を通信回線を介して結び、遠隔診断を可能にしたことを特徴とする発電設備の損傷診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は発電設備の機器の損傷を診断する装置に関する。
【０００２】
【従来の技術】
発電設備の機器の損傷診断は、ある定期的な時間間隔や運転回数間隔で機器を停止させ、部品の損傷状態を点検し、損傷状態がある基準値を越えたものについて必要な補修や交換を行う方法が一般的である。また、発電設備の機器に備え付けられた各種のセンサの出力を監視し、出力がある基準値を越えた場合に機器を停止し、損傷状態を点検する場合もある。
【０００３】
近年の計算機とインターネットに代表されるデータ通信の進歩にともない、機器の損傷診断において、機器の運転状態の機器情報を通信回線を介して遠隔監視施設に送信し、劣化を診断する装置が知られている（特開平11−3113号公報）。
【０００４】
【発明が解決しようとする課題】
発電設備はその運転履歴、運転条件、使用環境が異なれば、機器の損傷度合が変化する。例えば、ガスタービン発電設備などでは、近年の電力需要のピークシェイビングのための負荷変動や起動停止回数の増加などがある。負荷変動や起動停止回数の増加は、一般的には部品の熱疲労損傷を増加させる。また、発電設備の立地により、吸い込み空気温度や空気中に含まれる各種元素が異なり、部品の温度条件、荷重条件、環境条件が変化する。外気温が低くなれば、タービン翼に負荷される流体荷重は大きくなり、流体荷重による損傷は増加する。また、発電設備が臨海工業地帯であれば、海水に含まれるナトリウムやカリウム、大気中の煤煙に含まれる亜硫酸ガスや窒素酸化物などの腐食作用のある物質が部品に付着し、疲労損傷が増加する。
【０００５】
このように、運転履歴、運転条件、使用環境の変化を考慮した場合に、前述の損傷診断方法のように定期間隔の点検では信頼性が万全でない場合が想定される。逆に、立地条件の良い設備では、過度に余裕のある保守点検となる可能性がある。しかしながら、運転履歴、運転条件、使用環境の変化を運転員と保守員が記録し、損傷を診断する方法は極めて効率が悪い。
【０００６】
上記の特開平11−3113号公報などのように、通信回線を用いた機器の劣化診断装置を用いれば、運転履歴、運転条件、使用環境の変化を常時監視することが可能ではある。しかし、本例は損傷の数値解析モデルを有していないため、損傷を正確に診断できないことがある。
【０００７】
本発明の目的は、上記した従来技術の問題点に鑑み、運転履歴、運転条件、使用環境の変化をオンラインで監視し、その情報を機器の損傷の数値解析モデルに与え、オンラインで機器の損傷状態を監視し、機器の損傷を正確かつ迅速に算出し、機器の信頼性を確保するものである。
【０００８】
【課題を解決するための手段】
上記目的を達成するために、本発明は、発電設備に設けたセンサによるセンサ情報と該設備の運転に関する運転情報及びメンテナンス情報に基づいて、該発電設備の機器の損傷を数値解析モデルにより評価する診断装置において、前記センサ情報を前記数値解析モデルの境界条件に変換する境界条件算出手段を設けている。
【０００９】
これによれば、発電設備のセンサ情報を境界条件と該センサの関係式を用いて、損傷の数値解析モデルの境界条件に変換し、この境界条件を用いて数値解析モデルによる機器の損傷解析を実施するので、機器の損傷状態を正確に予測することができ、この結果をもとに機器の転換や交換および運用の仕方を決定できる。
【００１０】
また、本発明は、実験計画法を用いて、前記数値解析モデルの機器の損傷状態を評価する逆問題解析手段を設け、該手段による数値と実機器の損傷状態を評価する数値を比較し、前記損傷状態を評価する数値を統合する。すなわち、両者の差が最小となる境界条件を求め、該境界条件を設定した数値解析モデルにより損傷の数値解析を実施する。
【００１１】
これによれば、実際の機器の損傷と比較し、数値解析モデルを実機の損傷に一致するように境界条件を改良することができるので、機器の損傷状態をより正確に予測することができ、機器の信頼性を高めることができる。
【００１２】
また、本発明は、実験計画法を用いて、前記センサ情報と前記運転情報とメンテナンス情報を直交表に割付け、損傷の数値解析モデルを解析し、分散分析を行って、損傷に係わる評価項目の近似式を作成し、該近似式を用いて損傷診断を行う。これによれば，該近似式により損傷解析を行うことができるので，有限要素法のような数値モデルの解析を実施する場合に比べ、計算機の計算負荷を軽減し、迅速に損傷を評価することができる。
【００１３】
さらに、本発明は、前記センサ情報と前記運転情報を収集する計算機を設け、該計算機と上記の診断装置を通信回線を介して結ぶので、遠隔地での監視を行うことができる。
【００１４】
【発明の実施の形態】
以下、本発明の実施の形態について、図面を参照しながら説明する。図１は発電設備の損傷遠隔診断装置の一実施例を示す。発電設備は、圧縮機１と燃焼器２とタービン３から構成されるガスタービンと、該ガスタービンにより駆動される発電機４と、排熱回収ボイラ５と蒸気タービン６と蒸気タービン６により駆動される発電機７および復水器８から構成される。これら主要機器、すなわち、圧縮機１、燃焼器２、タービン３、発電機４，７、ボイラ５、蒸気タービン６、復水器８には、その運転状態を監視、制御するために、各種のセンサ９が取り付けられている。
【００１５】
これらのセンサ９はケーブルにより制御装置１０に接続され、センサ情報を制御装置１０に送信している。制御装置10は運転監視計算機11にケーブルを介して接続し、センサ情報と運転情報を運転監視計算機１１に送信する。運転監視計算機１１は発電設備側ファイヤーウォール１２、インターネット１３、遠隔監視施設側ファイヤーウォール１４を介して、遠隔損傷診断装置１５にセンサ情報と運転情報を送信する。メンテナンス情報は、遠隔診断装置１５にオンラインまたはオフラインで入力される。なお、制御装置１０と運転監視計算機１１は一体化されていてもよい。
【００１６】
図１には，遠隔損傷診断装置１５におけるデータの処理機能についても示している。送信されたセンサ情報と運転情報は入出力インターフェイス15aを介して、境界条件算出ステップ15bに送られ、センサ情報を数値解析モデルの境界条件に変換する。
【００１７】
ここで述べる境界条件とは、有限要素法の熱解析では、流体温度や熱伝達率などの熱条件であり、構造解析では、荷重、圧力、遠心力、加速度などの荷重条件である。その他の数値解析モデルにおいても、流体温度や熱伝達率などの熱条件と、荷重、圧力、遠心力、加速度などの荷重条件を指すものとする。また、境界条件は、数値解析モデルや損傷解析に変数として入力することのできる条件を指すものとし、詳細については後述する。
【００１８】
この境界条件は損傷の数値解析手段１５ｃに送られ、機器の損傷が算出される。数値解析手段１５ｃは、有限要素法解析手段１５ｃ１と損傷解析数値モデル１５ｃ２からなり、前者は境界条件が与えられると部品の温度，変形，応力，ひずみ等の物理量を算出する。後者は算出された部品の物理量から損傷量を算出する。損傷量はある決められた回数毎、たとえば、一起動停止あたりの損傷量を算出する。起動停止の判別は、タービンの回転数、燃料、発電量などの情報から容易に判断できる。
【００１９】
算出された損傷は、損傷評価手段１５ｄへ送られ、過去の運転の累積の損傷と加算する。累積の損傷があらかじめ定められた基準値を越えた場合には、入出力インターフェイス１５ａに警告を表示する。そして、電子メールなどの手段を用いて、運転監視計算機１１に損傷情報を送信する。なお、各手段において利用する変数は適宜データベース１５eより読み出され、各手段での結果も適宜、データベース１５ｅに書き込まれる。
【００２０】
図２に本実施例の処理フローを示す。ステップ１において、発電設備のセンサ情報を入手する。ステップ２において、発電設備の運転情報を入手する。センサ情報と運転情報の例は図３に示す。
【００２１】
ステップ３において、入力されたセンサ情報と運転情報から、有限要素法熱解析の境界条件である流体温度と熱伝達率を作成する。作成方法の詳細は、数３〜数１０の説明において述べる。ステップ４において、入力されたセンサ情報と運転情報から、有限要素法構造解析の境界条件である荷重条件を作成する。作成方法の詳細は、数１〜数２の説明において述べる。
【００２２】
ステップ５において、上記の境界条件を入力し、ステップ６において、有限要素法による熱解析と構造解析を行う。ステップ７において、有限要素法による熱解析と構造解析の解析結果として、損傷解析に用いる物理量を出力する。物理量としては、部材の温度、応力、ひずみ等がある。
【００２３】
ステップ８において、求められた物理量を損傷解析モデルに入力し、ステップ９において、損傷解析を実施する。損傷解析の詳細は、数１１の説明にて後述する。ステップ１０において、損傷結果を出力し、データベース１５ｅに書き込む。損傷解析結果としては疲労損傷率や、クリープ損傷率や、き裂長さ変形量など、部材の損傷を数値化したものである。
【００２４】
ステップ１１において、部材の損傷結果の累計をデータベース１５ｅより読み込み、ステップ１０において算出した損傷解析結果を積算し、データベース１５ｅに書き込む。ステップ１２において、データベースより損傷解析結果の積算、つまり、損傷解析結果の累計を読み込み、基準値と比較する。さらに、ステップ13において、運転の予定、たとえば、運転する時間や運転する回数や運転の形態を損傷解析に入力し、損傷を外挿する。このとき、有限要素法解析モデルに与える境界条件は、予め与えられる設計条件としたり、過去のセンサ情報と運転情報及びメンテナンス情報から算出する。ステップ１４において、損傷解析結果の出力を行う。
【００２５】
また、ステップ１２で、累積の損傷があらかじめ定められた基準値を越えた場合は、損傷が基準値を越えた警告を入出力インターフェイス15aに表示する。電子メールなどの手段を用いて、遠隔監視施設側ファイヤーウォール１４、インターネット１３、発電設備側ファイヤーウォール１２を経由して、運転監視計算機１１に損傷情報を送信する。損傷が基準を越えていなくても、定期的に機器の損傷情報を送信しても良く、損傷が基準を超える時期を予測して送信しても良い。
【００２６】
図２の例では、有限要素法解析モデルを用いて、その境界条件を、センサ情報と運転情報から算出する例を示した。本発明は、数値解析モデルが有限要素法解析モデルでない場合も含む。つまり、センサ情報データと運転情報データとメンテナンス情報データを数値解析モデルや損傷解析の変数として入力し、損傷を解析する場合も含むものである。
【００２７】
図３に、遠隔損傷診断装置１５の入出力インターフェイス15aに入力されるセンサ情報、運転情報及びメンテナンス情報を、ガスタービンを例として示す。ガスタービンには、排ガス温度ａ１、タービンのホイールスペース温度ａ２、圧縮機の吐出空気温度ａ３、吐出空気圧力ａ４、燃焼器の保炎器の温度ａ５、燃焼器の燃料流量ａ６、入口空気温度ａ７、入口空気圧力ａ８、入口空気湿度ａ９、入口可変翼開度ａ10、回転数ａ11、軸受振動ａ12、軸振動ａ13、軸受メタル温度ａ14等を検出するセンサ９が設けられている。また、試運転時や特殊計測時には、圧縮機各段落の圧力ａ15、空気温度ａ16、圧力変動ａ17、吸い込み空気のＮａセンサａ18、燃焼器の燃焼ガス圧力変動、温度変動、圧縮機翼ひずみａ19、圧縮機翼温度ａ20、タービン翼ひずみa21、タービン翼温度ａ22、燃焼器ひずみａ23、燃焼器温度ａ24、ケーシング温度a25、ケーシングひずみa26、ケーシングa27、ケーシング変位a28などを測定するセンサを設けることがあり，これらのセンサ情報を用いてもよい。
【００２８】
運転情報としては、起動停止回数ｂ１、燃焼時間ｂ２、トリップ回数（緊急停止回数）ｂ３、負荷変動の回数ｂ４、発電出力ｂ５、発電効率ｂ６、圧縮機効率b7等がある。
【００２９】
メンテナンス情報としては、定期検査時における損傷状況を数値化したデータで、例えば、損傷の寸法や面積や体積や重量およびその数、つまり、き裂長さｃ１や発生数ｃ２、減肉の重量ｃ３や体積ｃ４や面積c5とその数ｃ６、コーティングの剥離についてもその面積ｃ７とその数ｃ８、腐食についても面積ｃ９とその数ｃ10などがある。さらに、定期検査日c11、定期検査回数c12、圧縮機洗浄回数c13、圧縮機水洗浄日c14、圧縮機水洗浄水ｐＨc15も、メンテナンス情報として利用できる。
【００３０】
メンテナンス情報は入出力インタフェース15aを通じてデータベース15eに入力し、適宜利用する。また、発電設備に発電機器が複数台存在する場合は、その認証番号を送信する。なお、図３では各センサの出力値の変数名を便宜的に付与している。
【００３１】
入出力インターフェイス15aに入力された情報は、境界条件算出手段15bに渡され、数値解析のための境界条件が算出され、この境界条件が損傷の数値解析手段15cに渡される。例としてガスタービン動翼を説明する。
【００３２】
図４と図５にガスタービン動翼の有限要素法解析モデル15c1の全体図とその断面図を示す。ガスタービン動翼１６の有限要素法解析モデルは、図４（ａ）に示すように、動翼の形状を有限要素法解析節点と要素によりモデル化され、前縁16aや翼先端16ｂやプラッフォーム16ｃやシャンク16dなどを有している。ガスタービン高温部品は一般に冷却空気で冷却されているが、動翼も冷却されている。このため、有限要素法解析モデルも冷却孔17a〜17lに示す冷却孔がモデル化されている。
【００３３】
有限要素法解析モデル１５ｃ１に解析の境界条件を与えて、部品の運転状態の各種の物理量を算出する。機械荷重の境界条件としては、回転数に応じた遠心力や運転条件に応じた流体力を負荷する。熱荷重の境界条件としては、一般的に熱伝達率やガス温度を有限要素法解析モデルの熱交換のある部位、つまり、流体に曝される部位に与える。
【００３４】
図５を用いて、熱境界条件の与え方の一例を示す。図５には、図４に示した動翼の有限要素法解析モデルの翼の断面図を示す。熱境界は、翼の燃焼ガスに曝される面18a〜cと冷却ガスに曝される面19a〜lを区別して与える。場合によっては、翼の燃焼ガスに曝される面18a〜cと冷却ガスに曝される面19a〜lの個々の領域において分布を持たせたり、さらに個々の領域を分割して熱境界条件を有限要素法解析モデル１５ｃ１に与える。
【００３５】
境界条件算出手段15bにおける境界条件の算出方法を説明する。機械的な荷重条件として、回転部品に加わる遠心力Ｆ１は数１で与えられる。ここで、a11は図３で示した回転数であり、d11は定数である。
【００３６】
【数１】
0003788901.tif 000002
【００３７】
また、機械的な荷重条件として、流体力Ｆ２はいくつかのセンサ情報を変数とする数２で与えられる。
【００３８】
【数２】
0003788901.tif 000003
【００３９】
熱的な境界条件として、燃焼ガス側の熱伝達率Hgiはセンサ情報のいくつかを変数とする数３で与えられる。なお、ｉは分割した境界の番号である。
【００４０】
【数３】
0003788901.tif 000004
【００４１】
燃焼ガス側のガス温度Tgiはセンサ情報のいくつかを変数とする数４で与えられる。
【００４２】
【数４】
0003788901.tif 000005
【００４３】
冷却空気側の熱伝達率Hciはセンサ情報のいくつかを変数とする数５で与えられる。ｉは分割した境界の番号である。
【００４４】
【数５】
0003788901.tif 000006
【００４５】
冷却空気側のガス温度Tciはセンサ情報のいくつかを変数とする数６で与えられる。
【００４６】
【数６】
0003788901.tif 000007
【００４７】
さらに、数値流体解析を用いて得られた熱的境界条件を変数として入力して、熱境界条件を算出してもよい。数値流体解析により得られた燃焼ガスに曝される部位iの燃焼ガスの熱伝達率とガス温度をそれぞれHgci，Tgci、冷却空気に曝される部位ｊの冷却空気の熱伝達率とガス温度をそれぞれHｃcｊ，Tｃcｊとする。このときの燃焼ガスの熱伝達率Hgiとガス温度Tgiはセンサ情報のいくつかを変数とする数７、数８で与えられる。
【００４８】
【数７】
0003788901.tif 000008
【００４９】
【数８】
0003788901.tif 000009
【００５０】
冷却空気の熱伝達率Hcjとガス温度Tcjはセンサ情報のいくつかを変数とする数９、数１０で与えられる。
【００５１】
【数９】
0003788901.tif 000010
【００５２】
【数１０】
0003788901.tif 000011
【００５３】
上記した境界条件を示す数１〜数１０においては、適宜、図２に示した他のセンサ情報や運転情報あるいはメンテナンス情報を変数として加えることも可能である。
【００５４】
このように、遠隔損傷診断装置１５の境界条件算出手段15bにおいて、入出力インターフェイス15aから入力されたセンサ情報などは、数１〜数１０により境界条件に変換され、数値解析手段15cに渡される。
【００５５】
数値解析手段１５ｃでは、図４と図５に示す有限要素法解析が実施される。実施される解析には、定常温度解析や非定常温度解析などの熱解析と、静荷重解析や熱応力解析やクリープ解析や振動解析などの応力解析がある。数値解析ステップ１５ｃの有限要素法解析モデル15c1により、部品の損傷を評価するための部品の物理量、例えば、温度、変形、応力、ひずみ等が算出される。
【００５６】
図６は物理量の解析例で、動翼の運転時における応力ひずみと温度の履歴を示す。応力とひずみの変動範囲や部材の最高温度を用いて、損傷解析を行う。図に示すような温度，応力，ひずみのグラフを入出力インターフエース１５ａに出力してもよい。
【００５７】
損傷解析モデル15c2は得られた物理量を入力して損傷解析を実施する。疲労損傷解析としてはマイナー則、修正マイナー則などに代表される線形損傷則などがある。クリープ損傷解析としては寿命比則などが一般的である。また、き裂による損傷に関しては、破壊力学解析を損傷解析として用いてもよい。
【００５８】
損傷解析モデルは実機の損傷を統計処理し、寿命を支配する適切なモデルを選択する。部材の損傷Ｄは一般に数１１にて表わされる。
【００５９】
【数１１】
0003788901.tif 000012
【００６０】
ここで、ｉは部材番号、ｋは損傷の形態、σは応力、εはひずみ、Ｔは温度、Ｎは負荷回数、ｔは負荷時間である。損傷解析モデルは適宜、部材損傷を良く模擬できるモデルを選択し、実機の損傷の回帰モデルなどを用いてもよい。
【００６１】
図７〜図９に損傷解析の一例を示す。図７はガスタービン静翼の有限要素法構造解析結果である。解析結果として、（ａ）のような応力分布図（または、温度分布図）や、（ｂ）のような温度応力の関係が得られる。センサ情報等を用いて、境界条件を変更した解析結果を出力してもよい。
【００６２】
この部品の場合は、熱疲労によるき裂進展が部品の寿命を支配する損傷と考えられる。このため、得られた応力の変動範囲を用いて破壊力学解析を行う。破壊力学解析のモデルを図８に示す。数１２を起動停止回数Ｎにより積分し、き裂深さaを求め、き裂深さaがある基準に達したときに警報を出力する。
【００６３】
【数１２】
0003788901.tif 000013
【００６４】
ここで、 ΔJはＪ積分範囲であり図６に示した結果から算出される物理量、C0、mは材料定数を示す。C1はクリープの効果を表現する補正係数、C2は材料劣化を表現する補正係数である。運転モードの違いによる圧縮保持クリープの影響、補修による材料の劣化を係数Ｃ１、Ｃ２として考慮している。
【００６５】
図９は、き裂進展による損傷解析結果を示す。横軸は起動停止回数をある値で正規化したものである。損傷として、き裂深さaを損傷解析で算出しており、き裂深さaがある寸法に達したときに警告を出力する。
【００６６】
損傷解析として一般的なものとしては、前述した線形損傷則がある。これは、数１３に示した疲労損傷率Ｄｆと、数１４に示したクリープ損傷率Ｄｃがある。
【００６７】
【数１３】
0003788901.tif 000014
【００６８】
ここで、ｎｉは応力σiの繰返し数、Nｉは応力σiにおける寿命である。応力は数値解析より得られる。疲労損傷率Dｆを発電設備の起動停止毎に算出し，過去の疲労損傷率Dｆと合計し、Ｄｆの累積を算出する。
【００６９】
【数１４】
0003788901.tif 000015
【００７０】
ここで、ｔｉは温度Tmpi，応力σiの保持時間であり、Ｔｉは温度Tmpi，応力σiにおけるクリープ破断寿命である。温度Tmpi，応力σiは数値解析により得られ、ｔｉは運転情報データをもとに算出する。算出したＤｆとＤｃの和が、ある基準値を越えたときに警報を発する。
【００７１】
次に、算出された損傷率は、損傷評価手段１５ｄに送られる。損傷評価手段１５ｄでは、過去の累積損傷率をデータベースから読み出し、損傷解析手段１５ｃで算出した損傷率を適宜加算する。そして、損傷率の総和がある基準値を越えた場合、部品の損傷率に関する情報を入出力インターフェイス15aに表示する。また、損傷率に関する情報を通信回線を介して、発電設備の運転監視計算機11に送信する。さらに、部品の点検や補修や交換の推奨情報を運転監視計算機11に送信する。これにより、機器の損傷状態を正確にかつ迅速に予測することができ、機器の交換や運用の仕方を決めることができるので、機器の信頼性を高めることができる。
【００７２】
損傷率がある基準値を越えた場合に出力する情報としては、発電設備名、部品名、損傷率、損傷率の増加の予測、部品の点検や補修や交換の推奨情報や予定などである。この例では、運転の予定、つまり、運転形態、運転時間、起動停止回数などを予め入力しておき、損傷解析の累計から損傷の推移を出力する。
【００７３】
図１０は運転時間、起動停止回数と損傷率の関係の表示例を示す。損傷率がある基準値を越えない場合は、損傷がある基準に達するまでの運転時間や起動停止回数を表示する。これにより、部品の交換や点検などの計画を立てることができる。図１０では、設計条件による損傷の予測も同時に示している。これにより、発電設備の損傷の仕方が設計条件よりも厳しいか否かが容易に分かり、場合によっては点検したり、運転方法を変えたりする参考になる。
【００７４】
次に、本発明の第２の実施例を説明する。図１１は第２の実施例による発電設備の遠隔損傷診断装置を示す。図１に示した第一の実施例との相違は、損傷の数値解析モデル１５ｃに対し、実験計画法から得られた近似式を用いている点にある。
【００７５】
本実施例では、損傷の数値解析手段１５ｃに実験計画法の直交表作成手段１５ｃ3を用い、センサ情報と運転情報を変数として直交表に割付け、損傷解析を実施する。損傷解析手段は、有限要素法解析モデル１５ｃ１による解析と損傷数値解析モデル１５ｃ２による解析の両方または一方で構成されている。
【００７６】
まず、センサ情報と運転情報を変数ｘｉとして、直交表への割り付けを行う。直交表とは、実験計画法で使用される試験や解析における変数の変化のさせ方を表にしたものである。使用する直交表の種類や直交表へ割り付ける変数の選択は、データ入出力インターフェイス15aにより行われる。直交表の選択は、選択した変数ｘｉの数とその水準数および変数ｘｉが有する２因子交互作用の数と影響を勘案して行う。
【００７７】
直交表の種類としては、L8（27），L16（215），L32（231），L9（34），L27（313），L81（340），L12（211），L18（2137），L36（211312，313）などが良く知られている。水準の設定方法は、変数ｘｉの基本統計量がデータベース１５eに格納されているときは、この基本統計量を参考にして定める。
【００７８】
たとえば、統計解析により、設計変数ｘｉの頻度分布が正規分布で近似でき、その平均値ｍｉとその分散ｓｉ2が得られたとする。このときの水準設定例として、３水準の設定のときは数１５に示すように、設計変数ｘｉの平均値ｍｉと設計変数ｘｉの標準編差ｓｉと係数aより行う。係数aは０より大きく１０以下程度の範囲とする。
【００７９】
【数１５】
0003788901.tif 000016
【００８０】
基本統計量がデータベース15eに納められていない設計変数ｘｉについては、類似データまたは経験から平均値ｍｉ，分散ｓｉ2，係数ａを定める。
【００８１】
図１２に変数とその水準を示す。変数は実験計画法解析手段１５ｃ３で水準に割り付けられ、データベース15eに格納される。解析は直交表に従って実施される。
【００８２】
図１３の直交表に解析内容を示す。変数番号は、図１２の行に示した変数番号と対応しており、行の番号は解析の番号を示している。解析番号１の解析は、図１２に示した全ての変数の水準を水準１とした解析を実施するという意味である。また、解析番号２の解析は、変数番号１〜１２を水準２、変数番号１３を水準１として解析するという意味である。
【００８３】
次に、損傷の数値解析手段１５ｃにおいては、実験計画法解析ステップ１５ｃ3で割り付けられた変数ｘｉの直交表に従い、解析を実施する。この解析は、部品の運転条件を模擬した有限要素法モデル15c1による熱構造解析と、その結果として得られる物理量を用いた損傷数値解析15c2である。先述のように、損傷解析としては破壊力学を用いたき裂進展解析や線形損傷則による累積損傷率の解析などがある。
【００８４】
損傷の数値解析手段15cで得られた損傷Ｄはデータベースに格納される。次に、分散分析手段１５ｃ４によりデータベースの損傷Ｄについて分散分析を行い、各変数の効果を算出する。分散分析とは、解析結果に含まれる情報を、誤差の部分と本質的な部分、つまり、変数の効果とに分離し、解析結果から統計的に意味のある結論を導けるか否かを吟味する手段である。平方和、自由度、平均平方、Ｆ比、Ｐ値、寄与率などを算出する。このＦ比、Ｐ値により、変数ｘｉのうち有意である変数を選択することができる。
【００８５】
次に、分散分析手段15c4で得られた解析結果を用いて、近似式作成手段１５ｃ5により損傷Ｄの近似式を作成する。近似式は数１６にて表され、２次の応答曲面を示している。ｂは回帰により得られた係数、εは誤差である。ここで，応答曲面を作成するに当たっては、分散分析解析手段１５ｃ４で得られたＦ比を用い、ある設定した有意水準、例えば１％または５％の設計変数について近似式を作成する。近似式作成手段１５ｃ４での変数入力や近似式の表示は、データ入出力インターフェイス１５aにて行う。
【００８６】
【数１６】
0003788901.tif 000017
【００８７】
損傷Ｄの近似式にセンサ情報と運転情報が入力されると、部品の損傷が算出される。
【００８８】
図１４に、第２の実施例の処理フローを示す。ステップ２１において、変数を選択する。変数としては、センサ情報と運転情報である。ステップ２２において、変数ｘｉを直交表への割り付けを行う。直交表に従って、ステップ２３、ステップ２４、ステップ２５、ステップ２６を実施する。つまり、直交表の行の数の損傷解析を実施する。
【００８９】
ステップ２３において、有限要素法による熱解析と構造解析を行う。ステップ２４において、有限要素法による熱解析と構造解析の解析結果として、損傷解析に用いる物理量を出力する。物理量としては、部材の温度，応力，ひずみなどがある。ステップ２５において、求められた物理量を用を損傷解析モデルに入力し、ステップ２６において、損傷解析を実施する。
【００９０】
ステップ２７において、損傷解析結果について分散分析を行う。分散分析により、各変数の効果を算出する。また、平方和、自由度、平均平方、Ｆ比、Ｐ値、寄与率などを算出する。Ｆ比、Ｐ値により、変数ｘｉのうち有意である変数を選択することができる。有意水準としては一般的に１％や５％が用いられる。これにより、直交表に割り付けた設計変数の中から、損傷Ｄに影響を与える変数ｘｉを選び出し、影響の無い変数は取り除いて、変数を減らすことができる。
【００９１】
ステップ２８において、損傷の近似式を作成する。近似式を作成するに当たっては、ステップ２７で得られたＦ比を用いて、ある設定した有意水準、例えば１％または５％の設計変数について近似式を作成する。
【００９２】
ステップ２９において、発電設備のセンサ情報を入手する。ステップ３０において、発電設備の運転情報を入手する。ステップ３１において、作成した近似式を用いて損傷解析を実施する。ステップ３２において、損傷結果を出力し、データベース１５ｅに書き込む。損傷解析結果としては疲労損傷率や、クリープ損傷率や、き裂長さ変形量など、部材の損傷を数値化したものである。
【００９３】
ステップ３３において、部材の損傷結果の累計をデータベースより読み込み、ステップ３１において算出した損傷解析結果を積算し、データベース１５ｅに書き込む。ステップ３４において、データベースより損傷解析結果の積算、つまり、損傷解析結果の累計を読み込み、基準値と比較する。
【００９４】
さらに、ステップ３５においては、運転の予定、たとえば、運転する時間や運転する回数や運転の形態を損傷解析に入力し、損傷を外挿する。このとき、近似式に与えるセンサ情報と運転情報は、予め与えてある設計条件としたり、過去のセンサ情報と運転情報から算出する。ステップ３６において、損傷解析結果の出力を行う。その他は実施例１と同様である。
【００９５】
本実施例では、損傷解析により作成した近似式を用いて損傷を算出するので、第１の実施例に比べ、計算機の負荷を低減し、解析を高速化することができるので、機器の損傷診断が迅速に行える。
【００９６】
次に、本発明の第３の実施例を説明する。図１５に第３の実施例による遠隔損傷診断装置を示す。本実施例は、実験計画法を用いて機器の損傷と比較し、数値解析モデルと機器の損傷を合わせる逆問題解析手段２０を設けている。
【００９７】
逆問題解析手段２０は直交表作成手段20a、有限要素法解析手段20b、損傷解析手段20c、分散分析手段20d、境界条件作成手段20e（図１の15bと同じ）より構成されている。直交表作成手段20aでは、有限要素法解析手段20bで用いる解析モデルの境界条件を変数として水準に割付け、直交表に従って解析を実施する。
【００９８】
図１６は図４で示した動翼の熱境界条件を水準に割付けた例を示す。翼の特徴的な部位ごと、例えば、翼の前縁、翼先端、翼面腹側、翼面背側、冷却孔をそれぞれ高さ方向に領域分けし、熱境界条件の領域を割当てる。ここでは、数値流体解析などで得られた熱伝達率に対する係数を変数として、各々３水準に変化させて、有限要素法の熱解析を実施する。
【００９９】
図１７に直交表を示す。この直交表による解析の例では、３６ケースの有限要素法熱解析を実施するものであり、たとえば、解析番号１はすべての熱伝達率を数値流体解析で得られた値に係数１を乗じた熱伝達率とする解析である。
【０１００】
有限要素法熱解析手段20bで得られた物理量と実機の損傷を比較して、その残差が最小となる境界条件を同定する。たとえば、動翼の場合、実機の動翼の温度は金属強化組織の粗大化量で推定できるため、実機の翼の部位ごとの組織を観察し、各部位の温度を推定することができる。
【０１０１】
次に、推定した各部の温度と数値解析手段１５ｃの温度の結果が最小となる境界条件を求める。たとえば、数１７に示す応答ｙの式を作成する。ここで、ｋは翼の部位番号、Ｔｃｋは数値解析で得られた部位ｋの温度、Ｔｉｋは実機翼の観察で得られた部位ｋの温度である。
【０１０２】
【数１７】
0003788901.tif 000018
【０１０３】
直交解析の結果を分散分析手段20dにおいて分散分析を実施し、解析と実機の温度差の残差の平方和を最小とする熱境界を境界条件作成手段20eで作成する。
【０１０４】
図１８に変数、つまり、熱境界の水準と実機と有限要素法解析の温度の残差平方和ｙの関係を示す。たとえば、変数記号Ａ〜Ｅの熱伝達率を、残差平方和が小さくなる値に設定する。そして、本熱境界条件をデータベースに格納し、損傷解析手段１５ｃにおける有限要素法解析モデル１５ｃ１の境界条件に入力する。逆問題解析手段２０で得られた熱伝達率は、数７と数９で示した熱伝達率Hgci、Hccjとして、有限要素法解析モデル１５ｃ１の境界条件に反映される。
【０１０５】
本実施例においては、有限要素法解析モデル１５ｃ１の境界条件が、逆問題解析手段２０により実機と一致するように修正されている。このため，損傷解析手段１５ｃがより高精度化され、部品の信頼性を高めることができる。
【０１０６】
また、き裂長さなどの損傷記録がデータベース１５eに存在する場合は、これらの損傷について実機データと解析とが合うように境界条件を変更してもよい。その場合は、損傷解析手段１５ｃを直交解析にて行い、得られた損傷と実機損傷の残差平方和を最小とする境界条件を算出する。ここでいう損傷とは、実機の寿命を支配する主要な損傷を数値化したもので、き裂寸法、き裂数、減肉寸法、減肉数、ピット寸法、剥離寸法、剥離数などである。
【０１０７】
また、直交表作成手段20aで変数として用いるものは、解析の境界条件であればよく、熱伝達率、温度、荷重、境界の剛性、境界の減衰、流体力、材料の劣化度合などが考えられる。なお、直交解析で実施する解析として、図１５では有限要素法解析手段１５ｃ１と損傷解析手段１５ｃ２を示したが、これに限られず、直交解析可能な数学モデルであればよい。
【０１０８】
図１９に第3の実施例による逆問題解析を含む処理フローを示す。ステップ４１において、変数を選択する。変数は、有限要素法熱解析や構造解析の境界条件や数値解析モデルに入力する変数や、損傷解析モデルに入力する変数である。ステップ４２において、変数ｘｉを直交表への割り付ける。直交表に従って、ステップ４３〜ステップ４６を実施する。つまり、直交表の列の数の損傷解析を実施する。
【０１０９】
ステップ４７において、有限要素法解析で得られる結果である温度や応力及び、ひずみと実機の検査記録から推定した結果を比較し、それらの残差が最小となる境界条件を実働境界条件として同定する。つまり、数１８に示す解析結果と実機検査による推定値との誤差平方和ｙを算出し、これを最小化する境界条件を求める。
【０１１０】
【数１８】
0003788901.tif 000019
【０１１１】
ここで、kは比較する変数の番号、ｘａ(k)は解析で得られた変数、ｘｉ(k)は検査記録から得られた変数である。同定した実働境界条件を用いた有限要素法による熱構造解析により、実機運転状態での温度・応力分布を精度よく予測することができる。比較する変数を数１１に示した損傷Dとし、損傷解析モデルの境界条件を同定することも可能である．このときは、数１９のｙを算出する。
【０１１２】
【数１９】
0003788901.tif 000020
【０１１３】
ここで、kは比較する変数の番号、Dａ(k)は解析で得られた損傷、Ｄｉ(k)は検査記録から得られた損傷である。
【０１１４】
ステップ４９において、分散分析を実施し、数１８や数１９に示したｙを最小化する変数ｘｉを選択する。ステップ５０においてｙを最小化する変数ｘｉを決定する。ステップ５１において、決定した変数ｘｉを有限要素法解析や数値解析や損傷解析に与え、これを逆問題解析モデルとして、データベースに登録する。
【０１１５】
登録した逆問題解析モデルを用いて，図２に示したフローに基づいて損傷解析を行う。つまり、図１９のフローにおいて、有限要素法解析の境界条件ｘｉを同定した場合は、この境界条件を同定した有限要素法モデルを図１５の有限要素法モデル15c1とする。図１９のフローにおいて、損傷解析の変数ｘｉを同定した場合は、この変数ｘｉを同定した損傷解析モデルを、図１５の損傷解析モデル15c2とする。
【０１１６】
なお、本実施例においては、直交表を用いた分散分析により、数１７、数１８、数１９に示した誤差平方和ｙを最小化する発明について述べた。さらに、本発明は、ｙを最小化する手法として、最適化手法、例えば、線形計画法、非線形計画法、遺伝的アルゴリズム等を用いたものについても含むものとする。
【０１１７】
【発明の効果】
本発明によれば、発電設備のセンサ情報と運転情報を損傷解析モデルの境界条件に変換し、運転環境が異なる発電設備の実際の運転状態における損傷解析を行うことができるので、ガスタービンの信頼性を高めることができる。
【図面の簡単な説明】
【図１】本発明の第１の実施例を示す遠隔損傷診断装置のシステム構成図。
【図２】第１の実施例の処理フロー図。
【図３】第１の実施例で利用するセンサ情報、運転所情報及びメンテナンス情報の例示図。
【図４】第１の実施例で利用する数値解析モデルの構成図。
【図５】第１の実施例で利用する数値解析モデルの境界条件の説明図。
【図６】第１の実施例で利用する数値解析モデルの解析結果のグラフ。
【図７】ガスタービン静翼の有限要素法の解析結果を示す図。
【図８】破壊力学解析のモデル図。
【図９】き裂進展による損傷解析結果のグラフ。
【図１０】運転時間、起動停止回数に対する損傷結果、損傷予測のグラフ。
【図１１】本発明の第２の実施例を示す遠隔損傷診断装置のシステム構成図。
【図１２】第２の実施例における変数と水準の例示図。
【図１３】第２の実施例で利用する直交表の例示図。
【図１４】第２の実施例の処理フロー図。
【図１５】本発明の第３の実施例を示す遠隔損傷診断装置のシステム構成図。
【図１６】第３の実施例における変数と水準の例示図。
【図１７】第３の実施例で利用する直交表の例示図。
【図１８】第３の実施例で利用する分散分析結果のグラフ。
【図１９】第３の実施例の逆問題解析ステップの処理フロー図。
【符号の説明】
１…圧縮機、２…燃焼器、３…タービン、４…発電機、５…排熱回収ボイラ、６…蒸気タービン、７…発電機、８…復水器、９…センサ信号送信ケーブル、１０…制御装置、１１…運転監視計算機、１２…発電設備側ファイヤーウォール、１３…インターネット、１４…遠隔監視施設側ファイヤーウォール、１５…遠隔損傷診断装置、15a…入出力インターフェイス、15b…境界条件算出手段、15c…損傷数値解析手段、15c1…有限要素法解析モデル、15c2…損傷数値解析モデル、15c3…直交表作成手段、15c4…分散分析手段、15c5…近似式作成手段、15d…損傷評価手段、15e…データベース、１６…タービン動翼、16a…翼前縁、16b…翼先端、16c…プラットフォーム、16d…シャンク、17a〜171…翼冷却孔、18a…翼前縁熱境界、18b…翼背側熱境界、18c…翼腹側熱境界、19a〜19l…翼冷却孔熱境界、20…逆問題解析手段、20a…直交表作成手段、20b…有限要素法解析手段、20ｃ…損傷解析手段、20d…分散分析手段、20e…境界条件作成手段。
【図１】
0003788901.tif 000021
【図２】
0003788901.tif 000022
【図３】
0003788901.tif 000023
【図４】
0003788901.tif 000024
【図５】
0003788901.tif 000025
【図６】
0003788901.tif 000026
【図７】
0003788901.tif 000027
【図８】
0003788901.tif 000028
【図９】
0003788901.tif 000029
【図１０】
0003788901.tif 000030
【図１１】
0003788901.tif 000031
【図１２】
0003788901.tif 000032
【図１３】
0003788901.tif 000033
【図１４】
0003788901.tif 000034
【図１５】
0003788901.tif 000035
【図１６】
0003788901.tif 000036
【図１７】
0003788901.tif 000037
【図１８】
0003788901.tif 000038
【図１９】
0003788901.tif 000039