JPH08214458A - 放射状配電系統用潮流計算方法 - Google Patents
放射状配電系統用潮流計算方法Info
- Publication number
- JPH08214458A JPH08214458A JP3758895A JP3758895A JPH08214458A JP H08214458 A JPH08214458 A JP H08214458A JP 3758895 A JP3758895 A JP 3758895A JP 3758895 A JP3758895 A JP 3758895A JP H08214458 A JPH08214458 A JP H08214458A
- Authority
- JP
- Japan
- Prior art keywords
- feeder
- node
- state
- state variables
- distribution system
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims description 48
- 238000000034 method Methods 0.000 abstract description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
間で行なう。 【構成】 系統構成、システム容量、変電所電圧・位相
角、収束判定基準、負荷量等のデータをコンピュータに
入力するデータ入力ステップ(ステップ1)と、負荷量
の初期計算値をメインフィーダ及び分岐フィーダの末端
ノードから加算して各フィーダの先頭ノードにおける状
態変数の初期値を求める初期値計算ステップ(ステップ
2)と、前記状態変数を用いて各フィーダの各ノードに
おける状態量を先頭ノードから末端ノード方向へ逐次計
算する系統状態量の計算ステップ(ステップ3)と、各
フィーダの末端ノードにおけるエラー分だけ、各フィー
ダの先頭ノードの状態変数を修正する状態変数の修正ス
テップ(ステップ4)と、末端ノードにおける状態変数
を判定基準と比較して収束判定を行なう収束判定ステッ
プ(ステップ5)とを有する。
Description
算方法に関し、詳しくは、配電営業所の制御用コンピュ
ータにより、放射状配電系統における系統構成、システ
ム容量、負荷量、変電所データ等を用いて潮流計算を行
なう方法に関する。
算は、現在系統における負荷や発電力の変化、送電線の
停止等に際して電力系統の運用方法を検討する場合や、
将来系統において発電所、送電線、変電所の新増設計画
を立てる場合等に必要不可欠のものとなっている。この
潮流計算方法は、本来的に、ループ系統を含む一般的な
配電系統用に開発されたものであるが、主としてループ
系統である高圧配電系統への適用を目的としている。
ても従前の潮流計算方法をそのまま適用しており、その
場合には、放射状配電系統の特性上、収束特性が悪く、
計算に長時間を有するという問題があった。
れたもので、その目的とするところは、放射状配電系統
を対象として短時間で実行できるようにした放射状配電
系統用潮流計算方法を提供することにある。
するため、本発明は、放射状配電系統に関する種々のデ
ータを入力するデータ入力ステップと、状態変数の初期
値を求める初期値計算ステップと、前進計算としての系
統状態量の計算ステップと、後進計算としての状態変数
の修正ステップと、状態変数の収束判定ステップとを順
次実行することにより実現される。
構成、システム容量、変電所電圧・位相角、収束判定基
準、負荷量等のデータをコンピュータに入力し、初期値
計算ステップでは、負荷量の初期計算値をメインフィー
ダ及び分岐フィーダの末端ノードから加算して各フィー
ダの先頭ノードにおける状態変数の初期値を求める。ま
た、系統状態量の計算ステップでは、前記状態変数を用
いて各フィーダの各ノードにおける状態量を先頭ノード
から末端ノード方向へ逐次計算し、状態変数の修正ステ
ップでは、各フィーダの末端ノードにおけるエラー分だ
け、各フィーダの先頭ノードの状態変数を修正する。更
に、収束判定ステップでは、末端ノードにおける状態変
数を判定基準と比較して収束判定を行なう。なお、状態
変数には、例えば各フィーダに流れ込む電力が用いられ
る。
る。まず、本発明では、計算の対象を放射状配電系統に
特化し、状態変数を減らした前進計算及び後進計算から
なる潮流計算方式を用いる。ここで、状態変数はメイン
フィーダ及び各分岐フィーダの先頭ノードの状態量と
し、これらの状態変数を用いて、当該フィーダの各ノー
ドの系統状態を送電端から末端方向に向かって逐次計算
する(前進計算)。
端のノードにおける有効電力P、無効電力Q、系統電圧
V等の状態変数x0(f(x0))が与えられれば、その
下流側のノードの状態変数x1(f(x1))が求めら
れ、以下、末端のxnまで逐次、計算することができ
る。以下に、状態変数としての有効電力P(数式1)、
無効電力Q(数式2)及び系統電圧V(数式4)の詳細
な計算式を示す。なお、数式3は、数式4における系統
電圧Vk+1(=ek+1+jfk+1)とその有効成分ek+1及
び無効成分fk+1の計算式である。
の流出有効電力、Ploss,k:ノードkにおける損失電
力、Vk:ノードkの電圧、Qk:ノードkからの流出無
効電力、PLk:ノードkの有効電力負荷、QLk:ノード
kの無効電力負荷、rk:ノードkの線路抵抗、xk:ノ
ードkの線路インピーダンス、Yk1,Yk2:ノードkの
線路対地容量である。
フィーダの系統例であるが、上記各数式におけるQk'は
図2のQk'に相当している。更に、数式3のVk+1の式
における「*」は電流iの共役を示す。
なわちPn=Qn=0)という仮定から、状態変数を修正
することができる(後進計算)。初期値として与える先
頭ノードの初期状態量は推定値であるため、末端ノード
の電力量分だけ先頭ノードの状態量を修正し、末端ノー
ドの電力が零になるまで収束計算を行なう。一般に分岐
フィーダのある放射状系統では、フィーダの状態量によ
り各分岐フィーダの状態量も変更されるので、後進計算
は各分岐フィーダの末端ノードから先頭ノードに向けて
行なう。
している。図において、F1はメインフィーダ、F2,F
3は分岐フィーダ、N11〜N1n,N21〜N2n,N31〜N
3nはノードであり、N11,N21,N31が先頭ノード、N
1n,N2n,N3nが末端ノードである。
の上段の〜の順に各ノードでの状態量を計算してい
く。この際、状態変数を各分岐フィーダに流れ込む電力
とした場合、メインフィーダF1の2番目のノードN12
から分岐フィーダF2の先頭ノードN21に流れ込んだ電
力がこのフィーダF2における計算の基礎的な状態変数
となる。同様にして、フィーダF2の2番目のノードN
22から分岐フィーダF3の先頭ノードN31に流れ込んだ
電力がこのフィーダF3における計算の基礎的な状態変
数となる。一方、後進計算においては、各フィーダにつ
き図の下段の〜の順に状態変数を修正していく。
いことから、すべてのノード状態を状態変数とする従来
の方法に比べて、高速計算が可能になる。
ローチャートである。以下、順に説明する。 (1)データ入力 系統構成、システム容量(MVA,KVA)、変電所電
圧・位相角、収束判定基準、負荷量等のデータをコンピ
ュータに入力する(S1)。
行ない、末端ノードから加算した値を各分岐フィーダ
(例えば図3におけるフィーダF2)への流出電力(状
態変数の一つ)の初期値とする(S2)。
り、電源(先頭ノード)から末端ノードに向かって逐次
的に系統状態量を計算する(S3)。 (4)状態変数の修正(後進計算) 末端ノードから流出する電力はないという事実に基づい
て、状態変数を修正する。つまり、末端ノードにおける
エラー分だけ各分岐フィーダの先頭ノードの流入電力を
修正する。分岐フィーダの流入電力の修正量は、当該分
岐フィーダの修正量と当該分岐フィーダからの分岐フィ
ーダの修正量の総和となる(S4)。 (5)収束判定 末端ノードの電力の絶対値がある判定基準より小さくな
ったかどうかにより、収束判定を行なう(S5)。
て、207,354,528ノードシステムにつき潮流
計算時間を測定したシミュレーション結果を表1に示
す。
が変電所に接続されているシステムとした。一般に我が
国の配電系統は、変電所から複数のフィーダが出てお
り、各フィーダができるだけバランスするように負荷が
平均的に分散されている。従って、実系統もこのテスト
系統に非常に近い形態であると考えられる。
対象とした従来の潮流計算時間に対して1/2〜1/8
の値であり、本発明によれば大幅に計算時間を短縮でき
ることが確認された。なお、シミュレーションはトラン
スピュータ(約10MIPS)を用いて行なったもので
ある。
状配電系統のみを対象として、データ入力ステップ、初
期値計算ステップ、少ない状態変数を用いて先頭ノード
から末端ノード方向へ逐次計算する系統状態量の計算ス
テップ、末端ノードから先頭ノードに向かう状態変数の
修正ステップ、及び収束判定ステップを順次実行するこ
とにより、従来よりも短時間で潮流計算を行なうことが
できる。
示す図である。
である。
である。
Claims (2)
- 【請求項1】 放射状配電系統を対象とした潮流計算方
法において、 系統構成、システム容量、変電所電圧・位相角、収束判
定基準、負荷量等のデータをコンピュータに入力するデ
ータ入力ステップと、 負荷量の初期計算値をメインフィーダ及び分岐フィーダ
の末端ノードから加算して各フィーダの先頭ノードにお
ける状態変数の初期値を求める初期値計算ステップと、 前記状態変数を用いて各フィーダの各ノードにおける状
態量を先頭ノードから末端ノード方向へ逐次計算する系
統状態量の計算ステップと、 各フィーダの末端ノードにおけるエラー分だけ、各フィ
ーダの先頭ノードの状態変数を修正する状態変数の修正
ステップと、 末端ノードにおける状態変数を判定基準と比較して収束
判定を行なう収束判定ステップと、 を有することを特徴とする放射状配電系統用潮流計算方
法。 - 【請求項2】 状態変数が各フィーダに流れ込む電力で
ある請求項1記載の放射状配電系統用潮流計算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3758895A JP3312522B2 (ja) | 1995-02-02 | 1995-02-02 | 潮流計算装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3758895A JP3312522B2 (ja) | 1995-02-02 | 1995-02-02 | 潮流計算装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08214458A true JPH08214458A (ja) | 1996-08-20 |
JP3312522B2 JP3312522B2 (ja) | 2002-08-12 |
Family
ID=12501704
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP3758895A Expired - Lifetime JP3312522B2 (ja) | 1995-02-02 | 1995-02-02 | 潮流計算装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3312522B2 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100511839B1 (ko) * | 1997-02-24 | 2005-11-22 | 가부시끼가이샤 히다치 세이사꾸쇼 | 전력계통구성의 작성장치 |
JP2012005210A (ja) * | 2010-06-15 | 2012-01-05 | Hitachi Ltd | 配電系統潮流シミュレーション装置、配電系統潮流シミュレーション方法およびそのプログラム |
KR101221083B1 (ko) * | 2011-05-25 | 2013-01-11 | 주식회사 파워이십일 | 배전 시스템의 상태추정 방법 |
WO2014114107A1 (zh) * | 2013-01-22 | 2014-07-31 | 国家电网公司 | 一种分布式电源最优潮流优化方法 |
CN109980650A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种辐射型配电***潮流计算方法 |
CN111313425A (zh) * | 2020-01-15 | 2020-06-19 | 国网重庆市电力公司 | 基于变量空间最优选择的潮流模型线性化误差最小化方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5743195B2 (ja) * | 2011-04-15 | 2015-07-01 | 富士電機株式会社 | ループ状系統用潮流計算方法及び潮流計算装置 |
-
1995
- 1995-02-02 JP JP3758895A patent/JP3312522B2/ja not_active Expired - Lifetime
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100511839B1 (ko) * | 1997-02-24 | 2005-11-22 | 가부시끼가이샤 히다치 세이사꾸쇼 | 전력계통구성의 작성장치 |
JP2012005210A (ja) * | 2010-06-15 | 2012-01-05 | Hitachi Ltd | 配電系統潮流シミュレーション装置、配電系統潮流シミュレーション方法およびそのプログラム |
KR101221083B1 (ko) * | 2011-05-25 | 2013-01-11 | 주식회사 파워이십일 | 배전 시스템의 상태추정 방법 |
WO2014114107A1 (zh) * | 2013-01-22 | 2014-07-31 | 国家电网公司 | 一种分布式电源最优潮流优化方法 |
CN109980650A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种辐射型配电***潮流计算方法 |
CN109980650B (zh) * | 2019-05-21 | 2022-11-15 | 华北水利水电大学 | 一种辐射型配电***潮流计算方法 |
CN111313425A (zh) * | 2020-01-15 | 2020-06-19 | 国网重庆市电力公司 | 基于变量空间最优选择的潮流模型线性化误差最小化方法 |
CN111313425B (zh) * | 2020-01-15 | 2023-11-14 | 国网重庆市电力公司 | 基于变量空间最优选择的潮流模型线性化误差最小化方法 |
Also Published As
Publication number | Publication date |
---|---|
JP3312522B2 (ja) | 2002-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4101515B2 (ja) | 高電圧直流送電システムにおける有効電力制御 | |
Liu et al. | An efficient power flow algorithm for distribution systems with polynomial load | |
US8849614B2 (en) | System and method for monitoring and managing three-phase power flows in electrical transmission and distribution networks | |
Baran et al. | Distribution transformer models for branch current based feeder analysis | |
Kersting | The computation of neutral and dirt currents and power losses | |
Furlaneto et al. | Short circuit network equivalents of systems with inverter-based resources | |
JPH08214458A (ja) | 放射状配電系統用潮流計算方法 | |
Yan et al. | Analysis of radial distribution systems with embedded series FACTS devices using a fast line flow-based algorithm | |
Borzacchiello et al. | Unified formulation of a family of iterative solvers for power systems analysis | |
CN112865033B (zh) | 适应继电保护的全功率风电机组故障阻抗构建方法及*** | |
CN110133420A (zh) | 一种用于不平衡配电网的单相接地故障快速线性计算方法 | |
Verma et al. | Application of modified Gauss-Zbus iterations for solving the load flow problem in active distribution networks | |
CN113514731A (zh) | 一种含逆变型电源不平衡配电网的短路电流确定方法 | |
CN110208634B (zh) | 一种复杂电力***不对称短路电流直流分量获取方法 | |
CN112039061A (zh) | 基于电磁暂态仿真的潮流计算方法 | |
JP4483622B2 (ja) | 三相四線式負荷模擬装置 | |
Abdel-Akher et al. | An approach to determine a pair of power-flow solutions related to the voltage stability of unbalanced three-phase networks | |
CN110676877A (zh) | 一种孤岛微电网潮流检测方法 | |
CN115549093A (zh) | 一种新能源电力***在线建模与振荡分析的方法及*** | |
Medina et al. | A three phase load flow algorithm for shipboard power systems (SPS) | |
JP3719477B2 (ja) | 放射状配電系統における潮流計算方法 | |
CN108521143A (zh) | 一种计及不确定性的交直流混合微网并网鲁棒控制方法 | |
JP2018182982A (ja) | 潮流計算装置、潮流計算方法、及び潮流計算プログラム | |
JP2003143756A (ja) | 電力系統の縮約方法及び処理プログラム記録媒体 | |
Mok et al. | Power flow analysis for balanced and unbalanced radial distribution systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20020430 |
|
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 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090531 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090531 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100531 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110531 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110531 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120531 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120531 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130531 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130531 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140531 Year of fee payment: 12 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |