CN105487077B - 基于激光多点相干探测的水下发声目标位置估计方法及实现该方法的装置 - Google Patents
基于激光多点相干探测的水下发声目标位置估计方法及实现该方法的装置 Download PDFInfo
- Publication number
- CN105487077B CN105487077B CN201511021900.5A CN201511021900A CN105487077B CN 105487077 B CN105487077 B CN 105487077B CN 201511021900 A CN201511021900 A CN 201511021900A CN 105487077 B CN105487077 B CN 105487077B
- Authority
- CN
- China
- Prior art keywords
- msub
- water surface
- light
- laser
- mrow
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/86—Combinations of sonar systems with lidar systems; Combinations of sonar systems with systems not using wave reflection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
基于激光多点相干探测的水下发声目标位置估计方法及实现该方法的装置，涉及水下发声目标位置探测技术。为了解决“激光‑声”联合探测手段无法实现水下发声目标位置探测的问题。本发明首先搭建激光相干探测***，然后利用相干探测***对感兴趣的水域多个位置进行水表面声波探测，利用频谱分析、相位解调等方法实现各个点水表面声波波幅的测量，最后根据水表面声波波幅分布特征，波幅最大的位置视为水下声源中心的最优估计。本发明所述的方法和装置能够实现水下发声源位置的空中探测，可根据实际需求选择同步多点相干探测或异步多点相干探测两种方法，灵活性和机动性非常好，适用于水下发声目标探测，以及空对潜通信。
Description
技术领域
本发明涉及水下发声目标位置探测技术。
背景技术
包括水下测音器、声呐浮标及拖拽式声波***在内的声呐技术依然是当前水下发声目标探测的主流技术手段，例如利用声呐技术对MH370黑匣子的搜索，声呐技术是一种利用声波在水下的传播特性，通过电声转换和信息处理完成水下探测和通讯任务的电子设备，分为主动式和被动式两种类型。船载声呐技术由于机动性、灵活性以及隐蔽性较差，对水下发声目标的探测应用受到了诸多限制。
人们希望能够实现水下发声目标的遥测，于是提出了“激光-声”联合探测手段，目前的实验室研究显示了“激光-声”联合探测技术的潜力，激光相干探测法已被研究应用到水下声信号的提取中，取得了一定的成效，但目前“激光-声”联合探测手段仍然无法实现水下发声目标位置的探测。
发明内容
本发明的目的是为了解决“激光-声”联合探测手段无法实现水下发声目标位置探测的问题，提供两种基于激光多点相干探测的水下发声目标位置估计方法及实现上述两种方法的装置。
本发明所述的第一种基于激光多点相干探测的水下发声目标位置估计方法为同步多点探测法，该方法包括以下步骤：
步骤一、在水表面上任意选取四个点作为探测点，且四个探测点构成正方形；
步骤二、采用水下声信号激光相干探测法对四个探测点进行探测，分别得到四个探测点的探测信号Id(t)，即光电探测器接收到的信号，并对四个探测点的探测信号进行频谱分析或相位解调，得到四个探测点的水表面声波振幅值；
步骤三、判断四个探测点的水表面声波振幅值是否相同，如果相同，则将探测中心、即四个探测点构成的正方形的中心点O作为声源中心位置的最佳估计，如果不相同，则执行步骤四；
步骤四、根据振幅最大和次大探测点计算探测***光学探头移动的方向向量
其中，Amax1为四个探测点的水表面声波最大振幅值，Amax1对应的位置为探测点Smax1，Amax2为四个探测点的水表面声波次大振幅值，Amax2对应的位置为探测点Smax2，O为四个探测点构成的正方形的中心；
步骤五、按照步骤四得到的向量来移动四个光学探头，然后返回步骤二。
实现上述方法的装置包括激光器(1)、光隔离器(2)、1×4光纤耦合器(3)、四个光纤声光调制器(4)、四个光环形器(5)、四个探头(6)、四个光衰减器(7)、四个2×1耦合器(8)、4通道的光电探测器(9)和4通道的信号解调模块(10)；
激光器(1)发出的激光经过光隔离器(2)后进入1×4光纤耦合器(3)，1×4光纤耦合器(3)的四个输出端分别连接四个光纤声光调制器(4)的输入端；
每个光纤声光调制器(4)与一个光环形器(5)、一个光衰减器(7)、一个2×1耦合器(8)和一个探头(6)相对应，每个光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端，该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口，该光环形器(5)的二号端口连接探头(6)的光纤接口，该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端；
四个2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端，光电探测器(9)的输出端连接信号解调模块(10)的信号输入端，信号解调模块(10)的输出端用于连接上位机。
本发明所述的第二种基于激光多点相干探测的水下发声目标位置估计方法为异步多点探测方法，该方法为：采用水下声信号激光相干探测法对待探测水域进行逐点探测，得到各个探测点的探测信号Id(t)，即光电探测器接收到的信号，并对每个探测点的探测信号进行频谱分析或相位解调，得到每个探测点的水表面声波振幅值，将水表面声波振幅值最大的探测点的位置作为声源中心位置的最佳估计。
实现上述方法的装置包括激光器(1)、光隔离器(2)、光纤声光调制器(4)、光环形器(5)、探头(6)、光衰减器(7)、2×1耦合器(8)、光电探测器(9)和信号解调模块(10)；
激光器(1)发出的激光经过光隔离器(2)后进入光纤声光调制器(4)的输入端，光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端，该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口，该光环形器(5)的二号端口连接探头(6)的光纤接口，该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端，2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端，光电探测器(9)的输出端连接信号解调模块(10)的信号输入端，信号解调模块(10)的输出端用于连接上位机。
本发明针对水下发声源位置的空中探测这一问题提出了解决方案。本发明所述的两种方法可简述为如下过程：1)搭建激光相干探测***；2)利用相干探测***对感兴趣的水域多个位置进行水表面声波探测；3)利用频谱分析、相位解调等方法实现各个点水表面声波波幅的测量；4)根据水表面声波波幅分布特征，波幅最大的位置视为水下声源中心的最优估计。
根据基本多点相干探测的水下发声目标位置评估方法的基本原理，其实现方法可采用两种形式：1)同步多点探测；2)异步多点探测。
本发明所述的方法和装置能够实现水下发声源位置的空中探测，可根据实际需求选择同步多点相干探测或异步多点相干探测两种方法，灵活性和机动性非常好，适用于水下发声目标探测，以及空对潜通信。
附图说明
图1为实施方式一中的水表面声波的空间形态图；
图2为实施方式一中的水下声信号激光相干探测的基本原理示意图，11表示分光镜，12表示平面反射镜，13表示信号发生器，14表示功率放大器，15表示水下声源，16表示数据采集及处理***；
图3为实施方式一中四个探测点的初始位置分布图；
图4为实施方式一中声源中心定位示意图；
图5为实施方式一所述方法的流程图；
图6为实施方式二中的根据相干信号获得水表面声波信息的原理框图；
图7为实施方式五所述装置的结构示意图；
图8为实施方式七所述装置的结构示意图。
具体实施方式
具体实施方式一：结合图1至图5说明本实施方式，本实施方式所述的基于激光多点相干探测的水下发声目标位置估计方法为同步多点探测法，该方法为：
步骤一、在水表面上任意选取四个点作为探测点，且四个探测点构成正方形；
步骤二、采用水下声信号激光相干探测法对四个探测点进行探测，分别得到四个探测点的探测信号Id(t)，即光电探测器接收到的信号，并对四个探测点的探测信号进行频谱分析或相位解调，得到四个探测点的水表面声波振幅值；
步骤三、判断四个探测点的水表面声波振幅值是否相同，如果相同，则将探测中心、即四个探测点构成的正方形的中心点O作为声源中心位置的最佳估计，如果不相同，则执行步骤四；
步骤四、根据振幅最大和次大探测点计算探测***光学探头移动的方向向量
其中，Amax1为四个探测点的水表面声波最大振幅值，Amax1对应的位置为探测点Smax1，Amax2为四个探测点的水表面声波次大振幅值，Amax2对应的位置为探测点Smax2，O为四个探测点构成的正方形的中心；
步骤五、按照步骤四得到的向量来移动四个光学探头，然后返回步骤二。
水表面声波是由水下声辐射穿透水气界面而引起的一种横向微幅波，若以水表面声源扰动中心为坐标原点，水面法方向作为Z坐标轴建立坐标系，那么根据液体表面横向微幅波理论，水表面声波的三维行波解可由下式表示：
式中，η为位置(x,y)处水表面声波实时波幅，j为虚数符号，x为X轴坐标分量，y为Y轴坐标分量，ω为水表面声波的角频率(也即水下声信号的角频率)，e为自然对数的底。
由于水的粘滞作用，水表面横向微幅波的波幅不可能保持不变，也即水表面横向微幅波的范围不是无限的，水表面微幅波的振幅随着传播距离的增加而减小，即水表面微幅波区域的质点越靠近扰动源，其振幅越大。研究表明，水表面横向微幅波的波幅将以指数形式衰减，因此，可将上式修正为：
式中，τ为衰减系数。因此，由水下声辐射引起的水表面声波的空间形态如图1所示。整个波形呈现同心圆的结构，且中心幅度大，外沿幅度小。同心圆的圆心位置即为水下声辐射中心的所在位置。正是由于水表面横向微幅波的振幅衰减的特点，我们可以利用多点波幅探测的方法来实现水下目标发声位置的估计。
水下声信号激光相干探测的原理如图2所示，该图为典型的激光相干法探测水下声信号的光路图。激光器发出的激光束经声光调制器(AOM)衍射后，分成两束光，一束作为测量光，另一束作为参考光。测量光直接入射到被测水面，该水面的振动因为受到了水下声源的调制而包含有水下声源同频率的振动成分，位相被这些振动调制的反射光(实际表现为中心光强大四周光强小的光斑)入射至分光镜，并经分光镜反射；参考光经45度平面反射镜反射后，入射至分光镜，并经分光镜透射，在分光镜里内部，测量光与参考光汇合，发生了相干，相干光入射至光电探测器的探测面上。光电探测器接收相干信号，并将其转化为电信号，该电信号经过信号处理电路和数据采集及处理***处理后，被上位机读取，被上位机读取的相干信号做进一步处理可解调出水下声信号的特征。
同步多点相干探测法的4个探测点S1、S2、S3、S4的分布如图3所示，整体呈现正方形，其几何中心为O点。实践时，利用同步多点相干探测***对感兴趣的水域进行多点探测，得到4个探测点S1、S2、S3、S4处的水表面声波信息，选出水表面声波振幅最大和次大的两个探测位置，分别记为Smax1和Smax2，那么声源中心方向则位于振幅最大探测点Smax1和次大探测点Smax2与中心点O组成的张角范围内。以图3为例，水表面声波振幅最大的探测点为S1，次大探测点为S4，那么声源中心C相对于探测***的方向则位于向量和之间.为搜索水表面声波的最大振幅位置，探测***需按照方向移动，不断调整探测***的位置，使4个探测点探测到的水表面声波的振幅大小基本一致，如图4所示，此时认为探测中心点O即为声源中心的最优估计。
具体实施方式二：结合图6说明本实施方式，本实施方式是对实施方式一所述的一种基于激光多点相干探测的水下发声目标位置估计方法的进一步限定，根据实施方式一中的测量原理，步骤二中的探测信号Id(t)的获得方法为：
参考光的光振幅Eb随时间t的分布为：
Eb(t)＝Absin(ω0t+kzb+φ) (1)
式中，Ab为参考光光振幅Eb的幅值；
ω0为激光角频率；
k为激光波数；
zb为参考光的光程；
φ为激光初始相位。
测量光的光振幅Em随时间t的分布为：
Em(t)＝Amsin[(ω0+ωd)t+kzm+φ] (2)
式中，Am为测量光光振幅Eb的幅值；
zm为测量光的光程；
ωd为声光调制器调制频率。
参考光与测量光汇合后的合振幅为：
E(t)＝Eb(t)+Em(t)＝Am sin(ω0t+kzm+φ)+Ab sin[(ω0+ωd)t+kzb+φ] (3)
则相干光光强I(t)为：
光强由光电探测器接收，光电探测器不对高频项响应，即将Am 2cos2(ω0t+kzm+φ)、Ab 2cos2[(ω0+ωd)t+kzb+φ]、和AmAbcos[(2ω0+ωd)t+kzm+kzb+2φ]转变为直流信号，因此，去除直流分量后的光强信号Id表示为：
Id(t)＝KAmAbcos[k(zm-zb)]＝Adcos[ωdt+k(zm-zb)] (5)
式中，K为光电接收器的放大倍率，记Ad＝KAmAb，
由上式可知，光电探测器接收的信号是关于参考光和测量光光程差的调制信号，两束激光的光程差由两个臂的初始长度和水表面振动有关，其中水表面振动包含自然水表面的波动以及水下声信号激发的水表面声波两部分，因此，光程差zm-zb表示为：
zm-zb＝2(L1-L2)+2Ansin(2πfnt+φn)+2Assin(2πfst+φs) (6)
式中，Ansin(2πfnt+φn)为自然水表面的振动表达式(An为自然水表面振幅，fn为自然水表面振动频率，φn为自然水表面振动初相位)，Assin(2πfst+φs)为水下声源引起的水表面振动表达式(As为水表面声波振幅，fs为水表面声波振动频率，φs为水表面声波振动初相位)，L1为水面静止时测量光的初始光程，L2为水面静止时参考光的光程，那么，(L1-L2)为水面静止时两路光的初始光程差，记为ΔL，则光电探测器接收到的光强信号可由下式表述：
Id(t)＝Adcos{ωdt+k[2△L+2Ansin(2πfnt+φn)+2Assin(2πfst+φs)]} (7)
Assin(2πfst+φs)为由水下声源引起的水表面声波的信息。
由公式(7)可知，由光电探测器获得的探测信号的相位中包含了水下声源引起的水表面声波的信息：Assin(2πfst+φs)。采用常规的方法(例如频谱分析或相位解调等方法)对探测信号进行分析，即可得到水表面声波的振动幅值As和振动频率fs。可采用图6所示的原理得到水表面声波信息：Assin(2πfst+φs)。
具体实施方式三：本实施方式所述的基于激光多点相干探测的水下发声目标位置估计方法为异步多点探测方法，该方法为：采用水下声信号激光相干探测法对待探测水域进行逐点探测，得到各个探测点的探测信号Id(t)，即光电探测器接收到的信号，并对每个探测点的探测信号进行频谱分析或相位解调，得到每个探测点的水表面声波振幅值，将水表面声波振幅值最大的探测点的位置作为声源中心位置的最佳估计。
具体实施方式四：本实施方式是对实施方式三所述的基于激光多点相干探测的水下发声目标位置估计方法的进一步限定，探测信号Id(t)的获得方法为：
参考光的光振幅Eb随时间t的分布为：
Eb(t)＝Absin(ω0t+kzb+φ) (1)
式中，Ab为参考光光振幅Eb的幅值；
ω0为激光角频率；
k为激光波数；
zb为参考光的光程；
φ为激光初始相位。
测量光的光振幅Em随时间t的分布为：
Em(t)＝Amsin[(ω0+ωd)t+kzm+φ] (2)
式中，Am为测量光光振幅Eb的幅值；
zm为测量光的光程；
ωd为声光调制器调制频率。
参考光与测量光汇合后的合振幅为：
E(t)＝Eb(t)+Em(t)＝Am sin(ω0t+kzm+φ)+Absin[(ω0+ωd)t+kzb+φ] (3)
则相干光光强I(t)为：
光强由光电探测器接收，光电探测器将Am 2cos2(ω0t+kzm+φ)、Ab 2cos2[(ω0+ωd)t+kzb+φ]、和AmAbcos[(2ω0+ωd)t+kzm+kzb+2φ]转变为直流信号，因此，去除直流分量后的光强信号Id表示为：
Id(t)＝KAmAbcos[k(zm-zb)]＝Adcos[ωdt+k(zm-zb)] (5)
式中，K为光电接收器的放大倍率，记Ad＝KAmAb，
光程差zm-zb表示为：
zm-zb＝2(L1-L2)+2An sin(2πfnt+φn)+2As sin(2πfst+φs) (6)
式中，An sin(2πfnt+φn)为自然水表面的振动表达式(An为自然水表面振幅，fn为自然水表面振动频率，φn为自然水表面振动初相位)，As sin(2πfst+φs)为水下声源引起的水表面振动表达式(As为水表面声波振幅，fs为水表面声波振动频率，φs为水表面声波振动初相位)，L1为水面静止时测量光的初始光程，L2为水面静止时参考光的光程，那么，(L1-L2)为水面静止时两路光的初始光程差，记为ΔL，则光电探测器接收到的光强信号可由下式表述：
Id(t)＝Ad cos{ωdt+k[2△L+2An sin(2πfnt+φn)+2As sin(2πfst+φs)]} (7)
As sin(2πfst+φs)为由水下声源引起的水表面声波的信息。
具体实施方式五：结合图6和图7说明本实施方式，本实施方式是实现实施方式二所述的基于激光多点相干探测的水下发声目标位置估计方法的装置，所述装置包括激光器(1)、光隔离器(2)、1×4光纤耦合器(3)、四个光纤声光调制器(4)、四个光环形器(5)、四个探头(6)、四个光衰减器(7)、四个2×1耦合器(8)、4通道的光电探测器(9)和4通道的信号解调模块(10)；
激光器(1)发出的激光经过光隔离器(2)后进入1×4光纤耦合器(3)，1×4光纤耦合器(3)的四个输出端分别连接四个光纤声光调制器(4)的输入端；
每个光纤声光调制器(4)与一个光环形器(5)、一个光衰减器(7)、一个2×1耦合器(8)和一个探头(6)相对应，每个光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端，该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口，该光环形器(5)的二号端口连接探头(6)的光纤接口，该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端；
四个2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端，光电探测器(9)的输出端连接信号解调模块(10)的信号输入端，信号解调模块(10)的输出端用于连接上位机。
图7所示为实现同步多点探测法的装置的结构示意图：带尾纤的激光器(1)输出单频激光束，经过光隔离器(2)后输入到1×4光纤耦合器(3)内，该1×4光纤耦合器(3)按等光强输出4路激光，4路激光各自输入到光纤声光调制器(4)，光纤声光调制器(4)输出的0级光作为参考光，其输出的1级光作为探测光：4路探测光经光环形器(5)输出到4个探头(6)，4个探头(6)将探测激光出射到待探测水域，并接收水面返射回的激光束，即反射激光束，反射激光束的相位受到了水表面波动的调制，反射激光束经环形器后输入到2×1耦合器(8)；另一方面，参考光经过光衰减器(7)后也输入到2×1耦合器(8)，两束光汇合后产生了相干信号，4路相干信号由4通道的光电探测器(9)接收，光电探测器(9)将相干信号转换为电信号，进入4通道的信号解调模块(10)，信号解调模块(10)采用图6所示的原理对接收到的信号进行处理，得到水表面声波信息As sin(2πfst+φs)，将该信息发送至上位机，上位机即可获得4个探测点水表面声波的信息。
具体实施方式六：结合图7说明本实施方式，本实施方式是对实施方式五所述装置的进一步限定，本实施方式中，所述探头(6)包括多镜片光耦合器和定焦透镜，多镜片光耦合器的光纤接口作为探头(6)的光纤接口，探测光依次经过多镜片光耦合器和定焦透镜后入射到水表面，经水表面返回的激光依次经过定焦透镜和多镜片光耦合器后进入光环形器(5)的二号端口。
该光纤耦合器用于将探测激光出射到待探测水域，同时配合定焦透镜实现反射激光束的接收。
具体实施方式七：结合图8说明本实施方式，本实施方式是实现实施方式四所述的基于激光多点相干探测的水下发声目标位置估计方法的装置，所述装置包括激光器(1)、光隔离器(2)、光纤声光调制器(4)、光环形器(5)、探头(6)、光衰减器(7)、2×1耦合器(8)、光电探测器(9)和信号解调模块(10)；
激光器(1)发出的激光经过光隔离器(2)后进入光纤声光调制器(4)的输入端，光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端，该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口，该光环形器(5)的二号端口连接探头(6)的光纤接口，该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端，2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端，光电探测器(9)的输出端连接信号解调模块(10)的信号输入端，信号解调模块(10)的输出端用于连接上位机。
该装置利用一个光学探头对感兴趣的水域内进行逐点的水表面声波探测，以最大振幅处为声源中心位置的最优估计。
具体实施方式八：结合图8说明本实施方式，本实施方式是对实施方式七所述装置的进一步限定，本实施方式中，所述探头(6)包括多镜片光耦合器和定焦透镜，多镜片光耦合器的光纤接口作为探头(6)的光纤接口，探测光依次经过多镜片光耦合器和定焦透镜后入射到水表面，经水表面返回的激光依次经过定焦透镜和多镜片光耦合器后进入光环形器(5)的二号端口。
该光纤耦合器用于将探测激光出射到待探测水域，同时配合定焦透镜实现反射激光束的接收。
Claims (4)
1.基于激光多点相干探测的水下发声目标位置估计方法，其特征在于，该方法为：
步骤一、在水表面上任意选取四个点作为探测点，且四个探测点构成正方形；
步骤二、采用水下声信号激光相干探测法对四个探测点进行探测，分别得到四个探测点的探测信号Id(t)，即光电探测器接收到的信号，并对四个探测点的探测信号进行频谱分析或相位解调，得到四个探测点的水表面声波振幅值；
步骤三、判断四个探测点的水表面声波振幅值是否相同，如果相同，则将探测中心、即四个探测点构成的正方形的中心点O作为声源中心位置的最佳估计，如果不相同，则执行步骤四；
步骤四、根据振幅最大和次大探测点计算探测***光学探头移动的方向向量
其中，Amax1为四个探测点的水表面声波最大振幅值，Amax1对应的位置为探测点Smax1，Amax2为四个探测点的水表面声波次大振幅值，Amax2对应的位置为探测点Smax2，O为四个探测点构成的正方形的中心；
步骤五、按照步骤四得到的向量来移动四个光学探头，然后返回步骤二。
2.根据权利要求1所述的基于激光多点相干探测的水下发声目标位置估计方法，其特征在于，步骤二中的探测信号Id(t)的获得方法为：
参考光的光振幅Eb随时间t的分布为：
Eb(t)＝Absin(ω0t+kzb+φ) (1)
式中，Ab为参考光光振幅Eb的幅值；
ω0为激光角频率；
k为激光波数；
zb为参考光的光程；
φ为激光初始相位；
测量光的光振幅Em随时间t的分布为：
Em(t)＝Amsin[(ω0+ωd)t+kzm+φ] (2)
式中，Am为测量光光振幅Eb的幅值；
zm为测量光的光程；
ωd为声光调制器调制频率；
参考光与测量光汇合后的合振幅为：
E(t)＝Eb(t)+Em(t)＝Amsin(ω0t+kzm+φ)+Absin[(ω0+ωd)t+kzb+φ] (3)
则相干光光强I(t)为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>=</mo>
<msup>
<msub>
<mi>A</mi>
<mi>m</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>kz</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<msub>
<mi>A</mi>
<mi>b</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>kz</mi>
<mi>b</mi>
</msub>
<mo>+</mo>
<mi>&phi;</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>m</mi>
</msub>
<msub>
<mi>A</mi>
<mi>b</mi>
</msub>
<mi>cos</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>kz</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>kz</mi>
<mi>b</mi>
</msub>
<mo>+</mo>
<mn>2</mn>
<mi>&phi;</mi>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>m</mi>
</msub>
<msub>
<mi>A</mi>
<mi>b</mi>
</msub>
<mi>cos</mi>
<mo>&lsqb;</mo>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mi>t</mi>
<mo>+</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>b</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
光强由光电探测器接收，光电探测器将Am 2cos2(ω0t+kzm+φ)、Ab 2cos2[(ω0+ωd)t+kzb+φ]、和AmAbcos[(2ω0+ωd)t+kzm+kzb+2φ]转变为直流信号，因此，去除直流分量后的光强信号Id表示为：
Id(t)＝KAmAbcos[k(zm-zb)]＝Adcos[ωdt+k(zm-zb)] (5)
式中，K为光电接收器的放大倍率，记Ad＝KAmAb，
光程差zm-zb表示为：
zm-zb＝2(L1-L2)+2Ansin(2πfnt+φn)+2Assin(2πfst+φs) (6)
式中，Ansin(2πfnt+φn)为自然水表面的振动表达式，An为自然水表面振幅，fn为自然水表面振动频率，φn为自然水表面振动初相位，Assin(2πfst+φs)为水下声源引起的水表面振动表达式，As为水表面声波振幅，fs为水表面声波振动频率，φs为水表面声波振动初相位，L1为水面静止时测量光的初始光程，L2为水面静止时参考光的光程，那么，(L1-L2)为水面静止时两路光的初始光程差，记为ΔL，则光电探测器接收到的光强信号可由下式表述：
Id(t)＝Adcos{ωdt+k[2ΔL+2Ansin(2πfnt+φn)+2Assin(2πfst+φs)]} (7)
Assin(2πfst+φs)为由水下声源引起的水表面声波的信息。
3.实现权利要求2所述的基于激光多点相干探测的水下发声目标位置估计方法的装置，其特征在于，所述装置包括激光器(1)、光隔离器(2)、1×4光纤耦合器(3)、四个光纤声光调制器(4)、四个光环形器(5)、四个探头(6)、四个光衰减器(7)、四个2×1耦合器(8)、4通道的光电探测器(9)和4通道的信号解调模块(10)；
激光器(1)发出的激光经过光隔离器(2)后进入1×4光纤耦合器(3)，1×4光纤耦合器(3)的四个输出端分别连接四个光纤声光调制器(4)的输入端；
每个光纤声光调制器(4)与一个光环形器(5)、一个光衰减器(7)、一个2×1耦合器(8)和一个探头(6)相对应，每个光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端，该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口，该光环形器(5)的二号端口连接探头(6)的光纤接口，该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端；
四个2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端，光电探测器(9)的输出端连接信号解调模块(10)的信号输入端，信号解调模块(10)的输出端用于连接上位机。
4.根据权利要求3所述的实现基于激光多点相干探测的水下发声目标位置估计方法的装置，其特征在于，所述探头(6)包括多镜片光耦合器和定焦透镜，多镜片光耦合器的光纤接口作为探头(6)的光纤接口，探测光依次经过多镜片光耦合器和定焦透镜后入射到水表面，经水表面返回的激光依次经过定焦透镜和多镜片光耦合器后进入光环形器(5)的二号端口。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511021900.5A CN105487077B (zh) | 2015-12-29 | 2015-12-29 | 基于激光多点相干探测的水下发声目标位置估计方法及实现该方法的装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511021900.5A CN105487077B (zh) | 2015-12-29 | 2015-12-29 | 基于激光多点相干探测的水下发声目标位置估计方法及实现该方法的装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105487077A CN105487077A (zh) | 2016-04-13 |
| CN105487077B true CN105487077B (zh) | 2017-11-03 |
Family
ID=55674184
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511021900.5A Active CN105487077B (zh) | 2015-12-29 | 2015-12-29 | 基于激光多点相干探测的水下发声目标位置估计方法及实现该方法的装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105487077B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106546991A (zh) * | 2016-11-02 | 2017-03-29 | 哈尔滨工业大学（威海） | 水下发声目标空中探测装置及方法 |
| CN106595836B (zh) * | 2016-12-26 | 2019-04-26 | 哈尔滨工业大学 | 水下双声源频谱混叠情况下的频率提取方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6396770B1 (en) * | 1982-06-28 | 2002-05-28 | Bae Systems Information And Electronic Systems Integration Inc. | Steerable thermoacoustic array |
| CN101082671A (zh) * | 2007-06-27 | 2007-12-05 | 中国科学院上海光学精密机械研究所 | 相干激光遥感探测水下声信号的方法和装置 |
| CN104316160A (zh) * | 2014-11-20 | 2015-01-28 | 哈尔滨工业大学 | 基于小波脊的水下声信号瞬时频率解调方法 |
| CN104808208A (zh) * | 2015-04-16 | 2015-07-29 | 浙江大学 | 一种基于激光声源探测水下目标方位及尺寸的测量***及其测量方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8228760B2 (en) * | 2008-02-12 | 2012-07-24 | The United States Of America, As Represented By The Secretary Of The Navy | Airborne laser-acoustic mine detection system |
-
2015
- 2015-12-29 CN CN201511021900.5A patent/CN105487077B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6396770B1 (en) * | 1982-06-28 | 2002-05-28 | Bae Systems Information And Electronic Systems Integration Inc. | Steerable thermoacoustic array |
| CN101082671A (zh) * | 2007-06-27 | 2007-12-05 | 中国科学院上海光学精密机械研究所 | 相干激光遥感探测水下声信号的方法和装置 |
| CN104316160A (zh) * | 2014-11-20 | 2015-01-28 | 哈尔滨工业大学 | 基于小波脊的水下声信号瞬时频率解调方法 |
| CN104808208A (zh) * | 2015-04-16 | 2015-07-29 | 浙江大学 | 一种基于激光声源探测水下目标方位及尺寸的测量***及其测量方法 |
Non-Patent Citations (5)
| Title |
|---|
| Amplitude measurement of weak sinusoidal water surface acoustic wave using laser interferometer;Lieshan Zhang et al;《CHINESE OPTICS LETTERS》;20150910;第13卷(第9期);第1-5页 * |
| 基于小波脊的水表面声波激光干涉探测;张烈山 等;《光电子激光》;20150131;第26卷(第1期);第103-107页 * |
| 基于激光干涉的水表面声波探测技术研究;张晓琳;《中国博士学位论文全文数据库 信息科技辑》;20110815(第8期);第I135-22页 * |
| 水下声信号的激光干涉测量;张晓琳 等;《光学精密工程》;20100430;第18卷(第4期);第809-815页 * |
| 激光干涉法探测水下声信号;张晓琳 等;《大连海事大学学报》;20090228;第35卷(第1期);第53-56页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105487077A (zh) | 2016-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104808208B (zh) | 一种基于激光声源探测水下目标方位及尺寸的测量***及其测量方法 | |
| CN105092014B (zh) | 基于波束形成的分布式光纤声波探测装置及探测方法 | |
| CN104181505B (zh) | 一种基于近场源定位算法的多目标水声定位方法和*** | |
| CN108036852A (zh) | 一种光纤声传感器及多点声波检测装置 | |
| CN103954310B (zh) | 一种干涉型光纤传感器的大动态信号解调装置及解调方法 | |
| CN105606194A (zh) | 一种基于激光正交偏振干涉技术的水下声信号实时提取方法 | |
| CN102401691B (zh) | 全光纤激光多普勒三维测振仪 | |
| WO2020140869A1 (zh) | 一种扰动源多维空间定位***和方法 | |
| CN103047540A (zh) | 基于光纤传感的天然气管道泄漏监测光路*** | |
| CN108415067A (zh) | 一种基于微结构光纤分布式声波传感的地震波测量*** | |
| MX2011009600A (es) | Determinacion de categorias teoricas. | |
| CN106813766B (zh) | 声磁同测的分布式光纤传感*** | |
| CN105487077B (zh) | 基于激光多点相干探测的水下发声目标位置估计方法及实现该方法的装置 | |
| US9453900B2 (en) | Method and apparatus for three dimensional wavenumber-frequency analysis | |
| CN203929213U (zh) | 一种多点式激光多普勒测振*** | |
| Hwang et al. | Position estimation of sound source using three optical Mach-Zehnder acoustic sensor array | |
| Ballard et al. | Measurements and modeling of acoustic propagation in a scale model canyon | |
| CN203164417U (zh) | 基于多普勒效应的水下目标定位装置 | |
| Li et al. | Phase-shifted sensitivity calibration of fiber optic vector hydrophone based on heterodyne method | |
| CN203574661U (zh) | 一种基于光纤传感技术的水下对讲*** | |
| CN102997047A (zh) | 一种天然气管道泄漏光纤监测传感器的复用/解复用*** | |
| CN111128217B (zh) | 分布式多通道相干激光雷达语音侦听方法及装置 | |
| CN103226198A (zh) | 基于多普勒效应的水下目标定位方法和装置 | |
| Dahl et al. | Active intensity vortex and stagnation point singularities in a shallow underwater waveguide | |
| Pang et al. | Research of very low-frequency (VLF) hydroacoustic detection based on Fizeau interferometry |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |