Summary of the invention
The objective of the invention is difference at the dynamic frequency characteristic of subregion electrical network and interconnected network, system is rationally simplified to interconnected network, the off-load process is carried out reasonable layering, be implemented under the failure condition off-load fast and effectively, the invention provides a kind of low frequency deloading method that is used for the interconnected network frequency stabilization.
For realizing purpose of the present invention, the present invention adopts following proposal to be achieved:
A kind of low frequency deloading method that is used for the interconnected network frequency stabilization, its improvements are that described low frequency deloading method may further comprise the steps:
A, input interconnected network system parameters;
B, described interconnected network system is carried out subregion;
C, the people having the same aspiration and interest unit in the described subregion is carried out people having the same aspiration and interest equivalence;
D, the UFLS of subregion is adjusted;
E, described interconnected network system is carried out the dynamic frequency characteristic analysis;
F, determine load shedding position and load shedding quantity;
G, the interconnected network system of adjusting UFLS scheme.
A kind of optimized technical scheme provided by the invention is: in the described steps A, described interconnection network system parameter comprises electric network composition, load level, installed capacity, spinning reserve capacity, electric generator structure and generator parameter, load configuration and generator parameter, utilizing the trend computational tool to carry out trend calculates, the trend that obtains base regime is separated assumed load growing direction t=0.
Second kind of optimized technical scheme provided by the invention is: among the described step B, according to described interconnected network system parameters described interconnected network system is carried out subregion.
The third optimized technical scheme provided by the invention is: among the described step C, described people having the same aspiration and interest equivalence is that an equivalent unit participates in calculating, and the frequency characteristic of described equivalent unit is
The 4th kind of optimized technical scheme provided by the invention is: among the described step D, the UFLS of described subregion is realized with loadshedding equipment; Described loadshedding equipment adopts the UFLS pattern of " basic wheel+reserve wheel+urgent wheel ".
The 5th kind of optimized technical scheme provided by the invention is: in the described step e, the dynamic frequency characteristic of described interconnected network system is analyzed dynamic frequency characteristic Δ W (s)=S (MS of described interconnected network system on the basis that subregion is simplified
2+ DS+ ω
0J)
-1(Δ P
m(s)+L Δ P
L(s)).
The 6th kind of optimized technical scheme provided by the invention is: in the described step F, determine described load shedding position and load shedding quantity by the variation of voltage.
The 7th kind of optimized technical scheme provided by the invention is: among the described step G, determine load shedding position and the quantity interconnected network system UFLS scheme of adjusting according to the dynamic frequency characteristic of described interconnected network system and described change in voltage.
Compared with prior art, the beneficial effect that reaches of the present invention is:
A kind of low frequency deloading method that is used for the interconnected network frequency stabilization provided by the invention carries out people having the same aspiration and interest equivalence at the complexity of interconnected network UFLS to the subregion electrical network, analyzes the frequency characteristic of interconnected network on the basis of equivalence; At the difference of the dynamic frequency characteristic of subregion electrical network and interconnected network, adopt the method for layering off-load, the layering off-load is at first preferentially cutting corresponding load at the interconnected network subregion that is subjected to big disturbance; Secondly differentiate interconnected network load shedding amount and load shedding position automatically by the interconnected network system; And determine off-load quantity and off-load position according to the variation of voltage; Interconnected network is rationally simplified, the off-load process is carried out reasonable layering, realized under failure condition off-load fast and effectively.
Embodiment
Below in conjunction with accompanying drawing the specific embodiment of the present invention is described in further detail.
Fig. 1 is the flow chart that is used for the low frequency deloading method of interconnected network frequency stabilization of the present invention, and this low frequency deloading method may further comprise the steps:
A, input interconnected network system parameters;
B, the interconnected network system is carried out subregion;
C, the people having the same aspiration and interest unit in the subregion is carried out people having the same aspiration and interest equivalence;
D, the UFLS of subregion is adjusted;
E, the interconnected network system is carried out the dynamic frequency characteristic analysis;
F, determine load shedding position and load shedding quantity;
G, the interconnected network system of adjusting UFLS scheme.
Being implemented as follows of each step:
Because the factor differences such as electric network composition, load level, installed capacity, spinning reserve capacity, electric generator structure and parameter thereof, load configuration and parameter thereof between the electrical network of each department, thereby cause frequency characteristic difference.At first import the interconnected network system parameters, interconnected network is carried out subregion; The foundation of subregion is to utilize similar the classification and abbreviation of people having the same aspiration and interest group of planes frequency dynamic characteristic.These othernesses according to electrical network are divided into different subregions with electrical network, and in each subregion dynamic process, relative angle changes not too large generator between them, is called a people having the same aspiration and interest group of planes, they are merged into an equivalent unit participate in calculating.The instrument that subregion of the present invention adopts is that BPA power system analysis program means PSD-BPA carries out trend calculating, obtains the combination of a people having the same aspiration and interest group of planes.The equation of rotor motion of supposing j platform unit in the subregion is (1) formula:
The total h platform unit of the unification of setting up departments obtains (2) formula with all generator amature equation of motion additions that move units:
Because
Wherein
M
jInertia constant for unit j.
The inertia time constant that obtains equivalent generator is the inertia time constant sum of each generator reduction to reference power, that is:
S in the formula
jSpecified generated output for unit j; S
sSummation for the rated output power of all units of moving in the system.
So the equation of rotor motion of equivalent generator can be expressed as:
Be the active power perunit value of the power supply of remaining system,
Active power perunit value for the remaining system load.
Because
So formula can be written as
The frequency dynamic characteristic of interconnected network subregion electrical network can be obtained by (6) formula.The derivation of the frequency dynamic characteristic of whole interconnected network is as follows.The dynamic frequency characteristic of the interconnected network that different subregions are formed is drawn by following derivation.Derivation is based on following hypothesis:
(1) linearized system model;
(2) the meritorious trend of system is mainly determined by the phase angle of voltage in the network, and is little with its magnitude relation;
(3) classical generator model is not considered the effect of pressure regulation.
Obtaining the generator amature equation of motion is:
Network equation is:
In the formula: Δ P
G, Δ δ is the N dimensional vector;
Δ P
L, Δ θ is the K dimensional vector.
The H battle array satisfies following:
Get from formula (9) cancellation Δ θ according to formula (10):
ΔP
G=JΔ5+LΔP
L (11)
In the formula:
Each row element sum of J battle array is zero
Comprehensively (8), (9), (10) and (11) get:
In the formula:
M=diag[M
1…M
N] (19)
D=diag[D
1…D
K] (20)
Set up departments and original machine power disturbance Δ P takes place in the system
mWith load disturbance Δ P
L, carry out Laplace transformation and cancellation Δ δ gets:
ΔW(s)=S(MS
2+DS+ω
0J)
-1(ΔP
m(s)+LΔP
L(s)) (21)
Formula (21) is the Mathematical Modeling that adopts when the UFLS scheme is closed in complicated interconnected systems frequency dynamic process and school.
Above-mentioned is the derivation of complicated interconnected systems frequency dynamic, determines to determine in the frequency change process position and the quantity of UFLS below.The method that the present invention adopts is to determine the quantity and the position of load according to the variation of voltage.The electromotive force of generator is relevant with generator speed, and when system frequency reduced, the voltage of system also decreased; When generator voltage reduces along with the frequency reduction, make the idle minimizing of exerting oneself, then make system's reactive power vacancy bigger.Be subjected to the disturbance area the serious decline of voltage will occur,, determine the quantity and the position of load shedding by the variation of voltage in order to keep the constant of frequency.
Trend accounting equation commonly used is:
P
i=V
i∑
j∈iV
j(G
ijcosθ
ij+B
ijsinθ
ij)(i=1,2,…n) (22)
Q
1=V
i∑
j∈iV
j(G
ijsinθ
ij-B
ijcosθ
ij)(i=1,2,…n) (23)
In the formula: G
IjAnd B
IjFor the electricity of circuit is led and susceptance.
V
iAnd θ
iVoltage magnitude and phase angle for node i.
J ∈ i represents the node j that all are directly related with node i, comprises the situation of j=i.
In the formula:
Node power deviation equation was identical during node power increment equation and trend were calculated, but Δ P and Δ Q formula power deviation amount in node power deviation equation, the sign that they go to zero and finish as iteration in the Practical Calculation Chinese style, Δ P in the node power increment equation and Δ Q formula system power increment, the existence of this increment directly influences the variation of system frequency and voltage.
(24) formula is simplified: the reactance of each element is supposed i node active power amount of unbalance Δ P much larger than resistance in the electric power networks
iVariation only cause each node voltage phase angle θ
i, θ
jVariation, do not cause V
i, V
jVariation; The uneven Δ Q of reactive power
iVariation only cause each node V
i, V
jVariation, do not cause its phase angle θ
i, θ
jVariation, the N in the formula (24), J submatrix are omitted and it are reduced to:
Can get by (25) formula:
[ΔQ]=[0?L][ΔV/V] (26)
By formula (26) as can be known, the variation of node voltage will provide Useful Information, in order to determine to be subjected to the disturbance area, be necessary the decline in short-term of measuring voltage.
Will
Substitution formula (26):
Here only consider the variation of voltage in very short time at initial stage that disturbance takes place, because computing time is short, approximately think that reactive power flow does not also shift, the reactive power vacancy of i node only influences the voltage of this node.
The idle disturbance of node and the relation between the meritorious disturbance:
ΔP
i=K
iΔQ
i (31)
In the time of normal operation, each node power factor of system is generally 0.8~0.9, does not have unsteadyly too greatly, thinks K so can be similar to
1=K
2=... K
N, can get by top hypothesis:
ΔP=∑P
LS (35)
In the formula: Δ V
iVoltage decline minimum for node i;
N is all nodes that big change in voltage is arranged;
∑ P
LSBe the total load shedding power of system.
Obtain the cutting load amount of each comparatively serious node of voltage drop according to following formula, behind the corresponding load of each node excision, finally make system frequency be stabilized in one preferably on the level.
More than subregion frequency characteristic, interconnected systems frequency characteristic have been carried out detailed inference, also load shedding point and load shedding amount have been carried out quantizing to calculate.Next whole interconnected systems low frequency load shedding method is described, the load shedding method adopts layering load shedding method, by point and face, realizes the rapidity and the science of UFLS.
The layering cutting load is divided into the subregion ground floor and the interconnected systems second layer.
At the disturbance initial stage, because generated output withdraws from a large number, frequency sharply descends, and the result of calculation that low-frequency load reduction control system has little time to provide definite is in order to determine the position and the quantity of cutting load.The loadshedding equipment that be distributed in subregion this moment can adopt the UFLS pattern of " basic wheel+reserve wheel+urgent wheel " in time to move, and alleviates frequency and descends.Basic wheel moves and takes turns level, respectively takes turns particularly first run operating frequency of operating frequency, determines by common power shortage mode, promptly takes turns to be defined as according to the df/dt size and quickens the round that cutting load constitutes.When the basic wheel first round starts, quicken to cut basic wheel second and take turns and be defined as the urgent first round, take turns the first round substantially and quicken to cut basic wheel when starting second, third is taken turns and is defined as urgent second and takes turns, and the like.The operating frequency of reserve wheel should be not less than the starting frequency of taking turns the first round substantially, and takes certain time-delay, treats to move when system frequency is more stable, and frequency retrieval is arrived more than the limit value of permission, below is subregion load shedding scheme:
Discriminant when (1) frequency slowly descends is as follows:
F≤f
Qd, t 〉=t
QdLow-frequency start;
F≤f
Dz1, t 〉=t
Dz1Basic wheel first round action;
F≤f
Dz2, t 〉=t
Dz2Basic wheel second is taken turns action;
F≤f
Dz3, t 〉=t
Dz3The action of basic wheel third round;
F≤f
Dz4, t 〉=t
Dz4The action of basic wheel four-wheel.
(2) big, the frequency of the meritorious vacancy discriminant when very fast that descends
F≤f
Qd, t 〉=t
QdLow-frequency start;
f≤f
dz1,t≥t
dz1
Df/dt<(df/dt)
Dz1Basic wheel first round action;
(df/dt)
Dz1<df/dt<(df/dt)
Dz2Urgent first round action (cut substantially and take turns the first round, quicken to cut second simultaneously and take turns);
(df/dt)
Dz2<df/dt<(df/dt)
Dz3Urgent second takes turns action (cut substantially and take turns the first round, quicken to cut second and third simultaneously and take turns).
(3) recover underfrequency, the discriminant when frequency is hovered
F≤f
Qd, t 〉=t
QdLow-frequency start;
F≤f
Hb, t 〉=t
HbThe action of back wheel.
(4) prevent to take turns substantially the discriminant that the 3rd class " is crossed and cut "
The open basic wheel outlet in df/dt<0 (system frequency is in the decline process);
Outlet (system frequency is in uphill process) is taken turns in df/dt>0 locking substantially.
Wherein, f
Qd, t
QdBe the low-frequency start definite value and the time-delay definite value of this device, f
Dzi, t
DziBe the action definite value and the time-delay definite value of basic wheel i wheel, f
Hb, t
HbBe the operating value and the time-delay definite value of reserve wheel, (df/dt)
Dz1, (df/dt)
Dz2Be the urgent first round and the urgent second action definite value of taking turns, (df/dt)
Dz2Be slippage locking definite value.
In case big disturbance subregion low frequency load shedding equipment action alleviates frequency and descends, position and the quantity of sufficient time by low-frequency load reduction control system calculating load shedding is just arranged, below be total load shedding total amount and load shedding position:
If system has sufficient spinning reserve capacity (SR), the load shedding total amount is:
∑P
LS=ΔP
G-K
LΔf
ss-SR (37)
If do not have sufficient spinning reserve capacity in the system, the load shedding total amount is:
∑P
LS=ΔP
G-K
LΔf
ss (38)
∑ P in the formula
LSBe load shedding power; Δ P
GThe generated output that withdraws from for accident; K
LBe the Load Regulation effect coefficient.
In order to calculate the load shedding position, at first to determine the load shedding scope.The voltage of this node descends and only continues the very short time, but the voltage minimum point can go on record, and the load shedding power of node i can be expressed as:
EMS in the modern power systems (being called for short EMS) provides the data in formula (37), (38) and (39).
Low-frequency load reduction control system obtains corresponding data as required and calculates from the EMS system server, determine big disturbance area according to result of calculation, and sends the cutting load order to the low frequency load shedding equipment of this area, informs the quantity of cutting load.
Load shedding power also may appear at those and deducted in the subregion load shedding on the node of a large amount of loads this moment, in other words, whole load shedding scheme is applicable to the most serious power shortage and two kinds of situations of provincialism power shortage, finally makes system frequency be stabilized in one preferably on the level.
Should be noted that at last: only illustrate that in conjunction with the foregoing description technical scheme of the present invention is not intended to limit.Those of ordinary skill in the field are to be understood that: those skilled in the art can make amendment or are equal to replacement the specific embodiment of the present invention, but these modifications or change are all among the claim protection range that application is awaited the reply.