CN106972490B - The method and device of power transmission network Contingency analysis - Google Patents
The method and device of power transmission network Contingency analysis Download PDFInfo
- Publication number
- CN106972490B CN106972490B CN201710403029.8A CN201710403029A CN106972490B CN 106972490 B CN106972490 B CN 106972490B CN 201710403029 A CN201710403029 A CN 201710403029A CN 106972490 B CN106972490 B CN 106972490B
- Authority
- CN
- China
- Prior art keywords
- power
- power transmission
- transmission network
- distribution network
- node
- 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.)
- Active
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 181
- 238000004458 analytical method Methods 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 59
- 238000005259 measurement Methods 0.000 claims description 28
- 230000005611 electricity Effects 0.000 claims description 23
- 230000001351 cycling effect Effects 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 238000002347 injection Methods 0.000 claims 1
- 239000007924 injection Substances 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005520 electrodynamics Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
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/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
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of method and devices of power transmission network Contingency analysis.This method comprises: receiving the distribution network Equivalent Model that power distribution network control centre sends；The distribution network Equivalent Model is added in the Contingency analysis model of power transmission network, carries out Load flow calculation, wherein the Contingency analysis model includes multiple forecast accidents；According to the calculation of tidal current under each forecast accident, safety in operation of the power transmission network under the forecast accident is determined.The influence of distribution network change in topology in power transmission network accident is carried out approximation with distribution network equivalence by the present invention, is realized and is obtained more accurate power transmission network Contingency analysis result under few calculating and communication cost.
Description
Technical field
The present embodiments relate to Operation of Electric Systems and control technology more particularly to a kind of power transmission network Contingency analysis
Method and device.
Background technique
Influence of the existing power transmission network Contingency analysis method to power distribution network has only carried out very "ball-park" estimate: dividing
The influence of power distribution network injecting power at power distribution station low-voltage bus bar is only accounted in analysis, but ignores distribution network change in topology institute
The influence of caused power transfer.With the fast development of Information and Communication Technology, power distribution network, which has been assembled with, more and more to exist
Line automatic control device, this makes the topological structure of distribution network that may change in real time execution, continues to use existing defeated
Power grid Contingency analysis method may cause false alarm.One main cause of University of Southern California, U.S. power outage in 2011 is just
It is not accounting for the variation that distribution network may occur in operation in existing power transmission network Contingency analysis program, to make
At false alarm, the serious expansion of accident has been eventually led to.This shows to consider that distribution network is opened up in power transmission network Contingency analysis
The influence for flutterring variation has become extremely important.
In the prior art, global tidal current computing method can accurately embody distribution network change in topology to power transmission network trend
Influence, but calculate primary global trend the control centre of power transmission network and power distribution network needed to calculate separately and transmit electricity and match several times
Electric trend, and communicated several times.For practical power systems, it is contemplated that the number of accident is usually very huge, and pre- every time
Think to require to calculate primary global trend in accident, so for completion Contingency analysis, in the control of power transmission network and power distribution network
The heart must calculate separately over thousands of transmission and disttrbution trend and carry out over thousands of communication.So huge calculating and communication generation
Valence makes this method be difficult application on site in the Contingency analysis program of practical power systems.
Summary of the invention
In view of this, the embodiment of the present invention provides a kind of method and device of power transmission network Contingency analysis, few
Communication and calculate and obtain more accurate anticipation under cost and analyze result.
In a first aspect, the embodiment of the invention provides a kind of methods of power transmission network Contingency analysis, which comprises
Receive the distribution network Equivalent Model that power distribution network control centre sends；
The distribution network Equivalent Model is added in the Contingency analysis model of power transmission network, Load flow calculation is carried out,
Wherein, the Contingency analysis model includes multiple forecast accidents；
According to the calculation of tidal current under each forecast accident, operation peace of the power transmission network under the forecast accident is determined
Quan Xing.
Second aspect, the embodiment of the invention also provides a kind of device of power transmission network Contingency analysis, described device packets
It includes:
Equivalent Network model receiving module, the distribution network Equivalent Model sent for receiving power distribution network control centre；
Load flow calculation module, for the distribution network Equivalent Model to be added to the Contingency analysis model of power transmission network
In, carry out Load flow calculation, wherein the Contingency analysis model includes multiple forecast accidents；
Safety in operation determination module, for determining the transmission of electricity according to the calculation of tidal current under each forecast accident
Safety in operation of the net under the forecast accident.
The technical solution of the embodiment of the present invention, the distribution network equivalence mould sent by receiving power distribution network control centre
Distribution network Equivalent Model is added in the Contingency analysis model of power transmission network by type, Load flow calculation is carried out, according to each pre-
Think the calculation of tidal current under accident, determine safety in operation of the power transmission network under the forecast accident, due to distribution network equivalence
Model can approximation reflect distribution network change in topology, by the influence of the distribution network change in topology in power transmission network accident with matching
Electric network equivalence carries out approximation, realizes and obtains more accurate power transmission network forecast accident under few calculating and communication cost
Analyze result.
Detailed description of the invention
Fig. 1 is a kind of flow chart of the method for power transmission network Contingency analysis that the embodiment of the present invention one provides；
Fig. 2 is the distribution network Equivalent Model in the method for power transmission network Contingency analysis provided in an embodiment of the present invention
Schematic diagram；
Fig. 3 is a kind of flow chart of the method for power transmission network Contingency analysis provided by Embodiment 2 of the present invention；
Fig. 4 is a kind of flow chart of the method for power transmission network Contingency analysis that the embodiment of the present invention three provides；
Fig. 5 is a kind of structural schematic diagram of the device for power transmission network Contingency analysis that the embodiment of the present invention four provides.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining the present invention rather than limiting the invention.It also should be noted that in order to just
In description, only some but not all contents related to the present invention are shown in the drawings.
Embodiment one
Fig. 1 is a kind of flow chart of the method for power transmission network Contingency analysis that the embodiment of the present invention one provides, this implementation
Example is applicable to carry out power transmission network Contingency analysis according to the distribution network Equivalent Model for considering distribution network change in topology
The case where, this method can be executed by the device of power transmission network Contingency analysis, the device can by software and/or hardware Lai
It realizes, is generally configurable in power transmission network control centre, this method specifically comprises the following steps:
Step 110, the distribution network Equivalent Model that power distribution network control centre sends is received.
When distribution network topological structure changes, power distribution network control centre can according to Gaussian reduction or other
Method calculate distribution network Equivalent Model, by power distribution network node and load equivalence into boundary node, equivalence led for parallel connection
It receives and series admittance, and the distribution network Equivalent Model being calculated is sent to power transmission network control centre.Distribution network is equivalent
Model can approximation reflect distribution network change in topology.
Wherein, the distribution network Equivalent Model preferably includes at least one boundary node, and the boundary node is distribution
The node that net is connected with power transmission network, each boundary node passes through equivalent shunt admittance ground connection, if boundary node is no less than two,
Then connected between every two boundary node by equivalent series admittance.Fig. 2 is power transmission network anticipation thing provided in an embodiment of the present invention
Therefore the schematic diagram of the distribution network Equivalent Model in the method for analysis.As shown in Fig. 2, distribution network Equivalent Model includes boundary section
Point B1, boundary node B2..., boundary node Bk(k >=3), boundary node B1Pass through equivalent shunt admittance yB1Ground connection, boundary node
B2Pass through equivalent shunt admittance yB2Ground connection, boundary node BkPass through equivalent shunt admittance yBkGround connection, boundary node B1It is saved with boundary
Point B2Between be in series with equivalent series admittance yB12, boundary node B1With boundary node BkBetween be in series with equivalent series admittance yB1k,
Boundary node B2With boundary node BkBetween be in series with equivalent series admittance yB2k.Those skilled in the art will be seen that Fig. 2 is
One example, what is provided is to have the case where at least three boundary nodes, can also include one except in the case of in Fig. 2
Boundary node or two boundary nodes only need to be by the nodes and load equivalence in power distribution network in only one boundary node
For the shunt admittance of the boundary node, i.e., in distribution network Equivalent Model, which passes through equivalent shunt admittance ground connection；
When there are two boundary node, Equivalent Network model is similar with Fig. 2, and only a boundary node, i.e. two boundary nodes divide less
It is not grounded by respective equivalent shunt admittance, is connected between the two boundary nodes by equivalent series admittance.Wherein, equivalent
Shunt admittance is that the load of boundary node in power distribution network is carried out equivalence to obtain, and equivalent series admittance is by two in power distribution network
Load between boundary node carries out what equivalence obtained.
Step 120, the distribution network Equivalent Model is added in the Contingency analysis model of power transmission network, carries out tide
Stream calculation, wherein the Contingency analysis model includes multiple forecast accidents.
Wherein, Load flow calculation is electrodynamic noun, is first according to given electric network composition, parameter and generator, load etc.
The service condition of part determines the calculating of electric system each section steady-state operating condition parameter.Usually given service condition, which has, is
Each power supply and the power of load point, pivot point voltage, the voltage of equalization point and phase angle, running state parameters packet to be asked in system
Include power distribution, the power loss of network etc. of the voltage magnitude and phase angle and each branch of each bus nodes of power grid.Trend meter
It is the very important analytical calculation of electric system, to study the various problems proposed in systems organization and operation.To planning
In electric system, by Load flow calculation can with inspection institute propose Power System Planning scheme can meet the various methods of operation
Requirement；To running electric system, the change meeting of various load variations and network structure anticipated that by Load flow calculation
The safety of system will not be jeopardized, the voltage of all buses is whether within the scope of permission in system, various element (lines in system
Road, transformer etc.) whether will appear overload, and while being likely to occur overload, should take which precautionary measures etc. in advance.
By distribution network topological structure with distribution network Equivalent Model come approximate representation, in conjunction with the network topology knot of power transmission network
Structure carries out Load flow calculation to each forecast accident in the Contingency analysis model of power transmission network respectively, seeks in power transmission network
The transimission power of the voltage of each node and each transmission line of electricity.Forecast accident includes that the generator in power transmission network is out of service
Or transmission line of electricity is out of service, for actual electric system, it is contemplated that Accident analysis model includes at least two forecast accidents.
Step 130, according to the calculation of tidal current under each forecast accident, determine the power transmission network under the forecast accident
Safety in operation.
According to the calculation of tidal current under each forecast accident, the operating parameter and transmission of electricity of the node in power transmission network are judged
Whether the operating parameter of route meets the condition of setting.When meeting the condition of setting, power transmission network operational safety is determined；It is unsatisfactory for setting
When fixed condition, determining power transmission network, there are security risks, and provide corresponding safe early warning.
Wherein, according to the calculation of tidal current under each forecast accident, determine the power transmission network under the forecast accident
Safety in operation, it is optional to include:
According to the calculation of tidal current under each forecast accident, judge whether have in the power transmission network under the forecast accident
Whether node voltage is more than the default working voltage upper bound or default working voltage lower bound, and have the transimission power of transmission line of electricity super
Cross default operation upper power bound or default operation lower power bound；
If there is node voltage more than the default working voltage upper bound or default working voltage lower bound in the power transmission network, and/
Or, the transimission power for having transmission line of electricity is more than default operation upper power bound or default operation lower power bound, then power transmission network fortune is alerted
There is hidden danger of crossing the border in row.
Wherein, the operation upper bound that the working voltage upper bound is node voltage given in advance is preset, working voltage lower bound is preset
It is the operation lower bound of node voltage given in advance, default operation upper power bound is the transimission power of transmission line of electricity given in advance
The operation upper bound, default operation lower power bound is the operation lower bound of the transimission power of transmission line of electricity given in advance.
If there are at least one of following situations in power transmission network: having node in power transmission network in a forecast accident
Voltage is more than that default working voltage is previous or default working voltage lower bound, and the transimission power for having transmission line of electricity is more than default fortune
Row upper power bound or default operation lower power bound, then alert power transmission network operation and there is hidden danger of crossing the border, otherwise, it determines in the anticipation thing
Under therefore, power transmission network calculation of tidal current meets operating range given in advance.
The technical solution of the present embodiment, the distribution network Equivalent Model sent by receiving power distribution network control centre will
Distribution network Equivalent Model is added in the Contingency analysis model of power transmission network, Load flow calculation is carried out, according to each anticipation thing
Calculation of tidal current under therefore determines safety in operation of the power transmission network under the forecast accident, due to distribution network Equivalent Model
Can approximation reflect distribution network change in topology, by the influence power distribution network of the distribution network change in topology in power transmission network accident
Network equivalence carries out approximation, realizes and obtains more accurate power transmission network Contingency analysis under few calculating and communication cost
As a result.
Embodiment two
Fig. 3 is a kind of flow chart of the method for power transmission network Contingency analysis provided by Embodiment 2 of the present invention, this implementation
Example is optimized on the basis of the above embodiments, to " the distribution network Equivalent Model being added to the anticipation of power transmission network
In Accident analysis model, Load flow calculation is carried out " it is optimized, this method specifically comprises the following steps:
Step 310, the distribution network Equivalent Model that power distribution network control centre sends is received.
Wherein, the distribution network Equivalent Model includes at least one boundary node, the boundary node be power distribution network and
The node that power transmission network is connected, each boundary node passes through equivalent shunt admittance ground connection, if boundary node is no less than two, often
It is connected between two boundary nodes by equivalent series admittance.
Step 320, the current measurement voltage and current measured power of each boundary node are obtained.
Sensor can be set in boundary node, for acquiring the voltage and power of the current time boundary node, collects
Voltage be current measurement voltage, collected power is current measured power.Power distribution network control centre or power transmission network control
Center can directly be obtained from sensor in the current measurement voltage and current measured power for needing boundary node.
Step 330, it is led according to the current measurement voltage, current measured power, equivalent shunt admittance and equivalent series
It receives, calculates the net injecting power vector of the boundary node in the distribution network Equivalent Model.
According to the current measurement voltage of each boundary node, equivalent shunt admittance and equivalent series admittance, power distribution network is calculated
Equivalent power vector of the network Equivalent Model at boundary node, measurement power and equivalent power further according to each boundary node to
Measure the net injecting power vector for the boundary node being calculated in distribution network Equivalent Model.
Wherein, according to the current measurement voltage, current measured power, equivalent shunt admittance and equivalent series admittance,
The net injecting power vector of the boundary node in the distribution network Equivalent Model is calculated, it is optional to include:
According to the current measurement voltage, equivalent shunt admittance and equivalent series admittance, it is equivalent to calculate the distribution network
Equivalent power vector of the model at the boundary node；
According to the current measured power and equivalent power vector, the boundary section in the distribution network Equivalent Model is calculated
The net injecting power vector of point.
It is led first according to the equivalent parallel connection in the current measurement voltage and distribution network Equivalent Model of each boundary node
It receives with equivalent series admittance, the equivalence for calculating each boundary node introduced by distribution network Equivalent Model according to formula (1) is electric
Flow the equal currents vector of composition:
Wherein, VB1,0、VB2,0And VBk,0Respectively indicate boundary node B1Current measurement voltage, boundary node B2 it is current
Measure voltage and boundary node BkCurrent measurement voltage, ID1,0Indicate currently measurement voltage under at all boundary nodes
The dimension of the equal currents vector introduced by above-mentioned distribution network Equivalent Model, the equal currents vector is equal to of boundary node
Number.
After equal currents vector is calculated, further according to the current measurement electricity of each boundary node of equal currents vector sum
Pressure, is calculated equivalent power vector of the distribution network Equivalent Model at the boundary node according to formula (2):
Wherein, SD1,0Indicate currently measurement voltage under by above-mentioned distribution network Equivalent Model in above-mentioned boundary node
Equivalent power vector, symbol conj indicates each element in vector being converted to its corresponding conjugate complex number.
According to equivalent power vector sum current measured power, it is calculated in distribution network Equivalent Model according to formula (3)
Boundary node net injecting power vector:
SD2,0=SD,0-SD1,0 (3)
Wherein, SD2,0Indicate the net injecting power vector of boundary node, SD,0Indicate the measurement power S by boundary nodeB1,0、
SB2,0And SBk,0The column vector of composition, concrete form SD,0=trp ([SB1,0,SB2,0,…,SBk,0]), symbol trp is indicated to one
A vector carries out transposition operation.SB1,0Indicate boundary node B1Current measured power, SB2,0Indicate boundary node B2Current survey
Measure power, SB3,0Indicate boundary node B3Current measured power.Due to net injecting power vector by voltage change influenced compared with
Small, essentially definite value can be used directly in subsequent power flow equation.
Step 340, a forecast accident to be analyzed in power transmission network, the forecast accident as present analysis are obtained.
The forecast accident that can be analysed in advance finishes sequence, analyzes in sequence each forecast accident, i.e.,
First according to the sequence of forecast accident to be analyzed, one of those forecast accident to be analyzed is obtained, as the pre- of present analysis
Think accident.
Step 350, according to the distribution network Equivalent Model and net injecting power vector, the power transmission network is established current
Power flow equation under the forecast accident of analysis, and solve, obtain calculation of tidal current.
According to distribution network Equivalent Model and net injecting power vector, power transmission network is established under the forecast accident of present analysis
Power flow equation, such as formula (4):
Wherein, VTIndicate column vector composed by the voltage of each node of power transmission network, V under the forecast accident of present analysisB
Indicate the column vector composed by the voltage of the forecast accident lower boundary node of present analysis, STIndicate each node in power transmission network
Injecting power composition column vector, STTIndicate the column vector of the line power composition in power transmission network between every two node, STB
Indicate the column vector that the line power that boundary node is flowed to from power transmission network node forms, SBTIt indicates to flow to power transmission network from boundary node
The column vector of the line power composition of node, SBBIndicate the line power composition in boundary node between every two boundary node
Column vector, SD1It indicates in boundary node voltage VBUnder from equivalent power of the distribution network Equivalent Model at boundary node to
Amount, ID1It indicates in boundary node voltage VBUnder the equal currents introduced at all boundary nodes by distribution network Equivalent Model
Vector.VB1、VB2And VBkIndicate the boundary node B under the forecast accident of present analysis1、B2And BkVoltage, they meet relationship
Formula VB=trp ([VB1,VB2,…,VBk]).Therefore, above-mentioned vector STT、SBT、SBBAnd SD1It can be written as about VT、VBSecondary letter
Number, i.e., common " power equation " form in Power System Analysis textbook.
It solves above-mentioned power flow equation, obtains calculation of tidal current, i.e., the voltage of each node in power transmission network, then will transmission of electricity
The voltage of each node substitutes into the transimission power that the route in power transmission network between every two node is calculated in power equation in net.
Step 360, judge whether to have analyzed all forecast accidents to be analyzed, if it is not, then 340 are returned to step,
If so, thening follow the steps 370.
Judge whether not analyzed there are also forecast accident to be analyzed, if do not had there are also forecast accident to be analyzed
It is analyzed, then returns to step 340.
Step 370, according to the calculation of tidal current under each forecast accident, determine the power transmission network under the forecast accident
Safety in operation.
The technical solution of the present embodiment passes through the current measurement electricity according to boundary node on the basis of the above embodiments
Pressure and current measured power are calculated with the lesser net injecting power vector of voltage change, according to the net injecting power vector sum
Distribution network Equivalent Model establishes power flow equation of the power transmission network under the forecast accident of present analysis, need to only calculate primary net note
Enter vector power, can analyze to obtain the calculation of tidal current under all forecast accidents.In addition, matching in the embodiment of the present invention
Electric network Equivalent Model only just needs power distribution network control centre to be counted again when distribution network topological structure changes
It calculates, and distribution network topologies change number is not more than ten times according to actual electric network operating experience, in 24 hours, therefore, greatly
In most cases without recalculating distribution network Equivalent Model, reduce calculation amount compared with the existing technology, and reduce
Communication cost, it is contemplated that crash analysis result is more accurate, effectively avoids the false alarm being likely to occur in the prior art and false dismissal.
On the basis of the above embodiments, it according to the distribution network Equivalent Model and net injecting power vector, establishes
The power transmission network is also optional to include: before the power flow equation under the forecast accident of present analysis
If the forecast accident of present analysis is that generator is out of service, the output to the generator other than the generator
Power is adjusted, to meet Power Systems Constraints of Equilibrium in Contingency analysis, and according to generator adjusted
Output power calculates the injecting power of node where the generator in power transmission network；
If the forecast accident of present analysis is that route is out of service, according to the generator of node each in power transmission network
Output power and load power calculate the injecting power of the node；
According to the distribution network Equivalent Model and net injecting power vector, the power transmission network is established in the pre- of present analysis
Think the power flow equation under accident, comprising:
According to the injecting power of each node in the distribution network Equivalent Model, net injecting power vector sum power transmission network,
Establish power flow equation of the power transmission network under the forecast accident of present analysis.
Whether the anticipation for judging present analysis is that generator is out of service or route is out of service, is moved back if it is generator
It runs out, then adjusts the output power of remaining generator according to formula (5), to meet Power Systems in Contingency analysis
Constraints of Equilibrium.
Wherein, PGiIndicate generator out of service in preceding output power out of service, PGjAnd PGj,newIt respectively indicates
Remaining output power and output power adjusted for not having generator out of service to measure at current time, factor alphajTable
Show the generator power PG given in advance by power transmission network operations staffjRegulation coefficient.
According to the load power of node where the output power of generator adjusted and the generator, according to formula
(6) injecting power of the node is calculated:
ST,j=PGj,new-LDj (6)
Wherein, ST,jIndicate the injecting power of power transmission network arbitrary node j, LDjIndicate the load power of node j.
If the forecast accident of present analysis is that route is out of service, according to the output power in stage each in power transmission network
And load power, the injecting power of the node is calculated according to formula (7):
ST,j=PGj-LDj (7)
Wherein, ST,jIndicate the injecting power of power transmission network arbitrary node j, LDjIndicate the load power of node j, j indicates section
The serial number of point.
The injecting power in each stage in net injecting power vector sum power transmission network is substituted into power flow equation, that is, formula (4) to carry out
It solves, the voltage of the voltage of each node and boundary node in power transmission network, i.e. vector V is calculatedTWith vector VB, then will obtain
Power transmission network in each node voltage and boundary node voltage substitute into power equation, the anticipation in present analysis is calculated
The transimission power S of route under accidentTT、SBBAnd SBT, the voltage of each node, the voltage of boundary node and route in power transmission network
Transimission power, that is, calculation of tidal current.When generator is out of service, the output power of remaining generator is adjusted, is made
It meets Power Systems Constraints of Equilibrium in Contingency analysis, so that Contingency analysis result is more accurate.
Embodiment three
Fig. 4 is a kind of flow chart of the method for power transmission network Contingency analysis that the embodiment of the present invention three provides, this implementation
Example is a preferred embodiment on the basis of the above embodiments, and this method specifically comprises the following steps:
Step 410, power distribution network control centre calculates distribution network Equivalent Model according to Gaussian reduction, and by distribution network
Equivalent Model is sent to power transmission network control centre.
Distribution network Equivalent Model is as shown in Figure 2.The distribution network Equivalent Model includes at least one boundary node, institute
Stating boundary node is the node that power distribution network is connected with power transmission network, and each boundary node passes through equivalent shunt admittance ground connection, if side
Boundary's node is no less than two, then is connected between every two boundary node by equivalent series admittance.
Step 420, the distribution network Equivalent Model that power distribution network control centre sends is added to by power transmission network control centre
In the Contingency analysis model of power transmission network.
The current measurement voltage and current measured power for obtaining boundary node are being worked as according to formula (1) and formula (2) calculating
Equivalent power vector S of the distribution network Equivalent Model at boundary node under preceding measurement voltageD1,0, then counted according to formula (3)
Calculate the net injecting power vector of boundary node.
Step 430, it is assumed that power transmission network forecast accident to be analyzed is c, and power transmission network control centre judges forecast accident c for hair
Motor is out of service or route is out of service.
If forecast accident c is that generator is out of service, the output power of remaining generator is adjusted according to formula (5),
To meet Power Systems Constraints of Equilibrium in Contingency analysis；If forecast accident c is that route is out of service, do not have to adjust
It haircuts the output power of motor.
Step 440, net injecting power vector of the power transmission network control centre based on distribution network Equivalent Model and boundary node,
The power transmission network power flow equation at forecast accident c is established, and is solved.
The power transmission network power flow equation under forecast accident established is (4), in formula (4), vector STIn power transmission network
The injecting power S of arbitrary node jT,jIt is determined by formula (6) or formula (7).
Above-mentioned power transmission network power flow equation can be solved using Newton-Raphson method or Fast decomposition algorithms.
Step 450, whether in power transmission network have node voltage cross in advance given the operation upper bound or lower bound, and check if checking
Whether there is the transimission power of transmission line of electricity to cross given in advance the operation upper bound or lower bound.If there is such situation, then accuse
There is hidden danger of crossing the border in the alert operation of power transmission network at this time；Otherwise, show at above-mentioned forecast accident c, above-mentioned power transmission network power flow solutions are full
Foot operating range given in advance, system safety.
According to the solution of step 440, the node voltage V of power transmission network is obtainedTWith the voltage V of boundary nodeB, by power side
Journey calculates the transimission power S of the transmission line of electricity under above-mentioned forecast accident cTT、SBBAnd SBT, obtain the transmission of electricity under above-mentioned forecast accident c
Net power flow solutions.Check the transimission power of the node voltage and transmission line of electricity in the power transmission network power flow solutions.
Step 460, it checks whether and has had stepped through all forecast accidents to be analyzed, if it is, termination analysis；Otherwise
Step 430 is executed, analytical calculation is carried out to remaining forecast accident according to such step.
The technical solution of the present embodiment, take full advantage of distribution network Equivalent Model can approximation reflect that distribution network is opened up
The characteristics of flutterring variation carries out the influence of the distribution network change in topology in power transmission network forecast accident closely with distribution network equivalence
Seemingly, more accurate power transmission network Contingency analysis only need to be obtained with few calculating and communication cost as a result, providing
More accurate power grid security early warning, the final safety for improving operation of power networks.In addition, the distribution network in the present embodiment is equivalent
Model is only just recalculated when distribution network topological structure changes, and is run and passed through according to actual electric network
It tests, distribution network topologies change number is not more than ten times in 24 hours, therefore the step 410 of the present embodiment is in most of feelings
Without carrying out under condition, power transmission network control centre distribution network Equivalent Model at moment before is calculated since step 420,
In this case, compared to existing power transmission network Contingency analysis method, the present embodiment is not necessarily to additional communication cost.Needle
Electric system simulation to a variety of scales calculation shows that, the present embodiment is it is possible to prevente effectively from the false alarm being likely to occur in the prior art
And false dismissal.
Example IV
Fig. 5 is a kind of structural schematic diagram of the device for power transmission network Contingency analysis that the embodiment of the present invention four provides, this
The device of power transmission network Contingency analysis described in embodiment can be generally integrated in power transmission network control centre, as shown in figure 5, this
The device of power transmission network Contingency analysis described in embodiment includes: Equivalent Network model receiving module 510, Load flow calculation module
520 and safety in operation determination module 530.
Wherein, Equivalent Network model receiving module 510, for receiving the distribution network etc. that power distribution network control centre sends
It is worth model；
Load flow calculation module 520, for the distribution network Equivalent Model to be added to the Contingency analysis of power transmission network
In model, Load flow calculation is carried out, wherein the Contingency analysis model includes multiple forecast accidents；
Safety in operation determination module 530, for determining described defeated according to the calculation of tidal current under each forecast accident
Safety in operation of the power grid under the forecast accident.
Optionally, the distribution network Equivalent Model includes at least one boundary node, and the boundary node is power distribution network
The node being connected with power transmission network, each boundary node passes through equivalent shunt admittance ground connection, if boundary node is no less than two,
It is connected between every two boundary node by equivalent series admittance.
Optionally, the Load flow calculation module includes:
Measured value acquiring unit, for obtaining the current measurement voltage and current measured power of each boundary node；
Net injecting power computing unit, for according to the current measurement voltage, current measured power, equivalent shunt admittance
And equivalent series admittance, calculate the net injecting power vector of the boundary node in the distribution network Equivalent Model；
Forecast accident acquiring unit, for obtaining a forecast accident to be analyzed in power transmission network, as present analysis
Forecast accident；
Load flow calculation unit, for establishing described defeated according to the distribution network Equivalent Model and net injecting power vector
Power flow equation of the power grid under the forecast accident of present analysis, and solve, obtain calculation of tidal current；
Cycling element, for executing the operation for obtaining a forecast accident to be analyzed in power transmission network, until having analyzed institute
There is forecast accident to be analyzed.
Optionally, the Load flow calculation module further include:
Injecting power computing unit, for establishing according to the distribution network Equivalent Model and net injecting power vector
The power transmission network is before the power flow equation under the forecast accident of present analysis, if the forecast accident of present analysis is generator
It is out of service, then the output power of the generator other than the generator is adjusted, to meet electric power in Contingency analysis
System power Constraints of Equilibrium, and according to the output power of generator adjusted, calculate node where the generator in power transmission network
Injecting power；If the forecast accident of present analysis is that route is out of service, according to the power generation of node each in power transmission network
The output power and load power of machine, calculate the injecting power of the node；
The Load flow calculation unit is specifically used for:
According to the injecting power of each node in the distribution network Equivalent Model, net injecting power vector sum power transmission network,
Power flow equation of the power transmission network under the forecast accident of present analysis is established, and is solved, calculation of tidal current is obtained.
Optionally, the net injecting power computing unit is specifically used for:
According to the current measurement voltage, equivalent shunt admittance and equivalent series admittance, it is equivalent to calculate the distribution network
Equivalent power vector of the model at the boundary node；
According to the current measured power and equivalent power vector, the boundary section in the distribution network Equivalent Model is calculated
The net injecting power vector of point.
Optionally, the safety in operation determination module is specifically used for:
According to the calculation of tidal current under each forecast accident, judge whether have in the power transmission network under the forecast accident
Whether node voltage is more than the default working voltage upper bound or default working voltage lower bound, and have the transimission power of transmission line of electricity super
Cross default operation upper power bound or default operation lower power bound；
If there is node voltage more than the default working voltage upper bound or default working voltage lower bound in the power transmission network, and/
Or, the transimission power for having transmission line of electricity is more than default operation upper power bound or default operation lower power bound, then power transmission network fortune is alerted
There is hidden danger of crossing the border in row.
Power transmission network anticipation provided by any embodiment of the invention can be performed in the device of above-mentioned power transmission network Contingency analysis
The method of crash analysis has the corresponding functional module of execution method and beneficial effect.Not detailed description in the present embodiment
Technical detail, reference can be made to the method for the power transmission network Contingency analysis that any embodiment of that present invention provides.
Note that the above is only a better embodiment of the present invention and the applied technical principle.It will be appreciated by those skilled in the art that
The invention is not limited to the specific embodiments described herein, be able to carry out for a person skilled in the art it is various it is apparent variation,
It readjusts and substitutes without departing from protection scope of the present invention.Therefore, although being carried out by above embodiments to the present invention
It is described in further detail, but the present invention is not limited to the above embodiments only, without departing from the inventive concept, also
It may include more other equivalent embodiments, and the scope of the invention is determined by the scope of the appended claims.
Claims (8)
1. a kind of method of power transmission network Contingency analysis, which is characterized in that the described method includes:
Receive the distribution network Equivalent Model that power distribution network control centre sends；
The distribution network Equivalent Model is added in the Contingency analysis model of power transmission network, Load flow calculation is carried out, wherein
The Contingency analysis model includes multiple forecast accidents；
According to the calculation of tidal current under each forecast accident, operational safety of the power transmission network under the forecast accident is determined
Property；
The distribution network Equivalent Model includes at least one boundary node, and the boundary node is that power distribution network is connected with power transmission network
The node connect, each boundary node pass through equivalent shunt admittance ground connection, if boundary node is no less than two, every two boundary section
It is connected between point by equivalent series admittance；
The distribution network Equivalent Model is added in the Contingency analysis model of power transmission network, Load flow calculation is carried out, comprising:
Obtain the current measurement voltage and current measured power of each boundary node；
According to the current measurement voltage, current measured power, equivalent shunt admittance and equivalent series admittance, match described in calculating
The net injecting power vector of boundary node in electric network Equivalent Model；
Obtain a forecast accident to be analyzed in power transmission network, the forecast accident as present analysis；
According to the distribution network Equivalent Model and net injecting power vector, the power transmission network is established in the anticipation thing of present analysis
Power flow equation under therefore, and solve, obtain calculation of tidal current；
The operation for obtaining a forecast accident to be analyzed in power transmission network is executed, until having analyzed all anticipation things to be analyzed
Therefore.
2. the method according to claim 1, wherein according to the distribution network Equivalent Model and net injection function
Rate vector establishes the power transmission network before the power flow equation under the forecast accident of present analysis, further includes:
If the forecast accident of present analysis is that generator is out of service, to the output power of the generator other than the generator
It is adjusted, to meet Power Systems Constraints of Equilibrium in Contingency analysis, and according to the output of generator adjusted
Power calculates the injecting power of node where the generator in power transmission network；
If the forecast accident of present analysis is that route is out of service, according to the output of the generator of node each in power transmission network
Power and load power calculate the injecting power of the node；
According to the distribution network Equivalent Model and net injecting power vector, the power transmission network is established in the anticipation thing of present analysis
Power flow equation under therefore, comprising:
According to the injecting power of each node in the distribution network Equivalent Model, net injecting power vector sum power transmission network, establish
Power flow equation of the power transmission network under the forecast accident of present analysis.
3. according to the method described in claim 2, it is characterized in that, according to the current measurement voltage, current measured power, etc.
Be worth shunt admittance and equivalent series admittance, calculate the net injecting power of the boundary node in the distribution network Equivalent Model to
Amount, comprising:
According to the current measurement voltage, equivalent shunt admittance and equivalent series admittance, the distribution network Equivalent Model is calculated
Equivalent power vector at the boundary node；
According to the current measured power and equivalent power vector, the boundary node in the distribution network Equivalent Model is calculated
Net injecting power vector.
4. method according to claim 1 to 3, which is characterized in that according to the Load flow calculation knot under each forecast accident
Fruit determines safety in operation of the power transmission network under the forecast accident, comprising:
According to the calculation of tidal current under each forecast accident, judge whether there is node in the power transmission network under the forecast accident
Voltage is more than the default working voltage upper bound or default working voltage lower bound, and it is more than pre- for whether having the transimission power of transmission line of electricity
If running upper power bound or default operation lower power bound；
If there is node voltage more than the default working voltage upper bound or default working voltage lower bound in the power transmission network, and/or, have
The transimission power of transmission line of electricity is more than default operation upper power bound or default operation lower power bound, then alerts power transmission network operation and exist
It crosses the border hidden danger.
5. a kind of device of power transmission network Contingency analysis, which is characterized in that described device includes:
Equivalent Network model receiving module, the distribution network Equivalent Model sent for receiving power distribution network control centre；
Load flow calculation module, for the distribution network Equivalent Model to be added in the Contingency analysis model of power transmission network,
Carry out Load flow calculation, wherein the Contingency analysis model includes multiple forecast accidents；
Safety in operation determination module, for determining that the power transmission network exists according to the calculation of tidal current under each forecast accident
Safety in operation under the forecast accident；
The distribution network Equivalent Model includes at least one boundary node, and the boundary node is that power distribution network is connected with power transmission network
The node connect, each boundary node pass through equivalent shunt admittance ground connection, if boundary node is no less than two, every two boundary section
It is connected between point by equivalent series admittance；
The Load flow calculation module includes:
Measured value acquiring unit, for obtaining the current measurement voltage and current measured power of each boundary node；
Net injecting power computing unit, for according to the current measurement voltage, current measured power, equivalence shunt admittance and
Equivalent series admittance calculates the net injecting power vector of the boundary node in the distribution network Equivalent Model；
Forecast accident acquiring unit, the anticipation for obtaining a forecast accident to be analyzed in power transmission network, as present analysis
Accident；
Load flow calculation unit, for establishing the power transmission network according to the distribution network Equivalent Model and net injecting power vector
Power flow equation under the forecast accident of present analysis, and solve, obtain calculation of tidal current；
Cycling element is needed for executing the operation for obtaining a forecast accident to be analyzed in power transmission network until having analyzed
The forecast accident of analysis.
6. device according to claim 5, which is characterized in that the Load flow calculation module further include:
Injecting power computing unit, for according to the distribution network Equivalent Model and net injecting power vector, described in foundation
Power transmission network is before the power flow equation under the forecast accident of present analysis, if the forecast accident of present analysis exits for generator
Operation, then be adjusted the output power of the generator other than the generator, to meet electric system in Contingency analysis
Power-balance constraint, and according to the output power of generator adjusted, calculate the note of node where the generator in power transmission network
Enter power；If the forecast accident of present analysis is that route is out of service, according to the generator of node each in power transmission network
Output power and load power calculate the injecting power of the node；
The Load flow calculation unit is specifically used for:
According to the injecting power of each node in the distribution network Equivalent Model, net injecting power vector sum power transmission network, establish
Power flow equation of the power transmission network under the forecast accident of present analysis, and solve, obtain calculation of tidal current.
7. device according to claim 5, which is characterized in that the net injecting power computing unit is specifically used for:
According to the current measurement voltage, equivalent shunt admittance and equivalent series admittance, the distribution network Equivalent Model is calculated
Equivalent power vector at the boundary node；
According to the current measured power and equivalent power vector, the boundary node in the distribution network Equivalent Model is calculated
Net injecting power vector.
8. according to any device of claim 5-7, which is characterized in that the safety in operation determination module is specifically used
In:
According to the calculation of tidal current under each forecast accident, judge whether there is node in the power transmission network under the forecast accident
Voltage is more than the default working voltage upper bound or default working voltage lower bound, and it is more than pre- for whether having the transimission power of transmission line of electricity
If running upper power bound or default operation lower power bound；
If there is node voltage more than the default working voltage upper bound or default working voltage lower bound in the power transmission network, and/or, have
The transimission power of transmission line of electricity is more than default operation upper power bound or default operation lower power bound, then alerts power transmission network operation and exist
It crosses the border hidden danger.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710403029.8A CN106972490B (en) | 2017-06-01 | 2017-06-01 | The method and device of power transmission network Contingency analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710403029.8A CN106972490B (en) | 2017-06-01 | 2017-06-01 | The method and device of power transmission network Contingency analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106972490A CN106972490A (en) | 2017-07-21 |
| CN106972490B true CN106972490B (en) | 2019-10-22 |
Family
ID=59327462
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710403029.8A Active CN106972490B (en) | 2017-06-01 | 2017-06-01 | The method and device of power transmission network Contingency analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106972490B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1967960A (en) * | 2006-09-29 | 2007-05-23 | 天津大学 | Coherency generating sets of electric power system and polymeric method of asynchronous motor group parameter |
| CN101282039A (en) * | 2008-05-09 | 2008-10-08 | 天津大学 | Method of polymerization of induction motor load dynamic parameter taking account of slip and load factor |
| CN103441497A (en) * | 2013-09-05 | 2013-12-11 | 重庆大学 | Static equivalence method of external network in interconnected power grid |
| CN103997036A (en) * | 2014-05-16 | 2014-08-20 | 湖南智仁科技有限公司 | Motor aggregation method considering system stability and electromagnetic property |
| CN104766142A (en) * | 2015-04-21 | 2015-07-08 | 河海大学 | Transient stability constraint optimal power flow computation method based on EEAC and trace sensitivity |
| CN106099906A (en) * | 2016-05-13 | 2016-11-09 | 国电南瑞科技股份有限公司 | A kind of power system safety and stability online integrated campaign method |
-
2017
- 2017-06-01 CN CN201710403029.8A patent/CN106972490B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1967960A (en) * | 2006-09-29 | 2007-05-23 | 天津大学 | Coherency generating sets of electric power system and polymeric method of asynchronous motor group parameter |
| CN101282039A (en) * | 2008-05-09 | 2008-10-08 | 天津大学 | Method of polymerization of induction motor load dynamic parameter taking account of slip and load factor |
| CN103441497A (en) * | 2013-09-05 | 2013-12-11 | 重庆大学 | Static equivalence method of external network in interconnected power grid |
| CN103997036A (en) * | 2014-05-16 | 2014-08-20 | 湖南智仁科技有限公司 | Motor aggregation method considering system stability and electromagnetic property |
| CN104766142A (en) * | 2015-04-21 | 2015-07-08 | 河海大学 | Transient stability constraint optimal power flow computation method based on EEAC and trace sensitivity |
| CN106099906A (en) * | 2016-05-13 | 2016-11-09 | 国电南瑞科技股份有限公司 | A kind of power system safety and stability online integrated campaign method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106972490A (en) | 2017-07-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105186468B (en) | A kind of power distribution network directional over-current protection on-line tuning method and system of distributed power source access | |
| CN100590947C (en) | Large electric network on-line preventing control method based on static state and transient safety steady mode | |
| CN105576660B (en) | Local distribution network in-circuit emulation analytical equipment and system | |
| CN106469908A (en) | The static security synthesis aid decision-making method that a kind of multiple faults is coordinated | |
| US20140309952A1 (en) | Method for Analyzing Faults in Ungrounded Power Distribution Systems | |
| CN106058876B (en) | Consider the dynamic reactive Site Selection analysis method and system of Enhancement of Transient Voltage Stability | |
| CN102982230B (en) | A kind of short circuit current based on node impedance sensitivity exceeds standard aid decision-making method | |
| CN103279656B (en) | Based on the power distribution network net capability computing method of trend | |
| CN103971292B (en) | Fault accurate identification method based on real-time intra-station graph mapping | |
| CN106340907A (en) | Power system security and stability control strategy determining method and device | |
| CN106340897A (en) | Method for determining allowed capacity of distribution network distributed power supply | |
| EP2770600B1 (en) | Method and system for determining power consumption of loads | |
| CN102104252A (en) | Power system dynamic equivalence method suitable for electromagnetic transient analysis | |
| CN106781929A (en) | A kind of power distribution network relay protection training set and method | |
| CN106208041A (en) | A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis | |
| CN106096806A (en) | A kind of subjectively-based multi-infeed HVDC synchronizing point appraisal procedure integrated with objective weight | |
| CN104715144A (en) | Power distribution network closed loop power flow simulation algorithm based on real-time running data | |
| Li et al. | An integrated online dynamic security assessment system for improved situational awareness and economic operation | |
| Overbye et al. | Techniques for maintaining situational awareness during large-scale electric grid simulations | |
| CN106972490B (en) | The method and device of power transmission network Contingency analysis | |
| CN106022970A (en) | Active power distribution network measurement configuration method considering distributed power sources | |
| CN105911321A (en) | Electrical parameter transmission relationship-based T-type connection line electricity stealing identification for power distribution network | |
| CN109638871A (en) | Consider the extensive alternating current-direct current combined hybrid system major network division methods of wind power integration | |
| CN107219420B (en) | Simulation method and system for short-circuit current shunt coefficient of converter station | |
| CN103647284B (en) | A kind of voltage stabilization Forecasting Methodology of discontinuity surface problem when solving list |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221213 Address after: 518000 2nd floor, building a, Tsinghua campus, Shenzhen University Town, Xili street, Nanshan District, Shenzhen City, Guangdong Province Patentee after: Tsinghua Shenzhen International Graduate School Address before: 518000 Nanshan Zhiyuan 1001, Xue Yuan Avenue, Nanshan District, Shenzhen, Guangdong. Patentee before: TSINGHUA-BERKELEY SHENZHEN INSTITUTE |