🎓 Lesson 2
D2
Power Flow Equations: Deriving the Fundamental P-Q Mismatch Functions
Power flow equations tell us how much real (P) and reactive (Q) power must be supplied at each bus in a power system to keep voltages stable and flows balanced—like balancing a budget for electricity.
🎯 Learning Objectives
- ✓ Calculate real and reactive power mismatches (ΔP and ΔQ) at any bus given voltage phasors and admittance matrix entries
- ✓ Explain how the Jacobian matrix links voltage changes to power mismatches in Newton-Raphson load flow
- ✓ Apply polar-form power flow equations to formulate P-Q mismatches for a 3-bus radial system
- ✓ Analyze the impact of voltage angle vs. magnitude on active vs. reactive power flow sensitivity
📖 Why This Matters
In mining operations, reliable on-site power distribution is critical—for conveyor drives, ventilation fans, and automated drilling rigs. A single undetected power mismatch can cascade into voltage collapse, unplanned shutdowns, or equipment damage. Understanding P-Q mismatch functions isn’t just theory—it’s the first line of defense in designing robust mine electrical networks and validating protection coordination studies.
📘 Core Principles
Power flow begins with modeling each bus (node) using complex voltage V_i = |V_i|∠θ_i and complex power S_i = P_i + jQ_i. From nodal analysis, S_i = V_i · I_i*, and since I_i = Σ_j Y_ij V_j, substituting yields the fundamental P-Q expressions: P_i = Σ_j |V_i||V_j|(G_ij cosθ_ij + B_ij sinθ_ij), Q_i = Σ_j |V_i||V_j|(G_ij sinθ_ij − B_ij cosθ_ij). These equations are nonlinear due to trigonometric terms and voltage products—making them unsolvable analytically for large systems, hence the need for iterative methods like Newton-Raphson, where mismatches ΔP_i = P_i^spec − P_i^calc and ΔQ_i = Q_i^spec − Q_i^calc drive convergence.
📐 Key Calculation
The P-Q mismatch functions quantify the difference between scheduled (specified) and computed power injections. They serve as residuals in numerical solvers—and their structure determines convergence speed and robustness.
💡 Worked Example
Problem: For Bus 2 in a 3-bus system: V_1 = 1.05∠0° pu, V_2 = 1.00∠−0.0524 rad (−3°), V_3 = 0.98∠−0.0873 rad (−5°); Y_22 = 5.0 − j15.0 pu, Y_21 = −1.0 + j4.0 pu, Y_23 = −4.0 + j11.0 pu. Scheduled P₂ = 0.8 pu, Q₂ = 0.3 pu. Compute ΔP₂ and ΔQ₂.
1.
Step 1: Compute θ₂₁ = θ₂ − θ₁ = −0.0524 − 0 = −0.0524 rad; θ₂₃ = θ₂ − θ₃ = −0.0524 − (−0.0873) = 0.0349 rad
2.
Step 2: Evaluate P₂ = |V₂|²·G₂₂ + |V₂||V₁|(G₂₁ cosθ₂₁ + B₂₁ sinθ₂₁) + |V₂||V₃|(G₂₃ cosθ₂₃ + B₂₃ sinθ₂₃) = (1.00)²(5.0) + (1.00)(1.05)[(−1.0)cos(−0.0524) + (4.0)sin(−0.0524)] + (1.00)(0.98)[(−4.0)cos(0.0349) + (11.0)sin(0.0349)] ≈ 5.0 − 1.096 − 3.857 = 0.047 pu
3.
Step 3: Similarly compute Q₂ = |V₂|²·(−B₂₂) + |V₂||V₁|(G₂₁ sinθ₂₁ − B₂₁ cosθ₂₁) + |V₂||V₃|(G₂₃ sinθ₂₃ − B₂₃ cosθ₂₃) ≈ −(1.00)²(−15.0) + (1.00)(1.05)[(−1.0)(−0.0524) − (4.0)(0.9986)] + (1.00)(0.98)[(−4.0)(0.0349) − (11.0)(0.9994)] ≈ 15.0 − 4.152 − 10.921 = −0.073 pu
4.
Step 4: ΔP₂ = 0.8 − 0.047 = 0.753 pu; ΔQ₂ = 0.3 − (−0.073) = 0.373 pu
Answer:
ΔP₂ = 0.753 pu and ΔQ₂ = 0.373 pu — both significantly nonzero, indicating the initial voltage guess is far from solution and requires Newton-Raphson correction.
🏗️ Real-World Application
At the Escondida copper mine (Chile), a 2022 load flow study revealed persistent Q-mismatches > 0.25 pu at the 33 kV substation feeding the SAG mill drive. Field measurements confirmed leading VAR absorption due to oversized capacitor banks. By reformulating the Q-mismatch function to isolate shunt susceptance contributions and re-tuning reactive compensation, engineers reduced voltage deviation from ±4.2% to ±0.8%, improving motor torque stability and reducing harmonic resonance risk per IEEE Std 519-2022.