JP2018191482A - 最適潮流計算装置、最適潮流計算方法、および最適潮流計算プログラム - Google Patents
最適潮流計算装置、最適潮流計算方法、および最適潮流計算プログラム Download PDFInfo
- Publication number
- JP2018191482A JP2018191482A JP2017094459A JP2017094459A JP2018191482A JP 2018191482 A JP2018191482 A JP 2018191482A JP 2017094459 A JP2017094459 A JP 2017094459A JP 2017094459 A JP2017094459 A JP 2017094459A JP 2018191482 A JP2018191482 A JP 2018191482A
- Authority
- JP
- Japan
- Prior art keywords
- cross
- section
- constraint
- power system
- flow 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.)
- Granted
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
Description
以下では、まず、単断面の場合の最適潮流計算の定式化について、非特許文献1(電力系統標準モデルの普及・拡充調査専門委員会、「拡充系統モデル報告書(マニュアル)」、電気学会、2001年10月)の3.1.1 最適潮流計算の基礎的事項に基づいて説明する。次に、単断面を複数断面に拡張した場合の最適潮流計算の定式化について説明する。その次に、実施の形態1の最適潮流計算装置の装置構成と動作について説明する。
最適潮流計算は、種々の制約条件の下で、単一の目的関数、すなわち、
f(x,u,z) …(1)
が最小となるような変数x,u,zの最適値を求めることである。
KA×fA(x,u,z)+KB×fB(x,u,z) …(2)
を目的関数とすることもできる。
gA(x,u,z)=0 …(3)
gB(x,u,z)=0 …(4)
gC(x,u,z)≦0 …(5)
のように表わすことができる。
次に、上記の単断面の定式化を複数断面に拡張する方法について説明する。複数断面として、現時点の系統状態の他、電源の供給量および負荷の需要量の変化を考慮した将来の系統状態、系統故障などの理由により一部の系統が切り離された状態などを考えることができる。以下では、簡単のために、第1断面および第2断面からなる2つの断面の場合について説明するが、より多断面の場合も同様である。
f(x1,u1,z1;x2,u2,z2) …(6)
のように、各断面の変数を用いて目的関数が表される。なお、目的関数は全ての変数を用いる必要はなく、たとえば、第1断面の変数のみを用いて目的関数を表してもよい。
KA1×fA1(x1,u1,z1)+KB1×fB1(x1,u1,z1)
+KA2×fA2(x2,u2,z2)+KB2×fB(x2,u2,z2) …(7)
のように表される。ここで、第2断面よりも第1断面を重視する場合には、係数KA2およびKB2よりも係数KA1およびKB1を大きな値に設定してもよい。
gA1(x1,u1,z1)=0 …(8A)
gB1(x1,u1,z1)=0 …(9A)
gC1(x1,u1,z1)≦0 …(10A)
が定められる。同様に、第2断面について、上式(3)〜(5)にそれぞれ対応する、
gA2(x2,u2,z2)=0 …(8B)
gB2(x2,u2,z2)=0 …(9B)
gC2(x2,u2,z2)≦0 …(10B)
が定められる。上記の(8A)〜(10A)式および(8B)〜(10B)式は、断面ごとに少なくとも一部の変数の制約条件を規定する複数の第1の制約式である。
hA(x1p,x2p)=0 …(11A)
または hA(x1p,x2p)≦0 …(11B)
の制約式が追加される。上式(11A)は、第1断面の一部の電圧解x1pの値に応じて第2断面の対応する電圧解x2pの値が決まることを表す。上式(11B)は、第2断面の一部の電圧解x1pの値に応じて第2断面の対応する電圧解x2pの範囲が決まることを表す。
hB(u1p,u2p)=0 …(12A)
または hB(u1p,u2p)≦0 …(12B)
の制約式が追加される。第1断面の従属変数z1の一部をz1pとし、この一部の従属変数z1pにそれぞれ対応する第2断面の従属変数をz2pとすれば、
hC(z1p,z2p)=0 …(13A)
または hC(z1p,z2p)≦0 …(13B)
の制約式が追加される。
以下、具体例を挙げながら、複数断面の場合の最適潮流計算の定式化についてさらに説明する。
V(11)=V(21) …(14)
が考えられる。上式(14)において、V(11)は母線11の電圧を表し、V(21)は母線(21)の電圧を表す。
Tap(32)=Tap(42) …(15)
が考えられる。上式(15)において、Tap(32)は変圧器32のタップ比を表し、Tap(42)は変圧器42のタップ比を表す。
Q(34)=Q(44) …(16)
Q(35)=Q(45) …(17)
が考えられる。上記(16)において、Q(34)は電力用コンデンサ34による無効電力の投入量を表し、Q(44)は電力用コンデンサ44による無効電力の投入量を表す。上式(17)において、Q(35)は分路リアクトル35による無効電力の投入量を表し、Q(45)は分路リアクトル45による無効電力の投入量を表す。
F=a×P(31)×P(31)+b×P(31)+c …(18)
のように設定することができる。ただし、上式(18)は、発電コストを2次式で表した例である。
|V(12)−TargetV(12)|2 …(19)
のように表される。ただし、|…|2は大きさの2乗を表す。
|Q(34)−Qnow(34)|2 …(20)
のように表される。
F=Kc(31)×{a(31)×P(31)2+b(31)×P(31)+c(31)}
+Kc(41)×{a(41)×P(41)2+b(41)×P(41)+c(41)}
+Kv(12)×|V(12)−TargetV(12)|2+Kv(22)×|V(22)−TargetV(22)|2
+Kq(34)×|Q(34) − Qnow(34)|2+Kq(44)×|Q(44) −Qnow(44)|2
+Kq(35)×|Q(35) − Qnow(35)|2+Kq(45)×|Q(45) −Qnow(45)|2 …(21)
と表される。ただし、Kc(31)およびKc(41)は、それぞれ発電機31,41の発電コストに係る重み係数である。Kv(12)およびKv(22)は、それぞれ母線12,22の目標電圧からの差分の2乗に係る重み係数である。Kq(34),Kq(35),Kq(44),Kq(45)は、それぞれ調相設備34,35,44,45の現在値からの差分の2乗に係る重み係数である。a(31),b(31),c(31)は、発電機31の発電コスト係数であり、a(41),b(41),c(41)は、発電機41の発電コスト係数である(発電コストを2次式で表した例)。TargetV(12)およびTargetV(22)は、それぞれ母線12,22の目標電圧である。Qnow(34),Qnow(35),Qnow(44),Qnow(45)は、それぞれ調相設備34,35,44,45の無効電力出力の現在値を表す。Q(34),Q(35),Q(44),Q(45)は、それぞれ現在と異なる断面での調相設備34,35,44,45の無効電力出力を表す。
次に、実施の形態1による最適潮流計算装置100の構成および動作について説明する。実施の形態1に係る最適潮流計算装置100は、入力された電力および断面情報(前述のパラメータに対応する)に基づいて、変圧器変圧比、調相設備投入量、発電機出力、発電機端子電圧などの制御機器の適切な状態(前述の変数x,u,zに対応する)を出力する。
実施の形態1の最適潮流計算装置100によれば、複数の系統断面を一つの最適潮流計算問題として定式化して解を得ることが可能となる。したがって、現在状態からの改善内容を得るだけでなく、故障後および将来断面など、別の断面の系統状態の改善および悪化抑制を図ることができる。さらに、このようにして算出された適切な系統状態を用いることにより、現在状態に対して適切な制御対象の選択および制御量を算出することができる。これにより、電力系統監視制御分野および系統安定化分野において最適潮流計算結果を利用することができる。
実施の形態2では、実施の形態1の最適潮流計算装置100の機能に加えて、最適潮流計算においてパラメータとして用いられる断面情報のうち、現状状態以外の断面情報を自動的に作成するようにしたものである。これにより、故障後の断面、将来予測断面などを現在状態と同時に評価することが容易にできるので、現在状態の断面およびその他の断面を考慮した適切な系統状態を容易に得ることができる。以下、図面を参照して具体的に説明する。
Claims (9)
- 複数の制約式および目的関数において固定値として与えられる複数のパラメータに対応する電力系統の情報を、複数断面を構成する断面ごとに格納するデータ格納部と、
前記断面ごとの前記複数のパラメータに基づいて最適潮流計算を実行することによって、前記複数の制約式および前記目的関数に含まれる可変値である複数の変数の最適値を決定する最適潮流計算部とを備え、
前記複数の制約式は、
前記断面ごとに少なくとも一部の変数の関係を規定する複数の第1の制約式と、
異なる断面間で互いに対応する少なくとも一部の変数の関係を規定する1つまたは複数の第2の制約式とを含む、最適潮流計算装置。 - 前記断面ごとに前記複数のパラメータに対応する電力系統の情報の入力を受け付けるデータ入力部をさらに備える、請求項1に記載の最適潮流計算装置。
- 前記複数断面のうちの第1断面について、前記複数のパラメータに対応する電力系統の情報の入力を受け付けるデータ入力部と、
前記第1断面の電力系統の情報に基づいて前記複数断面のうちの前記第1断面と異なる他断面について、前記複数のパラメータに対応する電力系統の情報を生成する断面情報作成部とを備える、請求項1に記載の最適潮流計算装置。 - 前記複数断面の電力系統は、第1断面の電力系統と、故障によって前記第1断面の電力系統から一部が遮断された残りの電力系統である第2断面の電力系統とを含む、請求項1〜3のいずれか1項に記載の最適潮流計算装置。
- 前記第1断面の電力系統と前記第2断面の電力系統との両方に、第1母線と、前記第1母線に接続された同じ発電機とが含まれる場合、
前記複数の第2の制約式は、前記第1断面における前記第1母線の電圧と前記第2断面における前記第1母線の電圧とが等しいという制約式を含む、請求項4に記載の最適潮流計算装置。 - 前記第1断面の電力系統と前記第2断面の電力系統との両方に同じ変圧器が含まれる場合、
前記複数の第2の制約式は、前記第1断面における前記変圧器のタップ位置と前記第2断面における前記変圧器のタップ位置とが等しいという制約式を含む、請求項4または5に記載の最適潮流計算装置。 - 前記第1断面の電力系統と前記第2断面の電力系統との両方に同じ調相設備が含まれる場合、
前記複数の第2の制約式は、前記第1断面における前記調相設備による無効電力の投入量と前記第2断面における前記調相設備による無効電力の投入量とが等しいという制約式を含む、請求項4〜6のいずれか1項に記載の最適潮流計算装置。 - 複数の制約式および目的関数において固定値として与えられる複数のパラメータに対応する電力系統の情報を、複数断面を構成する断面ごとに記憶装置に格納するステップと、
プロセッサが、前記断面ごとの前記複数のパラメータに基づいて最適潮流計算を実行することによって、前記複数の制約式および前記目的関数に含まれる可変値である複数の変数の最適値を決定するステップとを備え、
前記複数の制約式は、
前記断面ごとに少なくとも一部の変数の関係を規定する複数の第1の制約式と、
異なる断面間で互いに対応する少なくとも一部の変数の関係を規定する1つまたは複数の第2の制約式とを含む、最適潮流計算方法。 - 複数の制約式および目的関数において固定値として与えられる複数のパラメータに対応する電力系統の情報を、複数断面を構成する断面ごとに格納するステップと、
前記断面ごとの前記複数のパラメータに基づいて最適潮流計算を実行することによって、前記複数の制約式および前記目的関数に含まれる可変値である複数の変数の最適値を決定するステップとをコンピュータに実行させるための最低潮流計算プログラムであって、
前記複数の制約式は、
前記断面ごとに少なくとも一部の変数の関係を規定する複数の第1の制約式と、
異なる断面間で互いに対応する少なくとも一部の変数の関係を規定する1つまたは複数の第2の制約式とを含む、最適潮流計算プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017094459A JP6785717B2 (ja) | 2017-05-11 | 2017-05-11 | 最適潮流計算装置、最適潮流計算方法、および最適潮流計算プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017094459A JP6785717B2 (ja) | 2017-05-11 | 2017-05-11 | 最適潮流計算装置、最適潮流計算方法、および最適潮流計算プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018191482A true JP2018191482A (ja) | 2018-11-29 |
JP6785717B2 JP6785717B2 (ja) | 2020-11-18 |
Family
ID=64480287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017094459A Active JP6785717B2 (ja) | 2017-05-11 | 2017-05-11 | 最適潮流計算装置、最適潮流計算方法、および最適潮流計算プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6785717B2 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2021083235A (ja) * | 2019-11-20 | 2021-05-27 | 株式会社日立製作所 | 電力需給計画装置 |
WO2021260788A1 (ja) | 2020-06-23 | 2021-12-30 | 三菱電機株式会社 | 調相設備の制御装置、および調相設備の制御方法 |
WO2022118495A1 (ja) | 2020-12-04 | 2022-06-09 | 株式会社日立製作所 | 電力系統運用計画作成支援装置および方法 |
JP2022532091A (ja) * | 2019-05-07 | 2022-07-13 | ヒタチ・エナジー・スウィツァーランド・アクチェンゲゼルシャフト | 動的予測、集約、および妥当性確認のための機器および方法 |
WO2022224730A1 (ja) * | 2021-04-21 | 2022-10-27 | 株式会社日立製作所 | 電力系統監視装置およびその監視方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11259450A (ja) * | 1998-03-09 | 1999-09-24 | Hitachi Ltd | 最適な出力決定方法および装置 |
JP2001339861A (ja) * | 2000-05-31 | 2001-12-07 | Rakunan Chin | 送電可能容量計算システム、最適潮流計算方法及び発電機の負荷配分を計算する方法並びにプログラム |
JP2010011735A (ja) * | 2004-12-14 | 2010-01-14 | Tokyo Electric Power Co Inc:The | 最適潮流計算方法及び最適潮流計算装置 |
JP2012080708A (ja) * | 2010-10-05 | 2012-04-19 | Tokyo Electric Power Co Inc:The | 系統運用策定装置 |
JP2016194849A (ja) * | 2015-04-01 | 2016-11-17 | 富士電機株式会社 | 潮流計算装置、潮流計算方法、及びプログラム |
-
2017
- 2017-05-11 JP JP2017094459A patent/JP6785717B2/ja active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11259450A (ja) * | 1998-03-09 | 1999-09-24 | Hitachi Ltd | 最適な出力決定方法および装置 |
JP2001339861A (ja) * | 2000-05-31 | 2001-12-07 | Rakunan Chin | 送電可能容量計算システム、最適潮流計算方法及び発電機の負荷配分を計算する方法並びにプログラム |
JP2010011735A (ja) * | 2004-12-14 | 2010-01-14 | Tokyo Electric Power Co Inc:The | 最適潮流計算方法及び最適潮流計算装置 |
JP2012080708A (ja) * | 2010-10-05 | 2012-04-19 | Tokyo Electric Power Co Inc:The | 系統運用策定装置 |
JP2016194849A (ja) * | 2015-04-01 | 2016-11-17 | 富士電機株式会社 | 潮流計算装置、潮流計算方法、及びプログラム |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2022532091A (ja) * | 2019-05-07 | 2022-07-13 | ヒタチ・エナジー・スウィツァーランド・アクチェンゲゼルシャフト | 動的予測、集約、および妥当性確認のための機器および方法 |
JP7198376B2 (ja) | 2019-05-07 | 2022-12-28 | ヒタチ・エナジー・スウィツァーランド・アクチェンゲゼルシャフト | 動的予測、集約、および妥当性確認のための機器および方法 |
JP2021083235A (ja) * | 2019-11-20 | 2021-05-27 | 株式会社日立製作所 | 電力需給計画装置 |
JP7245761B2 (ja) | 2019-11-20 | 2023-03-24 | 株式会社日立製作所 | 電力需給計画装置 |
WO2021260788A1 (ja) | 2020-06-23 | 2021-12-30 | 三菱電機株式会社 | 調相設備の制御装置、および調相設備の制御方法 |
JPWO2021260788A1 (ja) * | 2020-06-23 | 2021-12-30 | ||
JP7366263B2 (ja) | 2020-06-23 | 2023-10-20 | 三菱電機株式会社 | 調相設備の制御装置、および調相設備の制御方法 |
WO2022118495A1 (ja) | 2020-12-04 | 2022-06-09 | 株式会社日立製作所 | 電力系統運用計画作成支援装置および方法 |
JP7514748B2 (ja) | 2020-12-04 | 2024-07-11 | 株式会社日立製作所 | 電力系統運用計画作成支援装置および方法 |
WO2022224730A1 (ja) * | 2021-04-21 | 2022-10-27 | 株式会社日立製作所 | 電力系統監視装置およびその監視方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6785717B2 (ja) | 2020-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6785717B2 (ja) | 最適潮流計算装置、最適潮流計算方法、および最適潮流計算プログラム | |
Wang et al. | Coordination of multiple energy storage units in a low-voltage distribution network | |
Pandiarajan et al. | Fuzzy harmony search algorithm based optimal power flow for power system security enhancement | |
Rabiee et al. | Voltage security constrained multi-period optimal reactive power flow using benders and optimality condition decompositions | |
Mahmoud | Voltage stability analysis of radial distribution networks using catastrophe theory | |
Jayasankar et al. | Estimation of voltage stability index for power system employing artificial neural network technique and TCSC placement | |
Xu et al. | Case-based reasoning for coordinated voltage control on distribution networks | |
SE540926C2 (en) | Converter station power set point analysis system and method | |
JP6259783B2 (ja) | 送電設備計画支援システム及び方法 | |
Hajian et al. | A chance-constrained optimization approach for control of transmission voltages | |
Meah et al. | A new simplified adaptive control scheme for multi-terminal HVDC transmission systems | |
Mageshvaran et al. | Steady state load shedding to mitigate blackout in power systems using an improved harmony search algorithm | |
Nasrolahpour et al. | Optimal transmission congestion management by means of substation reconfiguration | |
CN114336634B (zh) | 一种电网***的潮流计算方法、装置和设备 | |
Wang et al. | A three-machine equivalent system frequency response model and its closed-form solution | |
Nie et al. | Wide-area damping power system control based on robust self-triggered model predictive control | |
Haileselassie et al. | Secondary control in multi-terminal VSC-HVDC transmission system | |
Mendonca et al. | Robust tuning of power system stabilisers to install in wind energy conversion systems | |
Liang et al. | Adaptive critic design based dynamic optimal power flow controller for a smart grid | |
CN111211583A (zh) | 用于混合ac/dc电网的网络*** | |
Vrakopoulou et al. | Optimal policy-based control of generation and HVDC lines in power systems under uncertainty | |
Ranamuka et al. | Mitigating tap changer limit cycles in modern electricity networks embedded with local generation units | |
Pohl et al. | Dc-sensitivities for impedance controllers in an agent-based power flow control system | |
Swain et al. | DSM for all‐day voltage profile improvement in a microgrid | |
Wang et al. | Frequency-based demand side response considering the discontinuity of the ToU tariff |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190723 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200623 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200630 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200908 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200917 |
|
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: 20200929 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20201027 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6785717 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |