特開 2023-3518 表層崩壊リスク評価装置、表層崩壊リスク評価方法及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023003518

(43)【公開日】2023-01-17

(54)【発明の名称】表層崩壊リスク評価装置、表層崩壊リスク評価方法及びプログラム

(51)【国際特許分類】

   G01W 1/00 20060101AFI20230110BHJP

   E02D 17/20 20060101ALI20230110BHJP

   G01D 21/00 20060101ALI20230110BHJP

   G08B 31/00 20060101ALI20230110BHJP

【ＦＩ】

G01W1/00 Z

E02D17/20 106

G01D21/00 D

G08B31/00 B

【審査請求】有

【請求項の数】6

【出願形態】ＯＬ

(21)【出願番号】P 2021104638

(22)【出願日】2021-06-24

(11)【特許番号】

(45)【特許公報発行日】2022-08-12

(71)【出願人】

【識別番号】520242399

【氏名又は名称】九州電力送配電株式会社

(71)【出願人】

【識別番号】501410779

【氏名又は名称】九州電技開発株式会社

(74)【代理人】

【識別番号】100136180

【弁理士】

【氏名又は名称】羽立 章二

(72)【発明者】

【氏名】小川 健

(72)【発明者】

【氏名】久保田 豊

(72)【発明者】

【氏名】堀本 さとみ

(72)【発明者】

【氏名】八嶋 俊雄

【テーマコード（参考）】

2D044

2F076

5C087

【Ｆターム（参考）】

2D044EA07

2F076BA01

2F076BA12

2F076BA18

2F076BB07

2F076BD14

2F076BD17

2F076BE03

2F076BE04

2F076BE06

2F076BE08

2F076BE09

2F076BE13

5C087AA02

5C087AA03

5C087AA04

5C087DD02

5C087EE08

(57)【要約】

【課題】 表層崩壊の危険性を評価することの精度を高めることに適した表層崩壊リスク評価装置等を提案する。
【解決手段】 表層崩壊リスク評価装置１において、地形評価部５は対象地点の地形を評価し、水流評価部７は周囲から水が集まる状況によって対象地点を評価する。危険度評価部９は、地形評価部５による地形評価に加えて、水流評価部７による谷あいや渓流などのように水が集まる個所による評価を加味して、分析対象地域における対象地点の斜面崩壊の危険性を評価する。
【選択図】 図１

【特許請求の範囲】

【請求項1】

分析対象地域における対象地点の斜面崩壊の危険性を評価する表層崩壊リスク評価装置であって、
基準水流量を記憶する基準水流量記憶部と、
前記分析対象地域において累積水流量が基準水流量以上である本流地点群を抽出する本流地点群抽出部と、
前記対象地点に対して前記本流地点群との距離を用いて水流評価値を計算する水流評価値計算部と、
前記水流評価値を用いて前記対象地点の表層崩壊の危険性を評価する危険度評価部を備える、表層崩壊リスク評価装置。

【請求項2】

前記基準水流量は、崩壊が生じた崩壊地点群に対して、前記崩壊地点群の各崩壊地点と前記本流地点群との距離に対する崩壊地点の個数の分布について、少なくとも、第１距離区分と、第１距離区分よりも遠い第２距離区分と、第２距離区分よりも遠い第３距離区分にそれぞれ属する崩壊地点の個数である第１個数と第２個数と第３個数について、第１個数が最も多く、第３個数が最も少なく、第２個数は第１個数よりも少なく第３個数よりも多い値である、請求項１記載の表層崩壊リスク評価装置。

【請求項3】

前記対象地点における地形を評価する地形評価値を計算する地形評価値計算部を備え、
前記危険度評価部は、前記地形評価値と前記水流評価値とを用いて危険評価値を計算して、前記危険評価値を用いて前記対象地点の表層崩壊の危険性を評価する、請求項１又は２に記載の表層崩壊リスク評価装置。

【請求項4】

前記地形評価値は、前記対象地点において、
縦断形状が凸であり、かつ、横断形状が凹である地形は、他の形状に比較して危険度が高いと評価される値であり、
地表の傾斜が大きいほど危険度が高いと評価される値であり、
前記水流評価値は、前記対象地点において、累積水流量が基準水流量以上である本流地点群に近いほど危険度が高いと評価される値である、請求項３記載の表層崩壊リスク評価装置。

【請求項5】

対象地点における斜面崩壊の危険性を評価する崩壊リスク評価方法であって、
崩壊リスク評価装置は、地形評価値計算部と、水流評価値計算部と、危険度評価部を備え、
前記地形評価値計算部及び前記水流評価値計算部が、それぞれ、地形評価値及び水流評価値を計算する評価値計算ステップと、
前記危険度評価部が、前記地形評価値及び前記水流評価値を用いて前記対象地点の表層崩壊の危険性を評価する危険度計算部を備え、
前記地形評価値は、前記対象地点において、
縦断形状が凸であり、かつ、横断形状が凹である地形は、他の形状に比較して危険度が高いと評価される値であり、
地表の傾斜が大きいほど危険度が高いと評価される値であり、
前記水流評価値は、前記対象地点において、累積水流量が基準水流量以上である本流地点群に近いほど危険度が高いと評価される値である、表層崩壊リスク評価方法。

【請求項6】

コンピュータを、請求項１から４のいずれかに記載の表層崩壊リスク評価装置として機能させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本願発明は、表層崩壊リスク評価装置、表層崩壊リスク評価方法及びプログラムに関し、特に、対象地点における表層崩壊の危険性を評価する表層崩壊リスク評価装置等に関する。

【背景技術】

【0002】

従来、降水量を用いて鉄塔付近の災害を予測することが知られていたところ、出願人は、地形情報を用いて危険斜面を抽出することを提案した（非特許文献１参照）。

【先行技術文献】

【非特許文献】

【0003】

【非特許文献1】九州電技開発株式会社，「電技だより 第65号 鉄塔敷地斜面災害予測システムの改良」，[online]，インターネット，＜URL：http://www.dengi.co.jp/dengidayori/no65/index.html＞

【発明の概要】

【発明が解決しようとする課題】

【0004】

非特許文献１記載の技術では、地形情報において、崩壊が生じた地形との類似性に着目するものであった。しかしながら、崩壊が生じた地形と同じ形状の地形であっても、崩壊の危険度までが同じとは限らなかった。そのため、より精度の高い危険度の判定が求められていた。

【0005】

そこで、本願発明は、表層崩壊の危険性を評価することの精度を高めることに適した表層崩壊リスク評価装置等を提案することを目的とする。

【課題を解決するための手段】

【0006】

本願発明の第１の観点は、分析対象地域における対象地点の斜面崩壊の危険性を評価する表層崩壊リスク評価装置であって、基準水流量を記憶する基準水流量記憶部と、前記分析対象地域において累積水流量が基準水流量以上である本流地点群を計算する本流地点群抽出部と、前記対象地点に対して前記本流地点群との距離を用いて水流評価値を計算する水流評価値計算部と、前記水流評価値を用いて前記対象地点の表層崩壊の危険性を評価する危険度評価部を備える。

【0007】

本願発明の第２の観点は、第１の観点の表層崩壊リスク評価装置であって、前記基準水流量は、崩壊が生じた崩壊地点群に対して、前記崩壊地点群の各崩壊地点と前記本流地点群との距離に対する崩壊地点の個数の分布について、少なくとも、第１距離区分と、第１距離区分よりも遠い第２距離区分と、第２距離区分よりも遠い第３距離区分にそれぞれ属する崩壊地点の個数である第１個数と第２個数と第３個数について、第１個数が最も多く、第３個数が最も少なく、第２個数は第１個数よりも少なく第３個数よりも多い値である。

【0008】

本願発明の第３の観点は、第１又は第２の観点の表層崩壊リスク評価装置であって、前記対象地点における地形を評価する地形評価値を計算する地形評価値計算部を備え、前記危険度評価部は、前記地形評価値と前記水流評価値とを用いて危険評価値を計算して、前記危険評価値を用いて前記対象地点の表層崩壊の危険性を評価する。

【0009】

本願発明の第４の観点は、第３の観点の表層崩壊リスク評価装置であって、前記地形評価値は、前記対象地点において、縦断形状が凸であり、かつ、横断形状が凹である地形は、他の形状に比較して危険度が高いと評価される値であり、地表の傾斜が大きいほど危険度が高いと評価される値であり、前記水流評価値は、前記対象地点において、累積水流量が基準水流量以上である本流地点群に近いほど危険度が高いと評価される値である。

【0010】

本願発明の第５の観点は、対象地点における斜面崩壊の危険性を評価する崩壊リスク評価方法であって、崩壊リスク評価装置は、地形評価値計算部と、水流評価値計算部と、危険度評価部を備え、前記地形評価値計算部及び前記水流評価値計算部が、それぞれ、地形評価値及び水流評価値を計算する評価値計算ステップと、前記危険度計算部が、前記地形評価値及び前記水流評価値を用いて前記対象地点の表層崩壊の危険性を評価する危険度計算部を備え、前記地形評価値は、前記対象地点において、縦断形状が凸であり、かつ、横断形状が凹である地形は、他の形状に比較して危険度が高いと評価される値であり、地表の傾斜が大きいほど危険度が高いと評価される値であり、前記水流評価値は、前記対象地点において、累積水流量が基準水流量以上である本流地点群に近いほど危険度が高いと評価される値である。

【0011】

本願発明の第６の観点は、コンピュータを、第１から第４のいずれかの観点の表層崩壊リスク評価装置として機能させるためのプログラムである。

【0012】

なお、本願発明を、本願発明の第６の観点のプログラムを記録するコンピュータ読み取り可能な記録媒体として捉えてもよい。

【発明の効果】

【0013】

非特許文献１記載の技術などでは、例えば各地点での降水量や地形のように、各地点を個別に分析する傾向にあった。出願人は、水流（例えば対象地域を区分したメッシュの間での水の流れなど）に着目し、実際に表層崩壊が生じた表層崩壊地点を分析した。各地点に近傍地点から流れる流量を累積した累積水流量に着目して、表層崩壊地点が、累積水流量が所定の値以上となる本流地点群に近いほど多く、遠いほど少なくなる関係を見出した。

【0014】

本願発明の各観点によれば、出願人の新たな知見に基づき、水流に着目して、対象地点に対して本流地点群との距離を用いて崩壊の危険性を評価することにより、表層崩壊の危険性の予測精度を高くすることができる。

【0015】

さらに、本願発明の第３の観点などにあるように、対象地点での地形の危険性を示す地形評価値を、水流評価値を用いて修正した危険評価値によって、水流評価値によって地形評価値による予測の精度を向上させることができる。

【図面の簡単な説明】

【0016】

【図1】本願発明の実施の形態に係る表層崩壊リスク評価装置１の（ａ）構成の一例を示すブロック図と、（ｂ）、（ｃ）及び（ｄ）動作の一例を示すフロー図である。

【図2】曲率を説明するための図である。

【図3】（ａ）、（ｂ）及び（ｃ）は、それぞれ、基準水流量が、２０、１００及び５０のときの、基準水流量の違いによる崩壊開始点の分布を示すグラフである。

【図4】基準水流量を２０とするときの、水流図を示す。

【図5】図３（ａ）において、崩壊開始点についての本流セル群との距離（ｍ）（横軸）と、データ数割合（％）との関係を示すグラフである。

【図6】地形評価値を水流評価値により補正した場合の一例を示す図である。

【発明を実施するための形態】

【0017】

以下では、図面を参照して、本願発明の実施例について説明する。なお、本願発明は、この実施例に限定されるものではない。

【実施例0018】

図１は、本願発明の実施の形態に係る表層崩壊リスク評価装置１の（ａ）構成の一例を示すブロック図と、（ｂ）、（ｃ）及び（ｄ）動作の一例を示すフロー図である。

【0019】

図１（ａ）を参照して、表層崩壊リスク評価装置１は、地形情報分析部３と、地形評価部５と、水流評価部７と、危険度評価部９（本願請求項の「危険度評価部」の一例）と、パラメータ記憶部１０を備える。

【0020】

地形情報分析部３とは、地理空間情報記憶部１１と、分割処理部１２と、形状解析部１３と、水文解析部１５を備える。

【0021】

地形評価部５とは、傾斜曲率計算部１７と、地形評価値計算部２１（本願請求項の「地形評価値計算部」の一例）を備える。

【0022】

水流評価部７は、基準水流量記憶部２３（本願請求項の「基準水流量記憶部」の一例）と、本流地点群抽出部２５（本願請求項の「本流地点群抽出部」の一例）と、水流評価値計算部２７（本願請求項の「水流評価値計算部」の一例）を備える。

【0023】

危険度評価部９は、危険評価値計算部２９と、危険度判定部３１を備える。

【0024】

パラメータ記憶部１０は、評価値などを計算するときに使用するパラメータを記憶する。

【0025】

地形情報分析部３において、地理空間情報記憶部１１は、地理空間情報を記憶する。地理空間情報は、例えば、空間上の特定の地点又は区域の位置を示す情報（位置情報）とそれに関連付けられた様々な事象に関する情報である。また、位置情報のみからなる情報であってもよい。地理空間情報には、地域における土地利用図、地質図、ハザードマップ等の主題図、都市計画図、地形図、地名情報、台帳情報、統計情報、空中写真、衛星画像などの情報がある。また、表層崩壊リスク評価装置１の利用者などが、現地で測定した情報であってもよい。

【0026】

分割処理部１２は、地理空間情報記憶部１１が記憶する地理空間情報から、危険性を評価する対象地点（本願請求項の「対象地点」の一例）を含む分析対象地域（本願請求項の「分析対象地域」の一例）を抽出して、分割情報を生成する。分割情報は、例えば、分析対象地域を一定間隔（例えば５ｍ間隔）で格子状に分割したメッシュデータである。以下では、分割で生じた各矩形地域を、セルという。

【0027】

形状解析部１３は、分割情報のそれぞれのセルに対して、地理空間情報での標高値データを対応づけて分割標高データとする。各セルに対応する標高値データは、例えば、地理空間情報での標高値データの一部又は全部であってもよく、これらの標高値データを統計処理したものであってもよい。

【0028】

水文解析部１５は、分割標高データのそれぞれの各セルに対して累積水流量を対応づけて地形情報とする。例えば、分割標高データを利用して、各セルから最も急な降下傾斜となる近傍セルへのフロー（流れ）の方向を求める。そして、各セルに対して、近傍セルから当該セルに流れるとする流量を累積して累積水流量を計算して地形情報とする。

【0029】

このように、地形情報は、例えば、分析対象地域を一定間隔で格子状に分割したメッシュデータで、各セルに対して地理空間情報での標高値データ及び累積水流量が対応づけられたものである。

【0030】

分割処理部１２、形状解析部１３及び水文解析部１５は、例えば、ＡｒｃＧＩＳ（ｅｓｒｉジャパン株式会社）などによって実現することができる。例えば、５ｍ間隔のＤＥＭデータを使用してラスターデータを作成し、小さな欠陥を取り除き、各セルから最も急な降下傾斜となる近傍セルへの流れの方向を示すラスターデータを生成して、各セルへ累積する流量のラスターデータを作成すればよい。

【0031】

図１（ｂ）は、地形評価部５の動作の一例を示すフロー図である。図２は、曲率を説明するための図である。図１（ｂ）及び図２を参照して、地形評価部５の動作を説明する。

【0032】

図１（ｂ）を参照して、傾斜曲率計算部１７は、地形情報を利用して、傾斜及び曲率を計算する（ステップＳＴＡ１）。地形評価値計算部２１は、地形評価値を計算する（ステップＳＴＡ２）。

【0033】

まず、傾斜について説明する。傾斜は、各セルの水平面からの傾きを示す地形量である。例えば、傾斜角は、対象セルとその近傍セルを比較し、それらの値の最大変化率を計算して求めることができる。この場合、原則として、対象セルとその８つの近傍セルの高度を比較したとき、距離に対するそれらの最大高度変化が対象セルからの最も急な降下傾斜となる。

【0034】

続いて、図２を参照して、曲率の計算について説明する。曲率は、斜面が尾根型か谷型かを表す数値である。縦断形状及び横断形状は、それぞれ、最大傾斜方向に平行な及び直交する方向の形状である。断面曲率及び平面曲率は、それぞれ、縦断方向及び横断方向の凹凸を示す曲率である。

【0035】

図２（ａ）は、断面曲率を説明するための図である。断面曲率は、斜面に平行であり、水の流れの加速と減速に影響する。図２（ａ）の左側（断面曲率が負の場合）の形状は、上方向に凸状であり、流れが減速する。図２（ａ）の中央（断面曲率が正の場合）の形状は、上方向に凹状であり、流れが加速する。図２（ｃ）の右側（断面曲率がゼロの場合）の形状は、表面が直線である。

【0036】

図２（ｂ）は、平面曲率を説明するための図である。断面曲率は、斜面に垂直であり、水の流れの収束と分岐に関係する。図２（ｂ）の左側（平面曲率が正の場合）の形状は、上方向に凸状であり、流れが分岐しやすい。図２（ｂ）の中央（平面曲率が負の場合）の形状は、上方向に凹状であり、流れが収束しやすい。図２（ｃ）の右側（平面曲率がゼロの場合）の形状は、表面が直線である。

【0037】

傾斜曲率計算部１７は、下記の式により、各セルの改合成曲率を計算する。ここで、ａは正の数であり、パラメータ記憶部１０に記憶されている。
改合成曲率＝平面曲率（凹が負）＋ａ×断面曲率（凸が負）

【0038】

なお、一般に、合成曲率は、凹型の傾向が強いように、「平面曲率（凹が負）－断面曲率（凸が負）」により計算される。そのため、上記の式により計算される値を「改合成曲率」という。

【0039】

なお、各セルの曲率の計算には、例えば、計算するセル及び近傍セルに含まれる標高値データを利用してもよい。

【0040】

実際に崩壊が生じた個所を調査した。以下では、崩壊が生じた個所を「崩壊箇所」といい、対比のために採用した、崩壊が生じていない個所を「健全個所」という。

【0041】

平面曲率については、健全個所では概ね－５～＋５の範囲にあり、崩壊個所では－１５～＋５であった。断面曲率については、健全個所では概ね－５～＋５の範囲にあり、崩壊個所では－５～＋１０であった。図２（ｃ）は、断面曲率と平面曲率を組み合わせたものである。崩壊箇所は、概ね、断面曲率が凸（負）で、かつ、平面曲率が凹（負）の個所であったことが確認された（図２（ｃ）の右上の欄を参照）。そこで、改合成曲率を斜面評価の指標の一つとして採用する。改合成曲率は、負の値が大きいほど崩壊が発生しやすい斜面という評価になる。

【0042】

傾斜について、健全個所では概ね３０°以下であるが、崩壊箇所では１０～５０°にばらついていた。傾斜と曲率との関係性を調べると、概ね、崩壊箇所の方が、傾斜角が大きい傾向がみられた。そこで、傾斜を、斜面評価の指標の一つとして採用する。傾斜は、値が大きいほど崩壊が発生しやすい斜面という評価になる。

【0043】

地形評価値計算部２１は、下記の式により、各セルの地形評価値を計算する。ここで、ｂは正の数であり、パラメータ記憶部１０に記憶されている。
地形評価値＝改合成曲率－ｂ×傾斜

【0044】

パラメータａ及びｂについては、健全個所及び崩壊箇所を用いて決定することができる。例えば、実際の崩壊箇所に対してａを変更して調べたところ、ａ＝１のときに改合成曲率が高くなるにつれて崩壊箇所数が多くなる傾向にあった。そのため、パラメータａ＝１を採用することができる。また、例えば、実際の崩壊箇所及び健全個所を調べると、ｂ＝０．３では、崩壊箇所のデータ数が－１０～－５より－５～０の方が多くなり、健全個所はすべて－１３以上となった。また、ｂ＝０．３５では、崩壊箇所に正の値がなく、健全個所はすべて－１５以上となった。ｂ＝０．４では、崩壊箇所に正の値がないが、健全個所が－２０～－１５と低い値にも及んでいた。この場合、ｂ＝０．３５を採用し、例えば、地形評価値が－１５以下で最も注意すべきとし、－１５よりも大きく－１０以下で次に注意すべきとし、－１０よりも大きく０以下で次に注意すべきとし、０よりも大きいと比較的安定していると評価する。

【0045】

図１（ｃ）は、水流評価部７の動作の一例を示すフロー図である。図３～図５は、水流評価値を説明するための図である。図１（ｃ）及び図３～図５を参照して、水流評価部７の動作を説明する。

【0046】

基準水流量記憶部２３は、基準水流量を記憶する。図１（ｃ）を参照して、本流地点群抽出部２５は、地形情報を利用して、累積水流量が基準水流量以上である本流セル群（本願請求項の「本流地点群」の一例）を抽出する（ステップＳＴＢ１）。水流評価値計算部２７は、各セルに対応して水流評価値を計算する（ステップＳＴＢ２）。水流評価値は、各セルと本流セル群との距離を用いて計算される値であり、例えば、距離が近いほど危険度が高く、遠いほど危険度が低いとされる値である。

【0047】

図３は、基準水流量の違いによる崩壊開始点の分布を示す。図３（ａ）、（ｂ）及び（ｃ）は、それぞれ、基準水流量が、２０、１００、及び、５０である。横軸は、崩壊開始点と本流セル群との距離（ｍ）である。縦軸は、データ数の割合（％）（各距離の崩壊開始点数について、崩壊開始点の全体数に対する割合）（四角のグラフ参照）と、各距離の崩壊開始点の地形評価値の平均（点）（三角のグラフ参照）である。図３を参照して、基準水流量について説明する。

【0048】

例えば、累積水流量が少ない個所は、周りから集まってくる水が少ない個所である。累積水流量が増えるにつれて、水が周りから集まる谷あいや渓流などが抽出できる。さらに累積水流量が増えると小さな河川（支流など）を形成する状況が抽出できる。累積水流量がさらに増えると、水がさらに集まって河川（主流など）を形成する状況を抽出できる。例えば、基準水流量が１００である場合は河川（本流など）を示し、基準水流量が５０である場合は小さな河川（支流など）を含めており、基準水流量が２０である場合は谷あいや渓流などをも含めていると考えられる。

【0049】

図３（ａ）を参照して、基準水流量が２０のとき、崩壊開始点が、本流セル群に近いほど多く、遠くなるにつれて少なくなる傾向がみられる（例えば、表１にあるように、本流からの距離が０ｍ（本願請求項の「第１距離区分」の一例）では４９個であり最も多く、本流からの距離が５０ｍよりも遠い（本願請求項の「第３距離区分」の一例）では１個であり最も少なく、本流からの距離が１０ｍよりも遠くて２０ｍ以下である場合（本願請求項の「第２距離区分」の一例）では２１個であり、これは４９個よりも少なく１個よりも多い。なお、本願発明では、崩壊開始点に限らず、例えば、最も崩壊箇所が大きいセルの個数などであってもよい。）。また、地形評価値について、本流セルに近いほど小さく、遠くなると大きくなる傾向がみられる。そのため、基準水流量は、２０を採用する。図４は、基準水流量を２０とするときの、水流図を示す。

【0050】

なお、図３（ｂ）及び（ｃ）を参照して、基準水流量が１００及び５０のときには、地形評価値について、本流セルに近いほど小さく、遠くなると大きくなる傾向がみられたものの、データ数の分布には明確な関連性が認められなかった。すなわち、一般には河川の水量が多い主流部分に着目しがちであるが、本願発明では、崩壊開始点の予測について基準水流量を採用することにより、一般的な分析では関心が及びにくい支流部分にも着目することができ、これによって崩壊開始点との関係性が明らかとなり、予測精度を向上させることができる。

【0051】

図５は、図３（ａ）において、崩壊開始点についての本流セル群との距離（ｍ）（横軸）と、データ数割合（％）との関係を示すグラフである。これを指数関数で近似したとき、近似式は、ｙ＝４８．１３８ｅ^-0.062xである。ここで、ｘは、本流セル群からの距離であり、ｙはセル数（割合％）である。

【0052】

表１は、近似曲線を使用して決定した水流評価値を示す。水流評価値計算部２７は、各セルに対応して、表１に従って、各セルの本流セル群からの距離に応じて水流評価値を求める。

【0053】

【表1】

【0054】

続いて、図１（ｄ）及び図６を参照して、危険度評価部９の動作を説明する。図６は、地形評価値を水流評価値により補正した場合の一例を示す図である。

【0055】

図１（ｄ）を参照して、危険評価値計算部２９は、各セルに対して、地形評価値を水流評価値により補正して、危険評価値を計算する（ステップＳＴＣ１）。危険度判定部３１は、各セルの危険評価値を利用して、各セルで表層崩壊が生じるリスクを判定する（ステップＳＴＣ２）。

【0056】

図６を参照して、危険評価値について説明する。図６（ａ）は、各セルの地形評価値の値を示す。図６（ａ）は、補正がない状態を示す。図６（ｂ）は、各セルの地形評価値に対して、水流評価値を乗除することにより補正した一例を示す。図６（ｃ）は、各セルの地形評価値に対して水流評価値を加減することにより補正した一例を示す。以下では、地形評価値を補正しない場合を「補正なし」という。地形評価値を水流評価値の乗除で補正する場合を「×補正」という。地形評価値を水流評価値の加減で補正する場合を「＋補正」ともいう。

【0057】

図６において危険度最小値が崩壊地内部で検出されたか否かを判定すると、図６（ａ）及び（ｃ）の補正なし及び＋補正では検出されたのに対し、図６（ｂ）の×補正では検出されなかった。同様にして合計２５箇所で判定したところ、補正なしの正当率は６５％、×補正の正答率も６５％であったのに対し、＋補正の正当率は７０％であった。そのため、＋補正が、最も精度よく検出できた。

【0058】

同様に、補正なし、×補正及び＋補正における全体の点数分布域を１０等分して１～１０ランク（１が最も危険）に区分し、崩壊開始点の平均ランクを求めた。補正なしは２．２１、×補正は３．０７であったのに対し、＋補正の正当率は２．１６であった。この場合、１に近いほど危険が高いと判定できている。そのため、＋補正が、最も精度よく検出できている。

【0059】

よって、危険評価値計算部２９は、例えば、次の式により危険評価値を計算する。ここで、水流評価値＝水流値／１０、水流値＝本流セル群（累積水流量≧２０）からの距離別崩壊数の近似式から求めた値、地形評価値＝改合成曲率－０．３５×傾斜、改合成曲率＝平面曲率（凹が負）＋断面曲率（凸が負）、である。なお、この例では、危険評価値は負であるほど危険と判定されるものであるため、水流評価値は減じる補正となる。
危険評価値＝地形評価値－水流評価値

【0060】

一般的には、危険評価値計算部２９は、パラメータ記憶部１０が記憶する正のパラメータａ、ｂ及びｃにより、次の式により危険評価値を計算することができる。ここで、水流評価値＝本流セル群（累積水流量≧２０）からの距離別崩壊数の近似式から求めた値、地形評価値＝改合成曲率－ｂ×傾斜、改合成曲率＝平面曲率（凹が負）＋ａ×断面曲率（凸が負）、である。なお、この例では、危険評価値は負であるほど危険と判定されるものであるため、水流評価値は減じる補正となる。
危険評価値＝地形評価値－ｃ×水流評価値

【0061】

なお、例えば、危険評価値計算部２９が、基準水流量並びにパラメータａ、ｂ及びｃを、実際に崩壊が生じた箇所を用いて決定してもよい。基準水流量は、例えば、崩壊開始点が本流セル群からの距離に応じて単調減少する傾向となるように決定する。このとき、指数関数などの所定の近似式を考慮して、近似式との誤差が最小となるようにしてもよい。パラメータａは、例えば、改合成曲率が高くなるにつれて崩壊箇所数が多くなる傾向がみられる値とする。パラメータｂは、崩壊箇所及び健全箇所の分布によって、崩壊箇所が小さくなって危険と判断され、かつ、健全箇所が大きくなって安全と判断される傾向となるように決定することができる。パラメータｃは、例えば、危険度最小値が崩壊地内部で検出された正答率が最も高くなるようにしたり、崩壊開始点の平均ランクが最も小さくなるようにしたりして、決定することができる。

【0062】

危険度判定部３１は、各セルの危険評価値について、小さい値であれば表層崩壊が生じるリスクが高いと判定し、大きい値であれば表層崩壊が生じるリスクが低いと判定する。

【符号の説明】

【0063】

１ 表層崩壊リスク評価装置、３ 地形情報分析部、５ 地形評価部、７ 水流評価部、９ 危険度評価部、１０ パラメータ記憶部、１１ 地理空間情報記憶部、１２ 分割処理部、１３ 形状解析部、１５ 水文解析部、１７ 傾斜曲率計算部、２１ 地形評価値計算部、２３ 基準水流量記憶部、２５ 本流地点群抽出部、２７ 水流評価値計算部、２９ 危険評価値計算部、３１ 危険度判定部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】