CN108899898A - A kind of reactive voltage partition method - Google Patents

A kind of reactive voltage partition method Download PDF

Info

Publication number
CN108899898A
CN108899898A CN201810764864.9A CN201810764864A CN108899898A CN 108899898 A CN108899898 A CN 108899898A CN 201810764864 A CN201810764864 A CN 201810764864A CN 108899898 A CN108899898 A CN 108899898A
Authority
CN
China
Prior art keywords
subregion
reactive
node
idle
oriented
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
CN201810764864.9A
Other languages
Chinese (zh)
Other versions
CN108899898B (en
Inventor
徐方维
舒勤
刘文通
钟俊
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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN201810764864.9A priority Critical patent/CN108899898B/en
Publication of CN108899898A publication Critical patent/CN108899898A/en
Application granted granted Critical
Publication of CN108899898B publication Critical patent/CN108899898B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30—Reactive power compensation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of reactive voltage partition methods, including：A, L-Q sensitivity matrix when network system is in critical state is calculated, and electrical distance is defined with this, suitable between class distance binding hierarchy clustering method is chosen and subregion internal layer is obtained to load bus progress subregion, and determine subregion number；B, the System Reactive Power trend under computing system a period of time discontinuity surface, the oriented idle topological diagram of weighting of building system, carries out subregion to topological network using the community discovery method based on side betweenness and obtains subregion outer layer；C, map overlay principle is formulated based on reactive reserve figureofmerit, and map overlay is carried out to subregion internal layer and subregion outer layer, obtain final division result.Its idle partitioned organization of gained is reasonable, it is ensured that has enough reactive reserves in subregion, achieves preferable subregion effect.

Description

A kind of reactive voltage partition method
Technical field
The present invention relates to electric network reactive-load voltage partitioning techniques fields, and in particular to a kind of reactive voltage partition method.
Background technique
How to carry out rationally effective reactive voltage subregion to power grid is the key that one of reactive power/voltage control scheme.With Power grid more bulky complex, also more difficult to the analysis monitoring during operation of power networks, traditional sub-area division is dependent on ground Region and personnel's operating experience are managed, can not reflect the actual operation conditions of power grid.In conjunction with traditional analysis and modern complexity Science carries out subregion to electric system, constructs reasonable electric network composition and is of great significance to safe operation of power system.
Currently, existing reactive voltage partition method is Voltage-Reactive Power sensitivity matrix by system, PQ decomposition method mostly The subregions such as B " matrix the electrical distance binding hierarchy clustering procedure, graph theoretic approach, simulated annealing and the Complex Networks Theory that define calculate Method is divided.
But the existing definition to electrical distance is only unilateral to the progress of the variables such as node voltage, impedance, power, electric current Analysis considers, can not reflect energy transmission size and distance between node simultaneously, in concentrated expression actual electric network operational process Shortcomings on electrical link.
In recent years, become one of the hot spot of current electric power research to power system network parsing by Complex Networks Theory, Now research shows that modern power systems have small world and uncalibrated visual servo characteristic, which determines single node failure in power grid System interlock will be caused to have a power failure, be suitable for solving the problems, such as sub-area division.But in complexity analysis network theory, with global search Mode computing module degree and then when determining the optimal network number of partitions, needs higher time complexity and space complexity, and Route weight defined in calculating process can not electrical link between concentrated expression node.
Summary of the invention
In order to solve the above-mentioned technical problem the present invention provides a kind of reactive voltage partition method.
The present invention is achieved through the following technical solutions：
A kind of reactive voltage partition method, including：
A, L-Q sensitivity matrix when network system is in critical state is calculated, and electrical distance is defined with this, by sentencing Disconnected comparison chooses between class distance binding hierarchy clustering method and obtains subregion internal layer to load bus progress subregion, and determines subregion Number；
B, the System Reactive Power trend under computing system a period of time discontinuity surface, the oriented idle topological diagram of weighting of building system, uses Community discovery method based on side betweenness carries out subregion to topological network and obtains subregion outer layer；
C, map overlay principle is formulated based on reactive reserve figureofmerit, and it is folded to carry out figure layer to subregion internal layer and subregion outer layer Add, obtains final division result.
This programme defines new electrical distance according to local voltage stability index, is constructed using Hierarchical Cluster Analysis Partition load internal layer；The idle transport network topology that Complex Networks Theory is constructed with the community discovery method based on side betweenness again Figure carries out subregion, constructs subregion outer layer；Finally, carrying out map overlay with corresponding division principle, sub-area division is realized.According to office Portion's voltage stability index proposes new electrical distance definition, and this definition explicit physical meaning more conducively determines number of partitions, adds The Decoupling Characteristics in strong region.The aeoplotropism of reactive power flow is considered during subregion, the oriented weighted topological figure of building is examined The line impedance influence idle to route is considered, has met modeling needs, reduce the calculation amount during subregion.It is folded by figure layer The method added compensates for corresponding to ectonexine subregion can not be to guaranteeing reactive balance and can not determine the defect of number of partitions, gained Idle partitioned organization is reasonable, it is ensured that has enough reactive reserves in subregion, achieves preferable subregion effect.
It can realize that the principle of the electrical link between reflection node is as follows using this programme：
Demand perfection differential to the L index of j node, the independent variable in L index have node active-power P, node reactive power Q, Node voltage amplitude V and node voltage phase place θ, wherein P=[P1,P2,L,Pn] T, Q=[Q1,Q2,L,Qn]T, V=[V1,V2,L, Vn]T, θ=[θ1,θ2,L,θn]T, can obtain its total differential equation is：
From the above equation, we can see that L index can comprehensive characterization node physical characteristic, the influence situation of L index can be indirect between each node Reflect the electrical link between node.
Load bus i is with electrical distance of the load bus j in n-dimensional space：
Wherein, dijElectrical distance matrix element between node i and node j, dij=-lg | αijαji|；αijFor system L Beacon decay matrix element,Node j is idle injection node, and node i is responsive node, The variation of node reactive power and the variation of L index are respectively indicated, the L index variation on j is
The oriented idle topological diagram of weighting of building system includes calculating the oriented most short electrical path of weighting
The oriented most short electrical path of weightingCalculation method be：
Wherein,For node viTo node vjOriented electrical path line set, XijFor line reactance, QijWhen being current Between trend section line reactive power, wij=Xij/QijFor route weight coefficient.Oriented weighting is calculated most using the above method Short electrical pathThe electrical path calculates the reactive power flow between node, and practical object is had more for line impedance Meaning is managed, the calculating of the electrical path is and then to obtain subregion internal layer to be applied to the community discovery method based on side betweenness.
The reactive reserve figureofmeritCalculation method be：
Wherein, QGqFor the idle power output of maximum of region q, QLqFor total load or burden without work in the q of region.
Compared with prior art, the present invention having the following advantages and benefits：
1, the present invention proposes new electrical distance according to local voltage stability index and defines, and this definition physical significance is bright Really, number of partitions, the Decoupling Characteristics of reinforced region are more conducively determined.The oriented of reactive power flow is considered during subregion Property, the oriented weighted topological figure of building considers the line impedance influence idle to route, meets modeling needs, reduces subregion Calculation amount in the process.Being compensated for by the method for map overlay can not be to guarantee reactive balance and nothing corresponding to ectonexine subregion Method determines the defect of number of partitions, and gained is idle, and partitioned organization is reasonable, it is ensured that has enough reactive reserves in subregion, achieves Preferable subregion effect.
Detailed description of the invention
Attached drawing described herein is used to provide to further understand the embodiment of the present invention, constitutes one of the application Point, do not constitute the restriction to the embodiment of the present invention.
Fig. 1 is the flow chart of this method.
Fig. 2 is IEEE-39 node system subregion internal layer schematic diagram.
Fig. 3 is IEEE-39 node system subregion outer layer schematic diagram.
Fig. 4 is 6 subregion schematic diagram of IEEE-39 node system.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below with reference to embodiment and attached drawing, to this Invention is described in further detail, and exemplary embodiment of the invention and its explanation for explaining only the invention, are not made For limitation of the invention.
Embodiment 1
A kind of reactive voltage partition method as shown in Figure 1, including：
A, L-Q sensitivity matrix when network system is in critical state is calculated, and electrical distance is defined with this, by sentencing Disconnected comparison chooses between class distance binding hierarchy clustering method and obtains subregion internal layer to load bus progress subregion, and determines subregion Number；
B, the System Reactive Power trend under computing system a period of time discontinuity surface, the oriented idle topological diagram of weighting of building system, uses Community discovery method based on side betweenness carries out subregion to topological network and obtains subregion outer layer；
C, map overlay principle is formulated based on reactive reserve figureofmerit, and it is folded to carry out figure layer to subregion internal layer and subregion outer layer Add, obtains final division result.
Embodiment 2
Principle based on the above embodiment, the present embodiment disclose a specific embodiment.
The present embodiment is described in detail this programme by taking IEEE-39 standard nodes system as an example.
L-Q sensitivity matrix when network system is in critical state is calculated, and electrical distance is defined with this, specifically, Electrical distance calculation method of the load bus i and load bus j in n-dimensional space be：
Wherein, dijElectrical distance matrix element between node i and node j, dij=-lg | αijαji|；αijFor system L Beacon decay matrix element,Node j is idle injection node, and node i is responsive node, The variation of node reactive power and the variation of L index are respectively indicated, the L index variation on j is
Electrical distance between computing system load bus carries out cluster point to load bus by Hierarchical clustering analysis algorithm Area obtains overloaded partition internal layer, and subregion internal layer schematic diagram is as shown in Figure 2.
The oriented idle topological diagram of weighting of building system includes calculating the oriented most short electrical path of weightingBased on side betweenness Community discovery method to obtain subregion internal layer as shown in Figure 3.
The oriented most short electrical path of weightingCalculation method be：
Wherein,For node viTo node vjOriented electrical path line set, XijFor line reactance, QijWhen being current Between trend section line reactive power, wij=Xij/QijFor route weight coefficient.
Calculate reactive reserve figureofmerit
Wherein, QGqFor the idle power output of maximum of region q, QLqFor total load or burden without work in the q of region.
Map overlay principle is formulated based on reactive reserve figureofmerit, and it is folded to carry out figure layer to subregion internal layer and subregion outer layer Add, obtains final division result and be detailed in Fig. 4.
In additive process, route 26-27,17-27,17-16 carrying reactive power flow be 67.1, -82.8,7.3, -41.6, 30.2, the reactive power flow of -40.8Mvar, critical circuits 17-27 carrying and the reactive power exchange of adjacent lines 26-27,17-16 carrying Gap is larger, carries out dividing the requirement for meeting interregional weak coupling to it, subregion 1 arrives subregion 6 in obtained division result Reactive reserve is followed successively by 53.4%, 21.6%, 43.6%, 30.0%, 64.6%, 20.7%, meets the requirement of reactive reserve.
Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects It is described in detail, it should be understood that being not intended to limit the present invention the foregoing is merely a specific embodiment of the invention Protection scope, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should all include Within protection scope of the present invention.

Claims (5)

1. a kind of reactive voltage partition method, which is characterized in that including：
A, L-Q sensitivity matrix when network system is in critical state is calculated, and electrical distance is defined with this, passes through judgement pair Subregion internal layer is obtained to load bus progress subregion than choosing between class distance binding hierarchy clustering method, and determines subregion number；
B, the System Reactive Power trend under computing system a period of time discontinuity surface, the idle topological diagram of the oriented weighting of building system, using being based on The community discovery method of side betweenness carries out subregion to topological network and obtains subregion outer layer；
C, map overlay principle is formulated based on reactive reserve figureofmerit, and map overlay is carried out to subregion internal layer and subregion outer layer, Obtain final division result.
2. a kind of reactive voltage partition method according to claim 1, which is characterized in that
Load bus i is with electrical distance of the load bus j in n-dimensional space：
Wherein, dijElectrical distance matrix element between node i and node j, dij=-lg | αijαji|；αijFor system L index Damping matrix element,Wherein load bus j is idle injection node, and load bus i is response section Point,The variation of node reactive power and the variation of L index are respectively indicated, the L index variation on j is
3. a kind of reactive voltage partition method according to claim 1, which is characterized in that the oriented weighting of building system is idle Topological diagram includes calculating the oriented most short electrical path of weighting
4. a kind of reactive voltage partition method according to claim 3, which is characterized in that the oriented weighting is most short electrical PathCalculation method be：
Wherein,For node viTo node vjOriented electrical path line set, XijFor line reactance, QijFor current time tide The line reactive power of flow section, wij=Xij/QijFor route weight coefficient.
5. a kind of reactive voltage partition method according to claim 1, which is characterized in that the reactive reserve figureofmerit Calculation method be：
Wherein, QGqFor the idle power output of maximum of region q, QLqFor total load or burden without work in the q of region.
CN201810764864.9A 2018-07-12 2018-07-12 A kind of reactive voltage partition method Active CN108899898B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810764864.9A CN108899898B (en) 2018-07-12 2018-07-12 A kind of reactive voltage partition method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810764864.9A CN108899898B (en) 2018-07-12 2018-07-12 A kind of reactive voltage partition method

Publications (2)

Publication Number Publication Date
CN108899898A true CN108899898A (en) 2018-11-27
CN108899898B CN108899898B (en) 2019-05-17

Family

ID=64349209

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810764864.9A Active CN108899898B (en) 2018-07-12 2018-07-12 A kind of reactive voltage partition method

Country Status (1)

Country Link
CN (1) CN108899898B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110189230A (en) * 2019-01-02 2019-08-30 国网冀北电力有限公司秦皇岛供电公司 A kind of construction method of the parsingization model of dynamic partition
CN111047453A (en) * 2019-12-04 2020-04-21 兰州交通大学 Detection method and device for decomposing large-scale social network community based on high-order tensor
CN111211577A (en) * 2019-11-22 2020-05-29 国网河南省电力公司电力科学研究院 Method for improving operation reliability of new energy in regional power grid subareas
CN111884220A (en) * 2020-07-31 2020-11-03 三峡大学 Electromagnetic looped network open-loop partition scheme generation method considering regional source-load balance
CN113224770A (en) * 2021-06-17 2021-08-06 华北电力大学（保定） Power grid reactive power partitioning method integrating electrical distance and edge betweenness
CN114512989A (en) * 2022-03-30 2022-05-17 国网浙江省电力有限公司湖州供电公司 Power grid partitioning method based on HEM sensitivity

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105186525A (en) * 2015-10-29 2015-12-23 山东大学 Reactive voltage control partitioning method under wind power integration
CN106451492A (en) * 2016-11-11 2017-02-22 华北电力大学 Partitioning method for reactive voltage of power distribution network by considering reactive power source response characteristics
CN107039984A (en) * 2017-06-19 2017-08-11 国网重庆市电力公司电力科学研究院 A kind of reactive voltage partition method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105186525A (en) * 2015-10-29 2015-12-23 山东大学 Reactive voltage control partitioning method under wind power integration
CN106451492A (en) * 2016-11-11 2017-02-22 华北电力大学 Partitioning method for reactive voltage of power distribution network by considering reactive power source response characteristics
CN107039984A (en) * 2017-06-19 2017-08-11 国网重庆市电力公司电力科学研究院 A kind of reactive voltage partition method

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110189230A (en) * 2019-01-02 2019-08-30 国网冀北电力有限公司秦皇岛供电公司 A kind of construction method of the parsingization model of dynamic partition
CN110189230B (en) * 2019-01-02 2023-06-16 国网冀北电力有限公司秦皇岛供电公司 Construction method of analytic model of dynamic partition
CN111211577A (en) * 2019-11-22 2020-05-29 国网河南省电力公司电力科学研究院 Method for improving operation reliability of new energy in regional power grid subareas
CN111211577B (en) * 2019-11-22 2022-04-05 国网河南省电力公司电力科学研究院 Method for improving operation reliability of new energy in regional power grid subareas
CN111047453A (en) * 2019-12-04 2020-04-21 兰州交通大学 Detection method and device for decomposing large-scale social network community based on high-order tensor
CN111884220A (en) * 2020-07-31 2020-11-03 三峡大学 Electromagnetic looped network open-loop partition scheme generation method considering regional source-load balance
CN111884220B (en) * 2020-07-31 2021-06-18 三峡大学 Electromagnetic looped network open-loop partition scheme generation method considering regional source-load balance
CN113224770A (en) * 2021-06-17 2021-08-06 华北电力大学（保定） Power grid reactive power partitioning method integrating electrical distance and edge betweenness
CN114512989A (en) * 2022-03-30 2022-05-17 国网浙江省电力有限公司湖州供电公司 Power grid partitioning method based on HEM sensitivity

Also Published As

Publication number Publication date
CN108899898B (en) 2019-05-17

Similar Documents

Publication Publication Date Title
CN108899898B (en) A kind of reactive voltage partition method
CN1327587C (en) Voltage control method based on soft partitions in power system
CN103902814B (en) Operation states of electric power system detection method based on dynamic partition
Xu et al. Structural and hierarchical partitioning of virtual microgrids in power distribution network
CN109149586B (en) Active power distribution network partition distributed voltage control method oriented to intelligent soft switch
CN105656198B (en) A kind of acquisition methods of powerline network redundant path strategy
CN102354981B (en) Distributed computation based voltage stability assessment method of sub-networks in interconnected power network
CN111817286B (en) Detection method for false data injection attack of direct current micro-grid cluster
CN107453928B (en) A kind of power communication network node different degree evaluation method and device
CN110783967A (en) Constraint aggregation-based virtual power plant output feasible region identification method and device
CN106208102A (en) A kind of active distribution network new distributed optimal reactive power based on Auxiliary Problem Principle
CN105184669A (en) 220kV urban ring network partitioning method based on node set GN splitting-up algorithm
CN103077480A (en) Safety checking method for power system
CN109100614A (en) A kind of transmission open acess system and method based on PMU device
CN104393590A (en) Intelligent identification method for transient instability mode of power grid
CN106897942A (en) A kind of power distribution network distributed parallel method for estimating state and device
CN108376983A (en) A kind of Power Network Status Estimation system and method based on boundary fusion
CN105656036B (en) Consider trend and the probability static security analysis method of sensitivity uniformity equivalence
Huang et al. Power flow calculation of distribution network with multiple energy routers
CN108667027A (en) Ac and dc systems power flow transfer search containing flexible direct current and quantitative analysis method
CN110707693B (en) Aggregate Kalman filtering dynamic state estimation method based on AMI full measurement point partition
Feng et al. Graph computing based distributed parallel power flow for AC/DC systems with improved initial estimate
CN110311372A (en) Sub-area division method based on spectral clustering
CN112103950B (en) Power grid partitioning method based on improved GN splitting algorithm
Sayeekumar et al. Graph theory and its applications in power systems-a review

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant