CN103050992A - Analyzing method for damping torque having small interfering stable influence on power system caused by wind power integration - Google Patents

Abstract

The invention discloses an analyzing method for a damping torque having a small interfering stable influence on a power system caused by wind power integration. The method comprises the following steps: performing depreciating treatment on a system linear matrix containing no fan set and a system linear matrix containing an appointed fan set, thereby obtaining a damping torque matrix provided by a generator before and after the fan is added, and then calculating a participating factor of each generator corresponding to an appointed vibration mode, thereby obtaining the influence of the accessing of a wind fan on the appointed vibration mode. The analyzing method can be utilized to clearly and detailedly give out how the appointed fan influences the electromechanical vibration mode of an appointed power system through each synchronous generator, so that a physical significance guidance is supplied to the small interfering vibrating stability analysis for the power system connected with the fan, the selection and the control design for a fan mounting place.

Description

Wind-electricity integration is on the damping torque analytical method of electric power system small interference stability impact

Technical field

A kind of wind power system of the present invention is specifically related to wind-electricity integration to the damping torque analytical method of electric power system small interference stability impact.

Background technology

The advantage that wind-powered electricity generation is pollution-free, renewable with it and the amount of holding is abundant, development is swift and violent in China.Along with wind-powered electricity generation research deeply and the progress of technology, the capacity of blower fan group constantly increases, operation of power networks is built up and dropped into to increasing big-and-middle-sized wind energy turbine set in succession.Because the continuous expansion of scale, wind energy turbine set fluctuates from simple local voltage on the impact of electrical network, the power quality problems such as harmonic pollution develop into electricity net safety stable, the frequency modulation peak regulation, the aspects such as economic dispatch, for example document " research of wind energy turbine set generate output confidence level " (" Proceedings of the CSEE " 10 phases in 2005), " external wind power generation guide rule and dynamic model brief introduction " (" electric power network technique " 12 phases in 2005), " High penetration of wind power inpower systems:an ISO perspective IEEE " and " impact of large Wind Farms on Power System Small Signal Stability and damping characteristic " (" electric power network technique " 13 phases in 2007) all sets forth to some extent.Therefore, the impact on power system safety and stability has important theory significance and practical value after the access of research wind energy turbine set.

China's wind-resources is abundanter, but be fit to the area of large-scale development wind-powered electricity generation generally all away from load center, thus may occur behind the large-scale wind power access electrical network that the line voltage level descends, the transmission of circuit overpower, system short circuit capacity increases and power system dynamic stability such as sexually revises at the series of problems.Document " containing combined iteration method and application that the electric power system tide of wind energy turbine set calculates " (" electric power network technique " 18 phases in 2005) and " determining the optimization method of electric power system wind power penetration limit under a kind of Static Security Constraints " (6 phases of " Proceedings of the CSEE " calendar year 2001) has been studied the static characteristic that the blower fan group accesses system behind the electrical network from aspects such as basic tidal current analysis, static security analysis; Aspect dynamic stability, research work mainly concentrates on the checking work that blower fan is set up mould and model, for example document " Dynamic Model of Wind Turbine Generator Sets research " (" Proceedings of the CSEE " 3 phases in 2005) and " large-sized speed-changing constant-frequency wind turbine generator group modeling and simulation " (" Proceedings of the CSEE " 6 phases in 2004).About the reciprocal effect of blower fan group and electric power system, the actual large electrical network of especially large capacity wind energy turbine set access still lacks systematic research at present on the impact of power system safety and stability operation.

The low-frequency oscillation that occurs in the electric power system derives from relatively waving between generator or the electric power generator group rotor, it is to occur local meritorious residue or shortage in system, relevant generator or electric power generator group therefore and acceleration and deceleration when reaching active balance, because the damping shortage causes the persistent oscillation of the active power that occurs in the system.The frequency range of vibration is generally between 0.2 ~ 2.5Hz, so be called low-frequency oscillation or electromechanical oscillations.In recent years low-frequency oscillation happens occasionally in China, has had a strong impact on power delivery and safe and stable operation between electrical network.Greatly develop new forms of energy especially wind energy be the direction of China Power development, along with the adding of large-scale wind power, it is complicated that power system operation constantly is tending towards, therefore very urgent on the research of Study of Power System Small Disturbance aspect impact for the wind-powered electricity generation access.

Summary of the invention

Goal of the invention: the object of the invention is in order to overcome the deficiencies in the prior art, provide a kind of and can access the access of blower fan to the effect on power system result, and obtain simultaneously more internal mechanism information, the calculating that provides quantitative for choosing of the selection in blower fan access place and controller and a kind of wind-electricity integration of analysis result are on the damping torque analytical method of electric power system small interference stability impact.

Technical scheme: wind-electricity integration of the present invention may further comprise the steps the damping torque analytical method of electric power system small interference stability impact:

(1) obtain the power system mesomeric state data by data acquisition and supervisory control system (SCADA system), EMS (EMS): generator terminal voltage, machine end are gained merit, bus is meritorious and bus is idle;

(2) substitution power system mesomeric state data in the system linearization matrix: generator terminal voltage, the machine end is meritorious, bus is meritorious and bus is idle;

(3) the power system static data that substitution traffic department provides in the system linearization matrix: electric power networks topological data, line impedance admittance data, transformer impedance no-load voltage ratio data;

(4) the generating set inherent data that substitution generating set production firm provides in the system linearization matrix: generator reactance data, excitation system data;

(5) the blower fan group inherent data that substitution blower fan group production firm provides in the system linearization matrix: the inner reactance data of blower fan, revolutional slip data, blower fan merit angular data;

(6) by step (2)-(4), do not contained the system linearization matrix of blower fan group:

$\mathrm{\Δ}{\stackrel{\·}{X}}_{\mathrm{without}}=A_{\mathrm{without}}\mathrm{\Δ}{X}_{\mathrm{without}}$

And:

$\left[\begin{array}{c}\mathrm{\Δ}\stackrel{\·}{\mathrm{\δ}}\\ \mathrm{\Δ}\stackrel{\·}{\mathrm{\ω}}\\ \mathrm{\Δ}\stackrel{\·}{Z}\end{array}\right]=\left[\begin{array}{ccc}0& {\mathrm{\ω}}_{0}I& 0\\ A_{21\mathrm{without}}& A_{22\mathrm{without}}& A_{23\mathrm{without}}\\ A_{31\mathrm{without}}& A_{32\mathrm{without}}& A_{33\mathrm{without}}\end{array}\right]\left[\begin{array}{c}\mathrm{\Δ\δ}\\ \mathrm{\Δ\ω}\\ \mathrm{\ΔZ}\end{array}\right]---\left(1\right)$

Wherein, X _WithoutFor not containing the system state variables of blower fan group, A _WithoutFor not containing the system linearization matrix of blower fan group, δ is generator's power and angle state variable vector, and ω is generator speed state variable vector, and Z is other state variable vector of synchronous motor, and Δ is linearized operator,

With

Be the differential operator of X, δ, ω and Ζ, I is the unit diagonal matrix, ω ₀Be rated angular velocity, $A_{21\mathrm{without}}=\∂\stackrel{\·}{\mathrm{\ω}}/\∂\mathrm{\δ},$ $A_{22\mathrm{without}}=\∂\stackrel{\·}{\mathrm{\ω}}/\∂\mathrm{\ω},$ $A_{23\mathrm{without}}=\∂\stackrel{\·}{\mathrm{\ω}}/\∂Z,$ $A_{31\mathrm{without}}=\∂\stackrel{\·}{Z}/\∂\mathrm{\δ},$ $A_{32\mathrm{without}}=\∂\stackrel{\·}{Z}/\∂\mathrm{\ω},$ $A_{33\mathrm{without}}=\∂\stackrel{\·}{Z}/\∂Z;$

(7) will not contain the system linearization matrix A of blower fan machine _WithoutCarry out depression of order and process, obtain

$\mathrm{\Δ}\stackrel{\·}{\mathrm{\ω}}=-M^{-1}B\left({\mathrm{\λ}}_i\right)\mathrm{\Δ\ω}$

Wherein, M is the diagonal matrix that comprises all synchronous generator inertia constant information, B (λ _i) be the matrix of N*N, its element B _Jk(λ _i) representing the damping torque that j platform generator provides to i Oscillatory mode shape by k platform generator, N is the number of units of generator in the system.λ _i=-ξ _i± j ω _iCharacteristic root for i Oscillatory mode shape of certain electromechanics of designated analysis in the system obtains according to formula (1):

B(λ _i)=-M((A _21without+A _23withoutA _zδwithout)A _δωwithout+(A _22without+A _23without)A _zωwithout) (2)

S=λ wherein _i=-ξ _i± j ω _i, $A_{\mathrm{z\δwithout}}=\∂Z/\∂\mathrm{\δ}={(\mathrm{sI}-A_{33\mathrm{without}})}^{-1}{A}_{31\mathrm{without}},$ $A_{\mathrm{z\ωwithout}}=\∂Z/\∂\mathrm{\ω}={(\mathrm{sI}-A_{33\mathrm{without}})}^{-1}{A}_{32\mathrm{without}},$ $A_{\mathrm{\δ\ωwithout}}=\∂\mathrm{\δ}/\∂\mathrm{\ω}={(\mathrm{sI}-(A_{11\mathrm{without}}+A_{13\mathrm{without}}{A}_{\mathrm{z\δwithout}}))}^{-1}(A_{12\mathrm{without}}+A_{13\mathrm{without}}{A}_{\mathrm{z\ωwithout}});$

Can get according to Fig. 1 and not contain blower fan, the damping torque matrix that provides for the generator of i Oscillatory mode shape is:

$T_{\mathrm{without}}=B\left({\mathrm{\λ}}_i\right)=\left[\begin{array}{ccccc}{B}_{11}\left({\mathrm{\λ}}_i\right)& ...& B_{1k}\left({\mathrm{\λ}}_i\right)& ...& B_{1N}\left({\mathrm{\λ}}_i\right)\\ ...& ...& ...& ...& ...\\ B_{j1}\left({\mathrm{\λ}}_i\right)& ...& B_{\mathrm{jk}}\left({\mathrm{\λ}}_i\right)& ...& B_{\mathrm{jN}}\left({\mathrm{\λ}}_i\right)\\ ...& ...& ...& ...& ...\\ B_{N1}\left({\mathrm{\λ}}_i\right)& ...& B_{\mathrm{Nk}}\left({\mathrm{\λ}}_i\right)& ...& B_{\mathrm{NN}}\left({\mathrm{\λ}}_i\right)\end{array}\right]---\left(3\right)$

(8) equally by step (2) ~ (5), obtain containing the system linearization matrix of specifying the blower fan group: $\mathrm{\Δ}{\stackrel{\·}{X}}_{\mathrm{with}}=A_{\mathrm{with}}\mathrm{\Δ}{X}_{\mathrm{with}}$

And:

$\left[\begin{array}{c}\mathrm{\&Delta;}\stackrel{\&CenterDot;}{\mathrm{\&delta;}}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{\mathrm{\&omega;}}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{Z}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{W}\end{array}\right]=\left[\begin{array}{cccc}0& {\mathrm{\&omega;}}_0\mathrm{<figure-callout\; id="0"\; label="I"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">I</figure-callout>}& 0& 0\\ A_{21\mathrm{with}}& A_{22\mathrm{with}}& A_{23\mathrm{with}}& A_{24\mathrm{with}}\\ A_{31\mathrm{with}}& A_{32\mathrm{with}}& A_{33\mathrm{with}}& A_{34\mathrm{with}}\\ A_{41\mathrm{with}}& A_{42\mathrm{with}}& A_{43\mathrm{with}}& A_{44\mathrm{with}}\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;\&delta;}\\ \mathrm{\&Delta;\&omega;}\\ \mathrm{\&Delta;Z}\\ \mathrm{\&Delta;W}\end{array}\right]---\left(4\right)$

Wherein, X _WithFor containing the system state variables of blower fan group, A _WithFor containing the system linearization matrix of blower fan group, W is the fan condition variable vector,

Be the differential operator of W, $A_{21\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;\mathrm{\&delta;},$ $A_{22\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;\mathrm{\&omega;},$ $A_{23\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;Z,$ $A_{24\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;W,$ $A_{31\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;\mathrm{\&delta;},$ $A_{32\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;\mathrm{\&omega;},$ $A_{33\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;Z,$ $A_{34\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;W,$ $A_{41\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;\mathrm{\&delta;},$ $A_{42\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;\mathrm{\&omega;},$ $A_{43\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;Z,$ $A_{44\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;W;$ Formula (4) can be as shown in Figure 2;

(9) will contain the system linearization matrix A of blower fan _WithCarry out depression of order and process, obtain the damping torque Matrix C (λ that generator provides _i), obtain according to formula (4):

C(λ _i)=-M((A _21with+(A _23with+A _24withA _wzwith)A _zδwith+A _24withA _wδwith)A _δωwith+(A _22with+(A _23with+A _24withA _wzwith)A _zωwith+A _24withA _wωwith) (5)

Wherein, s=λ _i=-ξ _i± j ω _i, $A_{\mathrm{w\&delta;with}}=\&PartialD;W/\&PartialD;\mathrm{\&delta;}={(\mathrm{sI}-A_{44\mathrm{with}})}^{-1}{A}_{41\mathrm{with}},$ $A_{\mathrm{w\&omega;with}}=\&PartialD;W/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-A_{44\mathrm{with}})}^{-1}{A}_{42\mathrm{with}},$ $A_{\mathrm{wzwith}}=\&PartialD;W/\&PartialD;Z={(\mathrm{sI}-A_{44\mathrm{with}})}^{-1}{A}_{43\mathrm{with}},$ $A_{\mathrm{z\&delta;with}}=\mathrm{\&delta;Z}/\&PartialD;\mathrm{\&delta;}={(\mathrm{sI}-A_{33\mathrm{with}}-A_{34\mathrm{with}}{A}_{\mathrm{wzwith}})}^{-1}(A_{31\mathrm{with}}+A_{34\mathrm{with}}{A}_{\mathrm{w\&delta;with}}),$ $A_{\mathrm{z\&omega;with}}=\mathrm{\&delta;Z}/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-A_{33\mathrm{with}}-A_{34\mathrm{with}}{A}_{\mathrm{wzwith}})}^{-1}(A_{32\mathrm{with}}+A_{34\mathrm{with}}{A}_{\mathrm{w\&omega;with}}),$ $A_{\mathrm{\&delta;\&omega;with}}=\&PartialD;\mathrm{\&delta;}/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-(A_{11\mathrm{with}}+(A_{13\mathrm{with}}+A_{14\mathrm{with}}{A}_{\mathrm{wzwith}})A_{\mathrm{z\&delta;with}}+A_{14\mathrm{with}}{A}_{\mathrm{w\&delta;with}}))}^{-1}$ $(A_{12\mathrm{with}}+(A_{31\mathrm{with}}+A_{14\mathrm{with}}{A}_{\mathrm{wzwith}})A_{\mathrm{z\&omega;with}}+A_{14\mathrm{with}}{A}_{\mathrm{w\&omega;with}});$

Can comprise blower fan according to Fig. 2, the damping torque matrix that provides for the generator of i Oscillatory mode shape is:

$T_{\mathrm{with}}=C\left({\mathrm{\&lambda;}}_i\right)=\left[\begin{array}{ccccc}{C}_{11}\left({\mathrm{\&lambda;}}_i\right)& ...& C_{1k}\left({\mathrm{\&lambda;}}_i\right)& ...& C_{1N}\left({\mathrm{\&lambda;}}_i\right)\\ ...& ...& ...& ...& ...\\ C_{j1}\left({\mathrm{\&lambda;}}_i\right)& ...& C_{\mathrm{jk}}\left({\mathrm{\&lambda;}}_i\right)& ...& C_{\mathrm{jN}}\left({\mathrm{\&lambda;}}_i\right)\\ ...& ...& ...& ...& ...\\ C_{N1}\left({\mathrm{\&lambda;}}_i\right)& ...& C_{\mathrm{Nk}}\left({\mathrm{\&lambda;}}_i\right)& ...& C_{\mathrm{NN}}\left({\mathrm{\&lambda;}}_i\right)\end{array}\right]---\left(6\right)$

(10) can obtain according to formula (5) and formula (6), behind the adding blower fan, the variable quantity matrix D (λ of the damping torque that generator provides _i).

$D\left({\mathrm{\&lambda;}}_i\right)=C\left({\mathrm{\&lambda;}}_i\right)-B\left({\mathrm{\&lambda;}}_i\right)=\left[\begin{array}{ccccc}{B}_{11}\left({\mathrm{\&lambda;}}_i\right)-C_{11}\left({\mathrm{\&lambda;}}_i\right)& ...& B_{1k}\left({\mathrm{\&lambda;}}_i\right)-C_{1k}\left({\mathrm{\&lambda;}}_i\right)& ...& B_{1N}\left({\mathrm{\&lambda;}}_i\right)-C_{1N}\left({\mathrm{\&lambda;}}_i\right)\\ ...& ...& ...& ...& ...\\ B_{j1}\left({\mathrm{\&lambda;}}_i\right)-C_{j1}\left({\mathrm{\&lambda;}}_i\right)& ...& B_{\mathrm{jk}}\left({\mathrm{\&lambda;}}_i\right)-C_{\mathrm{jk}}\left({\mathrm{\&lambda;}}_i\right)& ...& B_{\mathrm{jN}}\left({\mathrm{\&lambda;}}_i\right)-C_{\mathrm{jN}}\left({\mathrm{\&lambda;}}_i\right)\\ ...& ...& ...& ...& ...\\ B_{N1}\left({\mathrm{\&lambda;}}_i\right)-C_{N1}\left({\mathrm{\&lambda;}}_i\right)& ...& B_{\mathrm{Nk}}\left({\mathrm{\&lambda;}}_i\right)-C_{\mathrm{Nk}}\left({\mathrm{\&lambda;}}_i\right)& ...& B_{\mathrm{NN}}\left({\mathrm{\&lambda;}}_i\right)-C_{\mathrm{NN}}\left({\mathrm{\&lambda;}}_i\right)\end{array}\right]$

$=\left[\begin{array}{ccccc}{D}_{11}\left({\mathrm{\&lambda;}}_i\right)& ...& D_{1k}\left({\mathrm{\&lambda;}}_i\right)& ...& D_{1N}\left({\mathrm{\&lambda;}}_i\right)\\ ...& ...& ...& ...& ...\\ D_{j1}\left({\mathrm{\&lambda;}}_i\right)& ...& D_{\mathrm{jk}}\left({\mathrm{\&lambda;}}_i\right)& ...& D_{\mathrm{jN}}\left({\mathrm{\&lambda;}}_i\right)\\ ...& ...& ...& ...& ...\\ D_{N1}\left({\mathrm{\&lambda;}}_i\right)& ...& D_{\mathrm{Nk}}\left({\mathrm{\&lambda;}}_i\right)& ...& D_{\mathrm{NN}}\left({\mathrm{\&lambda;}}_i\right)\end{array}\right]---\left(7\right)$

(11) by formula (7) is reconstructed, off-diagonal element is put in order to diagonal element, formed new diagonal matrix D _New(λ _i):

$D_{\mathrm{new}}\left({\mathrm{\&lambda;}}_i\right)=\left[\begin{array}{ccccc}{D}_{11\mathrm{new}}\left({\mathrm{\&lambda;}}_i\right)& ...& 0& ...& 0\\ ...& ...& ...& ...& ...\\ 0& ...& D_{\mathrm{jjnew}}\left({\mathrm{\&lambda;}}_i\right)& ...& 0\\ ...& ...& ...& ...& ...\\ 0& ...& 0& ...& D_{\mathrm{NNnew}}\left({\mathrm{\&lambda;}}_i\right)\end{array}\right]---\left(8\right)$

Wherein, diagonal element D _Newj(λ _i), j=1,2 ... the damping torque that the blower fan that N represents to access provides by each generator;

(12) in the computing system each generator corresponding to the participation factors S that specifies Oscillatory mode shape _Ij:

$S_{\mathrm{ij}}=\frac{\&PartialD;{\mathrm{\&lambda;}}_i}{\&PartialD;\mathrm{\&Delta;}{D}_j}=w_i^T\frac{\&PartialD;}{\&PartialD;\mathrm{\&Delta;}{D}_j}(-M^{-1}(C\left({\mathrm{\&lambda;}}_i\right)+\mathrm{\&Delta;}{D}_j))\&times;v_i=-w_{\mathrm{ij}}\frac{1}{M_j}{v}_{\mathrm{ij}},j=\mathrm{1,2},...N---\left(9\right)$

Wherein, v _iFor corresponding to λ _iRight characteristic vector, v _IjBe v _iIn corresponding to △ ω _j, j=1,2 ... the component of N, j are j platform generator, w _iFor corresponding to λ _iLeft eigenvector, w _IjBe w _iIn corresponding to △ ω _j, j=1,2 ... the component of N;

Following formula shows: each generator is the sensitivity coefficient S that specifies the damping torque that Oscillatory mode shape provides to each generator blower fan in the computing system of formula (9) definition corresponding to the participation factors of specifying Oscillatory mode shape _Ij

(13) calculating of through type (8) and formula (9) can obtain the blower fan access to specifying the impact of Oscillatory mode shape:

$\mathrm{\&Delta;}{\mathrm{\&lambda;}}_i=\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}{S}_{\mathrm{ij}}{D}_{\mathrm{newj}}\left({\mathrm{\&lambda;}}_i\right)---\left(10\right)$

Wherein, △ λ _iFor specifying this variable of Oscillatory mode shape, can be used for arranging to specify the installation place of blower fan and choosing of controller;

(14) can oppositely prove the reasonability of formula (10) by following proof procedure:

With A _WithoutAnd A _WithAt keeping characteristics root λ _iPrerequisite under carry out depression of order and process, obtain damping torque matrix B (λ _i) and C (λ _i), therefore can obtain:

$w_i^T(-M^{-1}B)v_i={\mathrm{\&lambda;}}_{\mathrm{withouti}}---\left(11\right)$

$w_i^T(-M^{-1}C)v_i={\mathrm{\&lambda;}}_{\mathrm{withi}}---\left(12\right)$

Subtract each other and can get by formula (11) and formula (12), the variation of characteristic root behind the adding blower fan:

${\mathrm{\&lambda;}}_{\mathrm{withi}}-{\mathrm{\&lambda;}}_{\mathrm{withouti}}=\mathrm{\&Delta;}{\mathrm{\&lambda;}}_i=-w_i^T{M}^{-1}(C-B)v_i=-w_i^T{M}^{-1}D{v}_i$

$=\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}-w_{\mathrm{ij}}{D}_{\mathrm{newj}}{v}_{\mathrm{ij}}/M_j=\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}{S}_{\mathrm{ij}}{D}_{\mathrm{newj}}---\left(13\right)$

Wherein, λ _WithoutiBe the characteristic root of i Oscillatory mode shape of certain generator of not containing the blower fan group, λ _WithiCharacteristic root for i Oscillatory mode shape of certain generator of containing the blower fan group; By formula (10) as can be known: the diagonal matrix D of reconstruct _New(λ _i) in diagonal element measured the damping torque that blower fan provides to each generator, and participation factors S _IjThe damping torque that provides has been provided is and how to be converted into the impact of specifying electromechanical oscillations mode; So the physical significance of model analysis formula (10) is: specify blower fan to provide damping torque to every generator, its size is by diagonal matrix D _New(λ _i) in diagonal element D _Jjnew(λ _i), j=1,2 ... N tolerance; Multiply by again participation factors S _IjAfter, the damping torque of every generator acquisition just is converted into specifies blower fan by the impact of each generator on appointment electromechanical oscillations mode, and its size is by D _Newj(λ _i) S _IjTolerance, influence index △ λ _iBe exactly that blower fan passes through N platform generator to specifying the impact of electromechanical oscillations mode, be N item D _Jjnew(λ _i) S _Ij, j=1,2 ... the N sum; The physical significance of this model analysis formula specifies blower fan at first to every generator G1, G2 as shown in Figure 3 ... GN provides damping torque tolerance D _Newj(λ _i), j=1,2 ... N, then every generator G1, G2 ... GN passes through participation factors S with the damping torque that obtains _Ij, j=1,2 ... N transforms specifying the tolerance of motor oscillates Effect of Mode.

Beneficial effect: wind-electricity integration of the present invention is on the damping torque analytical method of electric power system small interference stability impact, can clearly at length provide and specify blower fan is the electric power system electromechanical oscillations mode that how affect appointment by each synchronous generator, thereby provides the physical significance guidance for the electric power system small interference oscillatory stability analysis of access blower fan, selection and the control design in assembling place.

Description of drawings

Fig. 1 when not containing blower fan, system linearity state equation schematic diagram;

Fig. 2 when comprising blower fan, system linearity state equation schematic diagram;

Fig. 3 is the physical significance schematic diagram of damping torque analytical method model analysis formula;

Fig. 4 is the system schematic that three machines, nine nodal analysis methods add the double-fed blower fan;

Fig. 5 is double-fed air-blower control model schematic diagram;

Fig. 6 is FJ type excitation system model schematic diagram;

Fig. 7 is the physical significance schematic diagram of damping torque analytical method model analysis breakdown;

Fig. 8 is a large grid wiring schematic diagram of actual measurement;

Fig. 9 is the physical significance schematic diagram of damping torque analytical method model analysis breakdown;

Figure 10 is that blower fan passes through each unit to the change amount result schematic diagram of modal characteristics root.

Embodiment

The below is elaborated to technical solution of the present invention, but protection scope of the present invention is not limited to described embodiment.

Embodiment: as shown in Figure 4, being the electric power system that three machines, nine nodal analysis methods add a double-fed blower fan, is incorporated into the power networks on the specific implementation process of the damping torque analytical method of power system small signal stability impact in the large-scale wind power field that proposes by electric power system shown in Figure 4 explanation the present invention; The object of model analysis is the appointment blower fan, and the content of model analysis is that the blower fan access is on specifying the impact of Oscillatory mode shape.

The parameter of the electric power system of three machines, nine nodal analysis methods is:

Synchronous generator:

Moment of inertia (p.u.): H ₁=23.6, H ₂=6.4, H ₃=3.01

Motor damping moment coefficient (p.u.): D ₁=D ₂=D ₃=0

D-axis transient state open circuit time constant (s): T _D01'=8.96s, T _D02'=6.0s, T _D03'=5.89s

The unsaturated synchronous reactance of d-axis (p.u.): x _D1=0.1460, x _D2=0.8958, x _D3=1.3130

Hand over the unsaturated synchronous reactance of axle (p.u.): x _Q1=0.0969, x _Q2=0.8645, x _Q3=1.2580

D-axis transient state reactance (p.u.): x' _D1=0.0608, x' _D2=0.1189, x' _D3=0.1813

Excitation system:

Voltage regulator gain (p.u.): K _A1=100, K _A2=100, K _A3=100

Exciter control system stabilizer gain (p.u.): K _F1=0.05, K _F2=0.05, K _F3=0.05

Voltage regulator time constant (s): T _A1=0.05s, T _A2=0.05s, T _A3=0.05s

T _B1=0.03s,T _B2=0.03s,T _B3=0.03s

T _C1=0.02s,T _C2=0.02s,T _C3=0.02s

Exciter control system stabilizer time constant (s): T _F1=0.7s, T _F2=0.7s, T _F3=0.7s

The double-fed blower fan:

Stator reactance (p.u.): x _s=2.9; Excitatory reactance (p.u.): x _m=2.6; Rotor reactance (p.u.): x _r=2.9;

Rotor resistance (p.u.): r _r=0.0013; Inertia constant (p.u.): T _j=3.4; Damping coefficient (p.u.): D=0;

Revolutional slip (p.u.): slip=-0.1

Transmission line reactance (p.u.): x _AB=0.246, x _AC=0.262, x _BC=0.1728, x _A1=x _B2=x _C3=0.01

Generator active power of output (100MW): P _Gen2=1.629; P _Gen3=0.85; G1 is the balancing machine of native system.

DFIG active power of output: P _DFIG=1.0

Load condition in the system: load active power (100MW): P _LoadA=2.25; P _LoadB=1.25; P _LoadC=1.00

Reactive load power (100MVA): Q _LoadA=0.70; Q _LoadB=0.50; Q _LoadC=0.55

1. the inearized model of double-fed blower fan

According to the document (S.Q.Bu that has delivered, W.Du, Z.Chen, L.Y.Xiao and H.F.Li, " ANovel ControlStrategy for the Doubly Fed Induction Generators to Improve Grid Fault Ride-ThroughCapacity ", 2 ^NdIEEE PES International Conference and Exhibition on Innovative Smart GridTechnologies (ISGT Europe) 2011, pp 1-7, Dec.2011), double-fed air-blower control model is as shown in Figure 5.

According to the document (S.Q.Bu that has delivered, W.Du, Z.Chen, L.Y.Xiao and H.F.Li, " A Novel ControlStrategy for the Doubly Fed Induction Generators to Improve Grid Fault Ride-ThroughCapacity ", 2 ^NdIEEE PES International Conference and Exhibition on Innovative Smart GridTechnologies (ISGT Europe) 2011, pp 1-7, Dec.2011; S.Q.Bu, W.Du, H.F.Wang, Z.Chen, L.Y.Xiao and H.F.Li, " Small-signal probabilistic stability of power systems considering thestochastic uncertainty of Grid-connected wind farm ", IET Conference on Renewable PowerGeneration (RPG 2011), pp 1-6, September 2011; F.Mei and B.Pal, " Modeling andSmall-signal Analysis of a Grid-connected Doubly-Fed Induction Generator ", IEEE PowerEngineering Society General Meeting 2005, pp 2101-2108, June 2005), the inearized model of double-fed blower fan can get, as the formula (14):

$\left\{\begin{array}{c}\frac{\mathrm{d\&Delta;}{X}_{\mathrm{DFIG}}}{\mathrm{dt}}=A_{\mathrm{DFIG}}\mathrm{\&Delta;}{X}_{\mathrm{DFIG}}+B_{\mathrm{DFIG}}\mathrm{\&Delta;}{V}_{\mathrm{DFIG}}\\ \mathrm{\&Delta;}{I}_{\mathrm{DFIG}}=C_{\mathrm{DFIG}}\mathrm{\&Delta;}{X}_{\mathrm{DFIG}}+D_{\mathrm{DFIG}}\mathrm{\&Delta;}{V}_{\mathrm{DFIG}}\end{array}\right.---\left(14\right)$

Wherein, Δ X _DFIG=[X _Dfig1X _Dfig2...] ^T, Δ I _DFIG=[I _Dfig1xI _Dfig1yI _Dfig2xI _Dfig2y...] ^T, Δ V _DFIG=[V _Dfig1xV _Dfig1yV _Dfig2xV _Dfig2y...] ^T, A _DFIG, B _DFIG, C _DFIGAnd D _DFIGBe diagonal matrix, the linearisation matrix of the element representation separate unit double-fed blower fan on its diagonal.

2. the inearized model of synchronous generator

Every synchronous generator has all been installed FJ type excitation in native system, as shown in Figure 6.According to document (Wang Xifan chief editor's " modern power systems analysis " and Ni Yixin, Chen Shousun, " theoritical analysis of dynamic power system " of Zhang Baolin work), the inearized model of synchronous generator can get, as the formula (15):

$\left\{\begin{array}{c}\frac{\mathrm{d\&Delta;}{X}_{\mathrm{Gen}}}{\mathrm{dt}}=A_{\mathrm{Gen}}\mathrm{\&Delta;}{X}_{\mathrm{Gen}}+B_{\mathrm{Gen}}\mathrm{\&Delta;}{V}_{\mathrm{Gen}}\\ \mathrm{\&Delta;}{I}_{\mathrm{Gen}}=C_{\mathrm{Gen}}\mathrm{\&Delta;}{X}_{\mathrm{Gen}}+D_{\mathrm{Gen}}\mathrm{\&Delta;}{V}_{\mathrm{Gen}}\end{array}\right.---\left(15\right)$

Wherein, Δ X _Gen=[X _Gen1X _Gen2...] ^T, Δ I _Gen=[I _Gen1xI _Gen1yI _Gen2xI _Gen2y...] ^T, Δ V _Gen=[V _Gen1xV _Gen1yV _Gen2xV _Gen2y...] ^T, A _Gen, B _Gen, C _GenAnd D _GenBe diagonal matrix, the linearisation matrix of the element representation separate unit synchronous generator on its diagonal.

3. system-wide inearized model

Can get system-wide inearized model in conjunction with the inearized model (formula (14)) of double-fed blower fan and the inearized model (formula (15)) of synchronous generator, as the formula (16):

$\left\{\begin{array}{c}\frac{\mathrm{d\&Delta;}{X}_{\mathrm{system}}}{\mathrm{dt}}=A_{\mathrm{system}}\mathrm{\&Delta;}{X}_{\mathrm{system}}+B_{\mathrm{system}}\mathrm{\&Delta;}{\mathrm{<figure-callout\; id="0"\; label="V\; system"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">V</figure-callout>}}_{\mathrm{system}}\\ 0=C_{\mathrm{system}}\mathrm{\&Delta;}{X}_{\mathrm{system}}+D_{\mathrm{system}}\mathrm{\&Delta;}{V}_{\mathrm{system}}\end{array}\right.---\left(16\right)$

Wherein, △ X _System=[△ X _GenΔ X _DFIG] ^T, △ V _System=[△ V _GenΔ V _DFIGΔ V _Rest] ^T, Δ V _RestExpression other node voltage vectors except synchronous generator and double-fed blower fan.

$A_{\mathrm{system}}=\left[\begin{array}{cc}{A}_{\mathrm{<figure-callout\; id="0"\; label="Gen"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">Gen</figure-callout>}}& 0\\ 0& A_{\mathrm{DFIG}}\end{array}\right]$

$B_{\mathrm{system}}=\left[\begin{array}{ccc}{\mathrm{<figure-callout\; id="0"\; label="B\; Gen"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">B</figure-callout>}}_{\mathrm{Gen}}& 0& 0\\ 0& {\mathrm{<figure-callout\; id="0"\; label="B\; DFIG"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">B</figure-callout>}}_{\mathrm{DFIG}}& 0\end{array}\right]$

$C_{\mathrm{system}}=\left[\begin{array}{cc}-{\mathrm{<figure-callout\; id="0"\; label="C\; Gen"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">C</figure-callout>}}_{\mathrm{Gen}}& 0\\ 0& -{\mathrm{<figure-callout\; id="0"\; label="C\; DFIG"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">C</figure-callout>}}_{\mathrm{DFIG}}\\ 0& 0\end{array}\right]$

$D_{\mathrm{system}}=\left[\begin{array}{ccc}{Y}_{\mathrm{GG}}-D_G& Y_{\mathrm{GD}}& Y_{\mathrm{GR}}\\ Y_{\mathrm{DG}}& Y_{\mathrm{DD}}-D_D& Y_{\mathrm{DR}}\\ Y_{\mathrm{RG}}& Y_{\mathrm{RD}}& Y_{\mathrm{RR}}\end{array}\right]$

Subscript G, D and R represent respectively synchronous generator connected node, blower fan connected node and residue node, and Y represents network admittance matrix.

4. modal analysis result and method

After the blower fan access, utilize formula (16) can make up system-wide linear condition equation, system has two electromechanic oscillation modes (0.4951+9.9663i) and (0.3358+7.4017i); Before blower fan does not access, at bus A place load is adjusted accordingly, do not changed with the trend that guarantees whole system.By the linearisation matrix as can be known, corresponding electromechanic oscillation mode be (0.4950+9.9660i) with (0.3385+7.3938i); Drawn by modal analysis method, before and after the blower fan access, the variable quantity of Oscillatory mode shape is-0.0001+0.0003i and 0.0027+0.0079i.

For second Oscillatory mode shape, utilize the rear impact on system of method research blower fan access of the present invention:

Before the blower fan access, the damping torque matrix B (λ that generator provides _i) be:

$B\left({\mathrm{\&lambda;}}_i\right)=\left[\begin{array}{ccc}-1.8864-107.16i& 0.37073+13.402i& 0.13774+3.6031i\\ -6.2103+200.79i& 5.6666-94.825i& -0.045335+11.080i\\ -8.4789+282.38i& -4.1318+58.935i& 5.7335-71.541i\end{array}\right]$

After the blower fan access, the damping torque Matrix C (λ that generator provides _i) be:

$C\left({\mathrm{\&lambda;}}_i\right)=\left[\begin{array}{ccc}-2.5485-109.12i& 0.26819+13.109i& 0.11612+3.5548i\\ -7.8113+196.45i& 5.4758-95.284i& -0.10052+10.949i\\ -10.152+277.24i& -4.3592+58.209i& 5.7029-71.574i\end{array}\right]$

After the blower fan access, the variable quantity matrix D (λ of the damping torque that generator provides _i) be:

$D\left({\mathrm{\&lambda;}}_i\right)=C\left({\mathrm{\&lambda;}}_i\right)-B\left({\mathrm{\&lambda;}}_i\right)=\left[\begin{array}{ccc}-0.6639-1.9600i& -0.1025-0.2930i& -0.0216-0.0483i\\ -1.6010-4.3400i& -0.1908-0.4590i& -0.0552-0.1310i\\ -1.6731-5.1400i& -0.2274-0.7260i& -0.0306-0.0330i\end{array}\right]$

With D (λ _i) diagonalization, the damping torque matrix D that the blower fan of access provides by each generator _New(λ _i) be:

$D_{\mathrm{new}}\left({\mathrm{\&lambda;}}_i\right)=\left[\begin{array}{ccc}-0.1290-0.3900\mathrm{<figure-callout\; id="0"\; label="i"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">i</figure-callout>}& 0& 0\\ 0& -0.0351-0.1010\mathrm{<figure-callout\; id="0"\; label="i"\; filenames="HDA00002557400400011.png,HDA00002557400400051.png"\; state="\{\{state\}\}">i</figure-callout>}& 0\\ 0& 0& -0.0165-0.0470i\end{array}\right]$

Each unit damping torque to the sensitivity of specifying Oscillatory mode shape is:

S _ij=[-0.0044+0.0005i -0.0530-0.0020i -0.0190+0.0004i] ^T

After the blower fan access, by each generator to being changed to that the oscillation mode characteristic root of research produces:

△λ _ij=[-0.0125-0.0322i 0.0127+0.0344i 0.0024+0.0057i] ^T

$\mathrm{\&Delta;}{\mathrm{\&lambda;}}_i=\underset{j=1}{\overset{3}{\mathrm{\&Sigma;}}}\mathrm{\&Delta;}{\mathrm{\&lambda;}}_{\mathrm{ij}-j}=0.0026+0.0079i$

The variable quantity of the described damping torque analytical method of present embodiment gained modal characteristics root is consistent with the model analysis acquired results, the physical interpretation that has provided modal analysis result based on the damping torque analytical method shown in Figure 7.

The damping torque analytical method that the present invention proposes, be easily understood, physical significance is clear, it is based upon damping torque that generator amature motion obtains, and this is notional, when being generalized to research blower fan access complicated multi-machine power system affected, the mechanism that can affect existing Oscillatory mode shape from the access of having analyzed in essence blower fan, physical significance is clear, can find the basic reason that causes existing Mode variation.

(1) distribution of damping torque and transmission: pass through the method, can answer the problem of " why ", namely by after observing the blower fan access, the characteristic root that each generator provides changes, and can judge clearly the variation of Oscillatory mode shape mainly by which platform unit is provided.The variation of blower fan Oscillatory mode shape that access causes is how to pass through each generator distribution and transmission.

(2) choosing of blower fan access place: blower fan can be accessed in the different location, drawing in difference access place by this analytical method is impact on system oscillation mode.In situation about can select, choose the less place access blower fan of impact; In the time can not selecting, should consider during access blower fan is taked control strategy.

The mode breakdown analytical method of using the present invention to propose is tested in China large system of a certain reality, the simplification winding diagram of this large system of reality as shown in Figure 8, test to as if system in the blower fan that accesses, purpose is to analyze the blower fan access to the impact of the appointment low-frequency oscillation mode in the system (code is SB).

SB pattern main manifestations be (Lyg+ zone, Ha+ zone, Sq+ zone, regional Xz+ zone Yc) to the power oscillation of (Cz+ zone, Zj+ zone, regional Nj+ zone Wx+ zone Sz), 13 equivalent blower fans (all blower fans that are installed in same wind energy turbine set will be a large capacity and equivalent blower fan with same dynamic characteristic by equivalence at this) connecting system will be arranged in this system.

When not accessing blower fan, the characteristic root that be can be calculated the SB pattern by modal analysis method is :-0.1588+4.9440i, behind 13 equivalent blower fan connecting systems, the characteristic root that can get the SB pattern is-0.1620+4.9577i, after adding blower fan, SB pattern feature root variable quantity is-0.0032+0.0137i.

According to the damping torque analysis theories, at first obtain before and after the adding blower fan unit damping torque contribute matrix B (λ _i), C (λ _i), obtain subsequently before and after the adding blower fan transformation matrices D (λ of unit damping contribution _i)=C (λ _i)-B (λ _i), its diagonalization is got blower fan by each unit damping contribution amount D _New(λ _i).Blower fan be multiply by each unit damping contribution by each unit damping contribution amount can obtain blower fan by the change amount of each unit to appointment modal characteristics root to the sensitivity of specifying Oscillatory mode shape, such as Fig. 9, shown in Figure 10, the variable quantity of SB pattern feature root is consistent with the result that modal analysis method obtains before and after the blower fan access that obtains with method shown in the present.

From the damping torque analytical method, the access that can see clearly blowing machine is how to distribute at the whole network on the impact of SB modal damping, thereby transmission forms, for further providing theoretical foundation to the air-blower control of blower fan connecting system and the research of synchronous motor control.

As mentioned above, although the specific preferred embodiment of reference has represented and has explained the present invention that it shall not be construed as the restriction to the present invention self.Under the spirit and scope of the present invention prerequisite that does not break away from the claims definition, can make in the form and details various variations to it.

Claims (2)

1. wind-electricity integration is characterized in that: comprise the steps: the damping torque analytical method of electric power system small interference stability impact

(1) obtain the power system mesomeric state data by data acquisition and supervisory control system (SCADA system), EMS (EMS): generator terminal voltage, machine end are gained merit, bus is meritorious and bus is idle;

(2) the described power system mesomeric state data of substitution in the system linearization matrix;

(3) the power system static data that substitution traffic department provides in the system linearization matrix: electric power networks topological data, line impedance admittance data, transformer impedance no-load voltage ratio data;

(4) the generating set inherent data that substitution generating set production firm provides in the system linearization matrix: generator reactance data, excitation system data;

(5) the blower fan group inherent data that substitution blower fan group production firm provides in the system linearization matrix: the inner reactance data of blower fan, revolutional slip data, blower fan merit angular data;

(6) by step (2) ~ (4), do not contained the system linearization matrix of blower fan group:

$\mathrm{\&Delta;}{\stackrel{\&CenterDot;}{X}}_{\mathrm{without}}=A_{\mathrm{without}}\mathrm{\&Delta;}{X}_{\mathrm{without}}$

And:

$\left[\begin{array}{c}\mathrm{\&Delta;}\stackrel{\&CenterDot;}{\mathrm{\&delta;}}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{\mathrm{\&omega;}}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{Z}\end{array}\right]=\left[\begin{array}{ccc}0& {\mathrm{\&omega;}}_{0}I& 0\\ A_{21\mathrm{without}}& A_{22\mathrm{without}}& A_{23\mathrm{without}}\\ A_{31\mathrm{without}}& A_{32\mathrm{without}}& A_{33\mathrm{without}}\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;\&delta;}\\ \mathrm{\&Delta;\&omega;}\\ \mathrm{\&Delta;Z}\end{array}\right]$

Wherein, X _WithoutFor not containing the system state variables of blower fan group, A _WithoutFor not containing the system linearization matrix of blower fan group, δ is generator's power and angle state variable vector, and ω is generator speed state variable vector, and Z is other state variable vector of synchronous motor, and Δ is linearized operator,

With

Be the differential operator of X, δ, ω and Ζ, I is the unit diagonal matrix, ω ₀Be rated angular velocity, $A_{21\mathrm{without}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;\mathrm{\&delta;},$ $A_{22\mathrm{without}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;\mathrm{\&omega;},$ $A_{23\mathrm{without}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;Z,$ $A_{31\mathrm{without}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;\mathrm{\&delta;},$ $A_{32\mathrm{without}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;\mathrm{\&omega;},$ $A_{33\mathrm{without}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;Z;$

(7) will not contain the system linearization matrix A of blower fan group _WithoutCarry out depression of order and process, obtain:

$\mathrm{\&Delta;}\stackrel{\&CenterDot;}{\mathrm{\&omega;}}=-M^{-1}B\left({\mathrm{\&lambda;}}_i\right)\mathrm{\&Delta;\&omega;}$

Wherein, M is the diagonal matrix that comprises all synchronous generator inertia constant information, B (λ _i) be the matrix of N*N, its element B _Jk(λ _i) representing the damping torque that j platform generator provides to i Oscillatory mode shape by k platform generator, N is the number of units of generator in the system, λ _i=-ξ _i± j ω _iCharacteristic root for i Oscillatory mode shape of certain electromechanics of designated analysis in the system obtains according to step (6):

B(λ _i)=-M((A _21without+A _23withoutA _zδwithout)A _δωwithout+(A _22without+A _23without)A _zωwithout)

S=λ wherein _i=-ξ _i± j ω _i, $A_{\mathrm{z\&delta;without}}=\&PartialD;Z/\&PartialD;\mathrm{\&delta;}={(\mathrm{sI}-A_{33\mathrm{without}})}^{-1}{A}_{31\mathrm{without}},$ $A_{\mathrm{z\&omega;without}}=\&PartialD;Z/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-A_{33\mathrm{without}})}^{-1}{A}_{32\mathrm{without}},$ $A_{\mathrm{\&delta;\&omega;without}}=\&PartialD;\mathrm{\&delta;}/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-(A_{11\mathrm{without}}+A_{13\mathrm{without}}{A}_{\mathrm{z\&delta;without}}))}^{-1}(A_{12\mathrm{without}}+A_{13\mathrm{without}}{A}_{\mathrm{z\&omega;without}});$

(8) same, by step (2) ~ (5), obtain containing the system linearization matrix of specifying the blower fan group:

$\mathrm{\&Delta;}{\stackrel{\&CenterDot;}{X}}_{\mathrm{with}}=A_{\mathrm{with}}\mathrm{\&Delta;}{X}_{\mathrm{with}}$

And:

$\left[\begin{array}{c}\mathrm{\&Delta;}\stackrel{\&CenterDot;}{\mathrm{\&delta;}}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{\mathrm{\&omega;}}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{Z}\\ \mathrm{\&Delta;}\stackrel{\&CenterDot;}{W}\end{array}\right]=\left[\begin{array}{cccc}0& {\mathrm{\&omega;}}_{0}I& 0& 0\\ A_{21\mathrm{with}}& A_{22\mathrm{with}}& A_{23\mathrm{with}}& A_{24\mathrm{with}}\\ A_{31\mathrm{with}}& A_{32\mathrm{with}}& A_{33\mathrm{with}}& A_{34\mathrm{with}}\\ A_{41\mathrm{with}}& A_{42\mathrm{with}}& A_{43\mathrm{with}}& A_{44\mathrm{with}}\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;\&delta;}\\ \mathrm{\&Delta;\&omega;}\\ \mathrm{\&Delta;Z}\\ \mathrm{\&Delta;W}\end{array}\right]$

Wherein, X _WithFor containing the system state variables of blower fan group, A _WithFor containing the system linearization matrix of blower fan group, W is the fan condition variable vector,

Be the differential operator of W, $A_{21\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;\mathrm{\&delta;},$ $A_{22\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;\mathrm{\&omega;},$ $A_{23\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;Z,$ $A_{24\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{\mathrm{\&omega;}}/\&PartialD;W,$ $A_{31\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;\mathrm{\&delta;},$ $A_{32\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;\mathrm{\&omega;},$ $A_{33\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;Z,$ $A_{34\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{Z}/\&PartialD;W,$ $A_{41\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;\mathrm{\&delta;},$ $A_{42\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;\mathrm{\&omega;},$ $A_{43\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;Z,$ $A_{44\mathrm{with}}=\&PartialD;\stackrel{\&CenterDot;}{W}/\&PartialD;W;$

(9) will contain the system linearization matrix A of blower fan group _WithCarry out depression of order and process, obtain the damping torque Matrix C (λ that generator provides _i), obtain according to step (8):

C(λ _i)=-M((A _21with+(A _23with+A _24withA _wzwith)A _zδwith+A _24withA _wδwith)A _δωwith

+(A _22with+(A _23with+A _24withA _wzwith)A _zωwith+A _24withA _wωwith)

Wherein, s=λ _i=-ξ _i± j ω _i, $A_{\mathrm{w\&delta;with}}=\&PartialD;W/\&PartialD;\mathrm{\&delta;}={(\mathrm{sI}-A_{44\mathrm{with}})}^{-1}{A}_{41\mathrm{with}},$ $A_{\mathrm{w\&omega;with}}=\&PartialD;W/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-A_{44\mathrm{with}})}^{-1}{A}_{42\mathrm{with}},$ $A_{\mathrm{wzwith}}=\&PartialD;W/\&PartialD;Z={(\mathrm{sI}-A_{44\mathrm{with}})}^{-1}{A}_{43\mathrm{with}},$ $A_{\mathrm{z\&delta;with}}=\mathrm{\&delta;Z}/\&PartialD;\mathrm{\&delta;}={(\mathrm{sI}-A_{33\mathrm{with}}-A_{34\mathrm{with}}{A}_{\mathrm{wzwith}})}^{-1}(A_{31\mathrm{with}}+A_{34\mathrm{with}}{A}_{\mathrm{w\&delta;with}}),$ $A_{\mathrm{z\&omega;with}}=\mathrm{\&delta;Z}/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-A_{33\mathrm{with}}-A_{34\mathrm{with}}{A}_{\mathrm{wzwith}})}^{-1}(A_{32\mathrm{with}}+A_{34\mathrm{with}}{A}_{\mathrm{w\&omega;with}}),$ $A_{\mathrm{\&delta;\&omega;with}}=\&PartialD;\mathrm{\&delta;}/\&PartialD;\mathrm{\&omega;}={(\mathrm{sI}-(A_{11\mathrm{with}}+(A_{13\mathrm{with}}+A_{14\mathrm{with}}{A}_{\mathrm{wzwith}})A_{\mathrm{z\&delta;with}}+A_{14\mathrm{with}}{A}_{\mathrm{w\&delta;with}}))}^{-1}$ $(A_{12\mathrm{with}}+(A_{31\mathrm{with}}+A_{14\mathrm{with}}{A}_{\mathrm{wzwith}})A_{\mathrm{z\&omega;with}}+A_{14\mathrm{with}}{A}_{\mathrm{w\&omega;with}});$

(10) can obtain according to step (7) and step (9), behind the adding blower fan, the matrix D (λ of the damping torque variable quantity that generator provides _i):

D(λ _i)=C(λ _i)-B(λ _i)；

(11) the damping torque variable quantity matrix D (λ by each generator is provided _i) be reconstructed, off-diagonal element is put in order to diagonal element, form new diagonal matrix D _New(λ _i), diagonal element D _Newj(λ _i), j=1,2 ... the damping torque that the blower fan that N represents to access provides by each generator;

(12) in the computing system each generator corresponding to specifying Oscillatory mode shape λ _iParticipation factors S _Ij:

$S_{\mathrm{ij}}=\frac{\&PartialD;{\mathrm{\&lambda;}}_i}{\&PartialD;\mathrm{\&Delta;}{D}_j}=w_i^T\frac{\&PartialD;}{\&PartialD;\mathrm{\&Delta;}{D}_j}(-M^{-1}(C\left({\mathrm{\&lambda;}}_i\right)+\mathrm{\&Delta;}{D}_j))\&times;v_i=-w_{\mathrm{ij}}\frac{1}{M_j}{v}_{\mathrm{ij}},j=\mathrm{1,2},...N$

Wherein, v _iFor corresponding to λ _iRight characteristic vector, v _IjBe v _iIn corresponding to △ ω _j, j=1,2 ... the component of N, j are j platform generator, w _iFor corresponding to λ _iLeft eigenvector, w _IjBe w _iIn corresponding to △ ω _j, j=1,2 ... the component of N;

(13) by the calculating of step (11) and step (12), can obtain the blower fan access to specifying the impact of Oscillatory mode shape:

$\mathrm{\&Delta;}{\mathrm{\&lambda;}}_i=\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}{S}_{\mathrm{ij}}{D}_{\mathrm{newj}}\left({\mathrm{\&lambda;}}_i\right)$

Wherein, _iFor specifying the change amount of Oscillatory mode shape, be used for arranging to specify the installation place of blower fan and choosing of controller.

2. be incorporated into the power networks on the damping torque analytical method of power system small signal stability impact in large-scale wind power according to claim 1 field, it is characterized in that: also comprise the reasonability that oppositely proves formula (10) by following proof procedure:

With A _WithoutAnd A _WithAt keeping characteristics root λ _iPrerequisite under carry out depression of order and process, obtain damping torque matrix B (λ _i) and C (λ _i), therefore obtain:

$w_i^T(-M^{-1}B)v_i={\mathrm{\&lambda;}}_{\mathrm{withouti}}$

$w_i^T(-M^{-1}C)v_i={\mathrm{\&lambda;}}_{\mathrm{withi}}$

Two formulas are subtracted each other, and add the variation of characteristic root behind the blower fan:

${\mathrm{\&lambda;}}_{\mathrm{withi}}-{\mathrm{\&lambda;}}_{\mathrm{withouti}}=\mathrm{\&Delta;}{\mathrm{\&lambda;}}_i=-w_i^T{M}^{-1}(C-B)v_i=-w_i^T{M}^{-1}D{v}_i$

$=\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}-w_{\mathrm{ij}}{D}_{\mathrm{newj}}{v}_{\mathrm{ij}}/M_j=\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}{S}_{\mathrm{ij}}{D}_{\mathrm{newj}}$

Wherein, λ _WithoutiBe the characteristic root of i Oscillatory mode shape of certain generator of not containing the blower fan group, λ _WithiCharacteristic root for i Oscillatory mode shape of certain generator of containing the blower fan group; Hence one can see that: the diagonal matrix D of reconstruct _New(λ _i) in diagonal element measured the damping torque that blower fan provides to each generator, and participation factors S _IjThe damping torque that provides has been provided is and how to be converted into the impact of specifying electromechanical oscillations mode; So the physical significance of the model analysis formula described in the step (13) is: specify blower fan to provide damping torque to every generator, its size is by diagonal matrix D _New(λ _i) in diagonal element D _Newj(λ _i), j=1,2 ... N tolerance; Multiply by again participation factors S _IjAfter, the damping torque of every generator acquisition just is converted into specifies blower fan by the impact of each generator on appointment electromechanical oscillations mode, and its size is by D _Newj(λ _i) S _IjTolerance, influence index △ λ _iBe exactly that blower fan passes through N platform generator to specifying the impact of electromechanical oscillations mode, be N item D _Jjnew(λ _i) S _Ij, j=1,2 ... the N sum.

Images