(57)【要約】【課題】Ｄ−Ｖｒｂ曲線における変曲点の抽出を自動化できるようにして地盤探査に伴う作業の簡略化を実現できるようにする。【解決手段】地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析装置であり、間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成した後、生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行い、選定した複数の変曲点に基づいて変曲点間の区間速度を算出する。特に、前記複数の変曲点の選定を、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの変化量ΔＶを因数として行う。【選択図】図１
2018077058.tif 000002
【特許請求の範囲】
【請求項１】
地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析装置であって、
該表面波探査解析装置は、間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成した後、生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行い、選定した複数の変曲点に基づいて変曲点間の区間速度を算出するものであり、
該表面波探査解析装置はまた、前記複数の変曲点の選定を、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの変化量ΔＶを因数として行うことを特徴とする表面波探査解析装置。
【請求項２】
前記表面波探査解析装置は、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分及び２階差分の少なくとも一方の計算によって前記変化量ΔＶを算出し、算出結果の符号が変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることを特徴とする請求項１に記載の表面波探査解析装置。
【請求項３】
前記表面波探査解析装置は、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分の計算によって前記変化量ΔＶを算出すると共に、深度データの１階差分の計算によって変化量ΔＤを算出した後、これらの算出結果から速度、深度についてデータの粗密度合いを算出し、算出した粗密度合いが速度、深度のそれぞれについてあらかじめ設定した値を上回るデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることを特徴とする請求項１に記載の表面波探査解析装置。
【請求項４】
地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析装置であって、
該表面波探査解析装置は、間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成した後、生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行い、選定した複数の変曲点に基づいて変曲点間の区間速度を算出するものであり、
該表面波探査解析装置はまた、前記複数の変曲点の選定を、前記入力信号ＳＡ、ＳＢの変化を因数として行うことを特徴とする表面波探査解析装置。
【請求項５】
前記表面波探査解析装置は、前記入力信号ＳＡ，ＳＢに対して単位周波数毎の加速度変化量ΔＬＡ、ΔＬＢを算出した後、算出した加速度変化量ΔＬＡとΔＬＢの差分を求めると共に、この差分の単位周波数あたりの変化量を求め、この変化量の符号が、“＋”から“−”に変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることを特徴とする請求項４に記載の表面波探査解析装置。
【請求項６】
地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析方法であって、
間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成するステップと、
生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行うステップと、
選定した複数の変曲点に基づいて変曲点間の区間速度を算出するステップと、を含み、
前記複数の変曲点の選定を行うステップは、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの変化量ΔＶを因数として行うことを特徴とする表面波探査解析方法。
【請求項７】
前記複数の変曲点の選定を行うステップにおいては、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分及び２階差分の少なくとも一方の計算によって前記変化量ΔＶを算出し、算出結果の符号が変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることを特徴とする請求項６に記載の表面波探査解析方法。
【請求項８】
前記複数の変曲点の選定を行うステップにおいては、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分の計算によって前記変化量ΔＶを算出すると共に、深度データの１階差分の計算によって変化量ΔＤを算出した後、これらの算出結果から速度、深度についてデータの粗密度合いを算出し、算出した粗密度合いが速度、深度のそれぞれについてあらかじめ設定した値を上回るデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることを特徴とする請求項６に記載の表面波探査解析方法。
【請求項９】
地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析方法であって、
間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成するステップと、
生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行うステップと、
選定した複数の変曲点に基づいて変曲点間の区間速度を算出するステップと、を含み、
前記複数の変曲点の選定を行うステップは、前記入力信号ＳＡ、ＳＢの変化を因数として行うことを特徴とする表面波探査解析方法。
【請求項１０】
前記複数の変曲点の選定を行うステップにおいては、前記入力信号ＳＡ，ＳＢに対して単位周波数毎の加速度変化量ΔＬＡ、ΔＬＢを算出した後、算出した加速度変化量ΔＬＡとΔＬＢの差分を求めると共に、この差分の単位周波数あたりの変化量を求め、この変化量の符号が、“＋”から“−”に変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることを特徴とする請求項９に記載の表面波探査解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地盤探査システムのための表面波探査解析方法及び表面波探査解析装置に関する。
【背景技術】
【０００２】
地盤探査の一手法として表面波探査法が知られている。表面波探査法は、地面の振動により地中を伝わる表面波（特にレイリー波）が、硬い地盤ほど速く伝わり、柔らかい地盤ほどゆっくりと伝わるという性質を持ち、振動する周波数が変わると伝わる深さが変化するという性質を利用している。表面波探査法では、起振機で地面に微弱な振動を与えて地中を伝わる表面波の速さを、起振現場に配置した少なくとも２つの検出器で検出し、これらの検出信号を用いて表面波の伝わる伝搬状況と速度を解析する。
【０００３】
このような表面波探査法を利用した地盤探査装置が特許文献１に記載されており、以下に簡単に説明する。
【０００４】
オペレータは、地盤探査を必要とする場所に起振機を設置すると共に、その近くの地盤上には間隔をおいて少なくとも２つの加速度検出器を設置する。起振機により地表面を起振周波数で上下に起振することにより、その周囲に表面波を発生させる。２つの加速度検出器からの検出信号は、Ａ（Analog）／Ｄ（Digital）変換等の信号処理機能を持つ計測部を通して加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）として出力される。計測部には、パーソナルコンピュータ等による解析装置が接続される。解析装置は、計測部からの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）を入力信号ＳＡ、ＳＢとして受け、入力信号ＳＡ、ＳＢに対してあらかじめ定められた解析プログラムに基づく信号処理を行い、伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出する。なお、伝搬平均速度は、通常、速度を表わすＶｒの上に平均を意味するバーを付して示されるが、ここでは、特に断りを入れない限り、表記の便宜上、バーを付さずにＶｒｂで示すこととする。
【０００５】
オペレータが、起振機による起振を、起振周波数を変化させながら繰り返し行うことにより、解析装置は深度Ｄ−伝搬平均速度Ｖｒｂ曲線（以下、Ｄ−Ｖｒｂ曲線と略称する）を生成し、それをモニターに表示する。Ｄ−Ｖｒｂ曲線の生成過程については後述する。
【０００６】
オペレータは、モニターに表示されたＤ−Ｖｒｂ曲線を確認し、Ｄ−Ｖｒｂ曲線上において地盤の物理的な性質が変化していると判定した箇所を変曲点（複数個）として選定する。変曲点の選定についてはマニュアル化されており、オペレータの選定作業はこのマニュアルに則って行われる。
【０００７】
すべての変曲点の選定後、選定された変曲点情報を基に変曲点間の区間速度が算出される。区間速度についても後述する。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００２−３４１０４８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
変曲点の選定についてはマニュアル化されているとは言え、オペレータの選定作業には熟練を必要とし、熟練したオペレータになるまでにはかなりの経験を必要とする。また、オペレータが変わると選定結果にもばらつきの生じることが避けられず、その結果、各種の算出結果にもばらつきが生じてしまう。
【００１０】
そこで、本発明の課題は、Ｄ−Ｖｒｂ曲線における変曲点の抽出を自動化できるようにして地盤探査に伴う作業の簡略化を実現できるようにすることにある。
【００１１】
本発明はまた、オペレータの違いによる計測結果のばらつきの解消と、計測精度の向上を実現できるようにすることにある。
【課題を解決するための手段】
【００１２】
本発明の第１の態様によれば、地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析装置であって、該表面波探査解析装置は、間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成した後、生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行い、選定した複数の変曲点に基づいて変曲点間の区間速度を算出するものであり、該表面波探査解析装置はまた、前記複数の変曲点の選定を、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの変化量ΔＶを因数として行うことを特徴とする表面波探査解析装置が提供される。
【００１３】
上記の第１の態様による表面波探査解析装置においては、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分及び２階差分の少なくとも一方の計算によって前記変化量ΔＶを算出し、算出結果の符号が変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることができる。
【００１４】
上記の第１の態様による表面波探査解析装置においてはまた、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分の計算によって前記変化量ΔＶを算出すると共に、深度データの１階差分の計算によって変化量ΔＤを算出した後、これらの算出結果から速度、深度についてデータの粗密度合いを算出し、算出した粗密度合いが速度、深度のそれぞれについてあらかじめ設定した値を上回るデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点としても良い。
【００１５】
本発明の第２の態様によれば、地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析装置であって、該表面波探査解析装置は、間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成した後、生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行い、選定した複数の変曲点に基づいて変曲点間の区間速度を算出するものであり、該表面波探査解析装置はまた、前記複数の変曲点の選定を、前記入力信号ＳＡ、ＳＢの変化を因数として行うことを特徴とする表面波探査解析装置が提供される。
【００１６】
上記の第２の態様による表面波探査解析装置においては、前記入力信号ＳＡ，ＳＢに対して単位周波数毎の加速度変化量ΔＬＡ、ΔＬＢを算出した後、算出した加速度変化量ΔＬＡとΔＬＢの差分を求めると共に、この差分の単位周波数あたりの変化量を求め、この変化量の符号が、“＋”から“−”に変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることができる。
【００１７】
本発明の第３の態様によれば、地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析方法であって、
間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成するステップと、
生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行うステップと、
選定した複数の変曲点に基づいて変曲点間の区間速度を算出するステップと、を含み、
前記複数の変曲点の選定を行うステップは、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの変化量ΔＶを因数として行うことを特徴とする表面波探査解析方法が提供される。
【００１８】
前記複数の変曲点の選定を行うステップにおいては、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分及び２階差分の少なくとも一方の計算によって前記変化量ΔＶを算出し、算出結果の符号が変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることができる。
【００１９】
前記複数の変曲点の選定を行うステップにおいてはまた、前記深度Ｄ−伝搬平均速度Ｖｒｂ曲線における速度データの１階差分の計算によって前記変化量ΔＶを算出すると共に、深度データの１階差分の計算によって変化量ΔＤを算出した後、これらの算出結果から速度、深度についてデータの粗密度合いを算出し、算出した粗密度合いが速度、深度のそれぞれについてあらかじめ設定した値を上回るデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点としても良い。
【００２０】
本発明の第４の態様によれば、地表面を上下に起振することにより、その周囲に発生する表面波を検出して、地盤探査を行う地盤探査システム用の表面波探査解析方法であって、
間隔Ｌをおいて起振現場に配置された少なくとも２つの加速度検出器の出力信号を、Ａ／Ｄ変換機能を持つ計測部を経由して入力信号ＳＡ、ＳＢとして受信し、受信した前記入力信号ＳＡ、ＳＢを処理して表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを算出し、前記伝搬平均速度Ｖｒｂ（ｆ）と前記深度Ｄ（ｆ）とを算出する処理を繰り返して深度Ｄ−伝搬平均速度Ｖｒｂ曲線を生成するステップと、
生成した深度Ｄ−伝搬平均速度Ｖｒｂ曲線に対して複数の変曲点の選定を行うステップと、
選定した複数の変曲点に基づいて変曲点間の区間速度を算出するステップと、を含み、
前記複数の変曲点の選定を行うステップは、前記入力信号ＳＡ、ＳＢの変化を因数として行うことを特徴とする表面波探査解析方法が提供される。
【００２１】
上記の第４の態様において前記複数の変曲点の選定を行うステップにおいては、前記入力信号ＳＡ，ＳＢに対して単位周波数毎の加速度変化量ΔＬＡ、ΔＬＢを算出した後、算出した加速度変化量ΔＬＡとΔＬＢの差分を求めると共に、この差分の単位周波数あたりの変化量を求め、この変化量の符号が、“＋”から“−”に変わるデータ位置を深度Ｄ、伝搬平均速度Ｖｒｂについて抽出して変曲点とすることができる。
【発明の効果】
【００２２】
本発明によれば、Ｄ−Ｖｒｂ曲線における変曲点の抽出を自動化できるようにしたことにより地盤探査に伴う作業の簡略化を実現することができる。
【００２３】
本発明によればまた、オペレータの違いによる計測結果のばらつきの解消と、計測精度の向上を実現することができる。
【図面の簡単な説明】
【００２４】
【図１】本発明が適用され得る、地盤探査装置の構成例を示した図である。
【図２】図１に示された起振機に与えられる電力の周波数バンドの一例（図ａ）と、図１に示された解析装置（パーソナルコンピュータ）で生成されるＤ−Ｖｒｂ曲線の一例（図ｂ）を示した図である。
【図３】本発明において表面波探査解析装置に格納される、図２（ｂ）に示されたＤ−Ｖｒｂ曲線のデータ構成の一例を示す。
【図４】本発明において表面波探査解析装置に格納される、ａｌｌデータのデータ構成の一例を示す。
【図５】本発明において表面波探査解析装置に格納される、ＦＬデータのデータ構成の一例を示す。
【図６】区間速度の算定例を説明するための図である。
【発明を実施するための形態】
【００２５】
本発明の実施形態について説明する前に、図１、図２を参照して、特許文献１に記載されている地盤探査装置について説明する。この地盤探査装置は、本発明にも適用可能である。
【００２６】
図１において、本地盤探査装置は、起振現場の地盤に設置される少なくとも２つの加速度検出器１１Ａ、１１Ｂと、これらの加速度検出器１１Ａ、１１Ｂからの加速度検出信号を受ける計測部１２と、を含む。計測部１２は、地震計部１２−１、Ａ／Ｄ変換部１２−２、通信部１２−３、発振部１２−４とから成る。地震計部１２−１は、ローパスフィルタ回路を内蔵し、アナログの加速度検出信号から加速度時系列信号を生成する。Ａ／Ｄ変換部１２−２は、地震計部１２−１からのアナログの加速度時系列信号をディジタルの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）に変換するためのもので、入力感度の自動調整機能を有する。通信部１２−３は、ディジタルの加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）を、計測部１２に接続される機器、例えばモニター付きのパーソナルコンピュータ（以下、ＰＣと略称する）１３に送信する。ＰＣ１３は加速度時系列信号Ａ（ｔ）、Ｂ（ｔ）を入力信号ＳＡ、ＳＢとして受け、内蔵メモリにあらかじめインストールされている解析プログラムソフトに基づいて入力信号ＳＡ、ＳＢの処理を行う機能を有し、地盤解析を行うための解析装置として機能する。ここで、（ｔ）は、時系列に沿って変化する任意の値を示す。
【００２７】
次に、動作について説明する。まず、探査場所に起振機１５、加速度検出器１１Ａ、１１Ｂを一直線上に設置する。加速度検出器１１Ａ、１１Ｂ間の距離をＬ（ｍ）とする。起振機１５を使って地表面を上下方向に起振することにより、起振機１５の周りに表面波を発生させる。地表面付近を伝搬する表面波（レイリー波）の上下振動を、加速度検出器１１Ａ、１１Ｂで検出する。加速度検出器１１Ａ、１１Ｂからの加速度検出信号は、地震計部１２−１のローパスフィルタ回路を通すことで、アナログの時系列信号となり、Ａ／Ｄ変換部１２−２に入力される。Ａ／Ｄ変換部１２−２でＡ／Ｄ変換された時系列信号Ａ（ｔ）、Ｂ（ｔ）は通信部１２−３からＰＣ１３へ転送される。
【００２８】
ＰＣ１３では、あらかじめ定められた解析プログラムに基づいて入力信号ＳＡ、ＳＢのパワースペクトルＧAA（ｆ）、ＧBB（ｆ）、クロススペクトルＧBA（ｆ）、伝達関数Ｈ（ｆ）、コヒーレンス関数γ2 （ｆ）等を計算する。これらのパワースペクトルＧAA（ｆ）、ＧBB（ｆ）、クロススペクトルＧBA（ｆ）、伝達関数Ｈ（ｆ）、コヒーレンス関数γ2 （ｆ）等は、内蔵のハードディスクに保存される。ＰＣ１３ではまた、伝達関数Ｈ（ｆ）より加速度検出器１１Ａ、１１Ｂ間の位相差Δθ（ｆ）を求め、続いてその時間差Δｔ（ｆ）を求める。ここで、（ｆ）は、周波数軸にそって変化する任意の値を示す。
【００２９】
ＰＣ１３では更に、時間差Δｔ（ｆ）と加速度検出器１１Ａ、１１Ｂ間の間隔Ｌより、表面波の伝搬平均速度Ｖｒｂ（ｆ）と深度Ｄ（ｆ）とを求める。
【００３０】
上記のパワースペクトルＧAA（ｆ）、ＧBB（ｆ）、クロススペクトルＧBA（ｆ）、伝達関数Ｈ（ｆ）、コヒーレンス関数γ2 （ｆ）等の計算過程については本発明の要旨ではなく、特許文献１に記載されているので省略するが、表面波の伝搬平均速度Ｖｒｂ（ｍ／ｓｅｃ）と深度Ｄ（ｍ）はそれぞれ、以下の式により算出する。 
Ｖｒｂ＝Ｌ／Δｔ（ｆ）＝２π×Ｆ×Ｌ／−Δθ（ｆ）
Ｄ＝λ／２＝Ｖｒｂ／２Ｆ＝π×Ｌ／−Δθ（ｆ）
但し、Ｆは励振（起振）信号の周波数、λは加速度検出信号の波長である。
【００３１】
上記の算出は所望のＤ−Ｖｒｂ曲線が得られるまで繰り返し行われる。すなわち、計測を１回行う毎に起振機１５に与える励振信号の周波数を変化させる。つまり、表面波の伝搬速度と逆数関係にある伝達関数Ｈ（ｆ）の位相差Δθ（ｆ）が、周波数別に計測される。次に、この位相差と周波数の関係から伝搬平均速度Ｖｒｂと深度Ｄとを計算し、繰り返し計測を行った結果としてモニター上にＤ−Ｖｒｂ曲線が生成、表示される。
【００３２】
図２（ａ）は、起振機１５に与えられる電力の周波数バンドを示し、周波数ｆ１〜ｆｎの領域を、ここではｆ１〜ｆｉ（バンドＢ１）、ｆｉ〜ｆｋ（バンドＢ２）、ｆｋ〜ｆｎ（バンドＢ３）の３つのバンドに分割している。この場合、図２（ｂ）に示すように、バンドＢ１での計測でＤ−Ｖｒｂ面には曲線Ｃ１が、バンドＢ２での計測で曲線Ｃ２が、バンドＢ３での計測で曲線Ｃ３がそれぞれ得られ、これらの曲線Ｃ１〜Ｃ３がモニター上に自動合成表示される。
【００３３】
図２（ｂ）には、Ｄ−Ｖｒｂ曲線に存在する変曲点のうち、以後の区間速度等の算出に必要とされる変曲点を○印で示している。以降の説明においては、特に断らない限り、Ｄ−Ｖｒｂ曲線に存在する変曲点のうち、以後の区間速度等の算出に必要とされる変曲点のみを変曲点と呼ぶこととする。
【００３４】
いずれにしても、オペレータは、モニター上の図２（ｂ）に示されるようなＤ−Ｖｒｂ曲線に対して変曲点の選定を行う。選定は、カーソルによる指示や、タッチペンによる指定で行われる。尚、変曲点は、Ｄ−Ｖｒｂ曲線が著しく変化している箇所が該当するが、計測する範囲において一番軟らかい地盤部分の速度を選択するため、Ｄ−Ｖｒｂ曲線の変化の左端、すなわち速度Ｖｒｂが小さい箇所を変曲点として選別する。
【００３５】
変曲点の選定が終了すると、ＰＣ１３は、選定された変曲点の情報を基に、区間速度等の算出を行う。
【００３６】
本発明による表面波探査解析装置は、ＰＣ１３と同様のパーソナルコンピュータで実現することができ、入力信号ＳＡ、ＳＢ（時系列信号Ａ（ｔ）、Ｂ（ｔ））の入力からＤ−Ｖｒｂ曲線の生成までは、ＰＣ１３に格納されている解析プログラムソフトと同じ機能を持つ解析プログラムソフトを搭載することで実現することができる。
【００３７】
次に、本発明の実施形態において扱われるデータの種類について説明する。
【００３８】
本発明による表面波探査解析装置は、図１で説明した少なくとも２つの加速度検出器１１Ａ、１１Ｂと、計測部１２と、電力増幅器１４と、起振機１５と共に組み合わせて構成される地盤探査システム用の表面波探査解析装置として適用され得る。それゆえ、本発明による表面波探査解析装置の実施形態においては、図１で説明したように、一定距離Ｌをおいて２箇所に設置された検出器からの検出信号を基に得られた伝搬平均速度Ｖｒｂ、深度Ｄ、励振信号の周波数Ｆの３種類のデータを、拡張子“＊＊．ＤＶ”のＤ−Ｖｒｂ曲線として内蔵のメモリに収納する。
【００３９】
図３はＤ−Ｖｒｂ曲線のデータ構成の一例を示す。Ｄ−Ｖｒｂ曲線のデータ構成は、データヘッダー部とデータ部とからなる。データヘッダー部にはソフト識別番号、計測番号、データ数、計測時速度軸最大値、等の変数が収納され、データ部には計算深度Ｄ（ｍ）、伝搬平均速度Ｖｒｂ（ｍ／ｓｅｃ）、周波数（Ｈｚ）の変数が収納される。
【００４０】
一方、計測に係るすべてのデータを、拡張子“＊＊．ａｌｌ”のａｌｌデータとして内蔵のメモリに収納している。
【００４１】
図４はａｌｌデータのデータ構成の一例を示す。ａｌｌデータのデータ構成もデータヘッダー部とデータ部とからなる。データヘッダー部にはソフト識別番号、計測番号、データ数、機器製造番号、等の変数が収納され、データ部には計算深度Ｄ（ｍ）、伝搬平均速度Ｖｒｂ（ｍ／ｓｅｃ）、周波数Ｆ（Ｈｚ）を含む１３種類の変数が収納される。
【００４２】
次に、本発明による表面波探査解析装置で実現される、表面波探査解析方法の特徴部分とも言える、Ｄ−Ｖｒｂ曲線における自動変曲点計算について説明する。
【００４３】
本発明における自動変曲点計算は、Ｄ−Ｖｒｂ曲線におけるデータの変化量ΔＶから変曲点を計算する第１の自動変曲点計算と、２箇所に設置された検出器からの検出信号に基づく時系列信号Ａ（ｔ）、Ｂ（ｔ）の変化を因数として変曲点を計算する第２の自動変曲点計算と、に大別される。第１、第２の自動変曲点計算のいずれも表面波探査解析装置によって実行することができる。
【００４４】
［第１の自動変曲点計算］
第１の自動変曲点計算は、前述したＤ−Ｖｒｂ曲線の生成後に実行される。第１の自動変曲点計算は、Ｄ−Ｖｒｂ曲線における速度データの変化ΔＶを因数とした自動変曲点計算であり、変化ΔＶの計算方法によって、さらに名称を“計算R”、“計算B”、“計算N”、“計算G”とする４種類の計算方法がある。
【００４５】
なお、Ｄ−Ｖｒｂ曲線は、前述のようにメモリに格納された拡張子“＊＊．ＤＶ”のＤＶデータに基づく。拡張子“＊＊．ＤＶ”のＤＶデータは、計算深度Ｄ、伝搬平均速度Ｖｒｂの２つの値が、計測時における励振信号の周波数Ｆの降順に格納されている。これらの計算深度、伝搬平均速度、周波数に関する計測データを、計測データＤＶ（Ｄ）ｉ、計測データＤＶ（Ｖ）ｉ、計測データＤＶ（Ｆ）ｉとし、総じて計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉとして表す。また、ｉは整数（ｉ＝１，２，・・・）であり、データの格納された順位を示す。
【００４６】
“計算R”は、変化ΔＶを１階差分中心）計算によって算定し、算定結果の符号が変わるデータ位置（伝搬平均速度、計算深度）を抽出し、変曲点とする。以下に詳しく説明する。
【００４７】
“計算R”算定内容
Ｄ−Ｖｒｂ曲線の計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉのうちの、計測データＤＶ（Ｖ）ｉの変化ΔＶに対して、１階差分を算定する。算定の際は、中心差分の計算を行い、条件式により、深度Ｄj、速度Ｖj、周波数Ｆjに関する変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊを判別する。なお、変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊは、深度Ｄｊ、速度Ｖｊ、周波数ｆｊによる変曲点データＦＬ（Ｄ）ｊ、ＦＬ（Ｖ）ｊ、ＦＬ（Ｆ）ｊでも表される。
【００４８】
つまり、後述する条件式を満たすiを用いて、ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉとなるデータを計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉの中から抽出し、これを変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊとする。また、ｊは整数（ｊ＝１，２，・・・）であり、周波数Ｆの降順となるデータ順位を示す。
【００４９】
１．１：速度データＤＶ（Ｖ）のｉ＋１番目とｉ−１番目の値の差分をΔＶｉとする。
ΔＶｉ＝ＤＶ（Ｖ）ｉ−１−ＤＶ（Ｖ）ｉ＋１
【００５０】
１．２：ＦＬ（Ｄ，Ｖ，Ｆ）ｊの有無を判別し、データ数を判別の上、
ＦＬ（Ｄ，Ｖ，Ｆ）１＝ＤＶ（Ｄ，Ｖ，Ｆ）１とする（初点の決定、および再設定）。
【００５１】
１．３：ΔＶｉとΔＶｉ＋２との条件により、１階差分の値の符号が、マイナスからプラスに転じる位置（ゼロクロス点）を抽出する。 
条件式１：（ΔＶｉ≧０∩ΔＶｉ＋１＜０）∪（ΔＶｉ＞０∩ΔＶｉ＋１≦０）
条件式１を満たす、ｉ（ｉ≧２）を選別する。 
条件式２：（ΔＶｉ≦０∩ΔＶｉ＋１＞０）∪（ΔＶｉ＜０∩ΔＶｉ＋１≧０）
条件式２を満たす、ｉ（ｉ≧２）を選別する。 
条件式３：条件式１あるいは条件式１を満たす、ｉ（ｉ≧２）を選別する。 
条件式１〜条件式３の選定は、基本は条件式２を用いるものとし、同一の現場においては、同一の条件式を適用するものとする。
【００５２】
１．４：選別したｉを代入して、ＤＶｉとする。この値に対して、以下の３つの条件を満たすか否かを確認する。(自動変曲点のデータ数、間隔の調整)
条件式４：自動変曲点の判断を行う深度Ｄｓ（任意で与える）に対して、
ＤＶ（Ｄ）ｉ≦ｋ×Ｄｓを満たす。但し、ｋは任意の値で通常は“１．１”とする。 
任意で与える削除(深度)レベルＤ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｄ）ｊ−１−ＤＶ（Ｄ）ｉ）≧Ｄ＿ＤＥを満たす。但し、Ｄ＿ＤＥは任意の値で通常は“０．１”とする。 
任意で与える削除(速度)レベルＶ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｖ）ｊ−１−ＤＶ（Ｖ）ｉ）≧Ｖ＿ＤＥを満たす。但し、Ｖ＿ＤＥは任意の値で通常は“４．０”とする。
【００５３】
１．５：上記１．４の条件を満たすＤＶｉにより、ＦＬｊを決定する。 
ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉ
【００５４】
以上のように、“計算R”によれば、これまでオペレータがＤ−Ｖｒｂ曲線に対して行っていた変曲点の選定と同様に、図２（ｂ）に○で示されるような、Ｄ−Ｖｒｂ曲線が著しく変化しかつ左端の変化箇所を、変曲点データＦＬ（Ｄ）ｊ、ＦＬ（Ｖ）ｊ、ＦＬ（Ｆ）ｊとして選別することができる。
【００５５】
“計算B”は、速度の変化ΔＶを２階差分（前進）計算によって算定し、算定結果の符号が変わるデータ位置（伝搬平均速度、計算深度）を抽出し、変曲点とする。以下に詳しく説明する。
【００５６】
“計算B”算定内容
Ｄ−Ｖｒｂ曲線のデータＤＶ（Ｄ，Ｖ，Ｆ）ｉのうちの、速度Ｖの変化に対して、２階差分を算定する。算定の際は、前進差分の計算を行い、条件式により変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊを算定する。
【００５７】
つまり、後述する条件式を満たすiを用いて、ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉとなるデータを計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉの中から抽出し、これを変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊとする。
【００５８】
２．１：速度データＤＶ（Ｖ）のｉ−２番目、ｉ−１番目およびｉ番目の速度の値による、２階差分をΔ２Ｖｉとする。 
Δ２Ｖｉ＝ＤＶ（Ｖ）ｉ−２−２ＤＶ（Ｖ）ｉ−１＋ＤＶ（Ｖ）ｉ
【００５９】
２．２：ＦＬ（Ｄ，Ｖ，Ｆ）ｊの有無を確認し、データ数を確認の上、
ＦＬ（Ｄ，Ｖ，Ｆ）１＝ＤＶ（Ｄ，Ｖ，Ｆ）１
とする（初点の決定、および再設定）。
【００６０】
２．３：Δ２ＶｉとΔ２Ｖｉ＋２との条件により、１階差分の値の符号がマイナスからプラスに転じる位置（ゼロクロス点）を抽出する。 
条件式１：（Δ２Ｖｉ≧０∩Δ２Ｖｉ＋１＜０）∪（Δ２Ｖｉ＞０∩Δ２Ｖｉ＋１≦０）
条件式を満たす、ｉ（ｉ≧２）を選別する。 
条件式２：（Δ２Ｖｉ≦０∩Δ２Ｖｉ＋１＞０）∪（Δ２Ｖｉ＜０∩Δ２Ｖｉ＋１≧０）
条件式を満たす、ｉ（ｉ≧２）を選別する。 
条件式３：条件式１あるいは条件式２を満たす、ｉ（ｉ≧２）を選別する。 
条件式１〜条件式３の選定は、基本は条件式２を用いるものとし、同一の現場においては、同一の条件式を適用するものとする。
【００６１】
２．４：選別したｉを代入して、ＤＶｉとする。この値に対して、以下の３つの条件を満たすか否かを確認する。(自動変曲点のデータ数、間隔の調整)
条件式４：自動変曲点の判断を行う深度Ｄｓ（任意で与える）に対して、
ＤＶ（Ｄ）ｉ≦ｋ×Ｄｓを満たす。但し、ｋは任意の値で通常は“１．１”とする。 
任意で与える削除(深度)レベルＤ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｄ）ｊ−１−ＤＶ（Ｄ）ｉ）≧Ｄ＿ＤＥを満たす。但し、Ｄ＿ＤＥは任意の値で通常は“０．１”とする。 
任意で与える削除(速度)レベルＶ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｖ）ｊ−１−ＤＶ（Ｖ）ｉ）≧Ｖ＿ＤＥを満たす。但し、Ｖ＿ＤＥは任意の値で通常は“４．０”とする。
【００６２】
２．５：上記２．４の条件を満たす、ＤＶｉにより、ＦＬｊを決定する。 
ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉ
【００６３】
以上のように、“計算B”によれば、これまでオペレータがＤ−Ｖｒｂ曲線に対して行っていた変曲点の選定と同様に、図２（ｂ）に○で示されるような、Ｄ−Ｖｒｂ曲線が著しく変化しかつ左端の変化箇所を、変曲点データＦＬ（Ｄ）ｊ、ＦＬ（Ｖ）ｊ、ＦＬ（Ｆ）ｊとして選別することができる。
【００６４】
“計算N”は、“R”/ “B”計算によって変化ΔＶを算定し、算定結果の符号が変わるデータ位置（伝搬平均速度、計算深度）を抽出し、変曲点とする。以下に詳しく説明する。
【００６５】
“計算N”算定内容
Ｄ−Ｖｒｂ曲線のデータＤＶ（Ｄ，Ｖ，Ｆ）ｉのうちの、速度Ｖの変化に対して、１階差分、２階差分を算定する。算定の際には１階差分は中立、２階差分は前進計算を行い、両計算の結果並びに、条件式により変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊを算定する。
【００６６】
つまり、後述する条件式を満たすiを用いて、ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉとなるデータを計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉの中から抽出し、これを変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊとする。
【００６７】
３．１：速度データＤＶ（Ｖ）の（ｉ−１）番目と（ｉ＋１）番目の速度Ｖの値の差分をΔＶｉとする。 
ΔＶｉ＝ＤＶ（Ｖ）ｉ−１−ＤＶ（Ｖ）ｉ＋１
【００６８】
３．２：変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊの有無を確認し、データ数を確認の上、
ＦＬ（Ｄ，Ｖ，Ｆ）１＝ＤＶ（Ｄ，Ｖ，Ｆ）１
とする（初点の決定、および再設定）。
【００６９】
３．３：差分ΔＶｉとΔＶｉ＋２との条件により、１階差分の値の符号が、マイナスからプラスに転じる位置（ゼロクロス点）を抽出し、その位置での深度、伝搬平均速度を抽出する。 
条件式１：（ΔＶｉ≧０∩ΔＶｉ＋１＜０）∪（ΔＶｉ＞０∩ΔＶｉ＋１≦０）
条件式を満たす、ｉ（ｉ≧２）を選別する。 
条件式２：（ΔＶｉ≦０∩ΔＶｉ＋１＞０）∪（ΔＶｉ＜０∩ΔＶｉ＋１≧０）
条件式を満たす、ｉ（ｉ≧２）を選別する。 
条件式３：条件式１あるいは条件式２を満たす、ｉ（ｉ≧２）を選別する。 
条件式１〜条件式３の選定は、基本は条件式１を用いるものとし、同一の現場においては、同一の条件式を適用するものとする。
【００７０】
３．４：上記３．３の条件を満たす、ＤＶｉにより、ＦＬｊを仮決定する。 
ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉ
【００７１】
３．５：上記３．４で得られた、変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊのうち深度Ｄと、任意で与える設定最大層厚ＤＴ＿Ｌとに対して、
ＦＬ（Ｄ）j＋１−ＦＬ（Ｄ）j＜ＤＴ＿Ｌを満たすjを抽出する。但し、ＤＴ＿Ｌは任意の値で通常は“３．０”とする。
【００７２】
３．６：ＦＬ（Ｄ，Ｖ，Ｆ）ｊのうち、周波数Ｆを用いて、ＦＬ（Ｆ）j＝ＤＶ（Ｆ）ｉならびにＦＬ（Ｆ）j+1＝ＤＶ（Ｆ）ｉとなるiを抽出し、各々をSi、Eiとする。
【００７３】
３．７：速度データＤＶ（Ｖ）のｉ−２番目、ｉ−１番目およびｉ番目の速度の値による、２階差分をΔ２Ｖiとする。 
Δ２Ｖi＝ＤＶ（Ｖ）i−２−２ＤＶ（Ｖ）i―１＋ＤＶ（Ｖ）i
但し、iは、SiからEiまでの整数とする。
【００７４】
３．８：Δ２ＶｉとΔ２Ｖｉ＋２との条件により、２階差分の値の符号が、マイナスからプラスに転じる位置（ゼロクロス点）を抽出し、その位置での深度、伝搬平均速度を抽出する。 
ここまでの計算N算定について要約すると、最初に１階差分による計算で、仮の変曲点を抽出し、次に、仮の変曲点の各点の深度の差が上記３．５で決められた３．０（ｍ）以内であるかを判定する。もし３．０ｍ以上空いた場合は、この３．０ｍ以上空いた範囲を抽出し、この範囲に対して、改めて２階差分の計算を実施する。 
条件式１：（Δ２Ｖｉ≧０∩Δ２Ｖｉ＋１＜０）∪（Δ２Ｖｉ＞０∩Δ２Ｖｉ＋１≦０）
条件式を満たす、ｉ（ｉ≧２）を選別する。 
条件式２：（Δ２Ｖｉ≦０∩Δ２Ｖｉ＋１＞０）∪（Δ２Ｖｉ＜０∩Δ２Ｖｉ＋１≧０）
条件式を満たす、ｉ（ｉ≧２）を選別する。 
条件式３：条件式１あるいは条件式２を満たす、ｉ（ｉ≧２）を選別する。 
条件式１〜条件式３の選定は、基本は条件式２を用いるものとし、同一の現場においては、同一の条件式を適用するものとする。
【００７５】
３．９：選別したｉを代入して、ＤＶｉとする。この値に対して、以下の３つの条件を満たすか否かを確認する。(自動変曲点のデータ数、間隔の調整)
条件式４：自動変曲点の判断を行う深度Ｄｓ（任意で与える）に対して、
ＤＶ（Ｄ）ｉ≦ｋ×Ｄｓを満たす。但し、ｋは任意の値で通常は“１．１”とする。 
任意で与える削除(深度)レベルＤ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｄ）ｊ−１−ＤＶ（Ｄ）ｉ）≧Ｄ＿ＤＥを満たす。但し、Ｄ＿ＤＥは任意の値で通常は“０．１”とする。 
任意で与える削除(速度)レベルＶ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｖ）ｊ−１−ＤＶ（Ｖ）ｉ）≧Ｖ＿ＤＥを満たす。但し、Ｖ＿ＤＥは任意の値で通常は“４．０”とする。
【００７６】
３．１０：上記３．９の条件を満たす、ＤＶｉにより、ＦＬｊｊを決定する。 
ＦＬ（Ｄ）ｊｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊｊ＝ＤＶ（Ｆ）ｉ
【００７７】
３．１１：上記３．４で得られたＦＬｊにＦＬｊｊを追加し、変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊとする。
【００７８】
以上のように、“計算N”によれば、これまでオペレータがＤ−Ｖｒｂ曲線に対して行っていた変曲点の選定と同様に、図２（ｂ）に○で示されるような、Ｄ−Ｖｒｂ曲線が著しく変化しかつ左端の変化箇所を、変曲点データＦＬ（Ｄ）ｊ、ＦＬ（Ｖ）ｊ、ＦＬ（Ｆ）ｊとして選別することができる。
【００７９】
“計算G”は、Ｄ−Ｖｒｂ曲線におけるデータ変化ΔＶ（伝搬平均速度値の変化）に加えてデータ変化ΔＤ（深度方向の変化）を同時に確認する。これらの値から、データの粗密の度合いＴＧを算定し、データの粗密の度合いＴＧの値があらかじめ設定した値ＧＲｌｅｖｅｌを上回るデータ位置（伝搬平均速度、計算深度）を抽出し、変曲点とする。
【００８０】
“計算G”算定内容
Ｄ−Ｖｒｂ曲線のデータＤＶ（Ｄ，Ｖ，Ｆ）ｉのうちの、速度Ｖの変化、並びに深度Ｄの変化を整理し、データの粗密の度合いの程度により変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊを算定する。
【００８１】
つまり、後述する条件式を満たすiを用いて、ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉとなるデータを計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉの中から抽出し、これを変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊとする。
【００８２】
４．１：Ｄ−Ｖｒｂ曲線のデータＤＶ（Ｄ，Ｖ，Ｆ）ｉの、速度Ｖ、深度Ｄのｉ−１番目およびｉ番目の値による変化量ΔＤ＿Ｖｉを以下の数式（１）により算出する。
【数１】
2018077058.tif 000003
但し、ΔＤ＿Ｖｉは速度、深度の変化量、Ｄｉはｉ番目の深度データ、Ｖｉはｉ番目の速度データ、ｉは整数（ｉ＝１，２，・・・）、ＶＴは計測時に設定した最大速度値であり、自動計算に際しては、この最大速度値より小さい速度データについて処理を行う。同様に、ＤＴは計測時に設定した最大深度値であり、自動計算に際しては、この最大深度値より小さい深度データについて処理を行う。
【００８３】
４．２：データの平均化に関する係数Ｔ＿Ｄｉを以下の数式（２）により設定する。
【数２】
2018077058.tif 000004
但し、Ｔ＿Ｄｉは平均化に関する係数で、Ｔ＿Ｄ１＝０、ＥＭは計算のための設定値（１〜９の整数で、通常４とする）である。
【００８４】
４．３：データの粗密の度合いを示す変数ＴＧｉを以下の数式（３）により設定する。
【数３】
2018077058.tif 000005
【００８５】
４．４：データの粗密度合いの程度を示す値として変数ＧＲＬを設定し、変数ＧＲＬと変数ＴＧｉとの条件を設け、これを満たすｉ（ｉ≧２）を選別する。但し、変数ＧＲＬはデータの粗密度合いの程度を示す値で、０．１〜１．５の値をとり、通常は０．５とする。 
条件式１：｜ＴＧｉ｜≧ＧＲＬ
条件式を満たす、ｉ（ｉ≧２）を選別する。
【００８６】
４．５：選別したｉを代入して、ＤＶｉとする。この値に対して、以下の３つの条件を満たすか否かを確認する。(自動変曲点のデータ数、間隔の調整)
条件式２：自動変曲点の判断を行う深度Ｄｓ（任意で与える）に対して、
ＤＶ（Ｄ）ｉ≦ｋ×Ｄｓを満たす。但し、ｋは任意の値で通常は“１．１”とする。 
任意で与える削除(深度)レベルＤ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｄ）ｊ−１−ＤＶ（Ｄ）ｉ）≧Ｄ＿ＤＥを満たす。 
任意で与える削除(速度)レベルＶ＿ＤＥに対して
ａｂｓ（ＦＬ（Ｖ）ｊ−１−ＤＶ（Ｖ）ｉ）≧Ｖ＿ＤＥを満たす。
【００８７】
４．６：上記４．５の条件を満たすＤＶｉにより、ＦＬｊを決定する。 
ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉ
【００８８】
以上のように、“計算G”によれば、これまでオペレータがＤ−Ｖｒｂ曲線に対して行っていた変曲点の選定と同様に、図２（ｂ）に○で示されるような、Ｄ−Ｖｒｂ曲線が著しく変化しかつ左端の変化箇所を、変曲点データＦＬ（Ｄ）ｊ、ＦＬ（Ｖ）ｊ、ＦＬ（Ｆ）ｊとして選別することができる。
【００８９】
［第２の自動変曲点計算］
第２の自動変曲点計算は、前述した計測部１２からの時系列信号Ａ（ｔ）、Ｂ（ｔ）を入力信号ＳＡ、ＳＢとして受けて実行される。すなわち、第２の自動変曲点計算は、２か所に設置された加速度検出器１１Ａ、１１Ｂからの検出信号を計測部１２経由で入力した入力信号ＳＡ、ＳＢの変化を因数とした自動変曲点計算である。加速度検出器１１Ａ、１１Ｂからの入力信号は、図４で説明したように、計測を実施した周波数に対する値として、“allデータ”に収録されている。加速度検出器１１Ａ、１１ＢからのＡｃｈ、Ｂｃｈ各々の入力信号ＳＡ、ＳＢに対して、単位周波数毎の加速度変化量ΔＬＡ、ΔＬＢを求め、さらに２つの加速度検出器１１Ａ、１１Ｂ間におけるΔＬＡならびにΔＬＢの差分を求める。この差分の単位周波数あたりの変化量を求め、この変化量の符号が、“＋”から“−”に変化するデータ位置を変曲点とする。以下に、第２の自動変曲点計算の算定内容について詳細に説明する。
【００９０】
第２の自動変曲点計算の算定内容
５．１：２つの加速度検出器１１Ａ、１１ＢによるＡｃｈ、Ｂｃｈ各々の入力信号ＳＡ、ＳＢに対して、単位周波数毎の加速度変化量ΔＬＡ、ΔＬＢを求め、さらに２つの加速度検出器１１Ａ、１１Ｂ間におけるΔＬＡならびにΔＬＢの差分を求める。この差分の単位周波数あたりの変化量を求め、変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊを算定する。
【００９１】
つまり、加速度変化量ΔＬＡならびにΔＬＢに基づく条件式を設け、この条件式を満たすＤＶ（Ｆ）ｉを用いて、ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ、ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ、ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉとなるデータを計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉの中から抽出し、これを変曲点データＦＬ（Ｄ，Ｖ，Ｆ）ｊとする。
【００９２】
５．２：加速度検出器１１Ａ、１１Ｂによる計測時に、Ａｃｈ、Ｂｃｈ各々で得られる周波数毎の最大加速度値を以下の数式（４）、（５）により変換し、デシベル値ＬＡＶａ、ＬＢＶａとする。
【数４】
2018077058.tif 000006
【数５】
2018077058.tif 000007
但し、ＬＡＶａはＡｃｈの振動加速度（ｄｂ）、ＳＡはＡｃｈの加速度（ｇａｌ）、ＬＢＶａはＢｃｈの振動加速度（ｄｂ）、ＳＢはＢｃｈの加速度（ｇａｌ）であり、いずれも周波数毎の振動加速度の最大値である。
【００９３】
ＳＡおよびＳＢは、計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉと同様に計測時にメモリに記録された拡張子＊＊．ａｌｌ”のａｌｌデータに格納されている。拡張子“＊＊．ＤＶ”のＤ−Ｖｒｂ曲線並びに拡張子“＊＊．ａｌｌ”のａｌｌデータは、計測に係る値をすべて同じ周波数変化に対して格納している。このことから、ＳＡおよびＳＢ、あるいはＬＡＶａおよびＬＢＶａは、計測データＤＶ（Ｄ，Ｖ，Ｆ）ｉと同様、整数ｉを用いて順位を表している（ｉ＝１，２，・・・）。
【００９４】
最後に、以下の数式（６）、（６‘）、（７）、（７’）により、ある点とその前後１点の３点の振動加速度ＬＡＶａ（i），ＬＡＶａ（i−１），ＬＡＶａ（i＋１）及び振動加速度ＬＢＶａ（i），ＬＢＶａ（i−１），ＬＢＶａ（i＋１）の平均化を行う。
【数６】
2018077058.tif 000008
ＬＡＶａ（i）＝ＬＡＶａ（i）‘                          （６’）
【数７】
2018077058.tif 000009
ＬＢＶａ（i）＝ＬＢＶａ（i）‘                          （７’）
【００９５】
５．３：振動加速度データのすべての打点に対し、以下の数式（８）、（９）により単位周波数毎の加速度変化量ΔＬＡ、ΔＬＢを求める。
【数８】
2018077058.tif 000010
【数９】
2018077058.tif 000011
但し、ΔＬＡ、ΔＬＢはＡｃｈ、Ｂｃｈの振動加速度データの１Ｈｚ当たりの変化量、ＤＶ（Ｆ）iは計測データＤＶのi番目の周波数、ＬＡＶａ（i）、ＬＢＶａ（i）は拡張子“＊＊．ａｌｌ”に格納されているＡｃｈ、Ｂｃｈのi番目の振動加速度（ｄｂ）を示す。
【００９６】
５．４：上記５．３で得られた、Ａｃｈ、Ｂｃｈ各々の振動加速度データの変化量の差ΔΔＬを求める。 
ΔΔＬ＝ΔＬＡ−ΔＬＢ
但し、ΔΔＬは１Ｈｚ当たりのＡｃｈの振動加速度データの変化量とＢｃｈの振動加速度データの変化量の差である。
【００９７】
最後に、以下の数式（１０）〜（１４）により、振動加速度データの変化量について、ある点とその前後１〜５点の平均化を１〜５回行う。平均化の内容については、境界面の自動判別の内容を確認しながら調整を行うことが望ましい。
【数１０】
2018077058.tif 000012
【数１１】
2018077058.tif 000013
【数１２】
2018077058.tif 000014
【数１３】
2018077058.tif 000015
【数１４】
2018077058.tif 000016
上記の数式（１０）〜数式（１４）に共通して、以下の数式（１５）が適用される。 
ΔΔＬ（i）＝ΔΔＬ’（i）    （１５）
【００９８】
数式（１０）、（１５）はある点とその前後１点の平均化の際の計算内容を示す。同様に、数式（１１）、（１５）はある点とその前後２点の平均化の際の計算内容、数式（１２）、（１５）はある点とその前後３点の平均化の際の計算内容、数式（１３）、（１５）はある点とその前後４点の平均化の際の計算内容、数式（１４）、（１５）はある点とその前後５点の平均化の際の計算内容をそれぞれ示す。数式（１０）〜（１４）は繰り返して行うことが可能であるが、５回を上限とする。
【００９９】
５．５：上記５．４で得られた振動加速度データの変化量の差ΔΔＬに対して以下の数式（１６）により単位周波数毎の変化量ΔΔΔＬを求める。
【数１５】
2018077058.tif 000017
但し、ＤＶ（Ｆ）iは計測データＤＶのi番目の周波数、ΔΔΔＬはΔΔＬデータの１Ｈｚ当たりの変化量、ΔΔＬ（i）はΔΔＬデータのi番目の値である。
【０１００】
５．６：上記５．５で得られたΔΔΔＬデータに対して境界面の判別を行う。 
ｋ＝ΔΔΔＬ（i）並びにｋｋ＝ΔΔΔＬ（i＋１）
条件（ｋ×ｋｋ＜０、且つｋｋ＜ｋ）（但し、ｋはΔΔΔＬデータのi番目、ｋｋはΔΔΔＬデータの（i＋１）番目の値を示す）を満たすiに対してＦＬ（Ｄ，Ｖ，ｆ）jを決定する。 
ＦＬ（Ｆ）ｊ＝ＤＶ（Ｆ）ｉ
ＦＬ（Ｄ）ｊ＝ＤＶ（Ｄ）ｉ
ＦＬ（Ｖ）ｊ＝ＤＶ（Ｖ）ｉ
【０１０１】
以上のように、“第２の自動変曲点計算”によれば、これまでオペレータがＤ−Ｖｒｂ曲線に対して行っていた、変曲点の選定と同様に、図２（ｂ）に○で示されるような、Ｄ−Ｖｒｂ曲線が著しく変化しかつ左端の変化箇所を、変曲点データＦＬ（Ｄ）ｊ、ＦＬ（Ｖ）ｊ、ＦＬ（Ｆ）ｊとして選別することができる。
【０１０２】
以上、第１、第２の自動変曲点計算について説明したが、上述した第１、第２の自動変曲点計算のいずれにおいても、図３で説明したように、変曲点の選定結果は、拡張子を“＊＊．ＦＬ”とし、ＦＬデータとして保存する。
【０１０３】
図５は、表面波探査解析装置に格納される、ＦＬデータのデータ構成の一例を示す。ＦＬデータのデータ構成もデータヘッダー部とデータ部とからなる。データヘッダー部にはソフト識別番号、計測番号、データ数等の変数が収納され、データ部には変曲点における計算深度Ｄ（ｍ）、伝搬平均速度Ｖ（ｍ／ｓ）、周波数Ｆ（Ｈｚ）の３つの値を収録する。
【０１０４】
ＦＬデータに収録されている伝搬平均速度Ｖｒｂは、地表面から計算深度Ｄまで、全体の地盤が振動した速度値を指す。ここで、変曲点によって区分した深度方向の各々範囲における速度値を区間速度として算定することで、深度方向における速度構造を算定することができる。
【０１０５】
図６は、区間速度の算定例を説明するための図である。図６（ａ）に示すように、Ｄ−Ｖｒｂ曲線から変曲点が抽出されると、地表面から第１、第２、第３、第４の４層の地盤が想定される。
【０１０６】
区間速度の算定については、例えば図２（ｂ）に示す第２層について言えば、第２層の上位の変曲点における伝搬平均速度Ｖｊ、計算深度Ｄｊ、並びに第２層の下位の変曲点における伝搬平均速度Ｖｊ＋１、計算深度Ｄｊ＋１により算定する。
【０１０７】
本実施形態では、区間速度を求める算定式を、区間速度計算式（１−１）、（１−２）、（１−３）として４種類用意している。これらの算定式の組み合わせ、条件式を用いて、探査場所における地盤の速度構造を求める。以下に、第２の区間速度の算定内容について詳細に説明する。
【０１０８】
区間速度算定内容
下記の区間速度計算式（１−１）〜（１−３）を用いて、得られた変曲点データから区間速度Ｖｒｊを計算する。ここでは、伝搬平均速度と区間速度を区別するために、これまでの説明において使用した表記を変更して、区間速度をＶｒｊと表記し、伝搬平均速度をＶの上にバーを付けて表記し、Ｖバーと呼ぶこととする。
【０１０９】
６．１：区間速度計算式（１−１）
【数１６】
2018077058.tif 000018
但し、ＶｊバーはＦＬデータのｊ番目の伝搬平均速度、ＤｊはＦＬデータのｊ番目の深度、jは整数（ｊ＝１，２，・・・）である。
【０１１０】
６．２：区間速度計算式（１−２）
【数１７】
2018077058.tif 000019
【０１１１】
６．３：区間速度計算式（１−３）
【数１８】
2018077058.tif 000020
【０１１２】
６．４：第１の区間速度計算方法（区間速度計算式（１−１）と下記の複数の条件式を用いる）
【数１９】
2018077058.tif 000021
【数２０】
2018077058.tif 000022
【数２１】
2018077058.tif 000023
【数２２】
2018077058.tif 000024
但し、ρＶは伝搬平均速度Ｖｒｂに対する補正値であり、後述する数式（１２−８）において、補正後の値CVが算出される。
【０１１３】
【数２３】
2018077058.tif 000025
【数２４】
2018077058.tif 000026
【数２５】
2018077058.tif 000027
【数２６】
2018077058.tif 000028
【数２７】
2018077058.tif 000029
【数２８】
2018077058.tif 000030
【数２９】
2018077058.tif 000031
【数３０】
2018077058.tif 000032
【数３１】
2018077058.tif 000033
但し、下記の条件式（１２−１１）の場合、区間速度計算式（１−２）を用いる。
【０１１４】
【数３２】
2018077058.tif 000034
【数３３】
2018077058.tif 000035
【数３４】
2018077058.tif 000036
【数３５】
2018077058.tif 000037
但し、下記の条件式（１２−１３）の場合、区間速度計算式（１−２）を用いる。
【０１１５】
【数３６】
2018077058.tif 000038
【数３７】
2018077058.tif 000039
【０１１６】
６．５：第２の区間速度計算方法（区間速度計算式（１−３）、（１−２）と下記の複数の条件式を用いる）
【数３８】
2018077058.tif 000040
【数３９】
2018077058.tif 000041
【数４０】
2018077058.tif 000042
【数４１】
2018077058.tif 000043
【０１１７】
６．６：第３の区間速度計算方法（区間速度計算式（１−１）、（１−２）と下記の複数の条件式を用いる）
【数４２】
2018077058.tif 000044
【数４３】
2018077058.tif 000045
【数４４】
2018077058.tif 000046
【数４５】
2018077058.tif 000047
【０１１８】
６．７：第４の区間速度計算方法（区間速度計算式（１−３）を用いる）
【数４６】
2018077058.tif 000048
以上のようにして、地表面から第１層、第２層、・・・、第j層の各層の地盤の区間速度が算定される。
【０１１９】
以上説明してきたように、本発明の実施形態によれば、Ｄ−Ｖｒｂ曲線における複数の変曲点を自動的に抽出することができるので、地盤探査に伴う作業の簡略化を実現することができるうえに、オペレータの違いによる計測結果のばらつきの解消と、計測精度の向上を実現することができる。そして、抽出した複数の変曲点によって区分される深度方向の複数の地盤層の各層における区間速度を算出することにより、探査地盤の深度方向における速度構造を算定することができる。
【産業上の利用可能性】
【０１２０】
本発明による表面波探査解析装置及び表面波探査解析方法は、地盤探査への適用に適している。
【符号の説明】
【０１２１】
１１Ａ、１１Ｂ    加速度検出器
１２    計測部
１３    パーソナルコンピュータ
【図１】
2018077058.tif 000049
【図２】
2018077058.tif 000050
【図３】
2018077058.tif 000051
【図４】
2018077058.tif 000052
【図５】
2018077058.tif 000053
【図６】
2018077058.tif 000054