JP2006106870A - Process management device - Google Patents

Process management device Download PDF

Info

Publication number
JP2006106870A
JP2006106870A JP2004289052A JP2004289052A JP2006106870A JP 2006106870 A JP2006106870 A JP 2006106870A JP 2004289052 A JP2004289052 A JP 2004289052A JP 2004289052 A JP2004289052 A JP 2004289052A JP 2006106870 A JP2006106870 A JP 2006106870A
Authority
JP
Japan
Prior art keywords
failure
pattern
contour line
analysis
data storage
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.)
Granted
Application number
JP2004289052A
Other languages
Japanese (ja)
Other versions
JP4434905B2 (en
Inventor
Kenzo Yonezawa
憲造 米沢
Tomio Yamada
富美夫 山田
Tomohito Takada
智史 高田
Nobutaka Nishimura
信孝 西村
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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP2004289052A priority Critical patent/JP4434905B2/en
Publication of JP2006106870A publication Critical patent/JP2006106870A/en
Application granted granted Critical
Publication of JP4434905B2 publication Critical patent/JP4434905B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing And Monitoring For Control Systems (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To efficiently discriminate whether or not a process is abnormal from the change of spectrums obtained by carrying out the wavelet-analysis of main component scores corresponding to process result data outputted from a process side. <P>SOLUTION: When the monitor control of an air conditioning device at each object building side is started, the main component analysis of process result data stored in a process result data storing means 7 is executed by a statistical analysis means 11, and the wavelet-analysis processing and contour line patterning processing of main component scores "Z<SB>i</SB>" obtained by the main component analysis is carried out by a process status deciding means 9 so that the area of a land can be calculated, and when this exceeds a preliminarily set threshold, it is decided that the process side is abnormal. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は、産業用や民生用の各種設備における異常や運転不能状況を予測したり、故障発生の原因を推定するプロセス管理装置に係わり、特にBEMS装置の数多くある機能のうち、故障検知、故障予知、故障診断を行うプロセス管理装置に関する。   The present invention relates to a process management apparatus for predicting abnormalities and inoperability in various industrial and consumer facilities, and estimating the cause of a failure. In particular, among many functions of a BEMS apparatus, failure detection, failure The present invention relates to a process management apparatus that performs prediction and failure diagnosis.

近年、コンピュータ処理技術の進歩に伴い、産業用や民生用の各種設備などを撮影して得られた画像データを処理して、これら各種設備の異常有無を判定し、設備の健全性を判定する方法が開発されている(特許文献1参照)。   In recent years, with the advancement of computer processing technology, image data obtained by photographing various industrial and consumer facilities is processed to determine the presence or absence of abnormalities in these various facilities and to determine the soundness of the facilities. A method has been developed (see Patent Document 1).

しかしながら、この方法では、画像処理によって監視対象の異常有無を判定していることから、ビル監視などに適用した場合、画像を通してビル各部の異常を把握できない限りは、ビル内各部における湿度データ、温度データ、弁開度データなどのプロセスデータが異常であっても、異常発生を検出できないというという問題があった。   However, in this method, the presence / absence of an abnormality to be monitored is determined by image processing. Therefore, when applied to building monitoring, etc., humidity data and temperature at each part in the building can be obtained unless the abnormality of each part of the building can be grasped through the image. Even if process data such as data and valve opening data is abnormal, there is a problem that the occurrence of abnormality cannot be detected.

このため、近年に建築の大規模化したビルでは、空調システムなどのビル設備に発生した異常の発見や処置が遅くなり、室内環境悪化やエネルギー浪費を引き起こす可能性が増大している。   For this reason, in a building with a large-scale building in recent years, the discovery and treatment of abnormalities occurring in building equipment such as an air conditioning system is slowed down, and the possibility of causing deterioration of the indoor environment and energy waste is increasing.

そこで、これらの問題を解決し、かつ高い省エネ効果が見込まれるシステムとして、BEMS(Bui1ding and Energy Management System:日本語ではビル・環境・エネルギー管理システム)装置が改めて今、注目されている。   Therefore, BEMS (Bui1ding and Energy Management System: Japanese Building, Environment, Energy Management System) devices are now attracting attention as a system that solves these problems and is expected to have a high energy saving effect.

BEMS装置は、ビル経営管理からエネルギー管理まで、ビルに関する全てのビル管理項目あるいは複数のビル管理項目を一括して管理するものであり、事務所や店舗など業務用ビルの省エネ対策として重要なものになっている。   BEMS equipment manages all building management items or multiple building management items from building management to energy management in a lump, and is important as an energy-saving measure for commercial buildings such as offices and stores. It has become.

この際、監視対象となるビル内に設けられたビル設備の異常有無を判定する方法として、ビル設備から出力されるプロセスデータをフーリェ変換してスペクトルを解析し、故障有無を判定する方法、またはシステム同定手法などを用いてビル設備の特性モデルを推定し、この特性モデルから出力されるプロセスデータと、ビル設備から出力されるプロセスデータとを比較し、故障有無を判定する方法などが使用される。
特開2001−311723号公報
At this time, as a method of determining the presence or absence of abnormality of the building equipment provided in the building to be monitored, a method of analyzing the spectrum by performing Fourier transform on the process data output from the building equipment, and determining the presence or absence of failure, or A method that estimates the characteristic model of a building facility using a system identification method, etc., compares the process data output from this characteristic model with the process data output from the building facility, and determines the presence or absence of a failure is used. The
JP 2001-31723 A

ところで、種々のビル設備の中で、暖房や冷房を行う空調システムは建物にとって欠かせないものになっており、ビル各部の温度や湿度、各空調機の消費電力量などを監視して最適な運転状況にすることは、極めて重要なことである。   By the way, in various building facilities, the air conditioning system that performs heating and cooling is indispensable for the building, and the optimum temperature is monitored by monitoring the temperature and humidity of each part of the building and the power consumption of each air conditioner. It is extremely important to get into the driving situation.

しかしながら、大規模なビルでは、出入り口ドアの開閉状況などにより、ビル各部の温度、湿度が急激に変化することがあり、フーリェ変換によって、時々刻々、変化するスペクトルを求めたり、システム同定手法によって、変動する特性モデルを求めたりすることができないという問題があった。   However, in large-scale buildings, the temperature and humidity of each part of the building may change suddenly depending on the opening and closing status of the entrance door, etc. The spectrum that changes from moment to moment due to Fourier transformation, or by system identification techniques, There is a problem that a characteristic model that fluctuates cannot be obtained.

また、このような大規模ビルでは、多数の空調機を備えた空調システムを各階毎に数台設置しなければならない。各空調機のいずれかに不具合が生じると、空調システムが正常な動作をしなくなり、居住者に快適な空調環境を提供できなくなる。さらに省エネ制御などが導入されていても、センサは通常、制御に最低限必要なものしか設置されず、故障及びその要因検知用のものは設置されないのが一般的である。   In such large-scale buildings, several air conditioning systems equipped with a large number of air conditioners must be installed on each floor. If any of the air conditioners malfunctions, the air conditioning system will not operate normally, and a comfortable air conditioning environment cannot be provided to the occupants. Further, even if energy saving control or the like is introduced, usually, only sensors necessary for the control are usually installed, and those for detecting a failure and its factor are generally not installed.

このため、設備機器の不具合が発生しても、これが発見されず、そのまま運転され、機器効率が大幅に低下した状態が継続されて、エネルギーが無駄に使われることが多い。   For this reason, even if a malfunction occurs in the equipment, it is not discovered and it is operated as it is, and the state in which the equipment efficiency is greatly reduced is continued and energy is often wasted.

そこで、近年、建物に設置された各設備の不具合を迅速に見つけるために、多数の建物と、センタとを通信回線で結び、センタに配置された監視員によって、各建物の設備を集中的に監視し、各建物側に管理装置が設置されていないときでも、センタ側で、各建物側に設けられた各設備の異常有無を監視できるようにした監視システムを導入することが多い。   Therefore, in recent years, in order to quickly find defects in each facility installed in a building, many buildings and the center are connected by a communication line, and the facilities in each building are concentrated by a supervisor placed in the center. In many cases, a monitoring system is installed that can monitor the presence or absence of abnormality of each facility provided on the building side on the center side even when a management device is not installed on each building side.

しかしながら、このような監視システムでは、センタ側に膨大な量のデータが集中することから、人手の監視による不具合の検知には限界があり、膨大な量のデータを自動的に解析して、効率良く故障の有無を判定することができる装置の開発が強く望まれていた。   However, in such a monitoring system, a huge amount of data is concentrated on the center side, so there is a limit to the detection of malfunctions due to manual monitoring, and a huge amount of data is automatically analyzed to improve efficiency. There has been a strong demand for the development of a device that can determine the presence or absence of a failure.

本発明は、上記の事情に鑑み、装置全体の構成を簡素化し、処理時間を短縮しつつ、プロセス側から出力されるプロセス実績データをウェーブレット解析して得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定でき、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減でき、保守費を低減できるプロセス管理装置を提供することを目的としている。   In view of the above circumstances, the present invention simplifies the configuration of the entire apparatus, shortens the processing time, and determines whether the process is normal from the change in spectrum obtained by wavelet analysis of process performance data output from the process side. The purpose is to provide a process management device that can efficiently estimate and estimate system performance, perform system performance verification, failure detection, failure prediction, and failure diagnosis, reduce energy consumption, and reduce maintenance costs. Yes.

上記の目的を達成するために本発明は、プロセスを時系列的に計測して得られる複数のプロセス実績データを蓄積するプロセス実績データ蓄積手段と、このプロセス実績データ蓄積手段に蓄積されている各プロセス実績データをウェーブレット解析するとともに、この解析結果を3次元の等高線パターン化して異常の有無を判定するプロセス状態判定手段とを備えたことを特徴としている。   In order to achieve the above object, the present invention provides a process result data storage means for storing a plurality of process result data obtained by measuring a process in time series, and each of the process result data storage means stored in the process result data storage means. It is characterized in that it includes a process state determination means for performing wavelet analysis of the process performance data and determining the presence or absence of abnormality by converting the analysis result into a three-dimensional contour pattern.

本発明によれば、装置全体の構成を簡素化して、処理時間を短縮しつつ、プロセス側から出力されるプロセス実績データをウェーブレット解析して得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定でき、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減でき、保守費を低減することができる。   According to the present invention, it is possible to efficiently determine whether a process is normal from a change in spectrum obtained by wavelet analysis of process performance data output from the process side while simplifying the configuration of the entire apparatus and reducing processing time. Discrimination and estimation can be performed, system performance verification, failure detection, failure prediction, and failure diagnosis can be performed with high accuracy, energy consumption can be reduced, and maintenance costs can be reduced.

《本発明で使用されるウェーブレット解析法の説明》
まず、本発明によるプロセス管理装置の具体的な説明に先立ち、本発明によるプロセス管理装置で使用されるウェーブレット解析法について説明する。
<< Description of Wavelet Analysis Method Used in the Present Invention >>
First, prior to specific description of the process management apparatus according to the present invention, a wavelet analysis method used in the process management apparatus according to the present invention will be described.

従来の解析法では、監視対象から出力されるプロセスデータなどの信号データ“x(t)”に対し、下記に示すフーリェ変換式を用いて、周波数特性などを解析する。
In the conventional analysis method, frequency characteristics and the like are analyzed by using the following Fourier transform formula for signal data “x (t)” such as process data output from the monitoring target.

これに対して本発明では、非定常信号を解析する手法として、近年、注目されている、下記に示す(連続)ウェーブレット変換式を用いて、周波数特性を解析する。
On the other hand, in the present invention, frequency characteristics are analyzed using the following (continuous) wavelet transform equation which has been attracting attention as a technique for analyzing unsteady signals in recent years.

ここで、“φ(・)”はウェーブレット解析で、マザーウェーブレットと呼ばれる基底関数であり、“φ((t−b)/a)”は、マザーウェーブレット“φ(t)”を“b”だけ平行移動し、“a”だけ、スケール(拡大縮小)したものである。   Here, “φ (•)” is a wavelet analysis and is a basis function called a mother wavelet, and “φ ((t−b) / a)” is a mother wavelet “φ (t)” that is only “b”. Translated and scaled (enlarged / reduced) by “a”.

この際、スケールパラメータ“a”に対応して、マザーウェーブレット“φ(t)”の幅が“a”倍になり、“1/a”が周波数に対応することから、スケールパラメータ“a”が大きいとき、基底関数が拡がり、またスケールパラメータ“a”が小さいとき、基底関数が狭まる。   At this time, the width of the mother wavelet “φ (t)” is “a” times corresponding to the scale parameter “a”, and “1 / a” corresponds to the frequency. When it is large, the basis function is expanded, and when the scale parameter “a” is small, the basis function is narrowed.

一般的に使用されいてるフーリェ変換式は、(連続)ウェーブレット変換式の基底関数を、
φ(t)=e−jωt …(3)
a=ω−1 …(4)
b=0 …(5)
とした場合に相当し、その基底関数は時間領域で無限の拡がりを持つ。
The commonly used Fourier transform is the basis function of the (continuous) wavelet transform,
φ (t) = e− jωt (3)
a = ω−1 (4)
b = 0 (5)
The basis function has an infinite extent in the time domain.

このため、フーリェ変換によって得られたスペクトルは、周波数軸のみの1次関数となって、信号の時間的情報が失われてしまい、観測データのどの部分の特徴かという時間依存性を判別することができない。   For this reason, the spectrum obtained by the Fourier transform becomes a linear function only on the frequency axis, the time information of the signal is lost, and the time dependency of which part of the observation data is characteristic is determined. I can't.

これに対し、ウェーブレット変換では、例えば下式に示すガボール関数と呼ばれる時間的に局所した関数を用いる。
In contrast, in the wavelet transform, for example, a temporally local function called a Gabor function shown in the following equation is used.

これにより、ウェーブレット変換によって得られたスペクトルは、周波数軸と、時間軸とに対する2次関数となり、信号の時間的情報を持つことから、観測データのどの部分の特徴かという時間依存性を判別することができる。   As a result, the spectrum obtained by the wavelet transform becomes a quadratic function with respect to the frequency axis and the time axis, and has temporal information of the signal, so that the time dependence of which part of the observation data is characteristic is determined. be able to.

なお、ウェーブレット変換理論に関する詳細は、例えば新井康平著「ウェーブレット解析の基礎理論」、森北出版株式会社(2000)などで、紹介されている。   Details on the wavelet transform theory are introduced in, for example, Kohei Arai “Basic Theory of Wavelet Analysis”, Morikita Publishing Co., Ltd. (2000).

しかしながら、このようなウェーブレット解析法を用いて、ビル内に設けられた各センサから出力されるプロセスデータを全て解析すると、解析結果が複雑になり過ぎて、故障の検知、予知、診断などが難しくなり、故障の発見が遅れてしまう恐れがある。   However, if all of the process data output from each sensor provided in the building is analyzed using such a wavelet analysis method, the analysis result becomes too complex, making it difficult to detect, predict, diagnose, etc. Therefore, there is a risk that the discovery of the failure will be delayed.

そこで、本発明では、ビル内に設けられる各センサの数が多いときには、主成分分析法を用いて、各センサから出力されるプロセスデータの主成分を抽出し、この抽出処理で得られた主成分をウェーブレット解析法で解析することにより、故障の早期発見、予知、診断などを実現している。   Therefore, in the present invention, when the number of sensors provided in the building is large, the principal component of the process data output from each sensor is extracted using the principal component analysis method, and the main component obtained by this extraction processing is extracted. By analyzing the components using the wavelet analysis method, early detection, prediction and diagnosis of failures are realized.

《第1の実施形態》
<実施形態の構成説明>
図1は本発明によるプロセス管理装置の第1の実施形態の全体構成を示すブロック図である。
<< First Embodiment >>
<Description of Configuration of Embodiment>
FIG. 1 is a block diagram showing the overall configuration of a first embodiment of a process management apparatus according to the present invention.

この図に示すプロセス管理装置1は、対象建物側(デパートビル2、事務所ビル3、病院ビル4等)に設置された各計測監視制御手段5とネットワーク6を介して接続され、各計測監視制御手段5で収集された検知結果(プロセス実績データ)を取り込んで蓄積するプロセス実績データ蓄積手段7と、主成分分析(PCA:Principa1 Components anlysis)を用いて、プロセス実績データ蓄積手段7に蓄積されているプロセス実績データを解析して主成分を生成する統計解析手段8と、統計解析手段8で得られた主成分をウェーブレット解析して監視制御対象になっている対象建物側に設置された各空調装置が正常か、異常かを判定するプロセス状態判定手段9と、このプロセス状態判定手段9で得られた等高線パターンなどの解析データなどを蓄積する解析データ蓄積手段10とを備えている。   The process management apparatus 1 shown in this figure is connected to each measurement monitoring control means 5 installed on the target building side (department building 2, office building 3, hospital building 4, etc.) via the network 6, and each measurement monitoring. The process result data storage means 7 that captures and accumulates the detection results (process result data) collected by the control means 5 and the principal component analysis (PCA: Principa1 Components anlysis) are stored in the process result data storage means 7. Statistical analysis means 8 for generating the principal component by analyzing the actual process performance data, and each installed on the target building side that is subject to monitoring control by wavelet analysis of the principal component obtained by the statistical analysis means 8 The process state determination means 9 for determining whether the air conditioner is normal or abnormal, and the analysis data such as the contour line pattern obtained by the process state determination means 9 are accumulated. Analysis data storage means 10.

プロセス実績データ蓄積手段7は、各デパートビル2、各事務所ビル3、各病院ビル4などの対象建物側に設置された各計測監視制御手段5から空調装置のエネルギー消費量、その他のプロセス値(外気温、室温、室内湿度等)などを含むプロセス実績データが出力されたとき、ネットワーク6を介してこれを取り込んで、プロセス状態判定手段9の判定結果などに基づき、正常時のプロセス実績データと、故障時のプロセス実績データとに分けて蓄積しながら、蓄積しているプロセス実績データを統計解析手段8に供給する。   The process result data storage means 7 includes the energy consumption of the air conditioner and other process values from each measurement monitoring control means 5 installed on the target building side such as each department store building 2, each office building 3, each hospital building 4, etc. When process performance data including (outside temperature, room temperature, indoor humidity, etc.) is output, it is fetched via the network 6, and based on the determination result of the process state determination means 9, the normal process performance data The accumulated process result data is supplied to the statistical analysis means 8 while being divided into process result data at the time of failure.

統計解析手段8は、各種建築物、負荷条件などの各パラメータを用いて、プロセスが正常なときのシミュレーションと各種故障が発生したときのシミュレーションを行うプロセスシミュレーション手段11と、プロセスシミュレーション手段11で得られたシミュレーション結果を蓄積するシミュレーションデータ蓄積手段12と、主成分分析を用いて、シミュレーションデータ蓄積手段12に蓄積されているシミュレーション結果、またはプロセス実績データ蓄積手段7に蓄積されているプロセス実績データを統計解析して、正常状態の特性値と、異状状態の特性値との差異を早期に、容易に抽出するのに必要な前段階の処理を行なう主成分分析手段13とを備えている。   The statistical analysis means 8 is obtained by the process simulation means 11 and the process simulation means 11 which perform simulation when the process is normal and simulation when various failures occur, using various parameters such as various buildings and load conditions. The simulation data storage means 12 for storing the simulation results and the simulation results stored in the simulation data storage means 12 or the process result data stored in the process result data storage means 7 using principal component analysis. Principal component analysis means 13 is provided that performs statistical analysis and performs a pre-stage process necessary for early and easy extraction of the difference between the characteristic value in the normal state and the characteristic value in the abnormal state.

そして、各対象建物側の空調装置などの監視制御を開始する前に、各種建築物、負荷条件などの各パラメータを用いて、プロセスが正常な状態であるときのシミュレーションと、各種故障状態のときのシミュレーションとを行わせてシミュレーション結果を蓄積し、このシミュレーション結果を主成分分析して、正常状態の特性値と、異状状態の特性値との差異を早期に、容易に抽出するのに必要なプロセス変数結合重みパラメータを求めるとともに、このプロセス変数結合重みパラメータを用いて、プロセスが正常な状態であるときのシミュレーション結果と、各種故障状態のときのシミュレーション結果とを主成分分析し、各主分析スコアをプロセス状態判定手段9に供給する。この後、各対象建物側の空調装置などの監視制御を行っているとき、プロセス変数結合重みパラメータを用いて、プロセス実績データ蓄積手段7を介して供給されるプロセス実績データを主成分分析し、主分析スコアをプロセス状態判定手段9に供給する。そして、プロセス管理装置1の運転期間が長くなり、プロセス実績データ蓄積手段7に正常時のプロセス実績データ、故障時のプロセス実績データが十分に蓄積された後は、これらのプロセス実績データを用いてプロセス変数結合重みパラメータを更新する。   Before starting monitoring and control of the air conditioners on the target building side, using various parameters such as various buildings and load conditions, simulations when the process is in a normal state and various failure states The simulation results are accumulated and the simulation results are accumulated, and the simulation results are analyzed by principal component analysis, which is necessary to quickly and easily extract the difference between the normal state characteristic values and the abnormal state characteristic values. In addition to obtaining process variable connection weight parameters, using this process variable connection weight parameter, principal component analysis of simulation results when the process is in a normal state and simulation results when various fault conditions occur, and each main analysis The score is supplied to the process state determination means 9. After this, when performing monitoring control of the air conditioners and the like on each target building side, the process result data supplied via the process result data storage means 7 is principally analyzed using the process variable coupling weight parameter, The main analysis score is supplied to the process state determination means 9. Then, after the operation period of the process management device 1 becomes long and the process result data storage means 7 has sufficiently stored the process result data at the normal time and the process result data at the time of failure, these process result data are used. Update the process variable binding weight parameter.

プロセス状態判定手段9は、前記(2)式に示す(連続)ウェーブレット変換式を用いて、統計解析手段8の主成分分析手段13で得られた各主分析スコアをウェーブレット解析するウェーブレット解析手段14と、ウェーブレット解析手段14で得られた解析結果を等高線パターン化して、解析データ蓄積手段10に蓄積する等高線パターン化手段15と、解析データ蓄積手段15に蓄積されている等高線パターンに基づき、海面に相当する振幅の高さを求めるとともに、この振幅の高さに基づき、等高線パターン化手段15で得られた等高線パターンを海面部分と陸地部分とに分けて陸地部分の面積を演算する陸地面積演算手段16と、解析データ蓄積手段10に蓄積されている等高線パターンなどに基づき、表示装置上に周波数軸、時間軸を持つ3次元の等高線パターンなどを表示しながら、陸地面積演算手段16で得られた陸地面積が予め設定されたしきい値を超えているとき、監視制御対象となっているプロセスが異常となっていると判定して、判定結果をプロセス実績データ蓄積手段7、解析データ蓄積手段10などに供給するとともに、表示装置上にプロセス異常を表示する正常/異常判定手段17とを備えている。   The process state determination means 9 uses the (continuous) wavelet transform equation shown in the above equation (2) to perform wavelet analysis on each main analysis score obtained by the principal component analysis means 13 of the statistical analysis means 8. The contour of the analysis result obtained by the wavelet analysis means 14 is converted into a contour pattern and stored in the analysis data storage means 10, and on the sea surface based on the contour line pattern stored in the analysis data storage means 15. A land area calculating means for calculating the corresponding amplitude height and dividing the contour line pattern obtained by the contour line patterning means 15 into a sea surface portion and a land portion based on the amplitude height and calculating the area of the land portion. 16 and the contour line pattern stored in the analysis data storage means 10, etc. When the land area obtained by the land area calculating means 16 exceeds a preset threshold value while displaying a three-dimensional contour line pattern having an axis, the process being monitored and controlled is abnormal. In addition to supplying the determination result to the process result data storage means 7, the analysis data storage means 10 and the like, a normal / abnormality determination means 17 for displaying a process abnormality on the display device is provided.

そして、各対象建物側の空調装置などの監視制御を開始する前に、統計解析手段8の主成分分析手段13で得られた各主分析スコアをウェーブレット解析して、解析結果を等高線パターン化し、解析データ蓄積手段10に蓄積した後、解析データ蓄積手段10に蓄積されている正常時の等高線パターン、異常時の等高線パターンから、海面に相当する振幅の高さを求める。この後、各対象建物側に設けられた空調装置などの監視制御を行っているとき、統計解析手段8の主成分分析手段13から出力されるプロセス実績データの主成分スコアをウェーブレット解析して、等高線パターンを求め、解析データ蓄積手段10に蓄積する。また、この動作と並行し、海面に相当する振幅の高さを用いて、等高線パターンに含まれる陸地部分(海面に相当する振幅の高さより大きい値の部分を)の面積を求め、この面積が予め設定されているしきい値を超えたとき、プロセスに異常が発生したと判定し、判定結果をプロセス実績データ蓄積手段7、解析データ蓄積手段10などに供給するとともに、表示装置上にプロセス異常を表示する。   And before starting monitoring control of the air conditioner etc. of each target building side, each main analysis score obtained by the principal component analysis means 13 of the statistical analysis means 8 is subjected to wavelet analysis, and the analysis result is converted into a contour pattern, After accumulating in the analysis data accumulating means 10, the height of the amplitude corresponding to the sea surface is obtained from the normal contour line pattern and the abnormal contour line pattern accumulated in the analysis data accumulating means 10. After that, when performing monitoring control such as an air conditioner provided on each target building side, wavelet analysis of the principal component score of the process result data output from the principal component analysis means 13 of the statistical analysis means 8 is performed, A contour line pattern is obtained and stored in the analysis data storage means 10. In parallel with this operation, the height of the amplitude corresponding to the sea surface is used to determine the area of the land portion (the portion having a value larger than the amplitude height corresponding to the sea surface) included in the contour pattern. When a preset threshold value is exceeded, it is determined that an abnormality has occurred in the process, the determination result is supplied to the process result data storage means 7, the analysis data storage means 10, etc., and the process abnormality is displayed on the display device. Is displayed.

解析データ蓄積手段10は、プロセス状態判定手段9の等高線パターン化手段15から新たな等高線パターンが出力される毎に、これを取り込み、1日分単位の時系列データとして蓄積しながら、プロセス状態判定手段9の陸地面積演算手段16から読み出し要求が出されたとき、蓄積している等高線パターンを読み出し、プロセス状態判定手段9の陸地面積演算手段16に供給する。また、プロセス状態判定手段9の正常/異常判定手段17から出力される正常/異常判定結果に基づき、プロセス状態判定手段9の等高線パターン化手段15から出力される等高線パターンを正常か、異常かに分類して、蓄積する。   The analysis data storage means 10 captures each time a new contour line pattern is output from the contour line patterning means 15 of the process state determination means 9 and stores it as time series data in units of one day while determining the process state. When a read request is issued from the land area calculation means 16 of the means 9, the accumulated contour line pattern is read and supplied to the land area calculation means 16 of the process state determination means 9. Whether the contour line pattern output from the contour line patterning unit 15 of the process state determining unit 9 is normal or abnormal is determined based on the normal / abnormal determination result output from the normal / abnormal determining unit 17 of the process state determining unit 9. Classify and accumulate.

<実施形態の動作説明>
次に、図1に示すブロック図、図2に示すフローチャートを参照しながら、各対象建物に設けられた空調装置などの監視制御を開始する前の動作と、各対象建物側の空調装置などの監視制御を行っているときの動作とに分けて、プロセス管理装置1の動作を説明する。
<Description of Operation of Embodiment>
Next, referring to the block diagram shown in FIG. 1 and the flowchart shown in FIG. 2, the operation before starting the monitoring control of the air conditioner provided in each target building, the air conditioner on the target building side, etc. The operation of the process management apparatus 1 will be described separately from the operation when monitoring control is performed.

[各対象建物側の空調装置などの監視制御を開始する前の動作]
まず、各対象建物側に設けられた空調装置などの監視制御を開始する前に、図2のフローチャートに示すシミュレーション処理が開始されて、統計解析手段8のプロセスシミュレーション手段11に入力されている各対象建物に関する各データ、各負荷条件データなどが時刻毎に整理される。
[Operation before starting monitoring and control of air conditioners, etc. at each target building]
First, before starting monitoring control of an air conditioner or the like provided on each target building side, the simulation process shown in the flowchart of FIG. 2 is started, and each of the processes input to the process simulation means 11 of the statistical analysis means 8 is started. Each data related to the target building, each load condition data, and the like are arranged for each time.

この後、プロセスシミュレーション手段11によって、予め入力されている対象建物側に関する各情報が読み出されて、建物側に設けられた部屋の大きさ、熱容量、壁面積、壁の熱貫流率等のデータが設定され(ステップS1)、各負荷条件データなどに基づき、空調負荷に関連する気象要素である外気温、外気湿度、日射量が事前に設定された後(ステップS2)、人の在室率変化に応じたシミュレーション時間分だけ、日射量を用いて窓ガラス面からの取得熱量、人の発熱量やOA等の電気機器、照明の発熱量などが予め計算されて、熱負荷データが設定される(ステップS3)。   After that, the process simulation means 11 reads information on the target building side that has been input in advance, and data such as the size of the room, the heat capacity, the wall area, and the thermal conductivity of the wall provided on the building side. Is set (step S1), and based on each load condition data and the like, the outside air temperature, the outside air humidity, and the amount of solar radiation, which are meteorological elements related to the air conditioning load, are set in advance (step S2), For the simulation time corresponding to the change, the amount of heat obtained from the window glass surface using the amount of solar radiation, the amount of heat generated by humans, electrical equipment such as OA, the amount of heat generated by lighting, etc. are calculated in advance, and heat load data is set. (Step S3).

この後、プロセスシミュレーション手段11によって、サンプリング制御周期(シミュレータ計算周期の整数倍)であるかどうかがチェックされ、サンプリング周期であるとき(ステップS4)、DDC(Direct Digita1 Contro11er)制御アルゴリズム方式で、室温設定値と室温との偏差、湿度設定値と室内湿度との偏差などに基づき、PIDパラメータ(速度形)が演算される(ステップS5)。   Thereafter, the process simulation means 11 checks whether or not the sampling control period is an integer multiple of the simulator calculation period, and when it is the sampling period (step S4), the room temperature is measured by the DDC (Direct Digita1 Contro11er) control algorithm method. Based on the deviation between the set value and room temperature, the difference between the humidity set value and room humidity, etc., the PID parameter (speed type) is calculated (step S5).

次いで、プロセスシミュレーション手段11によって、冷却コイル、加熱コイル、加湿器の順で熱交換器の理論計算式や、空気線図の計算式が用いられて、空調機内での空気の各状態点(温度、湿度、エンタルピーなど)が求められるとともに(ステップS6)、ビル内の部屋の壁が外壁(窓あり)、床、内壁、天井の部分に分けられて、伝熱方程式(伝熱モデル)で、熱の収支が記述され(ステップS7)、差分解法による室温と壁温との動特性変化が計算され、再循環空気と導入外気の混合計算により空調機入口での空気の状態点が求められる(ステップS8)。   Next, the process simulation means 11 uses the theoretical calculation formula of the heat exchanger and the calculation formula of the air diagram in the order of the cooling coil, the heating coil, and the humidifier, and each state point (temperature) of the air in the air conditioner , Humidity, enthalpy, etc.) (step S6), the wall of the room in the building is divided into the outer wall (with windows), the floor, the inner wall, and the ceiling, and the heat transfer equation (heat transfer model) The heat balance is described (step S7), the dynamic characteristic change between the room temperature and the wall temperature is calculated by the differential decomposition method, and the air state point at the air conditioner inlet is obtained by the mixing calculation of the recirculated air and the introduced outside air ( Step S8).

以下、プロセスシミュレーション手段11によって、時刻データが更新されながら、シミュレーション終了時刻になるまで、上述した計算が実行されて、プロセスが正常な状態であるときのシミュレーションと、各種故障状態のときのシミュレーションとが行われて、シミュレーションデータ蓄積手段12にプロセスが正常な状態であるときのプロセス実績データと、各種故障状態のときのプロセス実績データとが蓄積される(ステップS4〜S10)。   Hereinafter, the above-described calculation is executed by the process simulation means 11 until the simulation end time is reached while the time data is being updated, and the simulation when the process is in a normal state and the simulation when the various fault states occur. Is performed, and the process result data when the process is in a normal state and the process result data when various failure states are stored in the simulation data storage unit 12 (steps S4 to S10).

また、この動作と並行し、各対象建物側の空調装置などの監視制御を開始する前に、主成分分析手段13によって、シミュレーションデータ蓄積手段12に蓄積されているプロセスが正常なときのプロセス実績データと、各種故障が発生したときのプロセス実績データとが分析されて、正常状態の特性値と、異状状態の特性値との差異を早期に、容易に抽出するのに必要なプロセス変数結合重みパラメータが求められる。   Further, in parallel with this operation, the process results when the process accumulated in the simulation data accumulating unit 12 by the principal component analyzing unit 13 is normal before the monitoring control of the air conditioner etc. of each target building is started. Data and process performance data when various types of failures occur are analyzed, and process variable connection weights necessary to quickly and easily extract the difference between the normal state characteristic value and the abnormal state characteristic value A parameter is required.

例えば、プロセス実績データが変数の数がp個“X、X、……、x”、観測対象の数がn個の多変量データであり、この多変量データを用いて、p個より少ないm個の新しい変数“Z、Z、……、Z”を作る場合、下記に示す条件の下で、プロセス変数結合重みパラメータ“a11、a12、…、amp”が算出される。 For example, the process performance data is multivariate data with the number of variables “X 1 , X 2 ,..., X p ” and the number of observation objects is n, and using this multivariate data, p When creating fewer m new variables “Z 1 , Z 2 ,..., Z m ”, the process variable coupling weight parameters “a 11 , a 12 ,. Calculated.

Var(Z)≧Var(Z)≧…≧Var(Z) …(7)
i1 +ai2 +…十aip =1 …(8)
但し、i:i=1、2、…、m
Var(Z):分散
この際、新しい変数“Z”は新しい変数“Z”とは独立になるようにされる。また新しい変数“Z”は新しい変数“Z”および新しい変数“Z”とは独立になるようにされる。以下、新しい変数“Z”から新しい変数“Z”まで同様にされる。
Var (Z 1 ) ≧ Var (Z 2 ) ≧… ≧ Var (Z m ) (7)
a i1 2 + a i2 2 + ... 10 a ip 2 = 1 (8)
Where i: i = 1, 2,..., M
Var (Z i ): Variance At this time, the new variable “Z 1 ” is made independent of the new variable “Z 2 ”. The new variable “Z 3 ” is made independent of the new variable “Z 1 ” and the new variable “Z 2 ”. Hereinafter, the same process is performed from the new variable “Z 4 ” to the new variable “Z m ”.

このような条件の下で、プロセス変数結合重みパラメータ“a11、a12、…、amp”を求めることは、元の変数“X、X、…、X”の分散共分散行列の固有値と、固有ベクトルを計算することに帰着し、プロセス変数結合重みパラメータ“a11、a12、…、amp”は、固有ベクトルを求めることに他ならない。 Under such conditions, the process variable connection weight parameter "a 11, a 12, ... , a mp" to seek the original variables "X 1, X 2, ... , X p" covariance matrix of The process variable coupling weight parameters “a 11 , a 12 ,..., A mp ” are none other than finding the eigenvectors.

但し、元の変数“X、X、…、X”の単位と、新しい変数“Z、Z、……、Z”の単位とが異なることから、主成分分析の前に、元の変数“X、X、…、X”が正規化(標準化)される。 However, the original variable "X 1, X 2, ... , X p" and the units of the new variable "Z 1, Z 2, ...... , Z m" from that and units of different, in front of the principal component analysis , The original variables “X 1 , X 2 ,..., X p ” are normalized (standardized).

すなわち、k番目の変数の対象全データの平均値が“M”、標準偏差が“S”であるとき、下のように各変数“Xkj”が正規化される。 That is, when the average value of all target data of the kth variable is “M k ” and the standard deviation is “S k ”, each variable “X kj ” is normalized as follows.

kj→(Xkj−M)/S …(9)
但し、j:j=1〜n(nはデータの個数)
次いで、主成分分析手段によって、プロセス変数結合重みパラメータ“a11、a12、…、amp”が使用されて、次式に示す演算式で、正常な状態であるときプロセス実績データ(または、各種故障状態のときのプロセス実績データ)を構成している元の変数“X、X、……、x”が結合されて、主成分となる新しい変数“Z、Z、……、Z”が作成される。
X kj → (X kj −M k ) / S k (9)
However, j: j = 1 to n (n is the number of data)
Next, the process variable connection weight parameters “a 11 , a 12 ,..., A mp ” are used by the principal component analysis means, and the operation result data (or The original variables “X 1 , X 2 ,..., X p ” constituting the process result data at the time of various fault conditions are combined to form new variables “Z 1 , Z 2 ,. ..., Z m "is created.

(t)=a11・X(t)+a12・X(t)+…十a1p・X(t)…(10)
(t)=a21・X(t)+a22・X(t)+…十a2p・X(t)…(11)


(t)=am1・X(t)+am2・X(t)+…十amp・X(t)…(12)
但し、amp:プロセス変数結合重みパラメータ
t:時間
この場合、新しい変数(主成分スコア)“Z、Z、……、Z”は、次のような性質を持つものとする。
Z 1 (t) = a 11 · X 1 (t) + a 12 · X 2 (t) + ... Ten a 1p · X p (t) ... (10)
Z 2 (t) = a 21 · X 1 (t) + a 22 · X 2 (t) + ... Ten a 2p · X p (t) ... (11)
...
...
Z m (t) = a m1 · X 1 (t) + a m2 · X 2 (t) + ... 10 amp · X p (t) (12)
However, amp : process variable connection weight parameter t: time In this case, it is assumed that new variables (principal component scores) “Z 1 , Z 2 ,..., Z m ” have the following properties.

新しい変数“Z”は元の変数“X〜X”の情報が最大限集約されるようにする。また、新しい変数“Z”は元の変数“X〜X”の情報が新しい変数“Z”の次に出来るだけ多く集約されるようにする。以下、新しい変数“Z”から新しい変数“Z”まで、同様に考える。 The new variable “Z 1 ” allows the information of the original variables “X 1 to X p ” to be aggregated to the maximum extent. Further, the new variable “Z 2 ” is configured so that the information of the original variables “X 1 to X p ” is aggregated as much as possible after the new variable “Z 1 ”. Hereinafter, the same applies from the new variable “Z 3 ” to the new variable “Z m ”.

そして、主成分分析手段13によって、シミュレーションデータ蓄積手段12に蓄積されているプロセスが正常な状態であるときの全プロセス実績データに対する主成分分析処理と、各種故障状態のときの全プロセス実績データに対する主成分分析処理とが完了したとき、主成分分析で得られた正常時の主成分スコア“Z”と、故障時の主成分スコア“Z”とがプロセス状態判定手段9に供給される。 Then, the principal component analysis means 13 performs principal component analysis processing on all process result data when the process accumulated in the simulation data accumulation means 12 is in a normal state, and applies to all process result data when various failure states occur. when the principal component analysis process and is completed, the main component score "Z m" obtained at the time of normal in the principal component analysis, principal component scores "Z m" at the time of failure and is supplied to the process state judging means 9 .

また、上述した動作と並行し、プロセス状態判定手段9のウェーブレット解析手段14、等高線パターン化手段15によって、主成分分析手段13から出力される正常時の主成分スコア“Z”と、故障時の主成分スコア“Z”とがウェーブレット解析されるとともに、等高線パターン化されて、正常時の等高線パターンと、異常時の等高線パターンとが生成され、解析データ蓄積手段10に蓄積される。そして、ウェーブレット解析手段14、等高線パターン化手段15によって、主成分分析手段13から出力される全ての主成分スコア“Z”のウェーブレット解析処理、等高線パターン化処理が終了したとき、陸地面積演算手段16によって、解析データ蓄積手段10に蓄積されている正常時の等高線パターンと、異常時の等高線パターンとに基づき、海面に相当する振幅の高さを求められて、記憶される。 In parallel with the above-described operation, the normal component score “Z m ” output from the principal component analysis unit 13 by the wavelet analysis unit 14 and the contour line patterning unit 15 of the process state determination unit 9 and the failure time are displayed. The principal component score “Z m ” is subjected to wavelet analysis, and is also formed into a contour line pattern to generate a normal contour line pattern and an abnormal contour line pattern, which are stored in the analysis data storage means 10. When the wavelet analysis unit 14 and the contour line patterning unit 15 finish the wavelet analysis process and the contour line patterning process for all the principal component scores “Z m ” output from the principal component analysis unit 13, the land area calculation unit 16, the height of the amplitude corresponding to the sea level is obtained and stored based on the normal contour line pattern stored in the analysis data storage means 10 and the abnormal contour line pattern.

[各対象建物側の空調装置などの監視制御を行っているときの動作]
この後、プロセス管理装置1により、各対象建物側の空調装置などの監視制御が開始されると、プロセス実績データ蓄積手段7によって、デパートビル2、事務所ビル3、病院ビル4等の対象建物側に設置された各計測監視制御手段5から出力される各プロセス実績データ、例えば各空調プロセスのエネルギー消費量、その他のプロセス値(外気温、室温、室内湿度等)などを含むプロセス実績データが取り込まれ順次、記憶される。
[Operation when monitoring control of air conditioners etc. of each target building]
Thereafter, when the process control device 1 starts monitoring and controlling the air conditioners and the like of each target building, the process result data storage means 7 sets the target buildings such as the department store building 2, the office building 3, the hospital building 4, and the like. Each process result data output from each measurement monitoring control means 5 installed on the side, for example, process result data including energy consumption of each air conditioning process, other process values (outside temperature, room temperature, indoor humidity, etc.), etc. Captured and stored sequentially.

また、この動作と並行し、統計解析手段8の主成分分析手段13によって、各対象建物側の空調装置などの監視制御を開始する前に得られたプロセス変数結合重みパラメータ“a11、a12、…、amp”が用いられて、プロセス実績データ蓄積手段7に蓄積されているプロセス実績データに対し、前記(10)式〜(12)式に示す演算が行われて、主成分スコア“Z”が求められ、プロセス状態判定手段9に供給される。 Further, in parallel with this operation, the process variable combination weight parameters “a 11 , a 12 obtained by the principal component analysis unit 13 of the statistical analysis unit 8 before starting the monitoring control of the air conditioners and the like on each target building side are started. ,..., A mp ″ are used, and the operations shown in the above equations (10) to (12) are performed on the process result data stored in the process result data storage means 7 to obtain the principal component score “ Z m ″ is obtained and supplied to the process state determination means 9.

そして、プロセス状態判定手段9のウェーブレット解析手段14によって、前記(2)式に示す(連続)ウェーブレット変換式が用いられて、主成分スコア“Z”がウェーブレット解析された後、等高線パターン化手段15によって、ウェーブレット解析結果が等高線パターン化され、解析データ蓄積手段10に蓄積されるとともに、陸地面積演算手段16によって、各対象建物側の空調装置などの監視制御が開始される前に得られた海面の振幅を示す高さが用いられて、等高線パターン化手段15で得られた等高線パターンがレベル弁別され、陸地部分の面積(陸地面積)が求められる。 Then, the wavelet analysis unit 14 of the process state determination unit 9 uses the (continuous) wavelet transform equation shown in the above equation (2) to perform wavelet analysis of the principal component score “Z m ”, and then performs contour line patterning unit 15, the wavelet analysis result is contoured and accumulated in the analysis data accumulating unit 10, and obtained by the land area calculating unit 16 before the monitoring control of the air conditioner etc. of each target building is started. The height indicating the amplitude of the sea surface is used, and the contour line pattern obtained by the contour line patterning means 15 is subjected to level discrimination to obtain the land area (land area).

また、この動作と並行し、プロセス状態判定手段9の正常/異常判定手段17によって、解析データ蓄積手段10に蓄積されている等高線パターンが読み出されて、表示装置上に周波数軸、時間軸を持つ3次元の等高線パターンが表示されるとともに、陸地面積演算手段16から出力される陸地面積と、予め設定されているしきい値とが比較され、陸地面積演算手段16から出力される陸地面積が予め設定されているしきい値より小さいとき、プロセスが正常であると判定されて、判定結果が解析データ蓄積手段10などに供給され、等高線パターン化手段15から出力される今回の等高線パターンが正常時の等高線パターンとして、蓄積されるとともに、表示装置にプロセスが正常であることを示すメッセージなどが表示される。   In parallel with this operation, the normal / abnormal determination means 17 of the process state determination means 9 reads the contour line pattern stored in the analysis data storage means 10 and sets the frequency axis and time axis on the display device. A three-dimensional contour line pattern is displayed, and the land area output from the land area calculating means 16 is compared with a preset threshold value, and the land area output from the land area calculating means 16 is calculated. When the threshold value is smaller than a preset threshold value, it is determined that the process is normal, the determination result is supplied to the analysis data storage means 10 and the current contour line pattern output from the contour line patterning means 15 is normal. A time contour pattern is accumulated and a message indicating that the process is normal is displayed on the display device.

また、陸地面積演算手段16から出力される陸地面積が予め設定されているしきい値を超えているとき、正常/異常判定手段17によって、プロセスが異常であると判定されて、判定結果が解析データ蓄積手段10などに供給され、等高線パターン化手段15から出力される今回の等高線パターンが異常時の等高線パターンとして、蓄積されるとともに、表示装置にプロセスが異常であることを示すメッセージなどが表示される。   When the land area output from the land area calculation means 16 exceeds a preset threshold value, the normal / abnormality determination means 17 determines that the process is abnormal, and the determination result is analyzed. The current contour line pattern supplied to the data storage means 10 and the like and output from the contour line patterning means 15 is stored as a contour line pattern at the time of abnormality, and a message indicating that the process is abnormal is displayed on the display device. Is done.

そして、プロセス管理装置1の運転期間が長くなり、プロセス実績データ蓄積手段7に正常時のプロセス実績データ、故障時のプロセス実績データが十分に蓄積されると、主成分分析手段13によって、プロセス実績データ蓄積手段7に蓄積されている正常時のプロセス実績データ、故障時のプロセス実績データが分析されて、新たなプロセス変数結合重みパラメータ“a11、a12、…、amp”が求めら、プロセス変数結合重みパラメータ“a11、a12、…、amp”が更新される。 When the operation period of the process management apparatus 1 becomes long and the process result data storage unit 7 has enough process result data at normal time and process result data at failure, the process result data is stored by the principal component analysis unit 13. The normal process performance data and the failure process performance data stored in the data storage means 7 are analyzed, and new process variable coupling weight parameters “a 11 , a 12 ,..., A mp ” are obtained. The process variable connection weight parameter “a 11 , a 12 ,..., A mp ” is updated.

また、この動作と並行し、プロセス管理装置1の運転期間が長くなり、解析データ蓄積手段10に正常時の等高線パターンと、異常時の等高線パターンとが十分に蓄積されると、陸地面積演算手段16によって、解析データ蓄積手段10に蓄積されている正常時の等高線パターンと、異常時の等高線パターンとに基づき、海面に相当する振幅の高さが再度、求められて、新たな高さとして記憶される。   In parallel with this operation, when the operation period of the process management apparatus 1 becomes long, and the contour line pattern at the normal time and the contour line pattern at the abnormal time are sufficiently stored in the analysis data storage means 10, the land area calculating means 16, the height of the amplitude corresponding to the sea level is obtained again based on the normal contour line pattern stored in the analysis data storage means 10 and the abnormal contour line pattern, and stored as a new height. Is done.

<実施形態のシミュレーション結果>
次に、図3〜図6に示す説明図を用いて、プロセス管理装置1の正常/異常判定能力を確認したときのシミュレーション内容を説明する。なお、ウェーブレット変換ソフトとして、便宜的に音声解析ソフトを用いているので、図3〜図6に示す説明図では、横軸が秒単位になっているが、実単位は“×100時間”である。
<Simulation result of embodiment>
Next, simulation contents when the normality / abnormality determination capability of the process management apparatus 1 is confirmed will be described using the explanatory diagrams shown in FIGS. Since audio analysis software is used as wavelet conversion software for convenience, in the explanatory diagrams shown in FIGS. 3 to 6, the horizontal axis is in seconds, but the actual unit is “× 100 hours”. is there.

まず、対象プロセスとして、デパートビル2に設けられた空調プロセスを用い、またプロセス実績データの変数として、通常計測可能な下記の5つの変数(P=5)、すなわち現在、稼働している実際の空調プロセスなどに、新たにセンサを追加することなく、プロセス実績データを収集することができる変数を用いた。   First, the air conditioning process provided in the department store building 2 is used as the target process, and the following five variables (P = 5) that can be normally measured as the variables of the process performance data, that is, the actual operation currently in progress Variables that can collect process performance data without adding new sensors are used for air conditioning processes.

室温(計測ノイズ有) …(13)
冷水操作弁MV …(14)
室内湿度 …(15)
外気温度 …(16)
外気湿度 …(17)
そして、空調プロセスの運転時間が“8時30分”から“20時(計算周期:1分)”までという条件で、6日間、外気温度や室内熱負荷などの空調負荷条件を変えながら、1分単位で収集される冷房期間中のプロセス実績データのうち、6日分に対応する4,140個(6日×11.5時間×60)のプロセス実績データを用いて、プロセス変数結合重みパラメータ“a11、a12、…、amp”を求めた。
Room temperature (with measurement noise) (13)
Cold water operation valve MV (14)
Indoor humidity (15)
Outside temperature (16)
Outside air humidity (17)
Then, under the condition that the operation time of the air conditioning process is from “8:30” to “20:00 (calculation cycle: 1 minute)”, while changing the air conditioning load conditions such as the outside air temperature and the indoor heat load for 6 days, 1 Of the process performance data collected during the cooling period, the process variable combination weight parameter is used using 4,140 (6 days × 11.5 hours × 60) process performance data corresponding to 6 days. “A 11 , a 12 ,..., A mp ” was obtained.

次いで、正常時のプロセス実績データうち、1日分のプロセス実績データを選択し、上述したプロセス変数結合重みパラメータ“a11、a12、…、amp”を用いて、主成分分析を行ったとき、図3のトレンドグラフに示すように、変化する主成分スコア“Z”が得られ、この主成分スコア“Z”をウェーブレット解析したとき、図4の説明図に示すように、周波数軸、時間軸を持つ3次元の等高線パターン18と、2次元の等高線パターン19とが得られた。 Next, among the process result data at the normal time, the process result data for one day was selected, and the principal component analysis was performed using the process variable coupling weight parameters “a 11 , a 12 ,..., A mp ” described above. When the principal component score “Z i ” that changes is obtained as shown in the trend graph of FIG. 3 and this principal component score “Z i ” is subjected to wavelet analysis, as shown in the explanatory diagram of FIG. A three-dimensional contour line pattern 18 having an axis and a time axis and a two-dimensional contour line pattern 19 were obtained.

また、給気ファンの能力が低下したときのプロセス実績データを含む、1日分のプロセス実績データを選択し、上述したプロセス変数結合重みパラメータ“a11、a12、…、amp”を用いて主成分分析を行い、この主成分分析で得られた主成分スコア“Z”をウェーブレット解析したとき、図5の説明図に示すように、周波数軸、時間軸を持つ3次元の等高線パターン20と、2次元の等高線パターン21とが得られた。 Further, the process performance data for one day including the process performance data when the capacity of the air supply fan is reduced is selected, and the above-described process variable combination weight parameters “a 11 , a 12 ,..., A mp ” are used. When the principal component analysis is performed and the principal component score “Z i ” obtained by the principal component analysis is wavelet analyzed, a three-dimensional contour pattern having a frequency axis and a time axis as shown in FIG. 20 and a two-dimensional contour line pattern 21 were obtained.

また、冷却コイルの能力が低下したときのプロセス実績データを含む、1日分のプロセス実績データを選択し、上述したプロセス変数結合重みパラメータ“a11、a12、…、amp”を用いて主成分分析を行い、この主成分分析で得られた主成分スコア“Z”をウェーブレット解析したとき、図6の説明図に示すように、周波数軸、時間軸を持つ3次元の等高線パターン22と、2次元の等高線パターン23とが得られた。 Further, process performance data for one day including process performance data when the capacity of the cooling coil is reduced is selected, and the above-described process variable coupling weight parameters “a 11 , a 12 ,..., A mp ” are used. When the principal component analysis is performed and the principal component score “Z i ” obtained by the principal component analysis is wavelet analyzed, as shown in the explanatory diagram of FIG. 6, a three-dimensional contour pattern 22 having a frequency axis and a time axis is obtained. As a result, a two-dimensional contour line pattern 23 was obtained.

これらのシミュレーション結果から分かるように、図4に示す等高線パターン18の陸地面積、図5に示す等高線パターン20の陸地面積、図6に示す等高線パターン22の陸地面積が各々、異なることから、図4に示す等高線パターン18の陸地面積に対応する値をしきい値にすることにより、図4に示す正常時のプロセス実績データに対する3次元の等高線パターン18と異なる等高線パターン、例えば図5に示す3次元の等高線パターン20が得られたとき、あるいは図6に示す3次元の等高線パターン22が得られたとき、空調プロセスに異常が発生したと判断することができる。   As can be seen from these simulation results, the land area of the contour line pattern 18 shown in FIG. 4, the land area of the contour line pattern 20 shown in FIG. 5, and the land area of the contour line pattern 22 shown in FIG. By using a value corresponding to the land area of the contour line pattern 18 shown in FIG. 4 as a threshold value, a contour line pattern different from the three-dimensional contour line pattern 18 for the normal process performance data shown in FIG. 4, for example, the three-dimensional figure shown in FIG. When the contour line pattern 20 is obtained, or when the three-dimensional contour line pattern 22 shown in FIG. 6 is obtained, it can be determined that an abnormality has occurred in the air conditioning process.

また、表示装置に異常発生を示すメッセージとともに、図5に示す3次元の等高線パターン20が表示されたとき、波形形状から給気ファンの能力が低下したことを認識し、また図6に示す3次元の等高線パターン22が表示されたとき、波形形状から冷却コイルの能力が低下したことを認識することができる。   Further, when the three-dimensional contour line pattern 20 shown in FIG. 5 is displayed together with a message indicating the occurrence of an abnormality on the display device, it is recognized from the waveform shape that the ability of the air supply fan has decreased, and 3 shown in FIG. When the dimensional contour line pattern 22 is displayed, it can be recognized from the waveform shape that the capacity of the cooling coil has decreased.

<実施形態の効果>
このように、第1の実施形態では、各対象建物側の空調装置などの監視制御を開始したとき、プロセス変数結合重みパラメータ“a11、a12、…、amp”を使用して、プロセス実績データ蓄積手段7に蓄積されているプロセス実績データを主成分分析するとともに、この主成分分析で得られた主成分スコア“Z”をウェーブレット解析処理、等高線パターン化処理して、陸地部分の面積を求め、これが予め設定されているしきい値を超えたとき、プロセス側が異常になっていると判定するようにしているので、プロセス側から出力されるプロセス実績データを主成分スコアに変換して、ウェーブレット解析対象となるデータ種類数を抑制し、ウェーブレット解析処理、異常有無判定にかかる時間を短くしながら、ウェーブレット解析で得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定し、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減し、保守費を低減することができる。
<Effect of embodiment>
As described above, in the first embodiment, when the monitoring control of the air conditioner or the like on each target building side is started, the process variable combination weight parameters “a 11 , a 12 ,. The process result data stored in the result data storage means 7 is subjected to principal component analysis, and the principal component score “Z i ” obtained by the principal component analysis is subjected to wavelet analysis processing, contour line pattern processing, and When the area is obtained and this exceeds a preset threshold value, the process side is judged to be abnormal, so the process performance data output from the process side is converted to the principal component score. This reduces the number of data types subject to wavelet analysis, shortens the time required for wavelet analysis processing and abnormality determination, and reduces wavelet analysis. It is possible to efficiently discriminate and estimate whether the process is normal from the change in spectrum obtained by analysis, and to perform system performance verification, failure detection, failure prediction, and failure diagnosis with high accuracy, while reducing energy consumption and maintenance costs. Can be reduced.

また、第1の実施形態では、各対象建物側の空調装置などの監視制御を開始する前に、統計解析手段8のシミュレーション処理で得られた正常時のプロセス実績データ、異常時のプロセス実績データに基づき、プロセス変数結合重みパラメータ“a11、a12、…、amp”とともに、正常時の等高線パターン、異常時の等高線パターン、海面の振幅を示す高さとを求めるようにしているので、プロセス側から出力されるプロセス実績データの種類が多いときでも、プロセスの稼働を開始した時点から、ウェーブレット解析対象となるデータ種類数を抑制して、ウェーブレット解析処理、異常有無判定にかかる時間を短くしながら、ウェーブレット解析で得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定し、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減し、保守費を低減することができる。 Moreover, in 1st Embodiment, before starting monitoring control of the air conditioner etc. of each object building side, the process performance data at the time of normal obtained by the simulation process of the statistical analysis means 8, and the process performance data at the time of abnormality , And a process variable coupling weight parameter “a 11 , a 12 ,..., A mp ”, a normal contour line pattern, an abnormal contour line pattern, and a height indicating the amplitude of the sea surface are obtained. Even when there are many types of process result data output from the side, the number of data types subject to wavelet analysis can be reduced from the time the process starts to shorten the time required for wavelet analysis processing and abnormality determination. However, it is possible to efficiently discriminate and estimate whether the process is normal from the change in spectrum obtained by wavelet analysis, Performance verification and fault detection of the stem, failure prediction, along with a failure diagnosis can be performed with high accuracy, it is possible to reduce energy consumption, reduce maintenance costs.

なお、上述した際1の実施形態では、統計解析手段8によって、プロセス実績データ蓄積手段7に蓄積されているプロセス実績データを主成分分析して、データの種類数を低減しながら、プロセス状態判定手段9に供給し、ウェーブレット解析処理、等高線パターン化処理、陸地面積算出処理、正常/異常判定処理などを行うようにしているが、プロセス実績データのデータ種類数が少ないときには、統計解析手段8を削除し、プロセス実績データ蓄積手段7に蓄積されているプロセス実績データを直接、プロセス状態判定手段9に供給して、ウェーブレット解析処理、等高線パターン化処理、陸地面積算出処理、正常/異常判定処理などを行うようにしても良い。   In the first embodiment described above, the statistical analysis unit 8 performs principal component analysis on the process result data stored in the process result data storage unit 7 to reduce the number of types of data while determining the process state. The data is supplied to the means 9, and wavelet analysis processing, contour line patterning processing, land area calculation processing, normality / abnormality determination processing, and the like are performed. The process result data stored in the process result data storage unit 7 is directly supplied to the process state determination unit 9, and wavelet analysis processing, contour line patterning processing, land area calculation processing, normal / abnormality determination processing, etc. May be performed.

これにより、プロセス管理装置1全体の構成を簡素化し、処理時間を短縮しながら、プロセス側から出力されるプロセス実績データをウェーブレット解析して得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定し、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減し、保守費を低減することができる。   This simplifies the configuration of the entire process management device 1 and shortens the processing time, while efficiently determining whether the process is normal from the change in spectrum obtained by wavelet analysis of the process performance data output from the process side. It is possible to estimate, perform system performance verification, failure detection, failure prediction, and failure diagnosis with high accuracy, reduce energy consumption, and reduce maintenance costs.

《第2の実施形態》
図7は本発明によるプロセス管理装置の第2の実施形態の全体構成例を示すブロック図である。なお、この図において、図1の各部と同じ部分には、同じ符号が付してある。
<< Second Embodiment >>
FIG. 7 is a block diagram showing an example of the overall configuration of the second embodiment of the process management apparatus according to the present invention. In this figure, the same parts as those in FIG. 1 are denoted by the same reference numerals.

この図に示すプロセス管理装置31が図1に示すプロセス管理装置1と異なる点は、プロセス状態判定手段9から出力される等高線パターンを正常時の等高線パターン、異常時の等高線パターンに区別して蓄積する解析データ蓄積手段10に代えて、プロセス状態判定手段9から出力される等高線パターンを正常時の等高線パターン、故障内容別の等高線パターンに区別して蓄積するとともに、故障要因判定手段33から出力される故障内容別の陸地パターンを蓄積する故障知識データ蓄積手段32を設けるとともに、プロセス状態判定手段9でプラントが異常になっていると判定されたとき、プロセス状態判定手段9から出力される等高線パターンと故障知識データ蓄積手段32に蓄積されている故障内容別の陸地パターンとに基づき、どのシステム系統が異常になっているかを判定する故障要因判定手段33を設けたことである。   The process management apparatus 31 shown in this figure is different from the process management apparatus 1 shown in FIG. 1 in that the contour line patterns output from the process state determining means 9 are distinguished and stored as normal contour lines and abnormal contour lines. Instead of the analysis data accumulating means 10, the contour line pattern output from the process state determination means 9 is stored separately in a normal contour line pattern and a contour line pattern for each failure content, and the failure output from the failure factor determination means 33. The fault knowledge data storage means 32 for storing the land patterns according to the contents is provided, and when the process state determination means 9 determines that the plant is abnormal, the contour line pattern output from the process state determination means 9 and the failure Based on the land pattern for each failure content stored in the knowledge data storage means 32, System system is by providing the failure factor determination means 33 determines have become abnormal.

故障知識データ蓄積手段32は、プロセス状態判定手段9の正常/異常判定手段17から出力される正常/異常判定結果、故障要因判定手段33から出力される故障内容などに基づき、プロセス状態判定手段9の等高線パターン化手段15から出力される新たな等高線パターンを正常時の等高線パターンと、故障内容別の等高線パターンとに分類し、1日分単位の時系列データとして蓄積しながら、プロセス状態判定手段9の陸地面積演算手段16から読み出し要求が出されたとき、蓄積している正常時の等高線パターンと、故障内容別の等高線パターンとを読み出して、陸地面積演算手段16に供給する。また、故障要因判定手段33から新たな故障内容別の陸地パターン、例えば図8に示すように、A系統故障時に得られた等高線パターンをレベル弁別し、陸地に対応する升目を“1”、海に対する升目を“0”にしたA系統故障時の陸地パターン34、図9に示すように、B系統故障時に得られた等高線パターンをレベル弁別し、陸地に対応する升目を“1”、海に対する升目を“0”にしたB系統故障時の陸地パターン35、図10に示すように、C系統故障時に得られた等高線パターンをレベル弁別し、陸地に対応する升目を“1”、海に対する升目を“0”にしたC系統故障時の陸地パターン36などが出力される毎に、これを取り込んで蓄積し、故障要因判定手段33から読み出し要求が出されたとき、蓄積している故障内容別の陸地パターンを読み出して、故障要因判定手段33に供給する。   The failure knowledge data storage unit 32 is based on the normal / abnormal determination result output from the normal / abnormal determination unit 17 of the process state determination unit 9, the failure content output from the failure factor determination unit 33, and the like. The new contour line pattern output from the contour line patterning means 15 is classified into a normal contour line pattern and a contour line pattern for each failure content, and accumulated as time series data in units of one day, while the process state determining means When a read request is issued from the land area calculation means 16 of 9, the stored normal contour line patterns and the contour line patterns classified by failure contents are read and supplied to the land area calculation means 16. Further, the land pattern for each new failure content from the failure factor determination means 33, for example, the contour line pattern obtained at the time of system A failure is level-discriminated as shown in FIG. As shown in FIG. 9, the land pattern 34 at the time of system A failure with the grid for “0” is level discriminated as shown in FIG. 9, and the grid corresponding to the land is “1” for the sea The land pattern 35 at the time of failure of the B system with the grid set to “0”, and the contour line pattern obtained at the time of failure of the C system as shown in FIG. Each time the land pattern 36 at the time of C system failure is output with the value of “0” being stored, it is captured and stored, and when a read request is issued from the failure factor determination means 33, Land of It reads the pattern, and supplies the failure factor determination means 33.

故障要因判定手段33は、各対象建物側の空調装置などの監視制御を開始する前に、統計解析手段8のシミュレーション動作で、各故障内容別のプロセス実績データに対応する主成分スコア“Z”が生成され、プロセス状態判定手段9でウェーブレット解析、等高線パターン化されたときに得られる故障内容別の等高線パターン、海面の振幅を示す高さを取り込むとともに、海面の振幅を示す高さで故障内容別の等高線パターンをレベル弁別して、図8〜図10に示すような故障内容別の陸地パターン34〜36を生成し、故障知識データ蓄積手段32に蓄積する。 The failure factor determination means 33 performs the simulation operation of the statistical analysis means 8 before starting the monitoring control of the air conditioner etc. on each target building side, and the principal component score “Z i corresponding to the process result data for each failure content”. ”Is generated, the contour pattern for each failure content obtained when the process state determination means 9 performs wavelet analysis and contour patterning, the height indicating the amplitude of the sea surface is taken in, and the failure is detected at the height indicating the amplitude of the sea surface. The contour patterns according to contents are subjected to level discrimination to generate land patterns 34 to 36 according to fault contents as shown in FIGS. 8 to 10 and stored in the fault knowledge data storage means 32.

そして、各対象建物側の空調装置などの監視制御が行われ、プロセス状態判定手段9の正常/異常判定手段17でプロセス異常と判定されたとき、図11、図12に示すように、等高線パターン化手段15で得られた等高線パターン37の陸地部分に対応する升目に“1”を割り付け、海に対する升目に“0”を割り付けた陸地パターン38を作成した後、この陸地パターン38の各升目と、故障知識データ蓄積手段32に蓄積されている故障内容別の陸地パターン34〜36の各升目との積和を各々、求めた後、積和値が最大となる故障内容別の陸地パターンを選択して、この陸地パターンの故障内容と同じ故障が発生したと推定し、表示装置上に故障内容を表示するとともに、今回の故障内容を故障知識データ蓄積手段32に供給して、等高線パターン化手段15から出力された等高線パターン37を故障内容別の等高線パターンとして、蓄積する。   Then, when the monitoring control of the air conditioners and the like on each target building side is performed and the normal / abnormal determination unit 17 of the process state determination unit 9 determines that the process is abnormal, as shown in FIG. 11 and FIG. After assigning “1” to the grid corresponding to the land portion of the contour line pattern 37 obtained by the converting means 15 and creating the land pattern 38 assigned “0” to the sea, the land pattern 38 After obtaining the sum of products of each of the land patterns 34 to 36 for each failure content stored in the failure knowledge data storage means 32 and selecting the land pattern for each failure content having the maximum product sum value. Then, it is estimated that the same failure as the failure content of this land pattern has occurred, the failure content is displayed on the display device, and the current failure content is supplied to the failure knowledge data storage means 32, Contours pattern 37 as a fault condition by the contour pattern output from the high-line patterning means 15, accumulates.

この後、プロセス管理装置31の運転期間が長くなり、故障知識データ蓄積手段32に故障内容別の等高線パターンが十分に蓄積されたとき、陸地面積演算手段16で使用されている海の高さを取り込ませるとともに、故障知識データ蓄積手段32に蓄積されている故障内容別の等高線パターンを海の高さでレベル弁別して、陸地に対応する升目に“1”、海に対応する升目に“0”を割り当てた陸地パターンを作成し、故障知識データ蓄積手段32に蓄積されている故障内容別の陸地パターンを更新する。   Thereafter, when the operation period of the process management device 31 is extended and the contour line patterns for each failure content are sufficiently accumulated in the failure knowledge data accumulation means 32, the height of the sea used by the land area calculation means 16 is calculated. At the same time, the contour line pattern for each failure content stored in the failure knowledge data storage means 32 is discriminated by the level of the sea, and the grid corresponding to the land is “1”, and the grid corresponding to the sea is “0”. Is created, and the land pattern for each failure content stored in the failure knowledge data storage means 32 is updated.

このように、第2の実施形態では、故障内容別の等高線パターンに対応した故障内容別の陸地パターンを作成して、故障知識データ蓄積手段32に蓄積し、プロセス状態判定手段9でプロセスが異常であると判定されたとき、故障要因判定手段33によって、プロセス状態判定手段9から出力される等高線パターンを取り込ませて、陸地パターンを作成するとともに、この陸地パターンと故障知識データ蓄積手段32に蓄積されている故障内容別の陸地パターンとの類似性を判定し、この判定結果に基づき、故障内容を推定するようにしているので、プロセス側から出力されるプロセス実績データをウェーブレット解析して得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定し、故障発生時にはその要因を推定して、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減し、保守費を低減することができる。   As described above, in the second embodiment, a land pattern for each failure content corresponding to the contour line pattern for each failure content is created and stored in the failure knowledge data storage means 32, and the process state determination means 9 makes the process abnormal. Is determined, the failure factor determination means 33 takes in the contour line pattern output from the process state determination means 9 to create a land pattern and stores it in the land pattern and failure knowledge data storage means 32. Since the similarity with the land pattern for each failure content is determined and the failure content is estimated based on this determination result, it can be obtained by wavelet analysis of process performance data output from the process side Efficiently discriminate and estimate whether the process is normal from the change in the spectrum, and estimate the cause when a failure occurs, Performance verification and fault detection of the stem, failure prediction, along with a failure diagnosis can be performed with high accuracy, it is possible to reduce energy consumption, reduce maintenance costs.

また、第2の実施形態では、故障時の等高線パターンを格子状に分割し、陸地に対応する升目部分に“1”を割り当て、海に対応する升目に“0”を割り当てた故障内容別の陸地パターンを作成して、故障知識データ蓄積手段32に蓄積しておき、プロセス状態判定手段9でプロセスの異常と判定されたとき、故障要因判定手段33によって、プロセス状態判定手段9から出力される等高線パターンを格子状に分割し、陸地に対応する升目部分に“1”を割り当て、海に対応する升目に“0”を割り当てした陸地パターンを作成し、この陸地パターンの各升目と、故障知識データ蓄積手段32に蓄積されている故障内容別の陸地パターンの各升目との積和を取り、積和値が最大となる陸地パターンの故障内容に基づき、故障内容を推定するようにしているので、故障内容判定データの作成、故障内容判定処理を簡素化して、プロセス管理装置31の製造を容易にしながら、プロセス側から出力されるプロセス実績データをウェーブレット解析して得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定し、故障発生時には、その要因を素早く推定して、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減し、保守費を低減することができる。   Further, in the second embodiment, the contour line pattern at the time of failure is divided into a grid, and “1” is assigned to the grid portion corresponding to the land, and “0” is assigned to the grid corresponding to the ocean. A land pattern is created and stored in the failure knowledge data storage unit 32. When the process state determination unit 9 determines that the process is abnormal, the failure state determination unit 33 outputs the process pattern determination unit 9. Divide the contour line pattern into a grid, assign "1" to the grid corresponding to the land, and create a land pattern that assigns "0" to the sea corresponding to the sea. The product sum of each land pattern for each failure content stored in the data storage means 32 is summed, and the failure content is estimated based on the failure content of the land pattern with the maximum product sum value. Therefore, it is possible to simplify the creation of the failure content determination data and the failure content determination processing, and facilitate the manufacture of the process management device 31, while the process result data output from the process side is subjected to wavelet analysis. Efficiently discriminates and estimates whether the process is normal from changes, quickly estimates the cause of a failure, and enables accurate system performance verification, failure detection, failure prediction, and failure diagnosis, as well as energy consumption And maintenance costs can be reduced.

また、第2の実施形態では、各対象建物側の空調装置などの監視制御を開始する前に、統計解析手段8のシミュレーション動作で各故障内容別の主成分スコア“Z”を生成するとともに、プロセス状態判定手段9によって、主成分スコア“Z”をウェーブレット解析処理、等高線パターン化処理して、故障内容別の等高線パターン、海面の振幅を示す高さを作成した後、海面の振幅を示す高さで故障内容別の等高線パターンをレベル弁別して、故障内容別の陸地パターンを作成し、故障知識データ蓄積手段32に蓄積するようにしているので、シミュレーション処理で得られる異常時のプロセス実績データに基づき、故障内容推定処理で使用する故障内容別の陸地パターンを算出して、プロセスの稼働を開始した時点から、プロセス側から出力されるプロセス実績データをウェーブレット解析して得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定し、故障発生時には、その要因を推定して、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減し、保守費を低減することができる。 In the second embodiment, the principal component score “Z i ” for each failure content is generated by the simulation operation of the statistical analysis means 8 before starting the monitoring control of the air conditioner etc. of each target building. Then, the component state score “Z i ” is subjected to wavelet analysis processing and contour line patterning processing by the process state determination means 9 to create a contour line pattern for each failure content and a height indicating the amplitude of the sea surface. The contour pattern for each failure content is discriminated at the indicated height, and a land pattern for each failure content is created and stored in the failure knowledge data storage means 32. Based on the data, the land pattern for each failure content used in the failure content estimation process is calculated, and the process side Efficiently discriminate and estimate whether the process is normal from the changes in the spectrum obtained by wavelet analysis of the process performance data output from the system, and in the event of a failure, the cause is estimated to verify system performance, failure detection, and failure Prediction and failure diagnosis can be performed with high accuracy, energy consumption can be reduced, and maintenance costs can be reduced.

また、第2の実施形態では、プロセス管理装置31の運転期間が長くなり、故障知識データ蓄積手段32に故障内容別の等高線パターンが十分に蓄積されたとき、陸地面積演算手段16で使用されている海の高さを取り込ませるとともに、故障知識データ蓄積手段32に蓄積されている故障内容別の等高線パターンを海の高さでレベル弁別して、陸地に対応する升目に“1”、海に対応する升目に“0”を割り当てた陸地パターンを作成し、故障知識データ蓄積手段32に蓄積されている故障内容別の陸地パターンを更新するようにしているので、プロセス側から出力される異常時のプロセス実績データに基づき、故障内容推定処理で使用する故障内容別の陸地パターンを更新し、プラント側の状態変化などに応じて、故障内容別の陸地パターンを最適化しながら、プロセス側から出力されるプロセス実績データをウェーブレット解析して得られるスペクトルの変化からプロセスが正常かどうかを効率良く判別推定し、故障発生時には、その要因を推定して、システムの性能検証や故障検知、故障予知、故障診断を精度良く行うことができるとともに、エネルギー消費を削減し、保守費を低減することができる。   Further, in the second embodiment, when the operation period of the process management device 31 becomes long and the contour line pattern for each failure content is sufficiently accumulated in the failure knowledge data accumulation unit 32, it is used by the land area calculation unit 16. The level of the contour line for each failure content stored in the failure knowledge data storage means 32 is discriminated based on the level of the ocean, and the grid corresponding to the land is set to “1”, corresponding to the ocean. Since the land pattern assigned with “0” is created and the land pattern according to the failure content stored in the failure knowledge data storage means 32 is updated, the error pattern output from the process side is displayed. Based on the process performance data, the land pattern for each failure content used in the failure content estimation process is updated, and the land pattern for each failure content is changed according to the state change on the plant side. The process performance data output from the process side is optimized while efficiently determining whether or not the process is normal from the changes in the spectrum obtained by wavelet analysis. Performance verification, failure detection, failure prediction, and failure diagnosis can be performed with high accuracy, energy consumption can be reduced, and maintenance costs can be reduced.

《他の実施形態》
なお、本発明は、上記実施の形態に限定されるものではなく、その要旨を逸脱しない範囲で種々変形して実施できるものである。
<< Other embodiments >>
In addition, this invention is not limited to the said embodiment, It can implement variously in the range which does not deviate from the summary.

また、上述した実施形態では、具体例として各建物の空調プロセスを制御するBEMS装置の故障検知、故障予知、故障診断機能を用いているが、本発明が一般的なプロセスに適用できることは明らかである。   Further, in the above-described embodiment, the failure detection, failure prediction, and failure diagnosis functions of the BEMS device that controls the air conditioning process of each building are used as specific examples. However, it is obvious that the present invention can be applied to general processes. is there.

また、正常/異常の判定アルゴリズムとして、予め設定されたしきい値を逸脱したか否かにより判定するアルゴリズムを使用するようにしているが、瞬時値を用いて、正常/異常を判定するアルゴリズムを使用するようにしても良い。ただし、このような判定アルゴリズムを使用した場合“Noise Factor”によるデータ変動の影響で、誤判断の可能性があることから、ある時間幅のデータに対してしきい値を逸脱したデータの割合がどの程度かを評価し、正常・異常を判断するアルゴリズムと組み合わせて、ノイズに起因する誤判定のリスクを低減するようにしても良い。   In addition, as a normal / abnormal determination algorithm, an algorithm that is determined based on whether or not a preset threshold value is deviated is used. An algorithm for determining normal / abnormal using an instantaneous value is used. It may be used. However, if such a judgment algorithm is used, there is a possibility of misjudgment due to the influence of data fluctuation due to “Noise Factor”. The degree of misjudgment caused by noise may be reduced in combination with an algorithm for evaluating the degree of normality and abnormality.

また、上述した実施形態では、複数の建物に設置されている各空調設備などをリモートで監視、制御するBEMS装置を例にして、本発明を説明しているが、もちろん同様のアルゴリズムを建物単独で実施しても良い。この場合はリモートでなく、その建物に本発明によるプロセス管理装置が設置される。   In the above-described embodiment, the present invention has been described by taking as an example a BEMS device that remotely monitors and controls each air-conditioning facility installed in a plurality of buildings. You may carry out with. In this case, the process management apparatus according to the present invention is installed in the building, not remote.

本発明によるプロセス管理装置の第1の実施形態の全体構成例を示すブロック図である。It is a block diagram which shows the example of whole structure of 1st Embodiment of the process management apparatus by this invention. 図1に示すプロセスシミュレーション手段のシミュレーション動作例を示すフローチャートである。It is a flowchart which shows the simulation operation example of the process simulation means shown in FIG. 図1に示す統計解析手段をシミュレーション動作したときに得られる主成分スコアの一例を示すトレンドグラフである。It is a trend graph which shows an example of the principal component score obtained when the statistical analysis means shown in FIG. 図1に示すプロセス状態判定手段に正常時の主成分スコアを供給して、ウェーブレット解析、等高線パターン作成したときに得られた周波数軸、時間軸を持つ3次元の等高線パターンと、2次元の等高線パターンを示す説明図である。A normal component score is supplied to the process state determination means shown in FIG. 1 to obtain a wavelet analysis, a contour line pattern, a three-dimensional contour pattern having a frequency axis and a time axis, and a two-dimensional contour line. It is explanatory drawing which shows a pattern. 図1に示すプロセス状態判定手段に給気ファンの能力が低下したときの主成分スコアを供給して、ウェーブレット解析、等高線パターン作成したときに得られた周波数軸、時間軸を持つ3次元の等高線パターンと、2次元の等高線パターンを示す説明図である。The principal component score when the capacity of the air supply fan decreases is supplied to the process state judging means shown in FIG. 1, and the three-dimensional contour line having the frequency axis and the time axis obtained when the wavelet analysis and the contour line pattern are created. It is explanatory drawing which shows a pattern and a two-dimensional contour line pattern. 図1に示すプロセス状態判定手段に冷却コイルの能力が低下したときの主成分スコアを供給して、ウェーブレット解析、等高線パターン作成したときに得られた周波数軸、時間軸を持つ3次元の等高線パターンと、2次元の等高線パターンを示す説明図である。A three-dimensional contour pattern having a frequency axis and a time axis obtained by supplying a principal component score when the capacity of the cooling coil is reduced to the process state judging means shown in FIG. It is explanatory drawing which shows a two-dimensional contour line pattern. 本発明によるプロセス管理装置の第2の実施形態の全体構成例を示すブロック図である。It is a block diagram which shows the example of whole structure of 2nd Embodiment of the process management apparatus by this invention. 図7に示す故障要因判定手段によって、故障要因判定手段から出力されるA系統故障時の陸地パターン例を示す説明図である。It is explanatory drawing which shows the land pattern example at the time of the A system failure output from a failure factor determination means by the failure factor determination means shown in FIG. 図7に示す故障要因判定手段によって、故障要因判定手段から出力されるB系統故障時の陸地パターン例を示す説明図である。It is explanatory drawing which shows the land pattern example at the time of the B system failure output from a failure factor determination means by the failure factor determination means shown in FIG. 図7に示す故障要因判定手段によって、故障要因判定手段から出力されるC系統故障時の陸地パターン例を示す説明図である。It is explanatory drawing which shows the land pattern example at the time of C system failure output from a failure factor determination means by the failure factor determination means shown in FIG. 図7に示すプロセス状態判定手段から出力される等高線パターン例を示す説明図である。It is explanatory drawing which shows the example of the contour line pattern output from the process state determination means shown in FIG. 図7に示す故障要因判定手段で作成される陸地パターン例を示す説明図である。It is explanatory drawing which shows the land pattern example produced by the failure factor determination means shown in FIG.

符号の説明Explanation of symbols

1,31:プロセス管理装置
2:デパートビル
3:事務所ビル
4:病院ビル
5:計測監視制御手段
6:ネットワーク
7:プロセス実績データ蓄積手段
8:統計解析手段
9:プロセス状態判定手段
10:解析データ蓄積手段
11:プロセスシミュレーション手段
12:シミュレーションデータ蓄積手段
13:主成分分析手段
14:ウェーブレット解析手段
15:等高線パターン化手段
16:陸地面積演算手段
17:正常/異常判定手段
18,20,22:3次元の等高線パターン
19,21,23:2次元の等高線パターン
32:故障知識データ蓄積手段
33:故障要因判定手段
34:A系統故障時の陸地パターン
35:B系統故障時の陸地パターン
36:C系統故障時の陸地パターン
37:等高線パターン
38:陸地パターン
DESCRIPTION OF SYMBOLS 1,31: Process management apparatus 2: Department store building 3: Office building 4: Hospital building 5: Measurement monitoring control means 6: Network 7: Process performance data storage means 8: Statistical analysis means 9: Process state determination means 10: Analysis Data storage means 11: Process simulation means 12: Simulation data storage means 13: Principal component analysis means 14: Wavelet analysis means 15: Contour line patterning means 16: Land area calculation means 17: Normal / abnormality determination means 18, 20, 22: Three-dimensional contour line pattern 19, 21, 23: Two-dimensional contour line pattern 32: Failure knowledge data storage means 33: Failure factor determination means 34: Land pattern when system A fails 35: Land pattern when system B fails 36: C Land pattern at the time of system failure 37: Contour pattern 38: Land pattern N

Claims (7)

プロセスを時系列的に計測して得られる複数のプロセス実績データを蓄積するプロセス実績データ蓄積手段と、
このプロセス実績データ蓄積手段に蓄積されている各プロセス実績データをウェーブレット解析するとともに、この解析結果を3次元の等高線パターン化して異常の有無を判定するプロセス状態判定手段と、
を備えたことを特徴とするプロセス管理装置。
Process result data storage means for storing a plurality of process result data obtained by measuring processes in time series,
Wavelet analysis of each process performance data stored in the process performance data storage means, and a process state determination means for determining the presence or absence of abnormality by converting this analysis result into a three-dimensional contour line pattern;
A process management apparatus comprising:
請求項1に記載のプロセス管理装置において、
予め指定されたプロセス変数結合重みパラメータを用いて、前記プロセス実績データ蓄積手段に蓄積されている各プロセス実績データを主成分分析するとともに、この主成分分析で得られた主成分スコアを前記プロセス状態判定手段に供給する統計解析手段、
を備えたことを特徴とするプロセス管理装置。
The process management apparatus according to claim 1, wherein
The process result data stored in the process result data storage means is subjected to principal component analysis using a process variable combination weight parameter designated in advance, and the principal component score obtained by the principal component analysis is used as the process state. Statistical analysis means to be supplied to the judgment means;
A process management apparatus comprising:
請求項2に記載のプロセス管理装置において、
前記統計解析手段は、プロセスの正常状態と各種故障状態のシミュレーションを行うプロセスシミュレーション手段を備え、
シミュレーション処理で得られたプロセス実績データに基づき、主成分分析で使用するプロセス変数結合重みパラメータを算出した後、このプロセス変数結合重みパラメータを用いて、前記プロセス実績データ蓄積手段に蓄積されている各プロセス実績データを主成分分析するとともに、この主成分分析で得られた主成分スコアを前記プロセス状態判定手段に供給する、
ことを特徴とするプロセス管理装置。
The process management apparatus according to claim 2, wherein
The statistical analysis means includes a process simulation means for simulating the normal state of the process and various failure states,
After calculating the process variable coupling weight parameter used in the principal component analysis based on the process performance data obtained by the simulation process, each process variable data accumulated means is stored in the process performance data accumulation means using the process variable coupling weight parameter. The process result data is subjected to principal component analysis, and the principal component score obtained by the principal component analysis is supplied to the process state determination means.
A process management apparatus characterized by that.
請求項1、2、3のいずれかに記載のプロセス管理装置において、
故障内容別の等高線パターンに対応した故障内容別の陸地パターンを蓄積する故障知識データ蓄積手段と、
前記プロセス状態判定手段によって、プロセスが異常であると判定されたとき、前記プロセス状態判定手段から出力される等高線パターンから陸地パターンを作成するとともに、この陸地パターンと前記故障知識データ蓄積手段に蓄積されている故障内容別の陸地パターンとの類似性を判定し、この判定結果に基づき、故障内容を推定する故障要因判定手段と、
を備えたことを特徴とするプロセス管理装置。
In the process management device according to any one of claims 1, 2, and 3,
Failure knowledge data storage means for storing a land pattern by failure content corresponding to a contour line pattern by failure content;
When the process state determination means determines that the process is abnormal, a land pattern is created from the contour line pattern output from the process state determination means, and is stored in the land pattern and the failure knowledge data storage means. A failure factor determination means for determining similarity with a land pattern for each failure content and estimating the failure content based on the determination result;
A process management apparatus comprising:
請求項4に記載のプロセス管理装置において、
前記故障知識データ蓄積手段に蓄積される故障内容別の陸地パターン、前記故障要因判定手段によって生成される陸地パターンは、等高線パターンを格子状に分割し、陸地に対応する升目部分に“1”を割り当て、海に対応する升目に“0”を割り当てて作成され、
前記故障要因判定手段は、各陸地パターンの類似性を判定する際、前記プロセス状態判定手段から出力される等高線パターンに基づいて得られた陸地パターンの各升目と、前記故障知識データ蓄積手段に蓄積されている故障内容別の陸地パターンの各升目との積和を取り、積和値が最大となった陸地パターンの故障内容に基づき、故障内容を推定する、
ことを特徴とするプロセス管理装置。
The process management apparatus according to claim 4, wherein
The land pattern for each failure content stored in the failure knowledge data storage unit and the land pattern generated by the failure factor determination unit are obtained by dividing the contour line pattern into a grid pattern and setting “1” in the grid portion corresponding to the land. Assigned, created by assigning "0" to the cell corresponding to the sea,
When the failure factor determination means determines the similarity of each land pattern, each failure pattern of the land pattern obtained based on the contour line pattern output from the process state determination means is stored in the failure knowledge data storage means. Take the product sum with each square of the land pattern for each failure content, and estimate the failure content based on the failure content of the land pattern with the largest product sum value,
A process management apparatus characterized by that.
請求項4、5のいずれかに記載のプロセス管理装置において、
シミュレーションで得られた異常時のプロセス実績データを用いて故障内容別の等高線パターンを作成し、これら故障内容別の等高線パターンに基づき、故障内容別の陸地パターンを作成して、前記故障知識データ蓄積手段に蓄積する、
ことを特徴とするプロセス管理装置。
In the process management device according to any one of claims 4 and 5,
Using the process result data at the time of abnormality obtained by simulation, create contour lines for each failure content, create a land pattern for each failure content based on these contour patterns for each failure content, and store the failure knowledge data Accumulate in the means,
A process management apparatus characterized by that.
請求項4、5、6のいずれかに記載のプロセス管理装置において、
前記故障要因判定手段は、前記プロセス状態判定手段でプロセス異常と判定されたとき、故障内容を判定して、前記故障知識データ蓄積手段に蓄積されている等高線パターンのうち、プロセス異常と判定された等高線パターンに故障内容を付加し、前記故障知識データ蓄積手段に故障内容が付加された等高線パターンが十分に蓄積されたとき、この故障知識データ蓄積手段に蓄積されている故障内容が付加された等高線パターンから故障内容別の陸地パターンを作成して、前記故障知識データ蓄積手段に蓄積されている故障内容別の陸地パターンを更新する、
ことを特徴とするプロセス管理装置。
In the process management device according to any one of claims 4, 5, and 6,
The failure factor determination means determines the content of the failure when the process state determination means determines that the process is abnormal, and is determined to be process abnormality among the contour line patterns stored in the failure knowledge data storage means. Contour lines to which the fault contents stored in the fault knowledge data accumulating means are added when the fault contents are added to the contour pattern pattern and the contour line patterns with the fault contents are sufficiently accumulated in the fault knowledge data accumulating means. Create a land pattern for each failure content from the pattern, and update the land pattern for each failure content stored in the failure knowledge data storage means,
A process management apparatus characterized by that.
JP2004289052A 2004-09-30 2004-09-30 Process management device Expired - Fee Related JP4434905B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004289052A JP4434905B2 (en) 2004-09-30 2004-09-30 Process management device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004289052A JP4434905B2 (en) 2004-09-30 2004-09-30 Process management device

Publications (2)

Publication Number Publication Date
JP2006106870A true JP2006106870A (en) 2006-04-20
JP4434905B2 JP4434905B2 (en) 2010-03-17

Family

ID=36376576

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004289052A Expired - Fee Related JP4434905B2 (en) 2004-09-30 2004-09-30 Process management device

Country Status (1)

Country Link
JP (1) JP4434905B2 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008065821A (en) * 2006-09-06 2008-03-21 Fisher Rosemount Syst Inc Monitoring system and monitoring method
JP2009299934A (en) * 2008-06-10 2009-12-24 Samsung Electronics Co Ltd Failure diagnosing device using refrigerating cycle device
WO2011027607A1 (en) * 2009-09-07 2011-03-10 株式会社日立製作所 Anomaly detection and diagnostic method, anomaly detection and diagnostic system, and anomaly detection and diagnostic program
KR20170007408A (en) 2014-05-20 2017-01-18 도시바 미쓰비시덴키 산교시스템 가부시키가이샤 Manufacturing equipment diagnosis support system
US10436468B2 (en) 2014-03-24 2019-10-08 Nec Corporation Monitoring device, monitoring system, monitoring method, and non-transitory storage medium
WO2020090767A1 (en) * 2018-10-30 2020-05-07 国立研究開発法人宇宙航空研究開発機構 Abnormality diagnostic device, abnormality diagnostic method, and program
KR102261941B1 (en) * 2019-11-29 2021-06-08 한국생산기술연구원 Time series data fingerprinting and image deep learning analysis methods for product defect detection
KR20220059691A (en) * 2020-11-03 2022-05-10 한국생산기술연구원 Space filling curve transformation method and system for deep learning analysis of manufacturing time series data

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1078471A (en) * 1996-09-03 1998-03-24 Furukawa Electric Co Ltd:The Measuring method for partial discharge
JPH10258974A (en) * 1996-09-13 1998-09-29 Toshiba Corp Unsteady signal analyzer and medium recording unsteady signal analysis program
JPH10288639A (en) * 1997-02-13 1998-10-27 Toshiba Corp Apparatus and method for monitoring malfunction of electric apparatus
JP2000050488A (en) * 1998-07-28 2000-02-18 Kansai Tech Corp Method for discrimination cause of ground failure of high-voltage power distribution line
JP2003522577A (en) * 2000-02-18 2003-07-29 アーゴス インク Multivariate analysis of green and ultraviolet spectra of cell and tissue samples
JP2004234302A (en) * 2003-01-30 2004-08-19 Toshiba Corp Process management device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1078471A (en) * 1996-09-03 1998-03-24 Furukawa Electric Co Ltd:The Measuring method for partial discharge
JPH10258974A (en) * 1996-09-13 1998-09-29 Toshiba Corp Unsteady signal analyzer and medium recording unsteady signal analysis program
JPH10288639A (en) * 1997-02-13 1998-10-27 Toshiba Corp Apparatus and method for monitoring malfunction of electric apparatus
JP2000050488A (en) * 1998-07-28 2000-02-18 Kansai Tech Corp Method for discrimination cause of ground failure of high-voltage power distribution line
JP2003522577A (en) * 2000-02-18 2003-07-29 アーゴス インク Multivariate analysis of green and ultraviolet spectra of cell and tissue samples
JP2004234302A (en) * 2003-01-30 2004-08-19 Toshiba Corp Process management device

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008065821A (en) * 2006-09-06 2008-03-21 Fisher Rosemount Syst Inc Monitoring system and monitoring method
JP2009299934A (en) * 2008-06-10 2009-12-24 Samsung Electronics Co Ltd Failure diagnosing device using refrigerating cycle device
WO2011027607A1 (en) * 2009-09-07 2011-03-10 株式会社日立製作所 Anomaly detection and diagnostic method, anomaly detection and diagnostic system, and anomaly detection and diagnostic program
JP2011059790A (en) * 2009-09-07 2011-03-24 Hitachi Ltd Fault detection/diagnosis method, fault detection/diagnosis system, and fault detection/diagnosis program
US9483049B2 (en) 2009-09-07 2016-11-01 Hitachi, Ltd. Anomaly detection and diagnosis/prognosis method, anomaly detection and diagnosis/prognosis system, and anomaly detection and diagnosis/prognosis program
US10436468B2 (en) 2014-03-24 2019-10-08 Nec Corporation Monitoring device, monitoring system, monitoring method, and non-transitory storage medium
KR20170007408A (en) 2014-05-20 2017-01-18 도시바 미쓰비시덴키 산교시스템 가부시키가이샤 Manufacturing equipment diagnosis support system
US10996662B2 (en) 2014-05-20 2021-05-04 Toshiba Mitsubishi-Electric Industrial Systems Corporation Manufacturing equipment diagnosis support system
WO2020090767A1 (en) * 2018-10-30 2020-05-07 国立研究開発法人宇宙航空研究開発機構 Abnormality diagnostic device, abnormality diagnostic method, and program
US11692910B2 (en) 2018-10-30 2023-07-04 Japan Aerospace Exploration Agency Abnormality diagnostic device, abnormality diagnostic method, and program
KR102261941B1 (en) * 2019-11-29 2021-06-08 한국생산기술연구원 Time series data fingerprinting and image deep learning analysis methods for product defect detection
KR20220059691A (en) * 2020-11-03 2022-05-10 한국생산기술연구원 Space filling curve transformation method and system for deep learning analysis of manufacturing time series data
KR102558232B1 (en) 2020-11-03 2023-07-25 한국생산기술연구원 Space filling curve transformation method and system for deep learning analysis of manufacturing time series data

Also Published As

Publication number Publication date
JP4434905B2 (en) 2010-03-17

Similar Documents

Publication Publication Date Title
Hu et al. A statistical training data cleaning strategy for the PCA-based chiller sensor fault detection, diagnosis and data reconstruction method
Li et al. A data-driven strategy for detection and diagnosis of building chiller faults using linear discriminant analysis
Hu et al. Chiller sensor fault detection using a self-adaptive principal component analysis method
Wang et al. An online fault diagnosis tool of VAV terminals for building management and control systems
US9740545B2 (en) Equipment evaluation device, equipment evaluation method and non-transitory computer readable medium
Wang et al. Enhanced chiller fault detection using Bayesian network and principal component analysis
Yan et al. A sensor fault detection strategy for air handling units using cluster analysis
Wang et al. Online model-based fault detection and diagnosis strategy for VAV air handling units
Yu et al. A review of fault detection and diagnosis methodologies on air-handling units
Du et al. Fault detection and diagnosis for buildings and HVAC systems using combined neural networks and subtractive clustering analysis
Wang et al. A robust fault detection and diagnosis strategy for multiple faults of VAV air handling units
US9664400B2 (en) Automated technique of measuring room air change rates in HVAC system
Wang et al. A robust fault detection and diagnosis strategy for pressure-independent VAV terminals of real office buildings
Zhu et al. Fault diagnosis based operation risk evaluation for air conditioning systems in data centers
CN111047732A (en) Equipment abnormity diagnosis method and device based on energy consumption model and data interaction
Harmer et al. Using calibrated energy models for building commissioning and load prediction
O’Neill et al. Real time model-based energy diagnostics in buildings
Sun et al. Optimization of support vector regression model based on outlier detection methods for predicting electricity consumption of a public building WSHP system
WO2018112352A1 (en) Techniques of automated fault detection and related systems and methods
Yang et al. A hybrid model-based fault detection strategy for air handling unit sensors
Kang et al. Bayesian-Emulator based parameter identification for calibrating energy models for existing buildings
JP4434905B2 (en) Process management device
JP2005301582A (en) Process management device
Papadopoulos et al. Distributed diagnosis of sensor and actuator faults in air handling units in multi-zone buildings: a model-based approach
JP2004234302A (en) Process management device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070306

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090519

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090521

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090901

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091029

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20091124

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20091222

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

Free format text: PAYMENT UNTIL: 20130108

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees