物標探査装置、物標探査プログラム及び物標探査方法
【課題】音源から直接進行する進行波に混合した前方散乱波を用いて物標を正確に探査する。
【解決手段】伝搬空間内に音波を送波する音源と、前記伝搬空間内に存在する物標から前方に散乱する前方散乱波を受波する領域に配置されたトランスデューサアレイと、前記トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する減算処理手段と、前記減算処理手段で分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成するパッシブ位相共役処理手段と、前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する自己相関処理手段と、前記自己相関処理信号と前記共役性相関信号との相似性を判定する相関手段とを含む。
【解決手段】伝搬空間内に音波を送波する音源と、前記伝搬空間内に存在する物標から前方に散乱する前方散乱波を受波する領域に配置されたトランスデューサアレイと、前記トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する減算処理手段と、前記減算処理手段で分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成するパッシブ位相共役処理手段と、前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する自己相関処理手段と、前記自己相関処理信号と前記共役性相関信号との相似性を判定する相関手段とを含む。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、パッシブ位相共役処理を用いて物標を探査する物標探査装置、物標探査プログラム及び物標探査方法に関する。
【背景技術】
【０００２】
海中音響の分野において、信号を時間反転させる研究が増加している（非特許文献１〜３）。浅い海域において物標を探査する際、海中に送波した音波が海面や海底で反射し、これらの反射波が存在するために、物標の探査が困難になる場合がある。
【０００３】
近年、物標からの反射波に位相共役処理を施す探査法が注目されている（非特許文献４〜８）。これらの探査法は、物標から後方に散乱する後方散乱波を用いる。ここで、音波が物標に進入した場合について考える。音波が物標に進入すると、一部の音波が物標で散乱し、散乱波として音源側に逆行する。この散乱波を後方散乱波という。また、一部の音波が物標から前方に散乱する。この散乱波を前方散乱波という。
【非特許文献１】S.Kim, W. A. Kuperman, W. S. Hodgkiss, H. C. Song, G. F. Edelmann, and T. Akal, “Robusttime reversal focusing in the ocean,” J. Acoust. Soc. Am. 114, 145-157(2003)
【非特許文献２】S.C.Walker, P. Roux, and W.A. kuperman, “Focal depth shifting of a timereversal mirror in a rang-independent waveguide,” J. Acoust. Soc. Am. 118,1341-1347(2005)
【非特許文献３】S.C. Walker, W. A. Kuperman, and P. roux, “Active waveguide Green’s functionestimation with application to time-reversal focusing without a probe source ina rang-indepent waveguide,” J. Acoust. Soc. Am. 120, 2755-2763(2006)
【非特許文献４】G.Montaldo, M. Tanter, and M. Fink,” Revisiting iterative time reversalprocessing: Application to detection of multiple targets” J. Acoust. Soc. Am.115, 776-784(2004)
【非特許文献５】C.Prada, S. Manneville, D. Spoliansky, and M. Fink, “Decomposition of the timereversal operator: Detection and selective focusing on two scatterers,” J.Acoust. Soc. Am. 99, 2067-2076(1996)
【非特許文献６】N.Mordant, C. Prada, and M. Fink, “ Highly resolved detection and selectivefocusing in a waveguide using the D. O. R. T. method,” J. Acoust. Soc. Am. 105,2634-2642(1999)
【非特許文献７】D.H. chamberes, “Analysis of the time-reversal operator for scatterers of finitesize,” J. Acoust. Soc. Am. 112, 411-419(2002)
【非特許文献８】D.H. chamberes, A.K. Gautesen, “Time reversal for a single spherical scatterer,”J. Acoust. Soc. Am. 109, 2616-2624(2001)
【発明の開示】
【発明が解決しようとする課題】
【０００４】
非特許文献４〜８に開示された探査法では、物標から後方に散乱する後方散乱波を用いているが、小さい物標の場合や低周波の音波の場合、後方散乱波を受波する際、十分な後方散乱波を得られない場合がある。したがって、前記後方散乱波を用いた探査法では、物標を正確に探査できなくなる可能性が残されていた。
【０００５】
また、小さい物標の場合、一般的に前方散乱波のレベルが後方散乱波のレベルよりも高い。前記前方散乱波を受波する際、前方散乱波以外にも、音源から受波側に進行する進行波も、物標から離れた領域に存在する。前方散乱波は、音源からの音波が物標に進入して散乱したものであるから、前方散乱波が音源から直接進行する進行波に隠れてしまい、前方散乱波を用いて物標を正確に探査することができない場合がある。
【０００６】
本発明の目的は、音源から直接進行する進行波に混合した前方散乱波を用いて物標を正確に探査する物標探査装置、物標探査プログラム及び物標探査方法を提供することにある。
【課題を解決するための手段】
【０００７】
前記目的を達成するため、本発明に係る物標探査装置は、伝搬空間内に存在する物標を探査する物標探査装置であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波する音源と、
前記伝搬空間内に存在する物標から前方に散乱する前方散乱波を受波する領域に配置されたトランスデューサアレイと、
前記トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する減算処理手段と、
前記減算処理手段で分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成するパッシブ位相共役処理手段と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する自己相関処理手段と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する相関手段とを含むことを特徴とする。
【０００８】
以上の例では、本発明をハードウェアとしての物標探査装置として構築したが、これに限られるものではない。本発明は、物標を探索する方法として構築する、或いは物標探索装置の機能をコンピュータに実行させるソフトウェアとしてのプログラムとして構築してもよいものである。
【０００９】
本発明に係る物標探査プログラムは、伝搬空間内に存在する物標を探査する制御を行う物標探査プログラムであって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標の探査を制御するものであり、
コンピュータに、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する機能と、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成する機能と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する機能と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する機能とを実行させることを特徴とする。
【００１０】
本発明に係る物標探索方法は、伝搬空間内に存在する物標を探査する物標探査方法であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波し、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離し、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成し、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成し、
前記自己相関処理信号と前記共役性相関信号との相似性を判定することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、伝搬空間内に存在する物標に進入する音源からの音波の信号に対する自己相関関数と、物標から前方に散乱する前方散乱波に対するパッシブ位相共役結果とが一致することを利用することにより、伝搬空間内に存在する物標を正確に探査することができる。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の実施形態を図に基づいて説明する。
【００１３】
本発明の実施形態に係る物標探査装置は伝搬空間１内に存在する物標を探査する物標探査装置であって、図２に示すように、伝搬空間１内に存在する物標２に進入する音源３からの進行波４と、前記物標２から前方に散乱する前方散乱波５とを用いることにより、前記物標２を探査するものである。伝搬空間１として例えば海中を想定した場合、図２の１ａが海面であり、１ｂが海底であり、１が伝搬空間としての海中である。なお、音源３は海中１の定位置に設置する、或いは探査船から海中に吊下して設置するが、音源３の設置については特に限定されるものではない。図２において、音源３からの音波が物標２に入射した際に、音源３に向かう後方に生じる散乱波が後方散乱波５ａであり、トランスデューサアレイ６に向かう前方に生じる散乱波が前方散乱波５である。
【００１４】
本発明の実施形態は、物標２の位置に入射した音波の「自己相関関数」と、物標２から前方に散乱する前方散乱波５をトランスデューサアレイ６で受波してパッシブ位相共役処理を施した波形とが相似することを証明し、その証明結果に基づいて物標の探査を実現したことを特徴とするものである。
【００１５】
物標２の位置に入射した音波の「自己相関関数」と、物標２から前方に散乱する前方散乱波５をトランスデューサアレイ６で受波してパッシブ位相共役処理を施した波形とが相似することを検証する。
【００１６】
最初に位相共役について検証する。図２に示すように、例えば浅い海域において音源３とトランスデューサアレイ６とを対峙して設置する。音源３から送波される音波（進行波４）は、離れた地点に配置されたトランスデューサアレイ６により受波される。そして、トランスデューサアレイ６で受波した受波信号に対してパッシブ位相共役処理を行い、再びトランスデューサアレイ６から送波された場合、受波点における音波は、以下の式によって表される。
【数１】
式（１）は、文献 D. R. Jackson and D.R. Dowling, “Phase conjugation in underwater acoustics” J. Acoust. Soc. Am. 89, 171-181 (1991) で確立されている。
式（１）において、
【数２】
は音源３からトランスデューサアレイ６の各素子６ａへの伝搬に関するグリーン関数、
【数３】
はトランスデューサアレイ６の各素子６ａから受波点への伝搬に関するグリーン関数である。ｒは距離を、ｓ、ｎ、* は、それぞれ、音源、トランスデューサアレイの各素子の番号，複素共役を示している。
【００１７】
式（１）は周波数領域におけるものだが、以下の式では時間領域におけるものを示す。
【数４】
式（２）は、W. A. Kuperman, W. S. Hodgkiss, C. Song, T. Akal, C. Ferla and D. R. Jackson “Phase conjugation in the ocean; Experimental demonstration of an acoustic time-reversal mirror” J. Acoust. Soc. Am. 103, 25-40(1998) で確立されている。
式（２）において、Gω，ｔ，ｚはそれぞれ、グリーン関数，時間，深度である。Ｒは音源３からトランスデューサアレイ６までの水平距離、ｒはトランスデューサアレイ６から受波点との間の距離を示す。Ｓ（ω）は音源３から送波される音波の周波数スペクトル、ωは角周波数である。式（２）はタイムリバーサル波の音源近くの集束特性を調べるときにしばしば用いられる。
【００１８】
本発明者は、前記式（１），（２）が音源３の位置以外の受波点の音場にも適用可能であると考えている。このことについて考察する。すなわち、音源３の地点以外の任意の点の音場の位相共役について考察する。ここで、任意の点は、図２に示すように音源３とトランスデューサアレイ６との間にあるものと仮定する。
【００１９】
前記任意の点が新しい音源の位置であるものと考えられる。そして、新しい音源の駆動信号が元の音源３から供給される。したがって、元の音源３の周波数スペクトルＳ（ω）と、新しい音源の周波数スペクトルＳＲ（ω）の関係は以下のように表される。
【数５】
式（３）において、添字Ｒは任意の点をあらわす。
【００２０】
次に、式（２），（３）を用いて、音源３以外の任意の点の位相共役について考察する。一般的な位相共役波の場合、音源３から送波される音波パルスは図２に示す処理のため、トランスデューサアレイ６で受波され、その各素子６ａ毎に時間反転処理（タイムリバース処理）が行われる。その後、時間反転された信号（時間反転信号）が再び各素子６ａから送波され、音源３から送波される音波パルスと同様の音波パルスが音源３の位置に形成される。
【００２１】
本発明の実施形態において、上述した一般的な位相共役波と同様の特性が任意の地点で維持されているか否かを検証する。すなわち、音源３とトランスデューサアレイ６との間の任意の地点において、音源３から直接進行する音波（以下、直接波という）と時間反転信号との相似性を検証する。
【００２２】
シミュレーションの結果について考察する。先ず、任意の点における直接波と時間反転波の比較について考察する。
【００２３】
図２に示す処理において、海底１ｂまでの深さは１００ｍに設定し、音源３とトランスデューサアレイ６との間の距離は２．５ｋｍに設定してある。海水中の音速及び海水の密度はそれぞれ１５００ｍ／ｓ及び１０００ｋｇ／ｍ３である。そして、海底１ｂの堆積物における音速及び堆積物の密度はそれぞれ１６００ｍ／ｓ及び１５００ｋｇ／ｍ３である。
【００２４】
中心周波数５００Hz、パルス幅８サイクルのトーンバースト波が、深度５０ｍに設置された音源３から送波される。トーンバースト波の波形を図５（ａ）、そのスペクトルを図５（ｂ）にそれぞれ示す。
【００２５】
音波パルスの波形を詳しく調べるシミュレーションにおいてはワイドバンドパルスが用いられる。これらの音波パルスはトランスデューサアレイ６で受波され、その各素子６ａ毎に時間反転処理が施される。そして、前記時間反転処理した信号が各素子６ａから再送波される。再送波された音波は音源３とトランスデューサアレイ６との間（図２の任意の点）で受波される。
【００２６】
そして、時間反転信号と受波点で形成された駆動信号（音源３から供給される音波パルス）との相似性について検証する。音源３から１ｋｍの地点において、音源３から任意の点に直接進入する音波と、時間反転信号とを比較する。この場合、元の音源３から伝わる音波パルスはトーンバースト波である。図５（ｂ）に示すスペクトルは式（２）のＳ（ω）に相当する。
【００２７】
計算に用いる周波数スペクトルのバンド幅は３５０Ｈｚである。その理由は、周波数スペクトルのハンド幅を３５０Ｈｚに設定することが、任意の点における波形を詳細に表すからである。このバンド幅は、図５（ｂ）に示す周波数スペクトルの第３ローブまでのバンド幅に相当する。一方、新しい音源の周波数スペクトルは、元の音源３から任意の位置に直接進入する音波パルスのフーリェ変換により、容易に得られる。
【００２８】
新しい音源の周波数スペクトルＳR（ω）は、元の音源３の周波数スペクトルＳ（ω）と、式（３）により関連付けられる。新しい音源の位置は、元の音源３に対して距離１ｋｍ、深度５０ｍの地点である。元の音源３から新しい音源の位置に進入する音波パルスの波形を図６（ａ）、時間反転信号の波形を図６（ｂ）にそれぞれ示す。
【００２９】
図６（ｂ）に示す時間反転信号の進行方向は、図６（ａ）に示す直接波の進行方向と反対であるため、時間反転信号は時間を反転した座標で表してある。さらに、式（２）に示す周波数スペクトルＳ（ω）は、受波点に進入する音波を新しい音源と考えることによって、元の音源３からの音波パルスのフーリェ変換により直接求められる。図６（ａ），（ｂ）を比較すると、元の音源３から直接進入する音波の波形は時間反転信号の波形と一致する。
【００３０】
上述した検証結果から、位相共役性が音源３とトランスデューサアレイ６との間の任意の点においても維持されていることが分かる。
【００３１】
次に、音源３とトランスデューサアレイ６とを浅い海域に配置する。プローブ信号Ｐiとデータ信号Ｐdが音源３から送られ、音源３からの音波がトランスデューサアレイ６のｍ番目の素子で受波された場合（位置ｒm）、プローブ信号Ｐiとデータ信号Ｐdとの信号の相互相関関数は以下の式（４）によって表される。
【数６】
【００３２】
次に、トランスデューサアレイ６のすべての素子６ａの相互相関関数は以下の式（５）によって加えられる。
【数７】
Ｓ（ｔ）はプローブ信号の位相共役になることが知られている。この技術はパッシブ位相共役と呼ばれる。
なお、式（４），（５）は、D. Rouseff, D. R. Jackson, W. L. J. Fox, C. D. Jones, J. Ritcey and D. R. Dowling, “Underwater Acoustic Communication by Passive-Phase Conjugation; Theory and Experimental Results” IEEE J. Ocean, Eng. 26, 821-831 (2001) で確立されている。
【００３３】
音源３から進行する音波に対するパッシブ位相共役処理について考察する。図２のように、音源３とトランスデューサアレイ６とを浅い海域に配置する。音源３からトランスデューサアレイ６までの距離は３ｋｍである。物標２は音源３から距離１．５ｋｍ、深度５０ｍの地点に配置する。物標２の材料は鉄であり、高さと幅がともに２mの大きさのものである。
【００３４】
中心周波数５００Ｈｚ、パルス幅８サイクルのトーンバースト波が音源３から送波され、音波パルスが深度５０ｍに設置された音源３から送波される。音源３から送波された音波パルスは物標２に入射し、物標２から前方に散乱する前方散乱波５が発生する。
【００３５】
音源３から音波を送波すると、音源３からの進行波４と前方散乱波５とがトランスデューサアレイ６にほぼ同時に到達する。そして、トランスデューサアレイ６で受信した信号にパッシブ位相共役処理を施す。すなわち、トランスデューサアレイ６の各素子６ａで受波した信号に式（４）による相関処理を行う。さらに、トランスデューサアレイ６のすべての素子６ａの相互相関関数（式（５））を加える。この一連の相関処理がパッシブ位相共役処理である。
【００３６】
図７は、上述したパッシブ位相共役処理を経て得た進行波４と前方散乱波５との混合波のパッシブ位相共役信号を表している。図７から明らかなように、物標２による前方散乱波の影響は見られず、波形は図５（ａ）に示す音波パルスの自己相関関数の波形とほぼ一致する。パッシブ位相共役処理は、式（４）におけるプローブ信号Ｐｉとデータ信号Ｐｄとが同じであるから、混合波の自己相関関数となる。
【００３７】
一般に、音源３から送波される進行波４のレベルは前方散乱波５のレベルより大きく、前方散乱波５を検出するのは困難である。本発明の実施形態では、物標２が存在しないときに、予め音源３から音波（進行波４）を送波し、トランスデューサアレイ６の各素子６ａで受波した信号を保持しておく。その後、物標２が存在するときに音源３から音波（進行波４）を送波し、トランスデューサアレイ６の各素子６ａにより進行波４と前方散乱波５との混合波を受波し、混合波から進行波４を減じることにより、前方散乱波５を分離して顕在させる。
【００３８】
次に、物標２から前方に散乱した前方散乱波５に対するパッシブ位相共役処理について検証する。
【００３９】
本発明の実施形態では、物標２が存在するときにトランスデューサアレイで受波した混合波の信号から、物標２が存在しないときにトランスデューサアレイで受波した進行波４の信号を減ずる。そして、上述したパッシブ位相共役処理を、減じて顕在させた前方散乱波５に対して施す。前記減算した信号（前方散乱波５）に対して施すパッシブ位相共役処理は、式（４）におけるプローブ信号Ｐｉとデータ信号Ｐｄとが同じであるから、減算信号の自己相関関数となる。
【００４０】
物標２は距離１．５ｋｍ、深度５０mの地点に配置され、音波パルスが深度５０ｍに設置された音源３から送波される。トランスデューサアレイ６で受波した混合波から進行波４の成分を取り除いた後、パッシブ位相共役処理を施す。そのパッシブ位相共役処理を行った結果の波形を図８に示す。図８に示す波形は、音波パルスの繰り返し周期である１０秒を中心とした対称的な波形である。
【００４１】
以上のようにして得られた図８に示すパッシブ位相共役処理による波形が、物標２の位置を新しい音源の位置としてその音源から前方に放射された音波パルスにパッシブ位相共役処理を施した波形であることを確かめる必要がある。これを確かめるために、物標２が新しい音源の位置（任意の位置）にあることを確認する。この場合、前記新しい音源から放射される音波パルスは、トーンバースト波ではない。前記新しい音源が音波パルスを放射する駆動信号は、元の音源３から前記新しい音源の位置に直接進入する音波パルスである。前記新しい音源から放射されるパルスがトランスデューサ６で受波され、その受波信号にパッシブ位相共役処理を施す。
【００４２】
図９は、前記新しい音源から放射された信号にパッシブ位相共役処理を施した信号の波形を示している。図９に示す波形は図８に示す波形と実質的に一致している。図９に示す信号は、トランスデューサアレイ６で受波した混合波から進行波４を減算処理して得た物標２からの前方散乱波５のみであることが確認できた。
【００４３】
次に、物標２の位置とパッシブ位相共役処理との関係について考察する。伝搬空間１における伝搬環境は、これまで検証した際に記述したものと同じである。物標２の位置のみを変更しており、位置を変更した物標２から前方に散乱する前方散乱波５をトランスデューサアレイ６で受波し、その受波信号にパッシブ位相共役処理を施す。
【００４４】
図１０(ａ)，（ｂ），（ｃ）は、物標２の距離に関するパラメータとしてのパッシブ位相共役処理による波形を示している。物標２の距離（音源３に対する位置）は、１．３ｋｍ，１．６ｋｍ，１．７ｋｍである。また、物標２の深度はすべて５０ｍである。図１０から明らかなように、前記パッシブ位相共役処理による波形の変化は、物標２の距離によるところが大きいことがわかる。同様に、物標２の深度が変わると、前記パッシブ位相共役処理による波形は大きく変化する。図１１は、物標２の距離が１．５ｋｍであって、物標２の深度が２０ｍに変化した場合における前記パッシブ位相共役処理の波形を示している。図１０及び図１１からして、物標２の位置と前記パッシブ位相処理による波形とは、密接に関係していると考えられる。
【００４５】
以上の検証結果から、物標２から前方に散乱する前方散乱波５に対するパッシブ位相共役処理の結果は、音源３から送波されて物標２に進入する音波パルスと直接関連する。したがって、物標２から前方に散乱する前方散乱波５に対してパッシブ位相共役処理を施して得られる波形は、音源３から送波されて物標２の位置に進入する音波パルスの相関波形（自己相関関数）と同様である。一方、音源３を基準とする音場の任意の点に進入する音波パルスは伝搬環境により一意に決められる。すなわち、物標２から前方に散乱する前方散乱波５に対するパッシブ位相共役処理の結果は、物標２の位置に関する情報を保持する。
【００４６】
以上の関係はマッチドフィールド法として物標２の探査に用いることが可能である。物標２の音源３に対する距離が１．５ｋｍであり、物標２の深度が５０mにある場合、音源３を基準とする任意の点に進入する音波パルスの自己相関関数と、任意の点から前方に放射される音波パルスに施したパッシブ位相共役処理の結果との相関関数とは図１１に示すようになる。
【００４７】
図１２から明らかなように、前方散乱波５にパッシブ位相共役処理を行って生成した信号（前方散乱波の共役性相関信号）と、音源３を基準位置として伝搬空間内の任意の位置に存在する物標２の位置に入射する進行波４に自己相関処理を行って生成した信号（進行波の自己相関処理信号）との相関により得た相関値の最も高い位置Ｈが物標２の位置であることを明確に示している。
【００４８】
以上の検証結果を考察すると、音源３を基準とする伝搬空間の任意の位置から放射される音波パルスにパッシブ位相共役処理を施した結果と、前記任意の位置に音源３から進入する音波パルスの自己相関処理の結果とが一致することが証明できた。前記任意の位置に物標２が存在すれば、その前方散乱波は物標２に進入した音波パルスに対応する。したがって、物標２の位置に音源３から進入する進行波の自己相関処理の結果（進行波の自己相関処理信号）と、物標２から前方に散乱する前方散乱波５の信号にパッシブ位相共役処理を施した結果（前方散乱波の共役性相関信号）とが一致することとなり、進行波の自己相関処理信号と前方散乱波５の共役性相関処理信号との相関値が高く、任意の位置に物標２が存在しない場合の相関値が低くなる。
【００４９】
本発明の実施形態は、音源３から送波されて音源３を基準とした伝搬空間１の任意の位置に入射する音波（進行波４）のパッシブ位相共役処理の結果である共役性相関信号Ｓ２と、音源３から送波されて前記任意の位置と同一位置に入射する音波の自己相関処理の結果である進行波４の自己相関処理信号Ｓ１とが一致するか否かを監視することにより、音波が伝搬する伝搬空間１、例えば浅い海域における物標２を探査することを特徴とするものである。任意の位置に物標２が存在する場合に共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相関値が最も高くなるためである。以下、具体的に説明する。
【００５０】
本発明の実施形態に係る物標探査装置は図１に示すように、基本的な構成として、伝搬空間１内に存在する物標２に進入する音源からの進行波４と、前記物標２から前方に散乱する前方散乱波５とを用いることにより、前記物標２を探査するものである。そして、伝搬空間１内に音波を送波する音源３（図２参照）と、伝搬空間１内に存在する物標２から前方に散乱する前方散乱波５を受波する領域に配置されたトランスデューサアレイ６（図２参照）と、トランスデューサアレイ６に直接進行する進行波４と前方散乱波５との混合波から、トランスデューサアレイ６に直接進行した進行波４を減じることにより、前方散乱波５を分離する減算処理手段７と、減算処理手段７で分離した前方散乱波５にパッシブ位相共役処理を行うことにより、前方散乱波５の共役性相関信号Ｓ２を生成するパッシブ位相共役処理手段８と、音源３を基準位置として伝搬空間１内の任意の位置に入射する進行波４に自己相関処理を施すことにより、進行波４の自己相関処理信号Ｓ１を生成する自己相関処理手段９と、自己相関処理信号Ｓ１と共役性相関信号Ｓ２との相似性を判定する相関手段１０とを含むことを特徴とする。
【００５１】
本発明の実施形態では、伝搬空間１内に音波（進行波４）を音源３から送波し、トランスデューサアレイ６に直接進行する進行波４と前方散乱波５との混合波から、トランスデューサアレイ６に直接進行した進行波４を減じることにより、前方散乱波５を分離する。前方散乱波５を分離するにあたっては、物標２が存在しない伝搬空間１内を伝搬する進行波４をトランスデューサアレイ６で受波して進行波４の信号を予め取得しておき、次に前記混合波をトランスデューサアレイ６で受波し、前記混合波の信号から、トランスデューサアレイ６に直接進行した進行波４の信号を減じることにより、前方散乱波５を分離する。
【００５２】
次に、前記分離した前方散乱波５にパッシブ位相共役処理をパッシブ位相共役処理手段８で行うことにより、前方散乱波５の共役性相関信号Ｓ２を生成する。また、前記音源３を基準位置として伝搬空間１内の任意の位置に入射する進行波４に自己相関処理を自己相関処理手段９で施すことにより、進行波４の自己相関処理信号Ｓ１を生成する。そして、自己相関処理信号Ｓ１と共役性相関信号Ｓ２との相似性を相関手段１０で判定する。
【００５３】
なお、以上の説明では、本発明の実施形態に係る物標探査装置をハードウェアとして構築する場合について説明したが、これに限られるものではない。すなわち、本発明の実施形態に係る物標探査装置が実行する機能をプログラムとして構築し、そのプログラムをコンピュータに実行させることにより、本発明の実施形態に係る物標探査装置が実行する機能を実現するようにしてもよいものである。この場合、前記プログラムは記録媒体に記録され、その状態で商取引の対象となる。
【００５４】
本発明の実施形態に係る物標探査装置をさらに具体例を用いて詳細に説明する。
【００５５】
自己相関処理手段９は、環境情報生成部１１と、波形生成処理部１２と、メッシュ処理部１３と、自己相関処理部１４とを含んでいる。
【００５６】
環境情報生成部１１は、音源３から送波されてトランスデューサアレイ６で受波変換された校正信号を取り込み、不要な信号を除去すると共にＡ／Ｄ変換し、これを校正信号として用いる。さらに、環境情報生成部１１は、前記校正信号に基づいて実海域の定数、すなわち、音波ビームが伝搬する実際の海域の定数（海洋の音響定数）を演算する。具体的に説明する。
【００５７】
本発明の実施形態に係る物標探査装置は、例えば大陸棚のような比較的浅い海域で用いられる。その環境における音波伝搬に関する海洋の音響定数は、水深、海水の音速、海底堆積物（音速，密度）等のデータである。また、前記環境における音波パルスに関する海洋の音響定数は、音波パルスの中心周波数、音波パルスのスペクトル及び音波パルスのバンド幅等のデータである。海洋の音響定数を決定するための水深のデータは、測深機などの測定機器を用いることで直接計測により求められる。海洋の音響定数を決定するための海水の音速データは、水温の計測により、例えばメドイン（H.Medwin）の式などの既知の音速の式から求められる。
【００５８】
海洋の音響定数を決定するための海底堆積物のデータは、海図などから大凡の値を推定できるが、より正確な値は、既知の校正法により求められる。すなわち、物標が存在していると予想される探知海域に吊下した校正音源から送波された信号をトランスデューサアレイ６で受波し、その受波信号と予測した海底堆積物の定数とを変化させながら演算することで、実海域での堆積物のデータを取得する。ここに、海洋の音響定数とは、実際に音波が海中を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる環境因子を意味する。なお、海中での物標２を探索する場合を想定しているので、海洋音響定数を用いたが、これに限られるものではない。伝搬空間として海洋以外のものを用いる場合、前記音響定数は、実際に音波が伝搬空間を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる環境因子を意味することとなる。
【００５９】
環境情報生成部１１には、上述した音波伝搬及び音波パルスに関する実海域での海洋の音響定数を決定するためのデータが入力され、環境情報生成部１１は、前記入力された海洋の音響定数を決定するためのデータと、前記校正信号とに基づいて実海域の海洋の音響定数を決定する。環境情報生成部１１は、決定した実海域での海洋の音響定数を波形生成処理部１２に向けて出力する。なお、海中での物標２を探索する場合を想定しているので、海洋音響定数を用いたが、これに限られるものではない。伝搬空間として海洋以外のものを用いる場合、環境情報生成部１１は、実際に音波が伝搬空間を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる実際の音響定数を出力することとなる。
【００６０】
メッシュ処理部１３は図４に示すように、音源３とトランスデューサアレイ６との間の伝搬空間１内に物標２の位置を仮定し、その仮定した任意の位置情報を波形生成処理部１２に出力する。具体的に説明すると、伝搬空間１が例えば海中であるとすると、メッシュ処理部１３は図４に示すように、音源３とトランスデューサアレイ６との間の距離をｒとして任意の位置（距離）を水平座標上のｘ−ｙ面に設定する。さらに、メッシュ処理部１３は図４に示すように、海面１ａから海底１ｂまでの深度をｚとして任意の位置（深度）を垂直座標上のｙ−ｚ面に設定する。
【００６１】
そして、メッシュ処理部１３は、音源３に対する距離ｒを変更させて任意の位置をｘ−ｙ面内に順次移動させるとともに、音源３に対する深度ｚを変更させて任意の位置をｙ−ｚ面内に順次移動させ、その任意の位置情報（距離，深度）を波形生成処理部１２に出力する。
【００６２】
波形生成処理部１２は、環境情報生成部１１から出力される音響定数とメッシュ処理部１３から出力される位置情報とに基づいて、音源３から送波された音波が入射する点における新しい音源、すなわち音源３を基準とした音場における任意の位置における音波パルスを式（６）及び式（７）に基づいて求める。
【数８】
式（６）で求める音波パルスは音源３から任意の位置まで通常に伝搬した音波パルスである。式（６）において、Ｓ（ω）は音源３から送波される音波パルスの周波数スペクトルである。Ｐ（ｒ，ｚ，ω）は前記周波数スペクトルＳ（ω）に対応した任意の点における音圧であり、通常のノーマルモード法等により求める。ωｂはバンド幅である。
【００６３】
波形生成処理部１２は、例えば結合モード法を用いることにより、角周波数ωの音波の任意の点（ｒ，ｚ）における音圧Ｐ（ｒ，ｚ，ω）を式（７）に基づいて求め、これを式（６）に代入する。
【数９】
式（７）において、ｊは水平方向の区分領域（距離）を意味している。φ（ｚ，λｊ，ｍ）はモード関数である。
式（７）は、R. B. Evans, “A coupled mode solution for acoustic propagation in a wave guide with stepwise depth variations of a penetrable bottom.“ J. Acoust. Soc. Am., 74, 188-193(1983) によって確立されている。
【００６４】
波形生成処理部１２は、式（６），式（７）に基づいて計算した音源３を基準とする伝搬空間１内の任意の位置での音波（パルス）を自己相関処理部１４に出力する。
【００６５】
自己相関処理部１４は、波形生成処理１２から受け取った音波情報に基づいて、その音波信号に式（４）に基づく相関処理を行う。この場合、式（４）におけるプローブ信号Ｐｉとデータ信号Ｐｄとは同じであり、任意の位置に入射する進行波４を相関処理することは、自己相関処理の形態となる。したがって、自己相関処理部１４は、音源３を基準とする伝搬空間１内での任意の位置に入射する進行波４に自己相関処理を施すことにより、進行波４の自己相関処理信号Ｓ１を生成する。
【数１０】
【００６６】
自己相関処理部１４が出力する自己相関処理信号Ｓ１は、物標２が前記任意の位置に存在する場合、図８に示す波形の信号となる。
【００６７】
音源３から送波されて伝搬空間１内を伝搬した音波パルスがトランスデューサアレイ６で受波され、トランスデューサアレイ６は、受波音波を電気信号である受波信号に変換し、その受波信号を減算処理手段７に出力する。
【００６８】
減算処理手段７は、予め伝搬空間１内に物標２が存在しない場合に音源３からトランスデューサアレイ６に直接進行する進行波４の受波信号を記憶しておく。伝搬空間１内に物標２が存在する場合に前方散乱波５を含む進行波４がトランスデューサアレイ６に受波された場合、減算処理手段７は、進行波４と前方散乱波５との混合波から進行波４を減算することにより、前方散乱波５を分離し、その分離した前方散乱波５の受波信号をパッシブ位相共役処理手段８に出力する。
【００６９】
パッシブ位相共役処理手段８は、減算処理手段７で減算処理された減算信号を受け取ると、前記減算信号に式（４）及び式（５）に基づくパッシブ位相共役処理を施すことにより、前方散乱波５の共役性相関信号Ｓ２を生成する。パッシブ位相共役処理手段８は、前方散乱波５の共役性相関信号Ｓ２を相関手段１０に出力する。
【数１１】
【数１２】
【００７０】
パッシブ位相共役処理手段８は、前記任意の位置に物標２が存在する場合、図９に示すような波形の共役性相関信号Ｓ２を出力する。プローブ信号Ｐｉとデータ信号Ｐｄとは同一であるから、パッシブ位相共役処理手段８によるパッシブ位相共役処理は自己相関処理の形態となる。
【００７１】
相関手段１０は、相関処理部１６と、相関波形解析部１７とを含んでいる。
【００７２】
相関処理部１６は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、通常の相互相関式（８）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を判定する。相関処理部１６は、前記自己相関処理信号Ｓ１と前記共役性相関信号Ｓ２との相関で得た相関値が最も高い位置を物標位置として判定する。
【数１３】
式（８）において、Ｓ１（ｎ）は自己相関処理信号Ｓ１のデータ、Ｓ２（ｎ＋ｋ）は鏡や規制相関信号Ｓ２のデータである。
式（８）は、例えば、「ディジタル信号処理 江原義郎 著 P46-P47 東京電気大学出版局 １９９６年２月２０日発行」などに掲載されている相互相関関数である。
【００７３】
相関波形解析部１７は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、式（４）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を解析し、画像表示する。相関波形解析部１７が画像表示する一例を図１２に示す。
【００７４】
制御部１８は、自己相関処理手段９，パッシブ位相共役処理手段８，相関手段１０及びその他の構成要素を総合的に制御するものである。２０は、制御部１８の動作に必要な情報、及び演算に必要な作業領域を提供する記憶部、１９は、制御部１８を通して出力されるデータを外部に出力する指示部である。
【００７５】
次に、本発明の実施形態に係る物標探査装置を用いて伝搬空間、例えば海中１内に存在する物標２を探査する場合について説明する。
【００７６】
環境情報生成部１１は、音源３から送波されてトランスデューサアレイ６で受波変換された校正信号を取り込み、不要な信号を除去すると共にＡ／Ｄ変換し、これを校正信号として用い、前記校正信号に基づいて実海域の定数、すなわち、音波ビームが伝搬する実際の海域の定数（海洋の音響定数）を演算し、実際に音波が伝搬空間を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる実際の音響定数を波形生成処理部１２に出力する。
【００７７】
一方、メッシュ処理部１３は、音源３に対する距離ｒを変更させて任意の位置をｘ−ｙ面内に順次移動させるとともに、音源３に対する深度ｚを変更させて任意の位置をｙ−ｚ面内に順次移動させ、その任意の位置情報（距離，深度）を波形生成処理部１２に出力する。
【００７８】
波形生成処理部１２は、環境情報生成部１１から出力される音響定数とメッシュ処理部１３から出力される位置情報とに基づいて、音源３から送波された音波が入射する点における新しい音源、すなわち音源３を基準とした音場における任意の位置における音波パルスを式（６）及び式（７）に基づいて求め、伝搬空間１内の任意の位置での音波パルスを自己相関処理部１６に出力する。
【００７９】
自己相関処理部１４は、波形生成処理部１２から受け取った音波情報に基づいて、その音波信号に式（４）に基づく相関処理を行い、音源３を基準とする伝搬空間１内での任意の位置に入射する進行波４に自己相関処理を施すことにより、進行波４の自己相関処理信号Ｓ１を生成する（図３のステップ１００）。
【００８０】
自己相関処理部１４が出力する自己相関処理信号Ｓ１は、物標２が前記任意の位置に存在する場合、図８に示す波形の信号となる。
【００８１】
自己相関処理信号Ｓ１を生成する処理を行う際に、共役性相関信号Ｓ２を生成する処理を音源３からトーンバースト波を出力するタイミングで行う。
【００８２】
すなわち、音源３から送波されて伝搬空間１内を伝搬した音波パルスをトランスデューサアレイ６で受波し、トランスデューサアレイ６は、受波音波を電気信号である受波信号に変換し、その受波信号を減算処理手段７に出力する。
【００８３】
減算処理手段７は、予め伝搬空間１内に物標２が存在しない場合に音源３からトランスデューサアレイ６に直接進行する進行波４の受波信号を記憶しておく。伝搬空間１内に物標２が存在する場合に前方散乱波５を含む進行波４がトランスデューサアレイ６に受波された場合、減算処理手段７は、進行波４と前方散乱波５との混合波から進行波４を減算することにより、前方散乱波５を分離し、その分離した前方散乱波５の受波信号をパッシブ位相共役処理手段８に出力する。
【００８４】
パッシブ位相共役処理手段８は、減算処理手段７で減算処理された減算信号を受け取ると、前記減算信号に式（４）及び式（５）に基づくパッシブ位相共役処理を施すことにより、前方散乱波５の共役性相関信号Ｓ２を生成する。パッシブ位相共役処理手段８は、前方散乱波５の共役性相関信号Ｓ２を相関手段１０に出力する（図３のステップ２００）。
【００８５】
パッシブ位相共役処理手段８は、前記任意の位置に物標２が存在する場合、図９に示すような波形の共役性相関信号Ｓ２を出力する。プローブ信号Ｐｉとデータ信号Ｐｄとは同一であるから、パッシブ位相共役処理手段８によるパッシブ位相共役処理は自己相関処理の形態となる。
【００８６】
相関処理部１６は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、式（８）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を判定する（図３のステップ３００）。相関処理部１６が、自己相関処理信号Ｓ１と共役性相関信号Ｓ２との相関で得た相関値が最も高い位置を物標位置として判定する。
【００８７】
また、相関波形解析部１７は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、式（４）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を解析し、図１２に示すように画像表示する。
【００８８】
以上のように、本発明の実施形態によれば、伝搬空間内に存在する物標に進入する音源からの音波の信号に対する自己相関関数と、物標から前方に散乱する前方散乱波に対するパッシブ位相共役結果とが一致することを利用することにより、伝搬空間内に存在する物標を正確に探査することができる。
【００８９】
さらに、本発明の実施形態によれば、音源と基準として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施した進行波の自己相関処理信号と、減算処理で混合波から分離した前方散乱波にパッシブ位相共役処理を施した前方散乱波の共役性相関信号との相関を取ることにより、物標の位置を探査するため、物標の探査に要する時間を短縮することができる。何故ならば、共役性相関信号との相関を取る対象が、進行波に自己相関処理を施した自己相関処理信号であり、タイムリバーサル処理を施す必要がなく、時間を短縮することができるためである。
【産業上の利用可能性】
【００９０】
本発明によれば、広範囲の探知物を対象として、これらの探知物を容易に探知することができる。
【図面の簡単な説明】
【００９１】
【図１】本発明の実施形態に係る物標探査装置を示す構成図である。
【図２】本発明の実施形態における物標の探査処理を示す図である。
【図３】本発明の実施形態に係る物標探査装置を用いて物標の探査を行う場合を示すフローチャートである。
【図４】本発明の実施形態において、メッシュ処理部による任意の位置を仮定する様子を図示化して示す図である。
【図５】音源から送波するトーンバースト波を示す図である。
【図６】（ａ）は、元の音源から新しい音源に入射する音パルスの波形を示す図、（ｂ）は、時間反転した信号の波形を示す図である。
【図７】進行波と前方散乱波とが混合した混合波にパッシブ位相共役処理を施した信号の波形を示す図である。
【図８】音源から送波されて物標に入射する進行波にパッシブ位相共役処理を施した信号の波形を示す図である。
【図９】物標から前方に散乱した前方散乱波のみにパッシブ位相共役処理を施した信号の波形を示す図である。
【図１０】物標の範囲に関するパラメータとしてのパッシブ位相共役処理による波形を示す図である。
【図１１】物標の範囲を変更した場合におけるパッシブ位相共役処理の波形を示す図えある。
【図１２】マッチドフィールド法を適用して物標の位置を図式化して示す図である。
【符号の説明】
【００９２】
１ 伝搬空間
２ 物標
３ 音源
４ 進行波
５ 前方散乱波
６ トランスデューサアレイ
７ 減算処理手段
８ パッシブ位相共役処理手段
９ 自己相関処理手段
１０ 相関手段
Ｓ１ 自己相関処理手段
Ｓ２ 共役性相関信号
【技術分野】
【０００１】
本発明は、パッシブ位相共役処理を用いて物標を探査する物標探査装置、物標探査プログラム及び物標探査方法に関する。
【背景技術】
【０００２】
海中音響の分野において、信号を時間反転させる研究が増加している（非特許文献１〜３）。浅い海域において物標を探査する際、海中に送波した音波が海面や海底で反射し、これらの反射波が存在するために、物標の探査が困難になる場合がある。
【０００３】
近年、物標からの反射波に位相共役処理を施す探査法が注目されている（非特許文献４〜８）。これらの探査法は、物標から後方に散乱する後方散乱波を用いる。ここで、音波が物標に進入した場合について考える。音波が物標に進入すると、一部の音波が物標で散乱し、散乱波として音源側に逆行する。この散乱波を後方散乱波という。また、一部の音波が物標から前方に散乱する。この散乱波を前方散乱波という。
【非特許文献１】S.Kim, W. A. Kuperman, W. S. Hodgkiss, H. C. Song, G. F. Edelmann, and T. Akal, “Robusttime reversal focusing in the ocean,” J. Acoust. Soc. Am. 114, 145-157(2003)
【非特許文献２】S.C.Walker, P. Roux, and W.A. kuperman, “Focal depth shifting of a timereversal mirror in a rang-independent waveguide,” J. Acoust. Soc. Am. 118,1341-1347(2005)
【非特許文献３】S.C. Walker, W. A. Kuperman, and P. roux, “Active waveguide Green’s functionestimation with application to time-reversal focusing without a probe source ina rang-indepent waveguide,” J. Acoust. Soc. Am. 120, 2755-2763(2006)
【非特許文献４】G.Montaldo, M. Tanter, and M. Fink,” Revisiting iterative time reversalprocessing: Application to detection of multiple targets” J. Acoust. Soc. Am.115, 776-784(2004)
【非特許文献５】C.Prada, S. Manneville, D. Spoliansky, and M. Fink, “Decomposition of the timereversal operator: Detection and selective focusing on two scatterers,” J.Acoust. Soc. Am. 99, 2067-2076(1996)
【非特許文献６】N.Mordant, C. Prada, and M. Fink, “ Highly resolved detection and selectivefocusing in a waveguide using the D. O. R. T. method,” J. Acoust. Soc. Am. 105,2634-2642(1999)
【非特許文献７】D.H. chamberes, “Analysis of the time-reversal operator for scatterers of finitesize,” J. Acoust. Soc. Am. 112, 411-419(2002)
【非特許文献８】D.H. chamberes, A.K. Gautesen, “Time reversal for a single spherical scatterer,”J. Acoust. Soc. Am. 109, 2616-2624(2001)
【発明の開示】
【発明が解決しようとする課題】
【０００４】
非特許文献４〜８に開示された探査法では、物標から後方に散乱する後方散乱波を用いているが、小さい物標の場合や低周波の音波の場合、後方散乱波を受波する際、十分な後方散乱波を得られない場合がある。したがって、前記後方散乱波を用いた探査法では、物標を正確に探査できなくなる可能性が残されていた。
【０００５】
また、小さい物標の場合、一般的に前方散乱波のレベルが後方散乱波のレベルよりも高い。前記前方散乱波を受波する際、前方散乱波以外にも、音源から受波側に進行する進行波も、物標から離れた領域に存在する。前方散乱波は、音源からの音波が物標に進入して散乱したものであるから、前方散乱波が音源から直接進行する進行波に隠れてしまい、前方散乱波を用いて物標を正確に探査することができない場合がある。
【０００６】
本発明の目的は、音源から直接進行する進行波に混合した前方散乱波を用いて物標を正確に探査する物標探査装置、物標探査プログラム及び物標探査方法を提供することにある。
【課題を解決するための手段】
【０００７】
前記目的を達成するため、本発明に係る物標探査装置は、伝搬空間内に存在する物標を探査する物標探査装置であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波する音源と、
前記伝搬空間内に存在する物標から前方に散乱する前方散乱波を受波する領域に配置されたトランスデューサアレイと、
前記トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する減算処理手段と、
前記減算処理手段で分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成するパッシブ位相共役処理手段と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する自己相関処理手段と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する相関手段とを含むことを特徴とする。
【０００８】
以上の例では、本発明をハードウェアとしての物標探査装置として構築したが、これに限られるものではない。本発明は、物標を探索する方法として構築する、或いは物標探索装置の機能をコンピュータに実行させるソフトウェアとしてのプログラムとして構築してもよいものである。
【０００９】
本発明に係る物標探査プログラムは、伝搬空間内に存在する物標を探査する制御を行う物標探査プログラムであって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標の探査を制御するものであり、
コンピュータに、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する機能と、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成する機能と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する機能と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する機能とを実行させることを特徴とする。
【００１０】
本発明に係る物標探索方法は、伝搬空間内に存在する物標を探査する物標探査方法であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波し、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離し、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成し、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成し、
前記自己相関処理信号と前記共役性相関信号との相似性を判定することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、伝搬空間内に存在する物標に進入する音源からの音波の信号に対する自己相関関数と、物標から前方に散乱する前方散乱波に対するパッシブ位相共役結果とが一致することを利用することにより、伝搬空間内に存在する物標を正確に探査することができる。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の実施形態を図に基づいて説明する。
【００１３】
本発明の実施形態に係る物標探査装置は伝搬空間１内に存在する物標を探査する物標探査装置であって、図２に示すように、伝搬空間１内に存在する物標２に進入する音源３からの進行波４と、前記物標２から前方に散乱する前方散乱波５とを用いることにより、前記物標２を探査するものである。伝搬空間１として例えば海中を想定した場合、図２の１ａが海面であり、１ｂが海底であり、１が伝搬空間としての海中である。なお、音源３は海中１の定位置に設置する、或いは探査船から海中に吊下して設置するが、音源３の設置については特に限定されるものではない。図２において、音源３からの音波が物標２に入射した際に、音源３に向かう後方に生じる散乱波が後方散乱波５ａであり、トランスデューサアレイ６に向かう前方に生じる散乱波が前方散乱波５である。
【００１４】
本発明の実施形態は、物標２の位置に入射した音波の「自己相関関数」と、物標２から前方に散乱する前方散乱波５をトランスデューサアレイ６で受波してパッシブ位相共役処理を施した波形とが相似することを証明し、その証明結果に基づいて物標の探査を実現したことを特徴とするものである。
【００１５】
物標２の位置に入射した音波の「自己相関関数」と、物標２から前方に散乱する前方散乱波５をトランスデューサアレイ６で受波してパッシブ位相共役処理を施した波形とが相似することを検証する。
【００１６】
最初に位相共役について検証する。図２に示すように、例えば浅い海域において音源３とトランスデューサアレイ６とを対峙して設置する。音源３から送波される音波（進行波４）は、離れた地点に配置されたトランスデューサアレイ６により受波される。そして、トランスデューサアレイ６で受波した受波信号に対してパッシブ位相共役処理を行い、再びトランスデューサアレイ６から送波された場合、受波点における音波は、以下の式によって表される。
【数１】
式（１）は、文献 D. R. Jackson and D.R. Dowling, “Phase conjugation in underwater acoustics” J. Acoust. Soc. Am. 89, 171-181 (1991) で確立されている。
式（１）において、
【数２】
は音源３からトランスデューサアレイ６の各素子６ａへの伝搬に関するグリーン関数、
【数３】
はトランスデューサアレイ６の各素子６ａから受波点への伝搬に関するグリーン関数である。ｒは距離を、ｓ、ｎ、* は、それぞれ、音源、トランスデューサアレイの各素子の番号，複素共役を示している。
【００１７】
式（１）は周波数領域におけるものだが、以下の式では時間領域におけるものを示す。
【数４】
式（２）は、W. A. Kuperman, W. S. Hodgkiss, C. Song, T. Akal, C. Ferla and D. R. Jackson “Phase conjugation in the ocean; Experimental demonstration of an acoustic time-reversal mirror” J. Acoust. Soc. Am. 103, 25-40(1998) で確立されている。
式（２）において、Gω，ｔ，ｚはそれぞれ、グリーン関数，時間，深度である。Ｒは音源３からトランスデューサアレイ６までの水平距離、ｒはトランスデューサアレイ６から受波点との間の距離を示す。Ｓ（ω）は音源３から送波される音波の周波数スペクトル、ωは角周波数である。式（２）はタイムリバーサル波の音源近くの集束特性を調べるときにしばしば用いられる。
【００１８】
本発明者は、前記式（１），（２）が音源３の位置以外の受波点の音場にも適用可能であると考えている。このことについて考察する。すなわち、音源３の地点以外の任意の点の音場の位相共役について考察する。ここで、任意の点は、図２に示すように音源３とトランスデューサアレイ６との間にあるものと仮定する。
【００１９】
前記任意の点が新しい音源の位置であるものと考えられる。そして、新しい音源の駆動信号が元の音源３から供給される。したがって、元の音源３の周波数スペクトルＳ（ω）と、新しい音源の周波数スペクトルＳＲ（ω）の関係は以下のように表される。
【数５】
式（３）において、添字Ｒは任意の点をあらわす。
【００２０】
次に、式（２），（３）を用いて、音源３以外の任意の点の位相共役について考察する。一般的な位相共役波の場合、音源３から送波される音波パルスは図２に示す処理のため、トランスデューサアレイ６で受波され、その各素子６ａ毎に時間反転処理（タイムリバース処理）が行われる。その後、時間反転された信号（時間反転信号）が再び各素子６ａから送波され、音源３から送波される音波パルスと同様の音波パルスが音源３の位置に形成される。
【００２１】
本発明の実施形態において、上述した一般的な位相共役波と同様の特性が任意の地点で維持されているか否かを検証する。すなわち、音源３とトランスデューサアレイ６との間の任意の地点において、音源３から直接進行する音波（以下、直接波という）と時間反転信号との相似性を検証する。
【００２２】
シミュレーションの結果について考察する。先ず、任意の点における直接波と時間反転波の比較について考察する。
【００２３】
図２に示す処理において、海底１ｂまでの深さは１００ｍに設定し、音源３とトランスデューサアレイ６との間の距離は２．５ｋｍに設定してある。海水中の音速及び海水の密度はそれぞれ１５００ｍ／ｓ及び１０００ｋｇ／ｍ３である。そして、海底１ｂの堆積物における音速及び堆積物の密度はそれぞれ１６００ｍ／ｓ及び１５００ｋｇ／ｍ３である。
【００２４】
中心周波数５００Hz、パルス幅８サイクルのトーンバースト波が、深度５０ｍに設置された音源３から送波される。トーンバースト波の波形を図５（ａ）、そのスペクトルを図５（ｂ）にそれぞれ示す。
【００２５】
音波パルスの波形を詳しく調べるシミュレーションにおいてはワイドバンドパルスが用いられる。これらの音波パルスはトランスデューサアレイ６で受波され、その各素子６ａ毎に時間反転処理が施される。そして、前記時間反転処理した信号が各素子６ａから再送波される。再送波された音波は音源３とトランスデューサアレイ６との間（図２の任意の点）で受波される。
【００２６】
そして、時間反転信号と受波点で形成された駆動信号（音源３から供給される音波パルス）との相似性について検証する。音源３から１ｋｍの地点において、音源３から任意の点に直接進入する音波と、時間反転信号とを比較する。この場合、元の音源３から伝わる音波パルスはトーンバースト波である。図５（ｂ）に示すスペクトルは式（２）のＳ（ω）に相当する。
【００２７】
計算に用いる周波数スペクトルのバンド幅は３５０Ｈｚである。その理由は、周波数スペクトルのハンド幅を３５０Ｈｚに設定することが、任意の点における波形を詳細に表すからである。このバンド幅は、図５（ｂ）に示す周波数スペクトルの第３ローブまでのバンド幅に相当する。一方、新しい音源の周波数スペクトルは、元の音源３から任意の位置に直接進入する音波パルスのフーリェ変換により、容易に得られる。
【００２８】
新しい音源の周波数スペクトルＳR（ω）は、元の音源３の周波数スペクトルＳ（ω）と、式（３）により関連付けられる。新しい音源の位置は、元の音源３に対して距離１ｋｍ、深度５０ｍの地点である。元の音源３から新しい音源の位置に進入する音波パルスの波形を図６（ａ）、時間反転信号の波形を図６（ｂ）にそれぞれ示す。
【００２９】
図６（ｂ）に示す時間反転信号の進行方向は、図６（ａ）に示す直接波の進行方向と反対であるため、時間反転信号は時間を反転した座標で表してある。さらに、式（２）に示す周波数スペクトルＳ（ω）は、受波点に進入する音波を新しい音源と考えることによって、元の音源３からの音波パルスのフーリェ変換により直接求められる。図６（ａ），（ｂ）を比較すると、元の音源３から直接進入する音波の波形は時間反転信号の波形と一致する。
【００３０】
上述した検証結果から、位相共役性が音源３とトランスデューサアレイ６との間の任意の点においても維持されていることが分かる。
【００３１】
次に、音源３とトランスデューサアレイ６とを浅い海域に配置する。プローブ信号Ｐiとデータ信号Ｐdが音源３から送られ、音源３からの音波がトランスデューサアレイ６のｍ番目の素子で受波された場合（位置ｒm）、プローブ信号Ｐiとデータ信号Ｐdとの信号の相互相関関数は以下の式（４）によって表される。
【数６】
【００３２】
次に、トランスデューサアレイ６のすべての素子６ａの相互相関関数は以下の式（５）によって加えられる。
【数７】
Ｓ（ｔ）はプローブ信号の位相共役になることが知られている。この技術はパッシブ位相共役と呼ばれる。
なお、式（４），（５）は、D. Rouseff, D. R. Jackson, W. L. J. Fox, C. D. Jones, J. Ritcey and D. R. Dowling, “Underwater Acoustic Communication by Passive-Phase Conjugation; Theory and Experimental Results” IEEE J. Ocean, Eng. 26, 821-831 (2001) で確立されている。
【００３３】
音源３から進行する音波に対するパッシブ位相共役処理について考察する。図２のように、音源３とトランスデューサアレイ６とを浅い海域に配置する。音源３からトランスデューサアレイ６までの距離は３ｋｍである。物標２は音源３から距離１．５ｋｍ、深度５０ｍの地点に配置する。物標２の材料は鉄であり、高さと幅がともに２mの大きさのものである。
【００３４】
中心周波数５００Ｈｚ、パルス幅８サイクルのトーンバースト波が音源３から送波され、音波パルスが深度５０ｍに設置された音源３から送波される。音源３から送波された音波パルスは物標２に入射し、物標２から前方に散乱する前方散乱波５が発生する。
【００３５】
音源３から音波を送波すると、音源３からの進行波４と前方散乱波５とがトランスデューサアレイ６にほぼ同時に到達する。そして、トランスデューサアレイ６で受信した信号にパッシブ位相共役処理を施す。すなわち、トランスデューサアレイ６の各素子６ａで受波した信号に式（４）による相関処理を行う。さらに、トランスデューサアレイ６のすべての素子６ａの相互相関関数（式（５））を加える。この一連の相関処理がパッシブ位相共役処理である。
【００３６】
図７は、上述したパッシブ位相共役処理を経て得た進行波４と前方散乱波５との混合波のパッシブ位相共役信号を表している。図７から明らかなように、物標２による前方散乱波の影響は見られず、波形は図５（ａ）に示す音波パルスの自己相関関数の波形とほぼ一致する。パッシブ位相共役処理は、式（４）におけるプローブ信号Ｐｉとデータ信号Ｐｄとが同じであるから、混合波の自己相関関数となる。
【００３７】
一般に、音源３から送波される進行波４のレベルは前方散乱波５のレベルより大きく、前方散乱波５を検出するのは困難である。本発明の実施形態では、物標２が存在しないときに、予め音源３から音波（進行波４）を送波し、トランスデューサアレイ６の各素子６ａで受波した信号を保持しておく。その後、物標２が存在するときに音源３から音波（進行波４）を送波し、トランスデューサアレイ６の各素子６ａにより進行波４と前方散乱波５との混合波を受波し、混合波から進行波４を減じることにより、前方散乱波５を分離して顕在させる。
【００３８】
次に、物標２から前方に散乱した前方散乱波５に対するパッシブ位相共役処理について検証する。
【００３９】
本発明の実施形態では、物標２が存在するときにトランスデューサアレイで受波した混合波の信号から、物標２が存在しないときにトランスデューサアレイで受波した進行波４の信号を減ずる。そして、上述したパッシブ位相共役処理を、減じて顕在させた前方散乱波５に対して施す。前記減算した信号（前方散乱波５）に対して施すパッシブ位相共役処理は、式（４）におけるプローブ信号Ｐｉとデータ信号Ｐｄとが同じであるから、減算信号の自己相関関数となる。
【００４０】
物標２は距離１．５ｋｍ、深度５０mの地点に配置され、音波パルスが深度５０ｍに設置された音源３から送波される。トランスデューサアレイ６で受波した混合波から進行波４の成分を取り除いた後、パッシブ位相共役処理を施す。そのパッシブ位相共役処理を行った結果の波形を図８に示す。図８に示す波形は、音波パルスの繰り返し周期である１０秒を中心とした対称的な波形である。
【００４１】
以上のようにして得られた図８に示すパッシブ位相共役処理による波形が、物標２の位置を新しい音源の位置としてその音源から前方に放射された音波パルスにパッシブ位相共役処理を施した波形であることを確かめる必要がある。これを確かめるために、物標２が新しい音源の位置（任意の位置）にあることを確認する。この場合、前記新しい音源から放射される音波パルスは、トーンバースト波ではない。前記新しい音源が音波パルスを放射する駆動信号は、元の音源３から前記新しい音源の位置に直接進入する音波パルスである。前記新しい音源から放射されるパルスがトランスデューサ６で受波され、その受波信号にパッシブ位相共役処理を施す。
【００４２】
図９は、前記新しい音源から放射された信号にパッシブ位相共役処理を施した信号の波形を示している。図９に示す波形は図８に示す波形と実質的に一致している。図９に示す信号は、トランスデューサアレイ６で受波した混合波から進行波４を減算処理して得た物標２からの前方散乱波５のみであることが確認できた。
【００４３】
次に、物標２の位置とパッシブ位相共役処理との関係について考察する。伝搬空間１における伝搬環境は、これまで検証した際に記述したものと同じである。物標２の位置のみを変更しており、位置を変更した物標２から前方に散乱する前方散乱波５をトランスデューサアレイ６で受波し、その受波信号にパッシブ位相共役処理を施す。
【００４４】
図１０(ａ)，（ｂ），（ｃ）は、物標２の距離に関するパラメータとしてのパッシブ位相共役処理による波形を示している。物標２の距離（音源３に対する位置）は、１．３ｋｍ，１．６ｋｍ，１．７ｋｍである。また、物標２の深度はすべて５０ｍである。図１０から明らかなように、前記パッシブ位相共役処理による波形の変化は、物標２の距離によるところが大きいことがわかる。同様に、物標２の深度が変わると、前記パッシブ位相共役処理による波形は大きく変化する。図１１は、物標２の距離が１．５ｋｍであって、物標２の深度が２０ｍに変化した場合における前記パッシブ位相共役処理の波形を示している。図１０及び図１１からして、物標２の位置と前記パッシブ位相処理による波形とは、密接に関係していると考えられる。
【００４５】
以上の検証結果から、物標２から前方に散乱する前方散乱波５に対するパッシブ位相共役処理の結果は、音源３から送波されて物標２に進入する音波パルスと直接関連する。したがって、物標２から前方に散乱する前方散乱波５に対してパッシブ位相共役処理を施して得られる波形は、音源３から送波されて物標２の位置に進入する音波パルスの相関波形（自己相関関数）と同様である。一方、音源３を基準とする音場の任意の点に進入する音波パルスは伝搬環境により一意に決められる。すなわち、物標２から前方に散乱する前方散乱波５に対するパッシブ位相共役処理の結果は、物標２の位置に関する情報を保持する。
【００４６】
以上の関係はマッチドフィールド法として物標２の探査に用いることが可能である。物標２の音源３に対する距離が１．５ｋｍであり、物標２の深度が５０mにある場合、音源３を基準とする任意の点に進入する音波パルスの自己相関関数と、任意の点から前方に放射される音波パルスに施したパッシブ位相共役処理の結果との相関関数とは図１１に示すようになる。
【００４７】
図１２から明らかなように、前方散乱波５にパッシブ位相共役処理を行って生成した信号（前方散乱波の共役性相関信号）と、音源３を基準位置として伝搬空間内の任意の位置に存在する物標２の位置に入射する進行波４に自己相関処理を行って生成した信号（進行波の自己相関処理信号）との相関により得た相関値の最も高い位置Ｈが物標２の位置であることを明確に示している。
【００４８】
以上の検証結果を考察すると、音源３を基準とする伝搬空間の任意の位置から放射される音波パルスにパッシブ位相共役処理を施した結果と、前記任意の位置に音源３から進入する音波パルスの自己相関処理の結果とが一致することが証明できた。前記任意の位置に物標２が存在すれば、その前方散乱波は物標２に進入した音波パルスに対応する。したがって、物標２の位置に音源３から進入する進行波の自己相関処理の結果（進行波の自己相関処理信号）と、物標２から前方に散乱する前方散乱波５の信号にパッシブ位相共役処理を施した結果（前方散乱波の共役性相関信号）とが一致することとなり、進行波の自己相関処理信号と前方散乱波５の共役性相関処理信号との相関値が高く、任意の位置に物標２が存在しない場合の相関値が低くなる。
【００４９】
本発明の実施形態は、音源３から送波されて音源３を基準とした伝搬空間１の任意の位置に入射する音波（進行波４）のパッシブ位相共役処理の結果である共役性相関信号Ｓ２と、音源３から送波されて前記任意の位置と同一位置に入射する音波の自己相関処理の結果である進行波４の自己相関処理信号Ｓ１とが一致するか否かを監視することにより、音波が伝搬する伝搬空間１、例えば浅い海域における物標２を探査することを特徴とするものである。任意の位置に物標２が存在する場合に共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相関値が最も高くなるためである。以下、具体的に説明する。
【００５０】
本発明の実施形態に係る物標探査装置は図１に示すように、基本的な構成として、伝搬空間１内に存在する物標２に進入する音源からの進行波４と、前記物標２から前方に散乱する前方散乱波５とを用いることにより、前記物標２を探査するものである。そして、伝搬空間１内に音波を送波する音源３（図２参照）と、伝搬空間１内に存在する物標２から前方に散乱する前方散乱波５を受波する領域に配置されたトランスデューサアレイ６（図２参照）と、トランスデューサアレイ６に直接進行する進行波４と前方散乱波５との混合波から、トランスデューサアレイ６に直接進行した進行波４を減じることにより、前方散乱波５を分離する減算処理手段７と、減算処理手段７で分離した前方散乱波５にパッシブ位相共役処理を行うことにより、前方散乱波５の共役性相関信号Ｓ２を生成するパッシブ位相共役処理手段８と、音源３を基準位置として伝搬空間１内の任意の位置に入射する進行波４に自己相関処理を施すことにより、進行波４の自己相関処理信号Ｓ１を生成する自己相関処理手段９と、自己相関処理信号Ｓ１と共役性相関信号Ｓ２との相似性を判定する相関手段１０とを含むことを特徴とする。
【００５１】
本発明の実施形態では、伝搬空間１内に音波（進行波４）を音源３から送波し、トランスデューサアレイ６に直接進行する進行波４と前方散乱波５との混合波から、トランスデューサアレイ６に直接進行した進行波４を減じることにより、前方散乱波５を分離する。前方散乱波５を分離するにあたっては、物標２が存在しない伝搬空間１内を伝搬する進行波４をトランスデューサアレイ６で受波して進行波４の信号を予め取得しておき、次に前記混合波をトランスデューサアレイ６で受波し、前記混合波の信号から、トランスデューサアレイ６に直接進行した進行波４の信号を減じることにより、前方散乱波５を分離する。
【００５２】
次に、前記分離した前方散乱波５にパッシブ位相共役処理をパッシブ位相共役処理手段８で行うことにより、前方散乱波５の共役性相関信号Ｓ２を生成する。また、前記音源３を基準位置として伝搬空間１内の任意の位置に入射する進行波４に自己相関処理を自己相関処理手段９で施すことにより、進行波４の自己相関処理信号Ｓ１を生成する。そして、自己相関処理信号Ｓ１と共役性相関信号Ｓ２との相似性を相関手段１０で判定する。
【００５３】
なお、以上の説明では、本発明の実施形態に係る物標探査装置をハードウェアとして構築する場合について説明したが、これに限られるものではない。すなわち、本発明の実施形態に係る物標探査装置が実行する機能をプログラムとして構築し、そのプログラムをコンピュータに実行させることにより、本発明の実施形態に係る物標探査装置が実行する機能を実現するようにしてもよいものである。この場合、前記プログラムは記録媒体に記録され、その状態で商取引の対象となる。
【００５４】
本発明の実施形態に係る物標探査装置をさらに具体例を用いて詳細に説明する。
【００５５】
自己相関処理手段９は、環境情報生成部１１と、波形生成処理部１２と、メッシュ処理部１３と、自己相関処理部１４とを含んでいる。
【００５６】
環境情報生成部１１は、音源３から送波されてトランスデューサアレイ６で受波変換された校正信号を取り込み、不要な信号を除去すると共にＡ／Ｄ変換し、これを校正信号として用いる。さらに、環境情報生成部１１は、前記校正信号に基づいて実海域の定数、すなわち、音波ビームが伝搬する実際の海域の定数（海洋の音響定数）を演算する。具体的に説明する。
【００５７】
本発明の実施形態に係る物標探査装置は、例えば大陸棚のような比較的浅い海域で用いられる。その環境における音波伝搬に関する海洋の音響定数は、水深、海水の音速、海底堆積物（音速，密度）等のデータである。また、前記環境における音波パルスに関する海洋の音響定数は、音波パルスの中心周波数、音波パルスのスペクトル及び音波パルスのバンド幅等のデータである。海洋の音響定数を決定するための水深のデータは、測深機などの測定機器を用いることで直接計測により求められる。海洋の音響定数を決定するための海水の音速データは、水温の計測により、例えばメドイン（H.Medwin）の式などの既知の音速の式から求められる。
【００５８】
海洋の音響定数を決定するための海底堆積物のデータは、海図などから大凡の値を推定できるが、より正確な値は、既知の校正法により求められる。すなわち、物標が存在していると予想される探知海域に吊下した校正音源から送波された信号をトランスデューサアレイ６で受波し、その受波信号と予測した海底堆積物の定数とを変化させながら演算することで、実海域での堆積物のデータを取得する。ここに、海洋の音響定数とは、実際に音波が海中を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる環境因子を意味する。なお、海中での物標２を探索する場合を想定しているので、海洋音響定数を用いたが、これに限られるものではない。伝搬空間として海洋以外のものを用いる場合、前記音響定数は、実際に音波が伝搬空間を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる環境因子を意味することとなる。
【００５９】
環境情報生成部１１には、上述した音波伝搬及び音波パルスに関する実海域での海洋の音響定数を決定するためのデータが入力され、環境情報生成部１１は、前記入力された海洋の音響定数を決定するためのデータと、前記校正信号とに基づいて実海域の海洋の音響定数を決定する。環境情報生成部１１は、決定した実海域での海洋の音響定数を波形生成処理部１２に向けて出力する。なお、海中での物標２を探索する場合を想定しているので、海洋音響定数を用いたが、これに限られるものではない。伝搬空間として海洋以外のものを用いる場合、環境情報生成部１１は、実際に音波が伝搬空間を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる実際の音響定数を出力することとなる。
【００６０】
メッシュ処理部１３は図４に示すように、音源３とトランスデューサアレイ６との間の伝搬空間１内に物標２の位置を仮定し、その仮定した任意の位置情報を波形生成処理部１２に出力する。具体的に説明すると、伝搬空間１が例えば海中であるとすると、メッシュ処理部１３は図４に示すように、音源３とトランスデューサアレイ６との間の距離をｒとして任意の位置（距離）を水平座標上のｘ−ｙ面に設定する。さらに、メッシュ処理部１３は図４に示すように、海面１ａから海底１ｂまでの深度をｚとして任意の位置（深度）を垂直座標上のｙ−ｚ面に設定する。
【００６１】
そして、メッシュ処理部１３は、音源３に対する距離ｒを変更させて任意の位置をｘ−ｙ面内に順次移動させるとともに、音源３に対する深度ｚを変更させて任意の位置をｙ−ｚ面内に順次移動させ、その任意の位置情報（距離，深度）を波形生成処理部１２に出力する。
【００６２】
波形生成処理部１２は、環境情報生成部１１から出力される音響定数とメッシュ処理部１３から出力される位置情報とに基づいて、音源３から送波された音波が入射する点における新しい音源、すなわち音源３を基準とした音場における任意の位置における音波パルスを式（６）及び式（７）に基づいて求める。
【数８】
式（６）で求める音波パルスは音源３から任意の位置まで通常に伝搬した音波パルスである。式（６）において、Ｓ（ω）は音源３から送波される音波パルスの周波数スペクトルである。Ｐ（ｒ，ｚ，ω）は前記周波数スペクトルＳ（ω）に対応した任意の点における音圧であり、通常のノーマルモード法等により求める。ωｂはバンド幅である。
【００６３】
波形生成処理部１２は、例えば結合モード法を用いることにより、角周波数ωの音波の任意の点（ｒ，ｚ）における音圧Ｐ（ｒ，ｚ，ω）を式（７）に基づいて求め、これを式（６）に代入する。
【数９】
式（７）において、ｊは水平方向の区分領域（距離）を意味している。φ（ｚ，λｊ，ｍ）はモード関数である。
式（７）は、R. B. Evans, “A coupled mode solution for acoustic propagation in a wave guide with stepwise depth variations of a penetrable bottom.“ J. Acoust. Soc. Am., 74, 188-193(1983) によって確立されている。
【００６４】
波形生成処理部１２は、式（６），式（７）に基づいて計算した音源３を基準とする伝搬空間１内の任意の位置での音波（パルス）を自己相関処理部１４に出力する。
【００６５】
自己相関処理部１４は、波形生成処理１２から受け取った音波情報に基づいて、その音波信号に式（４）に基づく相関処理を行う。この場合、式（４）におけるプローブ信号Ｐｉとデータ信号Ｐｄとは同じであり、任意の位置に入射する進行波４を相関処理することは、自己相関処理の形態となる。したがって、自己相関処理部１４は、音源３を基準とする伝搬空間１内での任意の位置に入射する進行波４に自己相関処理を施すことにより、進行波４の自己相関処理信号Ｓ１を生成する。
【数１０】
【００６６】
自己相関処理部１４が出力する自己相関処理信号Ｓ１は、物標２が前記任意の位置に存在する場合、図８に示す波形の信号となる。
【００６７】
音源３から送波されて伝搬空間１内を伝搬した音波パルスがトランスデューサアレイ６で受波され、トランスデューサアレイ６は、受波音波を電気信号である受波信号に変換し、その受波信号を減算処理手段７に出力する。
【００６８】
減算処理手段７は、予め伝搬空間１内に物標２が存在しない場合に音源３からトランスデューサアレイ６に直接進行する進行波４の受波信号を記憶しておく。伝搬空間１内に物標２が存在する場合に前方散乱波５を含む進行波４がトランスデューサアレイ６に受波された場合、減算処理手段７は、進行波４と前方散乱波５との混合波から進行波４を減算することにより、前方散乱波５を分離し、その分離した前方散乱波５の受波信号をパッシブ位相共役処理手段８に出力する。
【００６９】
パッシブ位相共役処理手段８は、減算処理手段７で減算処理された減算信号を受け取ると、前記減算信号に式（４）及び式（５）に基づくパッシブ位相共役処理を施すことにより、前方散乱波５の共役性相関信号Ｓ２を生成する。パッシブ位相共役処理手段８は、前方散乱波５の共役性相関信号Ｓ２を相関手段１０に出力する。
【数１１】
【数１２】
【００７０】
パッシブ位相共役処理手段８は、前記任意の位置に物標２が存在する場合、図９に示すような波形の共役性相関信号Ｓ２を出力する。プローブ信号Ｐｉとデータ信号Ｐｄとは同一であるから、パッシブ位相共役処理手段８によるパッシブ位相共役処理は自己相関処理の形態となる。
【００７１】
相関手段１０は、相関処理部１６と、相関波形解析部１７とを含んでいる。
【００７２】
相関処理部１６は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、通常の相互相関式（８）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を判定する。相関処理部１６は、前記自己相関処理信号Ｓ１と前記共役性相関信号Ｓ２との相関で得た相関値が最も高い位置を物標位置として判定する。
【数１３】
式（８）において、Ｓ１（ｎ）は自己相関処理信号Ｓ１のデータ、Ｓ２（ｎ＋ｋ）は鏡や規制相関信号Ｓ２のデータである。
式（８）は、例えば、「ディジタル信号処理 江原義郎 著 P46-P47 東京電気大学出版局 １９９６年２月２０日発行」などに掲載されている相互相関関数である。
【００７３】
相関波形解析部１７は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、式（４）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を解析し、画像表示する。相関波形解析部１７が画像表示する一例を図１２に示す。
【００７４】
制御部１８は、自己相関処理手段９，パッシブ位相共役処理手段８，相関手段１０及びその他の構成要素を総合的に制御するものである。２０は、制御部１８の動作に必要な情報、及び演算に必要な作業領域を提供する記憶部、１９は、制御部１８を通して出力されるデータを外部に出力する指示部である。
【００７５】
次に、本発明の実施形態に係る物標探査装置を用いて伝搬空間、例えば海中１内に存在する物標２を探査する場合について説明する。
【００７６】
環境情報生成部１１は、音源３から送波されてトランスデューサアレイ６で受波変換された校正信号を取り込み、不要な信号を除去すると共にＡ／Ｄ変換し、これを校正信号として用い、前記校正信号に基づいて実海域の定数、すなわち、音波ビームが伝搬する実際の海域の定数（海洋の音響定数）を演算し、実際に音波が伝搬空間を伝搬する際に前記音波の伝搬に影響を及ぼす要因となる実際の音響定数を波形生成処理部１２に出力する。
【００７７】
一方、メッシュ処理部１３は、音源３に対する距離ｒを変更させて任意の位置をｘ−ｙ面内に順次移動させるとともに、音源３に対する深度ｚを変更させて任意の位置をｙ−ｚ面内に順次移動させ、その任意の位置情報（距離，深度）を波形生成処理部１２に出力する。
【００７８】
波形生成処理部１２は、環境情報生成部１１から出力される音響定数とメッシュ処理部１３から出力される位置情報とに基づいて、音源３から送波された音波が入射する点における新しい音源、すなわち音源３を基準とした音場における任意の位置における音波パルスを式（６）及び式（７）に基づいて求め、伝搬空間１内の任意の位置での音波パルスを自己相関処理部１６に出力する。
【００７９】
自己相関処理部１４は、波形生成処理部１２から受け取った音波情報に基づいて、その音波信号に式（４）に基づく相関処理を行い、音源３を基準とする伝搬空間１内での任意の位置に入射する進行波４に自己相関処理を施すことにより、進行波４の自己相関処理信号Ｓ１を生成する（図３のステップ１００）。
【００８０】
自己相関処理部１４が出力する自己相関処理信号Ｓ１は、物標２が前記任意の位置に存在する場合、図８に示す波形の信号となる。
【００８１】
自己相関処理信号Ｓ１を生成する処理を行う際に、共役性相関信号Ｓ２を生成する処理を音源３からトーンバースト波を出力するタイミングで行う。
【００８２】
すなわち、音源３から送波されて伝搬空間１内を伝搬した音波パルスをトランスデューサアレイ６で受波し、トランスデューサアレイ６は、受波音波を電気信号である受波信号に変換し、その受波信号を減算処理手段７に出力する。
【００８３】
減算処理手段７は、予め伝搬空間１内に物標２が存在しない場合に音源３からトランスデューサアレイ６に直接進行する進行波４の受波信号を記憶しておく。伝搬空間１内に物標２が存在する場合に前方散乱波５を含む進行波４がトランスデューサアレイ６に受波された場合、減算処理手段７は、進行波４と前方散乱波５との混合波から進行波４を減算することにより、前方散乱波５を分離し、その分離した前方散乱波５の受波信号をパッシブ位相共役処理手段８に出力する。
【００８４】
パッシブ位相共役処理手段８は、減算処理手段７で減算処理された減算信号を受け取ると、前記減算信号に式（４）及び式（５）に基づくパッシブ位相共役処理を施すことにより、前方散乱波５の共役性相関信号Ｓ２を生成する。パッシブ位相共役処理手段８は、前方散乱波５の共役性相関信号Ｓ２を相関手段１０に出力する（図３のステップ２００）。
【００８５】
パッシブ位相共役処理手段８は、前記任意の位置に物標２が存在する場合、図９に示すような波形の共役性相関信号Ｓ２を出力する。プローブ信号Ｐｉとデータ信号Ｐｄとは同一であるから、パッシブ位相共役処理手段８によるパッシブ位相共役処理は自己相関処理の形態となる。
【００８６】
相関処理部１６は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、式（８）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を判定する（図３のステップ３００）。相関処理部１６が、自己相関処理信号Ｓ１と共役性相関信号Ｓ２との相関で得た相関値が最も高い位置を物標位置として判定する。
【００８７】
また、相関波形解析部１７は、自己相関処理手段９から進行波４の自己相関処理信号Ｓ１と、パッシブ位相共役処理手段８から前方散乱波５の共役性相関信号Ｓ２とを受け取ると、式（４）に基づいて、共役性相関信号Ｓ２と自己相関処理信号Ｓ１との相似性を解析し、図１２に示すように画像表示する。
【００８８】
以上のように、本発明の実施形態によれば、伝搬空間内に存在する物標に進入する音源からの音波の信号に対する自己相関関数と、物標から前方に散乱する前方散乱波に対するパッシブ位相共役結果とが一致することを利用することにより、伝搬空間内に存在する物標を正確に探査することができる。
【００８９】
さらに、本発明の実施形態によれば、音源と基準として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施した進行波の自己相関処理信号と、減算処理で混合波から分離した前方散乱波にパッシブ位相共役処理を施した前方散乱波の共役性相関信号との相関を取ることにより、物標の位置を探査するため、物標の探査に要する時間を短縮することができる。何故ならば、共役性相関信号との相関を取る対象が、進行波に自己相関処理を施した自己相関処理信号であり、タイムリバーサル処理を施す必要がなく、時間を短縮することができるためである。
【産業上の利用可能性】
【００９０】
本発明によれば、広範囲の探知物を対象として、これらの探知物を容易に探知することができる。
【図面の簡単な説明】
【００９１】
【図１】本発明の実施形態に係る物標探査装置を示す構成図である。
【図２】本発明の実施形態における物標の探査処理を示す図である。
【図３】本発明の実施形態に係る物標探査装置を用いて物標の探査を行う場合を示すフローチャートである。
【図４】本発明の実施形態において、メッシュ処理部による任意の位置を仮定する様子を図示化して示す図である。
【図５】音源から送波するトーンバースト波を示す図である。
【図６】（ａ）は、元の音源から新しい音源に入射する音パルスの波形を示す図、（ｂ）は、時間反転した信号の波形を示す図である。
【図７】進行波と前方散乱波とが混合した混合波にパッシブ位相共役処理を施した信号の波形を示す図である。
【図８】音源から送波されて物標に入射する進行波にパッシブ位相共役処理を施した信号の波形を示す図である。
【図９】物標から前方に散乱した前方散乱波のみにパッシブ位相共役処理を施した信号の波形を示す図である。
【図１０】物標の範囲に関するパラメータとしてのパッシブ位相共役処理による波形を示す図である。
【図１１】物標の範囲を変更した場合におけるパッシブ位相共役処理の波形を示す図えある。
【図１２】マッチドフィールド法を適用して物標の位置を図式化して示す図である。
【符号の説明】
【００９２】
１ 伝搬空間
２ 物標
３ 音源
４ 進行波
５ 前方散乱波
６ トランスデューサアレイ
７ 減算処理手段
８ パッシブ位相共役処理手段
９ 自己相関処理手段
１０ 相関手段
Ｓ１ 自己相関処理手段
Ｓ２ 共役性相関信号
【特許請求の範囲】
【請求項１】
伝搬空間内に存在する物標を探査する物標探査装置であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波する音源と、
前記伝搬空間内に存在する物標から前方に散乱する前方散乱波を受波する領域に配置されたトランスデューサアレイと、
前記トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する減算処理手段と、
前記減算処理手段で分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成するパッシブ位相共役処理手段と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する自己相関処理手段と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する相関手段とを含むことを特徴とする物標探査装置。
【請求項２】
前記相関手段が、前記自己相関処理信号と前記共役性相関信号との相関で得た相関値が最も高い位置を物標位置として判定する請求項１に記載の物標探査装置。
【請求項３】
前記音源が、前記進行波としてパルス状のトーンバースト波を送波する請求項１に記載の物標探索装置。
【請求項４】
伝搬空間内に存在する物標を探査する制御を行う物標探査プログラムであって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標の探査を制御するものであり、
コンピュータに、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する機能と、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成する機能と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する機能と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する機能とを実行させることを特徴とする物標探査プログラム。
【請求項５】
前記コンピュータに、
前記自己相関処理信号と前記共役性相関信号との相関で得た相関値が最も高い位置を物標位置として判定する機能を実行させる請求項４に記載の物標探査プログラム。
【請求項６】
伝搬空間内に存在する物標を探査する物標探査方法であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波し、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離し、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成し、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成し、
前記自己相関処理信号と前記共役性相関信号との相似性を判定することを特徴とする物標探査方法。
【請求項７】
前記自己相関処理信号と前記共役性相関信号との相関で得た相関値が最も高い位置を物標位置として判定する請求項６に記載の物標探査方法。
【請求項８】
前記進行波としてパルス状のトーンバースト波を送波する請求項６に記載の物標探索方法。
【請求項１】
伝搬空間内に存在する物標を探査する物標探査装置であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波する音源と、
前記伝搬空間内に存在する物標から前方に散乱する前方散乱波を受波する領域に配置されたトランスデューサアレイと、
前記トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する減算処理手段と、
前記減算処理手段で分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成するパッシブ位相共役処理手段と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する自己相関処理手段と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する相関手段とを含むことを特徴とする物標探査装置。
【請求項２】
前記相関手段が、前記自己相関処理信号と前記共役性相関信号との相関で得た相関値が最も高い位置を物標位置として判定する請求項１に記載の物標探査装置。
【請求項３】
前記音源が、前記進行波としてパルス状のトーンバースト波を送波する請求項１に記載の物標探索装置。
【請求項４】
伝搬空間内に存在する物標を探査する制御を行う物標探査プログラムであって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標の探査を制御するものであり、
コンピュータに、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離する機能と、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成する機能と、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成する機能と、
前記自己相関処理信号と前記共役性相関信号との相似性を判定する機能とを実行させることを特徴とする物標探査プログラム。
【請求項５】
前記コンピュータに、
前記自己相関処理信号と前記共役性相関信号との相関で得た相関値が最も高い位置を物標位置として判定する機能を実行させる請求項４に記載の物標探査プログラム。
【請求項６】
伝搬空間内に存在する物標を探査する物標探査方法であって、
伝搬空間内に存在する物標に進入する音源からの進行波と、前記物標から前方に散乱する前方散乱波とを用いることにより、前記物標を探査するものであり、
伝搬空間内に音波を送波し、
トランスデューサアレイに直接進行する進行波と前記前方散乱波との混合波から、前記トランスデューサアレイに直接進行した進行波を減じることにより、前記前方散乱波を分離し、
前記分離した前記前方散乱波にパッシブ位相共役処理を行うことにより、前記前方散乱波の共役性相関信号を生成し、
前記音源を基準位置として伝搬空間内の任意の位置に入射する進行波に自己相関処理を施すことにより、前記進行波の自己相関処理信号を生成し、
前記自己相関処理信号と前記共役性相関信号との相似性を判定することを特徴とする物標探査方法。
【請求項７】
前記自己相関処理信号と前記共役性相関信号との相関で得た相関値が最も高い位置を物標位置として判定する請求項６に記載の物標探査方法。
【請求項８】
前記進行波としてパルス状のトーンバースト波を送波する請求項６に記載の物標探索方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【公開番号】特開２０１０−８２９５（Ｐ２０１０−８２９５Ａ）
【公開日】平成２２年１月１４日（２０１０．１．１４）
【国際特許分類】
【出願番号】特願２００８−１６９６３０（Ｐ２００８−１６９６３０）
【出願日】平成２０年６月２７日（２００８．６．２７）
【出願人】（０００００４２３７）日本電気株式会社 (19,353)
【出願人】（５０５３６３１７８）
【Ｆターム（参考）】
【公開日】平成２２年１月１４日（２０１０．１．１４）
【国際特許分類】
【出願日】平成２０年６月２７日（２００８．６．２７）
【出願人】（０００００４２３７）日本電気株式会社 (19,353)
【出願人】（５０５３６３１７８）
【Ｆターム（参考）】
[ Back to top ]