Part 1
Advanced Power-System Modeling & Per-Unit Analysis

01. Engineering Analysis သည် Model မှ စတင်သည်

Substation တစ်ခုမှာ 230 kV Bus၊ 115 kV Transmission Line၊ 13.8 kV Load နဲ့ Transformer တွေ ရှိတယ်လို့ စဉ်းစားကြည့်ပါ။ Voltage Level မတူတဲ့ Equipment တွေရဲ့ Impedance ကို ဘယ်လို ပေါင်းမလဲ။ Transformer Nameplate ပေါ်က 8% Impedance က Network Model ထဲမှာ ဘယ်တန်ဖိုးဖြစ်သွားမလဲ။ Software ထဲ Data ထည့်ပြီး ရရှိလာတဲ့ Result ကို ဘယ်လို စစ်ဆေးမလဲ။

ဒီမေးခွန်းတွေရဲ့ အစမှာ Physical System ကို လေ့လာလိုသော ပြဿနာနှင့် ကိုက်ညီသည့် Mathematical Model အဖြစ် ပြောင်းလဲခြင်း ရှိပါတယ်။ Model က Physical System ရဲ့ ဘယ်သဘောသဘာဝတွေကို ထိန်းထားသလဲ၊ ဘယ်အရာတွေကို ချန်ထားသလဲဆိုတာ သိမှ Result ကို မှန်ကန်စွာ အသုံးပြုနိုင်ပါတယ်။

လိုအပ်သော အခြေခံများမှာ Complex Numbers၊ RMS Phasors၊ Kirchhoff’s Laws၊ Three-Phase Power နှင့် Matrix Multiplication တို့ ဖြစ်ပါတယ်။ ဒီ Part ကို balanced, fundamental-frequency, positive-sequence steady-state modeling ဖြင့် စတင်ပါမယ်။ Unbalanced sequence networks ကို Part 2၊ power-flow solution ကို Part 3 မှာ ဆက်လက်တည်ဆောက်ပါမယ်။

02. Phasors၊ Power နှင့် Sign Conventions

ဒီစာမျက်နှာမှာ V နဲ့ I ကို RMS complex phasors အဖြစ် သုံးပါတယ်။ j = √(−1) ဖြစ်ပြီး * က complex conjugate ကို ဆိုလိုပါတယ်။ Balanced positive-sequence system အတွက် a-phase voltage နဲ့ current ကို အသုံးပြုလျှင်—

φ ဟာ corresponding phase voltage နဲ့ phase current ကြား angle ဖြစ်ပါတယ်။ Magnitude relation မှာ √3 ကို သုံးနိုင်ပေမယ့် actual V_ab phasor နဲ့ I_a ကို complex-power equation ထဲ တိုက်ရိုက်သုံးရင် 30° phase relationship ကိုပါ ထည့်စဉ်းစားရပါတယ်။ S = 3V_anI_a^* ကို အခြေခံထားခြင်းက ဒီရှုပ်ထွေးမှုကို ရှောင်ရှားပေးပါတယ်။

Voltage angle reference ကို a-phase positive-sequence reference အဖြစ် ဆက်တိုက်သုံးပါမယ်။ Line-to-line voltage magnitude ကို base အဖြစ် သတ်မှတ်ခြင်းက phasor angle reference ကို အလိုအလျောက် 30° ရွှေ့ခြင်း မဟုတ်ပါ။

03. Per-Unit System — Base ကို အရင်သတ်မှတ်ခြင်း

Per-unit ဟာ quantity ကို ရွေးထားတဲ့ reference scale နဲ့ normalize လုပ်ခြင်း ဖြစ်ပါတယ်။ Three-phase system အတွက် independent bases နှစ်ခုကို ရွေးပါမယ်—total three-phase S_B နဲ့ line-to-line RMS V_B,LL။

Z_B equation မှာ √3 မပါလာခြင်းကို phase quantities ဖြင့် ဆင်းသက်ကြည့်နိုင်ပါတယ်။

ဒီ bases တွေနဲ့ V_pu = V_an/V_B,ph၊ I_pu = I_a/I_B လို့ ရေးရင်—

Balanced condition မှာ V_pu ရဲ့ magnitude ဟာ |V_LL|/V_B,LL နဲ့ တူပါတယ်။ S_pu equation မှာ 3 သို့မဟုတ် √3 ထပ်မထည့်ပါ။

04. Base Conversion နှင့် Transformer Ratio

Physical impedance တစ်ခုကို base နှစ်မျိုးနဲ့ ဖော်ပြလည်း actual Ω တန်ဖိုးက မပြောင်းပါ။ ထို့ကြောင့်—

Voltage bases နှစ်ခုကို transformer ရဲ့ တူညီသော side ပေါ်မှာ နှိုင်းယှဉ်ရပါတယ်။ Percent impedance ကို အရင် 100 နဲ့စားပြီး per-unit ပြောင်းပါ။

Ideal Transformer ကို ဘယ်အခြေအနေမှာ Absorb လုပ်နိုင်သလဲ။

Nominal ratio a = V_H/V_L အတိုင်း voltage bases ရွေးလျှင် Z_H = a²Z_L နဲ့ Z_B,H = a²Z_B,L ဖြစ်တာကြောင့် Z_pu,H = Z_pu,L ဖြစ်ပါတယ်။ Nominal magnitude ratio ကို base system ထဲ absorb လုပ်နိုင်ပါတယ်။

05. Study Objective အလိုက် Equipment Model ရွေးချယ်ခြင်း

Transformer nameplate က |Z|% ပေးထားတယ်ဆိုရင် R နဲ့ X ကို သီးခြားသိပြီးသား မဟုတ်ပါ။ R/X data သို့မဟုတ် load-loss information ကို သင့်လျော်သော test condition နဲ့ ဆက်စပ်အသုံးပြုရပါတယ်။ ဒီ Part ရဲ့ lossless-transformer approximation ကို actual duty သို့မဟုတ် loss study တိုင်းအတွက် မသုံးရပါ။

Constant-PQ Load နှင့် Constant-Z Load

Constant-PQ load ကို fixed admittance အဖြစ် Y_bus ထဲ တိုက်ရိုက်ထည့်ခြင်းက load model ကို ပြောင်းလိုက်ခြင်းဖြစ်ပါတယ်။ V = 1 pu မှာ တူညီတဲ့ P/Q ရအောင် ချိန်ထားလည်း voltage ပြောင်းလျှင် response မတူတော့ပါ။

06. Transmission-Line Models နှင့် Line Charging

Per-phase series parameter z′ = r′ + jωL′ [Ω/km]၊ shunt parameter y′ = g′ + jωC′ [S/km]၊ line length ℓ [km] လို့ သတ်မှတ်ပါမယ်။ r′ ဟာ သက်ဆိုင်ရာ temperature/frequency မှ AC resistance ဖြစ်သင့်ပါတယ်။ L′ နဲ့ C′ မှာ conductor geometry၊ spacing၊ bundle နဲ့ earth effects တို့ ပါဝင်နိုင်ပါတယ်။

“ဒီ km ထက်ကျော်မှ long line” ဆိုတဲ့ rule တစ်ခုတည်းနဲ့ မဆုံးဖြတ်ပါ။ Cable charging၊ voltage level၊ frequency နဲ့ study objective ကလည်း အရေးပါပါတယ်။

Nominal-π ၏ Terminal Equations

I_s ကို sending terminal မှ line ထဲဝင်သော current၊ I_r ကို line မှ receiving load ဘက် ထွက်သော current လို့ သတ်မှတ်လျှင်—

Uniform distributed line အတွက် γ = √(z′y′), Z_c = √(z′/y′) လို့ သတ်မှတ်လျှင် A = D = cosh(γℓ), B = Z_csinh(γℓ), C = sinh(γℓ)/Z_c ဖြစ်ပါတယ်။ ဒါဟာ uniform-parameter sinusoidal model အတွင်း exact ဖြစ်ပြီး electromagnetic transient model အားလုံးကို အစားမထိုးပါ။

07. Common Training Network — Four-Bus Model

အောက်ပါ data တွေဟာ ဒီစာအုပ်အတွက် ရည်ရွယ်တည်ဆောက်ထားသော illustrative training data ဖြစ်ပါတယ်။ Actual utility study သို့မဟုတ် approved relay-setting data မဟုတ်ပါ။

Assumptions register: Balanced sinusoidal operation; transformer taps = 1 pu; zero modeled phase shift; magnetizing branches and transformer losses omitted; line conductance omitted; loads modeled as current injections rather than fixed shunts. Buses 2 and 3 have no additional load/generation.

ဒီ source boundary ကို short-circuit model အဖြစ် မယူဆပါ။ Fault study အတွက် upstream Thevenin/sequence impedances ကို ထည့်ရပါမယ်။ Power-flow model မှာ prescribed slack voltage သုံးခြင်းက upstream fault strength ကို သတ်မှတ်ပြီးသား မဟုတ်ပါ။

08. Y_bus တည်ဆောက်ခြင်း — Kirchhoff’s Current Law

Bus i မှ Bus k သို့ series current ကို I_ik = y_ik(V_i − V_k) လို့ ရေးနိုင်ပါတယ်။ Bus i မှ network ထဲဝင်သော current တွေကို ပေါင်းလျှင်—

ဒီ stamping rule က unity-tap branches အတွက် ဖြစ်ပါတယ်။ Parallel branches ရှိလျှင် admittance များ ပေါင်းထည့်ရပါတယ်။ Y_bus က generally sparse ဖြစ်လို့ network ကြီးလာလျှင် sparse linear algebra ကို အသုံးပြုရပါတယ်။

Off-Nominal / Phase-Shifting Transformer

From-side complex tap t = τe^jθ ကို V_internal = V_from/t လို့ define လုပ်ပြီး series admittance y၊ total branch charging jb_c ကို သုံးလျှင် branch stamp က—

t = 1 ဖြစ်လျှင် simple π stamp ကို ပြန်ရပါတယ်။ θ ≠ 0 ဖြစ်လျှင် generally Y_ft ≠ Y_tf ဖြစ်နိုင်ပါတယ်။ ဒီ tap convention နဲ့ branch model ကို MATPOWER User’s Manual, Section 3.2 — Branches တွင်လည်း ဖော်ပြထားပါတယ်။

09. Network Model မှ Operating Equations သို့

Y_bus ရရှိတာက network ရဲ့ voltage-current relationship ရရှိခြင်း ဖြစ်ပါတယ်။ Constant-power loads ပါတဲ့ network မှာ unknown bus voltages တွေကို သိပြီးသား မဟုတ်ပါ။

Unknown V ဟာ expression ရဲ့ နှစ်ဖက်လုံးမှာ ပါနေတာကြောင့် power-flow problem က nonlinear ဖြစ်ပါတယ်။ Bus 1 ကို reference/slack boundary လုပ်ပြီး Buses 2, 3 မှာ zero net power injection၊ Bus 4 မှာ S_inj,4 = −0.20 − j0.10 pu သတ်မှတ်နိုင်ပါတယ်။ အဲဒီ equations ကို Part 3 မှာ numerical solution ရှာပါမယ်။

10. Network Reduction — Boundary Behavior ကို ထိန်းထားခြင်း

Internal buses အားလုံးကို study တိုင်းမှာ ထိန်းထားဖို့ မလိုပါ။ Retained buses r နဲ့ eliminated buses e လို့ partition လုပ်ရင်—

Y_ee invertible ဖြစ်ပြီး eliminated buses မှာ independent injection I_e = 0 ဖြစ်လျှင် V_e ကို eliminate လုပ်နိုင်ပါတယ်။

ဒါကို Kron Reduction / Schur Complement လို့ ခေါ်ပါတယ်။ I_e ≠ 0 ဆိုရင် reduced equation မှာ injection term ကိုပါ ထိန်းရပါတယ်—

Reduction က သတ်မှတ်ထားသော linear model အတွင်း boundary voltage-current relation ကို ထိန်းပါတယ်။ Internal relay locations၊ local voltages သို့မဟုတ် nonlinear/dynamic behavior ကို အလိုအလျောက် ထိန်းပေးခြင်း မဟုတ်ပါ။ Software implementation မှာ explicit inverse ထုတ်ခြင်းထက် linear solves ကို အသုံးပြုသင့်ပါတယ်။

11. Engineering Verification နှင့် Common Mistakes

Calculation မှန်ကန်မှုကို independent checks များနဲ့ စစ်ရပါမယ်။ ဒီ training network အတွက် အောက်ပါ checks တွေက model structure ကို တိုက်ရိုက်စစ်ဆေးပေးပါတယ်။

1. Base consistency: Transformer နှစ်ဖက်စလုံးမှ converted per-unit impedance တူရမယ်။
2. Reciprocity: ဒီ unity-tap, non-phase-shifting example မှာ Y_bus = Y_bus^T ဖြစ်ရမယ်။ Complex symmetric ဖြစ်ခြင်းက Hermitian ဖြစ်ခြင်းနဲ့ မတူပါ။
3. Row sums: ဒီ network ရဲ့ row sums က [0, j0.0066125, j0.0066125, 0] ဖြစ်ရမယ်။ Shunt မရှိလျှင် row sums zero ဖြစ်မယ်။
4. Flat-voltage check: V = [1, 1, 1, 1]^T ဆိုရင် series currents zero ဖြစ်ပြီး line-charging currents ပဲ ကျန်ရမယ်။ ဒါဟာ model check ဖြစ်ပြီး loaded system ရဲ့ solution မဟုတ်ပါ။
5. Real losses: Line series current အတွက် P_loss,pu = r|I_series,pu|² ≥ 0 ဖြစ်ရမယ်။ Actual total three-phase loss = 3R|I_{series,actual}|² ဖြစ်ပါတယ်။ Charging ပါလျှင် terminal current နဲ့ series current မတူနိုင်ပါ။
6. Dimensional check: Ω → pu အတွက် Z_B နဲ့စား၊ siemens → pu အတွက် Z_B နဲ့မြှောက်။

Y_bus က အမြဲ invertible လို့ မယူဆပါ။ Shunt/reference treatment မပါသော connected passive series network ရဲ့ full nodal matrix က singular ဖြစ်နိုင်ပါတယ်။ Reference bus ကို သတ်မှတ်ပြီး သက်ဆိုင်ရာ reduced equations ဖြင့် ဖြေရှင်းရပါတယ်။

12. Practice Problems နှင့် Answer Checks

Answer ကို ကြည့်မီ assumption၊ base နဲ့ units ကို ကိုယ်တိုင် အရင်ရေးပါ။

1. Base calculation: 100 MVA, 34.5 kV system အတွက် I_B နှင့် Z_B ကို ရှာပါ။
2. Transformer conversion: 25 MVA, 115/13.8 kV, 10% impedance magnitude ရှိ transformer ကို 100 MVA base ပေါ် ပြောင်းပါ။ Resistance neglected လို့ ယူပါ။
3. Voltage-base change: 100 MVA, 115 kV base ပေါ် 0.10 pu ရှိ impedance ကို 100 MVA, 110 kV base ပေါ် ပြောင်းပါ။ Same physical side ကို သုံးပါ။
4. Model comparison: V = 1 pu မှာ 20 MW စားသုံးသော load ကို voltage 0.90 pu သို့ လျှော့လျှင် constant-P model နှင့် constant-Z model တို့၏ real power ကို နှိုင်းယှဉ်ပါ။
5. Y-bus change: Worked network ထဲက line charging ကို ဖြုတ်လိုက်လျှင် ဘယ် matrix entries တွေ ပြောင်းမလဲ။
6. Engineering judgment: Same MVA/kV/|Z|% ရှိ transformer နှစ်လုံးကို fault duty နဲ့ differential security studies အားလုံးမှာ identical လို့ ယူဆနိုင်သလား။ အကြောင်းပြပါ။

Part 1 မှ ရရှိသည့် Engineering Deliverable

ယခုအဆင့်မှာ base table၊ assumptions register၊ equipment model table၊ per-unit one-line diagram နှင့် verified Y_bus ရရှိထားပါတယ်။ ဒီအခြေခံကို Part 2 မှာ sequence networks အဖြစ် တိုးချဲ့ပြီး Part 3 မှာ operating voltages နဲ့ power flows ကို ဖြေရှင်းပါမယ်။

အဓိကသဘောတရား: Model တစ်ခု၏ တန်ဖိုးကို Equation ဘယ်လောက်ရှုပ်ထွေးသလဲဆိုတာထက်၊ မေးထားသော Engineering Question ကို လိုအပ်သည့် Accuracy ဖြင့် ဖြေနိုင်သလားဆိုတာဖြင့် အကဲဖြတ်ရပါတယ်။

References & Further Reading

အောက်ပါရင်းမြစ်များကို foundational reading နှင့် implementation reference အဖြစ် အသုံးပြုနိုင်ပါတယ်။ ဒီစာမျက်နှာရှိ training data၊ worked network နှင့် numerical examples များကို သီးခြားတည်ဆောက်တွက်ချက်ထားပါတယ်။

1. MIT OpenCourseWare — James L. Kirtley, Introduction to Electric Power Systems: Readings။ Polyphase Networks နှင့် network-analysis အခြေခံများကို ဆက်လက်လေ့လာရန်။
2. MATPOWER User’s Manual — Chapter 3, Modeling။ π branch model၊ complex tap နှင့် network formulation ကို စစ်ဆေးရန်။
3. MATPOWER — makeYbus source code။ Bus-admittance matrix တည်ဆောက်ပုံကို implementation နှင့် နှိုင်းယှဉ်ရန်။