CN106558880B - The aid decision-making method of branch overload problem in a kind of transient stability analysis - Google Patents
The aid decision-making method of branch overload problem in a kind of transient stability analysis Download PDFInfo
- Publication number
- CN106558880B CN106558880B CN201510639014.2A CN201510639014A CN106558880B CN 106558880 B CN106558880 B CN 106558880B CN 201510639014 A CN201510639014 A CN 201510639014A CN 106558880 B CN106558880 B CN 106558880B
- Authority
- CN
- China
- Prior art keywords
- branch
- unit
- overload
- load
- current
- 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
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000001052 transient effect Effects 0.000 title claims abstract description 39
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000000243 solution Substances 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 3
- 230000033228 biological regulation Effects 0.000 abstract description 4
- 230000005284 excitation Effects 0.000 abstract description 3
- 230000003068 static effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000005283 ground state Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Abstract
The present invention relates to a kind of aid decision-making methods of branch overload problem in transient stability analysis, comprising: whether branch overloads during judging transient stability analysis；If there is branch overload situations, inject unit virtual current one by one at unit or load bus, then the voltage of unit or load bus is calculated, and then computation overload branch current, to inject unit virtual current at different units or load bus, and unit and load are ranked up in the size that overload branch generates electric current；According to the overload quantity of overload branch and unit and load ordering scenario, active/idle adjustment amount of unit and load is determined；Transient stability analysis is checked.Technical solution provided by the invention solves the problems, such as that branch overload aid decision calculates in transient stability analysis of power system, due to considering the influence of the models such as excitation, speed regulation, load, so calculated result is more accurate credible compared with steady method.
Description
Technical field
The present invention relates to the aid decision-making methods in a kind of electric system, and in particular to branch in a kind of transient stability analysis
The aid decision-making method of overload problem.
Background technique
The safe and stable operation of holding electric system is the elementary object and basic principle of electric system, is transported in electric system
In the work such as professional etiquette is drawn, mode is analyzed, on-line early warning, it is necessary to control the object run mode of system and keep security and stability
In the range of.Thus, when branch overload occurs in ground state and forecast failure for system, it is necessary to be determined by branch overload auxiliary
Plan come adjust the method for operation with eliminate overload.
Traditional branch overload aid decision-making method is all based on stable state calculating, mostly uses DC power flow algorithm, sensitivity
Method, penalty method etc..Traditional branch overload aid decision is just for the mistake occurred in ground state trend or static security analysis
The method of load problem, method and static security analysis is directly associated, and for different static security analysis methods, needs to adopt
With corresponding aid decision-making method.Mainly include following several thinkings:
1) accurate tidal current computing method is used.
This method uses accurate power flow algorithm, is generally solved using Newton-Raphson approach.Static security analysis
In delete corresponding component models one by one, carry out Load flow calculation, and then analyze overload situations.It the case where for overloading, can
To be ranked up by the method for sensitivity to unit, load, then each group measure is combined and is verified, completed auxiliary
Help decision.
2) method of DC power flow is used.
This method is fairly simple, and calculating speed is fast, but calculates poor accuracy, so do not praised highly by everybody,
But one's early years computer technology and it is undeveloped when this method it is popular.
3) other methods
Also there is the branch overload aid decision-making method based on penalty method, but because computational accuracy is not high, also do not recognized by everybody
It can.Then make in addition, also having and calculating branch transfer distribution factor based on DC power flow algorithm for being ranked up to unit, load
The method verified with accurate trend method.This method combines the rapidity of DC power flow and the accuracy of accurate trend,
Have received the pretty good effect of comparison.
But above-mentioned all methods are all based only on steady-state model, there is no consider generator, excitation, speed regulation, bear
The transient characterisitics of the models such as lotus, direct current are not fully suitble to the overload problem occurred during transient stability analysis.
Summary of the invention
To solve above-mentioned deficiency of the prior art, the object of the present invention is to provide prop up to pass by a kind of transient stability analysis
The aid decision-making method of load problem is calculated this method solve branch overload aid decision in transient stability analysis of power system
Problem.
The purpose of the present invention is adopt the following technical solutions realization:
The present invention provides a kind of aid decision-making method of branch overload problem in transient stability analysis, and improvements exist
In the method includes the following steps:
Step 1: analyzing forecast failure using transient stability program PSD-BPA, analyzes the mistake of each branch after failure
Carry situation；Whether branch overloads during judging transient stability analysis using the power or current data of plateau, and sentences
The disconnected duration reaches 100 cycles or more；
Step 2: if there is branch overload situations, then using having been generated in transient stability program and carried out LU decomposition
The admittance matrix of (LU Decomposition, matrix decomposition), injects unit virtual current one by one at unit or load bus,
Then the voltage of unit or load bus, and then computation overload branch current are calculated, at different units or load bus
Unit virtual current is injected, and unit and load are ranked up in the size that overload branch generates electric current；
Step 3: according to overload branch overload quantity and step 2 in unit and load ordering scenario, determine unit and
Active/idle adjustment amount of load；
Step 4: transient stability analysis is checked: power flow stability data adjusted carry out Load flow calculation and transient stability point
Analysis judges whether to achieve the purpose that eliminate overload by the result of stability analysis, if do not eliminated, return to step 3 into
Row calculates, until maximum adjustment number is eliminated or reached to overload.
Further, in the step 1, the judgment method of plateau are as follows: take the power or current value shape of a period of time
At curve, then curve is fitted with PRONY method, amplitude, the frequency, damping information of curve are determined, when frequency is less than
When given value, i.e., it is believed that entering plateau, given value is adjusted according to actual electric network.
Further, the step 2 includes the following steps:
(1) if branch k overloads, unit virtual current is injected in unit or load bus i, forms equation group,
Matrix form is as follows:
Wherein: Y is admittance battle array, has been broken down into the form of Y=LDU, and L is the inferior triangular flap after decomposing, and D is to decompose
Diagonal element battle array afterwards, U are the upper triangular matrix after decomposing；For each node voltage column vector；For virtual Injection Current column vector,
It is 1 in the corresponding position upper value of node i, and other positions value is 0；Two end nodes of branch k are respectively respectively m and l；
(2) equation group in solution procedure (1), obtains column voltage vectorTake the voltage at the both ends branch kWithAnd branch
Road impedance Zk, calculate the electric current for then flowing through branch kAre as follows:
(3) i ∈ [1, n] is taken, carries out step (1) and (2) one by one, is found out pair for all units or load bus i
It answersAnd withSize be ranked up；Wherein n is the number of unit and load bus.
Further, in the step 3, the Adjustment principle of active/idle adjustment amount of unit and load includes:
1) it is adjusted according to the sequence of unit in step 2 and load；
2) target is adjusted are as follows:
Wherein: Δ PiFor the active power adjustment amount of node i, Δ QiFor the idle adjustment amount of node i, UiFor the voltage of node i,
Iover=Inow-Ilim it, InowFor the current value for adjusting preceding branch k, Ilim itFor the rated current of branch k, j is the imaginary number of plural number
Unit.
Further, in the step 4, maximum adjustment number is set as 5 times.
The excellent effect that technical solution provided by the invention has is:
(1) in transient stability analysis provided by the invention branch overload problem aid decision-making method, for transient stability
The branch overload problem occurred in analysis, creates integrated solution.Judged during multilayer output feedback network
It carries, calculated if overload by temporarily steady program and exports the virtual injection electricity of each adjustable node (generator and load bus) injection
The current value of branch is overloaded when stream, and power output and load adjustment are then carried out according to the sequence of above-mentioned current value, finally use transient stability
Program is checked.
(2) fully considered the transient characterisitics of the models such as generator, excitation, speed regulation, load, direct current, result relative to
Conventional method is more accurate, is the important foundation of aid decision-making method of the present invention.
Detailed description of the invention
Fig. 1 is the flow chart of the aid decision-making method of branch overload problem in transient stability analysis provided by the invention.
Specific embodiment
Specific embodiments of the present invention will be described in further detail with reference to the accompanying drawing.
The following description and drawings fully show specific embodiments of the present invention, to enable those skilled in the art to
Practice them.Other embodiments may include structure, logic, it is electrical, process and other change.Embodiment
Only represent possible variation.Unless explicitly requested, otherwise individual component and function are optional, and the sequence operated can be with
Variation.The part of some embodiments and feature can be included in or replace part and the feature of other embodiments.This hair
The range of bright embodiment includes equivalent obtained by the entire scope of claims and all of claims
Object.Herein, these embodiments of the invention can individually or generally be indicated that this is only with term " invention "
For convenience, and if in fact disclosing the invention more than one, the range for being not meant to automatically limit the application is to appoint
What single invention or inventive concept.
The present invention provides a kind of aid decision-making method of branch overload problem in transient stability analysis, and technical solution is overall
Thinking are as follows: branch overloading analysis is carried out to the system after failure based on transient stability program, if there is branch overload situations, then
Each unit and load are calculated to the current contribution of out-of-limit branch, to be ranked up to it, to measure using transient stability program
It after being combined, is checked using multilayer output feedback network, branch overload problem in transient stability analysis provided by the invention
Aid decision-making method flow chart as shown in Figure 1, the specific implementation steps are as follows:
Step 1: analyzing failure using transient stability program, analyzes the overload situations of each branch after failure.Judgement
The branch overload stage more stable using branch power after failure or electric current in transient process, and the duration judged is necessary
Reach 100 cycles or more.
Wherein: judge branch power or electric current more stable stage, it can be using the power or electric current for taking a period of time to put
It is worth forming curves, then curve is fitted with PRONY method, determines amplitude, the frequency, damping information of curve.Work as frequency
When less than given value, i.e., it is believed that entering the stabilization sub stage, given value can be adjusted according to actual electric network.
Step 2: if there is branch overload situations, then using having been generated in transient stability program and carried out LU decomposition
Admittance matrix, inject unit virtual current one by one at each unit or load, then calculate the voltage of each node, thus
The electric current for overloading branch at this time is calculated, then to inject unit virtual current at different units or load in overload branch
The size for the electric current that road generates is ranked up unit and load.Its detailed process is as follows:
(1) it if branch k (its two end node is respectively m and l) overloads, is injected in unit (or load) node i single
Position virtual current, then form equation group, matrix form is as follows:
Wherein: Y is admittance battle array, has been broken down into the form of Y=LDU, and L is the inferior triangular flap after decomposing, and D is to decompose
Diagonal element battle array afterwards, U are the upper triangular matrix after decomposing；For each node voltage column vector；For virtual Injection Current column vector,
It is only 1 in the corresponding position upper value of node i, and other positions value is 0.
(2) equation group in (1) is solved, column voltage vector is obtainedTake the voltage at the both ends branch kWithAnd branch roadlock
Anti- Zk, the electric current for then flowing through branch k can be calculatedAre as follows:
(3) i ∈ [1, n] is taken, (1) (2) step is carried out one by one, is found out for all i correspondingAnd with
Size be ranked up.Wherein n is the number of unit and load bus.
Step 3: according to the overload quantity I of branch koverAnd unit and load ordering scenario in step 2, determine unit and
Active/idle adjustment amount of load.Adjustment process needs to consider problem:
1) it is adjusted according to the sequence in step 2；
2) target is adjusted are as follows:
Wherein: Δ PiFor the active power adjustment amount of node i, Δ QiFor the idle adjustment amount of node i, UiFor the voltage of node i,
Iover=Inow-Ilim it, InowFor the current value for adjusting preceding branch k, Ilim itFor the rated current of branch k, j is the imaginary number of plural number
Unit.
Step 4: data adjusted carry out Load flow calculation, and then stability analysis is by the result judgement of stability analysis
It is no to achieve the purpose that eliminate overload, if do not eliminated, reenter step 3 and calculated, until overload is eliminated or reached
Maximum adjustment number.Maximum adjustment number can be usually set as 5 times.
Technical solution provided by the invention solves branch overload aid decision in transient stability analysis of power system and calculates
The problem of, propose the branch overload aid decision-making method based on transient stability analysis.This method utilizes replacing in circuit theory
For theorem and principle of stacking, first using the admittance battle array that LU decomposition has been formed and performed during multilayer output feedback network, calculate
The magnitude of current of branch is overloaded when each unit unitary current injects；Secondly, more being limited the quantity according to the electric current of overload branch, calculate each
The equivalent current of unit injects adjustment amount；Then adjustment amount is injected according to equivalent current and calculates corresponding power adjustment；Most
It is verified afterwards using transient stability program.Branch overload aid decision-making method based on transient stability analysis is encouraged due to considering
The influence of the models such as magnetic, speed regulation, load, so calculated result is more accurate credible compared with steady method.
The above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, although referring to above-described embodiment pair
The present invention is described in detail, those of ordinary skill in the art still can to a specific embodiment of the invention into
Row modification perhaps equivalent replacement these without departing from any modification of spirit and scope of the invention or equivalent replacement, applying
Within pending claims of the invention.
Claims (5)
1. the aid decision-making method of branch overload problem in a kind of transient stability analysis, which is characterized in that the method includes under
State step:
Step 1: forecast failure is analyzed with transient stability program PSD-BPA, and analyzes the overload situations of each branch after failure；Benefit
Whether branch overloads during judging transient stability analysis with the power of plateau or current data, and judges its duration
Reach 100 cycles or more；
Step 2: it is overloaded if there is branch, then utilizes the admittance square that LU decomposition has been generated and carried out in transient stability program
Battle array, injects unit virtual current one by one at unit or load bus, then calculates the voltage of unit or load bus, Jin Erji
Overload branch current is calculated, to inject unit virtual current at different units or load bus, generates electric current by overload branch
Size sort to unit and load；
Step 3: according to the unit and load ordering scenario in the overload quantity of overload branch and step 2, unit and load are determined
Active/idle adjustment amount；
Step 4: transient stability analysis is checked: power flow stability data adjusted carry out Load flow calculation and transient stability analysis, lead to
The result for crossing stability analysis judges whether to achieve the purpose that eliminate overload, if do not eliminated, returns to step 3 and counted
It calculates, until maximum adjustment number is eliminated or reached to overload.
2. aid decision-making method as described in claim 1, which is characterized in that in the step 1, the judgement side of plateau
Method includes: the power or current value forming curves for taking a period of time, is then fitted with PRONY method to curve, is determined bent
The amplitude of line, frequency, damping information, when frequency is less than given value, i.e., it is believed that entering plateau, given value is according to reality
Border power grid adjusts.
3. aid decision-making method as described in claim 1, which is characterized in that the step 2 includes the following steps:
(1) if branch k overloads, unit virtual current is injected in unit or load bus i, forms equation group, matrix
Form is as follows:
Wherein: Y is admittance battle array, has been broken down into the form of Y=LDU, and L is the inferior triangular flap after decomposing, and D is after decomposing
Diagonal element battle array, U are the upper triangular matrix after decomposing；For each node voltage column vector；For virtual Injection Current column vector,
The corresponding position upper value of node i is 1, and other positions value is 0；Two end nodes of branch k are respectively respectively m and l；
(2) equation group in solution procedure (1), obtains column voltage vectorTake the voltage at the both ends branch kWithAnd branch roadlock
Anti- Zk, calculate the electric current for then flowing through branch kAre as follows:
(3) i ∈ [1, n] is taken, step (1) and (2) is carried out one by one, is found out for all units or load bus i correspondingAnd withSize be ranked up；Wherein n is the number of unit and load bus.
4. aid decision-making method as described in claim 1, which is characterized in that in the step 3, unit and load it is active/
The Adjustment principle of idle adjustment amount includes:
1) it is adjusted according to the sequence of unit in step 2 and load；
2) target is adjusted are as follows:
Wherein: Δ PiFor the active power adjustment amount of node i, Δ QiFor the idle adjustment amount of node i, UiFor the voltage of node i, Iover=
Inow-Ilimit, InowFor the current value for adjusting preceding branch k, IlimitFor the rated current of branch k, j is the imaginary unit of plural number；The electric current of branch k is flowed through in expression.
5. aid decision-making method as described in claim 1, which is characterized in that in the step 4, maximum adjustment number is set as 5
It is secondary.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510639014.2A CN106558880B (en) | 2015-09-30 | 2015-09-30 | The aid decision-making method of branch overload problem in a kind of transient stability analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510639014.2A CN106558880B (en) | 2015-09-30 | 2015-09-30 | The aid decision-making method of branch overload problem in a kind of transient stability analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106558880A CN106558880A (en) | 2017-04-05 |
| CN106558880B true CN106558880B (en) | 2019-09-27 |
Family
ID=58417352
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510639014.2A Active CN106558880B (en) | 2015-09-30 | 2015-09-30 | The aid decision-making method of branch overload problem in a kind of transient stability analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106558880B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102214922A (en) * | 2011-06-27 | 2011-10-12 | 山东电力研究院 | Evaluation system of power network planning scheme |
| CN102968062A (en) * | 2012-11-22 | 2013-03-13 | 广东电网公司电力调度控制中心 | Method for implementing transient stability control batch computing based on PSD-BPA |
| CN103208804A (en) * | 2012-12-31 | 2013-07-17 | 武汉大学 | Branch voltage stability analysis based selection method for reactive power compensation point of microgrid |
| CN103606921A (en) * | 2013-12-10 | 2014-02-26 | 上海交通大学 | Transient voltage stabilization judgment method based on real-time generalized Thevenin equivalence |
-
2015
- 2015-09-30 CN CN201510639014.2A patent/CN106558880B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102214922A (en) * | 2011-06-27 | 2011-10-12 | 山东电力研究院 | Evaluation system of power network planning scheme |
| CN102968062A (en) * | 2012-11-22 | 2013-03-13 | 广东电网公司电力调度控制中心 | Method for implementing transient stability control batch computing based on PSD-BPA |
| CN103208804A (en) * | 2012-12-31 | 2013-07-17 | 武汉大学 | Branch voltage stability analysis based selection method for reactive power compensation point of microgrid |
| CN103606921A (en) * | 2013-12-10 | 2014-02-26 | 上海交通大学 | Transient voltage stabilization judgment method based on real-time generalized Thevenin equivalence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106558880A (en) | 2017-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103020853B (en) | Method for checking short-term trade plan safety | |
| CN106329523B (en) | Consider probabilistic active power distribution network intelligence Sofe Switch robust Optimization Modeling method | |
| Arboleya et al. | BFS algorithm for voltage-constrained meshed DC traction networks with nonsmooth voltage-dependent loads and generators | |
| CN105811452B (en) | Power distribution network distributed photovoltaic limit permeability analysis method and device | |
| CN107565582B (en) | The online method of discrimination of bulk power grid steady stability dominant pattern | |
| CN104156542B (en) | It is a kind of based on the active distribution system Simulation of stability method implicitly projected | |
| Hou et al. | Cluster computing-based trajectory sensitivity analysis application to the WECC system | |
| CN109066694A (en) | Multiple target tide optimization method containing the electric system of flow controller between line | |
| CN109462257A (en) | It is a kind of meter and multiple random variable Network Voltage Stability sensitivity discrimination method | |
| CN109494719A (en) | A kind of mesolow mixing power distribution network stratification impedance analysis method | |
| CN104716646A (en) | Node coupling degree analysis method based on injection currents | |
| CN109995023A (en) | The judgment method and system of stability of power system based on spectral clustering region division | |
| CN103678798B (en) | It is a kind of for the electro-magnetic transient real-time emulation method containing distributed power distribution network | |
| de Oliveira et al. | Decentralized three-phase distribution system static state estimation based on phasor measurement units | |
| CN105119279B (en) | A kind of distributed power source planing method and its system | |
| CN102769292B (en) | Method for enhancing voltage safety margin of wind farm | |
| Huang et al. | A practical method of transient stability analysis of stochastic power systems based on EEAC | |
| CN107994562A (en) | A kind of stability Design method of meter and the bavin storage type Vehicular power system of load characteristic | |
| CN106877309A (en) | Stability judgment method for direct-current micro-grid in master-slave control mode | |
| CN111682546B (en) | DC power flow improvement algorithm based on sensitivity analysis | |
| CN109638871A (en) | Consider the extensive alternating current-direct current combined hybrid system major network division methods of wind power integration | |
| Gupta et al. | Power system network equivalents: Key issues and challenges | |
| CN108988385A (en) | A kind of grid-connected optimal location choosing method of distributed generation resource | |
| CN106558880B (en) | The aid decision-making method of branch overload problem in a kind of transient stability analysis | |
| CN103384070A (en) | Power quality analytical calculation method based on PSASP |
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 |