(57)【要約】【課題】電圧制御機器の制御遅れを考慮した電圧・無効電力制御システムおよび電圧・無効電力制御方法を提供する。【解決手段】需要予測部４は、変電所で測定された有効電力および無効電力に基づいて、将来の変電所ごとの有効電力および無効電力を予測する。最適化部６は、メタヒューリスティック手法を用いて、予測された将来の変電所ごとの有効電力および無効電力に基づいて、各変電所の設備の最適制御量を探索する。【選択図】図３
2010136600.tif 000002
【特許請求の範囲】
【請求項１】
電力系統の電圧を適正範囲に維持するための制御を行う電圧無効電力制御システムであって、
変電所ごとの有効電力および無効電力を検出する検出部と、
前記測定した有効電力および無効電力に基づいて、将来の変電所ごとの有効電力および無効電力を予測する予測部と、
メタヒューリスティック手法を用いて、前記予測された将来の変電所ごとの有効電力および無効電力に基づいて、各変電所の設備の最適制御量を探索する最適化部と、
前記最適制御量に基づいて、前記各変電所の設備を制御する制御部とを備える、電圧無効電力制御システム。
【請求項２】
前記最適化部は、電力系統の有効電力ロス、変電所の変圧器のタップ位置の変化量、および変電所の調相設備の投入台数の変化量を変数として含む評価関数に基づいて前記最適制御量を探索する、請求項１に記載の電圧無効電力制御システム。
【請求項３】
前記最適化部は、パーティクル・スウォーム・オプティマイゼーション（ＰＳＯ）を用いて、前記最適制御量を探索する、請求項２記載の電圧無効電力制御システム。
【請求項４】
前記最適化部は、パーティクルの位置を表わす変数に、変電所の調相設備の投入台数、および変電所の変圧器のタップ位置を割当て、前記パーティクルの最適な位置の探索を行ない、
前記制御部は、前記最適な位置に基づいて、前記調相設備の投入台数および前記タップ位置を制御する、請求項３記載の電圧無効電力制御システム。
【請求項５】
前記最適化部は、タブーサーチ（ＴＳ）を用いて、前記最適制御量を探索する、請求項２記載の電圧無効電力制御システム。
【請求項６】
前記最適化部は、状態を表わす変数に、変電所の調相設備の投入台数、および変電所の変圧器のタップ位置を割当て、前記評価関数が最適な状態の探索を行ない、
前記制御部は、前記最適な状態に基づいて、前記調相設備の投入台数および前記タップ位置を制御する、請求項５記載の電圧無効電力制御システム。
【請求項７】
前記最適化部は、遺伝的アルゴリズム（ＧＡ）を用いて、前記最適制御量を探索する、請求項２記載の電圧無効電力制御システム。
【請求項８】
前記最適化部は、染色体を構成する遺伝子に、変電所の調相設備の投入台数、および変電所の変圧器のタップ位置を割当て、前記評価関数が最適な染色体の探索を行ない、
前記制御部は、前記最適な染色体に基づいて、前記調相設備の投入台数および前記タップ位置を制御する、請求項７記載の電圧無効電力制御システム。
【請求項９】
前記予測部は、自己回帰モデルに基づいて、前記有効電力の時間差分量および前記無効電力の時間差分量を予測する、請求項１〜８のいずれか１項に記載の電圧無効電力制御システム。
【請求項１０】
電力系統の電圧を適正範囲に維持するための制御を行う電圧無効電力制御方法であって、
変電所ごとの有効電力および無効電力を検出するステップと、
前記検出した有効電力および無効電力に基づいて、将来の変電所ごとの有効電力および無効電力を予測するステップと、
メタヒューリスティック手法を用いて、前記予測された将来の変電所ごとの有効電力および無効電力とに基づいて、各変電所の設備の最適制御量を算出するステップと、
前記最適制御量に基づいて、前記各変電所の設備を制御するステップとを備えた、電圧無効電力制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電圧無効電力制御システムおよび電圧無効電力制御方法に関し、特に、メタヒューリスティック法を用いた電圧無効電力制御システムおよび電圧無効電力制御方法に関する。
【背景技術】
【０００２】
電圧無効電力の自動制御（以下、Voltage Q (reactive power) Control；ＶＱＣ）は、電力系統の安定運用、適正電圧の維持、送電損失の低減機器の動作回数の最小化、および送電系統運用者の負担軽減の目的で導入されている。
【０００３】
近年では、中央ＶＱＣ方式が採用されている。この方式では、各制御装置に対して、中央給電指令所で計算された最適潮流計算の結果に基づく指令を行ない、系統全体の電圧を総合的に制御する。最適潮流計算では、膨大な組合せ最適化問題を解く必要があり、たとえば、特許文献１では、パーティクル・スウォーム・オプティマイゼーション（Particle Swarm Optimization：ＰＳＯ）を用いることが提案されている。特許文献１では、母線電圧の上下限、線路潮流の上限および電圧の安定性を制約条件とし、対象系統の電力損失最小化を目的関数としてＰＳＯによって、最適制御量を探索することが記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０００−１１６００３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１では、電圧制御機器の制御遅れが考慮されていない。すなわち、特許文献１では、上記のように最適制御量を決定した後、制御命令を出して電圧制御機器の制御動作が完了するまでの間に無視できない程度の時間を要する。そのため、決定した最適制御量が、実際に制御が完了する時点では最適なものでなくなってしまうことがある。
【０００６】
それゆえに、本発明の目的は、電圧制御機器の制御遅れを考慮した電圧無効電力制御システムおよび電圧無効電力制御方法を提供することである。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明は、電力系統の電圧を適正範囲に維持するための制御を行う電圧無効電力制御システムであって、変電所ごとの有効電力および無効電力を検出する検出部と、測定した有効電力および無効電力に基づいて、将来の変電所ごとの有効電力および無効電力を予測する予測部と、メタヒューリスティック手法を用いて、予測された将来の変電所ごとの有効電力および無効電力に基づいて、各変電所の設備の最適制御量を探索する最適化部と、最適制御量に基づいて、各変電所の設備を制御する制御部とを備える。
【０００８】
好ましくは、最適化部は、電力系統の有効電力ロス、変電所の変圧器のタップ位置の変化量、および変電所の調相設備の投入台数の変化量を変数として含む評価関数に基づいて最適制御量を探索する。
【０００９】
好ましくは、最適化部は、パーティクル・スウォーム・オプティマイゼーション（ＰＳＯ）を用いて、最適制御量を探索する。
【００１０】
好ましくは、最適化部は、パーティクルの位置を表わす変数に、変電所の調相設備の投入台数、および変電所の変圧器のタップ位置を割当て、パーティクルの最適な位置の探索を行ない、制御部は、最適な位置に基づいて、調相設備の投入台数およびタップ位置を制御する。
【００１１】
好ましくは、最適化部は、タブーサーチ（ＴＳ）を用いて、最適制御量を探索する。
好ましくは、最適化部は、状態を表わす変数に、変電所の調相設備の投入台数、および変電所の変圧器のタップ位置を割当て、評価関数が最適な状態の探索を行ない、制御部は、最適な状態に基づいて、調相設備の投入台数およびタップ位置を制御する。
【００１２】
好ましくは、最適化部は、遺伝的アルゴリズム（ＧＡ）を用いて、最適制御量を探索する。
【００１３】
好ましくは、最適化部は、染色体を構成する遺伝子に、変電所の調相設備の投入台数、および変電所の変圧器のタップ位置を割当て、評価関数が最適な染色体の探索を行ない、制御部は、最適な染色体に基づいて、調相設備の投入台数およびタップ位置を制御する。
【００１４】
好ましくは、予測部は、自己回帰モデルに基づいて、有効電力の時間差分量および無効電力の時間差分量を予測する。
【００１５】
また、本発明は、電力系統の電圧を適正範囲に維持するための制御を行う電圧無効電力制御方法であって、変電所ごとの有効電力および無効電力を検出するステップと、検出した有効電力および無効電力に基づいて、将来の変電所ごとの有効電力および無効電力を予測するステップと、メタヒューリスティック手法を用いて、予測された将来の変電所ごとの有効電力および無効電力とに基づいて、各変電所の設備の最適制御量を算出するステップと、最適制御量に基づいて、各変電所の設備を制御するステップとを備える。
【発明の効果】
【００１６】
本発明によれば、電圧制御機器の制御遅れを考慮した電圧無効電力制御を行なうことができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施形態の電圧無効電力制御システムの構成を表わす図である。
【図２】図１の変電所の構成を表わす図である。
【図３】図１のＶＱＣ装置の構成を表わす図である。
【図４】第１の実施形態の電圧無効電力制御システムの動作手順を表わすフローチャートである。
【図５】第２の実施形態の電圧無効電力制御システムの動作手順を表わすフローチャートである。
【図６】図５のステップＳ２１１の詳細な手順を表わすフローチャートである。
【図７】反復回数ｋの個体ｉの染色体Ｘikの例を表わす図である。
【図８】交差の例を説明するための図である。
【図９】突然変異の例を説明するための図である。
【図１０】第３の実施形態の電圧無効電力制御システムの動作手順を表わすフローチャートである。
【図１１】図１０のステップＳ４１０の詳細な手順を表わすフローチャートである。
【図１２】３つの手法による、ある変電所の母線電圧を表わす図である。
【図１３】３つの手法による、別の変電所の母線電圧を表わす図である。
【図１４】実績値を用いたシミュレーションによる電力系統の有効電力ロスと、３つの手法による電力系統の有効電力ロスを表わす図である。
【発明を実施するための形態】
【００１８】
以下、本発明の実施の形態について、図面を参照して説明する。
［第１の実施形態］
第１の実施形態は、ＰＳＯに基づいて、最適制御量を探索する電圧無効電力制御システムに関する。
【００１９】
（構成）
図１は、本発明の実施形態の電圧無効電力制御システムの構成を表わす図である。
【００２０】
図１を参照して、電圧無効電力制御システムは、変電所２−１〜２−Ｍからなる変電所群５と、発電機１２−１〜１２−Ｈからなる発電機群１１と、変電所２−１〜２−Ｍと無線または有線の通信回線７−１〜７−Ｍを介して接続されるＶＱＣ装置１から構成される。変電所２−１〜２−Ｍの各々は、隣接する変電所と連系線９−１〜９−２を介して接続される。
【００２１】
図２は、図１の変電所ｊ（ｊ＝１〜Ｍ）の構成を表わす図である。
図２を参照して、変電所ｊ（２−ｊ）は、通信部２２と、Ｐ検出部２４と、Ｑ検出部２６と、制御部２８と、タップ付き変圧器３０と、調相設備群３２とを備える。調相設備群３２は、変電所１〜Ｌにおいてのみ含まれるものとする。
【００２２】
タップ付き変圧器３０は、タップを含み、タップの位置を切り替えることで、電力系統の電圧と有効電力ロスが変化する。
【００２３】
変電所ｊ（ｊ＝１〜ＬＳ）において、調相設備群３２は、２０ＭＶＡ調相設備群３３−１と、３０ＭＶＡ調相設備群３３−２と、４０ＭＶＡ調相設備群３３−３とを含む。また、変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）において、調相設備群３２は、２０ＭＶＡ調相設備群３３−１と、３０ＭＶＡ調相設備群３３−２と、４０ＭＶＡ調相設備群３３−３と、６０ＭＶＡ調相設備群３３−４とを含む。
【００２４】
各調相設備群３３−１〜３３−４は、１台または複数台の調相設備を含む。調相設備の投入台数を切り替えることで、電力系統の電圧と有効電力ロスが変化する。
【００２５】
本発明の実施形態では、Ｍ個の変電所すべてにおいて、タップ付き変圧器３０のタップ位置が制御対象となることとし、Ｍ個の変電所のうちＬ個の変電所において、調相設備群３２の投入が制御対象となることとする。
【００２６】
Ｐ検出部２４は、時点ｎにおいて、変電所ｊの有効電力Ｐ（ｎ，ｊ）を検出する。
Ｑ検出部２６は、時点ｎにおいて、変電所ｊの無効電力Ｑ（ｎ，ｊ）を検出する。
【００２７】
制御部２８は、ＶＱＣ装置１から受信した最適制御量に基づいて、タップ付き変圧器３０のタップ位置と、２０ＭＶＡ、３０ＭＶＡ、４０ＭＶＡおよび６０ＭＶＡの調相設備の投入台数を制御する。
【００２８】
通信部２２は、検出した有効電力Ｐ（ｎ，ｊ）および無効電力ＱをＶＱＣ装置１へ通信回線ｊを通じて送信する。また、通信部２２は、ＶＱＣ装置１から通信回線ｊを通じて、変電所ｊの設備の最適制御量である、時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，ｊ）、時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，ｊ）、時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，ｊ）、時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ｊ）、およびタップ付き変圧器３０の時点ｎのタップ位置ＴＡＰ（ｎ，ｊ）を受信する。
【００２９】
図３は、図１のＶＱＣ装置１の構成を表わす図である。
図３を参照して、ＶＱＣ装置１は、通信部３と、需要予測部４と、最適化部６と、潮流計算部８とを備える。ＶＱＣ装置１は、たとえばＣＰＵとメモリを備えたコンピュータで実現できる。この場合、ＣＰＵが、需要予測、最適化計算、および潮流計算を行ない、計算結果および計算の途中結果をメモリに記憶する。
【００３０】
（需要予測部）
需要予測部４は、自己回帰モデル（Autoregressive model：ＡＲモデル)に基づいて、変電所ごとに、次の時点の有効電力予測量および無効電力予測量を算出する。
【００３１】
まず、有効電力予測量の自己回帰モデルについて説明する。
変電所ｊ（ｊ＝１〜Ｍ）の時点ｉの計測された有効電力をＰ（ｉ，ｊ）とする。式（１）に示すように、有効電力の１時点前後の時間差分量ＤＰ（ｉ、ｊ）を定義する。これは、有効電力は、平均値が大きく変動する非定常な時系列であるので、時間差分を取ることによって、定常な時系列に変換するためである。
【００３２】
需要予測部４は、式（２）に示す次数ｓのＡＲモデルに従って、変電所ｊの現在までに計測された有効電力の差分量を用いて、変電所ｊの次の時点ｎの有効電力の差分量の予測値ＤＰe（ｊ）を算出する。ここで、｛ａ1(j),ａ2(j),，ａs(j)｝は、変電所ｊの有効電力に関するＡＲ係数のセットである。
【００３３】
需要予測部４は、変電所ｊの次の時点の有効電力の差分量の予測値ＤＰe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された有効電力Ｐ（ｎ−１，ｊ）を用いて、式（３）にしたがって、変電所ｊの次の時点ｎの有効電力予測量Ｐe（ｊ）を算出する。
【００３４】
【数１】
2010136600.tif 000003
【００３５】
需要予測部４は、式（４）〜式（７）に従って、変電所ｊの現在までに計測された有効電力の差分量を用いて、ＡＲ係数｛ａ1(j),ａ2(j),，ａs(j)｝を算出する。
【００３６】
【数２】
2010136600.tif 000004
【００３７】
次に、無効電力予測量の自己回帰モデルについて説明する。
変電所ｊ（ｊ＝１〜Ｍ）の時点ｉの計測された無効電力をＱ（ｉ，ｊ）とする。式（８）に示すように、無効電力の１時点前後の時間差分量ＤＱ（ｉ、ｊ）を定義する。これは、無効電力は、平均値が大きく変動する非定常な時系列であるので、時間差分を取ることによって、定常な時系列に変換するためである。
【００３８】
需要予測部４は、式（９）に示す次数ｓのＡＲモデルに従って、変電所ｊの現在までに計測された無効電力の差分量を用いて、変電所ｊの次の時点ｎの無効電力の差分量の予測値ＤＱe（ｊ）を算出する。ここで、｛ｂ1(j),ｂ2(j),，ｂs(j)｝は、変電所ｊの無効電力に関するＡＲ係数のセットである。
【００３９】
需要予測部４は、変電所ｊの次の時点の無効電力の差分量の予測値ＤＱe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された無効電力Ｑ（ｎ−１，ｊ）を用いて、式（１０）にしたがって、変電所ｊの次の時点ｎの無効電力予測量Ｑe（ｊ）を算出する。
【００４０】
【数３】
2010136600.tif 000005
【００４１】
需要予測部４は、式（１１）〜式（１４）に従って、変電所ｊの現在までに計測された無効電力の差分量を用いて、ＡＲ係数｛ｂ1(j),ｂ2(j),，ｂs(j)｝を算出する。
【００４２】
【数４】
2010136600.tif 000006
【００４３】
（最適化部）
最適化部６は、ＰＳＯに基づいて、各変電所の調相設備群３２およびタップ付き変圧器３０の最適制御量を探索する。以下、最適制御量を探索する具体的な方法について説明する。
【００４４】
最適化部６は、式（１５）に示すように、エージェント（パーティクル）ｉの反復回数ｋの位置ベクトルＸikを設定する。エージェントｉは、エージェント１〜エージェントＩのＩ個とする。
【００４５】
【数５】
2010136600.tif 000007
【００４６】
式（１５）における、ＳＣik（ｊ）は、式（１６）および式（１７）で表わされる。式（１６）および式（１７）において、Ｎik＿ＳＣ２０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）におけるエージェントｉの反復回数ｋでの２０ＭＶＡの調相設備の投入台数を表わし、Ｎik＿ＳＣ３０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）におけるエージェントｉの反復回数ｋでの３０ＭＶＡの調相設備の投入台数を表わし、Ｎik＿ＳＣ４０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）におけるエージェントｉの反復回数ｋでの４０ＭＶＡの調相設備の投入台数を表わし、Ｎik＿ＳＣ６０（ｊ）は、変電所ｊ（ｊ＝ＬＳ＋１〜Ｌ）におけるエージェントｉの反復回数ｋでの６０ＭＶＡの調相設備の投入台数を表わす。
【００４７】
【数６】
2010136600.tif 000008
【００４８】
式（１５）におけるＴＡＰik（ｊ）は、変電所ｊ（ｊ＝１〜Ｍ）におけるエージェントｉの反復回数ｋでのタップ付き変圧器３０のタップ位置（中間位置を原点とする）を表わす。変電所ｊにおけるエージェントｉの反復回数ｋでのタップ付き変圧器３０のタップ比は、式（１８）で定義される。ｒｍ（ｊ）は、変電所ｊの変圧器のタップ位置が中間位置でのタップ比であり、ｒ０（ｊ）は、変電所ｊの変圧器のタップ位置が１つ変化するごとのタップ比の変化量である。ＴＡＰik（ｊ）は、タップ位置が中間位置のときに「０」となり、タップ比が増加する方向のタップ位置を正とし、タップ比が減少する方向のタップ位置を負とする。
【００４９】
【数７】
2010136600.tif 000009
【００５０】
式（１５）における要素数をＣＸ（＝３×ＬＳ＋４×（Ｌ−ＬＳ）＋Ｍ）とする。
最適化部６は、式（１９）に従って、エージェントｉの反復回数ｋの速度ベクトルＶikを算出する。ここで、Ｐｂｅｓｔik-1は、エージェントｉの反復回数（ｋ−１）での最良位置ベクトルである。Ｇｂｅｓｔk-1は、全エージェントの反復回数（ｋ−１）での最良位置ベクトルである。Ｃ１は反復回数ｋが増加するにつれて減少する変数であり、Ｃ２およびＣ３は定数である。ｒａｎｄは、０〜１の乱数値である。
【００５１】
【数８】
2010136600.tif 000010
【００５２】
最適化部６は、式（２０）に従って、エージェントｉの反復回数ｋの速度ベクトルＶikと、エージェントｉの反復回数（ｋ−１）の位置ベクトルＸik-1とに基づいて、エージェントｉの反復回数ｋの位置ベクトルＸikを算出する。
【００５３】
【数９】
2010136600.tif 000011
【００５４】
次に、最適化部６による、エージェントｉの反復回数ｋでの最良位置ベクトルＰｂｅｓｔikと、全エージェントの反復回数ｋでの最良位置ベクトルＧｂｅｓｔkの算出方法について説明する。
【００５５】
最適化部６は、式（２１）に従って、エージェントｉの反復回数ｋでのタップ位置の変化量ΔＴＡＰikを算出する。ここで、ＴＡＰ（ｎ−１，ｊ）は、変電所ｊ（ｊ＝１〜Ｍ）の時点（ｎ−１）でのタップ付き変圧器３０のタップ位置を表わす。つまり、ＴＡＰ（ｎ−１，ｊ）は、最適化部６で計算された時点（ｎ−１）での変電所ｊのタップ付き変圧器３０のタップ位置の最適値であって、これに基づいて、時点（ｎ−１）において、実際に変電所ｊのタップ付き変圧器３０のタップ位置が制御される。
【００５６】
【数１０】
2010136600.tif 000012
【００５７】
最適化部６は、式（２２）に従って、エージェントｉの反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを算出する。ここで、ここで、Ｎ＿ＳＣ２０（ｎ−１，ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）の時点（ｎ−１）での２０ＭＶＡの調相設備の投入台数を表わす。つまり、Ｎ＿ＳＣ２０（ｎ−１，ｊ）は、最適化部６で計算された時点（ｎ−１）での変電所ｊの２０ＭＶＡの調相設備の投入台数の最適値であって、これに基づいて、時点（ｎ−１）において、実際に変電所ｊの２０ＭＶＡの調相設備の投入台数が制御される。Ｎ＿ＳＣ３０（ｎ−１，ｊ）（ｊ＝１〜Ｌ）、Ｎ＿ＳＣ４０（ｎ−１，ｊ）（ｊ＝１〜Ｌ）、およびＮ＿ＳＣ６０（ｎ−１，ｊ）（ｊ＝（ＬＳ＋１）〜Ｌ）も、同様に、それぞれ変電所ｊの時点（ｎ−１）での３０ＭＶＡの調相設備の投入台数、４０ＭＶＡの調相設備の投入台数、６０ＭＶＡの調相設備の投入台数を表わす。
【００５８】
【数１１】
2010136600.tif 000013
【００５９】
最適化部６は、式（２３）に従って、潮流計算部８で算出されたエージェントｉの反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikと時点（ｎ−１）での電力系統の有効電力ロスＰｌｏｓｓ（ｎ−１）に基づいて、エージェントｉの反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossikを算出する。
【００６０】
【数１２】
2010136600.tif 000014
【００６１】
最適化部６は、式（２４）に示す評価関数に従って、エージェントｉの反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）を算出する。この評価関数は、エージェントｉの反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossik、エージェントｉの反復回数ｋでのタップ位置の変化量ΔＴＡＰik、およびエージェントｉの反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを変数として含む。この評価関数を最小化する方向に探索が行なわれる。
【００６２】
【数１３】
2010136600.tif 000015
【００６３】
最適化部６は、エージェントｉの反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）が、エージェントｉの反復回数（ｋ−１）での最良評価値ＥＶ（Ｐｂｅｓｔik-1）よりも小さい場合には、式（２５）で示されるように、エージェントｉの反復回数ｋでの最良位置ベクトルＰｂｅｓｔikの値をエージェントｉの反復回数ｋでの位置ベクトルＸikのベクトルの値と同一の値に設定し、式（２６）で示されるように、エージェントｉの反復回数ｋでの最良評価値ＥＶ（Ｐｂｅｓｔik）の値をエージェントｉの反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）と同一の値に設定する。
【００６４】
【数１４】
2010136600.tif 000016
【００６５】
最適化部６は、エージェントｉの反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）が、エージェントｉの反復回数（ｋ−１）での最良評価値ＥＶ（Ｐｂｅｓｔik-1）よりも等しいか大きい場合には、式（２７）で示されるように、エージェントｉの反復回数ｋでの最良位置ベクトルＰｂｅｓｔikの値をエージェントｉの反復回数（ｋ−１）での最良位置ベクトルＰｂｅｓｔik-1と同一の値に設定し、式（２８）で示されるようにエージェントｉの反復回数ｋでの最良評価値ＥＶ（Ｐｂｅｓｔik）の値をエージェントｉの反復回数（ｋ−１）での最良評価値ＥＶ（Ｐｂｅｓｔik-1）と同一の値に設定する。
【００６６】
【数１５】
2010136600.tif 000017
【００６７】
最適化部６は、全エージェント（ｉ＝１〜Ｉ）の最良評価値ＥＶ（Ｐｂｅｓｔik）の中で最小となるｉがｕであるときに、式（２９）に示すように、全エージェントの反復回数ｋでの最良位置ベクトルＧｂｅｓｔkの値をエージェントｕの反復回数ｋでの最良位置ベクトルＰｂｅｓｔukと同一の値に設定する。
【００６８】
【数１６】
2010136600.tif 000018
【００６９】
最適化部６は、式（３０）、（３１）、（３２）で表わされる次の時点ｎの最適制御量ＣＯＮ（ｎ）の値を、式（３４）で示されるように、最終の反復回数Ｋにおける全エージェントの最良位置ベクトルＧｂｅｓｔKと同一の値に設定する。ここで、Ｎ＿ＳＣ２０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における２０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ３０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における３０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ４０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における４０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ６０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）における６０ＭＶＡの調相設備の投入台数を表わす。また、ＴＡＰ（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｍ）におけるタップ付き変圧器３０のタップ位置（中間位置を原点とする）を表わす。
【００７０】
【数１７】
2010136600.tif 000019
【００７１】
（潮流計算部）
潮流計算部８は、各時点の最適制御量の探索中において、発電機Ｈをスイングノードとし、潮流計算を行なう。すなわち、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、次の時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、次の時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおけるエージェントｉの反復回数ｋでのタップ付き変圧器３０のタップ位置ＴＡＰik（１）〜ＴＡＰik（Ｍ）、変電所１〜変電所Ｌにおけるエージェントｉの反復回数ｋでの２０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ２０（１）〜Ｎik＿ＳＣ２０（Ｌ）、変電所１〜変電所Ｌにおけるエージェントｉの反復回数ｋでの３０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ３０（１）〜Ｎik＿ＳＣ３０（Ｌ）、変電所１〜変電所Ｌにおけるエージェントｉの反復回数ｋでの４０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ４０（１）〜Ｎik＿ＳＣ４０（Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおけるエージェントｉの反復回数ｋでの６０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ６０（ＬＳ＋１）〜Ｎik＿ＳＣ６０（Ｌ）に基づいて、潮流計算によってスイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。さらに、潮流計算部８は、式（３４）に従って、エージェントｉの反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikを算出する。潮流計算の手法は、たとえば、文献（「電力系統過渡解析論」、著者関根秦次、オーム社）などに記載された公知の手法と同等のものである。
【００７２】
【数１８】
2010136600.tif 000020
【００７３】
潮流計算部８は、時点ｎで最適制御量ＣＯＮ（ｎ）が決定された後、発電機Ｈをスイングノードとし、潮流計算を行なう。すなわち、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、次の時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、次の時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける次の時点ｎのタップ付き変圧器３０のタップ位置ＴＡＰ（ｎ，１）〜ＴＡＰ（ｎ，Ｍ）、変電所１〜変電所Ｌにおける次の時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，１）〜Ｎ＿ＳＣ２０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける次の時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，１）〜Ｎ＿ＳＣ３０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける次の時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，１）〜Ｎ＿ＳＣ４０（ｎ，Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける次の時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ＬＳ＋１）〜Ｎ＿ＳＣ６０（ｎ，Ｌ）に基づいて、潮流計算によって、スイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。潮流計算部８は、式（３５）に従って、次の時点ｎの電力系統の有効電力ロスＰｌｏｓｓ（ｎ）を算出する。
【００７４】
【数１９】
2010136600.tif 000021
【００７５】
（通信部）
通信部３は、各変電所ｊから通信回線ｊを通じて、検出された有効電力Ｐ（ｎ，ｊ）および無効電力Ｑ（ｎ，ｊ）を受信する。また、通信部３は、各変電所ｊへ通信回線ｊを通じて、変電所ｊの設備の最適制御量である、時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ｊ）（ｊ＝（ＬＳ＋１）〜Ｌの場合）、およびタップ付き変圧器３０の時点ｎのタップ位置ＴＡＰ（ｎ，ｊ）（ｊ＝１〜Ｍの場合）を受信する。
【００７６】
（動作）
次に、第１の実施形態における電圧無効電力制御の動作を説明する。
【００７７】
図４は、第１の実施形態の電圧無効電力制御システムの動作手順を表わすフローチャートである。
【００７８】
図４を参照して、変電所ｊ（ｊ＝１〜Ｍ）のＰ検出部２４は、時点ｎ（ｎ＝１〜Ｎ）の有効電力Ｐ（ｎ，ｊ）を検出する。また、変電所ｊのＱ検出部２６は、時点ｎ（ｎ＝１〜Ｎ）の無効電力Ｑ（ｎ，ｊ）を検出する。変電所ｊの通信部２２は、検出された有効電力Ｐ（ｎ，ｊ）および無効電力Ｑ（ｎ，ｊ）をＶＱＣ装置１の通信部３に送信する（ステップＳ１０１）。
【００７９】
次に、ＶＱＣ装置１において、時点ｎがＮ＋１に設定される（ステップＳ１０２）。
次に、変電所ｊ（ｊ＝１〜Ｍ）のＰ検出部２４は、時点ｎの有効電力Ｐ（ｎ，ｊ）を検出する。また、変電所ｊのＱ検出部２６は、時点ｎの無効電力Ｑ（ｎ，ｊ）を検出する。変電所ｊの通信部２２は、検出された有効電力Ｐ（ｎ，ｊ）および無効電力Ｑ（ｎ，ｊ）をＶＱＣ装置１の通信部３に送信する（ステップＳ１０３）
次に、ＶＱＣ装置１の需要予測部４は、式（４）〜式（７）に従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された有効電力の差分量を用いて、ＡＲ係数｛ａ1(j),ａ2(j),，ａs(j)｝を算出する。また、ＶＱＣ装置１の需要予測部４は、式（１１）〜式（１４）に従って、変電所ｊの現在までに計測された無効電力の差分量を用いて、ＡＲ係数｛ｂ1(j),ｂ2(j),，ｂs(j)｝を算出する（ステップＳ１０４）。
【００８０】
次に、ＶＱＣ装置１において、時点ｎが１だけインクリメントされる（ステップＳ１０５）。
【００８１】
次に、需要予測部４は、式（２）に示す次数ｓのＡＲモデルに従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された有効電力の差分量を用いて、変電所ｊの時点ｎの有効電力の差分量の予測値ＤＰe（ｊ）を算出する。需要予測部４は、変電所ｊの時点ｎの有効電力の差分量の予測値ＤＰe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された有効電力Ｐ（ｎ−１，ｊ）を用いて、式（３）にしたがって、変電所ｊの時点ｎの有効電力予測量Ｐe（ｊ）を算出する。
【００８２】
また、需要予測部４は、式（９）に示す次数ｓのＡＲモデルに従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された無効電力の差分量を用いて、変電所ｊの時点ｎの無効電力の差分量の予測値ＤＱe（ｊ）を算出する。需要予測部４は、変電所ｊの時点ｎの無効電力の差分量の予測値ＤＱe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された無効電力Ｑ（ｎ−１，ｊ）を用いて、式（１０）にしたがって、変電所ｊの時点ｎの有効電力予測量Ｑe（ｊ）を算出する（ステップＳ１０６）。
【００８３】
次に、ＶＱＣ装置１の最適化部６は、ｉ＝６〜Ｉについて、式（１５）〜（１８）で示すようなエージェントｉの反復回数１での位置ベクトルＸi1を擬似乱数で設定する。また、最適化部６は、エージェント１〜５の反復回数１での位置ベクトルＸp1（ｐ＝１〜５）の値を時点（ｎ−１）の全エージェントの反復回数Ｋでの最良位置ベクトルＧｂｅｓｔKと同一の値に設定する。また、最適化部６は、エージェントｉ（ｉ＝１〜Ｉ）の反復回数１での速度ベクトルＶi1を擬似乱数で設定する。さらに、最適化部６は、エージェントｉ（ｉ＝１〜Ｉ）の反復回数１での最良位置ベクトルＰｂｅｓｔi1の値を位置ベクトルＸi1の値と同一の値に設定する。さらに、最適化部６は、全エージェントの反復回数１での最良位置ベクトルＧｂｅｓｔ1の値をエージェントｉ（ｉ＝１〜Ｉ）の反復回数１での最良位置ベクトルＰｂｅｓｔi1のうちの最良のものと同一の値に設定する。最適化部６は、反復回数ｋ＝２に設定する（ステップＳ１０７）。
【００８４】
次に、最適化部６は、式（１９）に従って、ｉ＝１〜Ｉについて、エージェントｉの反復回数ｋでの速度ベクトルＶikを算出する（ステップＳ１０８）。
【００８５】
次に、最適化部６は、式（２０）に従って、ｉ＝１〜Ｉについて、エージェントｉの反復回数ｋでの位置ベクトルＸikを算出する（ステップＳ１０９）。
【００８６】
次に、最適化部６は、式（２１）に従って、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでのタップ位置の変化量ΔＴＡＰikを算出する。また、最適化部６は、式（２２）に従って、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを算出する（ステップＳ１１０）。
【００８７】
次に、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおけるエージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでのタップ付き変圧器３０のタップ位置ＴＡＰik（１）〜ＴＡＰik（Ｍ）、変電所１〜変電所Ｌにおけるエージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの２０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ２０（１）〜Ｎik＿ＳＣ２０（Ｌ）、変電所１〜変電所Ｌにおけるエージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの３０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ３０（１）〜Ｎik＿ＳＣ３０（Ｌ）、変電所１〜変電所Ｌにおけるエージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの４０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ４０（１）〜Ｎik＿ＳＣ４０（Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおけるエージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの６０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ６０（ＬＳ＋１）〜Ｎik＿ＳＣ６０（Ｌ）に基づいて、潮流計算によってスイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。さらに、潮流計算部８は、式（３４）に従って、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikを算出する。最適化部６は、式（２３）に従って、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossikを算出する（ステップＳ１１１）。
【００８８】
次に、最適化部６は、式（２４）に従って、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）を算出する。最適化部６は、エージェントｉの反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）が、エージェントｉの反復回数（ｋ−１）での最良評価値ＥＶ（Ｐｂｅｓｔik-1）よりも小さい場合には、式（２５）および（２６）で示されるように、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの最良位置ベクトルＰｂｅｓｔikの値をエージェントｉの反復回数ｋでの位置ベクトルＸikと同一の値に設定し、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの最良評価値ＥＶ（Ｐｂｅｓｔik）の値をエージェントｉの反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）と同一の値に設定する。最適化部６は、エージェントｉの反復回数ｋでの位置ベクトルＸikの評価値ＥＶ（Ｘik）が、エージェントｉの反復回数（ｋ−１）での最良評価値ＥＶ（Ｐｂｅｓｔik-1）よりも等しいか大きい場合には、式（２７）および（２８）に示されるように、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの最良位置ベクトルＰｂｅｓｔikの値をエージェントｉの反復回数（ｋ−１）での最良位置ベクトルＰｂｅｓｔik-1と同一の値に設定し、エージェントｉ（ｉ＝１〜Ｉ）の反復回数ｋでの最良評価値ＥＶ（Ｐｂｅｓｔik）の値をエージェントｉの反復回数（ｋ−１）での最良評価値ＥＶ（Ｐｂｅｓｔik-1）と同一の値に設定する（ステップＳ１１２）。
【００８９】
次に、最適化部６は、全エージェント（ｉ＝１〜Ｉ）の最良評価値ＥＶ（Ｐｂｅｓｔik）の中で最小となるｉがｕであるときに、式（２９）に示すように、全エージェントの反復回数ｋでの最良位置ベクトルＧｂｅｓｔkの値をエージェントｕの反復回数ｋでの最良位置ベクトルＰｂｅｓｔukと同一の値に設定する（ステップＳ１１３）。
【００９０】
次に、最適化部６は、ｋ＝Ｋであれば（ステップＳ１１４でＹＥＳ）、ステップＳ１１６に進み、ｋ＝Ｋでなければ（ステップＳ１１４でＮＯ）、ステップＳ１１５に移行して、反復回数ｋを１だけ増加した上で、ステップＳ１０８に戻る。
【００９１】
ステップＳ１１６において、最適化部６は、式（３０）、（３１）、（３２）で表わされる時点ｎの最適制御量ＣＯＮ（ｎ）の値を、式（３３）で示されるように、最終の反復回数Ｋにおける全エージェントの最良位置ベクトルＧｂｅｓｔKと同一の値に設定する。
【００９２】
さらに、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける時点ｎのタップ付き変圧器３０のタップ位置ＴＡＰ（ｎ，１）〜ＴＡＰ（ｎ，Ｍ）、変電所１〜変電所Ｌにおける時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，１）〜Ｎ＿ＳＣ２０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，１）〜Ｎ＿ＳＣ３０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，１）〜Ｎ＿ＳＣ４０（ｎ，Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ＬＳ＋１）〜Ｎ＿ＳＣ６０（ｎ，Ｌ）に基づいて、潮流計算によって、スイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。潮流計算部８は、式（３５）に従って、時点ｎの電力系統の有効電力ロスＰｌｏｓｓ（ｎ）を算出する。
【００９３】
さらに、通信部３は、変電所ｊに対して、最適制御量ＣＯＮ（ｎ）のうち、変電所ｊに関連する情報である、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，ｊ）、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，ｊ）、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，ｊ）、変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）の時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ｊ）、および変電所ｊ（ｊ＝１〜Ｍ）のタップ付き変圧器３０の時点ｎのタップ位置ＴＡＰ（ｎ，ｊ）を送信する（ステップＳ１１６）。
【００９４】
変電所ｊは、情報Ｎ＿ＳＣ２０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ３０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ４０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ６０（ｎ，ｊ）（ｊ＝（ＬＳ＋１）〜Ｌの場合）、およびＴＡＰ（ｎ，ｊ）（ｊ＝１〜Ｍの場合）を受信する。変電所の制御部２８は、これらの受信した情報に基づいて、タップ付き変圧器３０のタップ位置と、調相設備群３２の投入台数を制御する（ステップＳ１１７）。
【００９５】
（小括）
以上のように、本発明の実施形態の電圧無効電力制御システムによれば、将来の有効電力および無効電力を予測して、予測に基づいてＰＳＯによって最適制御量を探索するので、電圧制御機器の制御遅れを考慮した電圧無効電力制御を行なうことができる。
【００９６】
［第２の実施形態］
第２の実施形態は、タブーサーチ（Tabu Search:ＴＳ）に基づいて、最適制御量を探索する電圧無効電力制御システムに関する。
【００９７】
第２の実施形態の電圧無効電力制御システムは、最適化部６と潮流計算部８を除いて、第１の実施形態と同様である。以下では、最適化部６と潮流計算部８について説明する。
【００９８】
（最適化部）
最適化部６は、ＴＳに基づいて、各変電所の調相設備群３２およびタップ付き変圧器３０の最適制御量を探索する。以下、最適制御量を探索する具体的な方法について説明する。
【００９９】
最適化部６は、式（３６）に示すように、反復回数ｋの状態変数ベクトルＸkを設定する。
【０１００】
【数２０】
2010136600.tif 000022
【０１０１】
式（３６）におけるＳＣk（ｊ）は、式（３７）および式（３８）で表わされる。式（３７）および式（３８）において、Ｎk＿ＳＣ２０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）における反復回数ｋでの２０ＭＶＡの調相設備の投入台数を表わし、Ｎk＿ＳＣ３０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）における反復回数ｋでの３０ＭＶＡの調相設備の投入台数を表わし、Ｎk＿ＳＣ４０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）における反復回数ｋでの４０ＭＶＡの調相設備の投入台数を表わし、Ｎk＿ＳＣ６０（ｊ）は、変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）における反復回数ｋでの６０ＭＶＡの調相設備の投入台数を表わす。
【０１０２】
【数２１】
2010136600.tif 000023
【０１０３】
式（３６）におけるＴＡＰk（ｊ）は、変電所ｊ（ｊ＝１〜Ｍ）における反復回数ｋでのタップ付き変圧器３０のタップ位置（中間位置を原点とする）を表わす。変電所ｊにおける反復回数ｋでのタップ付き変圧器３０のタップ比は、式（３９）で定義される。ｒｍ（ｊ）は、変電所ｊの変圧器のタップ位置が中間位置でのタップ比であり、ｒ０（ｊ）は、変電所ｊの変圧器のタップ位置が１つ変化するごとのタップ比の変化量である。ＴＡＰk（ｊ）は、タップ位置が中間位置のときに「０」となり、タップ比が増加する方向のタップ位置を正とし、タップ比が減少する方向のタップ位置を負とする。
【０１０４】
【数２２】
2010136600.tif 000024
【０１０５】
式（３６）における要素数をＣＸ（＝３×ＬＳ＋４×（Ｌ−ＬＳ）＋Ｍ）とする。
次に、最適化部６による、次の反復回数（ｋ＋１）の状態変数ベクトルＸk+1および最良状態変数ベクトルＳk+1の算出方法について説明する。
【０１０６】
最適化部６は、式（４０）〜（４６）に従って、反復回数ｋでの状態変数ベクトルＸkの近傍ｉ（ｉ＝１〜ＣＸ）の状態変数ベクトルＲＸikを算出する。
【０１０７】
【数２３】
2010136600.tif 000025
【０１０８】
【数２４】
2010136600.tif 000026
【０１０９】
【数２５】
2010136600.tif 000027
【０１１０】
【数２６】
2010136600.tif 000028
【０１１１】
最適化部６は、、式（４７）に従って、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでのタップ位置の変化量ΔＴＡＰikを算出する。ここで、ＴＡＰ（ｎ−１，ｊ）は、変電所ｊ（ｊ＝１〜Ｍ）の時点（ｎ−１）でのタップ付き変圧器３０のタップ位置を表わす。つまり、ＴＡＰ（ｎ−１，ｊ）は、最適化部６で計算された時点（ｎ−１）での変電所ｊのタップ付き変圧器３０のタップ位置の最適値であって、これに基づいて、時点（ｎ−１）において、実際に変電所ｊのタップ付き変圧器３０のタップ位置が制御される。
【０１１２】
【数２７】
2010136600.tif 000029
【０１１３】
最適化部６は、式（４８）に従って、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを算出する。ここで、ここで、Ｎ＿ＳＣ２０（ｎ−１，ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）の時点（ｎ−１）での２０ＭＶＡの調相設備の投入台数を表わす。つまり、Ｎ＿ＳＣ２０（ｎ−１，ｊ）は、最適化部６で計算された時点（ｎ−１）での変電所ｊの２０ＭＶＡの調相設備の投入台数の最適値であって、これに基づいて、時点（ｎ−１）において、実際に変電所ｊの２０ＭＶＡの調相設備の投入台数が制御される。Ｎ＿ＳＣ３０（ｎ−１，ｊ）（ｊ＝１〜Ｌ）、Ｎ＿ＳＣ４０（ｎ−１，ｊ）（ｊ＝１〜Ｌ）、およびＮ＿ＳＣ６０（ｎ−１，ｊ）（ｊ＝（ＬＳ＋１）〜Ｌ）も、同様に、それぞれ変電所ｊの時点（ｎ−１）での３０ＭＶＡの調相設備の投入台数、４０ＭＶＡの調相設備の投入台数、６０ＭＶＡの調相設備の投入台数を表わす。
【０１１４】
【数２８】
2010136600.tif 000030
【０１１５】
最適化部６は、式（４９）に従って、潮流計算部８で算出された近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikと時点（ｎ−１）での電力系統の有効電力ロスＰｌｏｓｓ（ｎ−１）に基づいて、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossikを算出する。
【０１１６】
【数２９】
2010136600.tif 000031
【０１１７】
最適化部６は、式（５０）に示す評価関数に従って、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの状態変数ベクトルＲＸikの評価値ＥＶ（ＲＸik）を算出する。この評価関数は、近傍ｉの反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossik、近傍ｉの反復回数ｋでのタップ位置の変化量ΔＴＡＰik、および近傍ｉの反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを変数として含む。
【０１１８】
【数３０】
2010136600.tif 000032
【０１１９】
最適化部６は、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの状態変数ベクトルＲＸikの評価値ＥＶ（ＲＸik）のうちの最小値ＥＶ（ＲＸjk）を見つけて、これと反復回数ｋでの最良状態変数ベクトルＳkでの評価値ＥＶ（Ｓk）とを比較することによって、次の反復回数（ｋ＋１）での状態変数ベクトルＸk+1と最良状態変数ベクトルＳk+1の設定を行なうとともに、タブーリストの更新を行なう。
【０１２０】
（１）  状態変数ベクトルＸk+1の設定
最適化部６は、近傍ｊの状態変数ベクトルＲＸikがタブーリストに登録されていない場合には、式（５１）に従って状態変数ベクトルＸk+1を設定し、近傍ｊの状態変数ベクトルＲＸikがタブーリストに登録されている場合には、タブーリストに登録されていない近傍ｉ（ｉ＝１〜ＣＸ、ｉ≠ｊ）のうち評価値が最小である近傍ｕを特定し、式（５２）に従って状態変数ベクトルＸk+1を設定する。
【０１２１】
【数３１】
2010136600.tif 000033
【０１２２】
（２）  最良状態変数ベクトルＳk+1の設定
最適化部６は、評価値ＥＶ（ＲＸjk）が評価値ＥＶ（Ｓk）よりも小さい場合には、式（５３）に従って最良状態変数ベクトルＳk+1を設定し、評価値ＥＶ（ＲＸjk）が評価値ＥＶ（Ｓk）よりも等しいか大きい場合には、式（５４）に従って、最良状態変数ベクトルＳk+1を設定する。
【０１２３】
【数３２】
2010136600.tif 000034
【０１２４】
（３）  タブーリストの更新
最適化部６は、タブーリストに反復回数ｋでの状態変数ベクトルＸkを登録する。この際に、タブーリストがフルの場合には、最適化部６は、最も古く登録された状態変数ベクトルをタブーリストから削除した上で、状態変数ベクトルＸkを登録する。
【０１２５】
最適化部６は、式（５５）、（５６）、（５７）で表わされる次の時点ｎの最適制御量ＣＯＮ（ｎ）の値を、式（５８）で示されるように、最終の反復回数Ｋにおける最良状態変数ベクトルＳKと同一の値に設定する。ここで、Ｎ＿ＳＣ２０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における２０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ３０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における３０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ４０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における４０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ６０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）における６０ＭＶＡの調相設備の投入台数を表わす。また、ＴＡＰ（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｍ）におけるタップ付き変圧器３０のタップ位置（中間位置を原点とする）を表わす。
【０１２６】
【数３３】
2010136600.tif 000035
【０１２７】
（潮流計算部）
潮流計算部８は、各時点の最適制御量の探索中において、発電機Ｈをスイングノードとし、潮流計算を行なう。すなわち、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、次の時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、次の時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける近傍ｉの反復回数ｋでのタップ付き変圧器３０のタップ位置ＴＡＰik（１）〜ＴＡＰik（Ｍ）、変電所１〜変電所Ｌにおける近傍ｉの反復回数ｋでの２０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ２０（１）〜Ｎik＿ＳＣ２０（Ｌ）、変電所１〜変電所Ｌにおける近傍ｉの反復回数ｋでの３０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ３０（１）〜Ｎik＿ＳＣ３０（Ｌ）、変電所１〜変電所Ｌにおける近傍ｉの反復回数ｋでの４０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ４０（１）〜Ｎik＿ＳＣ４０（Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける近傍ｉの反復回数ｋでの６０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ６０（ＬＳ＋１）〜Ｎik＿ＳＣ６０（Ｌ）に基づいて、潮流計算によってスイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。さらに、潮流計算部８は、式（５９）に従って、近傍ｉの反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikを算出する。
【０１２８】
【数３４】
2010136600.tif 000036
【０１２９】
潮流計算部８は、時点ｎで最適制御量ＣＯＮ（ｎ）が決定された後、発電機Ｈをスイングノードとし、潮流計算を行なう。すなわち、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、次の時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、次の時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける次の時点ｎのタップ付き変圧器３０のタップ位置ＴＡＰ（ｎ，１）〜ＴＡＰ（ｎ，Ｍ）、変電所１〜変電所Ｌにおける次の時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，１）〜Ｎ＿ＳＣ２０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける次の時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，１）〜Ｎ＿ＳＣ３０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける次の時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，１）〜Ｎ＿ＳＣ４０（ｎ，Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける次の時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ＬＳ＋１）〜Ｎ＿ＳＣ６０（ｎ，Ｌ）に基づいて、潮流計算によって、スイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。潮流計算部８は、式（６０）に従って、次の時点ｎの電力系統の有効電力ロスＰｌｏｓｓ（ｎ）を算出する。
【０１３０】
【数３５】
2010136600.tif 000037
【０１３１】
（動作）
次に、第２の実施形態における電圧無効電力制御の動作を説明する。
【０１３２】
図５は、第２の実施形態の電圧無効電力制御システムの動作手順を表わすフローチャートである。
【０１３３】
図５を参照して、変電所ｊ（ｊ＝１〜Ｍ）のＰ検出部２４は、時点ｎ（ｎ＝１〜Ｎ）の有効電力Ｐ（ｎ，ｊ）を検出する。また、変電所ｊのＱ検出部２６は、時点ｎ（ｎ＝１〜Ｎ）の無効電力Ｑ（ｎ，ｊ）を検出する。変電所ｊの通信部２２は、検出された有効電力Ｐ（ｎ，ｊ）および無効電力Ｑ（ｎ，ｊ）をＶＱＣ装置１の通信部３に送信する（ステップＳ２０１）。
【０１３４】
次に、ＶＱＣ装置１において、時点ｎがＮ＋１に設定される（ステップＳ２０２）。
次に、変電所ｊ（ｊ＝１〜Ｍ）のＰ検出部２４は、時点ｎの有効電力Ｐ（ｎ，ｊ）を検出する。また、変電所ｊのＱ検出部２６は、時点ｎの無効電力Ｑ（ｎ，ｊ）を検出する。変電所ｊの通信部２２は、検出された有効電力Ｐ（ｎ，ｊ）および無効電力Ｑ（ｎ，ｊ）をＶＱＣ装置１の通信部３に送信する（ステップＳ２０３）
次に、ＶＱＣ装置１の需要予測部４は、式（４）〜式（７）に従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された有効電力の差分量を用いて、ＡＲ係数｛ａ1(j),ａ2(j),，ａs(j)｝を算出する。また、ＶＱＣ装置１の需要予測部４は、式（１１）〜式（１４）に従って、変電所ｊの現在までに計測された無効電力の差分量を用いて、ＡＲ係数｛ｂ1(j),ｂ2(j),，ｂs(j)｝を算出する（ステップＳ２０４）。
【０１３５】
次に、ＶＱＣ装置１において、時点ｎが１だけインクリメントされる（ステップＳ２０５）。
【０１３６】
次に、需要予測部４は、式（２）に示す次数ｓのＡＲモデルに従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された有効電力の差分量を用いて、変電所ｊの時点ｎの有効電力の差分量の予測値ＤＰe（ｊ）を算出する。需要予測部４は、変電所ｊの時点ｎの有効電力の差分量の予測値ＤＰe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された有効電力Ｐ（ｎ−１，ｊ）を用いて、式（３）にしたがって、変電所ｊの時点ｎの有効電力予測量Ｐe（ｊ）を算出する。
【０１３７】
また、需要予測部４は、式（９）に示す次数ｓのＡＲモデルに従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された無効電力の差分量を用いて、変電所ｊの時点ｎの無効電力の差分量の予測値ＤＱe（ｊ）を算出する。需要予測部４は、変電所ｊの時点ｎの無効電力の差分量の予測値ＤＱe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された無効電力Ｑ（ｎ−１，ｊ）を用いて、式（１０）にしたがって、変電所ｊの時点ｎの有効電力予測量Ｑe（ｊ）を算出する（ステップＳ２０６）。
【０１３８】
次に、ＶＱＣ装置１の最適化部６は、式（３６）〜（３９）で示すような反復回数１での状態変数ベクトルＸ1の値を時点（ｎ−１）での反復回数Ｋでの最良状態変数ベクトルＳKと同一の値に設定する。さらに、最適化部６は、反復回数１での最良状態変数ベクトルＳ1の値を状態変数ベクトルＸ1と同一の値に設定する。最適化部６は、反復回数ｋ＝１に設定する（ステップＳ２０７）。
【０１３９】
次に、最適化部６は、式（４０）〜（４６）に従って、反復回数ｋでの状態変数ベクトルＸkの近傍ｉ（ｉ＝１〜ＣＸ）の状態変数ベクトルＲＸikを算出する（ステップＳ２０８）。
【０１４０】
次に、最適化部６は、式（４７）に従って、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでのタップ位置の変化量ΔＴＡＰikを算出する。また、最適化部６は、式（４８）に従って、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを算出する（ステップＳ２０９）。
【０１４１】
次に、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでのタップ付き変圧器３０のタップ位置ＴＡＰik（１）〜ＴＡＰik（Ｍ）、変電所１〜変電所Ｌにおける近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの２０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ２０（１）〜Ｎik＿ＳＣ２０（Ｌ）、変電所１〜変電所Ｌにおける近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの３０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ３０（１）〜Ｎik＿ＳＣ３０（Ｌ）、変電所１〜変電所Ｌにおける近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの４０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ４０（１）〜Ｎik＿ＳＣ４０（Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの６０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ６０（ＬＳ＋１）〜Ｎik＿ＳＣ６０（Ｌ）に基づいて、潮流計算によってスイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。さらに、潮流計算部８は、式（５９）に従って、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikを算出する。最適化部６は、式（４９）に従って、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossikを算出する（ステップＳ２１０）。
【０１４２】
次に、最適化部６は、近傍ｉ（ｉ＝１〜ＣＸ）の反復回数ｋでの状態変数ベクトルＲＸikの評価値ＥＶ（ＲＸik）に従って、次の反復回数（ｋ＋１）での状態変数ベクトルＸk+1と最良状態変数ベクトルＳk+1の設定を行なうとともに、タブーリストの更新を行なう（ステップＳ２１１）。
【０１４３】
次に、最適化部６は、次の反復回数（ｋ＋１）がＫであれば（ステップＳ２１２でＹＥＳ）、ステップＳ２１４に進み、（ｋ＋１）がＫでなければ（ステップＳ２１２でＮＯ）、ステップＳ２１３に移行して、反復回数ｋを１だけ増加した上で、ステップＳ２０８に戻る。
【０１４４】
ステップＳ２１４において、最適化部６は、式（５５）、（５６）、（５７）で表わされる次の時点ｎの最適制御量ＣＯＮ（ｎ）の値を、式（５８）で示されるように、最終の反復回数Ｋにおける最良状態変数ベクトルＳKと同一の値に設定する。
【０１４５】
さらに、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける時点ｎのタップ付き変圧器３０のタップ位置ＴＡＰ（ｎ，１）〜ＴＡＰ（ｎ，Ｍ）、変電所１〜変電所Ｌにおける時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，１）〜Ｎ＿ＳＣ２０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，１）〜Ｎ＿ＳＣ３０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，１）〜Ｎ＿ＳＣ４０（ｎ，Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ＬＳ＋１）〜Ｎ＿ＳＣ６０（ｎ，Ｌ）に基づいて、潮流計算によって、スイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。潮流計算部８は、式（６０）に従って、時点ｎの電力系統の有効電力ロスＰｌｏｓｓ（ｎ）を算出する。
【０１４６】
さらに、通信部３は、変電所ｊに対して、最適制御量ＣＯＮ（ｎ）のうち、変電所ｊに関連する情報である、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，ｊ）、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，ｊ）、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，ｊ）、変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）の時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ｊ）、および変電所ｊ（ｊ＝１〜Ｍ）のタップ付き変圧器３０の時点ｎのタップ位置ＴＡＰ（ｎ，ｊ）を送信する（ステップＳ２１４）。
【０１４７】
変電所ｊは、情報Ｎ＿ＳＣ２０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ３０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ４０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ６０（ｎ，ｊ）（ｊ＝（ＬＳ＋１）〜Ｌの場合）、およびＴＡＰ（ｎ，ｊ）（ｊ＝１〜Ｍの場合）を受信する。変電所の制御部２８は、これらの受信した情報に基づいて、タップ付き変圧器３０のタップ位置と、調相設備群３２の投入台数を制御する（ステップＳ２１５）。
【０１４８】
（ステップＳ２１１の詳細）
図６は、図５のステップＳ２１１の詳細な手順を表わすフローチャートである。
【０１４９】
図６を参照して、最適化部６は、式（４０）〜（４６）に従って、反復回数ｋでの状態変数ベクトルＸkの近傍ｉ（ｉ＝１〜ＣＸ）の状態変数ベクトルＲＸikの評価値ＥＶ（ＲＸik）を算出する（ステップＳ３０１）。
【０１５０】
次に、最適化部６は、評価値ＥＶ（ＲＸik）（ｉ＝１〜ＣＸ）の中の最小値ＥＶ（ＲＸjk）を特定する（ステップＳ３０２）。
【０１５１】
最適化部６は、評価値ＥＶ（ＲＸjk）が、反復回数ｋでの最良状態変数ベクトルＳkの評価値ＥＶ（Ｓk）よりも小さい場合（ステップＳ３０３でＹＥＳ）には、式（５１）で示されるように、反復回数（ｋ＋１）での状態変数ベクトルＸk+1の値を状態変数ベクトルＲＸjkと同一の値に設定し（ステップＳ３０４）、式（５３）で示されるように、反復回数（ｋ＋１）での最良状態変数ベクトルＳk+1の値を状態変数ベクトルＲＸjkと同一の値に設定し（ステップＳ３０５）、状態変数ベクトルＸkをタブーリストに登録する（ステップ３０６）。
【０１５２】
また、最適化部６は、評価値ＥＶ（ＲＸjk）が、反復回数ｋでの最良状態変数ベクトルＳkの評価値ＥＶ（Ｓk）よりも等しいか大きい場合で（ステップＳ３０３でＮＯ）、かつ状態変数ベクトルＲＸjkがタブーリストに登録されていない場合（ステップＳ３０７でＹＥＳ）には、式（５１）で示されるように、反復回数（ｋ＋１）での状態変数ベクトルＸk+1の値を状態変数ベクトルＲＸjkと同一の値に設定し（ステップＳ３０８）、式（５４）で示されるように、反復回数（ｋ＋１）での最良状態変数ベクトルＳk+1の値を反復回数ｋでの最良状態変数ベクトルＳkと同一の値に設定し（ステップＳ３０９）、状態変数ベクトルＸkをタブーリストに登録する（ステップ３０６）。
【０１５３】
また、最適化部６は、評価値ＥＶ（ＲＸjk）が、反復回数ｋでの最良状態変数ベクトルＳkの評価値ＥＶ（Ｓk）よりも等しいか大きい場合で（ステップＳ３０３でＮＯ）、かつ状態変数ベクトルＲＸjkがタブーリストに登録されている場合（ステップＳ３０７でＮＯ）には、タブーリストに登録されていない近傍ｉ（ｉ＝１〜ＣＸ、ｉ≠ｊ）のうち評価値が最小である近傍ｕを特定し（ステップＳ３１０）、式（５２）で示されるように、反復回数（ｋ＋１）での状態変数ベクトルＸk+1の値を状態変数ベクトルＲＸukと同一の値に設定し（ステップＳ３１１）、式（５４）で示されるように、反復回数（ｋ＋１）での最良状態変数ベクトルＳk+1の値を反復回数ｋでの最良状態変数ベクトルＳkと同一の値に設定し（ステップＳ３０９）、状態変数ベクトルＸkをタブーリストに登録する（ステップ３０７）。
【０１５４】
（小括）
以上のように、本発明の実施形態の電圧無効電力制御システムによれば、将来の有効電力および無効電力を予測して、予測に基づいてＴＳによって最適制御量を探索するので、電圧制御機器の制御遅れを考慮した電圧無効電力制御を行なうことができる。
【０１５５】
［第３の実施形態］
第３の実施形態は、遺伝的アルゴリズム（Genetic Algorithm：ＧＡ）に基づいて、最適制御量を探索する電圧無効電力制御システムに関する。
【０１５６】
第３の実施形態の電圧無効電力制御システムは、最適化部６と潮流計算部８を除いて、第１の実施形態と同様である。以下では、最適化部６と潮流計算部８について説明する。
【０１５７】
（最適化部）
最適化部６は、ＧＡに基づいて、各変電所の調相設備群３２およびタップ付き変圧器３０の最適制御量を探索する。以下、最適制御量を探索する具体的な方法について説明する。
【０１５８】
最適化部６は、式（６１）に示すように、反復回数ｋの個体ｉの染色体Ｘikを設定する。個体ｉは、個体１〜個体ＩのＩ個とする。式（６１）〜式（６４）における各要素は、１つの遺伝子を表わす。
【０１５９】
【数３６】
2010136600.tif 000038
【０１６０】
式（６１）におけるＳＣik（ｊ）（ｊ＝１〜Ｌ）は、式（６２）および式（６３）で表わされる。式（６２）および式（６３）において、Ｎik＿ＳＣ２０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）における個体ｉの反復回数ｋでの２０ＭＶＡの調相設備の投入台数を表わし、Ｎik＿ＳＣ３０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）における個体ｉの反復回数ｋでの３０ＭＶＡの調相設備の投入台数を表わし、Ｎik＿ＳＣ４０（ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）における個体ｉの反復回数ｋでの４０ＭＶＡの調相設備の投入台数を表わし、Ｎik＿ＳＣ６０（ｊ）は、変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）における個体ｉの反復回数ｋでの６０ＭＶＡの調相設備の投入台数を表わす。
【０１６１】
式（６１）におけるＴＡＰik（ｊ）（ｊ＝１〜Ｍ）は、式（６４）で表わされる。式（６４）におけるＴＡＰik（ｊ）は、変電所ｊ（ｊ＝１〜Ｍ）における個体ｉの反復回数ｋでのタップ付き変圧器３０のタップ位置（中間位置を原点とする）を表わす。変電所ｊにおける個体ｉの反復回数ｋでのタップ付き変圧器３０のタップ比は、式（６４）で定義される。ｒｍ（ｊ）は、変電所ｊの変圧器のタップ位置が中間位置でのタップ比であり、ｒ０（ｊ）は、変電所ｊの変圧器のタップ位置が１つ変化するごとのタップ比の変化量である。ＴＡＰik（ｊ）は、タップ位置が中間位置のときに「０」となり、タップ比が増加する方向のタップ位置を正とし、タップ比が減少する方向のタップ位置を負とする。
【０１６２】
式（６１）における要素数をＣＸ（＝３×ＬＳ＋４×（Ｌ−ＬＳ）＋Ｍ）とする。
図７は、反復回数ｋの個体ｉの染色体Ｘikの例を表わす図である。
【０１６３】
図７を参照して、遺伝子Ｎik＿ＳＣ２０（１）には「５」が設定され、遺伝子Ｎik＿ＳＣ３０（１）には「３」が設定され、遺伝子Ｎik＿ＳＣ２０（２）には「４」が設定され、遺伝子ＴＡＰik（１）には「２」が設定され、遺伝子ＴＡＰik（２）には「７」が設定され、遺伝子ＴＡＰik（Ｍ）には「６」が設定されている。
【０１６４】
最適化部６は、次のようにして個体ｉの評価値を算出する。
最適化部６は、式（６５）に従って、個体ｉの反復回数ｋでのタップ位置の変化量ΔＴＡＰikを算出する。ここで、ＴＡＰ（ｎ−１，ｊ）は、変電所ｊ（ｊ＝１〜Ｍ）の時点（ｎ−１）でのタップ付き変圧器３０のタップ位置を表わす。つまり、ＴＡＰ（ｎ−１，ｊ）は、最適化部６で計算された時点（ｎ−１）での変電所ｊのタップ付き変圧器３０のタップ位置の最適値であって、これに基づいて、時点（ｎ−１）において、実際に変電所ｊのタップ付き変圧器３０のタップ位置が制御される。
【０１６５】
【数３７】
2010136600.tif 000039
【０１６６】
最適化部６は、式（６６）に従って、個体ｉの反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを算出する。ここで、ここで、Ｎ＿ＳＣ２０（ｎ−１，ｊ）は、変電所ｊ（ｊ＝１〜Ｌ）の時点（ｎ−１）での２０ＭＶＡの調相設備の投入台数を表わす。つまり、Ｎ＿ＳＣ２０（ｎ−１，ｊ）は、最適化部６で計算された時点（ｎ−１）での変電所ｊの２０ＭＶＡの調相設備の投入台数の最適値であって、これに基づいて、時点（ｎ−１）において、実際に変電所ｊの２０ＭＶＡの調相設備の投入台数が制御される。Ｎ＿ＳＣ３０（ｎ−１，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ４０（ｎ−１，ｊ）（ｊ＝１〜Ｌの場合）、およびＮ＿ＳＣ６０（ｎ−１，ｊ）（ｊ＝（ＬＳ＋１）〜Ｌの場合）も、同様に、それぞれ変電所ｊの時点（ｎ−１）での３０ＭＶＡの調相設備の投入台数、４０ＭＶＡの調相設備の投入台数、６０ＭＶＡの調相設備の投入台数を表わす。
【０１６７】
【数３８】
2010136600.tif 000040
【０１６８】
最適化部６は、式（６７）に従って、潮流計算部８で算出された個体ｉの反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikと時点（ｎ−１）での電力系統の有効電力ロスＰｌｏｓｓ（ｎ−１）に基づいて、個体ｉの反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossikを算出する。
【０１６９】
【数３９】
2010136600.tif 000041
【０１７０】
最適化部６は、式（６８）に示す評価関数に従って、個体ｉの染色体Ｘikの反復回数ｋでの評価値ＥＶ（Ｘik）を算出する。この評価関数は、個体ｉの染色体Ｘikの反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossik、個体ｉの染色体Ｘikの反復回数ｋでのタップ位置の変化量ΔＴＡＰik、および個体ｉの染色体Ｘikの反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを変数として含む。
【０１７１】
【数４０】
2010136600.tif 000042
【０１７２】
最適化部６は、次の順序で遺伝子操作を行なうことによって、反復回数（ｋ＋１）の個体ｉの染色体Ｘik+1を生成する。
【０１７３】
（１）  淘汰
最適化部６は、評価値が最大から第１番目〜第５番目の個体を評価値が最小である個体（エリート個体）で置換する。
【０１７４】
（２）  交差
最適化部６は、ＮＣＲ個（交差率×個体数Ｉ）の個体をランダムに選択して、選択した個体から（ＮＣＲ／２）個のペアを作成して、ランダムな位置でペアの２つの個体の染色体を交差させる。交差率は、たとえば０．６とする。
【０１７５】
図８は、交差の例を説明するための図である。図８（ａ）は、反復回数３の個体２および個体５の交差前の染色体を表わし、図８（ｂ）は、反復回数３の個体２および個体５の交差後の染色体を表わす。
【０１７６】
図８（ａ）および図８（ｂ）に示すように、遺伝子ＴＡＰ2(5)3（１）と遺伝子ＴＡＰ2(5)3（２）の間で交差が起こり、交差前の遺伝子ＴＡＰ23（２）〜ＴＡＰ23（Ｍ）の値が交差後の遺伝子ＴＡＰ53（２）〜ＴＡＰ53（Ｍ）の値となり、交差前の遺伝子ＴＡＰ53（２）〜ＴＡＰ53（Ｍ）の値が交差後の遺伝子ＴＡＰ23（２）〜ＴＡＰ23（Ｍ）の値となっている。
【０１７７】
（３）  突然変異
最適化部６は、ＮＶＲ個（突然変異率×個体数Ｉ）の個体をランダムに選択して、選択した各個体の染色体の１つの遺伝子をランダムに選択して、ランダムな値に変更する。突然変異率は、たとえば、０．３とする。
【０１７８】
図９は、突然変異の例を説明するための図である。図９（ａ）は、反復回数３の個体２の突然変位前の染色体を表わし、図９（ｂ）は、反復回数３の個体２の突然変異後の染色体を表わす。
【０１７９】
図９（ａ）および図９（ｂ）に示すように、遺伝子Ｎ23＿ＳＣ３０（２）で突然変異が起こり、交差前の遺伝子Ｎ23＿ＳＣ３０（２）の値が変更されている。
【０１８０】
最適化部６は、（３）の突然変位処理が終了した後、式（６９）に従って、反復回数（ｋ＋１）の個体ｉの染色体Ｘik+1を設定する。すなわち、（３）の突然変位の処理が終了した後の個体が次の世代（反復回数）の個体となる。
【０１８１】
また、最適化部６は、（３）の突然変位処理が終了した後で、評価値が最小である個体が第ｊ個体の場合には、式（７０）に従って、反復回数（ｋ＋１）のエリート個体の染色体Ｓk+1を設定する。
【０１８２】
【数４１】
2010136600.tif 000043
【０１８３】
最適化部６は、式（７１）、（７２）、（７３）で表わされる次の時点ｎの最適制御量ＣＯＮ（ｎ）の値を、式（７４）で示されるように、最終の反復回数Ｋにおける評価値が最小であるエリート個体の染色体ＳKと同一の値に設定する。ここで、Ｎ＿ＳＣ２０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における２０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ３０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における３０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ４０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｌ）における４０ＭＶＡの調相設備の投入台数を表わし、Ｎ＿ＳＣ６０（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）における６０ＭＶＡの調相設備の投入台数を表わす。また、ＴＡＰ（ｎ，ｊ）は、次の時点ｎの変電所ｊ（ｊ＝１〜Ｍ）におけるタップ付き変圧器３０のタップ位置（中間位置を原点とする）を表わす。
【０１８４】
【数４２】
2010136600.tif 000044
【０１８５】
（潮流計算部）
潮流計算部８は、各時点の最適制御量の探索中において、発電機Ｈをスイングノードとし、潮流計算を行なう。すなわち、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、次の時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、次の時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける個体ｉの反復回数ｋでのタップ付き変圧器３０のタップ位置ＴＡＰik（１）〜ＴＡＰik（Ｍ）、変電所１〜変電所Ｌにおける個体ｉの反復回数ｋでの２０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ２０（１）〜Ｎik＿ＳＣ２０（Ｌ）、変電所１〜変電所Ｌにおける個体ｉの反復回数ｋでの３０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ３０（１）〜Ｎik＿ＳＣ３０（Ｌ）、変電所１〜変電所Ｌにおける個体ｉの反復回数ｋでの４０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ４０（１）〜Ｎik＿ＳＣ４０（Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける個体ｉの反復回数ｋでの６０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ６０（ＬＳ＋１）〜Ｎik＿ＳＣ６０（Ｌ）に基づいて、潮流計算によってスイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。さらに、潮流計算部８は、式（７５）に従って、個体ｉの反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikを算出する。
【０１８６】
【数４３】
2010136600.tif 000045
【０１８７】
潮流計算部８は、時点ｎで最適制御量ＣＯＮ（ｎ）が決定された後、発電機Ｈをスイングノードとし、潮流計算を行なう。すなわち、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、次の時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、次の時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける次の時点ｎのタップ付き変圧器３０のタップ位置ＴＡＰ（ｎ，１）〜ＴＡＰ（ｎ，Ｍ）、変電所１〜変電所Ｌにおける次の時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，１）〜Ｎ＿ＳＣ２０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける次の時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，１）〜Ｎ＿ＳＣ３０（ｎ，Ｌ）、変電所１〜変電所ＬＭにおける次の時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，１）〜Ｎ＿ＳＣ４０（ｎ，Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける次の時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ＬＳ＋１）〜Ｎ＿ＳＣ６０（ｎ，Ｌ）に基づいて、潮流計算によって、スイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。潮流計算部８は、式（７６）に従って、次の時点ｎの電力系統の有効電力ロスＰｌｏｓｓ（ｎ）を算出する。
【０１８８】
【数４４】
2010136600.tif 000046
【０１８９】
（動作）
次に、第３の実施形態における電圧無効電力制御の動作を説明する。
【０１９０】
図１０は、第３の実施形態の電圧無効電力制御システムの動作手順を表わすフローチャートである。
【０１９１】
図１０を参照して、変電所ｊ（ｊ＝１〜Ｍ）のＰ検出部２４は、時点ｎ（ｎ＝１〜Ｎ）の有効電力Ｐ（ｎ，ｊ）を検出する。また、変電所ｊのＱ検出部２６は、時点ｎ（ｎ＝１〜Ｎ）の無効電力Ｑ（ｎ，ｊ）を検出する。変電所ｊの通信部２２は、検出された有効電力Ｐ（ｎ，ｊ）および無効電力Ｑ（ｎ，ｊ）をＶＱＣ装置１の通信部３に送信する（ステップＳ４０１）。
【０１９２】
次に、ＶＱＣ装置１において、時点ｎがＮ＋１に設定される（ステップＳ４０２）。
次に、変電所ｊ（ｊ＝１〜Ｍ）のＰ検出部２４は、時点ｎの有効電力Ｐ（ｎ，ｊ）を検出する。また、変電所ｊのＱ検出部２６は、時点ｎの無効電力Ｑ（ｎ，ｊ）を検出する。変電所ｊの通信部２２は、検出された有効電力Ｐ（ｎ，ｊ）および無効電力Ｑ（ｎ，ｊ）をＶＱＣ装置１の通信部３に送信する（ステップＳ４０３）
次に、ＶＱＣ装置１の需要予測部４は、式（４）〜式（７）に従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された有効電力の差分量を用いて、ＡＲ係数｛ａ1(j),ａ2(j),，ａs(j)｝を算出する。また、ＶＱＣ装置１の需要予測部４は、式（１１）〜式（１４）に従って、変電所ｊの現在までに計測された無効電力の差分量を用いて、ＡＲ係数｛ｂ1(j),ｂ2(j),，ｂs(j)｝を算出する（ステップＳ４０４）。
【０１９３】
次に、ＶＱＣ装置１において、時点ｎが１だけインクリメントされる（ステップＳ４０５）。
【０１９４】
次に、需要予測部４は、式（２）に示す次数ｓのＡＲモデルに従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された有効電力の差分量を用いて、変電所ｊの時点ｎの有効電力の差分量の予測値ＤＰe（ｊ）を算出する。需要予測部４は、変電所ｊの時点ｎの有効電力の差分量の予測値ＤＰe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された有効電力Ｐ（ｎ−１，ｊ）を用いて、式（３）にしたがって、変電所ｊの時点ｎの有効電力予測量Ｐe（ｊ）を算出する。
【０１９５】
また、需要予測部４は、式（９）に示す次数ｓのＡＲモデルに従って、変電所ｊ（ｊ＝１〜Ｍ）の現在までに計測された無効電力の差分量を用いて、変電所ｊの時点ｎの無効電力の差分量の予測値ＤＱe（ｊ）を算出する。需要予測部４は、変電所ｊの時点ｎの無効電力の差分量の予測値ＤＱe（ｊ）と、変電所ｊの時点（ｎ−１）に計測された無効電力Ｑ（ｎ−１，ｊ）を用いて、式（１０）にしたがって、変電所ｊの時点ｎの有効電力予測量Ｑe（ｊ）を算出する（ステップＳ４０６）。
【０１９６】
次に、ＶＱＣ装置１の最適化部６は、たとえばｉ＝６〜Ｉについて、式（６１）〜（６４）で示すような個体ｉの反復回数１での染色体Ｘi1を擬似乱数で設定する。また、最適化部６は、たとえば個体１〜個体５の反復回数１での染色体Ｘp1（ｐ＝１〜５）の値を時点（ｎ−１）の反復回数Ｋでのエリート個体の染色体ＳKと同一の値に設定する。さらに、最適化部６は、反復回数１でのエリート個体の染色体Ｓ1の値を個体ｉ（ｉ＝１〜Ｉ）の反復回数１での染色体Ｘi1のうちの評価値が最小である染色体と同一の値に設定する。最適化部６は、反復回数ｋ＝１に設定する（ステップＳ４０７）。
【０１９７】
次に、最適化部６は、式（６５）に従って、個体ｉの反復回数ｋでのタップ位置の変化量ΔＴＡＰikを算出する。また、最適化部６は、式（６６）に従って、個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの調相設備の投入台数の変化量ΔＳＣikを算出する（ステップＳ４０８）。
【０１９８】
次に、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでのタップ付き変圧器３０のタップ位置ＴＡＰik（１）〜ＴＡＰik（Ｍ）、変電所１〜変電所Ｌにおける個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの２０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ２０（１）〜Ｎik＿ＳＣ２０（Ｌ）、変電所１〜変電所Ｌにおける個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの３０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ３０（１）〜Ｎik＿ＳＣ３０（Ｌ）、変電所１〜変電所Ｌにおける個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの４０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ４０（１）〜Ｎik＿ＳＣ４０（Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの６０ＭＶＡの調相設備の投入台数Ｎik＿ＳＣ６０（ＬＳ＋１）〜Ｎik＿ＳＣ６０（Ｌ）に基づいて、潮流計算によってスイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。さらに、潮流計算部８は、式（７５）に従って、個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの電力系統の有効電力ロスＰｌｏｓｓikを算出する。最適化部６は、式（６７）に従って、個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの電力系統の有効電力ロスの変化量ΔＰlossikを算出する（ステップＳ４０９）。
【０１９９】
次に、最適化部６は、個体ｉの反復回数ｋでの染色体Ｘikの評価値ＥＶ（Ｘik）に従って、次の反復回数（ｋ＋１）での個体ｉ（ｉ＝１〜Ｉ）の染色体Ｘik+1とエリート個体の染色体Ｓk+1の設定を行なう（ステップＳ４１０）。
【０２００】
次に、最適化部６は、次の反復回数（ｋ＋１）がＫであれば（ステップＳ４１１でＹＥＳ）、ステップＳ４１３に進み、（ｋ＋１）がＫでなければ（ステップＳ４１１でＮＯ）、ステップＳ４１２に移行して、反復回数ｋを１だけ増加した上で、ステップＳ４０８に戻る。
【０２０１】
ステップＳ４１３において、最適化部６は、式（７１）、（７２）、（７３）で表わされる次の時点ｎの最適制御量ＣＯＮ（ｎ）の値を、式（７４）で示されるように、最終の反復回数Ｋにおける評価値が最小であるエリート個体の染色体ＳKと同一の値に設定する。
【０２０２】
さらに、潮流計算部８は、発電機１〜発電機（Ｈ−１）の有効電力ＰG（１）〜ＰG（Ｈ−１）、発電機１〜発電機Ｈの無効電力ＱG（１）〜ＱG（Ｈ）、時点ｎの変電所１〜変電所Ｍの有効電力予測量Ｐe（１）〜Ｐe（Ｍ）、時点ｎの変電所１〜変電所Ｍの無効電力予測量Ｑe（１）〜Ｑe（Ｍ）、変電所１〜変電所Ｍにおける時点ｎのタップ付き変圧器３０のタップ位置ＴＡＰ（ｎ，１）〜ＴＡＰ（ｎ，Ｍ）、変電所１〜変電所Ｌにおける時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，１）〜Ｎ＿ＳＣ２０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，１）〜Ｎ＿ＳＣ３０（ｎ，Ｌ）、変電所１〜変電所Ｌにおける時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，１）〜Ｎ＿ＳＣ４０（ｎ，Ｌ）、変電所（ＬＳ＋１）〜変電所Ｌにおける時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ＬＳ＋１）〜Ｎ＿ＳＣ６０（ｎ，Ｌ）に基づいて、潮流計算によって、スイングノードである発電機Ｈの有効電力ＰG（Ｈ）を求める。潮流計算部８は、式（７６）に従って、時点ｎの電力系統の有効電力ロスＰｌｏｓｓ（ｎ）を算出する。
【０２０３】
さらに、通信部３は、変電所ｊに対して、最適制御量ＣＯＮ（ｎ）のうち、変電所ｊに関連する情報である、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの２０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ２０（ｎ，ｊ）、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの３０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ３０（ｎ，ｊ）、変電所ｊ（ｊ＝１〜Ｌ）の時点ｎの４０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ４０（ｎ，ｊ）、変電所ｊ（ｊ＝（ＬＳ＋１）〜Ｌ）の時点ｎの６０ＭＶＡの調相設備の投入台数Ｎ＿ＳＣ６０（ｎ，ｊ）、および変電所ｊ（ｊ＝１〜Ｍ）のタップ付き変圧器３０の時点ｎのタップ位置ＴＡＰ（ｎ，ｊ）を送信する（ステップＳ４１３）。
【０２０４】
変電所ｊは、情報Ｎ＿ＳＣ２０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ３０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ４０（ｎ，ｊ）（ｊ＝１〜Ｌの場合）、Ｎ＿ＳＣ６０（ｎ，ｊ）（ｊ＝（ＬＳ＋１）〜Ｌの場合）、およびＴＡＰ（ｎ，ｊ）（ｊ＝１〜Ｍの場合）を受信する。変電所の制御部２８は、これらの受信した情報に基づいて、タップ付き変圧器３０のタップ位置と、調相設備群３２の投入台数を制御する（ステップＳ４１４）。
【０２０５】
（ステップＳ４１０の詳細）
図１１は、図１０のステップＳ４１０の詳細な手順を表わすフローチャートである。
【０２０６】
図１１を参照して、最適化部６は、式（６８）に従って、個体ｉ（ｉ＝１〜Ｉ）の反復回数ｋでの染色体Ｘikの評価値ＥＶ（Ｘik）を算出する（ステップＳ５０１）。
【０２０７】
次に、最適化部６は、個体ｉ（ｉ＝１〜Ｉ）について、評価値ＥＶ（Ｘik）が小さい順番に染色体Ｘikを並べる（ステップＳ５０２）。
【０２０８】
次に、最適化部６は、個体ｉ（ｉ＝１〜Ｉ）の染色体のうちの評価値ＥＶ（Ｘik）が最小の染色体Ｘjkを特定して、それをエリート個体とする（ステップＳ５０３）。
【０２０９】
次に、最適化部６は、染色体の淘汰処理を行なう。すなわち、最適化部６は、個体ｉ（ｉ＝１〜Ｉ）の染色体のうちの評価値ＥＶ（Ｘik）が最大から第１番目〜第５番目の５個の染色体を特定して、それらをエリート個体の染色体Ｘjkで置換する（ステップＳ５０４）。
【０２１０】
次に、最適化部６は、染色体の交差処理を行なう。すなわち、最適化部６は、ＮＣＲ個（交差率×個体数）の個体をランダムに選択して、選択した個体から（ＮＣＲ／２）個のペアを作成して、ランダムな位置でペアの２つの個体の染色体を交差させる（ステップＳ５０５）。
【０２１１】
次に、最適化部６は、染色体の突然変異処理を行なう。すなわち、最適化部６は、ＮＶＲ個（突然変異率×個体数）の個体をランダムに選択して、選択した各個体の染色体の１つの遺伝子をランダムに選択して、ランダムな値に変更する(ステップＳ５０６）。
【０２１２】
次に、最適化部６は、式（６９）に従って、反復回数（ｋ＋１）の個体ｉの染色体Ｘik+1の値を反復回数ｋの個体ｉの染色体Ｘikと同一の値に設定する（ステップＳ５０７）。
【０２１３】
次に、最適化部６は、個体ｉ（ｉ＝１〜Ｉ）の反復回数（ｋ＋１）での染色体Ｘik+1の評価値ＥＶ（Ｘik+1）を算出する。次に、最適化部６は、個体ｉ（ｉ＝１〜Ｉ）の染色体のうちの評価値ＥＶ（Ｘik+1）が最小の染色体Ｘjk+1を特定して、反復回数（ｋ＋１）のエリート個体の染色体Ｓk+1の値をその特定した染色体Ｘjk+1と同一の値に設定する（ステップＳ５０８）。
【０２１４】
（小括）
以上のように、本発明の実施形態の電圧無効電力制御システムによれば、将来の有効電力および無効電力を予測して、予測に基づいてＧＡによって最適制御量を探索するので、電圧制御機器の制御遅れを考慮した電圧無効電力制御を行なうことができる。
【０２１５】
［シミュレーション結果］
上記説明した本発明の第１〜第３の実施形態による電圧無効電力制御についてシミュレーションによって有効性を確かめた。
【０２１６】
制御時間間隔Δｔを１０分とし、式（２４）、（３０）、（６８）において、ｗ1を４×１０-4、ｗ2を７、ｗ3を５とした。式（１９）におけるＣ1を0.9から0.4の範囲で繰返し回数が増加するごとに減少させ、Ｃ2およびＣ3を２．０とした。
【０２１７】
図１２は、３つの手法による、ある変電所の母線電圧を表わす図である。図１３は、３つの手法による、別の変電所の母線電圧を表わす図である。図１２および図１３に示すように、いずれの変電所においても、母線電圧を制約条件（下限値以上、かつ上限値以下）を満たすように制御することができた。
【０２１８】
図１４は、実績値を用いたシミュレーションによる電力系統の有効電力ロスと、３つの手法による電力系統の有効電力ロスを表わす図である。図１４を参照して、３つの手法ででは、実績値を用いたシミュレーションに比べて、電力系統の有効電力ロスが約１％削減できた。
【０２１９】
（変形例）
本発明は、上記の実施形態に限定されるものではない。たとえば、以下のような変形例も当然に本発明に含まれる。
【０２２０】
（１）  変形例１
本発明の実施形態として、メタヒューリスティック手法の例として、ＰＳＯ、遺伝的アルゴリズム、タブーサーチについて説明したが、焼きなまし法などのような他のメタヒューリスティック手法を用いてもよい。
【０２２１】
（２）  変形例２
本発明の実施形態で説明した式（２１）、（２２）、（２３）に代えて、以下の式（７７）、（７８）、（７９）を用いることしてもよい。
【０２２２】
【数４５】
2010136600.tif 000047
【０２２３】
今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【符号の説明】
【０２２４】
１  ＶＱＣ装置、２−１〜２−ｍ  変電所１〜変電所Ｍ、３  通信部、４  需要予測部、５  変電所群、６  最適化部、７−１〜７−ｎ  通信回線、８  潮流計算部、９  連系線、１１  発電機群、１２〜１〜１２−Ｈ  発電機群、２２  通信部、２４  Ｐ検出部、２６  Ｑ検出部、２８  制御部、３０  タップ付き変圧器、３２  調相設備群、３３−１  ２０ＭＶＡ調相設備群、３３−２  ３０ＭＶＡ調相設備群、３３−３  ４０ＭＶＡ調相設備群、３３−４  ６０ＭＶＡ調相設備群。
【図１】
2010136600.tif 000048
【図２】
2010136600.tif 000049
【図３】
2010136600.tif 000050
【図４】
2010136600.tif 000051
【図５】
2010136600.tif 000052
【図６】
2010136600.tif 000053
【図７】
2010136600.tif 000054
【図８】
2010136600.tif 000055
【図９】
2010136600.tif 000056
【図１０】
2010136600.tif 000057
【図１１】
2010136600.tif 000058
【図１２】
2010136600.tif 000059
【図１３】
2010136600.tif 000060
【図１４】
2010136600.tif 000061