JP2001033548A - リアルタイム探査方法及び装置 - Google Patents

Abstract

(57)【要約】 【課題】 媒質中の物体の存在をリアルタイムで検出
でき、同時に、物体までの深さ検出も短時間で行えるリ
アルタイム探査方法及び装置を提供する。 【解決手段】 送受信部１６から移動方向前方への距離
Ｘ_i と同方向の垂直方向への距離Ｙ_j と送受信部１６が
放射する波動信号の媒質１中の伝搬速度Ｖ_k を座標とす
る３次元配列を予め作成し、送受信部１６を移動させな
がら、所定の時間間隔または移動距離毎に反射信号の反
射時間と信号強度を対応付けて保存し、３次元配列の各
ボクセル値を−Ｘ方向へ１ボクセル分シフトするととも
に、シフト後の各ボクセル値に対して、座標値Ｘ_i ，Ｙ
_j ，Ｖ_k で決定される反射時間に対応する保存された信
号強度を単純加算または加重加算し、当該加算処理の度
に、３次元配列の各（Ｘ_i ，Ｙ_j ）座標におけるＶ軸方
向の最大ボクセル値の何れかまたはその中の最大値が所
定閾値を超えたとき、物体２が存在すると判定する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、媒質の表面または
媒質中で送受信部を移動させながら、電磁波または音波
による波動信号を媒質中へ放射し、その媒質中に存在す
る物体からの反射信号を受信することにより、媒質中に
存在する物体を探査するリアルタイム探査方法及び装置
に関する。

【０００２】

【従来の技術】一般的なリアルタイム探査方法及び装置
としては、電磁波の反射を用いて地中にある埋設物また
は空洞を探査するリアルタイム探査方法及び装置があ
る。例えば、電子情報通信学会技術研究報告Ｖｏｌ．９
５，Ｎｏ．１９１、ｐｐ９〜１４「地中断面観測パター
ンに含まれる双曲線形状のリアルタイム認識方法」（１
９９５年７月）に開示されている方法（以下、従来方法
という。）が知られている。この従来方法の場合、図１
２に示すフローチャートに基づいて、アンテナが１ｃｍ
移動する毎にステップ３から８までの処理が繰り返し実
行される。ステップ４では反射信号を受信し、ステップ
５では反射信号の極大値と極小値を抽出し、各々を
「１」と「−１」に、それ以外を「０」とする。ステッ
プ６のセグメント化処理では、同一の値となった座標の
一定範囲内の近傍（例えば、反射時間で±２ｎｓ）にあ
る座標を同一のセグメントとして分類する。そして、ス
テップ７の加算処理では、ステップ６のセグメント化処
理で分類された、同一のセグメントに属する座標の組に
対して、反射時間の差が加算され、その加算値が所定の
値を超えた時に、ステップ８でブザーを鳴動させる。

【０００３】

【発明が解決しようとする課題】しかしながら、上記従
来方法では、ステップ６のセグメント化処理において、
同一の値となった座標の一定範囲内の近傍にある座標を
同一のセグメントとして分類したうえ、ステップ７の加
算処理において、それらを単純に加算しているだけであ
り、当該セグメントの傾き（双曲線形状の傾きに相当）
は全く考慮されていないため、電磁波の地中伝搬速度Ｖ
を正しく推定することができない。例えば、伝搬速度Ｖ
＝０．２ｃ₀ （比誘電率２５）の地中において、地表か
ら深さ１００ｃｍに埋設された物体からの電磁波の反射
時間は３３ｎｓであるが、平均的な地中の電磁波の伝搬
速度Ｖ＝０．３３ｃ₀ （比誘電率９）で深さを計算する
と、深さは１６５ｃｍと大きな誤差を生じてしまう。
尚、ｃ₀ は真空中の光速である（以下、本明細書中にお
いて同じ。）。従って、埋設物の深さに相当する送受信
アンテナからの正確な深さ方向の距離Ｙをリアルタイム
で求めることができないことが分かる。

【０００４】また、かかる深さ検出ができないという問
題を解決すべく、反射信号を保存する度に、前記移動方
向前方への前記送受信アンテナからの距離Ｘ、前記移動
方向に垂直な方向への前記送受信アンテナからの距離
Ｙ、波動信号の地中における伝搬速度Ｖ毎に加算の規則
を変更しながら加算処理を行っていては、処理時間が長
くなり、探査速度が遅くなってしまうという問題があ
る。加えて、ブザーを鳴動させるか否かの判定を行う際
に、伝搬速度が遅くなるような土壌（電磁波の減衰が大
きい粘土に相当）の場合に、反射信号強度が弱くなり、
埋設物を検知できないという問題が生じる。同様に、前
記移動方向に垂直な方向への前記送受信アンテナからの
距離Ｙが大きい場合も、距離が遠い分だけ反射信号強度
が弱くなるため、同様の問題が発生する。

【０００５】更に、上記従来方法を、媒質中で送受信ア
ンテナを回転させ且つその移動方向を中心軸として回転
移動させながら探査する方法に適用する場合、回転に対
する認識方法が配慮されていないため、埋設物が存在す
る方向を検出することができない。また、かかる回転す
る送受信アンテナを用いた探査に適応したリアルタイム
探査方法はなかった。

【０００６】本発明は上述のような従来技術の有する問
題点を解消するためになされたものであり、その目的
は、媒質中の物体の存在をリアルタイムで検出でき、同
時に、物体までの深さ検出も短時間で行えるリアルタイ
ム探査方法及び装置、更には、媒質中で送受信アンテナ
を回転移動させながら探査が行えるリアルタイム探査方
法及び装置を提供することを目的としている。

【０００７】

【課題を解決するための手段】この目的を達成するため
の本発明に係るリアルタイム探査方法の第一の特徴構成
は、特許請求の範囲の欄の請求項１に記載した通り、媒
質の表面または媒質中で送受信部を移動させながら、前
記送受信部から電磁波または音波による波動信号を前記
媒質中へ放射し、前記媒質中に存在する物体からの反射
信号を前記送受信部が受信することにより、前記媒質中
に存在する物体を探査するリアルタイム探査方法であっ
て、前記送受信部の所定個所の現在位置を基準とする前
記送受信部の移動方向前方への距離Ｘ_i と前記移動方向
に垂直な方向への距離Ｙ_j 及び前記波動信号の前記媒質
中における伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）と
する３次元配列を予め作成する配列作成工程と、前記送
受信部を移動させながら、所定の時間間隔または所定の
移動距離毎に前記反射信号の反射時間と信号強度を対応
付けて保存する信号保存工程と、前記信号保存工程が実
行される度に、前記３次元配列の各ボクセル値を−Ｘ方
向へ１ボクセル分シフトするとともに、前記シフト後の
各ボクセル値に対して、座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決定
される反射時間に対応する前記信号保存工程で保存され
た信号強度を単純加算または加重加算する信号加算工程
と、前記信号加算工程が実行される度に、前記３次元配
列の各（Ｘ_i ，Ｙ_j ）座標におけるＶ軸方向の最大ボク
セル値の何れかまたはその中の最大値が所定の閾値を超
えたときに、前記媒質中に物体が存在すると判定する物
体検出工程とを備えてなる点にある。

【０００８】尚、本明細書中、座標（Ｘ_i ，Ｙ_j ，Ｖ
_k ）とする３次元配列は、Ｘ_i ，Ｙ_j，Ｖ_k を夫々引数
とする３次元配列と同意である。ここで、ｉ，ｊ，ｋは
夫々座標値Ｘ_i ，Ｙ_j ，Ｖ_k の添字で例えば自然数であ
り、ｉ，ｊ，ｋを夫々特定することで、座標値Ｘ_i ，Ｙ
_j ，Ｖ_k が夫々特定され、更には、３次元配列の一つの
配列要素（ボクセルという）が特定される。また、Ｘ方
向またはＸ軸方向という場合、添字ｊ，ｋが共通するボ
クセル群の配列方向と同意である。＋Ｘ方向とは、添字
ｉの増加する方向で、−Ｘ方向とは、添字ｉの減少する
方向である。Ｙ及びＶについても同様である。

【０００９】同第二の特徴構成は、特許請求の範囲の欄
の請求項２に記載した通り、媒質中で送受信部をその移
動方向を中心軸として回転移動させながら、前記送受信
部から電磁波または音波による波動信号を前記媒質中へ
放射し、前記媒質中に存在する物体からの反射信号を前
記送受信部が受信することにより、前記媒質中に存在す
る物体を探査するリアルタイム探査方法であって、前記
送受信部の所定個所の現在位置を基準とする前記送受信
部の移動方向前方への距離Ｘ_i と前記移動方向に垂直な
方向への距離Ｙ_j 及び前記波動信号の前記媒質中におけ
る伝搬速度Ｖ_kを座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）とする３次
元配列を前記送受信部の所定の複数の回転角に対応させ
て予め複数作成する配列作成工程と、前記送受信部を回
転移動させながら、前記送受信部の回転角が前記所定の
回転角となる毎に前記反射信号の反射時間と信号強度を
対応付けて保存する信号保存工程と、前記信号保存工程
が実行される度に、前記複数の３次元配列の内の前記信
号保存工程実行時の回転角に対応する一の３次元配列の
各ボクセル値を−Ｘ方向へ１ボクセル分シフトするとと
もに、前記シフト後の各ボクセル値に対して、座標値Ｘ
_i ，Ｙ_j ，Ｖ_k で決定される反射時間に対応する前記信
号保存工程で保存された信号強度を単純加算または加重
加算する信号加算工程と、前記信号加算工程が実行され
る度に、前記一の３次元配列の各（Ｘ_i ，Ｙ_j ）座標に
おけるＶ軸方向の最大ボクセル値の何れかまたはその中
の最大値が所定の閾値を超えたときに、前記媒質中に物
体が存在すると判定する物体検出工程とを備えてなる点
にある。

【００１０】ここで、送受信部の回転角とは、送受信部
の基準回転姿勢からの回転軸芯回りの回転角である。ま
た、送受信部の所定の複数の回転角とは、例えば、基準
回転姿勢の回転角０°と、その後３０°刻みで３０°〜
３３０°の値をとる。尚、刻み幅（角度）は任意であ
り、また、必ずしも等間隔でなくてもよい。

【００１１】同第三の特徴構成は、特許請求の範囲の欄
の請求項３に記載した通り、上述の第二の特徴構成に加
えて、前記配列作成工程において、前記信号保存工程が
実行される各回転角に対応する前記複数の３次元配列以
外に、これらと同じ大きさの追加の３次元配列を作成
し、前記信号加算工程が実行される度に、前記複数の３
次元配列の各ボクセル毎の総和を、前記追加の３次元配
列において計算する配列加算工程を実行し、前記信号保
存工程実行時の回転角に対応する一の３次元配列を使用
する代わりに、前記配列加算工程で得られた前記追加の
３次元配列を使用して前記物体検出工程を実行する点に
ある。

【００１２】同第四の特徴構成は、特許請求の範囲の欄
の請求項４に記載した通り、上述の第二または第三の特
徴構成に加えて、前記物体検出工程において、前記所定
の閾値を超えた前記３次元配列のボクセル値の座標（Ｘ
_i ’，Ｙ_j ’，Ｖ_k ’）に対応する前記各回転角毎の前
記複数の３次元配列のボクセル値を比較し、そのボクセ
ル値が最大となる３次元配列に対応する前記回転角を前
記物体が存在する方位であると判定する方位検出工程を
実行する点にある。

【００１３】同第五の特徴構成は、特許請求の欄の請求
項５に記載した通り、上述の第一、第二、第三または第
四の特徴構成に加えて前記座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決
定される反射時間は、前記３次元配列と同じ大きさのル
ックアップテーブルに予め計算されて格納されているも
のを使用する点にある。尚、前記座標値Ｘ_i ，Ｙ_j ，Ｖ
_k で決定される反射時間は、前記信号加算工程において
使用される反射時間である。

【００１４】同第六の特徴構成は、特許請求の範囲の欄
の請求項６に記載した通り、上述の第一、第二、第三、
第四または第五の特徴構成に加えて、前記信号加算工程
において加重加算を実行する場合、加算するボクセルの
座標値Ｘ_i ，Ｙ_j ，Ｖ_k の内の少なくとも一つの座標値
に基づいて一義的に決定される振幅補正値を係数として
前記信号強度に乗じてから加算する点にある。

【００１５】同第七の特徴構成は、特許請求の範囲の欄
の請求項７に記載した通り、上述の第一、第二、第三、
第四または第五の特徴構成に加えて、前記信号加算工程
において単純加算を実行する場合、後続の前記物体検出
工程において、前記所定の閾値として、判定対象のボク
セルの座標値Ｘ_i ，Ｙ_j ，Ｖ_k の内の少なくとも一つの
座標値に基づいて一義的に決定される閾値を用いる点に
ある。

【００１６】同第八の特徴構成は、特許請求の範囲の欄
の請求項８に記載した通り、上述の第六または第七の特
徴構成に加えて、前記座標値Ｘ_i ，Ｙ_j ，Ｖ_k の内の少
なくとも一つの座標値に基づいて一義的に決定される値
は、前記３次元配列と同じ大きさのルックアップテーブ
ルに予め計算されて格納されているものを使用する点に
ある。

【００１７】前記座標値Ｘ_i ，Ｙ_j ，Ｖ_k の内の少なく
とも一つの座標値に基づいて一義的に決定される値と
は、上記第六の特徴構成を備える場合、即ち前記信号加
算工程において加重加算を実行する場合は、前記信号加
算工程において使用される前記振幅補正値のことであ
り、上記第七の特徴構成を備える場合、即ち前記信号加
算工程において単純加算を実行する場合は、前記物体検
出工程において使用される前記閾値のことである。

【００１８】この目的を達成するための本発明に係るリ
アルタイム探査装置の第一の特徴構成は、特許請求の範
囲の欄の請求項９に記載した通り、電磁波または音波に
よる波動信号を送受信可能な送受信部を備え、媒質の表
面または媒質中で前記送受信部を移動させながら、前記
波動信号を前記媒質中へ放射し、前記媒質中に存在する
物体からの反射信号を受信することにより、前記媒質中
に存在する物体を探査するリアルタイム探査装置であっ
て、前記送受信部の所定個所の現在位置を基準とする前
記送受信部の移動方向前方への距離Ｘ_i と前記移動方向
に垂直な方向への距離Ｙ_j 及び前記波動信号の前記媒質
中における伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）と
する３次元配列を作成する配列作成手段と、前記送受信
部を移動させながら、所定の時間間隔または所定の移動
距離毎に前記反射信号の反射時間と信号強度を対応付け
て保存する信号保存手段と、前記信号保存手段が前記保
存処理を実行する度に、前記３次元配列の各ボクセル値
を−Ｘ方向へ１ボクセル分シフトするとともに、前記シ
フト後の各ボクセル値に対して、座標値Ｘ_i ，Ｙ_j，Ｖ_k
で決定される反射時間に対応する前記信号保存手段が
保存した信号強度を単純加算または加重加算する信号加
算手段と、前記信号加算手段が前記加算処理を実行する
度に、前記３次元配列の各（Ｘ_i ，Ｙ_j ）座標における
Ｖ軸方向の最大ボクセル値の何れかまたはその中の最大
値が所定の閾値を超えたときに、前記媒質中に物体が存
在すると判定する物体検出手段とを備えてなる点にあ
る。

【００１９】同第二の特徴構成は、特許請求の範囲の欄
の請求項１０に記載した通り、電磁波または音波による
波動信号を送受信可能な送受信部を備え、媒質中で前記
送受信部をその移動方向を中心軸として回転移動させな
がら、前記波動信号を前記媒質中へ放射し、前記媒質中
に存在する物体からの反射信号を受信することにより、
前記媒質中に存在する物体を探査するリアルタイム探査
装置であって、前記送受信部の所定個所の現在位置を基
準とする前記送受信部の移動方向前方への距離Ｘ_i と前
記移動方向に垂直な方向への距離Ｙ_j 及び前記波動信号
の前記媒質中における伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ
_j ，Ｖ_k ）とする３次元配列を前記送受信部の所定の複
数の回転角に対応させて複数作成する配列作成手段と、
前記送受信部を回転移動させながら、前記送受信部の回
転角が前記所定の回転角となる毎に前記反射信号の反射
時間と信号強度を対応付けて保存する信号保存手段と、
前記信号保存手段が前記保存処理を実行する度に、前記
複数の３次元配列の内のその保存処理実行時の回転角に
対応する一の３次元配列の各ボクセル値を−Ｘ方向へ１
ボクセル分シフトするとともに、前記シフト後の各ボク
セル値に対して、座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決定される
反射時間に対応する前記信号保存手段が保存した信号強
度を単純加算または加重加算する信号加算手段と、前記
信号加算手段が前記加算処理を実行する度に、前記一の
３次元配列の各（Ｘ_i ，Ｙ_j ）座標におけるＶ軸方向の
最大ボクセル値の何れかまたはその中の最大値が所定の
閾値を超えたときに、前記媒質中に物体が存在すると判
定する物体検出手段とを備えてなる点にある。

【００２０】以下に上記特徴構成の作用及び効果を説明
する。上記リアルタイム探査方法に係る第一の特徴構成
によれば、探査対象の範囲において、且つ、必要となる
精度で、前記送受信部の所定個所、例えば送信アンテナ
と受信アンテナの中間点の現在位置を基準とする前記送
受信部の移動方向前方への距離Ｘ_i と前記移動方向に垂
直な方向への距離Ｙ_j 及び前記波動信号の前記媒質中に
おける伝搬速度Ｖ_k の三つのパラメータの範囲と刻みを
設定して、前記３次元配列を作成し、前記３次元配列に
上記した手順で座標値Ｘ_i ，Ｙ_j ，Ｖ_kで決定される反
射時間に対応する前記信号保存手段が保存した信号強度
を加算する加算処理を行うため、この加算処理により、
伝搬速度Ｖ_k が媒質中の伝搬速度と一致する場合には、
媒質中に物体が存在する場合に反射信号強度の最大値が
距離Ｘと反射時間ｔを座標とする２次元平面上で形成す
る双曲線に沿って、当該最大値を加算する処理に相当す
ることから、前記３次元配列のボクセル値が所定の閾値
を超えた時に前記媒質中に物体が存在するものと判定す
ることができる。逆に、媒質中に物体が存在しない場合
は、反射信号強度の最大値が距離Ｘと反射時間ｔを座標
とする前記２次元平面上で双曲線を形成しないため、設
定範囲内の伝搬速度Ｖ_k では加算処理毎に増加していく
ボクセル値が存在しないため、前記媒質中に物体が存在
するものとは判定されない。更に、前記閾値を超えたボ
クセル値に対応する座標値Ｙ_j ’より、埋設物の深さに
相当する前記移動方向に垂直な方向への前記送受信部か
ら物体までの距離Ｙをリアルタイムで求めることができ
る。

【００２１】同第二の特徴構成によれば、前記送受信部
の所定個所の現在位置を基準とする前記送受信部の移動
方向前方への距離Ｘ_i と前記移動方向に垂直な方向への
距離Ｙ_j 及び前記波動信号の前記媒質中における伝搬速
度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）とする３次元配列を
前記送受信部の所定の複数の回転角に対応させて予め複
数作成し、前記送受信部の回転角が前記所定の回転角と
なる毎に、上記第一の特徴構成と同様の信号保存工程、
信号加算工程、物体検出工程を実行することから、各回
転角について、上記第一の特徴構成と同様の作用効果を
発揮し、回転軸芯回りのどの方向に物体が存在していて
も、物体の存在をリアルタイムで検出できるとともに、
物体の回転軸芯からの離隔距離Ｙ（埋設物の深さに相
当）をリアルタイムで求めることができる。

【００２２】更に、同第三の特徴構成によれば、前記波
動信号の放射方向の中心から逸れたところに物体が存在
しても反射信号を受信することから、前記物体の存在す
る方位と一致する或いは近い回転角以外の回転角につい
ての信号加算工程においても物体の存在に対応して媒質
中の伝搬速度と一致する伝搬速度Ｖ_k のボクセルは加算
処理毎にボクセル値が増加するため、配列加算工程によ
り、全ての回転角に対応する３次元配列の処理結果を加
算することで、所謂Ｓ／Ｎ比の向上が図られ、物体の存
在の判定をより精度良く行うことができる。

【００２３】また、第四の特徴構成によれば、物体が存
在する方位もリアルタイムで検出することができる。

【００２４】第五の特徴構成によれば、前記信号加算工
程において、前記反射時間の計算を、その都度実行する
必要がなくなり、各工程の処理時間の短縮化が図れるの
で、探査速度を速くすることができる。ところで、物体
までの距離、特に回転軸芯からの離隔距離Ｙが遠い場合
には反射信号が減衰しているため、前方距離Ｘが徐々に
短くなるにもかかわらず、ボクセル値の加算が十分でな
い場合があり得る。また、伝搬速度が遅くなるような媒
質中（土壌）でも反射信号が減衰しているため、ボクセ
ル値の加算が十分でない場合があり得る。かかる場合
は、Ｓ／Ｎ比が低下することになり、物体の存在の判定
精度が低下する。

【００２５】しかし、第六の特徴構成によれば、反射強
度が減衰している場合でも、離隔距離Ｙ_j や伝搬速度Ｖ
_k に応じて振幅補正することで、物体までの距離が遠い
場合や伝搬速度が遅くなるような媒質中でも、物体の存
在の判定を精度良く行うことができる。

【００２６】また、第七の特徴構成によれば、物体まで
の距離が遠い場合や伝搬速度が遅くなるような媒質中で
も、その状況に応じて、前記物体検出工程において使用
する閾値を調整することができ、上記第六の特徴構成と
同様に物体の存在の判定を精度良く行うことができる。

【００２７】第八の特徴構成によれば、前記信号加算工
程において、前記振幅補正値の計算、或いは、前記物体
検出工程において前記閾値の計算を、その都度実行する
必要がなくなり、各工程の処理時間の短縮化が図れるの
で、探査速度を速くすることができる。

【００２８】上記リアルタイム探査装置に係る第一の特
徴構成は、各手段が対応する上記リアルタイム探査方法
に係る第一の特徴構成の各工程を実行可能であるため、
上記リアルタイム探査方法に係る第一の特徴構成と同様
の作用効果を奏することができる。

【００２９】また、上記リアルタイム探査装置に係る第
二の特徴構成は、各手段が対応する上記リアルタイム探
査方法に係る第二の特徴構成の各工程を実行可能である
ため、上記リアルタイム探査方法に係る第二の特徴構成
と同様の作用効果を奏することができる。

【００３０】

【発明の実施の形態】以下に本発明のリアルタイム探査
方法及び装置実施の形態を図面に基づいて説明する。図
１に示すように、媒質である土壌１にガス等の流体を配
送する鋼管などの物体２が埋設されており、送受信手段
である送受信機１０とデータ解析手段であるデータ解析
装置２０を備えたリアルタイム探査装置３の本体部３ａ
に接続された送受信部１６を付帯したボーリング装置の
回転先端部１７が地中に向かって回転移動しながら、前
記物体２の埋設位置を探査する。前記送受信部１６は送
信アンテナ１１と受信アンテナ１２からなり、前記送信
アンテナ１１と前記受信アンテナ１２は一定間隔で前記
回転推進先端部１７の先端部分に配置されている。前記
送受信機１０は例えば１００ＭＨｚ〜１ＧＨｚの図２
（１）に例示する単発のパルス信号を送信回路１３で発
生し、前記送信アンテナ１１より電磁波として土壌１に
放射する。

【００３１】前記回転推進先端部１７はその長手方向の
中心軸を回転軸芯として回転しながらその軸芯方向（図
１中のｘ方向）に向かって移動するため、前記回転推進
先端部１７の回転に伴って、前記送信アンテナ１１から
放射される電磁波の放射方向も前記回転軸芯回りに回転
する。

【００３２】前記送信アンテナ１１より放射された電磁
波の中の物体に入射した入射波４は物体２表面で反射散
乱し、その中の反射波５が前記受信アンテナ１２で受信
された後、前記送受信機１０の受信回路１４において、
図２（２）に例示するような反射信号として復調増幅さ
れる。前記送信アンテナ１１より放射され、受信アンテ
ナ１２で受信されるまでの時間差ΔＴは前記送受信部１
６の前記送信アンテナ１１と前記受信アンテナ１２の中
間点から物体２までの距離と土壌１の比誘電率εまたは
電磁波の伝搬速度ｖより一義的に決定される。

【００３３】前記送受信機１０に、前記受信回路１４の
増幅部の利得を前記時間差ΔＴに応じて変調する信号強
度変調手段１５を設け、前記時間差ΔＴが長くなるにつ
れて土壌１を伝搬する前記パルス信号の損失が大きくな
り、反射信号強度が減衰するのを振幅補正し、前記時間
差ΔＴ、つまりは反射時間ｔの増加に対して急激に減衰
しない反射信号強度分布を得ることができ、次段以降の
信号処理に必要な信号強度を確保する。具体的には、前
記信号強度変調手段１５は前記単発のパルス信号の送信
タイミングに同期して、前記時間差ΔＴの増加に伴い減
衰率を所定の変化率で自動的に低下させる減衰器で構成
され、前記増幅部の所定個所に挿入してある。

【００３４】データ解析装置２０はマイクロコンピュー
タや半導体メモリ等によって構成されるデータ処理部２
１と外部からの操作指示を入力するキーボード等の入力
部２２と各処理段階での画像データや出力結果を表示す
る陰極線管ディスプレイや液晶ディスプレイ等の表示部
２３と前記各処理段階での画像データや出力結果等を保
管格納する磁気ディスク等の外部補助記憶部２４から構
成されている。受信回路１４において反射信号は、前記
信号強度変調手段１５による振幅補正後に、波形のスム
ージング等の雑音除去処理やＡ／Ｄ変換処理等の前置処
理が施され、ディジタル信号として前記データ処理部２
１へ出力される。

【００３５】前記データ解析装置２０では前記ディジタ
ル化された反射信号より、前記送受信部１６の移動距離
ｘと前記反射波５の前記物体２からの反射時間ｔを座標
（ｘ，ｔ）とする２次元画像データが、前記送受信部１
６の移動に伴い逐次生成される。尚、前記送受信部１６
は、ｘ方向への移動に伴い回転するので、その回転角α
はα＝ｆ（ｘ）で表される移動距離ｘの関数となり、移
動距離ｘに応じて回転角αの方位に埋設された物体２か
らの反射波を受信することとなる。ここで、反射信号強
度を−１２７〜１２８の２５６階調で輝度表示し、図２
（２）に示すように、信号強度の正値を白（輝度大）、
信号強度の負値を黒（輝度小）、信号強度０を中間階調
として表示部２３に表示する。

【００３６】図３に、伝搬速度Ｖ＝０．４５ｃ₀ （比誘
電率５）の地中において、地表面から深さｘ＝１００ｃ
ｍ、ボーリング位置から水平距離ｙ＝２５ｃｍに埋設さ
れた物体を、前記送受信部１６の移動距離ｘが０．０６
ｃｍ刻みで地表面から１５０ｃｍ移動した時に探査して
得られた２次元画像例を示す。但し、図３では、高輝度
部分のみを黒色で表示し、低輝度及び中間階調部分は紙
面地色のままである。また、本実施形態の前記回転推進
先端部１７は移動距離ｘが０．０６ｃｍ進む毎に３０°
回転するので、０．７２ｃｍ移動した時に丁度１回転す
る。表１に本実施形態における移動距離ｘと回転角αの
関係の一部を示す。ところで、図３より、前記送受信部
１６が前記物体の存在する方位に向いている時に専ら反
射信号を受信するため、高輝度部分が断続的に破線状に
なっていることが分かる。

【００３７】

【表１】

【００３８】図４に示すように、前記データ処理部２１
は、制御部３０、配列作成手段３１、信号保存手段３
２、信号加算手段３３、配列加算手段３４、物体検出手
段３５、方位検出手段３６、２次元画像データ生成手段
３７、出力処理手段３８、内部データバス３９、制御・
アドレスバス４０等を備えて構成されている。

【００３９】前記配列作成手段３１、前記信号保存手段
３２、前記信号加算手段３３、前記配列加算手段３４、
前記物体検出手段３５及び前記方位検出手段３６は、夫
々後述する配列作成工程、信号保存工程、信号加算工
程、配列加算工程、物体検出工程及び方位検出工程を個
別に実行するための手段である。前記制御部３０は上記
各手段及びデータの入出力を管理制御する。前記２次元
画像データ生成手段３７は図３に例示したような２次元
画像データを生成し前記表示部２３に出力する。また、
前記出力処理手段３８は後述する種々の探査結果を前記
表示部２３へ出力したり、必要に応じて音声出力等する
ための処理手段である。尚、上述のように、前記データ
処理部２１は、マイクロコンピュータや半導体メモリ等
によって具体的に構成されるが、図４に示す各機能手段
はこれらマイクロコンピュータや半導体メモリ等の一部
または全部を使用して、内部データバス３９、制御・ア
ドレスバス４０によって有機的に結合され、所定のプロ
グラムを実行することで実現される。

【００４０】次に、前記データ処理部２１における典型
的なデータ処理手順及び各部の動作を、図３に示す探査
データの探査対象（媒質及び物体）を探査した実施例
（以下、第１実施例という。）について、図５に示すフ
ローチャートに基づいて説明する。

【００４１】図５に示す配列作成工程において、前記配
列作成手段３１が前記送受信部１６の前記中間点を基準
とする前記送受信部の移動方向前方への距離Ｘ_i と前記
移動方向に垂直な方向への距離Ｙ_j 及び前記波動信号の
前記媒質中における伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，
Ｖ_k ）とする３次元配列Ｄ_l を、前記送受信部１６の所
定の複数の回転角α_l （ｌは自然数）に対応させて複数
作成する。本実施形態では、上述の如く前記回転推進先
端部１７は移動距離ｘが０．０６ｃｍ増加する毎に３０
°回転するので、前記所定の複数の回転角α_l は、隣接
する回転角の間隔を３０°に設定し、０°〜３３０°の
１２通りとした。即ち、α ₁＝０°、α₂ ＝３０°、…
………、α₁₂＝３３０°となる。

【００４２】また、距離Ｘ_i とＹ_j の範囲を、２５ｃｍ
≦Ｘ_i ≦４９．４８ｃｍ、０ｃｍ≦Ｙ_j ≦３５ｃｍと
し、距離Ｘ_i とＹ_j の刻みを夫々０．７２ｃｍと２．５
ｃｍとした。尚、距離Ｘ_i の刻みを０．７２ｃｍとした
のは、前記送受信部１６の回転角αがある回転角α_l か
ら同じ回転角α_l まで１回転するのに要する移動距離ｘ
が０．７２ｃｍとなるからである。

【００４３】また、伝搬速度Ｖ_k は数１で定義され、ｋ
は１〜１１の自然数とした。

【００４４】

【数１】Ｖ_k ＝ｃ₀ ×１０^-0.04k-0.32

【００４５】従って、３次元配列Ｄ_l （Ｘ_i ，Ｙ_j ，Ｖ
_k ）の大きさは、３５×１５×１１となり、この大きさ
の３次元配列Ｄ_l が前記回転角α_l （ｌ＝１〜１２）に
対応して１２個作成される。

【００４６】更に、この配列作成工程において、後述す
る配列加算工程で使用する３次元配列Ｄ_l （ｌ＝１〜１
２）と同じ大きさの３次元配列Ｄ₁₃（Ｘ_i ，Ｙ_j ，Ｖ
_k ）を準備しておく。

【００４７】図５に示す信号保存工程において、前記信
号保存手段３２が、前記移動距離ｘが０．０６ｃｍ増加
し、前記送受信部１６の回転角αが前記回転角α_l と一
致する毎に逐次、図２（２）で表されるような物体から
の反射信号の前記受信回路１４においてディジタル化さ
れた所定の反射時間ｔの範囲内のものを、サンプリング
して保存する。本実施形態では、前記反射信号のサンプ
リングは、反射時間ｔが０〜２０ｎｓの範囲を０．０７
８ｎｓの時間間隔で行われ、１回の信号保存処理で２５
６点の信号強度が保存される。

【００４８】次に、前記信号保存工程が実行される度
に、図５に示す信号加算工程において、前記信号加算手
段３３が、前記複数の３次元配列Ｄ_l の内の前記信号保
存工程実行時の回転角α_l に対応する一の３次元配列の
各ボクセル値を−Ｘ方向へ１ボクセル分シフトするとと
もに、前記シフト後の各ボクセル値に対して、座標値Ｘ
_i ，Ｙ_j ，Ｖ_k で決定される反射時間に対応する前記信
号保存工程で保存された前記信号強度を加重加算する。

【００４９】上記シフト処理と加重加算は、数２及び数
３に示す要領で同時に処理される。尚、数２はｉ＝１〜
３４の各ボクセルに対する処理を、数３はｉ＝３５の各
ボクセルの処理を示しており、Ｂ（ｉ，ｊ，ｋ）は座標
（Ｘ_i ，Ｙ_j ，Ｖ_k ）のボクセル値を表している。ここ
で、数２の右辺第１項がシフト処理を、数２の右辺第２
項が加重加算を夫々表しており、数３の場合は、シフト
処理で０が代入されるため、右辺がシフト処理と加重加
算を表している。

【００５０】

【数２】Ｂ（ｉ，ｊ，ｋ）＝Ｂ（ｉ＋１，ｊ，ｋ）＋Ａ
_i,j,k ×ｒ（ｔ_i,j,k ）

【数３】 Ｂ（３５，ｊ，ｋ）＝Ａ_i,j,k ×ｒ（ｔ_35,j,k）

【００５１】数２及び数３において、ｒ（ｔ_i,j,k ）
は、前記信号保存工程において保存された前記反射信号
強度の反射時間ｔにおける信号強度である。また、反射
時間ｔ _i,j,k は座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決定される反
射時間で、数４で表される。

【００５２】

【数４】ｔ_i,j,k ＝２（Ｘ_i ² ＋Ｙ_j ² ）^1/2 ／Ｖ_k

【００５３】更に、Ａ_i,j,k は加重加算処理における係
数で、加算される前記反射信号強度ｒ（ｔ_i,j,k ）の振
幅補正値である。本実施形態では、Ａ_i,j,k は数５で表
され、伝搬速度が遅くなるような媒質中で反射信号強度
が減衰する場合においても、加算処理が適正に実行でき
るように対処している。従って、当該振幅補正値Ａ_i,
j,k は伝搬速度Ｖ_k のみの関数として表されている。

【００５４】

【数５】Ａ_i,j,k ＝Ｖ_k ^-1.5

【００５５】ところで、数２及び数３の演算処理におい
て、数４及び数５の計算を逐次実行しても構わないが、
各座標値Ｘ_i ，Ｙ_j ，Ｖ_k は、前記配列作成工程におい
て既に決定されているため、前記反射時間ｔ_i,j,k と前
記振幅補正値Ａ_i,j,k を予めルックアップテーブルとし
て計算したものを使用することで、加重加算処理の処理
時間を短縮化できる。

【００５６】次に、前記信号加算工程が実行される度
に、図５に示す配列加算工程において、前記配列加算手
段３４が、数６に示すように、前記回転角α_l （ｌ＝１
〜１２）に対応する１２個の３次元配列Ｄ_l （Ｘ_i ，Ｙ
_j ，Ｖ_k ）の各ボクセル値の総和を求め、前記他の一つ
の３次元配列Ｄ₁₃（Ｘ_i ，Ｙ_j ，Ｖ_k ）に代入する。

【００５７】

【数６】Ｄ₁₃（Ｘ_i ，Ｙ_j ，Ｖ_k ）＝Σ_l=1-12Ｄ_l （Ｘ
_i ，Ｙ_j ，Ｖ_k ）

【００５８】次に、図５に示す物体検出工程において、
前記物体検出手段３５が、前記配列加算工程で得られた
前記３次元配列Ｄ₁₃の全ボクセル値の最大値が所定の閾
値を超えた時に、前記媒質中に物体が存在すると判定す
る。

【００５９】図６に、前記閾値を６００００に設定した
場合の前記物体の検出結果を示す。図６において、縦軸
は前記物体までの推定離隔距離で距離Ｙ_j に相当し、横
軸は前記物体までの真の距離を示している。図中の◆印
は、前記媒質中に物体が存在すると判定されたときの物
体との位置関係を、前記３次元配列Ｄ₁₃の最大ボクセル
値を与える座標値Ｙ_j に基づいて表している。図６に示
すように、前記送受信部１６から前記物体までの真の距
離が約−４０ｃｍの地点（前記物体の約４０ｃｍ手前）
から前記物体の検出が開始しており、３次元配列Ｄ_l
（Ｘ_i ，Ｙ_j ，Ｖ_k ）の距離Ｘ_i の設定範囲内（２５ｃ
ｍ≦Ｘ_i≦４９．４８ｃｍ）で正しく検出されている。
また、推定された物体までの離隔距離も１２．５ｃｍ〜
３０ｃｍと、真の離隔距離２５ｃｍに対して妥当な値と
なっている。

【００６０】次に、図５に示す方位検出工程において、
前記方位検出手段３６が前記物体検出工程で前記所定の
閾値を超えた前記３次元配列Ｄ₁₃の最大ボクセル値の座
標（Ｘ_i ’，Ｙ_j ’，Ｖ_k ’）に対応する前記各回転角
α_l 毎の前記複数の３次元配列Ｄ_l のボクセル値Ｂ（Ｘ
_i ’，Ｙ_j ’，Ｖ_k ’）を比較し、そのボクセル値が最
大となる３次元配列Ｄ_l に対応する前記回転角α_l を前
記物体が存在する方位であると判定する。図３に示す探
査データに適用して得られた物体の方位（回転角α_l に
相当）毎の前記ボクセル値Ｂ（Ｘ_i ’，Ｙ_j ’，Ｖ
_k ’）の前記送受信部１６の移動距離ｘに伴う変化を図
７に示す。図７に示すように、全ての移動距離ｘにおい
て、回転角α₇ （＝１８０°）の時に、前記ボクセル値
が最大となっている。本第１実施例では、１８０°の方
位に前記物体が存在していることから妥当な結果となっ
ている。

【００６１】次に、本発明に係るリアルタイム探査方法
を同じ手順で他の媒質中に適用した第２実施例を示す。

【００６２】図８に、伝搬速度Ｖ＝０．２ｃ₀ （比誘電
率２５）の地中（媒質中）において、地表面から深さｘ
＝１００ｃｍ、ボーリング位置から水平距離ｙ＝２５ｃ
ｍに埋設された物体を、前記送受信部１６の移動距離ｘ
が０．０６ｃｍ刻みで地表面から１５０ｃｍ移動した時
に探査して得られた２次元画像例を示す。この媒質は伝
搬速度が遅くなるような土壌（電磁波の減衰が大きい粘
土に相当）であり、図８に示す画像より、反射信号強度
が弱くなっており、５０ｃｍ手前では反射信号が得られ
ていないことが分かる。

【００６３】図９に、第２実施例における物体検出工程
における、前記媒質中に埋設された物体の検出結果を示
す。図９は上記第１実施例の図６に対応するもので、縦
軸、横軸等の設定は同じである。

【００６４】図９に示すように、前記送受信部１６から
前記物体までの真の距離が約−２５ｃｍの地点（前記物
体の約２５ｃｍ手前）から物体の検出が開始しており、
３次元配列Ｄ_l （Ｘ_i ，Ｙ_j ，Ｖ_k ）の距離Ｘ_i の設定
範囲内（２５ｃｍ≦Ｘ_i ≦４９．４８ｃｍ）で検出され
ている。また、推定された物体までの離隔距離も２０ｃ
ｍ〜３５ｃｍと、真の離隔距離２５ｃｍに対して略妥当
な値となっている。

【００６５】このように、反射信号強度の減衰の大きな
媒質中にもかかわらず、物体の検出が可能となっている
のは、前記信号加算工程において、数５に示す前記振幅
補正値Ａ_i,j,k により振幅補正を行い、数２及び数３に
示す加重加算処理を行った効果である。従って、かかる
振幅補正を行わずに単純加算処理を行った場合は、閾値
を６００００に設定したのでは、検出不可能である。
尚、この閾値を下げると、伝搬速度が速く減衰が小さい
土壌の場合、物体までの距離がかなり遠い時にも検出さ
れてしまい、推定される物体までの離隔距離の算出誤差
が大きくなるという問題が発生する。

【００６６】更に、前記方位検出工程を実行した場合に
おける、前記物体の方位（回転角α _l に相当）毎の前記
ボクセル値Ｂ（Ｘ_i ’，Ｙ_j ’，Ｖ_k ’）の前記送受信
部１６の移動距離ｘに伴う変化を図１０に示す。図１０
に示すように、全ての移動距離ｘにおいて、回転角α₇
（＝１８０°）の時に、前記ボクセル値が最大となって
いる。本第２実施例では、１８０°の方位に埋設物が存
在していることから妥当な結果となっている。

【００６７】以下に他の実施形態を説明する。 〈１〉上記実施の形態において、前記配列加算工程と前
記物体検出工程を実行する代わりに、前記信号加算工程
が実行される度に、その加算処理が実行されたときの回
転角α_l に対応する３次元配列Ｄ_l （Ｘ_i ，Ｙ_j ，Ｖ
_k ）の各（Ｘ_i ，Ｙ _j ）座標におけるＶ軸方向の最大ボ
クセル値の何れかまたはその中の最大値が所定の閾値を
超えたときに、前記媒質中に物体が存在すると判定する
物体検出工程を実行しても構わない。

【００６８】〈２〉前記３次元配列Ｄ_l （Ｘ_i ，Ｙ_j ，
Ｖ_k ）の大きさや、距離Ｘ_i 、距離Ｙ _j 並びに伝搬速度
Ｖ_k の範囲及び刻み幅は、必ずしも上記実施の形態の値
に限定されるものではない。また、回転角α_l の総数や
設定値も、上記実施の形態の値に限定されるものではな
い。従って、前記３次元配列Ｄ_l の個数も上記実施の形
態の値に限定されるものではない。

【００６９】〈３〉上記実施の形態の前記信号加算工程
において、前記振幅補正値Ａ_i,j,k による振幅補正を行
わずに、数２及び数３における加重加算の係数Ａ_i,j,k
を１とする単純加算であっても構わない。媒質中の伝搬
速度が予測可能な場合等においては、上記したような振
幅補正を行う必要のない場合がある。また、かかる振幅
補正が必要な場合であっても、前記信号加算工程では単
純加算を行い、前記物体検出工程における前記所定の閾
値を、例えば、前記振幅補正値Ａ_i,j,k で前記定数の閾
値を除した値を使用するようにしてもよい。この場合、
この閾値を予めルックアップテーブル化するようにすれ
ば、判定時の計算時間が短縮できる。更に、前記振幅補
正値Ａ_i,j,k は必ずしも数５に示す補正式に依らなくて
も構わない。例えば、距離Ｙ_j 、或いは、距離Ｙ_j と伝
搬速度Ｖ_k に応じて変化するものであっても構わない。

【００７０】〈４〉上記実施の形態では、前記送受信部
１６はボーリング装置の回転推進先端部１７に付帯し、
媒質中を回転移動しながら、前記物体２の埋設位置を探
査する構成のものであったが、媒質表面を回転せずに移
動するものであっても構わない。また、この場合、前記
リアルタイム探査装置３の本体部３ａと別体で構成され
るのではなく、前記本体部と一体化していても構わな
い。前記送受信部１６が回転しない場合は、上記実施の
形態のリアルタイム探査方法は、図１１に示すように簡
略化することができる。

【００７１】また、図１１に示す配列作成工程では、３
次元配列は前記送受信部１６の回転角に対応させて複数
作成する必要がないため、３次元配列を一つだけを作成
する。また、３次元配列が一つだけであり、前記回転角
を検出する必要がないため、図５に示す前記配列加算工
程や前記方位検出工程は設ける必要がない。

【００７２】図５に示す前記信号保存工程が前記送受信
部１６の回転角αが前記回転角α_lと一致する毎に実行
されるのに対して、図１１に示す信号保存工程は、前記
送受信部１６の所定の移動時間間隔または所定の移動距
離毎に実行される点で異なる以外は、基本的に同じ処理
を行っている。また、信号加算工程は、常に同じ３次元
配列を使用する点で、図５に示す前記信号加算工程と異
なるが、加算処理自体は全く同じ処理を実行すればよ
い。更に、図１１に示す物体検出工程では、図５に示す
前記配列加算工程が不要であるため、当該工程で得られ
た前記他の一つの３次元配列Ｄ₁₃の代わりに、図１１に
示す信号加算工程で処理した３次元配列を使用すればよ
く、その他の実質的な処理は図５に示す物体検出工程と
同様に処理可能である。

【００７３】従って、前記送受信部１６が回転しない場
合は、上記実施の形態のリアルタイム探査装置３は、前
記配列加算手段３４と前記方位検出手段３６を必ずしも
備えている必要はない。また、前記配列加算手段３４と
前記方位検出手段３６を備え、前記配列作成手段３１、
前記信号保存手段３２、前記信号加算手段３３、及び、
前記物体検出手段３５が前記送受信部１６が回転する場
合と回転しない場合の両方に対応できるように構成され
ているのも好ましい。

【図面の簡単な説明】

【図１】本発明に係るリアルタイム探査装置のブロック
構成図

【図２】送信信号と受信信号の波形図

【図３】本発明に係るリアルタイム探査方法のデータ処
理手順に使用した２次元探査画像（第１実施例の探査対
象を表す）

【図４】本発明に係るリアルタイム探査装置のデータ解
析装置の機能ブロック図

【図５】本発明に係るリアルタイム探査方法の一実施の
形態のデータ処理手順を示すフローチャート

【図６】本発明に係るリアルタイム探査方法を図３に示
す探査対象に適用した物体検出工程における探査結果を
示す図

【図７】本発明に係るリアルタイム探査方法を図３に示
す探査対象に適用した方位検出工程における探査結果を
示す図

【図８】本発明に係るリアルタイム探査方法のデータ処
理手順に使用した２次元探査画像（第２実施例の探査対
象を表す）

【図９】本発明に係るリアルタイム探査方法を図８に示
す探査対象に適用した物体検出工程における探査結果を
示す図

【図１０】本発明に係るリアルタイム探査方法を図８に
示す探査対象に適用した方位検出工程における探査結果
を示す図

【図１１】本発明に係るリアルタイム探査方法の別実施
の形態のデータ処理手順を示すフローチャート

【図１２】従来技術のデータ処理手順を示すフローチャ
ート

【符号の説明】 １ 媒質 ２ 物体 ３ リアルタイム探査装置 ４ 入射波 ５ 反射波 １０ 送受信手段 １１ 送信アンテナ １２ 受信アンテナ １３ 送信回路 １４ 受信回路 １５ 信号強度変調手段 １６ 送受信部 １７ 回転推進先端部 ２０ データ解析装置 ２１ データ処理部 ２１ａ メモリ ２２ 入力部 ２３ 表示部 ２４ 外部補助記憶部 ３０ 制御部 ３１ 配列作成手段 ３２ 信号保存手段 ３３ 信号加算手段 ３４ 配列加算手段 ３５ 物体検出手段 ３６ 方位検出手段 ３７ ２次元画像データ生成手段 ３８ 出力処理手段

Claims (10)

【特許請求の範囲】
1. 【請求項１】 媒質の表面または媒質中で送受信部を移
   動させながら、前記送受信部から電磁波または音波によ
   る波動信号を前記媒質中へ放射し、前記媒質中に存在す
   る物体からの反射信号を前記送受信部が受信することに
   より、前記媒質中に存在する物体を探査するリアルタイ
   ム探査方法であって、 前記送受信部の所定個所の現在位置を基準とする前記送
   受信部の移動方向前方への距離Ｘ_i と前記移動方向に垂
   直な方向への距離Ｙ_j 及び前記波動信号の前記媒質中に
   おける伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）とする
   ３次元配列を予め作成する配列作成工程と、 前記送受信部を移動させながら、所定の時間間隔または
   所定の移動距離毎に前記反射信号の反射時間と信号強度
   を対応付けて保存する信号保存工程と、 前記信号保存工程が実行される度に、前記３次元配列の
   各ボクセル値を−Ｘ方向へ１ボクセル分シフトするとと
   もに、前記シフト後の各ボクセル値に対して、座標値Ｘ
   _i ，Ｙ_j ，Ｖ_k で決定される反射時間に対応する前記信
   号保存工程で保存された信号強度を単純加算または加重
   加算する信号加算工程と、 前記信号加算工程が実行される度に、前記３次元配列の
   各（Ｘ_i ，Ｙ_j ）座標におけるＶ軸方向の最大ボクセル
   値の何れかまたはその中の最大値が所定の閾値を超えた
   ときに、前記媒質中に物体が存在すると判定する物体検
   出工程とを備えてなるリアルタイム探査方法。
2. 【請求項２】 媒質中で送受信部をその移動方向を中心
   軸として回転移動させながら、前記送受信部から電磁波
   または音波による波動信号を前記媒質中へ放射し、前記
   媒質中に存在する物体からの反射信号を前記送受信部が
   受信することにより、前記媒質中に存在する物体を探査
   するリアルタイム探査方法であって、 前記送受信部の所定個所の現在位置を基準とする前記送
   受信部の移動方向前方への距離Ｘ_i と前記移動方向に垂
   直な方向への距離Ｙ_j 及び前記波動信号の前記媒質中に
   おける伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）とする
   ３次元配列を前記送受信部の所定の複数の回転角に対応
   させて予め複数作成する配列作成工程と、 前記送受信部を回転移動させながら、前記送受信部の回
   転角が前記所定の回転角となる毎に前記反射信号の反射
   時間と信号強度を対応付けて保存する信号保存工程と、 前記信号保存工程が実行される度に、前記複数の３次元
   配列の内の前記信号保存工程実行時の回転角に対応する
   一の３次元配列の各ボクセル値を−Ｘ方向へ１ボクセル
   分シフトするとともに、前記シフト後の各ボクセル値に
   対して、座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決定される反射時間
   に対応する前記信号保存工程で保存された信号強度を単
   純加算または加重加算する信号加算工程と、 前記信号加算工程が実行される度に、前記一の３次元配
   列の各（Ｘ_i ，Ｙ_j ）座標におけるＶ軸方向の最大ボク
   セル値の何れかまたはその中の最大値が所定の閾値を超
   えたときに、前記媒質中に物体が存在すると判定する物
   体検出工程とを備えてなるリアルタイム探査方法。
3. 【請求項３】 前記配列作成工程において、前記信号保
   存工程が実行される各回転角に対応する前記複数の３次
   元配列以外に、これらと同じ大きさの追加の３次元配列
   を作成し、 前記信号加算工程が実行される度に、前記複数の３次元
   配列の各ボクセル毎の総和を、前記追加の３次元配列に
   おいて計算する配列加算工程を実行し、 前記信号保存工程実行時の回転角に対応する一の３次元
   配列を使用する代わりに、前記配列加算工程で得られた
   前記追加の３次元配列を使用して前記物体検出工程を実
   行することを特徴とする請求項２記載のリアルタイム探
   査方法。
4. 【請求項４】 前記物体検出工程において、前記所定の
   閾値を超えた前記３次元配列のボクセル値の座標（Ｘ
   _i ’，Ｙ_j ’，Ｖ_k ’）に対応する前記各回転角毎の前
   記複数の３次元配列のボクセル値を比較し、そのボクセ
   ル値が最大となる３次元配列に対応する前記回転角を前
   記物体が存在する方位であると判定する方位検出工程を
   実行することを特徴とする請求項２または３記載のリア
   ルタイム探査方法。
5. 【請求項５】 前記座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決定され
   る反射時間は、前記３次元配列と同じ大きさのルックア
   ップテーブルに予め計算されて格納されているものを使
   用することを特徴とする請求項１、２、３または４記載
   のリアルタイム探査方法。
6. 【請求項６】 前記信号加算工程において加重加算を実
   行する場合、加算するボクセルの座標値Ｘ_i ，Ｙ_j ，Ｖ
   _k の内の少なくとも一つの座標値に基づいて一義的に決
   定される振幅補正値を係数として前記信号強度に乗じて
   から加算することを特徴とする請求項１、２、３、４ま
   たは５記載のリアルタイム探査方法。
7. 【請求項７】 前記信号加算工程において単純加算を実
   行する場合、後続の前記物体検出工程において、前記所
   定の閾値として、判定対象のボクセルの座標値Ｘ_i ，Ｙ
   _j ，Ｖ_k の内の少なくとも一つの座標値に基づいて一義
   的に決定される閾値を用いることを特徴とする請求項
   １、２、３、４または５記載のリアルタイム探査方法。
8. 【請求項８】 前記座標値Ｘ_i ，Ｙ_j ，Ｖ_k の内の少な
   くとも一つの座標値に基づいて一義的に決定される値
   は、前記３次元配列と同じ大きさのルックアップテーブ
   ルに予め計算されて格納されているものを使用すること
   を特徴とする請求項６または７記載のリアルタイム探査
   方法。
9. 【請求項９】 電磁波または音波による波動信号を送受
   信可能な送受信部を備え、媒質の表面または媒質中で前
   記送受信部を移動させながら、前記波動信号を前記媒質
   中へ放射し、前記媒質中に存在する物体からの反射信号
   を受信することにより、前記媒質中に存在する物体を探
   査するリアルタイム探査装置であって、 前記送受信部の所定個所の現在位置を基準とする前記送
   受信部の移動方向前方への距離Ｘ_i と前記移動方向に垂
   直な方向への距離Ｙ_j 及び前記波動信号の前記媒質中に
   おける伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）とする
   ３次元配列を作成する配列作成手段と、 前記送受信部を移動させながら、所定の時間間隔または
   所定の移動距離毎に前記反射信号の反射時間と信号強度
   を対応付けて保存する信号保存手段と、 前記信号保存手段が前記保存処理を実行する度に、前記
   ３次元配列の各ボクセル値を−Ｘ方向へ１ボクセル分シ
   フトするとともに、前記シフト後の各ボクセル値に対し
   て、座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決定される反射時間に対
   応する前記信号保存手段が保存した信号強度を単純加算
   または加重加算する信号加算手段と、 前記信号加算手段が前記加算処理を実行する度に、前記
   ３次元配列の各（Ｘ_i，Ｙ_j ）座標におけるＶ軸方向の
   最大ボクセル値の何れかまたはその中の最大値が所定の
   閾値を超えたときに、前記媒質中に物体が存在すると判
   定する物体検出手段とを備えてなるリアルタイム探査装
   置。
10. 【請求項１０】 電磁波または音波による波動信号を送
    受信可能な送受信部を備え、媒質中で前記送受信部をそ
    の移動方向を中心軸として回転移動させながら、前記波
    動信号を前記媒質中へ放射し、前記媒質中に存在する物
    体からの反射信号を受信することにより、前記媒質中に
    存在する物体を探査するリアルタイム探査装置であっ
    て、 前記送受信部の所定個所の現在位置を基準とする前記送
    受信部の移動方向前方への距離Ｘ_i と前記移動方向に垂
    直な方向への距離Ｙ_j 及び前記波動信号の前記媒質中に
    おける伝搬速度Ｖ_k を座標（Ｘ_i ，Ｙ_j ，Ｖ_k ）とする
    ３次元配列を前記送受信部の所定の複数の回転角に対応
    させて複数作成する配列作成手段と、 前記送受信部を回転移動させながら、前記送受信部の回
    転角が前記所定の回転角となる毎に前記反射信号の反射
    時間と信号強度を対応付けて保存する信号保存手段と、 前記信号保存手段が前記保存処理を実行する度に、前記
    複数の３次元配列の内のその保存処理実行時の回転角に
    対応する一の３次元配列の各ボクセル値を−Ｘ方向へ１
    ボクセル分シフトするとともに、前記シフト後の各ボク
    セル値に対して、座標値Ｘ_i ，Ｙ_j ，Ｖ_k で決定される
    反射時間に対応する前記信号保存手段が保存した信号強
    度を単純加算または加重加算する信号加算手段と、 前記信号加算手段が前記加算処理を実行する度に、前記
    一の３次元配列の各（Ｘ_i ，Ｙ_j ）座標におけるＶ軸方
    向の最大ボクセル値の何れかまたはその中の最大値が所
    定の閾値を超えたときに、前記媒質中に物体が存在する
    と判定する物体検出手段とを備えてなるリアルタイム探
    査装置。