JP5872441B2 - State estimation device and state estimation method - Google Patents

State estimation device and state estimation method Download PDF

Info

Publication number
JP5872441B2
JP5872441B2 JP2012256239A JP2012256239A JP5872441B2 JP 5872441 B2 JP5872441 B2 JP 5872441B2 JP 2012256239 A JP2012256239 A JP 2012256239A JP 2012256239 A JP2012256239 A JP 2012256239A JP 5872441 B2 JP5872441 B2 JP 5872441B2
Authority
JP
Japan
Prior art keywords
operation mode
harmonic
state estimation
power consumption
database
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.)
Active
Application number
JP2012256239A
Other languages
Japanese (ja)
Other versions
JP2014103830A (en
Inventor
保博 田口
保博 田口
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
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 JP2012256239A priority Critical patent/JP5872441B2/en
Publication of JP2014103830A publication Critical patent/JP2014103830A/en
Application granted granted Critical
Publication of JP5872441B2 publication Critical patent/JP5872441B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225Demand response systems, e.g. load shedding, peak shaving
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20End-user application control systems
    • Y04S20/222Demand response systems, e.g. load shedding, peak shaving

Landscapes

  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Description

本発明の実施形態は、電力を消費する機器の状態を推定する、ロードサーベイ技術に関する。   Embodiments described herein relate generally to a load survey technique for estimating a state of a device that consumes power.

省エネルギー(省エネ)意識の高まりや災害時におけるエネルギー自給のニーズなどを背景に、エネルギー管理システム(Energy Management System:EMS)が普及しつつある。EMSは電力の需要家(ビル、工場、家庭、店舗、事務所など)の有する需要家機器(以下、電気機器と総称する)の運転スケジュールを作成したり、デマンドレスポンスに対応したりする機能を備える。EMSによれば省エネはもとより、省コスト、CO2排出量の削減、再生可能エネルギーの導入などの多種多様なメリットがもたらされる。   Energy management systems (EMS) are becoming widespread against the background of increasing awareness of energy saving (energy saving) and the need for energy self-sufficiency in the event of a disaster. EMS has functions to create operation schedules for consumer equipment (hereinafter collectively referred to as electrical equipment) possessed by electricity consumers (buildings, factories, homes, stores, offices, etc.) and to respond to demand response. Prepare. According to EMS, not only energy saving but also various advantages such as cost saving, reduction of CO2 emission, introduction of renewable energy, etc. are brought about.

複数の電気機器を最適に運転させるためにはそれぞれの電気機器の現在の状態を知ることが必要である。しかし、電気機器に状態を通知する機能を持たせたり、状態を検知するセンサを個別に設けたりすることはコスト面などから得策でない。多数の電気機器を備えるビルや工場ではなおさらである。そこで、需要家内の配電系統に接続したセンサにより得た波形を高調波解析して稼働中の電気機器を推定する技術が提案されている。   In order to optimally operate a plurality of electrical devices, it is necessary to know the current state of each electrical device. However, it is not advantageous from the viewpoint of cost or the like to provide an electric device with a function of notifying the state or providing a sensor for detecting the state individually. This is especially true in buildings and factories with many electrical devices. In view of this, a technique has been proposed in which a waveform obtained by a sensor connected to a distribution system in a consumer is subjected to harmonic analysis to estimate an operating electrical device.

特許第4802129号公報Japanese Patent No. 4802129 特開2006−17456号公報JP 2006-17456 A 特許第3877269号公報Japanese Patent No. 3877269

提案されたいずれの技術も、稼働中の電気機器を推定すること、つまり電気機器のオン/オフを推定することを狙いとするものである。しかしながら近年の電気機器には単純なオン/オフだけでなく複数の運転モードを持つものがあり、これらの機器が混在すると既存の手法では対応が困難である。また、デマンドレスポンスに応じてオフする機器を決めようにも稼働中の電気機器の運転モードが分からなければ決めようがない。工場などで盲目的に電気機器を停止してしまうと膨大な損害がもたらされるケースもある。電気機器のオン/オフだけでなく、運転モードをも簡易に推定可能とする技術が要望されている。   All of the proposed technologies are aimed at estimating an electric device in operation, that is, estimating on / off of the electric device. However, some recent electrical devices have not only simple on / off but also a plurality of operation modes. When these devices are mixed, it is difficult to cope with existing methods. In addition, it is impossible to determine the device to be turned off according to the demand response unless the operation mode of the operating electric device is known. In some cases, blindly shutting down electrical equipment in factories can cause enormous damage. There is a demand for a technique that can easily estimate not only the on / off state of electrical equipment but also the operation mode.

目的は、電気機器の運転モードを簡易に推定可能な状態推定装置および状態推定方法を提供することにある。   An object is to provide a state estimation device and a state estimation method capable of easily estimating an operation mode of an electric device.

実施形態によれば、状態推定装置は、計算部と、データベースと、推定部とを具備する。計算部は、電気機器の接続される給電系統に生じる電力の波形を高調波解析して、有効電力高調波の次数ごとの高調波成分を算出する。データベースは、高調波の次数ごとの高調波含有率を電気機器の運転モードごとに記憶する。推定部は、算出された次数ごとの高調波成分から、データベースに記憶される高調波含有率に基づいて運転モードごとの消費電力を逆算し、その結果に基づいて電気機器の運転モードを推定する。   According to the embodiment, the state estimation device includes a calculation unit, a database, and an estimation unit. The calculation unit performs harmonic analysis on the waveform of the power generated in the power feeding system to which the electrical device is connected, and calculates a harmonic component for each order of the active power harmonic. A database memorize | stores the harmonic content rate for every order of a harmonic for every operation mode of an electric equipment. The estimation unit back-calculates the power consumption for each operation mode based on the harmonic content stored in the database from the calculated harmonic components for each order, and estimates the operation mode of the electric device based on the result. .

図1は、実施形態に係る状態推定装置を適用可能なエネルギー管理システムの一例を示す図である。FIG. 1 is a diagram illustrating an example of an energy management system to which the state estimation device according to the embodiment can be applied. 図2は、図1に示される状態推定装置10の一例を示す機能ブロック図である。FIG. 2 is a functional block diagram showing an example of the state estimation apparatus 10 shown in FIG. 図3は、測定装置20により取得される測定量の波形の一例を示す図である。FIG. 3 is a diagram illustrating an example of a waveform of a measurement amount acquired by the measurement apparatus 20. 図4は、計算部11の一例を示す機能ブロック図である。FIG. 4 is a functional block diagram illustrating an example of the calculation unit 11. 図5は、FFT処理部111から出力される高調波成分計算結果(実測高調波成分結果)の一例を示すグラフである。FIG. 5 is a graph showing an example of harmonic component calculation results (measured harmonic component results) output from the FFT processing unit 111. 図6は、運転モードM1における運転モード高調波成分の含有率の分布を示す図である。FIG. 6 is a diagram showing the distribution of the content ratio of the operation mode harmonic component in the operation mode M1. 図7は、運転モードM2における運転モード高調波成分の含有率の分布を示す図である。FIG. 7 is a diagram showing the distribution of the content rate of the operation mode harmonic component in the operation mode M2. 図8は、運転モードM3における運転モード高調波成分の含有率の分布を示す図である。FIG. 8 is a diagram showing the distribution of the content rate of the operation mode harmonic component in the operation mode M3. 図9は、運転モードM4における運転モード高調波成分の含有率の分布を示す図である。FIG. 9 is a diagram showing the distribution of the content rate of the operation mode harmonic component in the operation mode M4. 図10は、電気機器の運転モードの遷移の一例を示す図である。FIG. 10 is a diagram illustrating an example of the transition of the operation mode of the electric device. 図11は、推定部13における処理の一例を模式的に示す図である。FIG. 11 is a diagram schematically illustrating an example of processing in the estimation unit 13. 図12は、推定部13による処理を説明するための模式図である。FIG. 12 is a schematic diagram for explaining processing by the estimation unit 13. 図13は、推定部13における処理の他の例を模式的に示す図である。FIG. 13 is a diagram schematically illustrating another example of processing in the estimation unit 13. 図14は、推定部13における処理の他の例を模式的に示す図である。FIG. 14 is a diagram schematically illustrating another example of processing in the estimation unit 13. 図15は、推定部13における処理の他の例を模式的に示す図である。FIG. 15 is a diagram schematically illustrating another example of processing in the estimation unit 13. 図16は、実施形態に係る状態推定装置を適用可能なエネルギー管理システムの他の例を示す図である。FIG. 16 is a diagram illustrating another example of the energy management system to which the state estimation device according to the embodiment is applicable. 図17は、一つの装置の運転モードの推移の一例を示す模式図である。FIG. 17 is a schematic diagram illustrating an example of the transition of the operation mode of one device. 図18は、複数の装置の運転モードの推移の一例を示す模式図である。FIG. 18 is a schematic diagram illustrating an example of transition of operation modes of a plurality of apparatuses. 図19は、運転を急に停止した場合の損害を運転モードごとに示す模式図である。FIG. 19 is a schematic diagram showing damage when the operation is suddenly stopped for each operation mode. 図20は、全ての装置が停止した時の損害Pra、稼働台数n、および消費電力Pmをモードごとに配列して示す図である。FIG. 20 is a diagram showing the damage Pra, the number n of operations, and the power consumption Pm when all devices are stopped arranged for each mode. 図21は、実施形態に係る状態推定装置の他の例を示す機能ブロック図である。FIG. 21 is a functional block diagram illustrating another example of the state estimation device according to the embodiment. 図22は、結果表示部71により表示装置に表示される運転モード消費電力の一例を示す図である。FIG. 22 is a diagram illustrating an example of operation mode power consumption displayed on the display device by the result display unit 71. 図23は、各モードの継続時間を装置ごとに示すグラフの一例を示す図である。FIG. 23 is a diagram illustrating an example of a graph indicating the duration of each mode for each device.

図1は、実施形態に係る状態推定装置を適用可能なエネルギー管理システムの一例を示す図である。実施形態では需要家の一例としてビルを採りあげる。このほか病院、研究施設、学校、工場なども需要家になり得る。図1において、電力系統(power grid)からビル内の母線に引き込まれた電力は、母線からフィーダを介して例えば各階のローカル系統(local grid)40に分配される。   FIG. 1 is a diagram illustrating an example of an energy management system to which the state estimation device according to the embodiment can be applied. In the embodiment, a building is taken as an example of a consumer. In addition, hospitals, research facilities, schools, factories, etc. can be customers. In FIG. 1, power drawn from a power grid to a bus in the building is distributed from the bus to a local grid 40 on each floor via a feeder, for example.

電気機器30−1〜30−5はコンセントなどを介してローカル系統40に接続され、それぞれローカル系統40から電力の供給を受けて機能する。電気機器30−1〜30−5は、例えば、空調機(エアコン)、ノートPC、デスクトップPC、コピー/プリンタ、及び照明器具などである。工場などにおいては半導体製造機器なども電気機器に含まれる。   The electric devices 30-1 to 30-5 are connected to the local system 40 through an outlet or the like, and function by receiving power supplied from the local system 40, respectively. The electric devices 30-1 to 30-5 are, for example, an air conditioner (air conditioner), a notebook PC, a desktop PC, a copy / printer, and a lighting fixture. In factories and the like, semiconductor manufacturing equipment is also included in electrical equipment.

母線からローカル系統40に至るフィーダには測定装置20が接続される。測定装置20は、フィーダからローカル系統40に供給される電流、電圧、電力などの電気的物理量の値(以下、測定量と総称する)を取得する。取得された測定量は状態推定装置10に与えられる。   A measuring device 20 is connected to a feeder from the bus to the local system 40. The measuring device 20 acquires values of electrical physical quantities such as current, voltage and power supplied to the local system 40 from the feeder (hereinafter collectively referred to as measured quantities). The acquired measurement amount is given to the state estimation device 10.

状態推定装置10は測定量に基づいて、電気機器30−1〜30−5の状態や消費される電力負荷、有効電力の値などを推定する。実施形態では特に、状態推定装置10は、電気機器30−1〜30−5のそれぞれの運転モードと、当該運転モードにおける消費電力(以下、運転モード消費電力と称する)とを推定するとする。   The state estimation device 10 estimates the state of the electrical devices 30-1 to 30-5, the consumed power load, the value of active power, and the like based on the measured amount. In the embodiment, in particular, the state estimation device 10 estimates each operation mode of the electric devices 30-1 to 30-5 and power consumption in the operation mode (hereinafter referred to as operation mode power consumption).

状態推定装置10による推定の結果は、エネルギー管理サーバ50に渡される。エネルギー管理サーバ50は、状態推定装置10から取得した情報に基づいて、需要家の消費エネルギーのピークカット、ピークシフト、あるいはデマンドレスポンスなどの制御を実施する。   The result of estimation by the state estimation device 10 is passed to the energy management server 50. Based on the information acquired from the state estimation device 10, the energy management server 50 performs control such as peak cut, peak shift, or demand response of consumer energy consumption.

なお図1においてはビルの1階に至るフィーダに測定装置20を接続する例を示すが、この限りではない。例えば電力系統からの分岐線や母線に測定装置20を接続してビル全体での測定量を得るようにしても良い。また、状態推定装置10およびエネルギー管理サーバ50がビル監視室に設置される例を示すが、これも一例である。次に、複数の実施形態につき説明する。   In addition, although the example which connects the measuring apparatus 20 to the feeder which reaches the 1st floor of a building is shown in FIG. 1, it is not this limitation. For example, the measuring device 20 may be connected to a branch line or a bus line from the electric power system to obtain a measurement amount in the entire building. Moreover, although the example in which the state estimation apparatus 10 and the energy management server 50 are installed in a building monitoring room is shown, this is also an example. Next, a plurality of embodiments will be described.

[第1の実施形態]
図2は、図1に示される状態推定装置10の一例を示す機能ブロック図である。状態推定装置10は、例えば、CPU(Central Processing Unit)およびメモリを備えるハードウェアとしてのコンピュータと、CPUにより実行される命令を含むソフトウェアとしてのプログラムとが協働して実現される。
[First Embodiment]
FIG. 2 is a functional block diagram showing an example of the state estimation apparatus 10 shown in FIG. The state estimation device 10 is realized by, for example, a computer as hardware including a CPU (Central Processing Unit) and a memory and a program as software including instructions executed by the CPU.

状態推定装置10は、計算部11、保存部12、推定部13、およびデータベース14を備える。このうちデータベース14は、電気機器の運転モードごとの有効電力高調波成分を記憶する。つまりデータベース14は、高調波の次数ごとの高調波含有率(後述する)を電気機器の運転モードごとに記憶する。予めデータベース14に記憶された高調波成分と、計算部11により算出された高調波成分計算結果とを比較することで、電気機器の運転モードおよび運転モード消費電力を推定することが可能である。   The state estimation device 10 includes a calculation unit 11, a storage unit 12, an estimation unit 13, and a database 14. Among these, the database 14 memorize | stores the active power harmonic component for every operation mode of an electric equipment. That is, the database 14 stores the harmonic content (described later) for each harmonic order for each operation mode of the electric device. By comparing the harmonic component stored in the database 14 in advance with the harmonic component calculation result calculated by the calculation unit 11, it is possible to estimate the operation mode and the operation mode power consumption of the electric device.

計算部11は、測定装置20から与えられる測定量の波形を分析し、有効電力の高調波成分を計算する。つまり計算部11は、ローカル系統40に生じる電力の波形を高調波解析して、有効電力高調波の次数ごとの高調波成分を算出する。算出された高調波成分の計算結果(高調波成分計算結果)は、保存部12に保存されるとともに推定部13に渡される。   The calculation unit 11 analyzes the waveform of the measurement amount given from the measurement device 20 and calculates the harmonic component of the active power. That is, the calculation unit 11 performs harmonic analysis on the waveform of power generated in the local system 40 and calculates a harmonic component for each order of the active power harmonic. The calculated harmonic component calculation result (harmonic component calculation result) is stored in the storage unit 12 and passed to the estimation unit 13.

推定部13は、高調波成分計算結果に基づいてデータベース14を参照し、電気機器のそれぞれの運転モードと、運転モード消費電力とを推定する。つまり推定部13は、計算部11により算出された次数ごとの高調波成分から、データベース14に記憶される高調波含有率に基づいて運転モードごとの消費電力を逆算する。そしてその結果に基づいて、推定部13は電気機器の運転モードを推定する。   The estimation unit 13 refers to the database 14 based on the harmonic component calculation result, and estimates each operation mode and operation mode power consumption of the electrical device. That is, the estimation unit 13 back-calculates the power consumption for each operation mode based on the harmonic content stored in the database 14 from the harmonic components for each order calculated by the calculation unit 11. Based on the result, the estimation unit 13 estimates the operation mode of the electric device.

図3は、測定装置20により取得される測定量の波形の一例を示す図である。図3のグラフにおいて横軸を時間とし、縦軸を測定値とする。電圧V、電流I、および電力Pの各波形が示される。   FIG. 3 is a diagram illustrating an example of a waveform of a measurement amount acquired by the measurement apparatus 20. In the graph of FIG. 3, the horizontal axis is time, and the vertical axis is a measured value. Each waveform of voltage V, current I, and power P is shown.

電力の瞬時値をp、電圧の瞬時値をv、電流の瞬時値をiとすると式(1)の関係が成り立つ。

Figure 0005872441
When the instantaneous value of power is p, the instantaneous value of voltage is v, and the instantaneous value of current is i, the relationship of Equation (1) holds.
Figure 0005872441

式(1)から分かるように、測定量として電圧の瞬時値vと電流の瞬時値iとを得ることができれば、電力の瞬時値pを計算で求めることができる。   As can be seen from Equation (1), if the instantaneous voltage value v and the instantaneous current value i can be obtained as measured quantities, the instantaneous power value p can be obtained by calculation.

図4は、計算部11の一例を示す機能ブロック図である。計算部11はその処理機能としてFFT処理部111を備える。FFT処理部111は測定装置20から取得した測定量を高速フーリエ変換(Fast Fourier Transformation:FFT)し、測定量の波形に含まれる高調波成分を次数ごとに算出する。算出された高調波成分(高調波成分計算結果)は保存部12に記憶される。FFT処理は電力の瞬時値p、電圧の瞬時値v、電流の瞬時値iのいずれについても実施することが可能である。   FIG. 4 is a functional block diagram illustrating an example of the calculation unit 11. The calculation unit 11 includes an FFT processing unit 111 as its processing function. The FFT processing unit 111 performs fast Fourier transformation (FFT) on the measurement amount acquired from the measurement apparatus 20 and calculates a harmonic component included in the waveform of the measurement amount for each order. The calculated harmonic component (harmonic component calculation result) is stored in the storage unit 12. The FFT processing can be performed for any of the instantaneous value p of power, the instantaneous value v of voltage, and the instantaneous value i of current.

多くの電力機器では、運転時の電圧と電流との位相がほぼ等しいと看做せるので、各高調波成分の力率を力率=1と近似することが可能である。これに対し運転時の電圧と電流との位相差の大きい電気機器では次式(1)′のように高調波成分を補正する。
有効電力Pの高調波成分 = 電力Pの高調波成分×高調波成分の力率
= 電力Pの高調波成分×cos(θ) ・・・(1)′
θ:高調波成分の電圧と電流との位相差
力率:cos(θ)
図5は、FFT処理部111から出力される高調波成分計算結果(実測高調波成分結果)の一例を示すグラフである。図5において、横軸は高調波の次数を示し、縦軸は次数ごとの高調波成分の含有率(高調波含有率)を示す。図5には有効電力Pの高調波成分と電圧Vの高調波成分との双方が示される。高調波含有率は高調波の次数ごとに計算することができ、次式(2)で定義される。

Figure 0005872441
In many electric power devices, it can be considered that the phases of voltage and current during operation are almost equal, so that the power factor of each harmonic component can be approximated to power factor = 1. On the other hand, harmonic components are corrected as shown in the following equation (1) ′ in an electric device having a large phase difference between the voltage and current during operation.
Harmonic component of active power P = Harmonic component of power P x Power factor of harmonic component
= Harmonic component of power P × cos (θ) (1) ′
θ: Phase difference between voltage and current of harmonic components
Power factor: cos (θ)
FIG. 5 is a graph showing an example of harmonic component calculation results (measured harmonic component results) output from the FFT processing unit 111. In FIG. 5, the horizontal axis indicates the order of harmonics, and the vertical axis indicates the content of harmonic components (harmonic content) for each order. FIG. 5 shows both the harmonic component of the active power P and the harmonic component of the voltage V. The harmonic content rate can be calculated for each harmonic order, and is defined by the following equation (2).
Figure 0005872441

次に、図2に示されるデータベース14に記憶される内容につき説明する。図6〜図9に、データベース14に記憶される内容の一例を示す。これらは電気機器の一つの運転モードにおける高調波成分の次数ごとの含有率を示すグラフであり、横軸が次数を示し、縦軸が次数ごとの高調波含有率を示す。或る電気機器の運転モードをM1〜MNと表記して区別することにすると、図6〜図9は、それぞれ運転モードM1,M2,M3,M4における高調波含有率の分布を示す。運転モードごとの高調波含有率を総称して高調波含有率パターンと称する。   Next, the contents stored in the database 14 shown in FIG. 2 will be described. 6 to 9 show examples of contents stored in the database 14. These are graphs showing the content of harmonic components for each order in one operation mode of the electric device, with the horizontal axis showing the order and the vertical axis showing the harmonic content for each order. If the operation modes of a certain electrical device are distinguished from each other by expressing them as M1 to MN, FIGS. 6 to 9 show distributions of harmonic content in the operation modes M1, M2, M3, and M4, respectively. The harmonic content for each operation mode is collectively referred to as a harmonic content pattern.

図10は、電気機器の運転モードの遷移の一例を示す図である。或る電気機器は、起動直後の開始モードから過渡モード、定常モード1、定常モード2を経て終了モードに至る、5つの運転モードを有するとする。それぞれの運転モードM1,M2,M3,M4,M5における、有効電力で表した運転モード消費電力をPL1,PL2,PL3,PL4,PL5と表記する。
なお、高調波解析の手法はFFTに限らない。例えば、ウェーブレット変換などの手法も利用可能である。次に、上記構成における作用を説明する。
FIG. 10 is a diagram illustrating an example of the transition of the operation mode of the electric device. It is assumed that an electrical device has five operation modes from a start mode immediately after startup to a transition mode, a steady mode 1, a steady mode 2, and an end mode. In each of the operation modes M1, M2, M3, M4, and M5, the operation mode power consumption expressed as active power is expressed as PL1, PL2, PL3, PL4, and PL5.
Note that the method of harmonic analysis is not limited to FFT. For example, a technique such as wavelet transform can be used. Next, the operation of the above configuration will be described.

図11は、推定部13における処理の一例を模式的に示す図である。先ず、式(2)から算出されるi番目の次数の高調波含有率をαiとし、運転モード消費電力をPLとし、運転モード消費電力のi番目の次数の高調波の大きさをPL_iとすれば、式(3)の関係が成り立つ。

Figure 0005872441
FIG. 11 is a diagram schematically illustrating an example of processing in the estimation unit 13. First, let the i-th order harmonic content calculated from Equation (2) be αi, the operation mode power consumption be PL, and the i-th order harmonic magnitude of the operation mode power consumption be PL_i. For example, the relationship of Expression (3) is established.
Figure 0005872441

N次までの高調波を含む運転モードの運転モード消費電力PLは、式(3)を参照して式(4)により表すことができる。

Figure 0005872441
The operation mode power consumption PL of the operation mode including the harmonics up to the Nth order can be expressed by Expression (4) with reference to Expression (3).
Figure 0005872441

運転モードMjの運転モード消費電力をPLjとし、i番目の次数の高調波含有率をαijとすれば、式(5)の関係が成り立つ。

Figure 0005872441
When the operation mode power consumption of the operation mode Mj is PLj and the harmonic content of the i-th order is αij, the relationship of Expression (5) is established.
Figure 0005872441

さて、例えば図10における各運転モード消費電力PL1,PL2,PL3,PL4,PL5と、測定装置20を通過する有効電力Pとの間には、ロス分を無視できるとすると次式(6)の関係がある。

Figure 0005872441
Now, for example, assuming that the loss can be ignored between each operation mode power consumption PL1, PL2, PL3, PL4, PL5 in FIG. 10 and the active power P passing through the measuring device 20, the following equation (6) is satisfied. There is a relationship.
Figure 0005872441

式(6)は複数の電気機器の存在を考慮した式であるが、電気機器を1つとすれば式(6)は式(7)のようにも表わすことができる。

Figure 0005872441
Equation (6) is an equation that considers the presence of a plurality of electric devices. However, if one electric device is provided, equation (6) can also be expressed as equation (7).
Figure 0005872441

PがN次の高調波を含むと仮定し、i次の高調波成分をPiとすれば、式(8)の関係が成り立つ。

Figure 0005872441
Assuming that P includes an Nth-order harmonic and assuming that the i-th harmonic component is Pi, the relationship of Expression (8) is established.
Figure 0005872441

式(8)、式(4)、式(5)から、Pのi次の高調波と、運転モード消費電力のi次の高調波との間には式(9)の関係が成り立つ。

Figure 0005872441
From Expression (8), Expression (4), and Expression (5), the relationship of Expression (9) is established between the i-order harmonic of P and the i-order harmonic of the operation mode power consumption.
Figure 0005872441

式(9)を高調波の次数ごとにまとめ、高調波の次数を5次までで考えると、式(10)に示す行列式が得られる。

Figure 0005872441
When Expression (9) is summarized for each harmonic order and the harmonic orders are considered up to the fifth order, the determinant shown in Expression (10) is obtained.
Figure 0005872441

式(10)を運転モード消費電力PL1、PL2、PL3、PL4、PL5について書き直すと式(11)となる。

Figure 0005872441
When Expression (10) is rewritten for operation mode power consumption PL1, PL2, PL3, PL4, and PL5, Expression (11) is obtained.
Figure 0005872441

式(11)において、αijはデータベース14に記憶される量であり、つまり既知の量である。またP1、・・・、P5は測定装置20により測定された有効電力Pに含まれる高調波成分であるから(式(8))、式(11)に基づいて運転モード消費電力PL1、PL2、PL3、PL4、PL5を求めることができる。   In equation (11), αij is an amount stored in the database 14, that is, a known amount. Since P1,..., P5 are harmonic components included in the active power P measured by the measuring device 20 (Equation (8)), the operation mode power consumption PL1, PL2, PL3, PL4, and PL5 can be obtained.

すなわち図11において推定部13は、高調波成分P1〜P5を保存部12から取得し、高調波含有率αijをデータベース14から取得する。そして推定部13は、高調波成分P1〜P5から式(11)を用いて、高調波含有率αijに基づいて運転モード消費電力PL1、PL2、PL3、PL4、PL5を逆算する。   That is, in FIG. 11, the estimation unit 13 acquires the harmonic components P <b> 1 to P <b> 5 from the storage unit 12, and acquires the harmonic content rate αij from the database 14. And the estimation part 13 carries out back calculation of operation mode power consumption PL1, PL2, PL3, PL4, PL5 based on harmonic content rate alphaij using a harmonic component P1-P5 using Formula (11).

運転モード消費電力PL1、PL2、PL3、PL4、PL5は運転モードと一対一に対応するので、電気機器の運転モードを個別に推定することが可能になる。得られた運転モード消費電力PL1〜PL5および運転モードは、推定結果としてエネルギー管理サーバ50に渡される。   Since the operation mode power consumption PL1, PL2, PL3, PL4, and PL5 correspond to the operation mode on a one-to-one basis, it becomes possible to estimate the operation mode of the electric device individually. The obtained operation mode power consumption PL1 to PL5 and the operation mode are passed to the energy management server 50 as an estimation result.

図12は、推定部13による処理を説明するための模式図である。運転モードM1〜MNはそれぞれ固有の高調波パターンと、有効電力PL1〜PLNを持つ。運転モードM1〜MNの高調波パターンは式(10)のマトリクス[αij]における列ベクトルで示され、データベース14に予め記憶されている
一方、測定装置20により得られた測定量をFFT処理すると、測定点における有効電力の高調波成分(P1〜PN)が得られる。P1〜PNはPL1〜PLNに式(10)のマトリクス[αij]を乗算して得られるので、P1〜PNにマトリクス[αij]の逆行列を乗じればPL1〜PLNを逆算することが可能である。PL1〜PLNを求めることができれば、すなわち運転モードM1〜MNを特定することができる。
FIG. 12 is a schematic diagram for explaining processing by the estimation unit 13. Each of the operation modes M1 to MN has a unique harmonic pattern and active powers PL1 to PLN. The harmonic patterns of the operation modes M1 to MN are indicated by column vectors in the matrix [αij] of the equation (10) and are stored in advance in the database 14. On the other hand, when the measurement amount obtained by the measurement device 20 is FFT processed, Harmonic components (P1 to PN) of active power at the measurement point are obtained. Since P1 to PN are obtained by multiplying PL1 to PLN by the matrix [αij] of Expression (10), PL1 to PLN can be calculated backward by multiplying P1 to PN by the inverse matrix of the matrix [αij]. is there. If PL1 to PLN can be obtained, that is, the operation modes M1 to MN can be specified.

以上述べたように第1の実施形態では、測定した有効電力の波形をFFT処理により高調波解析し、次数ごとの高調波成分P1〜PNを算出する。また、運転モードごとの高調波パターンを予め求めておき、高調波含有率をデータベース14に記憶させる。そして、高調波含有率に基づいて、次数ごとの高調波成分P1〜PNから運転モードごとの有効電力PL1〜PLNを逆算することにより、電気機器の運転モードを推定するようにしている。これにより、一つの測定点(測定装置20)における測定結果から運転モードを電気機器の運転モードを計算して求めることができるようになり、電気機器の運転モードを簡易に推定することが可能となる。   As described above, in the first embodiment, the measured active power waveform is subjected to harmonic analysis by FFT processing, and the harmonic components P1 to PN for each order are calculated. Further, a harmonic pattern for each operation mode is obtained in advance, and the harmonic content is stored in the database 14. And based on the harmonic content rate, the operation mode of the electric equipment is estimated by calculating back the active powers PL1 to PLN for each operation mode from the harmonic components P1 to PN for each order. Accordingly, the operation mode can be obtained by calculating the operation mode of the electric device from the measurement result at one measurement point (measurement device 20), and the operation mode of the electric device can be easily estimated. Become.

つまり第1の実施形態によれば、ローカル系統40に接続される電気機器がビルの1つの階、あるいはビル全体にわたり広域に設置されていても、一つの測定点において得られた測定データに基づいて各電気機器の運転モードを推定することが可能になる。従って多数の計測装置や通信設備などを設けなくとも良いので、設備コストを抑えることが可能になる。また、運転モードを特定することにより、需要家におけるピークカットやピークシフトなどの需給制御もスマートに実施することが可能になる。   That is, according to the first embodiment, even if the electrical equipment connected to the local system 40 is installed in a wide area on one floor of the building or the entire building, it is based on the measurement data obtained at one measurement point. Thus, it becomes possible to estimate the operation mode of each electric device. Therefore, it is not necessary to provide a large number of measuring devices, communication facilities, and the like, so that the equipment cost can be reduced. In addition, by specifying the operation mode, it is possible to smartly implement supply and demand control such as peak cut and peak shift in the consumer.

さらには、デマンドレスポンスへのインテリジェントな対応が可能になる。つまり、電気機器には停止させることで損害を生じる運転モードがある。そこで、実施形態の手法により運転モードを推定し、デマンドレスポンスを要求された際には損害の少ない運転モードにある電気機器から停止させるようにすることで、コストインパクトを最小限に抑えることも可能になる。   Furthermore, intelligent response to demand response becomes possible. In other words, there is an operation mode in which electric devices cause damage when stopped. Therefore, it is possible to minimize the cost impact by estimating the operation mode using the method of the embodiment and stopping from the electrical equipment in the operation mode with less damage when demand response is requested. become.

[第2の実施形態]
図13は、推定部13における処理の他の例を模式的に示す図である。再び図6〜図9を参照すると、運転モードごとに、高調波成分の含有率(高調波パターン)に特徴のあることが分かる。例えば、図6の運転モードM1は次数5,7,11,13,17,19に特徴を持つ。図8の運転モードM3は次数3、5,7,9,11,13,15、17,19に特徴を持つ。図9の運転モードM4は次数3,5,7,9,11に特徴を持つ。
[Second Embodiment]
FIG. 13 is a diagram schematically illustrating another example of processing in the estimation unit 13. 6 to 9 again, it can be seen that there is a characteristic in the content rate (harmonic pattern) of the harmonic component for each operation mode. For example, the operation mode M1 in FIG. 6 is characterized by orders 5, 7, 11, 13, 17, and 19. The operation mode M3 in FIG. 8 is characterized by orders 3, 5, 7, 9, 11, 13, 15, 17, and 19. The operation mode M4 in FIG. 9 is characterized by orders 3, 5, 7, 9, and 11.

そこで、特徴のある次数を選択して、式(11)に代入することで、運転モードをより高い精度で推定することが可能になる。例えば、運転モードM1は17次、19次に、運転モードM3は3次、15次に、運転モードM4は、次数3,5,7に有意な特徴を示すので、例えば次数3,5,7,15,17を式(11)に当てはめると式(11)は式(12)のように示される。

Figure 0005872441
Therefore, it is possible to estimate the operation mode with higher accuracy by selecting a characteristic order and substituting it into the equation (11). For example, the operation mode M1 is significant in the 17th order, the 19th order, the operation mode M3 is the third order, the 15th order, and the operation mode M4 is significant in the orders 3, 5, and 7. , 15, and 17 are applied to the expression (11), the expression (11) is expressed as the expression (12).
Figure 0005872441

図13において、推定部13は、測定点の有効電力の高調波成分のうちP3、P5、P7、P15、P17を保存部12から取得し、データベース14から3次,5次,7次,15次,17次の高調波含有率を取得し、式(12)を用いて運転モード消費電力PL1、PL2、PL3、PL4、PL5を得ることが示される。   In FIG. 13, the estimation unit 13 acquires P3, P5, P7, P15, and P17 among the harmonic components of the active power at the measurement point from the storage unit 12, and obtains the third, fifth, seventh, and 15 from the database 14. It is shown that the next and 17th harmonic content is acquired and the operation mode power consumptions PL1, PL2, PL3, PL4, and PL5 are obtained using Expression (12).

図14は、推定部13における処理の他の例を模式的に示す図である。ここでは、電力ロスを考慮の対象とし、高調波成分ごとの電力ロスを推定して電力ロスによる誤差を補正する手法につき開示する。   FIG. 14 is a diagram schematically illustrating another example of processing in the estimation unit 13. Here, a method for correcting the error due to the power loss by estimating the power loss for each harmonic component with the power loss taken into consideration will be disclosed.

式(9)において、j番目の運転モード消費電力のi次の高調波に起因する電力ロスをPlosijと表記すると、式(13)、式(14)の関係が成り立つ。

Figure 0005872441
In Expression (9), when the power loss due to the i-th harmonic of the j-th operation mode power consumption is expressed as Plosij, the relations of Expression (13) and Expression (14) are established.
Figure 0005872441

式(13)において、rはj番目の運転モードにかかる抵抗を示す。式(14)において、QLijはj番目の運転モードの無効電力のi次の高調波成分を示す。Vijはj番目の運転モードの電圧のi次の高調波成分を示す。aij、bijは、j番目の運転モード消費電力のi次の有効電力高調波成分の近似式にかかる定数を示す。   In Equation (13), r represents the resistance applied to the jth operation mode. In Expression (14), QLij represents the i-th harmonic component of the reactive power in the j-th operation mode. Vij represents the i-th harmonic component of the voltage in the j-th operation mode. aij and bij represent constants relating to the approximate expression of the i-th order active power harmonic component of the j-th operation mode power consumption.

式(13)、式(14)の関係を式(10)に当てはめると、式(10)は式(15)と書き直せる。

Figure 0005872441
When the relationship between Expression (13) and Expression (14) is applied to Expression (10), Expression (10) can be rewritten as Expression (15).
Figure 0005872441

式(15)を運転モード消費電力PL1、PL2、PL3、PL4、PL5について書き直すと式(16)となる。

Figure 0005872441
When Expression (15) is rewritten with respect to operation mode power consumption PL1, PL2, PL3, PL4, and PL5, Expression (16) is obtained.
Figure 0005872441

式(16)において、γijは既知の量と看做すことができ、データベース14に予め記憶される。また、P1、・・・、P5は測定装置20により測定された有効電力Pに含まれる高調波成分であるから(式(8))、式(16)に基づいて運転モード消費電力PL1、PL2、PL3、PL4、PL5を求めることができる。   In Expression (16), γij can be regarded as a known amount, and is stored in the database 14 in advance. Since P1,..., P5 are harmonic components included in the active power P measured by the measuring device 20 (Equation (8)), the operation mode power consumption PL1, PL2 is based on Equation (16). , PL3, PL4, PL5 can be obtained.

すなわち図14において推定部13は、高調波成分P1〜P5を保存部12から取得し、高調波含有率αijを誤差成分で補正した係数γijをデータベース14から取得する。そして推定部13は、高調波成分P1〜P5から式(16)を用いて、補正した高調波含有率γijに基づいて運転モード消費電力PL1、PL2、PL3、PL4、PL5を逆算する。これにより電気機器の運転モードを個別に推定することが可能になる。   That is, in FIG. 14, the estimation unit 13 acquires the harmonic components P1 to P5 from the storage unit 12, and acquires the coefficient γij obtained by correcting the harmonic content rate αij with the error component from the database 14. And the estimation part 13 back-calculates operation mode power consumption PL1, PL2, PL3, PL4, PL5 based on the corrected harmonic content rate (gamma) ij using harmonic expression P1-P5 (16). Thereby, it becomes possible to estimate the operation mode of an electric equipment separately.

このように第2の実施形態では、高調波成分ごとの電力ロスを推定し、電力ロスによる誤差を補正して運転モードを推定することが可能になり、より実態に即した処理を実現することが可能になる。   As described above, in the second embodiment, it is possible to estimate the power loss for each harmonic component, correct the error due to the power loss, and estimate the operation mode, thereby realizing more realistic processing. Is possible.

[第3の実施形態]
図15は、推定部13における処理の他の例を模式的に示す図である。第1の実施形態においては、有効電力Pの高調波成分と、運転モード消費電力の高調波成分との間に式(9)、式(10)の関係があることを示した。式(10)に示される行列は正方行列であるが、実際には行と列の数が異なることがある。第3の実施形態ではこのような事態に対処する手法につき説明する。
[Third Embodiment]
FIG. 15 is a diagram schematically illustrating another example of processing in the estimation unit 13. In the first embodiment, it is shown that there is a relationship of Expression (9) and Expression (10) between the harmonic component of the active power P and the harmonic component of the operation mode power consumption. Although the matrix shown in Expression (10) is a square matrix, the number of rows and columns may actually be different. In the third embodiment, a method for dealing with such a situation will be described.

式(10)の左辺を実測高調波ベクトル、式(10)の行列における列ベクトルを運転モード高調波ベクトルとし、高調波の数をmとし、電気機器5の数をnとすれば、式(10)を式(17)のように近似式で表すことができる。

Figure 0005872441
If the left side of Equation (10) is the measured harmonic vector, the column vector in the matrix of Equation (10) is the operation mode harmonic vector, the number of harmonics is m, and the number of electrical devices 5 is n, then Equation (10) 10) can be expressed by an approximate expression as shown in Expression (17).
Figure 0005872441

式(17)の右辺から求まるベクトルを実測高調波近似値ベクトルP1′,P2′,…,Pm′と表記し、式(17)を式(18)のように表す。

Figure 0005872441
A vector obtained from the right side of Expression (17) is expressed as measured harmonic approximate value vectors P1 ′, P2 ′,..., Pm ′, and Expression (17) is expressed as Expression (18).
Figure 0005872441

第3の実施形態では、式(17)の左辺と式(18)の左辺との差(誤差)を最小にするベクトル(PL1,PL2,…,PLN)を求め、これを運転モード消費電力PL1、PL2、…、PLNとする。具体的には最小二乗法を用いて誤差を最小にする。   In the third embodiment, vectors (PL1, PL2,..., PLN) that minimize the difference (error) between the left side of Expression (17) and the left side of Expression (18) are obtained, and this is calculated as the operation mode power consumption PL1. , PL2,..., PLN. Specifically, the error is minimized using the least square method.

式(17)の左辺のベクトルと式(18)の左辺のベクトルとの差分ベクトルの、要素ごとの二乗(Sと表記する)を最小化すれば上記誤差を最小化できる。これを数式化すると式(19)を得る。

Figure 0005872441
The error can be minimized by minimizing the square of each element (denoted as S) of the difference vector between the vector on the left side of equation (17) and the vector on the left side of equation (18). When this is mathematically expressed, Equation (19) is obtained.
Figure 0005872441

以下、繁雑を避けるためにi=1からi=nまでの総和をΣで簡略表示する。   Hereinafter, in order to avoid complications, the sum from i = 1 to i = n is simply displayed as Σ.

式(19)を運転モード消費電力の各々で偏微分し、これらを0(零)にするベクトル(PL1,PL2,…,PLN)が運転モード消費電力を表す。式(20)は、式(19)をj番目の運転モード消費電力で偏微分した項を示す。

Figure 0005872441
A vector (PL1, PL2,..., PLN) that partially differentiates Equation (19) by each of the operation mode power consumption and sets these to 0 (zero) represents the operation mode power consumption. Equation (20) represents a term obtained by partial differentiation of Equation (19) with respect to the j-th operation mode power consumption.
Figure 0005872441

式(20)を展開すると式(21)を得る。

Figure 0005872441
When formula (20) is expanded, formula (21) is obtained.
Figure 0005872441

式(21)を0とし移項すると、行列式(22)を得る。

Figure 0005872441
Determining equation (21) as 0 and obtaining the determinant (22).
Figure 0005872441

式(22)右辺の行列は対称行列となり、その各要素はΣαki×αji(k列j行またはj列k行)、Σαji×αji(j列j行の対角成分)となる。[高調波次数の数≧運転モードの数]という条件が成り立てば、この対称行列は正則行列(Non-Singular Matrix)となり逆行列を計算することができる。すなわち、式(23)に基づいてベクトル(PL1,PL2,…,PLN)を求めることができ、従って運転モードおよび運転モード消費電力を計算することが可能である。

Figure 0005872441
The matrix on the right side of Expression (22) is a symmetric matrix, and each element thereof is Σαki × αji (k columns j rows or j columns k rows) and Σαji × αji (diagonal components of j columns j rows). If the condition [number of harmonic orders ≧ number of operation modes] is satisfied, this symmetric matrix becomes a regular matrix (Non-Singular Matrix), and the inverse matrix can be calculated. That is, the vector (PL1, PL2,..., PLN) can be obtained based on the equation (23), and therefore the operation mode and the operation mode power consumption can be calculated.
Figure 0005872441

すなわち図15において推定部13は、高調波成分P1〜Pmを保存部12から取得し、式(23)に基づいて運転モード消費電力PL1、PL2、・・・、PLnを算出する。これにより電気機器の運転モードを個別に推定することが可能になる。   That is, in FIG. 15, the estimation part 13 acquires the harmonic components P1-Pm from the preservation | save part 12, and calculates operation mode power consumption PL1, PL2, ..., PLn based on Formula (23). Thereby, it becomes possible to estimate the operation mode of an electric equipment separately.

[第4の実施形態]
図16は、実施形態に係る状態推定装置を適用可能なエネルギー管理システムの他の例を示す図である。図16においては、電気機器に装置の種別があることを明示する。図1においては電気機器30−1,30−2を装置Aとし、電気機器30−3,30−4を装置Bとし、電気機器30−5を装置Cとする。なお状態推定装置に符号100を付して示す。
[Fourth Embodiment]
FIG. 16 is a diagram illustrating another example of the energy management system to which the state estimation device according to the embodiment is applicable. In FIG. 16, it is clearly shown that there is a device type in the electrical equipment. In FIG. 1, the electrical devices 30-1 and 30-2 are referred to as a device A, the electrical devices 30-3 and 30-4 are referred to as a device B, and the electrical device 30-5 is referred to as a device C. In addition, the code | symbol 100 is attached | subjected and shown to a state estimation apparatus.

以下では説明を分かり易くするために、装置A、装置Bが2つの運転モードを有し、装置Cは1つの運転モードを有するとする。装置Aの2つの運転モードにおける運転モード消費電力を、それぞれPLA_1、PLA_2とする。装置Bの2つの運転モードにおける運転モード消費電力をそれぞれPLB_1、PLBA_2とする。装置Cの運転モード消費電力をPLC_1とする。これらの変数を用いて式(6)を式(24)のように表すことができる。

Figure 0005872441
In the following, for ease of explanation, it is assumed that apparatus A and apparatus B have two operation modes, and apparatus C has one operation mode. The operation mode power consumption in the two operation modes of the apparatus A is referred to as PLA_1 and PLA_2, respectively. The operation mode power consumption in the two operation modes of the apparatus B is assumed to be PLB_1 and PLBA_2, respectively. The operation mode power consumption of the apparatus C is set to PLC_1. Using these variables, Expression (6) can be expressed as Expression (24).
Figure 0005872441

第1の実施形態と同様の議論により、式(24)から式(25)のように、PLA_1、PLA_2、PLB_1、PLB_2、PLC_1を求めることができる。

Figure 0005872441
Through the same discussion as in the first embodiment, PLA_1, PLA_2, PLB_1, PLB_2, and PLC_1 can be obtained as in Expression (24) to Expression (25).
Figure 0005872441

電力ロスを考慮すると、式(16)を式(26)のように表すことができる。

Figure 0005872441
Considering the power loss, Expression (16) can be expressed as Expression (26).
Figure 0005872441

以上のように第4の実施形態によれば、電気機器に複数の種別があるケースでも、各装置における消費電力、および運転モードを特定することが可能になる。   As described above, according to the fourth embodiment, it is possible to specify the power consumption and the operation mode in each device even when there are a plurality of types of electrical equipment.

[第5の実施形態]
図17は、一つの装置の運転モードの推移の一例を示す模式図である。図17において横軸は時間を示し、縦軸は各運転モードにおける消費電力Pを示す。簡単のため3つの運転モードを想定し、モードM1を過渡モード、モードM2を定常モード、モードM3を終了モードとする。P1、P2、P3は、それぞれモードM1の装置1台あたりの消費電力、モードM2の装置2台あたりの消費電力、モードM3の装置1台あたりの消費電力を示す。
[Fifth Embodiment]
FIG. 17 is a schematic diagram illustrating an example of the transition of the operation mode of one device. In FIG. 17, the horizontal axis indicates time, and the vertical axis indicates power consumption P in each operation mode. For simplicity, three operation modes are assumed, mode M1 is a transient mode, mode M2 is a steady mode, and mode M3 is an end mode. P1, P2, and P3 respectively indicate power consumption per device in mode M1, power consumption per device in mode M2, and power consumption per device in mode M3.

図18は、複数の装置の運転モードの推移の一例を示す模式図である。つまり図18は図17に示す消費電力を全ての装置にわたり総和した値を示すグラフであり、測定点で測定されるモードM1、モードM2、モードM3の各消費電力を示す。当然ながら図18の縦軸スケールは図17の縦軸スケールよりも長い。図18におけるP1m、P2m、P3mは、それぞれモードM1の装置の消費電力の総和、モードM2の装置の消費電力の総和、モードM3の装置の消費電力の総和を示す。   FIG. 18 is a schematic diagram illustrating an example of transition of operation modes of a plurality of apparatuses. That is, FIG. 18 is a graph showing the sum of the power consumption shown in FIG. 17 over all the devices, and shows the power consumption in the modes M1, M2, and M3 measured at the measurement points. Of course, the vertical scale of FIG. 18 is longer than the vertical scale of FIG. P1m, P2m, and P3m in FIG. 18 indicate the total power consumption of the mode M1 device, the total power consumption of the mode M2 device, and the total power consumption of the mode M3 device, respectively.

各モードにおける装置1台の消費電力は装置の定格値として既知であり、データベース14に記憶されている。そこで、運転モードM1の装置稼働台数をn1、運転モードM2の装置稼働台数をn2、運転モードM3の装置稼働台数をn3とすると、式(27)が成り立つ。

Figure 0005872441
The power consumption of one device in each mode is known as the rated value of the device and is stored in the database 14. Therefore, when the number of operating devices in the operation mode M1 is n1, the number of operating devices in the operation mode M2 is n2, and the number of operating devices in the operation mode M3 is n3, the equation (27) is established.
Figure 0005872441

以上述べたように第5の実施形態では、高調波解析により得られた運転モードごとの消費電力を、装置1台あたりの消費電力で除算することで、運転モードごとの稼働台数を算出することが可能である。   As described above, in the fifth embodiment, the number of operating units for each operation mode is calculated by dividing the power consumption for each operation mode obtained by the harmonic analysis by the power consumption per device. Is possible.

[第6の実施形態]
図19は、運転を停止した場合の損害を運転モードごとに示す模式図である。図19のグラフの縦軸は、損害の大きさを表す指標Prを示す。図中、Pr1は、モードM1の装置1台が運転を停止した場合の損害を示す。Pr2は、モードM2の装置1台が運転を停止した場合の損害を示す。Pr3は、モードM3の装置1台が運転を停止した場合の損害を示す。データベース14は、運転モードごとに、稼動停止に伴う損害(リスク)を意味する指標Prを対応付けた情報を予め記憶する。
[Sixth Embodiment]
FIG. 19 is a schematic diagram showing damage when the operation is stopped for each operation mode. The vertical axis of the graph in FIG. 19 indicates an index Pr representing the magnitude of damage. In the figure, Pr1 indicates damage when one device in mode M1 stops operation. Pr2 indicates damage when one device in mode M2 stops operation. Pr3 indicates damage when one device in mode M3 stops operation. The database 14 stores in advance information that associates an index Pr that means damage (risk) associated with the operation stop for each operation mode.

図20は、全ての装置が停止した時の損害Pra、稼働台数n、および消費電力Pmをモードごとに配列して示す図である。Pra、n、およびPmの添え字1,2,3はモードM1,M2,M3に対応する。式(27)によりモードごとの稼働台数n1,n2,n3を知り得るので、Pr1,Pr2,Pr3に稼働台数n1,n2,n3を乗じればトータルでの損害Pr1a,Pr2a,Pr3aを求めることができる。P1m,P2m,P3mの求め方については既に説明したのでここでは省略する。   FIG. 20 is a diagram showing the damage Pra, the number n of operations, and the power consumption Pm when all devices are stopped arranged for each mode. Subscripts 1, 2, and 3 of Pra, n, and Pm correspond to modes M1, M2, and M3. Since the operating numbers n1, n2, and n3 for each mode can be obtained from the equation (27), the total damage Pr1a, Pr2a, and Pr3a can be obtained by multiplying the operating numbers n1, n2, and n3 by Pr1, Pr2, and Pr3. it can. Since the method for obtaining P1m, P2m, and P3m has already been described, it is omitted here.

図20によれば、停止時の損害はモードM2が最小で、モードM1、モードM3の順に大きくなることがわかる。逆に、稼動停止による負荷の低減効果は、モードM3、モードM1、モードM2の順に高いことが示される。推定部13は、推定した運転モードの電気機器を停止した場合のリスクも併せてエネルギー管理サーバ50に出力する。この結果を参照して、エネルギー管理サーバ50は、例えばデマンドレスポンスに対する、より有利な運転計画を立案することが可能になる。つまりエネルギー管理サーバ50は、推定された電気機器の運転モードと停止した際のリスクとに基づいて、稼動停止の可否の優先度を判定することができるようになる。   According to FIG. 20, it can be seen that the damage at the time of stoppage is the smallest in mode M2 and increases in the order of mode M1 and mode M3. On the contrary, the load reduction effect due to the operation stop is shown to be higher in the order of mode M3, mode M1, and mode M2. The estimation unit 13 also outputs to the energy management server 50 a risk when the electrical device in the estimated operation mode is stopped. With reference to this result, the energy management server 50 can make a more advantageous operation plan for a demand response, for example. That is, the energy management server 50 can determine the priority of whether or not the operation can be stopped based on the estimated operation mode of the electric device and the risk when the electric device is stopped.

例えば真夏に空調機器を停止すると、労働者の作業効率が低下して損害を生じることがある。あるいは工場における稼働中の製造機器が停止すると、製品の歩留まりが極端に悪化して大きな損害を生じるケースもある。第6の実施形態によればこのような事態に備えることが可能になる。   For example, if the air conditioner is stopped in midsummer, the work efficiency of workers may be reduced, causing damage. Alternatively, when the manufacturing equipment in operation in the factory stops, there is a case where the yield of the product is extremely deteriorated to cause great damage. According to the sixth embodiment, it is possible to prepare for such a situation.

[第7の実施形態]
図21は、状態推定装置の他の例を示す機能ブロック図である。図21においては図2と共通する部分に同じ符号を付して示し、ここでは異なる箇所についてのみ説明する。なお状態推定装置に符号70を付し、図2の構成と区別する。
状態推定装置70は、推定部13により推定した運転モード消費電力を、結果表示部71に渡す。結果表示部71は、ビル監視室などに設置される表示装置(図示せず)に運転モード消費電力をトレンドグラフなどの形態で、グラフィカルに表示する。結果表示部71は表示内容を例えば定期的に、あるいは運転モード消費電力の算出のたびごとに更新する。
[Seventh Embodiment]
FIG. 21 is a functional block diagram illustrating another example of the state estimation device. In FIG. 21, parts that are the same as those in FIG. 2 are given the same reference numerals, and only different parts will be described here. Note that the state estimation device is denoted by reference numeral 70 to distinguish it from the configuration of FIG.
The state estimation device 70 passes the operation mode power consumption estimated by the estimation unit 13 to the result display unit 71. The result display unit 71 graphically displays the operation mode power consumption in a form such as a trend graph on a display device (not shown) installed in a building monitoring room or the like. The result display unit 71 updates the display content, for example, regularly or every time the operation mode power consumption is calculated.

図22は、結果表示部71により表示装置に表示される運転モード消費電力の一例を示す図である。グラフの縦軸が運転モード消費電力を示し、横軸は時間(時刻)を示す。図22においては一日24時間にわたる結果を示す。図22においては運転モードM1、運転モードM2、運転モードM3、運転モードM4の各消費電力を一日24時間分プロットしたグラフが示される。   FIG. 22 is a diagram illustrating an example of operation mode power consumption displayed on the display device by the result display unit 71. The vertical axis of the graph represents the operation mode power consumption, and the horizontal axis represents time (time). In FIG. 22, the result over 24 hours a day is shown. FIG. 22 shows a graph in which each power consumption in the operation mode M1, the operation mode M2, the operation mode M3, and the operation mode M4 is plotted for 24 hours a day.

以上のように第7の実施形態では、運転モードの推定結果とその消費電力を視覚的に表示するようにしている。これによりユーザ(ビルの管理者など)はビル内における電気機器の運転モードとその消費電力を視覚的に把握できるようになり、電力需要のピークカットやピークシフトなどの需給制御に役立てることが可能になる。   As described above, in the seventh embodiment, the estimation result of the operation mode and the power consumption thereof are visually displayed. This allows users (such as building managers) to visually understand the operation mode and power consumption of electrical equipment in the building, which can be used for supply and demand control such as peak cut and peak shift of power demand. become.

[第8の実施形態]
図23は、各モードの継続時間を装置ごとに示すグラフの一例を示す図である。運転モードの継続時間は装置ごとに予め既知であるので、これを例えばテーブル化してデータベース14に記憶させておくようにする。推定部13は、装置の稼動停止に伴う損害と負荷削減との関係(図6:第6の実施形態)に加えて、図23に示される運転モードの継続時間をデータベース14から取得する。そして推定部13は、これらの情報に基づいて、推定された運転モードにある電気機器の負荷制御のタイミングを最適化する。
[Eighth Embodiment]
FIG. 23 is a diagram illustrating an example of a graph indicating the duration of each mode for each device. Since the duration of the operation mode is known in advance for each apparatus, this is, for example, tabulated and stored in the database 14. The estimation unit 13 acquires the duration of the operation mode shown in FIG. 23 from the database 14 in addition to the relationship between the damage caused by the operation stop of the apparatus and the load reduction (FIG. 6: sixth embodiment). And the estimation part 13 optimizes the timing of the load control of the electric equipment in the estimated operation mode based on such information.

つまり推定部13は、最大の効果を得ることの可能な装置の稼動停止のタイミングを、運転モードの継続時間を考慮して算出する。これにより、稼動停止に伴う損害と、負荷制御によるコスト抑制との双方のバランスを取った制御が可能となる。   That is, the estimation unit 13 calculates the operation stop timing of the apparatus that can obtain the maximum effect in consideration of the duration of the operation mode. As a result, it is possible to perform control that balances both damage due to operation stop and cost reduction by load control.

本発明のいくつかの実施形態を説明したが、これらの実施形態は例として提示するものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。   Although several embodiments of the present invention have been described, these embodiments are presented by way of example and are not intended to limit the scope of the invention. These novel embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the scope of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, and are included in the invention described in the claims and the equivalents thereof.

30−1〜30−5…電気機器、40…ローカル系統、20…測定装置、10,70,100…状態推定装置、50…エネルギー管理サーバ、11…計算部、12…保存部、13…推定部、14…データベース、71…結果表示部   30-1 to 30-5 ... electric equipment, 40 ... local system, 20 ... measuring device, 10, 70,100 ... state estimation device, 50 ... energy management server, 11 ... calculation unit, 12 ... storage unit, 13 ... estimation Part, 14 ... database, 71 ... result display part

Claims (10)

電気機器の接続される給電系統に生じる電力の波形を高調波解析して有効電力高調波の次数ごとの高調波成分を算出する計算部と、
前記高調波の次数ごとの高調波含有率を前記電気機器の運転モードごとに記憶するデータベースと、
前記算出された次数ごとの高調波成分から、前記データベースに記憶される高調波含有率に基づいて前記運転モードごとの消費電力を逆算し、その結果に基づいて前記電気機器の運転モードを推定する推定部とを具備する、状態推定装置。
A calculation unit that calculates a harmonic component for each order of active power harmonics by performing harmonic analysis on a power waveform generated in a power feeding system to which electrical equipment is connected;
A database for storing the harmonic content for each order of the harmonics for each operation mode of the electrical device;
From the calculated harmonic components for each order, the power consumption for each operation mode is calculated based on the harmonic content rate stored in the database, and the operation mode of the electrical device is estimated based on the result. A state estimation device comprising an estimation unit.
前記データベースは、運転モードごとの消費電力を前記電気機器の定格値として記憶し、
前記推定部は、前記推定された運転モードにおける前記電気機器の数を、前記逆算された消費電力と前記定格値とに基づいて推定する、請求項1に記載の状態推定装置。
The database stores the power consumption for each operation mode as the rated value of the electrical equipment,
The state estimation apparatus according to claim 1, wherein the estimation unit estimates the number of the electrical devices in the estimated operation mode based on the back-calculated power consumption and the rated value.
前記データベースは、稼動停止に伴うリスクを運転モードごとに対応付けた情報を予め記憶し、
前記推定部は、前記推定された運転モードの電気機器を停止した場合のリスクを出力する、請求項1に記載の状態推定装置。
The database previously stores information associating risks associated with operation suspensions for each operation mode,
The state estimation device according to claim 1, wherein the estimation unit outputs a risk when the electrical device in the estimated operation mode is stopped.
さらに、前記推定された運転モードを視覚的に表示する表示制御部を具備する、請求項1に記載の状態推定装置。   The state estimation apparatus according to claim 1, further comprising a display control unit that visually displays the estimated operation mode. 前記データベースは、前記運転モードごとの継続時間を記憶し、
前記推定部は、前記推定された運転モードの前記継続時間に基づいて、当該運転モードにある電気機器の負荷制御のタイミングを最適化する、請求項1に記載の状態推定装置。
The database stores a duration for each operation mode,
The state estimation device according to claim 1, wherein the estimation unit optimizes the timing of load control of an electric device in the operation mode based on the duration of the estimated operation mode.
コンピュータにより電気機器の状態を推定する状態推定方法であって、
前記コンピュータが、電気機器の接続される給電系統に生じる電力の波形を高調波解析し、
前記コンピュータが、有効電力高調波の次数ごとの高調波成分を算出し、
前記コンピュータが、前記高調波の次数ごとの高調波含有率を前記電気機器の運転モードごとにデータベースに記憶し、
前記コンピュータが、前記算出された次数ごとの高調波成分から、前記データベースに記憶される高調波含有率に基づいて前記運転モードごとの消費電力を逆算し、
前記コンピュータが、前記逆算された消費電力に基づいて前記電気機器の運転モードを推定する、状態推定方法。
A state estimation method for estimating the state of an electrical device by a computer,
The computer performs harmonic analysis of the waveform of power generated in a power feeding system to which electrical equipment is connected,
The computer calculates a harmonic component for each order of active power harmonics,
The computer stores the harmonic content for each harmonic order in a database for each operation mode of the electrical equipment,
The computer reversely calculates the power consumption for each operation mode based on the harmonic content stored in the database from the calculated harmonic components for each order.
A state estimation method in which the computer estimates an operation mode of the electric device based on the power consumption calculated in reverse.
前記コンピュータは、運転モードごとの消費電力を前記電気機器の定格値として前記データベースに記憶し、
前記推定することは、前記推定された運転モードにおける前記電気機器の数を、前記逆算された消費電力と前記定格値とに基づいて推定する、請求項6に記載の状態推定方法。
The computer stores power consumption for each operation mode in the database as a rated value of the electrical equipment,
The state estimation method according to claim 6, wherein the estimating includes estimating the number of the electric devices in the estimated operation mode based on the back-calculated power consumption and the rated value.
前記コンピュータは、稼動停止に伴うリスクを運転モードごとに対応付けた情報を予め前記データベースに記憶し、
前記推定することは、前記推定された運転モードの電気機器を停止した場合のリスクを出力する、請求項6に記載の状態推定方法。
The computer stores in advance in the database information that associates risks associated with the suspension of operation for each operation mode,
The state estimation method according to claim 6, wherein the estimating outputs a risk when the electrical device in the estimated operation mode is stopped.
前記コンピュータは、さらに、前記推定された運転モードを視覚的に表示する、請求項6に記載の状態推定方法。   The state estimation method according to claim 6, wherein the computer further visually displays the estimated operation mode. 前記コンピュータは、前記運転モードごとの継続時間を前記データベースに記憶し、
前記推定することは、前記推定された運転モードの前記継続時間に基づいて、当該運転モードにある電気機器の負荷制御のタイミングを最適化する、請求項6に記載の状態推定方法。
The computer stores the duration for each operation mode in the database,
The state estimation method according to claim 6, wherein the estimating includes optimizing a load control timing of an electric device in the operation mode based on the duration of the estimated operation mode.
JP2012256239A 2012-11-22 2012-11-22 State estimation device and state estimation method Active JP5872441B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012256239A JP5872441B2 (en) 2012-11-22 2012-11-22 State estimation device and state estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012256239A JP5872441B2 (en) 2012-11-22 2012-11-22 State estimation device and state estimation method

Publications (2)

Publication Number Publication Date
JP2014103830A JP2014103830A (en) 2014-06-05
JP5872441B2 true JP5872441B2 (en) 2016-03-01

Family

ID=51025878

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012256239A Active JP5872441B2 (en) 2012-11-22 2012-11-22 State estimation device and state estimation method

Country Status (1)

Country Link
JP (1) JP5872441B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104682386A (en) * 2015-03-10 2015-06-03 江苏省电力公司淮安供电公司 Control method for panoramic operation mode of regional power grid
CN112018760B (en) * 2020-08-27 2021-11-23 四川大学 Harmonic state estimation method based on asynchronous data

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6816078B2 (en) * 2000-04-12 2004-11-09 Central Research Institute Of Electric Power Industry System and method for estimating power consumption of electric apparatus, and abnormality alarm system utilizing the same
JP2002311068A (en) * 2001-04-10 2002-10-23 Yokogawa Electric Corp Logarithmic graph display device
JP2006158082A (en) * 2004-11-29 2006-06-15 Fuji Electric Systems Co Ltd Load operation status estimation device and method, and ground-fault generation source estimation device and method using therefor
JP4802129B2 (en) * 2007-03-19 2011-10-26 株式会社東芝 Power quality evaluation system
JP5127608B2 (en) * 2008-07-09 2013-01-23 株式会社東芝 Power quality evaluation system

Also Published As

Publication number Publication date
JP2014103830A (en) 2014-06-05

Similar Documents

Publication Publication Date Title
JP4705563B2 (en) Distribution system state estimation device, state estimation method and program thereof
Gu et al. Linearized voltage stability index for wide-area voltage monitoring and control
JP5957397B2 (en) State estimation device, state estimation method, and state estimation program
WO2016158659A1 (en) Power distribution monitoring control device
JP5517724B2 (en) Facility management apparatus and facility management method
KR20210055716A (en) Technology for data analysis from smart meters
KR20180070842A (en) Peak power detection and alarm system
JP2008236876A (en) Power quality evaluation system
JP5665250B2 (en) Power quality evaluation system
JP5872441B2 (en) State estimation device and state estimation method
JP5840541B2 (en) POWER LOAD ESTIMATION DEVICE, POWER LOAD ESTIMATION METHOD, AND POWER LOAD ESTIMATION PROGRAM
JP5818089B2 (en) Demand power control apparatus and power demand control method
CN110807539A (en) Fault equipment maintenance path selection method and device
JP2014153942A (en) Energy management support device, control method of energy management support device, control program and recording medium
WO2017155047A1 (en) Power distribution monitoring control system
JP2015225637A (en) Correlation analysis device, correlation analysis method, and program for correlation analysis
Bertho et al. Optimized power quality monitor placement based on a particle swarm optimization algorithm
JP2017103862A (en) Power estimating apparatus, power estimating method, and power estimating program
Jo et al. Scheduling of air-conditioner using occupancy prediction in a smart home/building environment
Kaur et al. On the use of corrective maintenance data for performance analysis of preventively maintained textile industry
JP6100188B2 (en) POWER LOAD ESTIMATION DEVICE, POWER LOAD ESTIMATION SYSTEM, POWER LOAD ESTIMATION METHOD, AND POWER LOAD ESTIMATION PROGRAM
JP2012017930A (en) Method of estimating used energy of apparatus in air conditioning system, and used energy estimating device incorporating the method
KR101743057B1 (en) Apparatus and method for determining stability of power system
JP6500712B2 (en) Power generation facility management device
Almehizia et al. Power system state estimation accuracy and observability evaluation with optimal PMU placement

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150303

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20151208

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: 20151215

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160113

R151 Written notification of patent or utility model registration

Ref document number: 5872441

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

R154 Certificate of patent or utility model (reissue)

Free format text: JAPANESE INTERMEDIATE CODE: R154