🎓 Lesson 3
D2
Bus Admittance Matrix Construction & Sparsity
The bus admittance matrix is a mathematical map that shows how easily electricity can flow between different points (buses) in a power system, like a wiring diagram written in numbers.
🎯 Learning Objectives
- ✓ Construct a Y-bus matrix from a given single-line diagram and line parameters
- ✓ Analyze sparsity patterns to identify and exploit computational efficiency in large-scale systems
- ✓ Explain how transformer taps and shunt elements affect Y-bus entries
- ✓ Apply Kron reduction to eliminate non-essential buses while preserving equivalent admittance relationships
📖 Why This Matters
In modern mining operations, reliable power delivery to haul trucks, crushers, and ventilation fans depends on stable grid integration—often via dedicated substations or microgrids. Accurately modeling how faults or load changes propagate requires solving load flow equations—and the Y-bus matrix is the essential 'engine' behind every commercial load flow solver. Its sparsity directly determines whether a 500-bus mine grid can be solved in milliseconds or minutes.
📘 Core Principles
The Y-bus is built from physical network components: transmission lines (modeled as π-sections), transformers (with tap ratios), and shunt elements (capacitors/reactors). Sparsity arises because each bus connects only to its immediate neighbors—so for a 1000-bus system with average degree 3–5, >99.5% of Y-bus entries are zero. This structural sparsity enables efficient storage (e.g., compressed sparse row format) and solution techniques (e.g., LU factorization with fill-in control). Understanding how topology maps to matrix structure is foundational for both modeling accuracy and real-time stability assessment.
📐 Y-Bus Assembly Rule
For an n-bus system, the Y-bus is assembled by summing contributions from each branch. Self-admittances accumulate all connected branch admittances; mutual admittances are negative branch admittances between directly connected buses. Transformer modeling requires incorporating off-nominal tap ratios (a ≠ 1) into off-diagonal terms and adjusting diagonal terms accordingly.
💡 Worked Example
Problem: A 3-bus system has: Bus 1–2: line with y₁₂ = 2.0 − j6.0 S; Bus 2–3: line with y₂₃ = 1.5 − j4.5 S; Bus 1–3: no direct connection. All shunt admittances neglected. Construct the 3×3 Y-bus.
1.
Step 1: Compute diagonal elements: Y₁₁ = y₁₂ = 2.0 − j6.0; Y₂₂ = y₁₂ + y₂₃ = (2.0 − j6.0) + (1.5 − j4.5) = 3.5 − j10.5; Y₃₃ = y₂₃ = 1.5 − j4.5
2.
Step 2: Compute off-diagonals: Y₁₂ = Y₂₁ = −y₁₂ = −2.0 + j6.0; Y₂₃ = Y₃₂ = −y₂₃ = −1.5 + j4.5; Y₁₃ = Y₃₁ = 0 (no direct link)
3.
Step 3: Assemble matrix and verify symmetry and zero-row-sum property for lossless case (here, losses included via conductance)
Answer:
Y-bus = [[2.0−j6.0, −2.0+j6.0, 0], [−2.0+j6.0, 3.5−j10.5, −1.5+j4.5], [0, −1.5+j4.5, 1.5−j4.5]]. Diagonal dominance holds; sparsity pattern matches network connectivity.
🏗️ Real-World Application
At Newmont’s Boddington Mine (Western Australia), a 132 kV internal distribution grid serving 12 major loads (crushers, conveyors, SAG mills) was modeled using a 27-bus Y-bus matrix. Engineers used sparse matrix solvers in PSS®E to perform daily contingency screening—identifying weak buses during peak loading. When a 33 kV feeder fault occurred, the pre-computed Y-bus enabled <80 ms voltage profile updates, triggering automated reactive compensation before process trips occurred.
📋 Case Connection
📋 Industrial Plant Power Design: Aluminum Smelter Load Flow Optimization
Severe voltage sag during anode changing cycles causing PLC trip cascades
📋 Hospital Power Systems: Emergency Generator Load Flow Validation
Voltage collapse observed during full transfer test due to motor inrush overwhelming AVR response