JP2006203985A - 下位系負荷モデルの定数算出方法及び算出プログラム - Google Patents

下位系負荷モデルの定数算出方法及び算出プログラム Download PDF

Info

Publication number
JP2006203985A
JP2006203985A JP2005011167A JP2005011167A JP2006203985A JP 2006203985 A JP2006203985 A JP 2006203985A JP 2005011167 A JP2005011167 A JP 2005011167A JP 2005011167 A JP2005011167 A JP 2005011167A JP 2006203985 A JP2006203985 A JP 2006203985A
Authority
JP
Japan
Prior art keywords
constant
value
positive
observation
negative
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2005011167A
Other languages
English (en)
Other versions
JP4370263B2 (ja
Inventor
Shintaro Komami
慎太郎 駒見
Norikazu Kanao
則一 金尾
Tomoyuki Ueda
智之 上田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hokuriku Electric Power Co
Original Assignee
Hokuriku Electric Power Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hokuriku Electric Power Co filed Critical Hokuriku Electric Power Co
Priority to JP2005011167A priority Critical patent/JP4370263B2/ja
Publication of JP2006203985A publication Critical patent/JP2006203985A/ja
Application granted granted Critical
Publication of JP4370263B2 publication Critical patent/JP4370263B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

【課題】 電力下位系統の負荷の物理的イメージの保存、パラメータの数の少ない負荷モデルの案出、並びにそのパラメータの算出。
【解決手段】 下位系統負荷を、定リアクタンスXLと可変コンダクタンスGLの直列回路と、定アドミタンスBCを並列接続してあるモデルと仮定。瞬時電圧低下時の観測データを利用し、観測・解析データと所定式によりBCを変数とするデータ適用式でXLを表し、代入値を利用してデータ適用式を解き、求めたXLを1次近似して1次関数の傾きの正負を判別。データ適用式の代入値を一定値ごと傾きが逆転する方向に増減させその度にデータ適用式を解き傾きを求める一連のルーチンを、傾きが正負逆転するまで繰り返して次関数の傾きを0に近づける。傾きが逆転した時点の代入値をBCとする。
【選択図】 図1

Description

本発明は、電力系統の下位系統の負荷を簡易モデルで模擬した場合に、そのモデルで使用する定数を算出する方法及びプログラムに関する。ここで下位系統とは、特別高圧(77、66kV)以下の系統を意味する。
電力系統において基幹系統の安定度解析を実施する場合には、発電機、送電線、および負荷などの模擬が安定度解析の結果に大きな影響を及ぼすので、これらは忠実にモデル化することが重要である。
従来、発電機、発電機制御、および系統制御に関するモデル化に大きな努力が払われてきた。そして、負荷のモデル化については、その重要性は認識されており、特に、電圧変動に対する負荷特性の重要性は古くから指摘されている。しかしながら、その実態を把握することの難しさから、十分な精度がある負荷モデルは現在まで存在しない。実態を把握することの難しさとは、配電線等のデータ量が膨大なことや、下位系統の詳細なデータを収集することが実際無理であることや、仮に下位系統の詳細なデータがあっても系統構成が変更されることがあること等を意味する。
この問題を解決するために送電線落雷時等の負荷応動を自動記録し、その分析から特性を把握する調査が進められている。また、上記記録データから統計処理による1次進み遅れ負荷動特性モデルの定数算定が下記式(1)(2)のように行われている(非特許文献1)。
Figure 2006203985
Figure 2006203985
P:負荷の有効電力、P0:負荷の無効電力初期値、V:負荷端の電圧、V0:負荷端の電圧初期値、α:有効電力電圧特性指数、αs:有効電力瞬時電圧特性指数、Tα:有効電力負荷時定数、Q:負荷の無効電力、Q0:負荷の無効電力初期値、β:無効電力電圧特性指数、βs:無効電力瞬時電圧特性指数、Tβ:無効電力負荷時定数
浅田、長尾他:「大幅電圧変動時の負荷特性」電気学会 電力技術研究会、1986年6月、PE-86-32
しかし、上述した従来の負荷モデルは、有効電力と無効電力を同一の形で扱っているので、物理的な対応が明らかでなく、それ故、解析結果における負荷応動の妥当性や現象解明が困難となる場合もある。また、モデルのパラメータの数が多すぎたり、不適切なパラメータ設定によって過度にその値がばらついたりするという問題点がある。
そこで、事故クリア直後の負荷の動特性は、主に誘導電動機によるものと考え、誘導電動機の数式モデルから導かれる簡易な動的負荷モデルを静的負荷モデルと組み合わせる方法が提案された(非特許文献2)。
浅田 他、「わが国の60Hz系統における負荷特性の実態分析」、電力中央研究所報告、平成12年4月、47-52P
この方法では図8に示すように、可変コンダクタンスGと可変サセプタンスBLと定コンデンサB1を並列接続したモデルとして負荷を考えている。しかし、このモデルを実系統で適用した場合、有効電力や無効電力の応答の傾向はあっているものの、実測値と一致しない(精度が悪い)という問題点があった。
これは、送配電線のインピーダンスZと、系統に存在するアドミタンスBoと、負荷のアドミタンスB1を、観測データに基づく観測・解析データ(観測点における有効電力P、無効電力Q、電圧V)や設備データから正確に推定することが困難であるところに原因がある。
また、負荷の一般的なイメージとして図2のように考えられることを想定すると、負荷損と漏れリアクタンスからなるインピーダンスを無視することは物理的なイメージからかけ離れることになる。
本発明は上記実情を考慮して創作されたもので、その目的とするところは、負荷の物理的イメージを保存することができると共に、パラメータの数をできる限り少なくできる負荷モデルを案出し、そのパラメータを、観測データから定数として算出することである。
請求項1の発明は、電力系統の観測点から下位系統の負荷を簡略化して、定リアクタンスXLと可変コンダクタンスGLからなる直列回路と、定アドミタンスBCを、並列接続してあるモデルと仮定した場合の下位系負荷モデルの定数算出方法である。
この負荷モデルは図1に示すように、下位系統の負荷そのもの(系統側のインピーダンス(送配電線)を除いた分)を、定リアクタンスXL背後に可変コンダクタンスGLを直列接続してあるものと考えることを前提とする。電気機器は電気エネルギー即ち有効電力を他のエネルギーに変換するが、変換場面は電気回路側から見れば、負荷の一般的イメージを表した図2に示すようにコンダクタンス(純抵抗)となる。このコンダクタンスは電気機器の負荷条件によって変化する。無効電力消費は直列リアクタンス(漏れリアクタンス)と並列リアクタンス(励磁リアクタンス)で行われる。これらのリアクタンスは一定である。電気機器の無効電力消費を補償するためには並列コンデンサ(補償コンデンサ)が用いられている。補償コンデンサの作用は励磁リアクタンスとの差し引きで考慮すればよい。励磁リアクタンスと補償コンデンサは電力系統から見れば単なる定アドミタンスである。ここで負荷の動特性を大きく作用するのは漏れリアクタンスと可変コンダクタンスである。コンダクタンスが可変ゆえに漏れリアクタンスの消費する無効電力も変化するからである。
上記理由を考慮すると観測点からの負荷は図3のように、送配電線インピーダンスZと、漏れリアクタンスXiと、可変コンダクタンスGを直列接続したものに、アドミタンスB(観測点以下のアドミタンス(励磁リアクタンスと補償コンデンサ)の差し引き合計を観測点に集約する。)を並列接続したものと見える。
そして、図3を等価モデルで表すと図4のように、定リアクタンスXLと可変コンダクタンスGLの直列回路に対して、定アドミタンスBCを並列接続したものとなる。ここで送配電線に含まれるR分は、可変コンダクタンスに含まれ、送配電線に含まれるX分は定リアクタンスに含まれるものとする。
そして、観測点で採取した観測データより、観測点における電圧V、観測点から負荷へ送電する有効電力P、無効電力Qが得られる。詳細に言えば、観測点(一次変電所または連系用変電所)に設置した観測装置によって、観測点における瞬時電圧低下時の電圧・電流の時系列データを採取する。そして、観測装置とネットワークで繋がれた観測データ解析装置では、観測装置から送られた電圧・電流の時系列データを利用して、有効電力、無効電力の時系列の解析データを計算して求める。
また、定リアクタンスXLと可変コンダクタンスGLとの直列回路における無効電力成分QL(t)、可変コンダクタンスGL、定リアクタンスXLは、それぞれ下記式(3)〜(5)によって表される。即ち、直列回路の無効電力成分QL(t)は観測点における無効電力Q(t)と定アドミタンスの無効電力成分から下記式(3)で表される。
Figure 2006203985
また、可変コンダクタンスは、観測点における有効電力P(t)、直列回路の無効電力成分QL(t)、観測点における電圧V(t)とから下記式(4)で表される。
Figure 2006203985
さらに、観測点における有効電力P(t)、直列回路の無効電力成分QL(t)、及び観測点における電圧V(t)とから定リアクタンスXLは下記式(5)で表される。
Figure 2006203985
これら式(3)〜(5)に従い、観測点から負荷側の定アドミタンスBC、定リアクタンスXL、及び可変コンダクタンスGLを計算できる。
まず、式(5)に式(3)を代入すれば、XL(t)はQ(t)、V(t)、P(t)、BCで表される。BC以外は観測・解析データより把握してあるので、XL(t)はBCを変数とする時系列の関数式(観測・解析データ適用式)で表され、定アドミタンスBCの大きさにより事故除去後(電圧回復後)の定リアクタンスXLは図6のように、右上がりや、一定、或いは右下がりとなる。ここで、送配電線のリアクタンスと負荷の漏れリアクタンスが一定であることから定リアクタンスXLは事故除去後一定であり、傾きがほぼ0の1次式、XL(t)=at+bで表されると考えられる。
従って、定アドミタンスBCを決定するには以下のようにする。観測・解析データ適用式の定アドミタンスBCに、適当な初期値(通常0から1の値)BCを代入して解き、その代入値から得られる定リアクタンスを1次近似して1次関数の傾きの正負を判別する正負判別ステップをまず行なう。その後に1次関数の傾きを0に近づける正負逆転ステップを行う。正負逆転ステップは、観測・解析データ適用式の代入値を、直前のステップで得られた傾きの正負が逆転する方向に一定値ごと増減させ、その度に観測・解析データ適用式を解き、1次近似して1次関数の傾き及び定常値を求める一連のルーチンを、傾きの正負が逆転するまで繰り返すものである。正負逆転ステップを連続して複数回行なう場合には二回目以降の正負逆転ステップでは、一連のルーチンで用いる一定値を、直前の正負逆転ステップで用いた値よりも絶対値で小さくすることによって1次関数の傾きを0に近づけるものとする。最後の正負逆転ステップの終了時点又は終了直前の代入値及び定常値を、それぞれ定アドミタンスBCの値及び定リアクタンスXLの値とする決定ステップを最後に行なう。
この後、可変コンダクタンスGLを求めるステップを行なっても良い。この場合は、電圧V、有効電力P、定アドミタンスBC、式(3)(4)に従って、可変コンダクタンスGLを求める。
請求項2の発明は、請求項1の発明をコンピュータで実行するためのプログラムである。即ち、電力系統の観測点から下位系統の負荷を簡略化して上述したようにモデル化し、観測点における瞬時電圧低下時の、電圧V・有効電力P・無効電力Qの時系列の観測・解析データを利用して負荷モデルの定数を求めるためにコンピュータに以下のステップを実行させる下位系負荷モデルの定数算出プログラムである。実行させるステップは、観測・解析データ適用式の定アドミタンスBCに代入する値を入力する入力フォームを出力装置に表示するステップと、入力フォームに入力された代入値を観測・解析データ適用式に代入して定リアクタンスXLを求め、求めた定リアクタンスXLを1次近似して1次関数の傾きの正負を判別する正負判別ステップと、1次関数の傾きを0に近づける正負逆転ステップと、定アドミタンスBCの値と定リアクタンスXLを決定する決定ステップである。
請求項1の発明は、定リアクタンスの値が事故除去後一定となると考えられることに基づいて、1次関数の傾きを0に近づけて定リアクタンスの値を求めているので、負荷の物理的イメージが保存される。また、定リアクタンスXLと可変コンダクタンスGLからなる直列回路と、定アドミタンスBCを、並列接続してあるモデルとして、下位系統の負荷を考えられるので、負荷のパラメータも少なく、そのパラメータである定数を求めやすい。さらに、観測・解析データに基づいて定アドミタンスBCの値を算出しているため、詳細な下位系統のデータがない場合、配電線などのデータ量が膨大で算出が困難な場合でも観測点から下位系統の負荷に関する定アドミタンスと定リアクタンスの値を求めることができる。
請求項2の発明のプログラムを実行すれば、請求項1の発明と同等の効果が得られる。
本発明の下位系負荷モデルの定数算出方法は、図4に示すように一次変電所または連系用変電所に設置した観測装置で、瞬時電圧低下時の観測データを自動的に採取し、その観測データを別地点の観測データ解析装置に送信し、観測データ解析装置で観測データを解析し、その観測・解析データを定数算出装置(下位系負荷モデルの定数算出プログラムがインストールされたコンピュータ)に送信して、定数算出装置で各種の定数を算出する。なお、算出された各種定数、観測・解析データ、並びに設備データより求めた電力系統簡易モデルの上位系インピーダンスは、電力系統安定度解析プログラムを実行する際に使用する。
コンピュータで定数算出プログラムを実行すると、図示しない各種手段が構築されて、以下のステップが図5のフローチャートに示すように順番に行われる。まず、入力フォーム表示手段によって、まず観測・解析データ適用式のBCをBiと置換し、観測・解析データ適用式に定アドミタンスの数値をBiとして代入するための入力フォームが記憶装置から読み込まれて出力装置に表示される。入力フォームに初期値としての一般的な値(0から1の値、例えば0)を作業者が入力装置から入力し、値を確定する操作を入力装置にすると、正負判別手段によって、入力フォームに入力された代入値を観測・解析データ適用式に入れて演算し定リアクタンスXLを求める。続いて、正負判別手段が、求めた定リアクタンスXLを最小二乗法によって1次近似し、その近似した1次関数の傾きa、定常値bを算出し、その後、傾きが負か否かを判別する。
負の場合は図の左側の正負逆転ステップに移行する。一番目の正負逆転ステップでは、正負逆転手段によってBiの代入値を一定値(0.1)増やし、正負判別ステップと同様の処理で1次関数の傾きa、定常値bを求め、傾きaが正か否かを判定し、判定結果が正になるまで、代入値を0.1増やし、同様にして傾きaが正か否かを判定する一連のルーチンを行なう。この一連のルーチンは、判定結果が正になるまで繰り返す。判定結果が正になった場合は二番目の正負逆転ステップに移行する。二番目の正負逆転ステップでは、一番目の正負逆転ステップと同様に行なうが、Biの代入値を0.01減らすこと、一連のルーチンを判定結果が負になるまで繰り返すことが異なる。判定結果が負になった場合は三番目の正負逆転ステップに移行する。三番目の正負逆転ステップでは二番目の正負逆転ステップと同様に行なうが、Biの代入値を0.001増やすこと、一連のルーチンを判定結果が正になるまで繰り返すことが異なる。判定結果が正になった場合は、決定ステップに移行する。決定ステップは、三番目の正負逆転ステップで判定結果が正になった場合のBiの値を定アドミタンスBCの値として決定(置換)すると共に、最後に求めた定常値bを定リアクタンスXLの値として決定(置換)し、これらの結果を出力装置に表示すると共に、記憶装置に保存し、プログラムを終了する。
一方、最初の正負判定ステップで傾きaが負でなかった場合は、図の右側の正負逆転ステップを同様にして三回行い、定アドミタンスBCを求める。また、最後に求めた定常値bが定リアクタンスXLとなる。
このようにして、正負逆転ステップを何回も連続して行う場合には、回を追うごとに、Biを増減する一定値を、直前のステップで用いた値よりも絶対値で小さくすることによって、1次関数の傾きaが徐々に0に近づく。なお、厳密に言えば、1次関数の傾きを判別したときに傾きが0の場合もありうるが、確率的には非常に低く、また、傾きが0の場合は次の正負逆転ステップによって、傾きが0から離れることになるが、一定値が小数点3桁(0.001)になるまでして定アドミタンスBCを求めているので、精度的には問題がないと考えられる。
なお、図5のフローチャートには記載していないが、最後の正負逆転ステップの後に、可変コンダクタンスGLを求めるステップを行なっても良い。この場合は、前述した式(3)を式(4)に代入して一纏めにすることによって、可変コンダクタンスGLが電圧V、有効電力P、定アドミタンスBCを変数とする関数で表される。従って、可変コンダクタンス算出手段がこの一纏めにした関数式に電圧V、有効電力P、定アドミタンスBCを代入することによって、可変コンダクタンスGLが算出され、これらの結果を出力装置に表示すると共に、記憶装置に保存する。ちなみに、可変コンダクタンスGLは図7に示すような応答となる。すなわち、瞬時電圧低下中、可変コンダクタンスは増大し、瞬低回復後、可変コンダクタンスは減少し定常値に落ち着く。
本発明は上記実施形態に限定されるものではない。例えば、三番目の正負逆転ステップでは、判定結果が正になった場合は正になる直前のBiの値をBCの値として決定すると共に、同様に正になる直前に求めた定常値bを定リアクタンスXLの値として決定しても良く、精度的には大差がないと考えられる。また、正負逆転ステップを1回だけ行って定アドミタンス等を求めても良いが、この場合、一定値には0.001以下の値(小数点3桁)を入れてやれば良く、但し計算に時間がかかる。
最終的な下位系統の負荷モデルを示す説明図である。 負荷の一般的イメージを示す説明図である。 一般的イメージを簡易化した負荷モデルを示す説明図である。 定数を算出する過程等を示す説明図である。 定数を算出する詳細なフローチャートである。 定リアクタンスXLの振る舞いを示すグラフである。 XL決定後算出される可変コンダクタンスGLの振る舞いを示すグラフである。 従来の負荷モデルを示す説明図である。
符号の説明
XL 定リアクタンス
GL 可変コンダクタンス
BC 定アドミタンス

Claims (2)

  1. 電力系統の観測点から下位系統の負荷を簡略化して、定リアクタンス(XL)と可変コンダクタンス(GL)からなる直列回路と、定アドミタンス(BC)を、並列接続してあるモデルと仮定し、
    観測点における無効電力Q(t)と定アドミタンスの無効電力成分から直列回路の無効電力成分QL(t)を下記式(3)で表し、
    Figure 2006203985
    観測点における有効電力P(t)、直列回路の無効電力成分Q(t)、及び観測点における電圧V (t)とから定リアクタンスを下記式(5)で表し、
    Figure 2006203985
    観測点における瞬時電圧低下時の、電圧(V)・有効電力(P)・無効電力(Q)の時系列の観測・解析データを利用し、事故除去後の観測・解析データと式(3)と式(5)により定アドミタンス(BC)を変数とする観測・解析データ適用式で定リアクタンス(XL)を表し、
    観測・解析データ適用式の定アドミタンス(BC)に仮の値を代入して解き、求めた定リアクタンス(XL)を1次近似して1次関数の傾きの正負を判別する正負判別ステップを行い、
    その後に1次関数の傾きを0に近づける正負逆転ステップを行い、
    正負逆転ステップは、観測・解析データ適用式の代入値を、直前のステップで得られた傾きの正負が逆転する方向に一定値ごと増減させ、その度に観測・解析データ適用式を解き、1次近似して1次関数の傾き及び定常値を求める一連のルーチンを、傾きの正負が逆転するまで繰り返すものであって、
    正負逆転ステップを連続して複数回行なう場合には二回目以降の正負逆転ステップでは、一連のルーチンで用いる一定値を、直前の正負逆転ステップで用いた値よりも絶対値で小さくすることによって1次関数の傾きを0に近づけ、
    最後の正負逆転ステップの終了時点又は終了直前の代入値及び定常値を、それぞれ定アドミタンス(BC)の値及び定リアクタンス(XL)の値とする決定ステップを行なうことを特徴とする下位系負荷モデルの定数算出方法。
  2. 電力系統の観測点から下位系統の負荷を簡略化してモデル化し、観測点における瞬時電圧低下時の、電圧(V)・有効電力(P)・無効電力(Q)の時系列の観測・解析データを利用して負荷モデルの定数を求めるためにコンピュータに以下のステップを実行させる下位系負荷モデルの定数算出プログラムであって、
    負荷モデルは、定リアクタンス(XL)と可変コンダクタンス(GL)からなる直列回路と、定アドミタンス(BC)を、並列接続すると共に、観測点における無効電力Q(t)と定アドミタンスの無効電力成分から直列回路の無効電力成分QL(t)を下記式(3)で表し、
    Figure 2006203985
    観測点における有効電力P(t)、直列回路の無効電力成分QL(t)、及び観測点における電圧V
    (t)とから定リアクタンスを下記式(5)で表し、
    Figure 2006203985
    事故除去後の観測・解析データと式(3)と式(5)により定アドミタンス(BC)を変数とする観測・解析データ適用式で定リアクタンス(XL)を表すモデルであって、
    実行させるステップは、
    観測・解析データ適用式の定アドミタンス(BC)に代入する値を入力する入力フォームを出力装置に表示するステップと、
    入力フォームに入力された代入値を観測・解析データ適用式に代入して定リアクタンス(XL)を求め、求めた定リアクタンス(XL)を1次近似して1次関数の傾きの正負を判別する正負判別ステップと、
    1次関数の傾きを0に近づける正負逆転ステップと、
    定アドミタンス(BC)の値を決定する決定ステップであって、
    正負逆転ステップは、観測・解析データ適用式の代入値を、直前のステップで得られた傾きの正負が逆転する方向に一定値ごと増減させ、その度に観測・解析データ適用式を解き、1次近似して1次関数の傾き及び定常値を求める一連のルーチンを、傾きの正負が逆転するまで繰り返すものであって、
    正負逆転ステップを連続して複数回行なう場合には、二回目以降の正負逆転ステップでは、一連のルーチンで用いる一定値を、直前の正負逆転ステップで用いた値よりも絶対値で小さくすることによって1次関数の傾きを0に近づけるものとし、
    決定ステップは、最後の正負逆転ステップの終了時点又は終了直前の代入値及び定常値を、それぞれ定アドミタンス(BC)の値及び定リアクタンス(XL)の値とするものである、ことを特徴とする下位系負荷モデルの定数算出プログラム。
JP2005011167A 2005-01-19 2005-01-19 下位系負荷モデルの定数算出方法及び算出プログラム Expired - Fee Related JP4370263B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005011167A JP4370263B2 (ja) 2005-01-19 2005-01-19 下位系負荷モデルの定数算出方法及び算出プログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005011167A JP4370263B2 (ja) 2005-01-19 2005-01-19 下位系負荷モデルの定数算出方法及び算出プログラム

Publications (2)

Publication Number Publication Date
JP2006203985A true JP2006203985A (ja) 2006-08-03
JP4370263B2 JP4370263B2 (ja) 2009-11-25

Family

ID=36961497

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005011167A Expired - Fee Related JP4370263B2 (ja) 2005-01-19 2005-01-19 下位系負荷モデルの定数算出方法及び算出プログラム

Country Status (1)

Country Link
JP (1) JP4370263B2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007221943A (ja) * 2006-02-17 2007-08-30 Hokuriku Electric Power Co Inc:The 電力系統における電圧安定性評価方法および装置ならびにプログラム
JP2008131829A (ja) * 2006-11-24 2008-06-05 Hokuriku Electric Power Co Inc:The 下位系統の負荷モデルの定数算出方法及びプログラム

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007221943A (ja) * 2006-02-17 2007-08-30 Hokuriku Electric Power Co Inc:The 電力系統における電圧安定性評価方法および装置ならびにプログラム
JP4584157B2 (ja) * 2006-02-17 2010-11-17 北陸電力株式会社 電力系統における電圧安定性評価方法および装置ならびにプログラム
JP2008131829A (ja) * 2006-11-24 2008-06-05 Hokuriku Electric Power Co Inc:The 下位系統の負荷モデルの定数算出方法及びプログラム

Also Published As

Publication number Publication date
JP4370263B2 (ja) 2009-11-25

Similar Documents

Publication Publication Date Title
WO2022080377A1 (ja) リチウムイオンバッテリーの寿命予測方法、放電容量維持率予測方法、寿命予測プログラム、放電容量維持率予測プログラム及び情報処理装置
CN107918103B (zh) 一种基于灰色粒子滤波的锂离子电池剩余寿命预测方法
CN110212592B (zh) 基于分段线性表达的火电机组负荷调节最大速率估计方法及系统
CN112288147B (zh) 一种BP-Adaboost强预测器预测发电机定子绝缘状态的方法
CN116756530B (zh) 用于新能源接入配电网的电网线损评估方法及系统
JP2006340552A (ja) 電力品質評価システムとその方法、およびプログラム
CN112816881A (zh) 电池压差异常检测方法、装置及计算机存储介质
CN108961094A (zh) 基于输电线路最小空气间隙在线量测的风偏故障预警方法
CN107834551A (zh) 一种基于支持向量机的配电网低电压预测方法
CN109378835A (zh) 基于互信息冗余性最优的大规模电力系统暂态稳定评估系统
CN104680017B (zh) 时变电力系统稳定性分析系统及方法
CN113011012A (zh) 一种基于Box-Cox变化的储能电池剩余寿命预测方法
CN116912589A (zh) 一种基于多模态融合的中压配电电缆状态评估方法
CN114970154A (zh) 一种暂态电压支撑能力量化评估指标构建方法
CN113904384B (zh) 一种基于梯度提升机的电网暂态稳定协调控制方法及系统
JP4370263B2 (ja) 下位系負荷モデルの定数算出方法及び算出プログラム
CN105514975B (zh) 一种光伏发电系统的能效预测方法
CN113300373A (zh) 一种基于prmse评价指标的稳定裕度值预测方法及装置
CN117394436A (zh) 基于随机场景分析的光伏承载能力评估方法及系统
CN110098610B (zh) 故障扰动下电力系统振荡主导模式的实时辨识方法及系统
CN108123436B (zh) 基于主成分分析和多元回归算法的电压越限预测模型
Kontis et al. Development of measurement-based load models for the dynamic simulation of distribution grids
CN116400266A (zh) 基于数字孪生模型的变压器故障检测方法、装置及介质
CN110943473A (zh) 一种基于广域测量系统和聚类理论的发电机同调识别方法
CN108920810A (zh) 基于大数据架构的电力仿真方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20071203

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090507

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090519

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090708

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090804

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090831

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120904

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130904

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees