(57)【要約】【課題】複数の地域系統を含む電力系統における負荷周波数制御の評価指標の公平性を効率よく簡便に評価できるようにする。【解決手段】情報処理装置１０に、複数の地域系統を含む系統モデルにおける各地域系統の系統パラメータを記憶させ、各地域系統に与えられる外乱をΔＰＤｉ（ｓ）、伝達関数をρｉｊ（ｓ）とした場合における地域要求量ΔＲｉ（ｓ）を定式化し、前記系統モデルについての評価指標ＣＰＳ１のＣＦ１の理論式を定式化し、この理論式の構成要素（Ｃij，Ｖｆij）をある地域系統にのみ外乱ΔＰＤｉ（ｓ）を与えた場合についてシミュレーションすることにより求め、求めた構成要素（Ｃij，Ｖｆij）と上記理論式とに基づいて、複数の前記地域系統に外乱を入力した場合におけるＣＦ１の値を求める。【選択図】図２
2008109810.tif 000002
【特許請求の範囲】
【請求項１】
複数の地域系統を含む電力系統における負荷周波数制御の評価指標の解析方法であって、
ＣＰＵ及びメモリを有する情報処理装置が、
前記電力系統をモデル化した系統モデルにおける前記各地域系統の系統パラメータを記憶するステップと、
前記各地域系統に与えられる外乱をΔＰＤｉ（ｓ）（i=1〜N，ｓはラプラス変数）、伝達関数をρｉｊ（ｓ）とした場合における地域要求量（Area Requirement）ΔＲｉ（ｓ）を、次式
2008109810.tif 000021
によって定義し、前記系統モデルについての評価指標ＣＰＳ１のＣＦ１の理論式を、ε１を基準値、Ｋを系統周波数特性定数、Ｐiを系統ｉの系統容量、αjを外乱の標準偏差の大きさを決める係数として次式
2008109810.tif 000022
で表した場合におけるＣij，Ｖｆijを、前記系統モデルについて、ある前記地域系統にのみ前記外乱ΔＰＤｉ（ｓ）を与えた場合についてシミュレーションすることにより求めるステップと、
求めた前記Ｃij，Ｖｆijと前記理論式とに基づいて、複数の前記地域系統に外乱を入力した場合における前記ＣＦ１の値を求めること、
を特徴とする負荷周波数制御の評価指標の解析方法。
【請求項２】
請求項１に記載の負荷周波数制御の評価指標の解析方法であって、
前記地域要求量ΔＲｉ（ｓ）は周波数偏差であること
を特徴とする負荷周波数制御の評価指標の解析方法。
【請求項３】
請求項１に記載の負荷周波数制御の評価指標の解析方法であって、
前記地域要求量ΔＲｉ（ｓ）は連系線潮流偏差であること
を特徴とする負荷周波数制御の評価指標の解析方法。
【請求項４】
請求項１に記載の負荷周波数制御の評価指標の解析方法であって、
前記系統パラメータには、前記各地域系統の容量、外乱の標準偏差、発電機の慣性定数、周波数特性定数、同期化係数、ＴＢＣパラメータ、速度調停率のうちの少なくともいずれかが含まれること
を特徴とする負荷周波数制御の評価指標の解析方法。
【請求項５】
請求項１に記載の負荷周波数制御の評価指標の解析方法であって、
前記シミュレーションにおいて、ある前記地域系統にのみ与えられる前記外乱ΔＰＤｉ（ｓ）は、標準白色外乱であること
を特徴とする負荷周波数制御の評価指標の解析方法。
【請求項６】
請求項１に記載の負荷周波数制御の評価指標の解析方法であって、
ガバナフリーとした場合における前記系統モデルと、ガバナロックとした場合における前記系統モデルの双方について前記ＣＦ１の値を求めること
を特徴とする負荷周波数制御の評価指標の解析方法。
【請求項７】
複数の地域系統を含む電力系統における負荷周波数制御の評価指標を解析する情報処理装置であって、
前記電力系統をモデル化した系統モデルにおける前記各地域系統の系統パラメータを記憶し、
前記各地域系統に与えられる外乱をΔＰＤｉ（ｓ）（i=1〜N，ｓはラプラス変数）、伝達関数をρｉｊ（ｓ）とした場合における地域要求量（Area Requirement）ΔＲｉ（ｓ）を、次式
2008109810.tif 000023
によって定義し、前記系統モデルについての評価指標ＣＰＳ１のＣＦ１の理論式を、ε１を基準値、Ｋを系統周波数特性定数、Ｐiを系統ｉの系統容量、αjを外乱の標準偏差の大きさを決める係数として次式
2008109810.tif 000024
で表した場合におけるＣij，Ｖｆijを、前記系統モデルについて、ある前記地域系統にのみ前記外乱ΔＰＤｉ（ｓ）を与えた場合についてシミュレーションすることにより求め、
求めた前記Ｃij，Ｖｆijと前記理論式とに基づいて、複数の前記地域系統に外乱を入力した場合における前記ＣＦ１の値を求めること、
を特徴とする情報処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の地域系統を含む電力系統における負荷周波数制御の評価指標の解析方法に関し、とくに負荷周波数制御の評価指標の公平性を効率よく簡便に評価できるようにするための技術に関する。
【背景技術】
【０００２】
近年の電力自由化に伴って、電力会社においては電力品質の維持が課題となっている。電力品質を決定する要素の一つとして、系統内の周波数を観測して系統内の発電量にフィードバックさせる、いわゆる負荷周波数制御（LFC：Load Frequency Control）がある。ここで電力系統は、通常は各電力会社間で連系運用がなされているため、各電力系統についての負荷周波数制御の性能評価が必要となるが、そのための評価指標として、従来より各地域の系統（以下、地域系統と称する。）の容量や各地域系統の動的な特性を反映したＡＲ標準偏差、ＡＣＥ標準偏差、ＭＡＣ、ＣＰＳなどが知られている。
【特許文献１】特開２００６−４２４５８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
ここで上記連系運用のもとでは、各地域系統の上記負荷周波数制御の性能を高い精度で公平に評価することが求められる。そして上記評価指標を用いた負荷周波数制御の性能評価を実施する場合には、上記各評価指標の公平性についても評価する必要がある。
【０００４】
しかしながら、各地域系統の連系によって負荷周波数制御が相互に影響を及ぼすため、一般に評価指標の公平性の評価は難しく、各評価指標の公平性を効率よく簡便に行う仕組みは存在しない。
【０００５】
本発明は、このような背景に鑑みてなされたもので、評価指標の公平性を効率よく簡便に評価することが可能な負荷周波数制御の評価指標の解析方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記目的を達成するための請求項１に記載の発明は、複数の地域系統を含む電力系統における負荷周波数制御の評価指標の解析方法であって、ＣＰＵ及びメモリを有する情報処理装置が、前記電力系統をモデル化した系統モデルにおける前記各地域系統の系統パラメータを記憶するステップと、前記各地域系統に与えられる外乱をΔＰＤｉ（ｓ）（i=1〜N，ｓはラプラス変数）、伝達関数をρｉｊ（ｓ）とした場合における地域要求量（Area Requirement）ΔＲｉ（ｓ）を、次式
2008109810.tif 000003
によって定義し、前記系統モデルについての評価指標ＣＰＳ１のＣＦ１の理論式を、ε１を基準値、Ｋを系統周波数特性定数、Ｐiを系統ｉの系統容量、αjを外乱の標準偏差の大きさを決める係数として次式
2008109810.tif 000004
で表した場合におけるＣij，Ｖｆijを、前記系統モデルについて、ある前記地域系統にのみ前記外乱ΔＰＤｉ（ｓ）を与えた場合についてシミュレーションすることにより求めるステップと、求めた前記Ｃij，Ｖｆijと前記理論式とに基づいて、複数の前記地域系統に外乱を入力した場合における前記ＣＦ１の値を求めることとする。
【０００７】
本発明では、ある地域系統にのみ外乱ΔＰＤｉ（ｓ）を与えてシミュレーションを行ってＣij，Ｖｆij（以下、基礎データと称する。）を求め、求めた前記Ｃij，Ｖｆijと前記理論式とに基づいて、複数の地域系統に外乱を入力した場合における評価指標ＣＰＳ１のＣＦ１の具体的な値を求めるようにしている。このため、評価指標ＣＰＳ１を算出しようとする電力系統が、複数の地域系統を含む複雑な電力系統である場合でも、シミュレーションの実施回数をできる限り少なくして、複数の地域系統に外乱を入力した場合におけるＣＦ１の値を理論的に求めることが可能である。そして、このようにして求めたＣＦ１の値を用いて負荷周波数制御の評価指標の性能を簡便に評価することができる。
【０００８】
また本発明のうち請求項２に記載の発明は、請求項１に記載の負荷周波数制御の評価指標の解析方法であって、前記地域要求量ΔＲｉ（ｓ）は周波数偏差であることとする。
【０００９】
また本発明のうち請求項３に記載の発明は、請求項１に記載の負荷周波数制御の評価指標の解析方法であって、前記地域要求量ΔＲｉ（ｓ）は連系線潮流偏差であることとする。
【００１０】
また本発明のうち請求項４に記載の発明は、請求項１に記載の負荷周波数制御の評価指標の解析方法であって、前記系統パラメータには、前記各地域系統の容量、外乱の標準偏差、発電機の慣性定数、周波数特性定数、同期化係数、ＴＢＣパラメータ、速度調停率のうちの少なくともいずれかが含まれることとする。
【００１１】
また本発明のうち請求項５に記載の発明は、請求項１に記載の負荷周波数制御の評価指標の解析方法であって、前記シミュレーションにおいて、ある前記地域系統にのみ与えられる前記外乱ΔＰＤｉ（ｓ）は、標準白色外乱であることとする。
【００１２】
また本発明のうち請求項６に記載の発明は、請求項１に記載の負荷周波数制御の評価指標の解析方法であって、ガバナフリーとした場合における前記系統モデルと、ガバナロックとした場合における前記系統モデルの双方について前記ＣＦ１の値を求めることとする。
【００１３】
また本発明のうち請求項７に記載の発明は、複数の地域系統を含む電力系統における負荷周波数制御の評価指標を解析する情報処理装置であって、前記電力系統をモデル化した系統モデルにおける前記各地域系統の系統パラメータを記憶し、前記各地域系統に与えられる外乱をΔＰＤｉ（ｓ）（i=1〜N，ｓはラプラス変数）、伝達関数をρｉｊ（ｓ）とした場合における地域要求量（Area Requirement）ΔＲｉ（ｓ）を、次式
2008109810.tif 000005
によって定義し、前記系統モデルについての評価指標ＣＰＳ１のＣＦ１の理論式を、ε１を基準値、Ｋを系統周波数特性定数、Ｐiを系統ｉの系統容量、αjを外乱の標準偏差の大きさを決める係数として次式
2008109810.tif 000006
で表した場合におけるＣij，Ｖｆijを、前記系統モデルについて、ある前記地域系統にのみ前記外乱ΔＰＤｉ（ｓ）を与えた場合についてシミュレーションすることにより求め、求めた前記Ｃij，Ｖｆijと前記理論式とに基づいて、複数の前記地域系統に外乱を入力した場合における前記ＣＦ１の値を求めることとする。
【発明の効果】
【００１４】
本発明によれば、複数の地域系統を含む電力系統における負荷周波数制御の評価指標の公平性を効率よく簡便に評価することができる。
【発明を実施するための最良の形態】
【００１５】
以下、本発明の一実施形態につき詳細に説明する。図１に本発明の一実施形態として説明する情報処理装置１０のハードウエア構成を示している。同図に示すように、情報処理装置１０は、パーソナルコンピュータやワークステーション、汎用機などのコンピュータであって、ＣＰＵ１１、ＲＡＭ・ＲＯＭなどのメモリ１２、ハードディスク１３、キーボードやマウスなどの入力装置１４、液晶ディスプレイやブラウン管ディスプレイなどの表示装置１５、ＬＡＮ（Local Area Network）などの通信ネットワークに接続するための機能を提供する通信インタフェース１６を備えている。
【００１６】
情報処理装置１０は、例えば異なる電力会社によって運用される複数の地域系統を含む電力系統における負荷周波数制御の評価指標のうち、ＮＥＲＣ(North America Electric Reliability Council：北米電力信頼度評議会)によって提案された評価指標ＣＰＳ１（Control Performance Standard 1）を求める機能を備えている。図２に評価指標ＣＰＳ１を解析するための情報処理装置１０が有する機能を示している。なお、これらの機能はＣＰＵ１１がメモリ１２に記憶されているプログラムを実行することによって実現される。
【００１７】
図２に示す各種機能のうち系統モデル記憶部２１は、上記評価指標ＣＰＳ１を解析する際の基本となる電力系統のモデル（以下、系統モデルと称する。）を記憶する。具体的には、系統モデル記憶部２１は、上記基本となる系統モデルとして図３Ａに示す３つの地域系統３１を含んで構成される系統モデル３０を記憶している。また系統モデル記憶部２１は、系統モデル３０を構成している各地域の系統モデルとして、各地域系統に共通する本質的な要素のみを考慮した、図３Ｂに示す系統モデル（以下、地域系統モデル３１０と称する。）を記憶している。なお、同図に示す地域系統モデル３１０において、ＰＥＬＤiは経済負荷配分指令値合計であり、ＰＤｉは負荷変動であり、ΔＰＴｉは連系線潮流
2008109810.tif 000007
であり、Ｔｉｊは同期化定数である。またＫＡＲｉは系統周波数特性定数、Δｆｉは周波数偏差であり、ΔＲｉは地域要求量であり、ＫＨｉはLFC比例ゲインであり、ＴＨｉはLFC積分時定数である。またＭｉは慣性定数であり、Ｄｉは系統負荷定数であり、Ｔgoviはガバナ時定数であり、Ｋgoviはガバナゲイン(=100/R)であり、Rは速度調定率である。
【００１８】
系統パラメータ記憶部２２は、各地域系統３１について設定されたパラメータ（以下、系統パラメータと称する。）を記憶する。図３Ｃは本実施形態において設定した系統パラメータである。
【００１９】
情報処理装置１０によって評価指標ＣＰＳ１を求める際には、上記系統モデル３０に加える外乱を変化させた各種のケース（Case）を設定し、設定した各ケースについて評価指標ＣＰＳ１を算出する。また情報処理装置１０によって評価指標ＣＰＳ１を効率よく簡便に算出すべく、情報処理装置１０は、ある地域系統３１にのみ外乱を与えたケース（以下、ベースケースと称する。）についてシミュレーションを行うことにより得られる評価指標ＣＰＳ１（以下、基礎データと称する。）と、後述のようにして導出される理論式とに基づいて、各地域系統３１に系統負荷の変動などの外乱が常時発生している実際のケースを想定したケースである、複数の地域系統３１に同時に外乱を与えたケース（以下、一般ケースと称する。）における評価指標ＣＰＳ１を算出する。
【００２０】
図２において、基礎データ生成部２３は、上記基礎データを生成する機能である。理論値算出部２４は、上記理論式によって評価指標ＣＰＳ１を算出する機能である。外乱生成部２５は、上記評価指標ＣＰＳ１の算出に際して地域系統３１に与える各種の外乱を生成する。
【００２１】
なお、情報処理装置１０は、上記評価指標ＣＰＳ１の算出に際し、ガバナの状態（ガバナフリー、ガバナロック）を考慮する。ガバナ状態設定部２６は、上記評価指標ＣＰＳ１の算出に際し、ガバナの状態を反映させる。ガバナ状態設定部２６によるガバナ状態の具体的な設定方法については後述する。
【００２２】
解析結果記憶部２７は、理論値算出部２４によって求められた評価指標ＣＰＳ１の値を記憶する。解析結果記憶部２７が記憶している値は適宜表示装置１５に出力されてユーザの利用に供される。
【００２３】
＝＝理論式の導出＝＝
次に、理論値算出部２４が基礎データに基づいて評価指標ＣＰＳ１を算出する際に用いられる前述した理論式について説明する。
【００２４】
＜評価指標ＣＰＳ１＞
ＣＰＳ１では、まず基準となる年を設定し、その1年間の周波数偏差の実績値により次式で定義される基準値ε１を決定する。  
2008109810.tif 000008
なお、系統の周波数の平均が基準周波数に等しい場合、ε１は基準周波数に対する標準偏差を表す。ＣＰＳ１では、核となる要素として、次のＣＦ１(Compliance Factor 1)が定義されている。  
2008109810.tif 000009
ここで、Ｂｉは、系統ｉの周波数バイアス値[MW/Hz]である。なお、上式における記号[]１は、一年間の平均を意味する。ＡＣＥｉは、系統ｉの発電電力の過不足量（調整残）であり、次式で定義される。
【００２５】
2008109810.tif 000010
【００２６】
＜伝達関数行列＞
本実施形態では、各地域系統モデル３１０に与える外乱をΔＰＤｉ（ｓ）（i=1〜N，ｓはラプラス変数）、伝達関数をρij（ｓ）とした場合における地域要求量（Area Requirement）ΔＲｉ（ｓ）を次式で定義する。  
2008109810.tif 000011
ここで図３Ｂに示した地域系統モデル３１０について、式（６）に基づき周波数偏差と連系線潮流偏差の伝達関数行列を求めると、次式（７）、（８）のようになる。  
2008109810.tif 000012
【００２７】
＜理論式＞
ＣＦ１の定義式（５）は、式（３）を代入すると次のような形に変形できる。  
2008109810.tif 000013
ここで分散をＶ[;]、平均をＥ[;]と表すと、上式（９）は次式（１０）で表される。  
2008109810.tif 000014
【００２８】
また定式化で導出した周波数偏差と連系線潮流偏差の伝達関数行列（７）、（８）を用い、Ｅ［Δｆi・ΔＰｎ］（系統iの周波数偏差と連系線潮流偏差の積の平均値）、Ｖ［Δｆi］（系統iの周波数偏差の分散）の予測式を求めると、式（１０）は結局、
2008109810.tif 000015
と表される。本実施形態において、理論値算出部２４は、式（１１）を理論式として求められるＣＦ１の平均値を評価指標ＣＰＳ１の理論値として求める。なお、Ｃijは系統ｊに標準外乱が入った時の系統iにおけるΔｆｉ・ΔＰiの平均値である。またＶｆijは系統ｊに標準外乱が入った時の系統ｉにおけるΔｆｉの分散であり、αjは外乱の標準偏差の大きさを決める係数（倍率）である。
【００２９】
ここでＣij，Ｖｆijは、上述した基礎データに対応し、基礎データ生成部２３は、ベースケースについてシミュレーションにより基礎データＣij，Ｖｆijを求める。図４Ａ乃至図４Ｃに、図３Ａに示す系統モデル３０についてのベースケースを示す。
【００３０】
なお、Ｎ系統が連系された系統モデルにはＮ個のベースケースが存在し、これらのＮ個すべてのベースケースについて基礎データを収集しなければ、上記理論式を用いて結果を予測することはできない。従って、シミュレーションは最低でＮ回（ベースケース数）必要となる。また１回のシミュレーションで系統数分（Ｎ個）の基礎データが得られるので、Ｎ回のシミュレーションによって上記の要素Ｃｉｊ，ＶｆｉｊのそれぞれについてＮ（ベースケース数）×Ｎ（系統数）個の基礎データが得られることになる。
【００３１】
＝＝ガバナ状態の設定方法＝＝
本実施形態ではガバナ（調速機）の特性を次式で表す。なお、ガバナフリー発電機は１次遅れ系で時定数が１．０[s]の特性を有しているものとする。  
2008109810.tif 000016
上式においてＫgoviは系統iのガバナゲイン [%MW/0.1Hz]である。ガバナゲインは次のようにして算出する。すなわち、基準周波数をFn[Hz]とすると全負荷遮断時の周波数変化ΔＦ（＝Ｆ0−Ｆn）は、
2008109810.tif 000017
となる。ただし、Ｒは速度調定率[%Hz/%MW]である。ガバナゲインは、0.1Hzの周波数変化に対して系統容量の何％出力変化するかを表すものであるから、全負荷遮断時の周波数変化により次式が求まる。  
2008109810.tif 000018
従って、式（１３）,（１４）より次式が導かれる。  
2008109810.tif 000019
ここではＲ＝4.5，Ｆｎ＝60であるので、これらを上式に代入すると、Ｋgovi＝3.73[%MW/0.1Hz]となる。従ってガバナフリー発電機の伝達関数は次式で表される。
【００３２】
2008109810.tif 000020
【００３３】
＝＝シミュレーション＝＝
以上の構成からなる情報処理装置１０によって、図４Ａ乃至図４Ｃに示す系統モデル３０およびパラメータを設定し、一般ケースについての評価指標ＣＰＳ１の理論値としてＣＦ１の平均値を求めた。まず図４Ａ乃至図４Ｃに示した系統モデル３０およびパラメータを設定し、次のベースケース（Case1〜Case3）についてシミュレーションを行って基礎データを生成した。  
Case1：系統1にのみ標準白色外乱を発生させた場合
Case2：系統2にのみ標準白色外乱を発生させた場合
Case3：系統3にのみ標準白色外乱を発生させた場合
【００３４】
なお、シミュレーションは全ケースにおいて3600×24×365[s]の１年間を対象とした。また系統構成の違いによる影響は非常に少ないと仮定し、シミュレーションはループ型モデルについてのみ行っている。またシミュレーションはガバナロック、ガバナフリーのそれぞれについて行った。
【００３５】
外乱はＬＦＣの制御対象であるフリンジ成分(数分〜数十分)としており、本実施形態では１分〜２０分周期の外乱をシミュレーションに用いた。外乱は負荷変動が正規分布に基づく白色外乱とし、帯域フィルタを用いて1分〜20分以外の周期外乱を除去し、外乱の標準偏差がγ√（Ｐｉ）となるように補正した。各地域系統の外乱は相関が無いものとした。
【００３６】
図５Ａ及び図５Ｂにガバナロックの場合の基礎データＣｉｊ，Ｖｆｉｊを、また図６Ａ，図６Ｂにガバナフリーの場合の基礎データＣｉｊ，Ｖｆｉｊを示している。なお、これらの図において、分散Ｖｆｉｊは標準偏差√（Ｖｆｉｊ） [Hz]として表記してある。
【００３７】
次に、前述した評価指標ＣＰＳ１の理論式（１１）を用い、以下の一般ケース（Case4〜Case7）についてＣＦ１の平均値を求めた。  
Case4：系統1，2に標準白色外乱を発生させた場合
Case5：系統1，3に標準白色外乱を発生させた場合
Case6：系統2，3に標準白色外乱を発生させた場合
Case7：系統1，2，3すべてに標準白色外乱を発生させた場合
【００３８】
結果を図７Ａ及び図７Ｂに示す。なお、図７Ａはガバナロックの場合であり、図７Ｂはガバナフリーの場合である。また上段はシミュレーションの結果であり、括弧内の値は理論値(予測値)である。
【００３９】
理論値をシミュレーションによる値と比較すると、多少の誤差はあるものの、ガバナロックとガバナフリーのどちらの場合についても数パーセントの誤差しかなく、理論式によって有意な値が得られることが確認できた。なお、誤差が生じた原因は、基礎データの収集に際して行ったシミュレーションが有限時間内についてのものであるため、白色外乱の平均が0にならないこと、シミュレーションの値に外乱の相関が現れていることなどが考えられる。これらの誤差はさらに長時間のシミュレーションを実施することにより限りなく0に近づけることができるものと考えられる。
【００４０】
以上のように、複数の地域系統からなる複雑な電力系統であっても、評価指標ＣＰＳ１について理論式によって高い精度で効率よく簡便に解析を行うことができることを検証できた。
【００４１】
＝＝評価指標ＣＰＳ１の性能評価＝＝
次に、Case1〜Case7の結果に基づき、評価指標ＣＰＳ１の性能評価を行った。各ケースについて求めたＣＦ１の平均値を図８Ａ及び図８Ｂに示す。図８Ａはガバナロックの場合であり、図８Ｂはガバナフリーの場合である。図中、外乱が与えられている系統には下線を付してある。
【００４２】
図８Ａのガバナロックの場合について、外乱が与えられている系統の値を表全体で比較するといずれもほぼ0.014[p.u]で一致していることがわかる。これは外乱が与えられている系統については評価指標ＣＰＳ１による評価が公平になされていることを意味する。一方、下線の引かれていない系統、すなわち外乱が与えられていない系統の値を表全体で比較すると、0.00002〜0.00009[p.u]の幅があるが、いずれの絶対値も非常に小さく、評価指標ＣＰＳ１の公平性の判断には殆ど影響を与えないことがわかる。
【００４３】
図８Ｂのガバナフリーの場合について、外乱が与えられている系統の値を表全体で比較すると、図８Ａに比べて値にばらつきがあることがわかる。これは発生させている外乱の大きさが同じであっても、個々の発生条件によって評価指標ＣＰＳ１の値が変わってくることを意味している。一方、下線の引かれていない系統、すなわち外乱が与えられていない系統の値を表全体で比較すると、同じケース内では値が一致しているが、異なるケース間では値が異なっており、評価指標ＣＰＳ１の公平性に問題があることがわかる。
【００４４】
以上に説明したように、本実施形態の情報処理装置１０は、ベースケースについてのみシミュレーションを行って基礎データを生成し、生成した基礎データと、理論式とによって一般ケースについての評価指標を取得することにより、計算負荷を抑えて複数の地域系統を含む電力系統における負荷周波数制御の評価指標を生成するので、評価指標の公平性を効率よく簡便に判断することができる。また本実施形態の情報処理装置１０は、より公平な判断が可能となるように評価指標を改善するための手段としても用いることができる。
【００４５】
以上の実施形態の説明は本発明の理解を容易にするためのものであり、本発明を限定するものではない。本発明はその趣旨を逸脱することなく、変更、改良され得ると共に本発明にはその等価物が含まれることは勿論である。
【図面の簡単な説明】
【００４６】
【図１】本発明の一実施形態として説明する情報処理装置１０のハードウエア構成を示す図である。
【図２】本発明の一実施形態として説明する情報処理装置１０が備える各種機能を示す図である。
【図３Ａ】本発明の一実施形態として説明する地域系統３１を含んで構成される系統モデル３０を示す図である。
【図３Ｂ】本発明の一実施形態として説明する地域系統モデル３１０を示す図である。
【図３Ｃ】本発明の一実施形態として説明する系統パラメータを示す図表である。
【図４Ａ】本発明の一実施形態として説明するベースケース（Case1）の一例を示す図である。
【図４Ｂ】本発明の一実施形態として説明するベースケース（Case2）の一例を示す図である。
【図４Ｃ】本発明の一実施形態として説明するベースケース（Case3）の一例を示す図である。
【図５Ａ】本発明の一実施形態として説明するガバナロックの場合の基礎データＣｉｊを示す図である。
【図５Ｂ】本発明の一実施形態として説明するガバナロックの場合の基礎データＶｆｉｊを示す図である。
【図６Ａ】本発明の一実施形態として説明するガバナフリーの場合の基礎データＣｉｊを示す図である。
【図６Ｂ】本発明の一実施形態として説明するガバナフリーの場合の基礎データＶｆｉｊを示す図である。
【図７Ａ】本発明の一実施形態として説明するガバナロックの場合のＣＦ１の平均値を示す図である。
【図７Ｂ】本発明の一実施形態として説明するガバナフリーの場合のＣＦ１の平均値を示す図である。
【図８Ａ】本発明の一実施形態として説明するガバナロックの場合のＣＦ１の平均値を示す図である。
【図８Ｂ】本発明の一実施形態として説明するガバナフリーの場合のＣＦ１の平均値を示す図である。
【符号の説明】
【００４７】
１０  情報処理装置
２１  系統モデル記憶部
２２  系統パラメータ記憶部
２３  基礎データ生成部
２４  理論値算出部
２５  外乱生成部
２６  ガバナ状態設定部
２７  解析結果記憶部
【図１】
2008109810.tif 000025
【図２】
2008109810.tif 000026
【図３Ａ】
2008109810.tif 000027
【図３Ｂ】
2008109810.tif 000028
【図３Ｃ】
2008109810.tif 000029
【図４Ａ】
2008109810.tif 000030
【図４Ｂ】
2008109810.tif 000031
【図４Ｃ】
2008109810.tif 000032
【図５Ａ】
2008109810.tif 000033
【図５Ｂ】
2008109810.tif 000034
【図６Ａ】
2008109810.tif 000035
【図６Ｂ】
2008109810.tif 000036
【図７Ａ】
2008109810.tif 000037
【図７Ｂ】
2008109810.tif 000038
【図８Ａ】
2008109810.tif 000039
【図８Ｂ】
2008109810.tif 000040