JPH11110367A - Automatic evaluation structure/parameter adjusting device for total state evaluation of plant - Google Patents

Automatic evaluation structure/parameter adjusting device for total state evaluation of plant

Info

Publication number
JPH11110367A
JPH11110367A JP9266198A JP26619897A JPH11110367A JP H11110367 A JPH11110367 A JP H11110367A JP 9266198 A JP9266198 A JP 9266198A JP 26619897 A JP26619897 A JP 26619897A JP H11110367 A JPH11110367 A JP H11110367A
Authority
JP
Japan
Prior art keywords
evaluation
fuzzy
state
parameter
data
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.)
Withdrawn
Application number
JP9266198A
Other languages
Japanese (ja)
Inventor
Yoshihiro Toriyama
佳洋 鳥山
Mitsuhiro Watanabe
光宏 渡辺
Makoto Uekusa
誠 植草
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.)
Fuji Electric Co Ltd
FFC Ltd
Original Assignee
Fuji Electric Co Ltd
FFC 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 Fuji Electric Co Ltd, FFC Ltd filed Critical Fuji Electric Co Ltd
Priority to JP9266198A priority Critical patent/JPH11110367A/en
Publication of JPH11110367A publication Critical patent/JPH11110367A/en
Withdrawn legal-status Critical Current

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

PROBLEM TO BE SOLVED: To reduce the burden of a plant operator unskilled in the setting of a fuzzy proposition by selecting and extracting an explanation parameter based on a significant explanation parameter index. SOLUTION: This device is provided with a 1st-order adjusting system (fuzzy proposition/membership function/fuzzy measure determining part) 10, 2nd-order adjusting system (G.A. type evaluating structure/parameter adjusting part) 20 and 3rd-order adjusting system (adjustment content collecting part) 30. In the 1st-order adjusting system 10, a target parameter and its total evaluation system are applied to a time sequential data base (DB) 101 of set plant data, and time sequential signal processing such as the arithmetic processing of a move average value/tendency value (change amount)} 102 is performed by a separately set signal processing parameter 111. Then, data to perform total state evaluation (evaluation object data) are separated by a threshold value corresponding to the target parameter and by considering only these separated evaluation object data, the explanation parameter is determined according to the significant explanation parameter index not to be affected by the relation between the explanation parameters in all the data.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、プラントの状態総
合評価における評価構造・パラメータの自動調整装置に
関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a device for automatically adjusting an evaluation structure and parameters in comprehensive evaluation of a plant state.

【0002】[0002]

【従来の技術】従来、プラントの状態総合評価に当た
り、プラントからの計測値情報はそのままの形で監視す
るのではなく、無駄時間の考慮、平滑化、傾向値の計算
等の信号処理を施してからプラントの状態評価を行って
いる。また、信号処理後のプラントデータを評価する方
法として、メンバシップ関数/ファジィ命題によるプラ
ント状態の言語的評価、つまり、「温度が高い」、「流
量値が上昇中」等の言語による状態評価が行われてい
る。ファジィ命題は認識すべきプラントの状態を説明す
る説明変数の挙動に関する知見であり、時系列データを
知見ごとの適合度に変換する働きを持つ。ファジィ命題
は、通常、複数の説明変数とそれに対するメンバシップ
関数により構成されている。更に、複数のファジィ命題
により言語的評価に対する適合度に変換されたプラント
状態は、ファジィ命題間の相互作用を考慮したファジィ
測度・積分により、複数の評価項目間の相互作用(相乗
的・相殺的な関係)を表現し、一つのプラント状態総合
評価値として統合した後、プラントのオペレータに通知
している。
2. Description of the Related Art Conventionally, in comprehensive evaluation of a plant state, measured value information from a plant is not monitored as it is, but signal processing such as consideration of dead time, smoothing, and calculation of a tendency value is performed. Is evaluating the state of the plant. As a method of evaluating the plant data after the signal processing, a linguistic evaluation of a plant state by a membership function / fuzzy proposition, that is, a state evaluation by a language such as “high temperature” or “flow rate value is increasing” is used. Is being done. The fuzzy proposition is knowledge about the behavior of an explanatory variable explaining the state of a plant to be recognized, and has a function of converting time-series data into a degree of fitness for each knowledge. A fuzzy proposition usually consists of a plurality of explanatory variables and their membership functions. Furthermore, the plant state converted into the degree of conformity to the linguistic evaluation by a plurality of fuzzy propositions is analyzed by a fuzzy measure / integration in consideration of the interaction between the fuzzy propositions, and the interaction between the plurality of evaluation items (synergistic and offsetting). Are integrated as one comprehensive evaluation value of the state of the plant, and then notified to the plant operator.

【0003】従来の技術では、無駄時間のシフト・平滑
化・傾向値の計算等の信号処理パラメータを、相互相関
分析による最適パラメータの決定による固定的な設定方
法により、あるいは、その後、プラントの状態評価精度
が向上するように、手作業で調整していた。また、説明
変数の選択及びファジィ命題の決定は、プラントの運転
員が直接、またはプラントの運転員から聞き取った情報
に基づいてエンジニアが行い、その都度、ファジィ測度
を自動的に再同定してプラントの状態総合評価の善し悪
しを評価し、評価精度向上のための調整作業を手作業に
より行っていた。
In the prior art, signal processing parameters such as shift of dead time, smoothing, calculation of a tendency value, and the like are fixed by a fixed setting method by determining an optimum parameter by cross-correlation analysis, or thereafter, the state of the plant is changed. Manual adjustments were made to improve the evaluation accuracy. In addition, the selection of explanatory variables and the determination of fuzzy propositions are performed directly by the plant operator or by an engineer based on information received from the plant operator, and each time, the fuzzy measure is automatically re-identified and the plant is re-identified. In order to evaluate the quality of the comprehensive evaluation of the condition, the adjustment work was performed manually to improve the evaluation accuracy.

【0004】[0004]

【発明が解決しようとする課題】従来においては、ファ
ジィ命題を人間が設定しているため、ファジィ命題の設
定に不慣れなプラント運転員にとっては大きな負担とな
っていた。また、プラントの挙動が複雑な場合には人間
が適切な知見を獲得することができず、ファジィ命題の
設定そのものが困難な場合もあった。一方、ファジィ命
題の設定を、目的変数と説明変数との間の相関関係に基
づいて有意な説明変数を選択し、その後、それらの説明
変数に対する挙動を捉えるメンバシップ関数を設定する
ことにより行う方式が考えられるが、この方式では、プ
ラントのように計測値が多い場合に、それらの変数が絡
み合って目的変数と各説明変数との間の明確な相関関係
を得ることが困難であり、定性的な特性を得ることさえ
危ぶまれる。まして定量的な特性を得ることはできな
い。しかし、異常状態のように変動の大きいデータに対
しては、定性的な関係を見いだすのは困難ではないた
め、本発明はこれに着目したものである。
Conventionally, human beings have set fuzzy propositions, which has been a heavy burden on plant operators who are not accustomed to setting fuzzy propositions. Also, when the behavior of the plant is complicated, humans cannot acquire appropriate knowledge, and setting the fuzzy proposition itself may be difficult. On the other hand, the fuzzy proposition is set by selecting significant explanatory variables based on the correlation between the objective variable and the explanatory variables, and then setting a membership function that captures the behavior for those explanatory variables. However, in this method, when there are many measured values such as in a plant, it is difficult to obtain a clear correlation between the objective variable and each explanatory variable due to the intertwining of those variables, and qualitatively You are even jeopardized to get the right properties. Furthermore, it is not possible to obtain quantitative characteristics. However, it is not difficult to find a qualitative relationship for data having a large fluctuation such as an abnormal state, and the present invention focuses on this.

【0005】そこで本発明の第1の課題は、設定した信
号処理パラメータのもとに信号処理加工を施した時系列
データに対して、多くの計測値が互いに絡み合うために
相関分析では明確な関係を見出せない場合でも、プラン
トの状態を説明する説明変数を適切に抽出し、その変数
の挙動を説明対象となるプラントの状態への適合度情報
に変換するファジィ命題・メンバシップ関数を決定し、
それらによる時系列データのファジィ命題適合度情報と
総合評価するべきプラント状態とに基づいてファジィ測
度を同定し、自動的にプラントの状態総合評価のための
評価構造・パラメータを決定して状態総合評価精度を計
算することにある。また、プラントの状態総合評価にお
いて、対象システムの複雑さから最適な信号処理パラメ
ータの設定を行うことはオペレータにとって一般に困難
であるほか、プロセスの挙動特性が時間的に変化するな
ど信号処理パラメータを挙動特性が変化するたびに調整
しなくてはならず、定周期で連続して評価構造・パラメ
ータ調整を行わざるを得ない運用系では、連続した時間
帯の時系列データによりそのたびに関連するプラントの
状態総合評価構造・パラメータの自動調整を行わなくて
はならない。更に、信号処理後の時系列データから説明
変数(ファジィ命題)を選択する際にも、精度向上のた
めの最適な組合せを選択する必要がある。また、信号処
理パラメータ・ファジィ命題の最適なものを選択するの
に、それらの組合せをすべて探索することは、計算量の
点から見て現実的ではない。そこで、本発明の第2の課
題は、有意義な時間内で探索可能な信号処理パラメータ
探索手段を付加することで、最適な信号処理パラメータ
と評価構造・パラメータを探索・調整し、プロセスの時
間的変化にも対応させることにある。
[0005] Therefore, a first object of the present invention is to provide a clear relationship in correlation analysis because many measured values are entangled with each other for time-series data subjected to signal processing based on set signal processing parameters. Even if it is not possible to find, properly extract explanatory variables that explain the state of the plant, determine a fuzzy proposition / membership function that converts the behavior of that variable into information on the degree of conformity to the state of the plant to be explained,
A fuzzy measure is identified based on the fuzzy proposition conformance information of the time series data and the plant state to be comprehensively evaluated, and the evaluation structure and parameters for the comprehensive state evaluation of the plant are automatically determined to comprehensively evaluate the state. Calculate accuracy. In addition, it is generally difficult for the operator to set the optimal signal processing parameters due to the complexity of the target system in the comprehensive evaluation of the plant status, and the behavior of the signal processing parameters, such as the process characteristics changing over time, is difficult. In operational systems, where adjustments must be made every time the characteristics change and evaluation structures and parameters must be adjusted continuously at regular intervals, the time series data for continuous time zones The automatic adjustment of the state comprehensive evaluation structure and parameters must be performed. Further, when selecting an explanatory variable (fuzzy proposition) from the time-series data after signal processing, it is necessary to select an optimal combination for improving accuracy. Further, it is not realistic to search all combinations of the signal processing parameters and fuzzy propositions in order to select the optimum ones from the viewpoint of the amount of calculation. Therefore, a second object of the present invention is to search and adjust an optimal signal processing parameter and an evaluation structure / parameter by adding a signal processing parameter search means capable of searching within a meaningful time, and to improve the time of the process. The idea is to adapt to changes.

【0006】[0006]

【課題を解決するための手段】上記課題を解決するた
め、請求項1記載の発明は、プラントの時系列データに
対し、目的変数とその総合評価方式を与え、別途設定し
た信号処理パラメータによる時系列的な信号処理を行っ
て、プラントの総合状態総合評価を与えるようなファジ
ィ命題・メンバシップ関数・ファジィ測度を評価構造・
パラメータとして同定すると共に、ファジィ命題ごとの
適合度演算、ファジィ積分演算、ファジィ積分値による
状態総合評価、状態総合評価精度計算を行う1次調整系
と、複数の時系列データを用い、個体としての前記信号
処理パラメータ及び評価構造・パラメータを逐次調整す
る2次調整系と、2次調整系により逐次調整された評価
構造・パラメータの中から評価精度の高い評価構造・パ
ラメータを生成する3次調整系と、を有する評価構造・
パラメータ自動調整装置において、前記1次調整系は、
プラントの時系列データに対し、設定した信号処理パラ
メータにより移動平均値及び変化量を演算する信号処理
手段と、目的変数に対するしきい値を用いて状態総合評
価対象データを分離する手段と、説明変数ごとに、異常
時データと全データの平均の差と、全データの分散とを
計算し、それらの比率を異常時データが全データから乖
離している度合いを示す説明変数有意指標として計算す
る手段と、この説明変数有意指標の大きい変数から順に
説明変数を選択する手段と、を備えたものである。すな
わち、請求項1記載の発明は、プラントの時系列データ
に対する信号処理手段と、目的変数に対するしきい値に
より状態総合評価を行うデータ(評価対象データ)を分
離し、この分離した評価対象データのみを考慮すること
により、全データにおける説明変数間の関係に左右され
ない説明変数有意指標に従って説明変数を決定する手段
を備え、前記説明変数有意指標により説明変数を選択、
抽出する説明変数決定手段に特徴を有する。
In order to solve the above-mentioned problems, the invention according to claim 1 provides an objective variable and an overall evaluation method thereof to a time-series data of a plant, and uses a signal processing parameter separately set. Performs sequential signal processing and evaluates fuzzy propositions, membership functions, and fuzzy measures that give a comprehensive evaluation of the overall state of the plant.
A primary adjustment system that identifies the parameters as well as performs the fitness calculation for each fuzzy proposition, fuzzy integration calculation, comprehensive evaluation of the state using the fuzzy integration value, and comprehensive calculation of the state evaluation accuracy, and a plurality of time-series data, A secondary adjustment system for sequentially adjusting the signal processing parameters and the evaluation structure / parameters; and a tertiary adjustment system for generating an evaluation structure / parameter with high evaluation accuracy from the evaluation structures / parameters sequentially adjusted by the secondary adjustment system. And an evaluation structure having
In the automatic parameter adjustment device, the primary adjustment system includes:
A signal processing means for calculating a moving average value and a change amount according to a set signal processing parameter with respect to a time series data of a plant; a means for separating state comprehensive evaluation target data using a threshold value for an objective variable; and an explanatory variable Means for calculating the difference between the mean of abnormal data and the average of all data and the variance of all data, and calculating the ratio thereof as an explanatory variable significant index indicating the degree of deviation of abnormal data from all data And means for selecting an explanatory variable in order from the variable having the largest explanatory variable significance index. That is, the invention according to claim 1 separates the signal processing means for the time series data of the plant and the data (evaluation target data) for performing the state comprehensive evaluation by the threshold value for the objective variable, and separates only the separated evaluation target data. In consideration of the above, comprising means for determining an explanatory variable according to an explanatory variable significance index not affected by the relationship between explanatory variables in all data, selecting an explanatory variable by the explanatory variable significance index,
The extraction variable determining means to be extracted is characterized.

【0007】請求項2記載の発明は、請求項1の発明に
おいて選択された説明変数に対し、説明変数有意指標の
符号からメンバシップ関数を決定する手段と、全データ
の平均値と最大値または最小値とを用いて、異常時デー
タをカバーし、異常時データの全データからの乖離が増
すにつれてメンバシップ関数の値が増加するようにメン
バシップ関数パラメータを決定する手段と、選択された
説明変数と前記メンバシップ関数とからファジィ命題を
決定する手段と、を備えたものである。すなわち、請求
項2記載の発明は、選択された説明変数に対し、プラン
トから得られる計測値データを評価対象状態への適合度
に変換するメンバシップ関数パラメータを、説明変数ご
との有意指標と統計量とに基づいて決定する手段によ
り、複数のファジィ命題・メンバシップ関数を自動的に
決定するファジィ命題・メンバシップ関数決定手段に特
徴を有する。
According to a second aspect of the present invention, for the explanatory variable selected in the first aspect of the present invention, means for determining a membership function from the sign of the explanatory variable significance index, an average value and a maximum value of all data or Means for using the minimum value to cover the abnormal data and determining the membership function parameters so that the value of the membership function increases as the deviation of the abnormal data from all data increases, and a selected explanation. Means for determining a fuzzy proposition from a variable and the membership function. In other words, according to the second aspect of the present invention, for the selected explanatory variable, a membership function parameter for converting measured value data obtained from the plant into a degree of conformity to the state to be evaluated is defined as a significance index and a statistical index for each explanatory variable. The fuzzy proposition / membership function determining means for automatically determining a plurality of fuzzy proposition / membership functions is characterized by the means for determining based on the quantity.

【0008】請求項3記載の発明は、請求項2により決
定されたメンバシップ関数及びファジィ命題を用い、プ
ラントの計測値データを評価対象状態への適合度に変換
する手段と、ファジィ測度積分により複数のファジィ命
題の適合度を一つのファジィ積分値に変換する手段と、
ファジィ命題の適合度と目的とするプラントの状態総合
評価値信号とによりファジィ測度を同定する手段と、前
記ファジィ積分値とこれに対するしきい値とからプラン
トの状態総合評価値信号を生成する手段と、プラントの
計測値データ内の実績値から状態総合評価の正誤判定を
行う手段と、正誤判定結果の中の正解数をカウントし、
状態総合評価精度を計算する手段と、を備えたものであ
る。すなわち、請求項3記載の発明は、請求項2により
決定されたメンバシップ関数及びファジィ命題を用いて
プラントの計測値データを評価対象状態への適合度に変
換する手段と、ファジィ測度積分により複数のファジィ
命題間の相互作用を考慮したファジィ命題適合度の統合
計算を行い、その値をプラントの状態総合評価値とし、
その値により単なる重み付け平均では評価・判定できな
いプロセス挙動の爆発的変化(相乗効果)や挙動が類似
した説明変数の評価レベルの抑制効果(相殺効果)を考
慮してプラントの状態総合評価を行う状態総合評価手段
と、ある時間帯のプラントデータによる状態総合評価成
功率を計算して評価構造及び評価パラメータを評価する
プラント状態総合評価精度計算手段とに特徴を有する。
According to a third aspect of the present invention, there is provided a means for converting measured value data of a plant into a degree of conformity to a state to be evaluated using the membership function and the fuzzy proposition determined in the second aspect, and a fuzzy measure integration. Means for converting the fitness of a plurality of fuzzy propositions into one fuzzy integral value;
Means for identifying a fuzzy measure based on the degree of conformity of the fuzzy proposition and the target state evaluation value signal of the plant; means for generating a state evaluation value signal of the plant from the fuzzy integration value and a threshold value for the fuzzy integration value; Means for performing correct / incorrect judgment of the state comprehensive evaluation from the actual value in the measured value data of the plant, and counting the number of correct answers in the correct / incorrect judgment result,
Means for calculating the state total evaluation accuracy. That is, a third aspect of the present invention provides a means for converting measured value data of a plant into a degree of conformity to a state to be evaluated using the membership function and the fuzzy proposition determined by the second aspect, and a plurality of means by fuzzy measure integration. Computes the fuzzy proposition goodness of fit taking into account the interaction between fuzzy propositions, and uses that value as the overall state evaluation value of the plant,
A state in which the plant state is comprehensively evaluated taking into account the explosive change in process behavior (synergistic effect) that cannot be evaluated and judged by simple weighted averaging based on that value (synergistic effect) and the effect of suppressing the evaluation level of the explanatory variables with similar behavior (cancelling effect). It is characterized by comprehensive evaluation means and plant state comprehensive evaluation accuracy calculating means for calculating a state comprehensive evaluation success rate based on plant data in a certain time zone to evaluate an evaluation structure and evaluation parameters.

【0009】請求項4記載の発明は、請求項1における
信号処理パラメータの組を乱数により発生する手段と、
請求項1〜請求項3の手段を用いて、発生した信号処理
パラメータによる信号処理の結果、発生するそれぞれ異
なった時系列データに対するファジィ命題・メンバシッ
プ関数・ファジィ測度(前述のようにこの三つの組を評
価構造・パラメータと呼び、これに信号処理パラメータ
を加えたものを個体と呼ぶことにする)の決定及び状態
総合評価精度の計算を行う手段と、設定したしきい値よ
り評価精度が低い個体を削除する手段と、残った個体に
対し、設定したしきい値より評価精度の低い個体に対し
てそのファジィ測度からファジィ命題間の相互作用係数
を求める手段と、相互作用係数が最小であるファジィ命
題を、請求項1における説明変数有意指標の次候補によ
るファジィ命題に置き換える手段(そのファジィ命題を
置換したファジィ命題と呼ぶ)と、削除した個体に対
し、信号処理パラメータの組を乱数により新規生成する
手段(そのファジィ命題を突然変異したファジィ命題と
呼ぶ)と、置換したファジィ命題及び突然変異したファ
ジィ命題を用いてファジィ測度を再同定し、プラントの
状態総合評価精度を再計算する手段と、請求項1〜請求
項3の手段を実施する際に用いた時系列データ(調整デ
ータ)を含み、それよりも設定数倍の時間帯の時系列デ
ータ(評価用データ)を設定する手段と、評価用データ
を用いて個体の状態評価精度を計算し、評価精度の向上
率が極小になったとき、または設定回数だけ再評価を繰
り返したときにこの反復計算を終了し、その時に評価精
度が最大である個体を評価構造・パラメータとして保存
する手段とを備えたものである。すなわち、請求項4記
載の発明は、状態判定成功率を向上させるための最適化
対象である信号処理パラメータ及びファジィ命題の最適
な組合せを探索するに当たり、前記請求項1〜3の発明
による一連の動作と、状態総合評価精度に関して最適な
信号処理パラメータ及びファジィ命題の組合せを遺伝的
アルゴリズム(GA)型方式により求める手段と、ファ
ジィ測度から計算するファジィ命題適合度間の相互作用
係数により類似した挙動を示す説明変数を削除し、状態
総合評価に有効な説明変数・ファジィ命題の選択頻度を
増大させ、最適な信号処理パラメータ・ファジィ命題の
相互作用係数によるGA型の信号処理パラメータ・ファ
ジィ命題自動調整手段に特徴を有する。
According to a fourth aspect of the present invention, there is provided means for generating a set of signal processing parameters according to the first aspect by random numbers,
As a result of the signal processing using the generated signal processing parameters, the fuzzy proposition, the membership function, and the fuzzy measure (as described above) for different time series data generated by using the means of claim 1 to claim 3 A set is called an evaluation structure / parameter, and the sum of the signal processing parameters and the signal processing parameters is called an individual), and a means for determining the overall evaluation accuracy of the state, and an evaluation accuracy lower than the set threshold value A means for deleting an individual, a means for obtaining an interaction coefficient between fuzzy propositions from a fuzzy measure for an individual having an evaluation accuracy lower than a set threshold value for a remaining individual, and a minimum interaction coefficient Means for replacing the fuzzy proposition with a fuzzy proposition by the next candidate of the explanatory variable significance index in claim 1 (fuzzy proposition with the fuzzy proposition replaced) And a means for newly generating a set of signal processing parameters by random numbers for the deleted individual (the fuzzy proposition is called a mutated fuzzy proposition). Means for re-identifying the fuzzy measure using the method and recalculating the overall state evaluation accuracy of the plant; and time-series data (adjustment data) used in implementing the means of claim 1 to claim 3. Means for setting time-series data (evaluation data) for several times the time zone, and calculating the state evaluation accuracy of the individual using the evaluation data, and when the rate of improvement in the evaluation accuracy is minimal, or When the re-evaluation is repeated the set number of times, the repetition calculation is terminated, and at this time, the individual having the highest evaluation accuracy is stored as an evaluation structure / parameter. That is, the invention according to claim 4 searches for an optimal combination of a signal processing parameter and a fuzzy proposition to be optimized to improve the success rate of state determination, and a series of the invention according to claims 1 to 3 described above. Means for finding the optimal combination of signal processing parameters and fuzzy propositions with respect to the operation and state evaluation accuracy by a genetic algorithm (GA) type method, and behavior similar to the interaction coefficient between the fuzzy proposition fitness calculated from the fuzzy measure Automatically adjusts GA-type signal processing parameters and fuzzy propositions by increasing the frequency of selection of explanatory variables and fuzzy propositions that are effective for comprehensive evaluation of states, and using optimal signal processing parameters and fuzzy proposition interaction coefficients The means are characterized.

【0010】請求項5記載の発明は、請求項4における
個体の削除処理において、前回調整時に保存した評価構
造・パラメータ(前回最大評価精度の個体)は削除しな
いように動作を修正する手段と、今回調整により、評価
精度が前回よりも設定率以上、向上した場合には削除す
るように動作を修正する手段とを備えたものである。す
なわち、請求項5記載の発明は、請求項4により複数時
刻帯の時系列データを逐次的に調整する過程において、
前回調整時に保存した信号処理パラメータ・ファジィ命
題を今回調整時において探索対象から除外しないように
する手段と、今回調整により状態総合評価精度が前回よ
りも設定率以上の向上がある場合には、前回調整時に保
存した信号処理パラメータ・ファジィ命題を除外する手
段とを有する信号処理パラメータ・ファジィ命題忘却
(過調整)防止手段に特徴を有する。
According to a fifth aspect of the present invention, in the individual deleting process according to the fourth aspect, means for correcting an operation so as not to delete the evaluation structure / parameter (an individual having the previous maximum evaluation accuracy) stored at the time of the previous adjustment; A means is provided for correcting the operation so that when the evaluation accuracy is improved by more than the set rate compared to the previous time, the operation is deleted. That is, according to a fifth aspect of the present invention, in the process of sequentially adjusting the time-series data of a plurality of time zones according to the fourth aspect,
Means for preventing the signal processing parameter / fuzzy proposition saved during the previous adjustment from being excluded from the search target at the time of this adjustment. A signal processing parameter / fuzzy proposition forgetting (over-adjustment) preventing means having a means for excluding the signal processing parameter / fuzzy proposition stored during the adjustment.

【0011】請求項6記載の発明は、請求項5に記載し
たプラントの状態総合評価における評価構造・パラメー
タ自動調整装置において、前記3次調整系は、複数の時
系列データにより評価構造・パラメータを連続的に調整
したときの各評価構造・パラメータを保存する手段と、
各評価構造・パラメータごとに各時系列データに対する
状態総合評価精度を計算する手段と、各時系列データに
対する状態総合評価精度の平均値を計算する手段と、そ
の平均値が最大である評価構造・パラメータを保存する
手段と、各時系列データごとに各評価構造・パラメータ
に対する状態総合評価精度を計算する手段と、各評価構
造・パラメータに対する状態総合評価精度の平均値を計
算する手段と、その平均値が最小である時系列データを
調整したときに得られた評価構造・パラメータを選択す
る手段と、選択した評価構造・パラメータから、ファジ
ィ測度の相互作用係数が最大であるファジィ命題を追加
ファジィ命題として保存する手段と、状態総合評価精度
の平均値が最大である評価構造・パラメータ中のファジ
ィ命題としてのマスタファジィ命題に前記追加ファジィ
命題を追加して新たな合成ファジィ命題を生成する手段
と、全ての時系列データに対して、前記マスタファジィ
命題及び合成ファジィ命題の状態総合評価精度を計算す
る手段と、合成ファジィ命題の状態総合評価精度がマス
タファジィ命題の状態総合評価精度よりも大きい場合
に、合成ファジィ命題をマスタファジィ命題とする手段
と、を備えたものである。すなわち、請求項6記載の発
明は、請求項4,5による状態総合評価構造・パラメー
タ(ファジィ命題・メンバシップ関数・ファジィ測度)
及び信号処理パラメータの逐次調整終了時において、各
時刻帯の時系列データによる評価構造・パラメータ調整
結果を保存する手段と、逐次調整を行った個々のデータ
に対し保存した全ての評価構造・パラメータにより再評
価(状態総合評価精度の計算)する手段と、全てのデー
タの状態総合評価精度が最大である評価構造・パラメー
タをマスタ評価構造・パラメータとして選択する手段
と、各時刻帯の時系列データの中で最も状態総合評価精
度が低い評価構造・パラメータをマスタ評価構造・パラ
メータに合成する手段とを備え、連続調整時の調整デー
タごとの過調整を防止し、全てのデータに対して状態総
合評価精度が高い評価構造・パラメータを出力する評価
構造・パラメータ自動調整結果再評価手段に特徴を有す
る。
According to a sixth aspect of the present invention, in the apparatus for automatically adjusting the evaluation structure and parameters in the comprehensive evaluation of the state of the plant according to the fifth aspect, the tertiary adjustment system converts the evaluation structure and parameters using a plurality of time-series data. Means for storing each evaluation structure / parameter when continuously adjusted,
Means for calculating the state comprehensive evaluation accuracy for each time-series data for each evaluation structure / parameter; means for calculating the average value of the state comprehensive evaluation accuracy for each time-series data; Means for storing parameters, means for calculating the state comprehensive evaluation accuracy for each evaluation structure / parameter for each time-series data, means for calculating the average value of the state comprehensive evaluation accuracy for each evaluation structure / parameter, and the average A means for selecting the evaluation structure and parameters obtained when adjusting the time-series data with the minimum value, and adding a fuzzy proposition with the largest interaction coefficient of the fuzzy measure from the selected evaluation structure and parameters As a fuzzy proposition in the evaluation structure / parameter with the maximum average value of the state comprehensive evaluation accuracy. Means for generating a new combined fuzzy proposition by adding the additional fuzzy proposition to the tafuzzy proposition, means for calculating the state comprehensive evaluation accuracy of the master fuzzy proposition and the combined fuzzy proposition for all time series data, Means for making the composite fuzzy proposition a master fuzzy proposition when the state comprehensive evaluation accuracy of the composite fuzzy proposition is greater than the state comprehensive evaluation accuracy of the master fuzzy proposition. That is, the invention according to claim 6 provides the state comprehensive evaluation structure and parameters (fuzzy proposition, membership function, fuzzy measure) according to claims 4 and 5.
Means for storing the evaluation structure / parameter adjustment result based on the time-series data of each time zone at the end of the sequential adjustment of the signal processing parameters, and all the evaluation structures / parameters stored for the individually adjusted individual data. Means for re-evaluation (calculation of state total evaluation accuracy), means for selecting an evaluation structure / parameter having the maximum state total evaluation accuracy of all data as a master evaluation structure / parameter, and a means for selecting time-series data of each time zone. A means for synthesizing the evaluation structure / parameter with the lowest state comprehensive evaluation accuracy to the master evaluation structure / parameter is provided to prevent over-adjustment for each adjustment data at the time of continuous adjustment, and comprehensive state evaluation for all data It is characterized by an evaluation structure / parameter automatic adjustment result re-evaluation means for outputting a highly accurate evaluation structure / parameter.

【0012】本発明は以下の作用をなす。プラントの状
態総合評価において、 (1)プラントの時系列データが、相関分析により全体
として目的変数に対する説明変数の抽出が困難な場合に
おいても、異常時データのみを考慮する説明変数有意指
標を求めることで、目的変数が異常時のように大きく変
動している場合に、各説明変数に対し、異常時データが
大きく乖離している場合を数値化し、プラントの異常を
説明する説明変数を自動的に抽出することができる。
The present invention has the following functions. In the comprehensive evaluation of the state of the plant, (1) To obtain an explanatory variable significance index that takes into account only abnormal data, even when it is difficult to extract the explanatory variable for the objective variable as a whole by correlation analysis of the time-series data of the plant. Then, when the target variable fluctuates greatly as in the case of an abnormality, for each explanatory variable, the case where the abnormal time data greatly deviates is quantified, and the explanatory variable explaining the plant abnormality is automatically set. Can be extracted.

【0013】(2)上記(1)により決定した説明変数
に対して、異常時データの乖離の度合いを表す有意指標
の極性を用いることで、メンバシップ関数の形を異常時
データと全体データとの乖離が増大するときに適合度も
増大する方向に決定することができる。また、説明変数
の平均値と最大値または最小値を用いることで、メンバ
シップ関数のパラメータを異常時データと全データとの
乖離がなくなる変数値の適合度を0に、異常時データと
全データとの乖離が最大になる変数値の適合度を1にす
るように、自動的に決定することができる。
(2) The polarity of the significance index indicating the degree of deviation of the abnormal data is used for the explanatory variable determined in the above (1), so that the form of the membership function can be changed between the abnormal data and the overall data. Can be determined in such a direction that the degree of matching also increases when the deviation of In addition, by using the average value and the maximum value or the minimum value of the explanatory variables, the parameter of the membership function is set to 0 so that the discrepancy between the abnormal data and all data is eliminated. Can be automatically determined so that the degree of conformity of the variable value that maximizes the difference between the two is set to 1.

【0014】(3)上記(1),(2)により決定した
ファジィ命題・メンバシップ関数により、ファジィ測度
を自動的に同定し、プラントの状態総合評価構造・パラ
メータを決定することができ、これらの評価構造・パラ
メータを用いて、説明変数へのファジィ命題間の相互作
用を考慮したプラントの状態総合評価を行い、更に正誤
判定を行うことで、プラントの状態総合評価精度の判定
を自動的に行うことができる。
(3) The fuzzy measure can be automatically identified by the fuzzy proposition and membership function determined by the above (1) and (2), and the overall state evaluation structure and parameters of the plant can be determined. By using the evaluation structure and parameters of the above, comprehensive evaluation of the state of the plant taking into account the interaction between fuzzy propositions to the explanatory variables is performed, and correctness / incorrectness is determined, thereby automatically determining the accuracy of the comprehensive evaluation of the state of the plant. It can be carried out.

【0015】(4)プラントの状態総合評価精度を向上
させるために、時系列データに対する信号処理パラメー
タの最適な組合せを、GA型探索により行うことで、効
率よく探索することができる。具体的には、まず、複数
の信号処理パラメータの組をランダムに発生させ、それ
らによる状態総合評価構造を決定し、評価精度を計算す
る。次に、それらの評価精度の大きさにより、それぞれ
の信号処理パラメータ及びファジィ命題、ファジィ測度
を削除/置換/突然変異といった改良作業を施して、新
たな個体に対する評価構造・パラメータの調整・評価精
度の計算を反復して実行する。また、ファジィ測度から
計算するファジィ命題適合度間の相互作用係数により、
探索に不必要な候補を探索対象から除外することで探索
効率を更に向上させることができる。
(4) In order to improve the overall evaluation accuracy of the state of the plant, the optimal combination of the signal processing parameters for the time series data can be efficiently searched by performing the GA type search. Specifically, first, a plurality of sets of signal processing parameters are randomly generated, a state comprehensive evaluation structure based on them is determined, and evaluation accuracy is calculated. Next, depending on the magnitude of the evaluation accuracy, the respective signal processing parameters, fuzzy propositions, and fuzzy measures are subjected to improvement work such as deletion / replacement / mutation to adjust the evaluation structure / parameter adjustment / evaluation accuracy for a new individual. Is repeatedly performed. Also, by the interaction coefficient between the fuzzy proposition fitness calculated from the fuzzy measure,
The search efficiency can be further improved by excluding the unnecessary candidates from the search target.

【0016】(5)ある時系列データによる、プラント
の状態総合評価構造・パラメータの逐次調整時の評価構
造・パラメータ探索過程において、前回調整時における
最大評価精度の評価構造・パラメータに対しては選択頻
度を制御することで、時系列データの変化に対して過調
整(オーバーフィッティング)を抑制することができ、
安定した評価構造・パラメータを獲得することができ
る。
(5) Based on certain time series data, in the overall state evaluation structure of the plant, the evaluation structure at the time of successive adjustment of parameters, and the parameter search process, the evaluation structure and parameters of the maximum evaluation accuracy at the previous adjustment are selected. By controlling the frequency, it is possible to suppress over-adjustment (over-fitting) for changes in time-series data,
A stable evaluation structure and parameters can be obtained.

【0017】(6)ある期間の複数時刻データによる評
価構造・パラメータ逐次調整結果を解析することで、最
も信頼性のある評価構造・パラメータ(マスタファジィ
命題)を選択することができる。更に、上記解析の結
果、最も調整困難な時系列データよりファジィ命題(追
加ファジィ命題)を抽出してマスタファジィ命題に追加
合成することができ、一層評価精度の高い評価構造・パ
ラメータを獲得することができる。
(6) The most reliable evaluation structure / parameter (master fuzzy proposition) can be selected by analyzing the evaluation structure / parameter sequential adjustment result based on a plurality of time data in a certain period. Furthermore, as a result of the above analysis, fuzzy propositions (additional fuzzy propositions) can be extracted from the time series data that is most difficult to adjust and can be additionally synthesized with the master fuzzy proposition, thereby obtaining evaluation structures and parameters with higher evaluation accuracy. Can be.

【0018】以上により、プラントの状態総合評価を行
う最適な状態総合評価構造・パラメータを自動的に決定
することができ、プラント特性の時間的変化に動的に対
応することが可能になる。
As described above, the optimum state comprehensive evaluation structure and parameters for performing the comprehensive state evaluation of the plant can be automatically determined, and it is possible to dynamically respond to the temporal change of the plant characteristics.

【0019】[0019]

【発明の実施の形態】以下、図に沿って本発明の実施形
態を説明する。ここでは、プラントの状態総合評価にお
ける、あるプラント状態を予測する評価構造・パラメー
タの自動調整の例(第1実施形態)と、プロセスのある
操作量(設定値)をオペレータに代わって出力する際の
評価構造・パラメータの自動調整の例(第2実施形態)
を想定して説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of the present invention will be described below with reference to the drawings. Here, an example of the automatic adjustment of the evaluation structure / parameter for predicting a certain plant state in the comprehensive evaluation of the state of the plant (first embodiment) and a case where a certain manipulated variable (set value) of the process is output on behalf of the operator Of automatic adjustment of evaluation structure / parameters of the second example (second embodiment)
It is assumed that

【0020】まず、図1は第1、第2実施形態の双方に
共通する本発明の全体構成を示す概念図であり、本発明
の評価構造・パラメータ自動調整装置は、1次調整系
(ファジィ命題・メンバシップ関数・ファジィ測度決定
部)10、2次調整系(G.A.型評価構造・パラメー
タ調整部)20及び3次調整系(調整内容まとめ部)3
0を備えている。1次調整系10では、設定したプラン
トデータの時系列DB(データベース)101に対し、
目的変数とその総合評価方式を与え、別途設定した信号
処理パラメータ111による時系列信号処理(移動平均
値・傾向値(変化量)の演算処理等)102を行い、プ
ラントの状態総合評価に用いるファジィ命題及びメンバ
シップ関数112、ファジィ測度113を同定して状態
総合評価構造・パラメータを獲得し、ファジィ命題ごと
の適合度演算処理114、ファジィ積分演算処理11
5、積分値による状態総合評価処理116、状態総合評
価精度計算処理117を行う。なお、ファジィ命題及び
メンバシップ関数112、ファジィ測度113は評価構
造・パラメータDB118に格納され、更に、この評価
構造・パラメータDB118と信号処理パラメータ11
1とが個体DB110に格納される。ここで、説明変数
決定処理103、メンバシップ関数パラメータ決定処理
104、ファジィ測度同定処理105の詳細については
後述する。
First, FIG. 1 is a conceptual diagram showing the overall configuration of the present invention which is common to both the first and second embodiments. The evaluation structure / parameter automatic adjustment apparatus of the present invention comprises a primary adjustment system (fuzzy system). Proposition, membership function, fuzzy measure determination unit) 10, secondary adjustment system (GA type evaluation structure / parameter adjustment unit) 20, and tertiary adjustment system (adjustment content summarization unit) 3
0 is provided. In the primary adjustment system 10, a time series DB (database) 101 of the set plant data is
An objective variable and its comprehensive evaluation method are given, time-series signal processing (operation processing of moving average value / trend value (change amount), etc.) 102 is performed using a separately set signal processing parameter 111, and fuzzy logic used for comprehensive evaluation of the state of the plant A proposition and membership function 112 and a fuzzy measure 113 are identified to obtain a state comprehensive evaluation structure and parameters, and a fitness calculation process 114 for each fuzzy proposition and a fuzzy integration calculation process 11
5. Perform state comprehensive evaluation processing 116 based on integral values and state comprehensive evaluation accuracy calculation processing 117. The fuzzy proposition and membership function 112 and the fuzzy measure 113 are stored in the evaluation structure / parameter DB 118, and the evaluation structure / parameter DB 118 and the signal processing parameter 11
1 is stored in the individual DB 110. Here, details of the explanatory variable determination processing 103, the membership function parameter determination processing 104, and the fuzzy measure identification processing 105 will be described later.

【0021】2次調整系20では、複数の時系列データ
を用い、1次調整系10の信号処理パラメータ111、
ファジィ命題・メンバシップ関数112の逐次調整を自
動的に行う。すなわちこの調整系20は、自動調整の開
始処理201、信号処理パラメータ111の組を乱数に
より複数組発生させる処理202、個体ごとの評価精度
の計算処理203、その計算の反復終了判定(停止)処
理204、計算した評価精度に基づく個体の淘汰処理
(パラメータの削除)205、評価精度が最大である個
体に対する前記淘汰の禁止処理206、1次調整系10
におけるファジィ命題の置換改良処理207、信号処理
パラメータの組の乱数発生による増殖処理208を有す
る。この調整系20におけるファジィ命題・メンバシッ
プ関数112の逐次調整には、GA型の最適パラメータ
探索方式を採用する。
The secondary adjustment system 20 uses a plurality of time-series data and uses the signal processing parameters 111 and
The sequential adjustment of the fuzzy proposition / membership function 112 is automatically performed. That is, the adjustment system 20 includes a process 201 for starting automatic adjustment, a process 202 for generating a plurality of sets of the signal processing parameters 111 using random numbers, a process 203 for calculating the evaluation accuracy for each individual, and a process for determining the end of the repetition (stop) of the calculation. 204, an individual selection process based on the calculated evaluation accuracy (deletion of parameters) 205, the selection prohibition process 206 for the individual having the highest evaluation accuracy 206, the primary adjustment system 10
And a multiplication process 208 for generating a random number of a set of signal processing parameters. The sequential adjustment of the fuzzy proposition / membership function 112 in the adjustment system 20 employs a GA type optimal parameter search method.

【0022】3次調整系30では、1次、2次調整系1
0,20を経て逐次調整した評価構造・パラメータか
ら、長期間の評価精度が最も大きいものをマスタファジ
ィ命題として抽出する。更に、マスタファジィ命題での
評価精度が低い時系列データに対し、その時系列データ
を調整した際に得られた評価構造・パラメータから追加
ファジィ命題(説明変数)を相乗効果最大命題として抽
出し、マスタファジィ命題と合成することにより、評価
精度の再計算、合成効果の確認の上で、全ての時系列デ
ータに対する評価精度の高い評価構造・パラメータを生
成する。すなわちこの調整系30は、調整まとめ開始処
理301、時系列DB101及び評価構造・パラメータ
DB118内のデータを取り込んで行う、過去に自動調
整した評価構造・パラメータの評価精度再計算処理30
2、マスタファジィ命題の抽出処理305、過調整構造
の抽出及び相乗効果最大命題(追加ファジィ命題)の抽
出処理306、マスタファジィ命題と相乗効果最大命題
との合成処理307、評価精度再計算及び合成効果の確
認処理308及び停止処理309を有する。
In the tertiary adjustment system 30, the primary and secondary adjustment systems 1
From the evaluation structures / parameters successively adjusted through 0 and 20, those having the longest evaluation accuracy are extracted as master fuzzy propositions. Furthermore, for time-series data having a low evaluation accuracy in the master fuzzy proposition, additional fuzzy propositions (explanatory variables) are extracted as synergistic maximum propositions from evaluation structures and parameters obtained when the time-series data is adjusted. By synthesizing with the fuzzy proposition, evaluation structures and parameters with high evaluation accuracy for all time-series data are generated after recalculating the evaluation accuracy and confirming the synthesis effect. That is, the adjustment system 30 performs the adjustment summary start processing 301, the evaluation structure / parameter evaluation accuracy recalculation processing 30 of the evaluation structure / parameters that have been automatically adjusted in the past, and performs the acquisition of data in the time series DB 101 and the evaluation structure / parameter DB 118.
2. Master fuzzy proposition extraction processing 305, over-adjustment structure extraction and synergy effect maximum proposition (additional fuzzy proposition) extraction processing 306, master fuzzy proposition and synergy effect maximum proposition synthesis processing 307, evaluation accuracy recalculation and synthesis An effect confirmation process 308 and a stop process 309 are provided.

【0023】次に、図2に基づいて、説明変数決定処理
103、メンバシップ関数パラメータ決定処理104の
全体を説明する。まず、プラントのセンサ情報を取り込
み、数式1に示すような一定量(T時刻分)のN変数の
情報をデータ1として、記憶装置に一旦、蓄積する。こ
こで、xvtは時刻tにおけるv番目の変数の値を表
す。また、xOtは時刻tにおける目的変数の値を表
す。
Next, the whole of the explanatory variable determination processing 103 and the membership function parameter determination processing 104 will be described with reference to FIG. First, sensor information of a plant is fetched, and information of a certain amount (T times) of N variables as shown in Expression 1 is temporarily stored as data 1 in a storage device. Here, xv t represents the value of the v-th variable at time t. XO t represents the value of the objective variable at time t.

【0024】[0024]

【数1】 (Equation 1)

【0025】次に、蓄積したデータに対して移動平均・
変化量算出の信号処理を行なう。信号処理により高周波
成分を取り除き、トレンドデータの傾向値を算出する。
具体的に、データ1に対し、各変数の時間方向に移動平
均処理(一般的であるので説明を省略する)を施したも
のと、同じく変化量算出(一定時間の変数値の変化量算
出)を行なったものとを並べたものを、数式2に示す如
く新たにデータAとおく。ここで、Nは、元の値の2倍
になっていることになる。また、目的変数xOtは移動
平均処理した変数のみとし、変化量を演算した目的変数
は一般のxvtと同じ扱いとする。
Next, a moving average and
The signal processing for calculating the amount of change is performed. High frequency components are removed by signal processing, and the trend value of the trend data is calculated.
Specifically, the data 1 is subjected to moving average processing (generally omitted for simplicity) in the time direction of each variable, and the change amount is calculated (calculation of the change amount of the variable value in a certain time). The result obtained by performing the above is arranged as data A as shown in Expression 2. Here, N is twice the original value. Also, the target variable xO t is only the variable subjected to the moving average processing, and the target variable for which the amount of change is calculated is treated the same as general xv t .

【0026】[0026]

【数2】 (Equation 2)

【0027】次に、信号処理を施したデータを説明変数
決定処理103に転送し、捉えるべき異常兆候を説明す
る有意な変数y1,y2,…,yv,…,yMを決定す
る。そして、決定した説明変数をメンバシップ関数パラ
メータ決定処理104に転送し、説明変数の変動を異常
兆候の度合に変換する、ファジィ命題のメンバシップ関
数パラメータを決定する。以上により、数式3に示すよ
うなファジィ命題を決定する。
Next, the signal-processed data is transferred to an explanatory variable determining process 103, and significant variables y1, y2,..., Yv,. Then, the determined explanatory variables are transferred to the membership function parameter determination processing 104, and the membership function parameters of the fuzzy proposition that convert the fluctuation of the explanatory variables into the degree of the abnormal sign are determined. Thus, a fuzzy proposition as shown in Expression 3 is determined.

【0028】[0028]

【数3】 (Equation 3)

【0029】ここで、メンバシップ関数SA,LAの種
類別は、メンバシップ関数パラメータ決定処理104に
おいて決定する。なお、図3(a),(b)にメンバシ
ップ関数の種類とパラメータを示す。
Here, the types of the membership functions SA and LA are determined in the membership function parameter determination processing 104. FIGS. 3A and 3B show the types and parameters of the membership functions.

【0030】次いで、図2〜図9を参照しながら、説明
変数決定処理103、メンバシップ関数パラメータ決定
処理104の詳細を説明する。図2において、信号処理
を施したデータAの中で、目的変数xOtの値が、設定
したしきい値Thを越えるSステップ前時刻のデータ
を、異常発生直前データ(便宜上、異常時データと呼
ぶ)としてデータAから分離し、これをデータBとす
る。またその件数をUとする。次に、各変数xv毎に、
全時刻のデータAの統計量と分離した異常発生直前のデ
ータBの平均値(mA,mB)・標準偏差σAを数式4によ
りそれぞれ算出する。
Next, the details of the explanatory variable determining process 103 and the membership function parameter determining process 104 will be described with reference to FIGS. 2, in the data A which has been subjected to signal processing, the value of the objective variable xO t is, the data of the S step before time exceeding the threshold value Th set, before an error occurred data (for convenience, the abnormal data The data A is separated from the data A, and is referred to as data B. Also, the number of cases is U. Next, for each variable xv,
The statistic of the data A at all times and the average value (m A , m B ) and the standard deviation σ A of the data B immediately before the occurrence of the separated error are calculated by Equation 4.

【0031】[0031]

【数4】 (Equation 4)

【0032】次に、各変数について説明変数有意指標I
expを数式5により求め、各変数を説明変数有意指標の
絶対値の降順に並び替える。
Next, the explanatory variable significance index I for each variable
exp is calculated by Expression 5, and the variables are sorted in descending order of the absolute value of the explanatory variable significance index.

【0033】[0033]

【数5】Iexp=(mB−mA)/σA I exp = (m B −m A ) / σ A

【0034】ここで、目的変数自身は本処理の対象から
除く。そして、並べ替えた変数の上位M件の変数を説明
変数として決定し、これを、 y1,y2,…,yv,…,yM とする。補足として、数式5による説明変数有意指標I
expは全時刻のデータすなわちデータAの散らばりに対
する、異常発生直前時刻のデータであるデータBが分離
している比率を表す。
Here, the target variable itself is excluded from the target of this processing. Then, the top M variables of the rearranged variables are determined as explanatory variables, and are defined as y1, y2,..., Yv,. As a supplement, the explanatory variable significance index I according to Equation 5
exp represents the ratio of data B, which is data immediately before the occurrence of an abnormality, to the distribution of data at all times, that is, data A, being separated.

【0035】図3に、メンバシップ関数の形状とパラメ
ータ[p1,p2]を示す。すなわち、メンバシップ関
数は図に示す帰属度関数であり、入力の値に応じて
[0,1]の値を出力する関数とする。更に、メンバシ
ップ関数の形状は、形についてはファジィラベルLA及
びSAによって、関数の折れ曲がり点はパラメータ[p
1,p2]によってそれぞれ決定される。以上のような
関数をメンバシップ関数とする。
FIG. 3 shows the shape of the membership function and the parameters [p1, p2]. That is, the membership function is a membership function shown in the figure, and is a function that outputs a value of [0, 1] according to an input value. Further, the shape of the membership function is determined by fuzzy labels LA and SA, and the bending point of the function is determined by the parameter [p
1, p2]. The above functions are referred to as membership functions.

【0036】図4において、説明変数決定処理103に
て計算した説明変数有意指標Iexpの符号により、説明
変数毎のメンバシップ関数の種類を以下のように決定す
る: Iexp<0の場合 メンバシップ関数SA(左上がりの
関数) Iexp>0の場合 メンバシップ関数LA(右上がりの
関数) その上で、メンバシップ関数のパラメータ[p1,p
2]を以下のように決定する: ・説明変数が変化量でない場合 Iexp<0の場合 SA[minA,mA] (図5参照) Iexp>0の場合 LA[mA,maxA] (図6参照) ・説明変数が傾向値(変化量)である場合: Iexp<0の場合 SA[mB,0] (図7参照) Iexp>0の場合 LA[0,mB] (図8参照) ここで、minAはデータAの最小値、maxAはデータ
Aの最大値、minBはデータBの最小値、maxBはデ
ータBの最大値である。
In FIG. 4, the type of the membership function for each explanatory variable is determined as follows based on the sign of the explanatory variable significance index I exp calculated in the explanatory variable determination processing 103: If I exp <0: Member Ship function SA (left-up function) If I exp > 0 Membership function LA (right-up function) Then, the parameters [p1, p
2] is determined as follows: - If the explanatory variables is not variation I exp <For 0 SA [min A, m A ] ( see FIG. 5) I exp> For 0 LA [m A, max A (See FIG. 6) • When the explanatory variable is a trend value (change amount): If I exp <0, SA [m B , 0] (see FIG. 7) If I exp > 0, LA [0, m B (See FIG. 8) Here, min A is the minimum value of data A, max A is the maximum value of data A, min B is the minimum value of data B, and max B is the maximum value of data B.

【0037】図9は、ファジィ命題ごとの適合度演算処
理114の動作ブロック図である。上述のようにして決
定されたファジィ命題とメンバシップ関数を用いて、信
号処理済みのプラントの計測値データに対し、計算対象
変数を選択してメンバシップ関数の値を計算し、ファジ
ィ命題適合度(評価対象状態への適合度)に変換する。
FIG. 9 is an operation block diagram of the fitness calculation processing 114 for each fuzzy proposition. Using the fuzzy proposition and the membership function determined as described above, for the measured value data of the signal-processed plant, the calculation target variable is selected and the value of the membership function is calculated. (Degree of conformity to the state to be evaluated).

【0038】次に、図10に基づいてファジィ測度同定
処理105につき説明する。まず、ファジィ測度及びフ
ァジィ積分(Choquet積分)について、その定義を行な
う。本実施形態では、有限集合上のファジィ測度のみを
扱う。そこで、数式6のようにおく。
Next, the fuzzy measure identification processing 105 will be described with reference to FIG. First, the fuzzy measure and the fuzzy integral (Choquet integral) are defined. In this embodiment, only a fuzzy measure on a finite set is handled. Therefore, it is set as in Equation 6.

【0039】[0039]

【数6】X={1,2,…,n}X = {1, 2,..., N}

【0040】定義1:ファジィ測度 ファジィ測度を数式7のように定義する。なお、数式7
において、2XはXのべき集合である。
Definition 1: Fuzzy measure A fuzzy measure is defined as in equation (7). Equation 7
Where 2 X is the power set of X.

【0041】[0041]

【数7】 (Equation 7)

【0042】定義2:ファジィ積分 関数f:X→Rのファジィ測度μに関するChoquet積分
は、数式8により定義される。ただし、x1,x2,…,
nはXの要素をFの値について昇順に並べ替えたもの
であり、 f(x1)≦f(x2)≦…≦f(xn) である。また、数式9を条件とする。
Definition 2: Fuzzy integral The function f: Choquet integral for the fuzzy measure μ of X → R is defined by equation (8). Where x 1 , x 2 , ...,
x n are those which sorts the elements of X in ascending order the values of the F, is f (x 1) ≦ f ( x 2) ≦ ... ≦ f (x n). Equation 9 is a condition.

【0043】[0043]

【数8】 (Equation 8)

【0044】[0044]

【数9】 (Equation 9)

【0045】ファジィ測度同定処理105は、図10に
示す如く、ファジィ命題適合度を入力とし、それを上記
定義2のファジィ積分処理を行なった結果得られるファ
ジィ積分値と、時系列データ内の目的変数との残差を最
小にするように、ファジィ測度を逐次的に修正する。収
束判定は、上記残差と設定しきい値との比較により行な
い、残差がしきい値以下になったとき、及びファジィ測
度修正回数が設定回数に達したときに反復修正及び残差
評価による収束判定・同定を終了する。
As shown in FIG. 10, the fuzzy measure identification processing 105 receives a fuzzy integral value obtained as a result of performing the fuzzy integration processing of the above-described definition 2 with the fuzzy proposition matching degree as an input, and an objective value in the time series data. Modify the fuzzy measure sequentially to minimize residuals with variables. The convergence determination is performed by comparing the above residual with a set threshold, and when the residual becomes equal to or less than the threshold, and when the number of times of fuzzy measure correction reaches the set number, it is determined by iterative correction and residual evaluation. The convergence judgment / identification ends.

【0046】図11はファジィ積分演算処理115の動
作を示すものであり、このファジィ積分演算処理は、入
力である複数のファジィ命題(例えば命題1,命題2)
の適合度(命題ごとの適合度演算処理114による)を
その値の降順に並び替え、命題適合度の階差数列を求
め、ファジィ測度μとの内積計算によりファジィ積分値
を獲得する。図12はファジィ測度μの例である。この
ファジィ測度は集合に対する関数であり、命題1の適合
度のファジィ測度と命題2の適合度のファジィ測度との
和が、命題1の適合度と命題2の適合度のファジィ測度
よりも小さく、命題1と命題2とを同時に見た場合の相
乗的な効果を表している。
FIG. 11 shows the operation of the fuzzy integral operation processing 115. In this fuzzy integral operation processing, a plurality of fuzzy propositions (eg, proposition 1 and proposition 2) which are inputs are used.
Are sorted in descending order of their values (by the suitability calculation process 114 for each proposition), a sequence of differences of the proposition suitability is obtained, and a fuzzy integral value is obtained by calculating an inner product with the fuzzy measure μ. FIG. 12 is an example of the fuzzy measure μ. This fuzzy measure is a function for the set, and the sum of the fuzzy measure of the fitness of Proposition 1 and the fuzzy measure of the fitness of Proposition 2 is smaller than the fuzzy measure of the fitness of Proposition 1 and the fitness of Proposition 2. This shows a synergistic effect when viewing Proposition 1 and Proposition 2 at the same time.

【0047】図13はファジィ積分の例を示している。
ファジィ積分は、ファジィ測度による相互作用を考慮し
た線形和演算(面積計算)であり、図12のファジィ測
度によるファジィ積分において、命題1及び命題2の適
合度の大小関係によって使用されるファジィ測度が異な
り、ファジィ積分値が説明変数間の相互作用により決定
されることを表している。図14はファジィ積分値によ
る状態総合評価116及び状態総合評価精度計算117
の動作ブロック図である。ファジィ積分値による状態総
合評価116は、ファジィ積分値とファジィ積分値に対
するしきい値との比較により行ない、それを状態総合評
価信号とすることで行なう。更に、時系列データベース
の目的変数データ(計測値データ)を用いて、実際のプ
ラント状態と状態総合評価信号との比較を行ない、状態
総合評価の正誤判定を行なう。判定結果の中から正解数
をカウントし、状態総合評価精度としてデータベースに
登録する。
FIG. 13 shows an example of fuzzy integration.
The fuzzy integral is a linear sum operation (area calculation) in consideration of the interaction by the fuzzy measure. In the fuzzy integral by the fuzzy measure shown in FIG. Differently, it indicates that the fuzzy integral value is determined by the interaction between the explanatory variables. FIG. 14 shows a state comprehensive evaluation 116 and a state comprehensive evaluation accuracy calculation 117 based on fuzzy integral values.
3 is an operation block diagram of FIG. The state comprehensive evaluation 116 based on the fuzzy integrated value is performed by comparing the fuzzy integrated value with a threshold value for the fuzzy integrated value, and using the result as a state comprehensive evaluation signal. Further, using the target variable data (measured value data) of the time series database, the actual plant state is compared with the state comprehensive evaluation signal to determine whether the overall state evaluation is correct or not. The number of correct answers is counted from the judgment result and registered in the database as the state total evaluation accuracy.

【0048】図15は、GA型信号処理パラメータ・フ
ァジィ命題自動調整処理の構成を示すものであり、前述
した図1における1次調整系10及び2次調整系20の
一部を抜き出したものに相当する。信号処理パラメータ
・ファジィ命題の自動調整が開始(処理201)される
と、信号処理のための信号処理パラメータの組(個体の
一部)を乱数により複数発生させる(処理202)。次
に、発生した個体毎に1次調整系10を動作させ、評価
構造・パラメータ(ファジィ命題とメンバシップ関数1
12・ファジィ測度113)を求める。次に、それらの
状態総合評価精度を計算して状態総合評価精度DB11
7Aに蓄積し、個体毎の評価精度を求める(処理20
3)。次いで、評価精度の値により、評価構造・パラメ
ータの削除を行なう(処理205)。これを淘汰処理と
呼び、削除された個体を死滅した個体と呼ぶ。次に、淘
汰処理後に残った個体に対して、ファジィ測度によるフ
ァジィ命題置換改良処理207を行ない、評価精度向上
の見込まれるファジィ命題を構成する。次に、1次調整
系10を再度動作させて、変更した個体に対する評価精
度の再計算を行なう(処理203,204)。ここまで
の、一連の処理を設定回数繰り返す(この一連の処理を
世代交代、繰り返すことを進化と呼ぶ)。更に、淘汰処
理205で死滅した個体に対しては、信号処理パラメー
タの組合せ乱数発生により、新しい個体として1次調整
系動作・評価精度の再計算対象とする。これを突然変異
処理208と呼ぶ。
FIG. 15 shows a configuration of the GA type signal processing parameter / fuzzy proposition automatic adjustment process, which is obtained by extracting a part of the primary adjustment system 10 and the secondary adjustment system 20 in FIG. Equivalent to. When the automatic adjustment of the signal processing parameter / fuzzy proposition is started (processing 201), a plurality of sets of signal processing parameters (part of an individual) for signal processing are generated by random numbers (processing 202). Next, the primary adjustment system 10 is operated for each generated individual, and the evaluation structure and parameters (fuzzy proposition and membership function 1
12. Obtain the fuzzy measure 113). Next, the state comprehensive evaluation accuracy is calculated to calculate the state comprehensive evaluation accuracy DB11.
7A and obtain the evaluation accuracy for each individual (processing 20
3). Next, the evaluation structure / parameter is deleted according to the value of the evaluation accuracy (process 205). This is called selection processing, and the deleted individual is called a dead individual. Next, a fuzzy proposition replacement improvement process 207 based on a fuzzy measure is performed on the individuals remaining after the selection process, thereby constructing a fuzzy proposition that is expected to improve evaluation accuracy. Next, the primary adjustment system 10 is operated again to recalculate the evaluation accuracy for the changed individual (processes 203 and 204). The series of processing up to this point is repeated a set number of times (this series of processing is called generational change and repeated is called evolution). Furthermore, the individual who died in the selection process 205 is subjected to re-calculation of the primary adjustment system operation / evaluation accuracy as a new individual by generating a combination random number of signal processing parameters. This is called a mutation process 208.

【0049】図16は、上述した個体の淘汰処理205
の動作ブロック図である。個体ごとの状態総合評価精度
を降順に並び替え、評価精度淘汰しきい値により、しき
い値以下の個体を削除する。この結果を個体DB110
に再登録する。図17は、ファジィ命題置換改良処理2
07の動作ブロック図である。死滅していない生存して
いる個体から状態総合評価精度の低い個体を選択する。
次に、選択した個体に対応するファジィ測度からファジ
ィ積分の入力であるファジィ命題の適合度間の相互作用
係数(Shapley値)を計算し、相互作用係数が最小のファ
ジィ命題(説明変数)を選択する。ここで、Shapley値と
は、ファジィ測度μ(μは全体集合X={x1,x2
…,xn}に対する2X上の集合関数で、x1,x2,…,
nはそれぞれ、ファジィ命題適合度である)から求め
られる値、 φ(μ)={φ(μ)(x1),φ(μ)(x2),…
…,φ(μ)(xn)} のことであり、数式10により計算される。なお、数式
11を条件とし、|・|は集合の要素数である。
FIG. 16 shows the individual selection process 205 described above.
3 is an operation block diagram of FIG. The state total evaluation accuracy for each individual is sorted in descending order, and individuals below the threshold are deleted according to the evaluation accuracy selection threshold. This result is stored in the individual DB 110
Register again. FIG. 17 shows fuzzy proposition replacement improvement processing 2
It is an operation block diagram of 07. Individuals with low overall state evaluation accuracy are selected from living individuals who have not died.
Next, from the fuzzy measure corresponding to the selected individual, calculate the interaction coefficient (Shapley value) between the fitness of the fuzzy proposition which is the input of the fuzzy integral, and select the fuzzy proposition (explanatory variable) with the smallest interaction coefficient I do. Here, the Shapley value is a fuzzy measure μ (μ is a whole set X = {x 1 , x 2 ,
, X n } is a set function on 2 X , x 1 , x 2 ,.
x n respectively, the values obtained from the fuzzy proposition is fit), φ (μ) = { φ (μ) (x 1), φ (μ) (x 2), ...
.., Φ (μ) (x n )}, and is calculated by Expression 10. Note that | · | is the number of elements in the set, on condition of Expression 11.

【0050】[0050]

【数10】 (Equation 10)

【0051】[0051]

【数11】 [Equation 11]

【0052】次に、上述した相互作用係数が最小である
ファジィ命題(説明変数)を、説明変数決定手段(請求
項1)における説明変数有意指標の次候補によるファジ
ィ命題に置き換える。置き換えたファジィ命題を置換さ
れたファジィ命題(新ファジィ命題)として個体DB1
10(評価構造・パラメータDB118)に再登録す
る。
Next, the fuzzy proposition (explanatory variable) having the minimum interaction coefficient is replaced with a fuzzy proposition based on the next candidate of the explanatory variable significance index in the explanatory variable determining means (claim 1). Individual DB1 as a replaced fuzzy proposition (new fuzzy proposition) with the replaced fuzzy proposition
10 (evaluation structure / parameter DB 118).

【0053】図18は、個体毎の評価精度の計算処理2
03及び反復終了判定処理204の動作ブロック図であ
る。まず、個体毎の評価精度の計算は、1次調整系10
で自動調整に用いた時刻のものを含む設定倍数の時系列
データにより行なう。次に、2次調整系20の反復終了
判定処理204において、反復回数が調整回数の上限に
達したとき、または評価精度向上率が極小になったとき
に、自動調整を終了する。最後に、調整の結果評価精度
が最大である個体を準最適な個体として、個体DB11
0(評価構造・パラメータDB118)に登録する。
FIG. 18 shows a calculation process 2 of evaluation accuracy for each individual.
FIG. 3 is an operation block diagram of a process No. 03 and a repetition end determination process 204. First, the calculation of the evaluation accuracy for each individual is performed by the primary adjustment system 10.
Is performed using time series data of a set multiple including the time used for automatic adjustment. Next, in the repetition end determination process 204 of the secondary adjustment system 20, when the number of repetitions reaches the upper limit of the number of adjustments, or when the evaluation accuracy improvement rate is minimized, the automatic adjustment ends. Finally, the individual with the maximum evaluation accuracy as a result of the adjustment is regarded as a sub-optimal individual, and the individual DB 11
0 (evaluation structure / parameter DB 118).

【0054】図19は、ファジィ命題忘却防止処理(評
価精度最大個体に対する淘汰禁止処理206)の動作ブ
ロック図である。図において、前回調整時系列データに
よる状態総合評価精度の中で、評価精度が最大である個
体(準最適な個体)を選択する。次に、今回調整時系列
データによる個体に対する淘汰処理205において、上
記選択個体に対しては淘汰処理を行なわないように淘汰
禁止指令を加える。ただし、今回調整により、状態総合
評価精度が前回より設定率以上の向上がある場合には、
本淘汰禁止処理を行なわないものとする。
FIG. 19 is an operation block diagram of the fuzzy proposition forgetting prevention processing (selection prohibition processing 206 for individuals with the highest evaluation accuracy). In the figure, an individual (sub-optimal individual) with the highest evaluation accuracy is selected from among the state comprehensive evaluation accuracy based on the previously adjusted time-series data. Next, in the selection process 205 for individuals based on the current adjusted time-series data, a selection prohibition command is added so that the selection process is not performed for the selected individual. However, if the state total evaluation accuracy has improved by more than the set rate from the previous time due to this adjustment,
This selection prohibition processing is not performed.

【0055】図20は、3次調整系30における過去自
動調整評価構造・パラメータの評価精度再計算処理30
2、マスタファジィ命題の抽出処理305の動作ブロッ
ク図である。図において、過去調整した時系列DB10
1とそれに調整した評価構造・パラメータDB108と
の全ての組合せの状態総合評価精度を再評価し、これを
評価精度DB117Aに格納する。次に、評価構造・パ
ラメータ毎に、各時系列DBの状態総合評価精度の平均
値を計算する。次いで、状態総合評価精度の平均値が最
大である評価構造からファジィ命題をマスタファジィ命
題として選択し、メンバシップ関数と共に一時格納す
る。
FIG. 20 shows a past automatic adjustment evaluation structure / parameter evaluation accuracy recalculation processing 30 in the tertiary adjustment system 30.
2 is an operation block diagram of a master fuzzy proposition extraction process 305. FIG. In the figure, the time series DB 10 adjusted in the past
The state comprehensive evaluation accuracy of all combinations of 1 and the adjusted evaluation structure / parameter DB 108 is reevaluated and stored in the evaluation accuracy DB 117A. Next, the average value of the state comprehensive evaluation accuracy of each time-series DB is calculated for each evaluation structure / parameter. Next, a fuzzy proposition is selected as a master fuzzy proposition from the evaluation structure in which the average value of the state comprehensive evaluation accuracy is the maximum, and is temporarily stored together with the membership function.

【0056】図21は、3次調整系30における過調整
構造の抽出処理及び相乗効果が最大である命題の抽出処
理306の動作ブロック図である。評価精度DB117
Aから、時系列DB101毎に、各評価構造・パラメー
タの評価精度の平均値を計算する。次に、評価精度の平
均値が最小である評価構造を選択する。次いで、選択し
た評価構造から相乗効果(相互作用)が最大であるファ
ジィ命題を追加ファジィ命題として選択し、メンバシッ
プ関数と共に一時格納する。
FIG. 21 is an operation block diagram of the extraction processing of the over-adjusted structure and the extraction processing 306 of the proposition with the maximum synergistic effect in the tertiary adjustment system 30. Evaluation accuracy DB117
From A, the average value of the evaluation accuracy of each evaluation structure / parameter is calculated for each time-series DB 101. Next, an evaluation structure with the minimum evaluation accuracy average value is selected. Next, a fuzzy proposition having the greatest synergy (interaction) is selected as an additional fuzzy proposition from the selected evaluation structure, and is temporarily stored together with the membership function.

【0057】図22は、あるプラント状態を予測する評
価構造・パラメータの自動調整の例である第1実施形態
の動作ブロック図である。図において、時系列データか
ら、異常予測対象となる目的変数を選択する。その際、
設定されたしきい値より上か下かによって目的変数に対
する異常状態の定義を行なう。次に、全ての説明変数に
対する、ラグタイム、平滑化・傾向値計算時間幅の候補
を設定し、これらを信号処理パラメータ候補とする。最
後に、目的変数と異常状態の定義、及び信号処理パラメ
ータを1〜3次調整系10〜30に転送し、前述した手
順により状態総合評価構造・パラメータの自動調整を行
なう。
FIG. 22 is an operation block diagram of the first embodiment which is an example of the automatic adjustment of the evaluation structure / parameter for predicting a certain plant state. In the figure, an objective variable to be subjected to abnormality prediction is selected from time-series data. that time,
The abnormal state of the objective variable is defined depending on whether the value is above or below the set threshold value. Next, lag time and smoothing / trend value calculation time width candidates are set for all explanatory variables, and these are set as signal processing parameter candidates. Finally, the definition of the objective variable and the abnormal state, and the signal processing parameters are transferred to the first to third order adjustment systems 10 to 30, and the state comprehensive evaluation structure and parameters are automatically adjusted according to the above-described procedure.

【0058】図23は、プロセスのある操作量(設定
値)をオペレータに代わって出力する際の評価構造・パ
ラメータの自動調整の例である第2実施形態の動作ブロ
ック図である。図において、時系列データから、全ての
説明変数に対する、ラグタイム、平滑化・傾向値計算時
間幅の候補を設定し、これらを信号処理パラメータ候補
とする。次に、オペレータ等によるプロセス操作量を状
態総合評価対象とし、この操作量を目的変数として時系
列データに追加する。最後に、オペレータ操作量及び信
号処理パラメータを1〜3次調整系10〜30に転送
し、前述した手順により状態総合評価構造・パラメータ
の自動調整を行なう。
FIG. 23 is an operation block diagram of the second embodiment which is an example of the automatic adjustment of the evaluation structure / parameter when outputting a certain operation amount (set value) of the process on behalf of the operator. In the figure, lag time and smoothing / trend value calculation time width candidates are set for all explanatory variables from time-series data, and these are set as signal processing parameter candidates. Next, a process operation amount by an operator or the like is set as a target for state comprehensive evaluation, and this operation amount is added to the time-series data as an objective variable. Finally, the operator's operation amount and signal processing parameters are transferred to the first to third order adjustment systems 10 to 30, and the state total evaluation structure and parameters are automatically adjusted according to the above-described procedure.

【0059】[0059]

【発明の効果】以上述べたように、本発明によれば、 (1)目的変数と各説明変数との間に非線形な特性があ
る場合でも、総合評価する状態を説明する説明変数を選
び出すことができる。 (2)プラントの状態総合評価を行なう際に、これまで
運転員の経験による知識をヒヤリング形式で獲得してい
た説明変数を、全時系列変数から自動的に抽出すること
ができる。 (3)プラント状態を説明する説明変数の挙動を捉え、
[0,1]の適合度に変換するファジィ命題のメンバシ
ップ関数の決定、及びファジィ積分によるプラントの状
態総合評価と評価精度の計算を自動的に行なうことがで
き、これまでプラントの運転員ないしその知識を獲得し
たエンジニアが行なっていた、プラント状態に対する知
見の抽出作業に伴う負荷を軽減することができる。 (4)信号処理パラメータの最適な組合せの探索をファ
ジィ命題間の相互作用を考慮したGA型選択により行な
うことで、効率良く実行できる。 (5)連続した複数時刻における評価構造・パラメータ
の逐次調整において、前回調整した評価構造・パラメー
タを考慮することで個々の時刻における過調整を抑制で
き、安定した評価構造・パラメータを獲得することがで
きる。 (6)ある一定期間の連続した逐次調整後、調整内容を
解析することで長い期間において評価精度が高く、有効
なマスタ評価構造・パラメータを獲得し、それに局所過
調整状態に陥った評価構造・パラメータを追加合成する
ことで、長期間における状態総合評価精度を更に向上さ
せることができる。
As described above, according to the present invention, (1) even when there is a non-linear characteristic between the objective variable and each explanatory variable, an explanatory variable explaining the state to be comprehensively evaluated is selected. Can be. (2) When performing a comprehensive evaluation of the state of a plant, it is possible to automatically extract, from all time-series variables, explanatory variables for which knowledge based on the experience of the operator has been acquired in the form of a hearing. (3) Capture the behavior of the explanatory variables that explain the plant state,
It is possible to automatically determine the membership function of the fuzzy proposition to be converted to the fitness of [0, 1], and automatically evaluate the overall state of the plant and calculate the evaluation accuracy by fuzzy integration. It is possible to reduce the load involved in extracting knowledge about the plant state, which was performed by an engineer who has acquired the knowledge. (4) The search for the optimal combination of the signal processing parameters is efficiently performed by performing the GA type selection in consideration of the interaction between the fuzzy propositions. (5) In the successive adjustment of the evaluation structure / parameters at a plurality of consecutive times, over-adjustment at each time can be suppressed by considering the evaluation structure / parameter adjusted last time, and a stable evaluation structure / parameter can be obtained. it can. (6) After continuous adjustment for a certain period of time, by analyzing the adjustment contents, the evaluation accuracy is high over a long period of time, an effective master evaluation structure / parameter is obtained, and the evaluation structure that falls into a local over-adjustment state is obtained. By additionally synthesizing the parameters, it is possible to further improve the overall state evaluation accuracy over a long period of time.

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

【図1】本発明の全体構成を示す概念図である。FIG. 1 is a conceptual diagram showing the overall configuration of the present invention.

【図2】説明変数決定処理を示す動作ブロック図であ
る。
FIG. 2 is an operation block diagram illustrating an explanatory variable determination process.

【図3】メンバシップ関数の種類とパラメータを示す図
である。
FIG. 3 is a diagram showing types and parameters of membership functions.

【図4】メンバシップ関数パラメータ決定処理の動作ブ
ロック図である。
FIG. 4 is an operation block diagram of a membership function parameter determination process.

【図5】メンバシップ関数の説明図である。FIG. 5 is an explanatory diagram of a membership function.

【図6】メンバシップ関数の説明図である。FIG. 6 is an explanatory diagram of a membership function.

【図7】メンバシップ関数の説明図である。FIG. 7 is an explanatory diagram of a membership function.

【図8】メンバシップ関数の説明図である。FIG. 8 is an explanatory diagram of a membership function.

【図9】ファジィ命題ごとの適合度演算処理の動作ブロ
ック図である。
FIG. 9 is an operation block diagram of a fitness calculation process for each fuzzy proposition.

【図10】ファジィ測度同定処理の動作ブロック図であ
る。
FIG. 10 is an operation block diagram of a fuzzy measure identification process.

【図11】ファジィ積分演算処理の動作ブロック図であ
る。
FIG. 11 is an operation block diagram of a fuzzy integration operation process.

【図12】ファジィ測度の例を示す図である。FIG. 12 is a diagram illustrating an example of a fuzzy measure.

【図13】ファジィ積分の例を示す図である。FIG. 13 is a diagram illustrating an example of fuzzy integration.

【図14】ファジィ積分による状態総合評価処理及び評
価精度計算処理の動作ブロック図である。
FIG. 14 is an operation block diagram of state comprehensive evaluation processing and evaluation accuracy calculation processing by fuzzy integration.

【図15】GA型信号処理パラメータ・ファジィ命題自
動調整手段の構成を示す図である。
FIG. 15 is a diagram showing a configuration of GA type signal processing parameter / fuzzy proposition automatic adjustment means.

【図16】個体の淘汰処理の動作ブロック図である。FIG. 16 is an operation block diagram of an individual selection process.

【図17】ファジィ命題置換改良処理の動作ブロック図
である。
FIG. 17 is an operation block diagram of a fuzzy proposition replacement improvement process.

【図18】個体ごとの評価精度の計算処理及び反復終了
判定処理の動作ブロック図である。
FIG. 18 is an operation block diagram of a calculation process of evaluation accuracy for each individual and a repetition end determination process.

【図19】評価精度最大個体に対する淘汰禁止処理の動
作ブロック図である。
FIG. 19 is an operation block diagram of a selection prohibition process for an individual with the highest evaluation accuracy.

【図20】過去自動調整評価構造・パラメータの評価精
度再計算処理及びマスタファジィ命題の抽出処理の動作
ブロック図である。
FIG. 20 is an operation block diagram of past automatic adjustment evaluation structure / parameter evaluation accuracy recalculation processing and master fuzzy proposition extraction processing.

【図21】過調整構造の抽出処理及び相乗効果が最大で
ある命題の抽出処理の動作ブロック図である。
FIG. 21 is an operation block diagram of the extraction processing of the over-adjusted structure and the extraction processing of the proposition with the maximum synergistic effect.

【図22】本発明の第1実施形態の動作ブロック図であ
る。
FIG. 22 is an operation block diagram of the first embodiment of the present invention.

【図23】本発明の第2実施形態の動作ブロック図であ
る。
FIG. 23 is an operation block diagram of the second embodiment of the present invention.

【符号の説明】[Explanation of symbols]

10 1次調整系 20 2次調整系 30 3次調整系 101 時系列DB 102 信号処理 103 説明変数決定処理 104 メンバシップ関数パラメータ決定処理 105 ファジィ測度同定処理 110 個体DB 111 信号処理パラメータ 112 ファジィ命題・メンバシップ関数 113 ファジィ測度 114 ファジィ命題ごとの適合度演算処理 115 ファジィ積分演算処理 116 積分値による状態総合評価処理 117 状態総合評価精度計算処理 117A 状態総合評価精度DB 118 評価構造・パラメータDB 201 自動調整開始処理 202 信号処理パラメータの組を乱数により複数組発
生させる処理 203 個体ごとの評価精度の計算処理 204 反復終了判定(停止)処理 205 評価精度に基づく個体の淘汰処理(パラメータ
の削除) 206 評価精度が最大である個体に対する淘汰禁止処
理 207 ファジィ命題の置換改良処理 208 信号処理パラメータの組の乱数発生による増殖
処理 301 調整まとめ開始処理 302 過去に自動調整した評価構造・パラメータの評
価精度再計算処理 305 マスタファジィ命題の抽出処理 306 過調整構造の抽出及び相乗効果最大命題の抽出
処理 307 マスタファジィ命題と相乗効果最大命題との合
成処理 308 評価精度再計算及び合成効果の確認処理 309 停止処理
DESCRIPTION OF SYMBOLS 10 Primary adjustment system 20 Secondary adjustment system 30 Tertiary adjustment system 101 Time series DB 102 Signal processing 103 Explanation variable determination processing 104 Membership function parameter determination processing 105 Fuzzy measure identification processing 110 Individual DB 111 Signal processing parameters 112 Fuzzy proposition Membership function 113 Fuzzy measure 114 Goodness-of-fit calculation processing for each fuzzy proposition 115 Fuzzy integration calculation processing 116 Comprehensive state evaluation processing based on integral values 117 Comprehensive state evaluation accuracy calculation processing 117A Comprehensive state evaluation accuracy DB 118 Evaluation structure / parameter DB 201 Automatic adjustment Start processing 202 Processing of generating a plurality of sets of signal processing parameters by random numbers 203 Calculation processing of evaluation accuracy for each individual 204 Repetition end determination (stop) processing 205 Selection processing of individual based on evaluation accuracy (parameter deletion) 206) Selection prohibition processing for individuals with the highest evaluation accuracy 207 Fuzzy proposition replacement improvement processing 208 Propagation processing by generating random numbers of signal processing parameter sets 301 Adjustment start processing 302 Evaluation accuracy of evaluation structures and parameters automatically adjusted in the past Recalculation processing 305 Extraction processing of master fuzzy proposition 306 Extraction of over-adjusted structure and extraction processing of maximum synergistic effect proposition 307 Synthesis processing of master fuzzy proposition and maximum synergy effect proposition 308 Evaluation accuracy recalculation and confirmation processing of synthesis effect 309 Stop processing

───────────────────────────────────────────────────── フロントページの続き (72)発明者 植草 誠 東京都日野市富士町1番地 富士ファコム 制御株式会社内 ──────────────────────────────────────────────────続 き Continued on the front page (72) Inventor Makoto Uekusa 1 Fuji-cho, Hino-shi, Tokyo Inside Fujifacom Control Co., Ltd.

Claims (6)

【特許請求の範囲】[Claims] 【請求項1】 プラントの時系列データに対し、目的変
数とその総合評価方式を与え、別途設定した信号処理パ
ラメータによる時系列的な信号処理を行って、プラント
の総合状態総合評価を与えるようなファジィ命題・メン
バシップ関数・ファジィ測度を評価構造・パラメータと
して同定すると共に、ファジィ命題ごとの適合度演算、
ファジィ積分演算、ファジィ積分値による状態総合評
価、状態総合評価精度計算を行う1次調整系と、 複数の時系列データを用い、個体としての前記信号処理
パラメータ及び評価構造・パラメータを逐次調整する2
次調整系と、 2次調整系により逐次調整された評価構造・パラメータ
の中から評価精度の高い評価構造・パラメータを生成す
る3次調整系と、 を有する評価構造・パラメータ自動調整装置において、 前記1次調整系は、 プラントの時系列データに対し、設定した信号処理パラ
メータにより移動平均値及び変化量を演算する信号処理
手段と、 目的変数に対するしきい値を用いて状態総合評価対象デ
ータを分離する手段と、 説明変数ごとに、異常時データと全データの平均の差
と、全データの分散とを計算し、それらの比率を異常時
データが全データから乖離している度合いを示す説明変
数有意指標として計算する手段と、 この説明変数有意指標の大きい変数から順に説明変数を
選択する手段と、 を備えたことを特徴とするプラントの状態総合評価にお
ける評価構造・パラメータ自動調整装置。
An objective variable and a comprehensive evaluation method are given to a time series data of a plant, and a time series signal processing is performed by a signal processing parameter set separately to give a comprehensive evaluation of a total state of the plant. Identify fuzzy propositions, membership functions, fuzzy measures as evaluation structures and parameters, calculate the fitness of each fuzzy proposition,
A primary adjustment system for performing a fuzzy integral operation, a state comprehensive evaluation based on a fuzzy integral value, and a state comprehensive evaluation accuracy calculation; and sequentially adjusting the signal processing parameters and the evaluation structure / parameter as an individual using a plurality of time-series data.
An evaluation structure / parameter automatic adjustment device comprising: a secondary adjustment system; and a tertiary adjustment system that generates an evaluation structure / parameter with high evaluation accuracy from evaluation structures / parameters sequentially adjusted by the secondary adjustment system. The primary adjustment system separates the time series data of the plant from the signal processing means that calculates the moving average value and the amount of change using the set signal processing parameters, and the state evaluation data using the threshold value for the objective variable. Means to calculate the difference between the mean of abnormal data and the average of all data, and the variance of all data for each explanatory variable, and calculate the ratio between them as an explanatory variable indicating the degree to which abnormal data deviates from all data. Means for calculating as a significant index, and means for selecting an explanatory variable in order from a variable having a large explanatory variable significant index, a plant state characterized by comprising: Evaluation in case evaluation structure parameter automatic adjusting apparatus.
【請求項2】 請求項1に記載したプラントの状態総合
評価における評価構造・パラメータ自動調整装置におい
て、 前記1次調整系は、 選択された説明変数に対し、説明変数有意指標の符号か
らメンバシップ関数を決定する手段と、 全データの平均値と最大値または最小値とを用いて、異
常時データをカバーし、異常時データの全データからの
乖離が増すにつれてメンバシップ関数の値が増加するよ
うにメンバシップ関数パラメータを決定する手段と、 選択された説明変数と前記メンバシップ関数とからファ
ジィ命題を決定する手段と、 を備えたことを特徴とするプラントの状態総合評価にお
ける評価構造・パラメータ自動調整装置。
2. The evaluation structure / parameter automatic adjustment apparatus in the comprehensive state evaluation of a plant according to claim 1, wherein the primary adjustment system determines a membership of the selected explanatory variable from a sign of an explanatory variable significance index. Using the means to determine the function and the average value and the maximum or minimum value of all data, cover the abnormal data, and the value of the membership function increases as the deviation of the abnormal data from all data increases Means for determining a membership function parameter, and means for determining a fuzzy proposition from a selected explanatory variable and the membership function as described above. Automatic adjustment device.
【請求項3】 請求項2に記載したプラントの状態総合
評価における評価構造・パラメータ自動調整装置におい
て、 前記1次調整系は、 決定されたメンバシップ関数及びファジィ命題を用い、
プラントの計測値データを評価対象状態への適合度に変
換する手段と、 ファジィ測度積分により複数のファジィ命題の適合度を
一つのファジィ積分値に変換する手段と、 ファジィ命題の適合度と目的とするプラントの状態総合
評価値信号とによりファジィ測度を同定する手段と、 前記ファジィ積分値とこれに対するしきい値とからプラ
ントの状態総合評価値信号を生成する手段と、 プラントの計測値データ内の実績値から状態総合評価の
正誤判定を行う手段と、 正誤判定結果の中の正解数をカウントし、状態総合評価
精度を計算する手段と、 を備えたことを特徴とするプラントの状態総合評価にお
ける評価構造・パラメータ自動調整装置。
3. The evaluation structure / parameter automatic adjustment apparatus in the comprehensive state evaluation of a plant according to claim 2, wherein the primary adjustment system uses a determined membership function and a fuzzy proposition.
Means for converting the measured value data of the plant into the degree of conformity to the state to be evaluated, means for converting the degree of conformity of a plurality of fuzzy propositions into one fuzzy integral value by fuzzy measure integration, degree of conformity and purpose of the fuzzy proposition Means for identifying a fuzzy measure based on a total state evaluation value signal of the plant to be performed; means for generating a total state evaluation value signal of the plant from the fuzzy integration value and a threshold value for the fuzzy integration value; Means for performing a correct / incorrect judgment of the state comprehensive evaluation from the actual value; and means for counting the number of correct answers in the result of the correct / incorrect judgment and calculating the state comprehensive evaluation accuracy. Evaluation structure / parameter automatic adjustment device.
【請求項4】 請求項3に記載したプラントの状態総合
評価における評価構造・パラメータ自動調整装置におい
て、 前記2次調整系は、 信号処理パラメータの組を乱数により発生する手段と、 発生した信号処理パラメータによる信号処理を前記信号
処理手段により行った結果、発生するそれぞれ異なった
時系列データに対する前記評価構造・パラメータの決定
及び状態総合評価精度の計算を行う手段と、 設定したしきい値より評価精度が低い前記個体を削除す
る手段と、 残った個体に対し、設定したしきい値より評価精度の低
い個体に対してそのファジィ測度からファジィ命題間の
相互作用係数を求める手段と、 相互作用係数が最小であるファジィ命題を、前記説明変
数有意指標の次候補によるファジィ命題に置換する手段
と、 削除された個体に対し、信号処理パラメータを乱数によ
り新規に生成する手段と、 置換されたファジィ命題と、信号処理パラメータの新規
生成により突然変異したファジィ命題とを用いてファジ
ィ測度を再同定し、プラントの状態総合評価精度を再計
算する手段と、 請求項1〜請求項3において用いた時系列データを調整
データとして含み、それよりも設定数倍の時間帯の時系
列データを評価用データとして設定する手段と、 前記評価用データを用いて個体の状態評価精度を計算
し、評価精度の向上率が極小になったとき、または設定
回数だけ再評価を繰り返したときにこの反復計算を終了
し、その時に評価精度が最大である個体を準最適な個体
として保存する手段と、 を備えたことを特徴とするプラントの状態総合評価にお
ける評価構造・パラメータ自動調整装置。
4. The evaluation structure / parameter automatic adjustment device in the comprehensive state evaluation of a plant according to claim 3, wherein the secondary adjustment system generates a set of signal processing parameters by random numbers, and the generated signal processing. Means for determining the evaluation structure / parameters for each different time-series data generated as a result of performing the signal processing by the parameters by the signal processing means and calculating the state total evaluation accuracy; Means for deleting the individual having a low value, for the remaining individuals, means for obtaining an interaction coefficient between fuzzy propositions from an individual whose evaluation accuracy is lower than a set threshold from its fuzzy measure, Means for replacing the smallest fuzzy proposition with a fuzzy proposition by the next candidate of the explanatory variable significance index; For the individual, re-identify the fuzzy measure using the means for newly generating signal processing parameters by random numbers, the replaced fuzzy proposition, and the fuzzy proposition mutated by new generation of signal processing parameters, and the state of the plant Means for recalculating the overall evaluation accuracy; means for including the time-series data used in claim 1 as adjustment data, and setting time-series data in a time zone several times larger than the adjustment data as evaluation data Calculate the state evaluation accuracy of the individual using the evaluation data, and when the rate of improvement of the evaluation accuracy is minimal, or when re-evaluation is repeated a set number of times, this repetition calculation ends, Means for storing an individual having the highest evaluation accuracy as a sub-optimal individual, and an evaluation structure and parameters in the comprehensive state evaluation of the plant, comprising: Data automatic adjustment device.
【請求項5】 請求項4に記載したプラントの状態総合
評価における評価構造・パラメータ自動調整装置におい
て、 前記2次調整系は、 前回調整時に保存した準最適な個体を削除しないように
動作を修正する手段と、 今回調整により、評価精度が前回よりも設定率以上、向
上した場合には削除するように動作を修正する手段と、 を備えたことを特徴とするプラントの状態総合評価にお
ける評価構造・パラメータ自動調整装置。
5. The evaluation structure / parameter automatic adjustment apparatus in the comprehensive state evaluation of a plant according to claim 4, wherein the secondary adjustment system corrects an operation so as not to delete a sub-optimal individual stored at the time of the previous adjustment. And a means for correcting the operation so that when the evaluation accuracy is improved by more than a set rate compared to the previous time, the operation is corrected so as to delete the evaluation accuracy. -Automatic parameter adjustment device.
【請求項6】 請求項5に記載したプラントの状態総合
評価における評価構造・パラメータ自動調整装置におい
て、 前記3次調整系は、 複数の時系列データにより評価構造・パラメータを連続
的に調整したときの各評価構造・パラメータを保存する
手段と、 各評価構造・パラメータごとに各時系列データに対する
状態総合評価精度を計算する手段と、 各時系列データに対する状態総合評価精度の平均値を計
算する手段と、 その平均値が最大である評価構造・パラメータを保存す
る手段と、 各時系列データごとに各評価構造・パラメータに対する
状態総合評価精度を計算する手段と、 各評価構造・パラメータに対する状態総合評価精度の平
均値を計算する手段と、 その平均値が最小である時系列データを調整したときに
得られた評価構造・パラメータを選択する手段と、 選択した評価構造・パラメータから、ファジィ測度の相
互作用係数が最大であるファジィ命題を追加ファジィ命
題として保存する手段と、 状態総合評価精度の平均値が最大である評価構造・パラ
メータ中のファジィ命題としてのマスタファジィ命題に
前記追加ファジィ命題を追加して新たな合成ファジィ命
題を生成する手段と、 全ての時系列データに対して、前記マスタファジィ命題
及び合成ファジィ命題の状態総合評価精度を計算する手
段と、 合成ファジィ命題の状態総合評価精度がマスタファジィ
命題の状態総合評価精度よりも大きい場合に、合成ファ
ジィ命題をマスタファジィ命題とする手段と、 を備えたことを特徴とするプラントの状態総合評価にお
ける評価構造・パラメータ自動調整装置。
6. The evaluation structure / parameter automatic adjustment apparatus according to claim 5, wherein the tertiary adjustment system adjusts the evaluation structure / parameter continuously using a plurality of time series data. Means for storing each evaluation structure / parameter of the above, means for calculating the state total evaluation accuracy for each time-series data for each evaluation structure / parameter, and means for calculating the average value of the state total evaluation accuracy for each time-series data Means for storing the evaluation structure / parameter having the maximum average value; means for calculating the state total evaluation accuracy for each evaluation structure / parameter for each time-series data; and state total evaluation for each evaluation structure / parameter. Means for calculating the average value of accuracy, and the evaluation structure obtained when adjusting the time series data with the minimum average value Means for selecting parameters, means for storing the fuzzy proposition with the largest interaction coefficient of the fuzzy measure as an additional fuzzy proposition from the selected evaluation structure and parameters, and evaluation structure for which the average value of the state comprehensive evaluation accuracy is the largest Means for adding the additional fuzzy proposition to the master fuzzy proposition as a fuzzy proposition in the parameters to generate a new composite fuzzy proposition, and for all time-series data, the state of the master fuzzy proposition and the composite fuzzy proposition Means for calculating the overall evaluation accuracy; and means for using the composite fuzzy proposition as a master fuzzy proposition when the state overall evaluation accuracy of the composite fuzzy proposition is greater than the state overall evaluation accuracy of the master fuzzy proposition. Evaluation system / parameter automatic adjustment device for comprehensive evaluation of plant condition.
JP9266198A 1997-09-30 1997-09-30 Automatic evaluation structure/parameter adjusting device for total state evaluation of plant Withdrawn JPH11110367A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP9266198A JPH11110367A (en) 1997-09-30 1997-09-30 Automatic evaluation structure/parameter adjusting device for total state evaluation of plant

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP9266198A JPH11110367A (en) 1997-09-30 1997-09-30 Automatic evaluation structure/parameter adjusting device for total state evaluation of plant

Publications (1)

Publication Number Publication Date
JPH11110367A true JPH11110367A (en) 1999-04-23

Family

ID=17427633

Family Applications (1)

Application Number Title Priority Date Filing Date
JP9266198A Withdrawn JPH11110367A (en) 1997-09-30 1997-09-30 Automatic evaluation structure/parameter adjusting device for total state evaluation of plant

Country Status (1)

Country Link
JP (1) JPH11110367A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012522307A (en) * 2009-03-31 2012-09-20 東京エレクトロン株式会社 Method and system for tool performance degradation and mismatch detection
CN111915171A (en) * 2020-07-21 2020-11-10 中国石油大学(华东) Comprehensive evaluation method for deep plugging and adjusting effect
EP3792750A1 (en) 2019-09-13 2021-03-17 Fujitsu Limited Information processing apparatus, information processing program, and information processing method

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012522307A (en) * 2009-03-31 2012-09-20 東京エレクトロン株式会社 Method and system for tool performance degradation and mismatch detection
EP3792750A1 (en) 2019-09-13 2021-03-17 Fujitsu Limited Information processing apparatus, information processing program, and information processing method
US11449360B2 (en) 2019-09-13 2022-09-20 Fujitsu Limited Information processing apparatus, non-transitory computer-readable storage medium for storing information processing program, and information processing method
CN111915171A (en) * 2020-07-21 2020-11-10 中国石油大学(华东) Comprehensive evaluation method for deep plugging and adjusting effect
CN111915171B (en) * 2020-07-21 2024-02-13 中国石油大学(华东) Comprehensive evaluation method for depth plugging effect

Similar Documents

Publication Publication Date Title
US20030033263A1 (en) Automated learning system
CN110705807A (en) Semiconductor production quality prediction method and system based on field data
CN110866832A (en) Risk control method, system, storage medium and computing device
CN112288137A (en) LSTM short-term load prediction method and device considering electricity price and Attention mechanism
CN111898637B (en) Feature selection algorithm based on Relieff-DDC
CN115221942A (en) Equipment defect prediction method and system based on time sequence fusion and neural network
JPH11110367A (en) Automatic evaluation structure/parameter adjusting device for total state evaluation of plant
US20220156573A1 (en) Machine Learning Engine Providing Trained Request Approval Decisions
CN112380273A (en) Data estimation method and device
US20090326830A1 (en) Interpreting dna
EP1849115B1 (en) Mass spectrometry precursor ion selection
JP7202757B1 (en) Information processing system, information processing method and program
Jelinek et al. Diagnostic with incomplete nominal/discrete data
AU2002353190B2 (en) Improvements in and relating to interpreting DNA
CN114820074A (en) Target user group prediction model construction method based on machine learning
CN109308565B (en) Crowd performance grade identification method and device, storage medium and computer equipment
Marwala et al. Neuro-rough models for modelling HIV
CN118193855B (en) Online response method and device and computer equipment
JP2004206167A (en) Case prediction device and method
CN114090795B (en) Two-way interaction intelligent service system and method based on equipment knowledge graph
Dube et al. Runtime Prediction of Machine Learning Algorithms in Automl Systems
CN117893216B (en) User behavior analysis method and system based on big data
Shi et al. How to combine CBR and RBR for diagnosing multiple medical disorder cases
KR20200021132A (en) Method for Processing Prediction by using Artificial Intelligence
CN117273467B (en) Multi-factor coupling-based industrial safety risk management and control method and system

Legal Events

Date Code Title Description
A300 Application deemed to be withdrawn because no request for examination was validly filed

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20041207