JP4195171B2 - 地盤構造推定方法 - Google Patents
地盤構造推定方法 Download PDFInfo
- Publication number
- JP4195171B2 JP4195171B2 JP2000128099A JP2000128099A JP4195171B2 JP 4195171 B2 JP4195171 B2 JP 4195171B2 JP 2000128099 A JP2000128099 A JP 2000128099A JP 2000128099 A JP2000128099 A JP 2000128099A JP 4195171 B2 JP4195171 B2 JP 4195171B2
- Authority
- JP
- Japan
- Prior art keywords
- ground
- type
- target
- sensor
- vibration detection
- 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.)
- Expired - Lifetime
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Description
【発明の属する技術分野】
本発明は、地盤の地下構造を推定する方法に係り、特に、予め推定システムにより得られた分散特性から振動検知センサの間隔を最適に設定し、効率よく地盤構造を推定することができる方法に関する。
【0002】
【従来の技術】
複数の振動検知センサを地表面に配置して、常に存在する微振動である常時微動を収録する微動アレイ観測では、複数の振動センサで地表面に伝わる表面波を計測している。表面波の波長と地盤を進行する深さとは密接な関係があり、波長が短いほど地表面付近を進行し、逆に波長が長いほど地表面から深層部まで進行することがわかっている。そのため、より深くまでの地盤構造を推定するためには、より長い波長を持つ波を検知する必要が生じる。
【0003】
このとき、地表面に配置した複数の振動センサでは、収録された波の時間的変化あるいはセンサ間の周波数の関係を算出しているため、センサ間隔が広いほど、より長い波長を持つ波を検知することができる。微動観測では、振動センサで収録された微動データを周波数解析により、周波数と位相速度の関係(以下、分散特性という)を導き、地盤構造に固有の特性を示す、この分散特性を用いることで地盤の層構造の推定が可能となる。
【0004】
従来、この種の微動観測による地盤構造の推定方法としては種々の文献があり、例えば振動センサの配置間隔と分散特性の関係を紹介したものがある。この場合は、地盤構造の推定を行う振動センサの配置間隔は、分散特性の関係のみにより決定され、決定された配置間隔で振動センサを配置して微動データを計測し、前記した方法により地盤の地下構造が推定される。
【0005】
【発明が解決しようとする課題】
ところで、前記した推定方法は、分散特性の関係によりセンサ間隔が決定されるが、地盤は多種多様のものがあるため、センサ間隔が対象地盤に最適な間隔とは異なることがある。このため、前記した推定方法では、対象地盤の正確な地下構造の推定が短時間で効率よく行えないという問題点があった。
【0006】
本発明は、このような問題に鑑みてなされたものであって、その目的とするところは、予め求められた地盤種別ごとのセンサ間隔と探査深度との関係式から、調査しようとする対象地盤と目標の探査深度に最適のセンサ間隔を決定し、このセンサ間隔で振動検知センサを配置して、短時間で効率的に、しかも地盤の地下構造を精度よく推定できる地盤構造推定方法を提供することにある。
【0007】
【課題を解決するための手段】
前記目的を達成すべく、本発明に係る地盤構造推定方法は、複数の振動検知センサと、前記複数の振動検知センサの出力を解析する解析装置とを備える推定システムにより、予め求めた地盤種別ごとのセンサ間隔と探査深度との関係式を使用して地盤の地下構造を推定する方法であって、該方法は、対象地盤の地盤種別を決める第1のステップと、対象地盤の目標深度を決める第2のステップと、決められた深度と前記関係式から決められた対象地盤の地盤種別に基づいてセンサ間隔を決める第3のステップと、前記第3のステップで決められたセンサ間隔により複数の振動検知センサを地盤上に設置して微動観測を行う第4のステップと、前記微動観測により得られたデータを前記解析装置に入力して分散特性を算出する第5のステップと、前記分散特性を基に逆解析により地盤の地下構造を推定する第6のステップとを備えることを特徴とする。
【0008】
また、本発明に係る地盤構造推定方法の好ましい具体的な態様としては、前記関係式は、目標の探査深度が同じ場合にはセンサ間隔は第1種地盤<第2種地盤<第3種地盤の関係となり、センサ間隔が同じ場合には目標の探査深度は第1種地盤>第2種地盤>第3種地盤の関係となることを特徴としている。
【0009】
このように構成された本発明の地盤構造推定方法によれば、予め求めた関係式から、地盤種別と目標の探査深度に最適なセンサ間隔を決定でき、そのセンサ間隔で振動検知センサを配置するため、短時間で効率よく、正確な地盤の地下構造を推定することができる。また、実際の微動アレイ観測の前に、振動検知センサの配置計画を容易とすることができる。前記関係式は、目標の探査深度が同じ場合、センサ間隔は第1種地盤<第2種地盤<第3種地盤の関係に設置すればよく、センサ間隔が同じ場合、目標の探査深度は第1種地盤>第2種地盤>第3種地盤の関係となり、探査深度の判断が容易となる。
【0010】
【発明の実施の形態】
以下、本発明に係る地盤構造推定方法を実施する推定システムの一実施形態を図面に基づき詳細に説明する。図1は、微動探査により地盤の地下構造を推定する推定システムの一実施形態のブロック図、図2は振動センサの配置図である。図1、2において、調査目的地の地表面には複数の振動検知センサ1〜4が設置されている。振動検知センサ1〜4は、図2に示すように正三角形上に振動検知センサ1、2、3が設置されると共に、振動検知センサ4が3つの振動検知センサの中心位置に設置され、センサアレイ5を構成している。
【0011】
センサアレイ5のサイズは、対象とする地盤あるいは、求めたい地盤構造の深さで決定され、目標深度が深いほどセンサアレイのサイズは大きくなるため、目標深度までの地下の構造を推定するためには、その深度に合わせたアレイサイズが必要となる。センサアレイ5のサイズは、センサアレイを構成する振動検知センサ1〜4の間隔a、bで定義される。振動検知センサ1〜4は常時微動を検出して計測するものであり、常時微動とは、海洋の振動、風や地殻の変動等の自然現象に起因する地球固有の振動の他、交通や工場等の人間活動に起因する人工の振動であり、周波数帯域は1〜30Hz程度である。
【0012】
振動検知センサ1〜4により計測された複数箇所の常時微動データはデータ収録装置6に入力され、その後、解析装置10に入力される。解析装置10にはデータ等を入力する入力装置11が接続され、解析装置10からの解析結果を出力する出力装置12が接続されている。出力装置12はプリンターや表示装置等が用いられる。推定システムは、複数の振動検知センサ1〜4から構成されるセンサアレイ5、そのデータ出力を解析する解析装置10、入力装置11及び出力装置12から構成される。
【0013】
解析装置10は、振動検知センサ1〜4からの微動観測による微動データよりレイリー波の分散特性を算出し、その分散特性から地盤の地下構造を推定するものである。地盤の地下構造の推定は、遺伝アルゴリズム法あるいは最小二乗法といった手法により逆解析を行う。
【0014】
ここで、建築基準法施行令で定められている地盤種別の分類について表1に示す(建築基準法施行令、告示、昭55建告第1798号第2)。地盤の工学的性質は、地層の形成過程によって最も大きく影響を受け、これらの性質ごとに地盤種別と呼ばれる分類がなされている。地盤種別は硬い地盤、軟らかい地盤、両者の中間のものと3種類に分類される。第1種地盤は硬く、主に山地あるいは岩盤がむき出しになったところ、第2種地盤は普通の平地あるいは台地、住宅が立ち並んでいるところ、第3種地盤は水田や河川のそば等、軟らかいところであり、このように地形的あるいは土地条件によって、おおよその地盤種別の区分は可能である。しかし、複雑な地形では人的な判断に大きく影響されるため、定量的に区分でき、誰にでも判断できる明確な方法が必要となる。後述する表2に記載されている地盤周期は、対象となる地盤において地震動の応答スペクトルの値が、ある周期近傍に渡って著しく卓越すると考えられる周期を示している。地盤周期(=卓越周期Tg)は、1点の常時微動の周期頻度曲線等から求められるため、容易に計測、解析することができる。微動データは様々な振動波が混合しているので、微動のような時間変化のデータはフーリエ変換により周波数(周期の逆数)と強度(頻度)の関係に変換できる。常時微動の周期頻度曲線は、こういった計算により導くことができる。
【0015】
【表1】
【0016】
また、建築基準法施行令による地盤種別と卓越周期の分類を表2に示す。(建築基準法施行令、通達、昭56住指発第96号)
【0017】
【表2】
【0018】
前記の如く構成された本実施形態の地盤構造の推定システムの動作について以下に説明する。先ず、予めセンサ間隔と探査深度との関係を求める。例えば、前記した推定システムにより、振動検知センサ1〜4は調査地盤の地表面に設置され、常時微動を計測する。振動検知センサ1〜4の測定データはデータ収録装置6に収録され、解析装置10に入力される。解析装置10は前記測定データから図3に示される位相速度の分散特性を算出する。そして、この分散特性より振動検知センサ1〜4の間隔ごとに最大の波長を図4のようにプロットして、センサ間隔(L)と最大波長(λ)との関係式を得る。前記の最大波長(λ)は、得られる分散特性の最大の位相速度÷その時の周波数で得ることができる。この図4に示される最大波長の近似曲線を引き、センサ間隔(L)と得られる最大波長(λ)の関係式を導く。この関係式を式1に示す。
λ=28.675×L0.709 (1)
【0019】
この後、PS検層等から得られた地盤層構造から理論的に分散特性を導き、工学的基盤の探査深度(D)と位相速度400m/s時の波長(λ)との関係を、図5に示すように地盤種別ごとに導き、式2、式3、式4に示される分散特性の波長(λ)と探査深度(D)との関係式を得る。
第1種地盤の場合 D=0.0532×λ1.369 (2)
第2種地盤の場合 D=0.1573×λ1.049 (3)
第3種地盤の場合 D=0.0036×λ1.549 (4)
【0020】
そして、式1と式2、式3、式4との関係から、センサ間隔(L)と探査深度(D)との間に、第1種地盤の場合、式5に示す関係を、第2種地盤の場合、式6に示す関係を、第3種地盤の場合、式7に示す関係を求め、これらの関係から図6に示す曲線a、b、cを求めた。
D=5.26×L0.971 (5)
D=5.32×L0.744 (6)
D=0.652×L1.10 (7)
【0021】
このようにして得られた図6に示される関係を用いて、調査地点における地盤の地下構造を推定する。なお、このようにして求められた関係は、近隣の地盤を調査するときに、何度でも活用することができる。図7に示す地盤構造を推定するフローチャートに基づいて説明する。図7において、先ず、第1ステップS1として、調査地点の対象地盤が、表1に示される地盤種別のどれに当たるかを調べる。全ての地盤は前記した3種類の地盤種別に分類することができ、ほとんどの場合、地形を見るだけで地盤種別を決めることができる。目視で判断しにくい場合は、卓越周期Tgを測定して表2を用いて決めることができる。
【0022】
そして、第2ステップS2として、対象地盤にて推定したい目標深度を決める。次に、第3ステップS3として、探査深度(D)とセンサ間隔(L)を示す図6の関係に基づいて、センサアレイ5のサイズを決める。例えば、対象地盤が第3種地盤で、目標深度が60mの場合、図6の探査深度の60mの横線と第3種地盤の曲線cとの交点よりセンサ間隔L=約60mを得ることができる。このように、事前にセンサ間隔を決めることができるため、センサの配置計画を容易に行うことができる。
【0023】
このようにして得られたセンサ間隔L=60mを、センサアレイ5の1辺として正三角形の各頂点および中心上に振動検知センサ1〜4を配置する。次いで、第4ステップS4として、前記のように得られた間隔のセンサアレイ5により微動観測を行う。この微動観測による測定データは、解析装置10に入力され、第5ステップS5として、周波数−波数スペクトル法(F−K法)あるいは空間自己相関法(SPAC法)といった解析手法で、調査地盤固有の周波数と位相速度との関係、つまり分散特性が算出される。そして、この分散特性を基に、第6ステップS6として、遺伝アルゴリズム法あるいは最小二乗法といった手法により逆解析を行い、すなわち、ある深さにおけるS波速度と層厚の少なくとも一方をパラメータとして、逆解析を行う(インバージョン法)ことで調査地点の地盤の地下構造を推定する。
【0024】
前記のようにして推定された地盤の地下構造は、予め求められた探査深度(D)とセンサ間隔(L)との関係から振動検知センサ1〜4を配置しているため、対象地盤と目標深度に最適なセンサ間隔となっており、地盤の地下構造を短時間で効率的に推定することができ、しかも精度の高い推定が可能となる。また、図6に示される関係は、近隣の地盤の地下構造を推定するときに、何度でも使用することができ、一度関係を求めておけば、再度効率よく、精度よく、地盤の地下構造を推定することができる。
【0025】
【発明の効果】
以上の説明から理解できるように、本発明の地盤構造推定方法は、予め求めたセンサ間隔と探査深度との関係から、対象地盤と目標の探査深度に最適なセンサ間隔を決めることができ、このセンサ間隔で配置したセンサアレイにより短時間で効率よく、しかも精度よく地盤の地下構造を推定することができる。また、微動観測の前に、振動検知センサの配置計画を容易に立てることができる。
【図面の簡単な説明】
【図1】本発明に係る地盤構造推定方法を実施する推定システムの一実施形態のブロック図。
【図2】図1の振動検知センサの配置図。
【図3】地盤構造の推定システムにより求められた分散特性図。
【図4】図3の分散特性図より各センサ間隔ごとに得られた最大の波長をプロットし、近似曲線を引いた特性図。
【図5】地盤層構造から理論的に分散特性を導き、探査深度と波長との関係を地盤種別ごとに導いた特性図。
【図6】図4と図5の関係から、センサ間隔と探査深度との関係を求めた関係図。
【図7】本発明に係る地盤構造推定方法の一実施形態のフローチャート。
【符号の説明】
1〜4 振動検知センサ
5 センサアレイ
10 解析装置
Claims (2)
- 複数の振動検知センサと、前記複数の振動検知センサの出力を解析する解析装置とを備える推定システムにより、予め求めた地盤種別ごとのセンサ間隔と探査深度との関係式を使用して地盤の地下構造を推定する方法であって、
該方法は、対象地盤の地盤種別を決める第1のステップと、対象地盤の目標の探査深度を決める第2のステップと、決められた探査深度と前記関係式から決められた対象地盤の地盤種別に基づいてセンサ間隔を決める第3のステップと、前記第3のステップで決められたセンサ間隔により複数の振動検知センサを地盤上に設置して微動観測を行う第4のステップと、前記微動観測により得られたデータを前記解析装置に入力して分散特性を算出する第5のステップと、前記分散特性を基に逆解析により地盤の地下構造を推定する第6のステップとを備えることを特徴とする地盤構造推定方法。 - 前記関係式は、目標の探査深度が同じ場合にはセンサ間隔は第1種地盤<第2種地盤<第3種地盤の関係となり、センサ間隔が同じ場合には目標の探査深度は第1種地盤>第2種地盤>第3種地盤の関係となることを特徴とする請求項1記載の地盤構造推定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000128099A JP4195171B2 (ja) | 2000-04-27 | 2000-04-27 | 地盤構造推定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000128099A JP4195171B2 (ja) | 2000-04-27 | 2000-04-27 | 地盤構造推定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2001305235A JP2001305235A (ja) | 2001-10-31 |
JP4195171B2 true JP4195171B2 (ja) | 2008-12-10 |
Family
ID=18637595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000128099A Expired - Lifetime JP4195171B2 (ja) | 2000-04-27 | 2000-04-27 | 地盤構造推定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4195171B2 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6876296B2 (ja) * | 2017-04-19 | 2021-05-26 | 公益財団法人鉄道総合技術研究所 | 支持地盤の深度測定方法 |
JP6886873B2 (ja) * | 2017-06-23 | 2021-06-16 | 戸田建設株式会社 | 地盤の表面波の位相速度推定方法 |
JP7181532B2 (ja) * | 2019-08-29 | 2022-12-01 | 公益財団法人鉄道総合技術研究所 | 支持地盤の深度推定方法 |
JP7103606B2 (ja) * | 2020-04-15 | 2022-07-20 | 株式会社ジオシス | 海底地盤の探査装置 |
-
2000
- 2000-04-27 JP JP2000128099A patent/JP4195171B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2001305235A (ja) | 2001-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Foti et al. | Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project | |
Molnar et al. | Application of microtremor horizontal-to-vertical spectral ratio (MHVSR) analysis for site characterization: State of the art | |
RU2106000C1 (ru) | Способ для оценки ламинарной структуры и других характеристик земли | |
Whiteley et al. | Landslide monitoring using seismic refraction tomography–The importance of incorporating topographic variations | |
CN103728659B (zh) | 一种提高地下岩溶探测精度的方法 | |
Anbazhagan et al. | Spatial variability of the depth of weathered and engineering bedrock using multichannel analysis of surface wave method | |
Coccia et al. | Application of Refraction Microtremor (ReMi) technique for determination of 1-D shear wave velocity in a landslide area | |
Gabàs et al. | Combination of geophysical methods to support urban geological mapping | |
CN105022097A (zh) | 一种土质边坡滑动面综合预报方法 | |
Farrugia et al. | Noninvasive techniques for site characterization of Alberta seismic stations based on shear‐wave velocity | |
JP6531934B2 (ja) | ハイブリッド表面波探査方法及びハイブリッド表面波探査システム | |
CN114236612A (zh) | 一种针对陆内水域地震背景噪声的信号处理方法 | |
JP4195171B2 (ja) | 地盤構造推定方法 | |
RU2346300C1 (ru) | Способ обнаружения возможности наступления катастрофических явлений | |
Durdağ et al. | Combined application of electrical resistivity tomography and multi-channel analysis of surface waves methods in the tunnel detection: A case study from Kocaeli University Campus Site, Turkey | |
CN113900145B (zh) | 一种基于微动信号的地下横波速度分布的推测方法和装置 | |
JP2001193046A (ja) | 地盤速度構造の推定方法及び推定システム | |
CN111221043B (zh) | 一种被动源电场法的工作参数优化方法 | |
Kamalb et al. | Comparison of Shear Wave Velocity Derived from PS Logging and MASW–A Case Study of Mymensingh Pourashava, Bangladesh. | |
Feng et al. | Rockfall localization from seismic polarization considering multiple triaxial geophones and frequency bands | |
RU2645790C1 (ru) | Способ определения границ субвертикальных протяженных объектов в геологической среде | |
Samyn et al. | Assessment of vulnerability to erosion: Digital mapping of a loess cover thickness and stiffness using spectral analysis of seismic surface-waves | |
JP2004138447A (ja) | 岩盤の物性評価方法 | |
JP2000321360A (ja) | 地下速度構造推定方法及び地下速度構造推定システム | |
Anbazhagan et al. | A Study on the Effectiveness of Various Geophysical Methods in Detecting Naturally Formed Cavities in Lateritic Deposit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070111 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080613 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080709 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080811 |
|
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: 20080903 |
|
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: 20080925 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 4195171 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: 20111003 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121003 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121003 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131003 Year of fee payment: 5 |
|
EXPY | Cancellation because of completion of term |