JP3617036B2 - 地盤探査装置及びそれに使用される解析プログラム - Google Patents
地盤探査装置及びそれに使用される解析プログラム Download PDFInfo
- Publication number
- JP3617036B2 JP3617036B2 JP2001151961A JP2001151961A JP3617036B2 JP 3617036 B2 JP3617036 B2 JP 3617036B2 JP 2001151961 A JP2001151961 A JP 2001151961A JP 2001151961 A JP2001151961 A JP 2001151961A JP 3617036 B2 JP3617036 B2 JP 3617036B2
- Authority
- JP
- Japan
- Prior art keywords
- sine wave
- wave
- ground
- ground exploration
- waves
- 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
- 238000005259 measurement Methods 0.000 claims description 49
- 230000001133 acceleration Effects 0.000 claims description 32
- 238000001228 spectrum Methods 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 21
- 230000005284 excitation Effects 0.000 claims description 18
- 238000012546 transfer Methods 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 9
- 230000035945 sensitivity Effects 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 5
- 230000010355 oscillation Effects 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 description 15
- 238000010586 diagram Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、表面波、特にレイリー波の特性を利用した地盤探査装置に関する。
【０００２】
【従来の技術】
土木建設の分野における構造物基礎、トンネル、法面等の地盤調査には物理探査法が多く活用されている。なかでも弾性波探査法は地盤の工学的性質を定量的に把握する手法として用いられ、その代表的な例として広範囲に把握する屈折波法、局部的に把握する直接波法（ボーリング孔等を利用する）が知られている。これらの両手法で取り扱う波は人工振源で生じる弾性波動のうち実体波と呼ばれるＰ波及びＳ波であり、地質条件等を考慮して両波が使い分けされている。
【０００３】
これらの実体波とは別に、表面波と呼ばれる波の存在が知られており、表面波は分散現象を伴い、その伝播速度は波長に依存する性質を持っていることも知られている。以下に、表面波の性質について簡単に説明する。
【０００４】
起振機で地表面に上下振動を与えると、実体波（Ｐ波、Ｓ波）以外に表面波と呼ばれる波が地盤中に発生する。このうち実体波は、振源から半球状に伝播していく。これに対して、表面波のうちレイリー波は、周波数に規制された高さ一定の円筒が径を広げていくように伝播していく。つまり、起振機より発生する実体波には指向性があり、粗密波（Ｐ波）は起振機直下、せん断波（Ｓ波）は４５度方向に伝わるエネルギーが最大で、水平方向には伝わりにくい。このため、地表面で大振幅を与えるのは唯一表面波のみであると言える。
【０００５】
実体波（Ｐ波、Ｓ波）の幾何学的減衰係数は、下方に対してはｒ−１（但し、ｒは距離）であるが、地表付近ではｒ−２である。
【０００６】
これに対して表面波の地表面での幾何学的減衰係数は、ｒ−０．５である。
【０００７】
表面波にはレイリー波の他にラブ波、Ｐ波に類似するＰ−Ｐモード波、板波等があるが、本発明ではレイリー波を扱う。
【０００８】
レイリー波速度とＰ波及びＳ波速度の関係は、図８に示す通りであり、これらの３種類の速度比はポアソン比νで異なるが、Ｓ波速度とレイリー波速度は近似した値を示す。図８において、縦軸はＶ／ＶＳ ＝Ｖ・（ρ／Ｇ）１／２ で表される値である。但し、ＶＳ はＳ波速度、ρは密度、Ｇは剛性率であり、ＶにはＰ波あるいはレイリー波速度が当てはまる。
【０００９】
レイリー波速度はＳ波速度に近似しており、土質工学的数値を良く表していることからＳ波速度とＮ値（標準貫入試験値）の相関式やｑａ値（許容支持力）との相関式を用いてそれぞれの土質工学的数値を推定することができる。
【００１０】
【発明が解決しようとする課題】
上記のような知見に基づいて、本発明は、表面波、特にレイリー波の分散性を利用し、非破壊で正確に迅速でしかも経済的に探査を行うことのできる地盤探査装置を提供しようとするものである。
【００１１】
【課題を解決するための手段】
本発明によれば、地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査装置において、間隔Ｌをおいて地盤上に配置された２つの加速度検出器と、前記２つの加速度検出器からの検出信号を受けて加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）を出力する地震計部を含む計測器と、前記時系列信号Ａ（ｔ）、Ｂ（ｔ）を受けてあらかじめ定められた解析プログラムに基づく信号処理を行う信号処理部であって、フーリエ変換を行ってパワースペクトル及びクロススペクトルを算出すると共に、算出されたパワースペクトル及びクロススペクトルを用いて伝達関数Ｈ（ｆ）を算出し、算出された伝達関数Ｈ（ｆ）より前記２つの加速度検出器間の位相差Δθ（ｆ）とその時間差Δｔ（ｆ）を算出し、更に算出された時間差Δｔ（ｆ）と前記間隔Ｌより表面波の伝播平均速度Ｖｒ（ｆ）と深度Ｄ（ｆ）とを算出する信号処理部とを備え、前記計測器は更に、地盤を起振する起振機に対して与える励振信号を発生する発振部を含み、前記発振部は、励振信号として、ステップサイン波、又はマルチサイン波・スウェプトサイン波を切り替えて発生することが可能であり、前記ステップサイン波は周波数の異なる複数種類の正弦波で規定され、前記マルチサイン波は周波数の異なる複数種類の正弦波を合成した合成波で規定され、前記スウェプトサイン波は周波数の異なる複数種類の正弦波を振幅が一定の状態で合成した合成波で規定され、前記マルチサイン波・スウェプトサイン波を使用する場合には合成される周波数領域ｆ １ 〜ｆ ｎ を複数のバンドに分割し、分割された各バンドにおいて前記マルチサイン波・スウェプトサイン波による計測を行うことを特徴とする地盤探査装置が提供される。
【００１４】
本地盤探査装置においては更に、前記信号処理部は、前記伝播平均速度Ｖｒ （ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝播平均速度Ｖｒ 曲線を生成する。
【００１５】
本地盤探査装置においては更に、前記計測器が更に、前記地震計部からのアナログの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）をディジタルの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）に変換して出力するＡ／Ｄ変換部を含み、該Ａ／Ｄ変換部は、入力レベルに応じて入力感度を自動調整する機能を有する。
【００１６】
本地盤探査装置においては更に、前記信号処理部は、前記励振信号として前記ステップサイン波が使用される場合には、ＤＦＴ（離散フーリエ変換）による処理を行い、前記励振信号として前記マルチサイン波・スウェプトサイン波が使用される場合には、ＦＦＴ（高速フーリエ変換）による処理を行う。
【００１９】
【発明の実施の形態】
図１を参照して、本発明の好ましい実施の形態による地盤探査装置の構成について説明する。図１において、本装置は、２つの加速度検出器１１Ａ、１１Ｂと、これらの加速度検出器１１Ａ、１１Ｂからの加速度検出信号を受ける計測器１２とを含む。計測器１２は、地震計部１２−１、Ａ／Ｄ変換部１２−２、通信部１２−３、発振部１２−４とから成る。地震計部１２−１は、ローパスフィルタ回路を内蔵し、アナログの加速度検出信号から加速度時系列信号を生成するためのものである。Ａ／Ｄ変換部１２−２は、地震計部１２−１からのアナログの加速度時系列信号をディジタルの加速度時系列信号に変換するためのもので、入力感度の自動調整機能を有する。通信部１２−３は、ディジタルの加速度時系列信号を計測器１２に接続される機器、例えばモニタ付きのパーソナルコンピュータ（以下、ＰＣと略称する）１３に送信するためのものである。ＰＣ１３は、あらかじめインストールされた解析プログラムソフトに基づいて信号処理を行い、地盤解析を行うためのものである。
【００２０】
特に、本発明ではこの信号処理にＦＦＴ（高速フーリエ変換）、ＤＦＴ（離散フーリエ変換）を採用しているが、これについては後述する。
【００２１】
発振部１２−４は、起振機１５に与える励振信号を発生するためのもので、この励振信号は電力増幅器１４で増幅されて起振機１５に与えられる。発振部１２−４は、励振信号として、サイン波信号、マルチサイン信号・スウェプトサイン信号を発生することができる。ここで、サイン波信号というのは、正弦波のことであり、周波数の異なる複数種類のサイン波信号を生成できる。マルチサイン信号というのは、図２（ａ）に示すような異なる周波数ｆ１ 〜ｆｎ のサイン波信号を様々な振幅を持つように合成したもので、一例をあげると図２（ｂ）のような合成波である。一方、スウェプトサイン信号というのは、異なる周波数ｆ１ 〜ｆｎ のサイン波信号を振幅が一定の状態になるように合成したもので、一例をあげると図２（ｃ）のような合成波、いわば周波数変調波である。なお、図２（ｂ）、図２（ｃ）は便宜上の合成波形であり、図２（ａ）のサイン信号の実際の合成波形を示すものでは無い。
【００２２】
後述する説明で明らかとなるように、本形態では、励振信号としてサイン波信号を用いる場合には信号処理は上記のＤＦＴにより行われる（以下では、これをステップサイン計測と呼ぶ）。一方、マルチサイン信号・スウェプトサイン信号を用いる場合には信号処理は上記のＦＦＴにより行われる（以下では、これをマルチサイン計測と呼ぶ）。ＦＦＴの場合、周波数は２ｎ で決まる固定値となり、ステップ幅も固定である。しかも、設定する周波数を低い方から高い方に変化させてゆく必要がある。これに対し、ＤＦＴの場合には、周波数を任意に設定することができ、設定する周波数も低い方から高い方、高い方から低い方への両方の変化が可能である。
【００２３】
以下に動作を説明する。まず、探査地点に起振機１５、加速度検出器１１Ａ、１１Ｂを一直線上に設置する。加速度検出器１１Ａ、１１Ｂ間の距離をＬ（ｍ）とする。起振機１５を使って地表面を上下方向に起振することにより、起振機１５の周りに表面波を発生させる。地表面付近を伝播する表面波（レイリー波）の上下振動を、加速度検出器１１Ａ、１１Ｂで検出する。加速度検出器１１Ａ、１１Ｂからの加速度検出信号は、地震計部１２−１のローパスフィルタ回路を通すことで、アナログの時系列信号Ａ（ｔ）、Ｂ（ｔ）となり、Ａ／Ｄ変換部１２−２に入力される。Ａ／Ｄ変換部１２−２でＡ／Ｄ変換された時系列信号Ａ（ｔ）、Ｂ（ｔ）は通信部１２−３からＰＣ１３へ転送される。
【００２４】
ＰＣ１３では、あらかじめ定められた解析プログラムに基づいて時系列信号Ａ（ｔ）、Ｂ（ｔ）のパワースペクトルＧＡＡ（ｆ）、ＧＢＢ（ｆ）、クロススペクトルＧＢＡ（ｆ）、伝達関数Ｈ（ｆ）、コヒーレンス関数γ２ （ｆ）等を計算する。これらのパワースペクトルＧＡＡ（ｆ）、ＧＢＢ（ｆ）、クロススペクトルＧＢＡ（ｆ）、伝達関数Ｈ（ｆ）、コヒーレンス関数γ２ （ｆ）等は、内蔵のハードディスクに保存される。ＰＣ１３ではまた、伝達関数Ｈ（ｆ）より加速度検出器１１Ａ、１１Ｂ間の位相差Δθ（ｆ）を求め、続いてその時間差Δｔ（ｆ）を求める。ＰＣ１３では更に、時間差Δｔ（ｆ）と加速度検出器１１Ａ、１１Ｂ間の間隔Ｌより、表面波の伝播平均速度Ｖｒ （ｆ）と深度Ｄ（ｆ）とを求める。そして、求められた伝播平均速度Ｖｒ （ｆ）と深度Ｄ（ｆ）とよりモニター上にＤ−Ｖｒ 分散曲線を表示させる。
【００２５】
図３をも参照して、解析プログラムに基づく、上記の表面波の伝播平均速度Ｖｒ （ｆ）と深度Ｄ（ｆ）の算出方法について更に詳しく説明する。
【００２６】
ＰＣ１３では、Ａ／Ｄ変換された時系列信号Ａ（ｔ）、Ｂ（ｔ）に対し、下記の数１、数２で表されるフーリエ変換を行って変換信号ＳＡ （ｆ）、ＳＢ （ｆ）を得た後（ステップＳ１）、これらの変換信号ＳＡ （ｆ）、ＳＢ （ｆ）を用いて下記の数３、数４、及び数５によりパワースペクトルＧＡＡ（ｆ）、ＧＢＢ（ｆ）及びクロススペクトルＧＢＡ（ｆ）を計算する（ステップＳ２）。
【００２７】
【数１】
【００２８】
【数２】
【００２９】
【数３】
【００３０】
【数４】
【００３１】
【数５】
次に、上記の各スペクトルを用いて、下記の数６により伝達関数Ｈ（ｆ）を、数７によりコヒーレンス関数γ２ （ｆ）をそれぞれ算出し（ステップＳ３、Ｓ４）、これらをＰＣ１３内蔵のハードディスクに保存する。
【００３２】
【数６】
【００３３】
【数７】
なお、コヒーレンス関数γ２ （ｆ）は、外乱ノイズの影響を見るために算出されるものである。すなわち、コヒーレンス関数γ２ （ｆ）は、０〜１の範囲の値をとるが、外乱ノイズレベルが大きくなってくるとその値が０に近づいてくる。そこで、コヒーレンス関数γ２ （ｆ）の値にしきい値を設定しておき、算出されたコヒーレンス関数γ２ （ｆ）の値がこのしきい値よりも小さくなった場合には算出された伝達関数Ｈ（ｆ）を採用しないようにされる。
【００３４】
続いて、伝達関数Ｈ（ｆ）の実数成分Ｒｅ（Ｈ（ｆ））と虚数成分Ｉｍ（Ｈ （ｆ））とを用いて、下記の数８により２つの信号間の位相差Δθ（ｆ）を求め、更に下記の数９により２つの信号間の時間遅れΔｔ（ｆ）に変換する（ステップＳ５）。
【００３５】
【数８】
【００３６】
【数９】
そして、時間遅れΔｔ（ｆ）と間隔Ｌとを用いて、下記の数１０、数１１により伝播平均速度Ｖｒ （ｍ／ｓｅｃ）と深度Ｄ（ｍ）とをそれぞれ計算する（ステップＳ６）。
【００３７】
【数１０】
【００３８】
【数１１】
なお、上記式中で、記号＊は複素共役を表す。
【００３９】
上記の算出は所望のＤ−Ｖｒ 曲線が得られるまで（ステップＳ７）繰り返し行われる。すなわち、計測を１回行う毎に起振機１５に与える励振信号の周波数を変化させる。つまり、表面波の伝播速度と逆数関係にある伝達関数Ｈ（ｆ）の位相差が、周波数別に計測される。次に、この位相差と周波数の関係から平均速度Ｖｒ と深度Ｄとを計算し、繰り返し計測を行った結果としてモニター上にＤ−Ｖｒ 分散曲線が表示される。
【００４０】
特に、ステップサイン計測の場合には、起振機１５にある周波数ｆｎ （図２ａ参照）の励振信号を加えて求められた伝播平均速度Ｖｒ （ｆｎ ）と深度Ｄ（ｆｎ ）で決まる値Ｐｎ をＤ−Ｖｒ 面上にプロットし、次に周波数ｆｎ−１ （図２ａ参照）の励振信号を加えて求められた伝播平均速度Ｖｒ （ｆｎ−１ ）と深度Ｄ（ｆｎ−１ ）で決まる値Ｐｎ−１ をＤ−Ｖｒ 面上にプロットする。これを繰り返し、最後に周波数ｆ１ （図２ａ参照）の励振信号を加えて求められた伝播平均速度Ｖｒ （ｆ１ ）と深度Ｄ（ｆ１ ）で決まる値Ｐ１ をＤ−Ｖｒ 面上にプロットすることで、図４に示すようなＤ−Ｖｒ 分散曲線が得られる。
【００４１】
一方、マルチサイン計測の場合には、起振機１５に、図２（ｂ）に示したようなマルチサイン信号あるいは図２（ｃ）に示したようなスウェプトサイン信号を加えると、合成波による振動に起因する時系列信号Ａ（ｔ）、Ｂ（ｔ）が得られる。これらの時系列信号は、図２（ａ）に示されるような周波数成分を含んでおり、これらの周波数成分をＦＦＴにより分離、抽出する。そして、抽出された周波数成分のそれぞれに対して上記のような処理を行うことで、図４に示すような値Ｐｎ 〜Ｐ１ が一度に得られる。
【００４２】
上記の説明で明らかなように、マルチサイン計測の場合には、Ｄ−Ｖｒ 分散曲線を得るまでの計測処理時間が短くて済む反面、合成波であるために振動のエネルギーが小さくなり、ノイズの影響を受け易くなる。また、図２（ｃ）に示したようなスウェプトサイン信号を印加する場合には、土中に存在する石やその他の物質により共振を起こすことがあるので、マルチサイン信号との切替えを行いながら計測を行う必要がある。
【００４３】
なお、振動エネルギーの弱小化の問題は、図５（ａ）に示すように、周波数ｆ１ 〜ｆｎ の領域を、例えばｆ１ 〜ｆｉ （バンドＢ１ ）、ｆｉ 〜ｆｋ （バンドＢ２ ）、ｆｋ 〜ｆｎ （バンドＢ３ ）の３つのバンドに分割してそれぞれのバンドにおいて合成波によるマルチサイン計測を行うことで解消できる。この場合、図５（ｂ）に示すように、バンドＢ１ での計測でＤ−Ｖｒ 面には曲線Ｃ１ が、バンドＢ２ での計測で曲線Ｃ２ が、バンドＢ３ での計測で曲線Ｃ３ がそれぞれ得られ、これらの曲線Ｃ１ 〜Ｃ３ は自動合成表示される。このようにするのは、合成される周波数信号が少ないほど、振動エネルギーは大きくなるからである。勿論、バンド数は２でも良い。
【００４４】
図６、図７はそれぞれ、ステップサイン計測、マルチサイン計測のモニターの表示例を示し、ここではＤ−Ｖｒ 分散曲線の他、コヒーレンス関数γ２ （ｆ）、伝達関数Ｈ（ｆ）のゲイン特性、伝達関数Ｈ（ｆ）の位相特性、リニアスペクトルがそれぞれマルチ画面形式で表示されている。但し、図６、図７のいずれにおいてもコヒーレンス関数γ２ （ｆ）はほぼ１に近い値がプロットされているので見えにくくなっている。また、リニアスペクトルというのは、パワースペクトルの平方根であり、表示画面で見る場合にはパワースペクトルよりもリニアスペクトルの方が見易いことを考慮したものである。
【００４５】
次に、ステップサイン計測と、マルチサイン計測について比較説明する。
【００４６】
ステップサイン計測では、前に述べたように、起振機１５の励振信号としてサイン波が使用され、解析プログラムによる解析ではＤＦＴが用いられる。ステップサイン計測の特徴は以下の通りである。
【００４７】
Ａ．ＤＦＴを用いると、任意のサイン波周波数を設定して計測を行うことができる。つまり、周波数バンド、分解能の制約が無く、任意の発振周波数を設定できる。周波数設定の制約が無ければ任意の発振周波数を設定してデータ収集を行うことができるので、モニターに表示されたＤ−Ｖｒ 分散曲線を観察しながら、計測値のプロットが疎となっている曲線領域について周波数を追加設定することで追加の計測を容易に行うことができる。
【００４８】
Ｂ．浅い方から深い深度にかけて探査を行うことができる。つまり、Ｄ−Ｖｒ 分散曲線を観察しながら、浅い方から深い深度にかけて探査を行い、目的の探査深度で計測を終了できる。
【００４９】
これは、本形態では周波数を低い方から高い方、高い方から低い方への両方向に掃引しながらの計測（ステップサイン計測）が可能であることによる。従って、周波数を高い方から低い方向に掃引しながら計測することにより、浅い方から深い深度に向かっての計測が可能となる。
【００５０】
Ｃ．以前計測した周波数と同じ周波数での計測が可能である。つまり、同じ地層構造の地点での周波数設定を簡略化できる。
【００５１】
これは、以前に計測した周波数設定順序をファイルリストとしてＰＣ１３のハードディスクに記憶でき、その記憶された周波数リストを元に以前計測した時と同じ周波数の設定順序で計測が可能であることによる。
【００５２】
Ｄ．自動入力レンジでの計測が可能である。つまり、入力レンジ変更の手間が省ける。
【００５３】
これは、本形態では、Ａ／Ｄ変換器１２−２に入力感度の自動調整機能を持たせていることによる。Ａ／Ｄ変換器１２−２では、入力レベルをあらかじめ設定した上下限値について監視し、下限値以下になると入力感度を上げ、上限値を越えると入力感度を下げる操作を自動的に行う。これにより、最適入力レベルで計測を行うことができる。但し、これは外来ノイズレベルの変動が大きい場所では不向きであり、計測時間が長くなる場合があり得る。そのような場合は手動レンジに切替えるようにする。
【００５４】
Ｅ．Ｄ−Ｖｒ 分散曲線に自動変換表示できる。つまり、毎回、伝達関数位相差曲線からＤ−Ｖｒ 変換を行うことが不要となり、手間が省け、しかも図６に示す如く最終結果であるＤ−Ｖｒ 分散曲線をリアルタイムで観測できる。
【００５５】
これは、本形態ではデータが更新されるたびにＤ−Ｖｒ 分散曲線が自動的に変換表示されるからである。
【００５６】
Ｆ．計測点の振動レベルを加速度（Ｇａｌ）で表示できる。つまり、図６の右下部分の画面領域に示されるように、実レベルの数値が読み取れる。図６で言えば、図１に示された起振機１５に近い方の加速度検出器１１ＡによるパワースペクトルＧａａ：３．４２６、起振機１５から遠い方の加速度検出器１１ＢによるパワースペクトルＧｂｂ：３．０６４が示されている。
【００５７】
Ｇ．短い時間での計測が可能である。つまり、周波数バンド、分解能の制約が無く、効率のよいデータ収集ができるため、計測時間が短くできる。
【００５８】
これは、本形態ではＤＦＴを使っているため任意の周波数を設定して信号処理が行えるからである。従って、接近した同じ様なデータを何度も収集する必要がなく、計測時間を短縮できる可能性がある。
【００５９】
Ｈ．大部分の計測点でステップサイン計測だけでカバーできる。
【００６０】
一方、マルチサイン計測では、前に述べたように、起振機１５の励振信号としてマルチサイン波、あるいはスウェプトサイン波が使用され、解析プログラムによる解析ではＦＦＴが用いられる。マルチサイン計測の特徴は以下の通りである。
【００６１】
ａ．図５で説明したように、最大３バンドの計測を自動的に合成できる。従って、計測後のＤ−Ｖｒ 曲線の合成作業が省ける。
【００６２】
これは、本形態では最大３バンドの計測データまで自動的に合成表示する機能を有するからである。そして、１バンドの周波数分解能が８００点になっており、きめ細かいＤ−Ｖｒ 分散曲線が計測でき、より精密な計測が出来る。
【００６３】
ｂ．周波数範囲を限定してマルチサイン波を周期性ランダム信号として発生できる。従って、バンド範囲を限定することにより、起振機の能力を向上できる。
【００６４】
これは、本形態では周波数範囲を限定できるために個々の周波数スペクトル値を大きくして周期性ランダム信号を発生できるからである。
【００６５】
ｃ．マルチサイン波、スウェプトサイン波の切替えが簡単に行える。
【００６６】
ｄ．ステップサイン計測と同様に、Ｄ−Ｖｒ 分散曲線に自動変換表示できる。
【００６７】
ｅ．ステップサイン計測と同様に、計測点の振動レベルを加速度（Ｇａｌ）で表示できる。
【００６８】
【発明の効果】
以上説明したように、本発明による地盤探査装置によれば、非破壊で正確に迅速でしかも経済的に探査を行うことのできる地盤探査装置を提供できる。
【図面の簡単な説明】
【図１】本発明の好ましい実施の形態による地盤探査装置の構成を示した図である。
【図２】本発明のマルチサイン計測の場合に発生されるマルチサイン波（ｂ）、スウェプトサイン波（ｃ）を説明するための波形図である。
【図３】本発明で使用される解析プログラムの動作を説明するためのフローチャート図である。
【図４】本発明のステップサイン計測におけるＤ−Ｖｒ 分散曲線の描画過程を説明するための図である。
【図５】本発明のマルチサイン計測において周波数バンドを３つに分割してＤ−Ｖｒ 分散曲線を得る場合の動作を説明するための図である。
【図６】本発明の地盤探査装置におけるモニターに表示される表示画面の一例を、ステップサイン計測の場合について示した図である。
【図７】本発明の地盤探査装置におけるモニターに表示される表示画面の一例を、マルチサイン計測の場合について示した図である。
【図８】半無限弾性体におけるＰ波、Ｓ波、レイリー波伝播速度とポアソン比とを関係を示した特性図である。
【符号の説明】
１１Ａ、１１Ｂ 加速度検出器
１２ 計測器
１３ パーソナルコンピュータ
【発明の属する技術分野】
本発明は、表面波、特にレイリー波の特性を利用した地盤探査装置に関する。
【０００２】
【従来の技術】
土木建設の分野における構造物基礎、トンネル、法面等の地盤調査には物理探査法が多く活用されている。なかでも弾性波探査法は地盤の工学的性質を定量的に把握する手法として用いられ、その代表的な例として広範囲に把握する屈折波法、局部的に把握する直接波法（ボーリング孔等を利用する）が知られている。これらの両手法で取り扱う波は人工振源で生じる弾性波動のうち実体波と呼ばれるＰ波及びＳ波であり、地質条件等を考慮して両波が使い分けされている。
【０００３】
これらの実体波とは別に、表面波と呼ばれる波の存在が知られており、表面波は分散現象を伴い、その伝播速度は波長に依存する性質を持っていることも知られている。以下に、表面波の性質について簡単に説明する。
【０００４】
起振機で地表面に上下振動を与えると、実体波（Ｐ波、Ｓ波）以外に表面波と呼ばれる波が地盤中に発生する。このうち実体波は、振源から半球状に伝播していく。これに対して、表面波のうちレイリー波は、周波数に規制された高さ一定の円筒が径を広げていくように伝播していく。つまり、起振機より発生する実体波には指向性があり、粗密波（Ｐ波）は起振機直下、せん断波（Ｓ波）は４５度方向に伝わるエネルギーが最大で、水平方向には伝わりにくい。このため、地表面で大振幅を与えるのは唯一表面波のみであると言える。
【０００５】
実体波（Ｐ波、Ｓ波）の幾何学的減衰係数は、下方に対してはｒ−１（但し、ｒは距離）であるが、地表付近ではｒ−２である。
【０００６】
これに対して表面波の地表面での幾何学的減衰係数は、ｒ−０．５である。
【０００７】
表面波にはレイリー波の他にラブ波、Ｐ波に類似するＰ−Ｐモード波、板波等があるが、本発明ではレイリー波を扱う。
【０００８】
レイリー波速度とＰ波及びＳ波速度の関係は、図８に示す通りであり、これらの３種類の速度比はポアソン比νで異なるが、Ｓ波速度とレイリー波速度は近似した値を示す。図８において、縦軸はＶ／ＶＳ ＝Ｖ・（ρ／Ｇ）１／２ で表される値である。但し、ＶＳ はＳ波速度、ρは密度、Ｇは剛性率であり、ＶにはＰ波あるいはレイリー波速度が当てはまる。
【０００９】
レイリー波速度はＳ波速度に近似しており、土質工学的数値を良く表していることからＳ波速度とＮ値（標準貫入試験値）の相関式やｑａ値（許容支持力）との相関式を用いてそれぞれの土質工学的数値を推定することができる。
【００１０】
【発明が解決しようとする課題】
上記のような知見に基づいて、本発明は、表面波、特にレイリー波の分散性を利用し、非破壊で正確に迅速でしかも経済的に探査を行うことのできる地盤探査装置を提供しようとするものである。
【００１１】
【課題を解決するための手段】
本発明によれば、地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査装置において、間隔Ｌをおいて地盤上に配置された２つの加速度検出器と、前記２つの加速度検出器からの検出信号を受けて加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）を出力する地震計部を含む計測器と、前記時系列信号Ａ（ｔ）、Ｂ（ｔ）を受けてあらかじめ定められた解析プログラムに基づく信号処理を行う信号処理部であって、フーリエ変換を行ってパワースペクトル及びクロススペクトルを算出すると共に、算出されたパワースペクトル及びクロススペクトルを用いて伝達関数Ｈ（ｆ）を算出し、算出された伝達関数Ｈ（ｆ）より前記２つの加速度検出器間の位相差Δθ（ｆ）とその時間差Δｔ（ｆ）を算出し、更に算出された時間差Δｔ（ｆ）と前記間隔Ｌより表面波の伝播平均速度Ｖｒ（ｆ）と深度Ｄ（ｆ）とを算出する信号処理部とを備え、前記計測器は更に、地盤を起振する起振機に対して与える励振信号を発生する発振部を含み、前記発振部は、励振信号として、ステップサイン波、又はマルチサイン波・スウェプトサイン波を切り替えて発生することが可能であり、前記ステップサイン波は周波数の異なる複数種類の正弦波で規定され、前記マルチサイン波は周波数の異なる複数種類の正弦波を合成した合成波で規定され、前記スウェプトサイン波は周波数の異なる複数種類の正弦波を振幅が一定の状態で合成した合成波で規定され、前記マルチサイン波・スウェプトサイン波を使用する場合には合成される周波数領域ｆ １ 〜ｆ ｎ を複数のバンドに分割し、分割された各バンドにおいて前記マルチサイン波・スウェプトサイン波による計測を行うことを特徴とする地盤探査装置が提供される。
【００１４】
本地盤探査装置においては更に、前記信号処理部は、前記伝播平均速度Ｖｒ （ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝播平均速度Ｖｒ 曲線を生成する。
【００１５】
本地盤探査装置においては更に、前記計測器が更に、前記地震計部からのアナログの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）をディジタルの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）に変換して出力するＡ／Ｄ変換部を含み、該Ａ／Ｄ変換部は、入力レベルに応じて入力感度を自動調整する機能を有する。
【００１６】
本地盤探査装置においては更に、前記信号処理部は、前記励振信号として前記ステップサイン波が使用される場合には、ＤＦＴ（離散フーリエ変換）による処理を行い、前記励振信号として前記マルチサイン波・スウェプトサイン波が使用される場合には、ＦＦＴ（高速フーリエ変換）による処理を行う。
【００１９】
【発明の実施の形態】
図１を参照して、本発明の好ましい実施の形態による地盤探査装置の構成について説明する。図１において、本装置は、２つの加速度検出器１１Ａ、１１Ｂと、これらの加速度検出器１１Ａ、１１Ｂからの加速度検出信号を受ける計測器１２とを含む。計測器１２は、地震計部１２−１、Ａ／Ｄ変換部１２−２、通信部１２−３、発振部１２−４とから成る。地震計部１２−１は、ローパスフィルタ回路を内蔵し、アナログの加速度検出信号から加速度時系列信号を生成するためのものである。Ａ／Ｄ変換部１２−２は、地震計部１２−１からのアナログの加速度時系列信号をディジタルの加速度時系列信号に変換するためのもので、入力感度の自動調整機能を有する。通信部１２−３は、ディジタルの加速度時系列信号を計測器１２に接続される機器、例えばモニタ付きのパーソナルコンピュータ（以下、ＰＣと略称する）１３に送信するためのものである。ＰＣ１３は、あらかじめインストールされた解析プログラムソフトに基づいて信号処理を行い、地盤解析を行うためのものである。
【００２０】
特に、本発明ではこの信号処理にＦＦＴ（高速フーリエ変換）、ＤＦＴ（離散フーリエ変換）を採用しているが、これについては後述する。
【００２１】
発振部１２−４は、起振機１５に与える励振信号を発生するためのもので、この励振信号は電力増幅器１４で増幅されて起振機１５に与えられる。発振部１２−４は、励振信号として、サイン波信号、マルチサイン信号・スウェプトサイン信号を発生することができる。ここで、サイン波信号というのは、正弦波のことであり、周波数の異なる複数種類のサイン波信号を生成できる。マルチサイン信号というのは、図２（ａ）に示すような異なる周波数ｆ１ 〜ｆｎ のサイン波信号を様々な振幅を持つように合成したもので、一例をあげると図２（ｂ）のような合成波である。一方、スウェプトサイン信号というのは、異なる周波数ｆ１ 〜ｆｎ のサイン波信号を振幅が一定の状態になるように合成したもので、一例をあげると図２（ｃ）のような合成波、いわば周波数変調波である。なお、図２（ｂ）、図２（ｃ）は便宜上の合成波形であり、図２（ａ）のサイン信号の実際の合成波形を示すものでは無い。
【００２２】
後述する説明で明らかとなるように、本形態では、励振信号としてサイン波信号を用いる場合には信号処理は上記のＤＦＴにより行われる（以下では、これをステップサイン計測と呼ぶ）。一方、マルチサイン信号・スウェプトサイン信号を用いる場合には信号処理は上記のＦＦＴにより行われる（以下では、これをマルチサイン計測と呼ぶ）。ＦＦＴの場合、周波数は２ｎ で決まる固定値となり、ステップ幅も固定である。しかも、設定する周波数を低い方から高い方に変化させてゆく必要がある。これに対し、ＤＦＴの場合には、周波数を任意に設定することができ、設定する周波数も低い方から高い方、高い方から低い方への両方の変化が可能である。
【００２３】
以下に動作を説明する。まず、探査地点に起振機１５、加速度検出器１１Ａ、１１Ｂを一直線上に設置する。加速度検出器１１Ａ、１１Ｂ間の距離をＬ（ｍ）とする。起振機１５を使って地表面を上下方向に起振することにより、起振機１５の周りに表面波を発生させる。地表面付近を伝播する表面波（レイリー波）の上下振動を、加速度検出器１１Ａ、１１Ｂで検出する。加速度検出器１１Ａ、１１Ｂからの加速度検出信号は、地震計部１２−１のローパスフィルタ回路を通すことで、アナログの時系列信号Ａ（ｔ）、Ｂ（ｔ）となり、Ａ／Ｄ変換部１２−２に入力される。Ａ／Ｄ変換部１２−２でＡ／Ｄ変換された時系列信号Ａ（ｔ）、Ｂ（ｔ）は通信部１２−３からＰＣ１３へ転送される。
【００２４】
ＰＣ１３では、あらかじめ定められた解析プログラムに基づいて時系列信号Ａ（ｔ）、Ｂ（ｔ）のパワースペクトルＧＡＡ（ｆ）、ＧＢＢ（ｆ）、クロススペクトルＧＢＡ（ｆ）、伝達関数Ｈ（ｆ）、コヒーレンス関数γ２ （ｆ）等を計算する。これらのパワースペクトルＧＡＡ（ｆ）、ＧＢＢ（ｆ）、クロススペクトルＧＢＡ（ｆ）、伝達関数Ｈ（ｆ）、コヒーレンス関数γ２ （ｆ）等は、内蔵のハードディスクに保存される。ＰＣ１３ではまた、伝達関数Ｈ（ｆ）より加速度検出器１１Ａ、１１Ｂ間の位相差Δθ（ｆ）を求め、続いてその時間差Δｔ（ｆ）を求める。ＰＣ１３では更に、時間差Δｔ（ｆ）と加速度検出器１１Ａ、１１Ｂ間の間隔Ｌより、表面波の伝播平均速度Ｖｒ （ｆ）と深度Ｄ（ｆ）とを求める。そして、求められた伝播平均速度Ｖｒ （ｆ）と深度Ｄ（ｆ）とよりモニター上にＤ−Ｖｒ 分散曲線を表示させる。
【００２５】
図３をも参照して、解析プログラムに基づく、上記の表面波の伝播平均速度Ｖｒ （ｆ）と深度Ｄ（ｆ）の算出方法について更に詳しく説明する。
【００２６】
ＰＣ１３では、Ａ／Ｄ変換された時系列信号Ａ（ｔ）、Ｂ（ｔ）に対し、下記の数１、数２で表されるフーリエ変換を行って変換信号ＳＡ （ｆ）、ＳＢ （ｆ）を得た後（ステップＳ１）、これらの変換信号ＳＡ （ｆ）、ＳＢ （ｆ）を用いて下記の数３、数４、及び数５によりパワースペクトルＧＡＡ（ｆ）、ＧＢＢ（ｆ）及びクロススペクトルＧＢＡ（ｆ）を計算する（ステップＳ２）。
【００２７】
【数１】
【数２】
【数３】
【数４】
【数５】
【００３２】
【数６】
【数７】
【００３４】
続いて、伝達関数Ｈ（ｆ）の実数成分Ｒｅ（Ｈ（ｆ））と虚数成分Ｉｍ（Ｈ （ｆ））とを用いて、下記の数８により２つの信号間の位相差Δθ（ｆ）を求め、更に下記の数９により２つの信号間の時間遅れΔｔ（ｆ）に変換する（ステップＳ５）。
【００３５】
【数８】
【数９】
【００３７】
【数１０】
【数１１】
【００３９】
上記の算出は所望のＤ−Ｖｒ 曲線が得られるまで（ステップＳ７）繰り返し行われる。すなわち、計測を１回行う毎に起振機１５に与える励振信号の周波数を変化させる。つまり、表面波の伝播速度と逆数関係にある伝達関数Ｈ（ｆ）の位相差が、周波数別に計測される。次に、この位相差と周波数の関係から平均速度Ｖｒ と深度Ｄとを計算し、繰り返し計測を行った結果としてモニター上にＤ−Ｖｒ 分散曲線が表示される。
【００４０】
特に、ステップサイン計測の場合には、起振機１５にある周波数ｆｎ （図２ａ参照）の励振信号を加えて求められた伝播平均速度Ｖｒ （ｆｎ ）と深度Ｄ（ｆｎ ）で決まる値Ｐｎ をＤ−Ｖｒ 面上にプロットし、次に周波数ｆｎ−１ （図２ａ参照）の励振信号を加えて求められた伝播平均速度Ｖｒ （ｆｎ−１ ）と深度Ｄ（ｆｎ−１ ）で決まる値Ｐｎ−１ をＤ−Ｖｒ 面上にプロットする。これを繰り返し、最後に周波数ｆ１ （図２ａ参照）の励振信号を加えて求められた伝播平均速度Ｖｒ （ｆ１ ）と深度Ｄ（ｆ１ ）で決まる値Ｐ１ をＤ−Ｖｒ 面上にプロットすることで、図４に示すようなＤ−Ｖｒ 分散曲線が得られる。
【００４１】
一方、マルチサイン計測の場合には、起振機１５に、図２（ｂ）に示したようなマルチサイン信号あるいは図２（ｃ）に示したようなスウェプトサイン信号を加えると、合成波による振動に起因する時系列信号Ａ（ｔ）、Ｂ（ｔ）が得られる。これらの時系列信号は、図２（ａ）に示されるような周波数成分を含んでおり、これらの周波数成分をＦＦＴにより分離、抽出する。そして、抽出された周波数成分のそれぞれに対して上記のような処理を行うことで、図４に示すような値Ｐｎ 〜Ｐ１ が一度に得られる。
【００４２】
上記の説明で明らかなように、マルチサイン計測の場合には、Ｄ−Ｖｒ 分散曲線を得るまでの計測処理時間が短くて済む反面、合成波であるために振動のエネルギーが小さくなり、ノイズの影響を受け易くなる。また、図２（ｃ）に示したようなスウェプトサイン信号を印加する場合には、土中に存在する石やその他の物質により共振を起こすことがあるので、マルチサイン信号との切替えを行いながら計測を行う必要がある。
【００４３】
なお、振動エネルギーの弱小化の問題は、図５（ａ）に示すように、周波数ｆ１ 〜ｆｎ の領域を、例えばｆ１ 〜ｆｉ （バンドＢ１ ）、ｆｉ 〜ｆｋ （バンドＢ２ ）、ｆｋ 〜ｆｎ （バンドＢ３ ）の３つのバンドに分割してそれぞれのバンドにおいて合成波によるマルチサイン計測を行うことで解消できる。この場合、図５（ｂ）に示すように、バンドＢ１ での計測でＤ−Ｖｒ 面には曲線Ｃ１ が、バンドＢ２ での計測で曲線Ｃ２ が、バンドＢ３ での計測で曲線Ｃ３ がそれぞれ得られ、これらの曲線Ｃ１ 〜Ｃ３ は自動合成表示される。このようにするのは、合成される周波数信号が少ないほど、振動エネルギーは大きくなるからである。勿論、バンド数は２でも良い。
【００４４】
図６、図７はそれぞれ、ステップサイン計測、マルチサイン計測のモニターの表示例を示し、ここではＤ−Ｖｒ 分散曲線の他、コヒーレンス関数γ２ （ｆ）、伝達関数Ｈ（ｆ）のゲイン特性、伝達関数Ｈ（ｆ）の位相特性、リニアスペクトルがそれぞれマルチ画面形式で表示されている。但し、図６、図７のいずれにおいてもコヒーレンス関数γ２ （ｆ）はほぼ１に近い値がプロットされているので見えにくくなっている。また、リニアスペクトルというのは、パワースペクトルの平方根であり、表示画面で見る場合にはパワースペクトルよりもリニアスペクトルの方が見易いことを考慮したものである。
【００４５】
次に、ステップサイン計測と、マルチサイン計測について比較説明する。
【００４６】
ステップサイン計測では、前に述べたように、起振機１５の励振信号としてサイン波が使用され、解析プログラムによる解析ではＤＦＴが用いられる。ステップサイン計測の特徴は以下の通りである。
【００４７】
Ａ．ＤＦＴを用いると、任意のサイン波周波数を設定して計測を行うことができる。つまり、周波数バンド、分解能の制約が無く、任意の発振周波数を設定できる。周波数設定の制約が無ければ任意の発振周波数を設定してデータ収集を行うことができるので、モニターに表示されたＤ−Ｖｒ 分散曲線を観察しながら、計測値のプロットが疎となっている曲線領域について周波数を追加設定することで追加の計測を容易に行うことができる。
【００４８】
Ｂ．浅い方から深い深度にかけて探査を行うことができる。つまり、Ｄ−Ｖｒ 分散曲線を観察しながら、浅い方から深い深度にかけて探査を行い、目的の探査深度で計測を終了できる。
【００４９】
これは、本形態では周波数を低い方から高い方、高い方から低い方への両方向に掃引しながらの計測（ステップサイン計測）が可能であることによる。従って、周波数を高い方から低い方向に掃引しながら計測することにより、浅い方から深い深度に向かっての計測が可能となる。
【００５０】
Ｃ．以前計測した周波数と同じ周波数での計測が可能である。つまり、同じ地層構造の地点での周波数設定を簡略化できる。
【００５１】
これは、以前に計測した周波数設定順序をファイルリストとしてＰＣ１３のハードディスクに記憶でき、その記憶された周波数リストを元に以前計測した時と同じ周波数の設定順序で計測が可能であることによる。
【００５２】
Ｄ．自動入力レンジでの計測が可能である。つまり、入力レンジ変更の手間が省ける。
【００５３】
これは、本形態では、Ａ／Ｄ変換器１２−２に入力感度の自動調整機能を持たせていることによる。Ａ／Ｄ変換器１２−２では、入力レベルをあらかじめ設定した上下限値について監視し、下限値以下になると入力感度を上げ、上限値を越えると入力感度を下げる操作を自動的に行う。これにより、最適入力レベルで計測を行うことができる。但し、これは外来ノイズレベルの変動が大きい場所では不向きであり、計測時間が長くなる場合があり得る。そのような場合は手動レンジに切替えるようにする。
【００５４】
Ｅ．Ｄ−Ｖｒ 分散曲線に自動変換表示できる。つまり、毎回、伝達関数位相差曲線からＤ−Ｖｒ 変換を行うことが不要となり、手間が省け、しかも図６に示す如く最終結果であるＤ−Ｖｒ 分散曲線をリアルタイムで観測できる。
【００５５】
これは、本形態ではデータが更新されるたびにＤ−Ｖｒ 分散曲線が自動的に変換表示されるからである。
【００５６】
Ｆ．計測点の振動レベルを加速度（Ｇａｌ）で表示できる。つまり、図６の右下部分の画面領域に示されるように、実レベルの数値が読み取れる。図６で言えば、図１に示された起振機１５に近い方の加速度検出器１１ＡによるパワースペクトルＧａａ：３．４２６、起振機１５から遠い方の加速度検出器１１ＢによるパワースペクトルＧｂｂ：３．０６４が示されている。
【００５７】
Ｇ．短い時間での計測が可能である。つまり、周波数バンド、分解能の制約が無く、効率のよいデータ収集ができるため、計測時間が短くできる。
【００５８】
これは、本形態ではＤＦＴを使っているため任意の周波数を設定して信号処理が行えるからである。従って、接近した同じ様なデータを何度も収集する必要がなく、計測時間を短縮できる可能性がある。
【００５９】
Ｈ．大部分の計測点でステップサイン計測だけでカバーできる。
【００６０】
一方、マルチサイン計測では、前に述べたように、起振機１５の励振信号としてマルチサイン波、あるいはスウェプトサイン波が使用され、解析プログラムによる解析ではＦＦＴが用いられる。マルチサイン計測の特徴は以下の通りである。
【００６１】
ａ．図５で説明したように、最大３バンドの計測を自動的に合成できる。従って、計測後のＤ−Ｖｒ 曲線の合成作業が省ける。
【００６２】
これは、本形態では最大３バンドの計測データまで自動的に合成表示する機能を有するからである。そして、１バンドの周波数分解能が８００点になっており、きめ細かいＤ−Ｖｒ 分散曲線が計測でき、より精密な計測が出来る。
【００６３】
ｂ．周波数範囲を限定してマルチサイン波を周期性ランダム信号として発生できる。従って、バンド範囲を限定することにより、起振機の能力を向上できる。
【００６４】
これは、本形態では周波数範囲を限定できるために個々の周波数スペクトル値を大きくして周期性ランダム信号を発生できるからである。
【００６５】
ｃ．マルチサイン波、スウェプトサイン波の切替えが簡単に行える。
【００６６】
ｄ．ステップサイン計測と同様に、Ｄ−Ｖｒ 分散曲線に自動変換表示できる。
【００６７】
ｅ．ステップサイン計測と同様に、計測点の振動レベルを加速度（Ｇａｌ）で表示できる。
【００６８】
【発明の効果】
以上説明したように、本発明による地盤探査装置によれば、非破壊で正確に迅速でしかも経済的に探査を行うことのできる地盤探査装置を提供できる。
【図面の簡単な説明】
【図１】本発明の好ましい実施の形態による地盤探査装置の構成を示した図である。
【図２】本発明のマルチサイン計測の場合に発生されるマルチサイン波（ｂ）、スウェプトサイン波（ｃ）を説明するための波形図である。
【図３】本発明で使用される解析プログラムの動作を説明するためのフローチャート図である。
【図４】本発明のステップサイン計測におけるＤ−Ｖｒ 分散曲線の描画過程を説明するための図である。
【図５】本発明のマルチサイン計測において周波数バンドを３つに分割してＤ−Ｖｒ 分散曲線を得る場合の動作を説明するための図である。
【図６】本発明の地盤探査装置におけるモニターに表示される表示画面の一例を、ステップサイン計測の場合について示した図である。
【図７】本発明の地盤探査装置におけるモニターに表示される表示画面の一例を、マルチサイン計測の場合について示した図である。
【図８】半無限弾性体におけるＰ波、Ｓ波、レイリー波伝播速度とポアソン比とを関係を示した特性図である。
【符号の説明】
１１Ａ、１１Ｂ 加速度検出器
１２ 計測器
１３ パーソナルコンピュータ
Claims (4)
- 地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査装置において、
間隔Ｌをおいて地盤上に配置された２つの加速度検出器と、
前記２つの加速度検出器からの検出信号を受けて加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）を出力する地震計部を含む計測器と、
前記時系列信号Ａ（ｔ）、Ｂ（ｔ）を受けてあらかじめ定められた解析プログラムに基づく信号処理を行う信号処理部であって、フーリエ変換を行ってパワースペクトル及びクロススペクトルを算出すると共に、算出されたパワースペクトル及びクロススペクトルを用いて伝達関数Ｈ（ｆ）を算出し、算出された伝達関数Ｈ（ｆ）より前記２つの加速度検出器間の位相差Δθ（ｆ）とその時間差Δｔ（ｆ）を算出し、更に算出された時間差Δｔ（ｆ）と前記間隔Ｌより表面波の伝播平均速度Ｖｒ（ｆ）と深度Ｄ（ｆ）とを算出する信号処理部とを備え、
前記計測器は更に、地盤を起振する起振機に対して与える励振信号を発生する発振部を含み、
前記発振部は、励振信号として、ステップサイン波、又はマルチサイン波・スウェプトサイン波を切り替えて発生することが可能であり、
前記ステップサイン波は周波数の異なる複数種類の正弦波で規定され、前記マルチサイン波は周波数の異なる複数種類の正弦波を合成した合成波で規定され、前記スウェプトサイン波は周波数の異なる複数種類の正弦波を振幅が一定の状態で合成した合成波で規定され、
前記マルチサイン波・スウェプトサイン波を使用する場合には合成される周波数領域ｆ １ 〜ｆ ｎ を複数のバンドに分割し、分割された各バンドにおいて前記マルチサイン波・スウェプトサイン波による計測を行うことを特徴とする地盤探査装置。 - 請求項１記載の地盤探査装置において、前記信号処理部は、前記伝播平均速度Ｖｒ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝播平均速度Ｖｒ曲線を生成することを特徴とする地盤探査装置。
- 請求項２記載の地盤探査装置において、前記計測器は更に、前記地震計部からのアナログの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）をディジタルの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）に変換して出力するＡ／Ｄ変換部を含み、該Ａ／Ｄ変換部は、入力レベルに応じて入力感度を自動調整する機能を有することを特徴とする地盤探査装置。
- 請求項３記載の地盤探査装置において、前記信号処理部は、前記励振信号として前記ステップサイン波が使用される場合には、ＤＦＴ（離散フーリエ変換）による処理を行い、前記励振信号として前記マルチサイン波・スウェプトサイン波が使用される場合には、ＦＦＴ（高速フーリエ変換）による処理を行うことを特徴とする地盤探査装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001151961A JP3617036B2 (ja) | 2001-05-22 | 2001-05-22 | 地盤探査装置及びそれに使用される解析プログラム |
| CNB021201188A CN1194232C (zh) | 2001-05-22 | 2002-05-17 | 地面分析系统 |
| US10/152,287 US6665617B2 (en) | 2001-05-22 | 2002-05-21 | Ground analyzing system and recording medium recording analysis program used therefor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001151961A JP3617036B2 (ja) | 2001-05-22 | 2001-05-22 | 地盤探査装置及びそれに使用される解析プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002341048A JP2002341048A (ja) | 2002-11-27 |
| JP3617036B2 true JP3617036B2 (ja) | 2005-02-02 |
Family
ID=18996705
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001151961A Expired - Lifetime JP3617036B2 (ja) | 2001-05-22 | 2001-05-22 | 地盤探査装置及びそれに使用される解析プログラム |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6665617B2 (ja) |
| JP (1) | JP3617036B2 (ja) |
| CN (1) | CN1194232C (ja) |
Families Citing this family (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102006009246B3 (de) * | 2006-02-28 | 2007-08-02 | GeoForschungsZentrum Potsdam Stiftung des öffentlichen Rechts | Verfahren und Vorrichtung zur Vorauserkundung beim Tunnelbau |
| JP4913544B2 (ja) * | 2006-10-30 | 2012-04-11 | 曙ブレーキ工業株式会社 | 衝突検知方法 |
| WO2008147809A1 (en) | 2007-05-24 | 2008-12-04 | Schlumberger Canada Limited | Near surface layer modeling |
| US8852128B2 (en) * | 2008-03-12 | 2014-10-07 | University Of Cincinnati | Computer system and method for assessing dynamic bone quality |
| WO2009120401A1 (en) * | 2008-03-28 | 2009-10-01 | Exxonmobil Upstream Research Company | Characterizing spatial variablility of surface waves in seismic processing |
| US8615362B2 (en) | 2008-10-10 | 2013-12-24 | Westerngeco L.L.C. | Near-surface geomorphological characterization based on remote sensing data |
| US8509027B2 (en) | 2008-11-26 | 2013-08-13 | Westerngeco L.L.C. | Continuous adaptive surface wave analysis for three-dimensional seismic data |
| MX2012003975A (es) | 2009-10-08 | 2012-06-25 | Geco Technology Bv | Interpretacion conjunta de las ondas de rayleigh y deteccion remota para geologia de superficie cercana. |
| CN102096093B (zh) * | 2010-11-29 | 2013-04-10 | 杨本才 | 一种利用微震点作为震源计算矿区地震波传播速度的方法 |
| TWI420516B (zh) | 2012-03-19 | 2013-12-21 | Wistron Corp | 硬碟讀寫率與振動源關係判斷方法及其系統 |
| US9239220B2 (en) | 2012-04-30 | 2016-01-19 | Conocophillips Company | Determination of near surface geophyscial properties by impulsive displacement events |
| WO2015106065A1 (en) * | 2014-01-10 | 2015-07-16 | Cgg Services (U.S.) Inc. | Device and method for mitigating cycle-skipping in full waveform inversion |
| US10061044B2 (en) * | 2015-02-10 | 2018-08-28 | Transcend Engineering and Technology, LLC | Systems, methods, and software for detecting the presence of subterranean tunnels and tunneling activity |
| CN105319577A (zh) * | 2015-10-20 | 2016-02-10 | 北京市水电物探研究所 | 一种天然源面波的勘探仪器及其使用方法 |
| CN106019372B (zh) * | 2016-05-26 | 2018-01-05 | 重庆大学 | 地震波传播速度获取方法、测深方法及地震波观测系统 |
| JP6843427B2 (ja) * | 2016-11-07 | 2021-03-17 | ビイック株式会社 | 表面波探査解析方法及び表面波探査解析装置 |
| JP6543298B2 (ja) * | 2016-11-29 | 2019-07-10 | パナソニックホームズ株式会社 | 振動波形データ間の同期ずれ量計算方法、及び、振動波形データ間の同期補正方法 |
| CN106768296B (zh) * | 2016-12-26 | 2019-04-19 | 中北大学 | 一种应用于地下挖掘系统的区分地面震动信号干扰的方法 |
| JP6884341B2 (ja) * | 2017-03-30 | 2021-06-09 | 公立大学法人大阪 | 土質判定方法 |
| EP3635185A1 (en) | 2017-06-09 | 2020-04-15 | Volvo Construction Equipment AB | An information system for a working machine |
| CN108254782B (zh) * | 2018-02-09 | 2019-11-05 | 中国地质大学（北京） | 一种边坡地震破坏失稳概率的获取方法及系统 |
| RU2691208C1 (ru) * | 2018-08-15 | 2019-06-11 | Федеральное государственное бюджетное учреждение науки Институт вычислительных технологий Сибирского отделения Российской академии наук (ИВТ СО РАН) | Способ оценки состояния контакта фундаментной плиты строящегося здания с грунтовым основанием |
| CN109871656B (zh) * | 2019-03-20 | 2023-04-25 | 东北大学 | 一种基于动态强度折减dda法的边坡稳定性分析系统 |
| CN110133665B (zh) * | 2019-05-16 | 2021-01-01 | 浙江大学 | 一种超声换能器基于包络线的多普勒测量方法 |
| WO2022241555A1 (en) * | 2021-05-19 | 2022-11-24 | University Of Manitoba | Systems and methods for in-situ characterization of permafrost sites |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4750157A (en) * | 1987-05-06 | 1988-06-07 | Standard Oil Production Company | Seismic vibrator earth impedance determination and compensation system |
| US4933911A (en) * | 1989-09-01 | 1990-06-12 | Amoco Corporation | Method for determining seismic velocities |
| JP2958362B2 (ja) * | 1990-04-28 | 1999-10-06 | 孝次 時松 | 地盤構造の計測解析判定方法 |
| US5243565A (en) * | 1991-08-13 | 1993-09-07 | University Of Miami | Method of measuring directional spectra of surface gravity waves |
| US5490062A (en) * | 1994-05-11 | 1996-02-06 | The Regents Of The University Of California | Real-time neural network earthquake profile predictor |
| US5642327A (en) * | 1994-06-08 | 1997-06-24 | Exxon Production Research Company | Method for creating a gain function for seismic data and method for processing seismic data |
| US6536553B1 (en) * | 2000-04-25 | 2003-03-25 | The United States Of America As Represented By The Secretary Of The Army | Method and apparatus using acoustic sensor for sub-surface object detection and visualization |
-
2001
- 2001-05-22 JP JP2001151961A patent/JP3617036B2/ja not_active Expired - Lifetime
-
2002
- 2002-05-17 CN CNB021201188A patent/CN1194232C/zh not_active Expired - Fee Related
- 2002-05-21 US US10/152,287 patent/US6665617B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US6665617B2 (en) | 2003-12-16 |
| JP2002341048A (ja) | 2002-11-27 |
| US20020193952A1 (en) | 2002-12-19 |
| CN1194232C (zh) | 2005-03-23 |
| CN1387051A (zh) | 2002-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3617036B2 (ja) | 地盤探査装置及びそれに使用される解析プログラム | |
| US11940580B2 (en) | Heterogeneous subsurface imaging systems and methods | |
| US20140104979A1 (en) | Ground-Penetrating Tunnel-Detecting Active Sonar | |
| US11841265B2 (en) | Heterogeneous subsurface imaging systems and methods | |
| CN106949953B (zh) | 连续原位测量岩土剪切波速测试仪及其测试方法 | |
| Chan | Variations of shear wave arrival time in unconfined soil specimens measured with bender elements | |
| Wang et al. | Estimating near‐surface rigidity from low‐frequency noise using collocated pressure and horizontal seismic data | |
| JP2862171B2 (ja) | 地層の物理特性の音響波を用いた非破壊測定方法 | |
| Niu et al. | Laboratory small-strain stiffness measurement using distributed acoustic sensing | |
| Filipiak et al. | Surface acoustic waves for the detection of small vibrations | |
| CN208477118U (zh) | 一种波速测试系统 | |
| CN113296149A (zh) | 一种基于微动探测技术的滑坡体地层分布情况探测方法 | |
| JP2952297B2 (ja) | 地盤の計測解析判定システム | |
| JP7237346B2 (ja) | 表面波探査解析装置及び表面波探査解析方法 | |
| JP4195171B2 (ja) | 地盤構造推定方法 | |
| JP2004138447A (ja) | 岩盤の物性評価方法 | |
| JP2006292415A (ja) | 空洞の探査方法 | |
| CN111551639B (zh) | 一种土壤弹性波色散曲线的测量方法 | |
| JP2000321360A (ja) | 地下速度構造推定方法及び地下速度構造推定システム | |
| JPH11248850A (ja) | 周波数−波数スペクトル解析方法 | |
| CN116972954B (zh) | 基于原位激振的岩土波速测量方法及装置 | |
| JP2000186319A (ja) | 地盤調査方法 | |
| Harrop et al. | Laser-Doppler vibrometer measurements of acoustic-to-seismic coupling in unconsolidated soils | |
| McNeese et al. | Investigation of piezoelectric bimorph bender transducers to generate and receive shear waves | |
| JP2011247737A (ja) | 音波を用いた探知方法、非接触音響探知システム、そのシステムで用いるプログラム、およびそのプログラムを記録した記録媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040618 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040804 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040930 |
|
| 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: 20041027 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20041029 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 3617036 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071119 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081119 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081119 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091119 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101119 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101119 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111119 Year of fee payment: 7 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121119 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121119 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131119 Year of fee payment: 9 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |