[go: up one dir, main page]

CN105487077B - Underwater sounding target position estimation method based on laser multipoint coherent detection and device for realizing method - Google Patents

Underwater sounding target position estimation method based on laser multipoint coherent detection and device for realizing method Download PDF

Info

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
detection
water surface
optical
laser
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
Application number
CN201511021900.5A
Other languages
Chinese (zh)
Other versions
CN105487077A (en
Inventor
张晓琳
张烈山
唐文彦
王军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology Shenzhen
Original Assignee
Harbin Institute of Technology Shenzhen
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology Shenzhen filed Critical Harbin Institute of Technology Shenzhen
Priority to CN201511021900.5A priority Critical patent/CN105487077B/en
Publication of CN105487077A publication Critical patent/CN105487077A/en
Application granted granted Critical
Publication of CN105487077B publication Critical patent/CN105487077B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/86Combinations of sonar systems with lidar systems; Combinations of sonar systems with systems not using wave reflection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/06Systems determining the position data of a target
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment 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)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

An underwater sounding target position estimation method based on laser multipoint coherent detection and a device for realizing the method relate to the underwater sounding target position detection technology. The underwater sounding target position detection method aims at solving the problem that the underwater sounding target position detection cannot be realized by a laser-sound combined detection means. The method comprises the steps of firstly building a laser coherent detection system, then carrying out water surface acoustic wave detection on a plurality of positions of an interested water area by using the coherent detection system, measuring the water surface acoustic wave amplitude of each point by using methods such as frequency spectrum analysis, phase demodulation and the like, and finally considering the position with the maximum amplitude as the optimal estimation of the center of an underwater sound source according to the distribution characteristics of the water surface acoustic wave amplitude. The method and the device can realize the aerial detection of the position of the underwater sound generating source, can select two methods of synchronous multipoint coherent detection or asynchronous multipoint coherent detection according to actual requirements, have very good flexibility and maneuverability, and are suitable for the detection of the underwater sound generating target and the aerial-to-submarine communication.

Description

基于激光多点相干探测的水下发声目标位置估计方法及实现 该方法的装置Method and implementation of underwater sound target position estimation based on laser multi-point coherent detection The means of the method

技术领域technical field

本发明涉及水下发声目标位置探测技术。The invention relates to an underwater sound target position detection technology.

背景技术Background technique

包括水下测音器、声呐浮标及拖拽式声波定位器在内的声呐技术依然是当前水下发声目标探测的主流技术手段,例如利用声呐技术对MH370黑匣子的搜索,声呐技术是一种利用声波在水下的传播特性,通过电声转换和信息处理完成水下探测和通讯任务的电子设备,分为主动式和被动式两种类型。船载声呐技术由于机动性、灵活性以及隐蔽性较差,对水下发声目标的探测应用受到了诸多限制。Sonar technology, including underwater sounders, sonobuoys and towed sonars, is still the mainstream technology for underwater sound target detection. For example, sonar technology is used to search for the black box of MH370. The propagation characteristics of sound waves underwater, electronic equipment that completes underwater detection and communication tasks through electro-acoustic conversion and information processing, is divided into two types: active and passive. Due to the poor mobility, flexibility and concealment of shipborne sonar technology, the application of underwater acoustic target detection is limited.

人们希望能够实现水下发声目标的遥测,于是提出了“激光-声”联合探测手段,目前的实验室研究显示了“激光-声”联合探测技术的潜力,激光相干探测法已被研究应用到水下声信号的提取中,取得了一定的成效,但目前“激光-声”联合探测手段仍然无法实现水下发声目标位置的探测。People hope to be able to realize the remote measurement of underwater sound targets, so the "laser-sound" joint detection method is proposed. The current laboratory research shows the potential of the "laser-sound" joint detection technology. The laser coherent detection method has been researched and applied to Some achievements have been made in the extraction of underwater acoustic signals, but the current "laser-acoustic" joint detection method is still unable to detect the location of underwater acoustic targets.

发明内容Contents of the invention

本发明的目的是为了解决“激光-声”联合探测手段无法实现水下发声目标位置探测的问题,提供两种基于激光多点相干探测的水下发声目标位置估计方法及实现上述两种方法的装置。The purpose of the present invention is to solve the problem that the "laser-acoustic" joint detection method cannot realize the position detection of the underwater sound target, and provide two methods for estimating the position of the underwater sound target based on laser multi-point coherent detection and the method for realizing the above two methods device.

本发明所述的第一种基于激光多点相干探测的水下发声目标位置估计方法为同步多点探测法,该方法包括以下步骤:The first underwater sound target position estimation method based on laser multi-point coherent detection described in the present invention is a synchronous multi-point detection method, which includes the following steps:

步骤一、在水表面上任意选取四个点作为探测点,且四个探测点构成正方形;Step 1. Randomly select four points on the water surface as detection points, and the four detection points form a square;

步骤二、采用水下声信号激光相干探测法对四个探测点进行探测,分别得到四个探测点的探测信号Id(t),即光电探测器接收到的信号,并对四个探测点的探测信号进行频谱分析或相位解调,得到四个探测点的水表面声波振幅值;Step 2, using the underwater acoustic signal laser coherent detection method to detect the four detection points, respectively obtain the detection signals I d (t) of the four detection points, that is, the signal received by the photoelectric detector, and detect the four detection points Spectrum analysis or phase demodulation is performed on the detection signal to obtain the amplitude values of the water surface acoustic waves at the four detection points;

步骤三、判断四个探测点的水表面声波振幅值是否相同,如果相同,则将探测中心、即四个探测点构成的正方形的中心点O作为声源中心位置的最佳估计,如果不相同,则执行步骤四;Step 3. Determine whether the amplitude values of the water surface acoustic waves at the four detection points are the same. If they are the same, use the detection center, that is, the center point O of the square formed by the four detection points, as the best estimate of the center position of the sound source. If they are not the same , execute step 4;

步骤四、根据振幅最大和次大探测点计算探测系统光学探头移动的方向向量 Step 4. Calculate the moving direction vector of the optical probe of the detection system according to the maximum and second maximum amplitude detection points

其中,Amax1为四个探测点的水表面声波最大振幅值,Amax1对应的位置为探测点Smax1,Amax2为四个探测点的水表面声波次大振幅值,Amax2对应的位置为探测点Smax2,O为四个探测点构成的正方形的中心;Among them, A max1 is the maximum amplitude value of the water surface acoustic wave at the four detection points, the position corresponding to A max1 is the detection point S max1 , A max2 is the second largest amplitude value of the water surface acoustic wave at the four detection points, and the corresponding position of A max2 is Detection point S max2 , O is the center of a square formed by four detection points;

步骤五、按照步骤四得到的向量来移动四个光学探头,然后返回步骤二。Step 5. According to the vector obtained in step 4 to move the four optical probes, then return to step two.

实现上述方法的装置包括激光器(1)、光隔离器(2)、1×4光纤耦合器(3)、四个光纤声光调制器(4)、四个光环形器(5)、四个探头(6)、四个光衰减器(7)、四个2×1耦合器(8)、4通道的光电探测器(9)和4通道的信号解调模块(10);The device for realizing the above method includes a laser (1), an optical isolator (2), a 1×4 fiber coupler (3), four fiber optic acousto-optic modulators (4), four optical circulators (5), four Probes (6), four optical attenuators (7), four 2×1 couplers (8), 4-channel photodetectors (9) and 4-channel signal demodulation modules (10);

激光器(1)发出的激光经过光隔离器(2)后进入1×4光纤耦合器(3),1×4光纤耦合器(3)的四个输出端分别连接四个光纤声光调制器(4)的输入端;The laser light emitted by the laser (1) enters the 1×4 fiber coupler (3) after passing through the optical isolator (2), and the four output ports of the 1×4 fiber coupler (3) are respectively connected to four fiber optic acousto-optic modulators ( 4) input terminal;

每个光纤声光调制器(4)与一个光环形器(5)、一个光衰减器(7)、一个2×1耦合器(8)和一个探头(6)相对应,每个光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端,该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口,该光环形器(5)的二号端口连接探头(6)的光纤接口,该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端;Each optical fiber acousto-optic modulator (4) corresponds to an optical circulator (5), an optical attenuator (7), a 2×1 coupler (8) and a probe (6), and each optical fiber acousto-optic An output end of the modulator (4) is connected to an input end of a 2×1 coupler (8) through an optical attenuator (7), and the other output end of the optical fiber acousto-optic modulator (4) is connected to an optical circulator (5 ), the No. 2 port of the optical circulator (5) is connected to the optical fiber interface of the probe (6), and the No. 3 port of the optical circulator (5) is connected to another input of the 2×1 coupler (8) end;

四个2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端,光电探测器(9)的输出端连接信号解调模块(10)的信号输入端,信号解调模块(10)的输出端用于连接上位机。The output ends of the four 2×1 couplers (8) are connected to the signal input ends of the photodetectors (9), and the output ends of the photodetectors (9) are connected to the signal input ends of the signal demodulation module (10), and the signal demodulation The output end of the module (10) is used for connecting to the upper computer.

本发明所述的第二种基于激光多点相干探测的水下发声目标位置估计方法为异步多点探测方法,该方法为:采用水下声信号激光相干探测法对待探测水域进行逐点探测,得到各个探测点的探测信号Id(t),即光电探测器接收到的信号,并对每个探测点的探测信号进行频谱分析或相位解调,得到每个探测点的水表面声波振幅值,将水表面声波振幅值最大的探测点的位置作为声源中心位置的最佳估计。The second underwater acoustic target position estimation method based on laser multi-point coherent detection described in the present invention is an asynchronous multi-point detection method. Obtain the detection signal I d (t) of each detection point, that is, the signal received by the photoelectric detector, and perform spectrum analysis or phase demodulation on the detection signal of each detection point to obtain the amplitude value of the water surface acoustic wave at each detection point , the position of the detection point with the largest amplitude value of the water surface acoustic wave is taken as the best estimate of the center position of the sound source.

实现上述方法的装置包括激光器(1)、光隔离器(2)、光纤声光调制器(4)、光环形器(5)、探头(6)、光衰减器(7)、2×1耦合器(8)、光电探测器(9)和信号解调模块(10);The device for realizing the above method includes a laser (1), an optical isolator (2), an optical fiber acousto-optic modulator (4), an optical circulator (5), a probe (6), an optical attenuator (7), and a 2×1 coupling device (8), photodetector (9) and signal demodulation module (10);

激光器(1)发出的激光经过光隔离器(2)后进入光纤声光调制器(4)的输入端,光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端,该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口,该光环形器(5)的二号端口连接探头(6)的光纤接口,该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端,2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端,光电探测器(9)的输出端连接信号解调模块(10)的信号输入端,信号解调模块(10)的输出端用于连接上位机。The laser light emitted by the laser (1) enters the input end of the optical fiber acousto-optic modulator (4) after passing through the optical isolator (2), and one output end of the optical fiber acousto-optic modulator (4) is connected to 2 × 1 One input end of the coupler (8), the other output end of the optical fiber acousto-optic modulator (4) is connected to the No. 1 port of the optical circulator (5), and the No. 2 port of the optical circulator (5) is connected to the probe (6), the No. 3 port of the optical circulator (5) is connected to another input end of the 2 × 1 coupler (8), and the output end of the 2 × 1 coupler (8) is connected to the photodetector (9 ), the output end of the photodetector (9) is connected to the signal input end of the signal demodulation module (10), and the output end of the signal demodulation module (10) is used to connect to the upper computer.

本发明针对水下发声源位置的空中探测这一问题提出了解决方案。本发明所述的两种方法可简述为如下过程:1)搭建激光相干探测系统;2)利用相干探测系统对感兴趣的水域多个位置进行水表面声波探测;3)利用频谱分析、相位解调等方法实现各个点水表面声波波幅的测量;4)根据水表面声波波幅分布特征,波幅最大的位置视为水下声源中心的最优估计。The present invention proposes a solution to the problem of aerial detection of the location of underwater sound sources. The two methods described in the present invention can be briefly described as the following process: 1) build a laser coherent detection system; 2) use the coherent detection system to detect water surface acoustic waves at multiple locations in the water area of interest; 3) use spectrum analysis, phase The method of demodulation realizes the measurement of the amplitude of the water surface acoustic wave at each point; 4) According to the distribution characteristics of the amplitude of the water surface acoustic wave, the position with the largest amplitude is regarded as the optimal estimate of the center of the underwater sound source.

根据基本多点相干探测的水下发声目标位置评估方法的基本原理,其实现方法可采用两种形式:1)同步多点探测;2)异步多点探测。According to the basic principle of the underwater acoustic target position assessment method of basic multipoint coherent detection, its implementation method can adopt two forms: 1) synchronous multipoint detection; 2) asynchronous multipoint detection.

本发明所述的方法和装置能够实现水下发声源位置的空中探测,可根据实际需求选择同步多点相干探测或异步多点相干探测两种方法,灵活性和机动性非常好,适用于水下发声目标探测,以及空对潜通信。The method and device of the present invention can realize the aerial detection of the location of the underwater sound source, and two methods of synchronous multi-point coherent detection or asynchronous multi-point coherent detection can be selected according to actual needs. The flexibility and maneuverability are very good, and it is suitable for underwater Down-emitting acoustic target detection, and air-to-submarine communications.

附图说明Description of drawings

图1为实施方式一中的水表面声波的空间形态图;Fig. 1 is the spatial shape diagram of the water surface acoustic wave in the first embodiment;

图2为实施方式一中的水下声信号激光相干探测的基本原理示意图,11表示分光镜,12表示平面反射镜,13表示信号发生器,14表示功率放大器,15表示水下声源,16表示数据采集及处理系统;2 is a schematic diagram of the basic principle of laser coherent detection of underwater acoustic signals in Embodiment 1. 11 represents a beam splitter, 12 represents a plane mirror, 13 represents a signal generator, 14 represents a power amplifier, 15 represents an underwater sound source, 16 Indicates the data acquisition and processing system;

图3为实施方式一中四个探测点的初始位置分布图;Fig. 3 is an initial position distribution diagram of four detection points in Embodiment 1;

图4为实施方式一中声源中心定位示意图;Fig. 4 is a schematic diagram of sound source center positioning in Embodiment 1;

图5为实施方式一所述方法的流程图;Fig. 5 is a flowchart of the method described in Embodiment 1;

图6为实施方式二中的根据相干信号获得水表面声波信息的原理框图;Fig. 6 is a functional block diagram of obtaining water surface acoustic wave information according to coherent signals in Embodiment 2;

图7为实施方式五所述装置的结构示意图;Fig. 7 is a schematic structural diagram of the device described in Embodiment 5;

图8为实施方式七所述装置的结构示意图。Fig. 8 is a schematic structural diagram of the device described in Embodiment 7.

具体实施方式detailed description

具体实施方式一:结合图1至图5说明本实施方式,本实施方式所述的基于激光多点相干探测的水下发声目标位置估计方法为同步多点探测法,该方法为:Specific Embodiment 1: This embodiment is described in conjunction with FIGS. 1 to 5. The method for estimating the position of an underwater sounding target based on laser multi-point coherent detection described in this embodiment is a synchronous multi-point detection method. The method is:

步骤一、在水表面上任意选取四个点作为探测点,且四个探测点构成正方形;Step 1. Randomly select four points on the water surface as detection points, and the four detection points form a square;

步骤二、采用水下声信号激光相干探测法对四个探测点进行探测,分别得到四个探测点的探测信号Id(t),即光电探测器接收到的信号,并对四个探测点的探测信号进行频谱分析或相位解调,得到四个探测点的水表面声波振幅值;Step 2, using the underwater acoustic signal laser coherent detection method to detect the four detection points, respectively obtain the detection signals I d (t) of the four detection points, that is, the signal received by the photoelectric detector, and detect the four detection points Spectrum analysis or phase demodulation is performed on the detection signal to obtain the amplitude values of the water surface acoustic waves at the four detection points;

步骤三、判断四个探测点的水表面声波振幅值是否相同,如果相同,则将探测中心、即四个探测点构成的正方形的中心点O作为声源中心位置的最佳估计,如果不相同,则执行步骤四;Step 3. Determine whether the amplitude values of the water surface acoustic waves at the four detection points are the same. If they are the same, use the detection center, that is, the center point O of the square formed by the four detection points, as the best estimate of the center position of the sound source. If they are not the same , execute step 4;

步骤四、根据振幅最大和次大探测点计算探测系统光学探头移动的方向向量 Step 4. Calculate the moving direction vector of the optical probe of the detection system according to the maximum and second maximum amplitude detection points

其中,Amax1为四个探测点的水表面声波最大振幅值,Amax1对应的位置为探测点Smax1,Amax2为四个探测点的水表面声波次大振幅值,Amax2对应的位置为探测点Smax2,O为四个探测点构成的正方形的中心;Among them, A max1 is the maximum amplitude value of the water surface acoustic wave at the four detection points, the position corresponding to A max1 is the detection point S max1 , A max2 is the second largest amplitude value of the water surface acoustic wave at the four detection points, and the position corresponding to A max2 is Detection point S max2 , O is the center of a square formed by four detection points;

步骤五、按照步骤四得到的向量来移动四个光学探头,然后返回步骤二。Step 5. According to the vector obtained in step 4 to move the four optical probes, then return to step two.

水表面声波是由水下声辐射穿透水气界面而引起的一种横向微幅波,若以水表面声源扰动中心为坐标原点,水面法方向作为Z坐标轴建立坐标系,那么根据液体表面横向微幅波理论,水表面声波的三维行波解可由下式表示:The water surface acoustic wave is a kind of transverse micro-amplitude wave caused by underwater acoustic radiation penetrating the water-air interface. If the center of the sound source disturbance on the water surface is taken as the coordinate origin and the normal direction of the water surface is used as the Z coordinate axis to establish a coordinate system, then according to the liquid According to the surface transverse micro-amplitude wave theory, the three-dimensional traveling wave solution of the water surface acoustic wave can be expressed by the following formula:

式中,η为位置(x,y)处水表面声波实时波幅,j为虚数符号,x为X轴坐标分量,y为Y轴坐标分量,ω为水表面声波的角频率(也即水下声信号的角频率),e为自然对数的底。In the formula, η is the real-time amplitude of the water surface acoustic wave at the position (x, y), j is the imaginary number symbol, x is the X-axis coordinate component, y is the Y-axis coordinate component, and ω is the angular frequency of the water surface acoustic wave (that is, the underwater angular frequency of the acoustic signal), and e is the base of the natural logarithm.

由于水的粘滞作用,水表面横向微幅波的波幅不可能保持不变,也即水表面横向微幅波的范围不是无限的,水表面微幅波的振幅随着传播距离的增加而减小,即水表面微幅波区域的质点越靠近扰动源,其振幅越大。研究表明,水表面横向微幅波的波幅将以指数形式衰减,因此,可将上式修正为:Due to the viscous effect of water, the amplitude of the transverse micro-amplitude wave on the water surface cannot remain unchanged, that is, the range of the transverse micro-amplitude wave on the water surface is not infinite, and the amplitude of the micro-amplitude wave on the water surface decreases with the increase of the propagation distance. Small, that is, the closer the particle in the micro-amplitude wave region of the water surface is to the disturbance source, the larger the amplitude will be. Studies have shown that the amplitude of the transverse micro-amplitude wave on the water surface will decay exponentially, so the above formula can be modified as:

式中,τ为衰减系数。因此,由水下声辐射引起的水表面声波的空间形态如图1所示。整个波形呈现同心圆的结构,且中心幅度大,外沿幅度小。同心圆的圆心位置即为水下声辐射中心的所在位置。正是由于水表面横向微幅波的振幅衰减的特点,我们可以利用多点波幅探测的方法来实现水下目标发声位置的估计。In the formula, τ is the attenuation coefficient. Therefore, the spatial shape of the water surface acoustic wave caused by underwater acoustic radiation is shown in Fig. 1. The entire waveform presents a concentric circle structure with a large center amplitude and a small outer edge amplitude. The center of the concentric circles is the location of the underwater sound radiation center. It is precisely because of the characteristics of the amplitude attenuation of the transverse micro-amplitude wave on the water surface that we can use the method of multi-point amplitude detection to realize the estimation of the sound position of the underwater target.

水下声信号激光相干探测的原理如图2所示,该图为典型的激光相干法探测水下声信号的光路图。激光器发出的激光束经声光调制器(AOM)衍射后,分成两束光,一束作为测量光,另一束作为参考光。测量光直接入射到被测水面,该水面的振动因为受到了水下声源的调制而包含有水下声源同频率的振动成分,位相被这些振动调制的反射光(实际表现为中心光强大四周光强小的光斑)入射至分光镜,并经分光镜反射;参考光经45度平面反射镜反射后,入射至分光镜,并经分光镜透射,在分光镜里内部,测量光与参考光汇合,发生了相干,相干光入射至光电探测器的探测面上。光电探测器接收相干信号,并将其转化为电信号,该电信号经过信号处理电路和数据采集及处理系统处理后,被上位机读取,被上位机读取的相干信号做进一步处理可解调出水下声信号的特征。The principle of laser coherent detection of underwater acoustic signals is shown in Figure 2, which is a typical optical path diagram of laser coherent detection of underwater acoustic signals. The laser beam emitted by the laser is diffracted by the acousto-optic modulator (AOM), and then split into two beams, one beam is used as the measurement beam, and the other beam is used as the reference beam. The measurement light is directly incident on the water surface to be measured. The vibration of the water surface contains vibration components of the same frequency as the underwater sound source because of the modulation of the underwater sound source. The light spot with small light intensity around it) enters the beam splitter and is reflected by the beam splitter; the reference light is reflected by the 45-degree flat mirror, enters the beam splitter, and is transmitted through the beam splitter. Inside the beam splitter, the measurement light and the reference light The light merges, coherence occurs, and the coherent light is incident on the detection surface of the photodetector. The photodetector receives the coherent signal and converts it into an electrical signal. After the electrical signal is processed by the signal processing circuit and the data acquisition and processing system, it is read by the host computer, and the coherent signal read by the host computer is further processed to be decipherable. Call out the characteristics of the underwater acoustic signal.

同步多点相干探测法的4个探测点S1、S2、S3、S4的分布如图3所示,整体呈现正方形,其几何中心为O点。实践时,利用同步多点相干探测系统对感兴趣的水域进行多点探测,得到4个探测点S1、S2、S3、S4处的水表面声波信息,选出水表面声波振幅最大和次大的两个探测位置,分别记为Smax1和Smax2,那么声源中心方向则位于振幅最大探测点Smax1和次大探测点Smax2与中心点O组成的张角范围内。以图3为例,水表面声波振幅最大的探测点为S1,次大探测点为S4,那么声源中心C相对于探测系统的方向则位于向量之间.为搜索水表面声波的最大振幅位置,探测系统需按照方向移动,不断调整探测系统的位置,使4个探测点探测到的水表面声波的振幅大小基本一致,如图4所示,此时认为探测中心点O即为声源中心的最优估计。The distribution of the four detection points S 1 , S 2 , S 3 , and S 4 of the synchronous multi-point coherent detection method is shown in Figure 3. The overall shape is a square, and its geometric center is point O. In practice, use the synchronous multi-point coherent detection system to conduct multi-point detection on the water area of interest, obtain the water surface acoustic wave information at the four detection points S 1 , S 2 , S 3 , and S 4 , and select the water surface acoustic wave with the largest amplitude and the second-largest two detection positions, respectively denoted as S max1 and S max2 , then the center direction of the sound source is within the angular range formed by the detection point S max1 with the largest amplitude, the second-largest detection point S max2 and the central point O. Taking Figure 3 as an example, the detection point with the largest amplitude of the acoustic wave on the water surface is S 1 , and the second largest detection point is S 4 , then the direction of the sound source center C relative to the detection system then in the vector with Between. In order to search for the maximum amplitude position of the water surface acoustic wave, the detection system needs to follow the direction Move and constantly adjust the position of the detection system so that the amplitudes of the water surface acoustic waves detected by the four detection points are basically the same, as shown in Figure 4. At this time, the detection center point O is considered to be the optimal estimate of the sound source center.

具体实施方式二:结合图6说明本实施方式,本实施方式是对实施方式一所述的一种基于激光多点相干探测的水下发声目标位置估计方法的进一步限定,根据实施方式一中的测量原理,步骤二中的探测信号Id(t)的获得方法为:Specific Embodiment 2: This embodiment is described in conjunction with FIG. 6. This embodiment is a further limitation of the underwater sound target position estimation method based on laser multi-point coherent detection described in Embodiment 1. According to Embodiment 1 Measuring principle, the acquisition method of the detection signal I d (t) in step 2 is:

参考光的光振幅Eb随时间t的分布为:The distribution of the light amplitude E b of the reference light with time t is:

Eb(t)=Absin(ω0t+kzb+φ) (1)E b (t)=A b sin(ω 0 t+kz b +φ) (1)

式中,Ab为参考光光振幅Eb的幅值;In the formula, A b is the amplitude of the reference light amplitude E b ;

ω0为激光角频率;ω 0 is the laser angular frequency;

k为激光波数;k is the laser wave number;

zb为参考光的光程;z b is the optical path of the reference light;

φ为激光初始相位。φ is the initial phase of the laser.

测量光的光振幅Em随时间t的分布为:The distribution of the light amplitude E m of the measured light over time t is:

Em(t)=Amsin[(ω0d)t+kzm+φ] (2)E m (t)=A m sin[(ω 0d )t+kz m +φ] (2)

式中,Am为测量光光振幅Eb的幅值;In the formula, A m is the amplitude of the measured light amplitude E b ;

zm为测量光的光程;z m is the optical path of the measuring light;

ωd为声光调制器调制频率。ω d is the modulation frequency of the AOM.

参考光与测量光汇合后的合振幅为:The combined amplitude of the reference light and measuring light is:

E(t)=Eb(t)+Em(t)=Am sin(ω0t+kzm+φ)+Ab sin[(ω0d)t+kzb+φ] (3)E(t)=E b (t)+E m (t)=A m sin(ω 0 t+kz m +φ)+A b sin[(ω 0d )t+kz b +φ] ( 3)

则相干光光强I(t)为:Then the coherent light intensity I(t) is:

光强由光电探测器接收,光电探测器不对高频项响应,即将Am 2cos20t+kzm+φ)、Ab 2cos2[(ω0d)t+kzb+φ]、和AmAbcos[(2ω0d)t+kzm+kzb+2φ]转变为直流信号,因此,去除直流分量后的光强信号Id表示为:The light intensity is received by the photodetector, which does not respond to high-frequency terms, that is, A m 2 cos 20 t+kz m +φ), A b 2 cos 2 [(ω 0d )t+kz b +φ], and A m A b cos[(2ω 0d )t+kz m +kz b +2φ] are transformed into DC signals, therefore, the light intensity signal I d after removing the DC component is expressed as:

Id(t)=KAmAbcos[k(zm-zb)]=Adcos[ωdt+k(zm-zb)] (5)I d (t)=KA m A b cos[k(z m -z b )]=A d cos[ω d t+k(z m -z b )] (5)

式中,K为光电接收器的放大倍率,记Ad=KAmAbIn the formula, K is the magnification of the photoelectric receiver, write A d =KA m A b ,

由上式可知,光电探测器接收的信号是关于参考光和测量光光程差的调制信号,两束激光的光程差由两个臂的初始长度和水表面振动有关,其中水表面振动包含自然水表面的波动以及水下声信号激发的水表面声波两部分,因此,光程差zm-zb表示为:It can be seen from the above formula that the signal received by the photodetector is a modulation signal about the optical path difference between the reference light and the measuring light. The optical path difference of the two laser beams is related to the initial length of the two arms and the vibration of the water surface, where the vibration of the water surface includes There are two parts of the natural water surface fluctuation and the water surface acoustic wave excited by the underwater acoustic signal. Therefore, the optical path difference z m -z b is expressed as:

zm-zb=2(L1-L2)+2Ansin(2πfnt+φn)+2Assin(2πfst+φs) (6)z m -z b =2(L 1 -L 2 )+2A n sin(2πf n t+φ n )+2A s sin(2πf s t+φ s ) (6)

式中,Ansin(2πfnt+φn)为自然水表面的振动表达式(An为自然水表面振幅,fn为自然水表面振动频率,φn为自然水表面振动初相位),Assin(2πfst+φs)为水下声源引起的水表面振动表达式(As为水表面声波振幅,fs为水表面声波振动频率,φs为水表面声波振动初相位),L1为水面静止时测量光的初始光程,L2为水面静止时参考光的光程,那么,(L1-L2)为水面静止时两路光的初始光程差,记为ΔL,则光电探测器接收到的光强信号可由下式表述:In the formula, A n sin(2πf n t+φ n ) is the vibration expression of the natural water surface (A n is the amplitude of the natural water surface, f n is the vibration frequency of the natural water surface, and φ n is the initial phase of the natural water surface vibration) , A s sin(2πf s t+φ s ) is the expression of water surface vibration caused by underwater sound source (A s is the amplitude of water surface acoustic wave, f s is the frequency of water surface acoustic wave vibration, φ s is the initial value of water surface acoustic wave vibration phase), L 1 is the initial optical path length of the measuring light when the water surface is stationary, L 2 is the optical path length of the reference light when the water surface is stationary, then (L 1 -L 2 ) is the initial optical path difference of the two paths of light when the water surface is stationary, Denoted as ΔL, the light intensity signal received by the photodetector can be expressed by the following formula:

Id(t)=Adcos{ωdt+k[2△L+2Ansin(2πfnt+φn)+2Assin(2πfst+φs)]} (7)I d (t)=A d cos{ω d t+k[2△L+2A n sin(2πf n t+φ n )+2A s sin(2πf s t+φ s )]} (7)

Assin(2πfst+φs)为由水下声源引起的水表面声波的信息。A s sin(2πf s t+φ s ) is the information of the water surface acoustic wave caused by the underwater sound source.

由公式(7)可知,由光电探测器获得的探测信号的相位中包含了水下声源引起的水表面声波的信息:Assin(2πfst+φs)。采用常规的方法(例如频谱分析或相位解调等方法)对探测信号进行分析,即可得到水表面声波的振动幅值As和振动频率fs。可采用图6所示的原理得到水表面声波信息:Assin(2πfst+φs)。It can be seen from the formula (7) that the phase of the detection signal obtained by the photodetector contains the information of the water surface acoustic wave caused by the underwater sound source: As sin(2πf s t+φ s ) . The vibration amplitude A s and vibration frequency f s of the water surface acoustic wave can be obtained by analyzing the detection signal by conventional methods (such as spectrum analysis or phase demodulation). The principle shown in Figure 6 can be used to obtain water surface acoustic wave information: A s sin(2πf s t+φ s ).

具体实施方式三:本实施方式所述的基于激光多点相干探测的水下发声目标位置估计方法为异步多点探测方法,该方法为:采用水下声信号激光相干探测法对待探测水域进行逐点探测,得到各个探测点的探测信号Id(t),即光电探测器接收到的信号,并对每个探测点的探测信号进行频谱分析或相位解调,得到每个探测点的水表面声波振幅值,将水表面声波振幅值最大的探测点的位置作为声源中心位置的最佳估计。Specific Embodiment Three: The method for estimating the position of an underwater sounding target based on laser multi-point coherent detection described in this embodiment is an asynchronous multi-point detection method. Point detection, to obtain the detection signal I d (t) of each detection point, that is, the signal received by the photoelectric detector, and perform spectrum analysis or phase demodulation on the detection signal of each detection point, and obtain the water surface of each detection point Acoustic amplitude value, the position of the detection point with the largest acoustic wave amplitude value on the water surface is used as the best estimate of the center position of the sound source.

具体实施方式四:本实施方式是对实施方式三所述的基于激光多点相干探测的水下发声目标位置估计方法的进一步限定,探测信号Id(t)的获得方法为:Specific embodiment four: this embodiment is the further limitation to the underwater sound target position estimation method based on laser multi-point coherent detection described in embodiment three, the acquisition method of detection signal I d (t) is:

参考光的光振幅Eb随时间t的分布为:The distribution of the light amplitude E b of the reference light with time t is:

Eb(t)=Absin(ω0t+kzb+φ) (1)E b (t)=A b sin(ω 0 t+kz b +φ) (1)

式中,Ab为参考光光振幅Eb的幅值;In the formula, A b is the amplitude of the reference light amplitude E b ;

ω0为激光角频率;ω 0 is the laser angular frequency;

k为激光波数;k is the laser wave number;

zb为参考光的光程;z b is the optical path of the reference light;

φ为激光初始相位。φ is the initial phase of the laser.

测量光的光振幅Em随时间t的分布为:The distribution of the light amplitude E m of the measured light over time t is:

Em(t)=Amsin[(ω0d)t+kzm+φ] (2)E m (t)=A m sin[(ω 0d )t+kz m +φ] (2)

式中,Am为测量光光振幅Eb的幅值;In the formula, A m is the amplitude of the measured light amplitude E b ;

zm为测量光的光程;z m is the optical path of the measuring light;

ωd为声光调制器调制频率。ω d is the modulation frequency of the AOM.

参考光与测量光汇合后的合振幅为:The combined amplitude of the reference light and measuring light is:

E(t)=Eb(t)+Em(t)=Am sin(ω0t+kzm+φ)+Absin[(ω0d)t+kzb+φ] (3)E(t)=E b (t)+E m (t)=A m sin(ω 0 t+kz m +φ)+A b sin[(ω 0d )t+kz b +φ] ( 3)

则相干光光强I(t)为:Then the coherent light intensity I(t) is:

光强由光电探测器接收,光电探测器将Am 2cos20t+kzm+φ)、Ab 2cos2[(ω0d)t+kzb+φ]、和AmAbcos[(2ω0d)t+kzm+kzb+2φ]转变为直流信号,因此,去除直流分量后的光强信号Id表示为:The light intensity is received by the photodetector, which converts A m 2 cos 20 t+kz m +φ), A b 2 cos 2 [(ω 0d )t+kz b +φ], and Am A b cos[(2ω 0d ) t+kz m + kz b +2φ] is transformed into a DC signal, therefore, the light intensity signal I d after removing the DC component is expressed as:

Id(t)=KAmAbcos[k(zm-zb)]=Adcos[ωdt+k(zm-zb)] (5)I d (t)=KA m A b cos[k(z m -z b )]=A d cos[ω d t+k(z m -z b )] (5)

式中,K为光电接收器的放大倍率,记Ad=KAmAbIn the formula, K is the magnification of the photoelectric receiver, write A d =KA m A b ,

光程差zm-zb表示为:The optical path difference z m -z b is expressed as:

zm-zb=2(L1-L2)+2An sin(2πfnt+φn)+2As sin(2πfst+φs) (6)z m -z b =2(L 1 -L 2 )+2A n sin(2πf n t+φ n )+2A s sin(2πf s t+φ s ) (6)

式中,An sin(2πfnt+φn)为自然水表面的振动表达式(An为自然水表面振幅,fn为自然水表面振动频率,φn为自然水表面振动初相位),As sin(2πfst+φs)为水下声源引起的水表面振动表达式(As为水表面声波振幅,fs为水表面声波振动频率,φs为水表面声波振动初相位),L1为水面静止时测量光的初始光程,L2为水面静止时参考光的光程,那么,(L1-L2)为水面静止时两路光的初始光程差,记为ΔL,则光电探测器接收到的光强信号可由下式表述:In the formula, A n sin(2πf n t+φ n ) is the vibration expression of the natural water surface (A n is the amplitude of the natural water surface, f n is the vibration frequency of the natural water surface, and φ n is the initial phase of the natural water surface vibration) , A s sin(2πf s t+φ s ) is the expression of water surface vibration caused by underwater sound source (A s is the amplitude of water surface acoustic wave, f s is the frequency of water surface acoustic wave vibration, φ s is the initial value of water surface acoustic wave vibration phase), L 1 is the initial optical path length of the measuring light when the water surface is stationary, L 2 is the optical path length of the reference light when the water surface is stationary, then (L 1 -L 2 ) is the initial optical path difference of the two paths of light when the water surface is stationary, Denoted as ΔL, the light intensity signal received by the photodetector can be expressed by the following formula:

Id(t)=Ad cos{ωdt+k[2△L+2An sin(2πfnt+φn)+2As sin(2πfst+φs)]} (7)I d (t)=A d cos{ω d t+k[2△L+2A n sin(2πf n t+φ n )+2A s sin(2πf s t+φ s )]} (7)

As sin(2πfst+φs)为由水下声源引起的水表面声波的信息。A s sin(2πf s t+φ s ) is the information of the water surface acoustic wave caused by the underwater sound source.

具体实施方式五:结合图6和图7说明本实施方式,本实施方式是实现实施方式二所述的基于激光多点相干探测的水下发声目标位置估计方法的装置,所述装置包括激光器(1)、光隔离器(2)、1×4光纤耦合器(3)、四个光纤声光调制器(4)、四个光环形器(5)、四个探头(6)、四个光衰减器(7)、四个2×1耦合器(8)、4通道的光电探测器(9)和4通道的信号解调模块(10);Specific embodiment five: This embodiment is described in conjunction with FIG. 6 and FIG. 7. This embodiment is a device for realizing the method for estimating the position of an underwater sounding target based on laser multi-point coherent detection described in Embodiment 2. The device includes a laser ( 1), optical isolator (2), 1×4 fiber coupler (3), four fiber optic acousto-optic modulators (4), four optical circulators (5), four probes (6), four optical Attenuator (7), four 2×1 couplers (8), 4-channel photodetector (9) and 4-channel signal demodulation module (10);

激光器(1)发出的激光经过光隔离器(2)后进入1×4光纤耦合器(3),1×4光纤耦合器(3)的四个输出端分别连接四个光纤声光调制器(4)的输入端;The laser light emitted by the laser (1) enters the 1×4 fiber coupler (3) after passing through the optical isolator (2), and the four output ports of the 1×4 fiber coupler (3) are respectively connected to four fiber optic acousto-optic modulators ( 4) input terminal;

每个光纤声光调制器(4)与一个光环形器(5)、一个光衰减器(7)、一个2×1耦合器(8)和一个探头(6)相对应,每个光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端,该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口,该光环形器(5)的二号端口连接探头(6)的光纤接口,该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端;Each optical fiber acousto-optic modulator (4) corresponds to an optical circulator (5), an optical attenuator (7), a 2×1 coupler (8) and a probe (6), and each optical fiber acousto-optic An output end of the modulator (4) is connected to an input end of a 2×1 coupler (8) through an optical attenuator (7), and the other output end of the optical fiber acousto-optic modulator (4) is connected to an optical circulator (5 ), the No. 2 port of the optical circulator (5) is connected to the optical fiber interface of the probe (6), and the No. 3 port of the optical circulator (5) is connected to another input of the 2×1 coupler (8) end;

四个2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端,光电探测器(9)的输出端连接信号解调模块(10)的信号输入端,信号解调模块(10)的输出端用于连接上位机。The output ends of the four 2×1 couplers (8) are connected to the signal input ends of the photodetectors (9), and the output ends of the photodetectors (9) are connected to the signal input ends of the signal demodulation module (10), and the signal demodulation The output end of the module (10) is used for connecting to the upper computer.

图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个探测点水表面声波的信息。Figure 7 is a schematic diagram of the structure of the device for realizing the synchronous multi-point detection method: the laser with a pigtail (1) outputs a single-frequency laser beam, which is input to the 1×4 fiber coupler (3) after passing through the optical isolator (2) Inside, the 1×4 fiber optic coupler (3) outputs 4 lasers at equal light intensity, and the 4 lasers are respectively input to the optical fiber acousto-optic modulator (4), and the 0-order light output by the optical fiber acousto-optic modulator (4) is used as The reference light, the first-class light output by it is used as the detection light: the 4-way detection light is output to the 4 probes (6) through the optical circulator (5), and the 4 probes (6) emit the detection laser light to the water area to be detected, and receive The laser beam reflected back from the water surface is the reflected laser beam. The phase of the reflected laser beam is modulated by the fluctuation of the water surface. The reflected laser beam is input to the 2×1 coupler (8) after passing through the circulator; After passing through the optical attenuator (7), it is also input to the 2×1 coupler (8), and the coherent signal is generated after the combination of the two beams of light. The 4-way coherent signal is received by the 4-channel photodetector (9), and the photodetector ( 9) convert the coherent signal into an electrical signal, enter the signal demodulation module (10) of 4 channels, and the signal demodulation module (10) uses the principle shown in Figure 6 to process the received signal to obtain the water surface acoustic wave information A s sin(2πf s t+φ s ), send this information to the host computer, and the host computer can obtain the information of the water surface acoustic waves at the 4 detection points.

具体实施方式六:结合图7说明本实施方式,本实施方式是对实施方式五所述装置的进一步限定,本实施方式中,所述探头(6)包括多镜片光耦合器和定焦透镜,多镜片光耦合器的光纤接口作为探头(6)的光纤接口,探测光依次经过多镜片光耦合器和定焦透镜后入射到水表面,经水表面返回的激光依次经过定焦透镜和多镜片光耦合器后进入光环形器(5)的二号端口。Embodiment 6: This embodiment is described in conjunction with FIG. 7. This embodiment is a further limitation to the device described in Embodiment 5. In this embodiment, the probe (6) includes a multi-mirror optical coupler and a fixed-focus lens. The optical fiber interface of the multi-mirror optical coupler is used as the optical fiber interface of the probe (6). The probe light is incident on the water surface after passing through the multi-mirror optical coupler and the fixed-focus lens in sequence, and the laser light returned through the water surface passes through the fixed-focus lens and the multi-lens lens in sequence After the optical coupler, enter the No. 2 port of the optical circulator (5).

该光纤耦合器用于将探测激光出射到待探测水域,同时配合定焦透镜实现反射激光束的接收。The fiber optic coupler is used to emit the detection laser to the water area to be detected, and at the same time cooperate with the fixed-focus lens to realize the reception of the reflected laser beam.

具体实施方式七:结合图8说明本实施方式,本实施方式是实现实施方式四所述的基于激光多点相干探测的水下发声目标位置估计方法的装置,所述装置包括激光器(1)、光隔离器(2)、光纤声光调制器(4)、光环形器(5)、探头(6)、光衰减器(7)、2×1耦合器(8)、光电探测器(9)和信号解调模块(10);Specific Embodiment Seven: This embodiment is described in conjunction with FIG. 8. This embodiment is a device for realizing the method for estimating the position of an underwater sounding target based on laser multi-point coherent detection described in Embodiment 4. The device includes a laser (1), Optical isolator (2), optical fiber acousto-optic modulator (4), optical circulator (5), probe (6), optical attenuator (7), 2×1 coupler (8), photodetector (9) And signal demodulation module (10);

激光器(1)发出的激光经过光隔离器(2)后进入光纤声光调制器(4)的输入端,光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端,该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口,该光环形器(5)的二号端口连接探头(6)的光纤接口,该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端,2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端,光电探测器(9)的输出端连接信号解调模块(10)的信号输入端,信号解调模块(10)的输出端用于连接上位机。The laser light emitted by the laser (1) enters the input end of the optical fiber acousto-optic modulator (4) after passing through the optical isolator (2), and one output end of the optical fiber acousto-optic modulator (4) is connected to 2 × 1 One input end of the coupler (8), the other output end of the optical fiber acousto-optic modulator (4) is connected to the No. 1 port of the optical circulator (5), and the No. 2 port of the optical circulator (5) is connected to the probe (6), the No. 3 port of the optical circulator (5) is connected to another input end of the 2 × 1 coupler (8), and the output end of the 2 × 1 coupler (8) is connected to the photodetector (9 ), the output end of the photodetector (9) is connected to the signal input end of the signal demodulation module (10), and the output end of the signal demodulation module (10) is used to connect to the upper computer.

该装置利用一个光学探头对感兴趣的水域内进行逐点的水表面声波探测,以最大振幅处为声源中心位置的最优估计。The device uses an optical probe to detect water surface acoustic waves point by point in the water area of interest, and takes the point of maximum amplitude as the optimal estimation of the center position of the sound source.

具体实施方式八:结合图8说明本实施方式,本实施方式是对实施方式七所述装置的进一步限定,本实施方式中,所述探头(6)包括多镜片光耦合器和定焦透镜,多镜片光耦合器的光纤接口作为探头(6)的光纤接口,探测光依次经过多镜片光耦合器和定焦透镜后入射到水表面,经水表面返回的激光依次经过定焦透镜和多镜片光耦合器后进入光环形器(5)的二号端口。Embodiment 8: This embodiment is described in conjunction with FIG. 8. This embodiment is a further limitation of the device described in Embodiment 7. In this embodiment, the probe (6) includes a multi-mirror optical coupler and a fixed-focus lens. The optical fiber interface of the multi-mirror optical coupler is used as the optical fiber interface of the probe (6). The probe light is incident on the water surface after passing through the multi-mirror optical coupler and the fixed-focus lens in sequence, and the laser light returned through the water surface passes through the fixed-focus lens and the multi-lens lens in sequence After the optical coupler, enter the No. 2 port of the optical circulator (5).

该光纤耦合器用于将探测激光出射到待探测水域,同时配合定焦透镜实现反射激光束的接收。The fiber optic coupler is used to emit the detection laser to the water area to be detected, and at the same time cooperate with the fixed-focus lens to realize the reception of the reflected laser beam.

Claims (4)

1.基于激光多点相干探测的水下发声目标位置估计方法,其特征在于,该方法为:1. The underwater sound target position estimation method based on laser multi-point coherent detection, it is characterized in that, this method is: 步骤一、在水表面上任意选取四个点作为探测点,且四个探测点构成正方形;Step 1. Randomly select four points on the water surface as detection points, and the four detection points form a square; 步骤二、采用水下声信号激光相干探测法对四个探测点进行探测,分别得到四个探测点的探测信号Id(t),即光电探测器接收到的信号,并对四个探测点的探测信号进行频谱分析或相位解调,得到四个探测点的水表面声波振幅值;Step 2, using the underwater acoustic signal laser coherent detection method to detect the four detection points, respectively obtain the detection signals I d (t) of the four detection points, that is, the signal received by the photoelectric detector, and detect the four detection points Spectrum analysis or phase demodulation is performed on the detection signal to obtain the amplitude values of the water surface acoustic waves at the four detection points; 步骤三、判断四个探测点的水表面声波振幅值是否相同,如果相同,则将探测中心、即四个探测点构成的正方形的中心点O作为声源中心位置的最佳估计,如果不相同,则执行步骤四;Step 3. Determine whether the amplitude values of the water surface acoustic waves at the four detection points are the same. If they are the same, use the detection center, that is, the center point O of the square formed by the four detection points, as the best estimate of the center position of the sound source. If they are not the same , execute step 4; 步骤四、根据振幅最大和次大探测点计算探测系统光学探头移动的方向向量 Step 4. Calculate the moving direction vector of the optical probe of the detection system according to the maximum and second maximum amplitude detection points 其中,Amax1为四个探测点的水表面声波最大振幅值,Amax1对应的位置为探测点Smax1,Amax2为四个探测点的水表面声波次大振幅值,Amax2对应的位置为探测点Smax2,O为四个探测点构成的正方形的中心;Among them, A max1 is the maximum amplitude value of the water surface acoustic wave at the four detection points, the position corresponding to A max1 is the detection point S max1 , A max2 is the second largest amplitude value of the water surface acoustic wave at the four detection points, and the position corresponding to A max2 is Detection point S max2 , O is the center of a square formed by four detection points; 步骤五、按照步骤四得到的向量来移动四个光学探头,然后返回步骤二。Step 5. According to the vector obtained in step 4 to move the four optical probes, then return to step two. 2.根据权利要求1所述的基于激光多点相干探测的水下发声目标位置估计方法,其特征在于,步骤二中的探测信号Id(t)的获得方法为:2. the underwater sound target position estimation method based on laser multipoint coherent detection according to claim 1, is characterized in that, the acquisition method of the detection signal I d (t) in step 2 is: 参考光的光振幅Eb随时间t的分布为:The distribution of the light amplitude E b of the reference light with time t is: Eb(t)=Absin(ω0t+kzb+φ) (1)E b (t)=A b sin(ω 0 t+kz b +φ) (1) 式中,Ab为参考光光振幅Eb的幅值;In the formula, A b is the amplitude of the reference light amplitude E b ; ω0为激光角频率;ω 0 is the laser angular frequency; k为激光波数;k is the laser wave number; zb为参考光的光程;z b is the optical path of the reference light; φ为激光初始相位;φ is the initial phase of the laser; 测量光的光振幅Em随时间t的分布为:The distribution of the light amplitude E m of the measured light over time t is: Em(t)=Amsin[(ω0d)t+kzm+φ] (2)E m (t)=A m sin[(ω 0d )t+kz m +φ] (2) 式中,Am为测量光光振幅Eb的幅值;In the formula, A m is the amplitude of the measured light amplitude E b ; zm为测量光的光程;z m is the optical path of the measuring light; ωd为声光调制器调制频率;ω d is the modulation frequency of the AOM; 参考光与测量光汇合后的合振幅为:The combined amplitude of the reference light and measuring light is: E(t)=Eb(t)+Em(t)=Amsin(ω0t+kzm+φ)+Absin[(ω0d)t+kzb+φ] (3)E(t)=E b (t)+E m (t)=A m sin(ω 0 t+kz m +φ)+A b sin[(ω 0d )t+kz b +φ] ( 3) 则相干光光强I(t)为:Then the coherent light intensity I(t) is: <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>&amp;omega;</mi> <mn>0</mn> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>kz</mi> <mi>m</mi> </msub> <mo>+</mo> <mi>&amp;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>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;omega;</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>&amp;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>&amp;phi;</mi> <mo>&amp;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>&amp;lsqb;</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&amp;omega;</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>&amp;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>&amp;phi;</mi> <mo>&amp;rsqb;</mo> <mo>+</mo> <msub> <mi>A</mi> <mi>m</mi> </msub> <msub> <mi>A</mi> <mi>b</mi> </msub> <mi>cos</mi> <mo>&amp;lsqb;</mo> <msub> <mi>&amp;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>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> <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>&amp;omega;</mi><mn>0</mn></msub><mi>t</mi><mo>+</mo><msub><mi>kz</mi><mi>m</mi></msub><mo>+</mo><mi>&amp;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>&amp;lsqb;</mo><mrow><mo>(</mo><msub><mi>&amp;omega;</mi><mn>0</mn></msub><mo>+</mo><msub><mi>&amp;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>&amp;phi;</mi><mo>&amp;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>&amp;lsqb;</mo><mrow><mo>(</mo><mn>2</mn><msub><mi>&amp;omega;</mi><mn>0</mn></msub><mo>+</mo><msub><mi>&amp;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>&amp;phi;</mi><mo>&amp;rsqb;</mo><mo>+</mo><msub><mi>A</mi><mi>m</mi></msub><msub><mi>A</mi><mi>b</mi></msub><mi>cos</mi><mo>&amp;lsqb;</mo><msub><mi>&amp;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>&amp;rsqb;</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow> 光强由光电探测器接收,光电探测器将Am 2cos20t+kzm+φ)、Ab 2cos2[(ω0d)t+kzb+φ]、和AmAbcos[(2ω0d)t+kzm+kzb+2φ]转变为直流信号,因此,去除直流分量后的光强信号Id表示为:The light intensity is received by the photodetector, which converts A m 2 cos 20 t+kz m +φ), A b 2 cos 2 [(ω 0d )t+kz b +φ], and Am A b cos[(2ω 0d ) t+kz m + kz b +2φ] is transformed into a DC signal, therefore, the light intensity signal I d after removing the DC component is expressed as: Id(t)=KAmAbcos[k(zm-zb)]=Adcos[ωdt+k(zm-zb)] (5)I d (t)=KA m A b cos[k(z m -z b )]=A d cos[ω d t+k(z m -z b )] (5) 式中,K为光电接收器的放大倍率,记Ad=KAmAbIn the formula, K is the magnification of the photoelectric receiver, write A d =KA m A b , 光程差zm-zb表示为:The optical path difference z m -z b is expressed as: zm-zb=2(L1-L2)+2Ansin(2πfnt+φn)+2Assin(2πfst+φs) (6)z m -z b =2(L 1 -L 2 )+2A n sin(2πf n t+φ n )+2A s sin(2πf s t+φ s ) (6) 式中,Ansin(2πfnt+φn)为自然水表面的振动表达式,An为自然水表面振幅,fn为自然水表面振动频率,φn为自然水表面振动初相位,Assin(2πfst+φs)为水下声源引起的水表面振动表达式,As为水表面声波振幅,fs为水表面声波振动频率,φs为水表面声波振动初相位,L1为水面静止时测量光的初始光程,L2为水面静止时参考光的光程,那么,(L1-L2)为水面静止时两路光的初始光程差,记为ΔL,则光电探测器接收到的光强信号可由下式表述:In the formula, A n sin(2πf n t+φ n ) is the vibration expression of the natural water surface, A n is the amplitude of the natural water surface, f n is the vibration frequency of the natural water surface, φ n is the initial phase of the natural water surface vibration, A s sin(2πf s t+φ s ) is the expression of the water surface vibration caused by the underwater sound source, A s is the amplitude of the water surface acoustic wave, f s is the vibration frequency of the water surface acoustic wave, and φ s is the initial phase of the water surface acoustic vibration , L 1 is the initial optical path length of the measurement light when the water surface is still, L 2 is the optical path length of the reference light when the water surface is still, then (L 1 -L 2 ) is the initial optical path difference of the two lights when the water surface is still, denoted as ΔL, the light intensity signal received by the photodetector can be expressed by the following formula: Id(t)=Adcos{ωdt+k[2ΔL+2Ansin(2πfnt+φn)+2Assin(2πfst+φs)]} (7)I d (t)=A d cos{ω d t+k[2ΔL+2A n sin(2πf n t+φ n )+2A s sin(2πf s t+φ s )]} (7) Assin(2πfst+φs)为由水下声源引起的水表面声波的信息。A s sin(2πf s t+φ s ) is the information of the water surface acoustic wave caused by the underwater sound source. 3.实现权利要求2所述的基于激光多点相干探测的水下发声目标位置估计方法的装置,其特征在于,所述装置包括激光器(1)、光隔离器(2)、1×4光纤耦合器(3)、四个光纤声光调制器(4)、四个光环形器(5)、四个探头(6)、四个光衰减器(7)、四个2×1耦合器(8)、4通道的光电探测器(9)和4通道的信号解调模块(10);3. realize the device of the underwater sound target position estimation method based on laser multipoint coherent detection claimed in claim 2, it is characterized in that, described device comprises laser device (1), optical isolator (2), 1 * 4 optical fiber Couplers (3), four fiber optic acousto-optic modulators (4), four optical circulators (5), four probes (6), four optical attenuators (7), four 2×1 couplers ( 8), 4-channel photodetector (9) and 4-channel signal demodulation module (10); 激光器(1)发出的激光经过光隔离器(2)后进入1×4光纤耦合器(3),1×4光纤耦合器(3)的四个输出端分别连接四个光纤声光调制器(4)的输入端;The laser light emitted by the laser (1) enters the 1×4 fiber coupler (3) after passing through the optical isolator (2), and the four output ports of the 1×4 fiber coupler (3) are respectively connected to four fiber optic acousto-optic modulators ( 4) input terminal; 每个光纤声光调制器(4)与一个光环形器(5)、一个光衰减器(7)、一个2×1耦合器(8)和一个探头(6)相对应,每个光纤声光调制器(4)的一个输出端通过光衰减器(7)连接2×1耦合器(8)的一个输入端,该光纤声光调制器(4)的另一个输出端连接光环形器(5)的一号端口,该光环形器(5)的二号端口连接探头(6)的光纤接口,该光环形器(5)的三号端口连接2×1耦合器(8)的另一个输入端;Each optical fiber acousto-optic modulator (4) corresponds to an optical circulator (5), an optical attenuator (7), a 2×1 coupler (8) and a probe (6), and each optical fiber acousto-optic An output end of the modulator (4) is connected to an input end of a 2×1 coupler (8) through an optical attenuator (7), and the other output end of the optical fiber acousto-optic modulator (4) is connected to an optical circulator (5 ), the No. 2 port of the optical circulator (5) is connected to the optical fiber interface of the probe (6), and the No. 3 port of the optical circulator (5) is connected to another input of the 2×1 coupler (8) end; 四个2×1耦合器(8)的输出端连接光电探测器(9)的信号输入端,光电探测器(9)的输出端连接信号解调模块(10)的信号输入端,信号解调模块(10)的输出端用于连接上位机。The output ends of the four 2×1 couplers (8) are connected to the signal input ends of the photodetectors (9), and the output ends of the photodetectors (9) are connected to the signal input ends of the signal demodulation module (10), and the signal demodulation The output end of the module (10) is used for connecting to the upper computer. 4.根据权利要求3所述的实现基于激光多点相干探测的水下发声目标位置估计方法的装置,其特征在于,所述探头(6)包括多镜片光耦合器和定焦透镜,多镜片光耦合器的光纤接口作为探头(6)的光纤接口,探测光依次经过多镜片光耦合器和定焦透镜后入射到水表面,经水表面返回的激光依次经过定焦透镜和多镜片光耦合器后进入光环形器(5)的二号端口。4. the device of realizing the underwater sound target position estimation method based on laser multi-point coherent detection according to claim 3, is characterized in that, described probe (6) comprises multi-mirror optical coupler and fixed-focus lens, and multi-mirror The optical fiber interface of the optical coupler is used as the optical fiber interface of the probe (6). The probe light is incident on the water surface after passing through the multi-mirror optical coupler and the fixed-focus lens in sequence. Enter the No. 2 port of the optical circulator (5) after the circulator.
CN201511021900.5A 2015-12-29 2015-12-29 Underwater sounding target position estimation method based on laser multipoint coherent detection and device for realizing method Active CN105487077B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511021900.5A CN105487077B (en) 2015-12-29 2015-12-29 Underwater sounding target position estimation method based on laser multipoint coherent detection and device for realizing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511021900.5A CN105487077B (en) 2015-12-29 2015-12-29 Underwater sounding target position estimation method based on laser multipoint coherent detection and device for realizing method

Publications (2)

Publication Number Publication Date
CN105487077A CN105487077A (en) 2016-04-13
CN105487077B true CN105487077B (en) 2017-11-03

Family

ID=55674184

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511021900.5A Active CN105487077B (en) 2015-12-29 2015-12-29 Underwater sounding target position estimation method based on laser multipoint coherent detection and device for realizing method

Country Status (1)

Country Link
CN (1) CN105487077B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106546991A (en) * 2016-11-02 2017-03-29 哈尔滨工业大学(威海) The device and method of audible target aerial reconnaissance under water
CN106595836B (en) * 2016-12-26 2019-04-26 哈尔滨工业大学 Frequency extraction method in the case of underwater dual sound source spectral aliasing

Citations (4)

* Cited by examiner, † Cited by third party
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 (en) * 2007-06-27 2007-12-05 中国科学院上海光学精密机械研究所 Method and device for detecting underwater acoustic signal by coherent laser remote sensing
CN104316160A (en) * 2014-11-20 2015-01-28 哈尔滨工业大学 Underwater sound signal instantaneous frequency demodulation method based on wavelet ridges
CN104808208A (en) * 2015-04-16 2015-07-29 浙江大学 Measurement system and method for detecting orientation and size of underwater target on basis of laser sound source

Family Cites Families (1)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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 (en) * 2007-06-27 2007-12-05 中国科学院上海光学精密机械研究所 Method and device for detecting underwater acoustic signal by coherent laser remote sensing
CN104316160A (en) * 2014-11-20 2015-01-28 哈尔滨工业大学 Underwater sound signal instantaneous frequency demodulation method based on wavelet ridges
CN104808208A (en) * 2015-04-16 2015-07-29 浙江大学 Measurement system and method for detecting orientation and size of underwater target on basis of laser sound source

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
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 (en) 2016-04-13

Similar Documents

Publication Publication Date Title
CN106500823B (en) Device for realizing high-sensitivity distributed acoustic wave sensing based on small-diameter multimode optical fiber
CN104808208B (en) Measurement system and method for detecting orientation and size of underwater target on basis of laser sound source
CN102879764B (en) A Method for Estimating the Azimuth of Underwater Sound Source
CN105092014B (en) Distribution type fiber-optic acoustic detection device and detection method based on Wave beam forming
JP2016524715A (en) Optical pulse compression reflector
CN110161458B (en) A Sound Source Localization System Based on Fiber-optic Michelson Interferometer
US12007259B2 (en) Multi-dimensional spatial positioning system and method for disturbance source
CN111443352A (en) Ocean remote sensing bathymetry method and device based on acousto-optic joint sensing
CN106813766B (en) Sound magnetic is the same as the distributed optical fiber sensing system surveyed
CN106546991A (en) The device and method of audible target aerial reconnaissance under water
CN103048662A (en) Three-beam all fiber coherent frequency-modulated continuous-wave laser radar
CN113640832A (en) A multi-beam coherent detection lidar
CN108318123A (en) A kind of near field sound chromatography test method of underwater acoustic array mutual radiation impedance
CN105487077B (en) Underwater sounding target position estimation method based on laser multipoint coherent detection and device for realizing method
CN207487813U (en) A kind of sound and image acquiring device based on wireless laser technology
CN105043524B (en) A kind of method for detecting vibration based on the imaging of laser speckle defocus
CN104819769B (en) A kind of vibration measurement device based on polarization singular point beam laser speckle
Jia et al. Direction of arrival estimation for single microelectromechanical systems vector hydrophone using modified wavelet packet de-noising
CN105372673A (en) Receiving and transmitting integrated inverse SAL (synthetic aperture ladar) system based on acousto-optic frequency shifter
CN105629259B (en) Orthoptic synthetic aperture laser imaging radar local oscillator strengthens optics plural numberization reception device
CN111366232A (en) All-fiber pulse laser Doppler vibrometer and signal processing method thereof
CN107941353B (en) Two-photon correlation-based coherent vortex topology charge measurement method
CN113702946A (en) Coaxial multi-field-of-view fusion linear frequency modulation continuous wave distance and speed measuring method and device
CN113432702B (en) A system and method for detecting acoustic signals across the cavitation layer in the ocean based on optical heterodyne
CN112147623A (en) Multi-region ranging method and system based on chaotic polarization radar

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