Q5Power System Engineering

Equal Area Criterion for a Fault Away from Line Ends

The equal area criterion is a graphical/energy-based method for assessing the transient (first-swing) stability of a single machine connected to an infinite bus (or an equivalent two-machine system reduced to this form) following a large disturbance, without needing to numerically integrate the full swing equation - it exploits the fact that the swing equation, multiplied through by d(delta)/dt and integrated, yields a statement in terms of areas under the power-angle curve, and the system remains stable (rotor angle delta oscillates but does not increase without bound) if and only if the 'accelerating' area (where mechanical power exceeds electrical power, causing the rotor to speed up) can be exactly balanced by an equal 'decelerating' area (where electrical power exceeds mechanical power, causing the rotor to slow back down) before delta reaches its maximum permissible value.

A 'fault away from line ends' refers to a fault location partway along a transmission line (rather than directly at the generator terminal bus or directly at the receiving-end/infinite bus), such that the network retains some nonzero transfer reactance (and hence nonzero power transfer capability, Pmax2 greater than zero) even during the fault period - this is in contrast to a fault directly at the sending-end bus of a single (non-parallel) line, which typically reduces the during-fault transfer power to exactly zero, since there is often no alternate path around such a fault location.

Three distinct power-angle curves are relevant: the pre-fault curve P1(delta)=Pmax1sin(delta) (using the pre-fault network's transfer reactance), the during-fault curve P2(delta)=Pmax2sin(delta) (using the reduced transfer reactance of the faulted network, computed via a star-delta network reduction technique that eliminates the fault point as an internal node, given that Pmax2 is nonzero for a fault away from line ends), and the post-fault curve P3(delta)=Pmax3*sin(delta) (using the network reactance after the faulted element is isolated by protective relaying and circuit breaker operation, generally different from - often intermediate between or in special cases equal to - the pre-fault reactance, depending on which specific line or element was faulted and subsequently removed).

Before the fault, the system operates at the initial angle delta0, where the mechanical input power Pm intersects the pre-fault curve: Pm=Pmax1*sin(delta0). When the fault occurs, the electrical power output suddenly drops to the (lower) during-fault curve P2(delta), so Pm now exceeds Pe (since P2(delta0) is less than Pm at the same angle delta0), causing the rotor to accelerate and delta to increase - this creates the accelerating area A1, the excess mechanical energy input integrated over the angle range from delta0 to the actual clearing angle deltac (whenever the fault is finally cleared by protective relaying and breaker operation).

Once the fault is cleared at angle deltac, the system transitions to the post-fault curve P3(delta), and if Pm is now less than P3(delta) at the clearing angle (which is generally the case immediately after clearing, since delta has increased beyond its equilibrium value), the rotor decelerates, and delta continues rising further but at a decreasing rate, tracing out the decelerating area A2 until delta reaches its maximum swing value delta_max, at which point d(delta)/dt returns to zero (momentarily) before delta swings back.

The system is stable for this particular fault-clearing scenario if the available decelerating area A2 (measured out to the maximum possible delta_max, the second intersection of the post-fault curve with Pm, or 180 degrees minus that intersection angle) is at least as large as the accelerating area A1 that has already accumulated by the actual clearing angle - if A2(available, up to delta_max) is less than A1, the system is unstable (delta increases without bound, and synchronism is lost) for that fault-clearing time.

Critical Clearing Angle and Critical Clearing Time

The critical clearing angle delta_cr is defined as the specific clearing angle at which the accelerating area A1 (from delta0 to delta_cr) exactly equals the maximum possible decelerating area A2 (from delta_cr to delta_max) - this represents the boundary (limiting) case: clearing the fault at any angle less than delta_cr (i.e., faster) guarantees stability with some margin, while clearing later than delta_cr results in instability. Setting A1=A2 (equal area condition) at the critical clearing angle:

Rearranging this equation to isolate cos(delta_cr) gives the standard closed-form expression for the critical clearing angle:

where delta0=arcsin(Pm/Pmax1) is the pre-fault operating angle and delta_max=pi-arcsin(Pm/Pmax3) is the maximum angle on the post-fault curve at which Pm is again matched (the second intersection point). Once delta_cr is known, the corresponding critical clearing time tcr (the maximum time the fault can be allowed to persist while still maintaining stability) is found by solving the swing equation during the fault period (using the during-fault power Pmax2) for the time at which delta first reaches delta_cr, starting from initial conditions delta(0)=delta0 and d(delta)/dt(0)=0:

For the special, simplified case where Pmax2=0 (fault directly at the sending end with no alternate path, giving zero power transfer during the fault - the classical, simpler equal-area problem, distinct from the 'fault away from line ends' case with Pmax2 nonzero examined in this question), this swing equation during the fault reduces to a constant-acceleration problem (Pm approximately constant, acting alone) with a simple closed-form solution: tcr=sqrt(4H(delta_cr-delta0)/(omega_sPm)). For the more general fault-away-from-line-ends case examined here (Pmax2 nonzero), the during-fault swing equation is nonlinear (since Pmax2*sin(delta) remains present) and generally must be solved by numerical step-by-step time-domain integration (such as the point-by-point or modified Euler/Runge-Kutta method) rather than by this simpler closed-form formula, to determine the actual critical clearing time corresponding to the critical clearing angle found from the equal-area criterion above.

In summary, the equal area criterion's key practical value lies precisely in its ability to determine the critical clearing angle and time without requiring full numerical time-domain simulation of the swing equation for every candidate clearing time, making it an efficient first-pass tool for protection engineers to determine the maximum permissible fault-clearing time for a given fault location and system operating condition, directly informing the required speed of protective relaying and circuit breaker operation for that specific contingency.

This combination of graphical/energy-based reasoning and the underlying differential-equation formulation illustrates why the equal area criterion has remained a standard teaching and quick-assessment tool in power system stability analysis for many decades, even as full numerical time-domain simulation has become the routine industry-standard method for detailed multi-machine stability studies.