(57)【要約】【課題】地すべり地表面の変位ベクトルを計測し、その結果から地すべり面を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することを課題とする。【解決手段】地表線等のデータ３を編集して解析断面形状を作成する断面モデル作成手段５と、地表に設置された変位測定手段６からの計測データを入力する計測データ入力手段９と、前記断面モデル作成手段５および計測データ入力手段９からのデータに基づいてすべり面形状を推定するすべり面推定手段１１とを有し、該すべり面推定手段１１は、ブロック毎に地すべり面を２次多項式で設定する高次多項式設定部２１と、ブロック毎に最小二乗法を適用して高次多項式のパラメータを求めるパラメータ算出部２３と、算出された夫々のブロックの２次多項式を連ねてすべり面形状を推定するすべり面推定部２５とを備えること。【選択図】図１
2008202994.tif 000002
【特許請求の範囲】
【請求項１】
地表の変位測定値から２次元鉛直断面における地すべり面形状を推定する推定方法において、
地すべり土塊を鉛直方向の分割線によって複数のブロックに分割するとともに、変位測定手段が各ブロックの地表に設置されるように分割し、該ブロック毎の地すべり面を高次多項式で設定し、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも１つの既知の境界条件に基づいて求め、算出されたブロック毎の前記高次多項式を連ねて地すべり土塊全体のすべり面形状を推定することを特徴とする地すべり面形状の推定方法。
【請求項２】
前記境界条件が分割線上の等高および等勾配条件であり、端部ブロックではさらに地表境界条件を含むことを特徴とする請求項１記載の地すべり面形状の推定方法。
【請求項３】
地表の変位測定値から２次元鉛直断面における地すべり面形状を推定する推定装置において、
地表線、地表変位計測点、境界点、および分割線位置のデータを編集して解析断面形状を作成する断面モデル作成手段と、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段と、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段とを有し、前記すべり面推定手段は、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも１つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部とを備えたことを特徴とする地すべり面形状の推定装置。
【請求項４】
前記パラメータ算出部では、前記境界条件を分割線上の等高および等勾配条件から、さらに端部ブロックでは地表境界条件から求めることを特徴とする請求項３記載の地すべり面形状の推定装置。
【請求項５】
地盤の地すべり面形状を推定するためのコンピュータを、
地表線、地表変位計測点、境界点、および分割線位置のデータを編集して鉛直断面形状を作成する断面モデル作成手段、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段として機能させ、さらに前記すべり面推定手段を、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも１つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部として機能させるための地すべり面形状の推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地表面の変位計測結果に基づいて地すべりのすべり面形状を推定する推定方法、推定装置、および推定プログラムに関する。
【背景技術】
【０００２】
一般に、建設現場で地すべりが発生した場合、対応の遅れが被災規模の拡大を招くことから、発生初期の緊急対応が重要である。そのため、地すべりの発生とその規模を迅速に把握する必要がある。
しかしながら、地すべりの発生メカニズムとその規模を把握する際においては、通常のボーリング調査等では一般に日数を要し、しかも地すべり地内での作業であるため安全確保の点からも慎重な対応が求められるなど問題な点が多い。
このため、地表面の変位計測結果に基づき地すべり面を推定することが行われてきた。地すべり面を推定する手法は種々知られており、すべり面の形状を円弧形状によって表す円弧すべり面法や、直線を複数連ねて表す多角形法や、高次の多項式例えば２次式で表す多項式法が知られている。
そして、例えば、非特許文献１には、各種の回転半径をもつ円弧の結合による複雑な地すべり面形状を推定する論文が示され、また、非特許文献２には、地表面の計測変位から地すべり面を推定して斜面の安定性を評価する評価手法についての論文が示されている。
【０００３】
また、特許文献１には、ＧＰＳ衛星からの電波を受けて、斜面変位を計測し該計測結果からすべり面およびすべり層を設定して、すべり安全率を求めて安定性を評価する防災監視システムが示されている。
【０００４】
【非特許文献１】「地すべり」Ｖｏｌ．２８，Ｎｏ．４，ＰＰ．６−１４，１９９２．「各種の回転半径をもつ円弧の結合による複雑な地下すべり面形状の推定」
【非特許文献２】「土と基礎」Ｖｏｌ．４９，Ｎｏ．７，ＰＰ．１０−１２，２００１．「計測変位を用いた斜面の安定性評価法」
【特許文献１】特許第３７３５２８０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
前記したように地すべり面の推定については、種々の手法が知られているが、主に地すべりの実態の把握を目的とし、計測器には、伸縮計や傾斜計が慣用的に用いられ、地すべり土塊や斜面の移動・変状を観察するに留まっていた。
しかし、近年の測量技術や計測技術の進歩に伴って、斜面上地表点の変位計測が経済的にも、精度的にも可能になると、より高度なシステムして、地表面変位計測結果に基づくすべり面形状の推定手法が必要となり、例えば、航空写真判読や、光波測距儀を用いたトータルステーション測量によって滑動土塊地表点の３次元変位時刻歴を求めて、その変位から地下に潜在する滑り面の形状を推定しようとするものや、特許文献１に示されているようなＧＰＳ衛星からの電波によって地表面変位を求めてすべり面形状を推定するもの、さらには３Ｄ（３次元）レーザスキャナによって地表面変位を計測してすべり面形状を推定するもの等が行われている。
このように、地すべりや斜面崩壊の予知を対象に斜面上地表点の変位計測が経済的にも可能となるともに、技術的にも精度の高いレベルでの測定が可能になるに伴って、より精度の高い地表面変位計測結果に基づくすべり面形状の推定手法が望まれている。
【０００６】
そこで、本発明は、このような背景に鑑みなされたものであり、地すべり地表面の変位ベクトルを計測し、その結果から地すべり面を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することを課題とする。
【課題を解決するための手段】
【０００７】
前記課題を解決するため、請求項１記載の発明は、地表の変位測定値から２次元鉛直断面における地すべり面形状を推定する推定方法において、地すべり土塊を鉛直方向の分割線によって複数のブロックに分割するとともに、変位測定手段が各ブロックの地表に設置されるように分割し、該ブロック毎の地すべり面を高次多項式で設定し、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも１つの既知の境界条件に基づいて求め、算出されたブロック毎の前記高次多項式を連ねて地すべり土塊全体のすべり面形状を推定することを特徴とする。
【０００８】
かかる請求項１記載の発明によれば、地すべり土塊を、地表の変位測定手段を含む領域毎に複数のブロックに分割して、ブロック毎に高次多項式ですべり面の形状を推定するため、１つの多項式で全体のすべり面形状を設定するのに比べて、複数の多項式で推定するため大規模で複雑なすべり面形状をより適確に推定することができる。
【０００９】
また、変位測定手段によって地すべり地表面変位ベクトルを計測し、その結果からすべり面形状を推定するため、地すべりの発生初期において、変状の規模を迅速に把握して、すべり機構の迅速な把握とその後の対策工の実施に役立てることができる。
【００１０】
また、地すべり機構把握のための通常のボーリング調査等では一般に日数を要し、しかも地すべり地内での作業であるため安全確保の点からも慎重な対応が求められていたが、本発明により地すべり地内での危険な作業を減らすことが可能になる。
【００１１】
また、請求項２記載の発明は、請求項１において、前記境界条件が分割線上の等高および等勾配条件であり、端部ブロックではさらに地表境界条件を含むことを特徴とする。
かかる請求項２記載の発明によれば、ブロック間の連続条件である分割線上の等高および等勾配条件を未知パラメータ方程式に組み込むことによって、ブロック毎の最小二乗法の適用によって得られる未知パラメータ方程式だけでは解を求めることができない場合でも、高次多項式のパラメータを算出できるようになる。
【００１２】
また、請求項３記載の発明は、地すべり面形状の推定装置にかかる発明であり、地表の変位測定値から２次元鉛直断面における地すべり面形状を推定する推定装置において、地表線、地表変位計測点、境界点、および分割線位置のデータを編集して解析断面形状を作成する断面モデル作成手段と、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段と、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段とを有し、前記すべり面推定手段は、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも１つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部とを備えたことを特徴とする。
【００１３】
かかる請求項３記載の発明によれば、請求項１と同様の作用効果を得ることができるとともに、大規模で複雑なすべり面形状をより適確に推定することができる装置を簡単な構成で達成することができる。
【００１４】
また、請求項４記載の発明は、請求項３において、前記パラメータ算出部では、前記境界条件を分割線上の等高および等勾配条件から、さらに端部ブロックでは地表境界条件から求めることを特徴とし、請求項２と同様に、ブロック間の連続条件である分割線線上の等高および等勾配条件を未知パラメータ方程式に組み込むことによって、ブロック毎の最小二乗法の適用によって得られる未知パラメータ方程式だけでは解を求めることができない場合でも、高次多項式のパラメータを算出できるようになる。
【００１５】
また、請求項５記載の発明は、地すべり面形状の推定プログラムにかかる発明であり、地盤の地すべり面形状を推定するためのコンピュータを、地表線、地表変位計測点、境界点、および分割線位置のデータを編集して鉛直断面形状を作成する断面モデル作成手段、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段として機能させ、さらに前記すべり面推定手段を、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも１つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部として機能させるための地すべり面形状の推定プログラムである。
【００１６】
かかる請求項５記載の発明によれば、プログラムをコンピータに読み取らせて、地すべり面の形状を推定させることができるので、地すべり面の形状推定をコンピュータに実行させて、大規模で複雑なすべり面形状をより迅速に、そして適確に得ることができる。
【発明の効果】
【００１７】
本発明によれば、地すべり地表面の変位ベクトルを計測し、その結果から大規模で複雑なすべり面形状を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することができる。
【発明を実施するための最良の形態】
【００１８】
以下、図面を参照して本発明の好適な実施の形態を例示的に詳しく説明する。但しこの実施の形態に記載されている構成部品の寸法、材質、形状、その相対的配置等は特に特定的な記載がない限りは、この発明の範囲をそれに限定する趣旨ではなく、単なる説明例に過ぎない。
【００１９】
参照する図面において、図１は本発明に係る地すべり面の推定方法が適用される地すべり面推定装置の全体構成を示すブロック図である。図２はすべり面座標に解析対象の断面地形図を示した説明図である。図３はブロック毎に分割したすべり面推定の説明図である。図４はすべり面推定手段における制御フローである。
【００２０】
地盤の地すべり面を推定する地すべり面推定装置１は、地表線、計測点、境界点等のデータ３の入力支援を行う断面モデル作成手段５、地表の変位測定手段６からの計測データ７を入力する計測データ入力手段９、前記断面モデル作成手段５および計測データ入力手段９からのデータに基づいてすべり面形状を推定するすべり面推定手段１１、さらに、すべり面推定結果のすべり面形状を画面表示するすべり面形状表示部１３を含んで構成される。
【００２１】
断面モデル作成手段５は、解析対象の地表面形状を作成する。地表面形状のデータを直接入力するか、または地表面形状が既にＣＡＤデータとして与えられている場合には該データを読み込む等によって、解析対象の断面地表形状のデータを入力するとともにデータベース１５に記憶する。また、斜面に複数個所設置される変位測定手段６の設置位置の座標データや、後述する境界条件としての地表面に現れるすべり発生位置の頂上境界点データ、またはすべり終了位置の脚部境界点データを入力するとともにデータベース１５に記憶する。
【００２２】
また、変位測定手段６の設置位置のデータが入力されることに伴って、該設置位置を１箇所含めて、地すべり土塊を水平方向に複数ブロックに分割するように鉛直方向に延びるブロック境界線（分割線）１７の位置情報が設定される。このブロック境界線１７の位置は、変位測定手段６の設置位置の中間点に自動的に設定されるようになっているが、作業者によって任意の位置座標を入力して設定することもできる。
【００２３】
計測データ入力手段９は、傾斜の地表面に設置された変位測定手段６からのデータが定期的に入力されるようにオンライン、または無線によって受信されるようになっている。なお、変位測定手段６はグローバルポジショニングシステム（ＧＰＳ）によるもの、トータルステーション、レーザスキャナ、または変位計等によるものいずれの機器を用いても良い。
【００２４】
すべり面推定手段１１は、前記断面モデル作成手段５で設定されたブロック境界線１７によって分割されたブロック毎に、地すべり面を高次多項式で設定する高次多項式設定部２１と、該ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも１つの既知の境界条件に基づいて求めるパラメータ算出部２３と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を形成するすべり面推定部２５とを有して構成されている。
【００２５】
高次多項式設定部２１では、二次多項式を採用してすべり面形状を推定するように構成されている。
パラメータ算出部２３では、すべり土塊がすべり面上を剛体として移動すると仮定しているので、地表計測点からすべり面までの深さＤが一定に保たれるから、異なる計測時点の計測座標を用いて、すべり面座標の多項式を表すことができる。そして、計測点の座標は初期座標に計測変位を加えた測定値であるため誤差を伴っており、その誤差を、最小二乗法を適用して解析し、二次多項式のパラメータを求める。このとき境界条件として、断面モデル作成手段５によって入力設定された境界点データ、さらにブロック境界線１７上の等高および等勾配条件を境界条件として算出される。このすべり面推定手段１１の詳細手順については後述する。
【００２６】
すべり面推定手段１１のすべり面推定部２５では、高次多項式設定部２１で設定した二次多項式にパラメータ算出部２３で求めたパラメータを用いて、ブロックごとの二次多項式求め、該求めた二次多項式を連ねて、ブロック全体のすべり面形状を推定する。
【００２７】
すべり面形状表示部１３は、すべり面推定手段１１によって推定されたすべり面形状を画面イメージとして出力し、断面モデル作成手段５で作成された解析対象の地表面形状とともにすべり面形状を画面表示部２７に表示する。
【００２８】
次に、すべり面推定手段１１における推定手順について、図２〜４を参照して説明する。
図４に示すフローチャートにおいて、ステップＳ１で開始すると、まずステップＳ２で断面モデル作成手段５によって作成された解析対象の断面モデルを入手する。次にステップＳ３では、断面モデル作成手段５によって設定されたブロック境界線１７によって、地すべり土塊を水平方向に複数ブロックＢ１、Ｂ２、Ｂ３（図３参照）に分割する。各ブロックには変位測定手段６が１つ設置されている。そして、ステップＳ４で、ブロック毎のすべり面形状を２次多項式ｚ＝ａｓ２＋ｂｓ＋ｃによって設定する。
なお、原理的には２次式以外に次の（１）式のような一般的な高次多項式を用いることができるが、後述する最小二乗法を適用する場合の数値計算上の精度を考慮して、本実施形態では２次多項式を採用した。
【００２９】
【数１】
2008202994.tif 000003
ここにｚ:すべり面を投影した鉛直断面内の鉛直座標
ｓ:        〃                  の水平座標
ａｋ:すべり面形状を表す多項式の係数
【００３０】
次に、ステップＳ５で異なる計測時点（ｔ＝ｔ１、ｔ２、ｔ３）の計測点座標を用いて、すべり面座標の多項式を作成する。
本すべり面推定法の基本仮定によると、すべり土塊はすべり面上を剛体的に移動すると仮定しているから、地表変位計測点の変位軌跡は、地下のすべり面形状を表していることになる。図２に示すように、今、一つの変位測定手段６が設置されている地表変位計測点（座標値（ｓｉ，ｚｉ））の変位軌跡に注目すると、変位を生じている期間中に地表計測点からすべり面までの深さＤは一定に保たれるから、異なる計測時点（ｔ＝ｔ１，ｔ２，ｔ３）の計測点座標を用いて、すべり面座標の多項式を次の（２）式のように表すことができる。なお、座標値（ｓｉ，ｚｉ）は、水平方向位置をｓｉ、鉛直方向位置をｚｉとする。またｎブロック中のｉ番目のブロック（測定点ｉ）における２次多項式を示すものとする。
【００３１】
【数２】
2008202994.tif 000004
【００３２】
次に、ステップＳ６で、式（２）の誤差式を作成する。
すべり面多項式のパラメータａｉ、ｂｉ、ｃｉ、Ｄｉに関するこのような方程式は、計測点数×計測時点数（時刻履歴）だけ作ることができるが、図２に示すように計測点の座標値（ｓｉ，ｚｉ）は初期座標に計測変位を加えた測定値であるため誤差を伴っている。前記（２）式の誤差をＶｉ１，Ｖｉ２，Ｖｉ３とすると、これらはそれぞれ式（３）のように表される。
【００３３】
【数３】
2008202994.tif 000005
【００３４】
次に、ステップＳ７で、すべり面多項式のパラメータａｉ、ｂｉ、ｃｉ、Ｄｉを求めるには、最小二乗法を適用して、これら誤差の二乗和を最小化する解を求める。すなわち、次の式（４）、式（５）の解を求めることになる。
Ｖは、ａ、ｂ、ｃ、Ｄの４変数の２次関数であり、しかもａ２  ，ｂ２，ｃ２，Ｄ２の係数は正であるからＶが、ａ、ｂ、ｃ、Ｄのそれぞれに関して最小値をもつ。
【００３５】
【数４】
2008202994.tif 000006
【００３６】
【数５】
2008202994.tif 000007
【００３７】
一般に、ｎ個の計測点に対して、ｍ個の計測時点座標（ｍ−１個の変位時刻歴）が存在する場合、ｎ×ｍ個の誤差は次の式（６）と表すことができる。
【００３８】
【数６】
2008202994.tif 000008
【００３９】
式（６）を｛Ｖ｝＝［Ｓ］｛Ａ｝−［Ｚ］と表すと最小二乗法では、次の式（７）を行うことになるから、たとえばａの偏微分式からは式（８）の方程式が導かれる。
【００４０】
【数７】
2008202994.tif 000009
【数８】
2008202994.tif 000010
【００４１】
同様の演算を行うと、（ｎ＋３）個の未知パラメータ｛Ａ｝に関する次式（９）の正規方程式が得られることになる。
【数９】
2008202994.tif 000011
【００４２】
次に、図４（ａ）のフローチャートから図４（ｂ）のフローチャートへ進み、ステップＳ８で、前記した１計測点の３時点（ｔ＝ｔ１、ｔ２、ｔ３）の計測点座標（Ｓｉ，Ｚｉ）が得られているケースにあてはめると、次の式（１０）なる方程式を得る。
【００４３】
【数１０】
2008202994.tif 000012
【００４４】
式（１０）は明らかに第３式と第４式の方程式が同一であり、この方程式は不定（係数マトリックス［Ｓ］が特異）となる。
このことは、計測点数や計測時点数が増加しても解消されることはなく、常に第３式は第４式以降と一次従属関係にある。すなわち、第４式以降の式の和が第３式と同一になる。
そこで、本発明では、ステップＳ９で示すように全体方程式の行列式が特異にならないように、各ブロックの正規方程式の第４式を取り除き、その代わりに、次の式（１１）〜（１４）に示す新たな境界条件による未知パラメータ方程式を組み込む。
【００４５】
式（１１）は図３に示すようなブロック境界線１７上の等高さ条件式である。例えば、図３の第２ブロックＢ２と第３ブロックＢ３との間のブロック境界線１７上の点Ｐにおいて等しい高さ（Ｚ）を有する。
【数１１】
2008202994.tif 000013
【００４６】
式（１２）、（１３）は端部ブロック場合の条件を示す。
図３に示すように第１ブロックＢ１の場合には、脚部境界点（ＳＨ、ＺＨ）の条件、また頂部に位置するブロックＢ３の場合には、頂部境界点（ＳＴ、ＺＴ）の条件を示している。
【数１２】
2008202994.tif 000014
【数１３】
2008202994.tif 000015
【００４７】
式（１４）は、ブロック境界線１７上における隣接多項式の等勾配条件を設定している。この条件は、ブロック境界線１７上での等高さ条件と同様に、ブロック間連成方程式となる。
【００４８】
【数１４】
2008202994.tif 000016
【００４９】
次に、ステップＳ１０で、式（１０）の第１式〜第３式の３つの式と、式（１１）、（１２）、（１３）、（１４）の式によって、ｉブロックにおける多項式パラメータ（ａｉ、ｂｉ、ｃｉ、Ｄｉ）を算出する。
なお、式（１１）、（１４）についてはブロック間の連成方程式であるため隣接するブロックの方程式が必要であり、また、式（１２）、（１３）については端部ブロックの地表境界点のデータが必要となる。このため結果的には全ブロックにおける未知パラメータ方程式を行列化して、ｎブロック全てのパラメータについて（ａｉ、ｂｉ、ｃｉ、Ｄｉ）×ｎ（ブロック数）を求める必要がある。
また、式（１４）の等勾配条件式については、特に用いなくても未知パラメータの解を求めることができるが、式（１４）があることによって、少ない計測時点であっても適切に解を求めることができるようになる。
【００５０】
さらに、次のステップＳ１１で、求められたブロック毎の２次多項式を連ねて全体ブロックの全体方程式としての２次多項式を作成する。
すなわち、第１ブロックＢ１の領域では、第１ブロックＢ１の多項式パラメータ（ａ１、ｂ１、ｃ１、Ｄ１）を用いた２次多項式が用いられ、第２ブロックＢ２の領域では、第２ブロックＢ２の多項式パラメータ（ａ２、ｂ２、ｃ２、Ｄ２）を用いた２次多項式が用いられ、第３ブロックＢ３の領域では、第３ブロックＢ３の多項式パラメータ（ａ３、ｂ３、ｃ３、Ｄ３）を用いた２次多項式が用いられて、全体として２次多項式によってすべり面形状が推定される。そしてステップＳ１１で終了する。
【００５１】
以上のようにブロック間の連続条件であるブロック境界線上の等高および等勾配条件を未知パラメータ方程式に組み込むことによって、すなわち、式（１１）、（１２）、（１３）、（１４）を組み込むことによって、ブロック毎の最小二乗法の適用によって得られる式（１０）の第１式〜第３式だけではパラメータを求めることができなかった場合でも、２次多項式のパラメータを算出できるようになる。
また、境界条件として、地表面表れている脚部、頂部の境界点データでなく、予め地中調査によって得られている地中境界データがあれば、それを未知パラメータ方程式に組み入れても良い。
さらに、本実施形態の例では、１計測点の３時点（ｔ＝ｔ１、ｔ２、ｔ３）の計測点座標（Ｓｉ，Ｚｉ）について説明したが、ブロック内で最低限１個の変位ベクトルが得られていれば解を得ることができる。
【００５２】
以上本実施の形態によれば、地すべり土塊を、地表の変位測定手段６を含む領域毎に複数のブロックＢ１、Ｂ２、Ｂ３に分割して、ブロック毎に２次多項式ですべり面の形状を設定するため、１つの多項式で全体のすべり面形状を設定するのに比べて、複数の多項式で推定するため大規模で複雑なすべり面形状をより迅速に適確に推定することができる。
【００５３】
また、変位測定手段６によって地すべり地表面変位ベクトルを計測し、その結果からすべり面形状を推定するため、地すべりの発生初期において、変状の規模を迅速に把握して、すべり機構の迅速な把握とその後の対策工の実施に役立つことができる。
【００５４】
また、地すべり機構把握のための通常のボーリング調査等では一般に日数を要し、しかも地すべり地内での作業であるため安全確保の点からも慎重な対応が求められていたが、本発明により地すべり地内での危険な作業を減らすことが可能になる。
【００５５】
また、本実施形態によれば、前記のすべり面推定手段１１を備えた地すべり面推定装置１によれば、ブロック毎に２次多項式ですべり面の形状を設定して大規模で複雑なすべり面形状をより適確に推定することができる装置を簡単な構成で達成することができる。
【００５６】
さらに、本実施形態によれば、プログラムをコンピータに読み取らせて、地すべり面の形状を推定させることができるので、地すべり面の形状推定をコンピュータに実行させて、大規模で複雑なすべり面形状をより迅速に、そして適確に得ることができる。
【産業上の利用可能性】
【００５７】
本発明によれば、地すべり地表面の変位ベクトルを計測し、その結果から大規模で複雑なすべり面形状を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することができるので、地すべり面形状を推定する装置への適用に際して有益である。
【図面の簡単な説明】
【００５８】
【図１】本発明に係る地すべり面の推定方法が適用される地すべり面推定装置の全体構成を示すブロック図である。
【図２】すべり面座標に解析対象の断面地形図を示した説明図である。
【図３】ブロック毎に分割したすべり面推定の説明図である。
【図４（ａ）】すべり面推定手段における制御フローチャートである。
【図４（ｂ）】すべり面推定手段における制御フローチャートである。
【符号の説明】
【００５９】
１              地すべり面推定装置
５              断面モデル作成手段
６              変位測定手段
９              計測データ入力部
１１            すべり面推定手段
１３            すべり面形状表示部
１５            データベース
１７            ブロック境界線（分割線）
２１            高次多項式設定部
２３            パラメータ算出部
２５            すべり面推定部管路部品
Ｂ１            第１ブロック
Ｂ２            第２ブロック
Ｂ３            第３ブロック
【図１】
2008202994.tif 000017
【図２】
2008202994.tif 000018
【図３】
2008202994.tif 000019
【図４（ａ）】
2008202994.tif 000020
【図４（ｂ）】
2008202994.tif 000021