JP2021027775A - 発電制御システム及び発電制御方法 - Google Patents
発電制御システム及び発電制御方法 Download PDFInfo
- Publication number
- JP2021027775A JP2021027775A JP2019146624A JP2019146624A JP2021027775A JP 2021027775 A JP2021027775 A JP 2021027775A JP 2019146624 A JP2019146624 A JP 2019146624A JP 2019146624 A JP2019146624 A JP 2019146624A JP 2021027775 A JP2021027775 A JP 2021027775A
- Authority
- JP
- Japan
- Prior art keywords
- power consumption
- power
- power generation
- time
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/20—Climate change mitigation technologies for sector-wide applications using renewable energy
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Control Of Electrical Variables (AREA)
Abstract
【課題】過剰に太陽電池の発電量を制限することを防止し、発電効率を向上させることができる太陽光発電システム及び太陽光発電制御方法を提供することを課題とする。【解決手段】太陽光発電システムは、太陽電池とパワーコンデショナを備えており、複数の直近の消費電力量の時間変化から、消費電力量の予測値を算出し、消費電力の予測値に基づいて指令値を求め、その指令値に従って、パワーコンディショナーが太陽電池の発電量を制御する。また、消費電力量の予測値は、３点以上の時刻の消費電力量から二次回帰関数により、算出することができる。【選択図】 図１
Description
本発明は、発電制御システム及び発電制御方法に関する。
従来、太陽光発電システムは、電力会社との売買契約に基づいて、その余剰電力は商用電力線に逆潮流させ電力会社に売電されていた。しかし、太陽光発電のような分散電源が増加するにともない、逆潮流による電力系統の電圧変動という弊害が生じることとなった。
そのため、現在では太陽光発電システムから電力会社への逆潮流を回避しなければならない場合があり、太陽電池が発電する電力量（発電量）を制御する必要が生じてきた。
太陽電池の発電量制御として、逆潮流を抑制するための制御方法が公知である。
そのため、現在では太陽光発電システムから電力会社への逆潮流を回避しなければならない場合があり、太陽電池が発電する電力量（発電量）を制御する必要が生じてきた。
太陽電池の発電量制御として、逆潮流を抑制するための制御方法が公知である。
現実の消費電力量は、時々刻々変化する一方で、パワーコンディショナーに太陽電池の発電量の制御にはタイムラグが発生する。そのため、発電量が消費電力量を超過しないように、消費電力量の変化を想定して、発電量の上限を設定する必要がある。すなわち、分単位や秒単位の短い時間で消費電力量が変化するため、現実的には、この変化を考慮して、過剰にマージン（余裕）を持たせた低い上限値での発電制御を行うことになる。その結果、発電効率が低下することになる。
上記課題を鑑み、本発明は、過剰に太陽電池の発電量を制限することを防止し、発電効率を向上させることができる太陽光発電システム及び太陽光発電制御方法を提供することを課題とする。
本発明に係る太陽光発電システムは、
太陽電池と、
パワーコンディショナーと、
制御装置とを備え、
前記制御装置は、消費電力量を取得し、回帰分析に基づいて得られた関数により複数の時刻の消費電力量から所定時間経過後の消費電力量の予測値を算出し、
前記予測値に基づき指令値を算出して、前記パワーコンディショナーに出力し、
前記パワーコンディショナーは、前記指令値に基づき、前記太陽電池の発電量を前記予測値以下に制御することを特徴とする。
太陽電池と、
パワーコンディショナーと、
制御装置とを備え、
前記制御装置は、消費電力量を取得し、回帰分析に基づいて得られた関数により複数の時刻の消費電力量から所定時間経過後の消費電力量の予測値を算出し、
前記予測値に基づき指令値を算出して、前記パワーコンディショナーに出力し、
前記パワーコンディショナーは、前記指令値に基づき、前記太陽電池の発電量を前記予測値以下に制御することを特徴とする。
このような構成とすることで、過剰に太陽電池の発電量を制限することを防止し、太陽電池の発電効率を向上させることが可能な太陽光発電システムを得ることができる。
なお、制御装置は、公知の技術（電力計）により計測された、負荷の消費電力を取得することができる。例えば、受変電設備において、負荷に供給する電力量を計測し、その計測値を制御装置に出力することにより、制御装置は消費電力量を取得できる。
なお、制御装置は、公知の技術（電力計）により計測された、負荷の消費電力を取得することができる。例えば、受変電設備において、負荷に供給する電力量を計測し、その計測値を制御装置に出力することにより、制御装置は消費電力量を取得できる。
また、本発明に係る太陽光発電システムは、
前記関数は、ｍ次多項式（ｍ自然数）であることを特徴とする。
前記関数は、ｍ次多項式（ｍ自然数）であることを特徴とする。
また、本発明に係る太陽光発電システムは、
前記関数は、２次多項式であることを特徴とする。
前記関数は、２次多項式であることを特徴とする。
また、本発明に係る太陽光発電システムは、
自然数ｎを３として、ｎ点の時刻と消費電力量の回帰分析から二次関数を求める第１のステップと、
ｎを１増加して、ｎ点の時刻と消費電力量の回帰分析から二次関数と相関係数を求める第２のステップと、
予め設定された閾値と前記相関係数とを比較する第３ステップとをこの順に実行し、
前記第３のステップにおいて前記相関係数が前記閾値以上であれば、前記第２のステップを実行し、前記第３のステップにおいて、前記相関係数が前記閾値より小さければ、ｎから１を減じたｎ点の時刻と消費電力量についての回帰分析により求めた二次関数から、消費電力量の前記予測値を算出することを特徴とする。
自然数ｎを３として、ｎ点の時刻と消費電力量の回帰分析から二次関数を求める第１のステップと、
ｎを１増加して、ｎ点の時刻と消費電力量の回帰分析から二次関数と相関係数を求める第２のステップと、
予め設定された閾値と前記相関係数とを比較する第３ステップとをこの順に実行し、
前記第３のステップにおいて前記相関係数が前記閾値以上であれば、前記第２のステップを実行し、前記第３のステップにおいて、前記相関係数が前記閾値より小さければ、ｎから１を減じたｎ点の時刻と消費電力量についての回帰分析により求めた二次関数から、消費電力量の前記予測値を算出することを特徴とする。
このような構成とすることで、消費電力量の計測データを用いてパワーコンディショナーに対する指令値を出力することができる。
また、消費電力量の予測のための関数が代数的に求められ、容易に消費電力量の予測が可能となる。
また、消費電力量の予測のための関数が代数的に求められ、容易に消費電力量の予測が可能となる。
また、本発明に係る太陽光発電システムは、
前記所定時間は、制御遅れ時間の２倍から１０倍であることを特徴とする。
前記所定時間は、制御遅れ時間の２倍から１０倍であることを特徴とする。
このような構成とすることで、指令値が、過剰にマージンを持つことを防止することができる。
また、本発明に係る太陽光発電システムは、
取得した消費電力量の代わりに、直近の消費電力量の移動平均を用いることを特徴とする。
取得した消費電力量の代わりに、直近の消費電力量の移動平均を用いることを特徴とする。
このような構成とすることで、消費電力量の短周期の変化が大きい場合でも、過剰な消費電力量の予測を防止することができる。
本発明に係る発電制御方法は、
太陽電池と、
パワーコンディショナーとを備えた太陽光発電システムにおいて、
複数の時刻の消費電力量から、所定時間経過後の消費電力量の予測値を算出する第１の工程と、
前記予測値を基に前記パワーコンディショナーに対する指令値を算出する第２の工程と、
前記パワーコンディショナーが前記指令値に従って前記太陽電池の発電量を前記予測値以下に制御する第３の工程とを有することを特徴とする。
太陽電池と、
パワーコンディショナーとを備えた太陽光発電システムにおいて、
複数の時刻の消費電力量から、所定時間経過後の消費電力量の予測値を算出する第１の工程と、
前記予測値を基に前記パワーコンディショナーに対する指令値を算出する第２の工程と、
前記パワーコンディショナーが前記指令値に従って前記太陽電池の発電量を前記予測値以下に制御する第３の工程とを有することを特徴とする。
このような発電制御方法とすることで、過剰に太陽電池の発電量を制限することを防止し、太陽電池の発電効率を向上させることができる。
また、本発明に係る発電制御方法は、
前記第１の工程は、
自然数ｎを３として、ｎ点の時刻と消費電力量の回帰分析から二次関数を求める第１のステップと、
ｎを１増加して、ｎ点の時刻と消費電力量の回帰分析から二次関数と相関係数を求める第２のステップと、
予め設定された閾値と前記相関係数とを比較する第３ステップとをこの順に実行し、
前記第３のステップにおいて前記相関係数が前記閾値以上であれば、前記第２のステップを実行し、前記第３のステップにおいて、前記相関係数が前記閾値より小さければ、ｎから１を減じたｎ点の時刻と消費電力量の回帰分析から求めた二次関数から消費電力量の前記予測値を算出することを含むことを特徴とする。
前記第１の工程は、
自然数ｎを３として、ｎ点の時刻と消費電力量の回帰分析から二次関数を求める第１のステップと、
ｎを１増加して、ｎ点の時刻と消費電力量の回帰分析から二次関数と相関係数を求める第２のステップと、
予め設定された閾値と前記相関係数とを比較する第３ステップとをこの順に実行し、
前記第３のステップにおいて前記相関係数が前記閾値以上であれば、前記第２のステップを実行し、前記第３のステップにおいて、前記相関係数が前記閾値より小さければ、ｎから１を減じたｎ点の時刻と消費電力量の回帰分析から求めた二次関数から消費電力量の前記予測値を算出することを含むことを特徴とする。
このような発電制御方法とすることで、消費電力量の計測データから消費電力量を予測することが可能となり、太陽電池の発電効率を向上させることができる。
また、本発明に係る発電制御方法は、取得した消費電力量の代わりに、直近の消費電力量の移動平均を用いることを特徴とする。
このような発電制御方法とすることで、消費電力量の短周期の変化が大きい場合でも、過剰な消費電力量の予測を防止して、太陽電池の発電効率を向上させることができる。
本発明によれば、過剰に太陽電池の発電量を制限することを防止し、発電効率を向上させることができる太陽光発電システム及び太陽光発電制御方法を得ることができる。
以下、図面を参照して本発明の実施形態について説明する。但し、以下の実施形態は、いずれも本発明の要旨の認定において限定的な解釈を与えるものではない。また、同一又は同種の部材については同じ参照符号を付して、説明を省略することがある。
（第１の実施形態）
図１は、本発明にかかる太陽光発電制御システム１（以下、単に発電制御システムと称す。）の一実施形態の主要構成を示す。
発電制御システム１は、太陽電池２、ＰＣＳ（パワーコンディショナ）３及び発電制御装置４（以下単に、制御装置と称す）を備え、太陽電池２により発電された電力を自家消費するシステムである。
なお、発電制御システム１は、受変電部５（配電盤）を含んでも良い。さらに発電制御システム１は、太陽電池２により発電された電力を蓄える蓄電池を備えてもよい。
図１は、本発明にかかる太陽光発電制御システム１（以下、単に発電制御システムと称す。）の一実施形態の主要構成を示す。
発電制御システム１は、太陽電池２、ＰＣＳ（パワーコンディショナ）３及び発電制御装置４（以下単に、制御装置と称す）を備え、太陽電池２により発電された電力を自家消費するシステムである。
なお、発電制御システム１は、受変電部５（配電盤）を含んでも良い。さらに発電制御システム１は、太陽電池２により発電された電力を蓄える蓄電池を備えてもよい。
ＰＣＳ３は太陽電池２から出力された直流電力を交流電力に変換するとともに太陽電池２の発電量（発電電力）を制御する。すなわち、ＰＣＳ３は、太陽電池２のＩ−Ｖ特性に従って太陽電池２の発電量を制御することができ、例えば公知のＭＰＰＴ法により発電量が最大となるように制御することができる。
ＰＣＳ３は、電力線によって受変電部５と電気的に接続されており、さらに受変電部５は、例えばエアコン等の１つ以上の電気設備である負荷６と電気的に接続されている。
従って、発電された電力は、ＰＣＳ３及び受変電部５を経由して、負荷６に送電され消費される。
従って、発電された電力は、ＰＣＳ３及び受変電部５を経由して、負荷６に送電され消費される。
また、受変電部５は、電力会社７（外部電源）等からの商用電力線に接続されており、商用電力線から電力供給を受け負荷６に電力を供給することができる。従って、受変電部５は、複数の電力源に対して電力線により電気的に接続されており、それらの複数の電力源から受電した電力を、電気的に接続されている負荷６へと電力の分配を行う。
さらに受変電部５は、電力源の供給電圧と負荷６の電気設備の規格電圧に差がある場合、電気設備の規格電圧に応じて変電を行う変圧器の機能を備えてもよい。
さらに受変電部５は、電力源の供給電圧と負荷６の電気設備の規格電圧に差がある場合、電気設備の規格電圧に応じて変電を行う変圧器の機能を備えてもよい。
制御装置４は、太陽電池２の発電量の上限値を設定するため、負荷６の消費電力量、又は受変電部５から負荷６への供給電力を計測する計測部を備える。
なお、受変電部５側に負荷６へ供給電力を計測する機能がある場合、制御装置４は、受変電部５から負荷への供給電力の値を取得すればよい。
また、電力の計測部は、受変電部５内に設けてもよく、制御装置４に内蔵して計測部のみを電力線に設置してもよく、別途独立して設けても良い。計測部は公知の電力計を用いることができる。例えば、電力会社等が消費電力量を測定する際に使用する電力計を用いることができる。
なお、消費電力量の取得方法は、上記に限定されず、公知の方法により取得することができ、例えば、電力会社等からの買電力量と太陽電池２の発電量の総和から算出してもよい。
なお、受変電部５側に負荷６へ供給電力を計測する機能がある場合、制御装置４は、受変電部５から負荷への供給電力の値を取得すればよい。
また、電力の計測部は、受変電部５内に設けてもよく、制御装置４に内蔵して計測部のみを電力線に設置してもよく、別途独立して設けても良い。計測部は公知の電力計を用いることができる。例えば、電力会社等が消費電力量を測定する際に使用する電力計を用いることができる。
なお、消費電力量の取得方法は、上記に限定されず、公知の方法により取得することができ、例えば、電力会社等からの買電力量と太陽電池２の発電量の総和から算出してもよい。
制御装置４は、消費電力量を超過しないように、計測した消費電力量に基づき、発電量の上限を算出し、ＰＣＳ３に対して、発電量の上限（又は上限を指定する値）を指令値として出力する。ＰＣＳ３は、指令値に従って、太陽電池２の発電量を、指定された上限以下となるように制御する。
図２は、消費電力量とＰＣＳ３に対する指令値の時間的推移を示すグラフである。図２において、点線は消費電力量、実線は指令値を示す。指令値は、例えば消費電力量の一次関数や消費電力量から所定量減ずることにより、消費電力量より低い値に設定されている。
図２（ａ）に観察されるように、消費電力量の変化は、相対的に緩やかに変化する曲線に短周期で変化する成分が重畳された挙動を示す。
図２は、消費電力量とＰＣＳ３に対する指令値の時間的推移を示すグラフである。図２において、点線は消費電力量、実線は指令値を示す。指令値は、例えば消費電力量の一次関数や消費電力量から所定量減ずることにより、消費電力量より低い値に設定されている。
図２（ａ）に観察されるように、消費電力量の変化は、相対的に緩やかに変化する曲線に短周期で変化する成分が重畳された挙動を示す。
図２（ａ）に示す例の場合、逆潮流を防止するため、いずれの時間においても指令値は消費電力量より低い値に設定されている。
しかし、現実にＰＣＳ３が太陽電池２の発電量を制御する場合には、有限の制御遅れ時間（応答時間）が存在し、消費電力量が計測された時刻と、太陽電池２の発電量制御が完了する時刻との間に、有限の時間経過が生じ、例えば数秒程度の時間的なずれ（タイムラグ）が生じる。
しかし、現実にＰＣＳ３が太陽電池２の発電量を制御する場合には、有限の制御遅れ時間（応答時間）が存在し、消費電力量が計測された時刻と、太陽電池２の発電量制御が完了する時刻との間に、有限の時間経過が生じ、例えば数秒程度の時間的なずれ（タイムラグ）が生じる。
例えば、消費電力量を時間ｔの関数ＣＷ（ｔ）とし、時刻ｔ１において、消費電力量ＣＷ（ｔ１）に基づいて、発電量の指令値ＰＷ（ｔ１）（＜ＣＷ（ｔ１））を確定し、ＰＣＳ３が太陽電池２の発電量がＰＷ（ｔ１）となるよう指令したとする。そして、有限の時間（Δｔ）経過後に太陽電池２の発電量がＰＷ（ｔ１）となるとする。
消費電力量は、時々刻々変化するため、時刻ｔ１+ΔｔにおいてはＣＷ１（ｔ１+Δｔ）となる。もし消費電力量が減少した場合、ＣＷ（ｔ１+Δｔ）＜ＣＷ（ｔ１）となる。
この場合、ＰＷ（ｔ１）＜ＣＷ（ｔ１）であるが、消費電力量の減少量ＣＷ（ｔ１）−ＣＷ（ｔ１+Δｔ）が、消費電力量と指令値とのマージンＣＷ（ｔ１）−ＰＷ（ｔ１）を上回ると、ＣＷ（ｔ１）−ＣＷ（ｔ１+Δｔ）＞ＣＷ（ｔ１）−ＰＷ（ｔ１）となる。制御遅れ時間により、時刻ｔ１＋Δｔにおいて太陽電池２の発電量は、ＰＷ（ｔ１）となり、消費電力量より大きくなる（ＰＷ（ｔ１）＞ＣＷ（ｔ１+Δｔ））。その結果、時刻ｔ１＋Δｔにおいて、逆潮流が発生することになる。
消費電力量は、時々刻々変化するため、時刻ｔ１+ΔｔにおいてはＣＷ１（ｔ１+Δｔ）となる。もし消費電力量が減少した場合、ＣＷ（ｔ１+Δｔ）＜ＣＷ（ｔ１）となる。
この場合、ＰＷ（ｔ１）＜ＣＷ（ｔ１）であるが、消費電力量の減少量ＣＷ（ｔ１）−ＣＷ（ｔ１+Δｔ）が、消費電力量と指令値とのマージンＣＷ（ｔ１）−ＰＷ（ｔ１）を上回ると、ＣＷ（ｔ１）−ＣＷ（ｔ１+Δｔ）＞ＣＷ（ｔ１）−ＰＷ（ｔ１）となる。制御遅れ時間により、時刻ｔ１＋Δｔにおいて太陽電池２の発電量は、ＰＷ（ｔ１）となり、消費電力量より大きくなる（ＰＷ（ｔ１）＞ＣＷ（ｔ１+Δｔ））。その結果、時刻ｔ１＋Δｔにおいて、逆潮流が発生することになる。
図２（ｂ）は、図２（ａ）の円で囲んだ領域近傍を拡大したグラフである。円で囲んだ領域は、消費電力量が減少する谷の部分が存在する。図２（ｂ）に示すように、時間Δｔのずれ（制御遅れ時間）が生じた場合には、上記のように指令値ＰＷ（ｔ１）が、消費電力量ＣＷ（ｔ１+Δｔ）を上回り、逆潮流が生じることになる。
このように、発電制御の制御遅れ時間と同程度の短周期での消費電力量の変化は、逆潮流を引き起こすリスクを高める。
このように、発電制御の制御遅れ時間と同程度の短周期での消費電力量の変化は、逆潮流を引き起こすリスクを高める。
そのため従来は、図２（ｃ）に示すように、例えば、１日の内の全ての時刻での消費電力量の短周期の変動幅の最大量（図中ΔＰ）を予め推定しておき、この最大量の変動が生じても逆潮流が生じないよう十分に大きなマージンを持たせて、ＰＣＳ３の指令値を設定する必要があった。
この場合、日によって消費電力量の挙動が変化し、変動の最大量も変わることがあるため、最大量は過去の実績等から推定する必要がある。
また、この最大量の変動が生じる時刻は不明であるため、全ての時刻に対して、この最大量の変動が生じても逆潮流が生じないよう、ＰＣＳ３の指令値を設定する必要がある。
従って、逆潮流のリスクを防止するため、多くの時刻において、消費電力量の短周期の変動は過大評価されることとなる。
太陽電池２の発電量は、このように設定されたマージンのために過剰に制限（低減）されてしまい、結果的に、太陽電池２の発電効率を低下させてしまうことがある。
この場合、日によって消費電力量の挙動が変化し、変動の最大量も変わることがあるため、最大量は過去の実績等から推定する必要がある。
また、この最大量の変動が生じる時刻は不明であるため、全ての時刻に対して、この最大量の変動が生じても逆潮流が生じないよう、ＰＣＳ３の指令値を設定する必要がある。
従って、逆潮流のリスクを防止するため、多くの時刻において、消費電力量の短周期の変動は過大評価されることとなる。
太陽電池２の発電量は、このように設定されたマージンのために過剰に制限（低減）されてしまい、結果的に、太陽電池２の発電効率を低下させてしまうことがある。
なお、実際には、発電制御の有限の制御遅れ時間は、（１）ＰＣＳ３自体の制御特性の他、（２）消費電力量の計測（電力計測器から制御端末への通信）、（３）指令値の計算（ただし、これは一瞬で完了するため無視できる）、（４）指令値の送信（計測端末からＰＣＳ３への通信）が合わさった値であるが、実測により求めることもできる。
従って、短周期の消費電力量の変化を予測し、太陽電池２の発電量の上限を決定することで、逆潮流を防止しつつ太陽電池２の発電効率を、従来より向上させることができる。
すなわち、制御遅れ時間に相当する時間の経過後の消費電力量を、直近の消費電力量のデータに対して多項式を用いた回帰分析により予測し、太陽電池２の発電量の上限を消費電力量の予測値以下に設定する。
具体的には、過去の数点の消費電力量の値から、将来の消費電力量を予測する。推定する時刻は、制御遅れ時間相当の短時間後の時刻であり、例えば制御遅れ時間の２倍から１０倍の時間が経過した時刻の値を予測する。このように短周期での消費電力量の予測をするものであり、従って、長周期の消費電力量の予測は不要である。
なお、上記のように、制御遅れ時間は、消費電力量の計測開始時から発電量の出力までであり、実測も可能である。
すなわち、制御遅れ時間に相当する時間の経過後の消費電力量を、直近の消費電力量のデータに対して多項式を用いた回帰分析により予測し、太陽電池２の発電量の上限を消費電力量の予測値以下に設定する。
具体的には、過去の数点の消費電力量の値から、将来の消費電力量を予測する。推定する時刻は、制御遅れ時間相当の短時間後の時刻であり、例えば制御遅れ時間の２倍から１０倍の時間が経過した時刻の値を予測する。このように短周期での消費電力量の予測をするものであり、従って、長周期の消費電力量の予測は不要である。
なお、上記のように、制御遅れ時間は、消費電力量の計測開始時から発電量の出力までであり、実測も可能である。
このように短周期での消費電力量の変動を予測して追随するため、従来のように、全時刻での短周期の変動幅の最大量を見込んで、消費電力量と指令値とのマージンを設定する必要がなくなり、予測値と指令値とのマージンは、変動幅の最大量（ΔＰ）より小さくすることができる。その結果、発電効率が向上する。
（消費電力量の予測１）
基本的に、その時点の消費電力量の動き、速度と加速度から次の値を予測する。
具体的には、直近の３点の時刻の消費電力量のデータより、消費電力量の時間に関しての回帰分析により以下の二次関数（以下、二次回帰関数と称す）を求める。
Ｙ＝Ａ＋ＢＸ＋ＣＸ２ （式１）
ここでＹは消費電力量、Ｘは時間である。
３点の時刻と消費電力量の計測値、すなわち時刻Ｘ１、Ｘ２、Ｘ３の消費電力量の値Ｙ１、Ｙ２、Ｙ３が得られれば、Ａ、Ｂ、Ｃの値は一意に求めることができる。
上記回帰関数により、時刻Ｘにおける消費電力量Ｙを算出することができる。
なお、上記のように短周期の変化を求めるものであり、予測する時刻Ｘは、制御遅れ時間に相当する時間程度経過した時刻である。なお、制御遅れ時間は、システム全体やＰＣＳ３により変動するため、時刻Ｘは可変に設定できる。
なお、二次方程式におけるＡ、Ｂ、Ｃを求めることは速度、加速度を求める事と同等であることは言うまでもない。
基本的に、その時点の消費電力量の動き、速度と加速度から次の値を予測する。
具体的には、直近の３点の時刻の消費電力量のデータより、消費電力量の時間に関しての回帰分析により以下の二次関数（以下、二次回帰関数と称す）を求める。
Ｙ＝Ａ＋ＢＸ＋ＣＸ２ （式１）
ここでＹは消費電力量、Ｘは時間である。
３点の時刻と消費電力量の計測値、すなわち時刻Ｘ１、Ｘ２、Ｘ３の消費電力量の値Ｙ１、Ｙ２、Ｙ３が得られれば、Ａ、Ｂ、Ｃの値は一意に求めることができる。
上記回帰関数により、時刻Ｘにおける消費電力量Ｙを算出することができる。
なお、上記のように短周期の変化を求めるものであり、予測する時刻Ｘは、制御遅れ時間に相当する時間程度経過した時刻である。なお、制御遅れ時間は、システム全体やＰＣＳ３により変動するため、時刻Ｘは可変に設定できる。
なお、二次方程式におけるＡ、Ｂ、Ｃを求めることは速度、加速度を求める事と同等であることは言うまでもない。
（消費電力量の予測２）
さらに一層、消費電力量データの過去の挙動を反映させて、消費電力量を予測することが可能である。
本実施形態においては、反映させる過去の消費電力量データ数を自動的に変更し、将来の消費電力量を予測することができる。具体的には、３点以上の時間に関する消費電力量の近似式を二次関数（二次回帰関数）として求め、その近似式の実測データに対する相関係数から、近似式の妥当性を判断し、妥当であると判断される最大のデータ数を採用する。
以下、本実施形態２について、詳細に説明する。
さらに一層、消費電力量データの過去の挙動を反映させて、消費電力量を予測することが可能である。
本実施形態においては、反映させる過去の消費電力量データ数を自動的に変更し、将来の消費電力量を予測することができる。具体的には、３点以上の時間に関する消費電力量の近似式を二次関数（二次回帰関数）として求め、その近似式の実測データに対する相関係数から、近似式の妥当性を判断し、妥当であると判断される最大のデータ数を採用する。
以下、本実施形態２について、詳細に説明する。
具体的には、以下の手順に従って、制御装置４の演算処理部は、二次関数（二次回帰関数）を求める。
ステップ０： 相関係数の閾値（Ｒｃ）を設定する。
ステップ１： 直近のｎ＝３点の時刻と消費電力量から二次回帰関数を算出する。
ステップ２： ｎを１増加して（ｎ＝ｎ＋１）、直近のｎ点の時刻と消費電力量から二次回帰関数を算出する。
ステップ３： 得られた二次回帰関数に対して相関係数を算出し、閾値Ｒｃと比較する。
ステップ４： 算出した相関係数が、閾値Ｒｃ以上であれば、処理継続と判断し、ステップ２に戻り、閾値Ｒｃより小さければ、処理完了と判断し、ｎを１減じて（ｎ＝ｎ−１）、直近のｎ点の時刻と消費電力量から算出した二次回帰関数を採用し、手順を終了する。
すなわち、相関係数ｒが閾値Ｒｃ以上であれば、ステップ２、３を繰り返し、相関係数ｒが閾値Ｒｃより小さければステップ４により、手順は終了する。
ステップ０： 相関係数の閾値（Ｒｃ）を設定する。
ステップ１： 直近のｎ＝３点の時刻と消費電力量から二次回帰関数を算出する。
ステップ２： ｎを１増加して（ｎ＝ｎ＋１）、直近のｎ点の時刻と消費電力量から二次回帰関数を算出する。
ステップ３： 得られた二次回帰関数に対して相関係数を算出し、閾値Ｒｃと比較する。
ステップ４： 算出した相関係数が、閾値Ｒｃ以上であれば、処理継続と判断し、ステップ２に戻り、閾値Ｒｃより小さければ、処理完了と判断し、ｎを１減じて（ｎ＝ｎ−１）、直近のｎ点の時刻と消費電力量から算出した二次回帰関数を採用し、手順を終了する。
すなわち、相関係数ｒが閾値Ｒｃ以上であれば、ステップ２、３を繰り返し、相関係数ｒが閾値Ｒｃより小さければステップ４により、手順は終了する。
上記ステップにより、制御装置４は、消費電力量の予測値を算出することができる。その後、その消費電力量の予測値に基づき、太陽電池２の発電量の上限値を算出し、ＰＣＳ３に指令値として出力し、ＰＣＳ３は、その上限値以下となるよう太陽電池２の発電量を制御する。
二次回帰関数（式１）は、直近の連続した計測時刻Ｘｉを用い、時刻Ｘｉのときの消費電力量をＹｉとし、データ数をｎ（自然数）とすると、
Ｃ＝（Ｓｘ２ｙＳｘｘ−ＳｘｙＳｘｘ２）／（ＳｘｘＳｘ２ｘ２−（Ｓｘｘ２）２）
Ｂ＝（ＳｘｙＳｘ２ｘ２−Ｓｘ２ｙＳｘｘ２）／（ＳｘｘＳｘ２ｘ２−（Ｓｘｘ２）２）
Ａ＝＜Ｙ＞−Ｂ＜Ｘ＞−Ｃ＜Ｘ２＞
となる。
ここで、
＜Ｘ＞＝ΣＸｉ／ｎ
＜ｙ＞＝ΣＹｉ／ｎ
＜Ｘ２＞＝ΣＸｉ ２／ｎ
Ｓｘｘ＝Σ（Ｘｉ−＜Ｘ＞）２／ｎ
＝ΣＸｉ ２／ｎ−＜Ｘ＞２
Ｓｘｙ＝Σ（Ｘｉ−＜Ｘ＞）（Ｙｉ−＜Ｙ＞）／ｎ
＝ΣＸｉＹｉ／ｎ−＜Ｘ＞＜Ｙ＞
Ｓｘｘ２＝Σ（Ｘｉ−＜Ｘ＞）（Ｘｉ ２−＜Ｘ２＞）／ｎ
＝ΣＸｉ ３／ｎ−＜Ｘ＞＜Ｘ２＞
Ｓｘ２ｘ２＝Σ（Ｘｉ ２−＜Ｘ２＞）２／ｎ
＝ΣＸｉ ４／ｎ−＜Ｘ２＞＜Ｘ２＞
Ｓｘ２ｙ＝Σ（Ｘｉ ２−＜Ｘ２＞）（Ｙｉ−＜Ｙ＞）／ｎ
＝ΣＸｉ ２Ｙｉ／ｎ−＜Ｘ２＞＜Ｙ＞
である。
なお、Σは、インデックス（ｉ）が１からｎまでの総和を示す。
また、相関係数ｒは、
ｒ＝√（１−Σ（Ｙｉ−（Ａ＋ＢＸｉ＋ＣＸｉ ２））２／（Σ（Ｙｉ−＜Ｙ＞）２））
で定義した。
なお、「√」は平方根を示す。
Ｃ＝（Ｓｘ２ｙＳｘｘ−ＳｘｙＳｘｘ２）／（ＳｘｘＳｘ２ｘ２−（Ｓｘｘ２）２）
Ｂ＝（ＳｘｙＳｘ２ｘ２−Ｓｘ２ｙＳｘｘ２）／（ＳｘｘＳｘ２ｘ２−（Ｓｘｘ２）２）
Ａ＝＜Ｙ＞−Ｂ＜Ｘ＞−Ｃ＜Ｘ２＞
となる。
ここで、
＜Ｘ＞＝ΣＸｉ／ｎ
＜ｙ＞＝ΣＹｉ／ｎ
＜Ｘ２＞＝ΣＸｉ ２／ｎ
Ｓｘｘ＝Σ（Ｘｉ−＜Ｘ＞）２／ｎ
＝ΣＸｉ ２／ｎ−＜Ｘ＞２
Ｓｘｙ＝Σ（Ｘｉ−＜Ｘ＞）（Ｙｉ−＜Ｙ＞）／ｎ
＝ΣＸｉＹｉ／ｎ−＜Ｘ＞＜Ｙ＞
Ｓｘｘ２＝Σ（Ｘｉ−＜Ｘ＞）（Ｘｉ ２−＜Ｘ２＞）／ｎ
＝ΣＸｉ ３／ｎ−＜Ｘ＞＜Ｘ２＞
Ｓｘ２ｘ２＝Σ（Ｘｉ ２−＜Ｘ２＞）２／ｎ
＝ΣＸｉ ４／ｎ−＜Ｘ２＞＜Ｘ２＞
Ｓｘ２ｙ＝Σ（Ｘｉ ２−＜Ｘ２＞）（Ｙｉ−＜Ｙ＞）／ｎ
＝ΣＸｉ ２Ｙｉ／ｎ−＜Ｘ２＞＜Ｙ＞
である。
なお、Σは、インデックス（ｉ）が１からｎまでの総和を示す。
また、相関係数ｒは、
ｒ＝√（１−Σ（Ｙｉ−（Ａ＋ＢＸｉ＋ＣＸｉ ２））２／（Σ（Ｙｉ−＜Ｙ＞）２））
で定義した。
なお、「√」は平方根を示す。
上記のように四則演算により容易に二次回帰関数及び相関係数を求めることができ、演算処理装置、例えばマイコン等により、上記手順を実行することができる。
従って、公知の電力計により消費電力量値を定期的に計測し、得られた計測値を順次記憶装置に保存し、演算処理装置は、記憶装置に保存した消費電力量値を、読み出して算出することができる。
なお、消費電力量は、上記に限らず公知の方法により取得すればよい。
従って、公知の電力計により消費電力量値を定期的に計測し、得られた計測値を順次記憶装置に保存し、演算処理装置は、記憶装置に保存した消費電力量値を、読み出して算出することができる。
なお、消費電力量は、上記に限らず公知の方法により取得すればよい。
以下図３を参照し、消費電力量の予測値を算出する具体例を示す。
以降の消費電力量の予測値を算出する各ステップにおいては、記憶装置に保存されているデータを、最新の現在の計測値から順に、過去の時刻の計測値を、演算処理装置が読み出すこととなる。
例えば、直近の５点の時刻Ｘ［秒］と消費電力量Ｙ［ｋＷ］の組合わせ（Ｘｉ、Ｙｉ）が、（０、６）、（−１、７）、（−２、５）、（−３、４）、（−４、５）であったとする。
なお、時刻Ｘは、１秒毎に測定を行い、現時点を０秒とし、過去の時刻を負の値で表現している。
なお、「直近」とは、最新の時刻とそれに連続した過去の時刻を意味し、上記例のように、各ステップで回帰分析に使用する計測値は、そのステップ実行時刻に対して直近の時刻の計測値である。
以降の消費電力量の予測値を算出する各ステップにおいては、記憶装置に保存されているデータを、最新の現在の計測値から順に、過去の時刻の計測値を、演算処理装置が読み出すこととなる。
例えば、直近の５点の時刻Ｘ［秒］と消費電力量Ｙ［ｋＷ］の組合わせ（Ｘｉ、Ｙｉ）が、（０、６）、（−１、７）、（−２、５）、（−３、４）、（−４、５）であったとする。
なお、時刻Ｘは、１秒毎に測定を行い、現時点を０秒とし、過去の時刻を負の値で表現している。
なお、「直近」とは、最新の時刻とそれに連続した過去の時刻を意味し、上記例のように、各ステップで回帰分析に使用する計測値は、そのステップ実行時刻に対して直近の時刻の計測値である。
ステップ０：
相関係数の閾値（判定値）Ｒｃを０．８に設定する。
このステップは、一度実行すればよく、以下のステップと異なり、各時刻において、消費電力の予測値を算出する毎に実行する必要はない。そのため、予め閾値（判定値）Ｒｃを記憶装置に保存しておけばよい。
以下、ステップ１から４は、太陽電池２を制御する各時刻に対して実行する。
相関係数の閾値（判定値）Ｒｃを０．８に設定する。
このステップは、一度実行すればよく、以下のステップと異なり、各時刻において、消費電力の予測値を算出する毎に実行する必要はない。そのため、予め閾値（判定値）Ｒｃを記憶装置に保存しておけばよい。
以下、ステップ１から４は、太陽電池２を制御する各時刻に対して実行する。
ステップ１：
図３（ａ）に示すように、直近のｎ＝３点の時刻、すなわち時刻Ｘ＝０、−１、−２［秒］における各消費電力量の値を用い、二次回帰関数を求める。
二次関数の各係数は、Ａ＝６、Ｂ＝−２．５、Ｃ＝−１．５となる。
例えば、１秒後（Ｘ＝１）の消費電力量の予測値は、Ｙ＝２．０［ｋＷ］となる。
なお、この場合、言うまでもなくｒ＝１である。
図３（ａ）に示すように、直近のｎ＝３点の時刻、すなわち時刻Ｘ＝０、−１、−２［秒］における各消費電力量の値を用い、二次回帰関数を求める。
二次関数の各係数は、Ａ＝６、Ｂ＝−２．５、Ｃ＝−１．５となる。
例えば、１秒後（Ｘ＝１）の消費電力量の予測値は、Ｙ＝２．０［ｋＷ］となる。
なお、この場合、言うまでもなくｒ＝１である。
ステップ２：
図３（ｂ）に示すように、直近のｎ＋１＝４点の時刻、すなわち時刻Ｘ＝０、−１、−２、−３［秒］における各消費電力量の値を用い、二次回帰関数を求める。
二次関数の各係数は、Ａ＝６．２、Ｂ＝−０．７、Ｃ＝−０．５となる。
例えば、１秒後（Ｘ＝１）の消費電力量の予測値はＹ＝４．０［ｋＷ］となる。
図３（ｂ）に示すように、直近のｎ＋１＝４点の時刻、すなわち時刻Ｘ＝０、−１、−２、−３［秒］における各消費電力量の値を用い、二次回帰関数を求める。
二次関数の各係数は、Ａ＝６．２、Ｂ＝−０．７、Ｃ＝−０．５となる。
例えば、１秒後（Ｘ＝１）の消費電力量の予測値はＹ＝４．０［ｋＷ］となる。
ステップ３：
得られた二次回帰関数に対して相関係数ｒを算出すると、ｒ＝０．９２となる。
得られた二次回帰関数に対して相関係数ｒを算出すると、ｒ＝０．９２となる。
ステップ４：
算出した相関係数ｒ（０．９２）が、閾値Ｒｃ（０．８）以上であるため、ステップ２に戻る。
算出した相関係数ｒ（０．９２）が、閾値Ｒｃ（０．８）以上であるため、ステップ２に戻る。
ステップ２：
図３（ｃ）に示すように、直近のｎ＋１＝５点の時刻、すなわち時刻Ｘ＝０、−１、−２、−３、−４［秒］における各消費電力量の値を用い、二次回帰関数を算出する。
二次関数の各係数は、Ａ＝６．５４、Ｂ＝−０．７９、Ｃ＝０．０７となる。
例えば、１秒後（Ｘ＝１）の消費電力量の予測値はＹ＝４．０［ｋＷ］となる。
図３（ｃ）に示すように、直近のｎ＋１＝５点の時刻、すなわち時刻Ｘ＝０、−１、−２、−３、−４［秒］における各消費電力量の値を用い、二次回帰関数を算出する。
二次関数の各係数は、Ａ＝６．５４、Ｂ＝−０．７９、Ｃ＝０．０７となる。
例えば、１秒後（Ｘ＝１）の消費電力量の予測値はＹ＝４．０［ｋＷ］となる。
ステップ３：
得られた二次回帰関数に対して相関係数ｒを算出すると、ｒ＝０．７０となる。
得られた二次回帰関数に対して相関係数ｒを算出すると、ｒ＝０．７０となる。
ステップ４：
算出した相関係数ｒ（０．７０）が、閾値Ｒｃ（０．８）より小さいため、直近のｎ−１＝４点の時刻に対する二次回帰関数を採用し、手順を終了する。
算出した相関係数ｒ（０．７０）が、閾値Ｒｃ（０．８）より小さいため、直近のｎ−１＝４点の時刻に対する二次回帰関数を採用し、手順を終了する。
従って、採用された二次回帰関数は、直近のｎ＝４点の実測データから得られた二次関数を採用し、二次関数の各係数は、Ａ＝６．２、Ｂ＝−０．７、Ｃ＝−０．５となる。
消費電力量の予測値はＹ＝４．０［ｋＷ］となる。
消費電力量の予測値はＹ＝４．０［ｋＷ］となる。
このように、複数の直近の過去及び現在の時刻と消費電力量との相関関係から、所定時間後の消費電力量の予測値を算出できる。ここで、所定時間とは、ＰＣＳ３の制御遅れ時間程度に設定されるが、制御遅れ時間の変動を考慮すると、好適には、例えば制御遅れ時間の２倍から１０倍である。
制御する太陽電池２の発電量の上限を、消費電力量の予測値以下に設定する。
例えば、消費電力量の予測値より一定数少ない値又は消費電力量の予測値の関数（比例、一次関数等）とし、例えば３．５［ｋＷ］とする。
例えば、消費電力量の予測値より一定数少ない値又は消費電力量の予測値の関数（比例、一次関数等）とし、例えば３．５［ｋＷ］とする。
制御装置４により求められた太陽電池２の発電量の上限は、ＰＣＳ３に出力され、ＰＣＳ３は、太陽電池２の発電量を制御する。消費電力量の予測値は、制御遅れ時間に相当する時間が経過した時刻の値であるため、有限の制御遅れ時間が経過した時刻において、消費電力量の予測値を上回らないように太陽電池２は制御されることとなる。
なお、上記のステップ０から４により、複数の予測値が得られる場合、最も多くのデータ点数から予測値を採用したが、上記ステップ０から４により求めた各二次回帰関数から算出された予測値のうち、最も低い値となる予測値を採用してもよい。
（第２の実施形態）
第１の実施形態においては、取得された消費電力量から直接回帰分析により消費電力量の予測値を算出した。本実施形態においては、取得された消費電力量の移動平均を算出し、算出された移動平均から回帰分析により消費電力量を算出する。
取得された消費電量の時間推移の短周期の変動が激しい場合、消費電力量の予測値が過剰に大きくなるか又は小さくなる、すなわち消費電力量を過大評価又は過小評価する場合がある。消費電力量を過大評価した場合、逆潮流のリスクが増大し、消費電力を過小評価した場合、太陽電池２の発電効率を低下することになる。
第１の実施形態においては、取得された消費電力量から直接回帰分析により消費電力量の予測値を算出した。本実施形態においては、取得された消費電力量の移動平均を算出し、算出された移動平均から回帰分析により消費電力量を算出する。
取得された消費電量の時間推移の短周期の変動が激しい場合、消費電力量の予測値が過剰に大きくなるか又は小さくなる、すなわち消費電力量を過大評価又は過小評価する場合がある。消費電力量を過大評価した場合、逆潮流のリスクが増大し、消費電力を過小評価した場合、太陽電池２の発電効率を低下することになる。
以下に示すように、消費電力量の移動平均に基づいて消費電力量の予測値を算出することで、過大評価又は過小評価を低減することができる。
図４は、第１の実施形態及び第２の実施形態により消費電力量の予測値の時間的推移例を比較して示す。
点線Ａは、実際に取得された消費電力量、破線Ｂは、消費電力量の値から算出した予測値、実線Ｃは、消費電力量の移動平均から回帰分析により算出した予測値を示す。
この場合、破線Ｂ及び実線Ｃは、それぞれ直近の３点のデータから回帰分析により２次関数を求め、予測値を算出した。
図４は、第１の実施形態及び第２の実施形態により消費電力量の予測値の時間的推移例を比較して示す。
点線Ａは、実際に取得された消費電力量、破線Ｂは、消費電力量の値から算出した予測値、実線Ｃは、消費電力量の移動平均から回帰分析により算出した予測値を示す。
この場合、破線Ｂ及び実線Ｃは、それぞれ直近の３点のデータから回帰分析により２次関数を求め、予測値を算出した。
破線Ｂは、点線Ａと比べ、短周期の変動量が大きいことが分かる。一方、実線Ｃは、破線Ｂと比べ、短周期の変動が抑えられ、過大評価又は過小評価が抑制されていることが分かる。
以下具体的に、制御装置４による、移動平均を用いた消費電力量の予測方法及び太陽電池２の発電量の制御方法について説明する。
なお、移動平均とは、時系列データにおいて直近の複数のデータを用いて平均値を計算する手法であるが、簡単のため、移動平均によって得られた平均値を「移動平均」と称することがある。
以下具体的に、制御装置４による、移動平均を用いた消費電力量の予測方法及び太陽電池２の発電量の制御方法について説明する。
なお、移動平均とは、時系列データにおいて直近の複数のデータを用いて平均値を計算する手法であるが、簡単のため、移動平均によって得られた平均値を「移動平均」と称することがある。
まず、消費電力量の移動平均を算出するための時間間隔を決定する。これは、その時間間隔に相当するデータ点数（以下このデータ点数をｋとする。）を決定することと同じである。例えば、実際に得られた消費電力量から、短周期で変動する平均的周期を求め、その０．５〜１倍の時間間隔とすることができる。
また、ＰＣＳ３のＭＰＰＴ法による発電量が最大となる条件の探索に要する時間（以下探索時間と称す）に相当する時間間隔で平均を行ってもよい。探索時間は、実際に使用しているＰＣＳ３に依存し、実測が可能である。特に探索時間に比べ消費電量の変動の周期が短い場合、探索時間に相当する時間（例えば探索時間の１倍から２倍の時間）の間隔で、消費電力量の移動平均を算出してもよい。
また、ＰＣＳ３のＭＰＰＴ法による発電量が最大となる条件の探索に要する時間（以下探索時間と称す）に相当する時間間隔で平均を行ってもよい。探索時間は、実際に使用しているＰＣＳ３に依存し、実測が可能である。特に探索時間に比べ消費電量の変動の周期が短い場合、探索時間に相当する時間（例えば探索時間の１倍から２倍の時間）の間隔で、消費電力量の移動平均を算出してもよい。
各時刻において、得られた消費電力量に対して、直近のｋ個、例えば３〜５個の消費電力量の移動平均を算出する。
なお、移動平均は、ｋ個の単純移動平均であってもよいが、ｋ個のデータに対して、それぞれに重みを付けた加重移動平均であってもよい。加重移動平均の場合、最新の消費電力量の重みを大きくし、線形加重移動平均や指数加重移動平均等を用いてもよい。
なお、移動平均は、ｋ個の単純移動平均であってもよいが、ｋ個のデータに対して、それぞれに重みを付けた加重移動平均であってもよい。加重移動平均の場合、最新の消費電力量の重みを大きくし、線形加重移動平均や指数加重移動平均等を用いてもよい。
このようにして、各時刻において移動平均を求めておき、記憶装置に保存する。
その後、記憶装置に保存した消費電力量の移動平均に対して、（消費電力量の予測１）又は（消費電力量の予測２）に記載の方法により、回帰分析に基づいて近似関数を導出し、得られた近似関数から予測値を算出する。
得られた予測値に対して、太陽電池２の発電量の上限を決定することができる。
その後、記憶装置に保存した消費電力量の移動平均に対して、（消費電力量の予測１）又は（消費電力量の予測２）に記載の方法により、回帰分析に基づいて近似関数を導出し、得られた近似関数から予測値を算出する。
得られた予測値に対して、太陽電池２の発電量の上限を決定することができる。
なお、太陽電池２の発電量の上限を決定するために使用する消費電力量の値として、得られた実際の消費電力量が増加傾向にある場合、第１の実施形態により算出した消費電力量の予測値は、過大評価となるリスクがあるため、第２の実施形態により算出した消費電力量の予測値を使用してもよい。
また、得られた実際の消費電力量が減少傾向にある場合、第１の実施形態により算出した消費電力量の予測値は過小評価の傾向があるが、逆潮流のリスク低減を優先して、第１の実施形態により算出した消費電力量の予測値を使用してもよい。
また、得られた実際の消費電力量、第１の実施形態により算出した消費電力量の予測値、第２の実施形態により算出した消費電力量の予測値を比較し、逆潮流のリスク低減を優先して、最も低い値を予測値として使用してもよい。
また、得られた実際の消費電力量が減少傾向にある場合、第１の実施形態により算出した消費電力量の予測値は過小評価の傾向があるが、逆潮流のリスク低減を優先して、第１の実施形態により算出した消費電力量の予測値を使用してもよい。
また、得られた実際の消費電力量、第１の実施形態により算出した消費電力量の予測値、第２の実施形態により算出した消費電力量の予測値を比較し、逆潮流のリスク低減を優先して、最も低い値を予測値として使用してもよい。
このようにして得られた太陽電池２の発電量の上限に基づき、制御装置４は、ＰＣＳ３に対して発電量の制御値を出力し、ＰＣＳ３が、太陽電池２の発電量を制御することができる。
なお、上記の各実施形態において、消費電力量の予測のために二次関数を用いたが、さらに高次の多項式を用いてもよい。例えば、ｍ次多項式を用いた回帰分析の場合、ステップ１においてｎを３ではなくｍ＋１とし、二次回帰関数をｍ次多項式とすればよい。この場合、最小二乗法によりｍ次多項式の係数を代数的に求めることができるため、四則演算可能なマイコン等によって容易にｍ次多項式を算出できる。
なお、さらに三角関数や他の関数（例えば直交関数）やガウス関数等を用いて予測値を計算してもよい。ただし、上記のようにｍ次多項式を用いることで、代数的にｍ次多項式の係数を求めることができ、高度な演算処理が可能なＣＰＵを必要とせず、安価なマイコンによってもシステムを構成することができる。
なお、さらに三角関数や他の関数（例えば直交関数）やガウス関数等を用いて予測値を計算してもよい。ただし、上記のようにｍ次多項式を用いることで、代数的にｍ次多項式の係数を求めることができ、高度な演算処理が可能なＣＰＵを必要とせず、安価なマイコンによってもシステムを構成することができる。
以上のように、太陽電池２の発電量の制御遅れ時間を考慮し、直近の消費電力量のデータから、所定時間経過後の消費電力量を予測値を算出し、算出された消費電力量を予測値を下回るように太陽電池２の発電量が制御される。そのため、消費電力量の短周期の変動の最大量（ΔＰ）を見積もる必要もなく、全ての時刻に対して過大にマージンを持たせてＰＣＳ３の指令値を設定することがない。
そのため、太陽電池２の発電量の上限を不要に低く設定することを防止でき、その結果、太陽電池２の発電量の効率を、向上させることができる。
そのため、太陽電池２の発電量の上限を不要に低く設定することを防止でき、その結果、太陽電池２の発電量の効率を、向上させることができる。
また、既存の太陽電池２及びＰＣＳ３に、上記制御を実行する制御装置４を組み込むことが可能であり、既存の太陽光発電システムを有効に活用し、発電効率を向上させることが可能である。
本発明によれば、逆潮流を防止しつつ太陽電池の発電効率を向上させることができ、特に既存の太陽光発電システムの発電効率を向上させることも可能であり、産業上の利用可能性は大きい。
１ 太陽光発電制御システム
２ 太陽電池
３ ＰＣＳ（パワーコンディショナ）
４ 発電制御装置
５ 受変電部（配電盤）
６ 負荷
７ 電力会社（外部電源）
２ 太陽電池
３ ＰＣＳ（パワーコンディショナ）
４ 発電制御装置
５ 受変電部（配電盤）
６ 負荷
７ 電力会社（外部電源）
Claims (9)
- 太陽電池と、
パワーコンディショナーと、
制御装置とを備え、
前記制御装置は、消費電力量を取得し、回帰分析に基づいて得られた関数により複数の時刻の消費電力量から所定時間経過後の消費電力量の予測値を算出し、
前記予測値に基づき指令値を算出して、前記パワーコンディショナーに出力し、
前記パワーコンディショナーは、前記指令値に基づき、前記太陽電池の発電量を前記予測値以下に制御することを特徴とする太陽光発電システム。 - 前記関数は、ｍ次多項式（ｍ自然数）であることを特徴とする請求項１記載の太陽光発電システム。
- 前記関数は、２次多項式であることを特徴とする請求項１又は２記載の太陽光発電システム。
- 自然数ｎを３として、ｎ点の時刻と消費電力量の回帰分析から二次関数を求める第１のステップと、
ｎを１増加して、ｎ点の時刻と消費電力量の回帰分析から二次関数と相関係数を求める第２のステップと、
予め設定された閾値と前記相関係数とを比較する第３ステップとをこの順に実行し、
前記第３のステップにおいて前記相関係数が前記閾値以上であれば、前記第２のステップを実行し、前記第３のステップにおいて、前記相関係数が前記閾値より小さければ、ｎから１を減じたｎ点の時刻と消費電力量についての回帰分析により求めた二次関数から、消費電力量の前記予測値を算出することを特徴とする請求項３記載の太陽光発電システム。 - 前記所定時間は、制御遅れ時間の２倍から１０倍であることを特徴とする請求項１乃至４のいずれか１項記載の太陽光発電システム。
- 取得した消費電力量の代わりに、直近の消費電力量の移動平均を用いることを特徴とする請求項１乃至５のいずれか１項記載の太陽光発電システム。
- 太陽電池と、
パワーコンディショナーとを備えた太陽光発電システムにおいて、
複数の時刻の取得された消費電力量から、所定時間経過後の消費電力量の予測値を算出する第１の工程と、
前記予測値を基に前記パワーコンディショナーに対する指令値を算出する第２の工程と、
前記パワーコンディショナーが前記指令値に従って前記太陽電池の発電量を前記予測値以下に制御する第３の工程とを有することを特徴とする発電制御方法。 - 前記第１の工程は、
自然数ｎを３として、ｎ点の時刻と消費電力量の回帰分析から二次関数を求める第１のステップと、
ｎを１増加して、ｎ点の時刻と消費電力量の回帰分析から二次関数と相関係数を求める第２のステップと、
予め設定された閾値と前記相関係数とを比較する第３ステップとをこの順に実行し、
前記第３のステップにおいて前記相関係数が前記閾値以上であれば、前記第２のステップを実行し、前記第３のステップにおいて、前記相関係数が前記閾値より小さければ、ｎから１を減じたｎ点の時刻と消費電力量の回帰分析から求めた二次関数から消費電力量の前記予測値を算出することを含むことを特徴とする請求項７記載の発電制御方法。 - 取得した消費電力量の代わりに、直近の消費電力量の移動平均を用いることを特徴とする請求項７又は８記載の発電制御方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019146624A JP6792272B1 (ja) | 2019-08-08 | 2019-08-08 | 発電制御システム及び発電制御方法 |
| JP2020179886A JP7366352B2 (ja) | 2019-08-08 | 2020-10-27 | 発電制御装置、発電制御方法、発電制御プログラム及び発電システム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019146624A JP6792272B1 (ja) | 2019-08-08 | 2019-08-08 | 発電制御システム及び発電制御方法 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020179886A Division JP7366352B2 (ja) | 2019-08-08 | 2020-10-27 | 発電制御装置、発電制御方法、発電制御プログラム及び発電システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP6792272B1 JP6792272B1 (ja) | 2020-11-25 |
| JP2021027775A true JP2021027775A (ja) | 2021-02-22 |
Family
ID=73455249
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019146624A Active JP6792272B1 (ja) | 2019-08-08 | 2019-08-08 | 発電制御システム及び発電制御方法 |
| JP2020179886A Active JP7366352B2 (ja) | 2019-08-08 | 2020-10-27 | 発電制御装置、発電制御方法、発電制御プログラム及び発電システム |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020179886A Active JP7366352B2 (ja) | 2019-08-08 | 2020-10-27 | 発電制御装置、発電制御方法、発電制御プログラム及び発電システム |
Country Status (1)
| Country | Link |
|---|---|
| JP (2) | JP6792272B1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022210768A1 (ja) * | 2021-03-31 | 2022-10-06 | 株式会社ラプラス・システム | 発電システムの制御端末及びその制御用プログラム、及び発電システムの製造方法 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7194767B2 (ja) * | 2021-03-26 | 2022-12-22 | 株式会社コンテック | 発電制御システム |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017038432A (ja) * | 2015-08-07 | 2017-02-16 | シャープ株式会社 | 制御装置、システムおよび制御方法 |
| JP6364567B1 (ja) * | 2018-03-09 | 2018-07-25 | 株式会社ラプラス・システム | 発電制御装置及びそれを用いた発電制御システム |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008066016A (ja) | 2006-09-05 | 2008-03-21 | Ebara Ballard Corp | 燃料電池システムの運転方法及び燃料電池システム |
| JP2011176948A (ja) | 2010-02-24 | 2011-09-08 | Panasonic Corp | 料金管理装置、料金管理方法及び料金管理システム |
| JP5799548B2 (ja) | 2011-03-31 | 2015-10-28 | 株式会社ノーリツ | 発電システム |
| JP6097947B2 (ja) | 2012-11-30 | 2017-03-22 | 株式会社Tbk | 発電システム及びその制御方法 |
-
2019
- 2019-08-08 JP JP2019146624A patent/JP6792272B1/ja active Active
-
2020
- 2020-10-27 JP JP2020179886A patent/JP7366352B2/ja active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017038432A (ja) * | 2015-08-07 | 2017-02-16 | シャープ株式会社 | 制御装置、システムおよび制御方法 |
| JP6364567B1 (ja) * | 2018-03-09 | 2018-07-25 | 株式会社ラプラス・システム | 発電制御装置及びそれを用いた発電制御システム |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022210768A1 (ja) * | 2021-03-31 | 2022-10-06 | 株式会社ラプラス・システム | 発電システムの制御端末及びその制御用プログラム、及び発電システムの製造方法 |
| JP7222502B1 (ja) * | 2021-03-31 | 2023-02-15 | 株式会社ラプラス・システム | 発電システムの制御端末及びその制御用プログラム、及び発電システムの製造方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2021027800A (ja) | 2021-02-22 |
| JP6792272B1 (ja) | 2020-11-25 |
| JP7366352B2 (ja) | 2023-10-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11502534B2 (en) | Electrical energy storage system with battery state-of-charge estimation | |
| JP7266587B2 (ja) | 再充電可能なバッテリを制御するための方法およびシステム | |
| US8527110B2 (en) | Charge/discharge control device and power generation system | |
| JP6792272B1 (ja) | 発電制御システム及び発電制御方法 | |
| KR102067830B1 (ko) | Ｅｓｓ 출력 제어 방법 | |
| JP5507669B2 (ja) | 電力供給システム、電力供給方法および電力供給システムの制御プログラム | |
| JP6255251B2 (ja) | 太陽光発電装置の出力推定方法および装置並びにこれを用いた電力系統監視装置 | |
| CN104636874B (zh) | 检测业务异常的方法及设备 | |
| WO2011040470A1 (ja) | 充放電制御装置および発電システム | |
| JP6548570B2 (ja) | 電力供給システム、電力供給システム用の制御装置およびプログラム | |
| JP2021027800A5 (ja) | ||
| JP2014033539A (ja) | 電力制御装置、方法、プログラム | |
| JP6782573B2 (ja) | 太陽光発電出力推定装置及び太陽光発電出力推定方法 | |
| JP6314543B2 (ja) | 蓄電池の充電状態推定方法及び装置 | |
| JP5709910B2 (ja) | 制御装置及び方法並びにプログラム、それを備えた自然エネルギー発電装置 | |
| JP2023016944A (ja) | 発電制御システム及び発電制御方法 | |
| JP6406391B2 (ja) | 発電制御装置および制御方法 | |
| JP5375948B2 (ja) | 冷却装置を有する電子装置および冷却プログラム | |
| JP6621704B2 (ja) | 発電システム | |
| JP2018143038A (ja) | 太陽光発電システム及び発電制御プログラム | |
| CN112943557B (zh) | 风电场、风力发电机组及其运行状态的预测方法和设备 | |
| JP6464848B2 (ja) | 出力平滑化装置及び出力平滑化方法 | |
| EP2961031B1 (en) | Frequency control method and frequency control system | |
| JP6750472B2 (ja) | 情報処理装置、情報処理システム、情報処理プログラム及び情報処理方法 | |
| WO2023054476A1 (ja) | 発電電力予測方法、発電電力予測装置及び太陽光発電システム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200522 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20200522 |
|
| A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20200817 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200820 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200907 |
|
| 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: 20200929 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20201028 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6792272 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |