論文概要：配管網における天然ガスの定常流れの解析

論文タイトル: 非線形微分代数方程式系の解法による配管網における天然ガスの定常流れ
著者: Shriram Srinivasan, Kaarthik Sundar (ロスアラモス国立研究所)

1. 問題の背景と課題

天然ガス配管網の定常状態流れのシミュレーションは、ネットワーク運用の最適化、配送能力の推定、災害計画などにおいて不可欠です。従来の研究では、重力および非線形慣性効果（対流加速度）を無視することで、支配方程式が非線形代数方程式系に簡略化され、解析的に積分可能な形（閉形式の第一積分）を得るアプローチが主流でした。

しかし、以下の理由から、これらの効果を無視したモデルには限界があります：

1. 重力の影響: 配管網は広範囲に及ぶため、山岳地帯や平地を横断し、配管が非水平になることが一般的です。勾配と流れの方向によって、重力は流れを加速または抵抗させ、圧力分布に大きな影響を与えます。
2. 非理想気体の状態方程式: 実際の天然ガスは理想気体とは異なり、圧力と密度の関係が非線形です。非理想気体の状態方程式（CNGA 式など）を用いる場合、慣性項や重力項を考慮すると、微分方程式の解析的な第一積分（閉形式解）が得られなくなります。
3. 既存研究の限界: 既存の研究の多くは、慣性や重力のいずれかを無視するか、単一パイプに限定されており、大規模ネットワークにおいて非理想気体・重力・慣性のすべてを考慮した定常流れのシミュレーション手法は存在しませんでした。

本研究の目的は、重力と慣性効果を無視せず、非理想気体の状態方程式を含む定常流れを大規模ネットワークでシミュレーションする手法を確立し、これらの効果の重要性を定量的に評価することです。

2. 提案手法と数学的枠組み

2.1 支配方程式

パイプ内の定常流れは、質量保存則と運動量保存則に基づき、以下の一次常微分方程式（ODE）と代数方程式の連立方程式（微分代数方程式系：DAE）として記述されます。

• 質量流量 $f$ は一定。
• 圧力 $p(x)$ の空間分布は、摩擦、重力、慣性項を含む ODE で記述されます。
• 非理想気体の状態方程式（$p = p(\rho)$）が用いられます。

2.2 数値解法：ニュートン・ラプソン法と前方感度方程式

第一積分が得られない場合、従来の代数方程式系への帰着は不可能です。そこで、以下の戦略を提案しています：

1. 残差関数の構築:
   各パイプの ODE を数値的に解き、入口条件（圧力 $p_i$ と流量 $f_{ij}$）から出口での予測圧力 $p_{ij}(L)$ を計算します。これを実際の出口圧力 $p_j$ と比較する残差関数 $F(p_i, p_j, f_{ij}) = 0$ を定義します。これにより、問題が有限次元の非線形方程式系（未知数はノード圧力とエッジ流量）のゼロ点探索問題として定式化されます。

2. ヤコビアン行列の評価（前方感度方程式）:
   ニュートン・ラプソン法を適用する際、残差関数のヤコビアン（偏微分係数）を計算する必要があります。これには、ODE の解に対する初期条件および流量の感度を求める「前方感度方程式（Forward Sensitivity ODEs）」を連立して解く手法を採用しました。これにより、微分方程式の構造を保持したまま効率的な反復計算が可能になります。

3. 変数変換によるスケーリング改善:
   ODE の数値的硬さ（stiffness）とスケーリングの問題を緩和するため、圧力 $p$ を $p^3$ などの変数 $\pi$ に変換する手法を提案しました。これにより、反復計算中の数値的不安定性を低減しています。

4. 2 点コロケーション法による初期値推定:
   ニュートン法が収束するための良い初期値を得るため、ODE を 2 点（入口と出口）でのみ満たすように近似する「2 点コロケーション法」を用いて粗い近似解を求め、これをニュートン反復の初期値として使用します。

3. 主要な成果と結果

3 つのケーススタディを通じて、提案手法の有効性と物理現象の特性を評価しました。

ケーススタディ 1: 既存データセット（GasLib）での検証

• 水平パイプ（重力無視）の理想気体モデルを用い、提案手法で得られた解を解析的な第一積分式と比較しました。
• 結果: コロケーション法による初期値推定が高精度であり、ODE 定式化との整合性が確認されました。また、この条件下では慣性効果が無視できることも確認されました。

ケーススタディ 2: 単一パイプにおける重力・慣性効果の分離評価

• ヤマル・ヨーロッパパイプラインのデータに基づき、勾配角を変化させて単一パイプの流れをシミュレーションしました。
• 慣性効果: どの勾配角、どの状態方程式（理想・非理想）においても、慣性項による圧力変化は極めて小さく（0.01% 未満）、実用上は完全に無視可能であることが結論付けられました。
• 重力効果: 勾配角の変化に対して圧力変化は非線形に増加し、特に傾斜が大きい場合や非理想気体の場合、その影響は顕著でした。重力を無視すると、出口圧力の推定に大きな誤差が生じます。

ケーススタディ 3: 大規模ネットワーク（米国北西部パイプライン網）への適用

• コロラド州からワシントン州にかけて広がる 9,228 km のネットワーク（977 ノード、1,033 エッジ）に提案手法を適用しました。
• 結果:
  • 約 80% のノードにおいて、重力を考慮することによる圧力差は 2% 未満でした。
  • しかし、一部のノードでは 7-8% の差が生じ、配管網が容量限界近くで運用されることを考慮すると、重力効果を考慮した高忠実度シミュレーションの必要性が確認されました。
  • 提案手法は、非理想気体と重力を考慮した大規模ネットワークの定常解を安定して計算できることを実証しました。

4. 結論と意義

1. 物理的知見: 天然ガス配管網の定常流れにおいて、慣性効果（対流加速度）は一般的に無視できるが、重力効果は配管の勾配に応じて重要であり、特に非理想気体の条件下では無視できないことを初めて大規模ネットワークで定量的に示しました。
2. 計算手法の革新: 解析的な第一積分が得られない非線形微分代数方程式系（DAE）に対しても、ニュートン・ラプソン法と前方感度方程式を組み合わせることで、効率的かつスケーラブルな数値解法を確立しました。
3. 汎用性: 提案された枠組みは天然ガスに限定されず、水素輸送、CO2 パイプライン、地域熱供給システムなど、非線形状態方程式と重力効果を持つ他の流体ネットワークにも適用可能です。

本研究は、複雑な物理現象を含む大規模ネットワークの定常状態解析において、従来の代数方程式への簡略化に依存せず、微分方程式の構造を保持したまま高精度なシミュレーションを行うための重要な基盤を提供しています。