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
Application number
JP3758895A
Other languages
English (en)
Other versions
JP3312522B2 (ja
Inventor
Yoshikazu Fukuyama
良和 福山
Shuichi Kikuchi
秀一 菊地
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.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric Co Ltd
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 Fuji Electric Co Ltd filed Critical Fuji Electric Co Ltd
Priority to JP3758895A priority Critical patent/JP3312522B2/ja
Publication of JPH08214458A publication Critical patent/JPH08214458A/ja
Application granted granted Critical
Publication of JP3312522B2 publication Critical patent/JP3312522B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

(57)【要約】 【目的】 放射状配電系統を対象とする潮流計算を短時
間で行なう。 【構成】 系統構成、システム容量、変電所電圧・位相
角、収束判定基準、負荷量等のデータをコンピュータに
入力するデータ入力ステップ(ステップ1)と、負荷量
の初期計算値をメインフィーダ及び分岐フィーダの末端
ノードから加算して各フィーダの先頭ノードにおける状
態変数の初期値を求める初期値計算ステップ(ステップ
2)と、前記状態変数を用いて各フィーダの各ノードに
おける状態量を先頭ノードから末端ノード方向へ逐次計
算する系統状態量の計算ステップ(ステップ3)と、各
フィーダの末端ノードにおけるエラー分だけ、各フィー
ダの先頭ノードの状態変数を修正する状態変数の修正ス
テップ(ステップ4)と、末端ノードにおける状態変数
を判定基準と比較して収束判定を行なう収束判定ステッ
プ(ステップ5)とを有する。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は放射状配電系統用潮流計
算方法に関し、詳しくは、配電営業所の制御用コンピュ
ータにより、放射状配電系統における系統構成、システ
ム容量、負荷量、変電所データ等を用いて潮流計算を行
なう方法に関する。
【0002】
【従来の技術及び発明が解決しようとする課題】潮流計
算は、現在系統における負荷や発電力の変化、送電線の
停止等に際して電力系統の運用方法を検討する場合や、
将来系統において発電所、送電線、変電所の新増設計画
を立てる場合等に必要不可欠のものとなっている。この
潮流計算方法は、本来的に、ループ系統を含む一般的な
配電系統用に開発されたものであるが、主としてループ
系統である高圧配電系統への適用を目的としている。
【0003】しかるに従来では、放射状配電系統に対し
ても従前の潮流計算方法をそのまま適用しており、その
場合には、放射状配電系統の特性上、収束特性が悪く、
計算に長時間を有するという問題があった。
【0004】本発明は上記問題点を解決するためになさ
れたもので、その目的とするところは、放射状配電系統
を対象として短時間で実行できるようにした放射状配電
系統用潮流計算方法を提供することにある。
【0005】
【課題を解決するための手段及び作用】上記目的を達成
するため、本発明は、放射状配電系統に関する種々のデ
ータを入力するデータ入力ステップと、状態変数の初期
値を求める初期値計算ステップと、前進計算としての系
統状態量の計算ステップと、後進計算としての状態変数
の修正ステップと、状態変数の収束判定ステップとを順
次実行することにより実現される。
【0006】すなわち、データ入力ステップでは、系統
構成、システム容量、変電所電圧・位相角、収束判定基
準、負荷量等のデータをコンピュータに入力し、初期値
計算ステップでは、負荷量の初期計算値をメインフィー
ダ及び分岐フィーダの末端ノードから加算して各フィー
ダの先頭ノードにおける状態変数の初期値を求める。ま
た、系統状態量の計算ステップでは、前記状態変数を用
いて各フィーダの各ノードにおける状態量を先頭ノード
から末端ノード方向へ逐次計算し、状態変数の修正ステ
ップでは、各フィーダの末端ノードにおけるエラー分だ
け、各フィーダの先頭ノードの状態変数を修正する。更
に、収束判定ステップでは、末端ノードにおける状態変
数を判定基準と比較して収束判定を行なう。なお、状態
変数には、例えば各フィーダに流れ込む電力が用いられ
る。
【0007】
【実施例】以下、図に沿って本発明の実施例を説明す
る。まず、本発明では、計算の対象を放射状配電系統に
特化し、状態変数を減らした前進計算及び後進計算から
なる潮流計算方式を用いる。ここで、状態変数はメイン
フィーダ及び各分岐フィーダの先頭ノードの状態量と
し、これらの状態変数を用いて、当該フィーダの各ノー
ドの系統状態を送電端から末端方向に向かって逐次計算
する(前進計算)。
【0008】図1は、1フィーダの系統例であり、送電
端のノードにおける有効電力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の計算式である。
【0009】
【数1】
【0010】
【数2】
【0011】
【数3】
【0012】
【数4】
【0013】これらの式において、Pk:ノードkから
の流出有効電力、Ploss,k:ノードkにおける損失電
力、Vk:ノードkの電圧、Qk:ノードkからの流出無
効電力、PLk:ノードkの有効電力負荷、QLk:ノード
kの無効電力負荷、rk:ノードkの線路抵抗、xk:ノ
ードkの線路インピーダンス、Yk1,Yk2:ノードkの
線路対地容量である。
【0014】また、図2は、線路対地容量を考慮した1
フィーダの系統例であるが、上記各数式におけるQk'は
図2のQk'に相当している。更に、数式3のVk+1の式
における「*」は電流iの共役を示す。
【0015】次に、末端ノードでは流出電力はない(す
なわちPn=Qn=0)という仮定から、状態変数を修正
することができる(後進計算)。初期値として与える先
頭ノードの初期状態量は推定値であるため、末端ノード
の電力量分だけ先頭ノードの状態量を修正し、末端ノー
ドの電力が零になるまで収束計算を行なう。一般に分岐
フィーダのある放射状系統では、フィーダの状態量によ
り各分岐フィーダの状態量も変更されるので、後進計算
は各分岐フィーダの末端ノードから先頭ノードに向けて
行なう。
【0016】図3は上記前進計算、後進計算の概念を示
している。図において、F1はメインフィーダ、F2,F
3は分岐フィーダ、N11〜N1n,N21〜N2n,N31〜N
3nはノードであり、N11,N21,N31が先頭ノード、N
1n,N2n,N3nが末端ノードである。
【0017】前進計算においては、各フィーダにつき図
の上段の〜の順に各ノードでの状態量を計算してい
く。この際、状態変数を各分岐フィーダに流れ込む電力
とした場合、メインフィーダF1の2番目のノードN12
から分岐フィーダF2の先頭ノードN21に流れ込んだ電
力がこのフィーダF2における計算の基礎的な状態変数
となる。同様にして、フィーダF2の2番目のノードN
22から分岐フィーダF3の先頭ノードN31に流れ込んだ
電力がこのフィーダF3における計算の基礎的な状態変
数となる。一方、後進計算においては、各フィーダにつ
き図の下段の〜の順に状態変数を修正していく。
【0018】本実施例によれば、状態変数が極端に少な
いことから、すべてのノード状態を状態変数とする従来
の方法に比べて、高速計算が可能になる。
【0019】図4はこの実施例の潮流計算手順を示すフ
ローチャートである。以下、順に説明する。 (1)データ入力 系統構成、システム容量(MVA,KVA)、変電所電
圧・位相角、収束判定基準、負荷量等のデータをコンピ
ュータに入力する(S1)。
【0020】(2)初期値計算 次に、各電圧の初期推定値を用いて負荷量の初期計算を
行ない、末端ノードから加算した値を各分岐フィーダ
(例えば図3におけるフィーダF2)への流出電力(状
態変数の一つ)の初期値とする(S2)。
【0021】(3)系統状態量の計算(前進計算) 求めた状態変数を用いて、前記数式1〜数式4等によ
り、電源(先頭ノード)から末端ノードに向かって逐次
的に系統状態量を計算する(S3)。 (4)状態変数の修正(後進計算) 末端ノードから流出する電力はないという事実に基づい
て、状態変数を修正する。つまり、末端ノードにおける
エラー分だけ各分岐フィーダの先頭ノードの流入電力を
修正する。分岐フィーダの流入電力の修正量は、当該分
岐フィーダの修正量と当該分岐フィーダからの分岐フィ
ーダの修正量の総和となる(S4)。 (5)収束判定 末端ノードの電力の絶対値がある判定基準より小さくな
ったかどうかにより、収束判定を行なう(S5)。
【0022】状態変数を各フィーダに流れ込む電力とし
て、207,354,528ノードシステムにつき潮流
計算時間を測定したシミュレーション結果を表1に示
す。
【0023】
【表1】
【0024】このテストシステムは三つの同じ部分系統
が変電所に接続されているシステムとした。一般に我が
国の配電系統は、変電所から複数のフィーダが出てお
り、各フィーダができるだけバランスするように負荷が
平均的に分散されている。従って、実系統もこのテスト
系統に非常に近い形態であると考えられる。
【0025】表1に示された計算時間は、ループ系統を
対象とした従来の潮流計算時間に対して1/2〜1/8
の値であり、本発明によれば大幅に計算時間を短縮でき
ることが確認された。なお、シミュレーションはトラン
スピュータ(約10MIPS)を用いて行なったもので
ある。
【0026】
【発明の効果】以上述べたように本発明によれば、放射
状配電系統のみを対象として、データ入力ステップ、初
期値計算ステップ、少ない状態変数を用いて先頭ノード
から末端ノード方向へ逐次計算する系統状態量の計算ス
テップ、末端ノードから先頭ノードに向かう状態変数の
修正ステップ、及び収束判定ステップを順次実行するこ
とにより、従来よりも短時間で潮流計算を行なうことが
できる。
【図面の簡単な説明】
【図1】1フィーダの系統例を示す図である。
【図2】線路対地容量を考慮した1フィーダの系統例を
示す図である。
【図3】実施例における前進計算及び後進計算の概念図
である。
【図4】実施例における計算手順を示すフローチャート
である。
【符号の説明】
1 メインフィーダ F2,F3 分岐フィーダ N11〜N1n,N21〜N2n,N31〜N3n ノード

Claims (2)

    【特許請求の範囲】
  1. 【請求項1】 放射状配電系統を対象とした潮流計算方
    法において、 系統構成、システム容量、変電所電圧・位相角、収束判
    定基準、負荷量等のデータをコンピュータに入力するデ
    ータ入力ステップと、 負荷量の初期計算値をメインフィーダ及び分岐フィーダ
    の末端ノードから加算して各フィーダの先頭ノードにお
    ける状態変数の初期値を求める初期値計算ステップと、 前記状態変数を用いて各フィーダの各ノードにおける状
    態量を先頭ノードから末端ノード方向へ逐次計算する系
    統状態量の計算ステップと、 各フィーダの末端ノードにおけるエラー分だけ、各フィ
    ーダの先頭ノードの状態変数を修正する状態変数の修正
    ステップと、 末端ノードにおける状態変数を判定基準と比較して収束
    判定を行なう収束判定ステップと、 を有することを特徴とする放射状配電系統用潮流計算方
    法。
  2. 【請求項2】 状態変数が各フィーダに流れ込む電力で
    ある請求項1記載の放射状配電系統用潮流計算方法。
JP3758895A 1995-02-02 1995-02-02 潮流計算装置 Expired - Lifetime JP3312522B2 (ja)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5743195B2 (ja) * 2011-04-15 2015-07-01 富士電機株式会社 ループ状系統用潮流計算方法及び潮流計算装置

Cited By (8)

* Cited by examiner, † Cited by third party
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