🎓 Lesson 14
D5
Q-V Curve Construction & Critical Point Identification
The Q-V curve shows how reactive power (Q) affects voltage (V) at a bus, helping engineers spot when the power system is about to lose voltage stability.
🎯 Learning Objectives
- ✓ Calculate Q-V curve points using Newton-Raphson load flow iterations for a specified test bus
- ✓ Identify the critical (nose) point on a Q-V curve and interpret its voltage and reactive power coordinates
- ✓ Analyze system vulnerability by comparing the distance from operating point to the nose point in per-unit voltage margin
- ✓ Explain the physical meaning of Q-V curve slope (dV/dQ) and its sign change near instability
📖 Why This Matters
Voltage collapse caused by reactive power shortage has triggered blackouts worldwide—including the 2003 Northeast US/Canada blackout. In mining operations, remote substations feeding crushing plants, conveyors, and ventilation fans are especially vulnerable during high-load shifts or cable faults. Understanding the Q-V curve lets you proactively identify weak buses, size VAR compensation (e.g., SVCs or capacitor banks), and avoid unplanned shutdowns that cost millions per hour in lost production.
📘 Core Principles
Voltage stability depends on the balance between reactive power supply and demand across network impedances. The Q-V curve emerges from the nonlinear power flow equations: as reactive demand Q increases at a bus, voltage V initially declines gradually—but beyond a critical point, small ΔQ causes large ΔV or even no solution, signaling instability. This 'nose' corresponds to the saddle-node bifurcation in the power flow Jacobian. The curve shape reflects system strength: stiff systems (strong grid connection, low R/X) yield flat, wide curves; weak systems (long feeders, high impedance) show sharp, narrow noses. Sensitivity analysis (e.g., modal participation of generators or capacitors) further reveals which resources most influence the critical point.
📐 Q-V Curve Point Generation
Each point on the Q-V curve is obtained by solving the load flow with a fixed P and swept Q value at the test bus. No single closed-form formula exists—but the iterative process relies on the reactive power balance equation at bus i: Q_i = V_i Σⱼ(V_j (G_ij sinθ_ij − B_ij cosθ_ij)). For teaching, we use continuation-based or parameterized Newton-Raphson methods; however, a practical approximation for initial estimation uses the sensitivity coefficient dV/dQ ≈ −1 / (2 × V_i × B_ii), where B_ii is the diagonal element of the bus susceptance matrix.
Reactive Power Sensitivity Approximation
dV/dQ ≈ −1 / (2 × V_i × B_ii)Linearized estimate of voltage change per unit reactive power change at bus i, derived from the power flow Jacobian's reactive power submatrix.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| dV/dQ | Voltage sensitivity to reactive power | pu V / pu Q | Rate of voltage change with respect to reactive power injection at bus i |
| V_i | Voltage magnitude at bus i | pu | Per-unit voltage at the test bus |
| B_ii | Diagonal element of bus susceptance matrix | pu | Negative value typically dominates reactive power flow; more negative = stronger reactive coupling |
Typical Ranges:
Strong transmission bus: −15 to −30 pu
Weak distribution bus (mine feeder): −3 to −10 pu
💡 Worked Example
Problem: At a mine substation bus (i), base MVA = 100, nominal voltage = 34.5 kV. Load flow yields B_ii = −8.2 pu (on 100 MVA base). Operating voltage is 0.97 pu. Estimate dV/dQ near this point.
1.
Step 1: Use sensitivity approximation: dV/dQ ≈ −1 / (2 × V_i × B_ii)
2.
Step 2: Plug in values: dV/dQ ≈ −1 / (2 × 0.97 × (−8.2)) = −1 / (−15.908) = +0.0629 pu V / pu Q
3.
Step 3: Interpret sign: Positive dV/dQ indicates local stability (voltage rises with increased Q injection); magnitude implies ~0.063 pu voltage rise per pu reactive power increase.
Answer:
The estimated sensitivity is +0.0629 pu V/pu Q — confirming stable local behavior. A sign reversal (negative dV/dQ) would warn of proximity to the nose point.
🏗️ Real-World Application
At the Escondida copper mine (Chile), engineers observed repeated undervoltage trips at the 66-kV SAG mill bus during monsoon season when auxiliary loads surged. Using ETAP v22.5, they constructed a Q-V curve sweeping Q from −0.2 to +0.5 pu at the bus (P fixed at 0.85 pu). The nose point occurred at Q = 0.38 pu, V = 0.812 pu — only 0.042 pu below normal 0.854 pu. This revealed insufficient reactive reserve. They installed a 25-MVAR STATCOM, shifting the nose to V = 0.841 pu and increasing margin by 69%, eliminating trips.