JP3923773B2 - プラント異常事象診断装置とその診断方法、ならびに記録媒体 - Google Patents
プラント異常事象診断装置とその診断方法、ならびに記録媒体 Download PDFInfo
- Publication number
- JP3923773B2 JP3923773B2 JP2001317304A JP2001317304A JP3923773B2 JP 3923773 B2 JP3923773 B2 JP 3923773B2 JP 2001317304 A JP2001317304 A JP 2001317304A JP 2001317304 A JP2001317304 A JP 2001317304A JP 3923773 B2 JP3923773 B2 JP 3923773B2
- Authority
- JP
- Japan
- Prior art keywords
- fluctuation
- data
- abnormal event
- chemical
- plant
- Prior art date
- Legal status (The legal status 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 status listed.)
- Expired - Lifetime
Links
Images
Landscapes
- Testing And Monitoring For Control Systems (AREA)
Description
【発明の属する技術分野】
本発明は、プラントの化学管理に関する異常事象を診断するプラント異常事象診断装置とその診断方法、ならびに記憶媒体に関する。
【0002】
【従来の技術】
発電プラントの化学管理においては、従来、運転員や化学管理員等が、収集されたデータ、例えば、化学管理に関するプラントデータ、その化学分析データおよびオンライン化学分析装置データ等を直接監視することにより、異常事象の発生を示唆するデータ変動を判断し、このデータ変動から異常事象の予兆や発生を見い出している。
【0003】
異常事象が発生した場含には、運転員や化学管理員等による異常事象発生原因の調査が実施される。この調査では、運転員や化学管理員等が、個人のもつ知識と経験に基づいて、収集されたデータの中から異常事象と関連があると考えるデータを選び出し、評価に適した形式に編集加工し、その結果をベースにして調査を行う。
【0004】
上述の異常事象の発生を示唆するデータ変動の判断および異常事象発生原因の調査は、例えば、図10に示すように構成された評価装置100上で行われる。
【0005】
図10において、データ収集手段2は、想定するプラント異常事象の影響度合を求める発電プラントの化学管理に関するプラントデータ、化学分析データおよびオンライン化学分析装置データなどを収集する。データ収集手段2により収集されたデータは、評価装置100の周辺機器として接続されている表示装置10に、表示手段6によってリアルタイムで表示される。また、トレンドグラフが、収集されたデータからトレンドグラフ処理手段101によって作成され、表示手段6によって表示装置10や印刷出力装置12に表示される。なお、符号11は操作員等からの各種の指示やデータを受け入れる入力装置である。
【0006】
異常事象の発生を示唆するプラントデータや化学分析データの変動の判断は、運転員や化学管理員等が、プラントデータ、水質化学分析データ、あるいはトレンドグラフを表示装置10の画面上で直接監視することによって行われる。
【0007】
データ収集手段2により収集されたデータのうち、予め運転員や化学管理員等によって選択されたデータが、化学データファイルF2に保存される。
【0008】
異状事象が発生した場合には、その原因を調査するために、運転員や化学管理員等が、化学データファイルF2の内容を、編集用のアプリケーションプログラムを使って評価に適した形式のデータに編集加工し、その結果を表示装置10の画面を通して参照することにより、異状事象発生時のデータを評価する。
【0009】
【発明が解決しようとする課題】
従来、発電プラントの化学管理において、異状事象の発生を示唆するデータ変動の判断、および異状事象の発生原因の調査には、多数の運転員や化学管理員等の労力を要しており、管理コストがかさむという問題があった。
【0010】
また、運転員や化学管理員等による評価基準は、個人の主観が介入することによって曖昧になりやすく、信頼性の高い評価結果を得難いという問題があった。
【0011】
さらに、評価基準の曖昧さは、プラントデータの初期の軽微な変動(異状事象の予兆など)を見過ごす結果を招き、異状を発生初期にくい止めることができないという問題があった。
【0012】
本発明は、このような課題を解決するためになされたもので、発電プラントの化学管理において、異常事象の予兆や発生に関する評価を、客観的な基準の下で客観的なツールを用いて行うことを可能とし、これにより、運転員や化学管理員等による判断や調査の労力を軽減して評価に費やされるコストの引き下げを図るとともに、評価の精度を向上させることのできるプラント異状事象診断装置とその診断方法ならびに記憶媒体を提供することを目的とする。
【0013】
【課題を解決するための手段】
上記目的を達成するために、請求項1の発明のプラント異常事象診断装置は、プラントの化学管理に関するデータを収集するデータ収集手段と、データ収集手段により収集されたデータからデータの変動を検知するための基準値を求める基準値計算手段と、基準値計算手段により求められた基準値により、データ収集手段により収集されたデータの変動を検知するプラントデータ変動検知手段と、想定するプラント異常事象における化学パラメータの変動を予測する化学パラメータ変動シミュレーション手段と、化学パラメータ変動シミュレーション手段により計算され、作成されたシミュレーション結果ファイルを用いて、基準値計算手段により求められた基準値により、化学パラメータの変動を検知するシミュレーション変動検知手段と、シミュレーション変動検知手段により作成されたシミュレーション変動検知結果ファイルを用いて、一定時間経過後のシミュレーション変動検知結果データを読み込み想定する異常事象の発生により変動が検知された化学パラメータのリストを作成し、この化学パラメータのリストの一つ一つの化学パラメータが想定する異常事象を肯定する度合いおよび想定する異常事象により変動が予測される化学パラメータの不検出時否定影響度、想定する異常事象により変動が予測されない化学パラメータの検出時否定影響度を与えて予兆診断マトリックスを作成する予兆診断マトリックス作成手段と、データ収集手段により収集されたプラントデータの変動検知結果ファイルおよび予兆診断マトリックスとから変動が検知された化学パラメータ、変動が検知されなかった化学パラメータのリストに基づいて、変動が検出されたパラメータの組み合わせによる異常事象を肯定する検出時肯定適合度、検出されなかった化学パラメータの組み合わせによる異常事象を否定する不検出時否定適合度、検出された化学パラメータの組み合わせによる異常事象を否定する検出時否定適合度を求め、この求められた適合度を組み合わせてプラントの異常事象を診断する異常事象診断手段と、異常事象診断手段により得られた異常事象の診断結果を表示する表示手段とを具備することを特徴とする。
【0014】
ここで、化学パラメータとは、シミュレーション手段において用いる場合のデータ、例えば水の導電率等をいう。
【0015】
本プラント異常事象診断装置は、プラントが正常に稼働している期間に収集されたプラントの化学管理に関するデータ、例えば、プロセスデータや水質化学分析データから、そのデータ変動を評価するための基準値を計算する。この基準値をベースに収集データの変動の有無を判定する。
【0016】
データ変動評価のための基準値は、例えば、プラントが正常に稼働している期間に収集されたデータに基づいて、予め決められた計算方法により計算されるので、人の主観を排除した客観性のある値として求められる。
【0017】
異常事象を診断するための予兆診断マトリックスは、想定し得るプラントの異常事象に対して、予め化学パラメータ変動シミュレータを作成し、このシミュレータを操作して、変動を予測し、先の基準値を用いたデータ変動検知と同様の手法により、化学パラメータの変動の有無を判定する。この判定結果に基づいて、想定し得るプラントの異常事象に対して、一定時間経過後に、変動が検知されると予想される化学パラメータの組み合わせが得られる。
【0018】
プラントの異常事象は、化学パラメータ変動シミュレータから得られたマトリックスと実際のデータ変動の判定結果とから、所定の診断基準に従って診断される。
【0019】
本プラント異常事象診断装置によれば、客観性の高い異常事象診断結果を得ることができる。また、異常事象の発生の判断のみならず、予兆レベルの異常事象発生の判断を診断結果として得ることもできるので、異常の早期発見が可能となる。さらに、データ変動検知結果と化学パラメータの予兆診断マトリックスとの関連性を評価しながら異常事象の診断を行うので、異常事象の原因の究明が客観的でしかも容易である。加えて、診断のための人的負荷が大幅に軽減され、管理コストを低減することができる。
【0020】
請求項2記載のプラント異常事象診断装置は、請求項1記載のプラント異常事象診断装置において、化学パラメータ変動シミュレーション手段が、プラント運転時に発生する可能性のある異常事象毎に、予め用意した化学パラメータの変化を予測することを特徴とする。これにより、化学パラメータを化学関連パラメータに拡張することができる。
【0021】
請求項3記載のプラント異常事象診断装置は、請求項2記載のプラント異常事象診断装置において、化学パラメータ変動シミュレーション手段が、化学関連パラメータの変動シミュレーションの結果を用いて、シミュレーション結果ファイルを作成することを特徴とする。
これにより、化学パラメータ用のシミュレーション結果ファイルは、化学関連パラメータ用に拡張することができる。
【0022】
請求項4記載のプラント異常事象診断装置は、請求項1記載のプラント異常事象診断装置において、基準値計算手段が、データ収集手段を用いて収集されたプラント運転データ、化学分析データおよびオンライン化学分析装置データを貯える化学データファイルを用いて、データ種別毎の移動平均値とその標準偏差基準値を求め、これら平均値と標準偏差とで定義される基準値をプラントデータ種別毎の基準値として基準値データファイルを作成することを特徴とする。
【0023】
データ変動検知のための基準値を計算する方法には、移動平均、相関、積算確率等が挙げられる。これら異なる方法で複数の基準値を計算しておき、その中から想定し得る異状事象毎に、シミュレーションによって予測されるデータの変動傾向に合った基準値を選択するように構成する。
これにより、より一層信頼性の高い異常事象診断結果を得ることができる。
【0024】
請求項5記載のプラント異常事象診断装置は、請求項1記載のプラント異常事象診断装置において、プラントデータ変動検知手段が、データ収集手段により収集される化学データファイルを対象にし、データの変動を基準値データファイルに含まれる各データ毎の基準値と照合する際に、データ毎に予め定めた期間の移動平均値と移動平均値の変化量および変動検知結果を格納したプラントデータ変動検知結果ファイルを作成することを特徴とする。
これにより、プラントデータ変動検知手段の精度を向上させることができる。
【0025】
請求項6記載のプラント異常事象診断装置は、請求項3記載のプラント異常事象診断装置において、シミュレーション変動検知手段が、シミュレーション結果ファイルを対象に、シミュレーション値の変動を基準値データファイルに含まれる各化学関連パラメータ毎の基準値と照合する際に、化学関連パラメータ毎に予め定めた期間の移動平均値と移動平均値の変化量および変動検知結果を格納したシミュレーション変動検知結果ファイルを作成することを特徴とする。
これにより、シミュレーション変動検知手段の精度を向上させることができる。
【0031】
【発明の実施の形態】
(第1の実施の形態)(請求項1乃至6対応)
以下、本発明の実施形態を図面に基づいて説明する。
【0032】
図1は、本発明の第1の実施の形態におけるプラント異常事象診断装置の機能構成を示すブロック図である(請求項1対応)。
【0033】
同図に示すように、本実施形態のプラント異常事象診断装置1は、プラントの化学管理に関するプラント運転データやその化学分析データ(以下、プラントデータと呼ぶ。)を収集し、化学データファイルF2として格納するデータ収集手段2と、化学データファイルF2内のプラントデータからデータ変動検知用の基準値を計算し基準値データベースD1に格納する基準値計算手段3と、基準値データベースD1内の基準値を用いて化学データファイルF2に格納されるデータの変動を検知し、その検知結果をプラントデータ変動検知結果ファイルF3として格納するプラントデータ変動検知手段4と、想定されるプラント異常事象における化学パラメータの変動を計算し、シミュレーション結果ファイルS2として格納する化学パラメータ変動シミュレーション手段7と、基準値データベースD1内の基準値を用いてシミュレーション結果ファイルS2に格納された計算結果の変動を検知し、その検知結果をシミュレーション変動検知結果ファイルS3として記憶装置に格納するシミュレーション変動検知手段8と、シミュレーション変動検知結果ファイルS3を用いてプラントの異常事象診断手段に用いる予兆診断マトリックスM1を作成する予兆診断マトリックス作成手段9と、収集されたプラントデータ変動検知結果ファイルF3と予兆診断マトリックスM1とから、所定の診断基準に従って前記プラントの異常事象を診断する異常事象診断手段5により得られた異常事象の診断結果を表示装置10や印刷出力装置12に表示する表示手段6と、操作員等からの各種の指示やデータを受け入れる入力装置11とから構成されている。
【0034】
図1に示すプラント異常事象診断装置1は、具体的には、パーソナルコンピュータなどの汎用のコンピュータシステムや、ハード・ソフト的に専用化されたコンピュータシステムにより実現することができる。
【0035】
図2は、図1に示した基準値計算手段3の動作を示すフローチャートである(請求項4対応)。
【0036】
データ収集手段2は、決められた周期で、データ変動・異常事象の診断に必要な種類のプラントデータを収集し、化学データファイルF2として記憶領域に書き込む。
【0037】
基準値計算手段3は、化学データファイルF2から各周期のプラントデータを読み込み、このプラントデータについて、例えば、データ種別毎の移動平均値とその標準偏差を計算し、これら平均値と標準偏差とで定義される基準値をプラントデータ種別毎の基準値として基準値データファイルD1に格納する。
【0038】
基準値計算手段3は、データの移動平均値とその標準偏差による基準値を計算したが、これに限らず、移動平均値の変化量、相関、積算確率等の別の計算方法による基準値を生成するものであってよい。
【0039】
また、基準値計算手段3は、移動平均、移動平均値の変化量、相関、積算確率等のいくつかの計算方法により複数の異なる基準値を計算するものであってよい(請求項5対応)。
【0040】
複数の基準値を有効に運用する方法としては、例えば、プラントデータの種別毎に、そのプラントデータの性質に合った計算方法による最適な基準値を選択する方法、プラントの状態に合った計算方法による基準値を選択する方法などがある。これにより、診断結果の信頼性の向上を期待できる。
【0041】
基準値を計算する周期は、基準値の種別毎にプラントの運転状況やデータの変動の仕方に応じて、適宜最適な周期に変更できるように構成することが望ましい。また、データ種別毎に異なる計算周期を設定できるようにしてもよい。
【0042】
同一周期に得られた複数のプラントデータについて、これらデータ間の相関に基づいてデータ変動の診断を行う機能を加えてもよい。例えば、決められた複数のプラントデータ間の相関に一定の条件を定義しておき、あるプラントデータが変動して前記相関の一定条件を満足しなくなったとき、有意なデータ変動であると判断するようにしてもよい。
【0043】
図3は、図1に示したプラントデータ変動検知手段4の動作を示すフローチャートである。
【0044】
プラントデータ変動検知手段4は、化学データファイルF2より最新のプラントデータを読み込み、読み込んだプラントデータと、図4に示す基準値データファイルD1内の該当する基準値、例えば、移動平均基準値上限値「移動平均値+3σ」(σ;標準偏差)とを比較する。この比較の結果、読み込んだ計算結果データが「移動平均値+3σ」値以上になった場合は有意な変動が得られたものと見なし、そのプラントデータの変動内容(変動データ)をプラントデータ変動検知結果ファイルF3へ書き込む。そのプラントデータが移動平均基準値の上限を逸脱している場合には「増加」、移動平均基準値の下減を逸脱している場合には「減少」というデータ変動内容をプラントデータ変動検知結果ファイルF3へ書き込む(請求項5対応)。
このようにしてプラントデータ変動検知結果ファイルF3に各種のプラントデータの変動データを蓄積して行く。
【0045】
図5は、図1に示した化学パラメータ変動シミュレーション手段7の動作を示すフローチャートである(請求項6対応)。
【0046】
化学パラメータ変動シミュレーション手段7は、プラント運転時に発生する可能性のある事象毎に、予め用意した化学パラメータのシミュレーション計算を行うもので、化学パラメータの初期値情報を呼び込み、それぞれの事象が発生した時に想定される化学パラメータの変動計算結果をシミュレーション結果ファイルS2に書き込む。
【0047】
図6は、図1に示したシミュレーション変動検知手段8の動作を示すフローチャートである。
【0048】
図3に示したプラントデータ変動検知手段4の動作を示すフローチャートと同様に、シミュレーション変動検知手段8は、シミュレーション結果ファイルS2よりシミュレーション計算結果データを読み込み、このシミュレーション計算結果データと、図4に示す基準値データファイルD1内の該当する基準値、例えば、移動平均基準値上限値「移動平均値+3σ」(σ;標準偏差)とを比較する。この比較の結果、読み込んだ計算結果データが「移動平均値+3σ」値以上になった場合は有意な変動が得られたものと見なし、以下のように計算結果データの変動内容をシミュレーション変動検知結果ファイルS3へ書き込む。例えば、計算結果データが移動平均基準値の上限を逸脱している場合には「増加」、移動平均基準値の下減を逸脱している場合には「減少」という計算結果データの変動内容をシミュレーション変動検知結果ファイルS3へ書き込む。
このようにしてシミュレーション変動検知結果ファイルS3に各種のシミュレーション計算結果の変動傾向を蓄積する。
【0049】
(第2の実施の形態)
図7は、図1に示した予兆診断マトリックス作成手段9の動作を示すフローチャートである。
【0050】
予兆診断マトリックス作成手段9は、事前に想定した異常事象に対し、シミュレーション変動検知結果ファイルS3より一定時間経過後のシミュレーション変動検知結果データを読み込み、想定する事象の発生により変動が検知された化学パラメータのリストを作成する。
【0051】
ここで、一つ一つの化学パラメータが異常事象を肯定する度合いを、検出時肯定影響度aRi(i=1〜m; m=異常事象毎のシミュレーション変動検知手段による変動検知データの数)とし、シミュレーション変動検知手段で変動が検知された化学パラメータに1、検知されない化学パラメータに0を与えるものとする。
【0052】
例えば、原子力プラントにおける異常事象の一つとして、「主復水器海水漏洩発生」を想定し、特定規模の海水漏洩が主復水器で発生した時のシミュレーション計算を化学パラメータ変動シミュレーション手段を用いて行う。
【0053】
シミュレーション変動検知手段を用いて、シミュレーション計算結果の変動検知を実施する場合、「主復水器海水漏洩発生」という異常事象が発生して、変動が検知される化学パラメータが、例えば、「ホットウェル導電率」、「復水導電率」であるとすると、これらのパラメータの肯定影響度は1と定義する。
【0054】
「想定した異常事象により変動が予測される化学パラメータ」の検出時肯定影響度aRiの一覧を図8(a)に示す。これを予兆診断マトリックスM1と呼び、予め作成する。
【0055】
なお、予兆診断マトリックスM1において、シミュレーション変動検知手段により変動が検知された化学パラメータの基準値となる移動平均値に対するシミュレーション結果は、変化比を考慮して、変動検知データが異常事象に影響を与える肯定影響度を、1、0とともにその中間値を使って表現してもよい。すなわち、変化比が5倍以上の場合、肯定影響度は0.8、2倍以上の場合、肯定影響度は0.4、1.5倍以上の場合、肯定影響度は0.05などとおく。
【0056】
上述した例において、例えば、原子力プラントにおける「主復水器海水漏洩発生」を診断するとき、「ホットウェル導電率」の肯定影響度を0.8、「復水導電率」の肯定影響度を0.4と表現してもよい。これにより、各異常事象への適合度の値に微妙な差を与えることができ、各異常事象の間でのランク分けが可能になる。
【0057】
想定する異常事象に対し、化学パラメータの変動が、シミュレーション変動検知手段によって検知されたにもかかわらず、実際のプラントデータでは変動が検知されない場合には、異常事象を否定する度合いを不検出時否定影響度uRi(i=1〜m;mは異常事象毎のシミュレーション変動検知手段による変動検知データの数)とし、シミュレーション変動検知手段によって変動が検知された化学パラメータに1、検知されない化学パラメータに0を与えものとする。想定した異常事象により変動が予測される化学パラメータの不検出時否定影響度uRiの一覧を図8(b)の予兆診断マトリックスM1に示す。これも、予め作成することができる。
【0058】
さらに、想定する異常事象に対応した、シミュレーション変動検知手段によって変動が検知されない化学パラメータに対して、プラントデータ変動検知手段により変動が検知された場合には、異常事象を否定する度合いを検出時否定影響度nRi(i=1〜n;nは異常事象毎のシミュレーション変動検知手段による変動無検知データの数)とし、シミュレーション変動検知手段によって変動が検知されない化学パラメータに1、検知された化学パラメータに0を与えるものとする。
【0059】
想定した異常事象により変動が予測されない化学パラメータの検出時否定影響度nRiの一覧を図8(c)の予兆診断マトリックスM1に示す。これも、予め作成することができる。
【0060】
なお、不検出時否定影響度や検出時否定影響度は、検出時肯定影響度と同様1、0とともにその中間値を使って表現してもよい。
【0061】
(第3の実施の形態)
まず、異常事象診断手段5が、化学パラメータについて作成した予兆診断マトリックスファイルM1に基づいて総合適合度を算出するプロセスを説明する。
【0062】
図9は、図1に示した異常事象診断手段5の動作を示すフローチャートである。異常事象診断手段5は、プラントデータ変動検知結果ファイルF3より各化学パラメータの変動判断結果を読み込む。次に、予兆診断マトリックスファイルM1を読み込み、変動が検知された化学パラメータ、変動が検知されなかった化学パラメータのリストに基づいて、変動が検出されたパラメータの組み合わせによる異常事象を肯定する検出時肯定適合度aR、検出されなかった化学パラメータの組み合わせによる異常事象を否定する不検出時否定適合度uR、検出された化学パラメータの組み合わせによる異常事象を否定する検出時否定適合度nRを下式により計算する。
【0063】
aR=1−(1−aR1)×(1−aR2)×(1−aR3)×・・・×(1−aRm)
uR=1−(1−uR1)×(1−uR2)×(1−uR3)×・・・×(1−nRm)
nR=1−(1−nR1)×(1−nR2)×(1−nR3)×・・・×(1−nRn)
【0064】
これらの結果に基づき、総合的に異常事象発生の可能性を診断する総合適合度tRを下式により計算する。
【0065】
tR=aR−α×{β×uR+(1−β)nR}
【0066】
なお、ここで示すα、βは不検出時適合度uR、検出時否定適合度nRにより、想定される異常事象を否定する場合の重みを設定するもので0〜1の間の数値で設定する。
【0067】
異常事象診断手段7は、この計算をすべての異常事象について行い、その結果を診断データファイルF4に保存する。
【0068】
例えば、前記予兆診断マトリックス手段M1により、異常事象の「主復水器海水漏洩発生」時に、化学パラメータの「ホットウェル導電率の増加」と「復水導電率の増加」が検出され、化学パラメータの「復水貯蔵タンク水導電率の増加」が検出されないと評価される場合、「主復水器海水漏洩発生」という一つの異常事象について診断される化学パラメータは「ホットウェル導電率」と「復水導電率」となる。この結果、それぞれの化学パラメータの検出時肯定適合度をaR1、aR2、不検出時否定適合度をuR1、uR2とし、それぞれに1を与えるものとする。また、「復水貯蔵タンク水導電率」は検出時否定適合度をnR1とし1を与えるものとする。
【0069】
次に、異常診断手段5が、プラントデータについて作成したプラントデータ変動検知結果ファイルF3に基づいて総合適合度を算出するプロセスを説明する。
【0070】
プラントデータ変動検知手段4により「ホットウェル導電率の増加」と「復水導電率の増加」が検出され、「復水貯蔵タンク水導電率の増加」が検出されない場合、「主復水器海水漏洩発生」という一つの異常事象についての検出時肯定適合度aRを求めると、
aR=1−(1−1)×(1−1)=1
となり、検出時肯定適合度は1となる。一方、不検出時否定適合度uRを求めると、
uR=1−(1−0)×(1−0)=0
となり、不検出時否定適合度は0となる。また、検出時否定適合度nRを求めると、
nR=1−(1−0)=0
となり、検出時否定適合度は0となる。
【0071】
以上をまとめると、総合適合度tRは下式により
tR=aR−α×{β×uR+(β−1)nR}=1−0=1
1となり、「主復水器海水漏洩発生」の異常事象の予兆または発生があると判断されることになる。
【0072】
次に、プラントデータ変動検知手段により、「ホットウェル導電率の増加」、「復水導電率の増加」、および「復水貯蔵タンク水導電率の増加」が検出された場合、「主復水器海水漏洩発生」という一つの異常事象についての検出時肯定適合度aRを求めると、
aR=1−(1−1)×(1−1)=1
となり、検出時肯定適合度は1となる。
【0073】
一方、不検出時否定適合度uRを求めると、
uR=1−(1−0)×(1−0)=0
となり、不検出時否定適合度は0となる。また、検出時否定適合度nRを求めると、
nR=1−(1−1)=1
となり、検出時否定適合度は1となる。
【0074】
以上をまとめると、総合適合度tRは下式によりtR=aR−α×{β×uR+(1−β)nR}=1−α×(1−β)
となり、仮に、α=0.5、β=0.5と設定した場合、総合適合度は0.75となり、「主復水器海水漏洩発生」の異常事象の予兆または発生が判断されるものの、その可能性(適合度)はやや低下することになる。
【0075】
また、「ホットウェル導電率の増加」と「復水導電率の増加」のいずれも検出されない場合、「主復水器海水漏洩発生」という一つの異常事象についての適合度aRを求めると、
aR=1−(1−0)×(1−0)
で求められ、適合度aRは0となり、「主復水器海水漏洩発生」は診断されないことになる。
【0076】
このように、不検出時否定適合度と検出時否定適合度とを、先に求めた検出時肯定適合度から減ずることにより、実際の専門家が無意識に行っているような現実度の高い診断が可能になる。
【0077】
さらに、個々の変動データ毎の異常事象への検出時肯定影響度aRiや不検出時否定適合度uRi、検出時否定適合度nRiの値は、運転員や化学管理員等による入力により、適宜、設定・変更することが可能である。
【0078】
表示手段6は、各異常事象に対する適合度の計算結果を、例えば、プラントデータの変動内容(増加,減少など)とともに表示装置10または印刷装置12より出カすることにより、異常事象の予兆や発生を運転員や化学管理員等に提供する。
【0079】
(第4の実施の形態)(請求項8対応)
以上説明したように、プラント異常事象診断装置の特徴的な手段は、例えば、ソフトディスク、CD−ROM等のコンピュータによって読み取り可能な記憶媒体にアプリケーションプログラムとして記憶することができる。
【0080】
すなわち、図1に示した本実施形態のプラント異常事象診断装置は、このような記憶媒体からアプリケーションプログラムを読み込んで実行するコンピュータの形態で実現することができる。
【0081】
【発明の効果】
本プラント異常事象診断装置によれば、客観性の高い異常事象診断結果を得ることができる。また、異常事象の発生の判断のみならず、予兆レベルの異常事象発生の判断を診断結果として得ることもできるので、異常の早期発見が可能となる。さらに、データ変動検知結果と化学パラメータの予兆診断マトリックスとの関連性を評価しながら異常事象の診断を行うので、異常事象の原因の究明が客観的でしかも容易である。加えて、診断のための人的負荷が大幅に軽減され、管理コストを低減することができる。
【図面の簡単な説明】
【図1】本発明の第1の実施の形態におけるプラント異常事象診断装置の構成を示す図である。
【図2】図1の基準値計算手段の動作を示すフローチャートである。
【図3】図1のプラントデータ変動検知手段の動作を示すフローチャートである。
【図4】基準値データベースの構成を示す図である。
【図5】図1の化学パラメータ変動シミュレーション手段の動作を示すフローチャートである。
【図6】図1のシミュレーション変動検知手段の動作を示すフローチャートである。
【図7】図1の予兆診断マトリックス作成手段の動作を示すフローチャートである。
【図8】(a)は、検出時肯定影響度を登録するための予兆診断マトリックスの一例を、(b)は、不検出時否定影響度を登録するための予兆診断マトリックスの一例を、(c)は、検出時否定影響度を登録するための予兆診断マトリックスの一例を示す図である。
【図9】図1の異常事象診断手段の動作を示すフローチャートである。
【図10】従来のプラント異常事象診断装置の構成を示す図である。
【符号の説明】
1…プラント異常事象診断装置、2…データ収集手段。3…基準値計算手段、4…プラントデータ変動検知手段、5…異常事象診断手段、6…表示手段、7…化学パラメータ変動シミュレーション手段、8…シミュレーション変動検知手段、9…予兆診断マトリックス作成手段、10…表示装置、11…入力装置、12…印刷出力装置
Claims (8)
- プラントの化学管理に関するデータを収集するデータ収集手段と、
前記データ収集手段により収集されたデータからデータの変動を検知するための基準値を求める基準値計算手段と、
前記基準値計算手段により求められた基準値により、前記データ収集手段により収集されたデータの変動を検知するプラントデータ変動検知手段と、
想定するプラント異常事象における化学パラメータの変動を予測する化学パラメータ変動シミュレーション手段と、
前記化学パラメータ変動シミュレーション手段により計算され、作成されたシミュレーション結果ファイルを用いて、前記基準値計算手段により求められた基準値により、化学パラメータの変動を検知するシミュレーション変動検知手段と、
前記シミュレーション変動検知手段により作成されたシミュレーション変動検知結果ファイルを用いて、一定時間経過後のシミュレーション変動検知結果データを読み込み前記想定する異常事象の発生により変動が検知された化学パラメータのリストを作成し、この化学パラメータのリストの一つ一つの化学パラメータが前記想定する異常事象を肯定する度合いおよび前記想定する異常事象により変動が予測される化学パラメータの不検出時否定影響度、前記想定する異常事象により変動が予測されない化学パラメータの検出時否定影響度を与えて予兆診断マトリックスを作成する予兆診断マトリックス作成手段と、
前記データ収集手段により収集されたプラントデータの変動検知結果ファイルおよび前記予兆診断マトリックスとから変動が検知された化学パラメータ、変動が検知されなかった化学パラメータのリストに基づいて、変動が検出されたパラメータの組み合わせによる異常事象を肯定する検出時肯定適合度、検出されなかった化学パラメータの組み合わせによる異常事象を否定する不検出時否定適合度、検出された化学パラメータの組み合わせによる異常事象を否定する検出時否定適合度を求め、この求められた適合度を組み合わせてプラントの異常事象を診断する異常事象診断手段と、
前記異常事象診断手段により得られた異常事象の診断結果を表示する表示手段とを具備することを特徴とするプラント異常事象診断装置。 - 前記化学パラメータ変動シミュレーション手段は、プラント運転時に発生する可能性のある異常事象毎に、予め用意した化学パラメータの変化を予測することを特徴とする請求項1記載のプラント異常事象診断装置。
- 前記化学パラメータ変動シミュレーション手段は、化学関連パラメータの変動シミュレーションの結果を用いて、シミュレーション結果ファイルを作成することを特徴とする請求項2記載のプラント異常事象診断装置。
- 前記基準値計算手段は、データ収集手段を用いて収集されたプラント運転データを貯える化学データファイルを用いて、データ種別毎の移動平均値とその標準偏差基準値を求め、これら平均値と標準偏差とで定義される基準値をプラントデータ種別毎の基準値として基準値データファイルを作成することを特徴とする請求項1記載のプラント異常事象診断装置。
- 前記プラントデータ変動検知手段は、前記データ収集手段により収集される前記化学データファイルを対象にし、データの変動を前記基準値データファイルに含まれる各データ毎の基準値と照合する際に、データ毎に予め定めた期間の移動平均値と移動平均値の変化量および変動検知結果を格納したプラントデータ変動検知結果ファイルを作成することを特徴とする請求項1記載のプラント異常事象診断装置。
- 前記シミュレーション変動検知手段は、前記シミュレーション結果ファイルを対象に、シミュレーション値の変動を前記基準値データファイルに含まれる各化学関連パラメータ毎の基準値と照合する際に、化学関連パラメータ毎に予め定めた期間の移動平均値と移動平均値の変化量および変動検知結果を格納したシミュレーション変動検知結果ファイルを作成することを特徴とする請求項3記載のプラント異常事象診断装置。
- プラントの化学管理に関するデータを収集するデータ収集工程と、
前記データ収集手段により収集されたデータからデータの変動を検知するための基準値を求める基準値計算工程と、
前記基準値計算工程により求められた基準値により、前記データ収集工程により収集されたデータの変動を検知するプラントデータ変動検知工程と、
想定するプラント異常事象における化学パラメータの変動を予測する化学パラメータ変動シミュレーション工程と、
前記化学パラメータ変動シミュレーション手段により計算され、作成されたシミュレーション結果ファイルを用いて、前記基準値計算工程により求められた基準値により、化学パラメータの変動を検知するシミュレーション変動検知工程と、
前記シミュレーション変動検知工程により作成されたシミュレーション変動検知結果ファイルを用いて、一定時間経過後のシミュレーション変動検知結果データを読み込み前記想定する異常事象の発生により変動が検知された化学パラメータのリストを作成し、この化学パラメータのリストの一つ一つの化学パラメータが前記想定する異常事象を肯定する度合いおよび前記想定する異常事象により変動が予測される化学パラメータの不検出時否定影響度、前記想定する異常事象により変動が予測されない化学パラメータの検出時否定影響度を与えて予兆診断マトリックスを作成する予兆診断マトリックス作成工程と、
前記データ収集工程により収集されたプラントデータの変動検知結果ファイルおよび前記予兆診断マトリックスとから変動が検知された化学パラメータ、変動が検知されなかった化学パラメータのリストに基づいて、変動が検出されたパラメータの組み合わせによる異常事象を肯定する検出時肯定適合度、検出されなかった化学パラメータの組み合わせによる異常事象を否定する不検出時否定適合度、検出された化学パラメータの組み合わせによる異常事象を否定する検出時否定適合度を求め、この求められた適合度を組み合わせてプラントの異常事象を診断する異常事象診断工程と、
前記異常事象診断工程により得られた異常事象の診断結果を表示する表示工程とを有することを特徴とするプラント異常事象診断方法。 - プラントの化学管理に関するデータを収集するデータ収集手順と、
前記データ収集手段により収集されたデータからデータの変動を検知するための基準値を求める基準値計算手順と、
前記基準値計算手順により求められた基準値により、前記データ収集手順により収集されたデータの変動を検知するプラントデータ変動検知手順と、
想定するプラント異常事象における化学パラメータの変動を予測する化学パラメータ変動シミュレーション手順と、
前記化学パラメータ変動シミュレーション手順により計算され、作成されたシミュレーション結果ファイルを用いて、前記基準値計算手段により求められた基準値により、化学パラメータの変動を検知するシミュレーション変動検知手順と、
前記シミュレーション変動検知手順により作成されたシミュレーション変動検知結果ファイルを用いて、一定時間経過後のシミュレーション変動検知結果データを読み込み前記想定する異常事象の発生により変動が検知された化学パラメータのリストを作成し、この化学パラメータのリストの一つ一つの化学パラメータが前記想定する異常事象を肯定する度合いおよび前記想定する異常事象により変動が予測される化学パラメータの不検出時否定影響度、前記想定する異常事象により変動が予測されない化学パラメータの検出時否定影響度を与えて予兆診断マトリックスを作成する予兆診断マトリックス作成手順と、
前記データ収集手順により収集されたプラントデータの変動検知結果ファイルおよび前記予兆診断マトリックスとから変動が検知された化学パラメータ、変動が検知されなかった化学パラメータのリストに基づいて、変動が検出されたパラメータの組み合わせによる異常事象を肯定する検出時肯定適合度、検出されなかった化学パラメータの組み合わせによる異常事象を否定する不検出時否定適合度、検出された化学パラメータの組み合わせによる異常事象を否定する検出時否定適合度を求め、この求められた適合度を組み合わせてプラントの異常事象を診断する異常事象診断手順と、
前記異常事象診断手順により得られた異常事象の診断結果を表示する表示手順とをコンピュータに実行させるプログラムを記録したコンピュータ読み取り可能な記録媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001317304A JP3923773B2 (ja) | 2001-10-15 | 2001-10-15 | プラント異常事象診断装置とその診断方法、ならびに記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001317304A JP3923773B2 (ja) | 2001-10-15 | 2001-10-15 | プラント異常事象診断装置とその診断方法、ならびに記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003122429A JP2003122429A (ja) | 2003-04-25 |
JP3923773B2 true JP3923773B2 (ja) | 2007-06-06 |
Family
ID=19135198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001317304A Expired - Lifetime JP3923773B2 (ja) | 2001-10-15 | 2001-10-15 | プラント異常事象診断装置とその診断方法、ならびに記録媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3923773B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI274436B (en) * | 2006-02-23 | 2007-02-21 | Iner Aec Executive Yuan | Method of supplying fuel to fuel cells |
JP5871157B2 (ja) * | 2011-10-03 | 2016-03-01 | 株式会社Ihi | 遠心圧縮設備のサージング防止方法 |
-
2001
- 2001-10-15 JP JP2001317304A patent/JP3923773B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2003122429A (ja) | 2003-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Which log level should developers choose for a new logging statement? | |
JP4361634B2 (ja) | 複数の診断関連情報を統合するためのシステム及び方法 | |
KR101348324B1 (ko) | 분석 시스템 및 보수 점검 지원 장치 | |
JP4873985B2 (ja) | 設備機器用故障診断装置 | |
EP0786725A1 (en) | Improvements in or relating to the prediction of product defects | |
JPH10510385A (ja) | ソフトウエア品質のアーキテクチャに基づく分析のための方法およびシステム | |
US7853546B2 (en) | Enhanced rule execution in expert systems | |
US20200143292A1 (en) | Signature enhancement for deviation measurement-based classification of a detected anomaly in an industrial asset | |
US20080275747A1 (en) | Incident/accident report analysis apparatus and method | |
Kitchenham et al. | Design metrics in practice | |
JP4973952B2 (ja) | プラント診断方法、プラント診断装置およびプラント診断用プログラム | |
JP3923773B2 (ja) | プラント異常事象診断装置とその診断方法、ならびに記録媒体 | |
JP2003315213A (ja) | 振動診断装置およびその方法 | |
Kleinbaum et al. | Assessing discriminatory performance of a binary logistic model: ROC curves | |
JPH086635A (ja) | プラント診断装置 | |
JPH0713617A (ja) | 不具合事象の原因推定方法 | |
Simon | A three-dimensional receiver operator characteristic surface diagnostic metric | |
JP3053902B2 (ja) | 異常原因診断方法 | |
CN116628564B (zh) | 用于检测发电机状态的模型训练方法和系统 | |
JPH07174617A (ja) | 回転機器健全性診断システム | |
JP7153585B2 (ja) | 異常原因推定方法、および、異常原因推定装置 | |
JPH04169849A (ja) | データ判別用しきい値設定方式 | |
CN117104453A (zh) | 一种智能舵机的能力等级确定方法、装置及设备 | |
CN115235749A (zh) | 一种机器设备故障诊断方法及装置 | |
JPH07146806A (ja) | テスト作業管理支援装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20040319 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040810 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061122 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061128 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070129 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070220 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070222 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 3923773 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100302 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110302 Year of fee payment: 4 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110302 Year of fee payment: 4 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120302 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130302 Year of fee payment: 6 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130302 Year of fee payment: 6 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130302 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140302 Year of fee payment: 7 |