🙋
電力潮流計算って何ですか？送電線に流れる電力が計算できるって聞いたけど。
🎓
大まかに言うと、発電所や変電所（母線）を結ぶ電力系統の中で、電圧や位相、電力の流れを求める計算だよ。例えば、新しい太陽光発電所をどこに接続するか検討する時、既存の系統が壊れないか調べるのに必須なんだ。このシミュレーターでは、上の「収束判定値ε」のスライダーを動かして、計算の厳しさを変えてみよう。値を小さくすると精度は上がるけど、計算回数が増えるよ。
🙋
「ガウス・ザイデル法」って何が特別なんですか？反復的に計算するということですか？
🎓
その通り！この方法は、最初に適当な電圧値を仮定して、母線の方程式を一つずつ順番に解き、解いたばかりの新しい値を使って次の母線を解いていくんだ。実務では小規模な系統や、初期設計の段階でよく使われるね。シミュレーターで「反復回数」が増えていく様子を見ると、電圧値が少しずつ真の値に近づいていくのが分かるよ。
🙋
母線に「スラック」「PV」「PQ」って種類があるみたいだけど、どう使い分けるんですか？
🎓
良い質問だね。スラック母線は系統全体の電力バランスを取る調整役で、電圧の基準になる。大きな発電所がこれにあたるよ。PV母線は電圧を一定に保ちたい発電所（例えば火力発電所）。PQ母線は消費地（工場や街）で、消費する電力が決まっている。このシミュレーターで「基準MVA」の値を変えると、全ての電力値がこの基準に対して「パーユニット」で表示される。実データを扱う時、この考え方は特に重要だ。
加速係数（1.5〜1.7程度）を導入して収束を促進できます。また、初期電圧推定値を系統の実運用値に近づける、または母線の種類（スラック・PV・PQ）を正しく設定することも重要です。反復回数上限を増やす前に、まずこれらの調整を試してください。
スラック母線は系統全体の電力収支を調整する基準母線で、電圧の大きさと位相角（通常0°）を固定します。設定時は、必ず1つの母線をスラックに指定し、その母線の有効電力・無効電力は計算結果として自動算出されるため入力不要です。
主な原因は、母線の電力指定値（Psch, Qsch）が実際の系統容量を超えているか、アドミタンス行列の入力ミスです。特に変圧器のタップ比や線路のインピーダンス値が現実的か確認してください。また、PV母線の電圧指定値が高すぎる場合も過大潮流の原因になります。
本ツールは3母線専用です。ただし、原理的にはガウス・ザイデル法のアルゴリズムを拡張すれば任意の母線数に対応可能ですが、本実装では母線数が固定されています。4母線以上の計算が必要な場合は、別の汎用潮流計算ツールをご利用ください。
系統計画・設計：新しい発電所（太陽光、風力など）や大規模工場を既存の送電網に接続する際、電圧低下や設備過負荷が起きないかを事前にシミュレーションします。ガウス・ザイデル法は初期検討の段階でよく用いられます。
系統運用：日々の電力需要の変化に応じて、どの発電所をどの出力で運転するか（経済負荷配分）、または停電復旧時の系統再構築手順を検討する際に、潮流計算結果が基盤データとなります。
スマートグリッド解析：需要家側に分散配置された電源（DER）や蓄電池の制御方策を評価します。従来とは電力の流れが双方向になるため、より高度な潮流解析が必要です。
CAEツール連携：PSS®EやDIgSILENT PowerFactoryといった商用の電力系統解析ソフトウェアの内部アルゴリズムの基礎となっています。大規模系統では計算速度の速いニュートン・ラフソン法などが使われますが、基本原理は同じです。
ガウス・ザイデル法を初めて使う時に陥りがちなポイントをいくつか挙げておくよ。まず「収束判定値εは小さければ小さいほど良い」という誤解。確かにεを1e-8とかにすれば精度は上がるけど、計算時間が爆増するし、実は計算機の丸め誤差の影響を受けて逆に収束しなくなることもあるんだ。実務では、電圧の精度を0.1%程度（ε=1e-3〜1e-4）に設定すれば十分なケースが多い。このシミュレータでεを変えながら反復回数がどう変わるか、体感してみて。
次に初期値の与え方。教科書では「定格電圧1.0∠0°で初期化」と書いてあるけど、実際の系統で大きな負荷がかかっているPQ母線などは、電圧が大きく低下していることもある。そんな時に1.0から始めると、収束までに時間がかかったり、最悪発散したりする。経験則として、PQ母線の初期電圧は少し低め（例えば0.95∠0°）に設定すると収束が早まることもあるよ。
最後に、この手法の限界を理解しておくこと。ガウス・ザイデル法は計算がシンプルでメモリ使用量が少ないのが利点だけど、収束性はシステムの「対角優位性」に大きく依存するんだ。つまり、送電線の抵抗(R)とリアクタンス(X)の比(R/X比)が大きい系統（例えば配電系統）では、更新式の分母にある$Y_{ii}$が小さくなりがちで、計算が不安定になることがある。実務では、このような系統には後述するニュートン・ラフソン法が必須になるんだ。