Solved Problems · Set 34

Load Frequency Control and AGC

Part 7 · Operation, Control and Economics — droop and load sharing, area stiffness, integral gain, tie-line bias and the area control error, from the first second after a load step to a settled interchange. Chapter 33 of the textbook.

Prof. Mithun Mondal 20 solved · 12 practice · 3 challenge GATE · ESE · University

Set 34 — Load Frequency Control and AGC

Twenty worked problems on the control loop that keeps generation equal to load second by second. Sets 27 and 28 decided which units run and at what output, on the assumption that demand was known; it never is, and the gap is closed by governors acting in seconds and by automatic generation control acting in minutes. Almost every question here is answered by one of three lines — \(\Delta f = -\Delta P_L/\beta\) for a lone area, \(\Delta f = -\Delta P_L/(\beta_1+\beta_2)\) with \(\Delta P_{12} = \beta_2\,\Delta f\) for two areas on droop, and \(\mathrm{ACE}_i = \Delta P_{\text{tie},i} + B_i\,\Delta f\) once the supplementary loop is in service. The work is in deciding which of the three the question describes and in getting every droop into megawatts per hertz before anything else happens.

Textbook Chapter 33 · 20 solved · 12 practice · 3 challenge · 12 MCQs

i Method Recap
  • Imbalance appears as a rate of change of frequency. The network stores no energy, so any mismatch is taken from the rotors: \(\dfrac{2H}{f_0}\dfrac{d(\Delta f)}{dt} = \Delta P_m - \Delta P_L - D\,\Delta f\). Setting \(t=0^+\), before any valve moves, gives \(\mathrm{RoCoF} = -f_0\,\Delta P_L/(2H)\) — the one quantity inertia alone decides.

  • Droop must be converted before use. A droop of \(R\%\) means \(R\%\) of nominal frequency drives the unit from no load to full load, so \(1/R = P_{\text{rated}}/\big[(R\%/100)f_0\big]\) MW/Hz. Nearly every wrong answer in this topic is a droop left in percent.

  • Droop is proportional action, so a residual error is structural. Each governor contributes \(-\Delta f/R_i\). If \(\Delta f\) returned to zero the extra generation would vanish with it, so the primary loop cannot remove its own error — and that is exactly what makes many governors able to share one frequency.

  • Regulations add because every unit sees the same frequency. \(1/R_{\text{eq}} = \sum_i 1/R_i\) over the responsive units only, and \(\beta = D + 1/R_{\text{eq}}\) in MW/Hz. A unit at its ceiling, or off governor control, contributes nothing.

  • The central result of primary control is \(\Delta f_{ss} = -\Delta P_L/\beta\), obtained either from the final-value theorem or in two lines from the steady-state power balance \(-\Delta f/R - D\,\Delta f = \Delta P_L\). Every lag contributes unity at \(s = 0\), so the dynamics drop out entirely.

  • Only an integrator removes the offset. A constant output from an integrator demands a zero input, so \(\Delta P_{\text{ref}} = -K_I\!\int\!\Delta f\,dt\) forces \(\Delta f = 0\) in any steady state that exists. A proportional supplementary term merely adds to \(\beta\) and never reaches zero.

  • A tie line is an integrator driven by a frequency difference: \(\Delta P_{12}(s) = \dfrac{2\pi T_{12}}{s}\big[\Delta F_1(s)-\Delta F_2(s)\big]\) with \(T_{12} = \dfrac{|V_1||V_2|}{X_{12}}\cos(\delta_1^0-\delta_2^0)\). Flow settles only when the two areas are at one frequency, which is what "synchronous interconnection" means.

  • Two areas on droop alone share by stiffness. \(\Delta f = -\Delta P_{L1}/(\beta_1+\beta_2)\) and \(\Delta P_{12} = -\beta_2\Delta P_{L1}/(\beta_1+\beta_2)\): the neighbour supports a disturbance it did not cause, which is wanted for seconds and unwanted after a minute.

  • The area control error resolves that. \(\mathrm{ACE}_i = \Delta P_{\text{tie},i} + B_i\,\Delta f\), with export positive. Setting \(\mathrm{ACE}_1 = \mathrm{ACE}_2 = 0\) and adding gives \((B_1+B_2)\Delta f = 0\), so any positive biases end at \(\Delta f = 0\), \(\Delta P_{12} = 0\). The bias decides the path, not the destination.

  • With \(B_i = \beta_i\) control is non-interactive: the disturbed area's ACE equals its own disturbance exactly and every other area's ACE is zero. The bias setting encodes the neighbour's expected response, so no communication between areas is needed.

  • Two allocations run in sequence. In the first seconds the increment divides by \(1/R_i\) — a technical split that ignores cost; over the following minutes AGC re-divides it by the participation factors \(\alpha_i = \big(d^2C_i/dP_i^2\big)^{-1}\big/\sum_j\big(d^2C_j/dP_j^2\big)^{-1}\) from Set 27. The transient state is physical, the final state economic.

Problem 1DrillDroop in MW/Hz

A 250 MW turbo-alternator has a speed regulation of 5% and is carrying 150 MW at exactly 50.00 Hz. The speed-changer setting is not touched.

  1. Find \(1/R\) in MW/Hz.
  2. The system frequency falls to 49.70 Hz. Find the new output.
  3. Find the frequency at which this characteristic would carry no load, and the frequency at which it would carry full load.
Solution

Convert the droop. A droop of 5% means that 5% of nominal frequency — that is, 2.5 Hz — moves the unit across its whole range of 250 MW:

\[ \frac{1}{R} = \frac{P_{\text{rated}}}{(R\%/100)\,f_0} = \frac{250}{0.05\times 50} = \frac{250}{2.5} = 100\ \text{MW/Hz} \]

Read it as a sensitivity: this machine offers 100 MW for every hertz the system frequency is depressed, until it runs out of capacity.

The new output. With \(\Delta P_{\text{ref}} = 0\) the governor characteristic is simply \(\Delta P = -\Delta f/R\):

\[ \Delta f = 49.70 - 50.00 = -0.30\ \text{Hz} \quad\Rightarrow\quad \Delta P = 100 \times 0.30 = 30\ \text{MW} \]
\[ P = 150 + 30 = 180\ \text{MW} \]

Note that the unit does not ask why the frequency fell, nor where. It responds to the local measurement of a global quantity, which is what allows primary control to work without any communication.

The ends of the characteristic. Full load requires a further 100 MW above the present 150 MW, and no load requires 150 MW less:

\[ \Delta P = +100 \Rightarrow \Delta f = -\frac{100}{100} = -1.00\ \text{Hz} \quad\Rightarrow\quad f_{\text{full load}} = 49.00\ \text{Hz} \]
\[ \Delta P = -150 \Rightarrow \Delta f = +1.50\ \text{Hz} \quad\Rightarrow\quad f_{\text{no load}} = 51.50\ \text{Hz} \]

The two ends are 2.50 Hz apart, which is the 5% of 50 Hz that the definition promised. The 150 MW operating point sits 60% of the way along, so the characteristic is not symmetric about 50 Hz — where it sits is decided entirely by the speed-changer setting.

The picture. Sliding the speed changer moves this line bodily, parallel to itself, without changing its slope:

f0 f0+df dP_ref df original setting after the speed changer moves slope = -R P (MW) f (Hz)
The droop characteristic. Frequency deviation moves the unit along the line; the speed-changer command moves the line itself.
The percentage droop is a statement about the whole range of the machine, not about its present output. A 250 MW unit and a 25 MW unit at the same 5% setting have regulations of 100 and 10 MW/Hz — a factor of ten apart — even though they look identical on a nameplate. That is the entire content of the conversion, and forgetting it is the commonest single error in this chapter.
Answer\(1/R = 100\) MW/Hz; output rises to 180 MW at 49.70 Hz; the characteristic runs from 51.50 Hz at no load to 49.00 Hz at full load
Problem 2Exam levelLoad Sharing

An isolated 50 Hz system is supplied by three units of 400 MW, 300 MW and 300 MW, all set to 4% droop and all on governor control with reserve available. The composite load has a damping constant of 12 MW/Hz. A block of 45 MW of load is switched on.

Find the settled frequency, the output change of each unit, and the megawatts contributed by load relief. Verify that the books balance.

Solution

Regulations, one unit at a time. Each uses its own rating with the common 4% and 50 Hz:

\[ \frac{1}{R_1} = \frac{400}{0.04\times50} = 200, \qquad \frac{1}{R_2} = \frac{1}{R_3} = \frac{300}{0.04\times50} = 150\ \text{MW/Hz} \]

They add, because every unit sees the same frequency. That is the whole justification — the three governors are not in series with anything, they are three parallel responses to one measurement:

\[ \frac{1}{R_{\text{eq}}} = 200 + 150 + 150 = 500\ \text{MW/Hz}, \qquad \beta = D + \frac{1}{R_{\text{eq}}} = 12 + 500 = 512\ \text{MW/Hz} \]

The frequency. One line:

\[ \Delta f = \frac{-\Delta P_L}{\beta} = \frac{-45}{512} = -0.08789\ \text{Hz} \qquad\Rightarrow\qquad f = 49.9121\ \text{Hz} \]

Under 0.09 Hz for a 45 MW step, on a system of 1000 MW of installed capacity. That is what a stiff area looks like.

The allocation. Each unit moves along its own droop line by \(-\Delta f/R_i\):

\[ \begin{array}{lcl} \Delta P_1 &=& 200 \times 0.08789 = 17.58\ \text{MW} \\ \Delta P_2 &=& 150 \times 0.08789 = 13.18\ \text{MW} \\ \Delta P_3 &=& 150 \times 0.08789 = 13.18\ \text{MW} \end{array} \]

The split 17.58 : 13.18 : 13.18 is exactly 400 : 300 : 300. Equal percentage droops produce sharing in proportion to rating, so every machine picks up the same 4.4% of its own capacity — the fair and safe outcome, and the reason grid codes specify a common droop band.

Load relief closes the books. The load itself is smaller at the lower frequency:

\[ -D\,\Delta f = 12 \times 0.08789 = 1.05\ \text{MW} \]
\[ 17.58 + 13.18 + 13.18 + 1.05 = 45.00\ \text{MW}\ \checkmark \]

Load damping supplies 2.3% of the correction here. It is small, but it is what remains when every governor has run out of reserve — Problem 5 shows what that looks like.

Two different sums have to agree, and checking them is the cheapest error trap in this subject. The governor increments must add to \(\Delta P_L + D\,\Delta f\), and the individual increments must be in the ratio of the individual \(1/R_i\). If either fails, a droop was left in percent somewhere.
Answer\(\beta = 512\) MW/Hz, \(f = \mathbf{49.912}\) Hz; 17.58, 13.18 and 13.18 MW from the three units, plus 1.05 MW of load relief
Problem 3DrillUnequal Droops

Two units run in parallel on an isolated 50 Hz system: unit A is 150 MW at 3% droop, unit B is 350 MW at 6% droop. Load damping may be neglected. The load rises by 40 MW.

  1. Find the new frequency and the share taken by each unit.
  2. Compare the sharing with the ratio of the ratings.
  3. Invert the question: what droop must be set on B so that the two units share strictly in proportion to their ratings?
Solution

Regulations. Different ratings and different percentages, so both must be converted separately:

\[ \frac{1}{R_A} = \frac{150}{0.03\times50} = \frac{150}{1.5} = 100\ \text{MW/Hz}, \qquad \frac{1}{R_B} = \frac{350}{0.06\times50} = \frac{350}{3.0} = 116.67\ \text{MW/Hz} \]

The small machine is nearly as stiff as the large one. Its 3% setting compensates almost entirely for its being less than half the size.

Frequency and shares, with \(D = 0\) so that \(\beta = 216.67\) MW/Hz:

\[ \Delta f = \frac{-40}{216.67} = -0.18462\ \text{Hz} \qquad\Rightarrow\qquad f = 49.8154\ \text{Hz} \]
\[ \Delta P_A = 100 \times 0.18462 = 18.46\ \text{MW}, \qquad \Delta P_B = 116.67 \times 0.18462 = 21.54\ \text{MW} \]

The comparison, which is the point of the problem:

\[ \begin{array}{lcc} & \text{share of the 40 MW} & \text{share of capacity} \\ \hline \text{Unit A} & 46.2\% & 30.0\% \\ \text{Unit B} & 53.8\% & 70.0\% \end{array} \]

Unit A picks up 12.3% of its own rating while unit B picks up 6.2% of its. The stiffer machine is worked more than twice as hard relative to its size, and would reach its ceiling first on a large disturbance.

The inversion. Proportional sharing requires the regulations themselves to be in the ratio of the ratings:

\[ \frac{1/R_B}{1/R_A} = \frac{350}{150} \quad\Rightarrow\quad \frac{1}{R_B} = 100 \times \frac{350}{150} = 233.33\ \text{MW/Hz} \]
\[ R_B\% = \frac{P_{\text{rated}}}{(1/R_B)\,f_0}\times 100 = \frac{350}{233.33\times 50}\times100 = 3.00\% \]

Equal to A's setting, as it had to be. The general statement is worth remembering: sharing is proportional to rating if and only if the percentage droops are equal, whatever the ratings are. Unequal percentages are used deliberately — a base-load nuclear unit given 8% droop moves less, an aeroderivative gas turbine given 3% moves more.

The two characteristics drawn together make the mechanism visible: one common frequency drop, two different horizontal excursions set by two different slopes.

50.000 49.815 UNIT A 150 MW, 3% UNIT B 350 MW, 6% 18.46 MW 21.54 MW unit output (MW) f (Hz)
One frequency, two slopes. The common vertical drop of 0.185 Hz produces different horizontal excursions, and the split is fixed by the slopes alone.
Percentage droop, not megawatts per hertz, is the setting an operator actually enters — and the two behave oppositely. Equal percentages give proportional sharing; equal MW/Hz would load a small machine to destruction. When a question quotes droops in percent it is telling you about fairness; when it quotes them in MW/Hz it is telling you about stiffness.
Answer\(f = 49.815\) Hz; A takes 18.46 MW (46.2%) against a 30% share of capacity; setting B to 3% droop restores proportional sharing
Problem 4Exam levelMeasured Stiffness

A control area of 6000 MW installed capacity is believed to have every unit on governor control at 4% droop, and its load damping is known to be 25 MW/Hz. A 250 MW generator trips; the recorded frequency settles at 49.83 Hz from an initial 50.00 Hz.

Find the area frequency response characteristic that the event implies, compare it with the value the data predicts, and state what fraction of the capacity was actually responding.

Solution

Invert the central result. A generator trip is a load increase as far as the surviving system is concerned, so \(\Delta P_L = +250\) MW:

\[ \beta_{\text{measured}} = \frac{-\Delta P_L}{\Delta f} = \frac{250}{0.17} = 1470.6\ \text{MW/Hz} \]

This is the only way \(\beta\) is ever really known. It is not computed from nameplates; it is measured from events, and every large interconnection maintains a running record of them.

What the data predicts, if the belief about the plant were true:

\[ \frac{1}{R} = \frac{6000}{0.04\times50} = 3000\ \text{MW/Hz}, \qquad \beta_{\text{expected}} = 25 + 3000 = 3025\ \text{MW/Hz} \]
\[ \Delta f_{\text{expected}} = \frac{-250}{3025} = -0.0826\ \text{Hz} \quad\text{against the measured } -0.17\ \text{Hz} \]

The system is half as stiff as the paperwork says. A factor of two is far too large for the damping constant to explain, so the discrepancy must be in the governor term.

Extract the governor response actually delivered. Damping is not in doubt:

\[ \frac{1}{R_{\text{actual}}} = \beta_{\text{measured}} - D = 1470.6 - 25 = 1445.6\ \text{MW/Hz} \]
\[ \frac{1445.6}{3000} = 0.482 \qquad\Rightarrow\qquad \text{48.2\% of the fleet responded} \]

Say it in megawatts, which is what an operator can act on. The capacity that behaved as 4%-droop plant is

\[ P_{\text{responsive}} = \frac{1}{R_{\text{actual}}}\times 0.04 \times 50 = 1445.6 \times 2 = 2891\ \text{MW} \]

So roughly 3100 MW of spinning plant contributed nothing. The usual reasons are all mundane and all common: units already at maximum output, units placed on a fixed megawatt setpoint for commercial reasons, units whose governor valves were fully open at the moment of the trip, and units whose governor deadband exceeded 0.17 Hz.

Why the distinction matters. Reserve on paper is not \(\beta\). Capacity that is spinning but blocked from responding appears in the reserve report and not in the frequency response, so an area can be fully reserved and still fall twice as far as expected. The quantity that must be tracked is responsive reserve, and this event has just measured it.

Every large system operator now runs exactly this calculation after every unit trip. The measured \(\beta\) from an event is compared with the declared governor settings, and the gap is the fleet's non-compliance. It has been closing steadily since 2016 because grid codes began to require it to be reported — and, more to the point, because the frequency bias settings of Problem 11 are supposed to equal \(\beta\), and a bias set from nameplate data would here be twice the truth.
Answer\(\beta_{\text{measured}} = \mathbf{1471}\) MW/Hz against 3025 predicted; only 2891 MW — 48.2% of the fleet — actually responded
Problem 5Exam levelLoad Damping

An area supplies 4000 MW at 50 Hz. Measurement shows that a 1% change in frequency produces a 1.5% change in load. The governors on the responsive plant give \(1/R = 2000\) MW/Hz.

  1. Find \(D\) and \(\beta\).
  2. A 400 MW generator trips. Find the settled frequency and the split of the 400 MW between governor action and load relief.
  3. What breaks if every responsive unit is already at its ceiling when the trip occurs?
Solution

Turn the percentages into a slope. One per cent of frequency is 0.5 Hz; 1.5% of the load is 60 MW:

\[ D = \frac{\partial P_L}{\partial f} = \frac{0.015\times 4000\ \text{MW}}{0.01 \times 50\ \text{Hz}} = \frac{60}{0.5} = 120\ \text{MW/Hz} \]

Positive, and it must be: motors driving fans and pumps slow with the system and demand less. It is a helpful term, and it is the only term that survives when the governors cannot move.

Area stiffness.

\[ \beta = D + \frac{1}{R} = 120 + 2000 = 2120\ \text{MW/Hz} \]

The settled frequency after losing 400 MW:

\[ \Delta f = \frac{-400}{2120} = -0.18868\ \text{Hz} \qquad\Rightarrow\qquad f = 49.811\ \text{Hz} \]

Where the 400 MW comes from:

\[ \begin{array}{lcl} \text{Governor response} & -\Delta f/R = 2000\times0.18868 &= 377.36\ \text{MW} \\ \text{Load relief} & -D\,\Delta f = 120\times0.18868 &= 22.64\ \text{MW} \\ \hline & & 400.00\ \text{MW}\ \checkmark \end{array} \]

Load relief supplies 5.7% here — larger than in Problem 2, because this area's damping is a larger fraction of its stiffness. In a real interconnection the fraction is usually between 2% and 8%.

The failure case. If every governor valve is already wide open, \(1/R\) drops out of \(\beta\) entirely and only the load's own sensitivity is left:

\[ \Delta f = \frac{-400}{D} = \frac{-400}{120} = -3.33\ \text{Hz} \qquad\Rightarrow\qquad f = 46.67\ \text{Hz} \]

Which is not a number any system reaches, because long before it does, under-frequency load shedding removes demand that the generators cannot supply. Working backwards, damping alone can arrest the fall at a 48.8 Hz relay setting only for a disturbance of \(120 \times 1.2 = 144\) MW; anything larger sheds load.

The design conclusion. \(\beta\) falls by a factor of nearly eighteen when the governor term vanishes, so the whole of an area's frequency performance rests on keeping some plant off its ceiling. That is why responsive spinning reserve is scheduled explicitly and why an area running every unit flat out at peak is running without primary control, however much capacity is synchronised.

The load damping constant is the floor under an area's frequency response, and it is a very low floor. Its usual role is to contribute a few per cent and to be neglected in rough work. Its real importance is what it says about the case where governors saturate: at that point \(\beta\) collapses to \(D\), the frequency excursion multiplies by \(\beta/D\), and load shedding rather than generation becomes the balancing mechanism.
Answer\(D = 120\), \(\beta = 2120\) MW/Hz; \(f = \mathbf{49.811}\) Hz with 377.4 MW from governors and 22.6 MW of relief. With governors saturated, \(\Delta f = -3.33\) Hz and load shedding operates
Problem 6DrillGenerator–Load Block

A control area has an inertia constant \(H = 4\) s on a 3000 MVA base and a load damping of \(D = 0.008\) pu MW/Hz on the same base; \(f_0 = 50\) Hz.

  1. Find \(K_p\) and \(T_p\) of the generator–load block, and \(D\) in MW/Hz.
  2. A 180 MW generator trips. Find the initial rate of change of frequency, and the time the frequency would take to fall to 49.50 Hz if nothing responded at all.
  3. Repeat the RoCoF if half the synchronous plant has been displaced by inverter-connected generation, halving the effective inertia.
Solution

The block parameters follow directly from their definitions:

\[ K_p = \frac{1}{D} = \frac{1}{0.008} = 125\ \text{Hz per pu MW}, \qquad T_p = \frac{2H}{f_0 D} = \frac{2\times4}{50\times0.008} = \frac{8}{0.4} = 20\ \text{s} \]
\[ D = 0.008 \times 3000 = 24\ \text{MW/Hz} \]

Both parameters have plain meanings. \(K_p = 1/D\) says that a permanent 1 pu MW surplus, left entirely to the load's own sensitivity, would raise the frequency until the load had grown by that same amount. \(T_p\) is the ratio of stored kinetic energy to damping, and 20 s is a very slow block — which is precisely why a governor is needed.

The first instant after the trip. Set \(\Delta P_m = 0\) and \(\Delta f = 0\) in the swing equation: no valve has moved and no damping has developed, so the whole imbalance comes out of the rotors:

\[ \Delta P_L = \frac{180}{3000} = 0.06\ \text{pu}, \qquad \mathrm{RoCoF} = \frac{-f_0\,\Delta P_L}{2H} = \frac{-50\times0.06}{2\times4} = -0.375\ \text{Hz/s} \]

Note what does not appear: droop, damping, turbine time constants. The first derivative of frequency is decided by inertia and the size of the event, and by nothing else.

Time to a threshold on inertia alone, holding the initial slope:

\[ t = \frac{0.50\ \text{Hz}}{0.375\ \text{Hz/s}} = 1.33\ \text{s} \]

A bound, not a prediction — the real trace bends upward as governors act, as Problem 7 shows. But it is the right bound for the protection engineer: it says how long the governors have.

Halve the inertia. Nothing else in the calculation changes:

\[ \mathrm{RoCoF} = \frac{-50\times0.06}{2\times2} = -0.75\ \text{Hz/s}, \qquad t = \frac{0.50}{0.75} = 0.67\ \text{s} \]

Twice as fast, and half the time available. The settled frequency, however, is untouched: \(\Delta f_{ss} = -\Delta P_L/\beta\) contains no \(H\) at all. Inertia sets the speed of the excursion; \(\beta\) sets its destination.

Why RoCoF became a planning quantity. Rate-of-change-of-frequency relays fitted to embedded generation trip at settings of the order of 0.5 to 1.0 Hz/s. On the original system a 180 MW loss produces 0.375 Hz/s and nothing trips; on the low-inertia system the same event produces 0.75 Hz/s and a tranche of distributed generation may disconnect, enlarging the disturbance while it is still developing.

Inertia and stiffness answer two different questions, and confusing them is the standard error. "How fast?" is answered by \(H\) through the RoCoF; "how far, in the end?" is answered by \(\beta\) through \(-\Delta P_L/\beta\). "How far, at worst?" — the nadir — needs both, plus the governor lags, and is the subject of Problem 15.
Answer\(K_p = 125\) Hz/pu MW, \(T_p = 20\) s, \(D = 24\) MW/Hz; \(\mathrm{RoCoF} = \mathbf{-0.375}\) Hz/s reaching 49.5 Hz in 1.33 s, and −0.75 Hz/s in 0.67 s at half the inertia
Problem 7Exam levelPrimary Transient

The area of Problem 6 has an equivalent droop of \(R = 3\) Hz per pu MW. Neglect the governor and turbine lags.

  1. Find \(\beta\) in pu MW/Hz and in MW/Hz, and the settled frequency after the 180 MW trip.
  2. Find the time constant of the closed primary loop and the time to reach 95% of the final deviation.
  3. Show that the initial slope of this response agrees with the RoCoF of Problem 6, and explain why it must.
Solution

Stiffness. Work in per unit and convert at the end:

\[ \frac{1}{R} = \frac{1}{3} = 0.3333\ \text{pu MW/Hz}, \qquad \beta = D + \frac{1}{R} = 0.008 + 0.3333 = 0.34133\ \text{pu MW/Hz} \]
\[ \beta = 0.34133 \times 3000 = 1024\ \text{MW/Hz}, \qquad \frac{1}{R} = 1000\ \text{MW/Hz} \]
\[ \Delta f_{ss} = \frac{-180}{1024} = -0.17578\ \text{Hz} \qquad\Rightarrow\qquad f = 49.824\ \text{Hz} \]

The closed-loop time constant. With the two fast lags removed, the loop reduces to a single lag whose time constant is divided by the loop gain:

\[ 1 + \frac{K_p}{R} = 1 + \frac{125}{3} = 42.667, \qquad T_p' = \frac{T_p}{1+K_p/R} = \frac{20}{42.667} = 0.469\ \text{s} \]
\[ \Delta f(t) = -\frac{\Delta P_L}{\beta}\left(1 - e^{-t/T_p'}\right) = -0.1758\left(1-e^{-t/0.469}\right)\ \text{Hz} \]

The governor loop has cut the open-loop 20 s down to under half a second — a factor of 42.7. A system left to its own load damping would drift for the better part of a minute; with droop it is settled in about a second and a half.

Times off the exponential:

\[ \begin{array}{lcl} 63\%\ \text{of the deviation} & t = T_p' &= 0.47\ \text{s} \\ 95\%\ \text{of the deviation} & t = 3T_p' &= 1.41\ \text{s} \\ \Delta f\ \text{at } t = 1\ \text{s} & & -0.155\ \text{Hz} \end{array} \]

The consistency check, which is the useful part of this problem. Differentiate the response at \(t = 0\):

\[ \left.\frac{d(\Delta f)}{dt}\right|_{0^+} = -\frac{\Delta P_L}{\beta}\cdot\frac{1}{T_p'} = \frac{-0.17578}{0.46875} = -0.375\ \text{Hz/s} \]

Exactly the RoCoF of Problem 6, computed there from \(H\) alone with no mention of \(R\) or \(D\).

Why it must agree. At \(t = 0^+\) the valve has not moved and \(\Delta f\) is still zero, so neither the droop term nor the damping term contributes anything, and the swing equation reduces to the inertial statement. Algebraically:

\[ \frac{\Delta P_L}{\beta\,T_p'} = \frac{\Delta P_L}{\beta}\cdot\frac{1+K_p/R}{T_p} = \frac{\Delta P_L}{\beta}\cdot\frac{\beta/D}{2H/(f_0 D)} = \frac{f_0\,\Delta P_L}{2H} \]

Both \(\beta\) and \(D\) cancel identically. Any closed-loop model that gets this wrong has an error in it, which makes the check worth two lines of anybody's time.

What the approximation costs. Restoring \(T_g\) and \(T_T\) makes the characteristic polynomial cubic, and the two fast lags supply enough phase lag for the frequency to undershoot its final value and recover — a dip of roughly twice the settled offset over five to ten seconds. The single-lag model gets the destination exactly right and the worst point badly wrong, which is the distinction Problem 15 turns on.

The droop loop is a gain of forty-two acting on a twenty-second plant, and every consequence follows from that one number. It divides the settling time by 42.7, divides the steady-state error by 42.7 relative to damping alone, and — because a loop of that gain closed around three lags is not automatically well damped — it is also the reason the real trace overshoots. Primary control is a fast, strong, deliberately imperfect loop.
Answer\(\beta = 1024\) MW/Hz, \(f = \mathbf{49.824}\) Hz; \(T_p' = 0.469\) s with 95% reached at 1.41 s; the initial slope is −0.375 Hz/s, identical to the RoCoF because at \(t=0^+\) neither droop nor damping has acted
Problem 8Exam levelIntegral Gain

The supplementary loop \(\Delta P_{\text{ref}} = -K_I\!\int\!\Delta f\,dt\) is now closed around the area of Problems 6 and 7.

  1. Write the characteristic equation and identify \(\omega_n\) and \(\zeta\).
  2. Find \(K_I\) for critical damping and for \(\zeta = 0.7\).
  3. With \(\zeta = 0.7\), find the peak frequency dip after the 180 MW trip and compare it with the permanent offset droop alone would have left.
Solution

Close the loop. The supplementary command enters at the same summing point as the droop feedback, so the two feedback paths add:

\[ \Delta F(s) = \frac{-K_p\,\Delta P_L(s)}{(1+sT_p) + K_p\!\left(\dfrac{1}{R} + \dfrac{K_I}{s}\right)} = \frac{-s\,K_p\,\Delta P_L(s)}{T_p s^2 + \left(1+\dfrac{K_p}{R}\right)s + K_pK_I} \]

The factor \(s\) in the numerator is the whole argument for integral action written as algebra: for a step \(\Delta P_L(s) = \Delta P_L/s\) that \(s\) cancels the pole at the origin, so the final-value theorem returns zero. The steady-state frequency error is gone for any \(K_I > 0\).

Standard second-order form. Dividing through by \(T_p\):

\[ \omega_n = \sqrt{\frac{K_pK_I}{T_p}} = \sqrt{\frac{125\,K_I}{20}} = \sqrt{6.25\,K_I}, \qquad 2\zeta\omega_n = \frac{1+K_p/R}{T_p} = \frac{42.667}{20} = 2.1333\ \text{s}^{-1} \]

Only \(\omega_n\) contains \(K_I\); the damping term is fixed by the primary loop. So raising the integral gain speeds the restoration and lowers \(\zeta\) at the same time — one knob, two opposed effects.

The two gains. From \(\omega_n = 2.1333/(2\zeta)\) and \(K_I = \omega_n^2 T_p/K_p = \omega_n^2/6.25\):

\[ \begin{array}{lccc} \zeta & \omega_n\ (\text{rad/s}) & \omega_n/2\pi\ (\text{Hz}) & K_I \\ \hline 1.0 & 1.0667 & 0.170 & 0.1820 \\ 0.7 & 1.5238 & 0.243 & 0.3715 \end{array} \]

Relaxing from critical damping to \(\zeta = 0.7\) doubles the gain and buys a faster return at the cost of a small overshoot past nominal.

The transient with \(\zeta = 0.7\). Inverting the transfer function for a step, with \(\Delta P_L = 0.06\) pu:

\[ \Delta f(t) = -\frac{K_p\,\Delta P_L}{T_p}\cdot\frac{e^{-\zeta\omega_n t}}{\omega_d}\sin\omega_d t, \qquad \frac{K_p\,\Delta P_L}{T_p} = \frac{125\times0.06}{20} = 0.375\ \text{Hz/s} \]
\[ \omega_d = \omega_n\sqrt{1-\zeta^2} = 1.5238\times0.7141 = 1.0882\ \text{rad/s} \]

The leading coefficient 0.375 Hz/s is again the RoCoF — the initial slope cannot be changed by any controller, only by inertia.

The peak. Differentiating, the extremum is at \(\omega_d t = \cos^{-1}\zeta\):

\[ t_{\max} = \frac{\cos^{-1}(0.7)}{1.0882} = \frac{0.7954}{1.0882} = 0.731\ \text{s} \]
\[ \Delta f_{\max} = -\frac{0.375}{1.0882}\,e^{-0.7\times1.5238\times0.731}\sin(0.7954) = -0.1129\ \text{Hz} \]
\[ f_{\min} = 49.887\ \text{Hz}, \qquad f_{\text{final}} = 50.000\ \text{Hz} \]

The comparison. Droop alone (Problem 7) left the system permanently at 49.824 Hz. With the supplementary loop the worst point is 49.887 Hz and the final point is nominal:

\[ \begin{array}{lcc} & \text{worst frequency} & \text{final frequency} \\ \hline \text{Droop only} & 49.824 & 49.824 \\ \text{Droop} + K_I\ (\zeta=0.7) & 49.887 & 50.000 \end{array} \]

Better on both counts, which looks like a free lunch and is not. The integrator demands the full 180 MW from the plant inside a couple of seconds — faster than a steam turbine may be loaded — and a loop this fast will excite the inter-area modes of Problem 9. Real settings are several times smaller, and the practical target is that the bulk of the restoration takes ten to thirty seconds.

The two loops are proportional-plus-integral control with the terms deliberately separated in speed and in place. The proportional term must be fast and distributed, so that many governors can share one frequency without fighting; the integral term must be slow and central, so that the error is finally removed without exciting anything. Merging them — one fast integrator per unit — is exactly the isochronous failure of Problem 14.
Answer\(K_I = 0.182\) for \(\zeta = 1\) and \(\mathbf{K_I = 0.372}\) for \(\zeta = 0.7\); the dip is 0.113 Hz at 0.73 s and the frequency returns to 50.000 Hz, against a permanent 49.824 Hz on droop alone
Problem 9DrillThe Tie Line

Two areas are joined by a tie line of reactance \(X_{12} = 1.5\) pu on a 1000 MVA base. Both terminal voltages are 1.0 pu and the steady-state angle across the tie is 20°.

  1. Find the synchronizing coefficient \(T_{12}\) in pu MW/rad, in MW/rad and in MW per degree.
  2. With \(H_1 = H_2 = 8\) s on that base, find the inter-area oscillation frequency.
  3. A second identical tie line is built in parallel. What happens to \(T_{12}\) and to the oscillation frequency?
Solution

The coefficient is the slope of the power-angle curve at the operating point — the same quantity that governed small-signal stability in Set 25:

\[ T_{12} = \frac{|V_1||V_2|}{X_{12}}\cos\!\big(\delta_1^0-\delta_2^0\big) = \frac{1.0\times1.0}{1.5}\cos 20^\circ = 0.6667\times0.93969 = 0.62646\ \text{pu MW/rad} \]
\[ T_{12} = 0.62646\times1000 = 626.5\ \text{MW/rad} = 626.5\times\frac{\pi}{180} = 10.93\ \text{MW per degree} \]

The last form is the operationally useful one. At the stated 20° the tie is carrying \((1000/1.5)\sin20^\circ = 228\) MW, and one further degree adds \((1000/1.5)(\sin21^\circ-\sin20^\circ) = 10.9\) MW — the linearised coefficient, evaluated exactly.

Why the tie behaves as an integrator. Angle is the integral of frequency, and with \(\delta\) in electrical radians and \(f\) in hertz, \(d\delta/dt = 2\pi\,\Delta f\):

\[ \Delta P_{12}(s) = \frac{2\pi T_{12}}{s}\big[\Delta F_1(s) - \Delta F_2(s)\big] \]

So the flow keeps changing for as long as the two areas are at different frequencies, and settles only when \(\Delta f_1 = \Delta f_2\). That single statement is what "synchronous interconnection" means, and it is why the two-area steady state of Problem 10 has one frequency and not two.

The inter-area mode. Two inertias swinging against one spring:

\[ \omega_{\text{osc}} = \sqrt{2\pi f_0 T_{12}\left(\frac{1}{2H_1}+\frac{1}{2H_2}\right)} = \sqrt{314.16\times0.62646\times\left(\tfrac{1}{16}+\tfrac{1}{16}\right)} \]
\[ = \sqrt{314.16\times0.62646\times0.125} = \sqrt{24.60} = 4.96\ \text{rad/s} \qquad\Rightarrow\qquad f_{\text{osc}} = \frac{4.96}{2\pi} = 0.789\ \text{Hz} \]

Within the 0.1–0.8 Hz band that characterises inter-area modes. They are far slower than the local machine modes of Set 25 because the effective inertias here are those of whole areas rather than of single rotors.

The second circuit. Two equal reactances in parallel halve \(X_{12}\), so:

\[ T_{12}' = 2\times0.62646 = 1.2529\ \text{pu MW/rad}, \qquad \omega_{\text{osc}}' = \sqrt2 \times 4.96 = 7.01\ \text{rad/s} \]
\[ f_{\text{osc}}' = 1.116\ \text{Hz} \]

A stiffer spring between unchanged masses: the mode rises by \(\sqrt2\). The interconnection is stronger in every steady-state sense — more transfer capability, a smaller angle for the same flow — and its oscillation is faster and, in general, better damped by the same power system stabilisers.

The consequence for control. Any supplementary controller must be slow compared with the inter-area mode, or it will excite it. At 0.79 Hz the mode period is 1.27 s, so an AGC loop with \(\omega_n = 1.52\) rad/s (0.24 Hz, from Problem 8) sits comfortably below it — but only by a factor of three, which is the practical reason AGC gains are set several times lower than the critically damped value.

A tie line has two quite different roles in this chapter and it is worth keeping them apart. In the steady-state problems it is a bookkeeping quantity — the flow that must return to schedule. In the dynamic problems it is a spring — the coupling that sets the inter-area mode and thereby the ceiling on every controller gain in both areas. \(T_{12}\) is the same number in both roles.
Answer\(T_{12} = 0.6265\) pu MW/rad \(= 626.5\) MW/rad \(= 10.93\) MW/degree; \(f_{\text{osc}} = \mathbf{0.789}\) Hz, rising to 1.116 Hz with a second parallel tie
Problem 10Exam levelTwo Areas on Droop

Two interconnected areas at 50 Hz. Area 1: 3000 MW of responsive capacity at 5% droop, \(D_1 = 30\) MW/Hz. Area 2: 5000 MW at 4% droop, \(D_2 = 45\) MW/Hz. The tie carries its scheduled flow. Supplementary control is out of service in both areas.

A 150 MW load increase occurs in area 1. Find the settled frequency, the change in tie-line flow with its sign, and the megawatts each area supplies. Compare with what area 1 would have suffered standing alone.

Solution

The two stiffnesses.

\[ \frac{1}{R_1} = \frac{3000}{0.05\times50} = 1200, \qquad \beta_1 = 1200 + 30 = 1230\ \text{MW/Hz} \]
\[ \frac{1}{R_2} = \frac{5000}{0.04\times50} = 2500, \qquad \beta_2 = 2500 + 45 = 2545\ \text{MW/Hz} \]

Set up the steady state. The tie has stopped changing, so the two areas are at one frequency. Write each area's power balance, taking \(\Delta P_{12}\) positive for flow out of area 1:

\[ \text{Area 1:}\quad -\frac{\Delta f}{R_1} - \Delta P_{L1} - \Delta P_{12} = D_1\Delta f \]
\[ \text{Area 2:}\quad -\frac{\Delta f}{R_2} + \Delta P_{12} = D_2\Delta f \]

The second equation alone gives \(\Delta P_{12} = \beta_2\,\Delta f\), which is the useful shortcut: area 2 has no disturbance of its own, so everything it does appears on the tie.

Solve.

\[ \Delta f = \frac{-\Delta P_{L1}}{\beta_1+\beta_2} = \frac{-150}{1230+2545} = \frac{-150}{3775} = -0.03974\ \text{Hz} \qquad\Rightarrow\qquad f = 49.960\ \text{Hz} \]
\[ \Delta P_{12} = \beta_2\,\Delta f = 2545\times(-0.03974) = -101.13\ \text{MW} \]

Negative, with the stated sign convention, means the flow is into area 1: area 2 is sending 101.13 MW to help a neighbour whose load it does not serve.

Who supplies what.

\[ \begin{array}{lccc} & \text{governors} & \text{load relief} & \text{total} \\ \hline \text{Area 1} & 1200(0.03974)=47.68 & 30(0.03974)=1.19 & 48.87\ \text{MW} \\ \text{Area 2} & 2500(0.03974)=99.34 & 45(0.03974)=1.79 & 101.13\ \text{MW} \\ \hline & & & 150.00\ \text{MW}\ \checkmark \end{array} \]

Area 1, where the load actually appeared, supplies only 32.6% of it. The split is \(\beta_1 : \beta_2\) and has nothing to do with which area is at fault.

Standing alone, area 1 would have had only its own stiffness:

\[ \Delta f_{\text{alone}} = \frac{-150}{1230} = -0.1220\ \text{Hz} \quad\text{against}\quad -0.0397\ \text{Hz interconnected} \]

Three times the excursion. This is the benefit that motivates interconnection in the first place, and it is delivered automatically by the physics, with no controller and no communication of any kind.

And the problem it creates. Area 2's generators have moved off their economic setpoints, area 2 is exporting 101 MW it never contracted to sell, and the tie flow no longer matches schedule. Governor support across a tie is exactly what one wants for the first few seconds and exactly what one does not want after a minute — which is the entire brief for the supplementary loop.

AREA 1 AREA 2 beta_1 = 1230 MW/Hz beta_2 = 2545 MW/Hz dP_L1 = +150 MW no disturbance supplies 48.87 MW supplies 101.13 MW X_12 101.13 MW dP_12 = -101.13 MW ACE_1 = dP_12 + B_1 df ACE_2 = -dP_12 + B_2 df = -150 MW = 0 one frequency: df = -0.0397 Hz
The two-area steady state on droop alone, with the area control errors that Problem 11 computes from it. Ideal bias makes the innocent area's error vanish.
On primary control the interconnection socialises every disturbance in proportion to stiffness, and it does so whether or not that is wanted. The frequency benefit is real and immediate; the commercial and dispatch consequences are unacceptable if left standing. Everything from Problem 11 onward exists to keep the first and undo the second.
Answer\(f = \mathbf{49.960}\) Hz; \(\Delta P_{12} = \mathbf{-101.13}\) MW, i.e. area 2 imports into area 1; area 1 supplies only 48.87 MW of its own 150 MW, and its excursion is a third of the 0.122 Hz it would have suffered alone
Problem 11DrillArea Control Error

For the state reached in Problem 10, both areas set their frequency bias equal to their own area frequency response characteristic: \(B_1 = 1230\), \(B_2 = 2545\) MW/Hz.

  1. Compute \(\mathrm{ACE}_1\) and \(\mathrm{ACE}_2\).
  2. Show that the only possible steady state with both integrators in service is \(\Delta f = 0\), \(\Delta P_{12} = 0\), and that this does not depend on the value of either bias.
Solution

Evaluate the two errors from the numbers already obtained, \(\Delta f = -0.03974\) Hz and \(\Delta P_{12} = -101.13\) MW, with export counted positive:

\[ \mathrm{ACE}_1 = \Delta P_{12} + B_1\Delta f = -101.13 + 1230(-0.03974) = -101.13 - 48.87 = -150.00\ \text{MW} \]
\[ \mathrm{ACE}_2 = -\Delta P_{12} + B_2\Delta f = +101.13 + 2545(-0.03974) = 101.13 - 101.13 = 0 \]

Read what those two numbers say. Area 1's controller sees exactly its own 150 MW deficit — not the 48.87 MW its governors happened to supply, and not the 101.13 MW that appeared on the tie, but the disturbance itself. Area 2's controller sees zero, and therefore does not move a single speed changer, even though its generators are at that moment producing an extra 101.13 MW.

That extra output is withdrawn automatically as area 1's integrator raises its own references and the frequency comes back to nominal. Area 2 lent support for the first seconds and is repaid without ever having been told what happened.

The general result behind it. With \(B_i = \beta_i\) and a disturbance \(\Delta P_{L1}\) in area 1:

\[ \mathrm{ACE}_1 = \frac{-\beta_2\Delta P_{L1}}{\beta_1+\beta_2} + \beta_1\!\left(\frac{-\Delta P_{L1}}{\beta_1+\beta_2}\right) = -\Delta P_{L1} \]
\[ \mathrm{ACE}_2 = \frac{+\beta_2\Delta P_{L1}}{\beta_1+\beta_2} + \beta_2\!\left(\frac{-\Delta P_{L1}}{\beta_1+\beta_2}\right) = 0 \]

Identically, for any \(\beta_1\) and \(\beta_2\). This is what "non-interactive control" means: the bias setting encodes the neighbour's expected response so exactly that the tie term and the frequency term cancel in every undisturbed area.

The equilibrium argument — and it needs no numbers at all. Each \(\Delta P_{\text{ref},i}\) is the output of an integrator driven by \(\mathrm{ACE}_i\). In a steady state every signal is constant; an integrator with a constant non-zero input produces a ramp, not a constant; so the only steady state that can exist has \(\mathrm{ACE}_1 = \mathrm{ACE}_2 = 0\):

\[ \Delta P_{12} + B_1\Delta f = 0, \qquad -\Delta P_{12} + B_2\Delta f = 0 \]
\[ \text{adding:}\quad (B_1+B_2)\,\Delta f = 0 \;\Longrightarrow\; \Delta f = 0 \;\Longrightarrow\; \Delta P_{12} = 0 \]

Both objectives are met for any positive pair of biases, and the biases have vanished from the answer. The final state does not depend on them at all.

So what does the bias decide? The path, not the destination. With \(B_i = \beta_i\) the innocent area is never disturbed and the whole correction is made by the area that caused it. With any other setting the innocent area's controller is given a non-zero error and starts moving plant for no reason, as Problem 12 shows in numbers.

The area control error is a single scalar carrying two objectives, and the trick is that a well-chosen weight makes it also carry the answer to "whose fault was it?". Neither term alone would do: frequency alone would make every area respond to every disturbance for ever, and tie flow alone would leave the frequency permanently offset if both areas were in balance with each other but jointly short of generation.
Answer\(\mathrm{ACE}_1 = \mathbf{-150}\) MW — exactly the disturbance — and \(\mathrm{ACE}_2 = \mathbf{0}\). The equilibrium \(\Delta f = 0\), \(\Delta P_{12} = 0\) holds for any positive biases
Problem 12Exam levelA Mis-Set Bias

The same two areas as Problem 10, but now a 250 MW load increase occurs in area 2. Supplementary control is in service.

  1. Find \(\Delta f\) and \(\Delta P_{12}\) immediately after primary control settles.
  2. Compute both ACEs with ideal biases.
  3. Repeat \(\mathrm{ACE}_1\) with \(B_1\) mistakenly set to 600 MW/Hz, and again with \(B_1 = 2000\) MW/Hz. Describe what area 1's controller does in each case.
Solution

The primary steady state, by symmetry with Problem 10 — the disturbed area is now the second one:

\[ \Delta f = \frac{-\Delta P_{L2}}{\beta_1+\beta_2} = \frac{-250}{3775} = -0.06623\ \text{Hz} \qquad\Rightarrow\qquad f = 49.934\ \text{Hz} \]
\[ \Delta P_{12} = -\beta_1\,\Delta f = -1230\times(-0.06623) = +81.46\ \text{MW} \]

Positive this time: area 1 exports 81.46 MW to help its neighbour. Area 2 supplies \(2545\times0.06623 = 168.54\) MW itself, and \(81.46 + 168.54 = 250\) MW. ✓

With ideal biases \(B_1 = 1230\), \(B_2 = 2545\) MW/Hz:

\[ \mathrm{ACE}_1 = 81.46 + 1230(-0.06623) = 81.46 - 81.46 = 0 \]
\[ \mathrm{ACE}_2 = -81.46 + 2545(-0.06623) = -81.46 - 168.54 = -250.00\ \text{MW} \]

The mirror image of Problem 11, as it must be. The area that caused the disturbance owns the whole of it; the innocent area is left alone.

Bias set too low: \(B_1 = 600\) MW/Hz, less than half its correct value.

\[ \mathrm{ACE}_1 = 81.46 + 600(-0.06623) = 81.46 - 39.74 = +41.72\ \text{MW} \]

A positive ACE is read as over-generation. Area 1's controller therefore starts winding its own references down by up to 41.7 MW while area 2 is still trying to make up a 250 MW shortfall. The two controllers work against each other, the frequency dips further before recovering, and both sets of turbines are moved further than necessary.

Bias set too high: \(B_1 = 2000\) MW/Hz.

\[ \mathrm{ACE}_1 = 81.46 + 2000(-0.06623) = 81.46 - 132.45 = -50.99\ \text{MW} \]

Now the error is negative, so area 1 raises its own generation by about 51 MW that nobody asked it for. Frequency is pushed back through nominal from below, area 2's integrator has to unwind some of its own correction, and the pair hunt briefly before settling.

Collect the three cases:

\[ \begin{array}{lccl} B_1\ (\text{MW/Hz}) & \mathrm{ACE}_1\ (\text{MW}) & \text{sign read as} & \text{area 1 does} \\ \hline 600 & +41.72 & \text{over-generation} & \text{reduces output, wrongly} \\ 1230 & 0 & \text{balance} & \text{nothing} \\ 2000 & -50.99 & \text{deficit} & \text{raises output, wrongly} \end{array} \]

Neither error is unstable. Problem 11 showed that \(\Delta f = 0\), \(\Delta P_{12} = 0\) is the only equilibrium for any positive bias, so all three cases end in the same place. What differs is the regulating energy spent getting there and the depth of the transient.

The asymmetry that matters in practice. Of the two errors, the under-set bias is worse. It makes the innocent area actively oppose the correction, so the frequency recovery is slowed at the moment it matters most. That is why grid codes specify a minimum bias — commonly a stated fraction of the measured \(\beta\) — rather than a tolerance band, and why the measurement of Problem 4 is repeated annually.

An area's bias is a statement about its own governors, and it must be re-measured, not re-declared. Problem 4's area believed \(\beta = 3025\) MW/Hz and actually had 1471. A bias set from the belief would be more than double the truth — the over-set case above, applied permanently, to every disturbance in the interconnection.
Answer\(f = 49.934\) Hz, \(\Delta P_{12} = +81.46\) MW; ideal bias gives \(\mathrm{ACE}_1 = 0\), \(\mathrm{ACE}_2 = -250\) MW. With \(B_1 = 600\), \(\mathrm{ACE}_1 = \mathbf{+41.7}\) MW and area 1 wrongly reduces generation; with \(B_1 = 2000\), \(\mathrm{ACE}_1 = \mathbf{-51.0}\) MW and it wrongly raises it
Problem 13Exam levelParticipation Factors

An area contains three units of 250, 200 and 150 MW, all on 4% droop, with load damping 10 MW/Hz. Their fuel-cost curves are

\[ C_1 = 0.0045P_1^2 + 5.2P_1 + 300, \quad C_2 = 0.0060P_2^2 + 4.8P_2 + 250, \quad C_3 = 0.0090P_3^2 + 6.0P_3 + 180 \]

in rupees per hour with \(P\) in MW. A 90 MW load increase occurs.

  1. Find the frequency and the output change of each unit after primary control settles.
  2. Find the participation factors and the final allocation once AGC has restored 50 Hz.
  3. State which units are unloaded and which are loaded further during the second stage.
Solution

Stage one — the droop allocation. Convert, add, divide:

\[ \frac{1}{R_1} = \frac{250}{2} = 125, \quad \frac{1}{R_2} = \frac{200}{2} = 100, \quad \frac{1}{R_3} = \frac{150}{2} = 75\ \text{MW/Hz} \]
\[ \sum\frac{1}{R_i} = 300, \qquad \beta = 10 + 300 = 310\ \text{MW/Hz}, \qquad \Delta f = \frac{-90}{310} = -0.29032\ \text{Hz} \]
\[ f = 49.710\ \text{Hz} \]

The primary shares:

\[ \Delta P_1 = 125(0.29032) = 36.29, \quad \Delta P_2 = 100(0.29032) = 29.03, \quad \Delta P_3 = 75(0.29032) = 21.77\ \text{MW} \]
\[ 36.29+29.03+21.77 = 87.10\ \text{MW}, \qquad -D\,\Delta f = 2.90\ \text{MW}, \qquad \text{total } 90.00\ \checkmark \]

A split of 125 : 100 : 75, which is 250 : 200 : 150 — the ratings, because the percentage droops are equal. Cost has played no part whatever.

Stage two — the participation factors. They come from the second derivatives of the cost curves, exactly as Set 27 derived them:

\[ \frac{d^2C_1}{dP_1^2} = 0.0090, \qquad \frac{d^2C_2}{dP_2^2} = 0.0120, \qquad \frac{d^2C_3}{dP_3^2} = 0.0180 \]
\[ \frac{1}{0.0090} = 111.11, \quad \frac{1}{0.0120} = 83.33, \quad \frac{1}{0.0180} = 55.56, \qquad \sum = 250.00 \]
\[ \alpha_1 = \frac{111.11}{250} = 0.4444, \qquad \alpha_2 = 0.3333, \qquad \alpha_3 = 0.2222 \]

The flattest incremental-cost curve takes the largest share, because a unit whose incremental cost rises slowly can absorb more megawatts before it becomes the expensive one.

The final allocation. AGC drives \(\Delta f\) to zero, so the droop terms all vanish and the load relief disappears with them. The units must supply the entire 90 MW from their reference settings:

\[ \Delta P_1 = 0.4444\times90 = 40.0, \quad \Delta P_2 = 0.3333\times90 = 30.0, \quad \Delta P_3 = 0.2222\times90 = 20.0\ \text{MW} \]
\[ \Delta P_{\text{ref},i} = \Delta P_i \quad\text{because } \Delta f = 0 \text{ removes the droop contribution entirely} \]

Compare the two stages, which is the object of the problem:

\[ \begin{array}{lccc} \text{Unit} & \text{after droop (MW)} & \text{after AGC (MW)} & \text{movement} \\ \hline 1 & 36.29 & 40.00 & +3.71 \\ 2 & 29.03 & 30.00 & +0.97 \\ 3 & 21.77 & 20.00 & -1.77 \\ \hline \text{sum} & 87.10 & 90.00 & +2.90 \end{array} \]

Unit 3 is unloaded during the second stage even though the system's total generation is rising, because it is the expensive machine at the margin. The extra 2.90 MW in the totals is the load relief being repaid as the frequency comes back to nominal — the load grows again on the way up.

The general shape. Two allocations run in sequence after every disturbance, and they answer different questions:

\[ \begin{array}{lll} \text{Seconds} & \text{divide by } 1/R_i & \text{technical: who can respond} \\ \text{Minutes} & \text{divide by } \alpha_i & \text{economic: who should} \\ \end{array} \]

The transient state is the physical one and the final state is the economic one. Units are unloaded and reloaded against each other while the total stays fixed, and that reshuffling is invisible on the frequency trace.

Droop and participation factors are two different answers to "how much of this does each unit take?", and both are right on their own time scale. Droop cannot be economic because a governor knows nothing about fuel prices; participation factors cannot be fast because an economic dispatch takes minutes to compute and turbines take minutes to move. Set 27 supplies the second; this chapter supplies the first; AGC is the handover between them.
AnswerDroop gives 36.29 / 29.03 / 21.77 MW at 49.710 Hz; AGC gives \(\alpha = 0.444/0.333/0.222\) and a final 40 / 30 / 20 MW at 50.000 Hz — unit 3 is unloaded by 1.77 MW during the changeover
Problem 14DrillIsochronous Control

Two 300 MW units share an isolated system. Each governor's speed reference is set by an instrument with an error of ±0.02 Hz.

  1. Show what happens if both governors are isochronous — that is, integral acting, holding speed exactly.
  2. With 4% droop on both instead, find the megawatt error in load sharing that the same 0.02 Hz reference error produces.
  3. State the general principle.
Solution

The isochronous case. An isochronous governor drives its valve until its own frequency error is zero, so its steady-state characteristic is a horizontal line at whatever frequency its instrument believes to be 50.00 Hz:

\[ \text{Unit A holds } 50.00\ \text{Hz}, \qquad \text{Unit B holds } 50.02\ \text{Hz} \]

These two demands cannot both be met. The system runs at one frequency, so at least one governor sees a non-zero error and keeps integrating.

Follow it through. Suppose the system settles anywhere below 50.02 Hz. Unit B sees a negative error and keeps raising its valve; unit A, seeing the frequency now above its own 50.00 Hz target, keeps lowering its valve. Neither stops:

\[ P_B \to 300\ \text{MW (its ceiling)}, \qquad P_A \to 0\ \text{MW (or reverse power, and trip)} \]

The load split is decided entirely by which instrument reads high, and by nothing else. A 0.02 Hz calibration difference — 0.04% — produces a 300 MW difference in output. Two integrators cannot share one controlled variable.

The droop case. A proportional characteristic converts the same reference error into a finite, small offset. With 4% on a 300 MW machine:

\[ \frac{1}{R} = \frac{300}{0.04\times50} = 150\ \text{MW/Hz} \]
\[ \Delta P_{\text{error}} = \frac{1}{R}\times 0.02 = 150\times0.02 = 3.0\ \text{MW} \]

Three megawatts out of three hundred — one per cent, and stable. The units share the load in the intended proportion plus a small bias, and neither runs away.

The scaling. The sharing error is \(\Delta f_{\text{instr}}/R\), so it grows as the droop is tightened:

\[ \begin{array}{lcc} \text{droop} & 1/R\ (\text{MW/Hz}) & \text{error from } 0.02\ \text{Hz} \\ \hline 8\% & 75 & 1.5\ \text{MW} \\ 4\% & 150 & 3.0\ \text{MW} \\ 2\% & 300 & 6.0\ \text{MW} \\ 0\% & \infty & \text{full range} \end{array} \]

Zero droop is the isochronous case, and the table shows it as the limit of a continuous trend rather than as a separate phenomenon. A stiffer unit is more responsive and more sensitive to reference error at the same time.

The principle, and it is not confined to power systems: a shared variable can be regulated by at most one integrator. Droop is the price of cooperation — it guarantees the permanent frequency offset of Problems 2 to 7, and in exchange it lets an unlimited number of governors act on one measurement without any of them needing to know the others exist.

Where isochronous control is still used: a single machine feeding an isolated load — a ship, an emergency generator, an island with one unit — where the sharing problem does not arise. It is also used deliberately on exactly one unit of a small island system, with every other unit on droop, so that one machine takes all the residual error. That is the same architecture as AGC, with the integrator placed at a machine rather than at a control centre.

The permanent frequency error of primary control is not a defect to be apologised for; it is the mechanism that makes primary control possible. Remove it and the system has no way of deciding how to divide a load increment between two willing machines. The error is then removed later, once, centrally, by the supplementary loop — which is a single integrator for the whole area, exactly as the principle demands.
AnswerIsochronous: one unit runs to full load and the other to zero, from a 0.02 Hz instrument error. With 4% droop the same error costs only \(150\times0.02 = \mathbf{3.0}\) MW — one per cent of rating
Problem 15ChallengeNadir vs Settling

An area carries 8000 MW with \(H = 5\) s on a 10 000 MVA base and \(\beta = 1600\) MW/Hz. An 800 MW unit trips. Model the governor response as a ramp: the mechanical power rises linearly from zero at the instant of the trip to the full 800 MW at \(t = T_r = 5\) s, and neglect load damping during the transient. Under-frequency load shedding is armed at 48.8 Hz.

  1. Derive the frequency nadir from this model and evaluate it.
  2. Compare it with the settled frequency. Does load shedding operate?
  3. Repeat with half the inertia, everything else unchanged, and comment.
Solution

Set up the swing equation with the ramp. Writing the area's stored energy as \(2HS/f_0\) with \(S\) the MVA base:

\[ \frac{2HS}{f_0}\frac{d(\Delta f)}{dt} = \Delta P_m(t) - \Delta P_L, \qquad \Delta P_m(t) = \Delta P_L\frac{t}{T_r}\ \ (t \le T_r) \]
\[ \frac{d(\Delta f)}{dt} = -K\left(1 - \frac{t}{T_r}\right), \qquad K \equiv \frac{f_0\,\Delta P_L}{2HS} = \mathrm{RoCoF} \]

Integrate. The deficit shrinks linearly and vanishes at \(T_r\), which is therefore where the frequency stops falling:

\[ \Delta f(t) = -K\left(t - \frac{t^2}{2T_r}\right) \qquad\Rightarrow\qquad \Delta f_{\text{nadir}} = \Delta f(T_r) = -\frac{K\,T_r}{2} \]

A compact and useful result: the nadir is the initial RoCoF multiplied by half the response time. Half, because the deficit is triangular — the area under a triangle, not a rectangle.

Evaluate for the stated area.

\[ K = \frac{50\times800}{2\times5\times10\,000} = \frac{40\,000}{100\,000} = 0.400\ \text{Hz/s} \]
\[ \Delta f_{\text{nadir}} = -0.400\times\frac{5}{2} = -1.00\ \text{Hz} \qquad\Rightarrow\qquad f_{\text{nadir}} = 49.00\ \text{Hz} \]

The settled frequency, which involves none of the same quantities:

\[ \Delta f_{ss} = \frac{-800}{1600} = -0.500\ \text{Hz} \qquad\Rightarrow\qquad f = 49.50\ \text{Hz} \]

The transient minimum is twice the permanent offset. Load shedding at 48.8 Hz does not operate — but the margin is 0.2 Hz, not the 0.7 Hz that the settling point alone would suggest.

Halve the inertia. Nothing in \(\beta\) changes, so the destination is untouched; only the path changes:

\[ K' = \frac{50\times800}{2\times2.5\times10\,000} = 0.800\ \text{Hz/s}, \qquad \Delta f_{\text{nadir}} = -0.800\times2.5 = -2.00\ \text{Hz} \]
\[ f_{\text{nadir}} = 48.00\ \text{Hz} \quad\text{against an unchanged } f_{ss} = 49.50\ \text{Hz} \]

Load shedding now operates and sheds a block of demand the system did not actually need to lose — the generators were always going to cover the deficit by \(t = 5\) s. The event is made worse by the protection, and the settled-state calculation gave no warning at all.

What would have to change to survive it. Requiring the nadir to stay above 48.8 Hz means \(K'T_r/2 \le 1.2\):

\[ T_r \le \frac{2\times1.2}{0.800} = 3.0\ \text{s} \]

The governors would have to deliver the whole 800 MW in three seconds instead of five. That is achievable for hydro and open-cycle gas, marginal for a non-reheat steam unit and out of the question for a reheat unit (Problem 16). The alternative levers are more inertia, a faster-than-governor resource such as a battery, or a lower shedding threshold.

The two traces make the point better than the algebra does. Same disturbance, same \(\beta\), same destination, different worst point:

50.0 49.5 48.8 48.0 0 5 10 15 20 time (s) f (Hz) H = 5 s nadir 49.0 H = 2.5 s nadir 48.0 settle 49.5 UFLS 48.8 Hz
Same loss, same area stiffness, same settling point — and one trace sheds load while the other does not. Inertia decides the depth of the dip, not its destination.
The examination asks for the settling point; the protection responds to the nadir; and the two are governed by different quantities. \(\beta\) alone sets where the frequency ends; \(H\), the governor and turbine lags, and \(\beta\) together set how far it dips on the way. Two systems with identical \(\beta\) can have nadirs a whole hertz apart, and only one of them keeps its load.
Answer\(\Delta f_{\text{nadir}} = -K T_r/2\); here 49.00 Hz against a settled 49.50 Hz, and load shedding does not operate. At half the inertia the nadir is 48.00 Hz and it does, with the settled frequency unchanged
Problem 16ChallengeThe Reheat Unit

A reheat steam turbine has \(G_T(s) = \dfrac{1+sK_rT_r}{(1+sT_T)(1+sT_r)}\) with \(K_r = 0.3\), \(T_r = 8\) s and \(T_T = 0.3\) s. A non-reheat unit on the same site has \(G_T(s) = 1/(1+0.3s)\).

Find the step response of the reheat unit, evaluate the fraction of its final response delivered at 1, 2 and 5 seconds, compare with the non-reheat unit, and state the consequence for the frequency nadir of Problem 15.

Solution

Expand in partial fractions. For a unit step, \(\Delta P_m(s) = G_T(s)/s\) with poles at \(0\), \(-1/T_T\) and \(-1/T_r\):

\[ \Delta P_m(t) = 1 + A\,e^{-t/T_T} + B\,e^{-t/T_r} \]
\[ A = \frac{1 - K_rT_r/T_T}{T_r/T_T - 1} = \frac{1-8}{25.667} = -0.27273, \qquad B = -\frac{1-K_r}{1-T_T/T_r} = -\frac{0.7}{0.9625} = -0.72727 \]

Check the sum: \(1 - 0.27273 - 0.72727 = 0\), so \(\Delta P_m(0) = 0\) as a strictly proper transfer function requires. The two residues are the two halves of the response, and their relative size is the whole story.

Read the structure. The fast term dies in about a second; after that the response is

\[ \Delta P_m(t) \approx 1 - 0.72727\,e^{-t/8} \qquad (t \gtrsim 1\ \text{s}) \]
\[ \text{prompt fraction} = 1 + B = 1 - \frac{1-K_r}{1-T_T/T_r} = 0.2727 \;\approx\; K_r \]

So the high-pressure cylinder delivers about 27% of the eventual response within a second, and the remaining 73% arrives with the reheater's 8 s time constant. The fraction \(K_r\) is exactly the share of the total power developed in the HP cylinder, which is what the parameter means physically.

Tabulate against the non-reheat unit:

\[ \begin{array}{lcc} t\ (\text{s}) & \text{reheat} & \text{non-reheat} \\ \hline 0.5 & 0.265 & 0.811 \\ 1 & 0.349 & 0.964 \\ 2 & 0.433 & 0.999 \\ 5 & 0.611 & 1.000 \\ 10 & 0.792 & 1.000 \\ 20 & 0.940 & 1.000 \end{array} \]
\[ \text{reheat reaches } 90\% \text{ at } t = -8\ln\!\left(\frac{0.10}{0.72727}\right) = 15.9\ \text{s} \]

The two curves on one axis:

1.0 0.5 K_r 0 0 5 10 15 20 time (s) dP_m (pu) non-reheat, T_T = 0.3 s reheat, T_r = 8 s
Step response of the two turbine models. The reheat unit's prompt contribution is limited to its high-pressure fraction; the rest arrives long after the frequency nadir has passed.

The consequence for the nadir. Problem 15's ramp model assumed the full response arrived by \(T_r = 5\) s. A fleet of reheat units delivers only 61% by then, so the effective ramp is much longer and the nadir correspondingly deeper. Crudely, replacing \(T_r = 5\) s by an effective 12 s in \(-KT_r/2\) multiplies the dip by 2.4.

The settling point is untouched — the reheat unit does eventually deliver its whole \(1/R\) contribution, so it appears in \(\beta\) at full value. It is only the nadir that suffers, which is exactly the quantity \(\beta\) does not describe.

The system-design conclusion. Two areas with identical \(\beta\) and identical inertia can have nadirs a hertz apart if one is rich in hydro and open-cycle gas and the other in reheat steam. Frequency-response markets have accordingly moved from paying for \(\beta\) to paying for delivery within a stated window — typically the megawatts sustained between 1 and 10 s. A reheat unit sells very little of that product, and a battery sells almost nothing else.

The reheater is a 5–10 second energy store in the middle of the control loop, and it is invisible in every steady-state calculation in this set. Its effect is to split the unit's governor response into a prompt fraction \(K_r\) that helps with the nadir and a delayed remainder that arrives only in time for the settling point. Both fractions count towards \(\beta\); only the first counts towards keeping the load-shedding relays quiet.
Answer\(\Delta P_m(t) = 1 - 0.2727e^{-t/0.3} - 0.7273e^{-t/8}\): 34.9% at 1 s, 43.3% at 2 s, 61.1% at 5 s against 96%, 99.9% and 100% for the non-reheat unit. The prompt fraction is \(K_r \approx 0.27\), and the nadir suffers while \(\beta\) does not
Problem 17ChallengeThe Hydro Governor

A hydro turbine is represented by \(G_T(s) = \dfrac{1-sT_w}{1+0.5\,sT_w}\) with a water starting time \(T_w = 1.0\) s.

  1. Find and interpret the response to a unit step increase in gate opening.
  2. Explain physically why the power initially moves the wrong way.
  3. Explain why a hydro governor must be given transient droop, and what that does to its contribution during the first seconds of a disturbance.
Solution

The step response. With \(T_w = 1\):

\[ \Delta P_m(s) = \frac{1-s}{s(1+0.5s)} = \frac{1}{s} - \frac{1.5}{1+0.5s} = \frac{1}{s} - \frac{3}{s+2} \]
\[ \Delta P_m(t) = 1 - 3e^{-2t/T_w} \]

Evaluate it, and the shape is startling:

\[ \begin{array}{lc} t\ (\text{s}) & \Delta P_m\ (\text{pu}) \\ \hline 0 & -2.00 \\ 0.25 & -0.82 \\ 0.549 & 0 \\ 1.0 & +0.59 \\ 2.0 & +0.945 \\ 3.0 & +0.993 \end{array} \]

Opening the gate to demand one more per-unit of power produces two per-unit less at the first instant. The power does not even change sign until \(t = \tfrac{T_w}{2}\ln 3 = 0.55\) s, and does not reach 90% of its target until 1.70 s.

The physics. The power developed is \(P \propto Hq\) — head times flow. Opening the gate raises the flow area instantly, but the water column in the penstock has mass and cannot accelerate instantly:

\[ \begin{array}{ll} \text{Instant } 0^+ & \text{gate area up, flow } q \text{ unchanged} \Rightarrow \text{velocity through the gate falls} \\ & \Rightarrow \text{head at the runner falls} \Rightarrow P \text{ falls} \\ \text{Over } \sim T_w & \text{the column accelerates, } q \text{ rises, head recovers} \Rightarrow P \text{ rises to target} \end{array} \]

The same effect works in reverse on closing: shutting a gate momentarily raises the power, and it is the reason a penstock needs a surge tank. The right-half-plane zero at \(s = +1/T_w\) is the transfer-function statement of a mechanical fact about water.

Why the governor must be detuned. A right-half-plane zero contributes phase lag while contributing gain. Close a fast proportional loop around it and the loop responds to the initial wrong-way power by demanding more gate, which makes the power fall further — a positive feedback that shows up as sustained hunting. The remedy is transient droop: a governor whose regulation is large immediately after a movement and relaxes to the permanent value as the water column settles:

\[ R(s) = R_p + (R_t - R_p)\frac{sT_R}{1+sT_R}, \qquad R_t \approx 5R_p,\ \ T_R \approx 5T_w \]

Typical values are a permanent droop of 4% and a transient droop of 20–40% with a reset time of several seconds. The governor is deliberately made five to ten times weaker on the time scale that matters for stability.

What that costs in frequency response. During the first second or two the machine behaves as though it had 25% droop rather than 4%, so its contribution to \(1/R\) is roughly a fifth of the nameplate value:

\[ \frac{1}{R_{\text{effective}}}\Big|_{t<2\ \text{s}} \approx \frac{1}{5}\cdot\frac{1}{R_p} \]

Which qualifies the folklore. Hydro plant is often called the fastest primary responder, and on a 10–30 s view it is: no thermal storage, no reheater, and full output available. On a 1–3 s view — the view the nadir takes — it is handicapped by its own water column and by the detuning that column forces on its governor.

The comparison worth carrying away. Over the first two seconds, a non-reheat steam unit delivers essentially all of its response, a hydro unit delivers a fraction of its response and starts by delivering the wrong sign, and a reheat unit delivers about 43% (Problem 16). Over sixty seconds all three deliver their full \(1/R\). Every one of them appears identically in \(\beta\).

The right-half-plane zero is the only element in this chapter that a controller cannot design around, only design for. Poles can be moved by feedback; a zero in the right half-plane cannot, and it puts a hard ceiling on the achievable bandwidth of any loop containing it. Hydro governing is the standard undergraduate example precisely because the physics behind the zero — water has mass — is so easy to state.
Answer\(\Delta P_m(t) = 1 - 3e^{-2t/T_w}\): it starts at −2 pu, crosses zero at 0.55 s and reaches 90% at 1.70 s. Transient droop of roughly \(5R_p\) is then required, cutting the unit's effective primary response to about a fifth in the first two seconds
Problem 18Exam levelTime Error

An interconnection runs at an average of 49.98 Hz for six hours because one area's AGC is out of service and its ACE is being held near zero by tie-flow alone.

  1. Find the accumulated synchronous-clock time error.
  2. Find how long the system must be run at 50.02 Hz to pay it back.
  3. If the tie flow was also 30 MW above schedule throughout, find the inadvertent interchange in MWh, and explain why \(\Delta f = 0\) alone does not prevent either error.
Solution

A synchronous clock counts cycles, so its error is the integral of the fractional frequency deviation:

\[ \Delta t = \int_0^{T}\frac{\Delta f}{f_0}\,d\tau = \frac{\Delta f}{f_0}\,T \]
\[ \Delta t = \frac{-0.02}{50}\times(6\times3600) = -4\times10^{-4}\times21\,600 = -8.64\ \text{s} \]

Every mains-timed clock in the interconnection is 8.64 seconds slow. The deviation was only 0.04% and would pass unremarked on any frequency trace; it is the integration over six hours that makes it visible.

Paying it back requires the opposite deviation for the same accumulated area:

\[ T_{\text{correction}} = \frac{|\Delta t|\,f_0}{\Delta f_{\text{corr}}} = \frac{8.64\times50}{0.02} = 21\,600\ \text{s} = 6\ \text{h} \]

Symmetric, because the correction offset is equal and opposite. Real time-error corrections use a smaller offset — a few millihertz — over a longer window, so that the correction is invisible to everything except the clocks.

The inadvertent interchange, which accumulates in exactly the same way but in energy:

\[ E_{\text{inadv}} = \int_0^{T}\Delta P_{\text{tie}}\,d\tau = 30\ \text{MW}\times6\ \text{h} = 180\ \text{MWh} \]

Energy delivered but not contracted. It is repaid by running the tie 30 MW below schedule for six hours, or 10 MW below for eighteen — payback in kind rather than in cash, which is the standard arrangement between interconnected utilities.

Why zero ACE does not prevent either. The supplementary loop drives \(\mathrm{ACE} \to 0\), which forces \(\Delta f \to 0\) and \(\Delta P_{\text{tie}} \to 0\) — but says nothing about their integrals:

\[ \Delta f \to 0 \quad\text{does not imply}\quad \int\Delta f\,dt \to 0 \]

Whatever error accumulated while the loop was correcting stays on the books for ever, and there is nothing in \(\mathrm{ACE}_i = \Delta P_{\text{tie},i} + B_i\Delta f\) that can remove it.

What is actually done. A slow outer correction is added to the ACE:

\[ \mathrm{ACE}_i' = \Delta P_{\text{tie},i} + B_i\,\Delta f + \underbrace{k_t\!\int\!\Delta f\,dt}_{\text{time-error}} + \underbrace{k_e\!\int\!\Delta P_{\text{tie},i}\,dt}_{\text{inadvertent}} \]

Two further integrators, deliberately made very slow — hours, not minutes — so that they never interact with the frequency-control loop. They are switched in by agreement between the interconnected areas, usually when the accumulated error passes a stated threshold such as ten seconds.

Three integrals are in play in this chapter and they are removed by three different mechanisms. The integral of ACE is removed by AGC in minutes. The integral of frequency error — clock time — is removed by an agreed correction over hours. The integral of tie deviation — inadvertent energy — is removed by scheduled payback over days. Each is one order of integration further out than the last, and each therefore needs a loop one order slower.
AnswerClocks lose \(\mathbf{8.64}\) s, repaid by 6 h at 50.02 Hz; inadvertent interchange is 180 MWh. Zero ACE fixes \(\Delta f\) and \(\Delta P_{\text{tie}}\) but never their integrals, so separate slow corrections are required
Problem 19Exam levelRamp Rates

An area's AGC computes a new set of setpoints every 4 s. Its regulating units are 400 MW and 300 MW machines whose permitted loading rate is 3% of rating per minute. At a given instant the area control error is −120 MW.

  1. Find the time the units need to deliver the correction.
  2. How many AGC cycles pass during that time?
  3. Reconcile this with the integral gain \(K_I = 0.372\) found in Problem 8, and state what limits the loop in practice.
Solution

The plant's ramp capability:

\[ 0.03\times400 = 12.0\ \text{MW/min}, \qquad 0.03\times300 = 9.0\ \text{MW/min}, \qquad \text{total } 21.0\ \text{MW/min} \]
\[ t = \frac{120}{21.0} = 5.71\ \text{min} = 343\ \text{s} \]

Nearly six minutes to answer a disturbance that primary control had arrested inside two seconds. That gap between the two time scales is not a defect of AGC; it is the reason the two loops exist separately.

The number of cycles:

\[ \frac{343\ \text{s}}{4\ \text{s per cycle}} = 86\ \text{cycles} \]

The controller issues 86 successive setpoints, each of which the plant only partially achieves before the next arrives. The AGC cycle is far faster than the plant, which is deliberate — it means the loop is always working from a fresh measurement and never waiting.

The reconciliation. Problem 8's gain of \(K_I = 0.372\) gave \(\omega_n = 1.52\) rad/s, a settling time of a few seconds. That controller would demand the whole correction inside about two seconds:

\[ \frac{120\ \text{MW}}{2\ \text{s}} = 60\ \text{MW/s} = 3600\ \text{MW/min} \quad\text{against}\quad 21\ \text{MW/min available} \]

A demand 170 times what the plant can deliver. The valves would saturate, the integrator would wind up against the limit, and the loop's actual behaviour would bear no relation to the linear design.

So what really limits \(K_I\)? Three ceilings, and the linear stability calculation is the loosest of them:

\[ \begin{array}{lll} 1 & \text{Plant ramp rate} & 21\ \text{MW/min here — the binding constraint} \\ 2 & \text{Inter-area modes} & \omega_n \ll \omega_{\text{osc}};\ 0.79\ \text{Hz in Problem 9} \\ 3 & \text{Damping } \zeta & \text{the } K_I \text{ of Problem 8 — the loosest} \end{array} \]

Setting the gain from criterion 3 alone gives a controller that is stable on paper and useless in service, because it spends every disturbance in saturation.

Design practice. The gain is set so that the loop's demand is comfortably inside the fleet's ramp capability — typically the bulk of the restoration in ten to thirty seconds for a small disturbance, with the loop simply rate-limited for a large one. Anti-windup on the integrator is mandatory, because without it the ACE integral accumulated during a five-minute ramp would drive a large overshoot the moment the plant catches up.

A design lever the problem hints at. If six minutes is too slow, the answer is not more gain but more regulating units. Adding a third 300 MW machine on regulation raises the fleet ramp to 30 MW/min and cuts the time to 4.0 min. Regulating capability is a scheduled quantity, procured like reserve, and it is what actually determines an area's AGC performance.

Every controller in this chapter is limited by something outside its own transfer function. The governor is limited by the turbine's thermal stress and the steam chest; the supplementary loop by the fleet's ramp rate and by the inter-area mode; the whole apparatus by how much plant is on regulation at all. A control gain chosen from a damping ratio alone has ignored all three.
Answer21 MW/min of ramp capability clears the 120 MW ACE in 5.71 min, across 86 AGC cycles. The \(K_I\) of Problem 8 would demand 3600 MW/min — 170 times the plant's capability — so ramp rate, not damping, is the binding constraint
Problem 20ChallengeA Low-Inertia Area

An area carries 12 000 MW with \(H = 4\) s on a 15 000 MVA base and \(\beta = 1800\) MW/Hz. Its largest single infeed is 900 MW. Take the ramp model of Problem 15 with \(T_r = 6\) s; load shedding is armed at 48.8 Hz.

  1. Find the RoCoF, the nadir and the settled frequency for the loss of the largest infeed today.
  2. Half the synchronous plant is displaced by inverter-connected wind and solar, halving the effective inertia while \(\beta\) is maintained by the remaining plant. Repeat.
  3. Size a battery delivering its full power within the first cycle so that the nadir is held at 49.0 Hz. State what the battery does not fix.
Solution

Today's system. Three quantities, three independent formulas:

\[ K = \frac{f_0\,\Delta P_L}{2HS} = \frac{50\times900}{2\times4\times15\,000} = \frac{45\,000}{120\,000} = 0.375\ \text{Hz/s} \]
\[ \Delta f_{\text{nadir}} = -\frac{K T_r}{2} = -0.375\times3 = -1.125\ \text{Hz} \quad\Rightarrow\quad f_{\text{nadir}} = 48.875\ \text{Hz} \]
\[ \Delta f_{ss} = \frac{-900}{1800} = -0.500\ \text{Hz} \quad\Rightarrow\quad f = 49.500\ \text{Hz} \]

The nadir clears 48.8 Hz by 75 mHz. That is a thin margin, and it is the margin the system is actually designed to, not the comfortable 0.7 Hz that the settling point suggests.

After the displacement. \(\beta\) is maintained by assumption, so the destination does not move; only the path does:

\[ K' = \frac{50\times900}{2\times2\times15\,000} = 0.750\ \text{Hz/s}, \qquad \Delta f_{\text{nadir}} = -0.750\times3 = -2.25\ \text{Hz} \]
\[ f_{\text{nadir}} = 47.75\ \text{Hz}, \qquad f_{ss} = 49.500\ \text{Hz unchanged} \]

Load shedding operates, and at 47.75 Hz so do the under-frequency protections of the generators themselves on many machines. A system that maintained every steady-state indicator has lost its ability to survive its own largest infeed.

Size the battery. Fast frequency response delivered essentially instantaneously reduces the deficit that inertia and governors must cover, so the nadir formula applies to \(\Delta P_L - P_b\):

\[ |\Delta f_{\text{nadir}}| = \frac{f_0\,(\Delta P_L - P_b)}{2H'S}\cdot\frac{T_r}{2} = \frac{50\,(900-P_b)}{60\,000}\times3 = \frac{900-P_b}{400} \]
\[ \text{Require } |\Delta f_{\text{nadir}}| = 1.00\ \text{Hz}: \qquad 900 - P_b = 400 \quad\Rightarrow\quad P_b = 500\ \text{MW} \]

Check: with 500 MW injected, \(K = 50\times400/60\,000 = 0.333\) Hz/s and the nadir is \(0.333\times3 = 1.00\) Hz, giving 49.00 Hz. ✓

The energy, which is the part that surprises people. The battery holds 500 MW for roughly the six seconds until the governors have taken over, then ramps out:

\[ E \approx 500\ \text{MW}\times\frac{6}{3600}\ \text{h} \approx 0.83\ \text{MWh}\ \text{of through-energy, say 5–10 MWh installed} \]

A very high power rating on a very small store. Frequency-response batteries are sized by megawatts and by C-rate, not by megawatt-hours, which is why they are physically much smaller than an energy-arbitrage installation of the same power.

What the battery does not fix. Three things, and each needs a different remedy:

\[ \begin{array}{lll} \text{RoCoF} & \text{still } 0.75\ \text{Hz/s at } t=0^+ & \text{needs inertia, or synchronous condensers} \\ \Delta f_{ss} & \text{still } -0.50\ \text{Hz} & \text{needs } \beta,\ \text{i.e. governors} \\ \text{Sustained energy} & \text{battery empties in minutes} & \text{needs plant on AGC} \end{array} \]

The first is the sharpest limitation. A battery responding within a cycle still responds after the frequency has begun to fall, so it cannot flatten the initial slope, and RoCoF relays on embedded generation see the full 0.75 Hz/s regardless. Only rotating mass — retained synchronous plant, synchronous condensers, or grid-forming converters emulating inertia — changes that number.

The three-lever summary for a low-inertia area, in the order the disturbance meets them:

\[ \begin{array}{llll} \text{Lever} & \text{acts on} & \text{time} & \text{quantity fixed} \\ \hline \text{Inertia / synchronous condensers} & \text{RoCoF} & 0\text{--}1\ \text{s} & 2HS \\ \text{Fast frequency response (battery)} & \text{nadir} & 0.1\text{--}10\ \text{s} & \text{effective } \Delta P_L \\ \text{Governor droop} & \text{settling point} & 2\text{--}30\ \text{s} & \beta \\ \text{AGC} & \text{return to nominal} & 30\ \text{s}\text{--}15\ \text{min} & \int\!\mathrm{ACE} \\ \end{array} \]
An area can hold every steady-state indicator constant and still lose the ability to survive its largest single event. \(\beta\) was maintained here and \(\Delta f_{ss}\) never moved; what changed was a quantity that appears in no steady-state equation in this set. That is why RoCoF and nadir have joined \(\beta\) as contracted, monitored, separately procured products, and why the frequency-control problem of the 2020s is a transient one rather than a steady-state one.
AnswerToday: 0.375 Hz/s, nadir 48.875 Hz, settling 49.500 Hz. At half inertia: 0.750 Hz/s, nadir 47.75 Hz — load shedding operates with the settling point unchanged. A 500 MW battery restores a 49.00 Hz nadir but changes neither the RoCoF nor the settling point
Practice

Practice Problems

Take \(f_0 = 50\) Hz throughout, and convert every droop to MW/Hz before doing anything else. Work each on paper before opening the answer; the answer is given so you can check yourself, the method deliberately is not.

  1. P1. A 400 MW unit with 4% droop is carrying 220 MW when the frequency falls to 49.60 Hz. Find its new output, the speed changer being untouched.

    Show answer
    \(1/R = 200\) MW/Hz, so \(\Delta P = 200\times0.4 = 80\) MW and the unit runs at 300 MW.
  2. P2. An area with \(\beta = 2400\) MW/Hz loses 600 MW of generation. Find the settled frequency.

    Show answer
    \(\Delta f = -600/2400 = -0.25\) Hz, so \(f = \mathbf{49.75}\) Hz.
  3. P3. Two units, 200 MW at 5% and 400 MW at 4%, share an isolated system with damping neglected. A 56 MW load increase occurs. Find \(\Delta f\) and each unit's increment.

    Show answer
    \(1/R = 80\) and \(200\) MW/Hz, sum 280. \(\Delta f = \mathbf{-0.20}\) Hz; increments 16 MW and 40 MW.
  4. P4. An area supplies 3000 MW. A 1% change in frequency produces a 2% change in load. Find \(D\).

    Show answer
    \(D = (0.02\times3000)/(0.01\times50) = 60/0.5 = \mathbf{120}\) MW/Hz.
  5. P5. An area has \(H = 6\) s on a 5000 MVA base. A 250 MW infeed trips. Find the initial RoCoF.

    Show answer
    \(-f_0\Delta P_L/(2HS) = -50(250)/(2\times6\times5000) = \mathbf{-0.208}\) Hz/s.
  6. P6. For an area with \(D = 0.012\) pu MW/Hz and \(H = 5\) s, find \(K_p\) and \(T_p\).

    Show answer
    \(K_p = 1/D = \mathbf{83.3}\) Hz/pu MW; \(T_p = 2H/(f_0D) = 10/0.6 = \mathbf{16.7}\) s.
  7. P7. A tie line has \(X_{12} = 0.8\) pu, \(V_1 = V_2 = 1.0\) pu and an angle difference of 25°. Find \(T_{12}\).

    Show answer
    \(T_{12} = \cos25^\circ/0.8 = 0.90631/0.8 = \mathbf{1.133}\) pu MW/rad.
  8. P8. Two areas have \(\beta_1 = 1000\) and \(\beta_2 = 2000\) MW/Hz. A 180 MW load increase occurs in area 2 with supplementary control out of service. Find \(\Delta f\), the tie flow change with sign, and each area's contribution.

    Show answer
    \(\Delta f = -180/3000 = \mathbf{-0.06}\) Hz. \(\Delta P_{12} = -\beta_1\Delta f = \mathbf{+60}\) MW, so area 1 exports 60 MW; area 2 supplies 120 MW.
  9. P9. For P8 with ideal biases, find both area control errors.

    Show answer
    \(\mathrm{ACE}_1 = 60 + 1000(-0.06) = \mathbf{0}\); \(\mathrm{ACE}_2 = -60 + 2000(-0.06) = \mathbf{-180}\) MW.
  10. P10. Repeat P9's \(\mathrm{ACE}_1\) with \(B_1\) mistakenly set to 500 MW/Hz, and say what area 1 does.

    Show answer
    \(60 + 500(-0.06) = \mathbf{+30}\) MW. Positive is read as surplus, so area 1 wrongly reduces generation while area 2 is short.
  11. P11. Three units have \(d^2C/dP^2 = 0.010,\ 0.016,\ 0.020\) Rs/MW²h. Find the participation factors and the split of an 85 MW AGC correction.

    Show answer
    Reciprocals 100, 62.5, 50, sum 212.5. \(\alpha = 0.4706,\ 0.2941,\ 0.2353\); split 40, 25 and 20 MW.
  12. P12. An interconnection runs at 49.99 Hz for 10 hours. Find the synchronous-clock error.

    Show answer
    \((-0.01/50)\times36\,000 = \mathbf{-7.2}\) s — every mains clock is 7.2 seconds slow.
Challenge

Challenge Problems

Three problems that need an idea rather than a formula — the awkward corners where a working assumption of this chapter quietly fails.

  1. C1 — Non-interactive control with many areas. An interconnection has \(n\) areas, each with its own \(\beta_i\), all mutually interconnected, and each setting \(B_i = \beta_i\). A disturbance \(\Delta P_{Lk}\) occurs in area \(k\). Prove that after primary control settles, \(\mathrm{ACE}_k = -\Delta P_{Lk}\) and \(\mathrm{ACE}_j = 0\) for every \(j \ne k\) — and identify the one assumption the proof needs that the two-area case hid.

    Show answer

    The one frequency. After the tie flows have stopped changing, every tie is an integrator with a settled output, so every pair of areas is at the same frequency: \(\Delta f_1 = \cdots = \Delta f_n = \Delta f\). This is the whole content of the interconnection being synchronous, and it is what makes an \(n\)-area problem no harder than a two-area one.

    Power balance of an undisturbed area \(j\). Its governors supply \(-\Delta f/R_j\) and its load relieves \(-D_j\Delta f\); with no local disturbance, all of that must leave on its ties:

    \[ \Delta P_{\text{tie},j} = -\frac{\Delta f}{R_j} - D_j\,\Delta f = -\beta_j\,\Delta f \]

    Note that this holds however many ties area \(j\) has and whatever their reactances: \(\Delta P_{\text{tie},j}\) is the net interchange deviation, and the internal distribution among its ties is irrelevant.

    Its ACE is then immediate:

    \[ \mathrm{ACE}_j = \Delta P_{\text{tie},j} + B_j\Delta f = -\beta_j\Delta f + \beta_j\Delta f = 0 \]

    The disturbed area follows from the fact that the ACEs must sum to the disturbance. Summing the balances of all \(n\) areas, every tie appears once with each sign and cancels:

    \[ \sum_i \beta_i\,\Delta f = -\Delta P_{Lk} \;\Longrightarrow\; \Delta f = \frac{-\Delta P_{Lk}}{\sum_i \beta_i} \]
    \[ \Delta P_{\text{tie},k} = -\beta_k\Delta f - \Delta P_{Lk} \;\Longrightarrow\; \mathrm{ACE}_k = -\beta_k\Delta f - \Delta P_{Lk} + \beta_k\Delta f = -\Delta P_{Lk} \]

    The hidden assumption is that the whole interconnection reaches one steady frequency — that is, that no tie is at an angle limit and none has tripped. If a tie saturates or opens, the interconnection splits into islands, each with its own frequency, and areas in the healthy island no longer see the disturbed area's \(\Delta f\) at all. Their ACEs are then non-zero for a reason no bias setting can correct, because the premise \(\Delta f_j = \Delta f_k\) has failed.

    The second assumption, invisible in the two-area case, is that each area's \(\Delta P_{\text{tie},i}\) is metered as a net sum over all its ties and with a consistent sign convention. An area with six interconnectors that mis-signs one of them has an ACE that is wrong by twice that tie's flow, permanently, and the error looks exactly like a mis-set bias — which is why tie metering is checked before bias settings are questioned.

  2. C2 — The offset that would not go away. An operator reports that after every large disturbance the frequency returns cleanly to 50.000 Hz, but the tie-line interchange settles about 40 MW above schedule and stays there for hours until manually corrected. Both areas' AGC is in service and both integrators are known to be working. Identify the possible causes, say which is most likely, and state how each would be confirmed.

    Show answer

    First, what the symptom rules out. Problem 11 showed that the only steady state with two working ACE integrators is \(\Delta f = 0\) and \(\Delta P_{12} = 0\), for any positive biases. A persistent tie offset with the frequency exactly at nominal is therefore not an equilibrium of the model at all. Something in the real loop differs from the model, and a mis-set bias is not a candidate — bias affects the path, never the destination.

    Cause 1: the tie-line term is missing from one ACE. If area 2's controller is running on frequency alone — \(\mathrm{ACE}_2 = B_2\Delta f\), the mode called flat frequency control — then its integrator is satisfied by \(\Delta f = 0\) whatever the tie is doing. Area 1's is then satisfied too, since its own ACE can be driven to zero by any combination. Every tie flow becomes an equilibrium. Confirm by reading the control mode of each area's AGC; flat frequency control is a legitimate mode, used by a small area inside a large interconnection, and it is frequently left selected by accident after a test.

    Cause 2: a metering error on the tie. If the interchange telemetry reads 40 MW low, the controller drives the measured deviation to zero and leaves a real 40 MW offset. The frequency is unaffected, because frequency is measured independently. Confirm by comparing the two areas' independent measurements of the same tie: they should sum to zero and will differ by 40 MW.

    Cause 3: a wrong schedule. The two areas' scheduled interchange values disagree by 40 MW — a data-entry error at the start of a delivery period, or a schedule change accepted by one party and not the other. Each controller is holding its own idea of schedule perfectly. Confirm by comparing the schedule files; this is why interchange schedules are confirmed electronically between control centres before each hour.

    Cause 4: integrator saturation or a deadband. If the ACE deadband is set at, say, ±50 MW to reduce plant movement, a 40 MW offset generates no control action at all and persists indefinitely. Confirm by reading the deadband setting; note that this cause predicts an offset that is always smaller than the deadband, which is a testable signature.

    Which is most likely. Cause 3, then cause 2. The symptom — appears after disturbances, persists for hours, cleared by manual action — is characteristic of a bookkeeping disagreement rather than a control fault, and both areas' controllers behaving perfectly is exactly what a schedule mismatch produces. The diagnostic that separates them in one step is to ask the neighbour what its ACE reads. If area 2 reports zero ACE at the same moment as area 1 reports zero ACE while the tie is 40 MW off, the two are working to different numbers, and the fault is in the schedule or the metering, not in the loops.

    The general lesson: a control loop with an integrator cannot leave a steady-state error in the quantity it actually measures. If an error persists, the loop is measuring something other than what you think it is.

  3. C3 — Governor deadband. Every governor in an area has a deadband of ±20 mHz: the valve does not move until \(|\Delta f|\) exceeds 0.02 Hz, after which it responds normally to the excess. The area has \(1/R = 1500\) MW/Hz and \(D = 30\) MW/Hz. Find the settled frequency for a 30 MW disturbance and for a 300 MW disturbance, compare with the deadband-free values, and explain the pattern.

    Show answer

    Model the deadband. Once the frequency has moved beyond the band, each governor responds to the excess deviation only, so the steady-state balance becomes

    \[ \frac{1}{R}\big(|\Delta f| - f_{db}\big) + D\,|\Delta f| = \Delta P_L \qquad (|\Delta f| > f_{db}) \]
    \[ |\Delta f| = \frac{\Delta P_L + f_{db}/R}{D + 1/R} = \frac{\Delta P_L + 1500(0.02)}{1530} = \frac{\Delta P_L + 30}{1530} \]

    The deadband adds a fixed 30 MW to the numerator — a constant penalty, independent of the size of the disturbance.

    The small disturbance. Without deadband, \(30/1530 = 0.0196\) Hz — which is inside the band, so on the ideal model the governors would never have moved at all. With the band:

    \[ |\Delta f| = \frac{30+30}{1530} = 0.0392\ \text{Hz}, \qquad \beta_{\text{eff}} = \frac{30}{0.0392} = 765\ \text{MW/Hz} \]

    Exactly double the deviation, exactly half the effective stiffness. The frequency has to travel through the whole band before any help arrives, and for a disturbance this size the band is as large as the response.

    The large disturbance.

    \[ |\Delta f| = \frac{300+30}{1530} = 0.2157\ \text{Hz} \quad\text{against}\quad 0.1961\ \text{Hz ideal}, \qquad \beta_{\text{eff}} = 1391\ \text{MW/Hz} \]

    Only 10% worse, and \(\beta_{\text{eff}}\) is 91% of nominal. The fixed 30 MW penalty is now a tenth of the disturbance instead of all of it.

    The pattern, and it is the point. Deadband degrades the response by a fixed amount in megawatts, so its proportional effect is enormous for small disturbances and negligible for large ones:

    \[ \begin{array}{lccc} \Delta P_L\ (\text{MW}) & |\Delta f|\ \text{ideal} & |\Delta f|\ \text{with band} & \beta_{\text{eff}}/\beta \\ \hline 30 & 0.0196 & 0.0392 & 0.50 \\ 300 & 0.1961 & 0.2157 & 0.91 \\ \end{array} \]

    Why this matters more than it looks. The frequency of an interconnection spends nearly all of its life responding to small, continuous load fluctuations, not to unit trips. Deadband therefore dominates ordinary operation — it is the main reason the frequency of a real system wanders within a band of a few tens of millihertz rather than sitting still — while being nearly invisible in the large-event studies from which \(\beta\) is measured. Problem 4's measured \(\beta\) came from a 250 MW event and is therefore close to the true value; the same area's response to a 20 MW fluctuation would be half as stiff.

    Why the deadband exists at all: without it, every governor valve on the system would be in continuous small-amplitude motion, wearing servo linkages and cycling turbine stress. It is a deliberate trade of frequency quality for plant life, and grid codes set a maximum — typically ±10 to ±15 mHz — rather than requiring zero.

Self-Test

Multiple-Choice Questions

  1. MCQ 1. A 500 MW unit with 4% droop on a 50 Hz system has \(1/R\) equal to:
    (a) 4 MW/Hz   (b) 125 MW/Hz   (c) 250 MW/Hz   (d) 500 MW/Hz

    Show answer
    (c)\(500/(0.04\times50) = 250\). Option (b) is what you get by dividing the rating by 4 instead of by \(0.04f_0\), which is the single commonest slip in this chapter. Problem 1.
  2. MCQ 2. The steady-state frequency deviation of an isolated area after a load step \(\Delta P_L\) is:
    (a) \(-\Delta P_L/D\)   (b) \(-R\,\Delta P_L\)   (c) \(-\Delta P_L/(D+1/R)\)   (d) \(-\Delta P_L/(D+R)\)

    Show answer
    (c). (a) is the case with every governor saturated — a real limit, not the general answer (Problem 5). (d) adds quantities with different units; \(D\) is MW/Hz and \(R\) is Hz/MW, and only \(1/R\) can be added to \(D\).
  3. MCQ 3. The initial rate of change of frequency after a generator trip depends on:
    (a) the area stiffness \(\beta\)   (b) the droop setting   (c) the inertia and the size of the trip   (d) the turbine time constant

    Show answer
    (c). At \(t = 0^+\) no valve has moved and \(\Delta f\) is still zero, so neither the droop term nor the damping term contributes anything. Problems 6 and 7.
  4. MCQ 4. The supplementary loop must be integral rather than proportional because:
    (a) it is faster   (b) only infinite d.c. gain gives zero steady-state error   (c) it is more stable   (d) proportional action would fight the governors

    Show answer
    (b). A proportional supplementary term \(-K_c\Delta f\) merely adds \(K_c\) to \(\beta\), giving a smaller offset but never zero. (a) is backwards — the integral loop is deliberately the slow one. Problem 8.
  5. MCQ 5. With \(B_i = \beta_i\) in both areas and a disturbance in area 1, \(\mathrm{ACE}_2\) equals:
    (a) \(-\Delta P_{L1}\)   (b) zero   (c) \(-\beta_2\Delta P_{L1}/(\beta_1+\beta_2)\)   (d) \(\beta_2\Delta f\)

    Show answer
    (b) — exactly zero, which is the design intent. (c) and (d) are each one of the two terms in area 2's ACE; the point is that they cancel. Problem 11.
  6. MCQ 6. The tie-line block relating \(\Delta P_{12}\) to the frequency deviations is:
    (a) \(2\pi T_{12}\,s\)   (b) \(2\pi T_{12}/s\)   (c) \(T_{12}/(2\pi s)\)   (d) \(T_{12}\)

    Show answer
    (b). Power depends on angle, and angle is the integral of frequency, so the tie is an integrator. (a) has it as a differentiator, which would make the flow settle instantly whatever the frequencies were. Problem 9.
  7. MCQ 7. Two areas on primary control alone, disturbance \(\Delta P_{L1}\) in area 1. The tie-flow change is:
    (a) \(-\beta_1\Delta P_{L1}/(\beta_1+\beta_2)\)   (b) \(-\beta_2\Delta P_{L1}/(\beta_1+\beta_2)\)   (c) \(-\Delta P_{L1}/2\)   (d) zero

    Show answer
    (b). The undisturbed area's stiffness appears in the numerator, because everything area 2 does has to travel down the tie. Swapping the subscripts, as in (a), is the standard error. Problem 10.
  8. MCQ 8. The AGC participation factor of a unit is proportional to:
    (a) its rating   (b) \(1/R_i\)   (c) the reciprocal of \(d^2C_i/dP_i^2\)   (d) its incremental cost

    Show answer
    (c). (b) is the droop allocation, which governs the first seconds and is a different split entirely — that difference is the whole content of Problem 13. (d) is what the dispatch equalises, not what it apportions by.
  9. MCQ 9. If every unit in an area is set to the same percentage droop, a load increment is shared in proportion to:
    (a) inertia   (b) rating   (c) incremental cost   (d) equally in MW

    Show answer
    (b), because \(1/R_i = P_{\text{rated},i}/(R\%\,f_0/100)\) is proportional to rating when \(R\%\) is common. Every machine then picks up the same fraction of its own capacity. Problems 2 and 3.
  10. MCQ 10. An area's inertia halves while \(\beta\) is held constant. Which quantity is unchanged after a generator trip?
    (a) the RoCoF   (b) the frequency nadir   (c) the settled frequency deviation   (d) all three change

    Show answer
    (c). \(\Delta f_{ss} = -\Delta P_L/\beta\) contains no \(H\) whatever. The RoCoF doubles and the nadir roughly doubles — which is exactly why a system can pass every steady-state check and still shed load. Problems 15 and 20.
  11. MCQ 11. Area 1's bias is set well below \(\beta_1\). A disturbance occurs in area 2. Area 1's controller:
    (a) does nothing   (b) reads a positive ACE and reduces its generation   (c) reads a negative ACE and raises it   (d) becomes unstable

    Show answer
    (b). The tie term (positive, since area 1 is exporting help) outweighs the under-weighted frequency term, and a positive ACE is read as over-generation. (d) is wrong and worth being clear about: the equilibrium is unchanged for any positive bias, so the error costs regulating energy, not stability. Problem 12.
  12. MCQ 12. With both supplementary loops in service and any positive bias settings, the final steady state has:
    (a) \(\Delta f = 0\) only   (b) \(\Delta P_{12} = 0\) only   (c) both zero   (d) values that depend on \(B_1\) and \(B_2\)

    Show answer
    (c). Setting both ACEs to zero and adding gives \((B_1+B_2)\Delta f = 0\), whence both. (d) is the intuitive answer and is wrong — bias decides the path, never the destination. Problem 11.
Reference

Key Formulas

StatementRelationNotes
Swing equation in frequency\(\dfrac{2H}{f_0}\dfrac{d(\Delta f)}{dt} = \Delta P_m - \Delta P_L - D\,\Delta f\)Per unit on the area base
Rate of change of frequency\(\mathrm{RoCoF} = -\dfrac{f_0\,\Delta P_L}{2H}\)At \(t=0^+\); independent of \(R\) and \(D\)
Generator–load block\(G_p(s) = \dfrac{K_p}{1+sT_p}\), \(K_p = 1/D\), \(T_p = 2H/(f_0D)\)\(T_p \approx 10\text{–}25\) s
Droop, percent to MW/Hz\(\dfrac{1}{R} = \dfrac{P_{\text{rated}}}{(R\%/100)\,f_0}\)Do this first, always
Governor characteristic\(\Delta P_g = \Delta P_{\text{ref}} - \Delta f/R\)Speed changer shifts the line
Composite regulation\(1/R_{\text{eq}} = \sum_i 1/R_i\)Responsive units only
Area frequency response\(\beta = D + 1/R_{\text{eq}}\) MW/HzMeasured from events, not nameplates
Primary steady state\(\Delta f_{ss} = -\Delta P_L/\beta\)Every lag is unity at \(s = 0\)
Primary loop time constant\(T_p' = \dfrac{T_p}{1+K_p/R}\)Typically 30–60 times shorter
Nadir, ramp model\(\Delta f_{\text{nadir}} = -\dfrac{\mathrm{RoCoF}\times T_r}{2}\)Half, because the deficit is triangular
Integral supplementary control\(\Delta P_{\text{ref}} = -K_I\!\int\!\Delta f\,dt\)Sign matters; \(K_I>0\)
Second-order parameters\(\omega_n = \sqrt{K_pK_I/T_p}\), \(2\zeta\omega_n = (1+K_p/R)/T_p\)Only \(\omega_n\) contains \(K_I\)
Turbine, non-reheat\(1/(1+sT_T)\), \(T_T = 0.2\text{–}0.5\) sGovernor lag \(T_g = 0.08\text{–}0.3\) s
Turbine, reheat\(\dfrac{1+sK_rT_r}{(1+sT_T)(1+sT_r)}\)Prompt fraction \(\approx K_r\); \(T_r = 5\text{–}10\) s
Turbine, hydro\(\dfrac{1-sT_w}{1+0.5sT_w}\)Starts at \(-2\) pu; needs transient droop
Synchronizing coefficient\(T_{12} = \dfrac{|V_1||V_2|}{X_{12}}\cos(\delta_1^0-\delta_2^0)\)MW per radian
Tie-line block\(\Delta P_{12}(s) = \dfrac{2\pi T_{12}}{s}\big[\Delta F_1 - \Delta F_2\big]\)An integrator, not a gain
Inter-area mode\(\omega_{\text{osc}} = \sqrt{2\pi f_0 T_{12}\left(\dfrac{1}{2H_1}+\dfrac{1}{2H_2}\right)}\)0.1–0.8 Hz; caps every gain
Two areas, primary only\(\Delta f = \dfrac{-\Delta P_{L1}}{\beta_1+\beta_2}\), \(\Delta P_{12} = \beta_2\Delta f\)Undisturbed area's \(\beta\) in the tie term
Area control error\(\mathrm{ACE}_i = \Delta P_{\text{tie},i} + B_i\,\Delta f\)Export positive; \(B_i\) in MW/Hz
Non-interactive bias\(B_i = \beta_i \Rightarrow \mathrm{ACE}_k = -\Delta P_{Lk},\ \mathrm{ACE}_{j\ne k} = 0\)Final state independent of \(B\)
Participation factors\(\alpha_i = \dfrac{1/(d^2C_i/dP_i^2)}{\sum_j 1/(d^2C_j/dP_j^2)}\)From Set 27; \(\sum\alpha_i = 1\)
Time error\(\Delta t = \dfrac{\Delta f}{f_0}\,T\)Zero ACE does not remove it
Diagnostics

Common Mistakes

  1. Using a droop in percent as though it were MW/Hz. The conversion \(1/R = P_{\text{rated}}/[(R\%/100)f_0]\) must be done before anything else, and it involves the machine's rating, which two units with the same percentage do not share — Problems 1 and 3.

  2. Writing \(\beta = D + R\). \(D\) is MW/Hz and \(R\) is Hz/MW; only \(1/R\) can be added to \(D\). A dimensional check catches it instantly — Problem 2.

  3. Including saturated or non-responsive units in \(1/R_{\text{eq}}\). Capacity that is spinning but blocked from moving contributes nothing to \(\beta\), and an area can be fully reserved and half as stiff as its paperwork claims — Problems 4 and 5.

  4. Expecting inertia to appear in the settled frequency. \(\Delta f_{ss} = -\Delta P_L/\beta\) contains no \(H\), and \(\mathrm{RoCoF} = -f_0\Delta P_L/2H\) contains no \(\beta\). Inertia sets the speed of the excursion, stiffness its destination — Problems 6 and 15.

  5. Quoting the settling point when the question is about protection. Load-shedding and generator under-frequency relays respond to the nadir, which can be twice the settled offset and is not predicted by \(\beta\) at all — Problems 15 and 20.

  6. Getting the sign of the integral loop backwards. A load increase makes \(\Delta f\) negative; \(\Delta P_{\text{ref}} = -K_I\!\int\!\Delta f\,dt\) with \(K_I > 0\) is therefore positive and raises generation. Reverse it and the loop drives the frequency away without limit — Problem 8.

  7. Putting the disturbed area's \(\beta\) in the tie-flow formula. \(\Delta P_{12} = -\beta_2\Delta P_{L1}/(\beta_1+\beta_2)\) uses the undisturbed area's stiffness, because everything that area contributes must travel down the tie — Problem 10.

  8. Believing that the bias setting changes the final state. Any positive pair of biases ends at \(\Delta f = 0\), \(\Delta P_{12} = 0\); the bias decides how much unnecessary plant movement happens on the way — Problems 11 and 12.

  9. Setting the frequency bias from declared governor settings. The bias should equal the measured \(\beta\). The area of Problem 4 would have set a bias more than twice its true response — Problems 4 and 12.

  10. Allocating the AGC correction by droop. Once \(\Delta f\) is back at zero every droop term has vanished; the final split is by participation factor, and units genuinely move in opposite directions during the handover — Problem 13.

  11. Choosing \(K_I\) from a damping ratio alone. The plant's ramp rate and the inter-area mode both bind long before the linear stability limit does, usually by two orders of magnitude — Problems 9 and 19.

  12. Assuming zero ACE clears the accumulated errors too. Driving \(\Delta f\) and \(\Delta P_{\text{tie}}\) to zero says nothing about their integrals, so clock time error and inadvertent interchange need separate, much slower loops — Problem 18.

  13. Treating all governor response as equally prompt. A reheat unit delivers about \(K_r\) of its response in the first second and a hydro unit briefly delivers the wrong sign, yet both appear in \(\beta\) at full value — Problems 16 and 17.

Looking Ahead

Frequency is one number for a whole interconnection, so one controller per area suffices and the entire chapter reduces to three algebraic lines plus a discussion of time scales. Voltage is not like that at all. Reactive power does not travel far, so the voltage at a bus is set mostly by the reactive balance in its own neighbourhood, and there are as many voltage problems as there are buses.

Chapter 34 takes up that second control task: excitation systems and the automatic voltage regulator, reactive compensation by capacitors, reactors and static var compensators, tap-changing transformers, and the reactive-power capability limits of a synchronous machine. The two loops are almost perfectly decoupled — active power with frequency here, reactive power with voltage there — which is what allowed this chapter to hold every voltage magnitude constant while linearising the tie-line power in Problem 9.