JP2022020291A - 埋設物探査装置及び埋設物探査方法 - Google Patents
埋設物探査装置及び埋設物探査方法 Download PDFInfo
- Publication number
- JP2022020291A JP2022020291A JP2020123696A JP2020123696A JP2022020291A JP 2022020291 A JP2022020291 A JP 2022020291A JP 2020123696 A JP2020123696 A JP 2020123696A JP 2020123696 A JP2020123696 A JP 2020123696A JP 2022020291 A JP2022020291 A JP 2022020291A
- Authority
- JP
- Japan
- Prior art keywords
- buried object
- signal
- buried
- vector
- exploration
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims description 16
- 239000013598 vector Substances 0.000 claims abstract description 60
- 239000011159 matrix material Substances 0.000 claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 30
- 238000010586 diagram Methods 0.000 description 14
- 238000012937 correction Methods 0.000 description 10
- 239000013535 sea water Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 241000257303 Hymenoptera Species 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
【課題】電波の発射がなされる媒質の電気的特性が変化する環境下においても安定した埋設物の探査が行える埋設物探査装置を提供することを目的とする。【解決手段】探査の際にセンサによるセンサ信号を取得し、信号ベクトルとして扱われる電波によって埋設されている埋設物の位置を推定する推定位置を、位置候補ベクトルから選択することで埋設物を探査する埋設物探査装置であって、センサ信号及び前記信号ベクトルを取得する信号取得部と、センサ信号に基づいて応答行列を更新して信号処理を行う信号処理部と、を備え、応答行列は、位置候補ベクトルと信号ベクトルとを関連付ける行列である、埋設物探査装置を提供する。【選択図】図1
Description
本発明は、埋設物探査装置及び埋設物探査方法に関する。
例えば非特許文献1には、電波を用いて海底下埋設物センシングする手法が開示されている。
電波を用いた海底下センシングシステムの開発、電子情報通信学会ソサイエティ大会、2019
しかしながら非特許文献1に記載されたような従来の手法では、海水の電気的特性が変化した際に探査装置が受信する受信信号の振幅及び位相の変化量に変化が生じ探査性能が劣化するという課題がある。
本発明の実施の形態の一態様は、電波が送信及び受信される媒質の電気的特性が変化する環境下においても安定した探査が行える埋設物探査装置を提供することを目的とする。
探査の際にセンサによるセンサ信号を取得し、信号ベクトルとして扱われる電波によって埋設されている埋設物の位置を推定する推定位置を、位置候補ベクトルから選択することで埋設物を探査する埋設物探査装置であって、センサ信号及び前記信号ベクトルを取得する信号取得部と、センサ信号に基づいて応答行列を更新して信号処理を行う信号処理部と、を備え、応答行列は、位置候補ベクトルと信号ベクトルとを関連付ける行列である、埋設物探査装置を提供する。
探査の際にセンサによるセンサ信号を取得し、信号ベクトルとして扱われる電波によって埋設されている埋設物の位置を推定する推定位置を、位置候補ベクトルから選択することで埋設物を探査する埋設物探査装置による埋設物探査方法であって、センサ信号及び信号ベクトルを取得する第1のステップと、センサ信号に基づいて応答行列を更新して信号処理を行う第2のステップと、を備え、応答行列は、位置候補ベクトルと信号ベクトルとを関連付ける行列である、埋設物探査方法を提供する。
本発明の実施の形態の一態様によれば、海水中から送信する電波を用いて、海底下に埋設されている物体を、海底面に接触することなく探査することが可能な装置を実現できる。
以下図面を用いて、本発明の実施の形態の一態様を詳述する。
(第1の実施の形態)
まず図1を用いて本実施の形態による埋設物探査装置4について説明する。図1は、本実施の形態による埋設物探査装置4の構成を示すブロック図である。埋設物探査装置4は、例えば海水などの損失媒質1の中から、砂などの損失媒質2の中に埋設されている埋設物3の位置を推定する推定位置
を算出することで埋設物3を探査する。具体的には例えば埋設物探査装置4は、後述の信号ベクトルyと埋設物の位置候補ベクトルxと応答行列A(σ,εr,h)とから推定位置
を算出する。
まず図1を用いて本実施の形態による埋設物探査装置4について説明する。図1は、本実施の形態による埋設物探査装置4の構成を示すブロック図である。埋設物探査装置4は、例えば海水などの損失媒質1の中から、砂などの損失媒質2の中に埋設されている埋設物3の位置を推定する推定位置
を算出することで埋設物3を探査する。具体的には例えば埋設物探査装置4は、後述の信号ベクトルyと埋設物の位置候補ベクトルxと応答行列A(σ,εr,h)とから推定位置
を算出する。
埋設物探査装置4には損失媒質1の導電率を計測する導電率計測計5、埋設物探査装置4から海底面2Aまでの距離hを測定する高度計6及び埋設物探査装置4の傾斜mを計測する傾斜計7などのセンサが接続される。
埋設物探査装置4は、導電率σと比誘電率εr、距離h及び傾斜mを含むセンサ信号を導電率計測計5、高度計6及び傾斜計7から取得する。導電率計測計5は、例えばCTD(Conductivity Temperature Depth profiler)計のような測器としてもよい。
埋設物探査装置4には送信アンテナ8と受信アンテナ9~12が接続されており、受信アンテナ9~12は、送信アンテナ8が送信した送信信号に対する信号y1~y4を要素とする信号ベクトルyとして受信する。
より一般的な表現として信号ベクトルyの要素数は、受信アンテナ9~12の個数であるn個とする。n=4の場合、信号ベクトルyは信号y1~y4を要素とするベクトルとして表すことができる。
信号ベクトルyは(1)式のように示すことができる。信号ベクトルyの要素である信号y1~y4は送信アンテナ8と受信アンテナ9~12との間で電波が伝搬したことによる振幅及び位相の変化を示す。
………(1)
送信アンテナ8及び受信アンテナ9~12が扱う信号は、例えば周波数が1MHzである電波とし、後述の図6に示すように、送信アンテナ8及び受信アンテナ9~12は、送信アンテナ8が受信アンテナ9~12のいずれとの距離が等しくなるように配置される。
本実施の形態において、海水中における電波の減衰量が1mあたり35dB程度と非常に大きい周波数帯域である1MHzの電波を使用することから、埋設物3の個数は探査可能な範囲には1つであるとする。これによって、位置候補ベクトルxにはスパース性を導入することができる。具体的には位置候補ベクトルxのうち、埋設物の個数に相当する非ゼロ要素の数は1つしかないと考えることができる。
埋設物探査装置4は、信号取得部13と、信号処理部14と、を備える。信号取得部13は、信号ベクトルy及びセンサ信号を取得する。信号処理部14は、取得した信号ベクトルy及びセンサ信号と、予め保存されている位置候補ベクトルx={x1、…、x130}及び後述の初期応答行列A(σ0,εr0,h0)と、を用いて信号処理を行う信号処理部14とを備える。
信号取得部13は、FPGA(Field-Programmable Gate Array)などの集積回路とし、信号処理部14は、CPU(Central Processing Unit)によって呼び出されるRAM(Random Access Memory)に記録されたプログラムとする。
例えば位置候補ベクトルx={x1、…、x130}及び初期応答行列A(σ0,εr0,h0)などはRAMに予め保存されているものとする。信号取得部13と信号処理部14とは例えばLAN(Local Area Network)経由で情報の送受信を行うものとする。
埋設物探査装置4は、埋設物3が位置候補ベクトルx={x1、…、x130}によって表される例えば130個のボクセル群15のいずれか1つの位置に埋設されているものと仮定することで、埋設物3を探査する。
より一般的な表現として位置候補ベクトルの要素である位置候補x1~xlは、ボクセル群15の個数であるl個とする。このとき位置候補ベクトルx={x1、…、xl}として表すことができる。
応答行列A(σ,εr,h)は、位置候補ベクトルxと信号ベクトルyとを関係づける行列であって、位置候補ベクトルxと信号ベクトルyとを用いて(3)式のように表すことができる。なお応答行列A(σ,εr,h)は、金属や誘電体といったような埋設物3の材質に応じて変更する必要がある。
………(3)
図2を用いて本実施の形態における埋設物探査処理について説明する。図2は、本実施の形態による埋設物探査処理の処理手順を示すフローチャートである。信号取得部13は、信号ベクトルyを受信アンテナ9~12から取得する(S1)。
次に信号取得部13は、センサ信号を導電率計測計5、高度計6及び傾斜計7から導電率σと比誘電率εr、距離h及び傾斜mを取得する(S2)。信号取得部13は、距離hと傾斜mとから、それぞれの送信アンテナ8及び受信アンテナ9~12と海底面2Aとの距離h1~h5を算出する。距離hは例えば送信アンテナ8から海底面2Aまでの距離とする。
次に信号処理部14は応答行列A(σ,εr,h)を更新する(S3)。応答行列A(σ,εr,h)は、導電率σ、比誘電率εr及び距離hに依存する例えば4行130列の行列であって(3)式で表すことができる。ここで、演算子「.*」は、行列の要素同士の積(要素積)を求めることを意味する。応答行列A(σ,εr,h)の行数は受信アンテナ9~12の数である4であって、応答行列A(σ,εr,h)の列数はボクセル群15の数である130とする。
………(4)
ここで初期応答行列A(σ0,εr0,h0)は、あらかじめ測定した初期導電率σ0、初期比誘電率εr0及び初期距離h0に基づいて算出される行列であって行数と列数は応答行列A(σ,εr,h)と同じとする。第1の補正用行列f(σ,εr,h)及び第2の補正用行列g(σ,εr,h)は、導電率σ、比誘電率εr及び距離hに依存する例えば4行130列の行列とする。
応答行列A(σ,εr,h)を更新するには初期応答行列A(σ0,εr0,h0)に対して第1の補正と第2の補正とを行う。具体的にはまず初期応答行列A(σ0,εr0,h0)と第1の補正用行列f(σ,εr,h)との要素積を求めることで初期応答行列A(σ0,εr0,h0)に対して第1の補正として積による補正を行う。次に第1の補正を行った初期応答行列A(σ0,εr0,h0)に、第2の補正用行列g(σ,εr,h)を足すことで第2の補正として和による補正を行う。
算出されるノルム和NSに対して2次ノルムの影響が大きい場合、推定位置
の数が多くなる。本実施の形態においては、送信アンテナ8から送信される電波の届く範囲を考えると、検出可能な埋設物数は多くないと考えられるため、推定位置
を算出するにあたっての1次ノルムの影響を大きくしている。このように1次ノルムによって埋設物数に制約を設けることで過学習を抑制している。
の数が多くなる。本実施の形態においては、送信アンテナ8から送信される電波の届く範囲を考えると、検出可能な埋設物数は多くないと考えられるため、推定位置
を算出するにあたっての1次ノルムの影響を大きくしている。このように1次ノルムによって埋設物数に制約を設けることで過学習を抑制している。
なお埋設物の位置を推定する際は、ノルム和NSからさらに以下の(6)式を用いて尤度Lを算出してもよい。係数Kは、尤度Lが取りうる値の範囲を制限するための値であって、例えば255とする。係数Kを255とすると、ノルム和NSが0の際に尤度Lは最大値である255となる。
………(6)
なおステップS1とステップS2及びステップS3とは、例えば並列して処理されるため、ステップS1,S2,S3の順でなくてもよい。例えば、ステップS2,S3,S1の順でもよいし、ステップS2,S1,S3の順でもよい。
次に図3を用いて本実施の形態の効果について説明を行う。図3は、本実施の形態による探査結果と従来の探査結果とを比較する図である。図3においてグラフ21は既存の探査結果を示しており、グラフ22は誤差がない場合を示しており、グラフ23は本実施の形態による探査結果を示している。
図3における探査の際に使用した送信アンテナ8及び受信アンテナ9~12は微小ダイポールアンテナとし、図3における探査の際の埋設物3の材質は金属とした。グラフ21及びグラフ23は、埋設物3の位置を相対座標の(0cm,-60cm,-10cm)の地点から(0cm,60cm,-10cm)の地点まで20cm間隔でy軸方向に移動させながら、それぞれの際における探査結果を示す。
埋設物3の埋設深さは10cmとし、高度計6で測定される距離hは50cmとした。横軸の真値は、相対座標におけるy軸の値とし、縦軸の推定値は、埋設物探査装置4によって算出された埋設物3が埋設されていると推定されるボクセル位置に対するy軸の値を示す。
グラフ21に示す従来の探査結果における誤差としてRMSE(Root Mean Square Error)を求めたところ6.43cmとなった。これに対してグラフ21に示す本実施の形態による探査結果としてRMSEを求めたところ2.92cmとなり、半分以下の誤差となった。
より具体的には、グラフ21の探査では導電率σ0として4.30をそのまま使用したが、グラフ22の探査ではセンサ情報によって計測された導電率σとして4.41を使用した。
以上のように本実施の形態においては、埋設物探査装置4は、例えばステップS3によって探査の際に応答行列を更新し、埋設物探査処理を行う際の状態に即したパラメータを使用することができる。
このため、埋設物探査装置4は、従来に比べて精度よく埋設物3を探査することが可能となる。なお本実施の形態においては、送信アンテナ8及び受信アンテナ9~12ごとに高度計6を設置するなどして、送信アンテナ8及び受信アンテナ9~12ごとの海底面2Aまでの高度を測定できれば、傾斜計7は不要となる。
また本実施の形態においては、送信アンテナを1つとし、受信アンテナを複数としたがこれに限らず、送信アンテナは複数でもよいし、受信アンテナは1つでもよい。また本実施の形態において信号取得部13は集積回路とし、信号処理部14はプログラムとしたが、これに限らず、信号取得部13と信号処理部14とは、ともにプログラムであってもよいし、ともに集積回路であってもよい。また信号取得部13をプログラムとし、信号処理部14を集積回路としてもよい。
信号取得部13と信号処理部14とは例えばLAN(Local Area Network)経由で情報の送受信を行うものとしたがこれに限らず、例えばデータバスなどを経由して情報の送受信を行ってもよい。
本実施の形態においては信号ベクトルyを得るために受信アンテナ9~12で観測される受信信号を積分することで、信号ベクトルyの信号雑音比であるSNR(Signal-Noise Ratio)を改善することが可能となる。
なお本実施の形態においては、1MHzの電波を使用する場合について述べたが、これに限らず100kHzの電波を使用してもよい。RMSEなどによってあらわされる誤差は大きくなってしまうものの、1MHzの電波を使用した場合と比べて埋設物3が深く埋設されている場合や海底面2Aからより高い高度から広い範囲を対象とする場合についても検出が可能となる。
まず、100kHzの電波によって、深い位置にある埋設物の探査や広い範囲の探査を実施し、埋設物が検出された場合には、1MHzの電波によって、より推定精度の高い探査を行ってもよい。
(第2の実施の形態)
第1の実施の形態においては、埋設物3が埋設されている位置を推定するまでとしたが、本実施の形態においては埋設物3が埋設されていると推定される位置を表示してもよい。以下では第1の実施の形態との違いに中心をおいて説明を行う。
第1の実施の形態においては、埋設物3が埋設されている位置を推定するまでとしたが、本実施の形態においては埋設物3が埋設されていると推定される位置を表示してもよい。以下では第1の実施の形態との違いに中心をおいて説明を行う。
図4,図5を用いて本実施の形態による埋設物探査装置30について説明する。図4は、本実施の形態による埋設物探査装置30の構成を示すブロック図である。図5は、本実施の形態による埋設物探査装置30が取得する画像35を示す図である。埋設物探査装置30には、画像35を撮像するカメラ33が接続されている。
埋設物探査装置30は、カメラ33が撮像した画像35を取得する画像取得部31と埋設物3が埋設されていると推定される推定位置を表示する推定位置表示部32とを備える。例えば、画像取得部31と推定位置表示部32とは、CPUによって呼び出されるRAMに記録されたプログラムとする。
図6を用いて推定位置表示部32が表示する推定位置について説明する。図6は、本実施の形態による埋設物探査装置30が算出する尤度Lを示す図である。図6に示す尤度マップ36は、画像35におけるそれぞれの個所の尤度分布を示す分布図である。
なお尤度マップ36には、送信アンテナ8と受信アンテナ9~12とがTXANTとANT1~4として表されていてもよい。領域36Aは、ボクセル群15の中で最も尤度Lが高い個所を示す。
推定位置表示部32は、尤度マップ36に基づいて例えば図7に示すような表示画像37を埋設物探査装置30に接続される図示せぬディスプレイなどに表示する。図7は、第2の実施の形態による埋設物探査装置が表示する画像を示す図である。本実施の形態による埋設物探査装置30は、表示画像37を提示することで、埋設物3が埋設されている位置を明確に使用者に提示することができる。
図8を用いて本実施の形態における埋設物探査処理について説明する。図8は本実施の形態による埋設物探査処理の処理手順を示すフローチャートである。ステップS1~S4までは第1の実施の形態と同様であるため説明を省略する。
ステップS1およびステップS4とともに、画像取得部31は、カメラ33から画像35を取得する(S10)。次に推定位置表示部32は、表示画像37を図示せぬディスプレイなどに表示する(S11)。
なおステップS1~ステップS4とステップS10とは、例えば並列して処理されるため、ステップS1~S4,S10の順でなくてもよい。例えば、ステップS1~S3,S10,S4の順でもよい。
なお本実施の形態においては、埋設物探査装置30は、表示画像37における領域36Aのようにボクセル1つ分の領域のみを埋設物3が埋設されているとして位置を表示してもよいがこれに限らない。例えば、図9のように複数のボクセルを表示して尤度Lをたとえば色分けして埋設されている可能性の高低を表示してもよい。
図9は、第2の実施の形態による埋設物探査装置30が表示する他の表示画像38を示す図である。本実施の形態においては、他の表示画像38のように画像35と尤度マップ36とを重ねて表示してもよい。
1,2……損失媒質、3……埋設物、4,30……埋設物探査装置、5……導電率計測計、6……高度計、7……傾斜計、8……送信アンテナ、9~12……受信アンテナ、13……信号取得部、14……信号処理部、15……ボクセル群、33……カメラ、31……画像取得部、32……推定位置表示部。
Claims (5)
- 探査の際にセンサによるセンサ信号を取得し、信号ベクトルとして扱われる電波によって埋設されている埋設物の位置を推定する推定位置を、位置候補ベクトルから選択することで前記埋設物を探査する埋設物探査装置であって、
前記センサ信号及び前記信号ベクトルを取得する信号取得部と、
前記センサ信号に基づいて応答行列を更新して信号処理を行う信号処理部と、を備え、
前記応答行列は、前記位置候補ベクトルと前記信号ベクトルとを関連付ける行列である、埋設物探査装置。 - 埋設物探査装置に接続されたカメラから画像を取得する画像取得部と、
前記埋設物が埋設されていると推定される推定位置を表示する推定位置表示部と、
をさらに備える、請求項1に記載の埋設物探査装置。 - 前記位置候補ベクトルの1次ノルムと前記位置候補ベクトルの2次ノルムとの和に基づく尤度に基づいて推定される位置を1つ表示する、
請求項2に記載の埋設物探査装置。 - 前記推定位置表示部は、前記位置候補ベクトルの1次ノルムと前記位置候補ベクトルの2次ノルムとの和に基づく尤度に基づいて埋設されている可能性の高低を表示する、
請求項2に記載の埋設物探査装置。 - 探査の際にセンサによるセンサ信号を取得し、信号ベクトルとして扱われる電波によって埋設されている埋設物の位置を推定する推定位置を、位置候補ベクトルから選択することで前記埋設物を探査する埋設物探査装置による埋設物探査方法であって、
前記センサ信号及び前記信号ベクトルを取得する第1のステップと、
前記センサ信号に基づいて応答行列を更新して信号処理を行う第2のステップと、を備え、
前記応答行列は、前記位置候補ベクトルと前記信号ベクトルとを関連付ける行列である、埋設物探査方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020123696A JP2022020291A (ja) | 2020-07-20 | 2020-07-20 | 埋設物探査装置及び埋設物探査方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020123696A JP2022020291A (ja) | 2020-07-20 | 2020-07-20 | 埋設物探査装置及び埋設物探査方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022020291A true JP2022020291A (ja) | 2022-02-01 |
Family
ID=80216018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020123696A Pending JP2022020291A (ja) | 2020-07-20 | 2020-07-20 | 埋設物探査装置及び埋設物探査方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2022020291A (ja) |
-
2020
- 2020-07-20 JP JP2020123696A patent/JP2022020291A/ja active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11067685B2 (en) | System, device and methods for imaging of objects using electromagnetic array | |
US9681269B2 (en) | Positioning accuracy using 3D building models | |
US7209752B2 (en) | Error estimate concerning a target device's location operable to move in a wireless environment | |
KR100938047B1 (ko) | 시스템 교정 방법, 위치 탐색 시스템 및 컴퓨터 | |
GB2559157A (en) | Apparatus, method and system for alignment of 3D datasets | |
US20100046327A1 (en) | Methods and Systems to Detect Changes in Multiple-Frequency Band Sonar Data | |
US10094940B2 (en) | Method for localizing a marine mammal in an underwater environment implemented by a PAM system, corresponding device, computer program product and non-transitory computer-readable carrier medium | |
US20210156991A1 (en) | System, device and method for imaging of objects using signal clustering | |
US20160091341A1 (en) | Method and apparatus for object localization | |
JP6396037B2 (ja) | データ解析装置及び方法 | |
WO2009055320A2 (en) | Method and apparatus for determining parameters for a parametric expression characterizing the phase of an acquired signal | |
JP6575002B2 (ja) | 位置推定装置、位置推定方法、及び位置推定プログラム | |
JP5691517B2 (ja) | 位置推定プログラム、位置推定装置、及び位置推定方法 | |
CN114270406A (zh) | 一种参数标定方法、装置及设备 | |
JP2022020291A (ja) | 埋設物探査装置及び埋設物探査方法 | |
JP7130943B2 (ja) | 電波環境推定装置および電波環境推定方法 | |
Rasol et al. | Analysis and calibration of ground penetrating radar shielded antennas | |
US9964640B2 (en) | Method for phase unwrapping using confidence-based rework | |
US20180275299A1 (en) | Acoustic Underground Infrastructure Locating and Mapping System | |
Afkhaminia et al. | Sparsity-based DOA estimation of 2-D rectangular array in the presence of gain and phase uncertainty | |
JP2007163271A (ja) | 地中レーダ画像処理法 | |
JP6819797B2 (ja) | 位置推定装置、位置推定方法とプログラム、並びに、位置推定システム | |
JP7315349B2 (ja) | 位置推定装置、及び位置推定方法 | |
KR101241926B1 (ko) | 위상차 오차 추정에 의한 3차원 인터페로미터 어레이용 레이돔의 적합성 판정방법 | |
KR20210045098A (ko) | Asf 맵 생성 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230607 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240221 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240312 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240425 |