CN117332602A

CN117332602A - Primary frequency modulation simulation method and device for wind driven generator

Info

Publication number: CN117332602A
Application number: CN202311353815.3A
Authority: CN
Inventors: 牛高远; 胡阳; 宋子秋; 房方; 刘吉臻
Original assignee: North China Electric Power University
Current assignee: North China Electric Power University
Priority date: 2023-10-18
Filing date: 2023-10-18
Publication date: 2024-01-02
Anticipated expiration: 2043-10-18
Also published as: CN117332602B

Abstract

The application provides a primary frequency modulation simulation method and device for a wind driven generator, wherein the method comprises the following steps: acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation; according to a relational expression of electromagnetic torque and electromagnetic torque variable of the wind driven generator corresponding to a preset wind speed interval, a double-mass model and a variable pitch model of the wind driven generator are combined to construct a preset state space expression corresponding to the preset wind speed interval, wherein the preset state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variable quantity in the preset wind speed interval; determining a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs; substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation so as to perform primary frequency modulation simulation on the wind driven generator.

Description

Primary frequency modulation simulation method and device for wind driven generator

Technical Field

The application relates to the technical field of wind power generation, in particular to a primary frequency modulation simulation method and device for a wind power generator.

Background

For a variable-speed wind turbine, the current modeling method for primary frequency modulation response of a fan comprises two types: a "black box model" based on input-output data driving and a "white box modeling" based on mechanism analysis. The black box model depends on the appointed input and output characteristics, a large amount of input and output data are needed, and high-precision representation of the dynamic characteristics of the unit structure can be realized through a machine learning algorithm. However, the method has poor interpretation of the physical meaning of the parameters and has strong dependence on the accuracy and the quantity of input data. The representative simulation software of the fans such as GH Bladed and FAST in the white box model contains a high-fidelity model of high-order nonlinear dynamics, but the model is too complex and is not suitable for the design of a controller, and most white box mechanism models of primary frequency modulation response of the fans based on small signal derivation have strict requirements on input signals of frequency modulation instructions and are not suitable for the condition that the machine parameters change under actual complex working conditions.

Disclosure of Invention

Accordingly, the present application is directed to at least a method and an apparatus for simulating primary frequency modulation of a wind driven generator, which convert power variation during primary frequency modulation of a wind driven generator into state space expressions corresponding to preset state parameters under different wind speed intervals, so as to convert the power variation into linear expressions corresponding to the preset state parameters, thereby calculating a preset state parameter according to the state space expression corresponding to the wind speed interval to which the current wind speed belongs when receiving a reference output power variation of the output power of the wind driven generator, and simulating operation of the wind driven generator according to the calculated preset state parameter to simulate primary frequency modulation of the wind driven generator, so as to solve the technical problems that in the prior art, a large amount of data is required to perform primary frequency modulation simulation and a constructed model is too complex, and achieve the technical effect of improving primary frequency modulation simulation efficiency of the wind driven generator.

The application mainly comprises the following aspects:

in a first aspect, an embodiment of the present application provides a primary frequency modulation simulation method for a wind turbine, where the method includes: acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation; according to a relational expression of electromagnetic torque and electromagnetic torque variable of the wind driven generator corresponding to a preset wind speed interval, a double-mass model and a variable pitch model of the wind driven generator are combined to construct a preset state space expression corresponding to the preset wind speed interval, wherein the preset state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variable quantity in the preset wind speed interval; determining a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs; substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation so as to perform primary frequency modulation simulation on the wind driven generator.

Optionally, the preset wind speed interval includes a first preset interval and a second preset interval, an upper limit value of the first preset interval is equal to a lower limit value of the second preset interval, the preset state space expression includes a first state space expression corresponding to the first preset interval and a second state space expression corresponding to the second preset interval, the target wind speed interval is one of the first preset interval and the second preset interval, the target state space expression is one of a first state space expression and a second state space expression, and when the preset wind speed interval is the first preset interval, the first state space expression is used for representing a relationship between a preset state variable of the wind power generator and the reference output power variation when only changing electromagnetic torque of the wind power generator; and when the preset wind speed interval is the second preset interval, the second state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variation when only the pitch angle of the wind driven generator is changed.

Optionally, the preset state space expression corresponding to the preset wind speed interval is constructed in the following manner: according to a fan blade aerodynamic torque formula, constructing a linear formula of aerodynamic torque and a plurality of preset parameters, wherein each preset parameter is a parameter in the fan blade aerodynamic torque formula and generates a variable parameter in the power generation process; according to a relational expression between the electromagnetic torque of the wind driven generator and the reference electromagnetic torque, combining the linear formula, a double-mass model and a variable pitch model of the wind driven generator, and constructing a standard state space expression of the wind driven generator, wherein the standard state space expression is used for representing a relational expression between preset state parameters of the wind driven generator and variation corresponding to the preset parameters and a relational expression between output power and initial output power of the wind driven generator, and the preset state parameters are parameters which generate variation in the power generation process; when only the electromagnetic torque of the wind driven generator is changed, the first state space expression is constructed through a relational expression of the electromagnetic torque and the reference output power variation, the standard state space expression and a pitch angle first-order inertia link expression in the pitch model when a pitch angle reference value is zero; when only the pitch angle of the wind driven generator is changed, the second state space expression is constructed by a relation of a pitch angle derivative and the reference output power variation, the standard state space expression, a relation of output power and initial output power and a relation expression between the electromagnetic torque and the reference electromagnetic torque when the electromagnetic torque variable is zero.

Optionally, a linear equation of aerodynamic torque with respect to a plurality of preset parameters is constructed by:

T _r ＝a×ω _r +b×β+c×V+d

wherein P is _r The wind turbine is characterized in that the wind turbine is wind turbine pneumatic power, ρ is air density, R is wind turbine rotor radius, V is wind speed and C _P (lambda, beta) is the wind energy utilization coefficient, lambda is the tip speed ratio of the fan rotor, beta is the pitch angle, omega _r Is the rotating speed of a fan rotor, T _r Is the aerodynamic torque omega of the wind driven generator _r Beta and V are preset parameters, a is omega _r And b is a coefficient corresponding to beta, c is a coefficient corresponding to V, and d is a constant term of the linear formula.

Optionally, a relation between the electromagnetic torque of the wind driven generator and the reference electromagnetic torque is:

wherein T is _g Electromagnetic torque for the generator rotor;a derivative of the electromagnetic torque of the generator rotor; τ _g Is a generator constant; t (T) _ref The reference electromagnetic torque of the wind driven generator;

the dual mass model is represented by the following formula:

wherein J is _r The rotational inertia of a fan rotor of the wind driven generator; j (J) _g The rotational inertia of a generator rotor of the wind driven generator; t (T) _s h _aft Equivalent intermediate shaft torque between the fan rotor and the generator rotor; t (T) _r The aerodynamic torque of the wind driven generator is used; t (T) _g Electromagnetic torque for the generator rotor; θ _r Is the angular displacement of the fan rotor; θ _g Is the angular displacement of the generator rotor; θ is the angular displacement difference between the fan rotor and the generator rotor; omega _r Is the angular velocity of the fan rotor;is the derivative of the angular velocity of the fan rotor; omega _g Is the angular velocity of the generator rotor; />Is the derivative of the angular speed of the generator rotor; n (N) _g Is the gear ratio of the gear box; a is the rigidity coefficient of the equivalent intermediate shaft; b is the damping coefficient of the equivalent intermediate shaft;

the pitch model is represented by the following formula:

ω _ref ＝K _w ×(P _g,0 +ΔP _ref )

wherein,the first-order pitch angle derivative of a variable pitch model is the first-order inertia link with amplitude limiting and speed limiting; τ _β The time constant of the first-order inertia link; beta _ref Is a pitch angle reference value; beta is the pitch angle; k (K) _P For the ratio of the first-order inertial linksCoefficients; omega _r Is the angular velocity of the fan rotor; omega _ref A reference angular velocity for the fan rotor; k (K) _I Integrating coefficients for the first-order inertial links; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; k (K) _W Is the proportionality coefficient between the angular speed and the output power of the fan rotor; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref A reference output power variation amount for the generator rotor;

the standard state space expression of the wind driven generator is as follows:

ΔP _g ＝P _g -P _g.0

wherein a is a coefficient corresponding to the angular speed of the fan rotor, b is a coefficient corresponding to the pitch angle, c is a coefficient corresponding to the wind speed, and d is a constant term of the linear formula; p (P) _g The output power of a generator rotor of the wind driven generator; ΔP _g The output power variation of the generator rotor of the wind driven generator;is the derivative, ω, of the angular displacement difference between the fan rotor and the generator rotor _r 、ω _g Theta, phi, beta and T _g Is a preset state parameter.

Optionally, the relation of the electromagnetic torque and the reference output power variation is:

wherein T is _g Electromagnetic torque for the generator rotor;a derivative of the electromagnetic torque of the generator rotor; τ _g Is a generator constant; t (T) _ref The reference electromagnetic torque of the wind driven generator; k is a proportionality coefficient between electromagnetic torque and output power; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref A reference output power variation amount for the generator rotor;

the relation between the output power of the wind driven generator and the electromagnetic matrix variable is as follows:

P _g ＝η×T _g ×ω _g ＝P _g.0 +η×T _g.0 ×Δω _g +η×ΔT _g ×ω _g.0

Δω _g ＝ω _g -ω _g.0

ΔT _g ＝T _g -T _g.0

wherein η is generator efficiency; omega _g Is the angular velocity of the generator rotor; t (T) _g.0 An initial electromagnetic torque for the generator rotor; Δω _g The angular velocity variation of the generator rotor; delta T _g The electromagnetic torque variation of the generator rotor; omega _g.0 An initial angular velocity of the generator rotor;

the first-order inertia link expression of the pitch angle in the variable pitch model when the pitch angle variable is zero is as follows:

wherein beta is _ref The pitch angle reference value is zero;the first-order pitch angle derivative of a variable pitch model is the first-order inertia link with amplitude limiting and speed limiting; τ _β The time constant of the first-order inertia link; beta is the pitch angle;

the first state space expression is:

wherein a is a coefficient corresponding to the angular speed of the fan rotor, b is a coefficient corresponding to the pitch angle, c is a coefficient corresponding to the wind speed, and d is a constant term of the linear formula; omega _r Is the angular velocity of the fan rotor;is the derivative of the angular velocity of the fan rotor; omega _g Is the angular velocity of the generator rotor; />Is the derivative of the angular speed of the generator rotor; θ is the angular displacement difference between the fan rotor and the generator rotor; />Is the derivative of the angular displacement difference between the fan rotor and the generator rotor; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; v (V) ₀ Is the initial wind speed; j (J) _r The rotational inertia of a fan rotor of the wind driven generator; j (J) _g The rotational inertia of a generator rotor of the wind driven generator; n (N) _g Is the gear ratio of the gear box; a is the rigidity coefficient of the equivalent intermediate shaft; and B is the damping coefficient of the equivalent intermediate shaft.

Optionally, the relation of the pitch angle derivative and the reference output power variation is:

wherein,the first-order pitch angle derivative of a variable pitch model is the first-order inertia link with amplitude limiting and speed limiting; τ _β The time constant of the first-order inertia link; beta is the pitch angle; k (K) _P The first-order inertia link is the proportionality coefficient; omega _r Is the angular velocity of the fan rotor; k (K) _I Integrating coefficients for the first-order inertial links; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; k (K) _W Is the proportionality coefficient between the angular speed and the output power of the fan rotor; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref A reference output power variation amount for the generator rotor;

the relation between the output power and the initial output power is as follows:

wherein P is _g The output power of a generator rotor of the wind driven generator; η is generator efficiency; omega _g Is the angular velocity of the generator rotor; p (P) _g.0 An initial output power of a generator rotor of the wind power generator; omega _g.0 An initial angular velocity of the generator rotor;

the relational expression between the electromagnetic torque and the reference electromagnetic torque when the electromagnetic torque variable is zero is:

wherein T is _g Electromagnetic torque for the generator rotor;a derivative of the electromagnetic torque of the generator rotor; τ _g Is a generator constant; t (T) _g.0 An initial electromagnetic torque for the generator rotor; omega _g.0 An initial angular velocity of the generator rotor;

the second state space expression is:

wherein a is a coefficient corresponding to the angular speed of the fan rotor, b is a coefficient corresponding to the pitch angle, c is a coefficient corresponding to the wind speed, and d is a constant term of the linear formula; omega _r Is the angular velocity of the fan rotor;is the derivative of the angular velocity of the fan rotor; omega _g Is the angular velocity of the generator rotor; />Is the derivative of the angular speed of the generator rotor; θ is the angular displacement difference between the fan rotor and the generator rotor; />Is the derivative of the angular displacement difference between the fan rotor and the generator rotor; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; v (V) ₀ Is the initial wind speed; j (J) _r The rotational inertia of a fan rotor of the wind driven generator; j (J) _g The rotational inertia of a generator rotor of the wind driven generator; n (N) _g For gear-box speed changeRatio of; a is the rigidity coefficient of the equivalent intermediate shaft; and B is the damping coefficient of the equivalent intermediate shaft.

In a second aspect, an embodiment of the present application further provides a primary frequency modulation simulation device for a wind driven generator, where the device includes: the acquisition module is used for acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation; the construction module is used for constructing a preset state space expression corresponding to a preset wind speed interval according to a relational expression corresponding to the electromagnetic torque and the electromagnetic torque variable of the wind driven generator in the preset wind speed interval and combining a double-mass model and a variable pitch model of the wind driven generator, wherein the preset state space expression is used for representing the relation between the preset state variable and the reference output power variable of the wind driven generator in the preset wind speed interval; the determining module is used for determining a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs; and the calculation module is used for substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation so as to perform primary frequency modulation simulation on the wind driven generator.

In a third aspect, embodiments of the present application further provide an electronic device, including: the wind turbine primary frequency modulation simulation method comprises a processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, the processor and the memory are communicated through the bus when the electronic device is running, and the machine-readable instructions are executed by the processor to perform the steps of the wind turbine primary frequency modulation simulation method in the first aspect or any possible implementation manner of the first aspect.

In a fourth aspect, the embodiments of the present application further provide a computer readable storage medium, where a computer program is stored, where the computer program is executed by a processor to perform the steps of the wind turbine primary frequency modulation simulation method described in the first aspect or any possible implementation manner of the first aspect.

The embodiment of the application provides a primary frequency modulation simulation method and device for a wind driven generator, wherein the method comprises the following steps: acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation; according to a relational expression of electromagnetic torque and electromagnetic torque variable of the wind driven generator corresponding to a preset wind speed interval, a double-mass model and a variable pitch model of the wind driven generator are combined to construct a preset state space expression corresponding to the preset wind speed interval, wherein the preset state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variable quantity in the preset wind speed interval; determining a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs; substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation to control the wind driven generator to carry out primary frequency modulation. The power variation during primary frequency modulation of the wind driven generator is converted into state space expressions corresponding to preset state parameters respectively in different wind speed intervals, so that the power variation is converted into linear expressions corresponding to the preset state parameters, when the reference output power variation of the output power of the wind driven generator is received, the preset state parameters are calculated according to the state space expressions corresponding to the wind speed intervals to which the current wind speed belongs, the wind driven generator is operated in a simulated mode according to the calculated preset state parameters, the primary frequency modulation of the wind driven generator is simulated, the technical problem that a large amount of data are needed to perform primary frequency modulation simulation in the prior art, and the built model is too complex is solved, and the technical effect of improving the primary frequency modulation simulation efficiency of the wind driven generator is achieved.

In order to make the above objects, features and advantages of the present application more comprehensible, preferred embodiments accompanied with figures are described in detail below.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered limiting the scope, and that other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.

Fig. 1 shows a flowchart of a primary frequency modulation simulation method for a wind driven generator according to an embodiment of the present application.

Fig. 2 shows a functional block diagram of a primary frequency modulation simulation device of a wind driven generator according to an embodiment of the present application.

Fig. 3 shows a schematic structural diagram of an electronic device according to an embodiment of the present application.

Detailed Description

For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more clear, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application, and it should be understood that the accompanying drawings in the present application are only for the purpose of illustration and description, and are not intended to limit the protection scope of the present application. In addition, it should be understood that the schematic drawings are not drawn to scale. A flowchart, as used in this application, illustrates operations implemented according to some embodiments of the present application. It should be appreciated that the operations of the flow diagrams may be implemented out of order and that steps without logical context may be performed in reverse order or concurrently. Moreover, one or more other operations may be added to the flow diagrams and one or more operations may be removed from the flow diagrams as directed by those skilled in the art.

In addition, the described embodiments are only some, but not all, of the embodiments of the present application. The components of the embodiments of the present application, which are generally described and illustrated in the figures herein, may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, as provided in the accompanying drawings, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present application without making any inventive effort, are intended to be within the scope of the present application.

In the prior art, how various data of a generator should be changed when the output power is changed is predicted by a black box model driven based on input-output data, and a large amount of data is required by a prediction mode; or the model of the prediction mode is complex, and the application difficulty is high, so that how to change various data of the generator when the output power changes is predicted through 'white box modeling' based on mechanism analysis.

Based on this, the embodiment of the application provides a wind driven generator primary frequency modulation simulation method and device, by converting the power variation of the wind driven generator during primary frequency modulation into state space expressions corresponding to preset state parameters respectively under different wind speed intervals, so as to convert the power variation into linear expressions corresponding to the preset state parameters, when the reference output power variation of the output power of the generator is received, a preset state parameter is calculated according to the state space expression corresponding to the wind speed interval to which the current wind speed belongs, so that the wind driven generator operates according to the calculated preset state parameter in a simulation manner to simulate the primary frequency modulation of the wind driven generator, thereby solving the technical problems that in the prior art, a large amount of data are required to perform primary frequency modulation simulation and a constructed model is too complex, and achieving the technical effect of improving the primary frequency modulation simulation efficiency of the wind driven generator, and the method comprises the following steps:

Referring to fig. 1, fig. 1 is a flowchart of a primary frequency modulation simulation method for a wind driven generator according to an embodiment of the present application. As shown in fig. 1, the primary frequency modulation simulation method for the wind driven generator provided by the embodiment of the application includes the following steps:

s101: and acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation.

Reference output power variation Δp _ref Is at the initial output power P of the wind driven generator _g.0 On the basis of the amount of change of (2), the output power P after the change _g ＝P _g.0 ±ΔP _ref . The reference output power variation refers to a theoretical conversion amount of the wind power generator, and the actual output power variation is equal to the reference output power variation.

In the process of acquiring and controlling wind driven generator to executeBefore the primary frequency modulation reference output power variation, the wind driven generator is in a steady state, and the initial fan rotor rotating speed of the wind driven generator is omega _r.0 The initial electromagnetic torque of the wind driven generator is T _g.0 Initial electromagnetic power P of output _g,0 And an initial electromagnetic power reference value P _ref,0 Equal wind speed at this time is the initial wind speed V ₀ 。

When the frequency of the power grid changes or the voltage output by the power grid is unstable, the output power of the wind driven generator needs to be regulated, and the reference output power variation delta P for executing primary frequency modulation is sent to the wind driven generator _ref And will tell that ΔP is increasing _ref Or decrease ΔP _ref To base the wind power generator on delta P _ref And adjusting to enable the wind driven generator to reenter a steady state.

S102: and according to a relational expression of electromagnetic torque and electromagnetic torque variable of the wind driven generator corresponding to a preset wind speed interval, combining a double-mass model and a variable pitch model of the wind driven generator, and constructing a preset state space expression corresponding to the preset wind speed interval.

The preset state space expression is used for representing the relation between a preset state variable of the wind driven generator and the reference output power variation in a preset wind speed interval.

The preset wind speed interval comprises a first preset interval and a second preset interval, the upper limit value of the first preset interval is equal to the lower limit value of the second preset interval, and the preset state space expression comprises a first state space expression corresponding to the first preset interval and a second state space expression corresponding to the second preset interval.

When the preset wind speed interval is the first preset interval, the first state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variation when only the electromagnetic torque of the wind driven generator is changed. And when the preset wind speed interval is the second preset interval, the second state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variation when only the pitch angle of the wind driven generator is changed.

Wherein the upper limit value of the first preset interval and the lower limit value of the second preset interval are rated wind speeds, and when the current wind speed of the wind driven generator is rated wind speed V _rated When the wind power generator output power can be changed by changing the electromagnetic torque or the pitch angle of the wind power generator, so that the rated wind speed can be classified into a first preset interval or a second preset interval.

The first preset interval is greater than the minimum operating wind speed V _min To less than rated wind speed V _rated The second preset interval is larger than or equal to the rated wind speed to smaller than the highest running wind speed V _max The left-right opening section of (2), i.e. the first preset section is (V) _min ，V _rated ) The second preset interval is [ V ] _rated ，V _max ) The method comprises the steps of carrying out a first treatment on the surface of the Or the first preset interval is larger than the lowest running wind speed V _min To less than or equal to rated wind speed V _rated The second preset section is larger than the rated wind speed to smaller than the highest running wind speed V _max The open interval of (a), i.e. the first preset interval is (V _min ，V _rated ]The second preset interval is (V) _rated ，V _max ). Rated wind speed V _rated The determination can be made according to the model of the wind power generator, and is generally set to 11.4 meters per second.

At a current wind speed less than or equal to the minimum operating wind speed V _min When or when the current wind speed is greater than or equal to the maximum operating wind speed V _max And when the wind driven generator stops working.

That is, when the current wind speed belongs to the first preset interval, the output power of the wind power generator is changed by the reference output power variation amount by changing the electromagnetic torque of the wind power generator; when the current wind speed belongs to a second preset interval, the output power of the wind driven generator is changed by the reference output power change amount by changing the pitch angle of the wind driven generator.

The preset state space expression corresponding to the preset wind speed interval is constructed by the following steps:

according to a fan blade aerodynamic torque formula, constructing a linear formula of aerodynamic torque and a plurality of preset parameters, wherein each preset parameter is a parameter in the fan blade aerodynamic torque formula and generates a variable parameter in the power generation process;

a linear equation for pneumatic torque in relation to a plurality of preset parameters is constructed by:

T _r ＝a×ω _r +b×β+c×V+d(4)

in the formulas (1) to (4), P _r The wind turbine is characterized in that the wind turbine is wind turbine pneumatic power, ρ is air density, R is wind turbine rotor radius, V is wind speed and C _P (lambda, beta) is the wind energy utilization coefficient, lambda is the tip speed ratio of the fan rotor, beta is the pitch angle, omega _r Is the rotating speed of a fan rotor, T _r Is the aerodynamic torque omega of the wind driven generator _r Beta and V are preset parameters, a is omega _r And b is a coefficient corresponding to beta, c is a coefficient corresponding to V, and d is a constant term of the linear formula.

Wherein C is _P (lambda, beta) is a high order nonlinear function between tip speed ratio and pitch angle. That is, as can be seen from the combination of the formula (3), T in the pneumatic system _r Is omega of _r The highly nonlinear functions of β and V are difficult to describe with simple mathematical models. Further, simplification can be achieved by linearization. The PWA model (Piece Wise Affine ) is selected, and the linear model is utilized to approximately represent the pneumatic power by dividing different operation conditionsNonlinear characteristics, and thus equation (4).

Different running conditions divided by the method are divided into a first preset interval or a second preset interval, and when the current wind speed belongs to the first preset interval, all preset state parameters when the output power changes the reference output power variation are predicted through a first state space expression corresponding to the first preset interval; when the current wind speed belongs to a second preset interval, predicting each preset state parameter when the output power changes by the reference output power variation through a second state space expression corresponding to the second preset interval.

And according to a relational expression between the electromagnetic torque of the wind driven generator and the reference electromagnetic torque, combining the linear formula, a double-mass model and a variable pitch model of the wind driven generator to construct a standard state space expression of the wind driven generator.

The standard state space expression is used for representing a relational expression between preset state parameters of the wind driven generator and variation amounts respectively corresponding to the preset parameters and a relational expression between output power and initial output power of the wind driven generator, and the preset state parameters are parameters which change in the power generation process.

Because the electromagnetic transient process of the generator is usually only millisecond, the variable frequency driving model formed by the generator and the frequency converter in the generator system can be accurately replaced by a first-order linear link to be approximate.

Further, a relational expression between the electromagnetic torque of the wind turbine and the reference electromagnetic torque is:

in the formula (5), T _g Electromagnetic torque for the generator rotor;a derivative of the electromagnetic torque of the generator rotor; τ _g Is a generator constant; t (T) _ref Reference electricity for wind power generatorMagnetic torque, i.e. in theory the electromagnetic torque of a wind power generator. And the electromagnetic torque of the wind power generator after the actual transformation is equal to the reference battery torque.

Due to the proportional relationship between electromagnetic torque and output power, and with reference to output power P _ref ＝ΔP _ref +P _g.0 Wherein ΔP _ref To reference the output power variation, P _g.0 Is the initial output power.

Further, the expression between the reference electromagnetic torque and the reference output power variation is:

T _ref ＝K×P _ref ＝K×ΔP _ref +K×P _g.0 (6)

in the formula (6), T _ref A reference electromagnetic torque for the generator rotor; k is a proportionality coefficient between electromagnetic torque and output power; p (P) _ref Is the reference output power; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref Is the reference output power variation of the generator rotor.

Substituting formula (6) into formula (5) to obtain a relational expression of the electromagnetic torque and the reference output power variation, wherein the relational expression is as follows:

in the formula (7), T _g Electromagnetic torque for the generator rotor;a derivative of the electromagnetic torque of the generator rotor; τ _g Is a generator constant; t (T) _ref The reference electromagnetic torque of the wind driven generator; k is a proportionality coefficient between electromagnetic torque and output power; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref Is the reference output power variation of the generator rotor.

The dual mass model is represented by the following formula:

in the formula (8), J _r The rotational inertia of a fan rotor of the wind driven generator; j (J) _g The rotational inertia of a generator rotor of the wind driven generator; t (T) _shaft Equivalent intermediate shaft torque between the fan rotor and the generator rotor; t (T) _r The aerodynamic torque of the wind driven generator is used; t (T) _g Electromagnetic torque for the generator rotor; θ _r Is the angular displacement of the fan rotor; θ _g Is the angular displacement of the generator rotor; θ is the angular displacement difference between the fan rotor and the generator rotor; omega _r Is the angular velocity of the fan rotor;is the derivative of the angular velocity of the fan rotor; omega _g Is the angular velocity of the generator rotor; />Is the derivative of the angular speed of the generator rotor; n (N) _g Is the gear ratio of the gear box; a is the rigidity coefficient of the equivalent intermediate shaft; and B is the damping coefficient of the equivalent intermediate shaft.

Further, derivative of θ can be obtained:

the double-mass model is characterized in that mechanical energy generated by a fan rotor of a wind driven generator is firstly transmitted to a low-speed shaft of a gear box, and then transmitted to a generator rotor part by a high-speed shaft after the speed change function of the gear box, so that the fan rotor is rotated to drive the generator rotor to move, and wind energy is converted into electric energy.

The pitch model is represented by the following formula:

equation (10),the first-order pitch angle derivative of a variable pitch model is the first-order inertia link with amplitude limiting and speed limiting; τ _β The time constant of the first-order inertia link; beta _ref Is a pitch angle reference value; beta is the pitch angle.

The expression between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor is:

β _ref ＝K _P ×(ω _r -ω _ref )+K _I ×∫(ω _r -ω _ref )＝K _P ×(ω _r -ω _ref )+K _I ×φ (11)

in the formula (11), K _P The first-order inertia link is the proportionality coefficient; omega _r Is the angular velocity of the fan rotor; omega _ref A reference angular velocity for the fan rotor; k (K) _I Integrating coefficients for the first-order inertial links;the derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; k (K) _W Is the proportionality coefficient between the angular velocity of the fan rotor and the output power.

That is, in order to write into a linear expression, Φ= ≡ (ω) _r -ω _ref ) And then obtain:

/>

the expression between the reference angular velocity and the reference output power variation is:

ω _ref ＝K _W ×(P _g,0 +ΔP _ref ) (13)

in formula (13), ω _ref A reference angular velocity for the fan rotor; k (K) _W Is the proportionality coefficient between the angular speed and the output power of the fan rotor; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref Is the reference output power variation of the generator rotor.

Substituting equations (11) and (13) into the expression between the first order pitch angle derivative of the variable pitch model and the reference output power variation in equation (10):

In the formula (14), K _P The first-order inertia link is the proportionality coefficient; omega _r Is the angular velocity of the fan rotor; omega _ref A reference angular velocity for the fan rotor; k (K) _I Integrating coefficients for the first-order inertial links;the derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; k (K) _W Is the proportionality coefficient between the angular speed and the output power of the fan rotor; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref Is the reference output power variation of the generator rotor.

And combining the formula (5), the formula (8), the formula (9), the formula (10) and the formula (12) to obtain the standard state space expression of the wind driven generator.

The standard state space expression of the wind driven generator is as follows:

in the formula (15), a is a coefficient corresponding to the angular speed of the fan rotor, b is a coefficient corresponding to the pitch angle, c is a coefficient corresponding to the wind speed, and d is a constant term of the linear formula; p (P) _g The output power of a generator rotor of the wind driven generator; ΔP _g The output power variation of the generator rotor of the wind driven generator;is the derivative of the angular displacement difference between the fan rotor and the generator rotor; omega _r 、ω _g Theta, phi, beta and T _g Is a preset state parameter. / >

That is, the derivative corresponding to each preset state parameter is written in the form of a state space expression so that the derivative corresponding to each preset state parameter is equal to the preset parameter (the preset parameter at this time is ω _ref 、β _ref 、T _ref ) Related to the following.

When only the electromagnetic torque of the wind driven generator is changed, the first state space expression is constructed through a relational expression of the electromagnetic torque and the reference output power variation, the standard state space expression and a pitch angle first-order inertia link expression in the pitch model when a pitch angle reference value is zero.

That is, when the current wind speed belongs to the first preset interval, only the electromagnetic torque of the wind driven generator is changed without changing the pitch angle, and then beta _ref Is 0 and is equal to the reference angular velocity omega of the fan rotor _ref Nor is it relevant.

At this time, the relation between the output power of the wind power generator and the electromagnetic matrix variable is as follows:

P _g ＝η×T _g ×ω _g ＝P _g.0 +η×T _g.0 ×Δω _g +η×ΔT _g ×ω _g.0 (16)

Δω _g ＝ω _g -ω _g.0 (17)

ΔT _g ＝T _g -T _g.0 (18)

in formulas (16) to (18), η is generator efficiency; omega _g Is the angular velocity of the generator rotor; t (T) _g.0 An initial electromagnetic torque for the generator rotor; Δω _g Is the angular velocity variation of the generator rotor and; delta T _g The electromagnetic torque variation of the generator rotor; omega _g.0 Is the initial angular velocity of the generator rotor.

P _g.0 ＝η×T _g.0 ×ω _g.0 (19)

Substituting equations (17) and (18) into equation (16) yields:

P _g ＝P _g.0 +η×T _g.0 ×ω _g +η×T _g ×ω _g.0 -2×η×T _g.0 ×ω _g.0

＝-P _g.0 +η×T _g.0 ×ω _g +η×T _g ×ω _g.0 (20)

further, ΔP _g ＝P _g -P _g.0 ＝-2×P _g.0 +η×T _g.0 ×ω _g +η×T _g ×ω _g.0 。

in the formula (21), beta _ref The pitch angle reference value is zero;the first-order pitch angle derivative of a variable pitch model is the first-order inertia link with amplitude limiting and speed limiting; τ _β The time constant of the first-order inertia link; beta is the pitch angle.

Further, the first state space expression may be obtained by combining equation (7), equation (15), equation (20), and equation (21):

in the formula (22), a is a coefficient corresponding to the angular speed of the fan rotor, b is a coefficient corresponding to the pitch angle, c is a coefficient corresponding to the wind speed, and d is a constant term of the linear formula; omega _r Is the angular velocity of the fan rotor;is the derivative of the angular velocity of the fan rotor; omega _g Is the angular velocity of the generator rotor; />Is the derivative of the angular speed of the generator rotor; θ is the angular displacement difference between the fan rotor and the generator rotor; />Is the derivative of the angular displacement difference between the fan rotor and the generator rotor; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; v (V) ₀ Is the initial wind speed; j (J) _r The rotational inertia of a fan rotor of the wind driven generator; j (J) _g The rotational inertia of a generator rotor of the wind driven generator; n (N) _g Is the gear ratio of the gear box; a is the rigidity coefficient of the equivalent intermediate shaft; and B is the damping coefficient of the equivalent intermediate shaft.

The initial wind speed may be understood as the wind speed at which the wind turbine is in a steady state condition before the reference output power variation for controlling the wind turbine to perform primary frequency modulation is obtained. When only the pitch angle of the wind driven generator is changed, the second state space expression is constructed by a relation of a pitch angle derivative and the reference output power variation, the standard state space expression, a relation of output power and initial electromagnetic torque, and a relation expression between the electromagnetic torque and the reference electromagnetic torque when the electromagnetic torque variable is zero.

That is, when the current wind speed belongs to the second preset interval, the change of the output power is achieved only by changing the pitch angle.

When the current wind speed belongs to the second preset interval, the wind driven generator is in a power limiting state, the electromagnetic torque keeps unchanged at the initial value, and further,

In the formula (23), P _g The output power of a generator rotor of the wind driven generator; η is generator efficiency; omega _g Is the angular velocity of the generator rotor; p (P) _g.0 An initial output power of a generator rotor of the wind power generator; omega _g.0 Is the initial angular velocity of the generator rotor.

in the formula (24), T _g Electromagnetic torque for the generator rotor;a derivative of the electromagnetic torque of the generator rotor; τ _g Is a generator constant; t (T) _g.0 An initial electromagnetic torque for the generator rotor; omega _g.0 Is the initial angular velocity of the generator rotor.

Combining equation (12), equation (13), equation (14), equation (23), and equation (24) to obtain the second state space expression:

/>

in the formula (25), a is a coefficient corresponding to the angular speed of the fan rotor, b is a coefficient corresponding to the pitch angle, c is a coefficient corresponding to the wind speed, and d is a constant term of the linear formula; omega _r Is the angular velocity of the fan rotor;is a fan rotorDerivative of angular velocity of the son; omega _g Is the angular velocity of the generator rotor; />Is the derivative of the angular speed of the generator rotor; θ is the angular displacement difference between the fan rotor and the generator rotor; / >Is the derivative of the angular displacement difference between the fan rotor and the generator rotor; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; v (V) ₀ Is the initial wind speed; j (J) _r The rotational inertia of a fan rotor of the wind driven generator; j (J) _g The rotational inertia of a generator rotor of the wind driven generator; n (N) _g Is the gear ratio of the gear box; a is the rigidity coefficient of the equivalent intermediate shaft; and B is the damping coefficient of the equivalent intermediate shaft.

S103: and determining a target state space expression corresponding to the target wind speed interval according to the target wind speed interval to which the current wind speed belongs.

The target wind speed interval is one of the first preset interval and the second preset interval, and the target state space expression is one of a first state space expression and a second state space expression.

The current wind speed may be understood as the wind speed when the reference output power variation amount for controlling the wind power generator to perform primary frequency modulation is acquired.

And when the current wind speed belongs to the first preset interval, the target wind speed interval is the first preset interval, and then the target state space expression is determined to be the first state space expression, namely, when only the electromagnetic torque of the wind driven generator is changed, the state space expression between the reference output power variation and each preset state variable is determined.

When the current wind speed belongs to a second preset interval, the target wind speed interval is the second preset interval, and then the target state space expression is determined to be the second state space expression, namely, when the pitch angle of the wind driven generator is only changed, the state space expression between the reference output power variation and each preset state variable is determined.

S104: substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation so as to perform primary frequency modulation simulation on the wind driven generator.

After the target state space expression is determined, the reference output power variation is substituted into the target state space expression, so that the value of each preset state variable is determined, and the wind driven generator is controlled to simulate according to the calculated preset state variable to simulate the primary frequency modulation of the wind driven generator.

By way of example, the first state space expression may also be reduced to:

wherein,x＝[ω _r ω _g θ φ β T _g ] ^T ，u＝ΔP _ref ，/>v＝[V ₀ P _g.0 ] ^T ，/>y＝ΔP _g ，/>F ₁ ＝-2×P _g.0 。

by way of example, the second state space expression may also be reduced to:

wherein,x＝[ω _r ω _g θφβT _g ] ^T ，u＝Δp _ref ，v＝[V ₀ p _g.0 ] ^T ，/>y＝ΔP _g ，/>F ₂ ＝-P _g.0 。

furthermore, the fact that the approximation degree of the low-order mechanism model to the complex nonlinearity is limited is considered at the same time, so that a neural network algorithm can be introduced to compensate the mechanism model deviation, the neural network model is used for predicting the first compensation and the second compensation, and a state space expression based on the neural network model is as follows:

In the formula (27), A _i 、B _i 、C _i 、D _i 、E _i 、F _i And f is a first compensation, and g is a second compensation corresponding to the first state space expression or the second state space expression.

That is, a "black box model" based on input-output data driving or a "white box modeling" based on mechanism analysis is predictedPredicted +.f. from the first state space expression or the second state space expression>Performing difference to obtain a first difference value; the predicted y of the black box model based on the input-output data driving or the white box modeling based on the mechanism analysis is differenced with the predicted y of the first state space expression or the second state space expression to obtain a second difference value, and each time the predicted y passes through the first state space expression or the second state space expression>And y is taken as sample data, and the corresponding first difference value and second difference value are taken as labels to train the neural network model, so that the neural network model can predict the first difference value and the second difference value, namely, the first compensation and the second compensation, and more accurate can be obtained without using a black box model or white box modeling>And y, thereby obtaining more accurate x= [ omega ] _r ω _g θφβT _g ] ^T 。

Because the system state quantity in the continuous state space equation can be observed, the joint identification of the wind turbine mechanism parameter and the data parameter in the PWA model in the state space equation can be completed based on the operation data in the actual working field of the wind turbine, that is, the coefficient such as the rigidity coefficient A of the equivalent intermediate shaft, the damping coefficient B of the equivalent intermediate shaft and the like can be inaccurate in value due to the service life and the like, so that the more accurate A, B can be conveniently obtained through joint identification, the calculated preset state variable is more accurate, and the sub-region transient characteristics of a plurality of sub-models approach the nonlinear transient characteristics of the system overall.

Based on the same application conception, the embodiment of the application also provides a wind driven generator primary frequency modulation simulation device corresponding to the wind driven generator primary frequency modulation simulation method provided by the embodiment, and because the principle of solving the problem by the device in the embodiment of the application is similar to that of the wind driven generator primary frequency modulation simulation method provided by the embodiment of the application, the implementation of the device can be referred to the implementation of the method, and the repetition is omitted.

Fig. 2 is a functional block diagram of a primary frequency modulation simulation device of a wind driven generator according to an embodiment of the present application. The wind power generator primary frequency modulation simulation device 10 includes: an acquisition module 101, a construction module 102, a determination module 103 and a calculation module 104.

An acquisition module 101, configured to acquire a reference output power variation amount for controlling the wind turbine to perform primary frequency modulation; the construction module 102 is configured to construct a preset state space expression corresponding to a preset wind speed interval according to a relational expression corresponding to an electromagnetic torque and an electromagnetic torque variable of the wind turbine in the preset wind speed interval, and by combining a dual-mass model and a variable pitch model of the wind turbine, where the preset state space expression is used for representing a relation between a preset state variable of the wind turbine and the reference output power variable in the preset wind speed interval; a determining module 103, configured to determine a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs; and the calculating module 104 is configured to substitute the reference output power variation into the target state space expression, calculate a preset state variable of primary frequency modulation, and perform primary frequency modulation simulation on the wind driven generator.

Based on the same application concept, referring to fig. 3, which is a schematic structural diagram of an electronic device provided in an embodiment of the present application, the electronic device 20 includes: the wind turbine generator system comprises a processor 201, a memory 202 and a bus 203, wherein the memory 202 stores machine-readable instructions executable by the processor 201, and when the electronic device 20 is running, the processor 201 and the memory 202 communicate through the bus 203, and the machine-readable instructions are executed by the processor 201 to perform the steps of the wind turbine generator system primary frequency modulation simulation method according to any one of the above embodiments.

In particular, the machine readable instructions, when executed by the processor 201, may perform the following: acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation; according to a relational expression of electromagnetic torque and electromagnetic torque variable of the wind driven generator corresponding to a preset wind speed interval, a double-mass model and a variable pitch model of the wind driven generator are combined to construct a preset state space expression corresponding to the preset wind speed interval, wherein the preset state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variable quantity in the preset wind speed interval; determining a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs; substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation so as to perform primary frequency modulation simulation on the wind driven generator.

Based on the same application conception, the embodiment of the application also provides a computer readable storage medium, wherein a computer program is stored on the computer readable storage medium, and the computer program executes the steps of the wind driven generator primary frequency modulation simulation method provided by the embodiment when being run by a processor.

Specifically, the storage medium can be a general storage medium, such as a mobile magnetic disk, a hard disk, and the like, when a computer program on the storage medium is run, the above-mentioned wind turbine primary frequency modulation simulation method can be executed, by converting the power variation of the wind turbine primary frequency modulation into state space expressions corresponding to preset state parameters respectively under different wind speed intervals, so as to convert the power variation into linear expressions corresponding to the preset state parameters, when the reference output power variation of the output power of the wind turbine is received, the preset state parameters are calculated according to the state space expressions corresponding to the wind speed intervals to which the current wind speed belongs, so that the wind turbine is simulated to operate according to the calculated preset state parameters, and the technical problems that in the prior art, a large amount of data are required to perform primary frequency modulation simulation and a built model is too complex are solved, the technical effects of improving the primary frequency modulation simulation efficiency of the wind turbine and simplifying the model are achieved.

It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described system and apparatus may refer to corresponding procedures in the foregoing method embodiments, which are not described herein again. In the several embodiments provided in this application, it should be understood that the disclosed systems, devices, and methods may be implemented in other manners. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be other manners of division in actual implementation, and for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some communication interface, device or unit indirect coupling or communication connection, which may be in electrical, mechanical or other form.

The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

In addition, each functional unit in each embodiment of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.

The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer readable storage medium executable by a processor. Based on such understanding, the technical solutions of the present application may be embodied in essence or a part contributing to the prior art or a part of the technical solutions, or in the form of a software product, which is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random access Memory (Random AccessMemory, RAM), a magnetic disk, or an optical disk, or other various media capable of storing program codes.

The foregoing is merely a specific embodiment of the present application, but the protection scope of the present application is not limited thereto, and any person skilled in the art can easily think about changes or substitutions within the technical scope of the present application, and the changes or substitutions are covered in the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims

1. A wind power generator primary frequency modulation simulation method, which is characterized by comprising the following steps:

acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation;

according to a relational expression of electromagnetic torque and electromagnetic torque variable of the wind driven generator corresponding to a preset wind speed interval, a double-mass model and a variable pitch model of the wind driven generator are combined to construct a preset state space expression corresponding to the preset wind speed interval, wherein the preset state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variable quantity in the preset wind speed interval;

determining a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs;

substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation so as to perform primary frequency modulation simulation on the wind driven generator.

2. The method of claim 1, wherein the preset wind speed interval includes a first preset interval and a second preset interval, an upper limit value of the first preset interval is equal to a lower limit value of the second preset interval, the preset state space expression includes a first state space expression corresponding to the first preset interval and a second state space expression corresponding to the second preset interval, the target wind speed interval is one of the first preset interval and the second preset interval, the target state space expression is one of the first state space expression and the second state space expression,

when the preset wind speed interval is the first preset interval, the first state space expression is used for representing the relation between a preset state variable of the wind driven generator and the reference output power variation when only the electromagnetic torque of the wind driven generator is changed;

and when the preset wind speed interval is the second preset interval, the second state space expression is used for representing the relation between the preset state variable of the wind driven generator and the reference output power variation when only the pitch angle of the wind driven generator is changed.

3. The method according to claim 2, wherein the preset state space expression corresponding to the preset wind speed interval is constructed by:

according to a relational expression between the electromagnetic torque of the wind driven generator and the reference electromagnetic torque, combining the linear formula, a double-mass model and a variable pitch model of the wind driven generator, and constructing a standard state space expression of the wind driven generator, wherein the standard state space expression is used for representing a relational expression between preset state parameters of the wind driven generator and variation corresponding to the preset parameters and a relational expression between output power and initial output power of the wind driven generator, and the preset state parameters are parameters which generate variation in the power generation process;

when only the electromagnetic torque of the wind driven generator is changed, the first state space expression is constructed through a relational expression of the electromagnetic torque and the reference output power variation, the standard state space expression and a pitch angle first-order inertia link expression in the pitch model when a pitch angle reference value is zero;

When only the pitch angle of the wind driven generator is changed, the second state space expression is constructed by a relation of a pitch angle derivative and the reference output power variation, the standard state space expression, a relation of output power and initial output power and a relation expression between the electromagnetic torque and the reference electromagnetic torque when the electromagnetic torque variable is zero.

4. A method according to claim 3, characterized in that the linear formula of the aerodynamic torque in relation to a plurality of preset parameters is constructed by:

T _r ＝a×ω _r +b×β+c×V+d

5. A method according to claim 3, wherein the relationship between the electromagnetic torque of the wind turbine and the reference electromagnetic torque is:

the dual mass model is represented by the following formula:

the pitch model is represented by the following formula:

β _ref ＝K _P ×(ω _r -ω _ref )+K _I ×∫(ω _r -ω _ref )

＝K _P ×(ω _r -ω _ref )+K _I ×φ

ω _ref ＝K _W ×(P _g,0 +ΔP _ref )

wherein,the first-order pitch angle derivative of a variable pitch model is the first-order inertia link with amplitude limiting and speed limiting; τ _β The time constant of the first-order inertia link; beta _ref Is a pitch angle reference value; beta is the pitch angle; k (K) _P The first-order inertia link is the proportionality coefficient; omega _r Is the angular velocity of the fan rotor; omega _ref A reference angular velocity for the fan rotor; k (K) _I For the first-order inertiaAn integral coefficient of the sexual link; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; k (K) _W Is the proportionality coefficient between the angular speed and the output power of the fan rotor; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref A reference output power variation amount for the generator rotor;

the standard state space expression of the wind driven generator is as follows:

ΔP _g ＝P _g -P _g.0

6. A method according to claim 3, wherein the electromagnetic torque is related to the reference output power variation by the relation:

P _g ＝η×T _g ×ω _g ＝P _g.0 +η×T _g.0 ×Δω _g +η×ΔT _g ×ω _g.0

Δω _g ＝ω _g -ω _g.0

ΔT _g ＝T _g -T _g.0

wherein beta is _ref The pitch angle reference value is zero;for pitch-changing modesThe first-order pitch angle derivative is used for the variable pitch model, and the variable pitch model is a first-order inertia link with amplitude limiting and speed limiting; τ _β The time constant of the first-order inertia link; beta is the pitch angle;

the first state space expression is:

Wherein a is a coefficient corresponding to the angular speed of the fan rotor, b is a coefficient corresponding to the pitch angle, c is a coefficient corresponding to the wind speed, and d is a constant term of the linear formula; omega _r Is the angular velocity of the fan rotor;is the derivative of the angular velocity of the fan rotor; omega _g Is the angular velocity of the generator rotor; />Is the derivative of the angular speed of the generator rotor; θ is the angular displacement difference between the fan rotor and the generator rotor; />Is the derivative of the angular displacement difference between the fan rotor and the generator rotor; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; v (V) ₀ Is the initial wind speed; j (J) _r The rotational inertia of a fan rotor of the wind driven generator; j (J) _g The rotational inertia of a generator rotor of the wind driven generator; n (N) _g Is the gear ratio of the gear box; a is the rigidity coefficient of the equivalent intermediate shaftThe method comprises the steps of carrying out a first treatment on the surface of the And B is the damping coefficient of the equivalent intermediate shaft.

7. A method according to claim 3, wherein the pitch angle derivative is related to the reference output power variation by:

wherein,the first-order pitch angle derivative of a variable pitch model is the first-order inertia link with amplitude limiting and speed limiting; sigma (sigma) _β The time constant of the first-order inertia link; beta is the pitch angle; k (K) _P The first-order inertia link is the proportionality coefficient; omega _r Is the angular velocity of the fan rotor; k (K) _I Integrating coefficients for the first-order inertial links; />The derivative of phi is the difference between the angular velocity of the fan rotor and the reference angular velocity of the fan rotor; k (K) _W Is the proportionality coefficient between the angular speed and the output power of the fan rotor; p (P) _g,0 An initial output power of a generator rotor of the wind power generator; ΔP _ref A reference output power variation amount for the generator rotor;

wherein P is _g The output power of a generator rotor of the wind driven generator; η is generator efficiency; omega _g Is the angular velocity of the generator rotor; p (P) _g.0 An initial output power of a generator rotor of the wind power generator; omega _g.0 Is an electric generatorAn initial angular velocity of the rotor;

the second state space expression is:

8. A wind turbine primary frequency modulation simulation device, the device comprising:

the acquisition module is used for acquiring a reference output power variation quantity for controlling the wind driven generator to execute primary frequency modulation;

the construction module is used for constructing a preset state space expression corresponding to a preset wind speed interval according to a relational expression corresponding to the electromagnetic torque and the electromagnetic torque variable of the wind driven generator in the preset wind speed interval and combining a double-mass model and a variable pitch model of the wind driven generator, wherein the preset state space expression is used for representing the relation between the preset state variable and the reference output power variable of the wind driven generator in the preset wind speed interval;

The determining module is used for determining a target state space expression corresponding to a target wind speed interval according to the target wind speed interval to which the current wind speed belongs;

and the calculation module is used for substituting the reference output power variation into the target state space expression, and calculating a preset state variable of primary frequency modulation so as to perform primary frequency modulation simulation on the wind driven generator.

9. An electronic device, comprising: a processor, a memory and a bus, said memory storing machine readable instructions executable by said processor, said processor and said memory communicating via said bus when the electronic device is operating, said machine readable instructions when executed by said processor performing the steps of the wind turbine primary frequency modulation simulation method according to any one of claims 1 to 7.

10. A computer readable storage medium, characterized in that the computer readable storage medium has stored thereon a computer program which, when executed by a processor, performs the steps of the wind turbine primary frequency modulation simulation method according to any of claims 1 to 7.