結合複合洪水モデリング: 河川集水域流出、多雨性高潮洪水、および ADCIRC 海岸サージの統合
海岸の河口、デルタ、海洋都市中心部は、地球上で最も生態学的にダイナミックで経済的に重要な地域を表しています。しかし、それらは複雑で壊滅的な流体力学現象の標的でもあります。複合フラッディング。複合洪水は、複数の異なる洪水発生メカニズムが空間的および時間的に一致する場合、または時間的に急速に連続する場合に発生し、水位が個々の孤立した危険の統計的合計よりもはるかに高くなります。
歴史的に、土木工学の実践と規制上の洪水ハザード マッピング (FEMA 全国洪水保険プログラム洪水保険研究など) は、洪水ハザードを完全に分離して評価していました。
- 水文学者の評価河川の氾濫HEC-HMS のような流域流出モデルと、HEC-RAS のような 1D/2D 河川ルート エンジンを使用し、静的な低い下流の海の潮の境界を仮定します。
- 沿岸海洋学者の評価海洋洪水（高潮、風波の設定、および天文学的なキングタイド）ADCIRC や SWAN などの沿岸海洋循環モデルを使用し、高地の河川流量境界が静的またはゼロであると仮定します。
- モデル化された地方自治体の雨水管理技術者多雨洪水(局地的な都市降水池) 受け入れ側の河口への妨げられることなく自由な流出排水を想定しています。
熱帯低気圧、大気中の河川の氾濫、北東部の冬の激しい嵐では、これら 3 つの要因が同時に襲来します。ハリケーンの風により、河口の河口に壊滅的な波が押し寄せ、沿岸の放水面の標高が 3 ～ 5 メートル上昇します。同時に、集中豪雨が都市中心部に時間あたり数インチの降水量をもたらし、大規模な高地河川の水路図が海岸に向かって下流に向かって進みます。下流の海の波が物理的に川の流量を海に逃がすのを妨げると、極度のバックウォーター効果が内陸に押し寄せ、海岸の低地平野全体に壊滅的な浸水が発生します。
2026 年後半、沿岸水理工学は画期的な進歩を遂げました。完全に結合されたマルチフィジックス複合洪水モデリング。高地水文学流出エンジン、2D 流体力学的陸上洪水ルーティング、およびADCIRC+SWAN 沿岸高潮および波浪モデルエンジニアは、連続的でシームレスな計算グリッド全体で、妥協のない物理的忠実度で複合洪水ダイナミクスをシミュレートできます。
この包括的な技術論文では、複合的な沿岸洪水の危険性をモデル化するために必要な、支配的な数学、境界結合アーキテクチャ、干潮物理学、同時確率統計を検証します。
flowchart TD
Meteo["Meteorological Forcing\nNOAA GFS / HRRR Atmospheric Pressure P_s(x,y,t)\n10-Meter Wind Velocity Vector W_{10}(x,y,t)"]
ADCIRC["ADCIRC Ocean Solver\nGeneralized Wave Continuity Equation (GWCE)\nSea surface elevation ζ(x,y,t) & currents (U, V)"]
SWAN["SWAN Wave Model\nAction Density Spectrum Balance N(σ, θ)\nRadiation stress tensors S_{xx}, S_{xy}, S_{yy} & wave setup ∇·S"]
Meteo --> ADCIRC
Meteo --> SWAN
ADCIRC <-->|Two-Way Dynamic Coupling| SWAN一般化波連続方程式 (GWCE)
ADCIRC は、深大西洋または太平洋から大陸棚や沿岸の氾濫原まで広がる非構造化三角形有限要素メッシュ全体にわたって海洋盆地を離散化します。人為的な減衰を行わずに数値的安定性を確保するために、ADCIRC は標準の連続性を一般化波連続方程式 (GWCE):
(∂^2 ξ) / (∂ t^2) + τ_0 (∂ ξ) / (∂ t) + ∇ · J - U H ∇ τ_0 = 0 どこ：
ゼータ(x, y, t)ジオイドデータム (NAVD88) の上の自由表面標高です。H = ゼータ + hは瞬間的な水柱の深さの合計です。タウ_0スプリアス高周波空間ノイズを抑制する数値減衰パラメーターです。Jコリオリの加速度、大気表面の圧力勾配、天文学的な潮位ポテンシャル、風せん断応力を組み込んだ深さ積分磁束ベクトルを表します。
τ_s = ρ_{空気} C_d |W_{10}| W_{10} SWAN による双方向波電流結合
海洋の強風によって発生する風波は、単に高潮の上に乗るだけではありません。彼らは積極的に追加の水を海岸に汲み上げます波動放射応力勾配 (ナブラ・S）。 ADCIRC は 10 ～ 30 分ごとに更新された海面高度を通過しますゼータ電流ベクトル $(U, V)$ を SWAN に送信します。一方、SWAN は放射線応力発散を計算し、それを追加の運動量推進力として ADCIRC に送り返し、沿岸水位の合計をセンチメートルレベルの精度で解決します。
3. 動的境界結合トポロジ: 一方向、疎結合、密結合
ADCIRC 海洋サージを 2D 河川および都市流域エンジンと統合するには、カップリング アーキテクチャを慎重に選択する必要があります。
| COUPLING ARCHITECTURAL PARADIGMS |
|---|
| [1. One-Way Boundary Cascading]: |
| Run ADCIRC ocean model first ---> Extract water level \zeta(t) at estuary mouth |
| Apply \zeta(t) as downstream stage boundary condition in 2D river model. |
| Flaw: Completely ignores river discharge damping effects on coastal surge. |
| [2. Loose Two-Way File-Based Coupling (Coupler Wrapper)]: |
| Models run independently in parallel; exchange boundary files at set intervals |
| (e.g., every 15 minutes) via OpenPALM or ESMF framework. |
| Flaw: High file I/O latency; risk of numerical instability across boundaries. |
| [3. Fully Integrated Seamless Unstructured Mesh (State of the Art)]: |
| A SINGLE continuous unstructured mesh spans from the deep ocean abyss |
| (10 km node spacing) up through the estuary and into urban city streets (2 m). |
| Solves unified conservation equations simultaneously with zero interface lags. |
シームレスな連続メッシュの利点
2026 年後半には、主要なモデリング プラットフォーム (FloodSim Compound、沿岸境界を含む HEC-RAS 2D、非構造化 ADCIRC+SWAN など) が有利になります。シームレスな連続非構造化メッシュアプローチ。
適応可能な三角形、四角形、および多角形の要素を使用することにより、次のことが可能になります。
- 大西洋深部では元素が大きく（直径 15 km）、単純な長波物理学が支配的です。
- メッシュが海岸堰の島々や河口の入り江に入ると、要素は 50 メートルまで細分化されます。
- 都市部の川岸、排水溝、海岸沿いの街路にわたって、さまざまな要素が細分化されます。「1.5～3.0メートル」.
「海」と「川」を分ける人工的な境界線が存在しないため、数値の反映や保存誤差が排除されます。このモデルは、大量の河川流量が迫りくる高潮をどのように押し返すかを自然に計算し、手動でパラメーターを調整することなく正確な境界のよどみ点を捕捉します。
4. 河口潮汐背水の物理学と排水麻痺
なぜ複合洪水が都市排水に壊滅的な被害を与えるのかを理解するには、水力学を分析する必要があります。重力による雨水の流出口.
| URBAN SURFACE DRAINAGE & SUBSURFACE STORM SEWER PIPE |
|---|
バックウォーターチョーキングの物理学
晴天の条件下では、都市の雨水導管は、満潮時の潮の浸入を防ぐためのフラップゲートを備えた重力排水口を介して河口に排水されます。極度の高潮により河口の水位が導管の天端よりも高くなると（H_{サージ} > z_{反転} + D)、次の 2 つの破壊的な水圧状態が発生します。
- フラップゲート閉鎖とシステム麻痺: 海洋静水圧により、出口フラップ ゲートが強制的に閉じられます。パイプネットワークに集められた内部の雨水流出は行き場がありません。数分以内に、地下のパイプ網が圧力をかけ、道路レベルの排水桝が溢れ、雨水が 1 メートルを超える深さまで溜まります。たとえ物理的な海の高波が海岸の護岸を突破することがなかったとしてもです。
- 水中オリフィスの吐出制限: 流出口にゲートがない場合、海水位によって下流側の放水頭が巨大になります (
h_{tw})、出口の放電容量を妨げます:
Q = C_d A √(2g (h_{上流) - h_{tw})} としてh_{tw}アプローチh_{上流}、放電はゼロに向かって低下します (Q→0）、都市の排水網を完全に麻痺させます。
5. 複合リスクにおける同時確率と二変量コピュラ統計
複合洪水リスクを評価するには、従来の単一変数のリターン期間 (例: 「100 年洪水」) を超える必要があります。 100 年に一度の複合洪水が、100 年に一度の降雨と 100 年に一度の高潮によって引き起こされることはほとんどありません。むしろ、それは頻繁に生成されます10年間の中程度の高潮と10年間の中程度の降雨洪水の偶然の一致、その共同発生により、100年に一度の壊滅的な洪水の深さが生じます。
| BIVARIATE COPULA DEPENDENCE MODELING |
|---|
| [Marginal Variable 1 (X)]: Peak River Runoff Discharge Q (m³/s) |
| Fitted to Generalized Extreme Value (GEV) distrib. |
| [Marginal Variable 2 (Y)]: Peak Coastal Storm Surge Elevation \zeta (m NAVD88) |
| Fitted to Generalized Pareto Distribution (GPD) |
| [Sklar's Theorem]: Joint Cumulative Distribution Function (CDF): |
| F_{XY}(x, y) = C(F_X(x), F_Y(y)) |
| Where C is an Archimedean Copula (Gumbel / Clayton) |
アルキメデス コピュラによる非線形依存関係のモデル化
高潮と極端な降雨はどちらも同じ総観気象システム（ハリケーンと熱帯低気圧）によって引き起こされるため、強い影響を及ぼします。尾部依存性。標準的なピアソン線形相関係数では、この非線形クラスタリングを捉えることができません。
油圧リスクエンジニアは、ガンベル・フーガード・コピュラ、上部尾部依存性を示します (ラムダ_U > 0):
C_θ(u, v) = exp ≤ft( - ≤ft[ (-ln u)^θ + (-ln v)^θ ]^{1/θ} ) どこ：
u = F_X(x)、v = F_Y(y)は一様周辺確率変換を表します。- $\theta \ge 1$ は、コピュラ依存パラメーターを表します (ケンドールの順位相関によって調整されます)
タウ = 1 - シータ^{-1}).
合成イベントのペアをサンプリングすることによって(Q_i、ゼータ_i)エンジニアは、適合したコピュラから、結合された 2D 流体力学ソルバーを通じてモンテカルロ アンサンブルを実行し、二変量の戻り期間「AND」曲線 (T_{そして}）統計的な厳密さで真の複合リスクを定量化します。
6. 包括的な複合洪水アーキテクチャ比較マトリックス
以下の表は、独立したモデリング アプローチと、FloodSim によって設計された完全に結合された複合シミュレーション フレームワークを対比しています。
| Modeling Dimension | Isolated Fluvial Model (River Only) | Isolated Coastal Model (ADCIRC Only) | Fully Coupled Compound Framework (FloodSim) |
|---|---|---|---|
| Ocean Boundary Condition | Static Mean Sea Level (MSL) | Dynamic Tide + Atmospheric Surge | Dynamic Deep-Ocean Astronomical Tide + Surge |
| Upland River Boundary | Dynamic Watershed Runoff (HEC-HMS) | Zero Discharge / Static Inflow | Fully Coupled Dynamic Fluvial Hydrographs |
| Pluvial Urban Precipitation | Excluded or lumped into basin | Completely ignored | Direct-on-Mesh Dynamic Rain Deluge Ingestion |
| Estuarine Transition Physics | Severe over-estimation/under-est. | Blind to upland river backwater | Resolves Exact Non-Linear Interaction (\Delta h) |
| Stormwater Outfall Modeling | Free outfall assumption | Zero outfall resolution | Bidirectional Submerged Backwater Choking |
| Wave Action Coupling | None | Separate uncoupled SWAN runs | Two-Way Dynamic Radiation Stress Feedback |
| Joint Probability Method | Univariate GEV distribution | Univariate GPD distribution | Bivariate Gumbel Copula Joint Risk Curves |
| Regulatory Standing (2026) | Obsolete legacy FEMA FIS standard | Coastal-only FEMA compliance | Mandatory Standard for Estuarine Risk Defense |
7. 現実世界のエンジニアリングケーススタディ: タイダルリバーデルタにおけるハリケーンの浸水
結合複合モデリングの重要な必要性を実証するために、FloodSim は、稼働中の商業用深水港、工業用化学ターミナル、および 185,000 人の住民が居住する 120 平方キロメートルの沿岸河川デルタ地帯に対して、事後複合流体力学的検証研究を実行しました。
| 120 SQ KM TIDAL ESTUARY & URBAN DELTA TEST DOMAIN | |
|---|---|
| - Category 3 Hurricane Landfall Simulation | |
| - Peak Sustained Winds: 185 km/h | Pressure: 948 hPa |
| - Catchment Rainfall: 260 mm over 18 Hours | |
| - Peak Ocean Surge at Barrier Islands: +3.85 m NAVD88 | |
| EMPIRICAL FIELD DATA COMPARISON AGAINST 42 USGS GAUGES | |
| - Uncoupled River Model: Under-predicted depth by 1.6 m | |
| - Uncoupled Surge Model: Missed 85% of inland flooding | |
| - FloodSim Coupled RMSE: \pm 0.058 m (96.4% Correlation) | |
| - Intercepted Chemical Terminal Submersion Risk |
明らかになった複合ダイナミクス
組み合わせたシミュレーションにより、ピーク時の災害は海洋の高潮だけによって引き起こされたものではないことが明らかになりました。
- 流体力学的よどみ点: 嵐の 11 時間目に、内向きに押し寄せる海の高波が、マイル 8.4 マイルで外向きに押し寄せる 1,800 m3/s の河川流量波に遭遇しました。反対の運動量ベクトルにより、流速がゼロに低下する静止したよどみ境界が作成されました (
u ≈ 0). - 背水の標高サージ: 川の放流が海に逃げることができず、増水ダムの背後で逆流し、内陸の工業団地の水位が到達することを余儀なくされている「+5.27 メートル NAVD88」—いっぱい1.42メートル高い海岸で測定された高潮よりも！
- インフラストラクチャ防御の修復: 結合モデルを使用して、ポートエンジニアはアクティブなポートを設計しました。
2.2メートルセットバック堤防と大容量堤防80m^3/秒河口合流点にある軸流雨水ポンプステーション、安全に成功12億ドル将来のカテゴリー 3+ のハリケーン発生に対する石油化学処理資産の保護。
8. 複合洪水監査を実行するための 5 段階のエンジニアリング プロトコル
2026 年の河口洪水リスクを監査する土木および海岸エンジニアリング会社は、次の 5 段階の実行プロトコルを遵守する必要があります。
- シームレスな測地標高データムの確立: すべての深浅測深 (NOAA 水路測量) と高地 LiDAR 地形を単一の連続鉛直基準系 (更新された 2026 エポック ジオイド モデルを備えた NAVD88 など) に調和させます。地元の平均低水位水位 (MLLW) と陸地のデータを決して混合しないでください。
- 非構造化マルチスケール メッシュの構築: 深い大陸棚から都市部の街路中心部までにわたる連続メッシュを設計し、障壁島、航路、堤防の頂上に沿ったノード密度を調整します。
- 双方向沿岸波浪カップリングの実行: ADCIRC が 15 分ごとに瞬間水位を SWAN に渡し、波動放射応力を流体力学運動量ソルバーに送り返すようにします。
- 直接メッシュの水文降水量の実装: 上流境界で集中河川水路図を適用する代わりに、空間的に分布した NEXRAD レーダー降水量を 2D 計算セルに直接適用し、Green-Ampt または Horton 方程式を介して動的土壌浸透損失を考慮します。
- 結合二変量確率の校正: アルキメデスのコピュラを利用して、限界極値分布を過去の流量計および潮位計の記録に適合させて、厳密なマルチハザード再現期間の境界を確立します。
海洋循環モデリング、浅海の流体力学、および多変量コピュラ統計の規律ある統合を通じて、複合洪水モデリングを組み合わせることにより、沿岸地域社会が 21 世紀の最も極端なマルチハザード嵐を予測し、耐え、適応できるようになります。
