JP3219116B2 - Error diagnosis method - Google Patents
Error diagnosis methodInfo
- Publication number
- JP3219116B2 JP3219116B2 JP33252293A JP33252293A JP3219116B2 JP 3219116 B2 JP3219116 B2 JP 3219116B2 JP 33252293 A JP33252293 A JP 33252293A JP 33252293 A JP33252293 A JP 33252293A JP 3219116 B2 JP3219116 B2 JP 3219116B2
- Authority
- JP
- Japan
- Prior art keywords
- section
- data
- function
- point
- sum
- 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
Links
Landscapes
- Testing And Monitoring For Control Systems (AREA)
- Complex Calculations (AREA)
Description
【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION
【0001】[0001]
【産業上の利用分野】本発明は、石油化学等の工業プロ
セスに設置されるプロセス制御システムに内蔵して使用
され、工業プロセスの異常診断あるいは故障診断等に用
いられる異常診断方法に関し、特にプロセスデータから
プロセスの正常と異常とを識別する境界値を容易に求め
る方法を実現するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an abnormality diagnosis method which is used in a process control system installed in an industrial process such as petrochemical or the like and is used for abnormality diagnosis or failure diagnosis of an industrial process. This realizes a method for easily obtaining a boundary value for discriminating between normal and abnormal processes from data.
【0002】[0002]
【従来の技術】種々の工業プロセスには図10に示すよ
うな分散型プロセス制御システムを設置することが多
い。対象となるプロセスPから温度、圧力、流量等のプ
ロセスデータをフイールドコントロールステーション
(FCS)111に入力して、これらのデータに制御演
算を施して、弁開度制御信号やリレー接点信号等を出力
する。2. Description of the Related Art Various industrial processes are often provided with a distributed process control system as shown in FIG. Process data such as temperature, pressure, and flow rate from the target process P is input to a field control station (FCS) 111, and these data are subjected to a control operation to output a valve opening control signal, a relay contact signal, and the like. I do.
【0003】オペレータステーション(OPS)112
では、オペレータが前記の入出力信号の監視やプロセス
の運転操作等を行う。同一出願人による特願平4−85
328異常診断方法は、上記の温度、圧力、流量等のプ
ロセスデータに許される上限及び下限値を経験により予
め求め、また上限値または下限値の隣合う2点間の値は
直線で近似して求め、上記のプロセスデータと照合して
異常診断を行うものである。The operator station (OPS) 112
Then, the operator monitors the input / output signals and operates the process. Japanese Patent Application No. 4-85 filed by the same applicant
In the 328 abnormality diagnosis method, the upper limit and the lower limit allowed for the process data such as the temperature, the pressure, and the flow rate are previously obtained by experience, and the value between two adjacent points of the upper limit or the lower limit is approximated by a straight line. Then, the abnormality diagnosis is performed by collating with the above process data.
【0004】更に汎用ワークステーション113を設置
して、AI手法により異常診断や故障診断を行う方法が
あり、システム稼働中に得られたデータに予め作成して
おいたプロダクションルール(IF THEN RUL
Eとも呼ばれる)を主体とした異常診断アルゴリズムを
適用して、プロセスやシステム内で発生した異常や故障
を発見するものである。Further, there is a method in which a general-purpose workstation 113 is installed and an abnormality diagnosis or a failure diagnosis is performed by the AI method. A production rule (IF THEN RUL) prepared in advance on data obtained while the system is operating.
E), and discovers abnormalities or failures that have occurred in the process or system by applying an abnormality diagnostic algorithm mainly based on E).
【0005】[0005]
【発明が解決しようとする課題】上記のようなAI手法
を用いた異常診断方式では幅広い異常診断が可能になっ
たが、異常診断アルゴリズムを構築する際、プロセスの
正常と異常との識別境界値を決定することが難しく、ま
たオペレータがAI手法を理解することも含めて多くの
時間を要するという問題があった。The abnormality diagnosis method using the AI method as described above enables a wide range of abnormality diagnosis. However, when constructing an abnormality diagnosis algorithm, a threshold value for discriminating between normal and abnormal processes is considered. Is difficult to determine, and it takes much time for the operator to understand the AI method.
【0006】本発明の目的は、上記の様な問題点を解決
するために、オペレータにとって理解しやすく、短い時
間に簡単にプロセスの正常と異常との識別境界値を決定
でき、しかも最適な基準に基づいてプロセスの正常と異
常を判断できる異常診断方法を実現することである。SUMMARY OF THE INVENTION It is an object of the present invention to solve the above-mentioned problems and to make it easy for an operator to understand, to easily determine a discrimination boundary value between a normal process and an abnormal process in a short time , and to determine an optimal reference value. Process based on normal and abnormal
It is an object of the present invention to realize an abnormality diagnosis method that can always determine .
【0007】[0007]
【課題を解決するための手段】上記のような目的を達成
するために本発明は、対象とするプロセスから順次計測
データを収集して、前記プロセスが異常であるか否かを
診断する異常診断方法において、 (a)収集保存しているプロセスの計測データを順次プ
ロットして、対象とするデータの範囲をL個の区間に等
分割する手順と、 (b)第1区間の始点を起点としてL個の区間につい
て、それぞれデータ各点からの距離との差の二乗和が最
少になる折れ線函数を求める手順と、 (c)第2区間の始点を起点としてL個の区間につい
て、それぞれデータ各点からの距離との差の二乗和が最
少になる折れ線函数を求め、以下同様にして第L区間の
始点を起点とした場合まで折れ線関数を求める手順と、 (d)第1区間の始点を起点とした場合にL個の区間に
ついてそれぞれ求めたデータ各点からの距離の差の二乗
和の合計をSUM1とし、同様にしてSUM2からSU
MLまで求め、SUM1からSUMLの中で最小値にな
るものについての折れ線函数を中央折れ線函数とする手
順と、 (e)上記のようにして求めた中央折れ線函数の各区間
毎に中央折れ線函数からの標準偏差σを求める手順と、 (f)上記標準偏差σのうち最大のものを共通標準偏差
値として、経験値係数αを乗じて、上限、下限の各折れ
線函数を作成する手順と、 (g)この上限、下限の各折れ線函数で挟まれる範囲に
存在しない計測データに基づいてプロセスが異常である
と診断して信号を出力する手順とを有することを特徴と
する異常診断方法である。In order to achieve the above object, the present invention sequentially collects measurement data from a target process and determines whether the process is abnormal. In the abnormality diagnosis method for diagnosing, (a) a procedure of sequentially plotting measurement data of a collected and stored process and equally dividing a range of target data into L sections; Starting from the starting point,
Then, a procedure for finding a polygonal function that minimizes the sum of squares of the difference from the distance from each point of the data, and (c) starting from the start point of the second section as the starting point
Then , a polygonal function that minimizes the sum of squares of the difference from the distance from each data point is obtained , and the same applies to the L-th section in the same manner.
A procedure for obtaining a polygonal line function until the start point is set as a starting point;
The square of the difference of the distance from each point
The sum of the sum is defined as SUM1, and similarly from SUM2 to SU
ML and the minimum value from SUM1 to SUML
The line function of the object as the central line function
Common standard deviation and procedures, the largest of (f) the standard deviation σ to obtain a forward, a standard deviation σ from the center fold line function for each section of the central fold line function obtained as (e) above (G) multiplying the empirical value coefficient α as the value to create each of the upper and lower broken line functions; and (g) the process is abnormal based on the measurement data that does not exist in the range between the upper and lower broken line functions. And outputting a signal after diagnosing the abnormality.
【0008】[0008]
【作用】この発明では、収集保存しているプロセスの計
測データを順次プロットして、対象とするデータの範囲
をL個の区間に等分割する。 L個の区間の中で、第1区
間の始点を起点としてL個の区間について、それぞれデ
ータ各点からの距離との差の二乗和が最少になる折れ線
函数を求める。 第2区間の始点を起点としてL個の区間
について、それぞれデータ各点からの距離との差の二乗
和が最少になる折れ線函数を求める。以下同様にして第
L区間の始点を起点とした場合まで折れ線関数を求め
る。 第1区間の始点を起点とした場合にL個の区間につ
いてそれぞれ求めたデータ各点からの距離の差の二乗和
の合計をSUM1とする。同様にしてSUM2からSU
MLまで求める。SUM1からSUMLの中で最小値に
なるものについての折れ線函数を中央折れ線函数とす
る。このようにして求めた中央折れ線函数は、対象とす
る全区間(L個の区間)のデータを最もよく近似した折
れ線函数である。 According to the present invention, the total of the collected and stored processes is calculated.
Measurement data is plotted sequentially, and the range of the target data
Is equally divided into L sections. First section of L sections
Starting from the start point between the L segments,
Line that minimizes the sum of squares of the difference from the distance from each point
Find a function. L sections starting from the start point of the second section
For each, the square of the difference from the distance from each data point
Find the linear function that minimizes the sum. And so on
Find the line function until the start point of the L section
You. If the starting point of the first section is the starting point,
Sum of squares of the difference in distance from each point
Is set to SUM1. Similarly, from SUM2 to SU
Find up to ML. From SUM1 to the minimum value in SUML
Assume that the line function of
You. The central line function obtained in this way is the target
Of the best approximation of the data of all sections (L sections)
Is a linear function.
【0009】前記の中央折れ線函数と当該中央折れ線函
数の作成に使用した計測データとの間の標準偏差を求
め、そのうちの最大のものを共通標準偏差値とし経験値
係数αを乗じて、各区間につきプロセスの異常限界であ
る上限、下限の各折れ線函数を作成しているので、対象
とする計測の範囲で診断は連続して行えるし、計測デー
タのバラツキによる誤診を少なくすることもできる。The above-described central polygonal function and the central polygonal function
Calculate the standard deviation between the measurement data used to create the numbers , and multiply the maximum value as the common standard deviation value by the empirical coefficient α, and calculate the upper and lower limits that are the process abnormal limits for each section. Since the line function is created, diagnosis can be performed continuously within the range of the target measurement, and erroneous diagnosis due to variation in measurement data can be reduced.
【0010】運転中の計測データがこの上限、下限を表
す各折れ線函数で挟まれる範囲に存在しない場合にプロ
セスが異常であると診断して信号を出力する。このプロ
セスの異常限界を表す上限、下限の各折れ線函数の作成
には予め一定の周期でサンプリングし、保存してある運
転データを使用するので変化に即応することができる。If the measured data during operation does not exist in the range between the broken line functions representing the upper and lower limits, the process is diagnosed as abnormal and a signal is output. Since the upper and lower broken line functions representing the abnormal limit of this process are created at regular intervals and sampling operation data is used, it is possible to respond immediately to changes.
【0011】[0011]
【実施例】以下図面を用いて本発明を説明する。図1は
本発明の方法を実現する機能ブロック図である。通常本
発明の機能はソフトウエア的に実現されるものである。
そのために汎用ワークステーション113はマンマシン
インタフェースとして使用する。BRIEF DESCRIPTION OF THE DRAWINGS FIG. FIG. 1 is a functional block diagram for realizing the method of the present invention. Usually, the functions of the present invention are realized by software.
Therefore, the general-purpose workstation 113 is used as a man-machine interface.
【0012】プロセスPからフイールドコントロールス
テーション(FCS)111を経由して温度、圧力その
他のプロセスデータを一定の周期で収集保存するデータ
収集部11、データ収集部11の動作を規定するデータ
収集ジェネレータ12、データ収集部11から診断の為
の知識ベースを作成する目的のデータを記憶させるプロ
セスデータ保存部13、対象とするデータの異常診断条
件を決定する異常監視範囲決定部14、ワークステーシ
ョンの専用画面を使用して、異常診断に必要とされるデ
ータ、タグナンバー、異常診断条件を定義して登録する
異常診断定義ファイル15、異常診断知識ベースを作成
する異常診断知識ベースジェネレータ16、異常診断の
ための条件が記述されている簡易異常診断知識ベース1
7、データ収集部11から送られたプロセスデータと簡
易異常診断知識ベース17の内容を比較して異常か否か
を診断する異常診断部18、異常診断の結果を表示する
異常診断表示部19とから構成される。A data collection unit 11 that collects and stores temperature, pressure, and other process data from the process P via a field control station (FCS) 111 at regular intervals, and a data collection generator 12 that defines the operation of the data collection unit 11 A process data storage unit 13 for storing target data for creating a knowledge base for diagnosis from the data collection unit 11, an abnormality monitoring range determining unit 14 for determining an abnormality diagnosis condition of target data, and a dedicated screen of a workstation. , An abnormality diagnosis definition file 15 for defining and registering data, tag numbers, and abnormality diagnosis conditions required for abnormality diagnosis, an abnormality diagnosis knowledge base generator 16 for creating an abnormality diagnosis knowledge base, and an abnormality diagnosis Simple abnormality diagnosis knowledge base 1 describing the conditions of
7. An abnormality diagnosis unit 18 that compares the process data sent from the data collection unit 11 with the contents of the simple abnormality diagnosis knowledge base 17 to determine whether there is an abnormality, and an abnormality diagnosis display unit 19 that displays the result of the abnormality diagnosis. Consists of
【0013】異常診断を開始すると、プロセスPからフ
イールドコントロールステーション(FCS)111を
経由して温度、圧力その他のプロセスデータを一定の周
期で収集しプロセスデータ保存部13へ保存したあと異
常診断部18へ処理を移す。異常診断部18はプロセス
データと簡易異常診断知識ベース17の内容を参照して
プロセスの異常を診断し、診断結果を異常診断表示部1
9によりワークステーション113の画面に表示する。When the abnormality diagnosis is started, temperature, pressure, and other process data from the process P via the field control station (FCS) 111 are collected at a fixed cycle, and stored in the process data storage unit 13. Move the process to The abnormality diagnosis unit 18 diagnoses a process abnormality by referring to the process data and the contents of the simple abnormality diagnosis knowledge base 17 and displays the diagnosis result in the abnormality diagnosis display unit 1.
9 is displayed on the screen of the workstation 113.
【0014】異常監視範囲決定部14が、対象とするデ
ータの異常診断条件を決定する動作は、プロセスデータ
保存部13に保存されているデータを基にして、後に説
明する中央折れ線函数31を求め、次に監視範囲を決定
することである。この結果により異常診断定義ファイル
15へ監視範囲を設定する。図2は本発明の処理動作の
STEPを示す。The operation of the abnormality monitoring range determination unit 14 to determine the abnormality diagnosis condition of the target data is performed by obtaining a central polygonal line function 31 described later based on the data stored in the process data storage unit 13. Then, the monitoring range is determined. The monitoring range is set in the abnormality diagnosis definition file 15 based on the result. FIG. 2 shows a STEP of the processing operation of the present invention.
【0015】STEP1〜4の、収集したデータから折
れ線函数を求める処理は特許請求の範囲(a)〜(c)
に、STEP5の、中央折れ線函数31を求める処理は
(d)に、STEP6の、標準偏差を求める処理は
(e)に、STEP7の、監視範囲を決定する処理は
(f)に、STEP8〜9の異常診断をして結果を出力
する処理は(g)にそれぞれ対応する。Claims (a) to (c) in which the processing for obtaining the polygonal line function from the collected data in STEPs 1 to 4 are described below.
The processing for obtaining the central polygonal line function 31 in STEP 5 is shown in (d), the processing for obtaining the standard deviation in STEP 6 is shown in (e), the processing for determining the monitoring range in STEP 7 is shown in (f), and the processing in STEPs 8-9. The process of diagnosing the abnormality and outputting the result corresponds to (g).
【0016】次に本発明にかかる異常診断方法の原理を
説明する。先ず最初に、プロセスの運転条件、例えば燃
料流量を増加させて各部の温度を上昇なせるなどの操
作、を変えてプロセスから所定の計測データを得る。図
3はプロセスから収集し保存してあった計測データを順
次プロットしたもので、対象とするデータの範囲X0〜
XLをL個の区間に等分割することを示す。この基本と
なる計測データについて下記に説明する中央折れ線函数
を求める。図3のX0〜XLはプロセスの運転条件の変動
範囲に相当する。所定の条件Xで運転中に計測値Yが上
記の中央折れ線函数から求めた監視範囲を離脱した場合
を異常と診断するものである。Next, the principle of the abnormality diagnosis method according to the present invention will be described. First, predetermined measurement data is obtained from the process by changing the operating conditions of the process, for example, the operation of increasing the temperature of each part by increasing the fuel flow rate. FIG. 3 is a graph in which measurement data collected and stored from the process is sequentially plotted, and the range of target data X 0 to
Indicating that equally divide the X L into L intervals. A central polygonal function described below is obtained for the basic measurement data. X 0 to X L in FIG. 3 corresponds to the variation range of the process operating conditions. An abnormality is diagnosed when the measured value Y deviates from the monitoring range obtained from the above-described central polygonal line function during operation under the predetermined condition X.
【0017】以下に最小二乗法の直線回帰分析の手法を
用いて、中央折れ線函数を求めるアルゴリズムを(1)
項に、監視範囲を求めるアルゴリズムを(2)項にST
EPを追って詳しく説明する。 (1)中央折れ線函数を表すアルゴリズムを求める。図
3はプロセスから収集し保存してあった計測データを順
次プロットしたもので、対象とするデータの範囲X0〜
XLをL個の区間に等分割することを示す。(STEP
2〜3) 運転条件Xのときの計測値がYである。XとYの関係は
範囲X0〜XLの小さい範囲では一次函数で表されるもの
と考えられる。全範囲はL等分してそれぞれの一次函数
を求める。An algorithm for obtaining a central polygonal line function using the method of linear regression analysis of the least squares method will be described below (1).
In the section, the algorithm for obtaining the monitoring range is described in the section (2), ST.
This will be described in detail following the EP. (1) Find an algorithm that represents the central polygonal line function. FIG. 3 is a graph in which measurement data collected and stored from the process is sequentially plotted, and the range of target data X 0 to
Indicating that equally divide the X L into L intervals. (STEP
2-3) The measured value under the operating condition X is Y. Relationship between X and Y is considered to be represented by a linear function is a small range of the range X 0 to X L. The entire range is divided into L equal parts to obtain the respective primary functions.
【0018】図4は対象とするデータの範囲X0〜XLの
中の第K区間を拡大して示し、その中にk1〜knのn個
のデータがあることを示す。この区間内のデータを折れ
線函数Y=aX+bで代表させると、データ各点とこの
折れ線函数31との距離の残差の二乗和(以下Sで表
す)が最少になる場合のa及びbで表したY=aX+b
(式1)がこれらのデータを最も正しく表すものとな
る。このa及びbは下記の式で求まるものである。 Y=aX+b 式1 各区間内にn個のデータがあるものとし、Σは区間毎に
i=1からi=nまでのn回の和を表す記号とする。[0018] Figure 4 is an enlarged view of the first K domain in a range X 0 to X L of the data of interest, indicating that there are n data k 1 to k n therein. If the data in this section is represented by a polygonal line function Y = aX + b, it is expressed by a and b when the sum of squares of the residual of the distance between each data point and the polygonal line function 31 is minimized (hereinafter represented by S). Y = aX + b
(Equation 1) represents these data most correctly. A and b are determined by the following equations. Y = aX + b Equation 1 Assume that there are n pieces of data in each section, and 記号 is a symbol representing the sum of n times from i = 1 to i = n for each section.
【0019】また接続点のX及びYの関係は図3のとお
りである。第K区間のa、b、Sをak、bk、Skで表
し計算には次の式を使用する。 ak=(ΣXiΣYi−nΣXiYi)/((ΣXi)2−nΣXi 2) 式2 bk=1/n×(ΣYi−akΣXi) 式3 Sk=Σ(Yi−(akXi+bk))2 式4 図5は、第K区間以降及び以前の区間では先に決定した
折れ線とつながることを制約条件として、区間内のデー
タ各点からの距離との残差の二乗和が最少になる合計L
個の折れ線函数を求める動作を示す。FIG. 3 shows the relationship between X and Y at the connection points. The a, b, and S of the K-th section are represented by a k , b k , and S k , and the following equation is used for calculation. a k = (ΣX i ΣY i -nΣX i Y i ) / ((ΣX i ) 2 -nΣX i 2 ) Formula 2 b k = 1 / n × (ΣY i -a k ΣX i) Formula 3 S k = Σ (Y i - (a k X i + b k)) 2 Equation 4 FIG. 5 shows that the total L that minimizes the sum of squares of the residual with the distance from each data point in the section under the constraint that connection to the polygonal line determined earlier in the section after the K-th section and before is used as a constraint.
The operation of obtaining the polygonal line functions will be described.
【0020】この制約条件があるため、第K区間の前後
の区間で適用する計算式は異なる。第K+1の区間から
a、b、Sの計算を始めるときは次の式を使用する。 ak+1=Σ((Xk−Xi)(Yk−Yi))/Σ(Xk−Xi)2 式5 bk+1=Yk−ak+1Xk 式6 Sk+1=Σ(Yi−(ak+1Xi+bk+1))2 式7 第K−1の区間からa、b、Sの計算を始めるときは次
の式を使用する。 ak-1=Σ((Xk-1−Xi)(YK-1−Yi))/Σ(Xk-1−Xi)2 式8 bk-1=Yk-1−ak-1Xk-1 式9 Sk-1=Σ(Yi−(ak-1Xi+bk-1))2 式10 図6は、第1区間、第2区間・・・第L区間の始点を起
点として各々L個の区間内でデータ各点からの距離との
残差の二乗和Sを求める折れ線函数のマトリクスを示
す。上記の式2〜式10に基づきこの各欄を計算する。
(STEP4) ここでSUMl’、al’、bl’、Sl’をl’=1
−Lとしたときに、SUM1は第1区間を始点として第
L区間迄の各区間を直線近似する折れ線函数を表す
a1、b1〜aL、bL及びS1〜SLを示す欄である。例え
ばLを6とし、第4区間を起点とすると、起点番号l’
は4となり、順に6までのa4〜a6、b4〜b6、S4〜
S6を求める。次にl’−1から逆順に1までのa3〜a
1、b3〜b1、S3〜S1を求める。このようにして、a1
〜a6、b1〜b6、S1〜S6を求める。このようにし
て、SUMLは第L区間の始点を起点として第L−1区
間の最後のデータ迄の各区間の折れ線を表すa、b、及
びSを示す欄である。Due to these constraints, the calculation formulas applied in the sections before and after the K-th section are different. The following equation is used to start the calculation of a, b, and S from the (K + 1) th section. a k + 1 = Σ ((X k −X i ) (Y k −Y i )) / Σ (X k −X i ) 2 Equation 5 b k + 1 = Y k −a k + 1 X k Equation 6 S k + 1 = Σ (Y i − (a k + 1 X i + b k + 1 )) 2 Equation 7 When starting the calculation of a, b, and S from the K−1 section, the following equation is used. . a k−1 = Σ ((X k−1 −X i ) (Y K−1 −Y i )) / Σ (X k−1 −X i ) 2 Equation 8 b k−1 = Y k−1 − a k−1 X k−1 Equation 9 S k−1 = Σ (Y i − (a k−1 X i + b k−1 )) 2 Equation 10 FIG. 6 shows a first section, a second section,. 13 shows a matrix of a polygonal function for calculating a sum of squares S of a difference from a distance from each data point in each of L sections starting from a start point of the L-th section. Each of these fields is calculated based on the above equations 2 to 10.
(STEP 4) Here, SUMl ', al', bl ', and Sl' are represented by l '= 1.
When a -L, SUM 1 shows the a 1, b 1 ~a L, b L and S 1 to S L representing a polygonal line function that linearly approximates the each section until the L section of the first section as the starting point Column. For example, if L is 6, and the fourth section is the starting point, the starting point number l ′
Is 4, and a 4 to a 6 , b 4 to b 6 , S 4 to
Determine the S 6. Next, a 3 -a from l′−1 to 1 in reverse order
1 , b 3 to b 1 and S 3 to S 1 are obtained. Thus, a 1
Aa 6 , b 1 b 6 and S 1 〜S 6 are obtained. In this manner, SUM L is a column indicating a, b, and S representing a polygonal line in each section from the start point of the L-th section to the last data of the L-1 section.
【0021】ここで第1区間を起点にしてL全区間につ
いて求めたS1〜SLの合計をSUM 1とする。同様にし
て第L区間までの各区間の始点を起点にもつそれぞれの
S1〜SLの合計を求めてSUM2〜SUMLとする。この
SUM1からSUMLの中で最少となる欄にある第1区間
のa1、b1〜第L区間のaL、bLをそれぞれの区間の折
れ線函数Y=aX+bのa及びbとして採用すれば対象
とする全区間のデータを最も良く近似すると考える。こ
れを中央折れ線函数と言う。(STEP5) それぞれの区間のこの折れ線函数をY1、Y2、・・YL
で表す。 (2)監視範囲を表すアルゴリズムを計算する。Here, starting from the first section, all L sections
S1~ SLSum of 1And Likewise
Of each section up to the L-th section
S1~ SLSUM for the sum ofTwo~ SUMLAnd this
SUM1To SUML1st section in the smallest column of
A1, B1~ A of the L-th sectionL, BLFor each section
If it is adopted as a and b of the linear function Y = aX + b
It is considered that the data of all the sections is best approximated. This
This is called the central polygonal function. (STEP5) This line function of each section is Y1, YTwo, ... YL
Expressed by (2) Calculate an algorithm representing the monitoring range.
【0022】図7に第1区間から第L区間までの各区間
毎にプロセスデータと上記中央折れ線函数からの標準偏
差を求める動作を第1区間から第L区間まで繰り返すこ
とを示す。(STEP6) 計算は次の式を使用する。 σk=[Σ(Yi−(akXi+bk))2−[Σ(Yi−(akXi+bk))]2 /n]/(n−1) の平方根 式11 図8は各区間毎にプロセスデータと中央折れ線函数から
求めた標準偏差σの中から最大のものを共通標準偏差値
σmaxとして選び、経験値係数αを乗じて、第1区間か
ら第L区間の区間ごとに上限、下限の各折れ線函数を作
成する手順を示す。FIG. 7 shows that the operation of obtaining the process data and the standard deviation from the central polygonal line function is repeated for each section from the first section to the L-th section from the first section to the L-th section. (STEP 6) The following formula is used for the calculation. σ k = [Σ (Y i − ( ak X i + b k )) 2 − [Σ (Y i − ( ak X i + b k ))] 2 / n] / (n−1) FIG. 8 shows that, for each section, the largest one among the standard deviations σ obtained from the process data and the central polygonal line function is selected as the common standard deviation value σ max , multiplied by the empirical coefficient α, and the first to the L-th section. The procedure for creating each of the upper and lower polygonal line functions for each section of will be described.
【0023】次の通り計算をする。 第1区間の監視範囲の上限折れ線函数YHIGH1=Y1+σ
max×α 〜第L区間の監視範囲の上限折れ線函数YHIGHL=YL+
σmax×α 第1監視範囲の下限折れ線函数YLOW1=Y1−σmax×α 〜第L監視範囲の下限折れ線函数YLOWL=YL−σmax×
α X軸の両端にある監視限界は異常監視範囲決定部14に
より予め設定する。The calculation is performed as follows. Upper limit line function Y HIGH1 = Y 1 + σ of the monitoring range of the first section
max × α to the upper limit polygonal function Y HIGHL = Y L +
sigma max × alpha lower polygonal line function Y of the first monitoring range LOW1 = Y 1 -σ max × α ~ the L monitoring range of the lower limit polygonal line function Y LOWL = Y L -σ max ×
The monitoring limits at both ends of the α X-axis are set in advance by the abnormality monitoring range determination unit 14.
【0024】上記のようにして得られた上限、下限の各
折れ線函数は異常診断定義ファイルに設定される。(S
TEP7〜8)図9は上記のようにして求めた中央折れ
線函数は31、上限折れ線函数は92、下限折れ線函数
は93で示している。計測データが上限折れ線函数92
と下限折れ線函数93で挟まれる範囲に存在しない場合
にプロセスは異常であると診断して信号を出力する。The upper and lower broken line functions obtained as described above are set in the abnormality diagnosis definition file. (S
(TEP 7 to 8) FIG. 9 shows the central polygonal function obtained as described above as 31, the upper limit polynomial function as 92, and the lower polynomial function as 93. Measurement data is upper limit line function 92
If the process does not exist in the range between the lower limit linear function 93 and the lower limit line function 93, the process is diagnosed as abnormal and a signal is output.
【0025】図10は本発明の異常診断方法を内蔵した
分散型プロセス制御システムの例を示す。対象となるプ
ロセスPから温度、圧力、流量等のプロセスデータをフ
イールドコントロールステーション(FCS)111に
入力して、これらのデータに制御演算を施して、弁開度
制御信号やリレー接点信号等を出力する。オペレータス
テーション(OPS)112では、オペレータが前記の
入出力信号の監視やプロセスの運転操作等を行う。FIG. 10 shows an example of a distributed process control system incorporating the abnormality diagnosis method of the present invention. Process data such as temperature, pressure, and flow rate from the target process P is input to a field control station (FCS) 111, and these data are subjected to a control operation to output a valve opening control signal, a relay contact signal, and the like. I do. In the operator station (OPS) 112, an operator monitors the input / output signals and operates the process.
【0026】汎用ワークステーション113は異常診断
や故障診断等特定の作業を行うための専用操作卓であ
る。The general-purpose workstation 113 is a dedicated console for performing specific operations such as abnormality diagnosis and failure diagnosis.
【0027】[0027]
【発明の効果】本発明によれば、運転中に保存したデー
タを用いて、プロセスが正常であるか否かの監視範囲を
決定するため、経年変化等による監視範囲の変更や運転
条件の変更に伴なって異常限界を変更するために上限、
下限の各折れ線函数を変更する要求にも速やかな対応が
可能である。また、対象とする全区間のデータを最もよ
く近似した中央折れ線函数をもとに判断をしているた
め、最適な基準に基づいてプロセスの正常と異常を判断
できる。 According to the present invention, a monitoring range for determining whether or not a process is normal is determined by using data stored during operation. The upper limit to change the abnormal limit with the
It is possible to promptly respond to a request to change each lower-order linear function. In addition, the data of all target sections
Judgment is made based on the approximate central line function
To determine whether the process is normal or abnormal based on optimal criteria
it can.
【0028】また経験を基にして監視範囲や異常限界を
合理的にかつ容易に決定することができる。保存データ
や監視範囲をワークステーションの専用画面へ表示する
ことにも適しているので視覚的把握にも便利である。以
上のように、異常診断及び故障診断を実現する方法とし
ては比較的簡単に、かつ短時間に結果を得られるのでシ
ステムの構築が簡便であり、保守性及び信頼性の向上に
も顕著な効果がある。The monitoring range and the abnormal limit can be rationally and easily determined based on experience. It is also suitable for displaying stored data and monitoring range on a dedicated screen of the workstation, so that it is convenient for visual grasp. As described above, as a method for realizing the abnormality diagnosis and the failure diagnosis, the result can be obtained relatively easily and in a short time, so that the system can be easily constructed, and the remarkable effect is also provided for the improvement of the maintainability and reliability. There is.
【図1】本発明の方法を実現する機能ブロック図であ
る。FIG. 1 is a functional block diagram for realizing the method of the present invention.
【図2】本発明の処理動作手順を示す。FIG. 2 shows a processing operation procedure of the present invention.
【図3】保存している対象プロセスデータをプロットし
た例を示す。FIG. 3 shows an example in which stored target process data is plotted.
【図4】図3の区間Kを拡大した説明図である。FIG. 4 is an explanatory diagram in which a section K in FIG. 3 is enlarged.
【図5】残差の二乗の和を求める動作を示す流れ図であ
る。FIG. 5 is a flowchart showing an operation for obtaining a sum of squares of residuals.
【図6】折れ線函数のマトリクスを示す図である。FIG. 6 is a diagram showing a matrix of polygonal functions.
【図7】標準偏差演算処理を示す流れ図である。FIG. 7 is a flowchart showing a standard deviation calculation process.
【図8】上下限折れ線函数の決定処理の流れ図である。FIG. 8 is a flowchart of a process for determining upper and lower limit polygonal line functions.
【図9】中央折れ線と上下限折れ線の関係を示す図であ
る。FIG. 9 is a diagram showing a relationship between a center polygonal line and upper and lower polygonal lines.
【図10】システム構成図である。FIG. 10 is a system configuration diagram.
11 データ収集部 12 データ収集ジェネレータ 13 プロセスデータ保存部 14 異常監視範囲決定部 15 異常診断定義ファイル 16 異常診断知識ベースジェネレ−タ 17 簡易異常診断知識ベース 18 異常診断部 19 異常診断表示部 31 プロセスデータを代表する折れ線函数 42 K区間の上限折れ線函数 43 K区間の下限折れ線函数 92 上限折れ線函数 93 下限折れ線函数 111 フイールドコントロールステーション(FC
S)、Pはプロセス 112 オペレータステーション(OPS) 113 汎用ワークステーションReference Signs List 11 data collection unit 12 data collection generator 13 process data storage unit 14 abnormality monitoring range determination unit 15 abnormality diagnosis definition file 16 abnormality diagnosis knowledge base generator 17 simple abnormality diagnosis knowledge base 18 abnormality diagnosis unit 19 abnormality diagnosis display unit 31 process data Line function that represents the equation 42 Upper limit line function in the K section 43 Lower limit line function in the K section 92 Upper limit line function 93 Lower limit line function 111 Field control station (FC
S) and P are processes 112 Operator station (OPS) 113 General-purpose workstation
───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G05B 23/02 G06F 17/17 ──────────────────────────────────────────────────続 き Continued on front page (58) Field surveyed (Int.Cl. 7 , DB name) G05B 23/02 G06F 17/17
Claims (1)
収集して、前記プロセスが異常であるか否かを診断する
異常診断方法において、 (a)収集保存しているプロセスの計測データを順次プ
ロットして、対象とするデータの範囲をL個の区間に等
分割する手順と、 (b)第1区間の始点を起点としてL個の区間につい
て、それぞれデータ各点からの距離との差の二乗和が最
少になる折れ線函数を求める手順と、 (c)第2区間の始点を起点としてL個の区間につい
て、それぞれデータ各点からの距離との差の二乗和が最
少になる折れ線函数を求め、以下同様にして第L区間の
始点を起点とした場合まで折れ線関数を求める手順と、 (d)第1区間の始点を起点とした場合にL個の区間に
ついてそれぞれ求めたデータ各点からの距離の差の二乗
和の合計をSUM1とし、同様にしてSUM2からSU
MLまで求め、SUM1からSUMLの中で最小値にな
るものについての折れ線函数を中央折れ線函数とする手
順と、 (e)上記のようにして求めた中央折れ線函数の各区間
毎に中央折れ線函数からの標準偏差σを求める手順と、 (f)上記標準偏差σのうち最大のものを共通標準偏差
値として、経験値係数αを乗じて、上限、下限の各折れ
線函数を作成する手順と、 (g)この上限、下限の各折れ線函数で挟まれる範囲に
存在しない計測データに基づいてプロセスが異常である
と診断して信号を出力する手順とを有することを特徴と
する異常診断方法。1. An abnormality diagnosis method for sequentially collecting measurement data from a target process and diagnosing whether the process is abnormal, comprising: (a) sequentially plotting measurement data of a collected and stored process; (B) dividing the range of the target data into L sections equally; and (b) starting from the start point of the first section as a starting point.
Then, a procedure for finding a polygonal function that minimizes the sum of squares of the difference from the distance from each point of the data, and (c) starting from the start point of the second section as the starting point
Then , a polygonal function that minimizes the sum of squares of the difference from the distance from each data point is obtained , and the same applies to the L-th section in the same manner.
A procedure for obtaining a polygonal line function until the start point is set as a starting point;
The square of the difference of the distance from each point
The sum of the sum is defined as SUM1, and similarly from SUM2 to SU
ML and the minimum value from SUM1 to SUML
The line function of the object as the central line function
Common standard deviation and procedures, the largest of (f) the standard deviation σ to obtain a forward, a standard deviation σ from the center fold line function for each section of the central fold line function obtained as (e) above (G) multiplying the empirical value coefficient α as the value to create each of the upper and lower broken line functions; and (g) the process is abnormal based on the measurement data that does not exist in the range between the upper and lower broken line functions. And outputting a signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP33252293A JP3219116B2 (en) | 1993-12-27 | 1993-12-27 | Error diagnosis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP33252293A JP3219116B2 (en) | 1993-12-27 | 1993-12-27 | Error diagnosis method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH07191740A JPH07191740A (en) | 1995-07-28 |
JP3219116B2 true JP3219116B2 (en) | 2001-10-15 |
Family
ID=18255869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP33252293A Expired - Fee Related JP3219116B2 (en) | 1993-12-27 | 1993-12-27 | Error diagnosis method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3219116B2 (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8271981B2 (en) * | 2006-12-12 | 2012-09-18 | International Business Machines Corporation | Detecting an extraordinary behavior |
JP2009002651A (en) * | 2008-10-06 | 2009-01-08 | Daikin Ind Ltd | Abnormality diagnosis system |
JP5458569B2 (en) * | 2008-12-25 | 2014-04-02 | 株式会社大林組 | Construction method of steel pipe concrete columns |
JP6802025B2 (en) * | 2016-09-30 | 2020-12-16 | メタウォーター株式会社 | Valve operation monitoring method and valve operation abnormality sign detection method |
JP6724847B2 (en) | 2017-03-31 | 2020-07-15 | オムロン株式会社 | Control device, control program, control system, and control method |
CN113128598B (en) * | 2021-04-22 | 2024-04-09 | 深信服科技股份有限公司 | Sensing data detection method, device, equipment and readable storage medium |
-
1993
- 1993-12-27 JP JP33252293A patent/JP3219116B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH07191740A (en) | 1995-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3043897B2 (en) | Plant operation support device | |
JP4873985B2 (en) | Failure diagnosis device for equipment | |
KR20170125237A (en) | Plant system, and fault detecting method thereof | |
JPH03137518A (en) | Diagnosing apparatus of fault | |
JP3219116B2 (en) | Error diagnosis method | |
JPH06309584A (en) | Plant operation support device | |
JPH0793018A (en) | Method and system for diagnosing operating state | |
JPH09288512A (en) | Plant state visualized system | |
JP2672576B2 (en) | Diagnosis support system for plants and equipment | |
JPH08202431A (en) | Operation procedure guide device | |
JPH06307897A (en) | Diagnostic device for plant dispersed hierarchical system | |
JP2000214924A (en) | Plant abnormality monitoring device and method for identifying abnormality occurrence place | |
JP2000056825A (en) | Plant monitor device | |
JP2645017B2 (en) | Plant diagnostic method and apparatus | |
JPH0731076B2 (en) | Abnormality diagnosis method for rotating machinery | |
JPH07174617A (en) | System for diagnosing soundness of rotary equipment | |
JPS5882358A (en) | Front service and diagnostic service processor for motorcar service factory | |
JP7000266B2 (en) | Analysis flow creation system and analysis flow creation method | |
JPS63313208A (en) | Plant diagnostic method | |
JPH04359640A (en) | Method of diagnosing abnormal cause | |
JPH0797436B2 (en) | Power plant transient data display | |
JPH0972596A (en) | Method for diagnosing air conditioning system | |
JPH06332529A (en) | Equipment fault diagnosing method and device therefor | |
JPH04299422A (en) | Abnormality diagnosing method | |
JPH02245698A (en) | Assisting device of plant abnormality diagnosis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070810 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080810 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080810 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090810 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100810 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100810 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110810 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120810 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120810 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130810 Year of fee payment: 12 |
|
LAPS | Cancellation because of no payment of annual fees |