Open access peer-reviewed chapter

# Wind Farms as Negative Loads and as Conventional Synchronous Generation – Modelling and Control

By Roberto Daniel Fernández, Pedro Eugenio Battaiotto and Ricardo Julián Mantz

Submitted: February 21st 2012Reviewed: October 3rd 2012Published: March 20th 2013

DOI: 10.5772/54089

## 1. Introduction

The concept of negative load [1,2] has been applied to wind generators to indicate their capability for delivering current meanwhile their voltage is imposed by the electrical system at the connection point. More recently, the same concept was applied for studying dispatch or spinning reserve considering that the total regulating power required at any moment depends on the sum the system load and the wind power which can counterbalance or increase load variations. In this way, aggregated variations must be investigated regarding wind power as negative load [3,4].

Traditionally, induction generators, squirrel cage and double fed (wound rotor) induction generators (DFIG) have been considered as current sources (or power ones) in power system analysis. Indeed, for analyzing power system stability in a linear frame, i.e. by small signal analysis, it is possible to find the power system eigenvalues and concluding about stability following the next steps:

1. Writing differential algebraic non linear model of the power system considering wind turbines wind induction machines as negative loads.

2. Linearizing the non linear model [2,5,6].

Checking the movement of the system eigenvalues when wind fixed speed generation is increased [7] or when different control strategies for active and reactive powers are applied to DFIGs wind farm [8,9].

Even when modelling wind generators or wind farms as current sources have shown through linear and non linear analysis [10] that wind farms can contribute to the power systems stability, it is important to consider that power systems have been developed from voltage synchronous generation, i.e. they (power systems) have not been developed with variable current sources as wind farms. For this reason, it is usually assumed that, as cited in [11], 'These plants (wind farm ones) exhibit static and dynamic characteristics that differ fundamentally from that of conventional generators. As a result, wind power plants do not fit the template for models of conventional generating facilities.'

However, considering wind farms as (variable) voltage sources could help not just for a better understanding of induction machines in an electrical grid but also for mimicking conventional voltage source behaviors with wind power plants.

In this way, this chapter developes the equivalent model of induction generators representing them as a voltage sources with a series impedance. As a consequence, aside from the variability of wind, it would allow to analyze wind energy generation as another voltage source in power systems and, then, it should be possible to introduce ''standard'' rules for conventional generation to non conventional ones.

The structure of this chapter is as follows. Firstly, the concept of ''general reference frame'' is introduced in order to analyze squirrel cage and double fed induction machines. The dynamics model of the induction machine is deduced by considering a common simplification about (fast) stator dynamics. Then, equivalent Thevenin models are presented according to internal voltages sources and as a which are functions of the active and reactive powers. Secondly, it is shown that it is possible controlling active and reactive wind generator powers. Finally, a Thevenin agregated model of the wind farm is proposed and some stability concerns of power systems are considered. In this way, Lyapunov theory is applied looking for demonstrating wind farms contribution to the power system stability considering wind farms as currents sources but also as voltage ones. In this last regard some wind farm control rules are derived from exploiting similarities between conventional generation and wind farms.

## 2. Induction machine equations in different reference frames

This section relies in [12] for presenting induction machine dynamics from a two axes general reference frame.

It is well known that electrical machines solve the problem of obtaining a rotating field by employing three windings sinusoidally distributed and separated by 120º (mechanical degrees) which are feeded by three sinusoidal stator 120° electrical degree phase shift. However, because of field distributions are the same along the third dimension (the machine shaft direction), these field distributions are analyzed in the plane where only are needed two linearly independent directions for characterizing any movement. The relationship between the three-phase (A, B and C for stator and a, b and c for rotor) and two phase voltages taking into account natural frames (fixed to the stator sD - sQ and fixed to the rotor rα - rβ for stator and rotor quantities, respectively) are:

[us0usDusQ]=23[1/21/21/211/21/203/23/2][usAusBusC]E1
[ur0urαurβ]=23[1/21/21/211/21/203/23/2][uraurburc]E2

By considering the two axes description, quantitative and qualitative analyses of induction machines can be simplified and also vector control concepts can be used. In this way, the contribution of vector control is based in controlling the induction machines active and reactive powers independently and/or controlling them as DC equivalent ones.

In a general reference frame, all induction machine variables are referred to a real axis known as direct axis xand to the quadrature axis yboth rotating at the reference frame speed ωg=dθg/dtas shown in Figure 1. In this figure, θgis the angle of the real axis xmeasure from sD.

Then, the current stator phasor defined in the general framework, is:

i¯sg=i¯sejθg=isx+jisy

where upper bar indicates phasor quantities.

Also,

u¯sg=u¯sejθg=usx+jusy

ψ¯sg=ψ¯sejθg=ψsx+jψsy

With u¯sand ψ¯sstator voltage and (linked) flux space phasors in the general reference frame.

At the rotor side, Figure 2 shows three frames, rotor (rαand rβ), stator (sDand sQ) and general (xand y) and their angles θr, 0 and θg, respectively.

Meanwhile, the current phasor in the rotor reference frame can be expressed as i¯r=|i¯r|ejαr, in the general reference frame is i¯rg=|i¯r|ejαrwith αr=αr(θgθr). Then,

i¯rg=|i¯r|ejαrej(θgθr)=i¯rej(θgθr)=irx+jiry.

Also, rotor voltage and (linked) flux space phasors in the same frame are:

u¯rg=u¯rej(θgθr)=urx+jury

ψ¯rg=ψ¯rej(θgθr)=ψrx+jψry.

Finally, induction machine phasor expressions are:

u¯rg=Rri¯rg+dψ¯rgdt+j(ωgωr)ψ¯rg,

u¯sg=Rsi¯sg+dψ¯sgdt+jωgψ¯sg

being Rsand Rrstator and rotor resistances, respectively, and u¯rg=00 when a squirrel cage machine is considered. Additionally, stator and rotor fluxes can be expressed in terms of current phasors, and the stator, rotor and magnetizing inductances (Ls, Lr, Lmrespectively)

ψ¯sg=Lsi¯sg+Lmi¯rg

ψ¯rg=Lri¯rg+Lmi¯sg.

Last four equations can be rewritten compactly in matrix form as follows:

[u¯sgu¯rg]=[Rs00Rr][i¯sgi¯rg]+ddt[LsLmLmLs][i¯sgi¯rg']+jωg[LsLmLmLr][i¯sgi¯rg]jωr[00LmLr][i¯sgi¯rg],E3

or by considering real and imaginary components:

[usxusyurxury]=[Rs+pLsωgLspLmωgLmωgLsRs+pLsωgLmpLmpLm(ωgωr)LmRr+pLr(ωgωr)Lr(ωgωr)LmpLm(ωgωr)LrRr+pLr][isxisyirxiry],E4

Where usx, usyand isx, isyare stator voltages and currents in the general reference frame. Identical considerations remain for rotor quantities. If ωg=00, it is obtained the 'conmutator model', but if it is employed ωg=ωsynthe expression [4] can be rewritten as:

[usxusyurxury]=[Rs+pLsω1LspLmω1Lmω1LsRs+pLsω1LmpLmpLmsω1LmRr+pLrsω1Lrsω1LmpLmsω1LrRr+pLr][isxisyirxiry],E5

with ωsynthe synchronous speed and ωsynωr=sωsynwith sthe slip.

## 3. Thevenin equivalent of asynchronous machines

### 3.1. General reference frame - Cartesian coordinates

As presented in Appendix A and beggining with [5,13,14]:

u¯sg=Rsi¯sg+jωgψ¯sgE6
u¯rg=Rri¯rg+dψ¯rgdt+j(ωgωr)ψ¯rg,E7

the Thevenin equivalent of an asynchronous machine in cartesian coordinates is:

u¯sg=Rsji¯sg+ωg(LsLm2Lr)i¯sg+u¯sgE8
du¯sgdt=ωgLmLrju¯rg+RrLm2Lr2ωgji¯sg(ωgωr)ju¯sgRrLru¯sgE9

Meanwhile expression (9) shows the internal voltage dynamics, expression (8) indicates how the stator current changes when u¯sgvaries (assuming that u¯sgis constant because of the induction machine is connected to an electrical grid where the voltage connection remains constant).

For clarity's sake, a typical qualitative analysis involves next steps which are carried out over a fixed speed (squirrel cage) wind generator:

1. At steady state dψ¯rgdt=00.

2. At time t=t1there is a rotor speed change (wind velocity changed).

3. From (7) with urg=00 a rotor flux variation appears (dψ¯rgdt0).

4. Rotor flux change produces an internal stator voltage change provided that u¯sg=jωgLmLrψ¯rg.

5. u¯sgdynamics evolves according to (9).

Active and reactive power are changed by modifying the internal voltage source. Indeed, provided that u¯sg=constant, i¯sgmust change and then the associated powers. As a consequence, in squirrel cage wind generators active and reactive powers change according to the wind velocity. i.e. they are uncontrollable from the wind generator point of view. On the other side, it is known that by modifying pitch blades the active power from squirrel cage wind generators can be regulated. In this case the reactive power is a consequence of the active power control.

### 3.2. Dynamic model in polar coordinates

The Thevenin equivalent, i.e. the internal voltage magnitude and its phase, allows analyzing and considering induction machines (fixed and variable speed wind generators) in an electrical perspective looking for integrating wind generators when studying stability issues of power systems.

Beginning with the Cartesian Cordinates:

du¯sg'dt=ddtusx2+usy2,

du¯sg'dt=1/2usx2+usy2(2usxdusxdt+2usydusydt),

From expressions (8) and (9) it is possible to obtain the internal voltage derivative and its phase:

du¯sg'dt=1u¯sg'[usx(ωgLmLruryRrLm2Lr2ωgisy)+usy(ωgLmLrurx+RrLm2Lr2ωgisx)RrLr(usx2+usy2)]E10

Because of tanδ=usyusx, the voltage phase derivative is

dδdt=1u¯sg'[usx(ωgLmLrurx+RrLm2Lr2ωgisx)usy(ωgLmLruryRrLm2Lr2ωgisy)(ωgωr)(usx2+usy2)]E11

### 3.3. Wind generator model considering active and reactive power delivered

By remembering classical expressions of active and reactive powers:

P=32Re(u¯i¯*)=32(uxix+uyiy)E12
Q=32Im(u¯i¯*)=32(uyixuxiy),E13

and considering expressions (10) and (11):

dU¯sgdt=1U¯sg[23RrLm2Lr2ωgQs+usx(ωgLmLrury)+usy(ωgLmLrurx)RrLr(usx2+usy2)]E14
dδdt=1U¯sg[23RrLm2Lr2ωgPs+usx(ωgLmLrurx)usy(ωgLmLrury)(ωgωr)(usx2+usy2)]E15

it is noted, as expected, that the internal voltage derivative is a function of reactive power while its phase depends on the active power. Note also, that active and reactive powers, Psand Qsrespectively, are internal power sources, i.e. they are not the wind generator output powers.

## 4. Grid flux reference frame

As it was previously indicated, vector control allows controlling machine behaviors in an easier way than others techniques. One of the advantages of properly choosing the reference frame position is that is possible to simplify analysis and control

Explaining vector control is not the scope of this chapter. About vector control of different kinds of machines see [12] and [15]

of electrical machines. In this way, virtual flux reference frame is chosen by considering a virtual flux from grid voltage [16,17]:

ψ¯g=usjωg=j|us|ejθgωg,E16

With θgthe voltage phase and ωgthe phasor speed. According to Figure 3 the proposed reference frame defines a virtual flux (ψ¯g)90ºaway from the voltage grid which belongs to the imaginary axis.

Power expressions, by virtue of the chosen reference frame, are:

P=32Re(u¯i¯*)=32(uxix+uyiy)=32(Uisy),

Q=32Im(u¯i¯*)=32(uyixuxiy)=32(Uisx),

where both expressions show the importance of choosing suitably the reference frame. Indeed it is really simple to control active and reactive power from controlling the stator currents in an independent way.

### 4.1. Squirrel cage machine as a load

According to the previously discussed, it is possible to obtain active and reactive power dynamic models of the induction motors by operating with the presented expressions. These models can be employed in power systems stability studies considering that induction motors are about 60%% of all loads [18]. In this way, it is better to begin with the cartesian coordinates already presented and consider a virtual flux reference frame for an induction machine motor. According to Appendix B:

P3U2Lm(ωgωr)B2RsLsωg3ULmA2RsLsdUdt+3U2LmA2RsLsωgdωgdtK01U2+K02U2ωg+Δωg+K1UdUdt+K2U2ωg+Δωgdωgdt

In both power expressions A=(LsLm2Lr)LrLmRr, B=LrLmRrRs, ωgis the line frequency, Δωgis the frequency deviation and ωris the rotor speed that, because of rotor inertia, remains practically constant in the temporal interval of interest.

## 5. Modelling of DFIG’s operated with vector control

Variable speed wind farms powered by double fed induction (wound rotor) generators (DFIGs) are the other power plants considered in this chapter. Figura 4 shows the main components of a DFIG wind turbine: the rotor, the mechanical transmission system, the doubly fed induction generator and the back to back converters with their respective controls. In general, converters C1 and C2 are operated in an independent way. Meanwhile C1 is operated via vector control driving active and reactive stator powers, C2 maintains the DC bus voltage constant. In subsynchronous speeds, the rotor of the DFIG machine consumes active power meanwhile at supersynchronous speeds delivers it. As a consequence when considering active power delivered by a DFIG wind generator it should be taking into account the active power in the rotor channel.

Because of DFIG machines control is made via C1 converter, all of this chapter considers only C1 control and avoids analyzing C2 even when its operation is similar.

Beginning with the cartesian model in the general reference frame:

u¯sg=Rsji¯sg+ωg(LsLm2Lr)i¯sg+u¯sg,E17
du¯sgdt=ωgLmLrju¯rg+RrLm2Lr2ωgji¯sg(ωgωr)ju¯sgRrLru¯sg,E18

it is noted, due to “j” operator, that there are unwanted coupling terms between rotor and stator circuits along yand xaxes. This coupling can be eliminated by utilizing an additional voltage component u¯rgdec=u¯rgdec1+u¯rgdec2in expression (18), where:

u¯rgdec1=RrLmLri¯sgE19
u¯rgdec2=LrLm(ωgωr)ωgu¯sgE20

Both values should be added to the rotor controllers. Note also that, according to (20)

u¯rgdec2=LrLm(ωgωr)ωgu¯sg=LrLmsu¯sgsu¯sg

with Lr=Lm+LlrLm, sthe slip and Llrthe leakage inductance. Last expression can be approximated, considering that z=Rs+jωg(LsLm2Lr)in (17) is small. Then, u¯sgUand u¯rglimits the variable speed range operation in modern DFIGs by approxiimatedly 30%slip when considering rotor converter size (nominal voltage=30%%Ugrid).

On the other side, even when the voltage feedforward can avoid the undesired coupling, this is not an optimum solution when expression (19) is seen under next considerations

• it is important to maximize the DFIG variable speed range operation, then some coupling can be tolerated,

• as presented later, looking for mimicking conventional synchronous plants operation and control, it can be useful to use the current input, eliminated by (19), for control purposes.

Figure 5 shows, on the left, a block diagram where undesired voltages are eliminated via the feedforward of the stator voltage and current; on the rigth, a simplified equivalent loop where only appears the rotor dynamics. Power references are transformed to rotor voltages ones by vector control [12] indicated as Kin Figure 5.

According to expression (17) and considering Rs0, it is possible to obtain stator currents and, then, delivered powers:

i¯sg=1(LsLm2Lr)ωg(u¯sgu¯sg)

P=32usyisy=32U(u¯syU)(LsLm2Lr)ωgE21
Q=32usyisx=32Uu¯sx(LsLm2Lr)ωgE22

This completes the machine model which is presented in Figure 6 where Δrepresents uncertainties and approximation errors (for example by neglecting Rsor those due to unmatched parameters between the model and the actual machine). After active and reactive powers are measured they are feedback to the controller which is usually a PI one. Note that in vector control (which is presented here by K’) there are more than one feedback loop. Indeed, aside from decoupling voltages (Figure 5) there are two rotor current loops with the externals ones for active and reactive powers control in Figure 6.

## 6. Static models of asynchronous machines

It is well known that some model simplifications, keeping a good compromise between behavior and results exactitude, allow a better qualitative understanding of different kind of processes. An example of this appears studying mechanical behaviors of a wind turbine in presence of a wind velocity change. In this case, and beginning from a steady state condition, electrical behavior can be considered significantly faster than the mechanical response. In this way, the dominant mode is the mechanical one and the electrical modes can be omitted (the stator transitory can be neglected as done in this chapter).

This section obtains the static model of induction machines derived from the dynamic one already presented. Obviously, there will always be time frames in which none of the models presented in this chapter is the best one to analyze a particular problem. In this way, if it is of importance to study the stator transitory of electrical machines, obviously, this dynamic can not be neglected as in this chapter was done. In any case, the importance of determining which model is more adequate to analyze a particular situation always lies in the physical knowledge underlying the problem.

### 6.1. Polar coordinates

Beginning with the dynamic model:

u¯sg=Rsji¯sg+ωg(LsLm2Lr)i¯sg+u¯sg

du¯sgdt=ωgLmLrju¯rg+RrLm2Lr2ωgji¯sg(ωgωr)ju¯sgRrLru¯sg=0,

then:

U¯sg=usx2+usy2 andδ=arctan(usyusx)

### 6.2. Static model as a function of the active and reactive stator powers delivered in the virtual flux reference frame

The voltage internal phase:

tanδ=usyusx

tanδ=RsQ+ωg(LsLm2Lr)P+32U2RsP+ωg(LsLm2Lr)Q

is highly dependent on grid voltage, on active power and, in a lesser way, on reactive power. Note that, if a DFIG is operated at unitary power factory, i.e. Q=0, then δ90°.

On the other side, the internal voltage resuls:

U¯sg=usx2+usy2=usycotδ2+1=usysinδ

U¯sg=Rs23Q+ωg(LsLm2Lr)23P+U2Usinδ

U¯sg=usx2+usy2=usxtanδ2+1=usxcosδ

U¯sg=Rs23P+ωg(LsLm2Lr)23QUcosδ

## 7. Wind farm control

Wind farms have become a visible component of interconnected power grids. In the beginnings of wind generation, when a low portion of the electrical power was delivered from this renewable energy, only simple engineering judgment were necessary to conclude about the negligible impact of wind generation on power systems. Nowadays, with wind farms, but also with high power wind generators in dispersed grids, approaching the output rating of conventional power plants, it is necessary understanding the way in which wind generation can impact and/or contribute to the power system stability.

This section presents different wind farms controls by considering them as current sources and voltages sources (Thevenin equivalent) in an effort for mimicking conventional generation.

The proposed wind farm control approaches are based in the named Lyapunov Theory which give place to linear and non linear wind farm controls. In these approaches, the DFIG capability for controlling active and reactive powers plays an important role in contributing to power system stability. Then, from here on, only DFIG wind farms are considered.

Additionally, Energy (Lyapunov) approach is not based on system linearization and the proposed analysis allows considering any wind farm in any power system in the same way avoiding transform every issue about integrating wind farms in power systems in a different problem.

As indicated, in order to contribute to the network stability, both active and reactive wind farm power controls of a wind farm are considered. Then, steady state controls (normal operating conditions) plus incremental corrections are proposed:

Qwf=QSC+ΔQ,

Pwf=PSC+ΔP.

with Pwf(Qwf)is the total active (reactive) power, PSC(QSC)the power reference given by a Supervisory Control [19,20] and ΔP(ΔQ)the wind farm correction which contributes to the power system stability.

Note that for both active and reactive corrective actions some power reserve is required. Indeed, about active power it will be expected a power reserve which is a function of the wind turbine operating point and about the reactive power correction, the total 'apparent' power of the DFIG machine will limit the corrective action.

### 7.1. Wind farm aggregated model

A complete model of a wind farm with a high number of wind generators, may lead to compute an excessive and impractical number of equations. The size of the wind farm model may be reduced by aggregating several wind turbines with similar incoming wind into a bigger turbine called aggregated turbine [14]. The mechanical and electric parameters per unit are preserved, and the nominal power is increased up to the sum of the nominal power of the whole set of turbines to obtain the parameters of the aggregated turbine. This procedure is employed in this chapter where the wind farm is modeled as one aggregated wind turbine. In this regard, the Thevenin equivalent can be obtained as in a classical problem of electrical systems, applying precisely Thevenin theorem, being the internal voltage and its phase calculated according to the impedance and voltage seen from the wind farm common connection point.

### 7.2. Wind farm control. Method of Lyapunov

Power system stability has been defined as that property of a power system that enables it to remain in a state of operating equilibrium under normal operating conditions and to regain an acceptable state of equilibrium after being subjected to a disturbance [5].

Lyapunov demonstrated that a nonlinear dynamic system:

x˙=f(x),f(0)=0,E23

around the equilibrium point x=00 is asymptotically stable if there exist a scalar function ν(x)>00 for (23) with derivative ν˙(x)<00. The last condition is relaxed to ν˙(x)00 provided that ν˙(x)=00 only vanish at x=00.

Lyapunov theory deals with dynamical systems without inputs. However, it is possible to employ Lyapunov theory in feedback design by making negative the Lyapunov derivative [21,22,23]. The incremental energy function νof a power system without wind farms, where conventional models of the synchronous generators and of load impedance are considered, is [23]:

ν=k=1NG(0.5Mkω˜k2PMkδ˜k)+k=1NL(PLkθ˜k+QLkVkdVk),E24

with

ω˜k=ωkωCOI,ωCOI=˙1MTΣk=1mMkωk,

δ˜k=δkδCOI,δCOI=˙1MTΣk=1mMkδk,

θ˜k=θkδCOI,MT=˙Σk=1mMk,

with NLand NGthe number of loads and generators, Mkthe machine inertia constant, ωkthe machine speed, δkthe angle of the voltage behind the transient reactance, θkthe angle at each bus, PMkthe mechanical power of the generators, PLkand QLkthe active and reactive load powers and Vkthe voltage at the connection point. The angles and speeds are measured with respect to the center of inertia (COI) reference frame (δCOIand ωCOI).

#### 7.2.1. Non-linear power control of wind farms as negative loads

To damp the electromechanical oscillations, i.e. frequency oscillations, the incremental energy function of the power system must decrease. The time derivative of this energy function considering wind farms as negative loads, i.e. acting through its active and reactive powers, yields:

ν˙=k=1NG(Mkω˜˙k2+PGkPMk)δ˜˙k+n=1NLPnθ˜˙n+n=1NLV˙nVnQnV˙wfVwfΔQΔPθ˜˙wf.E25

Then, the active and reactive powers (ΔPand ΔQ) must be chosen in order to allow the sufficient condition of the derivative of the incremental energy function. In equation (25), the expression in between parenthesis is zero (or lesser than zero) because of the generators power balance (equals the internal generator damping). The next two terms correspond to the power balance equations at the nodes and are zero. Then, looking for damping the electromechanical oscillations, the last two expressions must be less than zero.

About the active power two possibilities are choosen [10]:

ΔP=Kc1θ˜˙wf;withKc1>0,E26
ΔP=Kc2θ˜¨wf2θ˜˙wf;withKc2>0.E27

meanwhile the first expression which is the classical proportional frequency law, the second one is a kind of inertial response. Indeed, the second expression can be understood as a modification of (26) where Kc1is a variable gain Kc1=Kc2θ˜¨wf2which takes into account frequency derivative.

On the other hand, a non-linear control strategy of the wind farm through the reactive power ΔQcan be derived from expression (25) by considering that the wind farm emulates the behavior of a static var compensator:

ΔQ=buVwf2

V˙wfVwfΔQ=buVwfV˙wf=12buddtVwf2.

In this way, in order to keep the sufficient condition, some possibilities arise:

bu=KrdVwf2dt;withKr>0orE28
bu=Krsign(dVwf2dt);withKr>0,E29

Where buis the (equivalent) wind farm susceptance.

With respect to the active power control, note that the idea of a power reserve, as a percentage of the maximum available power, is very attractive from a point of view of the network stability and it is usually employed [19,20]. However, suppose that ΔP=00, the appropriate choice of the reactive power function (expressions (28) or (29)) implies that the energy function derivative (34) almost always decreases, i.e. any electromechanical oscillation is damped. Then, the wind farm reactive power contributes to damp the electromechanical oscillations of the power system. Being, in general, ΔP0, the reactive power function reinforces the wind farm contribution to the network stability.

Note that, with both active and reactive control laws, it is possible maximizing the use of the energy resource in order to contribute to damp the electromechanical oscillations by exploiting all the capabilities of the DFIG machines. In this way, it is possible designing a control law for the reactive power which takes into account the advantage of producing as much reactive power as possible considering the apparent power of the DFIGs.

#### 7.2.2. Non-linear active power control of wind farms as Thevenin equivalent

Due to it is expected that wind farms act as power plants [24], it is necessary to demostrate that wind farms behave as their equivalent synchronous generators (the conventional power plants) with proportional and derivative (inertial) frequency control laws. The equations representing the dynamic behavior of synchronous generators for the reduced model, are [5]:

δ˜˙k=ω˜k,E30
ω˜˙k=1Mk(PmkDkωkPGk).E31

where Mkis the inertia of the whole machine (synchronous generator plus prime mover), Pmkis the mechanical power produced by the prime mover, Dkis the component of internal friction of the generator and PGkis the electrical power injected in the network. Because of variables are in per unit, powers and torques are equal [5]. Mimicking the analysis for a wind farm with frequency control and inertial contribution yields [19]:

α˜˙k=Ωgk,E32
Ω˜˙gk=1Mgk(PtPwf),E33

withPtthe turbine power.

Considering the control as:

Pwf=PSC+ΔP,

ΔP=Kp(ω˜refω˜)+Kd(ω˜˙refω˜˙)=Kpω+Kdω˙,

where is included a PD control named “proportional and inertial” classical control laws, then:

ω˙=1Kd(Pt(Kpω˜+MΩ˜˙gk)PSC).E34

Note the similarity of this expression with (31) for synchronous generation. Then, the analysis focuses on addressing the control of the wind farm as in ν˙1which considers the derivative of the Energy Function of a synchronous generator in (34):

ν1=k=1NG(0.5Mkω˜k2+PGkθ˜kPMkθ˜k)ν˙1=k=1NG(Dω˜k)<0E35

The equivalent expression for a wind farm is:

ν1wf=0.5ωk2+PSCKdθwfPtKdθwfν˙1wf=ωkω˙k+PSCPtKdθ˙wf=(ω˙k+PSCPtKd)θ˙wf,

ν˙1wf=(MwfKdΩ˙gKpKdω)θ˙wf=MwfKdΩ˙gθ˙wfKpKdθ˙wf2,E36

being Kp>00 and Kd>00 it is verified the negative sign of the second term in the last expression. On the other side, in order to verify the negative sign of MwfKdΩ˙gθ˙wfit is necessary to analyse the wind farm convergence to an equilibrium point (e.p.), knowing that Mwf>0and considering that the wind farm (the aggregated turbine) is outside the equilibrium point.

Figure 7 presents the torque - speed curves of the aggregated wind turbine with the wind velocity as a parameter. The e.p., considering constant wind velocity, corresponds to nominal frequency at the wind farm connection point with constant shaft speed Ωe.p.of the aggregated turbine.

In order to verify the convergence to the equilibrium, two conditions outside the e.p. [19], which are consequence of electrical disturbances, will be analyzed. First, consider that because of a disturbance action, the aggregated wind turbine is operating at point A (Figure 7) with ΩA<Ωe.p., being the frequency θA<θe.p.. At that point, the wind farm generated power is higher than the nominal one, i.e. the wind farm is contributing to restore the frequency at the connection point.

When the disturbance disappears, the network returns to its normal configuration and the wind farm power is higher than that which maintains the power balance in the system. As the deviation of the frequency decreases, so does the wind farm generated power. As a consequence, the wind turbine torque decreases and the turbine speed experiences an increment until the speed reaches Ωp.e.. Thus, while the frequency reaches their nominal value, the wind turbine increasing its speed. Then, the (negative) sign of (36) is verified by considering:

θ˙wf=θ˙refθ˙A>0andE37
Ω˙g=Ωe.p.ΩAΔt>0.E38

If as a consequence of another disturbance the aggregated turbine is operating at point B (Figure 7), when the disturbance is removed the sign of (36) yields:

θ˙wf=θ˙refθ˙A<0andE39
Ω˙g=Ωe.p.ΩBΔt<0.E40

Expressions (38) and (40) verify the negative sign of (36).

#### 7.2.3. Non-linear reactive power control of wind farms from the Thevenin equivalent

In last subsection, by mimicking the reduced model of conventional synchronous generators with wind farms, it has been deduced that wind farms PD laws, i.e. the named frequency and inertial responses control approaches, allows to contribute to the power system stability. In order to deduce a wind farm reactive power control, it is necessary to include the named "Structure Preserving the Model" [21] or "the one axis model" [6] for the conventional synchronous generators. The dynamics of a kthsynchronous generator, respect to the COI is [21]:

δ˜˙k=ω˜k,E41
Mkω˜˙k=PmkDkωkPGkMkMTPCOIE42
TdokdEqkdt=xdkxdkxdkVn+kcos(δkθn+k)+EfdkxdkxdkEqkE43

being Tdokthe daxis transient open circuit time constant; Eqkthe qaxis voltage behind transient reactance, Efdkis the exciter voltage which is assumed constant (if the exciter control action is included in the generator model, at least one additional dynamic expression should be included) and xdkand xdkare daxis synchronous reactance and transient reactance, respectively.

An advantage of this model when compared to the classical one with two states (expressions (30) and (31)) is the possibility of including loads where the impedance are not constant. In order of making more clear the explanation it will not be included any wind farm in this step. According to [21] next terms are added to the already applied Lyapunov function:

ν2a=i=n+12n12xdin[Eqin2+Vi22EqinVicos(δinθi)],

ν2b=i=1nEfdiEqixdixdi,

ν2c=i=1nEqi22(xdixdi),

ν2d=12i=n+1n+Nl=n+1n+NBklVkVlcos(θiθl),

ν2e=i=n+12nxdinxdin4xdinxdin[Vi2Vi2cos(2(δinθi))].

Then, the derivative of the (new) Lyapunov function is:

ν˙=k=1NG(Mkω˜˙k2+PGkPMk)δ˜˙k+n=1NLPnθ˜˙n+n=1NLV˙nVnQn+ν˙2a+ν˙2b+ν˙2c+ν˙2d+ν˙2e,

ν˙=k=1NGDk(ω˜˙k2)k=1NGTdxdxdE˙'q2,

where:

(Pnθ˜˙n+ν2aθ˜nθ˜˙n+ν2dθ˜nθ˜˙n+ν2dθ˜nθ˜˙n)=(Pn+Pk)θ˜˙n=0(V˙nVnQn+ν2aVnV˙n+ν2dVnV˙n+ν2dVnV˙n)=(Qn+Qk)V˙nVn=0

i.e. the total active and reactive powers (consumed plus injected) is zero.

When it is considered the last state variable, the voltage behind transient reactance Eqk, in the derivative of the Lyapunov function:

ν˙3=ν2aEqdEqdt+ν2bEqdEqdt+ν2cEqdEqdt=k=1NGTdxdxdE˙'q2

Due to Td>0and xd>xdthe Lyapunov condition is met.

Mimicking the third state equation of the Structure Preserving Model with the aggregated wind turbine (with DFIGs) implies taking dynamical expressions from (18) and eliminating an undesired coupling:

du¯sgdt=ωgLmLrju¯rg+RrLm2Lr2ωgji¯sg(ωgωr)ju¯sgRrLru¯sgE44
Trdu¯sgdt=Lm2Lrωgji¯sg+ωgLmRrju¯rgu¯sgE45

According to virtual flux control, isycontrols active power meanwhile isxcan be used for reactive power control (Q=32usyisx). Operating with expression (45):

Trdusydt=Lm2Lrωgisx+ωgLmRrurxusy

Trdusydt=(LsLeq)ωgisx+ωgLmRrurxusy

and considering, expression (8), that Leq=(LsLm2Lr).

According to Figure 8, which presents the Thevenin equivalent and the correspoding phasor diagram, and considering expression (17) with Rs0, results:

i¯sx=1Leqωg(u¯sxu¯sx)=u¯sxLeqωg=Vsin(δθ)Leqωgi¯sx=Vcos[90(δθ)]Leqωg

Then,

Trdusydt=(LsLeq)LeqVcos[90(δθ)]+ωgLmRrurxusy,

1ωgRrLmTrdusydt=1ωgRrLm(LsLeq)LeqVcos[90(δθ)]+urx1ωgRrLmusy,

with δthe internal voltage phase, θthe voltage angle at the common connection point and Tr=Lr/Rr. Last expression is pretty similar to (43). Then, an equivalent Lyapunov analysis can be done which implies that it is possible to contribute to the power system stability as conventional generation does.

## 8. Tendencies in the study of wind farm contributions to system stability

It is worth to note that the objective of this chapter is not give a full list of new trends in wind farms control. The comments are focused on aspects in which, in the authors opinion, are expected some improves in the near future. This does not mean to avoid exploring other control techniques that have been successful in different fields of control systems.

In this way exploring another control Lyapunov functions will certainly be an important input not just for wind farm contributions in power systems but also for power systems stability in general. Indeed, acting as negative loads implied than wind farm control can be extrapolated to other kind of devices which can control active and/or reactive power independently, as FACTS. Other important subject is concerning with the load characterization. Then, including in the Lyapunov frame some loads models as the developed for the induction motor will allow to find new rules with the advantages of decentralized and local measures as the developed ones in the chapter.

However, there would be particular cases in which some rules based in local measures, as the frequency response, could be ineffective. This is presented in Figure 9 taking the topology of a power system from [5] and including a wind farm in different places. In part (b) of Figure 9, after a disturbance appears and because of the wind farm is in the middle of areas, wind farm frequency remains constant due to that measure is taken in the center of 'bouncing' between areas. Indeed, in Figure 9(b) wind farm is located in the COI of the system, then ω˜k=ωkωCOI=ωCOIωCOI=0.0

Note that the power system in Figure 9 can be analyzed as an academic subject of study due to actual symmetrical systems will always have differences which will be enough to allow that classical control laws contribute to the power system stability. In this cases, obviously, the Lyapunov control laws derived in this chapter could be less efective than other ones which take into account the nature of this electrical grid.

About grids with symmetrical topologies, Passivity theory [25,26,27] can allow to consider some other tools looking for contributing to the power system stability [28,29]. Because of Passivity considers finding a controller in such a way the dynamical system energy function takes the desired form ('energy shaping') and lately considering a power shaping approach, but also solving an energy function which qualifies as a Lyapunov function, Passivity theory implies an important step ahead in the power systems study. Some authors attempts about including wind farm control in the passivity frame can be found in [30,31].

## 9. Conclusions

This chapter can be divided in two parts, meanwhile the first one is devoted of exploring models and developing a Thevenin equivalent of squirrel cage induction generators (fixed wind generators) and DFIG ones (variable speed generators), after introducing an aggregated wind farm concept and considering a wind farm as another Thevenin equivalent, the second part of the chapter considers analyzing DFIG wind farm linear and non linear control.

The final objective of the chapter is to demonstrate that can be developed an equivalent behavior of wind farms to their counterparts, the synchronous conventional generators, by properly controlling wind farms. This point of view will help to promote wind farm integration but, at the same time, it will open new doors on contributing to the power systems stability via wind farms control.

About the first part, induction machines models were presented including an additional dynamic model of a squirrel cage induction machine as a load. Additionally, even when vector control was not explained, a control perspective was adopted in showing the way in that stator currents (and powers) can be controlled from rotor voltages in DFIG machines. Also, by considering active and reactive powers delivered for wind generators as input, a static model of induction machines was derived from its dynamical counterpart.

About controlling wind farms, it was possible maximizing the use of the energy resource in order to contribute to damp the electromechanical oscillations by exploiting all the capabilities of the DFIG machines. In this way, wind farms were operated under a Supervisory Control which imposed steady states power references and (added) corrective actions, under Lyapunov Theory which was shortly explained, were proposed as a complement of the references ones.

By considering wind farms as negative loads the Lyapunov frame showed that meanwhile wind farm proportional active power control law was the same as the classical proposed one, about the inertial effect a modification from that one was found. About reactive power, wind farms emulated the behavior of static VAR compensators giving place to a highly non linear control law.

By considering wind farms as voltage sources the Lyapunov frame demonstrated that classical active power laws, proportional and derivative with frequency, can contribute to the power system stability. Also, about reactive power, it was possible demonstrating that wind farms can contribute to damping electromechanical oscillations by controlling its internal (Thevenin) voltage via rotor voltage actions in DFIGs.

As a consequence of the Energy (Lyapunov) approach, the obtained control laws were not based on the linearization of the system. This assures a bigger domain of attraction of the wind farm contribution indicating that, even under severe disturbances, the proposed control laws will contribute to the power system stability.

It is important to note the local nature of the signals used for control laws which avoids any coordination with the rest of the system. Additionally, the damping is not dependent on the power flow direction and neither the kind of failures on the power system.

However, aside from Energy Functions included only some kind of loads and that line resistances must be neglected in the calculus of the functions, another limitation can be found from considering some power systems and wind farms locations where active (or reactive) power control laws can not be effective on contributing to the system stability.

Looking for solving some of the aforementioned drawbacks, future research in Lyapunov topics were indicated but also other perspective, which in fact calculates Energy Functions in control systems, as Passivity was proposed. In this way, the same border conditions remain for future control laws as the necessity of local uncoordinated measures, independency of power system topology, robustness, etc.

Finally, it is important to note some differences faced when looking for control laws in power systems. Electrical grids have distinctive characteristics when compared with other “conventional” control systems. In first place, as indicated, the necessity and importance of uncoordinated and robust control laws which contribute to the whole system stability but based from local measures. In second place, wind farm control laws must cover an important range of applications by considering different power systems and different places in which a wind farm can be connected. Indeed, control systems teach us that it is of fundamental importance knowing the "system model", however that model is highly changing and it has infinite equilibrium points when talking about power systems. All of these characteristics only emphasize that power systems are difficult to control and as technology evolves they are upgraded in such a way that some challenges are solved and, at the same time, new ones appear.

## Appendix A

Beginning with the expressions in the general reference frame and neglecting stator dynamics (dψ¯sg/dt=0)[5,13,14]:

u¯sg=Rsi¯sg+jωgψ¯sgE46
u¯rg=Rri¯rg+dψ¯rgdt+j(ωgωr)ψ¯rg,E47

with

ψ¯sg=Lsi¯sg+Lmi¯rg;ψ¯rg=Lri¯rg+Lmi¯sg

i¯rg=ψ¯rgLmi¯sgLr andψ¯sg=(LsLm2Lr)i¯sg+LmLrψ¯rg

which replaced in expression (46) gives

u¯sg=[Rs+jωg(LsLm2Lr)]i¯sg+jωgLmLrψ¯rgu¯sg=[Rs+jωg(LsLm2Lr)]i¯sg+u¯sg

u¯sg=[Rs+jωg(LsLm2Lr)]ji¯sg+u¯sgE48

where u¯sgis the internal voltage source and the associated impedance is [Rs+jωg(LsLm2Lr)]. Also, from expression (47):

u¯rg=Rri¯rg+dψ¯rgdt+j(ωgωr)ψ¯rg=Rrψ¯rgLmi¯sgLr+dψ¯rgdt+j(ωgωr)ψ¯rg,

dψ¯rgdtjωgLmLr=[u¯rg+RrLmLri¯sgRrψ¯rgLrj(ωgωr)ψ¯rg]jωgLmLr

du¯sgdt=ωgLmLrju¯rg+RrLm2Lr2ωgji¯sg(ωgωr)ju¯sgRrLru¯sgE49

Summarizing:

u¯sg=Rsji¯sg+ωg(LsLm2Lr)i¯sg+u¯sgE50
du¯sgdt=ωgLmLrju¯rg+RrLm2Lr2ωgji¯sg(ωgωr)ju¯sgRrLru¯sgE51

## Appendix B

P=32(usyisy),Q=32(usyisx).E52

Taking i¯sgand i¯rgfrom (46) and (47), respectively and consideringusy=U:

i¯sg=u¯sgjωgψ¯sgRs and ψ¯sg=Lsi¯sg+Lmi¯rgi¯sg=u¯sgjωg(Lsi¯sg+Lmi¯rg)Rs,

and ψ¯rg=Lri¯rg+Lmi¯sgi¯rg=dψ¯rgdtj(ωgωr)ψ¯rgRr.

Taking the last three expressions it is possible to get i¯rgas a function ofi¯sg, u¯sg, dωgdt, ωg, ωr. Then, replacing i¯rgin the second expression and from the first one it is possible to obtain i¯sgas a function of dωgdt, du¯sgdt, u¯sg. Finally, active and reactive powers become:

Due to Rs0when P>10HP[12] and considering low slip, some simplifications can be done:

PK01U2+K02U2ωg+Δωg+K1UdUdt+K2U2ωg+Δωgdωgdt

Q32U(ARsdUdt(ARsdωgdt+(ωgωr)BRs)Uωg)(ωgωr)A+

## Acknowledgments

This work was supported by National University of Patagonia San Juan Bosco.

## Notes

• Explaining vector control is not the scope of this chapter. About vector control of different kinds of machines see [12] and [15]

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2013 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Roberto Daniel Fernández, Pedro Eugenio Battaiotto and Ricardo Julián Mantz (March 20th 2013). Wind Farms as Negative Loads and as Conventional Synchronous Generation – Modelling and Control, Modeling and Control Aspects of Wind Power Systems, S. M. Muyeen, Ahmed Al-Durra and Hany M. Hasanien, IntechOpen, DOI: 10.5772/54089. Available from:

### More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

### Related Content

#### Modeling and Control Aspects of Wind Power Systems

Edited by S. M. Muyeen

Next chapter

#### An Integrated Power Flow Solution of Flexible AC Transmission Systems Containing Wind Energy Conversion Systems

By E. Barrios-Martinez, L.M. Castro, C.R. Fuerte-Esquivel and C. Angeles-Camacho

#### Wind Power

Edited by S. M. Muyeen

First chapter

#### Introduction

By S.M.Muyeen

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.