JP7808531B2 - 電力系統運用計画生成装置及び電力系統運用計画生成方法 - Google Patents

Description

本発明は、電力系統運用計画生成装置及び電力系統運用計画生成方法に関する。

電力系統では、温室効果ガスを排出しない再生可能エネルギー電源（以下、再エネと言う）の導入が進んでいる。代表的な再エネである太陽光発電（以下、ＰＶと言う）や風力発電は気象条件よって出力が変動する。

加えて、電気自動車（以下、ＥＶと言う）の普及を背景に、急速充電ステーションの普及や大容量化が進んでおり、この電力需要が電力系統に影響を与える可能性がある。こうした再エネ出力やＥＶ充電需要に代表される、不確実性を有するパラメータ（以下、不確実パラメータと言う）を考慮に入れた系統運用が求められている。

一方、系統運用には、前日等の事前の段階で決定すべき変数（以下、運用計画変数と言う）とリアルタイムに決定できる変数（以下、当日制御変数と言う）の２種類の変数が存在する。例えば、配電系統における運用計画変数としては、タップ調整機器（例えば、ＬＲＴ：Ｌｏａｄ Ｒａｔｉｏ ｃｏｎｔｒｏｌ Ｔｒａｎｓｆｏｍｅｒ、ＳＶＲ：Ｓｔｅｐ Ｖｏｌｔａｇｅ Ｒｅｇｕｌａｔｏｒ）タップ位置、自家用発電機や蓄電池、ＥＶなどの需要家所有の分散エネルギー源（以下、ＤＥＲと言う）の調整余力の予約量などが挙げられる。

また、当日制御変数としては、ＳＶＣ（Ｓｔａｔｉｃ Ｖａｒ Ｃｏｍｐｅｎｓａｔｏｒ）や、系統運用事業者が所有するＤＥＲ（例えば、系統用蓄電池）、予約した需要家所有ＤＥＲの調整余力の発動指令量などが挙げられる。

運用計画変数の決定方法としてまず考えられるのは、不確実パラメータの予測が当たると仮定し、系統運用ＫＰＩ（例えば、電圧許容範囲や線路容量の違反量、対応コスト）が最良となる値の組み合わせを選ぶ方法である。この方法では、予測が外れた場合は考慮されておらず、不確実パラメータの変動パターンによっては系統運用ＫＰＩの極端な悪化（例えば、電圧許容範囲や線路容量などの制約違反の発生、追加対応による極端なコスト増など）が発生する可能性がある。

そのため、不確実パラメータの予測が変動した場合も、系統運用ＫＰＩ（例えば、電圧許容範囲や線路容量の違反量、対応コスト）が極端に悪化しない（以下、「ロバストな」と表記する）ように運用計画変数を決定する必要がある。

ロバストな系統運用計画を立案する最も単純な方法は、不確実な変動に対し保守的に計画を立てる、例えば需要家所有ＤＥＲの調整余力の予約量を可能な限り大きくすることである。この方法では、確かにロバスト性が担保されるが、系統運用にかかるコストが増大するため、現実的ではない。

このことから、ロバストな系統運用計画を現実の系統運用に適用するには、予測外れに対するロバスト性を確保しつつ、系統運用コストが最小となるような系統運用計画の立案手法が求められる。

ロバストな系統運用を目的とした運用計画変数の決定技術の一つとして、特許文献１に記載される技術が知られている。特許文献１には、「電力配分装置、電力配分方法、およびプログラム」が記されている。

特許文献１では、各分散電源の過去の動作履歴に基づいて将来に出力できる電力の変動範囲を推定し、この変動範囲に基づいて予め定められた目的条件および制約条件を満たすことができる各分散電源の有効電力および無効電力を演算する。

特開２０２２－７７４５９号公報

特許文献１の最適化手法は、一般にロバスト最適化と呼ばれるものである。ロバスト最適化で求まる解における、目的関数（例えば、系統運用コスト）とロバスト性（例えば、不確実な変動に対する電力系統の安全性）のバランスは、不確実パラメータの変動範囲の大小により決定される。

即ち、変動範囲が大きければ、ロバスト性が増大する代わりに目的関数が悪化する。このことから、ロバスト最適化では変動範囲の適切な設定が肝要である。

これに対し、特許文献１における変動範囲の設定方法について、「予め定められた期間における分散電源の出力変化量の標準偏差に基づいて推定してよい」と記載されている。しかし、標準偏差に基づく推定では、不確実データの時間的・空間的変動特性を考慮せずに、変動範囲を設定してしまうことになる。

ここで、不確実データの時間的・空間的変動特性について説明する。出力範囲が０ｋＷから５００ｋＷの再生可能エネルギー機器（以下、再エネ機器）を例に考えると、ある時刻における出力が２００ｋＷのとき、次のタイムステップでの出力は１００ｋＷから３００ｋＷの範囲に収まる場合がある。また、近隣の再エネ機器の出力が２００ｋＷの場合、同時刻におけるこの再エネ機器の出力は１００ｋＷから３００ｋＷの範囲に収まる場合がある。

ところが、特許文献１ではこの時間的・空間的変動特性を考慮しないため、本来であれば出力が１００ｋＷから３００ｋＷの範囲に収まる場合も、０ｋＷから５００ｋＷを変動範囲として設定してしまい、結果として必要以上にロバスト性を重視した運用計画が出力され、系統運用コストが増大する可能性がある。

本発明の目的は、電力系統運用計画生成装置において、予測外れへのロバスト性と系統運用コスト最小化を両立した系統運用計画を求めることにある。

本発明の一態様の電力系統運用計画生成装置は、電力系統における不確実性パラメータを考慮に入れ、事前の段階で決定すべき運用計画変数とリアルタイムに決定できる当日制御変数を定め、前記電力系統の系統運用計画を生成する電力系統運用計画生成装置であって、前記電力系統運用計画生成装置は、プロセッサにより所定の演算を実行する演算装置を有し、前記演算装置は、前記プロセッサにより、不確実性を有する不確実データの過去実績とロバスト性パラメータに基づいて前記不確実データの時間当たりの変動量又は前記不確実データの空間当たりの変動量に関する不確実データ変動範囲を作成して前記不確実データ変動範囲の制約式を生成する不確実データ変動範囲作成部と、前記プロセッサにより、前記不確実データ変動範囲の前記制約式を参照して前記電力系統の前記系統運用計画を生成する系統運用計画生成部と、を有し、前記不確実データ変動範囲作成部は、前記プロセッサにより、前記不確実データを基に散布図を描き、前記不確実データの全てのデータを囲む面積最小の多角形の複数の頂点を求め、複数の前記頂点のうち一番左下の点と一番右上の点を結ぶ中心線を基準にして、複数の前記頂点を前記中心線の上側と下側で分け、前記中心線の前記上側と前記下側とで、前記不確実データ変動範囲の外枠の直線式の等号を不等号に変換した前記制約式を作成し、前記系統運用計画生成部は、前記プロセッサにより、前記制約式を参照して前記電力系統の前記系統運用計画を生成することを特徴とする。

本発明の一態様によれば、電力系統運用計画生成装置において、予測外れへのロバスト性と系統運用コスト最小化を両立した系統運用計画を求めることができる。

実施例１の電力系統運用計画生成装置の構成を示すブロック図である。 不確実データ過去実績の例を示す図である。 不確実データ変動範囲の例を示す図である。 系統運用計画の例を示す図である。 実施例１の処理フローの例を示す図である。 ロバスト性パラメータを用いた頂点移動の例を示す図である。 実施例２の電力系統運用計画生成装置の構成を示すブロック図である。 実施例２の処理フローの例を示す図である。 電力系統運用計画生成装置のハードウェア構成を示す図である。

以下に、電力系統運用計画生成装置および電力系統運用計画生成方法の実施例を図面に基づいて詳細に説明する。なお、この実施例によりこの発明が限定されるものではない。そして、各実施例は、処理内容を矛盾させない範囲で適宜組み合わせることが可能である。

実施例１は、適用事例を需要家所有ＤＥＲ、系統事業者所有ＤＥＲ、ＳＶＲ、ＬＲＴ、ＳＶＣを含む配電系統を対象とする。実施例１の電力系統運用計画生成装置は、不確実パラメータとして太陽光発電の出力を対象として、太陽光発電の出力の時間的変動特性を考慮しロバストな運用計画を策定する。

図１を参照して、実施例１の電力系統運用計画生成装置の構成について説明する。図１は、計算機装置を用いて実現される電力系統運用計画生成装置の演算部における処理内容を、その代表的な処理機能部として表記した図である。

図１に示すように、電力系統運用計画生成装置１０は、演算部２０と演算結果記憶部３０を有している。例えば、演算部２０は、プロセッサにより所定の演算を実行する演算装置である。ただし、運用計画策定装置１０は、図１に示した機能部以外にも、例えば、入力データを入力する入力部、他の端末と通信を行う通信インタフェース部などを有してもよい。

演算部２０は、不確実データ変動範囲作成部２１０と系統運用計画生成部２２０を有している。不確実データ変動範囲作成部２１０は、プロセッサにより、不確実性を有する不確実データの過去実績１１０とロバスト性パラメータ１２０に基づいて不確実データの時間当たりの変動量又は不確実データの空間当たりの変動量に関する不確実データ変動範囲３１０を作成して不確実データ変動範囲３１０の制約式情報２３０を生成する。系統運用計画生成部２２０は、プロセッサにより、制約式情報２３０を参照して電力系統の系統運用計画３２０を生成する。また、演算結果記憶部３０は情報記憶のためのデータベースなどを有しており、不確実データ変動範囲３１０と系統運用計画３２０と系統状態の将来予測３３０を記憶している。

演算部２０は、その処理遂行のために、外部からの入力として、不確実データ過去実績１１０とロバスト性パラメータ１２０と電力需要情報１３０と系統機器情報１４０と運用制約情報１５０と電力市場価格情報１６０と再エネ予測情報１７０と需要家所有ＤＥＲ情報１８０を得る。

この場合に使用する外部入力のうち不確実データ過去実績１１０を、図２を用いて説明する。

図２の上から１行目は各太陽光発電のＩＤ１１１を表し、２行目は各太陽光発電の設置場所の北緯を表し、３行目は各太陽光発電の設置場所の東経を表し、４行目以降の各行は各時間ステップにおける各太陽光発電の出力を表す。図２では例として、時間ステップの刻み幅を５分としている。

また、ロバスト性パラメータ１２０は数値データであり、ロバスト性パラメータ１２０をｐと置くと、例えば０≦ｐ≦１を満たす数値データである。

また、系統機器情報１４０は、例えば、各タップ調整機器の機器種別（例えば、ＬＲＴ、ＳＶＲなど）、設置バスＩＤ、タップ段数、１段毎のタップ幅、タップ基準位置、制御指令への応答速度、日当たりのタップ調整制限回数、設備コスト、設置コストなどを含む。

また、運用制約情報１５０は、系統運用の際に守るべき電圧の範囲や、潮流量の範囲である。

また、需要家所有ＤＥＲ情報１８０は、例えば、機器種別（例えば、自家発、電気自動車など）、設置バスＩＤ、設備容量、制御指令への応答速度などがある。また、上記需要家所有ＤＥＲ情報１８０は、個別に分散エネルギー源ＤＥＲ毎に入力しても、複数の分散エネルギー源ＤＥＲの合算値を入力してもよい。

図１に戻ると、不確実データ変動範囲作成部２１０は、外部からの入力データとして、不確実データ過去実績１１０とロバスト性パラメータ１２０を入手し、これらを基に不確実データの変動範囲とその制約式を生成し、不確実データ変動範囲３１０を演算結果記憶部３０に出力し、不確実データ変動範囲の制約式情報２３０を系統運用計画生成部２２０に出力する。

出力データのうち不確実データ変動範囲３１０を、図３を用いて説明する。

不確実データ変動範囲３１０は、時刻ｔにおける太陽光発電出力を横軸に取り、時刻ｔ＋１における太陽光発電出力を縦軸に取った場合におけるプロット３１１と変動範囲外枠３１２を含む。

プロット３１１は不確実データ過去実績１１０を基に散布される。例えば、１２時００分に２０ｋＷであった出力が、１２時０５分に４０ｋＷに変化した場合、横軸が２０ｋＷ、縦軸が４０ｋＷの点にプロット３１１が打たれる。

なお、図３では全時刻の過去実績を同一の散布図に描いているが、例えば朝、昼、夕方で散布図を分割しても良い。また、変動範囲外枠３１２は、プロット３１１の集合を全て囲む面積最小の多面体、即ち凸包を構成する辺である。

なお、図３では太陽光発電の出力の時間的変動特性を考慮するために、時刻ｔにおける太陽光発電出力を横軸に取り、時刻ｔ＋１における太陽光発電出力を縦軸に取ったが、空間的変動特性を考慮する場合は、ある太陽光発電の出力を横軸に取り、近隣の太陽光発電の同時刻における出力を縦軸に取れば良い。

また、不確実データ変動範囲の制約式情報２３０は、上記変動範囲外枠３１２の直線の式ｙ＝ａｘ＋ｂのうち、等号を不等号に変換したものである。ｘは図３の横軸の値、ｙは図３の縦軸の値、ａとｂは定数である。不等号の向きは、変動範囲外枠３１２がプロット３１１より上にある場合、ｙ≦ａｘ＋ｂである。また、変動範囲外枠３１２がプロット３１１より下にある場合、ｙ≧ａｘ＋ｂである。

再度図１に戻ると、系統運用計画生成部２２０は、外部からの入力データとして、電力需要情報１３０と系統機器情報１４０と運用制約情報１５０と電力市場価格情報１６０と再エネ予測情報１７０と需要家所有ＤＥＲ情報１８０を入手し、これらを基にロバスト最適化問題を解き、系統運用計画３２０と系統状態の将来予測３３０を演算結果記憶部３０に出力する。

出力データのうち、系統運用計画３２０を、図４を用いて説明する。

図４の各行は各時間ステップ３２１を表し、各列はタップ調節機器のＩＤ３２２および需要家所有ＤＥＲのＩＤ３２３を表す。各タップ調節機器のＩＤ３２２については、各時間ステップ３２１におけるタップ位置３２４を格納する。また、各需要家所有ＤＥＲのＩＤ３２３については、各時間ステップの調整余力の予約量２２５を格納する。調整余力の予約量３２５は需要増減双方の値を取り得るため、図４で示すように負の値も取り得る。

再度図１に戻ると、演算結果記憶部３０に記憶された情報は、計算結果表示部４１０を介して運用者・計画者４２０に表示される。また、演算結果記憶部３０に記憶された系統運用計画３２０に従い、コントローラ５１０から、需要家所有ＤＥＲ５２０、系統事業者所有ＤＥＲ５３０、ＳＶＲ５４０、ＬＲＴ５５０、ＳＶＣ５６０に制御指令が送られる。これにより、上記機器群が制御され、ロバスト性の確保と系統運用コストの減少を両立した系統運用が実現する。

続いて、実施例１の処理フローを説明する。図５に示す処理部２０におけるフローチャートに従い説明する。フローチャート内の処理のうち、処理２１１から処理２１６までは不確実データ変動範囲作成部２１０で実行される。また、処理２２１は系統運用計画生成部２２０で実行される。

処理２１１にて、不確実データ過去実績１１０を基に図３の散布図を描き、全プロットを囲む面積最小の多角形、即ち凸包の頂点を求める。これは、ギフト包装法などを適用することで求めることができる。

処理２１２にて、凸包頂点のうち一番左下の点と右上の点を結ぶ線（以下、中心線）を引く。処理２１３にて、上記中心線と各凸包頂点との距離ｄを計算する。

処理２１４にて、ロバスト性パラメータ１２０（以下、ｐとする）を上記距離ｄと掛け算した値（以下、ｄ’とする）を計算する。そして、凸包頂点を通る中心線の垂線上にあり、中心線との距離がｄ’になるよう、凸包頂点を中心線の垂線方向に移動する。以上の処理を全ての凸包頂点に対して行う。

処理２１４を図で示したのが図６である。図６における薄い灰色のプロット６０１は、図３におけるプロット３１１と同一である。プロット６０１のうち、濃い灰色のプロット６０２は、処理２１１で計算された凸包の頂点である。このうち、対象の凸包頂点６０３から、中央線６０４に向かい垂線６０５を引く。そして、中央線６０４との距離がｄ’になるよう、垂線６０５に沿って対象の凸包頂点６０３を移動する。図６では例として、ロバスト性パラメータｐが１．３のときにおける移動後の凸包頂点６０６を示している。

ここで、処理２１３と処理２１４にて、中心線と各凸包頂点との距離ｄを用いる理由を説明する。他の候補として、全プロットの中央の点を基準として、各凸包頂点との距離を用いる案もある。しかし、この場合、ロバスト性パラメータｐが１未満のときに、中央線の左下側付近のプロットおよび右上側付近のプロットが、更新後の凸包の外側になる可能性がある。中央線付近は、時刻ｔと時刻ｔ＋１の出力がほぼ同じエリアになるため、多くのプロットが集中する。この多くのプロットが集中するエリアを凸包の外側にしてしまうと、後述のロバスト最適化で考慮されない変動パターンとなり、結果として系統運用制約を違反する運用計画になる可能性がある。以上より、本発明では中心線と各凸包頂点との距離ｄを用いて、処理２１３と処理２１４を行っている。

処理２１５では、処理２１４で移動された後の凸包頂点同士を結ぶ直線の式ｙ＝ａｘ＋ｂを求める。そして処理２１６では、上記直線の式の等号を不等号に変換する。不等号の向きは、上記直線がプロット３１１より上にある場合、ｙ≦ａｘ＋ｂである。また、上記直線がプロット３１１より下にある場合、ｙ≧ａｘ＋ｂである。

このように、不確実データ変動範囲作成部２１０は、不確実データを基に散布図を描き、不確実データの全てのデータを囲む面積最小の多角形の複数の頂点を求める。

複数の頂点のうち一番左下の点と一番右上の点を結ぶ中心線を基準にして、複数の頂点を中心線の上側と下側で分ける。そして、中心線の上側と下側とで、不確実データ変動範囲の外枠の直線式の等号を不等号に変換した制約式を作成する。

中心線の上側の場合は、不確実データ変動範囲の外枠の直線式の等号を小なりイコールの不等号（≦）に変換する。中心線の下側の場合は、不確実データ変動範囲の外枠の直線式の等号を大なりイコールの不等号（≧）に変換する。

また、不確実データ変動範囲作成部２１０は、ロバスト性パラメータに基づいて、頂点と中心線との距離に基づいて不確実データ変動範囲の面積を拡大及び縮小して調節する。

このように、凸包を構成する線を制約式に変換して、ロバスト性パラメータにより凸包の範囲を調節して数理計画問題の制約式に追加して解く。

最後に、処理２２１にてロバスト最適化問題を定式化し、求解する。この最適化問題の目的関数として、系統運用コストや送電ロスなどの重み付け和として定義される系統運用ＫＰＩを採用する場合、目的関数は以下の（式１）で示す式になる。

Φ^ＫＰＩは系統運用系統運用ＫＰＩ、ｔは時間ステップ、Ｔはｔの集合即ち計画期間を表す。また、ｍｉｎとｍａｘの下付き文字は決定変数を表し、ｘは運用計画変数のベクトル（例えば、ＳＶＲのタップ位置制御量、需要家所有ＤＥＲの調整余力予約量）、ｙは当日制御変数のベクトル（例えば、ＳＶＣの出力、予約した調整余力の発動指令量）をそれぞれ表す。

また、上記ロバスト最適化問題の制約条件は、電圧や潮流量の許容範囲と、系統機器の運用可能範囲と、需要家ＤＥＲの調整可能量の上下限、不確実データ変動範囲作成部２１０で作成した不確実データ変動範囲の制約式情報２３０などである。

上記ロバスト最適化問題の特徴は、目的関数がｍｉｎおよびｍａｘの３層構造となっている点である。（式１）の外側から順に、系統運用ＫＰＩを最小化（最良化）する運用計画ｘを決定する部分、系統運用ＫＰＩを最大化（最悪化）する不確実パラメータｕの値（最過酷シナリオ）を決定する部分、系統運用ＫＰＩを最小化する当日制御ｙを決定する部分である。

上記３層構造のため、このロバスト最適化問題はそのまま解くことができない。そこで、ベンダーズ分解などの分解手法を用いて３つの最適化問題に分割して解く。分割後の最適化問題は混合整数計画（以下、ＭＩＬＰと言う）や、混合整数二次錐計画（以下、ＭＩＳＯＣＰと言う）などに帰着するため、ＧｕｒｏｂｉやＣＰＬＥＸ等の数理計画ソルバーで解くことができる。

このようにして、上記ロバスト最適化を解くことで、計画期間Ｔにおける系統運用ＫＰＩを最小（最良）にする運用計画ｘが得られる。

図７を参照して、実施例２の電力系統運用計画生成装置の構成について説明する。

実施例２を示す図７が実施例１を示す図１と異なる点は、処理部２０の中の不確実データ変動範囲７１０の入力データとして、ノイズデータ除去数７２０が追加されている点である。他の機能およびデータは図１と同一であるのでその説明は省略する。

ノイズデータ除去数７２０は、０以上の整数で表される値であり、不確実データの過去実績の中からノイズデータとして除去し、変動範囲作成に寄与しないデータの数を表す。なお、ノイズデータ除去数７２０は、不確実データの過去実績のデータ数の割合などで与えても良い。

続いて、図８を参照して、実施例２の処理フローを説明する。

フローチャート内の処理のうち、処理２１１および処理２１２から処理２２１まで図５に示す処理と同一であり、処理８１１と処理８１２と処理８１３が図５から追加された処理である。また、フローチャート開始時に数値データａを０に初期化した状態で保持しておく。

処理８１１では、上記ノイズデータ除去数７２０をｘとしたとき、ａ＜ｘを満たすかどうかをチェックする。満たす場合は処理８１２に移動し、満たさない場合は処理２１２に移動する。

処理８１２では、処理２１１で求めた頂点の数をｂとし、ａの値をａ＋ｂに更新する。

処理８１３では、ｂとｘ－ａの値の大小を比較し、小さい方をｃに代入する。これにより計算したｃは、不確実データから除去する頂点の数を示す。なお、除去する頂点の選び方は、ランダム選出や、最近傍点との距離の大きい順に選出などが考えられる。

処理８１３の後は再度処理２１１に移動する。以上の処理を繰り返すことで、不確実データから除去したデータ数ａがノイズデータ除去数ｘに達するまで、繰り返し凸包の頂点が除去される。

図９は、電力系統運用計画生成装置１０のハードウェア構成を示す図である。

図１に示す電力系統運用計画生成装置１０は、例えば、図９に示すように、入力装置９１０、出力装置９２０、メモリー９３０、記憶装置９４０及びＣＰＵ９５０を有する計算機で構成される。

また、図１、図７に示す機能である「～部」は、例えば、プロセッサ（ＣＰＵ９５０等）によりプログラムを実行することによりその「機能」が実現される。

例えば、図１、図７に示す不確実データ変動範囲作成部２１０、７１０は、プロセッサ（ＣＰＵ９５０等）によりプログラムを実行することにより不確実データ変動範囲作成機能が実現される。

また、図１、図７に示す系統運用計画生成部２２０は、プロセッサ（ＣＰＵ９５０等）によりプログラムを実行することにより系統運用計画生成機能が実現される。

上記実施例によれば、電力系統運用計画生成装置において、予測外れへのロバスト性と系統運用コスト最小化を両立した系統運用計画を求めることができる。

１０ 電力系統運用計画生成装置
２０ 演算部
３０ 演算結果記憶部
２１０ 不確実データ変動範囲作成部
２２０ 系統運用計画生成部
４１０ 計算結果表示部
５１０ コントローラ
５２０ 需要家所有ＤＥＲ
５３０ 系統事業者所有ＤＥＲ
５４０ ＳＶＲ
５５０ ＬＲＴ
５６０ ＳＶＣ
７１０ 不確実データ変動範囲作成部

Claims (6)

1. 電力系統における不確実性パラメータを考慮に入れ、事前の段階で決定すべき運用計画変数とリアルタイムに決定できる当日制御変数を定め、前記電力系統の系統運用計画を生成する電力系統運用計画生成装置であって、
   前記電力系統運用計画生成装置は、
   プロセッサにより所定の演算を実行する演算装置を有し、
   前記演算装置は、
   前記プロセッサにより、不確実性を有する不確実データの過去実績とロバスト性パラメータに基づいて前記不確実データの時間当たりの変動量又は前記不確実データの空間当たりの変動量に関する不確実データ変動範囲を作成して前記不確実データ変動範囲の制約式を生成する不確実データ変動範囲作成部と、
   前記プロセッサにより、前記不確実データ変動範囲の前記制約式を参照して前記電力系統の前記系統運用計画を生成する系統運用計画生成部と、を有し、
   前記不確実データ変動範囲作成部は、前記プロセッサにより、
   前記不確実データを基に散布図を描き、前記不確実データの全てのデータを囲む面積最小の多角形の複数の頂点を求め、
   複数の前記頂点のうち一番左下の点と一番右上の点を結ぶ中心線を基準にして、複数の前記頂点を前記中心線の上側と下側で分け、
   前記中心線の前記上側と前記下側とで、前記不確実データ変動範囲の外枠の直線式の等号を不等号に変換した前記制約式を作成し、
   前記系統運用計画生成部は、前記プロセッサにより、
   前記制約式を参照して前記電力系統の前記系統運用計画を生成することを特徴とする電力系統運用計画生成装置。
2. 前記不確実データ変動範囲作成部は、前記プロセッサにより、
   前記中心線の前記上側の場合は、前記不確実データ変動範囲の外枠の前記直線式の等号を小なりイコールの不等号に変換し、
   前記中心線の前記下側の場合は、前記不確実データ変動範囲の外枠の前記直線式の等号を大なりイコールの不等号に変換することを特徴とする請求項１に記載の電力系統運用計画生成装置。
3. 前記不確実データ変動範囲作成部は、
   前記プロセッサにより、前記ロバスト性パラメータに基づいて前記頂点と前記中心線との距離に基づいて前記不確実データ変動範囲の面積を拡大及び縮小して調節し、
   前記系統運用計画生成部は、
   前記プロセッサにより、ロバスト最適化問題を定式化して求解することを特徴とする請求項２に記載の電力系統運用計画生成装置。
4. 電力系統における不確実性パラメータを考慮に入れ、事前の段階で決定すべき運用計画変数とリアルタイムに決定できる当日制御変数を定め、前記電力系統の系統運用計画を生成する電力系統運用計画生成方法であって、
   プロセッサにより、不確実性を有する不確実データの過去実績とロバスト性パラメータに基づいて前記不確実データの時間当たりの変動量又は前記不確実データの空間当たりの変動量に関する不確実データ変動範囲を作成して前記不確実データ変動範囲の制約式を生成する不確実データ変動範囲作成ステップと、
   前記プロセッサにより、前記不確実データ変動範囲の前記制約式を参照して前記電力系統の前記系統運用計画を生成する系統運用計画生成ステップと、を有し、
   前記不確実データ変動範囲作成ステップは、前記プロセッサにより、
   前記不確実データを基に散布図を描き、前記不確実データの全てのデータを囲む面積最小の多角形の複数の頂点を求め、
   複数の前記頂点のうち一番左下の点と一番右上の点を結ぶ中心線を基準にして複数の前記頂点を前記中心線の上側と下側で分け、
   前記中心線の前記上側と前記下側とで、前記不確実データ変動範囲の外枠の直線式の等号を不等号に変換した前記制約式を作成し、
   前記系統運用計画生成ステップは、前記プロセッサにより、
   前記制約式を参照して前記電力系統の前記系統運用計画を生成することを特徴とする電力系統運用計画生成方法。
5. 前記不確実データ変動範囲作成ステップは、前記プロセッサにより、
   前記中心線の前記上側の場合は、前記不確実データ変動範囲の外枠の前記直線式の等号を小なりイコールの不等号に変換し、
   前記中心線の前記下側の場合は、前記不確実データ変動範囲の外枠の前記直線式の等号を大なりイコールの不等号に変換することを特徴とする請求項４に記載の電力系統運用計画生成方法。
6. 前記不確実データ変動範囲作成ステップは、前記プロセッサにより、
   前記ロバスト性パラメータに基づいて、前記頂点と前記中心線の距離に基づいて前記不確実データ変動範囲の面積を拡大及び縮小して調節し、
   前記系統運用計画生成ステップは、前記プロセッサにより、
   ロバスト最適化問題を定式化して求解することを特徴とする請求項４に記載の電力系統運用計画生成方法。