CN104701864A - Reactive power planning method using SVC as reactive power compensation equipment - Google Patents
Reactive power planning method using SVC as reactive power compensation equipment Download PDFInfo
- Publication number
- CN104701864A CN104701864A CN201510097259.7A CN201510097259A CN104701864A CN 104701864 A CN104701864 A CN 104701864A CN 201510097259 A CN201510097259 A CN 201510097259A CN 104701864 A CN104701864 A CN 104701864A
- Authority
- CN
- China
- Prior art keywords
- svc
- reactive power
- voltage
- node
- reactive
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
- H02J3/1821—Arrangements for adjusting, eliminating or compensating reactive power in networks using shunt compensators
- H02J3/1871—Methods for planning installation of shunt reactive power compensators
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/10—Flexible AC transmission systems [FACTS]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a reactive power planning method for an electrical power system. The reactive power planning method uses SVC and considers the transient process stability constraint of the system after large power system disturbance. The reactive power planning method using the SVC as reactive power compensation equipment includes that generating a reactive power compensation equipment selected-point index matrix through different voltage stabilization principles, and using fuzzy clustering to perform reactive power compensation device selected-point arrangement to obtain a reactive power compensation point distribution scheme suitable for the SVC running performance; through building a reactive power compensation optimization model based on different running modes and a dynamic model based on complex high-order system equipment, using a time domain simulation and optimization algorithm to obtain the optimal result of the SVC capacity configuration of the corresponding reactive power compensation point so as to realize SVC reactive power planning. According to a fuzzy clustering theory, a dynamic electrical power system theory, a time domain simulation technology and an interior point method optimization synthesis algorithm, the reactive power planning method realizes the reactive power planning which considers the transient process stability constraint after the system goes wrong, fully develops the dynamic response characteristic of the SVC and enables the reactive power compensation running efficiency of the electric power system to be effectively improved.
Description
Technical field
The present invention relates to a kind of SVC that adopts as the idle planing method of reactive-load compensation equipment, in order to improve power system planning and to run the stability of a system controlled.Belong to power system planning and operation field.
Background technology
Idle planning is by analyzing the idle characteristic under Power System Steady-state and dynamic operating condition, control the distribution of reactive power in the whole network and flowing, distribute installation site and the compensation capacity of reactive-load compensation equipment rationally, optimize voltage-regulation means, to improve system run all right, improve interregional tie-line power transmission, meet the safe operation of electrical network in Steady state and transient state situation and the demand of control voltage fluctuation.
Traditional idle planning acquiescence adopts switching type Shunt Capacitor Unit as reactive power compensator, and compensation equipment is a fundamental frequency equivalence capacitive reactance model.Idle planning application research contents comprises: the reconnaissance planning of (1) newly-increased reactive-load compensation equipment; (2) capacity planning of newly-increased reactive-load compensation equipment.
Along with being on the increase of power system reactive power compensation equipment, except original Shunt Capacitor Unit, the FACTS type reactive-load compensation equipment that SVC, STATCOM etc. have fast dynamic response characteristic is widely used gradually.
In the research of power system reactive power planning problem, for the actual utility of power system reactive power planning, the input of reactive-load compensation equipment is the voltage stability in order to improve system.For this reason, consider the idle planning problem of the transient process response of the dynamic devices such as generator, excitation system, SVC, in Transient Stability Constraints condition except equation of rotor motion and angle stability, also should consider the Voltage Stability Constraints in transient process.Except the constraints under stable situation, the constraints in transient process also comprises: (1) generator amature angle is relative to the range constraint at inertia center; (2) system each node transient process Voltage Stability Constraints; (3) the operation equation constraint of dynamic device.
Summary of the invention
Technical problem: for solving the problems of the technologies described above, the invention provides a kind of employing SVC simultaneously considering that steady operation and the transient process stability of a system retrain and carries out the idle planing method of power transmission network.
Technical scheme: a kind of SVC that adopts of the present invention as the idle planing method of reactive-load compensation equipment,
Concrete steps are as follows:
Step 1: based on system failure ranking method, calculates the transient process voltage stability index after most catastrophe failure respectively according to formula (1), (2), (3) under peak load, average load level:
System transient modelling Voltage Drop:
Wherein, V
sfor short-term stability process interior joint voltage lowest amplitude, V
0for node voltage initial magnitude;
Variation index in transient process:
Low-voltage duration index in transient process:
Wherein, wherein, t
clfor the fault clearance time, t
ffor the emulation termination time;
Step 2: each voltage stability index in the transient process under the peak load obtained according to step 1, average load level, adopt fuzzy clustering method, obtain respectively SVC under system peak load, average load level the node i location sets N of candidate's reactive power compensation
c;
Step 3: set up SVC capacity planning optimal model:
Select three rank or more generator model, single order excitation system model and mixing load model:
Three rank generator models:
Wherein, δ (t) is for q axle is relative to the angular displacement of synchronous rotating frame; ω
brotor rated speed and rotor t angular speed is respectively with ω (t); T
jwith T '
d0for time constant, T
mwith T
et () is respectively mechanical force moment and dynamo-electric moment; E
f, E '
qbe respectively stator excitation electromotive force and motor q axle transient internal voltage; X
dwith X '
dbe respectively d axle synchronous reactance and d axle transient reactance; I
dfor d shaft current;
Its set end voltage equation is:
Wherein, V
d, V
q(I
d, I
q) be respectively d shaft voltage and q shaft voltage (electric current); x
q, r
abe respectively q axle synchronous reactance and winding resistance;
Single order excitation system model:
Wherein, T
afor time constant; K
afor control coefrficient; V
refand V
mbe respectively node voltage reference value and amplitude; E
f0for excitation electric gesture initial value;
Target function is made for condition so that the idle input cost of SVC is minimum:
Wherein, C
fix_ifor idle investment fixed cost; C
i_v1for capacitive reactive power variable cost; C
i_v2for the idle variable cost of perception;
for node i capacitive SVC compensation capacity;
for the SVC compensation capacity of perception;
Static Security Constraints condition is
In constraints, P
gi, P
lithe generated power at node i place is exerted oneself and load is gained merit demand respectively; Q
gi, Q
lithe generator reactive at node i place is exerted oneself and reactive load demand respectively; Q
cifor the reactive compensation capacity of node i; N
bfor node set, N
gfor system generator set, N
cfor SVC compensates node set; G
ijwith B
ijbe respectively conductance and the susceptance of node admittance matrix; θ
ijfor the phase angle difference of node i and node j; V
i, V
imax, V
iminbe respectively node i voltage magnitude, voltage magnitude maximum and minimum value; P
imax, P
iminbe respectively generator node i meritorious output maximum and minimum value; Q
imax, Q
iminbe respectively the idle output maximum of generator node i and minimum value; Q
sVCi, Q
sVCimaxand Q
sVCiminbe respectively reactive power compensation both candidate nodes i reactive compensation capacity, reactive power compensation maximum and minimum value; S
ij, S
ijmaxbe respectively apparent power and the maximum thereof of circuit ij; y
ifor Boolean variable, y
i=1 represents that node i is candidate's reactive power compensation point;
In transient process, Voltage Stability Constraints comprises:
Short-term Voltage falls constraint:
Low-voltage duration constraints in Short-term Voltage Stability process: load bus voltage is at 0.75V
0≤ V
s(t)≤0.8V
0the duration of scope can not more than 20 cycles, i.e. 0.4s;
Reach voltage magnitude constraint: V after new steady-state process
i>=0.95V
0ii=1 ..., N
b;
Step 4: adopt hiding-trapezium integral method, by differential equation differencing method, changes algebraic step differentiation equation into by the differential equation in Optimized model; According to time-domain-simulation method and interior some optimization, according to target function and constraints, calculate the reactive compensation capacity of the idle planned position of each candidate of SVC under system peak load, average load operational mode respectively;
Step 5: the SVC reconnaissance factor obtained under different load level according to step 2,4, in conjunction with time-domain-simulation analysis, checks position and capacity that SVC finally needs compensation.
Fuzzy clustering algorithm in described step 2:
Based on each load bus index of transient process after catastrophe failure under two kinds of load levels that step 1 calculates, carry out data normalization;
Calculate the similarity r between each node i and j
ij, form fuzzy similarity matrix;
By the calculating of transitive closure, obtain fuzzy equivalent matrix;
By the element in fuzzy equivalent matrix, adopt fuzzy nonhierarchial clustering to be classified by all PQ nodes, accept horizontal threshold values λ by dynamic conditioning, λ is increased to gradually close to 1 from a less threshold values, obtains dynamic clustering result; Obtain some group nodes weak in all nodes;
Carry out F inspection, the newly-increased SVC candidate infield under selected catastrophe failure.
If the newly-increased reactive power compensation point position under two kinds of operating conditions is identical, then this newly-increased compensated position selected, and the reactive capability optimum results under complex optimum two kinds of operation levels.
If the newly-increased reactive power compensation point position under two kinds of operating conditions is different, then respectively on the basis of existing optimum results, optimize the planned capacity of newly-increased SVC after catastrophe failure under another operating condition.
When adopting TSC-TCR type SVC as planning apparatus, the principle that should be the 1/N of SVC total capacity according to the Capacity Selection of TCR carries out the configuration of TCR capacity.
Beneficial effect: compared with prior art, the inventive method has the following advantages:
The present invention is based on fuzzy clustering algorithm, time-domain-simulation analysis and interior point method optimum theory, adopt SVC alternatively reactive-load compensation equipment, consider its dynamic response characteristic, simultaneously, based on the dynamic action of three rank generator models and single order excitation system, consider above-mentioned constraints, invented under meeting the different operational mode of system, consider the idle planing method of steady operation and the constraint of the transient process stability of a system simultaneously.Carry out enforcement checking based on IEEE test macro, after considering fault under utilizing the quick features of response of SVC to achieve different running method, the power system reactive power of transient process scleronomic constraint is planned.Embodiment analysis shows, effective configuration SVC reactive power compensation, while guarantee economy, can realize the idle planning meeting different scleronomic constraint requirement, have theoretical and engineer applied value.
Accompanying drawing explanation
Fig. 1 is the idle planning flow chart of newly-increased SVC;
Fig. 2 is Dynamic Fuzzy Clustering Algorithm process schematic;
Fig. 3 is the membership function schematic diagram of variation index VDI;
Fig. 4 is transient process generator's power and angle situation of change schematic diagram under most catastrophe failure under peak load;
Fig. 5 is transient process load bus voltage magnitude situation of change schematic diagram under most catastrophe failure under peak load;
Fig. 6 is transient process system voltage recovery curve figure after compensating containing SVC under peak load;
Fig. 7 is the change curve at relative inertness center, generator amature angle in transient process after peak load level compensates containing SVC;
Fig. 8 is without SVC bucking-out system transient process load bus voltage change curve figure under average load level.
Embodiment
Below in conjunction with embodiment and Figure of description; the execution mode of technical solution of the present invention and concrete operating process are elaborated; but embodiment is only the preferred embodiment of the present invention; be noted that for those skilled in the art; under the premise without departing from the principles of the invention; can also make some improvement and equivalent replacement, these improve the claims in the present invention and are equal to the technical scheme after replacing, and all fall into protection scope of the present invention.
Fig. 1 is the idle planning flow chart of newly-increased SVC.As shown in Figure 1, the present invention relates to the idle planing method adopting SVC as reactive-load compensation equipment, concrete steps are as follows:
Step 1: based on system failure ranking method, calculates the transient process voltage stability index after most catastrophe failure respectively according to formula (1), (2), (3) under peak load, average load level:
System transient modelling Voltage Drop:
Wherein, V
sfor short-term stability process interior joint voltage lowest amplitude, V
0for node voltage initial magnitude;
Variation index in transient process:
Low-voltage duration index in transient process:
Wherein, wherein, t
clfor the fault clearance time, t
ffor the emulation termination time, the starting point divided by T is that this index does not have dimension, is a normalized index.
Step 2: each voltage stability index in the transient process under the peak load obtained according to step 1, average load level, adopt fuzzy clustering method, obtain respectively SVC under system peak load, average load level the node i location sets N of candidate's reactive power compensation
c;
Step 3: set up SVC capacity planning optimal model:
Select three rank or more generator model, single order excitation system model and mixing load model:
The generator amature equation of motion adopts third-order model:
Its voltage equation is
In formula, δ is the angular displacement of rotor q axle relative to synchronous rotating frame real axis x; Time constant M, T '
d0, D, T
afor famous value (unit s); Its dependent variable is per unit value.V
d, V
q, I
d, I
qfor the voltage and current phasor under d-q coordinate,
Single order excitation model:
V
mfor voltage magnitude.V
reffor reference voltage amplitude.T
e=E′
qI
q-(X′
d-X
q)I
dI
q。
According to the principle that state variable is not suddenlyd change, generator rotor angle and E ' can be obtained
q, E
f0and T
minitial value:
Transient process stable constraint after consideration system generation catastrophe failure, the target function of concrete optimization problem and constraints are:
(1), target function
Adopt SVC as reactive-load compensation equipment, consider that the newly-increased reactive power compensation object of planning function of transient process scleronomic constraint is chosen for idle investment cost minimum:
(2), Static Security Constraints condition
Static Security Constraints comprises system load flow equality constraint, voltage magnitude constraint, generated power, the constraint of idle bound, line transmission capacity-constrained, reactive-load compensation equipment capacity bound and reactive power compensation point restricted number.
Wherein, C
fix_ifor idle investment fixed cost; C
i_v1for capacitive reactive power variable cost; C
i_v2for the idle variable cost of perception;
for node i capacitive SVC compensation capacity;
for the SVC compensation capacity of perception; In constraints, P
gi, P
lithe generated power at node i place is exerted oneself and load is gained merit demand respectively; Q
gi, Q
lithe generator reactive at node i place is exerted oneself and reactive load demand respectively; Q
cifor the reactive compensation capacity of node i; N
bfor node set, N
gfor system generator set, N
cfor SVC compensates node set; G
ijwith B
ijbe respectively conductance and the susceptance of node admittance matrix; θ
ijfor the phase angle difference of node i and node j; V
i, V
imax, V
iminbe respectively node i voltage magnitude, voltage magnitude maximum and minimum value; P
imax, P
iminbe respectively generator node i meritorious output maximum and minimum value; Q
imax, Q
iminbe respectively the idle output maximum of generator node i and minimum value; Q
sVCi, Q
sVCimaxand Q
sVCiminbe respectively reactive power compensation both candidate nodes i reactive compensation capacity, reactive power compensation maximum and minimum value; S
ij, S
ijmaxbe respectively apparent power and the maximum thereof of circuit ij; y
ifor Boolean variable, y
i=1 represents that node i is candidate's reactive power compensation point;
(3), transient process scleronomic constraint
A. the center of inertia COI of selecting system is reference, and transient process angle stability constraint representation is:
Wherein,
set δ herein
maxbe 2, unit is rad.
On the other hand, for transient process Voltage-stabilizing Problems, need to consider to occur in the transient process of N-1 fault in system, the operational factor of system voltage within the limits prescribed.According to voltage time domain running orbit, NERC/WECC evaluation criterion, adopts following index:
B. transient process Voltage Drop constraint
First, according to the definition of system transient modelling Voltage Drop:
wherein, V
sfor transient process interior joint voltage lowest amplitude, V
0for node voltage initial magnitude.
Require load bus V
dip≤ 25%, non-load bus V
dip≤ 30%; Or load bus meets V within continuous 20 cycles
dip≤ 20%.
For load bus,
i.e. V
s(t)≤0.75V
0.Meanwhile, for preventing generator overexcitation, need to set voltage magnitude upper limit constraint V simultaneously
max.Therefore, transient process Voltage Drop constraint representation is:
C. transient process low-voltage duration constraints
Load bus voltage is at 0.75V
0≤ V
s(t)≤0.8V
0the duration of scope can not more than 20 cycles, i.e. 0.4s.
D. voltage magnitude constraint after new steady-state process is reached
V
i≥0.95V
0ii=1,...,N
b
From the Optimized model set up above, consider that the power system reactive power planning problem of multiple reactive-load compensation equipment is a Nonlinear Mixed Integer Programming Problem containing differential-algebraic equation.Dimension is solved, by PROBLEM DECOMPOSITION for newly-increased SVC reconnaissance and capacity optimize two processes for reducing.Meanwhile, consider different operational modes and catastrophe failure, the optimal result of reactive power compensation is determined in comprehensive analysis.Contingency screening and optimizing process adopt based on hiding-trapezium integral method time-domain-simulation, and integration step is 0.01s.
Step 4: adopt hiding-trapezium integral method, by differential equation differencing method, changes algebraic step differentiation equation into by the differential equation in Optimized model; According to time-domain-simulation method and interior some optimization, according to target function and constraints, calculate the reactive compensation capacity of each idle planned position of SVC under different load operational mode respectively;
Step 5: the SVC reconnaissance factor obtained under different load level according to step 2,4, in conjunction with time-domain-simulation analysis, checks position and capacity that SVC finally needs compensation.
1.2 times that system based on IEEE3 machine 9 node system, and is gained merit by the specific embodiment of the invention, load or burden without work rises to initial data.Correspondingly, the meritorious of generator is exerted oneself also according to the proportional increase of its original situation of exerting oneself 1.2 times.The average load level of system as system to be planned will be obtained.For this basic load operation operating mode, setting peak load be when load gain merit, reactive requirement level brings up to 1.2 times time situation.
In the specific embodiment of the invention, due under different load levels, the situation that most catastrophe failure occurs generally can be different, and the newly-increased compensation point position of SVC thus will be caused different, and added compensation capacity is different.Therefore, for distributing SVC compensation point position rationally, the present embodiment is based on two kinds of different operating conditions: under peak load level and average load level (average load level), carry out reconnaissance analysis respectively, consider the idle planning of SVC of transient process Voltage Stability Constraints to SVC.
In step 1, each transient process voltage stability index after most catastrophe failure is all based on phasor form, and each index describes the voltage stabilization situation comprising each load bus of power system network.
Before idle work optimization, that selects suitable reactive power compensation position effectively can reduce whole planning problem solves difficulty, and Reduction Computation time greatly.Consider SVC to be installed on load bus, therefore the present invention is only concerned about and sets up the relevant dynamic indicator of load bus.Based on domain simulation eur, the bus voltage stability index of step 1 is set up as follows:
The first: variation index in transient process;
Define according to system transient modelling Voltage Drop:
Wherein, V
sfor transient process node voltage lowest amplitude, V
0for node voltage initial magnitude.
According to V
dipdefinition, be eliminate step effect, adopt the S type function of expansion, definition variation index VDI (Voltage Deviation Index):
The second: transient process low-voltage duration index;
After calculating fault clearance, [t
cl, t
f] period T=t
f-t
clinterior each load bus voltage is lower than 0.8V
0accumulative total time, definition low-voltage duration index LVDI (Low Voltage Duration Index):
Wherein, t
clfor the fault clearance time, t
ffor the emulation termination time.
Adopt the set that obtains of fuzzy clustering algorithm respectively based on the disturbance under system different load level run pattern in step 2, belong to the result that many cover operational systems calculate respectively.Its computational process is:
Step 21: data normalization: because the Principles and ways of each Index Establishment is different, initial data aims of standardization eliminate the difference between different indexs in dimension and dimension, under making each index all be distributed in unified scale.Usually, data normalization is interval to [0,1].
Step 22: form fuzzy similarity matrix: based on the matrix after the normalization obtained above, can similarity r between computing node i and j
ij, form fuzzy similarity matrix.
Step 23: calculate fuzzy equivalent matrix: for making matrix meet reflectivity, symmetry and transitivity simultaneously, further by the calculating of transitive closure, fuzzy equivalent matrix can be obtained.
Step 24: Dynamic Fuzzy Clustering Algorithm: according to fuzzy equivalent matrix, can carry out the Dynamic Fuzzy Clustering Algorithm analysis of system.For different threshold values, the λ Level Matrix of system can be obtained, draw final cluster result.Dynamic clustering process as shown in Figure 2.
Step 25:F checks: according to different threshold values, and calculate F value, larger F value shows to there is more obvious difference between each grouping.
For selecting the reasonability of suitable value assessment classification, adopt F inspection
In formula: n
jit is the node number of jth group in fuzzy clustering result; R is the group number of fuzzy cluster analysis
F is that to obey the degree of freedom be the statistic that the F of (r-1, n-r) distributes.According to different threshold values, calculate F value, larger F value shows may there is more obvious difference between each grouping.
By early stage shunt capacitor program results and SVC program results carry out integrated decision-making analysis.If the SVC reconnaissance obtained under peak load level and under average load level is consistent, then using optimum results higher for SVC value in the two optimum results as final SVC program results.If the SVC dot position obtained under two kinds of load levels is different, then the SVC program results under peak load is first adopted to check the transient process stability under average load level.If the transient process stability under average load level can be maintained, then the very final program results of the SVC program results obtained under peak load level; Otherwise, when the horizontal SVC program results of peak load drops into, by the SVC candidate compensation point determined under average load level progressively little step-length increase progressively the mode of SVC capacity, check through time-domain-simulation, obtain final SVC program results.
The aspect considered in complex optimum interpretation of result in steps of 5 mainly comprises:
A. due under different load levels, the situation that most catastrophe failure occurs generally can be different, and the newly-increased compensation point position of SVC thus will be caused different, and added compensation capacity is different.Therefore, consider that the idle planning of SVC of transient process Voltage Stability Constraints need consider the problem broken down under peak load level and average load level two kinds of different running method.
If the newly-increased reactive power compensation point position b. under two kinds of operating conditions is identical, then this newly-increased compensated position selected, and the reactive capability optimum results under complex optimum two kinds of operation levels, generally, select capability value that optimum results under two kinds of operation levels is larger as optimum capacity configuration result;
If the newly-increased reactive power compensation point position c. under two kinds of operating conditions is different, then respectively on the basis of existing optimum results, optimize the planned capacity of newly-increased SVC after catastrophe failure under another operating condition.
D. finally in conjunction with time-domain-simulation analysis, position and capacity that SVC finally needs compensation is checked.
E., when adopting TSC-TCR type SVC as planning apparatus, the principle that should be the 1/N of SVC total capacity according to the Capacity Selection of TCR carries out the configuration of TCR capacity.
In the specific embodiment of the invention, under two kinds of different operation levels, setting failure mode is that in system, a certain bar circuit, when 0.1s, three-phase ground short trouble occurs, and fault continues 0.1s (5 power frequency periods), line disconnection afterwards, excision fault.
(1) the idle planning under peak load
In peak load level run situation, most catastrophe failure is that circuit 8 (node 8-node 9), when node 8 is in 0.1s, three-phase ground short circuit occurs, and cut-offs circuit excision fault after 0.1s.When dropping into without SVC, time-domain-simulation obtains in transient process that generator's power and angle curve is as Fig. 4, and load bus voltage change curve is as Fig. 5.
After 0.2s moment fault clearance, system voltage sharply declines, and voltage stability index destroys, and when 0.6s, system exciting voltage reaches the highest, causes voltage keeps to vibrate, and finally there occurs voltage collapse event.Voltage lose stable after, there is merit angle unstability.For this reason, need to consider that planning drops into the reactive-load compensation equipment with fast dynamic response ability, stablize as SVC etc. maintains system transient modelling process voltage.
(i) SVC reconnaissance
Based on the present embodiment situation, calculate the index that each load bus is corresponding, as shown in table 1.
The each load bus index of transient process after table 1 two kinds of catastrophe failures
Fig. 2 is specific embodiment of the invention Dynamic Fuzzy Clustering Algorithm process schematic.As shown in Figure 2, under peak load level, the Data Source of fuzzy clustering is VDI and LVDI of table 1.
Fuzzy clustering algorithm is through data normalization: because the Principles and ways of each Index Establishment is different, initial data aims of standardization eliminate the difference between different indexs in dimension and dimension, under making each index all be distributed in unified scale; Usually, data normalization is interval to [0,1];
Obtain the matrix after based on the normalization obtained above, can similarity r between computing node i and j
ij, form fuzzy similarity matrix; By the calculating of transitive closure, obtain fuzzy equivalent matrix; Fuzzy equivalent matrix is:
By
element in matrix, adopts fuzzy nonhierarchial clustering to be classified by all 6 PQ nodes, accepts horizontal threshold values λ by dynamic conditioning, λ is increased to gradually close to 1 from a less threshold values, obtains dynamic clustering result.Node 5, node 6 are classified as two the weakest group nodes earlier respectively.
Finally in conjunction with F inspection, selected catastrophe failure lower node 5, node 6 are newly-increased SVC candidate infield.
(ii) the capacity optimization of SVC
Target function adopts newly-increased idle input to minimize, and using the SVC optimum results do not considered in transient process Voltage Stability Constraints situation as running initial value, under peak load level, node 5, node 6 increase SVC susceptance capacity newly and be 0.Calculation optimization model obtains considering that the SVC idle work optimization result of transient process Voltage Stability Constraints is:
Need to set up SVC capacity susceptance 0.7687p.u. in node 5 transient process, inductive susceptance 0p.u..
Newly-increased SVC capacity susceptance 0.4883p.u. is needed, inductive susceptance 0p.u. in node 6 transient process.
Based on time-domain-simulation, consider catastrophe failure form mentioned above, obtain compensating rear system voltage recovery curve as shown in Figure 6 containing SVC, containing relative inertness center, generator amature angle change curve in system transient modelling process after SVC compensation as shown in Figure 7.Emulation continues 10s.Easily see that system meets angle stability constraint, and reach a new poised state.The frequency stability of transient process also obtains checking checking.
(2) the idle planning under average load
Under average load level, most catastrophe failure is that circuit 8 (node 8-node 9), when node 8 is in 0.1s, three-phase ground short circuit occurs, and cut-off circuit excision fault after 0.1s, time-domain-simulation continues 3s.Fig. 8 is each load bus change in voltage situation in transient process.In this transient process, occurred the situation that NERC/WECC voltage indexes destroys, but last voltage resume is to a more stable value, system is stabilized to a new state again.System does not lose angle stability in transient process.In addition, when arranging after other each circuits of system occur three phase short circuit fault separately, all do not occur the situation that NERC/WECC voltage indexes destroys, system merit angle also remains stable.
(i) SVC reconnaissance
Under average load level, Data Source is VDI and the LVDI index of table 1.Adopt similar thinking, the node 5 that is easy to get is regarded as the most weak node of voltage always and is separated with other node clusters.In conjunction with F inspection, the newly-increased SVC candidate infield obtained under average load level is node 5.
(ii) the capacity optimization of SVC
Under average load level, consider that the SVC idle work optimization result of transient process Voltage Stability Constraints is: need to set up SVC capacity susceptance 0.5615p.u. in node 5 transient process, inductive susceptance 0p.u..
After increasing SVC compensation, after system generation catastrophe failure, transient process voltage stability index and transient stability index all meet rule.Omit herein and compensate the curve of rear system parameters in transient process, analytical method is identical with the check method under peak load mentioned above.After optimizing, under fault, transient process index all meets indication range, and the verification of transient process frequency stabilization is passed through.
(3) unified plan result
The optimum results obtained when comprehensive different operation level and catastrophe failure, the result of study of SVC configuration is comprehensively analyzed, can determine that idle program results is:
Node 5 need increase configuration capacity susceptance 0.7687p.u., inductive susceptance 0p.u. newly;
Node 6 need increase configuration capacity susceptance 0.4883p.u., inductive susceptance 0p.u. newly.
When adopting TSC-TCR type SVC as planning apparatus, the principle that should be the 1/N of SVC total capacity according to the Capacity Selection of TCR carries out the configuration of TCR capacity.
The above is only the preferred embodiment of the present invention; be noted that for those skilled in the art; under the premise without departing from the principles of the invention, can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.
Claims (5)
1. adopt SVC as an idle planing method for reactive-load compensation equipment, it is characterized in that: concrete steps are as follows:
Step 1: based on system failure ranking method, calculates the transient process voltage stability index after most catastrophe failure respectively according to formula (1), (2), (3) under peak load, average load level:
System transient modelling Voltage Drop:
Wherein, V
sfor short-term stability process interior joint voltage lowest amplitude, V
0for node voltage initial magnitude;
Variation index in transient process:
Low-voltage duration index in transient process:
Wherein, wherein, t
clfor the fault clearance time, t
ffor the emulation termination time;
Step 2: each voltage stability index in the transient process under the peak load obtained according to step 1, average load level, adopt fuzzy clustering method, obtain respectively SVC under system peak load, average load level the node i location sets N of candidate's reactive power compensation
c;
Step 3: set up SVC capacity planning optimal model:
Select three rank or more generator model, single order excitation system model and mixing load model:
Three rank generator models:
Wherein, δ (t) is for q axle is relative to the angular displacement of synchronous rotating frame; ω
brotor rated speed and rotor t angular speed is respectively with ω (t); T
jwith T '
d0for time constant, T
mwith T
et () is respectively mechanical force moment and dynamo-electric moment; E
f, E '
qbe respectively stator excitation electromotive force and motor q axle transient internal voltage; X
dwith X '
dbe respectively d axle synchronous reactance and d axle transient reactance; I
dfor d shaft current;
Its set end voltage equation is:
Wherein, V
d, V
q(I
d, I
q) be respectively d shaft voltage and q shaft voltage (electric current); x
q, r
abe respectively q axle synchronous reactance and winding resistance;
Single order excitation system model:
Wherein, T
afor time constant; K
afor control coefrficient; V
refand V
mbe respectively node voltage reference value and amplitude; E
f0for excitation electric gesture initial value;
Target function is made for condition so that the idle input cost of SVC is minimum:
Wherein, C
fix_ifor idle investment fixed cost; C
i_v1for capacitive reactive power variable cost; C
i_v2for the idle variable cost of perception;
for node i capacitive SVC compensation capacity;
for the SVC compensation capacity of perception;
Static Security Constraints condition is
In constraints, P
gi, P
lithe generated power at node i place is exerted oneself and load is gained merit demand respectively; Q
gi, Q
lithe generator reactive at node i place is exerted oneself and reactive load demand respectively; Q
cifor the reactive compensation capacity of node i; N
bfor node set, N
gfor system generator set, N
cfor SVC compensates node set; G
ijwith B
ijbe respectively conductance and the susceptance of node admittance matrix; θ
ijfor the phase angle difference of node i and node j; V
i, V
imax, V
iminbe respectively node i voltage magnitude, voltage magnitude maximum and minimum value; P
imax, P
iminbe respectively generator node i meritorious output maximum and minimum value; Q
imax, Q
iminbe respectively the idle output maximum of generator node i and minimum value; Q
sVCi, Q
sVCimaxand Q
sVCiminbe respectively reactive power compensation both candidate nodes i reactive compensation capacity, reactive power compensation maximum and minimum value; S
ij, S
ijmaxbe respectively apparent power and the maximum thereof of circuit ij; y
ifor Boolean variable, y
i=1 represents that node i is candidate's reactive power compensation point;
In transient process, Voltage Stability Constraints comprises:
Short-term Voltage falls constraint:
Low-voltage duration constraints in Short-term Voltage Stability process: load bus voltage is at 0.75V
0≤ V
s(t)≤0.8V
0the duration of scope can not more than 20 cycles, i.e. 0.4s;
Reach voltage magnitude constraint: V after new steady-state process
i>=0.95V
0ii=1 ..., N
b;
Step 4: adopt hiding-trapezium integral method, by differential equation differencing method, changes algebraic step differentiation equation into by the differential equation in Optimized model; According to time-domain-simulation method and interior some optimization, according to target function and constraints, calculate the reactive compensation capacity of the idle planned position of each candidate of SVC under system peak load, average load operational mode respectively;
Step 5: the SVC reconnaissance factor obtained under different load level according to step 2,4, in conjunction with time-domain-simulation analysis, checks position and capacity that SVC finally needs compensation.
2. employing SVC according to claim 1 is as the idle planing method of reactive-load compensation equipment, it is characterized in that: the fuzzy clustering algorithm in described step 2:
Based on each load bus index of transient process after catastrophe failure under two kinds of load levels that step 1 calculates, carry out data normalization;
Calculate the similarity r between each node i and j
ij, form fuzzy similarity matrix;
By the calculating of transitive closure, obtain fuzzy equivalent matrix;
By the element in fuzzy equivalent matrix, adopt fuzzy nonhierarchial clustering to be classified by all PQ nodes, accept horizontal threshold values λ by dynamic conditioning, λ is increased to gradually close to 1 from a less threshold values, obtains dynamic clustering result; Obtain some group nodes weak in all nodes;
Carry out F inspection, the newly-increased SVC candidate infield under selected catastrophe failure.
3. employing SVC according to claim 1 is as the idle planing method of reactive-load compensation equipment, it is characterized in that: if the newly-increased reactive power compensation point position under two kinds of operating conditions is identical, then this newly-increased compensated position selected, and the reactive capability optimum results under complex optimum two kinds of operation levels.
4. employing SVC according to claim 1 is as the idle planing method of reactive-load compensation equipment, it is characterized in that: if the newly-increased reactive power compensation point position under two kinds of operating conditions is different, then respectively on the basis of existing optimum results, optimize the planned capacity of newly-increased SVC after catastrophe failure under another operating condition.
5. employing SVC according to claim 1 is as the idle planing method of reactive-load compensation equipment, it is characterized in that: when adopting TSC-TCR type SVC as planning apparatus, the principle that should be the 1/N of SVC total capacity according to the Capacity Selection of TCR carries out the configuration of TCR capacity.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510097259.7A CN104701864A (en) | 2015-03-04 | 2015-03-04 | Reactive power planning method using SVC as reactive power compensation equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510097259.7A CN104701864A (en) | 2015-03-04 | 2015-03-04 | Reactive power planning method using SVC as reactive power compensation equipment |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104701864A true CN104701864A (en) | 2015-06-10 |
Family
ID=53348734
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510097259.7A Pending CN104701864A (en) | 2015-03-04 | 2015-03-04 | Reactive power planning method using SVC as reactive power compensation equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104701864A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105939020A (en) * | 2016-05-14 | 2016-09-14 | 国家电网公司 | Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability |
CN109596979A (en) * | 2018-12-06 | 2019-04-09 | 浙江大学 | A kind of synchronous generator system dynamic state estimator method based on interior optimization algorithm |
CN110707754A (en) * | 2019-08-28 | 2020-01-17 | 广东工业大学 | Optimization method for water, wind and light power supply capacity configuration in micro-grid |
CN111598718A (en) * | 2020-04-03 | 2020-08-28 | 广东电网有限责任公司电力调度控制中心 | STATCOM configuration planning method, device and equipment |
CN112072662A (en) * | 2020-08-28 | 2020-12-11 | 武汉大学 | Method for avoiding inter-harmonic parallel resonance of one frequency band |
-
2015
- 2015-03-04 CN CN201510097259.7A patent/CN104701864A/en active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105939020A (en) * | 2016-05-14 | 2016-09-14 | 国家电网公司 | Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability |
CN109596979A (en) * | 2018-12-06 | 2019-04-09 | 浙江大学 | A kind of synchronous generator system dynamic state estimator method based on interior optimization algorithm |
CN109596979B (en) * | 2018-12-06 | 2020-04-17 | 浙江大学 | Synchronous power generation system dynamic state estimation method based on interior point optimization algorithm |
CN110707754A (en) * | 2019-08-28 | 2020-01-17 | 广东工业大学 | Optimization method for water, wind and light power supply capacity configuration in micro-grid |
CN110707754B (en) * | 2019-08-28 | 2022-11-22 | 广东工业大学 | Optimization method for water-wind-light power supply capacity configuration in micro-grid |
CN111598718A (en) * | 2020-04-03 | 2020-08-28 | 广东电网有限责任公司电力调度控制中心 | STATCOM configuration planning method, device and equipment |
CN111598718B (en) * | 2020-04-03 | 2021-11-30 | 广东电网有限责任公司电力调度控制中心 | STATCOM configuration planning method, device and equipment |
CN112072662A (en) * | 2020-08-28 | 2020-12-11 | 武汉大学 | Method for avoiding inter-harmonic parallel resonance of one frequency band |
CN112072662B (en) * | 2020-08-28 | 2022-02-22 | 武汉大学 | Method for avoiding inter-harmonic parallel resonance of one frequency band |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vahidnia et al. | Transient stability improvement through wide-area controlled SVCs | |
CN103036230B (en) | A kind of Dynamic Equivalence of the alternating current-direct current series-parallel connection bulk power grid based on engineer applied | |
CN104701864A (en) | Reactive power planning method using SVC as reactive power compensation equipment | |
Xia et al. | Direct nonlinear primal–dual interior-point method for transient stability constrained optimal power flow | |
Zhao et al. | Power system dynamics as primal-dual algorithm for optimal load control | |
Elliott et al. | Real power modulation strategies for transient stability control | |
Al-Hinai | Voltage collapse prediction for interconnected power systems | |
Jabr | Power flow based volt/var optimization under uncertainty | |
Sravanthi et al. | Critical clearing time and transient stability analysis of SCIG based wind farm with STATCOM | |
Wu et al. | Multiple DC coordinated suppression method for ultra-low frequency oscillations | |
Reddy et al. | Controlling Power Oscillations in Real and Reactive Power using Symmetrical HYBRID PFC (Power Flow Controller) | |
Kawabe et al. | Improvement of transient stability and short‐term voltage stability by rapid control of batteries on EHV network in power systems | |
Bamasak et al. | Tuning of Fractional Order PID controller for TCSC-based stabilizer using PSO for power system oscillation damping | |
Liu | Planning reactive power control for transmission enhancement | |
Amer et al. | SCG stability enhancement using STATCOM based-ANN controller | |
Prasad et al. | Transient stability enhancement of multi-machine power system using fuzzy controlled TCSC | |
CN105514999A (en) | Prevention and control model for power grid static voltage stabilization and algorithm thereof | |
Liu et al. | Robust wide‐area damping controller design for inter‐area oscillations with signals' delay | |
Calle | Transient stability constrained optimal power flow: Improved models and practical applications | |
Yuan et al. | Stochastic hybrid system modeling and state estimation of modern power systems under contingency | |
Pavlovsky et al. | Grids transfer capacity: calculation methodology and features | |
Jin et al. | Coordinated dynamic voltage stabilization based on model predictive control | |
Kawabe et al. | Stability boundary on PV plane for analysis of short-term voltage stability | |
Yi et al. | Transient stability model and preventive control based on phase shifting transformer in power system | |
CN103617361B (en) | A kind of 10kV power distribution network dynamic power flow on-line calculation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150610 |