(57)【要約】【課題】各対象地点において発生する斜面崩壊を正確に予測することができる斜面崩壊予測方法及び斜面崩壊予測装置を提供する。【解決手段】半減期をパラメータとする実効雨量の現時点までの時系列値を、対象地点の過去の降雨イベントの降雨量データを用いて算出し、長半減期、短半減期をパラメータとする第１及び第２実効雨量をＸ軸、Ｙ軸とするＸＹ座標平面に、第１実効雨量の時系列値の最大値を通りＹ軸と平行な第１境界線、第２実効雨量の時系列値の最大値を通りＸ軸と平行な第２境界線を設定して判定用ＸＹ座標平面を作成する。次に予測対象降雨イベントの降雨量データから、長半減期、短半減期をパラメータとする第３及び第４実効雨量の時系列値を、予測対象降雨イベントの進行と共に順次算出し、同じ時刻における第３実効雨量及び第４実効雨量の時系列値をＸ値、Ｙ値とする測定点が、判定用ＸＹ座標平面において第１境界線を超えた領域、又は第２境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定する。【選択図】図９
2015232537.tif 000002
【特許請求の範囲】
【請求項１】
対象地点における降雨イベントによる斜面崩壊の発生を予測する斜面崩壊予測方法であって、
a) ある時刻における蓄積雨量が半減するまでの時間である半減期をパラメータとして、過去の降雨イベントが前記蓄積雨量に及ぼす影響の持続性を表した雨量指標である実効雨量の現時点までの時系列値を、前記対象地点の過去の降雨イベントの降雨量データを用い、所定の長半減期及び該長半減期より短い短半減期についてそれぞれ算出し、
b) 前記長半減期をパラメータとする第１実効雨量をＸ軸、前記短半減期をパラメータとする第２実効雨量をＹ軸とするＸＹ座標平面に、前記第１実効雨量の時系列値の最大値を通り前記Ｙ軸と平行な第１境界線及び前記第２実効雨量の時系列値の最大値を通り前記Ｘ軸と平行な第２境界線を設定して判定用ＸＹ座標平面を作成し、
c) 斜面崩壊発生の予測対象となる降雨イベントの降雨量データから、前記長半減期をパラメータとする第３実効雨量の時系列値及び前記短半減期をパラメータとする第４実効雨量の時系列値を、前記予測対象降雨イベントの進行と共に順次算出し、
d) 同じ時刻における前記第３実効雨量の時系列値及び前記第４実効雨量の時系列値をそれぞれＸ値、Ｙ値とする点を測定点としたとき、該測定点が、前記判定用ＸＹ座標平面において前記第１境界線を超える領域に位置するとき、又は前記第２境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定することを特徴とする斜面崩壊予測方法。
【請求項２】
請求項１に記載の斜面崩壊予測方法において、
前記対象地点の過去の降雨イベントの降雨量データから求められた、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットし、
前記複数の点のうちある時刻の点を注目点としたとき、該注目点とＸ値が同じ点の中に前記注目点よりもＹ値が大きい点が存在しない場合、及び前記注目点よりもＸ値が大きい点の中に前記注目点とＹ値が同じ点及び前記注目点よりもＹ値が大きい点のいずれもが存在しない場合は、前記注目点を第３境界点として抽出する作業を繰り返し、
抽出された全ての第３境界点をＸ値が小さい順に繋いだ線を、前記第１実効雨量及び前記第２実効雨量の相乗効果による実効雨量の最大値を表す第３境界線として前記判定用ＸＹ座標平面に設定し、
前記判定用ＸＹ座標平面において、前記測定点が前記第３境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定することを特徴とする斜面崩壊予測方法。
【請求項３】
請求項２に記載の斜面崩壊予測方法において、
複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、各時刻における判別結果の集計から、斜面崩壊の発生を予測することを特徴とする斜面崩壊予測方法。
【請求項４】
請求項２に記載の斜面崩壊予測方法において、
複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、少なくとも１個の判定用ＸＹ座標平面において測定点が第２〜第５領域のいずれかに位置するとき、斜面崩壊の発生の可能性が高いと予測することを特徴とする斜面崩壊予測方法。
【請求項５】
請求項２に記載の斜面崩壊予測方法において、
複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、
その結果を表すシンボルを、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成することを特徴とする斜面崩壊予測方法。
【請求項６】
請求項１に記載の斜面崩壊予測方法において、
前記対象地点の過去の降雨イベントの降雨量データから求められた、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットし、
前記複数の点のうちある時刻の点を注目点としたとき、該注目点とＸ値が同じ点の中に前記注目点よりもＹ値が大きい点が存在しない場合、及び前記注目点よりもＸ値が大きい点の中に前記注目点とＹ値が同じ点及び前記注目点よりもＹ値が大きい点のいずれもが存在しない場合は、前記注目点を第３境界点として抽出する作業を繰り返し、
抽出された全ての第３境界点をＸ値が小さい順に繋いだ線を、前記第１実効雨量及び前記第２実効雨量の相乗効果による実効雨量の最大値を表す第３境界線として前記判定用ＸＹ座標平面に設定し、
前記判定用ＸＹ座標平面において、前記測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出することを特徴とする斜面崩壊予測方法。
【請求項７】
請求項６に記載の斜面崩壊予測方法において、
複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出し、
該単位時間あたりの降雨量を、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成することを特徴とする斜面崩壊予測方法。
【請求項８】
対象地点における降雨イベントによる斜面崩壊の発生を予測する斜面崩壊予測装置であって、
a) ある時刻における蓄積雨量が半減するまでの時間である半減期をパラメータとして、過去の降雨イベントが前記蓄積雨量に及ぼす影響の持続性を表した雨量指標である実効雨量の現時点までの時系列値を、前記対象地点の過去の降雨イベントの降雨量データを用い、所定の長半減期及び該長半減期より短い短半減期についてそれぞれ算出する、既往実効雨量算出手段と、
b) 前記長半減期をパラメータとする第１実効雨量をＸ軸、前記短半減期をパラメータとする第２実効雨量をＹ軸とするＸＹ座標平面に、前記第１実効雨量の時系列値の最大値を通り前記Ｙ軸と平行な第１境界線及び前記第２実効雨量の時系列値の最大値を通り前記Ｘ軸と平行な第２境界線を設定して判定用ＸＹ座標平面を作成する、判定用ＸＹ座標平面作成手段と、
c) 斜面崩壊発生の予測対象となる降雨イベントの降雨量データから、前記長半減期をパラメータとする第３実効雨量の時系列値及び前記短半減期をパラメータとする第４実効雨量の時系列値を、前記予測対象降雨イベントの進行と共に順次算出する、予測用実効雨量算出手段と、
d) 同じ時刻における前記第３実効雨量の時系列値及び前記第４実効雨量の時系列値をそれぞれＸ値、Ｙ値とする点を測定点としたとき、該測定点が、前記判定用ＸＹ座標平面において前記第１境界線を超える領域に位置するとき、又は前記第２境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定する判定手段と
を備えることを特徴とする斜面崩壊予測装置。
【請求項９】
請求項８に記載の斜面崩壊予測装置において、
前記判定用ＸＹ座標平面作成手段が、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットし、前記複数の点のうちある時刻の点を注目点としたとき、該注目点とＸ値が同じ点の中に前記注目点よりもＹ値が大きい点が存在しない場合、及び前記注目点よりもＸ値が大きい点の中に前記注目点とＹ値が同じ点及び前記注目点よりもＹ値が大きい点のいずれもが存在しない場合は、前記注目点を第３境界点として抽出する作業を繰り返し、抽出された全ての第３境界点をＸ値が小さい順に繋いだ線を、前記第１実効雨量及び前記第２実効雨量の相乗効果による実効雨量の最大値を表す第３境界線として前記判定用ＸＹ座標平面に設定し、
前記判定手段が、前記判定用ＸＹ座標平面において前記測定点が前記第３境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定することを特徴とする斜面崩壊予測装置。
【請求項１０】
請求項９に記載の斜面崩壊予測装置において、
前記既往実効雨量算出手段が、長さが異なる３種類以上の半減期について前記実効雨量の時系列値を算出し、
前記判定用ＸＹ座標平面作成手段が、前記３種類以上の半減期から２種類の半減期を選択することにより形成される、複数の、長半減期と短半減期の組み合わせについて、それぞれ前記第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
前記判定手段が、前記複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、各時刻における判別結果の集計から、斜面崩壊の発生を予測することを特徴とする斜面崩壊予測装置。
【請求項１１】
請求項９に記載の斜面崩壊予測装置において、
前記既往実効雨量算出手段が、長さが異なる３種類以上の半減期について前記実効雨量の時系列値を算出し、
前記判定用ＸＹ座標平面作成手段が、前記３種類以上の半減期から２種類の半減期を選択することにより形成される、複数の、長半減期と短半減期の組み合わせについて、それぞれ前記第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
前記判定手段が、前記複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、少なくとも１個の判定用ＸＹ座標平面において測定点が第２〜第５領域のいずれかに位置するとき、斜面崩壊の発生の可能性が高いと予測することを特徴とする斜面崩壊予測装置。
【請求項１２】
請求項９に記載の斜面崩壊予測装置において、
前記判定用ＸＹ座標平面作成手段が、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、
その結果を表すシンボルを、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成する二次元判定図作成手段を有することを特徴とする斜面崩壊予測装置。
【請求項１３】
請求項８に記載の斜面崩壊予測装置において、
前記判定用ＸＹ座標平面作成手段が、前記対象地点の過去の降雨イベントの降雨量データから求められた、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットし、前記複数の点のうちある時刻の点を注目点としたとき、該注目点とＸ値が同じ点の中に前記注目点よりもＹ値が大きい点が存在しない場合、及び前記注目点よりもＸ値が大きい点の中に前記注目点とＹ値が同じ点及び前記注目点よりもＹ値が大きい点のいずれもが存在しない場合は、前記注目点を第３境界点として抽出する作業を繰り返し、抽出された全ての第３境界点をＸ値が小さい順に繋いだ線を、前記第１実効雨量及び前記第２実効雨量の相乗効果による実効雨量の最大値を表す第３境界線として前記判定用ＸＹ座標平面に設定し、
前記判定用ＸＹ座標平面において、前記測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出する超過降雨量算出手段を備えることを特徴とする斜面崩壊予測装置。
【請求項１４】
請求項１３に記載の斜面崩壊予測装置において、
前記判定用ＸＹ座標平面作成手段が、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
前記超過雨量算出手段が、前記複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出し、
該単位時間あたりの降雨量を、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成する二次元判定図作成手段を有することを特徴とする斜面崩壊予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、降雨に伴って発生する崖崩れや山崩れなどの斜面崩壊の発生を予測する斜面崩壊予測方法及び斜面崩壊予測装置に関する。
【背景技術】
【０００２】
崖崩れや山崩れ等の斜面崩壊の発生は人命や財産等に甚大な被害をもたらす。このため、斜面崩壊の的確な予測方法の確立は、斜面崩壊を回避するための対策を講じたり、警戒・避難の時期を判断したりする上で重要である。
【０００３】
斜面崩壊は、降雨時間が長い場合や降雨強度が大きい場合に発生しやすい。ただし、降雨時間が短い場合でも降雨強度が非常に大きければ、あるいは、降雨強度が小さい場合でも降雨時間が長く、それまでの降雨総量が多ければ斜面崩壊が発生することがある。このように、斜面崩壊の発生は、降雨総量と降雨強度の両方の影響を受ける。従って、降雨総量の影響を評価するための雨量指標（長期的雨量指標）と、降雨強度の影響を評価するための雨量指標（短期的雨量指標）を併用して、斜面崩壊の発生を予測することが効果的である。
【０００４】
そこで、長期的雨量指標と短期的雨量指標を用いて降雨状況を表したスネーク曲線と呼ばれるグラフを用いた予測方法が一般的に知られている。スネーク曲線は、長期的雨量指標と短期的雨量指標をそれぞれ横軸、縦軸とするグラフに、各時刻における長期的雨量指標と短期的雨量指標を表す点をプロットし、それらを線で繋いだもので、降雨時系列による点の軌跡を追跡することにより斜面崩壊の発生を予測する。
【０００５】
図１に、長期的雨量指標として半減期７２時間の実効雨量（７２ｈ実効雨量）を、短期的雨量指標として半減期１．５時間の実効雨量（１．５ｈ実効雨量）を用いたスネーク曲線の例を示す。実効雨量は、ある地点における先行降雨がその後の該地点の蓄積雨量に及ぼす影響の持続性を、半減期をパラメータとして表現した雨量指標である。半減期をＭ[ｈ]とすると、時刻ｔ[ｈ] における実効雨量Ｘ(Ｍ，ｔ)[ｍｍ]は次式（１）で定義される。
【数１】
2015232537.tif 000003
ここで、Ｘ(Ｍ，ｔ−１)[ｍｍ]は１時間前の実効雨量、Ｒ(ｔ)[ｍｍ]は時刻ｔ−１〜ｔの間の雨量（時間雨量）である。減少係数α[ｈ-1]は、半減期Ｍを用いた次式（２）で求められる。
【数２】
2015232537.tif 000004
【０００６】
実効雨量の式（１）では、ある時刻に降った雨量Ｒが半減期Ｍ[ｈ]経過するとＲ／２となり、更に半減期Ｍ[ｈ]が経過するとＲ／４になるというように、指数関数的に降雨の影響が減少する。従って、半減期Ｍが短い場合は、降雨の影響は時間と共に急激に減少することになり、実効雨量は実質的に降雨強度を表す指標となる。一方、半減期が長い場合は、降雨の影響がより長期間持続することになり、半減期を無限大に設定すれば実効雨量は積算雨量に等しくなる。図２に、ある年の各月における時間雨量を示すグラフ（ａ）、及び、この時間雨量について、半減期Ｍ[ｈ]を１．５、７２、２１６、７２０、及び３０００に設定したときの実効雨量を表すグラフ（ｂ）〜（ｆ）を示す。図２においては、例えば半減期１．５[ｈ]の実効雨量は地面の表層に貯留される水分量の指標を表し、半減期７２[ｈ]の実効雨量は地面から深層に貯留される水分量の指標を表す。
【０００７】
スネーク曲線上の境界線ＣＬは、崩壊発生限界の雨量線を示す。スネーク曲線において、境界線ＣＬより右上側に点がプロットされた時点で斜面崩壊が発生すると判断される。従って、スネーク曲線を用いた斜面崩壊の予測では、縦軸及び横軸に用いる雨量指標（半減期）の適切な選択と、適切な境界線ＣＬの設定が重要となる。
【０００８】
スネーク曲線の横軸及び縦軸に用いられる長期的雨量指標及び短期的雨量指標として、以下の例が挙げられる。
第一の例は、斜面崩壊の予測のために長年用いられてきた例であり、長期的雨量指標として積算雨量を、短期的雨量指標として時間雨量を利用したものである（非特許文献１）。
第二の例は、斜面内部の雨水貯留量を、図３に示すようなタンクモデルで評価した土壌雨量指数を長期的雨量指標とする例である（非特許文献２、非特許文献３）。1993年に牧原及び平沢によって提案され、近年、広く用いられている。
第三の例は、降雨の総量と強度の影響を半減期により調節した実効雨量を指標とする例であり、半減期を異ならせた様々な実効雨量について多くの検討が加えられてきた。特に、1993年に開催された総合土砂災害対策検討会（建設省）で提案された、半減期７２時間の実効雨量を長期的雨量指標とし、半減期１．５時間の実効雨量を短期的雨量指標とする方法が普及している（非特許文献４）。
【０００９】
一方、境界線については、例えば図４に示すように、対象地域における既往の斜面崩壊発生実績及び斜面崩壊非発生実績に基づいて、「見逃し」と「空振り」を共にできるだけ少なくすることができる直線をスネーク曲線上に引き、この直線を境界線ＣＬとする方法が永年用いられてきた（非特許文献１）。近年では、ニューラルネットワークを用いた方法により、「見逃し」と「空振り」を最小にする非直線型の境界線ＣＬを引く方法（非特許文献４）が提案され、広く使用されるようになっている。
【００１０】
ここで、「見逃し」とは、実際は斜面崩壊が発生している事例を非発生と判断してしまうことを指し、「空振り」とは、実際は斜面崩壊が発生していない事例を斜面崩壊が発生すると判断してしまうことを指す。「見逃し」が起きると、警戒・避難の勧告や指示がなされていない状況で斜面崩壊が発生するため、人的・物的被害が拡大する。一方、「空振り」が繰り返されると、対象地域の住民の斜面崩壊に対する意識が次第に低下するため、実際には斜面崩壊が発生する場合でも、勧告や指示を受けても警戒したり、避難行動をとったりしなくなるおそれがある。
従って、斜面崩壊を予測する上では「見逃し」と「空振り」をいかに少なくできるかが重要となる。
【先行技術文献】
【非特許文献】
【００１１】
【非特許文献１】寺田秀樹・中谷洋明（2001）：土砂災害警戒避難基準雨量の設定手法，国土技術政策総合研究所資料第5号，国土交通省，58pp.
【非特許文献２】牧原康隆・平沢正信（1993）：斜面崩壊危険度予測におけるタンクモデルの精度，気象庁研究時報，Vol.45，No.2，p.35-70
【非特許文献３】国土交通省気象庁ホームページ、"土壌雨量指数",[online],平成26年5月2日検索],インターネット<URL:http://www.jma.go.jp/jma/kishou/know/bosai/dojoshisu.html>
【非特許文献４】倉本和正・鉄賀博己・東寛和・荒川雅生・中山弘隆・古川浩平（2001）：RBF ネットワークを用いた非線形がけ崩れ発生限界雨量線の設定に関する研究，土木学会論文集，No.672/VI-50，p.117-132
【非特許文献５】小杉賢一朗・藤本将光・山川陽祐・正岡直也・糸数哲・水山高久・木下篤彦（2013）：山体基岩内部の地下水位変動を解析するための実効雨量に基づく関数モデル，砂防学会誌，Vol.66，No.4，p.21-32
【発明の概要】
【発明が解決しようとする課題】
【００１２】
力学的視点に立つと、斜面崩壊は、降雨によって地下水位が上昇し、ある臨界レベルを超過したときに発生すると考えることができる。従って、地下水位が上昇する様子や臨界レベルを適切に表現することができれば、斜面崩壊の発生予測における「見逃し」や「空振り」を少なくすることができる。
【００１３】
地下水位が上昇する様子（上昇する量やタイミング）は、雨の降り方、つまり、降雨時間や単位時間あたりの降雨量等によって決まるが、同じように雨が降っても地点によって地下水位が上昇する様子が異なる。これは、植生、地形、土壌、及び岩石の性質（保水性や透水性）が地点ごとに異なるため、雨が浸透する速さ、上流域からの水の集まりやすさ、下流域への水の通過しやすさがばらつくからである。
同様に、斜面崩壊が発生する地下水位の臨界レベルも、地点の植生、地形、土壌や岩石の性質（間隙の量や破壊に対する強度）に依存するため、地点ごとに異なる。
【００１４】
従って、地下水位の上昇の様子を適切に表現することができる雨量指標や、斜面崩壊の発生及び非発生の境界を矛盾無く説明できる境界線ＣＬは地点毎に異なるはずである。しかしながら、従来は、地点に関係なく、一律の雨量指標を使用して地下水位の上昇の様子を表現し、既往の斜面崩壊の発生・非発生の実績に基づいて境界線ＣＬを設定しており、斜面崩壊発生の予測精度に限界があった。
【００１５】
本発明が解決しようとする課題は、各対象地点において発生する斜面崩壊を正確に予測することができる斜面崩壊予測方法及び斜面崩壊予測装置を提供することである。
【課題を解決するための手段】
【００１６】
上記課題を解決するために成された本発明に係る斜面崩壊予測方法は、
対象地点における降雨イベントによる斜面崩壊の発生を予測する方法であって、
a) ある時刻における蓄積雨量が半減するまでの時間である半減期をパラメータとして、過去の降雨イベントが前記蓄積雨量に及ぼす影響の持続性を表した雨量指標である実効雨量の現時点までの時系列値を、前記対象地点の過去の降雨イベントの降雨量データを用い、所定の長半減期及び該長半減期より短い短半減期についてそれぞれ算出し、
b) 前記長半減期をパラメータとする第１実効雨量をＸ軸、前記短半減期をパラメータとする第２実効雨量をＹ軸とするＸＹ座標平面に、前記第１実効雨量の時系列値の最大値を通り前記Ｙ軸と平行な第１境界線及び前記第２実効雨量の時系列値の最大値を通り前記Ｘ軸と平行な第２境界線を設定して判定用ＸＹ座標平面を作成し、
c) 斜面崩壊発生の予測対象となる降雨イベントの降雨量データから、前記長半減期をパラメータとする第３実効雨量の時系列値及び前記短半減期をパラメータとする第４実効雨量の時系列値を、前記予測対象降雨イベントの進行と共に順次算出し、
d) 同じ時刻における前記第３実効雨量の時系列値及び前記第４実効雨量の時系列値をそれぞれＸ値、Ｙ値とする点を測定点としたとき、該測定点が、前記判定用ＸＹ座標平面において前記第１境界線を超える領域に位置するとき、又は前記第２境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定することを特徴とする。
【００１７】
ここで、降雨量データには、例えば気象庁や国土交通省等の公的機関が提供するデータ、民間の気象情報提供サービス企業が提供するデータを利用することができる。また、実効雨量の時系列値とは、複数の時刻における実効雨量の値を、時間の経過に従って並べたものをいう。時系列値を求める時刻の間隔は等間隔でも良く、不等間隔でも良い。通常、上述の公的機関や民間企業による降雨量データの提供が一定時間毎（例えば１時間毎）に行われることを考慮すると、実効雨量の時系列値の算出も一定時間毎に行うことが好ましいが、降雨イベントの最中及びその前後は、実効雨量の時系列値を算出する間隔を短くし、ある降雨イベントと、その次の降雨イベントの間の降雨量がゼロの期間、或いは降雨イベントの最中であっても降雨量が微量である期間は、時系列値を求める間隔を長くしても良い。
【００１８】
上記の斜面崩壊予測方法においては、前記対象地点の過去の降雨イベントの降雨量データから求められた、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットし、
前記複数の点のうちある時刻の点を注目点としたとき、該注目点とＸ値が同じ点の中に前記注目点よりもＹ値が大きい点が存在しない場合、及び前記注目点よりもＸ値が大きい点の中に前記注目点とＹ値が同じ点及び前記注目点よりもＹ値が大きい点のいずれもが存在しない場合は、前記注目点を第３境界点として抽出する作業を繰り返し、
抽出された全ての第３境界点をＸ値が小さい順に繋いだ線を、前記第１実効雨量及び前記第２実効雨量の相乗効果による実効雨量の最大値を表す第３境界線として前記ＸＹ座標平面に設定し、
前記測定点が、前記第３境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定することが好ましい。
【００１９】
上記の方法において、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットしたものは、「スネーク曲線」に相当する。スネーク曲線は、斜面崩壊の発生の予測に一般的に用いられているグラフであり、上記方法では、スネーク曲線を構成する複数の点を使って第３境界線を設定した。このため、より適切な境界線に基づいて斜面崩壊の発生の可能性を予測することができる。
【００２０】
斜面崩壊の発生を正確に予測する上では、長半減期及び短半減期の選択が重要である。これら２種類の半減期の選択を誤ると、判定用ＸＹ座標平面に適切な境界線を設定することができず、また、測定点が降雨イベントの性質や対象地点の保水性・透水性を反映したものとならない。この結果、判定用ＸＹ座標平面における測定点の位置に基づいて斜面崩壊の発生の可能性を予測しても、その予測が「空振り」に終わったり、「見逃し」になったりするおそれがある。
【００２１】
そこで、本発明の斜面崩壊予測方法においては、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、各時刻における判別結果の集計から、斜面崩壊の発生を予測すると良い。
この方法によれば、斜面崩壊が発生しない状況であるにも関わらず、斜面崩壊が発生する可能性が高いと予測してしまう、「空振り」を減らすことができる。
【００２２】
また、本発明の斜面崩壊予測方法においては、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、少なくとも１個の判定用ＸＹ座標平面において測定点が第２〜第５領域のいずれかに位置するとき、斜面崩壊の発生の可能性が高いと予測することが好ましい。
この方法によれば、斜面崩壊の発生の可能性を見落としてしまう、「見逃し」を回避することができる。
【００２３】
さらに、本発明の斜面崩壊予測方法においては、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、
その結果を表すシンボルを、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成するようにしても良い。
本発明においては、長半減期及び短半減期の設定が重要であるが、前記二次元判定図を作成することにより、長半減期及び短半減期の組み合わせの違いによる判定結果の変化を視覚的に認識することができる。従って、二次元判定図を参考にして、適切な長半減期及び短半減期を設定することができる。
【００２４】
また、本発明の斜面崩壊予測方法は、上述の判定用ＸＹ座標平面において、測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出することを特徴とする。
この方法によれば、斜面崩壊の発生の可能性を定量的に予測することができる。
【００２５】
この場合、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出し、
該単位時間あたりの降雨量を、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成するようにしても良い。
この方法によれば、長半減期及び短半減期の組み合わせの違いによる斜面崩壊をもたらす降雨量の変化を視覚的に認識することができる。従って、二次元判定図を参考にして、適切な長半減期及び短半減期を設定することができる。
【００２６】
また、本発明に係る斜面崩壊発生予測装置は、対象地点における降雨イベントによる斜面崩壊の発生を予測する装置であって、
a) ある時刻における蓄積雨量が半減するまでの時間である半減期をパラメータとして、過去の降雨イベントが前記蓄積雨量に及ぼす影響の持続性を表した雨量指標である実効雨量の現時点までの時系列値を、前記対象地点の過去の降雨イベントの降雨量データを用い、所定の長半減期及び該長半減期より短い短半減期についてそれぞれ算出する、既往実効雨量算出手段と、
b) 前記長半減期をパラメータとする第１実効雨量をＸ軸、前記短半減期をパラメータとする第２実効雨量をＹ軸とするＸＹ座標平面に、前記第１実効雨量の時系列値の最大値を通り前記Ｙ軸と平行な第１境界線及び前記第２実効雨量の時系列値の最大値を通り前記Ｘ軸と平行な第２境界線を設定して判定用ＸＹ座標平面を作成する、判定用ＸＹ座標平面作成手段と、
c) 斜面崩壊発生の予測対象となる降雨イベントの降雨量データから、前記長半減期をパラメータとする第３実効雨量の時系列値及び前記短半減期をパラメータとする第４実効雨量の時系列値を、前記予測対象降雨イベントの進行と共に順次算出する、予測用実効雨量算出手段と、
d) 同じ時刻における前記第３実効雨量の時系列値及び前記第４実効雨量の時系列値をそれぞれＸ値、Ｙ値とする点を測定点としたとき、該測定点が、前記判定用ＸＹ座標平面において前記第１境界線を超える領域に位置するとき、又は前記第２境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定する判定手段と
を備えることを特徴とする。
【００２７】
上記の斜面崩壊予測装置においては、
前記判定用ＸＹ座標平面作成手段が、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットし、前記複数の点のうちある時刻の点を注目点としたとき、該注目点とＸ値が同じ点の中に前記注目点よりもＹ値が大きい点が存在しない場合、及び前記注目点よりもＸ値が大きい点の中に前記注目点とＹ値が同じ点及び前記注目点よりもＹ値が大きい点のいずれもが存在しない場合は、前記注目点を第３境界点として抽出する作業を繰り返し、抽出された全ての第３境界点をＸ値が小さい順に繋いだ線を、前記第１実効雨量及び前記第２実効雨量の相乗効果による実効雨量の最大値を表す第３境界線として前記判定用ＸＹ座標平面に設定し、
前記判定手段が、前記判定用ＸＹ座標平面において前記測定点が前記第３境界線を超える領域に位置するときに斜面崩壊が発生する可能性があると判定することが好ましい。
【００２８】
また、前記既往実効雨量算出手段が、長さが異なる３種類以上の半減期について前記実効雨量の時系列値を算出し、
前記判定用ＸＹ座標平面作成手段が、前記３種類以上の半減期から２種類の半減期を選択することにより形成される、複数の、長半減期と短半減期の組み合わせについて、それぞれ前記第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
前記判定手段が、前記複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、各時刻における判別結果の集計から、斜面崩壊の発生を予測するようにしても良い。
【００２９】
さらに、前記既往実効雨量算出手段が、長さが異なる３種類以上の半減期について前記実効雨量の時系列値を算出し、
前記判定用ＸＹ座標平面作成手段が、前記３種類以上の半減期から２種類の半減期を選択することにより形成される、複数の、長半減期と短半減期の組み合わせについて、それぞれ前記第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
前記判定手段が、前記複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、少なくとも１個の判定用ＸＹ座標平面において測定点が第２〜第５領域のいずれかに位置するとき、斜面崩壊の発生の可能性が高いと予測することも良い構成である。
【００３０】
さらに、上記の斜面崩壊予測装置において、
前記判定用ＸＹ座標平面作成手段が、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
これら複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が、前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第１領域、前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた第２領域、前記第１境界線を超過し且つ前記第２境界線よりも前記Ｘ軸側の第３領域、前記第２境界線を超過し且つ前記第１境界線よりも前記Ｙ軸側の第４領域、前記第１境界線及び前記第２境界線の両方を超過する第５領域のいずれの領域に位置するかを判別し、その結果を表すシンボルを、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成する二次元判定図作成手段を備えるようにしても良い。
【００３１】
また、前記判定用ＸＹ座標平面作成手段が、前記対象地点の過去の降雨イベントの降雨量データから求められた、同じ時刻の前記第１実効雨量の時系列値及び前記第２実効雨量の時系列値をそれぞれＸ値及びＹ値とする複数の点を前記判定用ＸＹ座標平面にプロットし、前記複数の点のうちある時刻の点を注目点としたとき、該注目点とＸ値が同じ点の中に前記注目点よりもＹ値が大きい点が存在しない場合、及び前記注目点よりもＸ値が大きい点の中に前記注目点とＹ値が同じ点及び前記注目点よりもＹ値が大きい点のいずれもが存在しない場合は、前記注目点を第３境界点として抽出する作業を繰り返し、抽出された全ての第３境界点をＸ値が小さい順に繋いだ線を、前記第１実効雨量及び前記第２実効雨量の相乗効果による実効雨量の最大値を表す第３境界線として前記判定用ＸＹ座標平面に設定し、
前記判定用ＸＹ座標平面において、前記測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出する超過降雨量算出手段を備えるようにしても良い。
前記第１〜第３境界線を超えることになる単位時間あたりの降雨量の大きさから、斜面崩壊の発生の可能性を定量的に予測することができる。
【００３２】
この場合、前記判定用ＸＹ座標平面作成手段が、複数の、長半減期及び短半減期の組み合わせについて、それぞれ第１〜第３境界線を設定した判定用ＸＹ座標平面を求め、
前記超過雨量算出手段が、前記複数の判定用ＸＹ座標平面の全てについて、対応する第３実効雨量の時系列値及び第４実効雨量の時系列値からなる測定点のうち直近の測定点が前記Ｘ軸と前記Ｙ軸と前記第１境界線と前記第２境界線と前記第３境界線とに囲まれた領域に位置するとき、次の測定点が前記第１〜第３境界線を超えることになる単位時間あたりの降雨量を算出し、
該単位時間あたりの降雨量を、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成する二次元判定図作成手段を設けると、長半減期及び短半減期の組み合わせの違いによる斜面崩壊をもたらす降雨量の変化を視覚的に認識することができ、上記二次元判定図を参考にして適切な長半減期及び短半減期を設定することができる。
【００３３】
なお、本発明において、斜面崩壊が発生する可能性の有無の判定に用いる測定点は、通常は直近の測定点であるが、これに限定されるものではない。例えば、第３実効雨量及び第４実効雨量の時系列値の取得間隔（測定間隔）が短い場合には、直近の測定点よりも１ないし複数回分だけ過去に遡った測定点を用いても良い。また、判定に用いる測定点は１個に限らず、複数個でも良い。
【発明の効果】
【００３４】
本発明に係る斜面崩壊予測方法及び装置では、対象地域の過去の降雨イベントの降雨量データに基づいて境界線が設定された判定用ＸＹ座標平面を作成し、斜面崩壊の発生を予測する対象となる降雨イベントの進行と共に得られる降雨量データから得られる測定点と前記判定用ＸＹ座標平面に設定された境界線との位置関係に基づいて、斜面崩壊の発生の可能性を判定する。従って、斜面崩壊の発生の予測に、対象地点の保水性や透水性、予測対象となる降雨イベントの性質等を反映させることができ、「見逃し」や「空振り」を回避した、正確な斜面崩壊予測を行うことができる。
【図面の簡単な説明】
【００３５】
【図１】斜面崩壊の発生予測に用いるスネーク曲線図の例を示す図。
【図２】時間雨量（ａ）と各種半減期の実効雨量（ｂ）〜（ｆ）を示すグラフ。
【図３】タンクモデルを用いた土壌雨量指数の説明図。
【図４】スネーク曲線図における境界線ＣＬの最適化を説明する図。
【図５】スネーク曲線図において、各時刻の点をその時刻に観測された地下水位の大きさに応じて塗り分けた例を示す図。
【図６】スネーク曲線における境界線ＣＬの定義を説明するための図。
【図７】判定用グラフを用いた予測対象となる降雨イベントのスネーク曲線の判定結果を示す図。
【図８】長半減期及び短半減期の組み合わせが異なる複数のスネーク曲線から求めた判定結果のシンボルを示す二次元判定図の一例を示す図。
【図９】本発明に係る斜面崩壊予測装置の概略構成図。
【図１０】本発明に係る斜面崩壊予測方法の実施例１に用いた庄原災害時の雨量データ。
【図１１】庄原災害時の二次元判定図による既往最大値超過の判定結果。
【図１２】本発明に係る斜面崩壊予測方法の実施例２を示すものであり、梨子沢災害の発生場所及びその付近の観測所を示す地図（国土地理院地図に加筆）。
【図１３】３箇所の観測所における観測データを示すグラフ。
【図１４】既往降雨イベント及び梨子沢災害の原因となった今回の降雨イベントの時間雨量及び積算雨量を示すグラフ。
【図１５】既往降雨イベント及び今回の降雨イベントの、横軸を土壌雨量指数、縦軸を時間雨量としたスネーク曲線と土壌雨量指数既往最大値と時間雨量既往最大値とＡＮ境界線によって得られる既往最大値を示す図。
【図１６】既往降雨イベント及び今回の降雨イベントの、横軸を半減期51hの実効雨量、縦軸を半減期0.9hの実効雨量としたスネーク曲線と半減期51h実効雨量既往最大値と半減期0.9h実効雨量既往最大値とＡＮ境界線によって得られる既往最大値を示す図。
【図１７】長半減期及び短半減期の組み合わせが異なる複数のスネーク曲線から求めた、測定点がＡＮ境界線を超えるまでに必要な１時間の降雨量を示す二次元判定図。
【図１８】現時点までの降雨量データに今後１時間の予測降雨量を加えて求めたスネーク曲線（ａ）及び、現時点までの降雨量データにＡＮ境界線を超過する最小の１時間降雨量を加えて求めたスネーク曲線（ｂ）。
【図１９】斜面崩壊予測システムの概念図。
【発明を実施するための形態】
【００３６】
斜面崩壊が発生する地下水位の臨界レベルを自然斜面において正確に予測するのは困難であることから、地下水位が既往最大値を超えて上昇しているときには、いつ斜面崩壊が発生してもおかしくない状態にあると考えるべきである。本発明では、このような考えに基づき、地下水位の既往最大値を基準として臨界レベルを設定することとした。これにより、斜面崩壊発生予測における「見逃し」を極力なくすことができると考えられる。
【００３７】
降雨中の地下水位の予測は、本発明者が2013年に提案した、基岩内部の地下水位の変動を解析するための関数モデルの考え方を活用した（非特許文献５）。
本発明者は、半減期の長い実効雨量と半減期の短い実効雨量をそれぞれ横軸、縦軸とするスネーク曲線において、各時刻の点をその時刻に観測された地下水位の大きさに応じて塗り分けたとき、横軸及び縦軸の実効雨量の半減期を最適化することにより、横軸及び縦軸の実効雨量の大きさに応じて点の色がきれいに塗り分けられることを、非特許文献５に示した。つまり、半減期の異なる２種類の実効雨量からスネーク曲線を作成したとき、２種類の半減期の組み合わせが最適であれば、地下水位上昇の様子を再現できる。
【００３８】
図５に、半減期７２時間の実効雨量を横軸に、半減期１．５時間の実効雨量を縦軸にしたスネーク曲線図（ａ）、及び半減期１０６１時間の実効雨量を横軸に、半減期７７時間の実効雨量を縦軸にしたスネーク曲線図（ｂ）の例を示す。同図では、各時刻の点の色をその時刻に観測された地下水位が低いものから順に塗り分けている。なお、ここでは、図５を白と黒の濃淡図で表しているが、実際は、茶色、赤色、黄色、黄緑色、青色に塗り分けている。同図（ａ）では、広い範囲に亘って異なる色の点が混在しているのに対して、同図（ｂ）では、複数の点が色ごとにまとまって存在しており、実効雨量と地下水位の間に相関関係があると推定される。
【００３９】
以上より、図５（ｂ）のスネーク曲線図上では、実効雨量が既往最大値を超過しているときは、地下水位も既往最大値を超過していると考えられ、このスネーク曲線図を用いて斜面崩壊の発生を予測する場合は、実効雨量の既往最大値を境界線ＣＬに設定することが、最も理に適っている。実効雨量の既往最大値は、本発明における、第１実効雨量の時系列値の最大値（第１境界線）及び第２実効雨量の時系列値の最大値（第２境界線）、並びに第１実効雨量と第２実効雨量の相乗効果による実効雨量の最大値（第３境界線）に相当する。第１〜第３境界線については後述する。
【００４０】
しかし、スネーク曲線図の縦軸及び横軸にとる実効雨量の半減期をどのような値に設定したときに、図５（ｂ）に示すようなスネーク曲線図を得ることができるか、すなわち、地下水位上昇の様子を最も良く表現できる半減期が何であるかは地点ごとに異なる。このため、適切な半減期の組み合わせを全ての対象地点において正確に予測することは困難である。
【００４１】
これに対して、横軸に設定する実効雨量の半減期（以下「長半減期」ともいう。）と縦軸に設定する実効雨量の半減期（以下「短半減期」ともいう）を種々の長さに変化させ、あらゆる半減期の組合せについて、地下水位の既往最大値の超過が起きているかどうかを判定すると、いかなる地下水位上昇パターンを示す地点であっても地下水位の既往最大値の超過を検出できる。地下水位の既往最大値では、少なくとも斜面崩壊は発生しておらず、それを超えたときに斜面崩壊が発生するか否かは不明であるため、上記のように、あらゆる半減期の組合せについて、地下水位の既往最大値の超過が起きているかどうかを判定することは、斜面崩壊の発生の予測の「空振り」や「見逃し」の回避につながる。
【００４２】
ただし、上述したように、長半減期と短半減期の最適な組み合わせを選択することができれば、１個の組み合わせについて、地下水位の既往最大値の超過が起きているかどうかを判定するだけで、「空振り」や「見逃し」を回避した斜面崩壊の予測を行うことができる。
また、長半減期と短半減期の全ての組み合わせについて、地下水位の既往最大値の超過が起きているかどうかを判定することは、現実的に不可能である。「空振り」や「見逃し」を少なくするためには、判定に用いる長半減期と短半減期の組み合わせの数が多い方が有効であるが、少なくとも２個の組み合わせについて、上述した判定を行うことにより、１個の組み合わせについて判定を行う場合よりも斜面崩壊の予測精度を上げることができる。
【００４３】
実効雨量の既往最大値、つまり地下水位の既往最大値は、あくまでも過去の降雨時における地下水位の最大値であり、斜面崩壊が発生するときの地下水位を表しているわけではない。このため、雨量の観測期間が短い間は、実効雨量が既往最大値を超過したことに基づいて斜面崩壊の発生を予測しても、「空振り」となる事態が頻繁に起こる可能性が高い。しかし、「空振り」が起きる度に実効雨量の既往最大値が更新され、境界線ＣＬが引き上げられるため、雨量の観測を継続することによって次第に「空振り」は減少する。
【００４４】
また、本発明で使用する雨量指標（2種類の実効雨量）は、従来技術と異なり、「地下水位上昇の様子を最も良く再現できる雨量指標」という明確な物理的意味を持っている（図５参照）。このため、実際の降雨時における地下水位の変動を計測することができれば、適切な半減期を特定できる。つまり、本発明では、対象地点での地下水位の観測結果を斜面崩壊の発生予測に直接反映できるという利点を持つ。このため、地下水位の観測データの蓄積により「空振り」を効率的に減らすことができる。
【００４５】
さらに、本発明では、実効雨量の既往最大値を超過したか否かを判断し、その判断結果に基づき、斜面崩壊の発生を予測する。このような判断は、多くの住民が常日頃から行っている「この程度の雨は、これまでも起きたことがあるから大丈夫だ」「今回のような大雨は、初めて起きるから危険だ」といった、自らの過去の経験に基づく斜面崩壊の危険性の判断に通じるものであるため、たとえ「空振り」が起きたとしても、従来の予測技術に比べて住民に許容され易く、警戒・避難に対する意識低下を招きにくい。
【００４６】
実効雨量が既往最大値を超過したか否かの判定は、例えば図６に示したスネーク曲線図を用いて行う。図６のスネーク曲線図は、本発明の判定用ＸＹ座標平面に相当する。
まず、対象地点において降雨イベントが発生する直前までの長期時間雨量データを用いて、半減期がそれぞれＭ１，Ｍ２（ただしＭ１≧Ｍ２とする）の2種類の実効雨量Ｘ(Ｍ１，ｔ)，Ｘ(Ｍ２，ｔ)を算定し、実効雨量Ｘ(Ｍ１，ｔ)を横軸に、実効雨量Ｘ(Ｍ２，ｔ)を縦軸にとってスネーク曲線図を描く。このスネーク曲線図で「既往最大値超過」は、以下の４つのケースに分類される。
【００４７】
まず、スネーク曲線図上に、実効雨量Ｘ(Ｍ１，ｔ) ＝Ｘｍａｘ(Ｍ１)（Ｘ(Ｍ１，ｔ)の既往最大値）と実効雨量Ｘ(Ｍ２，ｔ) ＝Ｘｍａｘ(Ｍ２)（Ｘ(Ｍ２，ｔ)の既往最大値）の２本の直線を引く。これら２本の直線が本発明の第１境界線及び第２境界線に相当する。これら２本の直線によって作られる４つの領域のうち領域Ｂは、実効雨量Ｘ(Ｍ１，ｔ)が単独で既往最大値を超過しているケースに相当する。領域Ｃは、実効雨量Ｘ(Ｍ２，ｔ)が単独で既往最大値を超過しているケースを表している。領域Ｄは、実効雨量Ｘ(Ｍ１，ｔ)と実効雨量Ｘ(Ｍ２，ｔ)の両方が、それぞれの既往最大値を超過しているケースである。
【００４８】
これに対して領域（Ｎ＋Ａ）は、実効雨量Ｘ(Ｍ１，ｔ)と実効雨量Ｘ(Ｍ２，ｔ)がいずれも既往最大値を超過していないケースであり、そのうち領域Ｎは実効雨量Ｘ(Ｍ１，ｔ)と実効雨量Ｘ(Ｍ２，ｔ)の組み合わせでみても既往範囲内にあるケースを、領域Ａは実効雨量Ｘ(Ｍ１，ｔ)と実効雨量Ｘ(Ｍ２，ｔ)の組合せでみると既往最大値を超過しているケースを表している。
【００４９】
領域Ａと領域Ｎを分離する境界線（以下「ＡＮ境界線」と呼ぶ）を定めるために、ＡＮ境界線上にある点（以下、ＡＮ境界点と呼ぶ）を次の方法で抽出する。すなわち、ある時刻におけるスネーク曲線上の点ｐに対し、点ｐとＸ(Ｍ１，ｔ)が等しい点の中に、点ｐよりもＸ(Ｍ２，ｔ)が大きいものが無く、且つ、点ｐよりもＸ(Ｍ１，ｔ)が大きい点の中に、点ｐとＸ(Ｍ２，ｔ)が等しい点、及び点ｐよりもＸ(Ｍ２，ｔ)が大きい点が無い場合に、点ｐをＡＮ境界点と判断する。（※請求項の表現に合わせました。）従って、最初に選んだ点ｐが本発明の注目点に相当し、ＡＮ境界点が第３境界点に相当する。
その上で、抽出された全てのＡＮ境界点を実効雨量Ｘ(Ｍ１，ｔ)が小さいものから順に並び換え、図６に示すように直線で結ぶ。これにより得られる階段状の線をＡＮ境界線とする。ＡＮ境界線が本発明の第３境界線に相当する。なお、並び換えた後の最初のＡＮ境界点の実効雨量Ｘ(Ｍ２，ｔ)は既往最大値Ｘｍａｘ(Ｍ２)に一致し、最後のＡＮ境界点の実効雨量Ｘ(Ｍ１，ｔ)は既往最大値Ｘｍａｘ(Ｍ１)に一致する。
【００５０】
以上の手順で領域Ｎおよび領域Ａ〜Ｄを定めた上で、当該イベント中のスネーク曲線図の軌跡を重ねることによって、各時刻の点がいずれの領域に属するかを調べる（図７）。その結果から、各時刻について、既往最大値超過の有無の判定、ならびに超過が起きている場合には、領域Ａ〜Ｄのいずれに位置するケースに該当するかの判定を行う。
【００５１】
以上説明した解析を、半減期Ｍ１及びＭ２を種々の値に設定して実施する。そして各時刻に得られた結果を図８の二次元判定図に表示する。二次元判定図は、横軸を半減期Ｍ１、縦軸を半減期Ｍ２とし、これら半減期Ｍ１，Ｍ２の組み合わせによるスネーク曲線図において、当該時刻の点が領域Ａ〜Ｄもしくは領域Ｎのいずれに存在しているかを、それぞれの領域に対応するシンボルで表した図である。従って、この二次元判定図は、請求項５及び１２に係る二次元判定図に相当する。なお、二次元判定図の対角線上の点は半減期Ｍ１ = 半減期Ｍ２のスネーク曲線図の結果、つまり、１種類の実効雨量を単独で評価した場合の、既往最大値超過の有無の判定結果を示している。この場合は、領域ＮもしくはＤのいずれかであり、領域Ａ〜Ｃにはならない。
【００５２】
斜面崩壊の形態には様々なものがある。地表面近くの厚さ0.5〜2ｍ程度の土層が崩れる「表層崩壊」と、山体を構成する岩石の層（基岩層）も含めて厚さ数十メートルに及ぶ大規模な崩れが生じる「深層崩壊」が両極端なものであるが、それらの中間的な形態を示す斜面崩壊も多く見られる。本発明の方法では、「あらゆる半減期を持つ２種類の実効雨量の組み合わせについて、既往最大値超過を判定する」ため、斜面崩壊の形態別に着目すべき雨量指標に違いがあることを定量的に示すことができると期待される。従って、その結果から降雨量データを基にして起こり得る斜面崩壊の形態を精度良く予測でき、よりきめの細かい警戒・避難体制を整えることができる。
【００５３】
また、本発明の方法では、斜面崩壊の発生予測だけでなく、地震や樹木の風倒被害、人為的な土地改変等に起因する斜面の変化が斜面崩壊に与える影響を評価できる。例えば、大きな地震を経験した地域では、斜面の強度が低下して過去に崩壊を引き起こさなかった規模の降雨でも崩壊が発生することが懸念され、境界線ＣＬの引き下げが検討されることが多々ある。本発明の手法に基づき、スネーク曲線図の半減期を連続的に変化させて降雨の長期的・短期的影響を漏らさず評価することによって、実際に既往最大値超過が発生しない状況下で斜面崩壊が起こるのか否かについて、客観的かつ定量的な検討を加えることが可能になる。
【００５４】
以上の本発明に係る斜面崩壊予測方法を実行するための斜面崩壊予測装置の概略図を図９に示す。斜面崩壊予測装置１０は、データ受信部１１、既往実効雨量算出部１２、判定用ＸＹ座標平面作成部１３、予測用実効雨量算出部１４、判定部１５、二次元判定図作成部１７、及びこれらを制御する制御部１６を備える。制御部１６には、プリンタやディスプレイなどの出力装置３０が接続されている。
【００５５】
制御部１６からの指令を受けてデータ受信部１１が降雨量データ提供装置２０から降雨量データを取得すると、制御部１６は該降雨量データを既往実効雨量算出部１２に送る。
降雨量データ提供装置２０は、例えば国土交通省の防災情報提供センターや観測所、気象庁のアメダス観測所等の過去の降雨量データが格納されたデータベースを備えた機関を指す。データ受信部１１は、降雨量データ提供装置２０から降雨量データを直接受信しても良く、インターネットを介してデータベースにアクセスすることにより降雨量データを取得しても良い。制御部１６を介してデータ受信部１１から降雨量データを受け取った既往実効雨量算出部１２は、予め設定された複数の半減期について、過去の降雨量データ及び新たに取得した降雨量データを用いて実効雨量の時系列値を算出する。
【００５６】
判定用ＸＹ座標平面作成部１３は、これら実効雨量の時系列値から複数の任意の２種類の半減期の組み合わせを抽出し、各組の２個の半減期のうち時間が長い方を長半減期、時間が短い方を短半減期としてスネーク曲線を求める。そして、スネーク曲線に基づいて第１〜第３の境界線を設定した判定用ＸＹ座標平面図を作成し、保存する。
判定用ＸＹ座標平面図は、新たな降雨量データを受信する毎に作成され、保存される。
【００５７】
一方、予測用実効雨量算出部１４は、新たに取得した降雨量データを用いて、予め設定された複数の半減期の実効雨量を算出する。そして、２種類の半減期の組み合わせについて求めた実効雨量の時系列値を示す点（測定点）を、対応する半減期の組み合わせの判定用ＸＹ座標平面図にプロットする。判定部１５は測定点の位置と第１〜第３境界線との位置関係から斜面崩壊の発生の可能性を予測する。判定部１５は、判定用ＸＹ座標平面作成部１３が作成した複数の判定用ＸＹ座標平面図について、斜面崩壊の発生の可能性を予測し、その結果を記憶する。また、測定点が第１〜第３境界線を越えていない場合は、境界線を越えることになる次の一定時間（例えば１時間）の（もしくは，「単位時間あたりの」）降雨量を予測する。従って、ここでは判定部１５が超過降雨量算出手段としても機能する。
【００５８】
二次元判定図作成部１７は、ＸＹ座標平面作成部１３が記憶する予測結果を表すシンボルを、横軸を長半減期、縦軸を短半減期とする二次元座標平面にプロットした二次元判定図を作成する。予測結果を表すシンボルとしては、例えば、斜面崩壊の可能性の高さに応じて大きさの異なる図形、色が異なる図形などを用いることができる。
次に、具体的な降雨量データを用いて解析を行った事例を説明する。
【実施例１】
【００５９】
（１）対象災害
広島県庄原市川北町では、2010年7月16日に発生した豪雨により斜面崩壊が集中的に発生し、土石流を伴って、死者1名、重傷者1名、家屋損壊38棟の被害をもたらした。崩壊した斜面は火山岩、堆積岩、火山灰由来のクロボク土層が混在する複雑な地質を呈し、スギ・ヒノキ人工林、広葉樹二次林を植生とするが、伐採跡地も含まれていた。斜面崩壊の大半は、幅10〜15ｍ、長さ20〜60ｍ、深さ0.5〜1.5ｍの規模で、表層崩壊に分類された。以下では、この災害を庄原災害と呼ぶ。
【００６０】
図１０に、庄原市の最寄りの広島県大戸観測所で観測された雨量データと、大雨警報発表時刻、斜面崩壊開始時刻、及び土砂災害警戒情報発表時刻を示す。図１０より、庄原災害を引き起こした降雨イベントは、降雨強度が大きい集中豪雨であると考えられるが、これに先行する5日間に累積雨量260〜270ｍｍ程の大雨があったことも災害発生に大きな影響を与えていたと思われる。図１０に示すように、広島地方気象台より庄原市に大雨警報が発表されたのは7月16日の16時39分であり、土砂災害警戒情報が発表されたのは18時10分であった。しかし、16時55分には斜面崩壊・土石流発生の目撃情報が寄せられており、17時前から多くの斜面崩壊が発生し始めていたと考えられる。
【００６１】
（２）解析に用いた雨量データ
図１０に示した広島県大戸観測所のデータは長期間に渡って遡って取得することができない。そこで、被災地から約7ｋｍ離れた地点にある気象庁アメダス庄原観測所（北緯34度51.6分，東経133度1.4分，標高300ｍ）の雨量データを用いて実効雨量既往最大値の計算を行った。庄原観測所の雨量データは、庄原災害から約34年7ヶ月前の1976年1月1日1時まで遡ることができる。実効雨量の計算は、初期値を０とした上で、データの先頭時刻から連続して行った。ただし既往最大値の評価は、災害が発生した年の1月1日午前0時から30年間遡った期間を対象とした。これは、実効雨量の計算結果に初期値が与える影響を取り除くことを目的とした処置である。
【００６２】
半減期Ｍ１，Ｍ２（ただしＭ１≧Ｍ２）を0.1ｈ〜3000ｈの範囲で、その対数値の変化量が一定になるようにして、計21通りに変化させた。ここで半減期の設定範囲（0.1ｈ〜3000ｈ）は，実効雨量が、降雨強度や積算雨量といった先行降雨の影響の持続性に関する両極端な指標も十分にカバーすることを意図して定めている。特に最大値（3000ｈ）に関しては、非特許文献５において、本発明者らが地下水位データの再現に用いた最大値（1400ｈ）を考慮し、その約2倍の値に設定している。半減期の設定間隔については、連続した傾向が得られることを確認した上で定めている。
【００６３】
（３）結果
庄原災害に関する二次元判定図を図１１に示す。降雨開始から1時間が経過した7月16日16時00分（積算雨量38ｍｍ；図１０参照）には、既に48.6ｈ≦Ｍ１≦ 1070ｈ、0.100ｈ≦Ｍ２≦ 0.469ｈの範囲で領域Ａが出現し、既往最大値超過が起きていたことがわかる（図１１ａ）。単独の実効雨量に着目した場合には、この時刻にいずれの実効雨量も既往最大値を超過していなかった。
72ｍｍの時間雨量が降った直後の 17時00分には、多くの半減期Ｍ１，Ｍ２の組合せで既往最大値超過が発生した（図１１ｂ）。ただし、この時刻には、136ｈ≦Ｍ２≦Ｍ１の範囲の組合せでは既往最大値超過は発生していなかった。さらに、63ｍｍの時間雨量が降った18時00分には、既往最大値超過は半減期Ｍ１、Ｍ２が大きな範囲にも拡大し（図１１ｃ）、降雨終了時の19時00分になると半減期Ｍ１が小さな範囲で縮小した（図１１ｄ）。
【００６４】
（４）考察
庄原災害時の降り始め1時間（15〜16時）の雨量は38ｍｍ（図１０）であり、既往最大値64ｍｍ（1998年8月22日に記録）と比べると、６割程度のあまり大きな値ではなかった。しかしながら、先行する5日間に累積271ｍｍの降雨が降っていた。二次元判定図（図１１）では、先行降雨の影響により48.6ｈ≦Ｍ１≦ 1070ｈの横軸実効雨量が大きく増加したことと、15〜16時の38ｍｍの時間雨量によって0.100ｈ≦Ｍ２≦ 0.469ｈの縦軸実効雨量が増加したことの組み合わせによって、既往最大値超過が16時00分の段階で発生したことを示すことができている。
【００６５】
一方、半減期72時間の実効雨量を長期的雨量指標、半減期1.5時間の実効雨量を短期的雨量指標に用いる従来技術の方法では、既往最大値超過が初めて判定されたのは17時00分であり（図示せず）、斜面崩壊の発生開始よりも遅れた。また従来技術のうち、土壌雨量指数を長期的雨量指標、時間雨量を短期的雨量指標に用いる方法で得られた結果を図１１の二次元判定図中のサブウィンドウ（各パネル中の実線の四角形）内に示した。この場合も、既往最大値超過が初めて判定されたのは17:00であり（図１１ｂのサブウィンドウ内のシンボル「●」参照）、斜面崩壊の発生開始よりも遅れることがわかる。
このことから、実際の斜面崩壊の発生開始よりも既往最大値超過の判定が遅れていた従来の判定方法と異なり、本発明に係る予測方法によれば、庄原災害発生の１時間前から危険度が増加している状況を予測でき、従来技術に比べて「見逃し」を少なく抑えることができる。
【実施例２】
【００６６】
（１）対象災害
長野県南木曽町では、台風８号の影響を受けた豪雨によって、2014年7月9日の17:40頃に梨子沢（なしざわ）と呼ばれる木曽川の支流で土石流が発生した。以下、この災害を「梨子沢災害」と呼ぶ。梨子沢災害が起こった当日は昼過ぎまで晴れており、雨の降り始めは土石流発生の約２時間前であった。このように、土石流は不意を突く形で発生しており、南木曽町に避難勧告が出されたのは、土石流発生から１０分後の17:50であった。また、土砂災害警戒情報が発表されたのは、さらに25分を経過した18:15であった。
【００６７】
図１２に、土石流が発生した梨子沢流域及びその付近の降雨観測所（気象庁の南木曽観測所、国土交通省の三留野観測所及び蘭観測所）の位置を示す。図１２から分かるように、南木曽観測所が土石流発生流域に最寄りの観測所である。
これらの３つの観測所における梨子沢災害当日の10分間雨量を図１３に示す。いずれの観測所でも7月9日の3:00以降は降雨が一旦休止していたが、15:30過ぎに再開し、短時間のうちに降雨強度が増したことがわかる。最大10分間雨量は蘭観測所で最も多く23ｍｍであり、次に多いのが南木曽観測所の17ｍｍであった。
【００６８】
15:30から土石流が発生する17:40までの積算雨量は、土石流発生流域に最寄りの南木曽観測所で84ｍｍに達したが、三留野観測所及び蘭観測所に比べれば少なかった。また、土石流発生までの１時間（16:40〜17:40）の雨量は南木曽観測所で70ｍｍに達したが、この値も蘭観測所の97ｍｍに比べると少なかった。このように、土石流発生流域に最寄りの南木曽観測所の雨量は、他の２つの観測所に比べて若干少ない傾向にあった。また気象庁の統計によると、南木曽観測所における既往最大１時間雨量は89ｍｍ、最大10分間雨量は18ｍｍであり、今回の降雨イベントにおいて既往最大値を上回る値を記録したわけではなかった。
【００６９】
（２）既往降雨イベントとの比較
（２−１）時間雨量、積算雨量、及び先行雨量の比較
梨子沢災害から３０年前の1984年まで遡り、その期間に発生した主な降雨イベントを抽出し、それら既往降雨イベントと、梨子沢災害をもたらした降雨イベントについて、南木曽観測所で観測された時間雨量（正時の値）と積算雨量を比較した。以下、抽出した既往降雨イベントを、その発生年を用いて「S60イベント」、「H16イベント」、「H12イベント」、「H18イベント」という。同様に、梨子沢災害をもたらした降雨イベントを「H26イベント」という。
図１４ａ〜１４ｄは、３つの既往降雨イベント（H16イベント、H12イベント、H18イベント）とH26イベントの時間雨量と積算雨量を示すグラフである。また、４つの既往イベント（S60イベント、H16イベント、H12イベント、H18イベント）と、H26イベントの諸元をまとめた表１を下記に示す。
【表１】
2015232537.tif 000005
【００７０】
なお、ここでは、24時間以上の無降雨期間によって降雨イベントを区切りっており、これによると、H26イベントは2014年7月6日の18:00に開始したことになる。図１４ｄ及び表１に示すように、H26イベントの、降雨開始から土石流発生までの積算雨量は194ｍｍであった。
【００７１】
図１４ａに示すH16イベントは、前述の既往最大１時間雨量（89ｍｍ）を観測したイベントであり、かなりの集中豪雨であったといえる。積算雨量は比較的少ないが、終盤にはH26イベントに匹敵する値にまで上昇していた。
図１４ｂに示すH12イベントは、時間雨量のピーク（最大値）は小さいものの比較的短時間にまとまった量の雨が降り、積算雨量はH26イベントよりも85ｍｍほど多くなった。
図１４ｃに示すH18イベントは、時間雨量のピークはH16イベントよりもさらに小さいが、降雨の継続時間が150時間と長く、積算雨量は非常に多くなった。
表１に示したS60イベントは、降雨継続時間が262時間と極めて長く、積算雨量はH26イベントの2.7倍に達したが、時間雨量のピークは小さかった。
表１から分かるように、各降雨イベントに先行する２日間の雨量（先行２日間雨量）は、H26イベントで最も多かったが、その他の先行雨量は、H26イベント以外の降雨イベントのいずれかが最大値を示した。
【００７２】
以上より、時間雨量、積算雨量、先行雨量の各々を指標とした場合に、H26イベントの特殊性は検出されず、H26イベントが既往降雨イベントを上回る異常な降雨イベントであったと結論づけることはできなかった。
【００７３】
（２−２）スネーク曲線を用いた比較
土砂災害警戒情報発表の基準には、横軸を土壌雨量指数、縦軸を時間雨量としたスネーク曲線が用いられることが多い。そこで、各降雨イベントのスネーク曲線を求めた。
まず、1984年1月1日の1:00 からH26イベント開始直前までの約30年間の時間雨量データ（南木曽観測所）を用いて土壌雨量指数を計算し、既往スネーク曲線を描いた（図１５ａ〜１５ｄ中の灰色丸印）。各図中の太い実線は既往スネーク曲線の最大値を結んだものであり、以下、「既往最大値線」という。この既往最大値線の右上の領域にスネーク曲線が入ると既往最大値を超過したことになる。
【００７４】
次に、既往降雨イベント（H16、H12、H18イベント）及びH26イベントのスネーク曲線を黒色丸印及び白色丸印で図１５ａ〜１５ｄ中に書き込んだ。既往降雨イベントのスネーク曲線はいずれもＡＮ境界線まで上昇しており、 H26イベント開始直前までの既往最大値の一部を形成する降雨イベントであったことが分かる。また、図１５ｄより、H26イベントのスネーク曲線は、土石流が発生した時刻になると既往最大値線を僅かに超えて既往最大値超過を示す領域に侵入していたことが分かる。
【００７５】
上述したように、土砂災害の警戒・避難には、図１５に示すスネーク曲線の他、横軸を半減期の長い実効雨量、縦軸を半減期の短い実効雨量とするスネーク曲線が活用されることも多い。そこで、図１５と同様の解析を、横軸を半減期51hの実効雨量、縦軸を半減期0.9hの実効雨量としたスネーク曲線を用いて行った。その結果を図１６に示す。図１６ｄをみると、H26イベントのスネーク曲線は土石流発生時には既往最大値線を顕著に超過しており、既往最大値超過が発生したことが明らかである。このような解析を、横軸実効雨量の半減期と縦軸実効雨量の半減期が異なる様々な組合せで行ったところ、横軸半減期が13h〜173h、縦軸半減期が0.1h〜2.0 hの範囲にあるときに、H26イベントの土石流発生時において既往最大値超過の発生が検出されることが分かった。
【００７６】
H26イベントでは、先行２日間の雨量が多く（表１）、さらにイベント開始以降の積算雨量が100ｍｍを超えたところに、既往最大値には及ばないが、かなりの強度の雨が降った（図１４ｄのイベント開始から72h経過付近）。上述した解析結果は、縦軸と横軸の雨量指標を適切に設定したスネーク曲線を用いることによって、このような特殊性を「既往最大値超過」として検出でき、その結果、既往降雨イベントを上回る異常な降雨イベントが可能であることを示唆するものであった。
【００７７】
（３）異常降雨イベント発生の予測
図１５ｄ及び図１６ｄでは、7月9日の17:00に既往範囲内にあったスネーク曲線が、その後、土石流発生までに降った49ｍｍの雨で既往最大値を超過した。このように、与えられた雨量を用いて既往最大値超過の有無を算出する代わりに、7月9日17:00の時点において、今後１時間にどれだけの雨が降ったら既往最大値を超過するかを計算した。その結果、図１５ｄの場合は47.1ｍｍ、図１６ｄの場合は35.6ｍｍの値が得られた。
【００７８】
このような計算を、横軸実効雨量の半減期（長半減期Ｍ１）と縦軸実効雨量の半減期（短半減期Ｍ２）を0.1ｈ〜3000ｈの範囲で変化させて作成したスネーク曲線を用いて行った。その結果を図１７に示す。図１７は、スネーク曲線から求められた既往最大値を超過する雨量を、横軸を「横軸実効雨量半減期(h)」、縦軸を「縦軸実効雨量半減期(h)」とする二次元座標上に示すと共に、雨量の大きさ毎に色分けしたグラフである。このグラフは、請求項７及び１４に係る二次元判定図に相当する。つまり、横軸実効雨量の半減期が40〜100hで、且つ縦軸実効雨量の半減期が1.3hの領域では既往最大値を超過する雨量は34.3ｍｍとなり、最も少なかった。実際には、これを上回る49ｍｍの雨が降った結果、土石流が発生した。
【００７９】
H26イベントで土石流発生の直接の誘因となった降雨（図１３ａ）のように、急激に発達する雲によりもたらされる集中豪雨の予報は難しく、今回の土石流の発生予測は困難であったのは事実である。ただし、（２−２）節で検討したように、H26イベントでは、先行雨量、ならびに降雨イベント開始以降の累積雨量の影響も大きかったと考えられる。このため、これらの影響を考慮した図１７の解析を行うことによって、17:00の時点で、34ｍｍという特別に大きくはない時間雨量によって、土石流をもたらす異常降雨イベントになることが予測できた可能性がある。
【００８０】
従って、降雨イベント中のある時点における、今後１時間の予測降雨量の横軸実効雨量及び縦軸実効雨量を計算し、その結果得られるスネーク曲線が既往最大値線を超えるか否かに基づき、土石流発生の可能性を予測することができる。予測降雨量は、気象データ（雨量、風向、風速、気温等）や天気図等から求めることができる。
具体的には、図１８（ａ）に示すように、17:00の時点における今後１時間の予測降雨量が49ｍｍのとき、その横軸実効雨量及び縦軸実効雨量の点（「18:00」で示す点）は既往最大値線を超えた領域に位置する。このことから、今後１時間の間に土石流が発生する危険性があると予測することができる。
【００８１】
また、今後１時間の降雨量を種々変化させたときのスネーク曲線を作成し、既往最大値線を超えたときの今後１時間の降雨量の中から最小値を求めて、この最小値と今後１時間の予測降雨量との比較から、土石流発生の可能性（土石流発生の危険性）を定量的に評価することも可能である。例えば、図１８（ｂ）に示すように、17:00の時点における今後１時間の降雨量が38ｍｍのときに、その横軸実効雨量及び縦軸実効雨量を示す点が既往最大値線上に位置する場合、この値（38ｍｍ）と予測降雨量（49ｍｍ）を比較し、今後１時間の間に土石流が発生する可能性の程度を判定することができる。例えば、既往最大値線上に位置するときの降雨量が予測降雨量の０〜５０％の場合は可能性が「非常に大きい」、５０〜８０％の場合は可能性が「大きい」、８０％〜１００％のときは可能性が「中位」、１００％〜のときは可能性が「小さい」といった判定が可能である。
【００８２】
さらに、上述した二次元判定図（図８）と本実施例で用いた二次元判定図（図１７）を用いて、土石流の発生を予測するシステムを構築することができる。例えば図１９（ａ）、（ｂ）は、ある地点における過去の降雨量データから得られた現時点における二種類の二次元判定図を示している。このシステムでは、これら二次元判定図が液晶ディスプレイ等に表示される。そして、液晶ディスプレイに表示された判定図のいずれかにおいてユーザが任意の点を選択（クリック）すると、該選択した点の長半減期Ｍ１と短半減期Ｍ２をそれぞれ横軸及び縦軸とするスネーク曲線が自動的に作成され（図１９（ｃ））、さらに、今後１時間の降雨量の違いによって得られる推定点が書き込まれる。図１９（ｃ）では、今後１時間の降雨量を0ｍｍ、20ｍｍ、38ｍｍ、40ｍｍ、49ｍｍ、60ｍｍの場合の推定点が表示された例を示す。
従って、ユーザは、推定点の位置から土石流の発生を引き起こす、今後１時間の降雨量を推定することができる。
【００８３】
このように、上記した第１実施例及び第２実施例はそれぞれ斜面崩壊予測方法として有用であるが、さらに、これらの実施例に示す方法を併用すれば、より一層斜面崩壊予測精度を向上させることができる。
【００８４】
本発明は直ちに実用化可能なものであり、斜面崩壊予測精度を向上させる効果も大きい。ただし、さらなる精度向上に向けて、以下の課題を挙げることができる。
【００８５】
実効雨量の既往最大値を算定する期間の設定方法の検討が必要である。上記実施例では、気象庁から公開され入手が容易なデータを用いることで、約30年間を解析の対象としたが、より長期間の降雨量データを利用すれば更に正確な既往最大値を算定することができ、崩壊予測における「空振り」を減らすことができる。その一方で、斜面の特性に経年変化があることを考慮すれば、むやみに過去に遡ることが必ずしも良いとは限らず、適切な期間があることも考えられる。より多くの災害事例を解析することで、知見を蓄積する必要がある。
【符号の説明】
【００８６】
１０…斜面崩壊予測装置
１１…データ受信部
１２…既往実効雨量算出部
１３…判定用ＸＹ座標平面作成部
１４…予測用実効雨量算出部
１５…判定部
１６…制御部
１７…二次元判定図作成部
２０…降雨量データ提供装置
３０…出力装置
【図１】
2015232537.tif 000006
【図２】
2015232537.tif 000007
【図４】
2015232537.tif 000008
【図６】
2015232537.tif 000009
【図７】
2015232537.tif 000010
【図８】
2015232537.tif 000011
【図９】
2015232537.tif 000012
【図１０】
2015232537.tif 000013
【図３】
2015232537.tif 000014
【図５】
2015232537.tif 000015
【図１１】
2015232537.tif 000016
【図１２】
2015232537.tif 000017
【図１３】
2015232537.tif 000018
【図１４】
2015232537.tif 000019
【図１５】
2015232537.tif 000020
【図１６】
2015232537.tif 000021
【図１７】
2015232537.tif 000022
【図１８】
2015232537.tif 000023
【図１９】
2015232537.tif 000024