CN104167735B - A kind of non-mechanism equivalent modeling method of wind energy turbine set and device - Google Patents

Abstract

The invention provides the non-mechanism equivalent modeling method of a kind of wind energy turbine set and device, with the overall external characteristic of wind energy turbine set for modeling object, wind speed is inputted for input variable respectively with wind farm grid-connected point voltage, wind energy turbine set, all with the response of the overall active power of wind energy turbine set and reactive power for output variable, by parameter identification method set up respectively wind energy turbine set response and grid-connected point voltage change between transfer function, and wind energy turbine set respond and wind speed change between transfer function.The non-mechanism equivalent model of complete wind energy turbine set is made up of the product of these two transfer functions.The method that the present invention proposes simplifies the process of wind energy turbine set Equivalent Modeling and required data volume, has agreed with again the final purpose of wind energy turbine set Equivalent Modeling well simultaneously, and namely setting up can the simple lower-order model of accurate description wind energy turbine set external characteristic.The non-mechanism equivalent model of the wind energy turbine set set up by the inventive method can be widely used in the simulation analysis containing wind-powered electricity generation electric power system, has good future in engineering applications.

Description

A kind of non-mechanism equivalent modeling method of wind energy turbine set and device

Technical field

The present invention relates to renewable energy power generation field, relate to a kind of non-mechanism equivalent modeling method of wind energy turbine set.

Background technology

Rapidly, the electric pressure that the scale of wind power integration electrical network constantly increases, accesses electrical network improves constantly in China's wind power generation development in recent years.The concentrated access of large-scale wind power can bring problems to electric power system, when studying wind power integration to the affecting of power system dynamic stability, need one can the wind energy turbine set model of accurate characterization wind energy turbine set dynamic characteristic.Digital Simulation is wind energy turbine set models applying affects problem on Electrical Power System Dynamic major way in research wind power integration.For a large-scale wind power field, if carry out detailed modeling to it, set up the model of each module in every typhoon group of motors and unit in wind energy turbine set, though ensure that the dynamic accuracy of wind energy turbine set, but simulation calculation is carried out in huge wind energy turbine set model access electric power system, not only amount of calculation is large, and is difficult to convergence.On the other hand, in the electric system simulation that accesses large-scale wind power is analyzed, primary study be the dynamic characteristic of wind energy turbine set entirety, instead of every platform unit dynamically concrete in wind energy turbine set.Therefore, Equivalent Modeling is carried out to wind energy turbine set, set up accurate and simple wind energy turbine set equivalent model and be necessary.

At present, wind energy turbine set Equivalent Modeling mainly adopts the mode of modelling by mechanism, has proposed much concrete method.The common feature of these methods is, must know the detailed model etc. of wind speed that Equivalent Modeling moment every typhoon power generator bears, the inner various electric equipment of wind energy turbine set, when these information are known, and the precision that the wind energy turbine set equivalent model of acquisition is all specifically good.But, from current engineering actual state, the wind speed that in wind energy turbine set, every typhoon power generator bears cannot be learnt by power grid operation enterprise, also the detailed model of the inner various electric equipment of wind energy turbine set is lacked, the information that can obtain often only has that the overall power of wind energy turbine set exports, the wind speed etc. of wind energy turbine set annex anemometer tower, therefore will carry out mechanism Equivalent Modeling to wind energy turbine set and be difficult to realize.

Summary of the invention

When carrying out wind energy turbine set Equivalent Modeling for prior art from grid side, effective information only has that the overall power of wind energy turbine set exports, the wind speed etc. of wind energy turbine set annex anemometer tower, so that the technical problem of mechanism Equivalent Modeling cannot be carried out, propose a kind of non-mechanism equivalent modeling method of wind energy turbine set.

Another object of the present invention also proposes a kind of non-mechanism Equivalent Modeling device of wind energy turbine set.

Above-mentioned purpose of the present invention is realized by the technical characteristic of independent claims, and dependent claims develops the technical characteristic of independent claims with alternative or favourable mode.

For reaching above-mentioned purpose, the technical solution adopted in the present invention is as follows:

A kind of non-mechanism equivalent modeling method of wind energy turbine set, use the wind farm grid-connected voltage of point and the wind speed of anemometer tower to be input variable, use the response of wind energy turbine set overall power to be output variable, set up the non-mechanism equivalent model of wind energy turbine set based on transfer function by discrimination method, this non-mechanism equivalent modeling method comprises the following steps:

Steps A: the equation setting up the non-mechanism equivalent model of wind energy turbine set, concrete steps are as follows:

Steps A-1: for the situation of wind energy turbine set overall dynamics response under electric network fault, now ignore the change of wind speed, set up the non-mechanism equivalent model of wind farm grid-connected point voltage and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pu}}\left(U\right)\\ f_{\mathrm{Pu}}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{\mathrm{Qu}}\left(U\right)\\ f_{\mathrm{Qu}}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U is wind farm grid-connected point voltage; H _pu(s) and H _qus () is respectively transfer function equation; g _puand g (U) _qu(U) algebraic equation that to be respectively with wind farm grid-connected point voltage U be independent variable; Transfer function equation H _pu(s) and H _qus () is expressed as follows respectively:

$\left\{\begin{array}{c}{H}_{\mathrm{Pu}}\left(s\right)=\frac{b_{\mathrm{pu}0}{s}^{\mathrm{pn}}+b_{\mathrm{pu}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pun}}}{s^{\mathrm{pm}}+a_{\mathrm{pu}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pum}}}\\ H_{\mathrm{Qu}}\left(s\right)=\frac{b_{\mathrm{qu}0}{s}^{\mathrm{qn}}+b_{\mathrm{qu}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qun}}}{s^{\mathrm{qm}}+a_{\mathrm{qu}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qum}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pu0, b _pu1..., b _pun, a _pu1..., a _pum, b _qu0, b _qu1..., b _qun, a _qu1..., a _qumfor constant coefficient; Algebraic equation g _puand g (U) _qu(U) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\\ g_{\mathrm{Qu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\end{array}\right.$

In formula, U, U ₀be respectively wind farm grid-connected point voltage and grid-connected point voltage steady-state value;

Steps A-2: for the situation of wind energy turbine set electrical power change under fluctuations in wind speed, set up the non-mechanism equivalent model of wind energy turbine set input wind speed and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pv}}\left(V\right)\\ f_{\mathrm{Pv}}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ Q=Q_0{f}_{\mathrm{Qv}}\left(V\right)\\ f_{\mathrm{Qv}}\left(V\right)=H_{\mathrm{Qv}}\left(s\right)g_{\mathrm{Qv}}\left(V\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; V is wind energy turbine set input wind speed; H _pv(s) and H _qvs () is respectively transfer function equation; g _pvand g (V) _qv(V) algebraic equation that to be respectively with wind energy turbine set input wind speed V be independent variable, this transfer function equation H _pv(s) and H _qvs () expression formula is:

$\left\{\begin{array}{c}{H}_{\mathrm{Pv}}\left(s\right)=\frac{b_{\mathrm{pv}0}{s}^{\mathrm{pn}}+b_{\mathrm{pv}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pvn}}}{s^{\mathrm{pm}}+a_{\mathrm{pv}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pvm}}}\\ H_{\mathrm{Qv}}\left(s\right)=\frac{b_{\mathrm{qv}0}{s}^{\mathrm{qn}}+b_{\mathrm{qv}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qvn}}}{s^{\mathrm{qm}}+a_{\mathrm{qv}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qvm}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pv0, b _pv1..., b _pvn, a _pv1..., a _pvm, b _qv0, b _qv1..., b _qvn, a _qv1..., a _qvmfor constant coefficient;

Algebraic equation g _pvand g (V) _qv(V) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\\ g_{\mathrm{Qv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\end{array}\right.$

In formula, V is wind energy turbine set input wind speed; V ₀for the steady-state value of wind energy turbine set input wind speed;

Steps A-3: for considering the situation of wind energy turbine set Whole Response under wind speed change and electric network fault simultaneously, set up with wind energy turbine set input wind speed and wind farm grid-connected point voltage for input variable, with the non-mechanism equivalent model that wind energy turbine set overall power is output variable, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{P\_v}\left(V\right)f_{P\_u}\left(U\right)\\ f_{P\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{P\_u}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{Q\_v}\left(V\right)f_{Q\_u}\left(U\right)\\ f_{Q\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{Q\_u}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.---(*)$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U, H _pu(s), H _qu(s), g _pu(U), g _qu(U) definition and expression formula are with step A-1, V, H _pv(s), H _qv(s), g _pv(V), g _qv(V) definition and expression formula are with step A-2;

Step B: the parameter of the non-mechanism equivalent model of identification wind energy turbine set, concrete steps are as follows:

Ant colony optimization algorithm identification abovementioned steps A-1 is adopted to set up equivalent model H _pu(s) and H _qus () and steps A-2 set up equivalent model H _pv(s) and H _qvthe parameter of (s);

Step C, parameter identification result according to abovementioned steps B, (*) formula in abovementioned steps A-3 that substitutes into obtains the complete non-mechanism equivalent model of wind energy turbine set.

In further embodiment, the realization of abovementioned steps 2 specifically comprises the following steps:

Step B-1: the non-mechanism equivalent model that corresponding steps A-1 is set up, according to the delta data of wind farm grid-connected point voltage and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pu(s) and H _quthe parameter of (s);

Step B-2: the non-mechanism equivalent model that corresponding steps A-2 is set up, according to the wind energy turbine set input delta data of wind speed and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pv(s) and H _qvthe parameter of (s).

According to improvement of the present invention, also propose a kind of non-mechanism Equivalent Modeling device of wind energy turbine set, comprising:

First module, for according to different sights, sets up non-mechanism equivalent model;

Second unit, for according to the aforementioned first module of ant colony optimization algorithm identification set up the model parameter of non-mechanism equivalent model; And

Unit the 3rd, the model parameter for aforementioned identification being obtained is updated to aforementioned (*) formula and obtains the complete non-mechanism equivalent model of wind energy turbine set.

In further embodiment, described first module comprises:

First module, for the situation according to wind energy turbine set overall dynamics response under electric network fault, set up the first non-mechanism equivalent model of wind farm grid-connected point voltage and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pu}}\left(U\right)\\ f_{\mathrm{Pu}}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{\mathrm{Qu}}\left(U\right)\\ f_{\mathrm{Qu}}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U is wind farm grid-connected point voltage; H _pu(s) and H _qus () is respectively transfer function equation; g _puand g (U) _qu(U) algebraic equation that to be respectively with wind farm grid-connected point voltage U be independent variable; Transfer function equation H _pu(s) and H _qus () is expressed as follows respectively:

$\left\{\begin{array}{c}{H}_{\mathrm{Pu}}\left(s\right)=\frac{b_{\mathrm{pu}0}{s}^{\mathrm{pn}}+b_{\mathrm{pu}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pun}}}{s^{\mathrm{pm}}+a_{\mathrm{pu}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pum}}}\\ H_{\mathrm{Qu}}\left(s\right)=\frac{b_{\mathrm{qu}0}{s}^{\mathrm{qn}}+b_{\mathrm{qu}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qun}}}{s^{\mathrm{qm}}+a_{\mathrm{qu}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qum}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pu0, b _pu1..., b _pun, a _pu1..., a _pum, b _qu0, b _qu1..., b _qun, a _qu1..., a _qumfor constant coefficient; Algebraic equation g _puand g (U) _qu(U) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\\ g_{\mathrm{Qu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\end{array}\right.$

In formula, U, U ₀be respectively wind farm grid-connected point voltage and grid-connected point voltage steady-state value;

Second module, for the situation according to wind energy turbine set electrical power change under fluctuations in wind speed, set up the second non-mechanism equivalent model of wind energy turbine set input wind speed and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pv}}\left(V\right)\\ f_{\mathrm{Pv}}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ Q=Q_0{f}_{\mathrm{Qv}}\left(V\right)\\ f_{\mathrm{Qv}}\left(V\right)=H_{\mathrm{Qv}}\left(s\right)g_{\mathrm{Qv}}\left(V\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; V is wind energy turbine set input wind speed; H _pv(s) and H _qvs () is respectively transfer function equation; g _pvand g (V) _qv(V) algebraic equation that to be respectively with wind energy turbine set input wind speed V be independent variable, this transfer function equation H _pv(s) and H _qvs () expression formula is:

$\left\{\begin{array}{c}{H}_{\mathrm{Pv}}\left(s\right)=\frac{b_{\mathrm{pv}0}{s}^{\mathrm{pn}}+b_{\mathrm{pv}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pvn}}}{s^{\mathrm{pm}}+a_{\mathrm{pv}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pvm}}}\\ H_{\mathrm{Qv}}\left(s\right)=\frac{b_{\mathrm{qv}0}{s}^{\mathrm{qn}}+b_{\mathrm{qv}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qvn}}}{s^{\mathrm{qm}}+a_{\mathrm{qv}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qvm}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pv0, b _pv1..., b _pvn, a _pv1..., a _pvm, b _qv0, b _qv1..., b _qvn, a _qv1..., a _qvmfor constant coefficient;

Algebraic equation g _pvand g (V) _qv(V) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\\ g_{\mathrm{Qv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\end{array}\right.$

In formula, V is wind energy turbine set input wind speed; V ₀for the steady-state value of wind energy turbine set input wind speed;

3rd module, for according to considering the situation of wind energy turbine set Whole Response under wind speed change and electric network fault simultaneously, set up with wind energy turbine set input wind speed and wind farm grid-connected point voltage for input variable, the non-mechanism equivalent model being output variable with wind energy turbine set overall power, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{P\_v}\left(V\right)f_{P\_u}\left(U\right)\\ f_{P\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{P\_u}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{Q\_v}\left(V\right)f_{Q\_u}\left(U\right)\\ f_{Q\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{Q\_u}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.---(*)$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output.

In further embodiment, described second unit comprises:

First recognition module, for the model parameter according to the aforementioned first non-mechanism equivalent model of ant colony optimization algorithm identification;

Second recognition module, for the model parameter according to the aforementioned second non-mechanism equivalent model of ant colony optimization algorithm identification.

From the above technical solution of the present invention shows that, the non-mechanism equivalent modeling method of wind energy turbine set that the present invention proposes and device, overcome current wind energy turbine set measurement information not enough, cause the deficiency being difficult to set up wind energy turbine set mechanism equivalent model, the interaction between the inner each modular character of wind turbine generator and the inner each factor of wind energy turbine set is have ignored in the non-modelling by mechanism method of invention proposition, thus simplify the process of wind energy turbine set Equivalent Modeling and required data volume, agreed with again the final purpose of wind energy turbine set Equivalent Modeling well simultaneously, namely setting up can the simple lower-order model of accurate description wind energy turbine set external characteristic.The non-mechanism equivalent model of the wind energy turbine set set up by the inventive method can be widely used in the simulation analysis containing wind-powered electricity generation electric power system, has good future in engineering applications.

Accompanying drawing explanation

Fig. 1 is the illustrative diagram of wind farm grid-connected point voltage change curve.

The active power that Fig. 2 a-2b is wind energy turbine set respectively under the such as grid-connected point voltage change of Fig. 1 and reactive power response curve schematic diagram.

Fig. 3 is the wind energy turbine set non-mechanism equivalent model parameter identification schematic flow sheet based on non-modelling by mechanism method.

Fig. 4 is the illustrative diagram of wind farm wind velocity change curve.

The active power that Fig. 5 a-5b is wind energy turbine set respectively under such as Fig. 4 wind speed change and reactive power response curve schematic diagram.

Fig. 6 is the complete structure schematic diagram of the non-mechanism equivalent model of wind energy turbine set.

Fig. 7 a-7b is active power of wind power field under grid-connected point voltage and wind speed change simultaneously and reactive power dynamic response curve schematic diagram respectively.

Fig. 8 a-8b is the active power of wind power field of Fig. 7 a-7b and the partial enlarged drawing of reactive power dynamic response curve respectively.

Fig. 9 is the illustrative diagram of the non-mechanism Equivalent Modeling device of an embodiment of the present invention wind energy turbine set.

Embodiment

In order to more understand technology contents of the present invention, institute's accompanying drawings is coordinated to be described as follows especially exemplified by specific embodiment.

According to preferred embodiment of the present invention, a kind of non-mechanism equivalent modeling method of wind energy turbine set, use the wind farm grid-connected voltage of point and the wind speed of anemometer tower to be input variable, use the response of wind energy turbine set overall power to be output variable, set up the non-mechanism equivalent model of wind energy turbine set based on transfer function by discrimination method, this non-mechanism equivalent modeling method comprises the following steps:

Steps A: the equation setting up the non-mechanism equivalent model of wind energy turbine set, concrete steps are as follows:

Steps A-1: for the situation of wind energy turbine set overall dynamics response under electric network fault, now ignore the change of wind speed, set up the non-mechanism equivalent model of wind farm grid-connected point voltage and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pu}}\left(U\right)\\ f_{\mathrm{Pu}}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{\mathrm{Qu}}\left(U\right)\\ f_{\mathrm{Qu}}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U is wind farm grid-connected point voltage; H _pu(s) and H _qus () is respectively transfer function equation; g _puand g (U) _qu(U) algebraic equation that to be respectively with wind farm grid-connected point voltage U be independent variable; Transfer function equation H _pu(s) and H _qus () is expressed as follows respectively:

$\left\{\begin{array}{c}{H}_{\mathrm{Pu}}\left(s\right)=\frac{b_{\mathrm{pu}0}{s}^{\mathrm{pn}}+b_{\mathrm{pu}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pun}}}{s^{\mathrm{pm}}+a_{\mathrm{pu}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pum}}}\\ H_{\mathrm{Qu}}\left(s\right)=\frac{b_{\mathrm{qu}0}{s}^{\mathrm{qn}}+b_{\mathrm{qu}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qun}}}{s^{\mathrm{qm}}+a_{\mathrm{qu}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qum}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pu0, b _pu1..., b _pun, a _pu1..., a _pum, b _qu0, b _qu1..., b _qun, a _qu1..., a _qumfor constant coefficient; Algebraic equation g _puand g (U) _qu(U) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\\ g_{\mathrm{Qu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\end{array}\right.$

In formula, U, U ₀be respectively wind farm grid-connected point voltage and grid-connected point voltage steady-state value;

Steps A-2: for the situation of wind energy turbine set electrical power change under fluctuations in wind speed, set up the non-mechanism equivalent model of wind energy turbine set input wind speed (i.e. the wind speed of anemometer tower) and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pv}}\left(V\right)\\ f_{\mathrm{Pv}}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ Q=Q_0{f}_{\mathrm{Qv}}\left(V\right)\\ f_{\mathrm{Qv}}\left(V\right)=H_{\mathrm{Qv}}\left(s\right)g_{\mathrm{Qv}}\left(V\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; V is wind energy turbine set input wind speed; H _pv(s) and H _qvs () is respectively transfer function equation; g _pvand g (V) _qv(V) algebraic equation that to be respectively with wind energy turbine set input wind speed V be independent variable, this transfer function equation H _pv(s) and H _qvs () expression formula is:

$\left\{\begin{array}{c}{H}_{\mathrm{Pv}}\left(s\right)=\frac{b_{\mathrm{pv}0}{s}^{\mathrm{pn}}+b_{\mathrm{pv}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pvn}}}{s^{\mathrm{pm}}+a_{\mathrm{pv}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pvm}}}\\ H_{\mathrm{Qv}}\left(s\right)=\frac{b_{\mathrm{qv}0}{s}^{\mathrm{qn}}+b_{\mathrm{qv}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qvn}}}{s^{\mathrm{qm}}+a_{\mathrm{qv}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qvm}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pv0, b _pv1..., b _pvn, a _pv1..., a _pvm, b _qv0, b _qv1..., b _qvn, a _qv1..., a _qvmfor constant coefficient;

Algebraic equation g _pvand g (V) _qv(V) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\\ g_{\mathrm{Qv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\end{array}\right.$

In formula, V is wind energy turbine set input wind speed; V ₀for the steady-state value of wind energy turbine set input wind speed;

Steps A-3: for considering the situation of wind energy turbine set Whole Response under wind speed change and electric network fault simultaneously, set up with wind energy turbine set input wind speed and wind farm grid-connected point voltage for input variable, with the non-mechanism equivalent model that wind energy turbine set overall power is output variable, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{P\_v}\left(V\right)f_{P\_u}\left(U\right)\\ f_{P\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{P\_u}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{Q\_v}\left(V\right)f_{Q\_u}\left(U\right)\\ f_{Q\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{Q\_u}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.---(*)$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U, H _pu(s), H _qu(s), g _pu(U), g _qu(U) definition and expression formula are with step A-1, V, H _pv(s), H _qv(s), g _pv(V), g _qv(V) definition and expression formula are with step A-2;

Step B: the parameter of the non-mechanism equivalent model of identification wind energy turbine set, concrete steps are as follows:

Ant colony optimization algorithm identification abovementioned steps A-1 is adopted to set up equivalent model H _pu(s) and H _qus () and steps A-2 set up equivalent model H _pv(s) and H _qvthe parameter of (s);

Step C, parameter identification result according to abovementioned steps B, (*) formula in abovementioned steps A-3 that substitutes into obtains the complete non-mechanism equivalent model of wind energy turbine set.

Equivalent modeling method in the present embodiment, establishes the equivalent model H of relation between wind farm grid-connected point voltage and Power Output for Wind Power Field with the form of transfer function _us () (comprises H _pu(s) and H _qu(s)), and wind energy turbine set inputs the equivalent model H of relation between wind speed and Power Output for Wind Power Field _vs () (comprises H _pv(s) and H _qv(s)), as a complete non-mechanism equivalent model (situation for studying wind energy turbine set Whole Response under wind speed change and electric network fault simultaneously), transfer function is wherein H _u(s) and H _vthe product of (s).

The method of the present embodiment, composition graphs 1 and above-mentioned steps, first the non-mechanism equivalent model of wind energy turbine set based on single input variable is set up, then based on decoupling zero modeling method, set up the non-mechanism equivalent model of wind energy turbine set of multivariable (wind energy turbine set input wind speed and Power Output for Wind Power Field), final what set up is a complete non-mechanism model of wind energy turbine set.

As preferred mode, in the present embodiment, the realization of abovementioned steps 2 specifically comprises the following steps:

Step B-1: the non-mechanism equivalent model that corresponding steps A-1 is set up, according to the delta data of wind farm grid-connected point voltage and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pu(s) and H _quthe parameter of (s);

Step B-2: the non-mechanism equivalent model that corresponding steps A-2 is set up, according to the wind energy turbine set input delta data of wind speed and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pv(s) and H _qvthe parameter of (s).

Shown in Fig. 1-Fig. 8 a, 8b, the enforcement of the non-mechanism equivalent modeling method of aforementioned wind energy turbine set is exemplarily described.

The equation of the non-mechanism equivalent model of wind energy turbine set is arranged in the steps A-1 of steps A, steps A-2 and steps A-3, repeats no more herein.

For the step B of identification model parameter, due to when Wind turbines terminal voltage and input wind speed change simultaneously, Wind turbines internal control system can make corresponding actions.Under normal circumstances, during the system failure Wind turbines dynamically want far faster than wind speed change cause dynamic.Then can think, when wind farm grid-connected point voltage and wind speed change, the change of the main response voltage of blower interior control system, because process time is short, can ignore the impact of wind speed change; During wind speed change, blower fan is in quasi-stable state process, control system can keep set end voltage substantially constant, can carry out decoupled identification to wind energy turbine set non-mechanism equivalent model under grid-connected point voltage change and wind speed change two kinds of situations thus, the final parameter identification result that merges can obtain the non-mechanism equivalent model of complete wind energy turbine set.

Step B-1: the non-mechanism equivalent model that corresponding steps A-1 is set up, according to the delta data of wind farm grid-connected point voltage and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pu(s) and H _quthe parameter of (s).

First, the wind energy turbine set of collection wind energy turbine set under grid-connected point voltage changes (example that Fig. 1 exemplarily gives wind farm grid-connected point voltage change curve) situation is meritorious, idle dynamic response curve (as shown in the solid line in Fig. 2 a, 2b), then parameter identification flow process is as shown in Figure 3 adopted, identification transfer function H _pu(s) and H _qus the parameter of (), the result obtained is as shown in table 1.The dynamic response of the non-mechanism equivalent model of the wind energy turbine set obtained is emulated as shown in phantom in Figure 2 by this identification result.

Parameter identification result under the grid-connected point voltage situation of change of table 1 and model relative error

Equivalent model	b 3	b 2	b 1	b 0	a 2	a 1	a 0	Model error
									P-U	-1.5067	4.3238	-5.7532	2.6132	36.01	40.81	2.6132	0.14％
Q-U	-6.8280	175.68	4.2096	1.3083	4.9014	0.6332	1.3083	2.00％

Step B-2: the non-mechanism equivalent model that corresponding steps A-2 is set up, according to the wind energy turbine set input delta data of wind speed and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pv(s) and H _qvthe parameter of (s).

First, the wind energy turbine set of collection wind energy turbine set in input wind speed change (Fig. 4 exemplarily gives an example of wind farm wind velocity change) situation is meritorious, idle dynamic response curve (as shown in the solid line in Fig. 5 a, 5b), then parameter identification flow process is as shown in Figure 3 adopted, identification transfer function H _pu(s) and H _qus the parameter of (), the result obtained is as shown in table 2.The dynamic response of the non-mechanism equivalent model of the wind energy turbine set obtained is emulated as shown in phantom in Figure 5 by this identification result.

Parameter identification result under table 2 wind speed situation of change and model relative error

Equivalent model	b 2	b 1	b 0	a 2	a 1	a 0	Relative error
								P-U	1.6826	4.1659	0.3545	19.50	10.46	0.3545	0.81％
Q-U	7.4130	2.0827	0.0301	16.10	0.7027	0.0301	3.71％

Step B-3: according to the parameter identification result of step B-1 and step B-2, obtains the complete non-mechanism equivalent model of wind energy turbine set as shown in (*) formula in steps A-3.

Step B-1 and step B-2 identification are obtained model coefficient and substitute into (*) formula, the non-mechanism equivalent model of wind energy turbine set as shown in Figure 6 can be obtained, the H in figure _us () comprises H _pu(s) and H _qu(s), the H in figure _vs () comprises H _pv(s) and H _qvs (), exporting y can be active-power P, also can be reactive power Q.

For Fig. 6 institute representation model, from transfer function basic theories, as function H _u(s) and H _vconstant coefficient a in (s) ₀and b ₀time equal, when input variable keeps constant, transfer function H _u(s) and H _vs (), without dynamic process, is maintained constant coefficient 1.Therefore, when inputting wind farm grid-connected point voltage u and constant stable state wind speed v, (v/v ₀) ³h _v(s)=1, without the dynamic response of wind speed, equivalent model only shows the dynamic characteristic of wind energy turbine set under change in voltage; When inputting the wind speed v of change and constant wind farm grid-connected point voltage u, (u/u ₀) ²h _v(s)=1, no-voltage dynamic response, equivalent model only shows the dynamic behaviour of the lower wind energy turbine set of wind speed change; And when wind farm grid-connected point voltage u and input wind speed v changes simultaneously, transfer function H _u(s) and H _vs the change of () difference response voltage u and wind speed v, it dynamically shows as H _u(s) and H _vthe product of (s), the Nonlinear Superposition of two kinds of dynamic characteristics that what namely model was complete is dynamically.

The ant colony optimization algorithm adopted in abovementioned steps adopts conventional steps of the prior art to realize.Particularly, in the present embodiment, this algorithm is utilized to carry out the step of parameter identification as follows:

In wind energy turbine set transfer function equivalent model, error target function is set up to be:

$E=\underset{k=1}{\overset{N}{\mathrm{\Σ}}}[Y\left(k\right)-Y_m\left(k\right){]}^T[Y\left(k\right)-Y_m\left(k\right)]$

In formula: N is sampling number; Y _mfor the meritorious of actual measurement or without work value (i.e. actual value); Y is calculate meritorious of equivalent model or without work value.

The wind energy turbine set equivalent model based on single input variable of aforementioned proposition is continuous time model, when carrying out parameter identification to it, because the inputoutput data gathered, process is discrete data, equivalent model need be done sliding-model control.

In order to improve equivalent model computational accuracy, the present embodiment adopts numerical integration algorithm to solve equivalent model.Numerical Integral Solving Method essence is, with certain discretization method, Continuous Variable Problems is converted into discrete variable problem.For in the non-mechanism equivalent model of wind energy turbine set voltage active power, model expression is:

$P=P_0\frac{b_{\mathrm{pu}0}{s}^{\mathrm{pn}}+b_{\mathrm{pu}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pun}}}{s^{\mathrm{pm}}+a_{\mathrm{pu}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pum}}}{\left(\frac{U}{U_0}\right)}^2---(3-12)$

Transfer function model is converted into state space equation form:

$\left\{\begin{array}{c}\stackrel{\·}{x}=\mathrm{Ax}+B{(U/U_0)}^2\\ P=\mathrm{Cx}+D{P}_0{(U/U_0)}^2\end{array}\right.---(3-13)$

In formula,

$A=\left[\begin{array}{ccccc}0& 1& ...& 0& 0\\ 0& 0& ...& 0& 0\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ 0& 0& ...& 0& 1\\ {-a}_{\mathrm{pu}1}& {-a}_{\mathrm{pu}2}& ...& {-a}_{\mathrm{pum}-1}& {-a}_{\mathrm{pum}}\end{array}\right]$

x＝[x ₁,x ₂,…,x _m] ^T

B＝[00…1] ^T，C＝[b _pu1-b _pu0a _pu1b _pu2-b _pu0a _pu2…b _pun-b _pu0a _pum]，D＝b _pu0

After providing equivalent model parameter and input variable data, adopt fourth-order Runge-Kutta method to carry out simulation calculation to state space equation, obtain equivalent model and export centrifugal pump Y.Fourth-order Runge-Kutta method in this example is ripe, does not repeat them here.

After obtaining equivalent model output centrifugal pump, the error target function value of corresponding equivalent model can be obtained.

Error target function E is the complicated function of parameter θ, is difficult to obtain analytic solutions, and solution space often exists multiple extreme value.Therefore, optimized algorithm must be very effective.The present embodiment adopts ant colony optimization algorithm.

In order to determine the order of equivalent model, when carrying out parameter identification to model, increase successively model order n (even n=1,2 ...), calculate its optimized parameter θ respectively _*and target function value E (n), when model order n increases, E (n) declines, and when n increases to the order needed for realistic model, E (n) starts to decline not obvious, then order corresponding when selecting E (n) to decline not obvious is as the best order of equivalent model.Equivalent model order is chosen as 1 ~ 5 rank, successively identification model parameter, and it is as shown in the table to calculate errors target function value under each order.

The comparison of error target function value under different order

Equivalent model order	Equivalent model number of parameters	P-U model E value	Q-U model E value
				1	3	4.5906	23.86
2	5	2.3432	5.4816
				3	7	0.7124	1.8621
4	9	0.7040	1.8206

5

11

0.7032

1.8181

As can be seen from table, when equivalent model order is less than 3, corresponding error target function value is larger; And after model order is greater than 3, the increase of order is little on the impact of E value, but time adding parameter identification of parameter to be identified and difficulty in equivalent model, therefore the best order of preference pattern is 3 rank.And then identification obtains wind energy turbine set non-mechanism equivalent model parameter.

For verifying the validity of this non-mechanism equivalent model, wind speed shown in point voltage change grid-connected shown in Fig. 1 and Fig. 4 is changed these two kinds of inputs superpose, the output of the non-mechanism equivalent model of wind energy turbine set is as shown in Fig. 7 a, 7b, the partial enlarged drawing in grid-connected point voltage change moment as shown in Fig. 8 a, 8b, and wind energy turbine set detailed model export between error as follows listed by.

	Global error (%)	The local error (%) of change in voltage
			Active power	0.82	0.59
Reactive power	4.20	3.13

Can find out, when inputting wind speed and changing separately, wind energy turbine set non-mechanism equivalent model has good dynamic accuracy, and meritorious, reactive power error is respectively 0.82% and 4.20%; In the local that the change of grid-connected point voltage and input wind speed change simultaneously, meritorious, idle error is respectively 0.59% and 3.13%, can meet error requirements completely, indicate previous embodiment of the present invention put forward the validity of the non-mechanism equivalent model of wind energy turbine set.

According to of the present invention open, as shown in Figure 9, a kind of non-mechanism Equivalent Modeling device of wind energy turbine set, comprising:

First module 100, for according to different sights, sets up non-mechanism equivalent model;

Second unit 200, for according to the aforementioned first module of ant colony optimization algorithm identification set up the model parameter of non-mechanism equivalent model; And

3rd unit 300, the model parameter for aforementioned identification being obtained is updated to aforementioned (*) formula and obtains the complete non-mechanism equivalent model of wind energy turbine set.

As shown in Figure 9, first module 100 comprises the first module 101, second module 102, the 3rd module 103, wherein:

First module 101, for the situation according to wind energy turbine set overall dynamics response under electric network fault, set up the first non-mechanism equivalent model of wind farm grid-connected point voltage and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pu}}\left(U\right)\\ f_{\mathrm{Pu}}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{\mathrm{Qu}}\left(U\right)\\ f_{\mathrm{Qu}}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U is wind farm grid-connected point voltage; H _pu(s) and H _qus () is respectively transfer function equation; g _puand g (U) _qu(U) algebraic equation that to be respectively with wind farm grid-connected point voltage U be independent variable; Transfer function equation H _pu(s) and H _qus () is expressed as follows respectively:

$\left\{\begin{array}{c}{H}_{\mathrm{Pu}}\left(s\right)=\frac{b_{\mathrm{pu}0}{s}^{\mathrm{pn}}+b_{\mathrm{pu}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pun}}}{s^{\mathrm{pm}}+a_{\mathrm{pu}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pum}}}\\ H_{\mathrm{Qu}}\left(s\right)=\frac{b_{\mathrm{qu}0}{s}^{\mathrm{qn}}+b_{\mathrm{qu}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qun}}}{s^{\mathrm{qm}}+a_{\mathrm{qu}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qum}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pu0, b _pu1..., b _pun, a _pu1..., a _pum, b _qu0, b _qu1..., b _qun, a _qu1..., a _qumfor constant coefficient; Algebraic equation g _puand g (U) _qu(U) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\\ g_{\mathrm{Qu}}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\end{array}\right.$

In formula, U, U ₀be respectively wind farm grid-connected point voltage and grid-connected point voltage steady-state value;

Second module 102, for the situation according to wind energy turbine set electrical power change under fluctuations in wind speed, set up the second non-mechanism equivalent model of wind energy turbine set input wind speed and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{\mathrm{Pv}}\left(V\right)\\ f_{\mathrm{Pv}}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ Q=Q_0{f}_{\mathrm{Qv}}\left(V\right)\\ f_{\mathrm{Qv}}\left(V\right)=H_{\mathrm{Qv}}\left(s\right)g_{\mathrm{Qv}}\left(V\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; V is wind energy turbine set input wind speed; H _pv(s) and H _qvs () is respectively transfer function equation; g _pvand g (V) _qv(V) algebraic equation that to be respectively with wind energy turbine set input wind speed V be independent variable, this transfer function equation H _pv(s) and H _qvs () expression formula is:

$\left\{\begin{array}{c}{H}_{\mathrm{Pv}}\left(s\right)=\frac{b_{\mathrm{pv}0}{s}^{\mathrm{pn}}+b_{\mathrm{pv}1}{s}^{\mathrm{pn}-1}+...+b_{\mathrm{pvn}}}{s^{\mathrm{pm}}+a_{\mathrm{pv}1}{s}^{\mathrm{pm}-1}+...+a_{\mathrm{pvm}}}\\ H_{\mathrm{Qv}}\left(s\right)=\frac{b_{\mathrm{qv}0}{s}^{\mathrm{qn}}+b_{\mathrm{qv}1}{s}^{\mathrm{qn}-1}+...+b_{\mathrm{qvn}}}{s^{\mathrm{qm}}+a_{\mathrm{qv}1}{s}^{\mathrm{qm}-1}+...+a_{\mathrm{qvm}}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pv0, b _pv1..., b _pvn, a _pv1..., a _pvm, b _qv0, b _qv1..., b _qvn, a _qv1..., a _qvmfor constant coefficient;

Algebraic equation g _pvand g (V) _qv(V) expression formula is:

$\left\{\begin{array}{c}{g}_{\mathrm{Pv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\\ g_{\mathrm{Qv}}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\end{array}\right.$

In formula, V is wind energy turbine set input wind speed; V ₀for the steady-state value of wind energy turbine set input wind speed;

3rd module 103, for according to considering the situation of wind energy turbine set Whole Response under wind speed change and electric network fault simultaneously, set up with wind energy turbine set input wind speed and wind farm grid-connected point voltage for input variable, the non-mechanism equivalent model being output variable with wind energy turbine set overall power, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{P\_v}\left(V\right)f_{P\_u}\left(U\right)\\ f_{P\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{P\_u}\left(U\right)=H_{\mathrm{Pu}}\left(s\right)g_{\mathrm{Pu}}\left(U\right)\\ Q=Q_0{f}_{Q\_v}\left(V\right)f_{Q\_u}\left(U\right)\\ f_{Q\_v}\left(V\right)=H_{\mathrm{Pv}}\left(s\right)g_{\mathrm{Pv}}\left(V\right)\\ f_{Q\_u}\left(U\right)=H_{\mathrm{Qu}}\left(s\right)g_{\mathrm{Qu}}\left(U\right)\end{array}\right.---(*)$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output.

As shown in Figure 9, second unit 200 comprises the first recognition module 201, second recognition module 202, wherein:

First recognition module 201, for the model parameter according to the aforementioned first non-mechanism equivalent model of ant colony optimization algorithm identification;

Second recognition module 202, for the model parameter according to the aforementioned second non-mechanism equivalent model of ant colony optimization algorithm identification.

With reference to the implementation process of the non-mechanism equivalent modeling method of aforementioned wind energy turbine set, the non-mechanism Equivalent Modeling device of wind energy turbine set of the present disclosure, the function of its first module 100, second unit 200, the 3rd unit 300 and the first module 101, second module 102, the 3rd module 103, first recognition module 201, second recognition module 202, effect detailed description in the above-described embodiments, does not repeat them here.

According to record of the present disclosure, those of ordinary skill in the art can go out corresponding program according to its content development and run in a computer system, thus realizes every step and the effect of the embodiment of aforementioned non-mechanism equivalent modeling method.

Although the present invention with preferred embodiment disclose as above, so itself and be not used to limit the present invention.Persond having ordinary knowledge in the technical field of the present invention, without departing from the spirit and scope of the present invention, when being used for a variety of modifications and variations.Therefore, protection scope of the present invention is when being as the criterion depending on those as defined in claim.

Claims (3)

1. the non-mechanism equivalent modeling method of a wind energy turbine set, it is characterized in that, use the wind farm grid-connected voltage of point and the wind speed of anemometer tower to be input variable, use the response of wind energy turbine set overall power to be output variable, set up the non-mechanism equivalent model of wind energy turbine set based on transfer function by discrimination method, this non-mechanism equivalent modeling method comprises the following steps:

Steps A: the equation setting up the non-mechanism equivalent model of wind energy turbine set, concrete steps are as follows:

Steps A-1: for the situation of wind energy turbine set overall dynamics response under electric network fault, now ignore the change of wind speed, set up the non-mechanism equivalent model of wind farm grid-connected point voltage and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{Pu}\left(U\right)\\ f_{Pu}\left(U\right)=H_{Pu}\left(s\right)g_{Pu}\left(U\right)\\ Q=Q_0{f}_{Qu}\left(U\right)\\ f_{Qu}\left(U\right)=H_{Qu}\left(s\right)g_{Qu}\left(U\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U is wind farm grid-connected point voltage; H _pu(s) and H _qus () is respectively transfer function equation; g _puand g (U) _qu(U) algebraic equation that to be respectively with wind farm grid-connected point voltage U be independent variable; Transfer function equation H _pu(s) and H _qus () is expressed as follows respectively:

$\left\{\begin{array}{c}{H}_{Pu}\left(s\right)=\frac{b_{pu0}{s}^{pn}+b_{pu1}{s}^{pn-1}+...+b_{pun}}{s^{pm}+a_{pu1}{s}^{pm-1}+...+a_{pum}}\\ H_{Qu}\left(s\right)=\frac{b_{qu0}{s}^{qn}+b_{qu1}{s}^{qn-1}+...+b_{qun}}{s^{qm}+a_{qu1}{s}^{qm-1}+...+a_{qum}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pu0, b _pu1..., b _pun, a _pu1..., a _pum, b _qu0, b _qu1..., b _qun, a _qu1..., a _qumfor constant coefficient; Algebraic equation g _puand g (U) _qu(U) expression formula is:

$\left\{\begin{array}{c}{g}_{Pu}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\\ g_{Qu}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\end{array}\right.$

In formula, U, U ₀be respectively wind farm grid-connected point voltage and grid-connected point voltage steady-state value;

Steps A-2: for the situation of wind energy turbine set electrical power change under fluctuations in wind speed, set up the non-mechanism equivalent model of wind energy turbine set input wind speed and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{Pv}\left(V\right)\\ f_{Pv}\left(V\right)=H_{Pv}\left(s\right)g_{Pv}\left(V\right)\\ Q=Q_0{f}_{Qv}\left(V\right)\\ f_{Qv}\left(V\right)=H_{Qv}\left(s\right)g_{Qv}\left(V\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; V is wind energy turbine set input wind speed; H _pv(s) and H _qvs () is respectively transfer function equation; g _pvand g (V) _qv(V) algebraic equation that to be respectively with wind energy turbine set input wind speed V be independent variable, this transfer function equation H _pv(s) and H _qvs () expression formula is:

$\left\{\begin{array}{c}{H}_{Pv}\left(s\right)=\frac{b_{pv0}{s}^{pn}+b_{pv1}{s}^{pn-1}+...+b_{pvn}}{s^{pm}+a_{pv1}{s}^{pm-1}+...+a_{pvm}}\\ H_{Qv}\left(s\right)=\frac{b_{qv0}{s}^{qn}+b_{qv1}{s}^{qn-1}+...+b_{qvn}}{s^{qm}+a_{qv1}{s}^{qm-1}+...+a_{qvm}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pv0, b _pv1..., b _pvn, a _pv1..., a _pvm, b _qv0, b _qv1..., b _qvn, a _qv1..., a _qvmfor constant coefficient;

Algebraic equation g _pvand g (V) _qv(V) expression formula is:

$\left\{\begin{array}{c}{g}_{Pv}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\\ g_{Qv}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\end{array}\right.$

In formula, V is wind energy turbine set input wind speed; V ₀for the steady-state value of wind energy turbine set input wind speed;

Steps A-3: for considering the situation of wind energy turbine set Whole Response under wind speed change and electric network fault simultaneously, set up with wind energy turbine set input wind speed and wind farm grid-connected point voltage for input variable, with the non-mechanism equivalent model that wind energy turbine set overall power is output variable, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{P\_v}\left(V\right)f_{P\_u}\left(U\right)\\ f_{P\_v}\left(V\right)=H_{Pv}\left(s\right)g_{Pv}\left(V\right)\\ f_{P\_u}\left(U\right)=H_{Pu}\left(s\right)g_{Pu}\left(U\right)\\ Q=Q_0{f}_{Q\_v}\left(V\right)f_{Q\_u}\left(U\right)\\ f_{Q\_v}\left(V\right)=H_{Pv}\left(s\right)g_{Pv}\left(V\right)\\ f_{Q\_u}\left(U\right)=H_{Qu}\left(s\right)g_{Qu}\left(U\right)\end{array}\right.---\left(1\right)$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U, H _pu(s), H _qu(s), g _pu(U), g _qu(U) definition and expression formula are with step A-1, V, H _pv(s), H _qv(s), g _pv(V), g _qv(V) definition and expression formula are with step A-2;

Step B: the parameter of the non-mechanism equivalent model of identification wind energy turbine set, concrete steps are as follows:

Adopt the transfer function equation H ignoring wind speed change that ant colony optimization algorithm identification abovementioned steps A-1 sets up _pu(s) and H _quthe transfer function equation H of s consideration fluctuations in wind speed that () and steps A-2 is set up _pv(s) and H _qvthe parameter of (s);

Step C, parameter identification result according to abovementioned steps B, (1) formula in abovementioned steps A-3 that substitutes into obtains the complete non-mechanism equivalent model of wind energy turbine set.

2. the non-mechanism equivalent modeling method of wind energy turbine set according to claim 1, it is characterized in that, the realization of abovementioned steps B specifically comprises the following steps:

Step B-1: the non-mechanism equivalent model that corresponding steps A-1 is set up, according to the delta data of wind farm grid-connected point voltage and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pu(s) and H _quthe parameter of (s);

Step B-2: the non-mechanism equivalent model that corresponding steps A-2 is set up, according to the wind energy turbine set input delta data of wind speed and the power response data of wind energy turbine set, adopts ant colony optimization algorithm identification H _pv(s) and H _qvthe parameter of (s).

3. a non-mechanism Equivalent Modeling device for wind energy turbine set, is characterized in that, comprising:

First module, for according to different sights, sets up non-mechanism equivalent model;

Second unit, for according to the aforementioned first module of ant colony optimization algorithm identification set up the model parameter of non-mechanism equivalent model; And

Unit the 3rd, the model parameter for aforementioned identification being obtained is updated in the expression formula of the non-mechanism equivalent model that first module is set up and obtains the complete non-mechanism equivalent model of wind energy turbine set;

Wherein:

Described first module comprises:

First module, for the situation according to wind energy turbine set overall dynamics response under electric network fault, set up the first non-mechanism equivalent model of wind farm grid-connected point voltage and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{Pu}\left(U\right)\\ f_{Pu}\left(U\right)=H_{Pu}\left(s\right)g_{Pu}\left(U\right)\\ Q=Q_0{f}_{Qu}\left(U\right)\\ f_{Qu}\left(U\right)=H_{Qu}\left(s\right)g_{Qu}\left(U\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; U is wind farm grid-connected point voltage; H _pu(s) and H _qus () is respectively transfer function equation; g _puand g (U) _qu(U) algebraic equation that to be respectively with wind farm grid-connected point voltage U be independent variable; Transfer function equation H _pu(s) and H _qus () is expressed as follows respectively:

$\left\{\begin{array}{c}{H}_{Pu}\left(s\right)=\frac{b_{pu0}{s}^{pn}+b_{pu1}{s}^{pn-1}+...+b_{pun}}{s^{pm}+a_{pu1}{s}^{pm-1}+...+a_{pum}}\\ H_{Qu}\left(s\right)=\frac{b_{qu0}{s}^{qn}+b_{qu1}{s}^{qn-1}+...+b_{qun}}{s^{qm}+a_{qu1}{s}^{qm-1}+...+a_{qum}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pu0, b _pu1..., b _pun, a _pu1..., a _pum, b _qu0, b _qu1..., b _qun, a _qu1..., a _qumfor constant coefficient; Algebraic equation g _puand g (U) _qu(U) expression formula is:

$\left\{\begin{array}{c}{g}_{Pu}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\\ g_{Qu}\left(U\right)={\left(\frac{U}{U_0}\right)}^2\end{array}\right.$

In formula, U, U ₀be respectively wind farm grid-connected point voltage and grid-connected point voltage steady-state value;

Second module, for the situation according to wind energy turbine set electrical power change under fluctuations in wind speed, set up the second non-mechanism equivalent model of wind energy turbine set input wind speed and power output, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{Pv}\left(V\right)\\ f_{Pv}\left(V\right)=H_{Pv}\left(s\right)g_{Pv}\left(V\right)\\ Q=Q_0{f}_{Qv}\left(V\right)\\ f_{Qv}\left(V\right)=H_{Qv}\left(s\right)g_{Qv}\left(V\right)\end{array}\right.$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output; V is wind energy turbine set input wind speed; H _pv(s) and H _qvs () is respectively transfer function equation; g _pvand g (V) _qv(V) algebraic equation that to be respectively with wind energy turbine set input wind speed V be independent variable, this transfer function equation H _pv(s) and H _qvs () expression formula is:

$\left\{\begin{array}{c}{H}_{Pv}\left(s\right)=\frac{b_{pv0}{s}^{pn}+b_{pv1}{s}^{pn-1}+...+b_{pvn}}{s^{pm}+a_{pv1}{s}^{pm-1}+...+a_{pvm}}\\ H_{Qv}\left(s\right)=\frac{b_{qv0}{s}^{qn}+b_{qv1}{s}^{qn-1}+...+b_{qvn}}{s^{qm}+a_{qv1}{s}^{qm-1}+...+a_{qvm}}\end{array}\right.$

In formula, pn≤pm, qn≤qm; b _pv0, b _pv1..., b _pvn, a _pv1..., a _pvm, b _qv0, b _qv1..., b _qvn, a _qv1..., a _qvmfor constant coefficient;

Algebraic equation g _pvand g (V) _qv(V) expression formula is:

$\left\{\begin{array}{c}{g}_{Pv}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\\ g_{Qv}\left(V\right)={\left(\frac{V}{V_0}\right)}^3\end{array}\right.$

In formula, V is wind energy turbine set input wind speed; V ₀for the steady-state value of wind energy turbine set input wind speed;

3rd module, for according to considering the situation of wind energy turbine set Whole Response under wind speed change and electric network fault simultaneously, set up with wind energy turbine set input wind speed and wind farm grid-connected point voltage for input variable, the non-mechanism equivalent model being output variable with wind energy turbine set overall power, its equation is:

$\left\{\begin{array}{c}P=P_0{f}_{P\_v}\left(V\right)f_{P\_u}\left(U\right)\\ f_{P\_v}\left(V\right)=H_{Pv}\left(s\right)g_{Pv}\left(V\right)\\ f_{P\_u}\left(U\right)=H_{Pu}\left(s\right)g_{Pu}\left(U\right)\\ Q=Q_0{f}_{Q\_v}\left(V\right)f_{Q\_u}\left(U\right)\\ f_{Q\_v}\left(V\right)=H_{Pv}\left(s\right)g_{Pv}\left(V\right)\\ f_{Q\_u}\left(U\right)=H_{Qu}\left(s\right)g_{Qu}\left(U\right)\end{array}\right.---(*)$

In formula, P, P ₀be respectively active power and the active power steady-state value of wind energy turbine set output; Q, Q ₀be respectively reactive power and the reactive power steady-state value of now wind energy turbine set output;

Described second unit comprises:

First recognition module, for the model parameter according to the aforementioned first non-mechanism equivalent model of ant colony optimization algorithm identification;

Second recognition module, for the model parameter according to the aforementioned second non-mechanism equivalent model of ant colony optimization algorithm identification.