iRIC を用いた土石流解析
エンジニアリング本部 防災・環境解析部 水圏解析グループ

田中 春樹

１. はじめに険性が高まってくると予想される。

降雨による斜面崩壊には、大きく分けて深層崩壊

本稿では、土石流・泥流を解析するソフトウェアのと表層崩壊の二種類ある。深層崩壊とは長期間の降ひとつとして、一般に無料で公開されている iRIC を雨により土壌中に雨水が蓄積し、基盤上までの土層用いた解析事例について紹介する。

が崩壊する現象である。一方、表層崩壊とは降雨強度が大きい場合に表層土が崩壊する現象であり、崩２. 土石流解析モデルの概要

壊と同時に表層に蓄積した水と土砂、泥が混ざった

２．１ iRIC とはiRIC（International River Interface Cooperative）ソ土石流・泥流が発生する。
近年の土石流・泥流による災害は、平成 25 年 10

フトウェアは、これまで USGS（アメリカ地質調査所）で月に東京都大島町、平成 26 年 8 月に広島県広島市

開発してきた MD_SWMS と（財）北海道河川防災研究で発生した。

センターで開発してきた RIC-Nays の機能を統合した

東京都大島町の例では、平成 25 年 10 月 16 日に河川流況・河床変動解析ソフトウェアであり、台風 26 号がもたらした湿った空気の影響で、最大

http://i-ric.org/ja/より一般に無料でダウンロードで100mm/h を超える猛烈な雨が降り続き、24 時間の降きる。

水量が 800mm を超える豪雨となり、斜面崩壊によるiRIC は、ポストプロセッサ、プリプロセッサ、ソルバ泥流が発生し、死者 35 名、行方不明者 4 名、家屋の３つの機能から構成されており、ソルバとして、二被害(全壊、半壊)86 戸となる大規模な土砂災害と次元または、三次元の流況解析や河床変動解析ななった 1）。

ど様々なソルバを搭載しており、GUI でプリプロセッ

広島県広島市の例では、平成 26 年 8 月 20 日にサ処理（メッシュ作成や地形、粗度、構造物データ作

時間 100mm を超える局地的な豪雨となり、安佐南区成など）からポストプロセッサ処理（平面分布図、コンと安佐北区では土石流（107 件）とがけ崩れ（59 件）

ター図、アニメーション）を一貫して実行することがでが同時多発的に発生し、死者 74 名、家屋被害（全壊、きる。

半壊）255 戸と甚大な人的被害をもたらした２）。

このうち土石流・泥流解析については後述する

今後想定される温暖化の影響により上述したよう

Morpho2DH ソルバを用いた。

な短時間の猛烈な豪雨の発生頻度が増加すると考えられ、表層崩壊による土石流・泥流が発生する危− 39 −

1/6

�

２．２ 解析モデル概要√�� � � �

Morpho2DH は一般曲線座標で境界適合座標を用いた非定常平面二次元の土石流・泥流モデルを主体とした解析ソルバであり、京都大学の竹林洋史
θ氏によって開発された。

��n � �=�∗ tan�� � �� �� ��n �� � � ��n ��√�� � � �（2）（3）

:流動方向の河床勾配

θｘ ：x 方向の河床勾配、θｙ ：y 方向の河床勾配

主な機能・特徴は以下の通りである。

θ� ：水及び土砂の混合物中の鉛直平均土砂濃① 斜面崩壊を初期条件とした土石流・泥流の流動・度�̅に関する流動方向の平衡河床勾配

堆積過程を表現可能なモデルである。
② 砂防ダムや貯砂ダム、家屋などの構造物を考慮できる。
③ 植生や植生高さを考慮できる。
④ 複数の崩壊地点を考慮できる。
２．３ 基礎式
水及び土砂の混合物の質量保存則 3)は以下の式

図 2 二相流モデルの概念図 5）で表される。

t

図 2 のように、水と土砂の混合物に対して河床近∂h ∂hu ∂hv E++=∂t ∂x ∂u c∗：時間、 h（1）

傍に層流域、その上に乱流域を有する二層流を考え、�̅ を一定とすると、以下の関係が得られる。

：土石流・泥流の流動深

u、v： x、y 方向における速度成分

tan �� ��∗ ： 静止堆積層の砂礫の堆積濃度
E ： 河床の浸食速度
斜面崩壊箇所v��⁄� � 1��̅ ��tan ����⁄� � 1��̅ � 1 �（4）�� ：土砂の内部摩擦角、�� ⁄�：層流層厚さの比

液体中の土砂の質量保存則 3)は下式で表される。

h����ｃｈ � ｃｈＵ �ｃｈＶ
�
��Ｅ��
��
��u（5）

運動量保存則 3)から以下の関係が得られる。
��� ���� �������
1 �� ���（6）��� ������� �� �� ����
��
��

yx

図 1 平面二次元モデルの概念図

ここで、以下の江頭らの式(2)を用い、�は流動方向の河床勾配式であることから式(3)の関係が求めら��� ���� �������
1 �� ������ ����
�（7）�� �� �� ��������ここに、� は重力、�� は河床位である。� は圧力であり、静水圧近似を用いる。�� は下式で表される。
�� � �� � ���̅ � �れる。

（8）

�� ：土砂等の密度、�：水の密度、�:土砂の密度− 40 −2/6

表 1 計算条件
�τbx � ��� � ��� ��� � � � ��対象地域√�� �� �

τby � ��� � ��� ��� � � � ���√�� � � �τ�� , τ�� : x 、y 方向における掃流力成分

広島県広島市可部東六丁目
メッシュサイズ：5m
メッシュ数（X方向：81、Y方向：111）
地形国土地理院 5mメッシュ(標高)3
斜面崩壊規模 崩壊可能土砂量は7.5m(X方向 5m、Ｙ方向 5m、渓床深さ 0.3m )浸食深さ0.1m
計算時間200秒
DT0.001（9）

計算メッシュ

（10）

fb は抵抗係数であり、以下の関係を用いる．

(乱流域)
fb �(層流域)

��� �
∝
8fb �

計算対象地域は平成 26 年 8 月豪雨により土石流（11）が発生した可部東六丁目とした。計算対象範囲には標高約 250m の山があり、その麓に家屋が存在しており、台川の左岸側の比較的地盤が高いところに家屋が密集している。

������� �� ��
4 � �� � ｃ��� ��ｋｆ�������ｃ� � （12）
�
�
25 �
�
� ��
�
�
�ｃ
�
�

斜面崩壊規模については、表 1 に示すように崩壊可能土砂量は崩壊範囲を計算メッシュ１つ分(25m2)とし、崩壊する渓床深さ 0.3m と仮定の条件とした。斜面崩壊箇所は、地図・空中写真閲覧サービス（国土ここに，��� は泥流の抵抗に関する係数であり、地理院）より平成 26 年 8 月豪雨（広島市北部）の航ある。ｋｆ =0.16、�� =0.0828、� は粒子の反発係数、�

浸食深さとは浸食が予想される渓床堆積土砂のユーザーが泥流の流れ方を見ながら設定する値では土砂の平均粒径である。
3)

河床位方程式 は以下の式で表される。
��ｂ�����
�∗

平均深さのことであり、本来、露岩調査を行い縦断的な基岩の連続性を考慮し設定する必要があるが、今（13）

回はモデルの特徴を把握するため、仮定条件として計算領域全域の平均浸食深さを 0.1m と低く設定しこれらの基礎式を用いて土石流・泥流の流動深、流速、堆積・浸食量を求める。

空写真を参考に設定した。

た。また、静止堆積濃度、土砂内部摩擦角、層流層厚の比などの土石流性質に関するパラメータはモデルデフォルト値とした。

２．４ 計算条件および検討ケース（２） 検討ケース（１） 計算条件
計算条件表を表 1 に示す。

2.2 のモデルの特徴に挙げた項目を把握するために表 2 に記す CASE1～CASE4 を検討ケースとした。
表 2 検討ケース
CASE1
CASE2
CASE3
CASE4− 41 −

3/6

家屋なし
家屋あり
対策施設あり, 家屋あり
斜面崩壊箇所複数あり

CASE1 は家屋、砂防ダムなどの障害物がない条２．５ 結果及び考察件である。

検討ケース毎の流動深（土石流の深さ）の時間変

CASE2 は家屋を考慮した条件である。家屋につい化および最終地盤変動量を図 4～8 に示す。ては iRIC の GUI にてポリゴンを作成して加えること

CASE1（家屋なし）においては、土石流発生直後ができるが、今回は国土基盤地図情報の建築物をは土石流の規模は非常に小さいが、流下とともに大シェープファイル形式に変換したものを iRIC のインきくなり、家屋に流入する時点で大きくなっていること

ポート機能より読み込み設定した。

がわかる。

CASE3 は、対策施設として砂防ダムを考慮した条

また、土石流発生から 60 秒後には家屋に到達件である。砂防ダムは、設置設定箇所を「8.20 土砂し、初期斜面崩壊深さ 0.3m のものが渓流の土石を災害 砂防・治山施設整備計画図」5)を参考に配置し、巻き込んで約 1m の流動深となったことがわかる。そ

堤体部標高を一律 E.L.108m（天端高 10m）に嵩上げの後は、家屋を考慮していないため土石流は扇形にした。

広く薄く流れており、地盤の低い河川に流れ込んで

CASE4 は、斜面崩壊箇所を複数個所設定した条いる。

件である。複数個所それぞれの斜面崩壊規模は同一とした。また、家屋、砂防ダムなどの障害物はないものとした。
図 3 に計算対象域と計算地形、家屋位置を示す。

図 4 流動深の時間変化（ＣＡＳＥ1）
CASE2（家屋あり）においては、土石流発生から60 秒後には家屋に到達し、土石流は家屋の上流側で堆積し流動深が 2m 近くせり上がっていることがわかる。また、家屋が障害物となり、土石流が道路や家
図 3 計算対象地域および地形条件

屋の間を流下していることがわかる。最終地盤変動量についても家屋なしと比較すると、最終的に家屋前面で堆積した土砂が多いことがわかる。

− 42 −

4/6

図 6 最終地盤変動量

図 5 流動深の時間変化（ＣＡＳＥ2）
CASE3（対策施設あり）においては、前述同様の土石流が発生すると仮定して、家屋上流側に砂防ダムを設定した。その結果、土石流発生から 60 秒後は砂防ダムを設定したことにより、家屋への土石流の流入はなく計算終了時までダムから越流することがなかった。
CASE4（複数斜面崩壊）においては、今回新たに

図 7 流動深の時間変化および地盤変動量（ＣＡＳＥ3）追加した斜面崩壊箇所からの土石流を支川とし、CASE4 以前の崩壊箇所からの土石流を本川とする。
土石流発生直後は 2 地点とも土石流の規模は非常に小さいが、流下とともに大きくなり、支川が 10～15秒後に本川と合流することで土石流が発達し、60 秒後の到達範囲も CASE1 と比較し広くなっており、土石流の流速が大きくなったことがわかる。その後も、本川と支川それぞれの渓谷からの土石流が家屋に流れ込み到達範囲も広くなることがわかる。
以上のように検討ケース毎に土石流の発生・流動・堆積過程の違いを視覚的に捉えることができた。
図 8 流動深の時間変化（CASE4）− 43 −

5/6

３. おわりにかあるが、出力としてアニメーションや Google Earth

土石流・泥流モデル（Morpho2DH）および iRIC の上への表示も容易に可能であることから、土石流のGUI が適用可能である事項、適用できない事項およ発生・流動・堆積過程を視覚的に捉えるツールとしてび課題点を表 3 に整理した。

は非常に優れたモデルであると考えられる。今後は実務上仕様を満たす計算条件、出力結果項目とモ

表 3 Morpho2D の適用可否項目の整理

デルの計算対象可否条件を照らし合わせて iRIC の適用可・家屋などの障害物の設定
・本川・支川の崩壊箇所の設定
・植生密度、植生高さの設定・計算結果項目は流動深、地盤変動量、流速など
・連続画像出力、アニメーション動画出力
・シェープファイルをインポート可能
・その他ソフト形式出力可（Google Earth用(.kml)、ESRIシェープファイル、CSVなど）
適用外 ・降雨による粒径の細かい土砂の浸食は非対象/課題点 ・流木や家屋流失計算は非対象
・本川・支川の崩壊開始時間の設定は不可
・最大流動深、到達時間の平面分布図は出力対象外
・ソルバーのソースは非公開

土石流・泥流モデル（Morpho2DH）を使用することを検討していきたいと考える。

＜参考文献＞iRIC の土石流・泥流モデルにおける適用外の項

1） 「平成 25 年台風第 26 号伊豆大島の土砂災害の概要」（国土交通省，平成 25 年 11 月 12 日,目を以下に示す。

本省発表資料）① 土石流による土砂が堆積したあとに降り続いた豪2） 「広島県で発生した土砂災害への対応状況」（国土交通省,平成 26 年 10 月 31 日）

雨による粒径の小さい土砂の浸食は対象として

3） 「土石流の数値シミュレーション」（日本流体力いない。
② 本川、支川の複数箇所での斜面崩壊は設定でき

学会数値流体力学部門 Web 会誌（第 12 巻，第るが、各々の発生時間設定をできず同時刻に発2 号,pp. 33-43）,2004 年,江頭進治，伊藤隆郭）4） 「伊豆大島で発生した泥流の平面二次元解析」

生させることしかできない。

（河川技術論文集(第 20 巻), 2014 年, 竹林洋

計算結果は時刻毎の流動深、変動量、速度など

史，江頭進治，藤田正治）

があるが、包絡値については出力対象外である
ため、別途 CSV 形式より値を抽出し包絡値を取

5） 「8.20 土砂災害 砂防・治山施設整備計画図」（平成 26 年 12 月 2 日,広島県）

得し、別途 GIS ソフトを用いて取り扱えるように処

6） 「火山噴火に起因した土砂災害予想区域図作理する手間がある。

成の手引き（案）」（平成 25 年 3 月,国土交通省）

土石流とは異なるが火山災害予想区域図作成にあたって 6)は最大流動深分布図、最大堆積深分布図、分布の時間変化を求められることから上記の課題は解決が望まれる。
③ iRIC の GUI 上の課題としては、結果出力時の凡例の小数点以下桁数設定ができないため、実務上の報告書に掲載する際は③と同様に別途 GISを用いて再描画する必要があると考えられる。
④ 公開されているソルバのソースは非公開であるため、ソルバ変更は不可能である。
上記のように適用外となる項目や課題点はいくつ− 44 −

6/6