JP2008202994A - 地すべり面形状の推定方法、推定装置および推定プログラム - Google Patents
地すべり面形状の推定方法、推定装置および推定プログラム Download PDFInfo
- Publication number
- JP2008202994A JP2008202994A JP2007037150A JP2007037150A JP2008202994A JP 2008202994 A JP2008202994 A JP 2008202994A JP 2007037150 A JP2007037150 A JP 2007037150A JP 2007037150 A JP2007037150 A JP 2007037150A JP 2008202994 A JP2008202994 A JP 2008202994A
- Authority
- JP
- Japan
- Prior art keywords
- block
- landslide
- ground
- estimating
- surface shape
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
【課題】地すべり地表面の変位ベクトルを計測し、その結果から地すべり面を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することを課題とする。
【解決手段】地表線等のデータ3を編集して解析断面形状を作成する断面モデル作成手段5と、地表に設置された変位測定手段6からの計測データを入力する計測データ入力手段9と、前記断面モデル作成手段5および計測データ入力手段9からのデータに基づいてすべり面形状を推定するすべり面推定手段11とを有し、該すべり面推定手段11は、ブロック毎に地すべり面を2次多項式で設定する高次多項式設定部21と、ブロック毎に最小二乗法を適用して高次多項式のパラメータを求めるパラメータ算出部23と、算出された夫々のブロックの2次多項式を連ねてすべり面形状を推定するすべり面推定部25とを備えること。
【選択図】図1
【解決手段】地表線等のデータ3を編集して解析断面形状を作成する断面モデル作成手段5と、地表に設置された変位測定手段6からの計測データを入力する計測データ入力手段9と、前記断面モデル作成手段5および計測データ入力手段9からのデータに基づいてすべり面形状を推定するすべり面推定手段11とを有し、該すべり面推定手段11は、ブロック毎に地すべり面を2次多項式で設定する高次多項式設定部21と、ブロック毎に最小二乗法を適用して高次多項式のパラメータを求めるパラメータ算出部23と、算出された夫々のブロックの2次多項式を連ねてすべり面形状を推定するすべり面推定部25とを備えること。
【選択図】図1
Description
本発明は、地表面の変位計測結果に基づいて地すべりのすべり面形状を推定する推定方法、推定装置、および推定プログラムに関する。
一般に、建設現場で地すべりが発生した場合、対応の遅れが被災規模の拡大を招くことから、発生初期の緊急対応が重要である。そのため、地すべりの発生とその規模を迅速に把握する必要がある。
しかしながら、地すべりの発生メカニズムとその規模を把握する際においては、通常のボーリング調査等では一般に日数を要し、しかも地すべり地内での作業であるため安全確保の点からも慎重な対応が求められるなど問題な点が多い。
このため、地表面の変位計測結果に基づき地すべり面を推定することが行われてきた。地すべり面を推定する手法は種々知られており、すべり面の形状を円弧形状によって表す円弧すべり面法や、直線を複数連ねて表す多角形法や、高次の多項式例えば2次式で表す多項式法が知られている。
そして、例えば、非特許文献1には、各種の回転半径をもつ円弧の結合による複雑な地すべり面形状を推定する論文が示され、また、非特許文献2には、地表面の計測変位から地すべり面を推定して斜面の安定性を評価する評価手法についての論文が示されている。
しかしながら、地すべりの発生メカニズムとその規模を把握する際においては、通常のボーリング調査等では一般に日数を要し、しかも地すべり地内での作業であるため安全確保の点からも慎重な対応が求められるなど問題な点が多い。
このため、地表面の変位計測結果に基づき地すべり面を推定することが行われてきた。地すべり面を推定する手法は種々知られており、すべり面の形状を円弧形状によって表す円弧すべり面法や、直線を複数連ねて表す多角形法や、高次の多項式例えば2次式で表す多項式法が知られている。
そして、例えば、非特許文献1には、各種の回転半径をもつ円弧の結合による複雑な地すべり面形状を推定する論文が示され、また、非特許文献2には、地表面の計測変位から地すべり面を推定して斜面の安定性を評価する評価手法についての論文が示されている。
また、特許文献1には、GPS衛星からの電波を受けて、斜面変位を計測し該計測結果からすべり面およびすべり層を設定して、すべり安全率を求めて安定性を評価する防災監視システムが示されている。
「地すべり」Vol.28,No.4,PP.6−14,1992.「各種の回転半径をもつ円弧の結合による複雑な地下すべり面形状の推定」
「土と基礎」Vol.49,No.7,PP.10−12,2001.「計測変位を用いた斜面の安定性評価法」
特許第3735280号公報
前記したように地すべり面の推定については、種々の手法が知られているが、主に地すべりの実態の把握を目的とし、計測器には、伸縮計や傾斜計が慣用的に用いられ、地すべり土塊や斜面の移動・変状を観察するに留まっていた。
しかし、近年の測量技術や計測技術の進歩に伴って、斜面上地表点の変位計測が経済的にも、精度的にも可能になると、より高度なシステムして、地表面変位計測結果に基づくすべり面形状の推定手法が必要となり、例えば、航空写真判読や、光波測距儀を用いたトータルステーション測量によって滑動土塊地表点の3次元変位時刻歴を求めて、その変位から地下に潜在する滑り面の形状を推定しようとするものや、特許文献1に示されているようなGPS衛星からの電波によって地表面変位を求めてすべり面形状を推定するもの、さらには3D(3次元)レーザスキャナによって地表面変位を計測してすべり面形状を推定するもの等が行われている。
このように、地すべりや斜面崩壊の予知を対象に斜面上地表点の変位計測が経済的にも可能となるともに、技術的にも精度の高いレベルでの測定が可能になるに伴って、より精度の高い地表面変位計測結果に基づくすべり面形状の推定手法が望まれている。
しかし、近年の測量技術や計測技術の進歩に伴って、斜面上地表点の変位計測が経済的にも、精度的にも可能になると、より高度なシステムして、地表面変位計測結果に基づくすべり面形状の推定手法が必要となり、例えば、航空写真判読や、光波測距儀を用いたトータルステーション測量によって滑動土塊地表点の3次元変位時刻歴を求めて、その変位から地下に潜在する滑り面の形状を推定しようとするものや、特許文献1に示されているようなGPS衛星からの電波によって地表面変位を求めてすべり面形状を推定するもの、さらには3D(3次元)レーザスキャナによって地表面変位を計測してすべり面形状を推定するもの等が行われている。
このように、地すべりや斜面崩壊の予知を対象に斜面上地表点の変位計測が経済的にも可能となるともに、技術的にも精度の高いレベルでの測定が可能になるに伴って、より精度の高い地表面変位計測結果に基づくすべり面形状の推定手法が望まれている。
そこで、本発明は、このような背景に鑑みなされたものであり、地すべり地表面の変位ベクトルを計測し、その結果から地すべり面を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することを課題とする。
前記課題を解決するため、請求項1記載の発明は、地表の変位測定値から2次元鉛直断面における地すべり面形状を推定する推定方法において、地すべり土塊を鉛直方向の分割線によって複数のブロックに分割するとともに、変位測定手段が各ブロックの地表に設置されるように分割し、該ブロック毎の地すべり面を高次多項式で設定し、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも1つの既知の境界条件に基づいて求め、算出されたブロック毎の前記高次多項式を連ねて地すべり土塊全体のすべり面形状を推定することを特徴とする。
かかる請求項1記載の発明によれば、地すべり土塊を、地表の変位測定手段を含む領域毎に複数のブロックに分割して、ブロック毎に高次多項式ですべり面の形状を推定するため、1つの多項式で全体のすべり面形状を設定するのに比べて、複数の多項式で推定するため大規模で複雑なすべり面形状をより適確に推定することができる。
また、変位測定手段によって地すべり地表面変位ベクトルを計測し、その結果からすべり面形状を推定するため、地すべりの発生初期において、変状の規模を迅速に把握して、すべり機構の迅速な把握とその後の対策工の実施に役立てることができる。
また、地すべり機構把握のための通常のボーリング調査等では一般に日数を要し、しかも地すべり地内での作業であるため安全確保の点からも慎重な対応が求められていたが、本発明により地すべり地内での危険な作業を減らすことが可能になる。
また、請求項2記載の発明は、請求項1において、前記境界条件が分割線上の等高および等勾配条件であり、端部ブロックではさらに地表境界条件を含むことを特徴とする。
かかる請求項2記載の発明によれば、ブロック間の連続条件である分割線上の等高および等勾配条件を未知パラメータ方程式に組み込むことによって、ブロック毎の最小二乗法の適用によって得られる未知パラメータ方程式だけでは解を求めることができない場合でも、高次多項式のパラメータを算出できるようになる。
かかる請求項2記載の発明によれば、ブロック間の連続条件である分割線上の等高および等勾配条件を未知パラメータ方程式に組み込むことによって、ブロック毎の最小二乗法の適用によって得られる未知パラメータ方程式だけでは解を求めることができない場合でも、高次多項式のパラメータを算出できるようになる。
また、請求項3記載の発明は、地すべり面形状の推定装置にかかる発明であり、地表の変位測定値から2次元鉛直断面における地すべり面形状を推定する推定装置において、地表線、地表変位計測点、境界点、および分割線位置のデータを編集して解析断面形状を作成する断面モデル作成手段と、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段と、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段とを有し、前記すべり面推定手段は、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも1つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部とを備えたことを特徴とする。
かかる請求項3記載の発明によれば、請求項1と同様の作用効果を得ることができるとともに、大規模で複雑なすべり面形状をより適確に推定することができる装置を簡単な構成で達成することができる。
また、請求項4記載の発明は、請求項3において、前記パラメータ算出部では、前記境界条件を分割線上の等高および等勾配条件から、さらに端部ブロックでは地表境界条件から求めることを特徴とし、請求項2と同様に、ブロック間の連続条件である分割線線上の等高および等勾配条件を未知パラメータ方程式に組み込むことによって、ブロック毎の最小二乗法の適用によって得られる未知パラメータ方程式だけでは解を求めることができない場合でも、高次多項式のパラメータを算出できるようになる。
また、請求項5記載の発明は、地すべり面形状の推定プログラムにかかる発明であり、地盤の地すべり面形状を推定するためのコンピュータを、地表線、地表変位計測点、境界点、および分割線位置のデータを編集して鉛直断面形状を作成する断面モデル作成手段、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段として機能させ、さらに前記すべり面推定手段を、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも1つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部として機能させるための地すべり面形状の推定プログラムである。
かかる請求項5記載の発明によれば、プログラムをコンピータに読み取らせて、地すべり面の形状を推定させることができるので、地すべり面の形状推定をコンピュータに実行させて、大規模で複雑なすべり面形状をより迅速に、そして適確に得ることができる。
本発明によれば、地すべり地表面の変位ベクトルを計測し、その結果から大規模で複雑なすべり面形状を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することができる。
以下、図面を参照して本発明の好適な実施の形態を例示的に詳しく説明する。但しこの実施の形態に記載されている構成部品の寸法、材質、形状、その相対的配置等は特に特定的な記載がない限りは、この発明の範囲をそれに限定する趣旨ではなく、単なる説明例に過ぎない。
参照する図面において、図1は本発明に係る地すべり面の推定方法が適用される地すべり面推定装置の全体構成を示すブロック図である。図2はすべり面座標に解析対象の断面地形図を示した説明図である。図3はブロック毎に分割したすべり面推定の説明図である。図4はすべり面推定手段における制御フローである。
地盤の地すべり面を推定する地すべり面推定装置1は、地表線、計測点、境界点等のデータ3の入力支援を行う断面モデル作成手段5、地表の変位測定手段6からの計測データ7を入力する計測データ入力手段9、前記断面モデル作成手段5および計測データ入力手段9からのデータに基づいてすべり面形状を推定するすべり面推定手段11、さらに、すべり面推定結果のすべり面形状を画面表示するすべり面形状表示部13を含んで構成される。
断面モデル作成手段5は、解析対象の地表面形状を作成する。地表面形状のデータを直接入力するか、または地表面形状が既にCADデータとして与えられている場合には該データを読み込む等によって、解析対象の断面地表形状のデータを入力するとともにデータベース15に記憶する。また、斜面に複数個所設置される変位測定手段6の設置位置の座標データや、後述する境界条件としての地表面に現れるすべり発生位置の頂上境界点データ、またはすべり終了位置の脚部境界点データを入力するとともにデータベース15に記憶する。
また、変位測定手段6の設置位置のデータが入力されることに伴って、該設置位置を1箇所含めて、地すべり土塊を水平方向に複数ブロックに分割するように鉛直方向に延びるブロック境界線(分割線)17の位置情報が設定される。このブロック境界線17の位置は、変位測定手段6の設置位置の中間点に自動的に設定されるようになっているが、作業者によって任意の位置座標を入力して設定することもできる。
計測データ入力手段9は、傾斜の地表面に設置された変位測定手段6からのデータが定期的に入力されるようにオンライン、または無線によって受信されるようになっている。なお、変位測定手段6はグローバルポジショニングシステム(GPS)によるもの、トータルステーション、レーザスキャナ、または変位計等によるものいずれの機器を用いても良い。
すべり面推定手段11は、前記断面モデル作成手段5で設定されたブロック境界線17によって分割されたブロック毎に、地すべり面を高次多項式で設定する高次多項式設定部21と、該ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも1つの既知の境界条件に基づいて求めるパラメータ算出部23と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を形成するすべり面推定部25とを有して構成されている。
高次多項式設定部21では、二次多項式を採用してすべり面形状を推定するように構成されている。
パラメータ算出部23では、すべり土塊がすべり面上を剛体として移動すると仮定しているので、地表計測点からすべり面までの深さDが一定に保たれるから、異なる計測時点の計測座標を用いて、すべり面座標の多項式を表すことができる。そして、計測点の座標は初期座標に計測変位を加えた測定値であるため誤差を伴っており、その誤差を、最小二乗法を適用して解析し、二次多項式のパラメータを求める。このとき境界条件として、断面モデル作成手段5によって入力設定された境界点データ、さらにブロック境界線17上の等高および等勾配条件を境界条件として算出される。このすべり面推定手段11の詳細手順については後述する。
パラメータ算出部23では、すべり土塊がすべり面上を剛体として移動すると仮定しているので、地表計測点からすべり面までの深さDが一定に保たれるから、異なる計測時点の計測座標を用いて、すべり面座標の多項式を表すことができる。そして、計測点の座標は初期座標に計測変位を加えた測定値であるため誤差を伴っており、その誤差を、最小二乗法を適用して解析し、二次多項式のパラメータを求める。このとき境界条件として、断面モデル作成手段5によって入力設定された境界点データ、さらにブロック境界線17上の等高および等勾配条件を境界条件として算出される。このすべり面推定手段11の詳細手順については後述する。
すべり面推定手段11のすべり面推定部25では、高次多項式設定部21で設定した二次多項式にパラメータ算出部23で求めたパラメータを用いて、ブロックごとの二次多項式求め、該求めた二次多項式を連ねて、ブロック全体のすべり面形状を推定する。
すべり面形状表示部13は、すべり面推定手段11によって推定されたすべり面形状を画面イメージとして出力し、断面モデル作成手段5で作成された解析対象の地表面形状とともにすべり面形状を画面表示部27に表示する。
次に、すべり面推定手段11における推定手順について、図2〜4を参照して説明する。
図4に示すフローチャートにおいて、ステップS1で開始すると、まずステップS2で断面モデル作成手段5によって作成された解析対象の断面モデルを入手する。次にステップS3では、断面モデル作成手段5によって設定されたブロック境界線17によって、地すべり土塊を水平方向に複数ブロックB1、B2、B3(図3参照)に分割する。各ブロックには変位測定手段6が1つ設置されている。そして、ステップS4で、ブロック毎のすべり面形状を2次多項式z=as2+bs+cによって設定する。
なお、原理的には2次式以外に次の(1)式のような一般的な高次多項式を用いることができるが、後述する最小二乗法を適用する場合の数値計算上の精度を考慮して、本実施形態では2次多項式を採用した。
図4に示すフローチャートにおいて、ステップS1で開始すると、まずステップS2で断面モデル作成手段5によって作成された解析対象の断面モデルを入手する。次にステップS3では、断面モデル作成手段5によって設定されたブロック境界線17によって、地すべり土塊を水平方向に複数ブロックB1、B2、B3(図3参照)に分割する。各ブロックには変位測定手段6が1つ設置されている。そして、ステップS4で、ブロック毎のすべり面形状を2次多項式z=as2+bs+cによって設定する。
なお、原理的には2次式以外に次の(1)式のような一般的な高次多項式を用いることができるが、後述する最小二乗法を適用する場合の数値計算上の精度を考慮して、本実施形態では2次多項式を採用した。
次に、ステップS5で異なる計測時点(t=t1、t2、t3)の計測点座標を用いて、すべり面座標の多項式を作成する。
本すべり面推定法の基本仮定によると、すべり土塊はすべり面上を剛体的に移動すると仮定しているから、地表変位計測点の変位軌跡は、地下のすべり面形状を表していることになる。図2に示すように、今、一つの変位測定手段6が設置されている地表変位計測点(座標値(si,zi))の変位軌跡に注目すると、変位を生じている期間中に地表計測点からすべり面までの深さDは一定に保たれるから、異なる計測時点(t=t1,t2,t3)の計測点座標を用いて、すべり面座標の多項式を次の(2)式のように表すことができる。なお、座標値(si,zi)は、水平方向位置をsi、鉛直方向位置をziとする。またnブロック中のi番目のブロック(測定点i)における2次多項式を示すものとする。
本すべり面推定法の基本仮定によると、すべり土塊はすべり面上を剛体的に移動すると仮定しているから、地表変位計測点の変位軌跡は、地下のすべり面形状を表していることになる。図2に示すように、今、一つの変位測定手段6が設置されている地表変位計測点(座標値(si,zi))の変位軌跡に注目すると、変位を生じている期間中に地表計測点からすべり面までの深さDは一定に保たれるから、異なる計測時点(t=t1,t2,t3)の計測点座標を用いて、すべり面座標の多項式を次の(2)式のように表すことができる。なお、座標値(si,zi)は、水平方向位置をsi、鉛直方向位置をziとする。またnブロック中のi番目のブロック(測定点i)における2次多項式を示すものとする。
次に、ステップS6で、式(2)の誤差式を作成する。
すべり面多項式のパラメータai、bi、ci、Diに関するこのような方程式は、計測点数×計測時点数(時刻履歴)だけ作ることができるが、図2に示すように計測点の座標値(si,zi)は初期座標に計測変位を加えた測定値であるため誤差を伴っている。前記(2)式の誤差をVi1,Vi2,Vi3とすると、これらはそれぞれ式(3)のように表される。
すべり面多項式のパラメータai、bi、ci、Diに関するこのような方程式は、計測点数×計測時点数(時刻履歴)だけ作ることができるが、図2に示すように計測点の座標値(si,zi)は初期座標に計測変位を加えた測定値であるため誤差を伴っている。前記(2)式の誤差をVi1,Vi2,Vi3とすると、これらはそれぞれ式(3)のように表される。
次に、ステップS7で、すべり面多項式のパラメータai、bi、ci、Diを求めるには、最小二乗法を適用して、これら誤差の二乗和を最小化する解を求める。すなわち、次の式(4)、式(5)の解を求めることになる。
Vは、a、b、c、Dの4変数の2次関数であり、しかもa2 ,b2,c2,D2の係数は正であるからVが、a、b、c、Dのそれぞれに関して最小値をもつ。
Vは、a、b、c、Dの4変数の2次関数であり、しかもa2 ,b2,c2,D2の係数は正であるからVが、a、b、c、Dのそれぞれに関して最小値をもつ。
一般に、n個の計測点に対して、m個の計測時点座標(m−1個の変位時刻歴)が存在する場合、n×m個の誤差は次の式(6)と表すことができる。
式(6)を{V}=[S]{A}−[Z]と表すと最小二乗法では、次の式(7)を行うことになるから、たとえばaの偏微分式からは式(8)の方程式が導かれる。
次に、図4(a)のフローチャートから図4(b)のフローチャートへ進み、ステップS8で、前記した1計測点の3時点(t=t1、t2、t3)の計測点座標(Si,Zi)が得られているケースにあてはめると、次の式(10)なる方程式を得る。
式(10)は明らかに第3式と第4式の方程式が同一であり、この方程式は不定(係数マトリックス[S]が特異)となる。
このことは、計測点数や計測時点数が増加しても解消されることはなく、常に第3式は第4式以降と一次従属関係にある。すなわち、第4式以降の式の和が第3式と同一になる。
そこで、本発明では、ステップS9で示すように全体方程式の行列式が特異にならないように、各ブロックの正規方程式の第4式を取り除き、その代わりに、次の式(11)〜(14)に示す新たな境界条件による未知パラメータ方程式を組み込む。
このことは、計測点数や計測時点数が増加しても解消されることはなく、常に第3式は第4式以降と一次従属関係にある。すなわち、第4式以降の式の和が第3式と同一になる。
そこで、本発明では、ステップS9で示すように全体方程式の行列式が特異にならないように、各ブロックの正規方程式の第4式を取り除き、その代わりに、次の式(11)〜(14)に示す新たな境界条件による未知パラメータ方程式を組み込む。
式(12)、(13)は端部ブロック場合の条件を示す。
図3に示すように第1ブロックB1の場合には、脚部境界点(SH、ZH)の条件、また頂部に位置するブロックB3の場合には、頂部境界点(ST、ZT)の条件を示している。
図3に示すように第1ブロックB1の場合には、脚部境界点(SH、ZH)の条件、また頂部に位置するブロックB3の場合には、頂部境界点(ST、ZT)の条件を示している。
式(14)は、ブロック境界線17上における隣接多項式の等勾配条件を設定している。この条件は、ブロック境界線17上での等高さ条件と同様に、ブロック間連成方程式となる。
次に、ステップS10で、式(10)の第1式〜第3式の3つの式と、式(11)、(12)、(13)、(14)の式によって、iブロックにおける多項式パラメータ(ai、bi、ci、Di)を算出する。
なお、式(11)、(14)についてはブロック間の連成方程式であるため隣接するブロックの方程式が必要であり、また、式(12)、(13)については端部ブロックの地表境界点のデータが必要となる。このため結果的には全ブロックにおける未知パラメータ方程式を行列化して、nブロック全てのパラメータについて(ai、bi、ci、Di)×n(ブロック数)を求める必要がある。
また、式(14)の等勾配条件式については、特に用いなくても未知パラメータの解を求めることができるが、式(14)があることによって、少ない計測時点であっても適切に解を求めることができるようになる。
なお、式(11)、(14)についてはブロック間の連成方程式であるため隣接するブロックの方程式が必要であり、また、式(12)、(13)については端部ブロックの地表境界点のデータが必要となる。このため結果的には全ブロックにおける未知パラメータ方程式を行列化して、nブロック全てのパラメータについて(ai、bi、ci、Di)×n(ブロック数)を求める必要がある。
また、式(14)の等勾配条件式については、特に用いなくても未知パラメータの解を求めることができるが、式(14)があることによって、少ない計測時点であっても適切に解を求めることができるようになる。
さらに、次のステップS11で、求められたブロック毎の2次多項式を連ねて全体ブロックの全体方程式としての2次多項式を作成する。
すなわち、第1ブロックB1の領域では、第1ブロックB1の多項式パラメータ(a1、b1、c1、D1)を用いた2次多項式が用いられ、第2ブロックB2の領域では、第2ブロックB2の多項式パラメータ(a2、b2、c2、D2)を用いた2次多項式が用いられ、第3ブロックB3の領域では、第3ブロックB3の多項式パラメータ(a3、b3、c3、D3)を用いた2次多項式が用いられて、全体として2次多項式によってすべり面形状が推定される。そしてステップS11で終了する。
すなわち、第1ブロックB1の領域では、第1ブロックB1の多項式パラメータ(a1、b1、c1、D1)を用いた2次多項式が用いられ、第2ブロックB2の領域では、第2ブロックB2の多項式パラメータ(a2、b2、c2、D2)を用いた2次多項式が用いられ、第3ブロックB3の領域では、第3ブロックB3の多項式パラメータ(a3、b3、c3、D3)を用いた2次多項式が用いられて、全体として2次多項式によってすべり面形状が推定される。そしてステップS11で終了する。
以上のようにブロック間の連続条件であるブロック境界線上の等高および等勾配条件を未知パラメータ方程式に組み込むことによって、すなわち、式(11)、(12)、(13)、(14)を組み込むことによって、ブロック毎の最小二乗法の適用によって得られる式(10)の第1式〜第3式だけではパラメータを求めることができなかった場合でも、2次多項式のパラメータを算出できるようになる。
また、境界条件として、地表面表れている脚部、頂部の境界点データでなく、予め地中調査によって得られている地中境界データがあれば、それを未知パラメータ方程式に組み入れても良い。
さらに、本実施形態の例では、1計測点の3時点(t=t1、t2、t3)の計測点座標(Si,Zi)について説明したが、ブロック内で最低限1個の変位ベクトルが得られていれば解を得ることができる。
また、境界条件として、地表面表れている脚部、頂部の境界点データでなく、予め地中調査によって得られている地中境界データがあれば、それを未知パラメータ方程式に組み入れても良い。
さらに、本実施形態の例では、1計測点の3時点(t=t1、t2、t3)の計測点座標(Si,Zi)について説明したが、ブロック内で最低限1個の変位ベクトルが得られていれば解を得ることができる。
以上本実施の形態によれば、地すべり土塊を、地表の変位測定手段6を含む領域毎に複数のブロックB1、B2、B3に分割して、ブロック毎に2次多項式ですべり面の形状を設定するため、1つの多項式で全体のすべり面形状を設定するのに比べて、複数の多項式で推定するため大規模で複雑なすべり面形状をより迅速に適確に推定することができる。
また、変位測定手段6によって地すべり地表面変位ベクトルを計測し、その結果からすべり面形状を推定するため、地すべりの発生初期において、変状の規模を迅速に把握して、すべり機構の迅速な把握とその後の対策工の実施に役立つことができる。
また、地すべり機構把握のための通常のボーリング調査等では一般に日数を要し、しかも地すべり地内での作業であるため安全確保の点からも慎重な対応が求められていたが、本発明により地すべり地内での危険な作業を減らすことが可能になる。
また、本実施形態によれば、前記のすべり面推定手段11を備えた地すべり面推定装置1によれば、ブロック毎に2次多項式ですべり面の形状を設定して大規模で複雑なすべり面形状をより適確に推定することができる装置を簡単な構成で達成することができる。
さらに、本実施形態によれば、プログラムをコンピータに読み取らせて、地すべり面の形状を推定させることができるので、地すべり面の形状推定をコンピュータに実行させて、大規模で複雑なすべり面形状をより迅速に、そして適確に得ることができる。
本発明によれば、地すべり地表面の変位ベクトルを計測し、その結果から大規模で複雑なすべり面形状を精度よく、かつ迅速に把握できる実用的な推定手法、推定装置、および推定プログラムを提供することができるので、地すべり面形状を推定する装置への適用に際して有益である。
1 地すべり面推定装置
5 断面モデル作成手段
6 変位測定手段
9 計測データ入力部
11 すべり面推定手段
13 すべり面形状表示部
15 データベース
17 ブロック境界線(分割線)
21 高次多項式設定部
23 パラメータ算出部
25 すべり面推定部管路部品
B1 第1ブロック
B2 第2ブロック
B3 第3ブロック
5 断面モデル作成手段
6 変位測定手段
9 計測データ入力部
11 すべり面推定手段
13 すべり面形状表示部
15 データベース
17 ブロック境界線(分割線)
21 高次多項式設定部
23 パラメータ算出部
25 すべり面推定部管路部品
B1 第1ブロック
B2 第2ブロック
B3 第3ブロック
Claims (5)
- 地表の変位測定値から2次元鉛直断面における地すべり面形状を推定する推定方法において、
地すべり土塊を鉛直方向の分割線によって複数のブロックに分割するとともに、変位測定手段が各ブロックの地表に設置されるように分割し、該ブロック毎の地すべり面を高次多項式で設定し、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも1つの既知の境界条件に基づいて求め、算出されたブロック毎の前記高次多項式を連ねて地すべり土塊全体のすべり面形状を推定することを特徴とする地すべり面形状の推定方法。 - 前記境界条件が分割線上の等高および等勾配条件であり、端部ブロックではさらに地表境界条件を含むことを特徴とする請求項1記載の地すべり面形状の推定方法。
- 地表の変位測定値から2次元鉛直断面における地すべり面形状を推定する推定装置において、
地表線、地表変位計測点、境界点、および分割線位置のデータを編集して解析断面形状を作成する断面モデル作成手段と、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段と、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段とを有し、前記すべり面推定手段は、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも1つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部とを備えたことを特徴とする地すべり面形状の推定装置。 - 前記パラメータ算出部では、前記境界条件を分割線上の等高および等勾配条件から、さらに端部ブロックでは地表境界条件から求めることを特徴とする請求項3記載の地すべり面形状の推定装置。
- 地盤の地すべり面形状を推定するためのコンピュータを、
地表線、地表変位計測点、境界点、および分割線位置のデータを編集して鉛直断面形状を作成する断面モデル作成手段、地表に設置された変位測定手段からの計測データを入力する計測データ入力手段、前記断面モデル作成手段および計測データ入力手段からのデータに基づいてすべり面形状を推定するすべり面推定手段として機能させ、さらに前記すべり面推定手段を、鉛直方向の分割線によって分割されたブロック毎に地すべり面を高次多項式で設定する高次多項式設定部と、ブロック毎に最小二乗法による正規方程式を組立てて前記ブロック毎の高次多項式のパラメータを少なくとも1つの既知の境界条件に基づいて求めるパラメータ算出部と、算出された夫々のブロックの高次多項式を連ねてすべり面形状を推定するすべり面推定部として機能させるための地すべり面形状の推定プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007037150A JP2008202994A (ja) | 2007-02-16 | 2007-02-16 | 地すべり面形状の推定方法、推定装置および推定プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007037150A JP2008202994A (ja) | 2007-02-16 | 2007-02-16 | 地すべり面形状の推定方法、推定装置および推定プログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2008202994A true JP2008202994A (ja) | 2008-09-04 |
Family
ID=39780681
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007037150A Pending JP2008202994A (ja) | 2007-02-16 | 2007-02-16 | 地すべり面形状の推定方法、推定装置および推定プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2008202994A (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018106596A (ja) * | 2016-12-28 | 2018-07-05 | ラピスセミコンダクタ株式会社 | 地表移動検出装置、無線タグ、地表移動検出方法及び災害救助支援システム |
CN111210502A (zh) * | 2019-11-25 | 2020-05-29 | 温州大学 | 一种预测海底管线遭受滑坡竖向冲击作用力的计算机方法 |
JP2020165746A (ja) * | 2019-03-29 | 2020-10-08 | 国際航業株式会社 | 地すべり面推定装置、及び地すべり面推定方法 |
CN112669572A (zh) * | 2020-12-17 | 2021-04-16 | 四方智能(武汉)控制技术有限公司 | 一种用于江河流域库岸智能巡查的无人船系统 |
CN114495434A (zh) * | 2022-02-08 | 2022-05-13 | 北京寒武智能科技有限公司 | 一种滑坡灾害临滑预测预报的方法 |
CN114925577A (zh) * | 2022-06-07 | 2022-08-19 | 中国地质科学院地质力学研究所 | 一种蠕动型滑坡体上砌体建筑物的易损性分析方法与系统 |
CN116399211A (zh) * | 2023-06-08 | 2023-07-07 | 湖南联智科技股份有限公司 | 一种基于两点北斗监测定位计算边坡面滑移的方法 |
-
2007
- 2007-02-16 JP JP2007037150A patent/JP2008202994A/ja active Pending
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018106596A (ja) * | 2016-12-28 | 2018-07-05 | ラピスセミコンダクタ株式会社 | 地表移動検出装置、無線タグ、地表移動検出方法及び災害救助支援システム |
JP2020165746A (ja) * | 2019-03-29 | 2020-10-08 | 国際航業株式会社 | 地すべり面推定装置、及び地すべり面推定方法 |
JP7170575B2 (ja) | 2019-03-29 | 2022-11-14 | 国際航業株式会社 | 地すべり面推定装置、及び地すべり面推定方法 |
CN111210502A (zh) * | 2019-11-25 | 2020-05-29 | 温州大学 | 一种预测海底管线遭受滑坡竖向冲击作用力的计算机方法 |
CN111210502B (zh) * | 2019-11-25 | 2023-11-07 | 温州大学 | 一种预测海底管线遭受滑坡竖向冲击作用力的计算方法 |
CN112669572A (zh) * | 2020-12-17 | 2021-04-16 | 四方智能(武汉)控制技术有限公司 | 一种用于江河流域库岸智能巡查的无人船系统 |
CN112669572B (zh) * | 2020-12-17 | 2023-07-04 | 四方智能(武汉)控制技术有限公司 | 一种用于江河流域库岸智能巡查的无人船系统 |
CN114495434A (zh) * | 2022-02-08 | 2022-05-13 | 北京寒武智能科技有限公司 | 一种滑坡灾害临滑预测预报的方法 |
CN114495434B (zh) * | 2022-02-08 | 2024-01-12 | 北京寒武智能科技有限公司 | 一种滑坡灾害临滑预测预报的方法 |
CN114925577A (zh) * | 2022-06-07 | 2022-08-19 | 中国地质科学院地质力学研究所 | 一种蠕动型滑坡体上砌体建筑物的易损性分析方法与系统 |
CN116399211A (zh) * | 2023-06-08 | 2023-07-07 | 湖南联智科技股份有限公司 | 一种基于两点北斗监测定位计算边坡面滑移的方法 |
CN116399211B (zh) * | 2023-06-08 | 2023-08-18 | 湖南联智科技股份有限公司 | 一种基于两点北斗监测定位计算边坡面滑移的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2008202994A (ja) | 地すべり面形状の推定方法、推定装置および推定プログラム | |
Argüelles-Fraga et al. | Measurement planning for circular cross-section tunnels using terrestrial laser scanning | |
JP6591131B1 (ja) | 構造物計測装置および構造物計測方法 | |
Jeong et al. | Applying unmanned aerial vehicle photogrammetry for measuring dimension of structural elements in traditional timber building | |
Riveiro et al. | A novel approach to evaluate masonry arch stability on the basis of limit analysis theory and non-destructive geometric characterization | |
JP6506215B2 (ja) | 地質境界面または断層面の予測方法 | |
JP6979364B2 (ja) | トンネル経路の地表面変位観測装置及びトンネル経路の地表面変位観測プログラム | |
RU2013129852A (ru) | Оценивание поверхностных данных | |
Górska et al. | Displacement monitoring and sensitivity analysis in the observational method | |
Tanoli et al. | 3D design modeling application in machine guidance system for earthwork operations | |
JP3867025B2 (ja) | トンネル管理図及びその作成方法並びにそのシステム | |
CN116663762A (zh) | 一种城市规划地下空间勘察测绘方法及系统 | |
JP7244336B2 (ja) | 山岳トンネルコンクリート厚測定方法および測定装置 | |
Fys et al. | Modeling of 3-D objects using geodetic and cartographic data and determining their volumes with an accuracy assessment | |
JP2006349579A (ja) | 変位情報測定システム及び方法並びにプログラム | |
Yang et al. | Laser-based intelligent perception method for tunnel invasion | |
Altan et al. | Photogrammetry and geographic information systems for quick assessment, documentation and analysis of earthquakes | |
JP3691454B2 (ja) | 地盤変位測位システム | |
JP2018132993A (ja) | 建物設計情報修正支援装置、建物設計情報修正支援方法、及びプログラム | |
Ekseth et al. | Improving the Quality of Ellipse of Uncertainty Calculations in Gyro Surveys to Reduce the Risk of Hazardous Events like Blowouts or Missing Potential Production through Incorrect Wellbore Placement | |
PRATALI MAFFEI et al. | The real in the virtual. The 3D model in the cultural heritage sector: The tip of the iceberg | |
Mahasuwanchai et al. | An Alternative Method for Long‐Term Monitoring of Thai Historic Pagodas Based on Terrestrial Laser Scanning Data: A Case Study of Wat Krachee in Ayutthaya | |
JP3867011B2 (ja) | トンネル現況展開図及びその作成方法 | |
Rodriguez-Gonzalvez et al. | Learning physical geodesy. Application case to geoid undulation computation | |
Klapa et al. | Inventory of various stages of construction using TLS technology |