CN106253313A - The determination method of power system node injecting power sensitivity coefficient - Google Patents
The determination method of power system node injecting power sensitivity coefficient Download PDFInfo
- Publication number
- CN106253313A CN106253313A CN201610733565.XA CN201610733565A CN106253313A CN 106253313 A CN106253313 A CN 106253313A CN 201610733565 A CN201610733565 A CN 201610733565A CN 106253313 A CN106253313 A CN 106253313A
- Authority
- CN
- China
- Prior art keywords
- power
- node
- generator
- formula
- voltage
- 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
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
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- 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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
The invention discloses a kind of determination method of power system node injecting power sensitivity coefficient, when ac and dc systems fault and emergency, the method coordinating between a plurality of direct current or direct current from different energy storage and Load Emergency to control.The generator electromagnetic power to be studied caused due to certain node injection unit power is changed to the sensitivity to this generator electromagnetic power of the node injecting power, which just can come it is determined that at node injecting power according to sensitivity, the size of injecting power, injecting power moment and persistent period, coordination between multiple node injecting powers, etc..The calculating process of node injecting power sensitivity coefficient is as follows: according to given electrical network, generates the whole network bus admittance matrix;Calculate initial trend, obtain each generator terminal voltage current amplitude phase angle initial value;Based on set end voltage electric current, calculate electromotor transient potential and generator's power and angle initial value;Meter sensitivity coefficient matrix, the transient stability of system when improve fault and emergency.
Description
Technical field
The invention belongs to electrical engineering technical field, relate to the merit under AC and DC power system operation troubles and the state of emergency
Rate is coordinated to control, and is specifically related to a kind of determination method of power system node injecting power sensitivity coefficient.
Background technology
Ultra-high/extra-high voltage HVDC Transmission Technology obtains wideling popularize in terms of China's transferring electricity from the west to the east and regional power grid interconnection answers
With, define some typical electrical web frames with distinct characteristics, if southwest Hydropower Base electric power is by three big extra-high voltages
Many direct currents delivery system that direct current is sent outside, northwester light thermoelectricity is sent outside through the alternating current-direct current mixing that extra-high voltage alternating current-direct current passage is sent outside and is
System, the Ac/dc Power Systems that Central China's multi-infeed DC/many direct currents are sent outside, East China many direct current channels feed-in receiving end electricity
Force system, the asynchronous interconnected network in south in planning, etc., additionally occur in that high-capacity direct circulation road receiving end is layered multiple spot feed-in
Form.Under electrical network normal operating condition, direct current transportation is responsible for the task of trans-regional large-capacity power conveying.Fault is with urgent
Under state, the quick controlled ability of direct current transportation can be that network stability control provides preferable means, can pass through multiple-circuit line
Coordination between passage controls, and improves the stability of electrical network.Utilize dc power emergent control can improve power system transient stability
Property, it is possible to provide control action for system is retracted synchronous regime from desynchronizing state.But choose which bar or which bar direct current
Passage is implemented to control, and controlled quentity controlled variable is how many, controls how the moment determines how to optimize if using many direct currents to coordinate this, etc., need
Study suitable theoretical method as control design case foundation.In addition to utilizing direct-current emergency power and controlling, energy storage and promptly bearing
Lotus controls to be all that available electric power system transient stability controls technological means, relates to the installation site of energy storage, choosing among these
Take which node load to be controlled, and control to measure the most suitable and each node injecting power control between how coordination etc.
Problem.
Summary of the invention
The invention provides a kind of determination method of power system node injecting power sensitivity coefficient, solve a plurality of directly
The technical problem coordinating to control between stream or direct current and multiple energy storage and Load Emergency.
The present invention is to solve above technical problem by the following technical programs:
The determination method of power system node injecting power sensitivity coefficient, utilizes dc power emergent control to improve system transient modelling
Stability, retracts synchronous regime for system from desynchronizing state and provides control action.Energy storage and emergency load control are available
Electric power system transient stability controls technological means.These technological means all show as the change of certain node injecting power, from net
Network equation can solve this injecting power effect to the impact of generator electromagnetic power in electrical network, i.e. so-called node injecting power spirit
The concept of sensitivity.Therefore, define owing to certain node injects the generator electromagnetic power change that unit power causes, inject for node
The power sensitivity to this generator electromagnetic power, it is possible to come it is determined that at which node injecting power according to sensitivity,
The size of injecting power, injecting power moment and persistent period, the coordination between multiple node injecting powers.
Specifically comprise the following steps that
Relational expression between the first step, grid nodes injecting power and each node voltage amplitude phase angle is:
;
Second step, node power equation is carried out linearisation, can obtain:
;
In formula: Δ P=[ΔP 1 ΔP 2 … ΔP n ]T, Δ Q=[ΔQ 1 ΔQ 2 … ΔQ n ]T, respectively node injects active power
With reactive power deviation;Δθ=[Δθ 1 Δθ 2 … Δθ n ]T, Δ V=[ΔV 1 ΔV 2 … ΔV n ]T, respectively node voltage phase
Angle and amplitude deviation;VD=diag[V 1 V 2 … V n ] be node voltage amplitude composition diagonal matrix;
3rd step, wheni≠jTime:
;
Wheni=jTime:
;
4th step, the formula inversion operation that second step is obtained:
;
In formula: S11, S21Represent that node injects active power to the whole network node voltage phase angle and the impact of amplitude respectively;S12, S22
Represent that node injects reactive power to the whole network node voltage phase angle and the impact of amplitude respectively;
5th step, generator electromagnetic power linearisation:
Assuming that the whole network has m platform synchronous generator, and the most numbered 1-m, all use classical model (transient potentialE′Permanent
Fixed), Qi ZhongiThe electromagnetic power of platform electromotor is represented by:
;
In formula:E Gi 'It is the transient potential of the electromotor,V Gi For generator terminal voltage amplitude, machine end node correspondence system saves
Pointi; i 'For the phase angle difference between electromotor built-in potential and set end voltage, equal to generator's power and angle and set end voltage phase angle it
Difference;Generator's power and angleWill not dash forward relative to the angular displacement of the X-axis of reference synchronization coordinate system XY, disturbance moment for generator amature
Become;Set end voltage phase angle i For set end voltage phasor relative to the phase place of the X-axis of reference synchronization coordinate system XY, transient stability analysis
In, network node voltage magnitude and phase place including Generator end are all undergone mutation;
It should be noted that the balance node chosen in Load flow calculation, it is intended that its voltage-phase is zero, represents initial steady state operating mode
Under, the voltage phasor just X-direction with XY reference coordinate of balance node overlaps, in the case of disturbance, and the voltage amplitude of this node
Value phase place also can change, and no longer overlaps with X-axis;
Injecting power disturbance issues motor electromagnetic changed power and can be expressed as:
;
In formula:
;
Above formula is written as following form:
;
In formula: M i1And M i2It is n unit row vector, wherein except theiIndividual element other elements outer are zero, have:
;
If system m platform electromotor row being write and stating equation, have:
;
In formula: Δ Pe=[ΔP e1 ΔP e2 … ΔP em ]T, change column vector for each generator electromagnetic power;
;
;
6th step, determine node injecting power sensitivity coefficient:
Can be obtained by above formula:
;
Represent the node injecting power change impact for each generator electromagnetic power;
Definition sensitivity matrix:
;
In formula: s-matrix the element representation impact on generator electromagnetic power of the node injecting power, referred to as electromotor electromagnetism merit
The rate sensitivity to node injecting power.
The present invention is based on node injecting power sensitivity coefficient, it is provided that a kind of in ac and dc systems fault and emergency
Time, the method coordinating between a plurality of direct current or direct current from different energy storage and Load Emergency to control.Use the method can improve event
The transient stability of system when barrier and emergency.
Detailed description of the invention
The determination method of power system node injecting power sensitivity coefficient, comprises the following steps:
Relational expression between the first step, grid nodes injecting power and each node voltage amplitude phase angle is:
;
Second step, node power equation is carried out linearisation, can obtain:
;
In formula: Δ P=[ΔP 1 ΔP 2 … ΔP n ]T, Δ Q=[ΔQ 1 ΔQ 2 … ΔQ n ]T, respectively node injects active power
With reactive power deviation;Δθ=[Δθ 1 Δθ 2 … Δθ n ]T, Δ V=[ΔV 1 ΔV 2 … ΔV n ]T, respectively node voltage phase
Angle and amplitude deviation;VD=diag[V 1 V 2 … V n ] be node voltage amplitude composition diagonal matrix;
3rd step, wheni≠jTime:
;
Wheni=jTime:
;
4th step, the formula inversion operation that second step is obtained:
;
In formula: S11, S21Represent that node injects active power to the whole network node voltage phase angle and the impact of amplitude respectively;S12, S22
Represent that node injects reactive power to the whole network node voltage phase angle and the impact of amplitude respectively;
5th step, generator electromagnetic power linearisation:
Assuming that the whole network has m platform synchronous generator, and the most numbered 1-m, all use classical model (transient potentialE′Permanent
Fixed), Qi ZhongiThe electromagnetic power of platform electromotor is represented by:
;
In formula:E Gi 'It is the transient potential of the electromotor,V Gi For generator terminal voltage amplitude, machine end node correspondence system saves
Pointi; i 'For the phase angle difference between electromotor built-in potential and set end voltage, equal to generator's power and angle and set end voltage phase angle it
Difference;Generator's power and angleWill not dash forward relative to the angular displacement of the X-axis of reference synchronization coordinate system XY, disturbance moment for generator amature
Become;Set end voltage phase angle i For set end voltage phasor relative to the phase place of the X-axis of reference synchronization coordinate system XY, transient stability analysis
In, network node voltage magnitude and phase place including Generator end are all undergone mutation;
It should be noted that the balance node chosen in Load flow calculation, it is intended that its voltage-phase is zero, represents initial steady state operating mode
Under, the voltage phasor just X-direction with XY reference coordinate of balance node overlaps, in the case of disturbance, and the voltage amplitude of this node
Value phase place also can change, and no longer overlaps with X-axis;
Injecting power disturbance issues motor electromagnetic changed power and can be expressed as:
;
In formula:
;
Above formula is written as following form:
;
In formula: M i1And M i2It is n unit row vector, wherein except theiIndividual element other elements outer are zero, have:
;
If system m platform electromotor row being write and stating equation, have:
;
In formula: Δ Pe=[ΔP e1 ΔP e2 … ΔP em ]T, change column vector for each generator electromagnetic power;
;
;
6th step, determine node injecting power sensitivity coefficient:
Can be obtained by above formula:
;
Represent the node injecting power change impact for each generator electromagnetic power;
Definition sensitivity matrix:
;
In formula: s-matrix the element representation impact on generator electromagnetic power of the node injecting power, referred to as electromotor electromagnetism merit
The rate sensitivity to node injecting power.
Claims (1)
1. a determination method for power system node injecting power sensitivity coefficient, comprises the following steps:
Relational expression between the first step, grid nodes injecting power and each node voltage amplitude phase angle is:
;
Second step, node power equation is carried out linearisation, can obtain:
;
In formula: Δ P=[ΔP 1 ΔP 2 … ΔP n ]T, Δ Q=[ΔQ 1 ΔQ 2 … ΔQ n ]T, respectively node injects active power
With reactive power deviation;Δθ=[Δθ 1 Δθ 2 … Δθ n ]T, Δ V=[ΔV 1 ΔV 2 … ΔV n ]T, respectively node voltage phase
Angle and amplitude deviation;VD=diag[V 1 V 2 … V n ] be node voltage amplitude composition diagonal matrix;
3rd step, wheni≠jTime:
;
Wheni=jTime:
;
4th step, the formula inversion operation that second step is obtained:
;
In formula: S11, S21Represent that node injects active power to the whole network node voltage phase angle and the impact of amplitude respectively;S12, S22Point
Biao Shi not inject reactive power to the whole network node voltage phase angle and the impact of amplitude by node;
5th step, generator electromagnetic power linearisation:
Assuming that the whole network has m platform synchronous generator, and the most numbered 1-m, all use classical model (transient potentialE′Permanent
Fixed), Qi ZhongiThe electromagnetic power of platform electromotor is represented by:
;
In formula:E Gi 'It is the transient potential of the electromotor,V Gi For generator terminal voltage amplitude, machine end node correspondence system saves
Pointi; i 'For the phase angle difference between electromotor built-in potential and set end voltage, equal to generator's power and angle and set end voltage phase angle it
Difference;Generator's power and angleWill not dash forward relative to the angular displacement of the X-axis of reference synchronization coordinate system XY, disturbance moment for generator amature
Become;Set end voltage phase angle i For set end voltage phasor relative to the phase place of the X-axis of reference synchronization coordinate system XY, transient stability analysis
In, network node voltage magnitude and phase place including Generator end are all undergone mutation;
It should be noted that the balance node chosen in Load flow calculation, it is intended that its voltage-phase is zero, represents initial steady state operating mode
Under, the voltage phasor just X-direction with XY reference coordinate of balance node overlaps, in the case of disturbance, and the voltage amplitude of this node
Value phase place also can change, and no longer overlaps with X-axis;
Injecting power disturbance issues motor electromagnetic changed power and can be expressed as:
;
In formula:
;
Above formula is written as following form:
;
In formula: M i1And M i2It is n unit row vector, wherein except theiIndividual element other elements outer are zero, have:
;
If system m platform electromotor row being write and stating equation, have:
;
In formula: Δ Pe=[ΔP e1 ΔP e2 … ΔP em ]T, change column vector for each generator electromagnetic power;
;
;
6th step, determine node injecting power sensitivity coefficient:
Can be obtained by above formula:
;
Represent the node injecting power change impact for each generator electromagnetic power;
Definition sensitivity matrix:
;
In formula: s-matrix the element representation impact on generator electromagnetic power of the node injecting power, referred to as electromotor electromagnetism merit
The rate sensitivity to node injecting power.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610733565.XA CN106253313B (en) | 2016-08-27 | 2016-08-27 | The determination method of electric system node injecting power sensitivity coefficient |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610733565.XA CN106253313B (en) | 2016-08-27 | 2016-08-27 | The determination method of electric system node injecting power sensitivity coefficient |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106253313A true CN106253313A (en) | 2016-12-21 |
CN106253313B CN106253313B (en) | 2019-04-16 |
Family
ID=57595564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610733565.XA Active CN106253313B (en) | 2016-08-27 | 2016-08-27 | The determination method of electric system node injecting power sensitivity coefficient |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106253313B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109256790A (en) * | 2018-10-22 | 2019-01-22 | 暨南大学 | A kind of energy-storage system configuration method, device and storage medium |
CN110474352A (en) * | 2019-07-05 | 2019-11-19 | 国家电网有限公司 | Operating condition evaluation method, device and its electronic equipment of generator motor |
CN110535136A (en) * | 2019-08-16 | 2019-12-03 | 广东工业大学 | A kind of AC power flow linear sensitivity analysis method |
CN110488193B (en) * | 2019-09-19 | 2020-10-02 | 广东电网有限责任公司 | Synchronous generator transient power angle solving method |
CN111864741A (en) * | 2020-07-27 | 2020-10-30 | 国网安徽省电力有限公司电力科学研究院 | Quantitative analysis method and system for influence of line parameter errors on power distribution |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102403724A (en) * | 2011-11-09 | 2012-04-04 | 深圳大学 | Symmetrical obtaining method for nodal voltage sensitivity in alternating-direct current parallel-serial power grid |
-
2016
- 2016-08-27 CN CN201610733565.XA patent/CN106253313B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102403724A (en) * | 2011-11-09 | 2012-04-04 | 深圳大学 | Symmetrical obtaining method for nodal voltage sensitivity in alternating-direct current parallel-serial power grid |
Non-Patent Citations (1)
Title |
---|
董飞飞等: "基于节点注入功率—电压变化比灵敏度的交换容量极限计算研究", 《陕西电力》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109256790A (en) * | 2018-10-22 | 2019-01-22 | 暨南大学 | A kind of energy-storage system configuration method, device and storage medium |
CN109256790B (en) * | 2018-10-22 | 2022-02-11 | 暨南大学 | Energy storage system configuration method and device and storage medium |
CN110474352A (en) * | 2019-07-05 | 2019-11-19 | 国家电网有限公司 | Operating condition evaluation method, device and its electronic equipment of generator motor |
CN110535136A (en) * | 2019-08-16 | 2019-12-03 | 广东工业大学 | A kind of AC power flow linear sensitivity analysis method |
CN110488193B (en) * | 2019-09-19 | 2020-10-02 | 广东电网有限责任公司 | Synchronous generator transient power angle solving method |
CN111864741A (en) * | 2020-07-27 | 2020-10-30 | 国网安徽省电力有限公司电力科学研究院 | Quantitative analysis method and system for influence of line parameter errors on power distribution |
Also Published As
Publication number | Publication date |
---|---|
CN106253313B (en) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106253313A (en) | The determination method of power system node injecting power sensitivity coefficient | |
CN106786724B (en) | A kind of control strategy of more times MMC-HVDC feed-ins pole weak grid | |
Wang et al. | Fault analysis for distribution networks with current-controlled three-phase inverter-interfaced distributed generators | |
Dobakhshari et al. | Control of microgrids: Aspects and prospects | |
CN104716886B (en) | based on virtual impedance voltage type converter virtual synchronous motor control method | |
Kim et al. | A novel three-phase four-leg inverter based load unbalance compensator for stand-alone microgrid | |
CN102570477B (en) | SVC and TCSC based on linearisation power flow algorithm model mixes planing method of layouting | |
CN110429632B (en) | Frequency consistency control method for asynchronous interconnected system containing double-loop flexible direct current two areas | |
Benchaib | Advanced control of AC/DC power networks: system of systems approach based on spatio-temporal scales | |
Kamel et al. | Novel and simple scheme for Micro-Grid protection by connecting its loads neutral points: A review on Micro-Grid protection techniques | |
Wang et al. | Voltage stability analysis of power systems with a large number of non-synchronous machine sources connected | |
Ndreko et al. | Impact of DC voltage control parameters on AC/DC system dynamics under faulted conditions | |
Hamoud et al. | Voltage sag and swell mitigation using D-STATCOM in renewable energy based distributed generation systems | |
Huang et al. | Transient stability analysis of Shanghai Power Grid with multiple HVDC links | |
Moreira et al. | Participation of multi-terminal HVDC grids in frequency regulation services | |
Li et al. | Research of nonlinear control strategy for VSC-HVDC system based on Lyapunov stability theory | |
Maleki et al. | Comparative study for improving damping oscillation of SMIB system with STATCOM and BESS using remote and local signal | |
Azbe et al. | Energy function for an interline power-flow controller | |
Feng et al. | Influence on transmission line relay protection under DFIG-based wind farm intergration | |
Jinhua et al. | Research on coordinated control technology for AC/DC hybrid system in China | |
Yi-xi et al. | An improved droop control strategy in the islanding operation | |
Mobarak et al. | Voltage Instability of Initiation Fault Duration as Influenced by Nodes Short Circuit Levels NSCL. | |
Wang et al. | The power flow calculation based on impedance specifications for low-voltage AC microgrid | |
Suvorov et al. | Analysis of the effect of wind power plants on the out-of-step protection operation | |
Lin et al. | Control Strategy for Low Frequency AC Transmission Under Single-Phase Outage Condition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |