CN105978003A - Power system additional wide-area damping controller design method giving consideration to time lag - Google Patents
Power system additional wide-area damping controller design method giving consideration to time lag Download PDFInfo
- Publication number
- CN105978003A CN105978003A CN201610489302.9A CN201610489302A CN105978003A CN 105978003 A CN105978003 A CN 105978003A CN 201610489302 A CN201610489302 A CN 201610489302A CN 105978003 A CN105978003 A CN 105978003A
- Authority
- CN
- China
- Prior art keywords
- power system
- time lag
- avr
- lag
- wadc
- 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.)
- Granted
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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a power system additional wide-area damping controller design method giving consideration to time lag, and the method comprises the steps: (1), enabling an oscillation mode with the weakest damping of a power system to serve as a target control mode, determining the optimal installation positions of AVR and FACTS equipment, further building a power system mode and a WADC model which comprise the AVR and FACTS equipment and do not give consideration to time lag; (2), introducing a time lag link on the basis of the power system mode and the WADC model which comprise the AVR and FACTS equipment and do not give consideration to time lag, and further building a closed loop time lag power system model giving consideration to a time lag factor; (3), building a WADC parameter-optimization mathematic model according to the closed loop time lag power system model; (4), employing a PSO algorithm, and solving the optimal parameters of WADC at step (3). The method gives consideration to the time lag factor in each iteration optimization process, employs the SOD algorithm to solve the key characteristic value of the right side of a closed loop time lag power system, so as to achieve a big purpose that the damping of the key characteristic value is the biggest.
Description
Technical field
The present invention relates to a kind of power system additional wide area damping control method for designing considering time lag.
Background technology
Mechanism that extensive interconnected network is complicated and severe service condition so that it is overall dynamic behaviour be difficult to analyze and
Control.The low-frequency oscillation repeatedly occurred in interconnected electric power system, threatens the safe and stable operation of electrical network, restricts region
The ability to transmit electricity on some major transmission line road interior and interval so that the initial target that trans-regional electric energy dispatches this Power System Interconnection is difficult
To complete.
Along with WAMS (Wide Area Measurement System, WAMS) technology development with become
Ripe, WAMS is used for extensive interconnected electric power system and forms closed loop control, be one of the direction of power system development.WAMS can
There is provided the information such as through-put power of generator's power and angle, rotating speed and interregional interconnection in real time, for design wide area damping control
With suppression system inter-area oscillations, it is provided that new data platform.Wide area signal based on WAMS, fully reflection system local oscillation
With the modal information of inter-area oscillations, provide it to damping controller and constitute closed loop control, it is possible to obtain superior control performance.
But, application Wide-area Measurement Information brings new problem also can to the control of system.Modern digital communication network is
The transmission of WAMS information provides strong physical support, although communication for information is convenient and swift, but, it is contemplated that transmission range is relatively
Far, Wide-area Measurement Information is transmitted in a communication network and be there is inevitable communication delay.When time lag is bigger, may result in damping
Controller action effect is deteriorated, even induction power system unstability.
At present, it is considered to the wide area stability contorting of time-delay can be largely classified into wide area robust PSS and controls and wide area damp
Control two big classes.Wide area robust PSS controls to primarily focus on robust control theory to be applied to power system, by means of robust control
Workbox processed can complete the design of single damping controller, the most existing sensitive based on loop molding, gain scheduling and mixing
The Wide-area Time-delay damping controller method for designing of the technology such as degree.But in the design process, it is respectively provided with not due to system and time lag
Definitiveness, needs system model and Time Delay are carried out depression of order and approximate processing.Model reduction and the error of approximation link, directly
Connect the performance affecting robust controller.
For the design of wide area damping control (Wide Area Damping Controller, WADC), mainly pin
To the power system stabilizer, PSS (Power System Stabilizer, PSS) in electrical network, automatic voltage regulator
(Automatic Voltage Regular, AVR) and typical flexible AC transmitting system (Flexible Alternative
Current Transmission Systems, FACTS) additional damping controller of equipment (SVC, TCSC etc.) is optimized and sets
Meter.Current research is mostly based on single, the power system of fixed time lag, system exists multiple, stochastic Time-Delay situation and grinds
Study carefully the most less.These researchs are the most all to use Method of Pole Placement design parameter, time lag factor is introduced Optimized model, or sets
Count suitable time lag compensation link, to improve robustness and the damping of controller, and have no special with power system key
Value indicative is calculated as the WADC on basis and optimizes the achievement in research designed.Calculate additionally, existing based on WAMS and Modified particle swarm optimization
The SVC additional damping controller Optimization Design of method, obtains the relevant information of oscillation mode by Prony analytical tool, should
By modified particle swarm optiziation optimal control parameter, but the method does not accounts for the time delay of Wide-area Measurement Information to controlling effect
Impact.
Summary of the invention
In order to solve the shortcoming of prior art, the present invention provides a kind of power system additional wide area damp control considering time lag
Device method for designing processed.When the method calculates the stability of the critical eigenvalue of real system and judgement system, take into full account communication
The time lag impact on system damping performance, can be greatly improved the damping of inter-area oscillation mode, strengthens the stability of time lag system.
For achieving the above object, the present invention is by the following technical solutions:
A kind of power system additional wide area damping control method for designing considering time lag, comprises the following steps:
Step (1): using oscillation mode the most weak for power system damping as target control pattern, determines that AVR and FACTS sets
Standby best position, further foundation do not consider time lag factor the electric power system model comprising AVR and FACTS equipment and
WADC model;
Step (2): do not considering the electric power system model comprising AVR and FACTS equipment and the WADC model of time lag factor
On the basis of, introduce Time Delay, and then set up the model of the closed loop time-lag power system considering time lag factor;
Step (3): according to the model of closed loop time-lag power system, build the mathematical model of WADC parameter optimization;
Step (4): use particle group optimizing (Particle Swarm Optimization, PSO) algorithm, solution procedure
(3) optimized parameter of WADC in;
In each iterative optimization procedure, it is considered to time lag factor, use SOD Algorithm for Solving closed loop time-lag power system
The critical eigenvalue on right side, to reach the damping of critical eigenvalue this target maximum.
In described step (1), use associating geometric measures method based on controllability and observability, determine AVR and FACTS equipment
Best position.
In described step (1), use associating geometric measures method based on controllability and observability, determine AVR and FACTS equipment
The detailed process of best position, including:
Step (1.1): the parameter of given AVR and FACTS equipment, sets some node of power system as AVR and FACTS
Standby alternative installation site;
Step (1.2): calculate containing installing the system load flow equation of AVR and FACTS equipment and the initial of state variable
Value, forms system linearity coefficient matrix;According to linearisation coefficient matrix, calculate the critical eigenvalue of power system and right
The left and right characteristic vector answered, and set the oscillation mode to be investigated of power system;
Step (1.3): under the oscillation mode set, with the reference value of AVR and FACTS equipment for input, with different
Feedback signal is output, calculates the associating geometric measures of the controllability and observability of AVR and FACTS equipment;
Step (1.4): compare the modulus value size of the associating geometric measures of controllability and observability, filter out associating geometric measures
The node of modulus value maximum is defined as the best position of AVR and FACTS equipment, using the control signal of its correspondence as additional
The feedback signal of WADC, input WADC constitutes closed loop control, thus increases the damping capacity of system.
In step (1.4), the address installed according to AVR and the FACTS equipment selected and feedback signal, set up following system
System linearized state-space model:
In formula: x is the state vector of system, u and y is respectively input vector and output vector;A, B, C, D are respectively system
State matrix, input matrix, output matrix and straight-through matrix.If the eigenvalue of matrix A is λi(i=1 ..., n), accordingly
Left and right eigenvectors matrix is: U=[u1,u2,…,un], V=[v1,v2,…,vn].Wherein, uiAnd viIt is respectively ith feature
The left and right characteristic vector of value.After standardization processing, U and V meets UHV=VHU=In, InFor n rank unit matrix.
System i-th oscillation mode λiControllability geometric measures mciWith controllability geometric measures moiCan be calculated by following formula:
In formula: bkKth for input matrix B arranges;clL row for output matrix C;α(ui,bk) it is input vector bkAnd a left side
Characteristic vector uiAngle;θ(cl T,vi) it is output vector clWith right characteristic vector viAngle;| | and | | | | represent respectively
Delivery value and European norm.
For oscillation mode λi, can control/the associating geometric measures of controllability is:
mcoi(k, l)=mci(k)moi(l)
In step (1.4), according to (1.3) calculated can the associating geometric measures of control/controllability determine AVR with
The best position of FACTS equipment.If mcoi(k, l) ≠ 0, then explanation can pass through ukAnd ylControl model λi。mcoi(k l) takes
U during maximumkAnd ylIt is maximally effective;Work as mcoi(k, when l) taking maximum, if ukAnd ylSignal take from the same area, shake
Swing pattern λiCan be controlled by this locality damping;Otherwise, if ukAnd ylSignal take from different regions, then need with wide-area control come
Damp corresponding oscillation mode.
In described step (1), the structure of additional WADC is similar to tradition PSS.It is made up of lead-lag link, belongs to dynamic
Compensator.Apply it to form closed loop control, the damping capacity of power system can be promoted.
In described step (2), time-lag power system model is as follows:
In formula: Δ x is the state variable of system.τ=[τ1,…,τi,…,τm]T, τi> 0 it is the transmission of i-th wide area signal
Time lag, i=1,2 ..., m, maximum of which time lag is expressed as τmax。It is the state matrix of system, for dense matrix;It is the hangover state matrix of system, for sparse matrix.Δ x (t) is the increment of t system state variables,
Δx(t-τi) it is t-τiThe increment of moment system state variables,Increment for t system state variables derivative.Δx
(0) it is the initial value (i.e. initial condition) of system state variables, and is abbreviated asThe feature of the time-lag power system that above formula represents
Equation is:
In formula: λ is characterized value, v is characterized the right characteristic vector that value is corresponding.
In described step (3), additional WADC Parametric optimization problem mathematical model is expressed as:
Max J, J=min{ ζi, i ∈ electromechanic oscillation mode set }
In formula: ζiFor the damping ratio of system i-th oscillation mode, cost function J is expressed as system whole electromechanical oscillations mould
Damping ratio minimum in formula, KdFor the amplification of the Control of Voltage link of AVR and FACTS equipment, KaFor additional damping controller
Amplification, T1~T4Time constant for lead-lag link;WithIt is respectively amplification KdHigher limit and
Lower limit;WithIt is respectively amplification KaHigher limit and lower limit;WithWhen being respectively
Between constant T1~T4Higher limit and lower limit.
In described step (4), use PSO algorithm, before the WADC optimized parameter in solution procedure (3), also include: arrange
Population number and maximum iteration time also initialize population.
In described step (4), according to spectral mapping theorem, the eigenvalue being positioned at left half complex plane of closed loop time-lag power system
Being mapped to Solution operator and be positioned at the unit circle of z-plane, the eigenvalue being positioned at right half complex plane is mapped in outside unit circle.
In described step (4), use equally spaced linear multistep method (Linear Multi-Step, LMS) that Solution operator is entered
Row processes, and according to eigenvalue λ and eigenvalue (composing) the μ corresponding relation of Solution operator of time-lag power system, tries to achieve time lag power train
The critical eigenvalue of system.
Solution operator T (h): X → X is defined as the initial condition (state) in the θ moment in space XBeing mapped to the h+ θ moment is
The linear operator of system state ψ.
In formula: h is transfer step-length, 0≤h≤τmax;xtThe part solution of etching system when=x (t+ θ) is t >=0.Mapped by spectrum
Theorem understands, and has following relation between spectrum μ and the eigenvalue λ of time lag system of Solution operator T (h):
In formula: σ () represents spectrum, represent eliminating.
By semigroup of operators theory, the Solution operator in Banach space is infinite dimensional.In order to calculate Solution operator
Eigenvalue, needs T (h) is carried out discretization.After discretization, obtain corresponding with Solution operator, a finite dimensional approximate matrix,
The critical eigenvalue of former time lag system is can be obtained by by calculating the eigenvalue of approximate matrix.
Based in SOD Algorithm for Solving time-lag power system critical eigenvalue in described step (4), use equally spaced LMS
Method carries out approximate processing to differential equations with delay (Delay Differential Equations, DDE) relevant Solution operator, can
Former problem is converted into the eigenvalue problem solving a canonical matrix.
Assume h≤τmax, Solution operator can be expressed as form:
In above formula, the first row of right-hand side expression is initial-value problem, and LMS method can be used to solve.
Use equally spaced some θj=jh, j=-N ..., 0, Solution operator is carried out discretization, N is more than or equal to τmax/h
Smallest positive integral, i.e.
The expression formula of LMS method is:
At mesh point, with ψjApproximation replaces ψ (θj), withApproximation replacesThen ψjWithMathematical relationship express
Formula:
In formula: TNDiscretization matrix for Solution operator T (h), it may be assumed that
TNLast block row Γ be the coefficient matrix of polynomial eigenvalue problem, be specifically represented by:
Especially, when system contains only a time lag and N=τ/h is integer, Γ can explicitly be expressed as follows:
Γ=[Γ0,0n×n(N-k-1),Γ1]
In each iterative process of PSO algorithm of described step (4), the search information of the most more new particle:
In formula: k is the number of times that current iteration calculates,D for particle i search speed ties up component, WithRepresent the upper and lower bound of search speed respectively,The d projecting to search volume for i-th particle ties up component;
pbestdOptimum position for i-th particle up till now projects to the d dimension component of search volume;gbestdFor current all grains
The optimum position of son search projects to the d of search volume and ties up component;c1And c2For aceleration pulse;Rand () and Rand ()
It is the random number between two 0 to 1;ω is weight coefficient.
For improving convergence rate, linear narrowing weight coefficient in each iterative process of PSO algorithm of step (4).
In formula: k is current iteration number of times;K is as ω=ωminTime algorithm iterations;ωmaxAnd ωminAt the beginning of Fen Bieweiing
The weight coefficient maximum arranged during beginningization algorithm and minima.
The invention have the benefit that
(1) the method for the present invention is for designing the additional WADC of AVR and FACTS equipment, interval can be greatly improved and shake
Swing the damping of pattern, strengthen the stability of time lag system;Calculating the critical eigenvalue of real system and stablizing of judgement system
During property, take into full account the communication delay impact on system damping performance;
(2) method of the present invention combines SOD and PSO both algorithms, and designed additional wide area damping control is same
Original local PSS is mutually coordinated, the most there is not the reciprocal action being unfavorable for system damping performance;
(3) method of the present invention is with critical eigenvalue for the additional wide area damping control of design, it is to avoid traditional
Wide area power system regulator based on POLE PLACEMENT USING and Robust Damping controller design method need to drop system model
Rank and time lag carries out the deficiency of approximate processing, purposiveness is strong, control device is direct, action effect is obvious, and the most any
Conservative.
Accompanying drawing explanation
Fig. 1 is time-lag power system schematic diagram;
Fig. 2 is the SVC stability analysis model considering additional wide-area damping control;
Fig. 3 is the TCSC stability analysis model installing additional wide-area damping control;
Fig. 4 is the calculation process of the controlled/ornamental associating geometric measures of target pattern;
Fig. 5 is additional WADC design flow diagram based on PSO algorithm;
Fig. 6 is time-lag power system based on the SOD algorithm additional WADC method for designing flow chart of the present invention.
Detailed description of the invention
The present invention will be further described with embodiment below in conjunction with the accompanying drawings:
As shown in Figure 1: in the modeling process of actual large-scale electrical power system, introduce Time Delay.Time-lag power system bag
Exporting time lag four part containing without time-lag power system, wide area Feedback Delays, additional WADC and wide area, the connection between each several part is closed
System is as shown in the figure.In FIG, yfFor the output without time-lag power system, ydfFor the wide area feedback signal after consideration Feedback Delays
And as the input of damping controller, ycFor the output of wide area damping control, ydcWide area damp delayed during for considering to export
The output of controller, also serves as inputting without the control of time-lag power system simultaneously.
Time lag system is positioned at the eigenvalue of left half complex plane and is mapped to Solution operator and is positioned at the unit circle of z-plane, and time lag
System is positioned at the eigenvalue of right half complex plane and is mapped in outside unit circle.Therefore, the eigenvalue of Solution operator is utilized, it is possible to solve
Corresponding time lag system critical eigenvalue.
As shown in Figures 2 and 3, the structure of additional WADC is similar with PSS, is made up of lead-lag link and belongs to dynamically benefit
Repay device.The state-space expression of additional WADC can be written as:
In formula: Δ xc=[Δ V1 ΔV2 ΔV3 ΔVIS]T;
In formula: Δ V1、ΔV2、ΔV3For the intermediate variable of WADC links, Δ VISOutput variable for WADC, it is provided that
Constituting feedback control to FACTS equipment, K is amplification, and T is time constant.
As shown in Figure 4, use associating geometric measures method based on controllability and observability, determine that AVR and FACTS equipment is
The detailed process of good installation site, including:
Step (1.1): the parameter of given AVR and FACTS equipment, sets some node of power system as AVR and FACTS
Standby alternative installation site;
Step (1.2): calculate containing installing the system load flow equation of AVR and FACTS equipment and the initial of state variable
Value, forms system linearity coefficient matrix;According to linearisation coefficient matrix, calculate the critical eigenvalue of power system and right
The left and right characteristic vector answered, and set the oscillation mode to be investigated of power system;
Step (1.3): under the oscillation mode set, with the reference value of AVR and FACTS equipment for input, with different
Feedback signal is output, calculates the associating geometric measures of the controllability and observability of AVR and FACTS equipment;
Step (1.4): compare the modulus value size of the associating geometric measures of controllability and observability, filter out associating geometric measures
The node of modulus value maximum is defined as the best position of AVR and FACTS equipment, using the control signal of its correspondence as additional
The feedback signal of WADC, input WADC constitutes closed loop control, thus increases the damping capacity of system.
As Fig. 5 shows, additional WADC based on PSO algorithm designs and specifically includes following steps:
Step (4.1): population number and maximum iteration time are set, and initialize population
Step (4.2): Dynamic simulation model;
Step (4.3): consider time lag factor, calculates the rightmost side critical eigenvalue of time-lag power system based on SOD algorithm
And damping ratio, design controller parameter;
Step (4.4): Population Regeneration, it may be judged whether reach the iterations upper limit or meet required precision, if it is not,
Then return step (4.2) to recalculate, if it is, output parameters optimization.
In described step (4.4), additional WADC Parametric optimization problem mathematical model is expressed as:
Max J, J=min{ ζi, i ∈ electromechanic oscillation mode set }
In formula: ζiFor the damping ratio of system i-th oscillation mode, cost function J is expressed as system whole electromechanical oscillations mould
Damping ratio minimum in formula, KdFor the amplification of the Control of Voltage link of AVR and FACTS equipment, KaFor additional damping controller
Amplification, T1~T4Time constant for lead-lag link;WithIt is respectively amplification KdHigher limit and
Lower limit;WithIt is respectively amplification KaHigher limit and lower limit;WithWhen being respectively
Between constant T1~T4Higher limit and lower limit.
In described step (4.4), in each iterative process of PSO algorithm, the search information of the most more new particle:
In formula: k is the number of times that current iteration calculates,D for particle i search speed ties up component, WithRepresent the upper and lower bound of search speed respectively,The d projecting to search volume for i-th particle ties up component;
pbestdOptimum position for i-th particle up till now projects to the d dimension component of search volume;gbestdFor current all grains
The optimum position of son search projects to the d of search volume and ties up component;c1And c2For aceleration pulse;Rand () and Rand ()
It is the random number between two 0 to 1;ω is weight coefficient.
For improving convergence rate, linear narrowing weight coefficient in each iterative process of PSO algorithm of step (4).
In formula: k is current iteration number of times;K is as ω=ωminTime algorithm iterations;ωmaxAnd ωminAt the beginning of Fen Bieweiing
The weight coefficient maximum arranged during beginningization algorithm and minima.
As shown in Figure 6, a kind of power system additional wide area damping control method for designing considering time lag of the present invention, bag
Include following steps:
Step (1): using oscillation mode the most weak for power system damping as target control pattern, determines that AVR and FACTS sets
Standby best position, further foundation do not consider time lag factor the electric power system model comprising AVR and FACTS equipment and
WADC model;
Step (2): do not considering the electric power system model comprising AVR and FACTS equipment and the WADC model of time lag factor
On the basis of, introduce Time Delay, and then set up the model of the closed loop time-lag power system considering time lag factor;
Step (3): according to the model of closed loop time-lag power system, build the mathematical model of WADC parameter optimization;
Step (4): use PSO algorithm, the optimized parameter of WADC in solution procedure (3).At each iterative optimization procedure
In, it is considered to time lag factor, use the critical eigenvalue of the SOD Algorithm for Solving closed loop time-lag power system rightmost side, to reach crucial
The damping of eigenvalue this target maximum.
In specific implementation process, in step (1), use controllability and observability associating geometric measures method, determine AVR and
The best position of FACTS equipment.
The address installed according to AVR and the FACTS equipment selected and feedback signal, set up following system linearization state empty
Between model:
In formula: x is the state vector of system, u and y is respectively input vector and output vector;A, B, C, D are respectively system
State matrix, input matrix, output matrix and straight-through matrix.If the eigenvalue of matrix A is λi(i=1 ..., n), accordingly
Left and right eigenvectors matrix is: U=[u1,u2,…,un], V=[v1,v2,…,vn].Wherein, uiAnd viIt is respectively ith feature
The left and right characteristic vector of value.After standardization processing, U and V meets UHV=VHU=In, InFor n rank unit matrix.
System i-th oscillation mode λiControllability geometric measures mciWith controllability geometric measures moiCan be calculated by following formula:
In formula: bkKth for input matrix B arranges;clL row for output matrix C;α(ui,bk) it is input vector bkAnd a left side
Characteristic vector uiAngle;θ(cl T,vi) it is output vector clWith right characteristic vector viAngle;| | and | | | | represent respectively
Delivery value and European norm.
For oscillation mode λi, can control/the associating geometric measures of controllability is:
mcoi(k, l)=mci(k)moi(l)
If mcoi(k, l) ≠ 0, then explanation can pass through ukAnd ylControl model λi。mcoi(k l) takes u during maximumkAnd ylIt is
Maximally effective;Work as mcoi(k, when l) taking maximum, if ukAnd ylSignal take from the same area, oscillation mode λiCan be by this locality
Damping controls;Otherwise, if ukAnd ylSignal take from different regions, then need to damp corresponding oscillation mode with wide-area control.
Although the detailed description of the invention of the present invention is described by the above-mentioned accompanying drawing that combines, but not the present invention is protected model
The restriction enclosed, one of ordinary skill in the art should be understood that on the basis of technical scheme, and those skilled in the art are not
Need to pay various amendments or deformation that creative work can make still within protection scope of the present invention.
Claims (10)
1. the power system additional wide area damping control method for designing considering time lag, it is characterised in that include following step
Rapid:
Step (1): using oscillation mode the most weak for power system damping as target control pattern, determine AVR and FACTS equipment
Best position, further foundation does not consider the electric power system model comprising AVR and FACTS equipment and the WADC of time lag factor
Model;
Step (2): do not considering the electric power system model comprising AVR and FACTS equipment and the base of WADC model of time lag factor
On plinth, introduce Time Delay, and then set up the model of the closed loop time-lag power system considering time lag factor;
Step (3): according to the model of closed loop time-lag power system, build the mathematical model of WADC parameter optimization;
Step (4): using PSO algorithm, in solution procedure (3), the optimized parameter of WADC, in each iterative optimization procedure, examines
Consider time lag factor, use the critical eigenvalue of the SOD Algorithm for Solving closed loop time-lag power system rightmost side, to reach critical eigenvalue
Damping this target maximum.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 1, it is special
Levy and be, in described step (1), use the associating geometric measures method of controllability and observability, determine that AVR and FACTS equipment is
The Optimal Feedback signal of good installation site and additional wide area damping control.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 2, it is special
Levy and be, in described step (1), determine the detailed process of the best position of AVR and FACTS equipment, including:
Step (1.1): the parameter of given AVR and FACTS equipment, using some node of power system as AVR and FACTS equipment
Alternative installation site;
Step (1.2): calculate containing system load flow equation and the initial value of state variable, the shape installing AVR and FACTS equipment
Become system linearity coefficient matrix;According to linearisation coefficient matrix, calculate the critical eigenvalue of power system and correspondence thereof
Left and right characteristic vector, and set the oscillation mode to be investigated of power system;
Step (1.3): under the oscillation mode set, with the reference value of AVR and FACTS equipment for input, with different feedbacks
Signal is output, calculates the associating geometric measures of the controllability and observability of AVR and FACTS equipment;
Step (1.4): compare the modulus value size of the associating geometric measures of controllability and observability, filter out the modulus value of associating geometric measures
Maximum node is defined as the best position of AVR and FACTS equipment, using the control signal of its correspondence as additional WADC's
Feedback signal, input WADC constitutes closed loop control, thus increases the damping capacity of system.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 1, it is special
Levying and be, additional wide area damping control is made up of lead-lag link, belongs to dynamic compensator.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 1, it is special
Levying and be, according to spectral mapping theorem, the eigenvalue being positioned at left half complex plane of closed loop time-lag power system is mapped to Solution operator position
In the unit circle of z-plane, the eigenvalue being positioned at right half complex plane is mapped in outside unit circle.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 5, it is special
Levy and be, utilize the eigenvalue of Solution operator, solve corresponding closed loop time-lag power system critical eigenvalue.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 1, it is special
Levy and be, in described step (4), use PSO intelligent algorithm, before the WADC optimized parameter in solution procedure (3), also include:
Population number and maximum iteration time are set and initialize population.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 6, it is special
Levy and be, in described step (4), use equally spaced LMS that Solution operator is processed, according to the eigenvalue of time-lag power system
With the eigenvalue corresponding relation of Solution operator, try to achieve the critical eigenvalue of time-lag power system.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 1, it is special
Levy and be, in each iterative process of PSO algorithm of described step (4), the search information of the most more new particle:
In formula: k is the number of times that current iteration calculates,D for particle i search speed ties up component, WithRepresent the upper and lower bound of search speed respectively;The d projecting to search volume for i-th particle ties up component;pbestd
Optimum position for i-th particle up till now projects to the d dimension component of search volume;gbestdFor current all particle search
Optimum position project to search volume d tie up component;c1And c2For aceleration pulse;Rand () and Rand () is two
Random number between individual 0 to 1;ω is weight coefficient.
A kind of power system additional wide area damping control method for designing considering time lag the most as claimed in claim 1, it is special
Levy and be, for improving convergence rate, linear narrowing weight coefficient in each iterative process of PSO algorithm of step (4).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610489302.9A CN105978003B (en) | 2016-06-27 | 2016-06-27 | A kind of additional wide area damping control design method of power system for considering time lag |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610489302.9A CN105978003B (en) | 2016-06-27 | 2016-06-27 | A kind of additional wide area damping control design method of power system for considering time lag |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105978003A true CN105978003A (en) | 2016-09-28 |
CN105978003B CN105978003B (en) | 2017-08-29 |
Family
ID=57019394
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610489302.9A Active CN105978003B (en) | 2016-06-27 | 2016-06-27 | A kind of additional wide area damping control design method of power system for considering time lag |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105978003B (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529013A (en) * | 2016-11-03 | 2017-03-22 | 国网宁夏电力公司电力科学研究院 | High-voltage direct-current additional frequency control method |
CN107346892A (en) * | 2017-04-24 | 2017-11-14 | 国家电网公司 | PSS design methods based on the input hierarchical control of wide area multi signal |
CN107910880A (en) * | 2017-11-29 | 2018-04-13 | 山东大学 | Wide area damping control optimal parameter setting method based on inter-area oscillation mode |
CN108808705A (en) * | 2018-07-13 | 2018-11-13 | 山东大学 | Time-lag power system electromechanic oscillation mode computational methods based on low order SOD-PS-II-R algorithms |
CN109149601A (en) * | 2018-09-05 | 2019-01-04 | 北京建筑大学 | Wide area damp of electrical power system controller parameter determines method and system |
CN110365009A (en) * | 2019-06-22 | 2019-10-22 | 东北电力大学 | A kind of wide area damp optimal control method suitable for interconnected electric power system |
CN111181167A (en) * | 2018-11-09 | 2020-05-19 | 国网内蒙古东部电力有限公司经济技术研究院 | Method for improving stability of power grid at transmitting end by using static var compensator |
CN112510700A (en) * | 2020-11-26 | 2021-03-16 | 电子科技大学 | Controller self-adaptive robust optimization method considering wind speed change |
CN117996741A (en) * | 2024-01-30 | 2024-05-07 | 山东大学 | Double-fed fan subsynchronous oscillation suppression method and system |
CN117996741B (en) * | 2024-01-30 | 2024-08-02 | 山东大学 | Double-fed fan subsynchronous oscillation suppression method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101465550A (en) * | 2007-12-21 | 2009-06-24 | 清华大学 | Sdelayed time processing and compensating system for wide domain damped control of electric power system |
EP2299555A1 (en) * | 2009-09-21 | 2011-03-23 | ABB Research Ltd. | Fault tolerant damping of electromechanical oscillations in power systems |
CN102801158A (en) * | 2012-07-31 | 2012-11-28 | 山东电力集团公司电力经济技术研究院 | Method for calculating time-lag electric power system eigenvalue and discriminating stability based on Pade approximation |
CN105468909A (en) * | 2015-11-24 | 2016-04-06 | 山东大学 | Time delay power system electromechanical oscillation mode computing method based on SOD-PS-R R algorithm |
-
2016
- 2016-06-27 CN CN201610489302.9A patent/CN105978003B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101465550A (en) * | 2007-12-21 | 2009-06-24 | 清华大学 | Sdelayed time processing and compensating system for wide domain damped control of electric power system |
EP2299555A1 (en) * | 2009-09-21 | 2011-03-23 | ABB Research Ltd. | Fault tolerant damping of electromechanical oscillations in power systems |
CN102801158A (en) * | 2012-07-31 | 2012-11-28 | 山东电力集团公司电力经济技术研究院 | Method for calculating time-lag electric power system eigenvalue and discriminating stability based on Pade approximation |
CN105468909A (en) * | 2015-11-24 | 2016-04-06 | 山东大学 | Time delay power system electromechanical oscillation mode computing method based on SOD-PS-R R algorithm |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529013B (en) * | 2016-11-03 | 2019-03-08 | 国网宁夏电力公司电力科学研究院 | A kind of high voltage direct current additional frequency control method |
CN106529013A (en) * | 2016-11-03 | 2017-03-22 | 国网宁夏电力公司电力科学研究院 | High-voltage direct-current additional frequency control method |
CN107346892A (en) * | 2017-04-24 | 2017-11-14 | 国家电网公司 | PSS design methods based on the input hierarchical control of wide area multi signal |
CN107346892B (en) * | 2017-04-24 | 2021-01-15 | 国家电网公司 | PSS design method based on wide area multi-signal input hierarchical control |
CN107910880A (en) * | 2017-11-29 | 2018-04-13 | 山东大学 | Wide area damping control optimal parameter setting method based on inter-area oscillation mode |
CN108808705A (en) * | 2018-07-13 | 2018-11-13 | 山东大学 | Time-lag power system electromechanic oscillation mode computational methods based on low order SOD-PS-II-R algorithms |
CN108808705B (en) * | 2018-07-13 | 2020-05-22 | 山东大学 | Time-lag power system electromechanical oscillation mode calculation method based on low-order SOD-PS-II-R algorithm |
CN109149601A (en) * | 2018-09-05 | 2019-01-04 | 北京建筑大学 | Wide area damp of electrical power system controller parameter determines method and system |
CN111181167A (en) * | 2018-11-09 | 2020-05-19 | 国网内蒙古东部电力有限公司经济技术研究院 | Method for improving stability of power grid at transmitting end by using static var compensator |
CN110365009A (en) * | 2019-06-22 | 2019-10-22 | 东北电力大学 | A kind of wide area damp optimal control method suitable for interconnected electric power system |
CN112510700A (en) * | 2020-11-26 | 2021-03-16 | 电子科技大学 | Controller self-adaptive robust optimization method considering wind speed change |
CN112510700B (en) * | 2020-11-26 | 2022-10-18 | 电子科技大学 | Controller self-adaptive robust optimization method considering wind speed change |
CN117996741A (en) * | 2024-01-30 | 2024-05-07 | 山东大学 | Double-fed fan subsynchronous oscillation suppression method and system |
CN117996741B (en) * | 2024-01-30 | 2024-08-02 | 山东大学 | Double-fed fan subsynchronous oscillation suppression method and system |
Also Published As
Publication number | Publication date |
---|---|
CN105978003B (en) | 2017-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105978003A (en) | Power system additional wide-area damping controller design method giving consideration to time lag | |
Li et al. | An efficient optimal control method for open-loop transient stability emergency control | |
CN104917174B (en) | Static voltage stability judgment method for active power distribution network | |
CN106374513B (en) | A kind of more microgrid dominant eigenvalues optimization methods based on leader-followers games | |
CN105184027A (en) | Power load modeling method based on interactive multi-model algorithm | |
CN104600697B (en) | Quasi-direct current optimal power flow method considering temperature influence | |
CN105468909A (en) | Time delay power system electromechanical oscillation mode computing method based on SOD-PS-R R algorithm | |
CN106532710A (en) | Microgrid power flow optimization method considering voltage stability constraint | |
CN106655226A (en) | Active power distribution network asymmetric operation optimization method based on intelligent soft open point | |
CN106655227A (en) | SOP-based active power distribution network feeder load balancing method | |
CN103401270A (en) | Method and device for determining accepting ability of medium and lower voltage distribution network to distributed power sources | |
CN106849078A (en) | A kind of power system mesomeric state frequency predication method of meter and dead-time voltage | |
CN103810646A (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
CN104794531A (en) | MDCOPF (modified direct current optimal power flow) method based on grid loss equivalent load model | |
CN105243432A (en) | NSGA-III-based UPFC (unified power flow controller) location and capacity multi-target configuration method | |
CN108054757A (en) | A kind of embedded idle and voltage N-1 Close loop security check methods | |
CN110611305B (en) | Photovoltaic access planning method considering distribution network voltage out-of-limit risk | |
CN101702521B (en) | State estimation method for electric power system considering influences of multi-balancing machine | |
CN106532767A (en) | Small-interfering stability assessment method for power electronic single feed electric power system | |
CN104410065B (en) | A kind of Multiobjective Decision Making Method of receiving end electrical network limiting short-circuit current | |
CN117592257A (en) | Construction method and system of digital twin model of gate control water network | |
CN109494748A (en) | Newton load flow calculation method based on node type and modified Jacobian matrix | |
CN103401271A (en) | Method and device for determining accepting ability of provincial power grid to distributed power sources | |
CN108808706A (en) | Time-lag power system electromechanic oscillation mode computational methods based on SOD-PS-II-R algorithms | |
CN105576653B (en) | A kind of 220kV sections power network power supply capacity optimization 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |