(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024166962
(43)【公開日】2024-11-29
(54)【発明の名称】需給計画装置および需給計画方法
(51)【国際特許分類】
H02J 3/00 20060101AFI20241122BHJP
G06Q 50/06 20240101ALI20241122BHJP
【ＦＩ】
H02J3/00 170
G06Q50/06
【審査請求】未請求
【請求項の数】11
【出願形態】ＯＬ
(21)【出願番号】P 2023083410
(22)【出願日】2023-05-19
(71)【出願人】
【識別番号】000005108
【氏名又は名称】株式会社日立製作所
(74)【代理人】
【識別番号】110001807
【氏名又は名称】弁理士法人磯野国際特許商標事務所
(72)【発明者】
【氏名】山崎 潤
(72)【発明者】
【氏名】高山 厚
(72)【発明者】
【氏名】石川 倫明
(72)【発明者】
【氏名】黒田 英佑
(72)【発明者】
【氏名】桐原 健太
(72)【発明者】
【氏名】今井 秀岳
(72)【発明者】
【氏名】板井 準
【テーマコード（参考）】
5G066
5L049
5L050
【Ｆターム（参考）】
5G066AA02
5G066AA03
5G066AE09
5G066HB01
5G066HB06
5L049CC06
5L050CC06
(57)【要約】
【課題】実際の系統運用状態を反映した運用制約条件を加味した網羅的な検証シナリオの生成と処理時間の低減を両立することができる需給計画装置を提供する。
【解決手段】各種の運用制約条件を満たした上で運用計画を策定する需給計画装置２００であって、過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、過去実績情報と過去制約情報を比較評価する過去実績評価部と、過去実績評価部の過去実績評価結果に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成部１７と、を備える。需給計画装置２００は、さらに、過去運用実績を反映した過去時点における運用制約を含む推定制約情報を算出し、算出された過去時点における推定制約情報に基づき、将来時点の推定制約情報を算出する制約条件推定部１５を備える。
【選択図】図１
【解決手段】各種の運用制約条件を満たした上で運用計画を策定する需給計画装置２００であって、過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、過去実績情報と過去制約情報を比較評価する過去実績評価部と、過去実績評価部の過去実績評価結果に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成部１７と、を備える。需給計画装置２００は、さらに、過去運用実績を反映した過去時点における運用制約を含む推定制約情報を算出し、算出された過去時点における推定制約情報に基づき、将来時点の推定制約情報を算出する制約条件推定部１５を備える。
【選択図】図１
【特許請求の範囲】
【請求項1】
各種の運用制約条件を満たした上で運用計画を策定する需給計画装置であって、
過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、前記過去実績情報と前記過去制約情報を比較評価する過去実績評価部と、
前記過去実績評価部の過去実績評価結果に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成部と、を備える
ことを特徴とする需給計画装置。
過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、前記過去実績情報と前記過去制約情報を比較評価する過去実績評価部と、
前記過去実績評価部の過去実績評価結果に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成部と、を備える
ことを特徴とする需給計画装置。
【請求項2】
前記過去実績情報は、指定された過去情報のデータ期間、過去情報の取得日時、過去情報の対象設備、過去情報の対象エリアのいずれか1つ以上の条件を満たす過去情報を含む
ことを特徴とする請求項１に記載の需給計画装置。
ことを特徴とする請求項１に記載の需給計画装置。
【請求項3】
前記過去実績評価部は、指定期間におけるｋＷｈ制約、指定期間における計画停止制約のいずれか１つ以上の制約項目を含む制約条件について前記過去実績情報と前記過去制約情報を比較評価する
ことを特徴とする請求項１に記載の需給計画装置。
ことを特徴とする請求項１に記載の需給計画装置。
【請求項4】
前記過去実績評価部は、前記過去実績情報から算出された平均値、上限値、下限値、度数分布のいずれか1つ以上を含む統計的な評価指標値を用いて前記過去実績情報と前記過去制約情報を比較評価する
ことを特徴とする請求項１に記載の需給計画装置。
ことを特徴とする請求項１に記載の需給計画装置。
【請求項5】
前記需給計画装置は、さらに、前記過去実績評価部において算出した過去実績情報の評価値と前記過去制約情報に基づき、前記過去運用実績を反映した過去時点における運用制約を含む推定制約情報を算出する制約条件推定部を備え、
前記シナリオ生成部は、前記制約条件推定部で算出された推定制約情報に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成する
ことを特徴とする請求項１に記載の需給計画装置。
前記シナリオ生成部は、前記制約条件推定部で算出された推定制約情報に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成する
ことを特徴とする請求項１に記載の需給計画装置。
【請求項6】
前記制約条件推定部は、算出された過去時点における推定制約情報に基づき、将来時点の推定制約情報を算出する
ことを特徴とする請求項５に記載の需給計画装置。
ことを特徴とする請求項５に記載の需給計画装置。
【請求項7】
前記制約条件推定部は、制約項目の過去実績情報との適合度を用いて、過去時点における運用制約を含む推定制約情報を算出する
ことを特徴とする請求項５に記載の需給計画装置。
ことを特徴とする請求項５に記載の需給計画装置。
【請求項8】
前記需給計画装置は、前記過去時点および前記将来時点の推定制約情報に基づき、各制約項目の設定値の増減における上下限および変化幅、制約の考慮対象となる期間の増減における開始時間、終了時間、および期間長の上下限および変化幅、制約の考慮対象となる期間の年度、季節、曜日、平休日種別、時刻のいずれか1つ以上を含む制約変更情報を生成する制約変更条件生成部をさらに備え、
前記シナリオ生成部は、前記推定制約情報及び前記制約変更情報に基づき、前記推定制約情報及び前記制約変更情報を加味した需給計画を策定する
ことを特徴とする請求項６に記載の需給計画装置。
前記シナリオ生成部は、前記推定制約情報及び前記制約変更情報に基づき、前記推定制約情報及び前記制約変更情報を加味した需給計画を策定する
ことを特徴とする請求項６に記載の需給計画装置。
【請求項9】
前記需給計画装置は、前記検証シナリオに含まれる制約情報に基づき、前記制約情報を加味した需給計画を策定する計画策定部をさらに備える
ことを特徴とする請求項６に記載の需給計画装置。
ことを特徴とする請求項６に記載の需給計画装置。
【請求項10】
各種の運用制約条件を満たした上で運用計画を策定する需給計画方法であって、
過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、前記過去実績情報と前記過去制約情報を比較評価する過去実績評価ステップと、
前記過去実績評価ステップにおいて算出した過去実績情報の評価値と前記過去制約情報に基づき、前記過去運用実績を反映した過去時点における運用制約を含む推定制約情報を算出する制約条件推定ステップと、
前記制約条件推定ステップで算出された推定制約情報に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成ステップと、
前記検証シナリオに含まれる制約情報に基づき、前記制約情報を加味した需給計画を策定する計画策定ステップと、を有する
ことを特徴とする需給計画方法。
過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、前記過去実績情報と前記過去制約情報を比較評価する過去実績評価ステップと、
前記過去実績評価ステップにおいて算出した過去実績情報の評価値と前記過去制約情報に基づき、前記過去運用実績を反映した過去時点における運用制約を含む推定制約情報を算出する制約条件推定ステップと、
前記制約条件推定ステップで算出された推定制約情報に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成ステップと、
前記検証シナリオに含まれる制約情報に基づき、前記制約情報を加味した需給計画を策定する計画策定ステップと、を有する
ことを特徴とする需給計画方法。
【請求項11】
前記需給計画方法は、さらに、前記推定制約情報に基づき、各制約項目の設定値の増減における上下限および変化幅、制約の考慮対象となる期間の増減における開始時間、終了時間、および期間長の上下限および変化幅、制約の考慮対象となる期間の年度、季節、曜日、平休日種別、時刻のいずれか1つ以上を含む制約変更情報を生成する制約変更条件生成ステップを有し、
前記シナリオ生成ステップにおいて、前記推定制約情報及び前記制約変更情報に基づき、前記推定制約情報及び前記制約変更情報を加味した需給計画を策定する
ことを特徴とする請求項１０に記載の需給計画方法。
前記シナリオ生成ステップにおいて、前記推定制約情報及び前記制約変更情報に基づき、前記推定制約情報及び前記制約変更情報を加味した需給計画を策定する
ことを特徴とする請求項１０に記載の需給計画方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、需給計画装置および需給計画方法に関する。
【背景技術】
【0002】
電力系統は、多くの発電機、負荷、送配電機器および制御装置から構成されている。電力系統を安定かつ低コストで運用するため、時々刻々と変動する電力需要に対して、電力の供給量のバランスを保つための需給計画の策定が必要となる。天候に依存して出力が大きく変動する再生可能エネルギー電源を大量に導入しつつ高効率な系統運用を実現するため、各国において市場主導型の電力系統運用への移行が進められている。
【0003】
電力系統運用を実現するための必要機能として、次世代の中央給電指令所におけるセキュリティ制約付き発電機起動停止計画（ＳＣＵＣ：Security Constrained Unit Commitment）およびセキュリティ制約付き経済配分（ＳＣＥＤ：Security Constrained Economic Dispatch）機能があり、広域メリットオーダーと各種運用制約を充足する発電起動停止計画および出力配分を策定することが可能である。前記のＳＣＵＣ／ＳＣＥＤ機能において発電計画を求めるための最適化問題として、目的関数は運用コスト最小化、遵守すべき運用制約は発電ユニット制約（出力上下限、出力変化速度、最小起動停止時間など）、系統制約（潮流制約など）が設定される。
【0004】
今般、燃料備蓄量不足や想定外の需要変動などを原因とした需給逼迫（ｋＷｈ制約）や、気象急変による再生可能エネルギーの長期出力低下（ＤＤＤ：Dark Doldrums、または風力発電タービンの凍結など）が発生している。また、再生可能エネルギー、分散型電源、ＥＶ、ＢＥＳＳの導入量増大に伴い、電力消費傾向・変動トレンドの事前想定が困難となり、かつ事故時の影響が複雑化していく傾向にある。
【0005】
前記２点を背景として、系統の長期安定運用を担う責任者（日本国内では一般送配電事業者）が長期的な安定運用への投資と説明責任を果たすため、月間・年間以上の長期制約を検討範囲とする計画策定、および各種の想定外変動を加味した複数の検証シナリオを考慮した計画策定による事前対策検証が必要となる。ここで検証シナリオとは、電力需給状況や設備運用特性などに基づき、電力系統の運用において満たすべき制約条件を定義した情報を指す。上記の検証シナリオにより定義される制約条件を満たしつつ、コスト最小化などの指標に基づく最適な需給計画が策定される。
【0006】
前記の背景の元、電力系統を運用する事業者が様々な生じ得る系統状態を想定した上で最適な電力系統の運用方法を決定するための需給計画技術が考案されている。例えば、特許文献１には、「電力系統の設備計画の策定のために評価される複数の評価対象設備から１つの設備を選択設備として選択し、電力系統の系統構成を示す系統構成情報を用いて選択設備の存在の有無を変更した系統構成を示す変更系統構成情報を作成する構成変更反映部、を備える。優先度算出装置は、変更系統構成情報と負荷量を用いて潮流計算を行う潮流計算部と、潮流計算の結果を用いて、電力系統を構成する設備であって選択設備以外の設備である他設備への影響の度合いを選択設備の更新の優先度として算出する優先度算出部と、を備え、優先度算出装置は、選択設備を変更することにより、複数の評価対象設備のそれぞれの優先度を算出する。」という記載がある。
【先行技術文献】
【特許文献】
【0007】
【特許文献1】特開２０２１－１９１１１２号公報
【発明の概要】
【発明が解決しようとする課題】
【0008】
近年、太陽光発電や風力発電に代表される再生可能エネルギーが系統に大量導入されつつある。それらの再生可能エネルギーは急峻かつ事前想定困難な発電出力変動を伴い、この結果として、電力需給状況が大きく変化する可能性がある。また、火力発電機の運用においては、天然ガスや石油などの資源産出国の政治情勢、燃料輸送や貯蔵施設におけるトラブル、燃料価格高騰・調達困難を起因としてｋＷｈ制約の観点で発電可能量に影響が生じる事象が発生している。さらに、気象予報から乖離した想定外の気象状況による電力需要増加に発電機の作業停止の時期が重なることにより、供給力不足による需給逼迫が発生する可能性がある。
【0009】
前記の系統運用状況を考慮し、網羅的なシナリオ検証に基づく事前対策とともに、実際の運用状態に伴い今後の系統運用において想定されるシナリオを随時更新する必要がある。ここで、機械的に検証シナリオを設定すると、目的に対してシナリオが過剰になることによる計画策定の計算時間増大、またはシナリオ不足により重要シナリオの検証漏れが生じる恐れがある。以上の観点より、検証目的および実際の運用状態を反映した過不足のない検証シナリオの生成・抽出が必要となる。
【0010】
特許文献１に記載の技術では、ある１つの設備の故障有無に対する影響を評価するため、複数設備故障の組合せや設備故障以外の過酷事象など、事前に規定されていない系統状態の評価が困難であり、または評価ケースの入力の労力が増大するため、網羅的なシナリオ検証の実現が困難である。また、シナリオ設定において実際の運用状態に基づく制約条件を反映しないため、系統運用状態の変化により検証が必要となるシナリオの検証漏れや、実情を反映しない不要なシナリオ検証に計算リソースを消費する恐れがある。
【0011】
本発明の目的は、実際の系統運用状態を反映した運用制約条件を加味した網羅的な検証シナリオの生成と処理時間の低減を両立することができる需給計画装置および需給計画方法を提供することにある。
【課題を解決するための手段】
【0012】
前記目的を達成するために、本発明の需給計画装置は、各種の運用制約条件を満たした上で運用計画を策定する需給計画装置であって、過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、前記過去実績情報と前記過去制約情報を比較評価する過去実績評価部と、前記過去実績評価部の過去実績評価結果に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成部と、を備えることを特徴とする。本発明のその他の態様については、後記する実施形態において説明する。
【発明の効果】
【0013】
本発明によれば、実際の系統運用状態を反映した運用制約条件を加味した網羅的な評価シナリオの生成と処理時間の低減を両立することができる。
【図面の簡単な説明】
【0014】
【図1】本実施形態に係る需給計画装置の機能構成の一例を示す図である。
【図2】本実施形態に係る需給計画システムのハードウェア構成の一例を示す図である。
【図3】検証シナリオ生成処理の全体処理の一例を示す図である。
【図4】制約情報の一例を示す図である。
【図5A】制約拡張情報の一例（制約対象設備の増減、制約量の増減を対象とする場合）を示す図である。
【図5B】制約拡張情報の一例（制約対象期間の増減を対象とする場合）を示す図である。
【図6】過去実績情報の一例を示す図である。
【図7A】過去計画情報の一例（複数の発電機に対する各時間インターバルの発電量）を示す図である。
【図7B】過去制約情報の一例（各時間インターバルにおいて各制約項目の制約条件が有効あるいは無効のいずれであるか）を示す図である。
【図8】過去実績情報の評価指標情報の一例を示す図である。
【図9】過去時点における推定制約情報の一例を示す図である。
【図10A】将来期間における推定制約情報の一例を示す図である。
【図11A】将来時点における制約変更情報を示す図であり、制約対象設備の増減、あるいは制約量の増減を対象とする場合の制約拡張情報の一例を示す図である。
【図11B】将来時点における制約変更情報を示す図であり、制約対象期間の増減を対象とする場合の制約拡張情報の一例を示す図である。
【図12】検証シナリオ情報の一例を示す図である。
【発明を実施するための形態】
【0015】
以下、本発明の実施形態について説明する。尚、下記はあくまでも実施の例に過ぎず、下記具体的内容に発明自体が限定されることを意図する趣旨ではない。
本発明の実施形態について、図面を参照しながら以下説明する。
本発明の実施形態について、図面を参照しながら以下説明する。
【0016】
図１は、本実施形態に係る需給計画装置２００の機能構成の一例を示す図である。需給計画装置２００は、需給計画決定部１０と、計画策定部１８と、情報記憶部３０と、を含む。需給計画決定部１０は、計測装置２１（図２参照）に接続されていても良い。計測装置２１は、複数の地点の計測装置２１ａ、２１ｂ、・・・の計測情報を計測する。情報記憶部３０は、制約情報ＤＢ３１、制約拡張情報ＤＢ３２、過去実績ＤＢ３３、過去計画ＤＢ３４を格納する。
【0017】
需給計画決定部１０は、制約情報取得部１１と、過去実績取得部１２と、過去計画取得部１３と、過去実績評価部１４と、制約条件推定部１５と、制約変更条件生成部１６と、シナリオ生成部１７と、を含む。
【0018】
制約情報取得部１１は、制約情報ＤＢ３１、制約拡張情報ＤＢ３２から、制約情報、制約拡張情報を取得する。過去実績取得部１２は、過去実績ＤＢ３３から、過去実績情報を取得する。過去計画取得部１３は、過去計画ＤＢ３４より過去の発電計画および過去の発電計画に関連する制約情報のいずれかを含む過去情報を取得する。
【0019】
なお、過去実績情報、過去計画情報として、各種発電機における発電出力、または送電線、変電所を含む各種系統設備における電力潮流のいずれか１つ以上を含むものとする。
【0020】
以下、発電機における発電出力を過去実績情報および過去計画情報として用いる実施例について記載する。
【0021】
過去実績評価部１４は、過去実績取得部１２において取得した過去実績情報に基づき、過去の発電実績に基づく制約情報およびその統計情報を評価値として算出する。過去実績評価部１４は、過去実績取得部１２において取得した過去実績情報および過去計画取得部１３において取得した過去計画情報に基づき、過去の発電実績に基づく制約情報およびその統計情報を算出し、過去計画情報に含まれる過去の発電計画および過去の発電計画に関連する制約情報との比較により差分を算出する。
【0022】
すなわち、過去実績評価部１４は、過去時点における運用実績の情報を含む過去実績情報および過去時点において設定された制約情報の情報を含む過去制約情報に基づき、過去実績情報と過去制約情報を比較評価する。過去実績評価部１４は、指定期間におけるｋＷｈ制約、指定期間における計画停止制約のいずれか１つ以上の制約項目を含む制約条件について過去実績情報と過去制約情報を比較評価する。
【0023】
制約条件推定部１５は、過去実績評価部１４において算出した過去実績情報の評価値と過去計画情報に含まれる過去の発電計画および過去の発電計画に関連する制約情報を用いて、過去の運用実績における制約情報を推定する。また、前記の過去の運用実績における制約情報の推定結果を用いて、将来における制約情報の推定値を算出する。
【0024】
すなわち、制約条件推定部１５は、過去実績評価部１４において算出された過去実績情報の評価値と過去制約情報に基づき、過去運用実績を反映した過去時点における運用制約を含む推定制約情報を算出する。制約条件推定部１５は、算出された過去時点における推定制約情報に基づき、将来時点の推定制約情報を算出する。
【0025】
制約変更条件生成部１６は、制約条件推定部１５において算出した将来における制約情報の推定値および制約拡張情報に基づき、制約変更情報を生成する。
【0026】
シナリオ生成部１７は、制約情報取得部１１が制約情報ＤＢ３１および制約拡張情報ＤＢ３２より取得した制約情報および制約拡張情報、制約条件推定部１５により算出した推定制約情報、および制約変更条件生成部１６により算出した制約変更情報に基づいて、制約条件項目ごとの設定値の組み合わせにより定義されるシナリオを生成する。
【0027】
計画策定部１８は、シナリオ生成部１７において生成された前記検証シナリオに含まれる制約情報を加味した需給計画を策定する。ここで、計画策定手法として、セキュリティ制約付き発電機起動停止計画（ＳＣＵＣ：Security Constrained Unit Commitment）およびセキュリティ制約付き経済配分（ＳＣＥＤ：Security Constrained Economic Dispatch）機能と呼ばれる既知の数理計画手法を適用し、目的関数として運用コスト最小化、遵守すべき運用制約として発電ユニット制約（出力上下限、出力変化速度、最小起動停止時間）、系統制約（潮流制約）を設定して良い。また、適用する目的関数および運用制約は左記に限定されるものではなく、現状の電力系統運用において適用される各種の目的関数および運用制約を設定して良い。
【0028】
図２は、本実施形態に係る需給計画システムのハードウェア構成の一例を示す図である。図２に示すように、本実施形態における需給計画システムは、発電機１０１と、変電所１０２と、調相機器１０３と、電力負荷１０４と、外部電力系統１０５と、計測装置２１ａ，２１ｂ，・・・と、情報通信ネットワーク１０８と、需給計画装置２００と、を含む。以下、発電機１０１と、変電所１０２と、調相機器１０３と、電力負荷１０４と、計測装置２１ａ，２１ｂ，・・・を含む電力系統を自系統と呼ぶことがある。
【0029】
発電機１０１は、電力を発生させる発電機であり、火力発電、水力発電、原子力発電、太陽光発電、風力発電、バイオマス発電、潮流発電を含む何れかの発電手法により電力を発生させる発電機である。
【0030】
発電機１０１ａは、電力系統において変電所１０２の高電圧側に設置される火力発電、水力発電、原子力発電などを含む大規模発電機であり、計測装置２１ａおよび情報通信ネットワーク１０８を通じて発電量を含む系統状態量を需給計画装置２００に送信する。また、発電機１０１ａは、計測装置２１ａおよび情報通信ネットワーク１０８を通じて需給計画装置２００より送信された制御指令情報を受信し、制御指令情報に応じて発電量を含む系統状態量を変化させる。
【0031】
発電機１０１ｂは、電力系統において変電所１０２の低電圧側に設置される太陽光発電、風力発電、コジェネレーションなどを含む中小規模発電機であり、計測装置２１ｂおよび情報通信ネットワーク１０８を通じて発電量を含む系統状態量を需給計画装置２００に送信する。
【0032】
変電所１０２は、電力系統内の送電線間に設置され、大規模発電機である発電機１０１ａが設置される高電圧側より送電される電力の電圧値を変更し、電力負荷１０４が設置されている低電圧側に送電する。変電所１０２には、電力コンデンサ、分路リアクトルなどの調相機器１０３が接続される。
【0033】
調相機器１０３は、電力系統内における無効電力を変化させることにより、電力系統内の電圧分布を制御する機器であり、電力コンデンサ、分路リアクトル、ＳＴＡＴＣＯＭ（Static Synchronous Compensator：自励式無効電力補償装置）、ＳＶＣ（Static Var Compensator：静止型無効電力補償装置）、などを含む。一部の調相機器１０３は、計測装置２１ｃおよび情報通信ネットワーク１０８を通じて需給計画装置２００より送信された制御指令情報を受信し、制御指令情報に応じて発電量を含む系統状態量を変化させる。
【0034】
電力負荷１０４は、電力を消費する電動機や照明器具などを含む施設であり、家庭、工場、ビルなどの施設を表す。
【0035】
外部電力系統１０５は、需給計画装置２００からの制御が不可能である外部の電力系統であり、連系線により自系統と接続されている。
【0036】
計測装置２１ａ，２１ｂ，・・・は、発電機１０１ａにおける発電量、調相機器１０３における調相量、送電線における潮流値、電圧値、などの系統状態量を計測するセンサを含み、計測した系統状態量を、情報通信ネットワーク１０８を通じて需給計画装置２００に送信する。
【0037】
情報通信ネットワーク１０８は、双方向にデータを伝送可能なネットワークである。情報通信ネットワーク１０８は、例えば、有線ネットワーク若しくは無線ネットワーク、又はそれらの組合せである。情報通信ネットワーク１０８は、いわゆるインターネットであっても良いし、専用線のネットワークであっても良い。
【0038】
需給計画装置２００は、図１に示す需給計画機能を実現するための装置である。需給計画装置２００は、計測装置２１ａ，２１ｂ，・・・において計測された系統状態量を、情報通信ネットワーク１０８を通じて受信する。また、需給計画装置２００は、送信された系統の系統状態量および内部に蓄積された情報を用いて算出した制御指令情報を、情報通信ネットワーク１０８を通じて計測装置２１に送信する。
【0039】
需給計画装置２００の内部構成として、ＣＰＵ（Central Processing Unit）２０１と、表示装置２０２と、通信手段２０３と、入力手段２０４と、メモリ２０５と、記憶装置２０６と、がバス線２１１に接続されている。ＣＰＵ２０１は、記憶装置２０６に格納された計算プログラムを実行して、系統状態の算出や、制御信号の生成、画像データの生成などを行い、画像データを表示装置２０２に表示させる。メモリ２０５は、表示用の画像データ、系統状態の算出結果データなどを一旦格納するメモリであり、例えば、ＲＡＭ（Random Access Memory）などである。通信手段２０３は、情報通信ネットワーク１０８を通じて計測装置２１から潮流値、電圧値などの系統状態量を取得する。制約情報取得部１１は、通信手段２０３を用いて、複数の計測装置２１で計測された時系列の計測情報を受信する。
【0040】
需給計画装置２００のユーザ（管理者）は、入力手段２０４の所定のインタフェースを通じて各種閾値などのパラメータを設定又は変更し、自系統の需給計画装置２００の動作を適切に設定できる。また、ユーザは、入力手段２０４の所定のインタフェースを通じて確認したいデータの種類を選択し、表示装置２０２に表示させることができる。
【0041】
記憶装置２０６は、各種プログラムおよびデータを保持する。記憶装置２０６は、例えば、ＨＤＤ（Hard Disk Drive）又はフラッシュメモリなどである。記憶装置２０６は、例えば、後述する各種機能を実現し得るプログラムおよびデータなどを保持する。記憶装置２０６に記憶されているプログラムおよびデータは、必要に応じてＣＰＵ２０１に読み出されて実行される。なお、記憶装置２０６は、情報記憶部３０を実現し、各種のデータベース（ＤＢ）を格納する。
【0042】
計測装置２１は、系統内の計測地点における計測情報を計測し、制約情報取得部１１に送信する。計測情報は、有効電力潮流、無効電力潮流、系統電圧、系統電流、電圧位相の中の少なくとも一つの種別を含む。計測装置２１は、複数の計測地点に夫々配置された計測装置２１ａ，２１ｂ，・・・，２１ｎを含む。計測装置２１は、例えばＰＭＵ（Power Monitoring Unit）などであり、電力に関する定期的な計測を行い、時系列の計測情報を送信する。
【0043】
図１に戻り、情報記憶部３０は、需給計画において考慮すべき制約情報を格納した制約情報ＤＢ３１、および制約情報を拡張して検証シナリオを生成するためのパラメータを格納した制約拡張情報ＤＢ３２を含む。
【0044】
制約情報ＤＢ３１には、各種制約項目が反映される対象設備を示す制約対象設備、各種設備の運用量の上限値を規定する制約上限量、各種設備の運用量の下限値を規定する制約下限量、各種設備の運用量の指定値を規定する制約指定量、各種制約項目が反映される期間を示す制約対象期間、の各種データ項目を含む制約情報が格納される。
【0045】
また、制約拡張情報ＤＢ３２には、各種制約上限量、制約下限量、または制約指定量のいずれかの値を一定の範囲内で変更可能とする際の数量増減時の下限値および上限値、左記の数量増減時の下限値および上限値の範囲内で制約値を変動させる場合の刻み幅を規定する数量増減時の変化刻み幅、制約対象期間を一定の範囲内で変更可能とする際の期間増減時の開始時間および終了時間の下限値および上限値、左記の期間増減時の下限値および上限値の範囲内で制約値を変動させる場合の刻み幅を規定する期間増減時の変化刻み幅、各種制約項目が反映される期間の種別（時刻、昼間・夜間、平日・休日、曜日、月など）を示す制約対象期間種別、の各種データ項目を含む制約拡張情報が格納される。
【0046】
ここで、制約情報の一例として、図４に示す制約項目２に関連する制約情報では、制約対象設備として送電線Ａ－１およびＡ－２、制約上限量として３００［ＭＷ］、制約対象期間としては各時刻が設定されており、制約下限量および制約指定量は規定されていないことを表している。
【0047】
また、制約拡張情報の一例として、図５Ａの制約上限量について説明すると、元となる制約情報の５０％から２００％の範囲内において制約量の１０％を刻み幅として設定可能とする制約拡張情報が示されている。
【0048】
ここで、図４に示す制約項目２に示す制約情報をベース制約とする場合、前記の制約拡張情報に基づき、３００［ＭＷ］×５０［％］＝１５０［ＭＷ］を下限、３００［ＭＷ］×２００［％］＝６００［ＭＷ］を上限、３００［ＭＷ］×１０［％］＝３０［ＭＷ］を刻み幅として、当該の制約対象設備における制約上限量を１５０、１８０、…、５７０、６００［ＭＷ］のいずれかに設定された検証シナリオが作成される。
【0049】
また、制約拡張情報の別の一例として、制約対象期間について図５Ｂに示す制約拡張情報を用いて、制約期間を２０２３年４月から２０３５年３月までとする制約情報をベース制約とする場合、制約対象期間の開始時間としては２０２３年４月から２０２８年４月までの１ヶ月刻み、終了時間としては２０３２年３月から２０３５年３月までの１ヶ月刻み、制約対象期間種別としては全日、休日のみ、平日の０時から７時および１９時から２３時、のいずれかを設定された検証シナリオが作成される。
【0050】
＜検証シナリオ生成処理Ｓ３００＞
需給計画装置２００による検証シナリオ生成処理Ｓ３００について説明する。
図３は、検証シナリオ生成処理Ｓ３００の全体処理を示す。本処理の実施目的の１つは、過去の計画策定時に設定された運用制約条件の実績値または推定値、および過去の計画策定時点から本機能を使用する現時点までの間に取得された運用実績値（発電機における発電出力値、送電線などの送配電設備における潮流値、発電機や変圧器などの送配電設備における制御指令値情報など）に基づき、過去の計画策定時に設定された運用制約条件の適用状況を実績値との比較により評価することである。
需給計画装置２００による検証シナリオ生成処理Ｓ３００について説明する。
図３は、検証シナリオ生成処理Ｓ３００の全体処理を示す。本処理の実施目的の１つは、過去の計画策定時に設定された運用制約条件の実績値または推定値、および過去の計画策定時点から本機能を使用する現時点までの間に取得された運用実績値（発電機における発電出力値、送電線などの送配電設備における潮流値、発電機や変圧器などの送配電設備における制御指令値情報など）に基づき、過去の計画策定時に設定された運用制約条件の適用状況を実績値との比較により評価することである。
【0051】
また、前記の評価の情報に基づき、現時点より将来の運用制約条件を推定し、推定された将来の運用制約条件に基づき将来の運用計画を算出することも本処理の目的の１つである。
【0052】
本処理が必要となる背景として、発送分離などの事業制度変更を経て電力事業の自由化が進む中で、送配電事業者が安定な電力系統運用を実現するため、多様な事業者の運用状態を限定された公開情報に基づき推定する必要が生じることがあげられる。
【0053】
（ステップＳ３０１）
まず、制約情報取得部１１は、需給計画の策定において考慮が必要となる制約情報を情報記憶部３０に含まれる制約情報ＤＢ３１より取得する（ステップＳ３０１）。ここで、制約情報の取得方法として、使用者による分析ツール・ソフトウェアを用いた手動入力、制約情報が格納されたファイルインポート、他システムとのＤＢ連携のいずれかを採用して良い。
まず、制約情報取得部１１は、需給計画の策定において考慮が必要となる制約情報を情報記憶部３０に含まれる制約情報ＤＢ３１より取得する（ステップＳ３０１）。ここで、制約情報の取得方法として、使用者による分析ツール・ソフトウェアを用いた手動入力、制約情報が格納されたファイルインポート、他システムとのＤＢ連携のいずれかを採用して良い。
【0054】
（制約情報）
図４は、制約情報３１Ａの一例を示す図である。図４には、制約情報ＤＢ３１に格納される制約情報３１Ａの一例を示す。図４の各行で表される各制約項目に対して、制約対象設備３１１、制約上限量３１２、制約下限量３１３、制約指定量３１４、制約対象期間３１５の設定値が格納される。ここで、制約対象設備３１１と制約対象期間３１５については、いずれの制約項目についても設定が必要であり、制約上限量３１２、制約下限量３１３、制約指定量３１４については、そのいずれか１つ以上の設定が必要である。
図４は、制約情報３１Ａの一例を示す図である。図４には、制約情報ＤＢ３１に格納される制約情報３１Ａの一例を示す。図４の各行で表される各制約項目に対して、制約対象設備３１１、制約上限量３１２、制約下限量３１３、制約指定量３１４、制約対象期間３１５の設定値が格納される。ここで、制約対象設備３１１と制約対象期間３１５については、いずれの制約項目についても設定が必要であり、制約上限量３１２、制約下限量３１３、制約指定量３１４については、そのいずれか１つ以上の設定が必要である。
【0055】
制約対象設備３１１は、当該の制約が適用される対象設備が設定される。設定される設備種別として、発電機、送電線、変電所、電力負荷、連系点ノードのいずれか、あるいはそれら設備の組合せを設定して良い。
【0056】
制約上限量３１２は、当該の制約項目において、指定された設定値以上の値を保つ必要があることを意味する。制約上限量として、有効電力［ｋＷ］または［ＭＷ］、無効電力［ｋｖａｒ］または［Ｍｖａｒ］、有効電力量［ｋＷｈ］または［ＭＷｈ］、系統電圧［Ｖ］または［ｋＶ］、力率［％］、電流位相あるいは電圧位相［ｄｅｇ］または［ｒａｄ］を設定して良い。
【0057】
制約下限量３１３は、当該の制約項目において、指定された設定値以下の値を保つ必要があることを意味する。制約下限量として、有効電力［ｋＷ］または［ＭＷ］、無効電力［ｋｖａｒ］または［Ｍｖａｒ］、有効電力量［ｋＷｈ］または［ＭＷｈ］、系統電圧［Ｖ］または［ｋＶ］、力率［％］、電流位相あるいは電圧位相［ｄｅｇ］または［ｒａｄ］を設定して良い。
【0058】
制約指定量３１４は、当該の制約項目において、指定された設定値と同値を保つ必要があることを意味する。制約指定量として、有効電力［ｋＷ］または［ＭＷ］、無効電力［ｋｖａｒ］または［Ｍｖａｒ］、有効電力量［ｋＷｈ］または［ＭＷｈ］、系統電圧［Ｖ］または［ｋＶ］、力率［％］、電流位相あるいは電圧位相［ｄｅｇ］または［ｒａｄ］を設定して良い。
【0059】
制約対象期間３１５は、制約の遵守可否を判定するための期間が設定される。制約対象期間として、ある複数の時間インターバルを含む期間を設定して良い。また、個々の時間インターバルに対して個別の制約値が設定される制約項目に対して、制約対象期間として個々の時間インターバルを設定して良い。また、制約対象期間として、ある一定の条件を満たす期間を設定して良い。左記の条件として、年度、季節、曜日、平休日種別、時刻のいずれかを用いて、制約対象期間の条件を設定して良い。
【0060】
（ステップＳ３０２）
次に、制約情報取得部１１は、検証シナリオ生成において制約情報の変更に用いられるパラメータが格納された制約拡張情報を情報記憶部３０に含まれる制約拡張情報ＤＢより取得する（ステップＳ３０２）。
次に、制約情報取得部１１は、検証シナリオ生成において制約情報の変更に用いられるパラメータが格納された制約拡張情報を情報記憶部３０に含まれる制約拡張情報ＤＢより取得する（ステップＳ３０２）。
【0061】
ここで、制約拡張情報の取得方法として、使用者による分析ツール・ソフトウェアを用いた手動入力、制約拡張情報が格納されたファイルインポート、他システムとのＤＢ連携のいずれかを採用して良い。
【0062】
（制約拡張情報）
図５Ａは、制約拡張情報３２Ａの一例（制約対象設備の増減、制約量の増減を対象とする場合）を示す図である。図５Ｂは、制約拡張情報３２Ｂの一例（制約対象期間の増減を対象とする場合）を示す図である。図５Ａ、図５Ｂは、制約拡張情報ＤＢ３２に格納される制約拡張情報の一例を示す。
図５Ａは、制約拡張情報３２Ａの一例（制約対象設備の増減、制約量の増減を対象とする場合）を示す図である。図５Ｂは、制約拡張情報３２Ｂの一例（制約対象期間の増減を対象とする場合）を示す図である。図５Ａ、図５Ｂは、制約拡張情報ＤＢ３２に格納される制約拡張情報の一例を示す。
【0063】
図５Ａは、制約対象設備３２１の増減、あるいは制約量（制約上限量３２２、制約下限量３２３）の増減を対象とする場合の制約拡張情報３２Ａの一例を示している。制約対象設備３２１の増減について、検証シナリオとして設定可能とする各制約項目の対象とする系統設備の組合せ条件を指定して良い。図５Ａでは、ベースとする制約条件における当該制約項目の対象となる系統設備の設備数に対する割合を用いて、制約対象設備として設定可能とする設備数の上下限を指定する例を示す。
【0064】
制約量の増減について、図４に示す制約上限量３１２、制約下限量３１３、制約指定量３１４のいずれかに対して、検証シナリオとして設定可能とする制約量の上下限および変動幅を指定して良い。図５Ａでは、ベースとする制約条件における当該制約項目の制約量に対する割合を用いて、制約量として設定可能とする範囲の上下限を指定する例を示す。
【0065】
具体的には、制約対象設備３２１は、ベース制約の設備数の５０％から対象系統の全連系設備であるのに対し、制約上限量３２２は、ベース制約の制約量の５０％からベース制約の制約量の２００％の範囲であり、制約下限量３２３は、ベース制約の制約量の８０％からベース制約の制約量の１２０％の範囲であることが示されている。
【0066】
図５Ｂには、制約対象期間３２４の増減を対象とする場合の制約拡張情報３２Ｂの一例を示している。制約対象期間３２４の増減について、検証シナリオとして設定可能とする各制約項目の適用対象期間を指定して良い。図５Ｂでは、ベースとする制約条件における当該制約項目の適用対象期間に対して、当該期間の開始時間、終了時間の上下限値、変化幅、および期間種別を用いて、検証シナリオとして設定可能とする制約対象期間を指定する例を示す。ここで、期間種別として、年度、季節、曜日、平休日種別、時刻を設定して良い。
【0067】
（ステップＳ３０３）
次に、過去実績取得部１２は、各種発電機における発電出力、または送電線、変電所を含む各種系統設備における電力潮流のいずれか１つ以上を含む過去の電力系統運用実績から構成される過去実績情報を情報記憶部３０に含まれる過去実績ＤＢ３３より取得する（ステップＳ３０３）。ここで、過去実績情報の取得方法として、使用者による分析ツール・ソフトウェアを用いた手動入力、制約拡張情報が格納されたファイルインポート、他システムとのＤＢ連携のいずれかを採用して良い。
次に、過去実績取得部１２は、各種発電機における発電出力、または送電線、変電所を含む各種系統設備における電力潮流のいずれか１つ以上を含む過去の電力系統運用実績から構成される過去実績情報を情報記憶部３０に含まれる過去実績ＤＢ３３より取得する（ステップＳ３０３）。ここで、過去実績情報の取得方法として、使用者による分析ツール・ソフトウェアを用いた手動入力、制約拡張情報が格納されたファイルインポート、他システムとのＤＢ連携のいずれかを採用して良い。
【0068】
図６は、過去実績情報３３Ａの一例を示す図である。図６には、過去実績ＤＢ３３に格納される過去実績情報３３Ａの一例を図６に示す。過去実績情報３３Ａは、指定された過去情報のデータ期間、過去情報の取得日時、過去情報の対象設備、過去情報の対象エリアのいずれか1つ以上の条件を満たす過去情報を含む。図６には、需給計画に含まれる各系統設備の過去実績情報３３Ａの一例として、複数の発電機（発電機ａ－１，ａ―２，ｂ－１，ｂ－２）に対する各時間インターバルの発電量を示している。過去実績情報に含まれる系統設備の種別として、発電機以外の各種設備、すなわち送電線、変電所、電力負荷、連系点ノードにおける電力潮流量を過去計画情報として取得して良い。
【0069】
（ステップＳ３０４）
次に、過去計画取得部１３は、各種発電機における発電出力、または送電線、変電所を含む各種系統設備における電力潮流のいずれか１つ以上を含む過去時点において策定した電力系統運用計画から構成される過去計画情報を情報記憶部３０に含まれる過去計画ＤＢ３４より取得する（Ｓ３０４）。ここで、過去計画情報の取得方法として、使用者による分析ツール・ソフトウェアを用いた手動入力、制約拡張情報が格納されたファイルインポート、他システムとのＤＢ連携のいずれかを採用して良い。
次に、過去計画取得部１３は、各種発電機における発電出力、または送電線、変電所を含む各種系統設備における電力潮流のいずれか１つ以上を含む過去時点において策定した電力系統運用計画から構成される過去計画情報を情報記憶部３０に含まれる過去計画ＤＢ３４より取得する（Ｓ３０４）。ここで、過去計画情報の取得方法として、使用者による分析ツール・ソフトウェアを用いた手動入力、制約拡張情報が格納されたファイルインポート、他システムとのＤＢ連携のいずれかを採用して良い。
【0070】
図７Ａは、過去計画情報３４Ａの一例（複数の発電機に対する各時間インターバルの発電量）を示す図である。図７Ｂは、過去計画情報３４Ｂの一例（各時間インターバルにおいて各制約項目の制約条件が有効あるいは無効のいずれであるか）を示す図である。図７Ａ、図７Ｂは、過去計画ＤＢ３４に格納される過去計画情報の一例を示す。
【0071】
図７Ａには、需給計画に含まれる各系統設備の過去計画情報３４Ａの一例として、複数の発電機（発電機ａ－１，ａ―２，ｂ－１，ｂ－２）に対する各時間インターバルの発電量を示している。過去計画情報に含まれる系統設備の種別として、発電機以外の各種設備、すなわち送電線、変電所、電力負荷、連系点ノードにおける電力潮流量を過去計画情報として取得して良い。
【0072】
図７Ｂには、需給計画策定において加味された過去計画情報３４Ｂの一例として、各時間インターバルにおいて各制約項目の制約条件（発電機グループＡのｋＷｈ制約、発電機ｂ－１の計画停止制約、発電機ｂ－２の計画停止制約、発電機グルーブＢのｋＷｈ制約）が有効あるいは無効のいずれであるかを示している。図において、チェックが付記されているものが有効を意味する。
【0073】
ここで、過去制約情報として、制約条件の有効無効状態に加え、各制約項目において設定された制約量、すなわち制約上限量、制約下限量、制約指定量を取得して良い。また、複数の系統設備に関わる制約項目については、各時間インターバルにおいて、いずれの系統設備が各制約項目の対象となっているかを過去制約情報として取得して良い。
【0074】
図７Ｂに示す発電機グループＡのｋＷｈ制約の過去制約情報を用いて、図４に示す制約項目１として反映する場合、発電機グループＡが図４の制約対象設備に、制約条件が有効となっている時間コマの情報が図４の制約対象期間に、それぞれ設定される。また、図７Ｂに記載された情報に加えて、発電機グループＡのｋＷｈ制約の上限値が過去制約情報に含まれる場合は、左記情報に基づき図４の制約上限量３１２が設定される。
【0075】
（ステップＳ３０５）
次に、過去実績評価部１４は、過去実績情報に基づく評価指標を算出する（ステップＳ３０５）。ここで、評価指標は、ある一定期間の過去実績情報に基づき、平均値、最大値、最小値、度数分布（ヒストグラム）を例とした統計的指標値として算出される。
次に、過去実績評価部１４は、過去実績情報に基づく評価指標を算出する（ステップＳ３０５）。ここで、評価指標は、ある一定期間の過去実績情報に基づき、平均値、最大値、最小値、度数分布（ヒストグラム）を例とした統計的指標値として算出される。
【0076】
図８は、過去実績情報３３Ａの評価指標情報３５Ａの一例を示す図である。図８は、過去実績評価部１４において算出される評価指標の一例を示す。図８に示す過去実績の評価指標情報３５Ａの例では、評価指標の種別として、各時刻の発電機の出力を表す発電出力３５１、発電機が起動してから次に停止するまでの時間間隔を表す起動時間長３５２、発電機が停止してから次に起動するまでの時間間隔を表す停止時間長３５３、ある一定期間において発電機が起動または停止した回数を表す起動停止回数３５４、を算出した例を示す。
【0077】
これらの評価指標種別それぞれについて、平均値、最大値、最小値、度数分布（ヒストグラム）を例とした統計的指標値が過去実績情報３３Ａより算出される。ここで、各統計的指標値を算出する期間については、事前設定された固定値、本装置を使用するユーザによる手動設定などの手段により設定される。
【0078】
過去実績評価部１４は、過去実績情報から算出された平均値、上限値、下限値、度数分布のいずれか1つ以上を含む統計的な評価指標値を用いて過去実績情報と過去制約情報を比較評価する。
【0079】
（ステップＳ３０６）
次に、ステップＳ３０５において過去実績評価部１４により算出された過去実績情報３３Ａの評価指標に基づき、制約条件推定部１５は、過去の電力系統運用において適用された制約条件を推定する（ステップＳ３０６）。
次に、ステップＳ３０５において過去実績評価部１４により算出された過去実績情報３３Ａの評価指標に基づき、制約条件推定部１５は、過去の電力系統運用において適用された制約条件を推定する（ステップＳ３０６）。
【0080】
ここで、過去実績情報３３Ａの評価指標および過去計画情報３４Ａ，３４Ｂに含まれる過去制約情報３１ＡＰ（制約情報３１Ａのうち過去制約情報）（図９参照）を比較評価することにより、過去の電力系統運用において適用された制約条件を推定して良い。
【0081】
このとき、過去制約情報３１ＡＰ（図９参照）において規定される各制約項目について、過去実績情報３３Ａの評価指標との適合度を算出する。図４の例に示す制約項目１は、制約の対象設備が発電機ａ－１、発電機ａ－２、発電機ｂ－１の発電出力であり、制約の内容として制約下限量３１３が設定され、制約対象期間３１５が複数の時間刻みを含む一定期間（図４の例では１．５年）である。このような制約項目の過去実績情報３３Ａとの適合度を算出する場合、制約対象期間３１５と同等の期間における発電量［Ｗｈ］の実績値ＥＮＲを過去実績情報３３Ａより抽出し、過去制約情報３１ＡＰにおける電力量下限ＥＮｌｏｗと比較する。このとき、本制約項目における過去実績情報３３Ａとの適合度ｒ１を（１）式より算出して良い。ここで、αｌｏｗは事前設定される固定パラメータとする。
【数1】
【0082】
同様に、制約項目として電力量上限ＥＮｈｉｇｈが過去制約情報において規定される場合、本制約項目における過去実績情報との適合度ｒ２を（２）式より算出して良い。ここで、αｈｉｇｈは事前設定される固定パラメータとする。
【数2】
【0083】
また、図４の例に示す制約項目２は、制約の対象設備が送電線Ａ－１、送電線Ａ－２の電力潮流であり、制約の内容として電力潮流上限が設定され、制約対象期間が各時刻である。このような制約項目の過去実績情報との適合度を算出する場合、ある期間における電力潮流の平均値ＰＮＲａｖｇを過去実績情報の評価指標より抽出し、過去制約情報における潮流上限ＰＮｈｉｇｈと比較する。このとき、本制約項目における過去実績情報との適合度ｒ３を（３）式より算出して良い。ここで、βｈｉｇｈは事前設定される固定パラメータとする。
【数3】
【0084】
また、本制約項目に関する適合度算出の別手法として、ある期間における電力潮流の度数分布ｆ(ＰＮＲ)を過去実績情報の評価指標より抽出し、過去制約情報における潮流上限ＰＮｈｉｇｈと比較する。このとき、本制約項目における過去実績情報との適合度ｒ４を（４）式より算出して良い。ここで、ｆ％(ＰＮＲ)は、度数分布ｆ(ＰＮＲ)より算出される相対度数、βｈｉｇｈは事前設定される固定パラメータとする。
【数4】
【0085】
また、図４の例に示す制約項目３は、制約の対象設備が発電機ａ－１の発電出力であり、制約の内容として発電出力指定値が設定され、制約対象期間が各時刻である。このような制約項目の過去実績情報との適合度を算出する場合、ある期間における発電出力の過去実績値ＰＲ(ｔ)を過去実績情報より抽出し、過去制約情報における出力指定値Ｐｆｉｘと比較する。このとき、本制約項目における過去実績情報との適合度ｒ５を（５）式より算出して良い。ここで、γは事前設定される固定パラメータ、g(ＰＲ(ｔ))は時刻ｔにおける適合度、ｎは適合度を評価する過去実績情報の時間幅あるいは時間コマ数、とする。
【数5】
【0086】
以上の事例に示す通り、過去制約情報において制約として規定される制約項目のそれぞれについて適合度ｒｎを算出する。
【0087】
次に、前記において算出された適合度ｒｎを用いて推定制約情報３１Ｅを算出する。
図９は、過去時点における推定制約情報３１Ｅの一例を示す図である。図９には、過去実績情報３３Ａの評価指標および過去計画情報３４Ａ，３４Ｂに含まれる過去制約情報３１ＡＰにより推定される推定制約情報３１Ｅを示す。
図９は、過去時点における推定制約情報３１Ｅの一例を示す図である。図９には、過去実績情報３３Ａの評価指標および過去計画情報３４Ａ，３４Ｂに含まれる過去制約情報３１ＡＰにより推定される推定制約情報３１Ｅを示す。
【0088】
図９の例では、過去制約情報３１ＡＰとして制約項目１、制約項目２、制約項目３が規定され、適合度の判定閾値として０．５が指定された例を示す。このとき、制約項目２については、過去実績情報より算出された適合度が０．３であり、判定閾値を下回っている。判定結果より、過去制約情報に格納された制約上限量３００［ＭＷ］に対して情報更新が必要と判定され、図８に示す過去実績情報３３Ａの評価指標、例えば過去の一定期間における対象設備の運用実績の平均値である４００［ＭＷ］を新たな制約上限量として設定することを表す。
【0089】
ここで、制約項目１および制約項目３については、適合度が判定閾値を上回っているため、推定制約情報は過去制約情報と同等の内容になる。また、制約項目２については、適合度が判定閾値を下回っているため、過去の計画策定時点において想定した運用制約とは異なる運用制約が適用された影響が過去実績に表れていると推測される。
【0090】
このような場合は、推定制約情報３１Ｅとして、過去制約情報３１ＡＰと異なる値が設定される。ここで、制約項目２は制約上限に関する制約項目であるため、過去実績情報の評価指標における一定期間における上限値を推定制約情報として設定して良い。同様に、制約下限に関する制約項目について過去制約情報と異なる推定制約情報を設定する場合は、過去実績情報の評価指標における一定期間における下限値を推定制約情報として設定して良い。また、制約指定量に関する制約項目について過去制約情報と異なる推定制約情報を設定する場合は、過去実績情報の評価指標における一定期間における平均値または中央値を推定制約情報として設定して良い。
【0091】
（ステップＳ３０７）
次に、制約条件推定部１５は、将来の電力系統運用において適用される将来推定制約情報を算出する（ステップＳ３０７）。ここで、過去計画策定時点、および過去計画策定時点より未来の１つ以上の時点における推定制約情報を用いて、将来の推定制約情報を算出して良い。
次に、制約条件推定部１５は、将来の電力系統運用において適用される将来推定制約情報を算出する（ステップＳ３０７）。ここで、過去計画策定時点、および過去計画策定時点より未来の１つ以上の時点における推定制約情報を用いて、将来の推定制約情報を算出して良い。
【0092】
図１０Ａは、将来期間における推定制約情報３１Ｅの一例を示す図である。図１０Ｂは、図１０Ａに示す過去推定制約情報３１ＥＰ（太い破線）および将来推定制約情報３１ＥＦ（太い実線）を時系列グラフ上に表示した一例を示す図である。
【0093】
図１０Ａには、制約条件推定部１５において算出される将来推定制約情報３１ＥＦの一例を示す。図１０Ａの例では、制約項目１、制約項目２、制約項目３、の各制約項目について、過去期間１、過去期間２、過去期間３の各期間において過去実績情報を用いて過去推定制約情報３１ＥＰを算出した結果を示す。さらに、それらの過去期間における過去推定制約情報３１ＥＰを用いて、将来期間４、将来期間５における将来推定制約情報３１ＥＦを算出した例を示す。
【0094】
ここで、過去期間における推定制約情報３１Ｅを用いて、時系列方向の制約条件値の変更を自己回帰（ＡＲ）、移動平均（ＭＡ）、自己回帰移動平均（ＡＲＭＡ）などの線形時系列解析手法を用いてモデル化し、将来期間における将来推定制約情報を算出して良い。また、電力需要、電力市場取引、送電容量など電力系統運用に影響を及ぼす外部要因との相関関係を多変量線形回帰、非線形回帰などの手法を用いてモデル化し、上記の外部要因の将来予測値を上記のモデルの入力とすることで、将来期間における将来推定制約情報を算出して良い。
【0095】
（ステップＳ３０８）
次に、制約変更条件生成部１６は、制約拡張情報３２Ａ，３２Ｂ、過去実績情報３３Ａ、過去時点における過去推定制約情報３１ＥＰ、および将来推定制約情報３１ＥＦを用いて、制約変更情報３１Ｍを算出する（ステップＳ３０８）。
次に、制約変更条件生成部１６は、制約拡張情報３２Ａ，３２Ｂ、過去実績情報３３Ａ、過去時点における過去推定制約情報３１ＥＰ、および将来推定制約情報３１ＥＦを用いて、制約変更情報３１Ｍを算出する（ステップＳ３０８）。
【0096】
図１１Ａは、将来時点における制約変更情報３１Ｍを示す図であり、制約対象設備の増減、あるいは制約量の増減を対象とする場合の制約拡張情報３２ＡＭの一例を示す図である。図１１Ｂは、将来時点における制約変更情報３１Ｍを示す図であり、制約対象期間の増減を対象とする場合の制約拡張情報３２ＢＭの一例を示す図である。図１１Ａは、図５Ａの制約拡張情報３２Ａに対し、制約変更条件生成部１６において算出される制約変更情報３１Ｍの例である。また、図１１Ｂは、図５Ｂの制約拡張情報３２Ｂに対し、制約変更条件生成部１６において算出される制約変更情報３１Ｍの例である。なお、制約変更情報３１Ｍには、過去制約情報からの変更点が明記されている。
【0097】
図１１Ａの制約対象設備の増減について、検証シナリオとして設定可能とする各制約項目の対象とする系統設備の組合せ条件を指定して良い。図１１Ａでは、ベースとする制約条件における当該制約項目の対象となる系統設備の設備数に対する割合を用いて、制約対象設備として設定可能とする設備数の上下限を指定する例を示す。制約量の増減について、制約上限量、制約下限量、制約指定量のいずれかに対して、検証シナリオとして設定可能とする制約量の上下限および変動幅を指定して良い。
【0098】
図１１Ｂの制約対象期間の増減について、検証シナリオとして設定可能とする各制約項目の適用対象期間を指定して良い。図１１Ｂでは、ベースとする制約条件における当該制約項目の適用対象期間に対して、当該期間の開始時間、終了時間の上下限値、変化幅、および期間種別を用いて、検証シナリオとして設定可能とする制約対象期間を指定する例を示す。ここで、期間種別として、年度、季節、曜日、平休日種別、時刻を設定して良い。
【0099】
過去時点における計画策定に用いられた制約拡張情報に対し、過去実績情報より統計的な指標を算出し、その算出結果を過去制約拡張情報に反映して良い。このとき、ｋＷｈ制約に関連する制約項目については、過去計画情報を用いて一定期間の発電エネルギー積算値を算出し、複数の期間あるいは条件に対して左記のエネルギー積算値をそれぞれ算出し、取得した複数期間のエネルギー積算値の集合に対する統計情報、すなわち平均値、中央値、分散、標準偏差、事前規定した閾値に基づく分位点におけるサンプル値、のいずれかをｋＷｈ制約における制約上限または制約下限として、過去制約情報および過去制約拡張情報を設定して良い。
【0100】
同様に、計画停止に関連する制約については過去計画情報から事前規定した設備の停止期間を抽出し、上記と同様の積算値を用いた統計情報に基づき、作業停止制約における制約上限または制約下限として、過去制約情報および過去制約拡張情報を設定して良い。上記以外の制約項目について、各種制約に関連する過去計画情報の特徴量を算出し、統計分析を適用することにより、過去制約情報および過去制約拡張情報を設定して良い。
【0101】
また、過去時点における推定制約情報、および将来推定制約情報を用いて、過去制約拡張情報を設定して良い。一例として、将来時点における推定制約情報に含まれる制約上限値が過去時点において想定された過去制約情報の設定値よりも大きな値を取る場合、制約上限の緩和が将来的に想定される。このようなケースにおいて、図１１Ａの例に示す「数量増減時の上限」の「制約上限量」の設定値２００％をさらに増加させることにより、将来時点において制約上限が想定より緩和される条件を含む検証シナリオを設定することができる。
【0102】
制約変更条件生成部１６は、過去時点および将来時点の推定制約情報に基づき、各制約項目の設定値の増減における上下限および変化幅、制約の考慮対象となる期間の増減における開始時間、終了時間、および期間長の上下限および変化幅、制約の考慮対象となる期間の年度、季節、曜日、平休日種別、時刻のいずれか1つ以上を含む制約変更情報３１Ｍを生成する。
【0103】
（ステップＳ３０９）
次に、シナリオ生成部１７は、推定制約情報と制約変更情報に基づき、検証シナリオ情報３７を生成する（ステップＳ３０９）。ここで、推定制約情報３１Ｅに含まれるそれぞれの制約項目について、ステップＳ３０８において算出した制約変更情報３１Ｍに含まれるそれぞれの制約項目の変動範囲、すなわち変動上限、変動下限、変動幅に基づき生じ得る制約設定値の組合せを算出し、選択した全ての制約項目について左記制約設定値の組合せのいずれかを選択し、左記設定値の総ての組合せ条件を規定することで検証シナリオを生成して良い。
次に、シナリオ生成部１７は、推定制約情報と制約変更情報に基づき、検証シナリオ情報３７を生成する（ステップＳ３０９）。ここで、推定制約情報３１Ｅに含まれるそれぞれの制約項目について、ステップＳ３０８において算出した制約変更情報３１Ｍに含まれるそれぞれの制約項目の変動範囲、すなわち変動上限、変動下限、変動幅に基づき生じ得る制約設定値の組合せを算出し、選択した全ての制約項目について左記制約設定値の組合せのいずれかを選択し、左記設定値の総ての組合せ条件を規定することで検証シナリオを生成して良い。
【0104】
図１２は、検証シナリオ情報３７の一例を示す図である。図１２は、シナリオ生成部１７により生成される検証シナリオ情報３７の一例を示す。図１２の例では、各行に示される個々の検証シナリオ（シナリオａ―０、シナリオａ―１、シナリオｂ－０）に対して、各制約項目３７１（制約項目１），３７２（制約項目２）の適用有無、各制約項目の制約対象設備、各制約項目の制約量、各制約項目の制約対象期間、およびシナリオ生成条件を検証シナリオ情報として格納している。
【0105】
図１２の例では、図１０に示す制約項目であるｋＷｈ制約および計画停止制約を用いて、制約適用有無、制約対象設備、制約量、制約対象期間を図１１に示す制約拡張情報に基づく範囲内で設定し、各制約項目の制約条件を組合せることで検証シナリオを生成する例を示す。ここで、シナリオ生成条件３７３として、制約情報３１Ａおよび制約拡張情報３２Ａ、３２Ｂに基づき、ベースとなる制約条件に対してどのような操作を加えて検証シナリオにおける制約条件を設定したか、の情報を格納して良い。図１２の例では、ベースの制約条件における制約量の増減、および複数の制約項目を図１１に示す制約拡張情報に基づく範囲内で設定した条件のいずれを組み合わせて当該検証シナリオにおける制約条件を設定したか、の情報をシナリオ生成条件として格納している。
【0106】
（ステップＳ３１０）
計画策定部１８は、ステップＳ３０９において生成された検証シナリオを用いて、検証シナリオに基づく発電計画を算出する（ステップＳ３１０）。このとき、検証シナリオにより規定された制約条件を反映した最適化問題を以下の例に示す通り定式化し、最適化問題の計算結果に基づき、発電機の起動停止状態および発電出力を含む発電計画を算出して良い。ここで、下記の制約条件を満たした上で、目的関数を最小化する問題として最適化問題が定義される。
計画策定部１８は、ステップＳ３０９において生成された検証シナリオを用いて、検証シナリオに基づく発電計画を算出する（ステップＳ３１０）。このとき、検証シナリオにより規定された制約条件を反映した最適化問題を以下の例に示す通り定式化し、最適化問題の計算結果に基づき、発電機の起動停止状態および発電出力を含む発電計画を算出して良い。ここで、下記の制約条件を満たした上で、目的関数を最小化する問題として最適化問題が定義される。
【0107】
ここで目的関数の一例は、以下の（６）式で表される。ここで、数式中の各記号は、ｔ：時刻を表すインデックス、ｉ：発電機を表すインデックス、Ｃｉ(ｐ)：発電機ｉの発電出力ｐに対する運用コスト、ＳＵＣｉ：発電機ｉの起動コスト、Ｄｉ(ｒ)：発電機ｉの調整力拠出ｒに対する運用コスト、ＰＣ：制約違反に伴うペナルティコスト、をそれぞれ表す。
【数6】
【0108】
制約条件は以下のいずれか１つ以上を含む。（ａ）需給バランス制約：割当対象の電力需要が発電出力の合計と一致する、（ｂ）発電機出力上下限制約：各発電機が出力下限から出力上限の範囲内の出力値を取る、（ｃ）発電機出力変化速度制約：異なる時刻間の各発電機の出力値の変化幅が出力変化速度上限以内となる、（ｄ）最小起動時間・最小停止時間制約：発電機の連続した起動時間または停止時間が規定された最小時間以上となるよう運用する、（ｅ）指定起動・停止期間：運用者が指定した期間において発電機を強制的に停止または運転する、（ｆ）運転予備力・必要調整力制約：実運用時と需給計画との電力需要の誤差を補正する発電機余力の必要量を確保する。
【0109】
ここで、検証シナリオにおいて規定されない制約条件についても、最適化問題における制約条件として設定して良い。このとき、過去計画情報より一定の期間を抽出し、当該期間において設定されている制約情報の中で、左記検証シナリオにおいて規定されている以外の制約条件を抽出し、その設定値を用いて最適化問題の制約条件を設定して良い。
【0110】
本実施形態の需給計画装置２００および需給計画方法は、次の特徴を有する。
（１）各種の運用制約条件を満たした上で運用計画を策定する需給計画装置２００であって、過去時点における運用実績の情報を含む過去実績情報３３Ａおよび過去時点において設定された制約情報の情報を含む過去制約情報３１ＡＰに基づき、過去実績情報３３Ａと過去制約情報３１ＡＰを比較評価する過去実績評価部１４と、過去実績評価部１４の過去実績評価結果に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成部１７と、を備える。これにより、実際の系統運用状態を反映した運用制約条件を加味した網羅的な検証シナリオの生成と処理時間の低減を両立することができる。
（１）各種の運用制約条件を満たした上で運用計画を策定する需給計画装置２００であって、過去時点における運用実績の情報を含む過去実績情報３３Ａおよび過去時点において設定された制約情報の情報を含む過去制約情報３１ＡＰに基づき、過去実績情報３３Ａと過去制約情報３１ＡＰを比較評価する過去実績評価部１４と、過去実績評価部１４の過去実績評価結果に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成部１７と、を備える。これにより、実際の系統運用状態を反映した運用制約条件を加味した網羅的な検証シナリオの生成と処理時間の低減を両立することができる。
【0111】
すなわち、供給力不足による需給逼迫などの事前対策検証のため、各種の想定外変動を加味した複数シナリオを考慮した計画策定が必要となる。このとき、系統運用状況を考慮した網羅的なシナリオ検証に基づく事前対策とともに、実際の運用状態に伴い今後の系統運用において想定されるシナリオを随時更新する必要がある。ここで、機械的に検証シナリオを設定すると、目的に対してシナリオが過剰になることによる計画策定の計算時間増大、またはシナリオ不足により重要シナリオの検証漏れが生じる恐れがある。このため、実際の系統運用状態を反映した運用制約条件を加味した発電機起動停止計画・出力配分を策定することにより、網羅的な評価シナリオの生成と処理時間の低減を両立することができる。
【0112】
（２）前記（１）において、過去実績情報３３Ａは、指定された過去情報のデータ期間、過去情報の取得日時、過去情報の対象設備、過去情報の対象エリアのいずれか1つ以上の条件を満たす過去情報を含む。
【0113】
（３）前記（１）において、過去実績評価部１４は、指定期間におけるｋＷｈ制約、指定期間における計画停止制約のいずれか１つ以上の制約項目を含む制約条件について過去実績情報３３Ａと過去制約情報３１ＡＰを比較評価する。
【0114】
（４）前記（１）において、過去実績評価部１４は、過去実績情報３３Ａから算出された平均値、上限値、下限値、度数分布のいずれか1つ以上を含む統計的な評価指標値を用いて過去実績情報３３Ａと過去制約情報３１ＡＰを比較評価する。
【0115】
（５）前記（１）において、需給計画装置２００は、さらに、過去実績評価部１４において算出した過去実績情報３３Ａの評価値と過去制約情報３１ＡＰに基づき、過去運用実績を反映した過去時点における運用制約を含む推定制約情報３１Ｅを算出する制約条件推定部１５を備え、シナリオ生成部１７は、制約条件推定部１５で算出された推定制約情報３１Ｅに基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成する。
【0116】
（６）前記（５）において、制約条件推定部１５は、算出された過去時点における推定制約情報３１Ｅ（（過去推定制約情報３１ＥＰ）に基づき、将来時点の推定制約情報（将来推定制約情報３１ＥＦ）を算出する。
【0117】
（７）前記（５）において、制約条件推定部１５は、制約項目の過去実績情報との適合度を用いて、過去時点における運用制約を含む推定制約情報３１Ｅを算出する。
【0118】
（８）前記（６）において、需給計画装置２００は、過去時点および将来時点の推定制約情報３１Ｅに基づき、各制約項目の設定値の増減における上下限および変化幅、制約の考慮対象となる期間の増減における開始時間、終了時間、および期間長の上下限および変化幅、制約の考慮対象となる期間の年度、季節、曜日、平休日種別、時刻のいずれか1つ以上を含む制約変更情報３１Ｍを生成する制約変更条件生成部１６をさらに備え、シナリオ生成部１７は、推定制約情報３１Ｅ及び制約変更情報３１Ｍに基づき、推定制約情報３１Ｅ及び制約変更情報３１Ｍを加味した需給計画を策定する。
【0119】
（９）前記（６）において、需給計画装置２００は、検証シナリオに含まれる制約情報に基づき、制約情報を加味した需給計画を策定する計画策定部１８をさらに備える。
【0120】
（１０）各種の運用制約条件を満たした上で運用計画を策定する需給計画方法であって、
過去時点における運用実績の情報を含む過去実績情報３３Ａおよび過去時点において設定された制約情報の情報を含む過去制約情報３１ＡＰに基づき、過去実績情報３３Ａと過去制約情報３１ＡＰを比較評価する過去実績評価ステップ(ステップＳ３０６参照)と、過去実績評価ステップにおいて算出した過去実績情報３３Ａの評価値と過去制約情報３１ＡＰに基づき、過去運用実績を反映した過去時点における運用制約を含む推定制約情報３１Ｅを算出する制約条件推定ステップ（ステップＳ３０７参照）と、制約条件推定ステップで算出された推定制約情報に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成ステップ（ステップＳ３０９参照）と、検証シナリオに含まれる制約情報に基づき、前記制約情報を加味した需給計画を策定する計画策定ステップ（ステップＳ３１０参照）と、を有する。これにより、実際の系統運用状態を反映した運用制約条件を加味した網羅的な検証シナリオの生成と処理時間の低減を両立することができる。
過去時点における運用実績の情報を含む過去実績情報３３Ａおよび過去時点において設定された制約情報の情報を含む過去制約情報３１ＡＰに基づき、過去実績情報３３Ａと過去制約情報３１ＡＰを比較評価する過去実績評価ステップ(ステップＳ３０６参照)と、過去実績評価ステップにおいて算出した過去実績情報３３Ａの評価値と過去制約情報３１ＡＰに基づき、過去運用実績を反映した過去時点における運用制約を含む推定制約情報３１Ｅを算出する制約条件推定ステップ（ステップＳ３０７参照）と、制約条件推定ステップで算出された推定制約情報に基づき、各種の運用制約条件の組み合わせとして規定される検証シナリオを生成するシナリオ生成ステップ（ステップＳ３０９参照）と、検証シナリオに含まれる制約情報に基づき、前記制約情報を加味した需給計画を策定する計画策定ステップ（ステップＳ３１０参照）と、を有する。これにより、実際の系統運用状態を反映した運用制約条件を加味した網羅的な検証シナリオの生成と処理時間の低減を両立することができる。
【0121】
（１１）前記（１０）において、需給計画方法は、さらに、推定制約情報に基づき、各制約項目の設定値の増減における上下限および変化幅、制約の考慮対象となる期間の増減における開始時間、終了時間、および期間長の上下限および変化幅、制約の考慮対象となる期間の年度、季節、曜日、平休日種別、時刻のいずれか1つ以上を含む制約変更情報を生成する制約変更条件生成ステップ（ステップＳ３０８参照）を有し、シナリオ生成ステップにおいて、推定制約情報３１Ｅ及び制約変更情報３１Ｍに基づき、推定制約情報３１Ｅ及び制約変更情報３１Ｍを加味した需給計画を策定する。
【符号の説明】
【0122】
１０ 需給計画決定部
１１ 制約情報取得部
１２ 過去実績取得部
１３ 過去計画取得部
１４ 過去実績評価部
１５ 制約条件推定部
１６ 制約変更条件生成部
１７ シナリオ生成部
１８ 計画策定部
２１ 計測装置
３０ 情報記憶部
３１ 制約情報ＤＢ
３１Ａ 制約情報
３１ＡＰ 過去制約情報
３１Ｅ 推定制約情報
３１ＥＦ 将来推定制約情報
３１ＥＰ 過去推定制約情報
３１Ｍ 制約変更情報
３２ 制約拡張情報ＤＢ
３２Ａ，３２Ｂ 制約拡張情報
３２ＡＭ、３２ＢＭ 制約拡張情報
３３ 過去実績ＤＢ
３３Ａ 過去実績情報
３４ 過去計画ＤＢ
３４Ａ，３４Ｂ 過去計画情報
３５Ａ 評価指標情報
３７ 検証シナリオ情報
１０１ 発電機
１０２ 変電所
１０３ 調相機器
１０４ 電力負荷
１０５ 外部電力系統
１０８ 情報通信ネットワーク
２００ 需給計画装置
１１ 制約情報取得部
１２ 過去実績取得部
１３ 過去計画取得部
１４ 過去実績評価部
１５ 制約条件推定部
１６ 制約変更条件生成部
１７ シナリオ生成部
１８ 計画策定部
２１ 計測装置
３０ 情報記憶部
３１ 制約情報ＤＢ
３１Ａ 制約情報
３１ＡＰ 過去制約情報
３１Ｅ 推定制約情報
３１ＥＦ 将来推定制約情報
３１ＥＰ 過去推定制約情報
３１Ｍ 制約変更情報
３２ 制約拡張情報ＤＢ
３２Ａ，３２Ｂ 制約拡張情報
３２ＡＭ、３２ＢＭ 制約拡張情報
３３ 過去実績ＤＢ
３３Ａ 過去実績情報
３４ 過去計画ＤＢ
３４Ａ，３４Ｂ 過去計画情報
３５Ａ 評価指標情報
３７ 検証シナリオ情報
１０１ 発電機
１０２ 変電所
１０３ 調相機器
１０４ 電力負荷
１０５ 外部電力系統
１０８ 情報通信ネットワーク
２００ 需給計画装置
【図1】
【図2】
【図3】
【図4】
【図5A】
【図5B】
【図6】
【図7A】
【図7B】
【図8】
【図9】
【図10A】
【図10B】
【図11A】
【図11B】
【図12】