JP2020095451A - Parameter estimation device and method - Google Patents

Parameter estimation device and method Download PDF

Info

Publication number
JP2020095451A
JP2020095451A JP2018232579A JP2018232579A JP2020095451A JP 2020095451 A JP2020095451 A JP 2020095451A JP 2018232579 A JP2018232579 A JP 2018232579A JP 2018232579 A JP2018232579 A JP 2018232579A JP 2020095451 A JP2020095451 A JP 2020095451A
Authority
JP
Japan
Prior art keywords
parameter
likelihood
target parameter
estimated
target
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.)
Pending
Application number
JP2018232579A
Other languages
Japanese (ja)
Inventor
民圭 曹
Minkyu Jo
民圭 曹
良和 石井
Yoshikazu Ishii
良和 石井
新治 村田
Shinji Murata
新治 村田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2018232579A priority Critical patent/JP2020095451A/en
Publication of JP2020095451A publication Critical patent/JP2020095451A/en
Pending legal-status Critical Current

Links

Images

Abstract

To improve accuracy of estimating an unknown parameter related to a facility.SOLUTION: A parameter estimation device 1 estimating a parameter related to a facility 3 acquires an observed value from a sensor 2 measuring a physical amount related to the facility, and acquires a facility model formula to calculate an observed value based on a target parameter (41), calculates likelihood of an estimated observed value, which is obtained by substituting the estimated value of the target parameter into the facility model formula, to the observed value by a first likelihood function, calculates likelihood of the estimated value of the target parameter based on the first likelihood function and a second likelihood function for obtaining probability that the estimated value of the target parameter is within a predetermined range (42), estimates the target parameter based on the calculated likelihood (43), and evaluates the estimated target parameter according to a predetermined evaluation criterion (44).SELECTED DRAWING: Figure 1

Description

本発明は、パラメータ推定装置および方法に関する。 The present invention relates to a parameter estimation device and method.

エネルギを効率的に使用するため、ビルや商店施設にBEMS(Building Energy Management System)が普及しつつある。BEMSは、建物の熱源設備、空調設備、電力設備等のユーティリティ設備を管理する。BEMSは、ユーティリティ設備にセンサおよび通信装置を付けて計測し、見える化し、管理することで、設備の保全および設備の最適化を実行する。 BEMS (Building Energy Management System) is becoming widespread in buildings and commercial facilities in order to use energy efficiently. BEMS manages utility equipment such as heat source equipment, air conditioning equipment, and electric power equipment of a building. BEMS carries out maintenance of facilities and optimization of facilities by attaching sensors and communication devices to utility facilities, measuring, visualizing and managing them.

さらに近年では、CEMS(Community Energy Management System)の重要性が高まっている。CEMSは、複数のBEMSを管理し、地域全体のエネルギを管理する。BEMSおよびCEMSの主要な機能は、設備状態の見える化である。設備状態とは、消費電力だけではなく、効率、流量、温度、圧力など、BEMSおよびCEMSが管理するユーティリティ設備の動作状態を表す値である。しかし、設備の持つすべての状態を知るために、多種多様なセンサをユーティリティ設備に設けることは、高コストであり、現実的ではない。 Furthermore, in recent years, the importance of CEMS (Community Energy Management System) is increasing. CEMS manages several BEMS and manages the energy of the whole area. The main function of BEMS and CEMS is visualization of equipment status. The equipment state is a value representing not only power consumption but also an operation state of utility equipment managed by BEMS and CEMS, such as efficiency, flow rate, temperature, and pressure. However, it is expensive and impractical to provide various kinds of sensors in the utility equipment in order to know all the states of the equipment.

そこで、一つの状態を検知するセンサデータに基づいて、他の状態の値を推定することが考えられる。すなわち、BEMSまたはCEMSの管理下にある計測データを用いて、未知の状態を推定することにより、検知対象の状態毎にセンサを取り付けずに、設備状態を管理することができる。 Therefore, it is conceivable to estimate the value of another state based on the sensor data for detecting one state. That is, by estimating the unknown state by using the measurement data under the control of BEMS or CEMS, the equipment state can be managed without attaching a sensor for each state of the detection target.

ところで、ユーティリティ設備の状態を推定するためには、パラメータを推定する必要がある。しかし、未知のパラメータを推定する場合、物理的な値ではない値に推定値が収束する場合がある。未知のパラメータの推定精度を向上する技術として、特許文献1の技術が提案されている By the way, in order to estimate the state of utility equipment, it is necessary to estimate parameters. However, when estimating an unknown parameter, the estimated value may converge to a value that is not a physical value. As a technique for improving the estimation accuracy of unknown parameters, the technique of Patent Document 1 has been proposed.

特許文献1では、「分電盤接続されている複数の電気機器による消費電力の合計値を用いて、前記複数の電気機器の稼動状態を確率モデルによりモデル化したときのモデルパラメータを推定するパラメータ推定ステップを含み、前記パラメータ推定ステップでは、前記複数の電気機器それぞれの動作傾向から事前知識として予め定めることのできる電力データの特徴に基づいて、尤度関数によって算出される尤度が最大になるようなモデルパラメータを推定し、前記確率モデルは、Factorial HMM(Hidden Markov Model)であり、前記尤度は、前記Factorial HMMによってモデル化された前記消費電力の合計値のパターンの、実測した前記消費電力の合計値に対する確からしさを示す値である。」という記載がある。 In patent document 1, "a parameter for estimating a model parameter when a working model of a plurality of electric devices is modeled by a probabilistic model using a total value of power consumption by a plurality of electric devices connected to a distribution board. In the parameter estimation step, including the estimation step, the likelihood calculated by the likelihood function is maximized based on the characteristics of the power data that can be predetermined as prior knowledge from the operation tendency of each of the plurality of electric devices. Such a model parameter is estimated, the probabilistic model is a Factorial HMM (Hidden Markov Model), and the likelihood is the actually measured consumption of the total power consumption pattern modeled by the Factorial HMM. It is a value indicating the certainty with respect to the total value of electric power."

特開2017−49221号公報JP, 2017-49221, A

しかし、特許文献1では、Factorial HMM以外のモデルを用いる場合に限定されてしまい、他のモデルを用いて未知パラメータを推定することは難しい。 However, Patent Document 1 is limited to the case of using a model other than the Physical HMM, and it is difficult to estimate the unknown parameter using another model.

ところで、ユーティリティ設備モデルの多数のパラメータを容易に推定するためには、観測値の種類の数が多い必要がある。観測値とは、センサなどの計測装置により実測された値である。しかし多種類の観測値を得るためには、上述のように、センサなどの計測装置をユーティリティ設備に取り付ける必要があり、コスト上の問題がある。 By the way, in order to easily estimate many parameters of the utility equipment model, it is necessary to have a large number of types of observation values. The observed value is a value actually measured by a measuring device such as a sensor. However, in order to obtain many kinds of observation values, it is necessary to attach a measuring device such as a sensor to the utility equipment as described above, which causes a cost problem.

さらに、観測値の種類の数が少ない場合は、いわゆる劣問題となるため、推定値が収束しない可能性がある。したがって、限られた種類の観測値を用いて、真値に近い推定値を取得する手法が必要である。 Furthermore, when the number of types of observed values is small, it becomes a so-called inferior problem, and thus the estimated values may not converge. Therefore, there is a need for a method of obtaining an estimated value close to the true value by using a limited number of observation values.

本発明は、上記課題に鑑みてなされたもので、その目的は、設備に関する未知のパラメータの推定精度を向上できるようにしたパラメータ推定装置および方法を提供することにある。 The present invention has been made in view of the above problems, and an object thereof is to provide a parameter estimation device and method capable of improving the estimation accuracy of an unknown parameter regarding equipment.

上記課題を解決すべく、本発明に従うパラメータ推定装置は、設備に関するパラメータを推定する装置であって、設備に関する物理量を測定するセンサから観測値を取得し、観測値を対象パラメータに基づいて算出する設備モデル式を取得し、設備モデル式に対象パラメータの推定値を代入して得られる推定観測値の観測値に対する尤度を第1尤度関数により算出し、対象パラメータの推定値があらかじめ用意された所定範囲内に属する確率を求める第2尤度関数と第1尤度関数とに基づいて、対象パラメータの推定値の尤度を算出し、算出された尤度に基づいて対象パラメータを推定し、推定された対象パラメータを所定の評価基準にしたがって評価する。 In order to solve the above-mentioned problems, a parameter estimation device according to the present invention is a device for estimating a parameter related to equipment, and acquires an observed value from a sensor that measures a physical quantity related to the equipment, and calculates the observed value based on a target parameter. Obtain the equipment model formula and calculate the likelihood of the estimated observation value obtained by substituting the estimated value of the target parameter into the equipment model formula by the first likelihood function, and the estimated value of the target parameter is prepared in advance. The likelihood of the estimated value of the target parameter is calculated based on the second likelihood function and the first likelihood function for obtaining the probability of belonging to the predetermined range, and the target parameter is estimated based on the calculated likelihood. , The estimated target parameter is evaluated according to a predetermined evaluation criterion.

本発明によれば、第1尤度関数と第2尤度関数とに基づいて算出される対象パラメータの推定値の尤度から、対象パラメータを推定し、評価することができる。 According to the present invention, the target parameter can be estimated and evaluated from the likelihood of the estimated value of the target parameter calculated based on the first likelihood function and the second likelihood function.

パラメータ推定装置を含むシステム全体の構成図である。It is a block diagram of the whole system including a parameter estimation device. パラメータ推定装置の構成図である。It is a block diagram of a parameter estimation apparatus. パラメータ推定処理を示すフローチャートである。It is a flowchart which shows a parameter estimation process. 設備モデル式を決定する処理を示すフローチャートである。It is a flowchart which shows the process which determines an equipment model type|formula. 対象パラメータの適正範囲を設定する処理を示すフローチャートである。It is a flow chart which shows the processing which sets up the proper range of a target parameter. 対象パラメータの適正範囲の設定例を示す説明図である。It is explanatory drawing which shows the example of setting of the suitable range of a target parameter. 対象パラメータの適正範囲の他の設定例を示す説明図である。It is explanatory drawing which shows the other example of setting of the suitable range of a target parameter. 対象パラメータの適正範囲のさらに別の設定例を示す説明図である。It is explanatory drawing which shows another example of setting of the appropriate range of a target parameter. 第2実施例に係り、適正範囲の設定処理を示すフローチャートである。9 is a flow chart showing a setting process of an appropriate range according to the second embodiment.

以下、図面に基づいて、本発明の実施の形態を説明する。本実施形態に係るパラメータ推定装置では、対象パラメータ(以下、パラメータまたは未知パラメータと呼ぶ場合がある)の推定値が所定範囲での分布密度が高くなるように、尤度を設定して学習する。 Hereinafter, embodiments of the present invention will be described with reference to the drawings. In the parameter estimation device according to the present embodiment, the likelihood is set and learning is performed so that the estimated value of the target parameter (hereinafter, sometimes referred to as parameter or unknown parameter) has a high distribution density in a predetermined range.

本実施形態では、未知パラメータの推定値を確率的に自動算出し、評価する。パラメータの推定値が、所定範囲において尤度が高くなる所定の分布に従うと仮定することにより、対象パラメータの良否を判定する。 In this embodiment, the estimated value of the unknown parameter is stochastically automatically calculated and evaluated. The quality of the target parameter is determined by assuming that the estimated value of the parameter follows a predetermined distribution in which the likelihood increases in a predetermined range.

本実施形態によれば、パラメータの適正範囲に基づく尤度関数を設定し、観測値の推定値の尤度関数と適正範囲に基づく尤度関数とを用いてパラメータを推定し、パラメータの良否を判断することができる。したがって、本実施形態では、ユーティリティ設備における観測値の種類の数よりも多い数のパラメータを推定する場合に、真値に近い推定値を算出することができる。 According to the present embodiment, the likelihood function based on the appropriate range of the parameter is set, the parameter is estimated using the likelihood function of the estimated value of the observed value and the likelihood function based on the appropriate range, and the quality of the parameter is determined. You can judge. Therefore, in the present embodiment, when estimating a larger number of parameters than the number of types of observed values in utility equipment, it is possible to calculate an estimated value close to a true value.

本実施形態に係るパラメータ推定装置およびパラメータ推定方法は、例えば、以下のように表現することもできる。 The parameter estimation device and the parameter estimation method according to the present embodiment can also be expressed as follows, for example.

「表現1」ユーティリティ設備のパラメータ推定装置において、パラメータの適正範囲に基づき、尤度関数を設定し、パラメータを推定し、パラメータを判定するパラメータ推定装置。 "Expression 1" A parameter estimating device for utility equipment, which sets a likelihood function, estimates a parameter, and determines a parameter based on an appropriate range of the parameter.

「表現2」ユーティリティ設備のパラメータ推定方法において、パラメータの適正範囲に基づき、尤度関数を設定し、パラメータを推定し、パラメータを判定するパラメータ推定方法。 “Expression 2” In the parameter estimation method for utility equipment, a parameter estimation method that sets a likelihood function, estimates the parameter, and determines the parameter based on the appropriate range of the parameter.

「表現3」表現2記載のパラメータ推定方法において、一つの連続した適正範囲に基づく尤度設定をするステップを含むパラメータ推定方法。 “Expression 3” The parameter estimation method according to Expression 2, including the step of performing likelihood setting based on one continuous appropriate range.

「表現4」表現2記載のパラメータ推定方法において、複数の連続区間を有する適正範囲に基づく尤度設定をするステップを含むパラメータ推定方法。 “Expression 4” The parameter estimation method according to Expression 2, including a step of performing likelihood setting based on an appropriate range having a plurality of continuous sections.

「表現5」表現2記載のパラメータ推定方法において、推定パラメータに対して、一定時間の変化量の適正範囲に基づく尤度設定をするステップを含むパラメータ推定方法。 Expression 5” The parameter estimation method described in Expression 2, including a step of performing likelihood setting on the estimated parameter based on an appropriate range of the amount of change for a certain period of time.

「表現6」表現2記載のパラメータ推定方法において、設備モデル式を記憶部から検索し、推定するパラメータを決定するステップを含むパラメータ推定方法。 “Expression 6” The parameter estimation method according to Expression 2, including the step of searching the storage unit for the equipment model formula and determining the parameter to be estimated.

「表現7」ユーティリティ設備のパラメータを推定する方法において、複数のパラメータの適正範囲に基づき、複数の尤度関数の組合せ生成して設定し、パラメータを推定し、パラメータを判定するパラメータ推定方法。 [Expression 7] A method of estimating parameters of utility equipment, which is a method of estimating parameters by estimating and estimating parameters by generating and setting combinations of a plurality of likelihood functions based on appropriate ranges of the parameters.

「表現8」表現2記載のパラメータ推定方法において、各尤度設定の組合せで、適正範囲を用いてないパラメータに対しても、パラメータ判定ステップで、適正範囲設定ステップで設定した適正範囲を適用し、良否を判定するパラメータ推定方法。 [Expression 8] In the parameter estimation method described in Expression 2, even in a combination of likelihood settings, even if a parameter does not use an appropriate range, the appropriate range set in the appropriate range setting step is applied in the parameter determination step. , A parameter estimation method for determining pass/fail.

図1〜図8を用いて第1実施例を説明する。本実施例では、ビルディング内のユーティリティ設備についての未知パラメータを推定する場合を例に挙げて説明するが、これに限らず、例えば、工場の生産設備、商業施設の冷蔵設備、遊戯設備、エレベータ、エスカレータ、食品配送設備等の他の施設の他の設備にも適用することができる。 A first embodiment will be described with reference to FIGS. In the present embodiment, the case of estimating the unknown parameter for the utility equipment in the building will be described as an example, but the present invention is not limited to this. For example, factory production equipment, refrigeration equipment for commercial facilities, play equipment, elevators, It can also be applied to other facilities of other facilities such as escalator and food delivery facility.

<全体の構成> <Overall structure>

図1は、パラメータ推定装置1を含むシステムの全体構成を示す。パラメータ推定装置1は、例えば、BEMS4とセンサ2とに通信ネットワークCNを介して接続される。センサ2は、ユーティリティ設備(以下、設備)3と接続されており、計測情報(センサデータ、観測値とも呼ぶ)を、BEMS4またはパラメータ推定装置1へ送信する。 FIG. 1 shows the overall configuration of a system including the parameter estimation device 1. The parameter estimation device 1 is connected to, for example, the BEMS 4 and the sensor 2 via the communication network CN. The sensor 2 is connected to the utility equipment (hereinafter, equipment) 3 and transmits measurement information (also referred to as sensor data or observation value) to the BEMS 4 or the parameter estimation device 1.

設備3は、例えば、ボイラ設備、空調装置、ガスタービン発電装置、電力関連機器等である。センサ2は、設備3の各種物理量を計測し、その計測結果である計測情報を出力する装置である。センサ2は、例えば、消費電力、効率、流量、温度、圧力、液位、質量、音、光量などを計測する。色彩、味、匂いといった心理物理量を計測するセンサであってもよい。 The facility 3 is, for example, a boiler facility, an air conditioner, a gas turbine power generator, a power-related device, or the like. The sensor 2 is a device that measures various physical quantities of the equipment 3 and outputs measurement information that is the measurement result. The sensor 2 measures, for example, power consumption, efficiency, flow rate, temperature, pressure, liquid level, mass, sound, light amount, and the like. It may be a sensor that measures psychophysical quantities such as color, taste, and smell.

センサ2で計測されたデータ(計測情報)は、通信ネットワークCNを介してパラメータ推定装置1へ直接的に送信されてもよいし、BEMS4などの装置と通信ネットワークCNとを介してパラメータ推定装置1へ送信されてもよい。さらに、センサ2は、定期的にデータを送信してもよいし、パラメータ推定装置1またはBEMS4からの要求に応じてデータを送信してもよい。あるいは、センサ2の計測したデータを記録媒体に格納し、記録媒体からパラメータ推定装置1へ読み込ませてもよい。 The data (measurement information) measured by the sensor 2 may be directly transmitted to the parameter estimation device 1 via the communication network CN, or may be transmitted via a device such as the BEMS 4 and the communication network CN. May be sent to. Furthermore, the sensor 2 may transmit data regularly, or may transmit data in response to a request from the parameter estimation device 1 or the BEMS 4. Alternatively, the data measured by the sensor 2 may be stored in a recording medium and read into the parameter estimation device 1 from the recording medium.

なお、本実施例では、BEMS4を例に挙げるが、これに限らず、CEMSなどのユーティリティ設備3を管理するシステムであればよい。 In the present embodiment, the BEMS 4 is taken as an example, but the present invention is not limited to this, and any system that manages the utility equipment 3 such as CEMS may be used.

パラメータ推定装置1は、例えばCPU(マイクロプロセッサ)、メモリ、通信機能、表示機能、入力機能(いずれも図示せず)などを備えた装置である。パラメータ推定装置1は、パーソナルコンピュータを利用して構成することもできるし、専用コンピュータとして構成することもできる。 The parameter estimation device 1 is a device including, for example, a CPU (microprocessor), a memory, a communication function, a display function, an input function (none of which is shown). The parameter estimation device 1 can be configured using a personal computer or a dedicated computer.

パラメータ推定装置1の機能構成の詳細は、図2と共に後述する。ここでは、パラメータ推定装置1の概要を先に説明する。パラメータ推定装置1は、設備モデル式決定部41と、適正範囲設定部42と、パラメータ推定部43と、パラメータ判定部44といった機能を有している。 Details of the functional configuration of the parameter estimation device 1 will be described later with reference to FIG. Here, the outline of the parameter estimation device 1 will be described first. The parameter estimation device 1 has functions such as an equipment model formula determination unit 41, an appropriate range setting unit 42, a parameter estimation unit 43, and a parameter determination unit 44.

図2は、パラメータ推定装置1の機能構成を示す。パラメータ推定装置1は、例えば、取得部10と、インタフェース部20と、記憶部30と、演算部40とを含む。これら各機能10〜40は、メモリに格納された所定のコンピュータプログラム(不図示)をマイクロプロセッサが読み込んで実行することにより、実現される。機能10〜40のいずれか少なくとも一部をハードウェア回路として構成してもよい。 FIG. 2 shows a functional configuration of the parameter estimation device 1. The parameter estimation device 1 includes, for example, an acquisition unit 10, an interface unit 20, a storage unit 30, and a calculation unit 40. Each of these functions 10 to 40 is realized by the microprocessor reading and executing a predetermined computer program (not shown) stored in the memory. At least a part of any one of the functions 10 to 40 may be configured as a hardware circuit.

インタフェース部20は、例えば、入力部21と出力部22とを有する。演算部40は、例えば、設備モデル式決定部41と、適正範囲設定部42と、パラメータ推定部43と、パラメータ判定部44とを備える。 The interface unit 20 has, for example, an input unit 21 and an output unit 22. The calculation unit 40 includes, for example, an equipment model formula determination unit 41, an appropriate range setting unit 42, a parameter estimation unit 43, and a parameter determination unit 44.

取得部10は、設備3のセンサ2、BEMS4(あるいはCEMS)、インターネット上の情報源(不図示)からデータ(情報とも呼ぶ)を取得する。取得されたデータは、記憶部30に格納される。 The acquisition unit 10 acquires data (also referred to as information) from the sensor 2 of the equipment 3, the BEMS 4 (or CEMS), and an information source (not shown) on the Internet. The acquired data is stored in the storage unit 30.

インタフェース部20の入力部21は、パラメータ推定装置1に対してユーザが情報を入力するための機能である。ユーザは、例えばキーボードまたはマウスなどの入力装置(不図示)を用いて、設備3などの情報を入力部21へ入力する。キーボードなどに代えて、あるいはキーボードなどと共に、音声認識装置などの入力装置を用いてもよい。 The input unit 21 of the interface unit 20 is a function for the user to input information to the parameter estimation device 1. The user uses the input device (not shown) such as a keyboard or a mouse to input information about the facility 3 and the like to the input unit 21. An input device such as a voice recognition device may be used instead of or in addition to the keyboard.

入力部21に入力された設備3などの情報は、記憶部30に転送されて記憶される。入力部21に入力された情報の一部または全部を、出力部22により表示出力させることもできる。さらに、入力部21に入力された情報は、記憶部30を介して演算部40に送られる。 Information of the equipment 3 and the like input to the input unit 21 is transferred to and stored in the storage unit 30. A part or all of the information input to the input unit 21 can be displayed and output by the output unit 22. Further, the information input to the input unit 21 is sent to the calculation unit 40 via the storage unit 30.

入力される情報としては、例えば、設備3の仕様、設備3の設置場所、設備3と他の設備との接続関係、設備3について推定するパラメータがある。入力情報の例は、これらに限定されない。 The input information includes, for example, the specifications of the equipment 3, the installation location of the equipment 3, the connection relationship between the equipment 3 and other equipment, and parameters for estimating the equipment 3. Examples of input information are not limited to these.

出力部22は、ユーザに対して情報を提供する機能である。出力部22は、例えば、ディスプレイ、プリンタ、音声合成装置などの出力装置(不図示)を用いて、パラメータの推定結果などを表示する。出力部22は、入力部21へ入力された情報を表示したり、演算部40から記憶部30に格納された情報を読み出して表示したり、記憶部30に記憶されている情報を読み出して表示したりすることができる。 The output unit 22 has a function of providing information to the user. The output unit 22 displays the parameter estimation result and the like using an output device (not shown) such as a display, a printer, and a voice synthesizer. The output unit 22 displays the information input to the input unit 21, reads and displays the information stored in the storage unit 30 from the calculation unit 40, and reads and displays the information stored in the storage unit 30. You can

記憶部30は、取得部10と、インタフェース部20と、演算部40との間で、情報を送受信する。記憶部30は、入力された情報を格納すると共に、要求された情報を取り出して要求元へ送信する。記憶部30には、設備3の仕様情報などのほかに、後述する設備モデル式と適正範囲とパラメータ初期値などが格納される。 The storage unit 30 transmits/receives information among the acquisition unit 10, the interface unit 20, and the calculation unit 40. The storage unit 30 stores the input information, extracts the requested information, and transmits it to the request source. In addition to the specification information of the equipment 3, the storage unit 30 stores an equipment model formula (to be described later), an appropriate range, parameter initial values, and the like.

設備モデル式とは、取得部10から取得するデータとインタフェース部20から取得する推定パラメータとの関係を示す式である。設備モデル式は、取得部10から取得されるデータと推定パラメータとの関係を表現可能な複数の式として構成することができ、物理式を含んでもよい。「所定範囲」としての適正範囲とは、設備種類、設備仕様、設置環境といったタグ(種類情報)が設定された、各パラメータの適正範囲である。 The equipment model formula is a formula indicating the relationship between the data acquired from the acquisition unit 10 and the estimation parameter acquired from the interface unit 20. The equipment model formula can be configured as a plurality of formulas that can represent the relationship between the data acquired from the acquisition unit 10 and the estimated parameter, and may include a physical formula. The appropriate range as the “predetermined range” is an appropriate range of each parameter in which tags (type information) such as equipment type, equipment specifications, and installation environment are set.

演算部40は、設備モデル式決定部41と、適正範囲設定部42と、パラメータ推定部43と、パラメータ判定部44とを含む。演算部40は、インタフェース部20と記憶部30とから演算に必要な情報を取得して対象パラメータを推定し、インタフェース部20に推定結果を送信する。対象パラメータの推定結果は、出力部22に表示される。 The calculation unit 40 includes an equipment model formula determination unit 41, an appropriate range setting unit 42, a parameter estimation unit 43, and a parameter determination unit 44. The calculation unit 40 acquires information necessary for calculation from the interface unit 20 and the storage unit 30, estimates the target parameter, and transmits the estimation result to the interface unit 20. The estimation result of the target parameter is displayed on the output unit 22.

パラメータの推定結果を含む情報を、設備導入のシミュレーションシステムなどの外部装置(不図示)へ提供してもよい。パラメータ推定装置1から外部装置へ有線または無線により情報を転送してもよいし、記録媒体に格納して外部装置へ提供してもよい。 Information including the parameter estimation result may be provided to an external device (not shown) such as a facility introduction simulation system. Information may be transferred from the parameter estimation device 1 to an external device by wire or wirelessly, or may be stored in a recording medium and provided to the external device.

<パラメータ推定装置1の全体の処理の概要> <Outline of Overall Processing of Parameter Estimation Device 1>

図3は、パラメータ推定装置1の全体処理を示すフローチャートである。ステップS12〜S15の詳細は後述する。 FIG. 3 is a flowchart showing the overall processing of the parameter estimation device 1. Details of steps S12 to S15 will be described later.

パラメータ推定装置1は、インタフェース部20と記憶部30とからデータを取得する(S11)。パラメータ推定装置1は、ステップS11で取得されたデータに基づいて、設備モデル式を決定する(S12)。パラメータ推定装置1は、ステップS11で取得されたデータに基づいて、適正範囲を設定する(S13)。 The parameter estimation device 1 acquires data from the interface unit 20 and the storage unit 30 (S11). The parameter estimation device 1 determines an equipment model formula based on the data acquired in step S11 (S12). The parameter estimation device 1 sets an appropriate range based on the data acquired in step S11 (S13).

ここで、ステップS12とステップS13とは、並列で処理されてもいい。ステップS12を先に処理し、ステップS13を後で処理してもよい。ステップS12とステップS13の処理の順番は問わない。 Here, step S12 and step S13 may be processed in parallel. Step S12 may be processed first, and step S13 may be processed later. The order of the processes of step S12 and step S13 does not matter.

パラメータ推定装置1は、ステップS12で決定された設備モデル式とステップS13で設定された適正範囲とに基づいて、対象パラメータを推定する(S14)。 The parameter estimation device 1 estimates the target parameter based on the equipment model formula determined in step S12 and the appropriate range set in step S13 (S14).

パラメータ推定装置1は、ステップS14で算出されたパラメータの推定値を評価する(S15)。評価とは、例えば、算出されたパラメータの推定値の良否を判定することである。 The parameter estimation device 1 evaluates the estimated value of the parameter calculated in step S14 (S15). The evaluation is, for example, to judge the quality of the estimated value of the calculated parameter.

パラメータ推定装置1は、ステップS15で「否」と判定されると、ステップS14へ戻ってパラメータ推定を再度行う。一方、パラメータ推定装置1は、ステップS15で「良」と判定された場合、ステップS14で算出されたパラメータの推定値を出力する(S16)。 When the parameter estimation device 1 determines “NO” in step S15, the parameter estimation device 1 returns to step S14 and performs parameter estimation again. On the other hand, when the parameter estimation device 1 determines "good" in step S15, the parameter estimation device 1 outputs the estimated value of the parameter calculated in step S14 (S16).

<設備モデル式決定部41の処理> <Process of Equipment Model Formula Determination Unit 41>

図4は、図3中にステップS12で示された、設備モデル式決定部41での処理の詳細を示すフローチャートである。 FIG. 4 is a flowchart showing details of the processing in the equipment model formula determination unit 41 shown in step S12 in FIG.

設備モデル式決定部41は、図3のステップS11で取得されたデータに基づいて、記憶部30に格納されている設備モデル式の中から演算に使用する設備モデル式を一つ以上検索する(S121)。 The equipment model formula determination unit 41 searches for one or more facility model formulas to be used for calculation from the facility model formulas stored in the storage unit 30 based on the data acquired in step S11 of FIG. 3 ( S121).

記憶部30に格納されている設備モデル式には、例えば、「設備種類」、「式の左辺のパラメータ」、「式の右辺のパラメータ」といったタグが付いている。そこで、設備モデル式決定部41は、タグに基づいて設備モデル式を検索する。設備モデル式の検索には、クラスタリングを用いればいいが、この限りではない。 The equipment model formulas stored in the storage unit 30 are tagged with, for example, “equipment type”, “parameter on left side of formula”, and “parameter on right side of formula”. Therefore, the equipment model formula determination unit 41 searches for the equipment model formula based on the tag. Clustering may be used to retrieve the equipment model formula, but this is not the only option.

設備モデル式決定部41は、ステップS121で検索された設備モデル式のパラメータを確認し、推定対象のパラメータを決める(S122)。以下、ステップS121とステップS122の一例を、蒸気ボイラのガス消費量を計算する式を用いて説明する。 The equipment model formula determination unit 41 confirms the parameters of the equipment model formula retrieved in step S121 and determines parameters to be estimated (S122). Hereinafter, an example of step S121 and step S122 will be described using an equation for calculating the gas consumption of the steam boiler.

G=(H+Hb)/Hrout×Grout・・・式(1) G=(H+Hb)/Hrout×Grout (1)

H=f(S,Es,Ew) ・・・式(2) H=f(S, Es, Ew)... Formula (2)

Hb=f(S,Ek,Ew,Br)・・・式(3) Hb=f(S, Ek, Ew, Br)... Formula (3)

Es=f(P) ・・・式(4) Es=f(P)... Formula (4)

Ew=f(Tw) ・・・式(5) Ew=f(Tw)... Formula (5)

Ek=f(Tk) ・・・式(6) Ek=f(Tk) Equation (6)

Tk=g(P) ・・・式(7) Tk=g(P)... Formula (7)

ここで、「G」はガス消費量、「H」は熱量、「Hb」はブロー損失熱量、「Hrout」は定格熱出力、「Grout」は定格ガス消費量、「S」は蒸気量、「Es」は蒸気エンタルピ、「Ew」は給水エンタルピ、「Ek」は缶水エンタルピ、「Br」はブロー係数、「P」は使用圧力、「Tw」は給水温度、「Tk」は缶水温度である。 Here, "G" is the gas consumption amount, "H" is the heat amount, "Hb" is the blow loss heat amount, "Hrout" is the rated heat output, "Grout" is the rated gas consumption amount, "S" is the steam amount, " "Es" is steam enthalpy, "Ew" is water supply enthalpy, "Ek" is can water enthalpy, "Br" is blow coefficient, "P" is working pressure, "Tw" is water supply temperature, and "Tk" is water temperature. is there.

各式は、物理モデル式または近似式である。例えば、給水エンタルピEwの式(5)は給水温度の二次近似式であり、ガス消費量Gの式(1)は物理モデル式である。 Each formula is a physical model formula or an approximate formula. For example, the formula (5) of the water supply enthalpy Ew is a quadratic approximation formula of the water supply temperature, and the formula (1) of the gas consumption G is a physical model formula.

ステップS11において、インタフェース部20から、推定パラメータとしてブロー係数Brを、設備種類として蒸気ボイラをそれぞれ取得するものとする。さらに、ステップS11では、取得部10と入力部21とから、蒸気量S、使用圧力P、給水温度Tw、定格熱出力Hrout、定格ガス消費量Groutを取得するものとする。この場合、ステップS121では、記憶部30の設備モデル式を検索することにより、式(1)、式(2)、式(3)、式(4)、式(5)、式(6)、式(7)を取得する。 In step S11, the blow coefficient Br is acquired as the estimation parameter and the steam boiler is acquired as the equipment type from the interface unit 20. Furthermore, in step S11, the amount of steam S, the working pressure P, the feed water temperature Tw, the rated heat output Hrout, and the rated gas consumption amount Grout are acquired from the acquisition unit 10 and the input unit 21. In this case, in step S121, by searching the equipment model formula in the storage unit 30, the formula (1), formula (2), formula (3), formula (4), formula (5), formula (6), Expression (7) is acquired.

設備モデル式決定部41は、ステップS11において推定パラメータとして取得されたブロー係数Br以外に推定が必要なパラメータを決める(S122)。ここでは、ステップS11において、ブロー係数Br以外のすべてのパラメータを取得できたため、推定パラメータはブロー係数Brと決定する。これに対し、もし蒸気量Sが取得部10から取得できない場合は、蒸気量Sも推定パラメータであると決定される。 The equipment model formula determination unit 41 determines parameters that need to be estimated in addition to the blow coefficient Br acquired as the estimation parameter in step S11 (S122). Here, in step S11, all parameters other than the blow coefficient Br have been acquired, so the estimated parameter is determined to be the blow coefficient Br. On the other hand, if the steam amount S cannot be acquired from the acquisition unit 10, the steam amount S is also determined to be an estimation parameter.

<適正範囲設定部42の処理> <Process of Appropriate Range Setting Unit 42>

図5は、図3中にステップS13で示された、適正範囲設定部42での処理の詳細を示すフローチャートである。 FIG. 5 is a flowchart showing details of the processing in the proper range setting unit 42 shown in step S13 in FIG.

適正範囲設定部42は、設備モデル式決定部41により決定された推定パラメータの適正範囲を記憶部30から検索し、検索された適正範囲を推定対象パラメータに設定する(S131)。すなわち、適正範囲設定部42は、設備種類と推定対象パラメータの名称とをタグとして、記憶部30が格納している適正範囲の中から検索し、各推定対象パラメータに対する適正範囲を設定する(S131)。 The appropriate range setting unit 42 searches the storage unit 30 for the appropriate range of the estimation parameter determined by the equipment model formula determination unit 41, and sets the retrieved appropriate range as the estimation target parameter (S131). That is, the appropriate range setting unit 42 searches the appropriate range stored in the storage unit 30 using the equipment type and the name of the estimation target parameter as tags, and sets the appropriate range for each estimation target parameter (S131). ).

適正範囲の検索には、例えばクラスタリングを用いることができるが、検索手法は問わない。適正範囲を見つけることのできない推定対象パラメータについては、適正範囲を設定しない。 For example, clustering can be used to search for the appropriate range, but the searching method does not matter. The appropriate range is not set for the estimation target parameter for which the appropriate range cannot be found.

図6〜図8に示すように、適正範囲には、様々なパターンがある。図6は、例えば、給水温度の適正範囲が「20℃〜40℃」と一つの連続した範囲の場合である。図7は、給水温度の適正範囲が「10℃〜20℃」と「30℃〜40℃」との、複数の連続した区間を有する場合である。 As shown in FIGS. 6 to 8, there are various patterns in the proper range. FIG. 6 shows a case where the proper range of the feed water temperature is, for example, one continuous range of “20° C. to 40° C.”. FIG. 7 shows a case where the proper range of the supply water temperature has a plurality of continuous sections of “10° C. to 20° C.” and “30° C. to 40° C.”.

図8は、ある一定時間(t1−t2)において温度変化をするパラメータの適正範囲を表す。点Aは、時刻t1における温度T1である。点B1は、時刻t2において取り得る最低温度T2である。点B2は、時刻t2において取り得る最高温度T3である。 FIG. 8 shows an appropriate range of parameters that change the temperature during a certain fixed time (t1-t2). Point A is the temperature T1 at time t1. The point B1 is the lowest temperature T2 that can be taken at the time t2. The point B2 is the maximum temperature T3 that can be taken at the time t2.

図8の例において、ある一定時間(t1とt2の差分)による温度変化範囲は、ΔT1からΔT2の間である。このように、時間差分(t1−t2)によるパラメータの変化量ΔTを適正範囲として設定することも可能である。 In the example of FIG. 8, the temperature change range for a certain fixed time (difference between t1 and t2) is between ΔT1 and ΔT2. In this way, it is also possible to set the parameter variation amount ΔT due to the time difference (t1−t2) as an appropriate range.

図5に戻る。適正範囲設定部42は、ステップS131で設定された適性範囲を用いて、尤度関数を設定する(S132)。尤度関数の設定方法を、蒸気ボイラを例にあげて詳細に説明する。 Returning to FIG. The appropriate range setting unit 42 sets the likelihood function using the appropriate range set in step S131 (S132). A method of setting the likelihood function will be described in detail by taking a steam boiler as an example.

設備モデル式決定部41により決定された設備モデル式が、式(1)、式(2)、式(3)、式(4)、式(5)、式(6)、式(7)であるとし、さらに、ガス消費量Gの計測値が記憶部30に格納されていると仮定する。 The equipment model formulas determined by the equipment model formula determination unit 41 are the formula (1), the formula (2), the formula (3), the formula (4), the formula (5), the formula (6), and the formula (7). In addition, it is assumed that the measured value of the gas consumption G is stored in the storage unit 30.

式(1)、式(2)、式(3)、式(4)、式(5)、式(6)、式(7)の各パラメータに値を入力することにより、ガス消費量Gの推定値(推定観測値)を算出することができる。算出されたガス流量Gの推定値の尤度は、尤度関数により計算できる。例えば、計測値と推定値との誤差である観測誤差が正規分布であるとき、尤度関数は以下の式(8)で求められる。式(8)において、「X」は確率関数、「μ」は平均、「σ」は分散である。ステップS132の説明として、図5にも式(8)を示す。 By inputting values to the parameters of the formula (1), the formula (2), the formula (3), the formula (4), the formula (5), the formula (6), and the formula (7), the gas consumption G An estimated value (estimated observed value) can be calculated. The likelihood of the calculated estimated value of the gas flow rate G can be calculated by a likelihood function. For example, when the observation error, which is the error between the measured value and the estimated value, has a normal distribution, the likelihood function is calculated by the following equation (8). In Expression (8), “X” is a probability function, “μ” is an average, and “σ 2 ”is a variance. As an explanation of step S132, equation (8) is also shown in FIG.

f(X:μ,σ)=(1/(2πσ1/2)exp(−(X−μ)/2σ
・・・式(8)
f(X:μ,σ 2 )=(1/(2πσ 2 ) ½ )exp(−(X−μ) 2 /2σ 2 ).
...Formula (8)

適正範囲設定部42は、ステップS131で設定されたパラメータの適正範囲を用いて、尤度関数を更新する。 The appropriate range setting unit 42 updates the likelihood function using the appropriate range of the parameters set in step S131.

例えば、給水温度が推定対象パラメータであり、ステップS131で適正範囲を「20℃〜40℃」と設定した場合、パラメータの尤度関数として一様分布を適用する。ここでは一様分布を用いるが、一様分布以外の分布も適用できる。給水温度の推定値が20℃から40℃の範囲内の場合、その尤度は「1/(40−20)」であり、それ以外の場合は尤度が「0」である。 For example, when the water supply temperature is the estimation target parameter and the appropriate range is set to “20° C. to 40° C.” in step S131, the uniform distribution is applied as the likelihood function of the parameter. A uniform distribution is used here, but a distribution other than the uniform distribution can also be applied. When the estimated value of the water supply temperature is in the range of 20°C to 40°C, the likelihood is "1/(40-20)", and in other cases, the likelihood is "0".

適正範囲設定部42は、パラメータの尤度関数を式(8)に示す尤度関数に乗算することにより、尤度関数を設定する。ここでは図6の適正範囲に基づく尤度関数の設定を説明したが、図7または図8に示した適正範囲に基づく尤度関数設定も、図6の場合の尤度関数設定を例として応用することが可能である。推定対象パラメータが複数であり、それぞれの適正範囲がある場合は、式(8)の尤度関数にそれぞれのパラメータの尤度関数を乗算することにより、尤度関数を設定すればよい。 The appropriate range setting unit 42 sets the likelihood function by multiplying the likelihood function shown in Expression (8) by the likelihood function of the parameter. Here, the setting of the likelihood function based on the appropriate range in FIG. 6 has been described, but the likelihood function setting based on the appropriate range shown in FIG. 7 or 8 is also applied by taking the likelihood function setting in the case of FIG. 6 as an example. It is possible to When there are a plurality of estimation target parameters and there are respective appropriate ranges, the likelihood function may be set by multiplying the likelihood function of equation (8) by the likelihood function of each parameter.

<パラメータ推定部43の処理> <Process of Parameter Estimating Unit 43>

パラメータ推定部43は、設備モデル式決定部41により決定された設備モデル式と適正範囲設定部42により設定された推定対象パラメータの適正範囲とに基づく尤度関数を用いることにより、パラメータを推定する。パラメータ推定手法としては、例えば粒子フィルタ法や、マルコフ連鎖モンテカルロ法を用いることができる。 The parameter estimation unit 43 estimates a parameter by using a likelihood function based on the equipment model formula determined by the equipment model formula determination unit 41 and the appropriate range of the estimation target parameter set by the appropriate range setting unit 42. . As the parameter estimation method, for example, a particle filter method or a Markov chain Monte Carlo method can be used.

<パラメータ判定部44の処理> <Process of Parameter Determination Unit 44>

パラメータ判定部44は、パラメータ推定部43により算出されたパラメータの推定値の良否を判定する。良否の判定基準としては、例えば、推定値の収束度、適正範囲、観測誤差などがある。 The parameter determination unit 44 determines the quality of the estimated value of the parameter calculated by the parameter estimation unit 43. Examples of the quality criterion include the degree of convergence of the estimated value, the appropriate range, and the observation error.

推定値の収束度には、推定ステップごとの推定値の分散を用いればよい。これに代えて、標準偏差などを用いてもよい。推定値の収束度に基づく良否判定では、予め設定した閾値を基準に判定される。適正範囲に基づく良否判定では、推定値が適正範囲内に収束した場合に「良」と判定し、収束しない場合に「否」と判定する。観測誤差に基づく良否判定では、予め設定された閾値を基準に判定される。推定値の収束度、適正範囲、観測誤差の複数の判定基準(評価基準)に基づく良否判定において、いずれか一つでも「否」と判定された場合には、パラメータ推定部43において再度推定処理を行う(S14)。 The variance of the estimated value for each estimation step may be used as the degree of convergence of the estimated value. Instead of this, a standard deviation or the like may be used. In the pass/fail judgment based on the degree of convergence of the estimated value, the judgment is made based on a preset threshold value. In the pass/fail determination based on the appropriate range, it is determined as “good” when the estimated value converges within the appropriate range, and is determined as “fail” when the estimated value does not converge. In the quality judgment based on the observation error, the judgment is made based on a preset threshold value. In the quality judgment based on a plurality of judgment criteria (evaluation criteria) of the convergence of the estimated value, the appropriate range, and the observation error, if any one of the judgments is “No”, the parameter estimation unit 43 performs the estimation processing again. Is performed (S14).

<出力部22の処理> <Processing of output unit 22>

出力部22は、入力部21、記憶部30、および演算部40からの情報を受信して、出力装置に表示する。出力部22は、演算部40のパラメータ判定部44による判定基準毎の判定結果と、推定値と、適正範囲とをそれぞれ出力装置に表示する。 The output unit 22 receives the information from the input unit 21, the storage unit 30, and the calculation unit 40, and displays it on the output device. The output unit 22 displays the determination result of each determination criterion by the parameter determination unit 44 of the calculation unit 40, the estimated value, and the appropriate range on the output device.

このように構成される本実施例によれば、推定対象パラメータを効率よく正確に推定することができる。本実施例によれば、観測値の種類の数よりも多いパラメータを推定する場合であっても、適正範囲内の推定値の尤度が高くなるように設定された尤度関数を用いることにより、各推定対象パラメータの推定値を適正範囲内に収束させることができ、真値に近い推定値を効率的に取得することができる。 According to this example configured as described above, the estimation target parameter can be estimated efficiently and accurately. According to the present embodiment, even when estimating more parameters than the number of types of observed values, by using the likelihood function set so that the likelihood of estimated values within the appropriate range is high. The estimated values of the estimation target parameters can be converged within the proper range, and the estimated value close to the true value can be efficiently acquired.

なお、本実施例では、適正範囲設定部42において様々なパターンの適正範囲を用いることにより、推定対象パラメータについてより多くの情報を得ることができる。したがって、この場合は、より一層、真値に近い推定値を効率的に取得することができる。 In addition, in the present embodiment, by using the appropriate ranges of various patterns in the appropriate range setting unit 42, it is possible to obtain more information about the estimation target parameter. Therefore, in this case, the estimated value closer to the true value can be acquired more efficiently.

図9を用いて第2実施例を説明する。本実施例は、第1実施例の変形例であるため、第1実施例との相違を中心に説明する。 A second embodiment will be described with reference to FIG. Since this embodiment is a modification of the first embodiment, the description will focus on the differences from the first embodiment.

第1実施例では、適正範囲設定部42において尤度関数を設定する場合、すべての推定対象パラメータの尤度関数を乗算することにより、全体としての尤度関数を設定する場合を説明した。しかし、もしパラメータ判定部44での判定結果が「否」であった場合、パラメータ推定部43で再度推定することになるが、再推定してもその判定結果が「否」となってしまい、推定処理(S14)と判定処理(S15)とが繰り返される可能性が考えられる。 In the first embodiment, when setting the likelihood function in the appropriate range setting unit 42, the case where the likelihood function as a whole is set by multiplying the likelihood functions of all estimation target parameters has been described. However, if the determination result by the parameter determination unit 44 is “No”, the parameter estimation unit 43 will re-estimate the result, but even if it is re-estimated, the determination result will be “No”, The estimation process (S14) and the determination process (S15) may be repeated.

図9は、本実施例における適正範囲設定部42の処理S13Aを示す。図9に示すフローチャートと図5で述べたフローチャートとを比較すると、本実施例ではステップS132Aが第1実施例のステップS132と異なる。そこで、この異なるステップS132Aについて説明する。 FIG. 9 shows processing S13A of the appropriate range setting unit 42 in this embodiment. Comparing the flowchart shown in FIG. 9 with the flowchart described in FIG. 5, step S132A in this embodiment is different from step S132 in the first embodiment. Therefore, the different step S132A will be described.

適正範囲設定部42は、ステップS131で設定された適正範囲に基づいて尤度関数を設定する場合に、各推定対象パラメータの適正範囲の組合せを計算して、尤度関数を設定する(S132A)。 When setting the likelihood function based on the appropriate range set in step S131, the appropriate range setting unit 42 calculates the combination of the appropriate ranges of the estimation target parameters and sets the likelihood function (S132A). ..

例えば、推定対象のパラメータがパラメータPaとパラメータPbとの2つある場合、尤度設定の組合せは、パラメータPaの適正範囲による尤度設定、パラメータPbの適正範囲による尤度設定、パラメータPaの適正範囲とパラメータPbの適正範囲とによる尤度設定の、合計3つある。 For example, when there are two parameters to be estimated, that is, the parameter Pa and the parameter Pb, the combinations of the likelihood settings include the likelihood setting based on the appropriate range of the parameter Pa, the likelihood setting based on the appropriate range of the parameter Pb, and the appropriateness of the parameter Pa. There are a total of three likelihood settings based on the range and the appropriate range of the parameter Pb.

そして、本実施例のパラメータ推定部43は、それぞれの組合せの尤度設定を用いて、パラメータ推定処理を行う。本実施例のパラメータ判定部44は、それぞれの組合せの尤度設定によるパラメータ推定値の良否を判定する。 Then, the parameter estimation unit 43 of the present embodiment performs the parameter estimation process using the likelihood setting of each combination. The parameter determination unit 44 of this embodiment determines the quality of the parameter estimation value based on the likelihood setting of each combination.

パラメータ判定部44は、第1実施例で述べた判定基準を用いることができる。その判定基準の中で、適正範囲の設定されていない推定対象パラメータに対しても、ステップS131で設定された適正範囲を適用して、良否を判定する。 The parameter determination unit 44 can use the determination criteria described in the first embodiment. Among the criteria, the appropriate range set in step S131 is applied to the estimation target parameter for which the appropriate range is not set, and the quality is determined.

例えば、パラメータPaの適正範囲だけに基づいて尤度関数を設定する組合せの場合、パラメータ判定部44は、ステップS131で設定されたパラメータPbの適正範囲に基づいて、パラメータPbの推定値の良否を判定する。 For example, in the case of the combination in which the likelihood function is set only based on the proper range of the parameter Pa, the parameter determination unit 44 determines whether the estimated value of the parameter Pb is good or bad based on the proper range of the parameter Pb set at step S131. judge.

このように構成される本実施例も第1実施例と同様の作用効果を奏する。さらに、本実施例では、複数の尤度設定の組合せを用いて推定するため、パラメータ判定部44で「良」と判定されるパラメータ推定値を取得する確率を高めることができる。 This embodiment, which is configured in this way, also exhibits the same effects as the first embodiment. Furthermore, in the present embodiment, since the estimation is performed using a combination of a plurality of likelihood settings, it is possible to increase the probability of acquiring the parameter estimation value that is determined as “good” by the parameter determination unit 44.

なお、本発明は、上述した実施形態に限定されない。当業者であれば、本発明の範囲内で、種々の追加や変更等を行うことができる。上述の実施形態において、添付図面に図示した構成例に限定されない。本発明の目的を達成する範囲内で、実施形態の構成や処理方法は適宜変更することが可能である。 The present invention is not limited to the above embodiment. Those skilled in the art can make various additions and changes within the scope of the present invention. The above embodiment is not limited to the configuration example illustrated in the accompanying drawings. The configuration and the processing method of the embodiment can be appropriately changed within the scope of achieving the object of the present invention.

また、本発明の各構成要素は、任意に取捨選択することができ、取捨選択した構成を具備する発明も本発明に含まれる。さらに特許請求の範囲に記載された構成は、特許請求の範囲で明示している組合せ以外にも組み合わせることができる。 Further, each constituent element of the present invention can be arbitrarily selected, and an invention having a selected configuration is also included in the present invention. Furthermore, the configurations described in the claims can be combined in addition to the combinations specified in the claims.

1:パラメータ推定装置、2:センサ、3:設備、41:設備モデル式決定部、42:適正範囲設定部、41:パラメータ推定部、42:パラメータ判定部 1: Parameter estimation device, 2: Sensor, 3: Equipment, 41: Equipment model formula determination unit, 42: Appropriate range setting unit, 41: Parameter estimation unit, 42: Parameter determination unit

Claims (12)

設備に関するパラメータを推定する装置であって、
前記設備に関する物理量を測定するセンサから観測値を取得し、
前記観測値を対象パラメータに基づいて算出する設備モデル式を取得し、
前記設備モデル式に前記対象パラメータの推定値を代入して得られる推定観測値の前記観測値に対する尤度を第1尤度関数により算出し、
前記対象パラメータの推定値があらかじめ用意された所定範囲内に属する確率を求める第2尤度関数と前記第1尤度関数とに基づいて、前記対象パラメータの推定値の尤度を算出し、
前記算出された尤度に基づいて前記対象パラメータを推定し、
前記推定された対象パラメータを所定の評価基準にしたがって評価する、
パラメータ推定装置。
A device for estimating parameters relating to equipment,
Obtaining an observation value from a sensor that measures a physical quantity related to the equipment,
Acquire the equipment model formula to calculate the observed value based on the target parameter,
The likelihood for the observed value of the estimated observed value obtained by substituting the estimated value of the target parameter into the equipment model formula is calculated by the first likelihood function,
The likelihood of the estimated value of the target parameter is calculated based on the second likelihood function and the first likelihood function for obtaining the probability that the estimated value of the target parameter belongs to a predetermined range prepared in advance,
Estimating the target parameter based on the calculated likelihood,
Evaluating the estimated target parameter according to a predetermined evaluation criterion,
Parameter estimation device.
前記対象パラメータの評価結果を出力する、
パラメータ推定装置。
Output the evaluation result of the target parameter,
Parameter estimation device.
前記対象パラメータの推定値の前記所定範囲における分布密度が前記所定範囲外での分布密度よりも大きくなるように前記対象パラメータの尤度が設定されている、
請求項1に記載のパラメータ推定装置。
The likelihood of the target parameter is set such that the distribution density of the estimated value of the target parameter in the predetermined range is larger than the distribution density outside the predetermined range,
The parameter estimation device according to claim 1.
前記設備モデル式は、設備の種類に応じてあらかじめ用意されている、
請求項1に記載のパラメータ推定装置。
The equipment model formula is prepared in advance according to the type of equipment,
The parameter estimation device according to claim 1.
前記センサから得られる観測値の数よりも前記対象パラメータの数が大きく設定されている、
請求項1に記載のパラメータ推定装置。
The number of the target parameter is set to be larger than the number of observation values obtained from the sensor,
The parameter estimation device according to claim 1.
前記設備モデル式には、異なる複数の対象パラメータが含まれている、
請求項1に記載のパラメータ推定装置。
The equipment model formula includes a plurality of different target parameters,
The parameter estimation device according to claim 1.
前記所定範囲は、一つの連続した範囲である、
請求項1に記載のパラメータ推定装置。
The predetermined range is one continuous range,
The parameter estimation device according to claim 1.
前記所定範囲は、複数の連続した範囲からなる、
請求項1に記載のパラメータ推定装置。
The predetermined range includes a plurality of continuous ranges,
The parameter estimation device according to claim 1.
前記所定範囲は、あらかじめ設定される単位時間で変化する、
請求項1に記載のパラメータ推定装置。
The predetermined range changes in a preset unit time,
The parameter estimation device according to claim 1.
前記所定範囲は、複数の対象パラメータそれぞれに関する前記所定範囲の組合せにしたがって前記第2尤度関数と前記第1尤度関数とに基づいて、前記複数の対象パラメータの推定値の尤度をそれぞれ算出する、
請求項1に記載のパラメータ推定装置。
In the predetermined range, the likelihoods of the estimated values of the plurality of target parameters are calculated based on the second likelihood function and the first likelihood function according to a combination of the predetermined ranges for each of the plurality of target parameters. To do
The parameter estimation device according to claim 1.
前記推定された複数の対象パラメータを所定の評価基準にしたがって評価する場合、前記所定範囲の組合せに含まれていない所定範囲を持つ対象パラメータに対しても当該対象パラメータに関する前記所定範囲を設定して評価する、
請求項10に記載のパラメータ推定装置。
In the case of evaluating the estimated plurality of target parameters according to a predetermined evaluation criterion, the predetermined range regarding the target parameter is set for a target parameter having a predetermined range not included in the combination of the predetermined ranges. evaluate,
The parameter estimation device according to claim 10.
設備に関するパラメータを計算機により推定する方法であって、
前記計算機は、
前記設備に関する物理量を測定するセンサから観測値を取得し、
前記観測値を対象パラメータに基づいて算出する設備モデル式を取得し、
前記設備モデル式に前記対象パラメータの推定値を代入して得られる推定観測値の前記観測値に対する尤度を第1尤度関数により算出し、
前記対象パラメータの推定値があらかじめ用意された所定範囲内に属する確率を求める第2尤度関数と前記第1尤度関数とに基づいて、前記対象パラメータの推定値の尤度を算出し、
前記算出された尤度に基づいて前記対象パラメータを推定し、
前記推定された対象パラメータを所定の評価基準にしたがって評価する、
パラメータ推定方法。
A method for estimating parameters relating to equipment by a computer,
The calculator is
Obtaining an observation value from a sensor that measures a physical quantity related to the equipment,
Acquire the equipment model formula to calculate the observed value based on the target parameter,
The likelihood for the observed value of the estimated observed value obtained by substituting the estimated value of the target parameter into the equipment model formula is calculated by the first likelihood function,
The likelihood of the estimated value of the target parameter is calculated based on the second likelihood function and the first likelihood function for obtaining the probability that the estimated value of the target parameter belongs to a predetermined range prepared in advance,
Estimating the target parameter based on the calculated likelihood,
Evaluating the estimated target parameter according to a predetermined evaluation criterion,
Parameter estimation method.
JP2018232579A 2018-12-12 2018-12-12 Parameter estimation device and method Pending JP2020095451A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018232579A JP2020095451A (en) 2018-12-12 2018-12-12 Parameter estimation device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018232579A JP2020095451A (en) 2018-12-12 2018-12-12 Parameter estimation device and method

Publications (1)

Publication Number Publication Date
JP2020095451A true JP2020095451A (en) 2020-06-18

Family

ID=71086418

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018232579A Pending JP2020095451A (en) 2018-12-12 2018-12-12 Parameter estimation device and method

Country Status (1)

Country Link
JP (1) JP2020095451A (en)

Similar Documents

Publication Publication Date Title
US10261485B2 (en) Systems and methods for detecting changes in energy usage in a building
JP2016177682A (en) Facility evaluation device, facility evaluation method and computer program
US6996508B1 (en) System and method for remote retrofit identification of energy consumption systems and components
JP7006766B2 (en) Piping diagnostic method, piping diagnostic device, and piping diagnostic system
CN109840671B (en) Operation improvement effect calculation device, operation improvement effect calculation method, and recording medium
Zhang et al. Multivariate and multimodal wind distribution model based on kernel density estimation
JP6862130B2 (en) Anomaly detection device, anomaly detection method, and program
JP6181079B2 (en) Computational fluid dynamics system and method of use
KR20180137513A (en) Monitoring device, monitoring method and program of target device
Monari et al. Characterization of an airflow network model by sensitivity analysis: parameter screening, fixing, prioritizing and mapping
JP2017010232A (en) Plant diagnostic device and plant diagnostic method
US10817630B2 (en) Apparatus and method for analyzing buildings
CN117664281B (en) Ultrasonic water meter fault detection and automatic calibration method and system based on Internet of Things
JP2021076309A (en) Estimation method, estimation device, program, and generation method of learned model
CN106575312B (en) Method for detecting anomalies in a distribution network, in particular a water distribution network
US20160275407A1 (en) Diagnostic device, estimation method, non-transitory computer readable medium, and diagnostic system
WO2016136391A1 (en) Fault point locating device and method, electric power system monitoring system, and facility planning support system
CN111948534B (en) Generator state early warning method and system
KR20140109596A (en) Apparatus method for predicting and correcting meter reading data in a range not reading
JP2020095451A (en) Parameter estimation device and method
TW202042079A (en) Water quality monitoring system and control method capable of effectively reducing the manpower of operating equipment and achieving automatic operation
JP2020035107A (en) Plant data analysis system
CN111443602A (en) Hidden Markov-based valve health degree estimation and life prediction method
Remmen et al. Refinement of Dynamic Non-Residential Building Archetypes Using Measurement Data and Bayesian Calibration
US20180276322A1 (en) Failure diagnosis apparatus, monitoring apparatus, failure diagnosis method and recording medium