JP5602273B1

JP5602273B1 - 電力系統シミュレータならびに電力系統のシミュレーションプログラムおよびシミュレーション方法

Info

Publication number: JP5602273B1
Application number: JP2013088511A
Authority: JP
Inventors: 裕介福原; 学倉本; 道雄尾野本
Original assignee: Chugoku Electric Power Co Inc
Current assignee: Chugoku Electric Power Co Inc
Priority date: 2013-04-19
Filing date: 2013-04-19
Publication date: 2014-10-08
Anticipated expiration: 2033-04-19
Also published as: JP2014212650A

Abstract

【課題】最過酷事故の事故条件を選定することにより、検討すべき事故条件を絞り込み、検討時間を短縮する。
【解決手段】複数の事故条件から１つをシミュレーション対象として順次選択する条件選択部と、シミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定する初期値設定部と、シミュレーション対象の事故条件に対して、潮流量をパラメータとして潮流計算を行う潮流計算部と、潮流計算部の計算結果に基づいて、電力系統の過渡安定度の判定を行う安定度判定部と、複数の事故条件のそれぞれに対する安定度判定部の判定結果に基づいて、最過酷事故の事故条件を選定する最過酷事故選定部と、を有し、初期値設定部は、現在のシミュレーション対象の事故条件に対する安定度判定部の判定結果に基づいて、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定する。
【選択図】図１

Description

本発明は、電力系統シミュレータ、ならびに電力系統のシミュレーションプログラムおよびシミュレーション方法に関する。

電力系統で事故が発生した場合の影響を確認するため、電力系統シミュレータを用いて対象の系統条件および事故条件における潮流変化をシミュレートして過渡安定度の判定を行うことにより、安定的に流し得る潮流量の安定限界の探索が行われている（例えば特許文献１を参照）。電力系統シミュレータのこれらの機能を実現するためのシミュレーションプログラムとしては、Ｌ法と呼ばれる潮流計算プログラムや、Ｙ法と呼ばれる電力系統安定度解析プログラム（電力系統動特性過渡安定度解析プログラム）などが従来から広く知られている。

このように、対象となる事故条件をシミュレートして、事故発生後も電力系統を適正に維持することができる場合には安定と判定し、適正に維持することができない場合には不安定と判定している。また、安定と不安定との境界である安定限界を用いて電力系統を安定に運用している。

特開２０１１−１１５００３号公報

上記のようなシミュレーションプログラムを用いて多くの事故条件について検討する必要がある場合、検討者は、事故条件ごとにコンピュータを操作しつつ検討したり、事故条件を一括で設定したうえでシミュレーションを行って検討したりすることとなる。しかしながら、前者の方法では、検討者がコンピュータの操作を長時間行う必要がある一方、後者の方法では、一括で設定した事故条件が不足しないよう条件数を多くする必要があるため、シミュレーション時間が長くなる。そのため、これらの方法は、緊急で検討を行って短時間で結果を求める必要がある場合には不向きである。

前述した課題を解決する主たる本発明は、電力系統における複数の事故条件をシミュレートして、前記複数の事故条件のうち、前記電力系統に安定的に流し得る潮流量の安定限界が最小となる最過酷事故の事故条件を選定する電力系統シミュレータであって、前記複数の事故条件から１つをシミュレーション対象として順次選択する条件選択部と、前記シミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定する初期値設定部と、前記シミュレーション対象の事故条件に対して、潮流量をパラメータとして潮流計算を行う潮流計算部と、前記潮流計算部の計算結果に基づいて、前記電力系統の過渡安定度の判定を行う安定度判定部と、前記複数の事故条件のそれぞれに対する前記安定度判定部の判定結果に基づいて、前記最過酷事故の事故条件を選定する最過酷事故選定部と、を有し、前記初期値設定部は、現在のシミュレーション対象の事故条件に対する前記安定度判定部の判定結果に基づいて、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定することを特徴とする電力系統シミュレータである。

本発明の他の特徴については、添付図面及び本明細書の記載により明らかとなる。

本発明によれば、すべての事故条件のうち、電力系統に安定的に流し得る潮流量の安定限界が最小となる最過酷事故の事故条件を選定することにより、検討すべき事故条件を絞り込み、検討時間を短縮することができる。

本発明の一実施形態における電力系統シミュレータの構成を示すブロック図である。本発明の一実施形態における電力系統シミュレータの各部の機能をコンピュータに実現させるためのプログラムの動作を説明するフローチャートである。本発明の一実施形態における電力系統シミュレータの各部の機能をコンピュータに実現させるためのプログラムの動作を説明するフローチャートである。本発明の一実施形態における電力系統のシミュレーション方法による最過酷事故の選定方法を説明する図である。従来の電力系統のシミュレーション方法による最過酷事故の選定方法を説明する図である。

本明細書および添付図面の記載により、少なくとも以下の事項が明らかとなる。

＝＝＝電力系統シミュレータの構成＝＝＝
以下、図１を参照して、本発明の一実施形態における電力系統シミュレータの構成について説明する。

図１に示されている電力系統シミュレータ１０は、事故ケース選択部１０１、初期値設定部１０２、潮流計算部１０３、安定度判定部１０４、最過酷事故選定部１０５、入力部１０６、出力部１０７、記憶部１０８、およびバス１０９含んで構成されている。また、事故ケース選択部１０１、潮流計算部１０３、最過酷事故選定部１０５、入力部１０６、出力部１０７、および記憶部１０８は、バス１０９を介して互いに接続されている。なお、電力系統シミュレータ１０の機能は、入力部１０６、出力部１０７、記憶部１０８、およびバス１０９を備えるコンピュータ１００によって実現することができる。

＝＝＝電力系統のシミュレーション方法＝＝＝
以下、図２ないし図４を適宜参照して、本実施形態における電力系統シミュレータによる電力系統のシミュレーション方法について説明する。なお、以下の説明においては、電力系統における事故ケース１ないし事故ケースＫ（Ｋ個の事故条件）を当該順序でシミュレートして、最過酷事故ケース（最過酷事故の事故条件）を選定するものとする。ここで、最過酷事故ケースとは、事故ケース１ないし事故ケースＫのうち、対象の電力系統に安定的に流し得る潮流量の安定限界が最小となる事故ケースである。

前述したように、電力系統シミュレータ１０の機能は、コンピュータ１００によって実現することができる。例えば、コンピュータ１００にシミュレーションプログラムを実行させることによって、事故ケース選択部１０１、初期値設定部１０２、潮流計算部１０３、安定度判定部１０４、および最過酷事故選定部１０５に相当する処理を実行することができる。図２および図３は、これら各部に相当する処理をコンピュータに実行させるためのプログラムの動作を示している。

図２は、最初のシミュレーション対象である事故ケース１に対する動作を示している。

プログラムの実行が開始されると、まず、シミュレーションデータを作成するとともに、事故ケース１ないし事故ケースＫを設定する（Ｓ１１）。シミュレーションデータは、電力系統の状態を示す系統情報に基づいて作成され、この系統情報は、例えば、入力部１０６を介してオンラインで取得され、さらにバス１０９を介して記憶部１０８に記憶されている。一方、事故ケース１ないし事故ケースＫは、例えば、記憶部１０８に記憶されている事故ケースからＫ個を選択してもよく、また、キーボードやマウスなどの入力部１０７を用いて入力してもよい。

次に、事故ケースｋ（＝１）に対して、事故ケース選択部１０１（条件選択部）に相当する条件選択処理Ｓ１２を実行し、シミュレーション対象として選択する。なお、このとき、事故ケース１を暫定的に最過酷事故ケースに設定する。また、潮流計算におけるパラメータとしての潮流量を所定の初期値に設定する（Ｓ１３）。

Ｓ１４ないしＳ１６の処理は、潮流計算部１０３および安定度判定部１０４に相当する機能（潮流計算処理および安定度判定処理）を実現する。すなわち、シミュレーションを実行し、潮流量の初期値に対して潮流計算を行うとともに、潮流計算結果に基づいて電力系統の過渡安定度（安定または不安定）の判定を行う（Ｓ１４）。そして、安定限界が求まるまでの間（Ｓ１５：ＮＯ）、潮流量を変更しながら（Ｓ１６）、シミュレーションを実行する（Ｓ１４）。

例えば、潮流量の初期値に対して過渡安定度の判定を行い、安定と判定した場合には、不安定と判定されるまでの間、所定のステップ幅で潮流量を順次増加させることによって安定限界を求めることができる。一方、潮流量の初期値に対して不安定と判定した場合には、安定と判定されるまでの間、所定のステップ幅で潮流量を順次減少させることによって安定限界を求めることができる。また、潮流量の初期値からバイナリサーチ（二分探索）を行って安定限界の探索を行ってもよい。バイナリサーチでは、潮流量を変更するステップ幅を順次１／２倍ずつ小さくしつつ、安定と判定した場合には潮流量を増加させ、不安定と判定した場合には潮流量を減少させる。

安定限界が求まると（Ｓ１５：ＹＥＳ）、初期値設定部１０２に相当する初期値設定処理Ｓ１７を実行し、安定度判定処理の判定結果に基づいて、潮流量の初期値を更新する。より具体的には、事故ケース１（現在のシミュレーション対象）に対する安定度判定処理によって不安定と判定された最小の潮流量を、事故ケース２（次のシミュレーション対象）をシミュレートする際の潮流量の初期値として設定する。

ここで、図２に示した事故ケース１に対する動作の具体例を図４に示す。図４においては、安定との判定結果を○で示し、不安定との判定結果を×で示しており、安定限界では安定と判断されるものとする。

図４の事故ケース１において、潮流量の初期値（潮流計算：１回目）に対して安定と判定している。そのため、所定のステップ幅で潮流量を順次増加させ（潮流計算：２〜４回目）、２回目および３回目では安定と判定し、４回目で初めて不安定と判定して、安定限界を求めている。そして、不安定と判定された最小の潮流量（潮流計算：４回目）を、新たな潮流量の初期値として設定している。

図３は、２番目以降のシミュレーション対象である事故ケース２ないし事故ケースＫに対する動作を示している。

まず、ｋをインクリメントしたうえで（Ｓ２１）、条件選択処理Ｓ２２を実行し、事故ケースｋ（＞１）をシミュレーション対象として選択する。また、潮流量を、初期値設定処理Ｓ１７によって設定された初期値に設定する（Ｓ２３）。

Ｓ２４ないしＳ２６の処理は、潮流計算部１０３および安定度判定部１０４に相当する機能（潮流計算処理および安定度判定処理）を実現する。すなわち、シミュレーションを実行し、潮流量の初期値に対して潮流計算を行うとともに、過渡安定度の判定を行い（Ｓ２４）、不安定と判定した場合（Ｓ２５：ＮＯ）には、安定と判定されるまでの間、所定のステップ幅で潮流量を順次減少させて（Ｓ２６）、安定限界を求める。一方、潮流量の初期値に対して安定と判定した場合（Ｓ２５：ＹＥＳ）には、事故ケースｋに対してそれ以上のシミュレーションを実行せず、安定限界の探索も行わない。

次に、ｋ＝Ｋの終了条件が成立していなければ（Ｓ２７：ＮＯ）、Ｓ２８の判定を行う。ここで、潮流量の初期値に対して安定と判定し（Ｓ２５：ＹＥＳ）、安定限界の探索を行わなかった場合には、潮流量が初期値のままである（Ｓ２８：ＹＥＳ）ため、そのままＳ２１に戻る。

一方、潮流量の初期値に対して不安定と判定し（Ｓ２５：ＮＯ）、安定限界を求めた場合には、潮流量が初期値から減少している（Ｓ２８：ＮＯ）ため、初期値設定処理Ｓ２９を実行し、安定度判定処理の判定結果に基づいて、潮流量の初期値を更新する。より具体的には、事故ケースｋ（現在のシミュレーション対象）に対する安定度判定処理によって不安定と判定された最小の潮流量を、事故ケースｋ＋１（次のシミュレーション対象）をシミュレートする際の潮流量の初期値として設定して、Ｓ２１に戻る。なお、このとき、現在の事故ケースｋを暫定的に最過酷事故ケースに設定する。

以上の動作を、シミュレーション対象の事故ケースを順次選択しながら事故ケースＫまで繰り返し、ｋ＝Ｋの終了条件が成立すると（Ｓ２７：ＹＥＳ）処理を終了する。このとき、安定度判定処理によって最後に安定限界を求めた事故ケース、すなわち、Ｓ２９において最後に更新された暫定的な最過酷事故ケースを、最終的に最過酷事故ケースとして選定する。

ここで、図３に示した事故ケース２ないし事故ケースＫに対する動作の具体例を図４に示す。

図４の事故ケース２において、事故ケース１の初期値設定処理Ｓ１７によって設定された潮流量の初期値（潮流計算：１回目）に対して安定と判定している。そのため、事故ケース２に対してそれ以上のシミュレーションを実行せず、安定限界の探索も行っていない。また、潮流量の初期値を更新することなく、事故ケース３に移行する。

図４の事故ケース３において、事故ケース２の場合と同じ潮流量の初期値（潮流計算：１回目）に対して不安定と判定している。そのため、所定のステップ幅で潮流量を順次減少させ（潮流計算：２〜３回目）、２回目では安定と判定し、３回目で初めて不安定と判定して、安定限界を求めている。そして、不安定と判定された最小の潮流量（潮流計算：３回目）を、新たな潮流量の初期値として設定している。

そして、事故ケース４ないし事故ケースＫのいずれにおいても、事故ケース３の初期値設定処理Ｓ２９によって設定された潮流量の初期値（潮流計算：１回目）に対して安定と判定した場合には、事故ケース３が最終的に最過酷事故ケースとして選定されることとなる。

このように、本実施形態の電力系統のシミュレーション方法では、事故ケースｋ（現在のシミュレーション対象）に対する安定度判定結果に基づいて、事故ケースｋ＋１（次のシミュレーション対象）をシミュレートする際の潮流量の初期値を設定している。そのため、別の事故ケースの安定限界付近から効率的に安定限界の探索を行うことができ、さらに、各事故ケースに対する安定度判定結果に基づいて、最過酷事故ケースを選定することにより、検討すべき事故条件を絞り込み、検討時間を短縮することができる。

また、より具体的には、事故ケースｋに対する安定度判定処理によって不安定と判定された最小の潮流量を、事故ケースｋ＋１をシミュレートする際の潮流量の初期値として設定している。そのため、初期値設定処理によって設定された潮流量の初期値に対して安定と判定した場合には、この事故ケースに対して安定限界の探索を行わないことにより、シミュレーション回数を削減して、より短時間で最過酷事故ケースを選定することができる。

ここで、比較のため、従来の電力系統のシミュレーション方法による最過酷事故の選定方法を図５に示す。図５に示すように、事故ケース２ないし事故ケースＫにおいても常に安定限界の探索を行う場合、事故ケース２や事故ケースＫにおいて、本実施形態の電力系統のシミュレーション方法よりもシミュレーション回数が多くなっている。そのため、本実施形態の電力系統のシミュレーション方法は、緊急で検討を行って短時間で結果を求める必要がある場合により適している。

前述したように、電力系統シミュレータ１０において、事故ケース１ないし事故ケースＫ（複数の事故条件）から１つをシミュレーション対象として順次選択しながら、潮流量をパラメータとして潮流計算を行うとともに、潮流計算結果に基づいて電力系統の過渡安定度の判定を行い、各事故ケースに対する安定度判定結果に基づいて最過酷事故ケースを選定することによって、検討すべき事故条件を絞り込み、検討時間を短縮することができる。

また、事故ケースｋ（現在のシミュレーション対象）に対する安定度判定結果に基づいて、事故ケースｋ＋１（次のシミュレーション対象）をシミュレートする際の潮流量の初期値を設定することによって、別の事故ケースの安定限界付近から効率的に安定限界の探索を行うことができる。より具体的には、事故ケースｋに対する安定度判定処理によって不安定と判定された最小の潮流量を、事故ケースｋ＋１をシミュレートする際の潮流量の初期値として設定し、初期値設定処理によって設定された潮流量の初期値に対して安定と判定した場合には、当該事故ケースに対して安定限界の探索を行わないことによって、シミュレーション回数を削減して、より短時間で最過酷事故ケースを選定することができる。

また、系統情報や事故ケースが記憶されている記憶部１０８をさらに備えることによって、系統情報をオンラインで取得して記憶部１０８に記憶させておき、これを読み出してシミュレーションデータを作成することができ、記憶部１０８に記憶されている事故ケースからＫ個を読み出して事故ケース１ないし事故ケースＫを設定することができる。

また、コンピュータに、電力系統シミュレータ１０の各部に相当する処理を実行させるためのプログラムにおいて、事故ケース１ないし事故ケースＫから１つをシミュレーション対象として順次選択しながら、潮流量をパラメータとして潮流計算を行うとともに、潮流計算結果に基づいて電力系統の過渡安定度の判定を行い、各事故ケースに対する安定度判定結果に基づいて最過酷事故ケースを選定することによって、検討すべき事故条件を絞り込み、検討時間を短縮することができる。

また、事故ケース１ないし事故ケースＫから１つをシミュレーション対象として順次選択しながら、潮流量をパラメータとして潮流計算を行うとともに、潮流計算結果に基づいて電力系統の過渡安定度の判定を行い、各事故ケースに対する安定度判定結果に基づいて最過酷事故ケースを選定することによって、検討すべき事故条件を絞り込み、検討時間を短縮することができる。

なお、上記実施形態は、本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明は、その趣旨を逸脱することなく、変更、改良され得るとともに、本発明にはその等価物も含まれる。

１０電力系統シミュレータ
１００コンピュータ
１０１事故ケース選択部
１０２初期値設定部
１０３潮流計算部
１０４安定度判定部
１０５最過酷事故選定部
１０６入力部
１０７出力部
１０８記憶部
１０９バス

Claims

電力系統における複数の事故条件をシミュレートして、前記複数の事故条件のうち、前記電力系統に安定的に流し得る潮流量の安定限界が最小となる最過酷事故の事故条件を選定する電力系統シミュレータであって、
前記複数の事故条件から１つをシミュレーション対象として順次選択する条件選択部と、
前記シミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定する初期値設定部と、
前記シミュレーション対象の事故条件に対して、潮流量をパラメータとして潮流計算を行う潮流計算部と、
前記潮流計算部の計算結果に基づいて、前記電力系統の過渡安定度の判定を行う安定度判定部と、
前記複数の事故条件のそれぞれに対する前記安定度判定部の判定結果に基づいて、前記最過酷事故の事故条件を選定する最過酷事故選定部と、
を有し、
前記初期値設定部は、現在のシミュレーション対象の事故条件に対する前記安定度判定部の判定結果に基づいて、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定することを特徴とする電力系統シミュレータ。
請求項１に記載の電力系統シミュレータであって、
前記条件選択部によって最初に前記シミュレーション対象として選択された事故条件に対しては、
前記安定度判定部は、前記電力系統の過渡安定度の判定を行って前記安定限界を求め、
前記初期値設定部は、前記安定度判定部によって不安定と判定された最小の潮流量を、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値として設定し、
前記条件選択部によって２番目以降に前記シミュレーション対象として選択された事故条件に対しては、
前記安定度判定部は、前記初期値設定部によって設定された潮流量の初期値に対する前記電力系統の過渡安定度の判定を行い、不安定と判定した場合には、順次減少させた潮流量に対する前記電力系統の過渡安定度の判定を行って前記安定限界を求め、
前記初期値設定部は、前記安定度判定部が前記安定限界を求めた場合には、前記安定度判定部によって不安定と判定された最小の潮流量を、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値として設定し、
前記最過酷事故選定部は、前記安定度判定部が最後に前記安定限界を求めた事故条件を前記最過酷事故の事故条件として選定することを特徴とする電力系統シミュレータ。
請求項１または請求項２に記載の電力系統シミュレータであって、
前記電力系統の状態を示す系統情報と前記複数の事故条件とが記憶されている記憶部をさらに有することを特徴とする電力系統シミュレータ。
コンピュータに、
電力系統における複数の事故条件から１つをシミュレーション対象として順次選択する条件選択処理と、
前記シミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定する初期値設定処理と、
前記シミュレーション対象の事故条件に対して、潮流量をパラメータとして潮流計算を行う潮流計算処理と、
前記潮流計算処理の計算結果に基づいて、前記電力系統の過渡安定度の判定を行う安定度判定処理と、
前記複数の事故条件のそれぞれに対する前記安定度判定処理の判定結果に基づいて、前記複数の事故条件のうち、前記電力系統に安定的に流し得る潮流量の安定限界が最小となる最過酷事故の事故条件を選定する最過酷事故選定処理と、
を実行させ、
前記初期値設定処理は、現在のシミュレーション対象の事故条件に対する前記安定度判定処理の判定結果に基づいて、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定することを特徴とする電力系統のシミュレーションプログラム。
請求項４に記載の電力系統のシミュレーションプログラムであって、
前記条件選択処理によって最初に前記シミュレーション対象として選択された事故条件に対しては、
前記安定度判定処理は、前記電力系統の過渡安定度の判定を行って前記安定限界を求め、
前記初期値設定処理は、前記安定度判定処理によって不安定と判定された最小の潮流量を、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値として設定し、
前記条件選択処理によって２番目以降に前記シミュレーション対象として選択された事故条件に対しては、
前記安定度判定処理は、前記初期値設定処理によって設定された潮流量の初期値に対する前記電力系統の過渡安定度の判定を行い、不安定と判定した場合には、順次減少させた潮流量に対する前記電力系統の過渡安定度の判定を行って前記安定限界を求め、
前記初期値設定処理は、前記安定度判定処理が前記安定限界を求めた場合には、前記安定度判定処理によって不安定と判定された最小の潮流量を、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値として設定し、
前記最過酷事故選定処理は、前記安定度判定処理が最後に前記安定限界を求めた事故条件を前記最過酷事故の事故条件として選定することを特徴とする電力系統のシミュレーションプログラム。
電力系統における複数の事故条件をシミュレートして、前記複数の事故条件のうち、前記電力系統に安定的に流し得る潮流量の安定限界が最小となる最過酷事故の事故条件を選定する電力系統のシミュレーション方法であって、
前記複数の事故条件から１つをシミュレーション対象として順次選択する条件選択ステップと、
前記シミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定する初期値設定ステップと、
前記シミュレーション対象の事故条件に対して、潮流量をパラメータとして潮流計算を行う潮流計算ステップと、
前記潮流計算ステップの計算結果に基づいて、前記電力系統の過渡安定度の判定を行う安定度判定ステップと、
前記複数の事故条件のそれぞれに対する前記安定度判定ステップの判定結果に基づいて、前記最過酷事故の事故条件を選定する最過酷事故選定ステップと、
を有し、
前記初期値設定ステップは、現在のシミュレーション対象の事故条件に対する前記安定度判定ステップの判定結果に基づいて、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値を設定することを特徴とする電力系統のシミュレーション方法。
請求項６に記載の電力系統のシミュレーション方法であって、
前記条件選択ステップによって最初に前記シミュレーション対象として選択された事故条件に対しては、
前記安定度判定ステップは、前記電力系統の過渡安定度の判定を行って前記安定限界を求め、
前記初期値設定ステップは、前記安定度判定ステップによって不安定と判定された最小の潮流量を、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値として設定し、
前記条件選択ステップによって２番目以降に前記シミュレーション対象として選択された事故条件に対しては、
前記安定度判定ステップは、前記初期値設定ステップによって設定された潮流量の初期値に対する前記電力系統の過渡安定度の判定を行い、不安定と判定した場合には、順次減少させた潮流量に対する前記電力系統の過渡安定度の判定を行って前記安定限界を求め、
前記初期値設定ステップは、前記安定度判定ステップが前記安定限界を求めた場合には、前記安定度判定ステップによって不安定と判定された最小の潮流量を、次のシミュレーション対象の事故条件をシミュレートする際の潮流量の初期値として設定し、
前記最過酷事故選定ステップは、前記安定度判定ステップが最後に前記安定限界を求めた事故条件を前記最過酷事故の事故条件として選定することを特徴とする電力系統のシミュレーション方法。