千曲川の破堤による氾濫地区について,速報として再現計算をしてみました。 計算条件は得られる情報で以下のとおりにしました。
- 地形条件は国土地理院の5mメッシュ
- 氾濫条件は,破堤幅70m,越流水深1mで,本間の越流公式(完全越流)で計算すると,約100m3/sのため,連続式の右辺の同量を加えた。つまり,フラックス0の水塊がある状態。計算は,10 x 10mのメッシュで計算したため,0.2cmづつ5メッシュに分割して設定。
右側のツールバーのwheel zoom等で拡大,縮小できます。
計算結果
水深と流速ベクトルを図化しました。 右側のツールバーのwheel zoom等で拡大,縮小できます。
- 破堤30分後
車両基地のほうに向かって一直線です。相対的に低いようですね。
- 破堤1時間後
盛土があるのか南北方向に広がります。
- 破堤3時間後
車両基地はアウトですね。どんどん広がります。
- 破堤10時間後
浅川の左岸側まで広がりました。地形条件がおかしい?
地理院の情報と比較すると浸水範囲が狭いので,氾濫流量が多いか氾濫時間が長いのかどっちかだと思います。引き続き検証してみます。
浸水深の情報は地理院より公開されておりますので,それも参照下さい。 https://www1.gsi.go.jp/geowww/201910/shinsui/01_shinsui_chikuma.pdf
その他
- 数値計算は,現象を理解するためには有効な手段の一つなので,今後の河道計画,都市計画の立案,見直しの基礎資料になるかと思います。
- ご依頼いただければ,他地域も計算します。
ソースコード
MITライセンスですべてオープンにしております。ご活用下さい。
標高設定のためモジュールは以下の参照