CN105305504B - Transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method - Google Patents
Transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method Download PDFInfo
- Publication number
- CN105305504B CN105305504B CN201510856649.8A CN201510856649A CN105305504B CN 105305504 B CN105305504 B CN 105305504B CN 201510856649 A CN201510856649 A CN 201510856649A CN 105305504 B CN105305504 B CN 105305504B
- Authority
- CN
- China
- Prior art keywords
- branch road
- unit
- branch
- mrow
- power
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000005540 biological transmission Effects 0.000 title claims abstract description 9
- 230000035945 sensitivity Effects 0.000 claims abstract description 28
- 238000012546 transfer Methods 0.000 claims abstract description 13
- 230000005283 ground state Effects 0.000 claims abstract description 6
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 230000005611 electricity Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention belongs to Power System and its Automation technical field, and in particular to a kind of transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method.Characterized in that, first under normal operating mode, a kind of typical load pattern is taken, forms ground state trend initial data；The maximum operation branch road of active power is found out, disconnects its power flow transfer of the branch road to other branch roads；The most easily out-of-limit branch road found out in other branch roads；Calculating is sent by sensitivity of each unit to most easily out-of-limit branch road in end regions, is calculated adjustment amount and is reversely adjusted to branch road l after most easily out-of-limit branch road overload, overloadeNonoverload；Calculate now all Branch Power Flow summations, as maximum transfer power.All controllable units are ranked up by the present invention according to level of sensitivity, select sensitivity maximum during each pairing adjustment and minimum unit is matched, as a result more rationally reliable independent of balancing machine during adjustment.
Description
Technical field
The invention belongs to Power System and its Automation technical field, and in particular to one kind is adjusted based on equal and opposite quantities in pairs
The transmission cross-section limit method for solving of method.
Background technology
Modern interconnected network is in large scale and complicated, and strong and weak different, the conveying energy of part circuit is contacted between net
Power is inevitably restricted by the hot steady limit or the dynamic steady limit, is calculated for convenience of analysis and operational management, operations staff are led to
The contact section often power network being divided between some subsystems and system.The transmission of electricity limit of each contact section is grasped exactly, can
On the premise of safe operation of power system is ensured, farthest to meet the workload demand of all subregion.
The analysis of power transmitting capability and calculate all be all the time power system research important subject, existing meter
Calculation method is broadly divided into two classes：Deterministic Methods and probabilistic approach.Deterministic Methods mainly have direct solving method, repeated power flow
Method, Continuation Method, optimal load flow method, wherein, direct solving method causes circuit overload directly to increase interregional power,
Cut down survey region by half again and exchange power, its algorithm is simple but precision is not high, it is difficult to meets practical application request, rear several sides
Although method improves accuracy, but needed because solution procedure includes equation parameter, nonlinear equation processing and calculating process
Iteration etc. is wanted, causes calculating speed excessively slow, it is actual engineering can not to be applied to；Probabilistic approach mainly has stochastic programming, enumerated
Method and Monte Carlo Analogue Method etc., although these methods effectively meter and various uncertain factors, because it is dependent on a large amount of
Stochastic production sample, the calculating of large-scale power grid can not be adapted to.
The content of the invention
To solve the above problems, the invention discloses a kind of transmission cross-section Solving Limit based on equal and opposite quantities in pairs adjusting method
Solution method, is concretely comprised the following steps：
Step 1 takes a kind of typical load pattern under normal operating mode, forms ground state trend initial data：J-th strip branch
Road ljActive power be Pj, ljMaximum allowable power beJ=1,2 ...；
Step 2 finds out the maximum operation branch road of active power, is set to lk, disconnected branches lkThen equivalent to generation most serious event
Barrier, its power flow transfer is to other branch roads, in lkBranch road l in the state of disconnectionjActive power beNow all branch roads are not
Overload；
The most easily out-of-limit branch road that step 3 is found out in other branch roads, i.e. disconnected branches lkOther front and rear each branch roads have
The difference of work(powerMinimum branch road, is set to branch road le, then branch road leRequired adjustment amount is Δ Pe；
Step 4 is calculated and sent by each unit in end regions to branch road leSensitivity, and by sensitivity descending arrangement form join
Adjustment form is examined, sensitivity identical is matched somebody with somebody the first unit and last position unit according to the variable capacity descending sort of unit
Pair and adjustment, ith pairing adjustment amount be Δ Pi,
In formula, Se1For the sensitivity of the first unit, SemFor the sensitivity of last position unit,For the maximum of the first unit
Output limit value,For the EIAJ limit value of last position unit；
Step 5 judges present branch leAdjustment amountWhether Δ P is exceedede, wherein n is that present branch has matched
The logarithm of unit, if the branch road l not less than ifeDo not overload, the first unit reaches with the power that there must be a unit in last position unit
Limit value, it is removed from sequencing table, and continue again to be matched new the first unit and last position unit, return to step 4；
If the branch road l more than ifeIt has been overloaded that, present branch leAdjustment amount Δ P will reduce Δ PiTo branch road leNonoverload；
Step 6 calculates now all Branch Power Flow summations, as maximum transfer power.
Beneficial effect：
The present invention proposes that a kind of improved equal and opposite quantities in pairs adjusting method applies to solve maximum transfer power.According to sensitive
All controllable units are ranked up by degree size, and each pairing selects sensitivity maximum when adjusting and minimum unit is matched,
It is as a result more rationally reliable independent of balancing machine during adjustment.
Brief description of the drawings
Fig. 1 is the flow chart of the method for the invention；
Fig. 2 is the reference adjustment sequence list in the present invention；
Fig. 3 is IEEE39 nodes modular system wiring diagram in the embodiment of the present invention；
In figure：M is all units sum.
Embodiment
With reference to Figure of description and one group of specific embodiment, the present invention will be further described.
Fig. 1 is the flow chart of the method for the invention.The present invention's concretely comprises the following steps：
Step 1 takes a kind of typical load pattern under normal operating mode, forms ground state trend initial data：J-th strip branch
Road ljActive power be Pj, ljMaximum allowable power beJ=1,2 ...；
Step 2 finds out the maximum operation branch road of active power, is set to lk, disconnected branches lkThen equivalent to generation most serious event
Barrier, its power flow transfer is to other branch roads, in lkBranch road l in the state of disconnectionjActive power beNow all branch roads are not
Overload；
The most easily out-of-limit branch road that step 3 is found out in other branch roads, i.e. disconnected branches lkOther front and rear each branch roads have
The difference of work(powerMinimum branch road, is set to branch road le, then branch road leRequired adjustment amount is Δ Pe；
Step 4 is calculated and sent by each unit in end regions to branch road leSensitivity, and by sensitivity descending arrangement form join
Adjustment form is examined, sensitivity identical is matched somebody with somebody the first unit and last position unit according to the variable capacity descending sort of unit
Pair and adjustment, ith pairing adjustment amount be Δ Pi,
In formula, Se1For the sensitivity of the first unit, SemFor the sensitivity of last position unit,For the maximum of the first unit
Output limit value,For the EIAJ limit value of last position unit；
Step 5 judges present branch leAdjustment amountWhether Δ P is exceedede, wherein n is that present branch has matched
The logarithm of unit, if the branch road l not less than ifeDo not overload, the first unit reaches with the power that there must be a unit in last position unit
Limit value, it is removed from sequencing table, and continue again to be matched new the first unit and last position unit, return to step 4；
If the branch road l more than ifeIt has been overloaded that, present branch leAdjustment amount Δ P will reduce Δ PiTo branch road leNonoverload；
Step 6 calculates now all Branch Power Flow summations, as maximum transfer power.
Equal and opposite quantities in pairs method is a kind of method based on sensitivity, and its implication is：When adjusting system active power, it is
The node that the node that each increase is contributed finds a reduction output is paired, and vice versa；The node of pairing is active
The plus-minus amount of power is equal.In this way, the unbalanced power brought during conventional sensitivity class method adjustment unit output can be avoided
All undertaken by balancing machine, cause to balance the out-of-limit situation of acc power.Meanwhile in engineering in practice in the absence of the general of balancing machine
Read, be only concerned the plus-minus output situation of specific unit during adjustment during scheduling, and system power balances after ensureing adjustment, therefore reversely
Equivalent matched pair technique meets the scheduling rule of engineering in practice, and once adjusts a pair of nodes, can improve the efficiency of algorithm.
All controllable units are ranked up according to level of sensitivity, each pairing selects sensitivity maximum when adjusting and minimum
Unit matched, it is as a result more rationally reliable independent of balancing machine during adjustment.Detailed process is as follows：
If the active power that branch road l need to be adjusted is Δ Pl, active controllable unit number is m, according to its spirit to branch road road l
Sensitivity carries out descending order arrangement, forms reference adjustment sequence list as shown in Figure 2, and the big one end of sensitivity is head end,
The other end is end, is ranked up when sensitivity is identical according to the variable capacity of unit.
It is one group that the first unit and last position unit are selected during pairing.The first unit sensitivity is Sl1, increase output, last position machine
Group sensitivity is Slm, reduce and contribute, two units plus-minus power generating value is equal, then the first unit is to branch road l sensibility in practice：
S′l1=Sl1-Slm (1)
S′l1＞ 0, above formula can be seen that S 'l1Size and balancing machine selection it is not related.
The adjustment amount of unit is：
And in real system, the minimax output limit value of unit is once considered as during pairing adjustment process adjusting, therefore,
Actual adjustment amount is the minimum value in following three：
Obtain Δ PiAfterwards, the active power that need to be adjusted is corrected, if once pairing adjustment fails to meet to require, the first unit
Or must there is a power of the assembling unit to reach limit value in last position unit, it is removed from sequencing table, and continues new the first machine
Group with last position unit match and adjust again, until branch road l adjustment amount meet demand.
Using IEEE39 nodes modular system as research example, model and computational methods above are verified and analyzed.
The system includes 39 nodes, 10 generators, 12 transformers, and 34 circuits can be divided into 3 regions, and balance nodes are
Bus31.Reference voltage is 100kV, reference power 100MVA.System wiring figure and region division such as Fig. 3.
Region 1 is that region 3 is feeding section, region 1 and area by electric region it can be seen from ground state trend operation result
The discontinuity surface of domain 3 includes 3 transmission lines of electricity, respectively Bus16-19, Bus21-22, Bus16-24, zoning 1 to region 3
Maximum transfer power.Power flow simulation adjustment is carried out using BPA.
The branch road Bus21-22 that trend is maximum under the basic method of operation is disconnected, its active power is transferred to other two branch
Road.Ground state trend and section tidal current after disconnected branches Bus21-22 are as shown in table 1 below：
Table 1
As can be seen from Table 1, branch road Bus16-19 is most easily out-of-limit branch road, Δ P16-19=45.4MW.Calculate adjustable machine
Sensitivity of the group to Bus16-19 and by descending arrangement form with reference to adjustment table, it is as shown in table 2 below.
Table 2
The maximum B of Unit us34 of sensitivity it can be seen from upper table 2, the minimum unit of sensitivity is Bus38, therefore the
The unit for once carrying out reversely pairing adjustment is Bus34 and Bus38.Increase node Bus34 output, its variable capacity is
102MW, node Bus38 output being reduced, its variable capacity is 830MW, and the adjustment amount calculated according to formula (7) is 48.5561MW,
It is 48.5561MW that formula (8), which obtains final adjustment amount,.After adjusting for the first time, branch road Bus16-19 active is 501.7MW, out-of-limit.Machine
Group equal and opposite quantities in pairs adjustment result is as shown in table 3 below.
Table 3
According to above-mentioned calculating process, Δ P=2MW is set, continues equal and opposite quantities in pairs adjustment, but purpose is so that branch
Road Bus16-19 is not again out-of-limit, and now the active sum of branch road is maximum transfer power.The active tide of branch road after second of adjustment
Stream is as shown in below table 4.
Table 4
As can be seen from Table 4, the maximum transfer power in region 1 to region 3 is 1189.3MW.Regional internet has become electricity
The inexorable trend of Force system development, turns into that to influence bulk power grid development important between Power System Interconnection rear region the problem of electrical grid transmission ability
Factor.Can be rapidly and accurately to maximum transfer power under N-1 Static Security Constraints with equal and opposite quantities in pairs adjusting method is improved
The problem of solved, the fast and accurately feature of this method can apply to engineering in practice.
Claims (1)
1. a kind of transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method, it is characterised in that concretely comprise the following steps：
Step 1 takes a kind of typical load pattern under normal operating mode, forms ground state trend initial data：J-th strip branch road lj
Active power be Pj, ljMaximum allowable power beJ=1,2 ...；
Step 2 finds out the maximum operation branch road of active power, is set to lk, disconnected branches lkThen equivalent to occur most serious failure, its
Power flow transfer is to other branch roads, in lkBranch road l in the state of disconnectionjActive power beThe now equal nonoverload of all branch roads；
The most easily out-of-limit branch road that step 3 is found out in other branch roads, i.e. disconnected branches lkThe active power of other front and rear each branch roads
DifferenceMinimum branch road, is set to branch road le, then branch road leRequired adjustment amount is Δ Pe；
Step 4 is calculated and sent by each unit in end regions to branch road leSensitivity, and by sensitivity descending arrangement form with reference to adjustment
The first unit and last position unit are matched and adjusted according to the variable capacity descending sort of unit by form, sensitivity identical
Whole, the adjustment amount of ith pairing is Δ Pi,
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>min</mi>
<mo>{</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>e</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>e</mi>
<mi>m</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>&Delta;P</mi>
<mn>1</mn>
<mi>max</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&Delta;P</mi>
<mi>m</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>}</mo>
<mo>,</mo>
</mrow>
In formula, Se1For the sensitivity of the first unit, SemFor the sensitivity of last position unit, Δ P1 maxFor the EIAJ of the first unit
Limit value,For the EIAJ limit value of last position unit；
Step 5 judges present branch leAdjustment amountWhether Δ P is exceedede, wherein n is that present branch has matched unit
Logarithm, if the branch road l not less than ifeNot overloading, the first unit reaches limit value with the power that must have a unit in last position unit,
It is removed from sequencing table, and continues again to be matched new the first unit and last position unit, return to step 4；It is if super
Cross then branch road leIt has been overloaded that, present branch leAdjustment amount Δ P will reduce Δ PiTo branch road leNonoverload so that branch road leTrend weight
It is new not out-of-limit；
Step 6 calculates now all Branch Power Flow summations, as maximum transfer power.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510856649.8A CN105305504B (en) | 2015-11-30 | 2015-11-30 | Transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510856649.8A CN105305504B (en) | 2015-11-30 | 2015-11-30 | Transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105305504A CN105305504A (en) | 2016-02-03 |
| CN105305504B true CN105305504B (en) | 2018-01-23 |
Family
ID=55202431
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510856649.8A Expired - Fee Related CN105305504B (en) | 2015-11-30 | 2015-11-30 | Transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105305504B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107294103B (en) * | 2017-07-24 | 2020-01-14 | 广东工业大学 | Section flow control method and device |
-
2015
- 2015-11-30 CN CN201510856649.8A patent/CN105305504B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN105305504A (en) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103311926B (en) | Power system cascading failure simulation method based on unified power flow controller | |
| CN106655177B (en) | Distributed generation resource maximum access capability calculation method based on extension Second-order cone programming | |
| CN108023364B (en) | Power distribution network distributed generation resource maximum access capability calculation method based on convex difference planning | |
| CN101882167A (en) | Wind power station equivalent modeling method of large-scale wind power concentration access power grid | |
| CN107221945A (en) | A kind of UHVDC Transmission Lines forecast failure aid decision-making method and device | |
| CN102403716A (en) | Dynamic equalizing method for multi-infeed alternating/direct-current power grid | |
| Tanti et al. | Optimal placement of custom power devices in power system network to mitigate voltage sag under faults | |
| Beigvand et al. | Voltage stability analysis in radial smart distribution grids | |
| Xia et al. | Enhanced particle swarm optimisation applied for transient angle and voltage constrained discrete optimal power flow with flexible AC transmission system | |
| Dixit et al. | Optimal placement of SVC for minimizing power loss and improving voltage profile using GA | |
| CN107425527A (en) | A kind of THE UPFC static security prevention and control method | |
| Vural | Interior point‐based slack‐bus free‐power flow solution for balanced islanded microgrids | |
| CN105305504B (en) | Transmission cross-section limit method for solving based on equal and opposite quantities in pairs adjusting method | |
| CN103280797A (en) | Day-ahead static security correction method | |
| Zhang et al. | Impact of dynamic load models on transient stability-constrained optimal power flow | |
| Venkatesh et al. | Investigations on hybrid line stability ranking index with polynomial load modeling for power system security | |
| Hussein et al. | Employment of PSO algorithm to improve the neural network technique for radial distribution system state estimation | |
| Avudayappan et al. | Congestion management in deregulated power system using hybrid cat-firefly algorithm with TCSC and SVC FACTS devices | |
| CN115940157A (en) | Method, device and equipment for automatically generating load flow scene of stability control strategy checking task | |
| Yang et al. | Optimal usage of transmission capacity with FACTS devices in the presence of wind generation: A two-stage approach | |
| Adebayo et al. | Static voltage stability enhancement using FACTS controller | |
| CN102420426B (en) | Power grid partial loss reduction method based on injection power sensitivity | |
| Muthulakshmi et al. | The proper location and sizing of multiple distributed generators for maximizing voltage stability using PSO | |
| CN111211581B (en) | Method and device for determining new energy acceptance capacity in power system | |
| CN103280841A (en) | On-grid point confirmation method for intermittent energy sources |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180123 Termination date: 20201130 |