JP5159695B2 - 配電系統状態推定方法及び配電系統状態推定装置 - Google Patents
配電系統状態推定方法及び配電系統状態推定装置 Download PDFInfo
- Publication number
- JP5159695B2 JP5159695B2 JP2009114882A JP2009114882A JP5159695B2 JP 5159695 B2 JP5159695 B2 JP 5159695B2 JP 2009114882 A JP2009114882 A JP 2009114882A JP 2009114882 A JP2009114882 A JP 2009114882A JP 5159695 B2 JP5159695 B2 JP 5159695B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- node
- power
- distribution system
- calculation
- 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.)
- Expired - Fee Related
Links
Description
まず、特許文献1の手法の場合、潮流計算を用いているため、計算時間が膨大となり、オンラインで高速に状態推定を行うのには不向きである。また、負荷の力率を一定と仮定しているため、状態推定の誤差が大きくなるおそれもある。
この方法によれば、ノードごとの偏差及び偏差の合計という2つの観点から、各ノードの利用率及び力率の妥当性を評価することができる。また、所定の係数を乗算することにより、評価値に対する各項の重み付けを調整することができる。
図1は、配電系統状態推定方法を実現する配電系統状態推定システム10の構成例を示す図である。
配電系統の状態推定手法における最適化問題を定式化した目的関数を、次の式1a及び式1bに示す。そして、本最適化問題の制約条件として、式2aによる力率cosθの制約条件及び式2bによる負荷7の利用率ruの制約条件が定義される。ここで、i(=1〜n)は、ノードの番号であり、電圧、通過有効電力及び通過無効電力の測定値が取得可能なノードのみに関して各値を有し、それ以外のノードの各値はゼロとする。すなわち、式1におけるΣの計算は、電圧、通過有効電力及び通過無効電力が測定可能なノードについての和を示す。なお、2乗されている計算項は、絶対値に置き換えてもよい。
図4は、配電系統の状態推定の処理例を示すフローチャートである。本処理は、配電系統状態推定システム10の情報処理装置100において、処理部101が、メモリ102又は記憶部104に記憶された設備データ及び通信部103により受信された計測データから、各ノードの有効電力、無効電力及び電圧を計算するものであり、最適化の手法としてメタヒューリスティクスの手法の一つであるタブサーチを適用する。以下に、その処理の概要を説明する。
続いて、各処理の詳細を説明する。
===ノード毎の負荷P、Qの算出===
各ノードの設備容量が予め分かっていることを前提条件として、各ノードの負荷の利用率及び力率を暫定的に設定することにより、設備容量、利用率及び力率から各ノードの有効電力及び無効電力を算出することができる。具体的には、以下の式3及び式4による。
各ノードの有効電力及び無効電力が求められた場合、ブランチの通過有効電力及び通過無効電力を算出することができる。具体的には、以下の式5及び式6に示すように、ブランチi−1〜i(ノードi)の通過電力P、Qは、ノードi〜ノードnの各ノードのP、Qをそれぞれ合計する。
===ノードiに流れ込む電流とノードiの電圧の算出===
ノードi−1〜i間の電圧降下を計算し、変電所送出し電圧から逐次引き算を行う。具体的には、以下の式7、式8及び式9による。
以上から、測定可能ノードに関して、測定不可ノードの利用率及び力率を変化させながら、式3〜式9を用いて算出した通過有効電力、通過無効電力及び電圧降下の値を式1の目的関数f(x)に代入して評価値を求め、図4のS405及びS406に示すように、探索処理ごとにそれまでの最高評価値と比較し、更新する。これにより、目的関数f(x)の評価値が最小になる測定不可ノードの利用率及び力率のパターンを特定する。そして、その利用率及び力率のパターンから、測定不可ノードの電圧、有効電力及び無効電力を計算し、出力する。
図5は、現在負荷から近傍負荷を求める方法について説明する図である。現在負荷である利用率rui及び力率pfi(iはノード番号)に対して、近傍負荷として利用率rui±aiΔru及び力率pfi±biΔpfを設定する。係数ai及びbiは、1、2、3、4のいずれかの数値であり、正負の符号は乱数により決定される。Δruは、利用率の変化単位であり、例えば、現在負荷の利用率ruiの0.25%が用いられる。Δpfは、力率の変化単位であり、例えば、現在負荷の力率pfiの0.2%が用いられる。
図4に示した配電系統の状態推定の処理において、現在状態(現在負荷)を遷移させる方法について説明する。
(1)予め設定した現在状態(各ノードの負荷の利用率及び力率)に対して100通りのパターン(近傍負荷)を作成する。
(2)作成した100通りのパターンについて、評価値(目的関数の計算結果)を算出する。
(3)最も評価の高い(すなわち、評価値の小さい)パターンを現在状態に置き換える。
(4)置き換えた現在状態を元にして、さらに100通りのパターンを作成し、(2)に戻る。
評価した全てのパターンに対して最高評価の値(評価値の最小値)及びパターン(利用率と力率)を記憶しておき、設定した回数の探索[(2)〜(4)の処理]を終了した際に解として算出する。探索ごとの最高評価の値のうち、さらに最も小さい評価値及びそのパターンを最適解として記憶するようにしてもよい。
次に、図8〜図15を参照しながら、上記の配電系統状態推定方法を適用した実施例について説明する。特に、上記の状態推定方法による最適化の結果(計算値)と、潮流計算の結果(一部を真値、すなわち、計測値として用いる)とを比較する。
図10は、潮流計算結果を計測値とする地点である3点を示す。ノード番号7、13、19の各ノードの電圧及びノードの区間6−7、12−13、18−19の各ブランチの通過有効電力と通過無効電力の潮流計算結果が計測値として示されている。
図11は、各ノードの電圧、有効電力、無効電力及び通過電流について、潮流計算による真値と、状態推定結果との比較を示す。なお、状態推定処理における探索回数は、1000回である。
図12は、最適解の評価値の収束状況、電圧誤差の発生確率及び通過電流誤差の発生確率を示す。探索回数は、100回以下で十分と考えられる。
図13は、潮流計算結果を計測値とする地点である4点を示す。ノード番号4、10、13、19の各ノードの電圧及びノードの区間3−4、9−10、12−13、18−19の各ブランチの通過有効電力と通過無効電力の潮流計算結果が計測値として示されている。
図14は、各ノードの電圧、有効電力、無効電力及び通過電流について、潮流計算による真値と、状態推定結果とを比較を示す。計測地点を4点とした方が、3点とした方より、真値と、状態推定結果との差が小さくなっている。なお、状態推定処理における探索回数は、200回である。
図15は、最適解の評価値の収束状況、電圧誤差の発生確率及び通過電流誤差の発生確率を示す。探索回数は、100回以下で十分と考えられる。
以上、本発明を実施するための形態について説明したが、上記実施の形態は本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明はその趣旨を逸脱することなく変更、改良され得るとともに、本発明にはその等価物も含まれる。例えば、上記実施の形態では、最適解探索手法としてタブサーチを示したが、式1に示す目的関数f(x)の評価値が最小になるような各ノードの利用率及び力率を特定できるのであれば他の手法でもよい。例えば、遺伝的アルゴリズム(Genetic Algorithm:GA)や焼きなまし法を用いるようにしてもよい。
2 計測機能付き開閉器
3 開閉器
4 受電点
5 SVR
6 連系点
7 負荷
8 通信装置
10 配電系統状態推定システム
100 情報処理装置(配電系統状態推定装置)
101 処理部
104 記憶部
CB 配電用遮断器
G 分散型電源
Claims (10)
- コンピュータにより、配電系統におけるノードの状態を推定する配電系統状態推定方法であって、
前記ノードは、電圧、通過有効電力及び通過無効電力を測定可能な測定可能ノードと、それらを測定できない測定不可ノードとを含み、
前記コンピュータは、
前記配電系統のうち、1の配電線に関する定数である配電線データ、当該配電線につながる各ノードの設備容量、利用率及び力率に基づいて、各ノードの電圧計算値、有効電力計算値及び無効電力計算値、通過有効電力計算値及び通過無効電力計算値を算出する計算式を記憶するステップと、
前記測定可能ノードの電圧、通過有効電力及び通過無効電力の測定値と計算値を入力パラメータとし、各測定値及び各計算値のうち、対応する測定値と、計算値との差に関する評価値を算出する目的関数を記憶するステップと、
前記配電線データ、前記設備容量、前記測定可能ノードの電圧測定値、通過有効電力測定値及び通過無効電力測定値を取得するステップと、
各測定不可ノードの利用率及び力率を変化させながら、前記計算式を用いて前記測定可能ノードの電圧計算値、有効電力計算値及び無効電力計算値を求め、前記目的関数を用いて前記評価値を求めることにより、当該評価値が最小になるような各測定不可ノードの利用率及び力率の組合せを特定する最適化ステップと、
特定した前記利用率及び力率の組合せに基づいて、前記計算式を用いて各測定不可ノードの電圧計算値、有効電力計算値及び無効電力計算値を求め、状態推定の結果として出力するステップと、
を実行することを特徴とする配電系統状態推定方法。 - 請求項1に記載の配電系統状態推定方法であって、
前記計算式は、前記設備容量と、前記利用率とを乗算することにより、前記有効電力計算値を算出し、当該有効電力計算値を加算することにより、前記通過有効電力計算値を算出する第1の計算式を含む
ことを特徴とする配電系統状態推定方法。 - 請求項2に記載の配電系統状態推定方法であって、
前記計算式は、前記力率からリアクタンス率を求め、前記有効電力計算値に前記リアクタンス率を乗算し、当該乗算値を前記力率で除算することにより、前記無効電力計算値を算出し、当該無効電力計算値を加算することにより前記通過無効電力計算値を算出する第2の計算式を含む
ことを特徴とする配電系統状態推定方法。 - 請求項3に記載の配電系統状態推定方法であって、
前記配電線データは、各ブランチの抵抗、リアクタンス及び亘長を含み、
前記計算式は、前記通過有効電力計算値、前記通過無効電力計算値及び前記配電線データから、前記電圧計算値を算出する第3の計算式を含む
ことを特徴とする配電系統状態推定方法。 - 請求項1に記載の配電系統状態推定方法であって、
前記目的関数は、各測定値及び各計算値のうち、対応する測定値と、計算値との差の絶対値の合計値又はその合計値に所定の係数を乗算した値と、当該差の合計値の絶対値又はその絶対値に所定の係数を乗算した値とを加算した値を前記評価値として算出する
ことを特徴とする配電系統状態推定方法。 - 請求項1に記載の配電系統状態推定方法であって、
前記最適化ステップは、タブサーチの手法を用いる
ことを特徴とする配電系統状態推定方法。 - 請求項6に記載の配電系統状態推定方法であって、
前記最適化ステップは、
(1)現在状態として各ノードの利用率及び力率を記憶し、当該現在状態から所定個数の近傍状態を設定し、記憶する。
(2)所定個数の各近傍状態について、前記目的関数を用いて前記評価値を求める。
(3)所定個数の各近傍状態のうち、前記評価値が最も小さい前記近傍状態に前記現在状態を更新し、当該近傍状態をタブーリストに追加する。
(4)当該近傍状態の評価値が最適解の評価値より小さい場合、当該近傍状態に前記最適解を更新する。
(5)(3)で更新した前記現在状態から所定個数の近傍状態を設定し、記憶し、(2)に戻る。
(6)所定の条件を満たすまで(2)〜(5)の探索処理を繰り返す。
の手順によることを特徴とする配電系統状態推定方法。 - 請求項7に記載の配電系統状態推定方法であって、
前記現在状態から所定個数の近傍状態を設定するには、
前記現在状態である各ノードの利用率ごとに、所定の変化単位を設定し、当該変化単位に1〜4のいずれかの数値を乗算し、その乗算値に乱数により決定した正負の符号を付与し、当該符号を付与した乗算値を前記利用率ごとに加算し、近傍状態とし、
前記現在状態である各ノードの力率ごとに、所定の変化単位を設定し、当該変化単位に1〜4のいずれかの数値を乗算し、その乗算値に乱数により決定した正負の符号を付与し、当該符号を付与した乗算値を前記力率ごとに加算し、近傍状態とする
ことを特徴とする配電系統状態推定方法。 - 請求項7に記載の配電系統状態推定方法であって、
前記所定の条件は、前記探索処理を所定の回数繰り返すこと又は前記評価値が所定の範囲内であることである
ことを特徴とする配電系統状態推定方法。 - 配電系統におけるノードの状態を推定する配電系統状態推定装置であって、
前記ノードは、電圧、通過有効電力及び通過無効電力を測定可能な測定可能ノードと、それらを測定できない測定不可ノードとを含み、
前記配電系統のうち、1の配電線に関する定数である配電線データ、当該配電線につながる各ノードの設備容量、利用率及び力率に基づいて、各ノードの電圧計算値、有効電力計算値及び無効電力計算値、通過有効電力計算値及び通過無効電力計算値を算出する計算式を記憶する手段と、
前記測定可能ノードの電圧、通過有効電力及び通過無効電力の測定値と計算値を入力パラメータとし、各測定値及び各計算値のうち、対応する測定値と、計算値との差に関する評価値を算出する目的関数を記憶する手段と、
前記配電線データ、前記設備容量、前記測定可能ノードの電圧測定値、通過有効電力測定値及び通過無効電力測定値を取得する手段と、
各測定不可ノードの利用率及び力率を変化させながら、前記計算式を用いて前記測定可能ノードの電圧計算値、有効電力計算値及び無効電力計算値を求め、前記目的関数を用いて前記評価値を求めることにより、当該評価値が最小になるような各測定不可ノードの利用率及び力率の組合せを特定する最適化手段と、
特定した前記利用率及び力率の組合せに基づいて、前記計算式を用いて各測定不可ノードの電圧計算値、有効電力計算値及び無効電力計算値を求め、状態推定の結果として出力する手段と、
を備えることを特徴とする配電系統状態推定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009114882A JP5159695B2 (ja) | 2009-05-11 | 2009-05-11 | 配電系統状態推定方法及び配電系統状態推定装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009114882A JP5159695B2 (ja) | 2009-05-11 | 2009-05-11 | 配電系統状態推定方法及び配電系統状態推定装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010263754A JP2010263754A (ja) | 2010-11-18 |
JP5159695B2 true JP5159695B2 (ja) | 2013-03-06 |
Family
ID=43361366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009114882A Expired - Fee Related JP5159695B2 (ja) | 2009-05-11 | 2009-05-11 | 配電系統状態推定方法及び配電系統状態推定装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5159695B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101824357B1 (ko) * | 2017-03-28 | 2018-01-31 | 엘에스산전 주식회사 | 무효전력 모델값 산출 방법 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6080349B2 (ja) | 2010-11-26 | 2017-02-15 | キヤノン株式会社 | 光学部材および撮像装置 |
US8648499B2 (en) * | 2011-01-27 | 2014-02-11 | General Electric Company | Systems, methods, and apparatus for accelerating volt/VAR load flow optimization |
KR101221083B1 (ko) * | 2011-05-25 | 2013-01-11 | 주식회사 파워이십일 | 배전 시스템의 상태추정 방법 |
JP5372125B2 (ja) * | 2011-12-08 | 2013-12-18 | 中国電力株式会社 | 配電系統制御装置、配電系統制御装置の制御方法及びプログラム |
JP5372124B2 (ja) * | 2011-12-08 | 2013-12-18 | 中国電力株式会社 | 配電系統管理支援装置、配電系統管理支援装置の制御方法及びプログラム |
JP5949232B2 (ja) * | 2012-07-05 | 2016-07-06 | 富士電機株式会社 | 電力系統監視システム、その推定値の異常判定方法、そのプログラム |
CN104022503B (zh) * | 2014-06-18 | 2017-01-04 | 中国科学院自动化研究所 | 一种带有储能设备的智能微电网电能优化控制方法 |
CN104834221A (zh) * | 2015-05-25 | 2015-08-12 | 中国科学院自动化研究所 | 一种基于可变误差的非线性系统自适应最优控制方法 |
JP6801457B2 (ja) * | 2017-01-04 | 2020-12-16 | 富士電機株式会社 | 状態推定装置、状態推定方法、及び状態推定プログラム |
CN110991896B (zh) * | 2019-12-04 | 2022-08-12 | 南京理工大学 | 主动配电网电压运行状态的点-区-网多颗粒度评估方法 |
CN110929900B (zh) * | 2019-12-17 | 2022-08-12 | 南京理工大学 | 一种电网智能检修管理系统及方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004064969A (ja) * | 2002-07-31 | 2004-02-26 | Kansai Electric Power Co Inc:The | 配電線の区間負荷算出方法、配電線の負荷監視制御システム及びプログラム |
JP4718943B2 (ja) * | 2005-09-15 | 2011-07-06 | 財団法人電力中央研究所 | 配電系統の負荷分布推定方法、装置及びプログラム、並びに電圧推定方法、装置及びプログラム |
JP4705563B2 (ja) * | 2006-12-20 | 2011-06-22 | 株式会社日立製作所 | 配電系統の状態推定装置、状態推定方法及びそのプログラム |
-
2009
- 2009-05-11 JP JP2009114882A patent/JP5159695B2/ja not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101824357B1 (ko) * | 2017-03-28 | 2018-01-31 | 엘에스산전 주식회사 | 무효전력 모델값 산출 방법 |
Also Published As
Publication number | Publication date |
---|---|
JP2010263754A (ja) | 2010-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5159695B2 (ja) | 配電系統状態推定方法及び配電系統状態推定装置 | |
Chen et al. | Robust capacity assessment of distributed generation in unbalanced distribution networks incorporating ANM techniques | |
Li et al. | Flexible transmission expansion planning associated with large‐scale wind farms integration considering demand response | |
Zhang et al. | Probabilistic evaluation of available load supply capability for distribution system | |
Zhang et al. | Probabilistic analysis for optimal power flow under uncertainty | |
Wang et al. | Decision tree based online stability assessment scheme for power systems with renewable generations | |
JP2008154418A (ja) | 配電系統の状態推定装置、状態推定方法及びそのプログラム | |
Chen et al. | Available transfer capability evaluation in a deregulated electricity market considering correlated wind power | |
El‐Araby et al. | Reactive power reserve management tool for voltage stability enhancement | |
US20150168465A1 (en) | Method and apparatus for electric power system distribution state estimations | |
Venzke et al. | Chance-constrained AC optimal power flow integrating HVDC lines and controllability | |
Hong et al. | Calculating line losses in smart grid: A new rule of thumb | |
Saberi et al. | A heuristic benders-decomposition-based algorithm for transient stability constrained optimal power flow | |
Hojjat et al. | Chance‐constrained programming approach to stochastic congestion management considering system uncertainties | |
Aigner et al. | Robust approximation of chance constrained DC optimal power flow under decision-dependent uncertainty | |
Madadi et al. | Integrated transmission expansion and PMU planning considering dynamic thermal rating in uncertain environment | |
Kumar et al. | Multi‐objective design of advanced power distribution networks using restricted‐population‐based multi‐objective seeker‐optimisation‐algorithm and fuzzy‐operator | |
US20140074305A1 (en) | Method and device for producing a state signal | |
Alberdi et al. | Overhead line ampacity forecasting and a methodology for assessing risk and line capacity utilization | |
Lyu et al. | Impacts of wind power integration on generation dispatch in power systems | |
Seta et al. | Distribution system planning with representation of uncertainties based on interval analysis | |
Menke et al. | Heuristic monitoring method for sparsely measured distribution grids | |
Alberdi et al. | Overhead line ampacity forecasting with a focus on safety | |
KR101983807B1 (ko) | 리덕션 기법을 이용한 송전손실계수 제공 방법 및 그 장치 | |
Olfatinezhad et al. | Flexible two‐stage robust model for moving the transmission and reactive power infrastructures expansion planning towards power system integration of renewables |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20101215 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120131 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120228 |
|
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: 20121204 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20121211 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5159695 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20151221 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20151221 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
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 |