🧑🎓
「潮流計算」って何ですか？電力の流れを計算するのは分かるけど、具体的に何が求まるんですか？
🎓
ざっくり言うと、発電所で作った電気が、どのくらいの電圧で、どちらの方向に、どれだけ流れているかを計算するんだ。例えば、上のシミュレーターで「バス3 負荷 P」を増やすと、バス3の電圧が下がったり、バス1からバス3への電力の流れが増えたりするのがリアルタイムで見えるよ。
🧑🎓
え、そうなんですか！じゃあ、送電線の「リアクタンス X」を変えると、何が変わるんですか？
🎓
実務では、送電線が長くなるとこのリアクタンスが大きくなるんだ。シミュレーターで「線路1-3 リアクタンス X」のスライダーを大きくしてみて。そうすると、バス1からバス3への電力の流れが制限されて、別の経路（例えばバス1→2→3）に回り込む様子が観察できるよ。これが「潮流ルート」の変化だ。
🧑🎓
なるほど！でも、バス1の電圧はなぜ1.0で固定なんですか？他のバスみたいに変えられないのはなぜ？
🎓
良いところに気づいたね。バス1は「スラック母線」と呼ばれる、システム全体の電圧と周波数の基準なんだ。発電機の出力を変えたり負荷を変えたりすると、系統全体の電力の過不足が生じる。その調整役をこのスラック母線が一手に引き受けるんだ。だから固定されているんだよ。シミュレーターで負荷を増やすと、バス1からの送電電力が自動で増えてバランスを取るのが見えるはずだ。
潮流計算の基本は、各母線（バス）で電力の収支が合う（キルヒホッフの法則）という条件から立てられる非線形連立方程式を解くことです。送電線の電力潮流は、両端の母線電圧の大きさ（V）、位相差（δ）、および送電線のリアクタンス（X）で決まります。
$$P_{ij}= \dfrac{V_i V_j \sin(\delta_i - \delta_j)}{X_{ij}}$$
$P_{ij}$: 母線iから母線jへ流れる有効電力、$V_i, V_j$: 母線i, jの電圧、$\delta_i, \delta_j$: 母線i, jの電圧位相角、$X_{ij}$: 母線i-j間の送電線リアクタンス。有効電力の流れは、主に両端の電圧の位相差によって決まります。
一方、無効電力の潮流は、電圧の大きさそのものの差に大きく影響されます。無効電力は系統内の電圧レベルを維持するために重要な役割を果たします。
$$Q_{ij}= \dfrac{V_i^2 - V_i V_j \cos(\delta_i - \delta_j)}{X_{ij}}$$
$Q_{ij}$: 母線iから母線jへ流れる無効電力。電圧が高い方の母線から低い方の母線へ無効電力が流れる傾向があります。シミュレーターで負荷の無効電力Qを増やすと、そのバスの電圧が下がるのはこのためです。
よくある誤解と注意点
まず、「リアクタンスXを小さくすれば、いつでも送電能力が上がる」というのは誤解です。確かに、送電線のリアクタンスが小さいほど、同じ電圧差で流せる電力は大きくなります。しかし、実際の送電線では、リアクタンスを小さくする（例えば太い導体を使う）と、今度は線路の静電容量（充電電流）の影響が無視できなくなります。特に長距離送電では、この充電電流によって軽負荷時に受電端の電圧が送電端より高くなる「フェランチ効果」が発生し、電圧制御が難しくなります。シミュレーターで「線路1-3」を架空線から地中ケーブルに変えたつもりで、リアクタンスXを大幅に小さくし、かつ静電容量Bを大きく設定してみてください。負荷が少ない時に、バス3の電圧がバス1を上回る現象を確認できるはずです。
次に、母線電圧は「1.0」が常に最適と思いがちですが、実務では系統全体の電圧プロファイルを考慮します。例えば、送電端の電圧を1.05[p.u.]、中間点を1.02、需要端を0.98のように、段階的に下げて設定することがあります。これは、無効電力の流れを最小化して送電損失を減らすためです。シミュレーターで全ての母線電圧を1.0に固定して負荷を変えると、無効電力の潮流が大きくなり、結果として送電損失が増加する様子を観察できます。
最後に、潮流計算の結果は「静的なスナップショット」であることを理解しましょう。この計算は、ある特定の運転状態（例えば、午後2時のピーク時）での安定状態を解いています。しかし、実際の系統では、太陽光発電の出力急変や大規模電動機の起動などによる「動的な現象」が常に発生しています。潮流計算で電圧が許容範囲内だからといって、瞬時的な電圧降下（サグ）が起きない保証にはなりません。動的現象の解析には、別の「過渡安定度解析」ツールが必要です。
関連する工学分野
この潮流計算の技術は、電力系統の「血管」の状態を診断するものと言えます。これと深く連携するのが、「系統保護・リレー工学」です。潮流計算で求められる各送電線の電流値は、過電流リレーの設定値の基礎データになります。例えば、シミュレーターで「線路1-2」の潮流が最大500Aと分かれば、リレーは600A以上でトリップするように設定します。また、潮流の方向は方向性リレーの動作判定に不可欠です。
もう一つの重要な分野は「電力システム経済学」です。ここでは、潮流計算の結果である送電損失が直接、コストに跳ね返ります。発電コストが安い石炭火力が遠隔地にあっても、そこから需要地まで送電する際の損失分を加味すると、割高になることがあります。この損失を考慮した上で、どの発電機をどの出力で動かすのが最も経済的かを決定する「経済負荷配分」の計算は、潮流計算を内包した最適化問題です。シミュレーターでバス1の発電コストを安く、バス2の発電コストを高く設定し、最適化アルゴリズムを走らせると、送電損失分も考慮して発電機出力が自動調整される、といったイメージです。
さらに応用範囲を広げると、「電力電子工学」とも密接に関わります。近年、系統に導入されるFACTS（柔軟交流送電システム）装置、例えばSTATCOM（無効電力補償装置）やUPFC（統合電力流制御装置）は、潮流計算のパラメータ（リアクタンスや電圧）を能動的に変化させることで、電力の流れを制御します。シミュレーター上で、ある母線に「理想的な無効電力注入源」を追加し、電圧を一定に保つ操作をしてみてください。それがSTATCOMの基本的な役割です。
発展的な学習のために
まず次の一歩は、「ニュートン・ラフソン法」の理解です。潮流計算の核心は、あの非線形連立方程式をどう解くかです。手順としては：1) 適当な初期値（電圧=1.0、位相角=0）から始め、2) 現在の値で計算した電力と目標値の差（ミスマッチ）を求め、3) その差をゼロに近づけるように変数を修正する、という反復計算を繰り返します。この「修正量」を求めるために、システムの感度を表すヤコビアン行列というものを計算し、連立一次方程式を解きます。数式で書くと、修正方程式は $$\begin{bmatrix} \Delta P \\ \Delta Q \end{bmatrix} = [J] \begin{bmatrix} \Delta \delta \\ \Delta V \end{bmatrix}$$ の形になります。[J]がヤコビアン行列です。このアルゴリズムの流れを追うことで、なぜシミュレーターがパラメータを変えると「即座に」新しい解が得られるのか（実は内部で高速に反復計算している）、その仕組みが腹落ちします。
次に、「直流潮流法」を学ぶと視界が開けます。これは、電圧変動や無効電力の影響を一旦無視し、有効電力の流れのみに注目した超簡略化モデルです。先ほどの有効電力の式で、位相差が十分小さい時 $\sin(\delta_i - \delta_j) \approx \delta_i - \delta_j$、電圧を1.0と近似すると、$$P_{ij} \approx \dfrac{\delta_i - \delta_j}{X_{ij}}$$ という線形の関係式が得られます。これを使うと、大規模系統の大まかな潮流分布を一瞬で計算でき、系統計画の初期検討や市場分析で重宝されます。シミュレーターの結果と比較してみると、その近似の精度を体感できます。
これらの基礎を押さえたら、実際の大規模系統データを扱ってみることをお勧めします。例えば、IEEEが提供する標準テストシステム（IEEE 14バス系統、30バス系統など）には、発電機、負荷、送電線の詳細なデータが揃っています。NovaSolverのようなツールやオープンソースソフトウェア（MATPOWER等）でこれらのデータを読み込み、自分でケーススタディ（「あの送電線が故障したらどうなるか」など）を実行することで、教科書の知識が生きた実践力に変わります。