Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH07113995B2 - Inference method and device - Google Patents
[go: Go Back, main page]

JPH07113995B2 - Inference method and device - Google Patents

Inference method and device

Info

Publication number
JPH07113995B2
JPH07113995B2 JP1004526A JP452689A JPH07113995B2 JP H07113995 B2 JPH07113995 B2 JP H07113995B2 JP 1004526 A JP1004526 A JP 1004526A JP 452689 A JP452689 A JP 452689A JP H07113995 B2 JPH07113995 B2 JP H07113995B2
Authority
JP
Japan
Prior art keywords
state
degree
event
inference
calculation
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 - Fee Related
Application number
JP1004526A
Other languages
Japanese (ja)
Other versions
JPH02183895A (en
Inventor
明 鍛治
武一 丸山
順一 鈴木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP1004526A priority Critical patent/JPH07113995B2/en
Publication of JPH02183895A publication Critical patent/JPH02183895A/en
Publication of JPH07113995B2 publication Critical patent/JPH07113995B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Devices For Executing Special Programs (AREA)
  • Alarm Systems (AREA)

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は推論方法および装置に係り、特に異常の波及の
予測に配慮した、プラントの運転支援に好適な推論方法
および装置に関する。
The present invention relates to an inference method and apparatus, and more particularly to an inference method and apparatus suitable for plant operation support in consideration of prediction of anomaly spread.

〔従来の技術〕[Conventional technology]

従来、プラントの運転支援に用いられている推論装置
は、特開昭62−75720号公報に記載の発明のように、異
常伝播を表わす情報により推論を行うもので、各異常事
象に実測値を反映することは行っていない。
Conventionally, the inference device used to support the operation of a plant, as in the invention described in Japanese Patent Laid-Open No. 62-75720, makes inference based on information representing anomalous propagation, and an actual measurement value is set for each abnormal event. It has not been reflected.

〔発明が解決しようとする課題〕[Problems to be Solved by the Invention]

上記従来技術では、プラントの異常発生時、異常原因か
らその波及を推論するが、プラントの状態が異常の波及
を反映して時々刻々に変化していく点を考慮しておら
ず、精度が充分でなかった。即ち、異常の発火点から異
常の伝播情報に基づいて推論しただけでは、異常の波及
経路上の各事象の実際の状態が適切に反映されず、精度
のよい異常波及の予測ができない。
In the above-mentioned conventional technology, when an abnormality occurs in the plant, the ripples are inferred from the cause of the abnormality, but the fact that the state of the plant changes momentarily reflecting the ripples of the abnormality is not considered, and the accuracy is sufficient. It wasn't. That is, only by inferring from the firing point of the abnormality based on the propagation information of the abnormality, the actual state of each event on the propagation path of the abnormality cannot be properly reflected, and accurate prediction of the abnormality spread cannot be performed.

本発明の目的は、異常が発生し、その波及経路上のプラ
ント状態が変動している場合でも、発生した異常の伝播
によるプラント状態の変化を精度よく推論するにある。
An object of the present invention is to accurately infer a change in the plant state due to the propagation of the anomaly even when the anomaly occurs and the plant state on the transmission path is fluctuating.

〔課題を解決するための手段〕 上記の課題は、プラント事象の状態量を測定し、測定さ
れた状態量が当該プラントの他の事象の状態量におよぼ
す波及効果を推論する方法において、前記プラント事象
の状態量の測定が時系列的に連続して行われ、該状態量
の変化が、該状態量の変化よりも時間的にあとの前記他
の事象の状態量に及ぼす波及効果の推論が予め定められ
た演算のためのデータに基づいて連続して行なわれる推
論方法により達成される。
[Means for Solving the Problem] The above problem is to measure the state quantity of a plant event, and infer the ripple effect that the measured state quantity has on the state quantities of other events of the plant, in the plant. The state quantity of the event is continuously measured in time series, and the reason why the change of the state quantity has a ripple effect on the state quantity of the other event that is later in time than the change of the state quantity can be inferred. This is achieved by an inference method that is continuously performed based on data for a predetermined calculation.

また、プラント事象が他の事象に波及する経路上の事象
の状態量が連続して測定され、推論された前記プラント
事象が他の事象に及ぼす波及効果が、前記波及経路上の
事象の状態量に基づいて補正される請求項1に記載の推
論方法としてもよい。
In addition, the state quantity of the event on the path where the plant event spreads to other events is continuously measured, and the ripple effect that the inferred plant event has on other events is the state quantity of the event on the spread path. The inference method according to claim 1 that is corrected based on

さらに、推論された波及効果に基づいてプラント事象の
状態をしめす数値が推論され、測定されたプラント状態
量から算出されたプラント事象の状態をしめす数値と前
記推論された数値との差を減少させるように予め定めら
れた演算のためのデータが変更される手順を備えている
請求項1および2に記載の推論方法としてもよい。
Further, a numerical value indicating the state of the plant event is deduced based on the inferred ripple effect, and the difference between the numerical value indicating the state of the plant event calculated from the measured plant state quantity and the inferred numerical value is reduced. The inference method according to claim 1 or 2 may further include a procedure for changing the data for the predetermined calculation.

上記の課題は、また、プラントの下位事象の状態が上位
事象に及ぼす効果を推論する装置に、プラント事象の状
態量を連続して記憶するプラントデータ記憶手段と、該
記憶手段に接続して設けられ前記状態量を事象の状態度
をしめす数値に変換する状態度演算手段と、該状態度演
算手段に接続して設けられ前記下位事象の状態度をしめ
す数値に基づいて該下位事象の状態が上位事象に及ぼす
波及効果を状態度をしめす数値として出力する状態度波
及推論手段と、前記状態度演算手段に接続して設けられ
状態度をしめす数値を記憶する状態度記憶手段と、該状
態度記憶手段と前記状態度波及推論手段とに接続して設
けられ状態度波及推論手段が出力する数値を状態度記憶
手段が記憶している数値に基づいて補正する状態度合成
手段と、を備えることによっても達成される。
The above-mentioned problem is also provided in a device that infers the effect of the state of a lower event of a plant on an upper event, a plant data storage unit that continuously stores the state quantity of a plant event, and a device connected to the storage unit. State degree calculating means for converting the state quantity into a numerical value indicating the degree of state of the event, and the state of the lower order event based on the numerical value indicating the degree of state of the lower order event provided connected to the state degree calculating means. A state degree influence reasoning means for outputting a ripple effect exerted on an upper event as a numerical value indicating the degree of state, a state degree storing means connected to the state degree calculating means for storing a numerical value indicating the state degree, and the state degree State degree synthesis means for correcting the numerical value output by the state degree influence inference means, based on the numerical value stored in the state degree storage means, and connected to the storage means and the state degree influence inference means. It is also achieved by the.

また、状態度演算手段および状態度波及推論手段に接続
して設けられ前記各手段に演算のためのデータを供給す
る知識ベースと、状態度演算手段および状態度波及推論
手段に接続して設けられ状態度演算手段が出力する状態
度の数値と状態度波及推論手段が出力する状態度の数値
との差を演算する状態度誤差検出手段と、該誤差検出手
段および前記知識ベースに接続して設けられ前記差を減
少させるように前記知識ベースに記憶された演算のため
のデータを書き換える学習装置と、を備えている請求項
4に記載の推論装置としてもよい。
Further, a knowledge base provided to be connected to the state degree calculating means and the state degree influence inference means and for supplying data for calculation to each of the means, and provided to be connected to the state degree calculating means and the state degree influence inference means State degree error detecting means for calculating the difference between the state degree numerical value output by the state degree calculating means and the state degree numerical value output by the state degree influence inference means; and a state degree error detecting means connected to the error detecting means and the knowledge base. And a learning device that rewrites the data for the operation stored in the knowledge base so as to reduce the difference, and the inference device according to claim 4.

さらに、学習装置が、状態度誤差検出手段に接続されプ
ラント事象の状態度を時系列に記憶する来歴状態度記憶
手段と、該来歴状態度記憶手段に接続され誤差検出手段
が検出する差が零になる演算のためのデータを算出する
最適配分演算手段と、該最適配分演算手段に接続して設
けられ演算のためのデータを時系列に記憶する来歴演算
データ記憶手段と、該来歴演算データ記憶手段および前
記最適配分演算手段に接続して設けられ最適配分演算手
段が出力する演算のためのデータと来歴演算データ記憶
手段が出力する演算のためのデータとから平均化された
演算のためのデータを知識ベースに出力する演算データ
平均化手段と、を備えている請求項5に記載の推論装置
としてもよい。
Further, the learning device is connected to the degree-of-state error detecting means and stores the degree of state of the plant event in time series, and the difference detected by the error detecting means connected to the history-state degree storing means is zero. Optimal distribution calculation means for calculating data for calculation, history calculation data storage means connected to the optimum distribution calculation means for storing data for calculation in time series, and history calculation data storage Means and data for arithmetic operation provided by the optimal allocation arithmetic means and output by the optimal allocation arithmetic means, and data for arithmetic operation averaged from data for arithmetic operation output by the history arithmetic data storage means. And an arithmetic data averaging means for outputting to the knowledge base.

〔作用〕[Action]

プラント事象の状態量が連続して測定され、測定された
事象(下位事象)の状態量が他の事象(上位事象)に及
ぼす波及効果が、連続的に推論されるので、下位事象の
状態量の変動に基づく、上位事象への波及効果の大きさ
が常に下位事象の状態量を反映して推論される。
The state quantity of a plant event is continuously measured, and the ripple effect of the measured state quantity of an event (lower event) on another event (upper event) is continuously inferred. It is inferred that the size of the ripple effect on the higher-ranking event based on the fluctuation of the above always reflects the state quantity of the lower-ranking event.

また、下位事象の状態量が上位事象の状態量に、その間
に介在するひとつ以上の事象を経て影響するとき、その
間の事象の状態量も連続的に測定され、その値に基づい
て、下位事象が上位事象に及ぼす波及効果の推論量が補
正される。
Also, when the state quantity of a lower event affects the state quantity of a higher event through one or more intervening events, the state quantity of the event during that time is also continuously measured, and based on that value, the lower event The inference amount of the ripple effect that the above has on the upper event is corrected.

さらに、プラント事象の実際の状態と、推論された状態
とが比較され、推論された状態が実際の状態に近ずくよ
うに、推論の演算のデータを変更する手順が設けられる
ので、推論が繰り返えされると、次第に前記演算のデー
タがプラント事象の特性に合った値に近ずく。
Furthermore, a procedure is provided to compare the actual state of the plant event with the inferred state, and to change the data of the inference operation so that the inferred state approaches the actual state. When returned, the data of the operation gradually approaches a value that matches the characteristics of the plant event.

〔実施例〕〔Example〕

以下、本発明を適用した推論装置の一実施例を第1図に
より説明する。本実施例の推論装置10は、監視対象であ
るプラント1の各部の事象の状態量を示すセンサーから
の信号を、プロセス入力装置2を介して周期的に取り込
み、工学単位に変換した後格納するプラントデータ記憶
手段であるプラントデータベース3と、運転員のプラン
ト運転に係るノーハウを記憶した知識ベース4と、該知
識ベース4および前記プラントデータベース3に接続さ
れた推論処理装置5と、該推論処理装置5に接続された
推論結果記憶メモリ6とを備えている。プラントデータ
ベース3は、従って最大取り込み周期の遅れでプラント
の現在状態を反映しており、本装置は、プラントデータ
ベース3のプラント現在値を常に反映しながら知識ベー
スに従って推論処理を行う。推論処理の結果は、推論結
果記憶メモリ6に格納され、該推論結果記憶メモリ6に
接続されたCRT表示器7を介して運転員に知らされる。
An embodiment of the inference apparatus to which the present invention is applied will be described below with reference to FIG. The inference apparatus 10 of the present embodiment periodically takes in a signal from a sensor indicating the state quantity of an event in each part of the plant 1 to be monitored via the process input device 2, converts it into an engineering unit, and stores it. A plant database 3, which is a plant data storage unit, a knowledge base 4 that stores know-how related to plant operation of an operator, an inference processing device 5 connected to the knowledge base 4 and the plant database 3, and the inference processing device. And an inference result storage memory 6 connected to 5. Therefore, the plant database 3 reflects the current state of the plant with a delay of the maximum fetch period, and the present apparatus performs inference processing according to the knowledge base while always reflecting the plant current value of the plant database 3. The result of the inference processing is stored in the inference result storage memory 6 and notified to the operator via the CRT display 7 connected to the inference result storage memory 6.

次に本装置における推論処理を第2図により説明する。
本装置は、プロセスの状態をプロセス入力装置を介して
プラントデータベース3に入力されるプロセス状態信号
に基づいて常時監視し、プラント状態に何らかの異常が
検知されたら、その異常検知時点を起点(発火点)とし
て、知識ベース4に格納されたデータに基づいて、発生
した異常の波及を推論する。以下、プラントの状態をそ
の状態が異常かそうでないかを示す異常度で表現する。
推論処理装置5は、プラントデータベース3に接続され
た状態度演算手段であるプロセス異常度演算装置5A、該
プロセス異常度演算装置5Aに接続されたプロセス異常度
判定装置5B、該プロセス異常度判定装置5Bに接続された
状態度波及推論手段である異常度波及推論装置5C、該異
常度波及推論装置5Cに接続された状態度合成手段である
異常度合成装置5D、前記プロセス異常度演算装置5Aと前
記異常度合成装置5Dに接続された状態度記憶手段である
プロセス異常度記憶装置5Eを備えており、前記プロセス
異常度演算装置5A、異常度波及推論装置5Cおよび異常度
合成装置5Dは、それぞれ前記知識ベース4に接続され、
異常度合成装置5Dには推論結果記憶メモリである異常波
及経路情報記憶装置6が接続されている。
Next, the inference process in this device will be described with reference to FIG.
This device constantly monitors the process state based on the process state signal input to the plant database 3 via the process input device, and if any abnormality is detected in the plant state, the abnormality detection time point is set as the starting point (ignition point). ), Based on the data stored in the knowledge base 4, the spread of the abnormality that has occurred is inferred. Hereinafter, the state of the plant is expressed by an abnormality degree indicating whether the state is abnormal or not.
The inference processing device 5 includes a process abnormality degree calculating device 5A which is a state degree calculating means connected to the plant database 3, a process abnormality degree determining device 5B connected to the process abnormality degree calculating device 5A, and the process abnormality degree determining device. Anomaly degree spread inference device 5C which is a state degree spread reasoning means connected to 5B, an abnormality degree synthesis device 5D which is a state degree synthesis means connected to the abnormality degree spread reasoning device 5C, and the process abnormality degree computing device 5A The abnormality degree synthesis device 5D is provided with a process abnormality degree storage device 5E that is a state degree storage means connected to the abnormality degree synthesis device 5D, and the process abnormality degree calculation device 5A, the abnormality degree spread reasoning device 5C, and the abnormality degree synthesis device 5D, respectively. Connected to the knowledge base 4,
The abnormal degree synthesis device 5D is connected to an abnormal spread path information storage device 6 which is a memory for storing inference results.

プロセス異常度演算装置5Aは、運転員、プラントメーカ
のノーハウとして知識ベース4に格納された、プロセス
量と異常度の関係を表わすメンバーシップ関数により、
プロセス量を、0から1で正規化された異常度に変換す
る。例えば、横軸にプロセス量Q、縦軸に異常度Eをと
った第3図に示されるように、プロセス量が、250℃の
場合には、それに対応したプロセス異常度は、知識ベー
ス4に格納されているメンバーシップ関数Mを用いて0.
8(80%)と求められる。メンバーシップ関数は、プロ
セス量の特性に合った形で運転員等のノーハウに基づい
て定められる。プロセス異常度判定装置5Bは、プロセス
異常度演算装置5Aで演算されたプロセス異常度が、該プ
ロセスが異常であるかどうかを判定するしきい値を越え
ているかどうかを判定する。第4図に異常波及を推論す
るための知識例を示す。図ではA5なる異常事象が発生す
るための原因をA4からA1までさかのぼって示している。
ここで、A1〜A5(以下Aiという)は、プロセス量データ
Biから変換された異常度αiを有する。プロセス異常度
判定装置5Bで、ある事象例えばA2が異常と判定された
時、異常度波及推論装置5Cは、第4図のルートに従って
A2の異常がA3の異常度α3にどう波及するかを推論す
る。この場合、A2の異常度α2から推論されるA3の異常
度と、プロセス量B3から演算されるA3の異常度が異常度
合成装置5Dにより合成され、この結果に基づいて、A3
予測異常度Zが推論される。
The process abnormality degree computing device 5A uses the membership function representing the relationship between the process amount and the abnormality degree, which is stored in the knowledge base 4 as the know-how of the operator and the plant maker.
The process amount is converted into an abnormal degree normalized from 0 to 1. For example, as shown in FIG. 3 in which the horizontal axis represents the process amount Q and the vertical axis represents the abnormality degree E, when the process amount is 250 ° C., the process abnormality degree corresponding to it is stored in the knowledge base 4. 0 using the stored membership function M.
It is required to be 8 (80%). The membership function is determined based on the know-how of the operator etc. in a form that matches the characteristics of the process amount. The process abnormality degree determining device 5B determines whether the process abnormality degree calculated by the process abnormality degree calculating device 5A exceeds a threshold value for determining whether the process is abnormal. FIG. 4 shows an example of knowledge for inferring an abnormal influence. The figure shows retroactively cause for A5 becomes abnormal event is generated from the A 4 to A 1.
Here, A 1 to A 5 (hereinafter referred to as A i ) are process amount data
It has anomaly degree α i converted from B i . When a certain event, for example, A 2 is determined to be abnormal by the process abnormality degree determining device 5B, the abnormality degree influence inferring device 5C follows the route of FIG.
Infer how the anomaly of A 2 affects the anomaly degree α 3 of A 3 . In this case, the abnormality degree of A 3 inferred from the abnormality degree α 2 of A 2 and the abnormality degree of A 3 calculated from the process amount B 3 are synthesized by the abnormality degree synthesizing device 5D, and based on this result, A A predicted anomaly degree Z of 3 is inferred.

第5図はこの処理の流れをフローチャートで示したもの
で、最終段の事象まで推論が行われたら、その結果が異
常度波及経路情報記憶装置5Eに格納され、ガイドメッセ
ージ等としてCRT表示器7に表示される。
FIG. 5 is a flow chart showing the flow of this processing. When the event at the final stage is inferred, the result is stored in the abnormality degree transmission route information storage device 5E and is displayed as a guide message or the like on the CRT display 7 Is displayed in.

事象A1から事象A5迄の異常の伝播する様子およびそれを
推論により予測する様子を第6A〜6C図に示す。まず、時
刻t1において、A1の異常が検知され、この異常の波及に
よるA2〜A5の異常度が推論される。時刻t1における状態
を示す第6A図では、A1の異常が発生したばかりで、まだ
A2以降に波及していないため、A2〜A5のプロセス異常度
α(図の△)は正常に近い値を示している。一方、推論
された異常度R(図の○)によれば、A2,A3の異常の発
生が予測されている。合成異常度(図の×)は、それぞ
れの事象の現在の状態を示すプロセス異常度△と前段の
推論された異常度による波及効果とが合成された値で、
これが当該段の予測異常度Zとなる。
Figures 6A-6C show the propagation of anomalies from event A 1 to event A 5 and how they are predicted by inference. First, at time t 1 , the abnormality of A 1 is detected, and the degree of abnormality of A 2 to A 5 due to the spread of this abnormality is inferred. In Figure 6A, which shows the condition at time t 1 , the anomaly at A 1 has just occurred and is still
Since it does not spread to A 2 and later, the process abnormality degree α (Δ in the figure) of A 2 to A 5 shows a value close to normal. On the other hand, according to the inferred abnormality degree R (◯ in the figure), the occurrence of the abnormality of A 2 and A 3 is predicted. The composite anomaly degree (x in the figure) is a value obtained by combining the process anomaly degree Δ indicating the current state of each event and the ripple effect due to the anomaly degree inferred in the previous stage,
This is the predicted abnormality degree Z of the stage.

時刻t2を示す第6B図になると、A2に異常が波及し、A2
プロセス異常度△が推論異常度○よりも高くなってい
る。A3の異常度の推論には、A2のプロセス異常度とA2
推論異常度が合成された、合成異常度Z(図の×)が用
いられ、推論の精度が高められている。異常度の合成
は、本実施例ではプロセス異常度と推論異常度の最大を
とったが、両者の加重平均などの方法もある。
When the Product 6B showing the time t 2, the abnormality is spread to A 2, the process abnormality of the A 2 △ is higher than ○ inference error probability. The inference of the abnormality of the A 3, process error probability inference degree of abnormality A 2 of A 2 is synthesized, (× in the figure) Synthesis abnormality degree Z is used, the accuracy of the inference is enhanced. In the present embodiment, the maximum of the process abnormality degree and the inference abnormality degree is taken as the synthesis of the abnormality degree, but there is a method such as a weighted average of both.

時刻t3を示す第6C図では、A1の異常が、A3以降の事象に
波及してきている。
In FIG. 6C showing the time t 3 , the abnormality of A 1 is spreading to the events after A 3 .

次に異常波及推論に学習装置を付加し、より正確に波及
予測を行う仕組みの例を第7図に示す。この例では、第
1図に示した装置に加えて、プロセス異常度演算装置5A
及び異常度波及推論装置5Cに接続された状態度誤差検出
手段である異常度誤差検出装置5Hと、該異常度誤差検出
装置5Hと知識ベース4に接続された学習装置5Gとが、設
けられている。これは、オンライン推論とは別に、推論
の知識ベースを、推論結果とプロセス異常度により補正
していくもので、異常度誤差検出装置5Hでは、プラント
の異常現象が整定した時点で、プロセス異常度演算装置
5Aが出力するプロセス異常度αと、異常度波及推論装置
5Cが出力する推論異常度Rの差(誤差という)が演算さ
れ、得られた結果が学習装置5Gに送られる。学習装置5G
は異常度誤差検出装置5Hから送られてくる誤差を最少と
するように、知識ベース4に格納された異常波及伝播係
数βを変更する。ここで、変更された異常波及伝播係数
βに従って推論が行われ、再度誤差が評価されるという
学習が継続され、プロセスの特性に合った異常波及伝播
係数βが定められる。
Next, FIG. 7 shows an example of a mechanism in which a learning device is added to the extraordinary spread inference to more accurately predict the spread. In this example, in addition to the device shown in FIG.
And an abnormality degree error detecting device 5H which is state degree error detecting means connected to the abnormality degree influence inferring device 5C, and a learning device 5G connected to the abnormality degree error detecting device 5H and the knowledge base 4. There is. In addition to online inference, this is to correct the inference knowledge base based on the inference result and process anomaly level.In the anomaly level error detection device 5H, the process anomaly level is set when the plant abnormal phenomenon is settled. Arithmetic unit
Process anomaly degree α output by 5A and anomaly degree spread reasoning device
The difference (called an error) between the inference abnormality degrees R output by 5C is calculated, and the obtained result is sent to the learning device 5G. Learning device 5G
Changes the abnormal propagation coefficient β stored in the knowledge base 4 so as to minimize the error sent from the abnormality degree error detection device 5H. Here, the inference is performed according to the changed abnormal propagation coefficient β, the learning that the error is evaluated again is continued, and the abnormal propagation coefficient β that matches the characteristics of the process is determined.

以下に、上記推論の詳細を、火力発電プラントを例にと
り、第8図を参照して説明する。第8図はプラントに発
生する異常を予測するための知識を、フォールト・トリ
ー(FT)として因果関係を表現したものである。例えば
事象“ケーシング熱的変形”は、事象“負荷変化率
大”、“蒸気条件急変”、“ヒータ温度特高”のいずれ
かを原因として発生することを意味している。また、各
事象内には、事象のプロセス異常度Yを導くための関係
式と、自事象から上位事象(ここでは事象間の位置関係
を表わす表現として、フォールト・トリーの頂点方向、
すなわち、図の左方向を上位、逆方向を下位とする)へ
の波及の程度を示す異常波及伝播係数βが定義される。
The details of the above inference will be described below with reference to FIG. 8 by taking a thermal power plant as an example. FIG. 8 shows the causal relationship as the fault tree (FT), which is the knowledge for predicting the abnormality that occurs in the plant. For example, the event "casing thermal deformation" means that it occurs due to any of the events "large load change rate", "steam condition sudden change", and "heater temperature extraordinary". Further, in each event, a relational expression for deriving the process abnormality degree Y of the event and an upper event from the self event (here, as an expression representing the positional relationship between the events, the direction of the top of the fault tree,
That is, the abnormal ripple propagation coefficient β is defined, which indicates the degree of ripple to the left in the figure as the upper side and the opposite direction as the lower side.

事象のプロセス異常度Yは、その事象が意味する状態を
示したもので、0から1の間の連続量で表される。例え
ばY=0とはその事象が異常である程度は0、すなわち
正常であることを表わし、逆にY=1であれば、その事
象は全く異常であることを表している。事象のプロセス
量異常度はプロセス入力値から、そのプロセス量によっ
て表現される状態に基いて定められた関数(メンバシッ
プ関数)に従って求められる。これらの異常度を持った
事象群により、上位方向に推論(前向推論)が行われ、
最終段の観測事象である振動の発生する程度(異常度)
が予測される。
The process abnormality degree Y of an event indicates a state that the event means, and is represented by a continuous amount between 0 and 1. For example, Y = 0 means that the event is abnormal to some extent 0, that is, normal, and if Y = 1, it means that the event is totally abnormal. The process amount abnormality degree of an event is obtained from a process input value according to a function (membership function) defined based on the state represented by the process amount. Reasoning (forward reasoning) is performed in the upper direction by the event group having these abnormalities,
Extent of occurrence of vibration, which is the last observation event (abnormality)
Is predicted.

事象“ケーシング熱的変形”の異常度を、その下位事象
群、“負荷変化率大”、“蒸気条件急変”、“ヒータ温
度特高”の異常度に基づいて推論するために、これら事
象群の関係が、プロダクション・ルール(以下ルールと
いう)、として表わされる。“負荷変化率大”から“ケ
ーシング熱的変形”が発生する程度が“負荷変化率大”
の異常度から100%伝達するとすれば、この関係は、以
下のルールとして表わされる。
In order to infer the anomalous degree of the event "casing thermal deformation" based on the anomaly levels of its subordinate event groups, "large load change rate", "sudden change in steam condition", and "heater temperature high", these event groups Is expressed as a production rule (hereinafter referred to as a rule). The degree of occurrence of "casing thermal deformation" from "large load change rate" is "large load change rate"
Assuming that 100% is transmitted from the anomaly degree of, this relationship is expressed as the following rule.

もし、“負荷変化率大”ならば、“ケーシング熱的
変形”が起こる。
If the "load change rate is large", "casing thermal deformation" occurs.

同様に“蒸気条件急変”から“ケーシング熱的変形”が
発生する程度、“ヒータ温度特高”から“ケーシング熱
的変形”が発生する程度は、以下のようなルールで表わ
される。
Similarly, the degree of occurrence of "casing thermal deformation" due to "sudden change of steam condition" and the extent of occurrence of "casing thermal deformation" due to "heater temperature high" are expressed by the following rules.

もし、“蒸気条件急変”ならば、“ケーシング熱的
変形”が起こる。
If "steam condition sudden change", "casing thermal deformation" occurs.

もし、“ヒータ温度特高”ならば、“ケーシング熱
的変形”が起こる。
If "heater temperature is very high", "casing thermal deformation" occurs.

ここで、複数の下位事象が上位事象に波及するとき、そ
の波及の度合はコンバイン関数cfにより表わされる。コ
ンバイン関数cf(x1x2)は、入力x1,x2が整数である場
合には、cf(x1,x2)=x1+x2−x1*x2で表わされる。
例えば a“負荷変化率大”の異常度Y1=0.3 b“蒸気条件急変”の異常度Y2=0.4 c“ヒータ温度特高”の異常度Y3=0.5 とすると、“ケーシング熱的変形”の異常度Y4は、〜
のルールとa〜dの異常度とから、下記のように推論
される。
Here, when a plurality of lower events spread to upper events, the degree of the spread is represented by a combine function cf. The combine function cf (x 1 x 2 ) is represented by cf (x 1 , x 2 ) = x 1 + x 2 −x 1 * x 2 when the inputs x 1 and x 2 are integers.
For example, if the "amount of load change" abnormality is Y 1 = 0.3 b "steam condition sudden change" Y 2 = 0.4 c "heater temperature extraordinary" abnormality Y 3 = 0.5, "casing thermal deformation""Abnormality level Y 4 is ~
It is inferred as follows from the rule of 1 and the abnormalities of a to d.

i)ルール実行 Y4′=cf(Y4,Y1)=cf(0,0.3) …(1) =0.3−0*0.3=0.3 ii)ルール実行 Y4″=cf(Y4′,Y2)=cf(0.3,0.4) …(2) =0.3+0.4−0.3*0.4=0.58 iii)ルール実行 Y4=cf(Y4″,Y3)=cf(0.58,0.5) …(3) =0.58+0.5−0.58*0.5=0.79 上述のように、“負荷変化率大”の異常度が0.3、“蒸
気条件急変”の異常度が0.4、“ヒータ温度特高”の異
常度が0.5であれば、これら下位事象の異常の波及によ
り、最終的に“ケーシング熱的変形”が異常となる程度
は、79%と推論される。尚、上記(1),(2),
(3)式で、cf関数に上位事象の異常度Y4,Y4′,Y4
が入っているが、これは単に計算ロジック上のことであ
る。
i) Rule execution Y 4 ′ = cf (Y 4 , Y 1 ) = cf (0,0.3) (1) = 0.3−0 * 0.3 = 0.3 ii) Rule execution Y 4 ″ = cf (Y 4 ′, Y) 2) = cf (0.3,0.4) ... (2) = 0.3 + 0.4-0.3 * 0.4 = 0.58 iii) rule execution Y 4 = cf (Y 4 " , Y 3) = cf (0.58,0.5) ... (3 ) = 0.58 + 0.5-0.58 * 0.5 = 0.79 As mentioned above, the abnormality of "Large load change rate" is 0.3, the abnormality of "Rapid change of steam condition" is 0.4, and the abnormality of "heater temperature high" is If it is 0.5, it is inferred that the extent to which "casing thermal deformation" becomes abnormal is 79% due to the spread of anomalies of these subordinate events. In addition, the above (1), (2),
In Eq. (3), the cf function contains the abnormalities Y 4 , Y 4 ′, Y 4 ″ of the upper events.
Is included, but this is just in the calculation logic.

次に第9図により、推論処理装置の動作を詳細に説明す
る。
Next, the operation of the inference processing apparatus will be described in detail with reference to FIG.

(1)各事象のプロセス量が、各事象の異常を定義した
異常度定義式に基いて、あるスキャン速度で監視され
る。プロセス量は異常度定義式に定められたメンバシッ
プ関数に基いて、0から1の範囲のプロセス異常度αj
に変換される。
(1) The process amount of each event is monitored at a certain scan speed based on an abnormality degree definition formula that defines the abnormality of each event. The process amount is based on the membership function defined in the abnormality definition formula, and the process abnormality in the range of 0 to 1 α j
Is converted to.

(2)異常度定義式に基いて算出される各事象のプロセ
ス異常度αjが、異常判定レベル(しきい値)THを超え
た場合、異常事象発生と判定され、異常と判定された事
象を発火点として上位事象への異常波及推論が開始され
る。
(2) If the process abnormality degree α j of each event calculated based on the abnormality degree definition formula exceeds the abnormality determination level (threshold value) TH, it is determined that an abnormal event has occurred, and the event is determined to be abnormal. With the ignition point as the starting point, anomalous influence reasoning on higher-order events is started.

(3)異常が下位事象から上位事象に伝播,波及するに
は時間遅れを伴う場合が多く、予測に遅れを生む原因と
なるが、実際に当該事象が異常となる前に、下位事象群
の異常推論値を用いて予測を行うことにより、時間遅れ
が補われる。従って事象jの推論結果に基づく予測異常
度Zjは、下位事象群からの異常波及を推論した推論異常
度Rjと、プロセス量から定義されるプロセス異常度αj
とが合成されたもので、 Zj=cpf(αj,Rj) cpf:合成関数 である。
(3) Propagation and propagation of anomalies from lower-level events to higher-level events often involves a time delay, which causes a delay in prediction. However, before the event actually becomes abnormal, the lower-level event group The time delay is compensated by making a prediction using the anomaly inference value. Therefore, the predictive abnormality degree Z j based on the inference result of the event j is the inference abnormality degree R j that infers the abnormality spread from the lower event group and the process abnormality degree α j defined by the process amount.
And are combined, and Z j = cpf (α j , R j ) cpf: a combined function.

(4)事象jから上位事象kへの異常度の伝播量Xjは、
事象jの予測異常度Zjに異常波及伝播係数βjを乗じて
算出される。すなわち、 Xj=βj*Zj である。但し、異常の波及が整定した状態では、Xj=β
j*αjと考えてよい。
(4) The propagation amount X j of the degree of abnormality from the event j to the upper event k is
It is calculated by multiplying the predicted abnormal degree Z j of the event j by the abnormal ripple propagation coefficient β j . That is, X j = β j * Z j . However, when the ripple of the anomaly is settled, X j = β
You can think of it as j * α j .

(5)上位事象kへ影響する下位事象が複数(J1〜Jn
ある場合、それぞれの事象の異常度の伝播量Xj1〜Xjn
組み合わされて、事象kの推論異常度Rkとなる。Rkは下
位事象群J1,……Jnと上位事象kとの関係を示すコンバ
イン関数cfを用いて、 Rk=cf(Xj1,Xj2,…,Xjn) で示される。ここでコンバイン関数cfは、上位/下位事
象間の関係により、下記のように定義される。
(5) Lower events that affect the order event k multiple (J 1 through J n)
In some cases, the propagation amounts X j1 to X jn of the abnormalities of the respective events are combined to form the inference abnormality R k of the event k. R k is represented by R k = cf (X j1 , X j2 , ..., X jn ) by using a combine function cf showing the relationship between the lower event group J 1 , ... J n and the upper event k. Here, the combine function cf is defined as follows according to the relationship between upper / lower events.

(i)OR:cf(x1,x2)=1−(1−x1)*(1−x2) (ii)AND:cf(x1,x2)=x1*x2 事象kの予測異常度Zkは、事象kのプロセス量から定義
されるプロセス異常度αkと、前記推論異常度Rkから、
合成関数cpfにより、 Zk=cpf(αk,Rk) で示される。
(I) OR: cf (x 1, x 2) = 1- (1-x 1) * (1-x 2) (ii) AND: cf (x 1, x 2) = x 1 * x 2 events k the prediction error probability Z k of the process abnormality degree alpha k, which is defined from the process of event k, from the inference error probability R k,
The composite function cpf, represented by Z k = cpf (α k, R k).

事象kから更に上位の事象への異常の伝播量Xkは、異常
度波及伝播係数βkを用いて、 Xk=βk*Zk で示され、上位事象をLとし、事象kがk1〜knあるとす
れば、事象Lの推論異常度RLは、コンバイン関数cfを用
い、 RL=cf(Xk1,Xk2,…Xkn) で示される。さらに事象Lの予測異常度ZLは、合成関数
cpfを用い、 ZL=cpf(αL,RL) で示される。
Propagation quantity X k abnormal to events further higher from the event k, using the abnormal degree spread propagation coefficient beta k, indicated by X k = beta k * Z k, the order event is L, event k is k if there 1 to k n, the inference error probability R L of events L uses a combined function cf, R L = cf (X k1, X k2, ... X kn) represented by. Furthermore, the predictive anomaly degree Z L of event L is the composite function
It is shown by Z L = cpf (α L , R L ) using cpf.

(6)上記(3)〜(5)を繰り返すことにより、フォ
ールト・トリーでの波及推論が行われる。
(6) By repeating the above (3) to (5), the ripple reasoning in the fault tree is performed.

(7)異常波及路の決定には、発火点から生じた各波及
路のランク計算が行われ、ランク値(RNK)の高いもの
が上位波及路とされる。ランク計算は、予測異常度Zを
同一波及路上で平均をとることにより行われる。
(7) To determine the abnormal spread route, the rank calculation of each spread route generated from the ignition point is performed, and the one with the higher rank value (RNK) is determined as the upper spread route. The rank calculation is performed by averaging the predictive abnormalities Z on the same propagation route.

n:発火点より最終段の事象までの同一波及路上の事象数 次に、学習装置5Gの構成の実施例を第10図に示す。異常
度誤差検出装置5Hに接続して、来歴状態度記憶手段であ
る来歴異常度記憶装置5GAが設けられ、該来歴異常度記
憶装置5GAに最適配分演算手段である最適配分演算装置5
GBが接続されている。最適配分演算装置5GBには来歴演
算データ記憶手段である来歴伝播係数記憶装置5GCが接
続され、該来歴伝播係数記憶装置5GCの出力側には分散
処理装置5GDが接続されている。さらに、前記最適配分
演算装置5GBの出力側および来歴伝播係数記憶装置5GCに
接続して演算データ平均化手段である伝播係数平均化装
置5GEが設けられている。該伝播係数平均化装置5GEの出
力側は、異常波及経路情報記憶装置5E、知識ベース4お
よび来歴伝播係数記憶装置5GCに接続されている。この
構成の学習装置5Gを備えた第7図に示される推論装置に
おいて、プロセス異常度演算装置5Aで求められた事象j
のプロセス異常度αjは、異常度誤差検出装置5Hを経由
して、学習装置5G内の来歴異常度記憶装置5GAに送ら
れ、記憶される。最適配分演算装置5GBは、来歴伝播係
数記憶装置5GCに記憶されている事象ごとの時系列異常
度波及伝播係数βjeと、来歴異常度記憶装置5GAに時系
列に記憶されているプロセス異常度αjを読みこみ、こ
れに基づいて最適な中間伝播係数βjmを定め、伝播係数
平均化装置5GEへ出力する。伝播係数平均化装置5GEは、
この中間伝播係数βjmが入力されると、対応する時系列
異常度波及伝播係数βjeを来歴伝播係数記憶装置5GCか
ら読み出し、入力された中間伝播係数βjmと時系列異常
度波及伝播係数βjeとの平均値βjm/2+βje/2を求め、
これを今回の学習による伝播係数βjaとして異常波及経
路情報記憶装置5E、知識ベース4および来歴伝播係数記
憶装置5GCへ出力する。来歴伝播係数記憶装置5GCへ入力
された伝播係数βjaは最新の時系列異常度波及伝播係数
βjeとして記憶され、次回の学習処理に用いられ、知識
ベース4に送られたものは、オンライン推論のデータと
して用いられる。これら時系列異常度波及伝播係数βje
は、プロセスの特性に合ったものであれば、分散はゼロ
になり、従って時系列異常度波及伝播係数βjeの分散は
未知量によるものとなり、該伝播係数βjeの確かさを示
すものとなる。分散処理装置5GDは、来歴伝播係数記憶
装置に入力された伝播係数βjaの分散の度合を分散値と
して求め、この分散値を推論における「あいまい量」と
して出力/表示することにより、より精度の高い推論結
果が運転員に提示される。
n: the number of events on the same propagation route from the ignition point to the final stage event Next, FIG. 10 shows an example of the configuration of the learning device 5G. Connected to the abnormality degree error detection device 5H, a history abnormality degree storage device 5GA which is history state degree storage means is provided, and the history abnormality degree storage device 5GA has an optimum distribution calculation device 5 which is optimum distribution calculation means.
GB is connected. The optimal distribution arithmetic device 5GB is connected with a history propagation coefficient storage device 5GC, which is history calculation data storage means, and a distributed processing device 5GD is connected to the output side of the history propagation coefficient storage device 5GC. Further, a propagation coefficient averaging device 5GE, which is a calculation data averaging means, is provided, which is connected to the output side of the optimum distribution computing device 5GB and the history propagation coefficient storage device 5GC. The output side of the propagation coefficient averaging device 5GE is connected to the abnormal propagation route information storage device 5E, the knowledge base 4, and the history propagation coefficient storage device 5GC. In the inference device shown in FIG. 7 equipped with the learning device 5G having this configuration, the event j determined by the process abnormality degree computing device 5A
The process abnormality degree αj is transmitted to the history abnormality degree storage device 5GA in the learning device 5G via the abnormality degree error detection device 5H and is stored therein. The optimal distribution arithmetic unit 5GB consists of the time-series abnormality degree propagation coefficient β je for each event stored in the history propagation coefficient storage device 5GC and the process abnormality degree α stored in time series in the history abnormality coefficient storage device 5GA. It reads j , determines the optimum intermediate propagation coefficient β jm based on this, and outputs it to the propagation coefficient averaging device 5GE. The transmission coefficient averaging device 5GE is
When this intermediate propagation coefficient β jm is input, the corresponding time-series anomaly degree propagation coefficient β je is read from the history propagation coefficient storage device 5GC, and the input intermediate propagation coefficient β jm and the time-series anomaly degree propagation coefficient β j are input. Find the average value β jm / 2 + β je / 2 with je ,
This is output to the abnormal propagation path information storage device 5E, the knowledge base 4, and the history propagation coefficient storage device 5GC as the propagation coefficient β ja by this learning. The propagation coefficient β ja input to the history propagation coefficient storage device 5GC is stored as the latest time-series anomaly spread propagation coefficient β je , used for the next learning process, and sent to the knowledge base 4 for online inference. Used as data. These time-series anomaly degree spread propagation coefficient β je
Is zero if it matches the characteristics of the process, and therefore the variance of the time series anomalous degree propagation coefficient β je is due to an unknown quantity, and indicates the certainty of the coefficient β je. Become. The distributed processing device 5GD obtains the degree of dispersion of the propagation coefficient β ja input to the history propagation coefficient storage device as a dispersion value, and outputs / displays this dispersion value as an “ambiguous amount” in inference, thereby providing a higher accuracy. High inference results are presented to the operator.

本学習装置の処理フローを第11図に示す。ステップ1で
は、異常現象が整定した時点での時系列t0,t1,…tm
おいて、下位事象J1〜Jnから伝播波及される上位事象K
の推論異常度とのプロセス異常度誤差が演算される。α
j1t0,αj2t0,…αjnt0は時刻t0における下位事象J1
Jnのプロセス異常度で、αkt0は時刻t0における上位事
象Kのプロセス異常度である。βj1〜βjnは、事象J1
Jnの前記異常度波及伝播係数βjeであり、この伝播係数
と時系列t0〜tnにおけるプロセス異常度αj1t0,…,α
j2t0,…,αjnt0,…,αj1tn,…,αjntn,を用いて
cf関数で演算された事象Kの推論異常度と実際の事象K
の時刻t0〜tmにおけるプロセス異常度αkt0,…αktm
の誤差(すなわち推論誤差)を0とする連立方程式がβ
j1〜βjnを未知数として立てられる。αj1t0,αj2t0
…,αjnt0,…,αj1tm,…,αjntm,およびαkt0
…αktmはプロセス異常度で時刻t0〜tmで既知である。
The processing flow of this learning device is shown in FIG. In step 1, in the time series t 0 , t 1 , ... T m at the time when the abnormal phenomenon is settled, the upper event K propagated from the lower events J 1 to J n is propagated.
The process error degree error with the inference error degree is calculated. α
j1 t 0 , α j2 t 0 , ... α jn t 0 are lower events J 1 ~ at time t 0
In the process abnormality degree of J n , α k t 0 is the process abnormality degree of the upper event K at time t 0 . β j1 to β jn are event J 1 to
The anomaly degree transmission coefficient β je of J n , and the process abnormality degree α j1 t 0 , ..., α in the time series t 0 to t n .
Using j2 t 0 ,…, α jn t 0 ,…, α j1 t n ,…, α jn t n ,
Reasoning anomaly degree of event K calculated by cf function and actual event K
Of time t 0 ~t process in the m degree of abnormality α k t 0, ... α k t m,
The system of equations that makes the error of
It is possible to set j1 to β jn as unknowns. α j1 t 0 , α j2 t 0 ,
…, Α jn t 0 ,…, α j1 t m ,…, α jn t m , and α k t 0 ,
... α k t m is a process abnormality degree and is known from time t 0 to t m .

ステップ2ではステップ1で立てられた連立方程式の解
として、(m−n+1)組の異常度波及伝播係数が求め
られる。
In step 2, as the solution of the simultaneous equations established in step 1, (m-n + 1) sets of anomalous degree propagation factors are obtained.

ステップ3では、この(m−n+1)組の異常度波及伝
播係数群から、各組の変更量Eが求められ、ステップ4
で変更量Eが最小となる異常度波及伝播係数(βj1′,
βj2′,…,βjn′)が1組求められる。この
(βj1′,βj2′,…,βjn′)が前述の中間伝播係
数,βjmである。
In step 3, the change amount E of each group is obtained from the (m-n + 1) group of the degree-of-abnormality propagation factors, and step 4
Anomalousness propagation coefficient (β j1 ′,
One set of β j2 ′, ..., β jn ′) is obtained. This (β j1 ′, β j2 ′, ..., β jn ′) is the aforementioned intermediate propagation coefficient, β jm .

ステップ5で、ステップ4で求められた中間伝播係数
と、来歴伝播係数記憶装置5GCに時系列に記憶されてい
る異常度波及伝播係数の平均値が算出され、得られた値
(βj1,βj2,…,βjn)が学習による伝播係数βja
して、異常度波及経路情報記憶装置、来歴伝播係数記憶
装置、および知識ベース4へ送られる。以上のプロセス
により得られた伝播係数βjaが、オンライン推論の伝播
係数として用いられる。
In step 5, the average value of the intermediate propagation coefficient obtained in step 4 and the anomalous degree propagation coefficient stored in the history propagation coefficient storage device 5GC in time series is calculated, and the obtained values (β j1 , β j2 , ..., β jn ) are sent to the abnormality degree propagation route information storage device, the history propagation coefficient storage device, and the knowledge base 4 as the propagation coefficient β ja by learning. The propagation coefficient β ja obtained by the above process is used as the propagation coefficient for online inference.

〔発明の効果〕〔The invention's effect〕

本発明によれば、推論にプラント事象の推論時の実際の
状態量を用いているので、プラントの実際の状態が反映
された精度の高い推論値が得られ、予測されるプラント
状態に対応することを可能とする効果がある。また、学
習機能を備えて、プラントの特性にあった知識ベースを
構築し、推論の信頼性を高める効果がある。
According to the present invention, since the actual state quantity at the time of inference of a plant event is used for inference, a highly accurate inference value that reflects the actual state of the plant is obtained and corresponds to the predicted plant state. There is an effect that makes it possible. In addition, it has a learning function and builds a knowledge base that matches the characteristics of the plant, and has the effect of increasing the reliability of inference.

【図面の簡単な説明】[Brief description of drawings]

第1図は本発明を適用した推論装置の全体構成を示す概
念図、第2図は第1図に示された推論装置の推論処理部
分の詳細構成を示すブロック図、第3図はメンバシップ
関数の例を示すグラフ、第4図は事象の異常が波及・伝
播するようすを概念的に表わすフロー図、第5図は推論
処理手順を示すフローチャート、第6A〜6C図は事象の異
常の波及・伝播と予知のようすを示す概念図、第7図は
学習機能を備えた本発明の実施例の構成を示すブロック
図、第8図はフォールト・トリー図、第9図はフォール
ト・トリーによる異常波及予知を示すフロー図、第10図
は学習装置の詳細構成例を示すブロック図であり、第11
図は学習機能の演算手順を示すフロー図である。 3…プラントデータ記憶手段(プラントデータベー
ス)、4…知識ベース、5…推論処理装置、5A…状態度
演算手段(プロセス異常度演算装置)、5C…状態度波及
推論手段(異常度波及推論装置)、5D…状態度合成手段
(異常度合成装置)、5E…状態度記憶手段(プロセス異
常度記憶装置)、5G…学習装置、5GA…来歴状態度記憶
手段(来歴異常度記憶装置)、5GB…最適配分演算手段
(最適配分演算装置)、5GC…来歴演算データ記憶手段
(来歴伝播係数演算装置)、5GE…演算データ平均化手
段(伝播係数平均化装置)、5H…状態度誤差検出手段
(異常度誤差検出装置)。
FIG. 1 is a conceptual diagram showing the overall configuration of an inference apparatus to which the present invention is applied, FIG. 2 is a block diagram showing the detailed configuration of the inference processing section of the inference apparatus shown in FIG. 1, and FIG. 3 is membership. Fig. 4 is a graph showing an example of a function, Fig. 4 is a flow chart conceptually showing how event anomalies spread and propagate, Fig. 5 is a flow chart showing inference processing procedures, and Figs. 6A to 6C are event abnormal spreads. -Conceptual diagram showing the state of propagation and prediction, Fig. 7 is a block diagram showing the configuration of the embodiment of the present invention having a learning function, Fig. 8 is a fault tree diagram, and Fig. 9 is an abnormality due to fault tree FIG. 10 is a flow chart showing the spread prediction, and FIG. 10 is a block diagram showing a detailed configuration example of the learning device.
The figure is a flowchart showing the calculation procedure of the learning function. 3 ... Plant data storage means (plant database), 4 ... Knowledge base, 5 ... Inference processing device, 5A ... State degree calculation means (process abnormality degree calculation device), 5C ... State degree influence inference means (abnormality influence influence inference device) , 5D ... State degree synthesis means (abnormality degree synthesis apparatus), 5E ... State degree storage means (process abnormality degree storage apparatus), 5G ... Learning apparatus, 5GA ... History state degree storage means (history abnormality degree storage apparatus), 5GB ... Optimal distribution calculation means (optimal distribution calculation device), 5GC ... History calculation data storage means (history propagation coefficient calculation device), 5GE ... Calculation data averaging means (propagation coefficient averaging device), 5H ... State degree error detection means (abnormal) Error detector).

Claims (6)

【特許請求の範囲】[Claims] 【請求項1】プラント事象の状態量を測定し、測定され
た状態量が当該プラントの他の事象の状態量に及ぼす波
及効果を推論する方法において、前記プラント事象の状
態量の測定が時系列的に連続して行われ、該状態量の変
化が、該状態量の変化よりも時間的にあとの前記他の事
象の状態量に及ぼす波及効果の推論が、予め定められた
演算のためのデータに基づいて連続して行なわれること
を特徴とする推論方法。
1. A method of measuring a state quantity of a plant event and inferring a ripple effect of the measured state quantity on a state quantity of another event of the plant, wherein the state quantity of the plant event is measured in time series. For the predetermined calculation, the inference of the ripple effect that the change of the state quantity exerts on the state quantity of the other event temporally after the change of the state quantity is performed for a predetermined calculation. An inference method characterized by being performed continuously based on data.
【請求項2】プラント事象が他の事象に波及する経路上
の事象の状態量が連続して測定され、推論された前記プ
ラント事象が他の事象に及ぼす波及効果が、前記波及経
路上の事象の状態量に基づいて補正されることを特徴と
する請求項1に記載の推論方法。
2. The state quantity of an event on a path where a plant event propagates to another event is continuously measured, and the inferred ripple effect that the plant event has on another event is the event on the ripple path. The inference method according to claim 1, wherein the inference method is corrected based on the state quantity.
【請求項3】推論された波及効果に基づいてプラント事
象の状態をしめす数値が推論され、測定されたプラント
状態量から算出されたプラント事象の状態をしめす数値
と前記推論された数値との差を減少させるように予め定
められた演算のためのデータが変更される手順を備えて
いることを特徴とする請求項1もしくは2に記載の推論
方法。
3. A numerical value indicating the state of a plant event is deduced based on the inferred ripple effect, and the difference between the numerical value indicating the state of the plant event calculated from the measured plant state quantity and the inferred numerical value. 3. The inference method according to claim 1, further comprising a step of changing data for a predetermined operation so as to reduce
【請求項4】プラントの下位事象の状態が上位事象に及
ぼす効果を推論する装置において、プラント事象の状態
量を連続して記憶するプラントデータ記憶手段と、該記
憶手段に接続して設けられ前記状態量を事象の状態度を
しめす数値に変換する状態度演算手段と、該状態度演算
手段に接続して設けられ前記下位事象の状態度をしめす
数値に基づいて該下位事象の状態が上位事象に及ぼす波
及効果を状態度をしめす数値として出力する状態度波及
推論手段と、前記状態度演算手段に接続して設けられ状
態度をしめす数値を記憶する状態度記憶手段と、該状態
度記憶手段と前記状態度波及推論手段とに接続して設け
られ状態度波及推論手段が出力する数値を状態度記憶手
段が記憶している数値に基づいて補正する状態度合成手
段と、を備えていることを特徴とする推論装置。
4. An apparatus for inferring the effect of a lower event state of a plant on an upper event, comprising plant data storage means for continuously storing the state quantity of a plant event, and a plant data storage means connected to the storage means. A state degree calculating means for converting a state quantity into a numerical value indicating the degree of state of an event, and a state of the lower event based on a numerical value indicating the degree of state of the lower order event, which is provided in connection with the state degree calculating means. And a state degree memory means for storing the numerical value indicating the state degree, which is connected to the state degree calculating means, and outputs the influence effect on the state degree as a numerical value indicating the degree of state. And a state degree synthesizing means which is provided in connection with the state degree influence inference means and corrects the numerical value output by the state degree influence inference means based on the numerical value stored in the state degree storage means. Inference and wherein the.
【請求項5】状態度演算手段および状態度波及推論手段
に接続して設けられ前記各手段に演算のためのデータを
供給する知識ベースと、状態度演算手段および状態度波
及推論手段に接続して設けられ状態度演算手段が出力す
る状態度の数値と状態度波及推論手段が出力する状態度
の数値との差を演算する状態度誤差検出手段と、該誤差
検出手段および前記知識ベースに接続して設けられ前記
差を減少させるように前記知識ベースに記憶された演算
のためのデータを書き換える学習装置と、を備えている
ことを特徴とする請求項4に記載の推論装置。
5. A knowledge base, which is provided in connection with the state degree calculation means and the state degree spread inference means, and supplies data for calculation to each of the means, and is connected to the state degree calculation means and the state degree spread inference means. State degree error detecting means for calculating the difference between the state degree numerical value output by the state degree calculating means and the state degree numerical value output by the state degree influence inference means, and connected to the error detecting means and the knowledge base. And a learning device that rewrites data for arithmetic operations stored in the knowledge base so as to reduce the difference, and the inference device according to claim 4.
【請求項6】学習装置が、状態度誤差検出手段に接続さ
れプラント事象の状態度を時系列に記憶する来歴状態度
記憶手段と、該来歴状態度記憶手段に接続され誤差検出
手段が検出する差が零になる演算のためのデータを算出
する最適配分演算手段と、該最適配分演算手段に接続し
て設けられ演算のためのデータを時系列に記憶する来歴
演算データ記憶手段と、該来歴演算データ記憶手段およ
び前記最適配分演算手段に接続して設けられ最適配分演
算手段が出力する演算のためのデータと来歴演算データ
記憶手段が出力する演算のためのデータとから平均化さ
れた演算のためのデータを知識ベースに出力する演算デ
ータ平均化手段と、を備えていることを特徴とする請求
項5に記載の推論装置。
6. A learning device is connected to the degree-of-state error detecting means and stores the degree of state of a plant event in time series, and the error detecting means is connected to the history-state degree storing means to detect it. Optimal distribution calculation means for calculating data for calculation in which the difference becomes zero, history calculation data storage means connected to the optimum distribution calculation means for storing data for calculation in time series, and the history Of the calculation averaged from the data for calculation output by the optimum distribution calculation means provided connected to the calculation data storage means and the optimum distribution calculation means and the data for calculation output by the history calculation data storage means 6. The inference apparatus according to claim 5, further comprising: an arithmetic data averaging unit that outputs the data for the calculation to a knowledge base.
JP1004526A 1989-01-11 1989-01-11 Inference method and device Expired - Fee Related JPH07113995B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP1004526A JPH07113995B2 (en) 1989-01-11 1989-01-11 Inference method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP1004526A JPH07113995B2 (en) 1989-01-11 1989-01-11 Inference method and device

Publications (2)

Publication Number Publication Date
JPH02183895A JPH02183895A (en) 1990-07-18
JPH07113995B2 true JPH07113995B2 (en) 1995-12-06

Family

ID=11586492

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1004526A Expired - Fee Related JPH07113995B2 (en) 1989-01-11 1989-01-11 Inference method and device

Country Status (1)

Country Link
JP (1) JPH07113995B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10095813B2 (en) * 2013-11-18 2018-10-09 The Boeing Company Safety analysis of a complex system using component-oriented fault trees

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5966715A (en) * 1982-10-08 1984-04-16 Hitachi Ltd Plant failure spread prediction method
JPS6190210A (en) * 1984-10-11 1986-05-08 Nippon Atom Ind Group Co Ltd Plant abnormality diagnosis device
JPS63314609A (en) * 1987-06-17 1988-12-22 Mitsubishi Electric Corp Check system for equipment affected by trouble
JPH073642B2 (en) * 1987-06-19 1995-01-18 山武ハネウエル株式会社 Process abnormality diagnosis method

Also Published As

Publication number Publication date
JPH02183895A (en) 1990-07-18

Similar Documents

Publication Publication Date Title
US4630189A (en) System for determining abnormal plant operation based on whiteness indexes
EP0030459B2 (en) System for monitoring steam condenser performance
JP6658250B2 (en) Error diagnosis method, error diagnosis device, and error diagnosis program
JP5164954B2 (en) Device diagnostic method and device diagnostic device
CN101676559B (en) Method of alarm mask generation and condition monitoring of wind turbines
US20120296606A1 (en) Method, computer program, and system for performing interpolation on sensor data for high system availability
JP6702297B2 (en) Abnormal state diagnosis method and abnormal state diagnosis device
JPH03263227A (en) Knowledge data base processing system and expert system
JP2020047078A (en) Data processing method, data processing device, and data processing program
KR102735790B1 (en) A system for measuring and controlling water level of tank and a method thereof
JPH11161327A (en) Process abnormality diagnosis method and apparatus
JP7355108B2 (en) Prediction method, prediction device, recording medium
JPH07113995B2 (en) Inference method and device
JPH10301619A (en) Fuel cell power plant diagnostic equipment
JP2650914B2 (en) Process abnormality diagnosis device
JP2702206B2 (en) Operating status monitor for hydraulic machinery
JPH1115517A (en) Plant condition monitoring system
JPH07181292A (en) Plant operation support device
JP7850617B2 (en) Equipment deterioration diagnostic device and method
JP2540727B2 (en) Plant parameter abnormal amount search method and apparatus
JP6719612B1 (en) Outlier factor estimation support device, outlier factor estimation support method, and outlier factor estimation support program
KR20250040892A (en) Diagnostic device, diagnostic system, and diagnostic method for process processing equipment
WO2026070761A1 (en) Computer program, information processing method, and information processing device
JPH06103354B2 (en) Nuclear power plant diagnostics
JPS5952792A (en) plant monitoring equipment

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071206

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081206

Year of fee payment: 13

LAPS Cancellation because of no payment of annual fees