JP2001193046A - 地盤速度構造の推定方法及び推定システム - Google Patents
地盤速度構造の推定方法及び推定システムInfo
- Publication number
- JP2001193046A JP2001193046A JP2000003589A JP2000003589A JP2001193046A JP 2001193046 A JP2001193046 A JP 2001193046A JP 2000003589 A JP2000003589 A JP 2000003589A JP 2000003589 A JP2000003589 A JP 2000003589A JP 2001193046 A JP2001193046 A JP 2001193046A
- Authority
- JP
- Japan
- Prior art keywords
- ground
- theoretical
- velocity
- dominant
- destination
- 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.)
- Granted
Links
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
(57)【要約】
【課題】 短時間に、正確に地盤速度構造を推定する方
法、及び推定システムを提供する。 【解決手段】 調査目的地にて複数の振動センサを有す
る微動観測装置により地盤の微動探査を行い、その分散
特性から参考速度モデルを初期値にして逆解析により地
盤速度構造を推定する方法であり、前記調査目的地周辺
の複数の地盤データから複数の理論卓越周期を算出し、
前記調査目的地にて計測した観測卓越周期と前記複数の
理論卓越周期とを比較して最も近似している理論卓越周
期を有する地盤データを選択し、この選択された地盤デ
ータを参考速度モデルとして逆解析を行い、前記調査目
的地の地盤速度構造を推定することを特徴とする。
法、及び推定システムを提供する。 【解決手段】 調査目的地にて複数の振動センサを有す
る微動観測装置により地盤の微動探査を行い、その分散
特性から参考速度モデルを初期値にして逆解析により地
盤速度構造を推定する方法であり、前記調査目的地周辺
の複数の地盤データから複数の理論卓越周期を算出し、
前記調査目的地にて計測した観測卓越周期と前記複数の
理論卓越周期とを比較して最も近似している理論卓越周
期を有する地盤データを選択し、この選択された地盤デ
ータを参考速度モデルとして逆解析を行い、前記調査目
的地の地盤速度構造を推定することを特徴とする。
Description
【０００１】
【発明の属する技術分野】本発明は、レイリー波の分散
特性から地盤の速度構造を求める逆解析法（インバージ
ョン法）を用いて、調査目的地の地盤の速度構造を推定
する方法および推定システムに関する。
特性から地盤の速度構造を求める逆解析法（インバージ
ョン法）を用いて、調査目的地の地盤の速度構造を推定
する方法および推定システムに関する。
【０００２】
【従来の技術】地下構造を調査する方法としては、ボー
リング孔を用いる方法や、反射法のような地震探査法が
知られている。また、これらの方法に加えて、最近では
レイリー波、ラブ波等の表面波を利用する方法が提案さ
れ、この方法では、比較的深層部を調査する場合には常
時微動を用い、宅地等の浅層部を調査する場合には起震
機等を用いてその上下成分であるレイリー波を観測し、
周波数−波数スペクトル法や空間自己相関法により周波
数と位相速度との関係（以下、観測分散特性という）を
得て、地下の地盤速度構造から計算される理論上の周波
数と位相速度との関係（以下、理論分散特性という）
と、前記の観測分散特性とを比較して地盤モデルの修正
を繰り返すことにより、地下の地盤速度構造を推定する
インバージョン法が用いられることが多い。
リング孔を用いる方法や、反射法のような地震探査法が
知られている。また、これらの方法に加えて、最近では
レイリー波、ラブ波等の表面波を利用する方法が提案さ
れ、この方法では、比較的深層部を調査する場合には常
時微動を用い、宅地等の浅層部を調査する場合には起震
機等を用いてその上下成分であるレイリー波を観測し、
周波数−波数スペクトル法や空間自己相関法により周波
数と位相速度との関係（以下、観測分散特性という）を
得て、地下の地盤速度構造から計算される理論上の周波
数と位相速度との関係（以下、理論分散特性という）
と、前記の観測分散特性とを比較して地盤モデルの修正
を繰り返すことにより、地下の地盤速度構造を推定する
インバージョン法が用いられることが多い。
【０００３】インバージョン法には最小２乗法がよく用
いられるが、最近では遺伝アルゴリズムを用いる方法が
知られている。このインバージョン法に用いられるパラ
メータは、Ｐ波速度、Ｓ波速度、密度、層厚の４種類が
あるが、理論分散特性に大きく影響を与えるのはＳ波速
度と層厚であるので、この２種類をパラメータとして計
算させることが一般的である。
いられるが、最近では遺伝アルゴリズムを用いる方法が
知られている。このインバージョン法に用いられるパラ
メータは、Ｐ波速度、Ｓ波速度、密度、層厚の４種類が
あるが、理論分散特性に大きく影響を与えるのはＳ波速
度と層厚であるので、この２種類をパラメータとして計
算させることが一般的である。
【０００４】また、インバージョン法では観測分散特性
と理論分散特性を比較する必要があり、その一例が「物
理探査」第５０巻、第２号（１９９７）９３−１０６頁
の「短周期微動のアレイ観測による盛岡市域の地価速度
構造の推定」に記載されている。この従来例では、非線
形最小２乗法により地下の地盤速度構造を推定してお
り、インバージョン法の初期モデル（初期地下速度構造
モデル）を設定する場合に、観測分散特性に基づいて表
層のＳ波速度の概略と下層の速度のおおよその範囲を推
定した上で、既存の地図等から試行錯誤で設定してい
る。
と理論分散特性を比較する必要があり、その一例が「物
理探査」第５０巻、第２号（１９９７）９３−１０６頁
の「短周期微動のアレイ観測による盛岡市域の地価速度
構造の推定」に記載されている。この従来例では、非線
形最小２乗法により地下の地盤速度構造を推定してお
り、インバージョン法の初期モデル（初期地下速度構造
モデル）を設定する場合に、観測分散特性に基づいて表
層のＳ波速度の概略と下層の速度のおおよその範囲を推
定した上で、既存の地図等から試行錯誤で設定してい
る。
【０００５】
【発明が解決しようとする課題】しかしながら、物理探
査学会第９5回学術講演会論文集（１９９６）７６−８
０頁の「平滑化制約付き非線形最小２乗法によるラブ波
位相速度のインバージョン」にも記載されているよう
に、非線形最小２乗法を用いるインバージョン法は、方
程式が悪循環になり、解の安定性と一意性の問題が発生
する。したがって前記した従来の方法では、初期モデル
（初期地下速度構造モデル）を試行錯誤で設定するの
で、観測分散特性と理論分散特性が略一致して得られる
解が必ずしも正しいとは限らず、場合によっては結果を
得るまで何回もトライする必要があり、このため多大な
時間を要するという問題がある。また、既存の地図等か
ら情報を得るにしても、具体的にどの情報を用いるかは
不明であり、特に調査地が複雑な地形内にあるときはど
の情報を採用するかが不明確であり、採用した情報によ
り結果が異なることもあった。
査学会第９5回学術講演会論文集（１９９６）７６−８
０頁の「平滑化制約付き非線形最小２乗法によるラブ波
位相速度のインバージョン」にも記載されているよう
に、非線形最小２乗法を用いるインバージョン法は、方
程式が悪循環になり、解の安定性と一意性の問題が発生
する。したがって前記した従来の方法では、初期モデル
（初期地下速度構造モデル）を試行錯誤で設定するの
で、観測分散特性と理論分散特性が略一致して得られる
解が必ずしも正しいとは限らず、場合によっては結果を
得るまで何回もトライする必要があり、このため多大な
時間を要するという問題がある。また、既存の地図等か
ら情報を得るにしても、具体的にどの情報を用いるかは
不明であり、特に調査地が複雑な地形内にあるときはど
の情報を採用するかが不明確であり、採用した情報によ
り結果が異なることもあった。
【０００６】本発明は、このような問題に鑑みてなされ
たものであって、その目的とするところは、短時間で、
正確に地盤の速度構造を推定できる方法を提供すること
にある。また、前記した方法により短時間で、正確に地
盤の速度構造を推定できるシステムを提供することを目
的とする。
たものであって、その目的とするところは、短時間で、
正確に地盤の速度構造を推定できる方法を提供すること
にある。また、前記した方法により短時間で、正確に地
盤の速度構造を推定できるシステムを提供することを目
的とする。
【０００７】
【課題を解決するための手段】前記目的を達成すべく、
本発明に係る地盤速度構造の推定方法は、調査目的地に
て複数の振動センサを有する微動観測装置により地盤の
微動探査を行い、その分散特性から参考速度モデルを初
期値にして逆解析により地盤速度構造を推定する方法で
あり、前記調査目的地周辺の複数の地盤データから複数
の理論卓越周期を算出し、前記調査目的地にて計測した
観測卓越周期と前記複数の理論卓越周期とを比較して最
も近似している理論卓越周期を有する地盤データを選択
し、この選択された地盤データを参考速度モデルとして
逆解析を行い、前記調査目的地の地盤速度構造を推定す
ることを特徴とする。この構成によれば、地盤データベ
ースから算出した複数の理論卓越周期のうち、観測卓越
周期と最も近似している理論卓越周期を有する地盤デー
タを選んで、この地盤データを初期値として逆解析を行
なうため、短時間で正確に地盤の速度構造を推定するこ
とができる。
本発明に係る地盤速度構造の推定方法は、調査目的地に
て複数の振動センサを有する微動観測装置により地盤の
微動探査を行い、その分散特性から参考速度モデルを初
期値にして逆解析により地盤速度構造を推定する方法で
あり、前記調査目的地周辺の複数の地盤データから複数
の理論卓越周期を算出し、前記調査目的地にて計測した
観測卓越周期と前記複数の理論卓越周期とを比較して最
も近似している理論卓越周期を有する地盤データを選択
し、この選択された地盤データを参考速度モデルとして
逆解析を行い、前記調査目的地の地盤速度構造を推定す
ることを特徴とする。この構成によれば、地盤データベ
ースから算出した複数の理論卓越周期のうち、観測卓越
周期と最も近似している理論卓越周期を有する地盤デー
タを選んで、この地盤データを初期値として逆解析を行
なうため、短時間で正確に地盤の速度構造を推定するこ
とができる。
【０００８】また、本発明に係る地盤速度構造の推定方
法の好ましい具体的な態様としては、前記観測卓越周期
及び理論卓越周期の範囲は、０．１秒から１秒の範囲で
あることを特徴としている。この構成によれば、前記の
範囲は実際に建築物の耐震設計をする場合に、S波速度
で４００ｍ／ｓを越える工学的基盤に相当するため、そ
の範囲で効率的に比較することができ、短時間で正確に
地盤速度構造を推定することができる。
法の好ましい具体的な態様としては、前記観測卓越周期
及び理論卓越周期の範囲は、０．１秒から１秒の範囲で
あることを特徴としている。この構成によれば、前記の
範囲は実際に建築物の耐震設計をする場合に、S波速度
で４００ｍ／ｓを越える工学的基盤に相当するため、そ
の範囲で効率的に比較することができ、短時間で正確に
地盤速度構造を推定することができる。
【０００９】本発明に係る地盤速度構造推定システム
は、調査目的地にて微動探査を行なう複数の振動センサ
と、該複数の振動センサが接続された解析装置と、該解
析装置に入力する入力装置と、前記解析装置に接続され
た出力装置と、前記解析装置に接続された地盤データベ
ースとを備え、前記解析装置は、前記振動センサからの
微動データから観測分散特性を算出する分散特性算出手
段と、前記入力装置により入力された複数の地点におけ
る地盤データベースの地盤データからそれぞれの地点に
おける複数の理論卓越周期を算出する卓越周期算出手段
と、前記調査目的地の観測卓越周期を計測する計測手段
と、逆解析により地盤の速度構造を算出して前記出力装
置へ出力する速度構造算出手段とを有し、前記複数の理
論卓越周期と前記調査目的地の観測卓越周期とを比較
し、最も近似している理論卓越周期を有する地盤データ
を選択して前記速度構造算出手段へ参考速度モデルとし
て入力し、前記速度構造算出手段において前記参考速度
モデルを初期値にして逆解析を行い、前記調査目的地の
地盤速度構造を算出することを特徴としている。この構
成によれば、前記した推定方法により、調査目的地の地
盤速度構造を短時間に、正確に推定することができる。
は、調査目的地にて微動探査を行なう複数の振動センサ
と、該複数の振動センサが接続された解析装置と、該解
析装置に入力する入力装置と、前記解析装置に接続され
た出力装置と、前記解析装置に接続された地盤データベ
ースとを備え、前記解析装置は、前記振動センサからの
微動データから観測分散特性を算出する分散特性算出手
段と、前記入力装置により入力された複数の地点におけ
る地盤データベースの地盤データからそれぞれの地点に
おける複数の理論卓越周期を算出する卓越周期算出手段
と、前記調査目的地の観測卓越周期を計測する計測手段
と、逆解析により地盤の速度構造を算出して前記出力装
置へ出力する速度構造算出手段とを有し、前記複数の理
論卓越周期と前記調査目的地の観測卓越周期とを比較
し、最も近似している理論卓越周期を有する地盤データ
を選択して前記速度構造算出手段へ参考速度モデルとし
て入力し、前記速度構造算出手段において前記参考速度
モデルを初期値にして逆解析を行い、前記調査目的地の
地盤速度構造を算出することを特徴としている。この構
成によれば、前記した推定方法により、調査目的地の地
盤速度構造を短時間に、正確に推定することができる。
【００１０】
【発明の実施の形態】以下、本発明に係る地盤速度構造
推定システムの一実施形態を図面に基づき詳細に説明す
る。図１は、本発明に係る地盤速度構造の推定システム
の一実施形態のブロック図、図２は調査目的地Ｐと、周
辺の地盤データベースから選定した３地点Ａ、Ｂ、Ｃの
位置関係を示す説明図である。図１〜２において、微動
観測により地盤振動探査を行なう推定システム１０につ
いて説明する。調査目的地Ｐ上において、４個の振動セ
ンサ１１〜１４を用意し、これらの振動センサ１１、１
２、１３が半径Ｒの円周上の３点に等間隔に設置される
と共に、振動センサ１４が円周の中心に設置され、セン
サアレイ１５を構成している。複数の振動センサは調査
対象地盤の常時微動を検出するために配置され、表面波
の鉛直成分を同時に検出するものである。センサアレイ
１５のサイズは、調査地盤、求めたい深さにより決定さ
れるが、ここでは半径Rを３ｍから５０ｍまでの範囲で
適宜に変えて計測する。なお、常時微動とは１Hｚ〜数
１０Hｚまでのレイリー波である。
推定システムの一実施形態を図面に基づき詳細に説明す
る。図１は、本発明に係る地盤速度構造の推定システム
の一実施形態のブロック図、図２は調査目的地Ｐと、周
辺の地盤データベースから選定した３地点Ａ、Ｂ、Ｃの
位置関係を示す説明図である。図１〜２において、微動
観測により地盤振動探査を行なう推定システム１０につ
いて説明する。調査目的地Ｐ上において、４個の振動セ
ンサ１１〜１４を用意し、これらの振動センサ１１、１
２、１３が半径Ｒの円周上の３点に等間隔に設置される
と共に、振動センサ１４が円周の中心に設置され、セン
サアレイ１５を構成している。複数の振動センサは調査
対象地盤の常時微動を検出するために配置され、表面波
の鉛直成分を同時に検出するものである。センサアレイ
１５のサイズは、調査地盤、求めたい深さにより決定さ
れるが、ここでは半径Rを３ｍから５０ｍまでの範囲で
適宜に変えて計測する。なお、常時微動とは１Hｚ〜数
１０Hｚまでのレイリー波である。
【００１１】４個の振動センサ１１〜１４により検出さ
れた複数箇所の常時微動データは解析装置２０に入力さ
れ、解析装置２０にはデータ等を入力する入力装置２１
が接続され、解析装置２０からの解析結果を出力する出
力装置２２が接続されている。解析装置２０には地盤デ
ータベース２３が接続され、地盤データを取り込むこと
ができる。微動観測装置は、複数の振動センサ１１〜１
４から構成されるセンサアレイ１５、その出力を解析す
る解析装置２０、入力装置２１及び出力装置２２から構
成される。
れた複数箇所の常時微動データは解析装置２０に入力さ
れ、解析装置２０にはデータ等を入力する入力装置２１
が接続され、解析装置２０からの解析結果を出力する出
力装置２２が接続されている。解析装置２０には地盤デ
ータベース２３が接続され、地盤データを取り込むこと
ができる。微動観測装置は、複数の振動センサ１１〜１
４から構成されるセンサアレイ１５、その出力を解析す
る解析装置２０、入力装置２１及び出力装置２２から構
成される。
【００１２】解析装置２０は、振動センサ１１〜１４か
らの微動データよりレイリー波の分散特性を算出する分
散特性算出手段３１と、入力装置２１により入力された
地盤データベース２３の複数の地盤データから、それぞ
れの理論卓越周期を算出する卓越周期算出手段３２、調
査目的地における観測卓越周期を計測する計測手段３
３、逆解析により地盤の速度構造を算出して出力装置２
２へ出力する速度構造算出手段３４とを有している。そ
して、複数の理論卓越周期と調査目的地の観測卓越周期
とを比較し、最も近似している理論卓越周期を有する地
盤データを選択して入力装置２１により参考速度モデル
として入力し、速度構造算出手段３４により選択された
地盤データを初期値として逆解析を行い、調査目的地Ｐ
の地盤速度構造を算出して推定し、出力装置２２より出
力する。なお、卓越周期とは、地盤に特有の周波数特性
の内、ピークの周期のことをいう。
らの微動データよりレイリー波の分散特性を算出する分
散特性算出手段３１と、入力装置２１により入力された
地盤データベース２３の複数の地盤データから、それぞ
れの理論卓越周期を算出する卓越周期算出手段３２、調
査目的地における観測卓越周期を計測する計測手段３
３、逆解析により地盤の速度構造を算出して出力装置２
２へ出力する速度構造算出手段３４とを有している。そ
して、複数の理論卓越周期と調査目的地の観測卓越周期
とを比較し、最も近似している理論卓越周期を有する地
盤データを選択して入力装置２１により参考速度モデル
として入力し、速度構造算出手段３４により選択された
地盤データを初期値として逆解析を行い、調査目的地Ｐ
の地盤速度構造を算出して推定し、出力装置２２より出
力する。なお、卓越周期とは、地盤に特有の周波数特性
の内、ピークの周期のことをいう。
【００１３】振動センサ１１〜１４は、複数個用いて常
時微動の表面波の鉛直成分であるレイリー波を同時に検
出するものであり、センサアレイ１５の半径Ｒは調査対
象の地盤、求めたい深さにより決定される。このとき１
つのセンサアレイ１５で得られる分散特性は周波数範囲
が狭いため、センサアレイ１５の分散特性を合成して幅
広い周波数範囲で分散特性を作成する。実際には、セン
サアレイ１５の半径を変えて、順次計測を行なう。
時微動の表面波の鉛直成分であるレイリー波を同時に検
出するものであり、センサアレイ１５の半径Ｒは調査対
象の地盤、求めたい深さにより決定される。このとき１
つのセンサアレイ１５で得られる分散特性は周波数範囲
が狭いため、センサアレイ１５の分散特性を合成して幅
広い周波数範囲で分散特性を作成する。実際には、セン
サアレイ１５の半径を変えて、順次計測を行なう。
【００１４】解析装置２０は地盤の速度構造を逆解析法
（インバージョン法）により算出する前記の速度構造算
出手段３４を備え、入力装置２１を介して解析が行われ
る。この速度構造算出手段３４は、地盤データベース２
３などを参考にして初期値のデータベースモデル（参考
速度モデル）を決定し、解析装置２０の分散特性算出手
段３１にて微動観測により求めた観測分散特性と合致す
るように参考速度モデルを修正する逆解析を行い、当該
地盤の速度構造を算出して推定し、出力装置２２に結果
を出力するものである。
（インバージョン法）により算出する前記の速度構造算
出手段３４を備え、入力装置２１を介して解析が行われ
る。この速度構造算出手段３４は、地盤データベース２
３などを参考にして初期値のデータベースモデル（参考
速度モデル）を決定し、解析装置２０の分散特性算出手
段３１にて微動観測により求めた観測分散特性と合致す
るように参考速度モデルを修正する逆解析を行い、当該
地盤の速度構造を算出して推定し、出力装置２２に結果
を出力するものである。
【００１５】図２において、地盤データベース２３から
選定した３地点Ａ、Ｂ、Ｃは、調査目的地Ｐの周辺に位
置しており、３地点は調査目的地から２００ｍから３０
０ｍ程度離れて位置しており、これら３地点は例えばボ
ーリング調査により柱状図が得られている。そして、こ
れらの柱状図データは、地盤データベース２３から解析
装置２０に入力される。また、調査目的地Ｐに前記した
センサアレイ１５が設置される。
選定した３地点Ａ、Ｂ、Ｃは、調査目的地Ｐの周辺に位
置しており、３地点は調査目的地から２００ｍから３０
０ｍ程度離れて位置しており、これら３地点は例えばボ
ーリング調査により柱状図が得られている。そして、こ
れらの柱状図データは、地盤データベース２３から解析
装置２０に入力される。また、調査目的地Ｐに前記した
センサアレイ１５が設置される。
【００１６】前記の如く構成された本実施形態の地盤速
度構造の推定システム１０を使用して、地盤速度構造を
推定する方法について以下に説明する。振動センサ１１
〜１４により検出された複数箇所の常時微動データは解
析装置２０の分散特性算出手段３１により解析され、調
査地盤のレイリー波の分散特性が得られる。一つの振動
センサから得られる分散特性は周波数範囲が狭いので、
複数の振動センサ１１〜１４から得られた複数の分散特
性を合成し、幅広い周波数範囲の分散特性を作成する。
度構造の推定システム１０を使用して、地盤速度構造を
推定する方法について以下に説明する。振動センサ１１
〜１４により検出された複数箇所の常時微動データは解
析装置２０の分散特性算出手段３１により解析され、調
査地盤のレイリー波の分散特性が得られる。一つの振動
センサから得られる分散特性は周波数範囲が狭いので、
複数の振動センサ１１〜１４から得られた複数の分散特
性を合成し、幅広い周波数範囲の分散特性を作成する。
【００１７】地盤データベース２３を用いて、図３に示
す調査目的地Pの周辺の３地点Ａ、Ｂ、Ｃを入力装置２
１により入力し、３地点の比較的入手しやすいボーリン
グ柱状図を基に、解析装置２０によりＳ波速度を算出す
る。ここでは、Ｓ波速度の算出式である式（１）に従っ
て、N値、地表面からの深さH、時代係数K1、土質係数K2
から換算する。
す調査目的地Pの周辺の３地点Ａ、Ｂ、Ｃを入力装置２
１により入力し、３地点の比較的入手しやすいボーリン
グ柱状図を基に、解析装置２０によりＳ波速度を算出す
る。ここでは、Ｓ波速度の算出式である式（１）に従っ
て、N値、地表面からの深さH、時代係数K1、土質係数K2
から換算する。
【００１８】
【数１】 Ｖs＝６８．９×Ｎ0.173×Ｈ0.195×Ｋ1×Ｋ2 （１） ここで、Vｓ：Ｓ波速度、N：N値、H：地表面からの深さ（ｍ） K1：時代係数（沖積：1.000、洪積：1.306） K2：土質係数（粘土：1.000、砂：1.085、礫：1.189）である。
【００１９】このようにして、Ｓ波速度を算出し、図４
に示す地盤速度構造図を求め出力装置２２により出力す
る。そして、この地盤速度構造から、卓越周期算出手段
３２により、Ｓ波の重複反射理論に従って地盤の増幅率
（基盤面を１としたときの、地表面での波の振動の大き
さ）を各周波数で求め、図５に示す３地点の理論上の卓
越周期Ｆａ、Ｆｂ、Ｆｃを算出する。
に示す地盤速度構造図を求め出力装置２２により出力す
る。そして、この地盤速度構造から、卓越周期算出手段
３２により、Ｓ波の重複反射理論に従って地盤の増幅率
（基盤面を１としたときの、地表面での波の振動の大き
さ）を各周波数で求め、図５に示す３地点の理論上の卓
越周期Ｆａ、Ｆｂ、Ｆｃを算出する。
【００２０】この重複反射理論は、水平方向に半無限大
に広がる複数の水平層からなり基盤層として半無限地盤
を有する地盤において、水平変位のみを考え、表層とあ
る１層の振幅の関係を導き出すものであり、全ての層の
振幅は計算でき、歪、加速度も変位から求めることがで
きるものである。なお、ここでは理論上の卓越周期を求
めているが、３地点Ａ、Ｂ、Ｃにおいて微動計を用いて
振動を検出し、フーリエ変換する手順を踏むことによっ
て実測データとしてのパワースペクトルから卓越周期を
求めてもよい。
に広がる複数の水平層からなり基盤層として半無限地盤
を有する地盤において、水平変位のみを考え、表層とあ
る１層の振幅の関係を導き出すものであり、全ての層の
振幅は計算でき、歪、加速度も変位から求めることがで
きるものである。なお、ここでは理論上の卓越周期を求
めているが、３地点Ａ、Ｂ、Ｃにおいて微動計を用いて
振動を検出し、フーリエ変換する手順を踏むことによっ
て実測データとしてのパワースペクトルから卓越周期を
求めてもよい。
【００２１】調査目的地Ｐにおいて、地表面で水平方向
の常時微動を計測手段３３により計測し、そのパワース
ペクトルを求めることにより観測卓越周期Ｆｐを計測す
る。観測卓越周期は、基本的には直交する水平方向（例
えば、N−Ｓ、Ｅ−Ｗ方向）で地表面にて振動を計測
し、そのスペクトルの相乗平均で得られる。なお、前記
した観測卓越周期の計測方法の他に、水平方向と鉛直方
向の比のピークを卓越する方法で計測してもよく、水平
方向だけを計測してそのパワースペクトルから卓越周期
を計測してもよい。このようにして得られた調査目的地
Ｐにおける観測卓越周期Fｐを図６に示す。
の常時微動を計測手段３３により計測し、そのパワース
ペクトルを求めることにより観測卓越周期Ｆｐを計測す
る。観測卓越周期は、基本的には直交する水平方向（例
えば、N−Ｓ、Ｅ−Ｗ方向）で地表面にて振動を計測
し、そのスペクトルの相乗平均で得られる。なお、前記
した観測卓越周期の計測方法の他に、水平方向と鉛直方
向の比のピークを卓越する方法で計測してもよく、水平
方向だけを計測してそのパワースペクトルから卓越周期
を計測してもよい。このようにして得られた調査目的地
Ｐにおける観測卓越周期Fｐを図６に示す。
【００２２】つぎに、調査目的地Ｐの観測卓越周期Ｆｐ
と、地盤データベース上の３地点の理論上の卓越周期Ｆ
ａ、Ｆｂ、Ｆｃとを比較し、最も近似しているものを選
択し、選択された理論上の卓越周期Ｆｃを有する地盤デ
ータを参考速度モデルとして採用する。この参考速度モ
デルを入力装置２１により入力し、これと前記したレイ
リー波の分散特性に基づいて、インバージョン法を用い
て速度構造算出手段３４により、調査目的地Ｐにおける
地盤速度構造を算出して推定し、出力装置２２より出力
する。
と、地盤データベース上の３地点の理論上の卓越周期Ｆ
ａ、Ｆｂ、Ｆｃとを比較し、最も近似しているものを選
択し、選択された理論上の卓越周期Ｆｃを有する地盤デ
ータを参考速度モデルとして採用する。この参考速度モ
デルを入力装置２１により入力し、これと前記したレイ
リー波の分散特性に基づいて、インバージョン法を用い
て速度構造算出手段３４により、調査目的地Ｐにおける
地盤速度構造を算出して推定し、出力装置２２より出力
する。
【００２３】この場合、微動を計測した調査目的地Ｐと
地盤データベース２３の地盤における速度構造は大きく
変わらないので、解析装置２０は採用された参考速度モ
デルに基づいて、インバージョン法の探索範囲を限定し
て調査目的地Ｐの速度構造を算出して推定する。例えば
土質データの各層におけるＳ波速度と層厚を、参考速度
モデルの０．５から２倍程度に限定することにより、短
時間で正確に地盤の速度構造を推定することができる。
地盤データベース２３の地盤における速度構造は大きく
変わらないので、解析装置２０は採用された参考速度モ
デルに基づいて、インバージョン法の探索範囲を限定し
て調査目的地Ｐの速度構造を算出して推定する。例えば
土質データの各層におけるＳ波速度と層厚を、参考速度
モデルの０．５から２倍程度に限定することにより、短
時間で正確に地盤の速度構造を推定することができる。
【００２４】また、堆積等の状況によっては、調査目的
地Ｐと地盤データベースの３地点Ａ、Ｂ、Ｃの土質は大
きく変わらず、このためＳ波速度は近い値と予想される
が、層厚は大きく変わる可能性もある。この場合は、例
えば各層のＳ波速度は参考速度モデルの０．５から２倍
程度に限定し、層厚は参考速度モデルの０．３から３倍
程度に限定し、層厚の限定範囲をＳ波速度より広くする
ことにより、短時間で正確に地盤の速度構造を推定する
ことができる。このような判定は、オペレータが調査目
的地Ｐと地盤データベースの３地点間の距離、地質状況
に応じて総合的に行なう。
地Ｐと地盤データベースの３地点Ａ、Ｂ、Ｃの土質は大
きく変わらず、このためＳ波速度は近い値と予想される
が、層厚は大きく変わる可能性もある。この場合は、例
えば各層のＳ波速度は参考速度モデルの０．５から２倍
程度に限定し、層厚は参考速度モデルの０．３から３倍
程度に限定し、層厚の限定範囲をＳ波速度より広くする
ことにより、短時間で正確に地盤の速度構造を推定する
ことができる。このような判定は、オペレータが調査目
的地Ｐと地盤データベースの３地点間の距離、地質状況
に応じて総合的に行なう。
【００２５】よく知られているように地盤構造は地域に
よって多種多様であり、同一地域内で僅かな距離しか離
れていない地点でも、全く異なる構造であることも珍し
くない。しかし、ある地盤速度構造が決まれば、伝達し
やすい波の周波数が決定されるため、それぞれの地盤特
有の周波数特性が定義され、そのピークである卓越周期
が存在する。
よって多種多様であり、同一地域内で僅かな距離しか離
れていない地点でも、全く異なる構造であることも珍し
くない。しかし、ある地盤速度構造が決まれば、伝達し
やすい波の周波数が決定されるため、それぞれの地盤特
有の周波数特性が定義され、そのピークである卓越周期
が存在する。
【００２６】本発明に係る地盤速度構造の推定方法によ
れば、調査対象地盤におけるレイリー波の分散特性を求
めると共に、調査目的地の近傍の地盤に関する地盤デー
タベースから、卓越周期を基に最適な地盤データを選択
してインバージョン法の初期速度構造モデル（参考速度
モデル）を求め、これらの調査対象地点におけるレイリ
ー波の分散特性と参考速度モデルに基づいてインバージ
ョン法により調査目的地の地盤速度構造を算出して推定
するので、短時間で正確に地盤の速度構造を推定するこ
とができる。特に、建築物の耐震設計を行なう場合に、
考慮すべきS波速度で４００ｍ／ｓを越える工学的基盤
までを調査対象深度とするときには、その卓越周期が
０．１秒から１秒の範囲となるため、その範囲で比較す
ることは効率的であり、かつ妥当である。なお、前記し
た実施の形態では、調査目的地の周辺のサンプル地を３
地点とした例を示したが、これに限られるものでなく、
さらに多くのサンプル地から理論卓越周期を算出して比
較するようにしてもよい。
れば、調査対象地盤におけるレイリー波の分散特性を求
めると共に、調査目的地の近傍の地盤に関する地盤デー
タベースから、卓越周期を基に最適な地盤データを選択
してインバージョン法の初期速度構造モデル（参考速度
モデル）を求め、これらの調査対象地点におけるレイリ
ー波の分散特性と参考速度モデルに基づいてインバージ
ョン法により調査目的地の地盤速度構造を算出して推定
するので、短時間で正確に地盤の速度構造を推定するこ
とができる。特に、建築物の耐震設計を行なう場合に、
考慮すべきS波速度で４００ｍ／ｓを越える工学的基盤
までを調査対象深度とするときには、その卓越周期が
０．１秒から１秒の範囲となるため、その範囲で比較す
ることは効率的であり、かつ妥当である。なお、前記し
た実施の形態では、調査目的地の周辺のサンプル地を３
地点とした例を示したが、これに限られるものでなく、
さらに多くのサンプル地から理論卓越周期を算出して比
較するようにしてもよい。
【００２７】
【発明の効果】以上の説明から理解できるように、本発
明の地盤速度構造の推定方法によれば、調査目的地から
計測した観測卓越周期と、周辺の地盤データベースから
算出した複数の理論卓越周期とを比較し、最も近似して
いる理論卓越周期を有する地盤データを選択し、選択さ
れた地盤データを参考速度モデルとして逆解析を行なっ
て地盤速度構造を算出して推定するため、地下の地盤速
度構造を短時間で、正確に推定することができる。
明の地盤速度構造の推定方法によれば、調査目的地から
計測した観測卓越周期と、周辺の地盤データベースから
算出した複数の理論卓越周期とを比較し、最も近似して
いる理論卓越周期を有する地盤データを選択し、選択さ
れた地盤データを参考速度モデルとして逆解析を行なっ
て地盤速度構造を算出して推定するため、地下の地盤速
度構造を短時間で、正確に推定することができる。
【００２８】また、前記した地盤速度構造の測定方法に
おいて、卓越周期の範囲は、０．１秒から１秒の範囲に
することにより、地下の地盤速度構造をさらに短時間
に、効率よく推定することができる。本発明の地盤速度
構造の推定システムによれば、前記した推定方法によ
り、調査目的地の地盤速度構造を短時間に、正確に推定
することができる。
おいて、卓越周期の範囲は、０．１秒から１秒の範囲に
することにより、地下の地盤速度構造をさらに短時間
に、効率よく推定することができる。本発明の地盤速度
構造の推定システムによれば、前記した推定方法によ
り、調査目的地の地盤速度構造を短時間に、正確に推定
することができる。
【図１】本発明に係る地盤速度構造推定システムの一実
施形態のブロック図。
施形態のブロック図。
【図２】地盤データベースより選んだ調査目的地周辺の
３地点の位置関係を示す平面図。
３地点の位置関係を示す平面図。
【図３】図２の３地点の柱状図。
【図４】図２の３地点の地盤速度構造図。
【図５】図２の３地点の理論卓越周期を示す特性図。
【図６】調査目的地の観測卓越周期を示す特性図。
【図７】選択された理論卓越周期の横軸を広い範囲とし
た特性図。
た特性図。
１０ 推定システム １１、１２、１３、１４ 振動センサ １５ センサアレイ ２０ 解析装置 ２１ 入力装置 ２２ 出力装置 ２３ 地盤データベース ３１ 分散特性算出手段 ３２ 卓越周期算出手段 ３３ 計測手段 ３４ 速度構造算出手段
Claims (3)
- 【請求項１】 調査目的地にて複数の振動センサを有す
る微動観測装置により地盤の微動探査を行い、その分散
特性から参考速度モデルを初期値にして逆解析により地
盤速度構造を推定する方法であり、 前記調査目的地周辺の複数の地盤データから複数の理論
卓越周期を算出し、前記調査目的地にて計測した観測卓
越周期と前記複数の理論卓越周期とを比較して最も近似
している理論卓越周期を有する地盤データを選択し、こ
の選択された地盤データを参考速度モデルとして逆解析
を行い、前記調査目的地の地盤速度構造を推定すること
を特徴とする地盤速度構造の推定方法。 - 【請求項２】 前記観測卓越周期及び理論卓越周期の範
囲は、０．１秒から１秒の範囲であることを特徴とする
請求項１記載の地盤速度構造の推定方法。 - 【請求項３】 調査目的地にて微動探査を行なう複数の
振動センサと、該複数の振動センサが接続された解析装
置と、該解析装置に入力する入力装置と、前記解析装置
に接続された出力装置と、前記解析装置に接続された地
盤データベースとを備え、 前記解析装置は、前記振動センサからの微動データから
観測分散特性を算出する分散特性算出手段と、前記入力
装置により入力された複数の地点における前記地盤デー
タベースの地盤データからそれぞれの地点における複数
の理論卓越周期を算出する卓越周期算出手段と、前記調
査目的地の観測卓越周期を計測する計測手段と、逆解析
により地盤の速度構造を算出して前記出力装置へ出力す
る速度構造算出手段とを有し、 前記複数の理論卓越周期と前記調査目的地の観測卓越周
期とを比較し、最も近似している理論卓越周期を有する
地盤データを選択して前記速度構造算出手段へ参考速度
モデルとして入力し、前記速度構造算出手段において前
記参考速度モデルを初期値にして逆解析を行い前記調査
目的地の地盤速度構造を算出することを特徴とする地盤
速度構造の推定システム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000003589A JP4203199B2 (ja) | 2000-01-12 | 2000-01-12 | 地盤速度構造の推定方法及び推定システム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000003589A JP4203199B2 (ja) | 2000-01-12 | 2000-01-12 | 地盤速度構造の推定方法及び推定システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001193046A true JP2001193046A (ja) | 2001-07-17 |
| JP4203199B2 JP4203199B2 (ja) | 2008-12-24 |
Family
ID=18532507
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000003589A Expired - Lifetime JP4203199B2 (ja) | 2000-01-12 | 2000-01-12 | 地盤速度構造の推定方法及び推定システム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4203199B2 (ja) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103261916A (zh) * | 2010-09-30 | 2013-08-21 | 韩国电力公社 | 地震的累积绝对速度推定装置 |
| JP2017193841A (ja) * | 2016-04-19 | 2017-10-26 | 株式会社大林組 | 地盤推定方法 |
| JP2018178618A (ja) * | 2017-04-19 | 2018-11-15 | 公益財団法人鉄道総合技術研究所 | 支持地盤の深度測定方法 |
| JP2019007808A (ja) * | 2017-06-23 | 2019-01-17 | 戸田建設株式会社 | 地盤の表面波の位相速度推定方法 |
| JP2020133272A (ja) * | 2019-02-20 | 2020-08-31 | 株式会社大林組 | 地盤推定方法 |
| JP2021032055A (ja) * | 2019-08-29 | 2021-03-01 | 公益財団法人鉄道総合技術研究所 | 支持地盤の深度推定方法 |
| CN112882098A (zh) * | 2021-01-20 | 2021-06-01 | 长江勘测规划设计研究有限责任公司 | 一种天然源面波速度计算及地形校正方法 |
| JP2021169945A (ja) * | 2020-04-15 | 2021-10-28 | 株式会社ジオシス | 海底地盤の探査装置 |
| CN114185093A (zh) * | 2021-12-07 | 2022-03-15 | 中国石油大学(北京) | 一种基于瑞雷面波反演的近地表速度模型建立方法及装置 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2021031863A (ja) * | 2019-08-19 | 2021-03-01 | 積水化学工業株式会社 | 地盤調査方法 |
-
2000
- 2000-01-12 JP JP2000003589A patent/JP4203199B2/ja not_active Expired - Lifetime
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103261916B (zh) * | 2010-09-30 | 2016-02-10 | 韩国电力公社 | 地震的累积绝对速度推定装置 |
| CN103261916A (zh) * | 2010-09-30 | 2013-08-21 | 韩国电力公社 | 地震的累积绝对速度推定装置 |
| JP2017193841A (ja) * | 2016-04-19 | 2017-10-26 | 株式会社大林組 | 地盤推定方法 |
| JP2018178618A (ja) * | 2017-04-19 | 2018-11-15 | 公益財団法人鉄道総合技術研究所 | 支持地盤の深度測定方法 |
| JP2019007808A (ja) * | 2017-06-23 | 2019-01-17 | 戸田建設株式会社 | 地盤の表面波の位相速度推定方法 |
| JP7183855B2 (ja) | 2019-02-20 | 2022-12-06 | 株式会社大林組 | 地盤推定方法 |
| JP2020133272A (ja) * | 2019-02-20 | 2020-08-31 | 株式会社大林組 | 地盤推定方法 |
| JP2021032055A (ja) * | 2019-08-29 | 2021-03-01 | 公益財団法人鉄道総合技術研究所 | 支持地盤の深度推定方法 |
| JP7181532B2 (ja) | 2019-08-29 | 2022-12-01 | 公益財団法人鉄道総合技術研究所 | 支持地盤の深度推定方法 |
| JP2021169945A (ja) * | 2020-04-15 | 2021-10-28 | 株式会社ジオシス | 海底地盤の探査装置 |
| JP7103606B2 (ja) | 2020-04-15 | 2022-07-20 | 株式会社ジオシス | 海底地盤の探査装置 |
| CN112882098B (zh) * | 2021-01-20 | 2022-03-11 | 长江勘测规划设计研究有限责任公司 | 一种天然源面波速度计算及地形校正方法 |
| CN112882098A (zh) * | 2021-01-20 | 2021-06-01 | 长江勘测规划设计研究有限责任公司 | 一种天然源面波速度计算及地形校正方法 |
| CN114185093A (zh) * | 2021-12-07 | 2022-03-15 | 中国石油大学(北京) | 一种基于瑞雷面波反演的近地表速度模型建立方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4203199B2 (ja) | 2008-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7764572B2 (en) | Methods and systems for acoustic waveform processing | |
| Asten et al. | Optimised array design for microtremor array studies applied to site classification; comparison of results with SCPT logs | |
| US7590491B2 (en) | Signal integration measure for seismic data | |
| EP2166378A2 (en) | Interferometric Directional Balancing | |
| JP2001193046A (ja) | 地盤速度構造の推定方法及び推定システム | |
| Lin et al. | A study on issues relating to testing of soils and pavements by surface wave methods | |
| US20140142854A1 (en) | Method for locating a microseismic event | |
| Drijkoningen | The usefulness of geophone ground-coupling experiments to seismic data | |
| Yoon et al. | Combined active-passive surface wave measurements for near-surface site characterization | |
| JP4263840B2 (ja) | 地盤速度構造の推定方法 | |
| RU2357078C2 (ru) | Способ определения траектории скважины | |
| García-Jerez et al. | Determination of elastic properties of shallow sedimentary deposits applying a spatial autocorrelation method | |
| JP2000321360A (ja) | 地下速度構造推定方法及び地下速度構造推定システム | |
| Rosenblad et al. | Comparative study of refraction microtremor (ReMi) and active source methods for developing low-frequency surface wave dispersion curves | |
| Foti et al. | Application of FK analysis of surface waves for geotechnical characterization | |
| Rosenblad et al. | Performance of active and passive methods for measuring low-frequency surface wave dispersion curves | |
| JP4195171B2 (ja) | 地盤構造推定方法 | |
| JPH06294793A (ja) | 地層の物理特性の音響波を用いた非破壊測定方法 | |
| TOKIMATSlf et al. | Three-dimensional soil profiling in Kobe area using microtremors | |
| JP2001107349A (ja) | 地盤速度構造の推定方法および推定システム | |
| Gohil et al. | Pavement Testing Using Nondestructive MASW Approach | |
| JPH11248850A (ja) | 周波数−波数スペクトル解析方法 | |
| Kavand et al. | Determination of shear wave velocity profile of sedimentary deposits in Bam city (southeast of Iran) using microtremor measurements | |
| Moffat et al. | Comparación del promedio de la velocidad de onda de corte en los primeros 30 m usando ensayos downhole, MASW y bender elements | |
| Bawadi et al. | Relationship between Shear Wave Velocity and SPT-N Value for Residual Soils |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061018 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080325 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080409 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080917 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20081010 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 4203199 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111017 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121017 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121017 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131017 Year of fee payment: 5 |
|
| EXPY | Cancellation because of completion of term |