確率論的ダム決壊および堤防決壊シミュレーション: 動的地盤侵食モデリングと流体力学的波のルーティング

ダム、貯水池、洪水防止堤防システムは、現代社会において最も重要な民間インフラ資産の一部です。農業灌漑用の水の供給を確保し、再生可能な水力発電を生成し、歴史的な氾濫原に住む何百万もの住民を浸水から守ります。しかし、土製のダムや堤防が壊滅的に決壊すると、数分で膨大な量の位置エネルギーが放出され、下流のコミュニティ、橋、産業回廊全体を壊滅させる可能性のある暴力的で壊滅的な洪水波が発生します。

歴史的に、ダム決壊リスク評価は、一次元 (1D) 水力経路と組み合わせた静的で経験的な決壊公式 (古典的なフレーリッヒ方程式、マクドナルド-ラングリッジ-モノポリス方程式、フォン トゥーン-ジレット方程式など) に依存していました。これらのレガシー フレームワークでは次のことが行われます。

• エンジニアは、履歴回帰散布図に基づいて、固定の破損形状 (台形の破損形状、底部の幅、側面の傾斜) と固定の破損の開発時間 (たとえば、1.5 時間) を仮定しました。
• 決壊流出は、2D の氾濫原の谷の広がり、土砂輸送の嵩密度、峡谷の壁からの波の反射を完全に無視した 1D 断面を通って下流に送られました。

この静的な経験的アプローチには、危険な工学的限界があります。土の堤防は、任意の幾何学的スケジュールに沿って崩壊することはありません。複雑で非線形な場合は失敗します地盤侵食力学: 進行性の後方越流侵食、内部配管導管のせん断侵食、および回転斜面安定性スリップサークルの崩壊。さらに、静的な仮定では、地質工学的パラメータの不確実性（土壌の圧縮凝集力、粒度分布（D_{50})、臨界せん断応力(tau_c）、貯水池の水位の変動。

2026 年後半には、規制上の安全基準 (連邦エネルギー規制委員会 (FERC) および米国陸軍工兵隊 (USACE) によって義務付けられています) により、確率的物理ベースのダム決壊および堤防決壊シミュレーション。一時的な地盤工学的ブリーチ侵食の物理学を、GPU で高速化された 2 次元 (2D) 浅水方程式ソルバーおよびモンテカルロ アンサンブル不確かさの定量化と直接組み合わせることで、ダムの安全技術者は、妥協のない数学的精度で下流の浸水の危険性をマッピングできます。

この包括的な技術論文では、最先端のダム決壊評価を実行するために必要な、地盤工学的侵食物理学、堆積物を含む非ニュートン波の流体力学、モンテカルロアンサンブルサンプリング、および緊急避難ハザードマッププロトコルを検討しています。

1. 地盤破壊力学:越水と配管内部侵食

土堤ダムは主に、水圧越水または内部地盤工学的配管という 2 つの物理的メカニズムのいずれかによって破壊されます。

過剰せん断応力浸食方程式

物理ベースの破壊モデル (USDA-ARS WinDAM や DL-Breach など) では、破壊境界に沿った土壌の剥離は次の条件によって支配されます。過剰せん断応力の配合:

ε_r = k_d (τ_e - τ_c)^a

どこ：

• イプシロン_r土壌剥離侵食速度 (m/s または cm/hr) を表します。
• k_dは土壌浸食係数 (cm^3/N・s）、現場ジェット浸食試験（JET）または水中ピンホール分散試験によって決定されます。
• tau_eは、流出境界に流れる水によって及ぼされる有効な流体力学的せん断応力です。

τ_e = γ_w R_h S_f

どこガンマwは水の単位重量、R_hは突破水路の水力半径であり、S_fはエネルギーの傾きです。

• tau_cは、圧縮された盛土土の臨界せん断応力 (Pa) であり、これを下回ると剥離が発生しません。
• $a$ は経験的な指数であり、通常は $a = 1.0$ と見なされます。

流水が発達中の切れ込みを通って加速すると、tau_e500 Pa を超えることもあり、土壌の臨界せん断強度を大幅に上回ります (tau_c ≈ 5 - 20 Pa) - 指数関数的な体積浸食速度を促進します。

2. 動的流出ハイドログラフ計算：幅広堰とオリフィスの水力学

ブリーチノッチが広く深くなると、瞬間的な流出が発生します。Q_{違反}(t)の間で一時的に切り替わります浸漬オリフィス流（初期の内部配管時）および自由表面の幅の広い堰の流れ（一度ブレーク屋根が崩壊すると）。

flowchart TD
    Pool["Reservoir Pool Level: H_res(t)"] --> Geom["Instantaneous Breach Notch Geometry<br/>• B_w(t), Z_s, H_b(t)"]
    Geom -->|Internal Piping Phase| Pipe["Pressurized Orifice Flow<br/>• Q = C_d · A_pipe(t) · √(2g · ΔH)"]
    Geom -->|Open Breach Crest Phase| Weir["Broad-Crested Weir Flow<br/>• Q = C_w · B_w(t) · (H_res - Z_b)^1.5 + C_s · Z_s · (H_res - Z_b)^2.5"]
    Pipe --> Wave["Downstream 2D Hydrodynamic Flood Wave Routing"]
    Weir --> Wave

複合台形堰の形成

開いたブリーチトレンチが確立されると、長方形と三角形の広頭堰を組み合わせた方程式によって流出量が計算されます。

Q_{違反}(t) = C_{w1} B_w(t) ≤ft( H_{res}(t) - Z_b(t) )^{3/2} + C_{w2} Z_s ≤ft( H_{res}(t) - Z_b(t) )^{5/2}

どこ：

• B_w(t)拡張するブリーチノッチの瞬間的な底幅 (m) です。
• Z_b(t)ブリーチノッチの瞬間的な底面高さであり、岩盤または非浸食性基礎層に到達するまで下方に浸食されます。
• Z_sは台形の切れ目の側面の傾斜です (Z_s = (デルタ x) / (デルタ z)、通常、粘着性の粘土の場合は $0.5$ から $1.0$、非粘着性の砂の場合は最大 $1.5$）。
• C_{w1}、C_{w2}は、浸水放水効果のために校正された堰流量係数です。

同時に貯水池の標高もH_{レス}(t)リザーバーステージとストレージの保存方程式に従って動的に低下します。

frac{dH_{res}}{dt} = frac{Q_{流入}(t) - Q_{漏洩}(t) - Q_{放水路}(t)}{A_{res}(H_{res})}

どこA_{レス}(H_{レス})標高における貯水池の表面積です。H_{レス}.

3. 流体力学的波のルーティング: 2D 非ニュートンダム破壊波

ダムの決壊波は、きれいな透明な水ではありません。それは、水、粉砕された堤防の砂利、浮遊した堆積物、そして根こそぎにされた木々が混ざり合った、激しく高密度のスラリーです。

急峻な山の峡谷や狭い川の谷では、堆積物の体積濃度（C_v) $20\% - 45\%$ を超える可能性があり、洪水の波が非ニュートン超濃縮スラリー.

非ニュートンビンガム塑性流体力学

>

見かけのせん断応力の計算式:

τ = τ_y + µ_B · (∂u / ∂z)

>

パラメータ:

• τ_y: 堆積物スラリーの降伏応力 (Pa) — 流体運動を開始するには克服する必要があります

• μ_B：ビンガム動塑性粘度（Pa・s）

• かさ密度: ρ_slurry = ρ_w · (1 - C_v) + ρ_s · C_v

2D 衝撃捕捉流体力学ソルバー

数値的不安定や非物理的な拡散を発生させずに、この暴力的で高速の波を下流に送るには、次のようにします。

• モデルは 2D 浅水方程式を利用します。ビンガムプラスチックまたはハーシェル・バルクレーの流れ抵抗のソースターム (タウ_{bx}、タウ_{by}).
• ソルバーは、急峻な問題を解決するために二次 TVD (Total variation Diminishing) ショック捕捉スキームを強制します。ダムブレイク波ボアを超える速度で進む12～20m/秒下流の谷を下ります。
• 運動量方程式は、波が川の急な曲がり角を曲がる際の横方向の遠心力を解析し、峡谷の壁に衝突する超高波の波の遡上を正確に予測します。

4. 確率的モンテカルロアンサンブルフレームワーク: 地質工学的不確実性の定量化

法医学ダムの安全工学では、正確な地質工学的特性が 100% 確実に分かることはありません。土盛りの圧縮品質、含水率、凝集性、および内部均質性は、直線の長さによって異なります。

現代のエンジニアは、非現実的で保守的、または危険なほど楽観的な単一の「最悪の」決定論的シミュレーションを実行するのではなく、モンテカルロ確率アンサンブル.

flowchart TD
  Params["Input Parameter Probability Distributions\n• Reservoir Head: Normal(μ = Full Pool, σ = 1.2m)\n• Soil Erodibility k_d: Log-Normal(μ = 10 cm³/N-s, σ = 0.4)\n• Critical Shear Stress τ_c: Weibull(λ = 8.5 Pa, k = 2.1)\n• Cohesion c': Normal(μ = 25 kPa, σ = 5 kPa)\n• Breach Width B_w: Beta distribution (Froehlich)"]
  LHS["Latin Hypercube Sampling (LHS)\nGenerates 500 Stratified Scenarios"]
  GPU["Automated GPU Solver Cluster\n500 Coupled 2D Simulations (4.2 GPU Hours)"]
  Hazard["Probabilistic Hazard Outputs\nP10, P50, P90 Inundation Maps & Arrival Time Bands"]

  Params --> LHS --> GPU --> Hazard

確率論的な洪水浸水指標

システムは、500 回の実行アンサンブルから、空間的に分布した超過確率を抽出します。

• P90 浸水マップ: 信頼水準 $90\%$ の保守的な洪水境界。民間インフラ保護設計 (原子力発電所、病院複合施設、化学薬品貯蔵庫) に使用されます。
• P50 浸水マップ: 予想されるシナリオの中央値。自治体の区域設定、保険料率の評価、および主要な避難ルートに利用されます。
• 波の到達時間の不確実性: 特定の下流人口密集地への高速波面の可能な限り最も早い到達を定量化します (T_{到着、P10}）、緊急事態管理者に避難許可期間の防御可能な安全マージンを提供します。

5. 包括的なダム破壊方法の比較マトリックス

次の表は、従来の経験的手法と、最新の物理ベースの確率的 2D ダム決壊シミュレーション フレームワークを対比させたものです。

6. 現実世界のエンジニアリングケーススタディ: 危険性の高い土製ダムの安全性監査

重要な民間インフラ資産に関するこの確率論的枠組みを検証するために、FloodSim は、住民 95,000 人が暮らす都市部の谷の 6 マイル上流に位置する、32,000 エーカーフィート (3,950 万立方メートル) の飲料水および治水池を貯水する高さ 85 フィート (26 メートル) の土堤ダムの緊急行動計画 (EAP) ダム決壊監査を実行しました。

生命の安全に関する重要な発見:

1. 危険な「ラグタイム」幻想: 従来の 1D モデルでは、最初の洪水波が下流の 6 マイルを市内中心部まで移動するのに 65 分かかると予測していました。対照的に、高解像度 2D GPU ソルバーは、渓谷の狭まりと本流の河川回廊を流れる超臨界流を考慮して、高速水力掘削が都市中心部に衝突することを明らかにしました。わずか38分。もし緊急事態管理者が従来のモデルに依存していたら、避難指示は 30 分近く遅れて届いていたでしょう。
2. 州間橋のスーパー洗車: 2D シミュレーションは、超過流速をマッピングしました。6.2m/秒州間高速道路の橋の橋脚で、波が到達してから 12 分以内に橋脚の基礎を破壊する底部せん断応力を計算します。
3. 緊急行動計画の再調整：市の緊急事態管理局は、正確な確率的到着曲線を備え、4 つの避難サイレンの位置を変更し、高速道路の自動ゲート閉鎖を確立し、認定された 22 分間の避難許可プロトコルを確立し、完全な公共生活の安全を確保しました。

7. 下流域の避難危険分類: オーストラリアの危険脆弱性曲線

公共の生命安全リスクを評価するには、生の洪水深 ($h$) と速度 ($v$) を人間の安定性の危険カテゴリーに変換する必要があります。秒速 3.0 メートルで流れる深さ 0.4 メートルの浅い水流は、流体力学的抗力によって大人の歩行者や自動車が道路から押し流されるため、高さ 1.0 メートルの水たまりよりもはるかに致死的です。

最新の緊急行動計画では、オーストラリアの緊急事態管理ハザード曲線 (H1 ～ H6):

<!--ec:block {"type":"table","hasHeaderRow":true,"rows":[{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"危険性分類カテゴリ (H1 - H6)"}],"isHeader":true}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"総流体力学強度: D \\times V = h \\times \\sqrt{u^2 + v^2} (m²/s)"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"スパン","テキスト":"""}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"[H1: すべての人に安全]: D \\times V \\le 0.4 \m²/s, h \\le 0.3m"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"安全に歩くことができる。子供や高齢者も安全です。"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"スパン","テキスト":"""}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"[H2: 小型には危険]: D \\times V \\le 0.6 \m²/s, h \\le 0.5m"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"小さなお子様や体の弱いお年寄りには安全ではありません。"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"スパン","テキスト":"""}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"[H3: 車両にとって危険]: D \\times V \\le 0.8 \m²/s, h \\le 0.6m"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"セダンや軽自動車はトラクションと浮きを失います。"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"スパン","テキスト":"""}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"[H4: すべての人にとって安全ではない]: D \\times V \\le 1.0 \m²/s, h \\le 1.0m"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"大人は足を払いました。すべての車両が浮いています。"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"スパン","テキスト":"""}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"スパン","text":"[H5: 構造的損傷]: D \\times V \\le 4.0 \m²/s, h \\le 4.0m"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"補強されていない石積みの建物は損傷を受けます。"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"スパン","テキスト":"""}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"type":"span","text":"[H6: 壊滅的な崩壊]: D \\times V > 4.0 \m²/s"}]}]},{"type":"tableRow","cells":[{"type":"tableCell","content":[{"_type":"span","text":"住宅建物の完全な構造破壊"}]}]}]} -->

動的 GIS レイヤーをエクスポートすることで、H_1を通してH_6決壊したダムの下流で 5 分ごとに分類することで、緊急事態管理者はどの避難道路が安全に走行可能であり、どの時点で致命的な死の罠に移行するかを特定できます。

8. ダムの安全技術者のための規制遵守プロトコル

2026 年にダム崩壊の危険性評価を実施するダムの安全技術者、地盤工学コンサルタント、地方自治体の公益事業責任者は、次の 5 段階のコンプライアンス チェックリストに従う必要があります。

1. 現場噴流浸食試験 (JET) の実施: 浸食係数を決して推定しないでください (k_d）教科書の土壌分類表だけから。標準化された ASTM D5852 / D6665 テストを実行して、真の過剰せん断パラメータを校正します。
2. 動的地盤破壊拡大のモデル化: 静的な台形の仮定を破棄します。後方ヘッドカット移動、側斜面の回転崩壊、貯留層プールのドローダウンをモデル化する結合された物理ベースの侵食アルゴリズムを採用します。
3. GPU アクセラレーションによる 2D ショック キャプチャ ソルバーの導入: 高解像度の LiDAR から派生した 2D 非構造化メッシュ全体でルート違反流出波を解析し、重要な下流の谷の地形、道路の堤防、暗渠が確実に水圧調整されるようにします。
4. 地盤工学パラメータの不確実性を定量化する: 層化されたラテン ハイパーキューブ パラメーター分布にわたって少なくとも 250 回のモンテカルロ アンサンブル実行を実行し、防御可能な P10、P50、および P90 の危険コンターを生成します。
5. 深さ速度曲線による下流の避難ルートの分類: 流体力学的出力を標準に変換しますH_1 - H_6人命と車両の危険性を分類し、緊急事態管理責任者に人命を救う実用的な情報を提供します。

地質工学的土壌力学、GPU で高速化された非ニュートン浅水流体力学、確率的モンテカルロ不確実性定量化の規律ある統合を通じて、現代のダム決壊シミュレーションは、土木技術者が想像できる最も壊滅的な破壊シナリオから社会を守ることができるようにします。