CN108493921A - Thermal stability security domain fast construction method based on key node injecting power space - Google Patents

Thermal stability security domain fast construction method based on key node injecting power space Download PDF

Info

Publication number
CN108493921A
CN108493921A CN201810133357.5A CN201810133357A CN108493921A CN 108493921 A CN108493921 A CN 108493921A CN 201810133357 A CN201810133357 A CN 201810133357A CN 108493921 A CN108493921 A CN 108493921A
Authority
CN
China
Prior art keywords
node
circuit
security domain
key node
thermal stability
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
CN201810133357.5A
Other languages
Chinese (zh)
Other versions
CN108493921B (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.)
Tianjin University
China Southern Power Grid Co Ltd
Original Assignee
Tianjin University
China Southern Power Grid 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 Tianjin University, China Southern Power Grid Co Ltd filed Critical Tianjin University
Publication of CN108493921A publication Critical patent/CN108493921A/en
Application granted granted Critical
Publication of CN108493921B publication Critical patent/CN108493921B/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/01—Arrangements for reducing harmonics or ripples
    • 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48—Controlling the sharing of the in-phase component
    • 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/40—Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of thermal stability security domain fast construction methods based on key node injecting power space, including：State variable is set；Determine control variable；Build thermal stability security domain.The present invention proposes the concept of key node, being defined on the thermal stability security domain dimensionality reduction of injecting power spatially at only by the simplification security domain of key node spatial description；Hyperplane characteristic based on security domain boundaries is proposed a kind of quick calculation method of safe domain coefficient according to the trend constraint absorbing boundary equation of circuit.Thus obtained thermal stability security domain construction method overcomes traditional simulation fitting process and calculates time length, is not suitable for the shortcomings that application on site, can effectively improve its effect in the big system of electric power in LINE REAL TIME MONITORING, defence and control.Safety margin information and real-time control measure of the system in current point of operation can further be obtained using method proposed by the present invention, directive function is played to operation of power networks scheduling, there is very strong practical value.

Description

Thermal stability security domain fast construction method based on key node injecting power space
Technical field
The present invention relates to a kind of fast construction methods of thermal stability security domain, and in particular to one kind is injected based on key node The thermal stability security domain fast construction method of power space.
Background technology
In recent years, large-scale blackout continuously emerges caused by being destroyed due to security and stability in the big system of electric power, shows There is an urgent need to strong means for electric system to realize its online actual time safety monitoring, defence and control.For a long time, electric power The analysis method of security of system and stability typically presses given scenario or the method for operation in one or more of failure modes Under, the conclusion that system is safe or dangerous, stable or unstable is obtained by simulation calculation.Such methods are known as point by point method, although It is still played an important role in Power System Analysis, but is difficult to propose the overall evaluation to the operating status of electric system. The method of security domain is the new method that grows up on the basis of point by point method, it considers a problem from the angle in domain, description Be it is whole can safe and stable operation region, it is possible to provide safety margin and optimum control information can make electric system real online When regulation measure it is more scientific, more effective.At present in engineer application, traditional thermal stability security domain absorbing boundary equation is usually by counting After value simulation calculation obtains the largely Critical operating points that are suitably distributed, then formed by least square fitting mathematically.This Although kind of a method can describe security domain boundaries and precision is higher, since it needs a large amount of numerical simulations to calculate, to Bring heavy computation burden, spent time is longer, can only off-line calculation use online.
Invention content
The technical problem to be solved in the present invention is：Absorbing boundary equation and security domain boundaries based on Line Flow constraint surpass Flatness of the response defines the key node of injecting power spatially, provides a kind of heat based on key node injecting power spatially Stablize security domain fast construction method, overcoming security domain emulation approximating method to calculate, the time is long, is not suitable for lacking for application on site Point improves practical value of the thermal stability security domain in the big system of electric power in LINE REAL TIME MONITORING, defence and control.
In order to achieve the above-mentioned object of the invention, a kind of thermostabilization based on key node injecting power space proposed by the present invention Security domain fast construction method, it is characterized in that：Include the following steps：
Step 1, setting state variable：
It screens in actual electric network there are effective power flow is out-of-limit or the circuit of heavily loaded problem, setting is all, and there are effective power flows to get over The circuit of limit or heavily loaded problem is state variable set F=[circuit 1, circuit 2 ..., circuit k ..., circuit m], first in set F The total number of element is m；
Step 2 determines control variable, including：
1. calculating effective power flow sensitivity of all node active power outputs variations to circuit in set F, base in actual electric network In DC flow model, effective power flow sensitivity S of the node i to circuit kk-iCalculation formula is as follows
In formula (1), a, b are two endpoints of circuit k, XaiFor the mutual impedance between node a and node i, XbiFor node b Mutual impedance between node i, xkFor the reactance of circuit k；
The absolute value of sensitivity is bigger to indicate that node active power output variation is bigger on the influence of the effective power flow of circuit, thus Key node of the node of sensitivity order of magnitude sequence preceding ten as circuit k is selected, set A is denoted as1；
2. calculating tunability of all nodes to circuit in set F, tunability P of the node i to circuit kadjIt calculates public Formula is as follows：
In formula (2), Pi.min、Pi.max、Pi.nowThe active power lower limit of node i, the active power upper limit are indicated respectively and are worked as Preceding active injection power；
Tunability is bigger to indicate that the ability of node mitigation Line Flow load is bigger, thus selects tunability size Key node of the node of sequence preceding ten as circuit k, is denoted as set A2；
3. taking set A1With set A2Union as key node set i.e. control variables collection, A=A1∪A2；
Step 3, structure thermal stability security domain：
According to the state variable set F set in the control variables collection A and step 1 determined in step 2, to build Dimensionality reduction thermal stability security domain based on key node injecting power space：
In formula (3), nAFor the dimension in key node injecting power space, αk-iIt is that i-th of key node is safe at k-th The hyperplane coefficient on domain boundary, PiFor the injecting power of i-th of key node, ckIndicate k-th of safety of remaining non-key node pair The influence coefficient of domain equation；
Wherein, hyperplane factor alphak-iCalculation formula be：
In formula (4), Uk.maxAnd Ik.maxThe respectively reference voltage and rated current of circuit k；
Influence coefficient ckCalculation formula be：
In formula (5), Pk.nowAnd Pk.maxThe current active power and active power limit of circuit k, P are indicated respectivelyi.nowTo work as The active power injection rate of preceding method of operation lower node i.
The present invention carries on the basis of traditional thermal stability security domain is in node injecting power spatial model according to application scenario The concept for having gone out key node, be defined on the thermal stability security domain dimensionality reduction of injecting power spatially at only by key node vector The simplification security domain of description；Hyperplane characteristic based on security domain boundaries is released according to the trend constraint absorbing boundary equation of circuit and is pacified The quick calculation method of gamut boundary hyperplane coefficient.The thus obtained thermostabilization based on key node injecting power spatially Security domain overcomes security domain emulation approximating method and calculates time length, is not suitable for the shortcomings that application on site, can effectively improve Effect of the thermal stability security domain in the big system of electric power in LINE REAL TIME MONITORING, defence and control.Using side proposed by the present invention Method can further obtain safety margin information and real-time control measure of the system in current point of operation, be dispatched to operation of power networks To directive function, there is very strong practical value.
Specific implementation mode
A kind of thermal stability security domain fast construction method based on key node injecting power space proposed by the present invention, Mentality of designing is to define the dimensionality reduction thermal stability security domain based on key node injecting power space, is indicated as follows
In above formula, nAFor the dimension of key node Injection Space, m is the sum of security domain boundaries hyperplane, αk-iIt is i-th A key node is in the hyperplane coefficient of k-th of security domain boundaries, PiFor the injecting power of i-th of key node, ckIndicate remaining The influence coefficient of k-th of safe domain equation of non-key node pair.
Based on defined above, this method is as follows：
Step 1, setting state variable：
There are the equipment such as the circuits of thermostabilizations such as effective power flow is out-of-limit or heavily loaded as research in screening actual electric network State variable set F=[circuit 1, circuit 2 ..., circuit k ..., circuit m] is arranged in object by taking circuit as an example, element in set F Total number be m；
Step 2 determines control variable：
Alleged control variable refers to the section being affected to the effective power flow of circuit in state variable set F in the present invention Point, is defined as key node.The determination method for controlling variable, that is, key node is as follows
1. calculating effective power flow sensitivity of all node active power outputs variations to circuit in set F, base in actual electric network In DC flow model, effective power flow sensitivity S of the node i to circuit kk-iCalculation formula is as follows
In formula (1), a, b are two endpoints of circuit k, XaiFor the mutual impedance between node a and node i, XbiFor node b Mutual impedance between node i, xkFor the reactance of circuit k；
The absolute value of sensitivity is bigger to indicate that node active power output variation is bigger on the influence of the effective power flow of circuit, thus Key node of the node of sensitivity order of magnitude sequence preceding ten as circuit k is selected, set A is denoted as1；
2. calculating tunability of all nodes to circuit in set F, tunability P of the node i to circuit kadjIt calculates public Formula is as follows：
In formula (2), Pi.min、Pi.max、Pi.nowThe active power lower limit of node i, the active power upper limit are indicated respectively and are worked as Preceding active injection power；
Tunability is bigger to indicate that the ability of node mitigation Line Flow load is bigger, thus selects tunability size Key node of the node of sequence preceding ten as circuit k, is denoted as set A2；
3. taking set A1With set A2Union as key node set i.e. control variables collection, A=A1∪A2；
Step 3, structure thermal stability security domain：
According to the state variable set F set in the control variables collection A and step 1 determined in step 2, to build Dimensionality reduction thermal stability security domain based on key node injecting power space：
In formula (3), hyperplane factor alphak-iCalculation formula be：
In formula (4), Uk.maxAnd Ik.maxThe respectively reference voltage and rated current of circuit k；
Influence coefficient ckCalculation formula be：
In formula (5), Pk.nowAnd Pk.maxThe current active power and active power limit of circuit k, P are indicated respectivelyi.nowTo work as The active power injection rate of preceding method of operation lower node i.

Claims (1)

1. a kind of thermal stability security domain fast construction method based on key node injecting power space, it is characterized in that：Including with Lower step：
Step 1, setting state variable：
Screen in actual electric network there are effective power flow is out-of-limit or the circuit of heavily loaded problem, setting it is all there are effective power flow it is out-of-limit or The circuit of heavily loaded problem is state variable set F=[circuit 1, circuit 2 ..., circuit k ..., circuit m], element in set F Total number is m；
Step 2 determines control variable, including：
1. effective power flow sensitivity of all node active power outputs variations to circuit in set F in actual electric network is calculated, based on straight Flow tide model, effective power flow sensitivity S of the node i to circuit kk-iCalculation formula is as follows
In formula (1), a, b are two endpoints of circuit k, XaiFor the mutual impedance between node a and node i, XbiFor node b and node Mutual impedance between i, xkFor the reactance of circuit k；
The absolute value of sensitivity is bigger to indicate that node active power output variation is bigger on the influence of the effective power flow of circuit, thus selects Key node of the node of sensitivity order of magnitude sequence preceding ten as circuit k, is denoted as set A1；
2. calculating tunability of all nodes to circuit in set F, tunability P of the node i to circuit kadjCalculation formula is such as Under：
In formula (2), Pi.min、Pi.max、Pi.nowThe active power lower limit of node i, the active power upper limit are indicated respectively and are currently had Work(injecting power；
Tunability is bigger to indicate that the ability of node mitigation Line Flow load is bigger, thus selects the sequence of tunability size Key node of preceding ten node as circuit k, is denoted as set A2；
3. taking set A1With set A2Union as key node set i.e. control variables collection, A=A1∪A2；
Step 3, structure thermal stability security domain：
According to the state variable set F set in the control variables collection A and step 1 determined in step 2, to which structure is based on The dimensionality reduction thermal stability security domain in key node injecting power space：
In formula (3), nAFor the dimension in key node injecting power space, αk-iIt is i-th of key node in k-th of security domain boundaries Hyperplane coefficient, PiFor the injecting power of i-th of key node, ckIndicate k-th of safe domain equation of remaining non-key node pair Influence coefficient；
Wherein, hyperplane factor alphak-iCalculation formula be：
In formula (4), Uk.maxAnd Ik.maxThe respectively reference voltage and rated current of circuit k；
Influence coefficient ckCalculation formula be：
In formula (5), Pk.nowAnd Pk.maxThe current active power and active power limit of circuit k, P are indicated respectivelyi.nowCurrently to transport The active power injection rate of line mode lower node i.
CN201810133357.5A 2018-01-19 2018-02-08 Method for quickly constructing thermal stability security domain based on key node injection power space Active CN108493921B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201810057632 2018-01-19
CN201810057632X 2018-01-19

Publications (2)

Publication Number Publication Date
CN108493921A true CN108493921A (en) 2018-09-04
CN108493921B CN108493921B (en) 2021-04-23

Family

ID=63340070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810133357.5A Active CN108493921B (en) 2018-01-19 2018-02-08 Method for quickly constructing thermal stability security domain based on key node injection power space

Country Status (1)

Country Link
CN (1) CN108493921B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109193813A (en) * 2018-09-30 2019-01-11 中国南方电网有限责任公司 A kind of electric system thermal stability security domain radar drawing drawing method
CN109217339A (en) * 2018-10-25 2019-01-15 国网天津市电力公司 A kind of construction method in the Static Voltage Security domain based on PMU configuration
CN109285089A (en) * 2018-10-25 2019-01-29 天津大学 A kind of screening technique of electric system thermostabilization safety-critical unit
CN109542972A (en) * 2018-09-29 2019-03-29 天津大学 A kind of hyperspace thermal stability security domain representation method based on correlation model

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101281637A (en) * 2008-05-09 2008-10-08 天津大学 Electric power system optimizing swim and real time pricing method based on hyperplane form safety field periphery
CN103413043A (en) * 2013-08-09 2013-11-27 中国南方电网有限责任公司 Solving method for power system multi-dimensional space cross section thermal stability limit boundary
CN104008279A (en) * 2014-05-13 2014-08-27 南京邮电大学 Method for solving power network static security domain
US20150355655A1 (en) * 2014-06-06 2015-12-10 Shanghai Jiao Tong University Method for optimizing the flexible constraints of an electric power system
CN106655201A (en) * 2016-09-20 2017-05-10 天津大学 Security domain-based safe optimization and control method for electric power thermal stability

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101281637A (en) * 2008-05-09 2008-10-08 天津大学 Electric power system optimizing swim and real time pricing method based on hyperplane form safety field periphery
CN103413043A (en) * 2013-08-09 2013-11-27 中国南方电网有限责任公司 Solving method for power system multi-dimensional space cross section thermal stability limit boundary
CN104008279A (en) * 2014-05-13 2014-08-27 南京邮电大学 Method for solving power network static security domain
US20150355655A1 (en) * 2014-06-06 2015-12-10 Shanghai Jiao Tong University Method for optimizing the flexible constraints of an electric power system
CN106655201A (en) * 2016-09-20 2017-05-10 天津大学 Security domain-based safe optimization and control method for electric power thermal stability

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109542972A (en) * 2018-09-29 2019-03-29 天津大学 A kind of hyperspace thermal stability security domain representation method based on correlation model
CN109542972B (en) * 2018-09-29 2023-01-03 天津大学 Multi-dimensional space thermal stability security domain representation method based on correlation model
CN109193813A (en) * 2018-09-30 2019-01-11 中国南方电网有限责任公司 A kind of electric system thermal stability security domain radar drawing drawing method
CN109193813B (en) * 2018-09-30 2021-11-02 中国南方电网有限责任公司 Radar map drawing method for thermal stability security domain of electric power system
CN109217339A (en) * 2018-10-25 2019-01-15 国网天津市电力公司 A kind of construction method in the Static Voltage Security domain based on PMU configuration
CN109285089A (en) * 2018-10-25 2019-01-29 天津大学 A kind of screening technique of electric system thermostabilization safety-critical unit
CN109217339B (en) * 2018-10-25 2022-03-11 国网天津市电力公司 Construction method of static voltage security domain based on PMU configuration

Also Published As

Publication number Publication date
CN108493921B (en) 2021-04-23

Similar Documents

Publication Publication Date Title
CN108493921A (en) Thermal stability security domain fast construction method based on key node injecting power space
Nazir et al. A two-stage chance constrained volt/var control scheme for active distribution networks with nodal power uncertainties
CN100590947C (en) Large electric network on-line preventing control method based on static state and transient safety steady mode
Lian et al. Distributed resilient optimal current sharing control for an islanded DC microgrid under DoS attacks
JP6839769B2 (en) Medium pressure power distribution circuit How to evaluate the safety of network closed loop operation
WO2016197484A1 (en) Optimal configuration method for voltage sag monitoring node
CN104934971B (en) Dynamic section control method based on power flow transfer ratio
EP3375062A1 (en) Hierarchical robust model predictive voltage and var control with coordination and optimization of autonomous der voltage control
Bao et al. Field verification of frequency control by energy-intensive loads for isolated power systems with high penetration of wind power
CN102664417B (en) Control method and device for secondary voltage
CN104281057A (en) Composite PID fuzzy control method applied to transformer cooling system
CN104767199B (en) Directional control method for cut set cross section load flow in electric power system
Selim et al. Fast quasi‐static time‐series analysis and reactive power control of unbalanced distribution systems
CN103928939B (en) Direct-current commutation failure impacts the computational methods of lower alternating current interconnection tie power fluctuation peak
CN107451743B (en) Active power distribution network island operation method considering source load operation characteristics
Ming et al. Mixed H 2/H∞ control for nonlinear closed-loop Stackelberg games with application to power systems
CN103715687B (en) The congested real-time control method of a kind of active distribution network branch power
CN108711850A (en) The determination method of safe Loop Closing Operation is carried out to medium voltage distribution network
Zhang et al. Real-time optimal voltage control using measurement-based aggregate load model
CN106159957B (en) A kind of optimization method exchanging weak section steady-state overvoltage Pre-control measures
WO2017107358A1 (en) Method for controlling power flow of chained stable tie-lines of power grid, and server
CN110071499B (en) Method and system for rapidly judging safety loop closing of power distribution network
CN108233384B (en) Method for stabilizing load side voltage under network attack based on Petri network
KR101848993B1 (en) Dispatcher Training Simulator
CN103248052B (en) Saturated switching control method for three-phase parallel active power filter

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