JP2013226005A - 配電網電力フロー分析システム及び方法 - Google Patents

配電網電力フロー分析システム及び方法 Download PDF

Info

Publication number
JP2013226005A
JP2013226005A JP2012097484A JP2012097484A JP2013226005A JP 2013226005 A JP2013226005 A JP 2013226005A JP 2012097484 A JP2012097484 A JP 2012097484A JP 2012097484 A JP2012097484 A JP 2012097484A JP 2013226005 A JP2013226005 A JP 2013226005A
Authority
JP
Japan
Prior art keywords
matrix
node
row
distribution network
relationship
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
JP2012097484A
Other languages
English (en)
Other versions
JP5426716B2 (ja
Inventor
Ting Jia Ou
庭嘉 歐
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.)
Institute of Nuclear Energy Research
Original Assignee
Institute of Nuclear Energy Research
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 Institute of Nuclear Energy Research filed Critical Institute of Nuclear Energy Research
Priority to JP2012097484A priority Critical patent/JP5426716B2/ja
Priority to EP12166226.6A priority patent/EP2660738A1/en
Priority to CA2775682A priority patent/CA2775682A1/en
Priority to US13/459,336 priority patent/US20130289905A1/en
Publication of JP2013226005A publication Critical patent/JP2013226005A/ja
Application granted granted Critical
Publication of JP5426716B2 publication Critical patent/JP5426716B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/32—Circuit design at the digital level
    • G06F30/33—Design verification, e.g. functional simulation or model checking
    • G06F30/3323—Design verification, e.g. functional simulation or model checking using formal methods, e.g. equivalence checking or property checking
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20—End-user application control systems
    • Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
    • Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Evolutionary Computation (AREA)
  • Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Geometry (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

【課題】配電網電力フロー分析システム及び方法の提供。
【解決手段】本発明の配電網電力フロー分析システム及び方法は、第１関係マトリクスと第２関係マトリクスを利用し、配電網の電力フローを分析できる。該第１関係マトリクスはノード注入電流と分岐電流の関係とされる。第２関係マトリクスはノード電圧差と分岐電流の関係とされる。本発明のシステムと方法はまた、ノード或いはインピーダンス或いは並列回路が新たに増設される状況に応用可能で、其の他の周知の方法と比較すると、本発明のシステムと方法は、配電網電力フローの計算に対して良好な強健性と実行速度を有し、且つ必要なメモリ空間は比較的少ない。
【選択図】図１１

Description

本発明は一種の配電網電力フロー分析システム及び方法に関する。
図１を参照されたい。それは周知のバスの３相インピーダンス模型表示図である。それはバス０からバスｋの３相インピーダンス模型であり、Ｚマトリクスは式（１）のように、４×４マトリクスとして表示される。
Figure 2013226005
さらに式（２）のように、クロンのリダクション法（Ｋｒｏｎ’ｓ Ｒｅｄｕｃｔｉｏｎ）により、３×３マトリクスに減らす。
Figure 2013226005
これにより、バス０からバスｋまでの電圧と電流の関係式はＺマトリクスにより式（３）のように表される。
Figure 2013226005
上の式を電圧差△Ｖで、式（４）のように一般式として表示することができる。
Figure 2013226005
周知の配電網電力フロー分析システム及び方法は、複雑な計算を必要とし、且つ新たなノード或いはインピーダンスが加わる時、インピーダンス（Ｚ）マトリクスの変動が相当に大きく、計算もまた相当に煩雑となり、実行速度が下がり、及び、大きなメモリ空間が必要となり、且つ強健性もまた不良である。
これにより、一種の新規で且つ進歩性を有する配電網電力フロー分析システム及び方法を提供することで、上述の問題を解決する必要がある。
本発明は一種の配電網電力フロー分析システムを提供し、それは、第１関係マトリクス構築装置と第２関係マトリクス構築装置及び配電網電力フロー分析装置を包含する。
該第１関係マトリクス構築装置は第１関係マトリクスを構築するのに用いられ、それはノード（ｎｏｄｅ）注入電流マトリクスと分岐電流マトリクスの関係とされ、そのうちノード注入電流マトリクスは複数の注入電流とされ、分岐電流マトリクスは複数のノード間の分岐電流とされる。
該第２関係マトリクス構築装置は第２関係マトリクスを構築するのに用いられ、それはノード電圧差マトリクスと分岐電流マトリクスの関係とされ、そのうち、ノード電圧差マトリクスは参考点と其の他のノードの電圧差とされる。
該配電網電力フロー分析装置は該第１関係マトリクスと該第２関係マトリクスに基づき、配電網電力フローを分析する。
本発明はまた、配電網電力フロー分析方法を提供する。それは以下のステップを包含する。
（ａ）第１関係マトリクスを構築する。該第１関係マトリクスは、ノード（ｎｏｄｅ）注入電流マトリクスと分岐電流マトリクスの関係とされ、そのうちノード注入電流マトリクスは複数の注入電流とされ、分岐電流マトリクスは複数のノード間の電流とされる。
（ｂ）第２関係マトリクスを構築する。該第２関係マトリクスはノード電圧差マトリクスと分岐電流マトリクスの関係とされ、そのうち、ノード電圧差マトリクスは参考ノードと其の他のノードの電圧差とされる。
（ｃ）該第１関係マトリクスと該第２関係マトリクスに基づき、配電網電力フローを分析する。
本発明のシステム及び方法は、第１関係マトリクスと第２関係マトリクスを利用し、配電網電力フローを分析し、且つ新たにノード或いはインピーダンス或いは並列回路が増設される状況に応用可能で、其の他の周知の方法と比較すると、本発明のシステムと方法は、配電網電力フローの計算に対して良好な強健性と実行速度を有し、且つ必要なメモリ空間は比較的少ない。
周知のノードの３相インピーダンスモデル表示図である。 本発明の実施例の５つのノードを具えた配電網電力フロー分析システムの表示図である。 本発明の第１の新ノード追加実施例の配電網電力フロー分析システムの表示図である。 本発明の第２の新ノード追加実施例の配電網電力フロー分析システムの表示図である。 本発明の第３の新インピーダンス追加実施例の配電網電力フロー分析システムの表示図である。 本発明の第４の新ノードと新インピーダンス追加実施例の配電網電力フロー分析システムの表示図である。 本発明の配電網電力フロー分析方法のフローチャートである。 本発明によるシミュレーションの配電網電力フロー分析システム表示図である。 本発明と周知の方法の実行時間比較図である。 本発明と周知の方法の反復次数比較図である。 本発明の配電網電力フロー分析システムの回路表示図である。
図２を参照されたい。それは本発明の実施例の５つのノードを具えた配電網の表示図である。図７は本発明の配電網電力フロー分析方法のフローチャートである。図１１は本発明の配電網電力フロー分析システムの回路図である。図２、図７及び図１１を組み合わせて、本発明のシステムと方法を説明する。本発明の配電網電力フロー分析システム１１は、第１関係マトリクス構築装置１１０、第２関係マトリクス構築装置１２０及び配電網電力フロー分析装置１５０を包含する。該配電網電力フロー分析装置１５０は本分析方法を使用するプログラマブル装置、たとえばスマートメーターとされ得る。ステップＳ７１を参照されたい。該第１関係マトリクス構築装置１１０は第１関係マトリクスを構築するのに用いられ、該第１関係マトリクスはノード注入電流マトリクスと分岐電流マトリクスの関係とされ、そのうち、ノード注入電流マトリクスは複数の注入電流とされ、分岐電流マトリクスは複数のノード間の電流とされる。
等価電流注入法を利用し、ノードｉの電流はその第ｋ次の反復の注入パワー及び電圧によりその等価電流値を、以下の式（５）のように表示できる。
Figure 2013226005
 そのうち、Ｖi k、Ｉi k はノードｉの第ｋ次反復の電圧及び電流値である。
本発明の実施例において、図２のノード（ｎｏｄｅ）注入電流Ｉと分岐電流Ｂの関係は式（６）のように表示される。
Figure 2013226005
これにより、分岐電流Ｂと注入電流Ｉの関係式は、第１関係マトリクスＢIにより式（７ａ）のように表示される。
Figure 2013226005
式（７ａ）は式（７ｂ）のように一般式で表示され得る。
Figure 2013226005
 そのうち、ノード注入電流マトリクス［Ｉ］は複数の注入電流とされ、分岐電流マトリクス［Ｂ］は複数のノード間の分岐電流とされる。これにより、［ＢI］はノード注入電流マトリクス［Ｉ］と分岐電流マトリクス［Ｂ］の第１関係マトリクスとされ且つ僅かに０と１を包含する上三角マトリクスとされる。
ステップＳ７２を参照されたい。該第２関係マトリクス構築装置１２０は、第２関係マトリクスを構築するのに用いられ、該第２関係マトリクスは、ノード電圧差マトリクスと分岐電流マトリクスの関係とされ、そのうち、ノード電圧差マトリクスは参考ノードと其の他のノードの電圧差とされる。
図２を参照されたい。ノード（ｎｏｄｅ）１、２、３の電圧は式（８）のように書き換えられる。
Figure 2013226005
 そのうち、Ｖiはノードｉの電圧とされ、Ｚijはノードｉ、ｊ間のインピーダンスとされる。
上の式（８）中のＶ3は、式（９）のように改められる。
Figure 2013226005
式（９）より、ノードの電圧差は分岐電流、分岐インピーダンスの関数として表示される。これによりノード電圧差マトリクスは式（１０ａ）のように表示される。
Figure 2013226005
式（１０ａ）の△Ｖは式（１０ｂ）のように一般式で表示され得る。
Figure 2013226005
 そのうち、ノード電圧差マトリクス［△Ｖ］は参考ノード電圧Ｖ0と其の他のノードの電圧差とされ、［ＺV-BC］はノード電圧差マトリクス［△Ｖ］と分岐電流マトリクス［Ｂ］の第２関係マトリクスとされる。該第２関係マトリクスはノード間のインピーダンスとされ、且つ下三角マトリクスとされる。
ステップＳ７５を参照されたい。該配電網電力フロー分析装置１５０は第１関係マトリクス［ＢI］及び該第２関係マトリクス［ＺV-BC］に基づき、配電網電力フローを分析する。これにより、式（７ｂ）を式（１０ｂ）に代入して式（１１）が得られる。
Figure 2013226005
Figure 2013226005
 式（５）及び式（１２）の解を求めて、配電網電力フローを分析できる。
図３を参照されたい。それは本発明の第１のノード（ｎｏｄｅ）追加実施例の配電網電力フロー分析システムの表示図である。図２及び図３を参照されたい。本発明の第１のノード追加実施例は、新ノード６が該参考ノード０に接続される。且つ新インピーダンスＺnewが新ノード６と参考ノード０の間に設けられる。
図３、図７及び図１１を参照されたい。ステップＳ７３において、新たにノード或いはインピーダンスが追加されたかを判断し、イエスであれば、ステップＳ７４において第１関係マトリクスと第２関係マトリクスを修正する。
本発明の配電網電力フロー分析システム１１は、また第１関係マトリクス修正装置１３０及び第２関係マトリクス修正装置１４０を包含し、それらは新たにノード或いはインピーダンスが追加された時、該第１関係マトリクス及び該第２関係マトリクスを修正する。
上述の第１ノード追加実施例において、該第１関係マトリクス修正装置１３０は第１関係マトリクスに一行及び一列を追加し、新たに追加された行列の対角線位置は１とされ、其の他は０とされ、以下の式（１３ａ）のように表示される。
Figure 2013226005
該第２関係マトリクス修正装置１４０は第２関係マトリクスに一行及び一列を追加し、新たに追加された行列の対角線位置は新たに追加されたインピーダンスとされ、其の他は０とされ、以下の式（１３ｂ）のように表示される。
Figure 2013226005
式（１３ｂ）は式（１３ｃ）のように一般式で表示され得る。
Figure 2013226005
式（１３ｃ）に式（１３ａ）を代入して式（１３ｄ）が得られる。
Figure 2013226005
図４を参照されたい。図４は本発明の第２のノード追加実施例の配電網電力フロー分析システムの表示図である。図２を図４に合わせて参照されたい。本発明の第２のノード追加実施例によると、新規のノード６が第ｋ個のノード（本実施例では第５個のノード５）に接続され、且つ新インピーダンスＺnewが新ノード６と第５個のノード５の間に設けられる。
図４及び図１１を合わせて参照されたい。該第１関係マトリクス修正装置１３０は、該第１関係マトリクスに１行と１列を追加し、並びに該第１関係マトリクスの第ｋ行（本実施例では第５行とされる）の値をこの新たに追加された行に複製する。この新たに追加された行列の対角線位置は１とされ、其の他は０とされ、式（１４ａ）のように示される。
Figure 2013226005
 そのうち、ｃｏｌ．（ｋ）はもとの第１関係マトリクスの第ｋ行（本実施例では第５行）の値である。
該第２関係マトリクス修正装置１４０は該第２関係マトリクスに１行及び１列を追加し、並びに該第２関係マトリクスの第ｋ列（本実施例では第５列とされる）の値をこの新たに追加された列に複製する。この新たに追加された行列の対角線位置は新インピーダンスとされ、其の他は０とされ、式（１４ｂ）のように示される。
Figure 2013226005
式（１４ｂ）は式（１４ｃ）のように一般式で表示され得る。
Figure 2013226005
 そのうち、ｒｏｗ．（ｋ）はもとの第２関係マトリクスの第ｋ列（本実施例では第５列）の値とされる。
式（１４ｃ）に式（１４ａ）を代入して式（１４ｄ）が得られる。
Figure 2013226005
図５を参照されたい。それは本発明の第３のインピーダンス追加実施例の配電網電力フロー分析システムの表示図である。図２を図５に合わせて参照されたい。本発明の第３のインピーダンス追加実施例によると、新インピーダンスＺnewが第ｉ個のノード（本実施例では第４個のノード４）と第ｊ個のノード（本実施例では第５個のノード５）の間に増加される。
図６及び図１１を参照されたい。該第１関係マトリクス修正装置１３０は、該第１関係マトリクスに１行及び１列を追加し、並びに該第１関係マトリクスの第ｉ行（本実施例では第４行）の値より第ｊ行（本実施例では第５行）の値を減じて追加された行に置く。この追加された行列の対角線位置は１とされ、其の他は０とされる。本実施例では、新インピーダンスＺnewが第４ノード４と第５ノード５の間に追加され、第４ノード４と第５ノード５の注入電流は、式（１５ａ）のように表示される。
Figure 2013226005
第１関係マトリクスは式（１５ｂ）のように修正される。
Figure 2013226005
Ｂnewを移項すると、式（１５ｂ）は式（１５ｃ）のように表される。
Figure 2013226005
式（１５ｃ）は一般式（１５ｄ）のように表される。
Figure 2013226005
 そのうち、ｃｏｌ．（ｉ−ｊ）はもとの第１関係マトリクスの第ｉ行（本実施例では第４行）の値より第ｊ行（本実施例では第５行）の値を減じたものである。
該第２関係マトリクス修正装置１４０は該第２関係マトリクスに１行と１列を追加し、並びに該第２関係マトリクスの第ｉ列の値より第ｊ列の値を減じてこの追加された列に置き、追加された行列の対角線位置は該新インピーダンスとされ、其の他は０とされる。
ＫＶＬをこの新インピーダンスが追加された回路に応用し、式（１６ａ）のように表示する。
Figure 2013226005
式（１６ａ）と式（１０ａ）を合併し、ノード電圧差マトリクスは以下の式（１６ｂ）のように表される。
Figure 2013226005
式（１６ｂ）は一般式（１６ｃ）のように表示できる。
Figure 2013226005
 そのうち、ｒｏｗ．（ｉ−ｊ）はもとの第２関係マトリクスの第ｉ列（本実施例では第４列）の値より第ｊ列（本実施例では第５列）の値を減じたものである。
式（１５ｃ）に式（１６ｃ）を代入して式（１６ｄ）が得られる。
Figure 2013226005
 新たな列（ｒｏｗ）は、もとのＺの第ｉ列の値より第ｊ列の値を減じたものであり、Ｚnewが対角線位置に置かれる。そのうち、
Figure 2013226005
反復プロセスが式（５）、式（１６ｄ）、式（１６ｅ）を用いて、予め調整された電圧差△ｖにいたるまで実行される。
図６を参照されたい。それは、本発明の第４のノードとインピーダンス追加実施例の配電網電力フロー分析システムの表示図である。図２を図６に合わせて参照されたい。本発明の第４のノードとインピーダンス追加実施例によると、新ノード６が第ｉ個のノード（本実施例では第２個のノード２）と第ｊ個のノード（本実施例では第５個のノード５）の間に追加され、且つ新インピーダンスＺnewが、該新ノード６と該第ｉ個のノード（本実施例では第２個のノード２）の間に設けられる。
図５及び図１１を合わせて参照されたい。該第１関係マトリクス修正装置１３０は、該第１関係マトリクスに２行及び２列を追加し、並びに第１関係マトリクスの第ｉ行（本実施例では第２行）の値を追加された第１行に複製し、且つ該第１関係マトリクスの追加された第１行の値より第ｊ行（本実施例では第５行）の値を減じて追加された第２行に置き、追加された行列の対角線位置は１とされ、其の他は０とされ、式（１７ａ）のように表示される。
Figure 2013226005
式（１７ａ）は一般式（１７ｂ）のように表示され得る。
Figure 2013226005
 そのうち、ｃｏｌ．（ｉ）はもとの第１関係マトリクスの第ｉ行（本実施例では第２行）の値であり、ｃｏｌ．（ｋ−ｊ）は、第１関係マトリクスの追加された第１行の値より第ｊ行（本実施例では第５行）の値を減じたものである。
該第２関係マトリクス修正装置１４０は、該第２関係マトリクスを２行及び２列追加し、並びに該第２関係マトリクスの第ｉ列（本実施例では第２列）の値を追加された第１列に複製し、且つ該第２関係マトリクスの追加された第１列の値より第ｊ列（本実施例では第５列）の値を減じて追加された第２列に置き、追加された第１行及び追加された第１列の対角線位置は新インピーダンスとされ、其の他は０とされる。ＫＶＬをこの追加された回路に応用すると、式（１８ａ）のように表示される。
Figure 2013226005
式（１８ａ）と式（１０ａ）を合併すると、ノード電圧差マトリクスは式（１８ｂ）のように表示される。
Figure 2013226005
式（１８ｂ）は一般式（１８ｃ）のように表示される。
Figure 2013226005
 そのうち、ｒｏｗ．（ｉ）はもとの第２関係マトリクスの第ｉ列（本実施例では第２列）の値である。ｒｏｗ．（ｋ−ｊ）は、第２関係マトリクスの追加された第１列の値より第ｊ列（本実施例では第５列）の値を減じた値である。
式（１７ｂ）を式（１８ｃ）に代入して式（１８ｄ）が得られる。
Figure 2013226005
Figure 2013226005
反復プロセスは、予調整電圧差△Ｖにいたるまで、式（５）、式（１８ｄ）及び式（１８ｅ）により実行される。
式（１８ｄ）は式（１９）のように表示され得る。
Figure 2013226005
クロンのリダクション法（Ｋｒｏｎ’ｓ Ｒｅｄｕｃｔｉｏｎ）を利用し、ノード電圧差マトリクスは以下の式（２０）のように表される。
Figure 2013226005
図８を参照されたい。図８は本発明の配電網システムシミュレーション図である。それは、８個のノードを具えている。本発明は、周知のガウスＺマトリクス法（Ｇａｕｓｓｉｍｐｌｉｃｉｔ Ｚ−ｍａｔｒｉｘ）、周知のニュートン・ラフソン法及び本発明の方法でシミュレーション比較を行った。
図９を参照されたい。図９は本発明と周知の方法の実行時間比較図である。そのうち、第１の周知の方法は周知のガウスＺマトリクス法とされ、第２の周知の方法は、ニュートン・ラフソン法とされる。図９から分かるように、本発明の実行時間（ｎｏｒｍａｌｉｚｅｄ ｅｘｅｃｕｔｉｏｎ ｔｉｍｅ，ＮＥＴ）は周知の方法より遥かに短く、且つノード数が多くなるほど、本発明の方法はより効率的且つスピーディーとなる。
図１０を参照されたい。図１０は本発明と周知の方法の反復回数比較図である。
本発明のシステム及び方法は、該第１関係マトリクス及び該第２関係マトリクスを利用し、配電網電力フローを分析でき、且つ新たにノード或いはインピーダンスが追加される状況に応用可能であり、其の他の周知の方法と比較すると、本発明のシステム及び方法は、配電網電力フローの計算に対して良好な強健性と実行速度を有し、且つ必要なメモリ空間が比較的少ない。
以上述べたことは、本発明の実施例にすぎず、本発明の実施の範囲を限定するものではなく、本発明の特許請求の範囲に基づきなし得る同等の変化と修飾は、いずれも本発明の権利のカバーする範囲内に属するものとする。
１１ 配電網電力フロー分析システム
１１０ 第１関係マトリクス構築装置
１２０ 第２関係マトリクス構築装置
１３０ 第１関係マトリクス修正装置
１４０ 第２関係マトリクス修正装置
１５０ 配電網電力フロー分析装置

Claims (16)

  1. 配電網電力フロー分析システムにおいて、
    第１関係マトリクス構築装置であって、第１関係マトリクスを構築するのに用いられ、該第１関係マトリクスはノード注入電流マトリクスと分岐電流マトリクスの関係とされ、該ノード注入電流マトリクスは複数の注入電流とされ、該分岐電流マトリクスは複数のノード間の分岐電流とされる、上記第１関係マトリクス構築装置と、
    第２関係マトリクス構築装置であって、第２関係マトリクスを構築するのに用いられ、該第２関係マトリクスはノード電圧差マトリクスと分岐電流マトリクスの関係とされ、該ノード電圧差マトリクスは参考ノードと其の他のノードの電圧差とされる、上記第２関係マトリクス構築装置と、
    配電網電力フロー分析装置であって、該第１関係マトリクスと該第２関係マトリクスに基づき、配電網電力フローを分析する、上記配電網電力フロー分析装置と、
    を包含することを特徴とする、配電網電力フロー分析システム。
  2. 請求項１記載の配電網電力フロー分析システムにおいて、該システムの装置は、配電網電力フロー分析方法を使用するプログラマブル装置であるスマートメーターを包含することを特徴とする、配電網電力フロー分析システム。
  3. 請求項１記載の配電網電力フロー分析システムにおいて、該第１関係マトリクスは上三角マトリクスとされ、該第２関係マトリクスは下三角マトリクスとされることを特徴とする、配電網電力フロー分析システム。
  4. 請求項１記載の配電網電力フロー分析システムにおいて、ノード或いはインピーダンスが新たに追加される時に、該第１関係マトリクスと該第２関係マトリクスを修正するのに用いられる第１関係マトリクス修正装置と第２関係マトリクス修正装置をさらに包含することを特徴とする、配電網電力フロー分析システム。
  5. 請求項４記載の配電網電力フロー分析システムにおいて、追加された新ノードが参考ノードに接続され、且つ新インピーダンスが新ノードと該参考ノードの間に接続されると、該第１関係マトリクス修正装置は該第１関係マトリクスに１行及び１列を追加し、この追加された行列の対角線位置を１とし、其の他は０とし、該第２関係マトリクス修正装置は該第２関係マトリクスに１行及び１列を追加し、この追加された行列の対角線位置を該新インピーダンスとし、其の他を０とすることを特徴とする、配電網電力フロー分析システム。
  6. 請求項４記載の配電網電力フロー分析システムにおいて、追加された新ノードが第ｋ個のノードに接続され、且つ新インピーダンスが該新ノードと該第ｋ個のノードの間に接続される時、該第１関係マトリクス修正装置は該第１関係マトリクスに１行及び１列を追加し、並びに該第１関係マトリクスの第ｋ行の値をこの新たに追加された行に複製し、該追加された行の対角線位置を０とし、其の他を０とし、該第２関係マトリクス修正装置は該第２関係マトリクスに１行及び１列を追加し、並びに該第２関係マトリクスの第ｋ列の値をこの新たに追加された列に複製し、この新たに追加された行列の対角線位置を該新インピーダンスとし、其の他を０とすることを特徴とする、配電網電力フロー分析システム。
  7. 請求項４記載の配電網電力フロー分析システムにおいて、新インピーダンスが第ｉ個のノードと第ｊ個のノードの間に追加される時、該第１関係マトリクス修正装置は該第１関係マトリクスに１行及び１列を追加し、並びに該第１関係マトリクスの第ｉ行の値より第ｊ行の値を減じて、この新たに追加された行に置き、この新たに追加された行列の対角線位置を１とし、其の他は０とし、該第２関係マトリクス修正装置は該第２関係マトリクスに１行及び１列を追加し、並びに該第２関係マトリクスの第ｉ列の値をこの新たに追加された列に置き、この新たに追加された行列の対角線位置を該新インピーダンスとし、其の他を０とすることを特徴とする、配電網電力フロー分析システム。
  8. 請求項４記載の配電網電力フロー分析システムにおいて、新ノードが第ｉ個のノードと第ｊ個のノードの間に追加され、且つ新インピーダンスが該新ノードと該ｉ個のノードの間に設けられる時、該第１関係マトリクス修正装置は該第１関係マトリクスに２行及び２列を追加し、並びに該第１関係マトリクスの第ｉ行の値をこの新たに追加された第１行に複製し、且つ該第１関係マトリクスの新たに追加された第１行の値より第ｊ行の値を減じてこの新たに追加された第２行に置き、この新たに追加された行列の対角線位置を１とし、其の他を０とし、該第２関係マトリクス修正装置は、該第２関係マトリクスに２行及び２列を追加し、並びに該第２関係マトリクスの第ｉ列の値をこの新たに追加された列に複製し、且つ該第２関係マトリクスの新たに追加された第１列の値より第ｊ列の値を減じてこの新たに追加された第２列に置き、該第２関係マトリクスの新たに追加された第１行と新たに追加された第１列の対角線位置を該新インピーダンスとし、其の他は０とすることを特徴とする、配電網電力フロー分析システム。
  9. 配電網電力フロー分析方法において、
    （ａ）第１関係マトリクスを構築するステップであって、該第１関係マトリクスは、ノード注入電流マトリクスと分岐電流マトリクスの関係とされ、該ノード注入電流マトリクスは複数の注入電流とされ、該分岐電流マトリクスは複数のノード間の電流とされるステップ、
    （ｂ）第２関係マトリクスを構築するステップであって、該第２関係マトリクスはノード電圧差マトリクスと分岐電流マトリクスの関係とされ、該ノード電圧差マトリクスは参考ノードと其の他のノードの電圧差とされるステップ、
    （ｃ）該第１関係マトリクスと該第２関係マトリクスに基づき、配電網電力フローを分析するステップ、
    以上を包含することを特徴とする、配電網電力フロー分析方法。
  10. 請求項９記載の配電網電力フロー分析方法において、該第１関係マトリクスは上三角マトリクスとされることを特徴とする、配電網電力フロー分析方法。
  11. 請求項９記載の配電網電力フロー分析方法において、該第２関係マトリクスはノード間インピーダンスとされて下三角マトリクスとされることを特徴とする、配電網電力フロー分析方法。
  12. 請求項９記載の配電網電力フロー分析方法において、ノード或いはインピーダンス或いは並列回路が新たに追加される時に、該第１関係マトリクスと該第２関係マトリクスを修正する第１関係マトリクス修正ステップと第２関係マトリクス修正ステップをさらに包含することを特徴とする、配電網電力フロー分析方法。
  13. 請求項１２記載の配電網電力フロー分析方法において、追加された新ノードが参考ノードに接続され、且つ新インピーダンスが新ノードと該参考ノードの間に接続されると、該第１関係マトリクス修正ステップで、装置は該第１関係マトリクスに１行及び１列を追加し、この追加された行列の対角線位置を１とし、其の他は０とし、該第２関係マトリクス修正ステップで、該第２関係マトリクスに１行及び１列を追加し、この追加された行列の対角線位置を該新インピーダンスとし、其の他を０とすることを特徴とする、配電網電力フロー分析方法。
  14. 請求項１２記載の配電網電力フロー分析方法において、追加された新ノードが第ｋ個のノードに接続され、且つ新インピーダンスが該新ノードと該第ｋ個のノードの間に接続される時、該第１関係マトリクス修正ステップで、該第１関係マトリクスに１行及び１列を追加し、並びに該第１関係マトリクスの第ｋ行の値をこの新たに追加された行に複製し、該追加された行の対角線位置を０とし、其の他を０とし、該第２関係マトリクス修正ステップで、該第２関係マトリクスに１行及び１列を追加し、並びに該第２関係マトリクスの第ｋ列の値をこの新たに追加された列に複製し、この新たに追加された行列の対角線位置を該新インピーダンスとし、其の他を０とすることを特徴とする、配電網電力フロー分析方法。
  15. 請求項１２記載の配電網電力フロー分析方法において、新インピーダンスが第ｉ個のノードと第ｊ個のノードの間に追加される時、該第１関係マトリクス修正ステップで、該第１関係マトリクスに１行及び１列を追加し、並びに該第１関係マトリクスの第ｉ行の値より第ｊ行の値を減じて、この新たに追加された行に置き、この新たに追加された行列の対角線位置を１とし、其の他は０とし、該第２関係マトリクス修正ステップで、該第２関係マトリクスに１行及び１列を追加し、並びに該第２関係マトリクスの第ｉ列の値より第ｊ列の値を減じて、この新たに追加された列に置き、この新たに追加された行列の対角線位置を該新インピーダンスとし、其の他を０とすることを特徴とする、配電網電力フロー分析方法。
  16. 請求項１２記載の配電網電力フロー分析方法において、新ノードが第ｉ個のノードと第ｊ個のノードの間に追加され、且つ新インピーダンスが該新ノードと該ｉ個のノードの間に設けられる時、該第１関係マトリクス修正ステップで、該第１関係マトリクスに２行及び２列を追加し、並びに該第１関係マトリクスの第ｉ行の値をこの新たに追加された第１行に複製し、且つ該第１関係マトリクスの新たに追加された第１行の値より第ｊ行の値を減じてこの新たに追加された第２行に置き、この新たに追加された行列の対角線位置を１とし、其の他を０とし、該第２関係マトリクス修正ステップで、該第２関係マトリクスに２行及び２列を追加し、並びに該第２関係マトリクスの第ｉ列の値をこの新たに追加された列に複製し、且つ該第２関係マトリクスの新たに追加された第１列の値より第ｊ列の値を減じてこの新たに追加された第２列に置き、該第２関係マトリクスの新たに追加された第１行と新たに追加された第１列の対角線位置を該新インピーダンスとし、其の他は０とすることを特徴とする、配電網電力フロー分析方法。
JP2012097484A 2012-04-23 2012-04-23 配電網電力フロー分析システム及び方法 Expired - Fee Related JP5426716B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2012097484A JP5426716B2 (ja) 2012-04-23 2012-04-23 配電網電力フロー分析システム及び方法
EP12166226.6A EP2660738A1 (en) 2012-04-23 2012-04-30 Distribution power flow analysis system and method
CA2775682A CA2775682A1 (en) 2012-04-23 2012-04-30 Distribution power flow analyisis system and method
US13/459,336 US20130289905A1 (en) 2012-04-23 2012-04-30 Distribution power flow analysis system and method

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2012097484A JP5426716B2 (ja) 2012-04-23 2012-04-23 配電網電力フロー分析システム及び方法
EP12166226.6A EP2660738A1 (en) 2012-04-23 2012-04-30 Distribution power flow analysis system and method
CA2775682A CA2775682A1 (en) 2012-04-23 2012-04-30 Distribution power flow analyisis system and method
US13/459,336 US20130289905A1 (en) 2012-04-23 2012-04-30 Distribution power flow analysis system and method

Publications (2)

Publication Number Publication Date
JP2013226005A true JP2013226005A (ja) 2013-10-31
JP5426716B2 JP5426716B2 (ja) 2014-02-26

Family

ID=67146603

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012097484A Expired - Fee Related JP5426716B2 (ja) 2012-04-23 2012-04-23 配電網電力フロー分析システム及び方法

Country Status (4)

Country Link
US (1) US20130289905A1 (ja)
EP (1) EP2660738A1 (ja)
JP (1) JP5426716B2 (ja)
CA (1) CA2775682A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104484234A (zh) * 2014-11-21 2015-04-01 中国电力科学研究院 一种基于gpu的多波前潮流计算方法和系统
CN108292845A (zh) * 2015-10-16 2018-07-17 超级电力研究所 用于高压直流电网的互连设备
CN110991816A (zh) * 2019-11-13 2020-04-10 广州供电局有限公司 一流配电网的建设水平监测方法和装置
CN115995853A (zh) * 2023-02-21 2023-04-21 国网山东省电力公司临沂供电公司 一种基于功率控制的交直流配电系统过程管理方法及系统

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104092210B (zh) * 2014-07-14 2016-09-21 国家电网公司 一种具备n-2开断快速扫描功能的在线静态安全分析方法
CN104899396B (zh) * 2015-06-19 2017-08-08 大连海事大学 一种修正系数矩阵的快速分解法潮流计算方法
CN106130025B (zh) * 2016-06-24 2018-05-08 西南石油大学 一种基于回路分析法计算三相配电网潮流的线性方法
CN106451456B (zh) * 2016-08-30 2019-01-29 邓宏伟 一种基于环网电力系统的直算方法
CN106451556B (zh) * 2016-11-24 2021-10-12 中国电力科学研究院有限公司 一种确定配电网馈线接纳分布式电源能力的方法和装置
CN107665184B (zh) * 2017-11-16 2020-05-22 大连海事大学 基于关联矩阵运算的潮流计算导纳矩阵计算方法
CN107944682B (zh) * 2017-11-16 2020-04-07 大连海事大学 基于Matlab矩阵运算的潮流计算导纳矩阵计算方法
CN109980651B (zh) * 2019-05-21 2022-11-15 华北水利水电大学 一种环网型配电系统潮流计算方法
US11169187B2 (en) 2019-06-28 2021-11-09 King Fahd University Of Petroleum And Minerals Zig zag based load flow method and system for extended radial distribution systems
CN110460059B (zh) * 2019-07-26 2023-02-03 国网山东省电力公司泰安供电公司 一种电力系统潮流计算方法
CN110544941B (zh) * 2019-09-05 2022-03-29 国网黑龙江省电力有限公司电力科学研究院 一种电气节点源流能量分配关系计算方法
CN112364586B (zh) * 2020-11-11 2024-05-14 南方电网科学研究院有限责任公司 直流换流器的直流回路阻抗计算方法、装置和计算机设备

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5568399A (en) * 1995-01-31 1996-10-22 Puget Consultants Inc. Method and apparatus for power outage determination using distribution system information
US8024051B2 (en) * 2009-02-24 2011-09-20 Oracle America, Inc. Parallel power grid analysis
US8810253B2 (en) * 2010-09-20 2014-08-19 Thomas L. Marzetta Characterization of electrical power distribution systems using characterization matrices
US8712711B2 (en) * 2011-07-21 2014-04-29 Cisco Technology, Inc. Identification of electrical grid phase information for end-points in a grid network

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104484234A (zh) * 2014-11-21 2015-04-01 中国电力科学研究院 一种基于gpu的多波前潮流计算方法和系统
CN108292845A (zh) * 2015-10-16 2018-07-17 超级电力研究所 用于高压直流电网的互连设备
CN108292845B (zh) * 2015-10-16 2021-10-22 超级电力研究所 用于高压直流电网的互连设备
CN110991816A (zh) * 2019-11-13 2020-04-10 广州供电局有限公司 一流配电网的建设水平监测方法和装置
CN115995853A (zh) * 2023-02-21 2023-04-21 国网山东省电力公司临沂供电公司 一种基于功率控制的交直流配电系统过程管理方法及系统
CN115995853B (zh) * 2023-02-21 2023-08-18 国网山东省电力公司临沂供电公司 一种基于功率控制的交直流配电系统过程管理方法及系统

Also Published As

Publication number Publication date
EP2660738A1 (en) 2013-11-06
CA2775682A1 (en) 2013-10-30
JP5426716B2 (ja) 2014-02-26
US20130289905A1 (en) 2013-10-31

Similar Documents

Publication Publication Date Title
JP5426716B2 (ja) 配電網電力フロー分析システム及び方法
CN104049175B (zh) 具有多个多相母线的电网中的故障分析
Li et al. Distribution grid impedance & topology estimation with limited or no micro-PMUs
CN111181164B (zh) 一种改进的主从分裂输配协同潮流计算方法及系统
Lisboa et al. A fast power flow method for radial networks with linear storage and no matrix inversions
Zhang et al. A projection iterative algorithm boundary element method for the Signorini problem
Grisales-Noreña et al. Power flow methods used in AC distribution networks: An analysis of convergence and processing times in radial and meshed grid configurations
Murari et al. Graph‐theoretic based approach for the load‐flow solution of three‐phase distribution network in the presence of distributed generations
Yao et al. Contingency analysis based on partitioned and parallel holomorphic embedding
US20090132975A1 (en) Circuit Splitting in Analysis of Circuits at Transistor Level
Li et al. A loop-analysis theory based power flow method and its linear formulation for low-voltage DC grid
Zhang Improved GNN method with finite-time convergence for time-varying Lyapunov equation
US8645883B2 (en) Integrated circuit simulation using fundamental and derivative circuit runs
WO2018188393A1 (zh) 柔性直流电网潮流线性化确定方法及系统、计算机存储介质
CN114204564A (zh) 一种含逆变型新能源电网的短路电流计算方法及装置
CN105447214A (zh) 器件参数的确定方法和装置
CN111339624B (zh) 基于psasp和emtp/atp短路电流直流分量计算方法
Linan et al. Energy and helicity fluxes in line-tied eruptive simulations
CN106529089B (zh) 用于含小阻抗支路电网的补偿法快速分解法潮流计算方法
CN115130415A (zh) 供电网络瞬态仿真方法、装置、设备和介质
CN105260317B (zh) 一种测试用例的选取方法
Paszyńska Graph-grammar greedy algorithm for reutilization of partial LU factorization over 3D tetrahedral grids
CN107609274B (zh) 基于传输线与级别调度法的二维静磁场并行有限元方法
TWI421705B (zh) 分散式電力潮流分析系統及方法
Babu et al. An efficient power flow method for distribution system studies under various load models

Legal Events

Date Code Title Description
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: 20131119

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131128

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees