CN114563075B - Separation method of deep sea sound field multi-path arrival structure based on single-vector hydrophone - Google Patents
Separation method of deep sea sound field multi-path arrival structure based on single-vector hydrophone Download PDFInfo
- Publication number
- CN114563075B CN114563075B CN202210168409.9A CN202210168409A CN114563075B CN 114563075 B CN114563075 B CN 114563075B CN 202210168409 A CN202210168409 A CN 202210168409A CN 114563075 B CN114563075 B CN 114563075B
- Authority
- CN
- China
- Prior art keywords
- component
- sea
- wave
- vector hydrophone
- horizontal
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
-
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
-
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
-
- 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
Abstract
The invention relates to the technical field of underwater sound detection and sonar, and relates to a separation method of a deep sea sound field multi-path arrival structure based on a single-vector hydrophone. According to the separation method of the multi-path reaching structure of the deep sea sound field based on the single vector hydrophone, which is provided by the technical scheme of the invention, a broadband signal radiated by an underwater sound source with the depth larger than the distribution position is collected through the single vector hydrophone distributed near the deep sea surface, and a scalar sound pressure signal and a particle vibration velocity signal in a time domain are recorded; the method provided by the invention has small calculation amount and low requirement on hardware, and can obtain a pure time domain wave field with a single path.
Description
Technical Field
The invention relates to the technical field of underwater sound detection and sonar, and relates to a separation method of a deep sea sound field multi-path arrival structure based on a single-vector hydrophone.
Background
The acoustic speed characteristics of the ocean surface and sea surface undulations have a significant effect on acoustic propagation. On the one hand, the strong surface wind and sea-air action of the ocean in winter can form a surface acoustic channel on the sea surface, most of the sound ray energy which is larger than the cut-off frequency of the acoustic channel can be bound in the acoustic channel, and the sound ray energy is less contacted with the sea bottom, so that the propagation loss of the propagation mode is reduced. On the other hand, in winter, the wave and swell on the sea surface can significantly affect the reflection and scattering of sound lines on the sea surface under the action of the wind wave, thereby affecting the propagation loss of the sound field, and the influence of the surface sound channel sound propagation by the shallow sea swell [ J ] in physical science and news, 2021,70 (05): 195-204.) is pointed out that the propagation loss at 70km can be increased by 10dB more in consideration of the swell. The wave field recorded in the marine environment is the comprehensive result of superposition of sound rays of different paths in a sound field, the information of the sea surface in the deep sea is carried by sea surface reflection waves, and if the sea surface reflection paths can be separated independently, the influence of sea surface sound velocity change or sea surface fluctuation on sound field phase, transmission loss and other sound transmission rules can be researched independently. In addition, for communication signals, effective signals of the communication signals are carried by direct waves, the multipath paths can cause interference on the signals carried by the direct waves, and the separation of the multipath paths is beneficial to direct wave extraction.
In water acoustics, a beam forming method is adopted based on a vertical array to separate sound rays with different arrival angles. The beamforming is a Conventional pitch angle or azimuth angle estimation method, and is a method for estimating a target pitch angle or azimuth angle by using a cross-spectral density matrix (CSDM), and common methods include a Conventional Beamforming (CBF) method and an adaptive beamforming method.
However, when the beam forming method is used to separate sound rays with different arrival angles, the method has high requirements on the vertical aperture of the receiving array, and the method has large calculation amount. In geophysics, the deep sea bottom reflection information and the multiple underwater reflections are usually separated by calculating the sea bottom reflection coefficient. No case of separation of multi-path arrival structures in water based on a single vector detector has been seen at present.
Disclosure of Invention
The invention aims to make up for the blank of multi-path separation in water based on a single vector detector in the existing scheme and solve the problems of large calculation amount and large difficulty in offshore construction layout of a wave field separation method based on an array, thereby providing an up-going wave field and down-going wave field separation method of a deep sea sound field by using a single vector hydrophone. The invention separates different paths of sound rays in a deep sea sound field based on different polar response characteristics of a vertical vibration velocity detector to an uplink wave field and a downlink wave field and a single vector hydrophone so as to solve the problem of multi-path separation of two most main underwater sound rays in deep sea.
In order to solve the technical problems, the separation method of the multi-path arrival structure of the deep sea sound field based on the single-vector hydrophone, provided by the technical scheme of the invention, is characterized in that a vector hydrophone submerged buoy system is arranged near the sea surface, a sound pressure component and a particle vibration velocity component of an underwater sound source radiation signal are received, and the separation of sound rays of different paths is realized by utilizing the characteristic that a vibration velocity detector has opposite polar response to a wave field with opposite arrival angles and by adopting a mode of adding or subtracting different components; the separation steps of the up-going wave field and the down-going wave field of the sound field are as follows:
step 1): collecting broadband signals radiated by an underwater sound source with the depth larger than the deployment position through a single vector hydrophone deployed near the sea surface in the deep sea, and recording scalar sound pressure signals p (t) of a time domain and components V of particle vibration velocity signals in the vertical z direction z (t) component V of the particle velocity signal in the horizontal x-direction x (t) and the component V of the particle velocity signal in the horizontal y-direction y (t), wherein the x, y and z directions are mutually perpendicular;
step 2): obtaining the sound pressure signal p (t) and the components V of the particle vibration velocity signal in the directions of x, y and z by using fast Fourier transform x (t)、V y (t) and V z (t) corresponding spectra P (f), V x (f)、V y (f) And V z (f) And calculating the acoustic energy flow I in the x, y and z directions x 、I y And I z ；
Step 3): using the acoustic energy flow I in the x and y directions x And I y To obtain the azimuth angle of the underwater sound source
Using the component V of the particle velocity signal in the horizontal x direction x (t) the particle velocity signal is in the horizontal y-directionComponent V y (t) and azimuth of underwater sound source
Radial vibration velocity component V of the mesa time domain r ；
Step 4): simulating by using a ray method BELLHOP to respectively obtain the amplitudes of the direct wave and the sea surface reflected wave on a horizontal component and a vertical component;
step 5): calculating a first amplitude ratio alpha of a vertical component and a horizontal component of a direct wave 1 And a second amplitude ratio alpha of the vertical component to the horizontal component of the sea surface reflected wave 2 ：
Wherein, amp z1 Representing the amplitude, amp, of the direct wave in the vertical component x1 Representing the amplitude of the direct wave in the horizontal component, amp z2 Representing the amplitude, amp, of the sea surface reflection wave in a vertical component x2 Representing the amplitude of the sea surface reflection waves in a horizontal component;
step 6): respectively acquiring a downlink wave field and an uplink wave field of the deep sea sound field to achieve the purpose of separating the uplink wave field and the downlink wave field of the deep sea sound field; wherein the first amplitude ratio alpha is obtained by passing through 1 Corrected radial vibration velocity component V of time domain r Component V of vibration velocity perpendicular to time domain z (t) subtracting to obtain a down-going wave field V of the sound field down ：
V down ＝α 1 *V r -V z (t)
By passing through said second amplitude ratio alpha 2 Corrected time domain radial vibration velocity component V r Component V of vibration velocity perpendicular to z (t) adding to obtain the up-going wavefield V of the sound field up ：
V up ＝α 2 *V r +V z (t)；
Wherein, V r Is the radial vibration velocity component of the time domain.
As an improvement of the method, the sound energy flow I in the x, y and z directions is calculated in the step 2) x 、I y And I z The method comprises the following specific steps:
acoustic power flow in x direction:
acoustic energy flow in the y-direction:
acoustic energy flow in z direction:
wherein the superscripts denote complex conjugate operators, symbols
Representing the real part of the fetched data, P (f) i ) Is a frequency point f i Acoustic pressure signal spectrum of (V) x (f i ) Is a particle vibration velocity signal in the x direction at a frequency point f i Frequency spectrum of (V) y (f i ) Is a particle vibration velocity signal in the x direction at a frequency point f i Spectrum of (a), i =1,2,3 1 And f L Respectively, the minimum and maximum frequencies involved in the operation.
As an improvement of the above method, the step 3) specifically includes:
using acoustic energy flow I in x and y directions x And I y Calculating azimuth of underwater sound source
Using the component V of the time domain particle velocity signal in the horizontal x direction x (t) the particle vibration velocity signal is in waterComponent V in the flat y-direction y (t) and underwater sound source azimuth
Radial vibration velocity component V of synthetic time domain r ：
As an improvement of the method, the deep sea has a sea depth H ranging from 1000 to 6000m.
As an improvement of the method, the distribution depth of the vector hydrophone ranges from 0m to 300m.
As an improvement of the method, the horizontal distance range between the underwater sound source and the vector hydrophone is the arrival range of the direct wave.
As a modification of the above method, the depth of the underwater sound source is 300-600m.
The beneficial effects of the invention are: the up-going wave field or the down-going wave field can be obtained by simple operation of the radial vibration velocity component and the vertical vibration velocity component recorded by the single vector hydrophone subsurface buoy distributed near the sea surface, and the purpose of wave field separation is realized.
Drawings
FIG. 1 is a schematic diagram of propagation paths of direct waves and sea surface reflected waves radiated by a deep sea underwater sound source;
FIG. 2 simulates the acoustic velocity profile;
fig. 3 (a) is a waveform diagram of a time domain sound pressure signal p (t) received by a vector hydrophone;
FIG. 3 (b) shows a component V of a time-domain particle velocity signal received by a vector hydrophone in the vertical z-direction z (t) waveform diagrams;
FIG. 3 (c) shows a component V of a time-domain particle velocity signal received by a vector hydrophone in the horizontal x-direction x (t) waveform diagrams;
FIG. 3 (d) is a vector hydrophoneComponent V of the received time domain particle velocity signal in the horizontal y-direction y (t) waveform diagrams;
FIG. 4 illustrates the arrival angles of the direct wave and the sea surface reflection path obtained by Bellhop simulation;
FIG. 5 (a) shows an up-going wave field obtained by adding the horizontal vibration velocity and the vertical vibration velocity components at a receiving-transmitting distance of 2.2693 km;
FIG. 5 (b) shows the up-going wave field obtained by adding the horizontal vibration velocity and the vertical vibration velocity components when the receiving-transmitting distance is 3.5161 km;
FIG. 6 (a) shows a downlink wave field obtained by subtracting the horizontal vibration velocity component and the vertical vibration velocity component at a receiving-transmitting distance of 2.2693 km;
FIG. 6 (b) shows a downlink wave field obtained by subtracting the horizontal vibration velocity component and the vertical vibration velocity component at a transmit-receive distance of 3.5161 km;
fig. 7 shows correlation coefficients of the horizontal component of the direct wave obtained by combining the wave fields and the simulation.
Detailed Description
The technical scheme provided by the invention is further illustrated by combining the following embodiments.
The invention provides a separation method of a multi-path arrival structure of a deep sea sound field based on a single vector hydrophone. Firstly, a vector hydrophone subsurface system arranged near the deep sea surface receives broadband sound pressure and three-component particle vibration velocity signals radiated by an underwater sound source, the amplitudes of direct waves and sea surface reflected waves on vertical components and horizontal components are obtained by BELLHOP simulation, the horizontal components and the vertical components are corrected, and the corrected horizontal components and the corrected vertical components are added to obtain uplink wave signals; and carrying out subtraction operation on the two signals to obtain a downlink wave signal. The invention realizes the separation of the uplink wave field and the downlink wave field by the wave field addition or subtraction operation by utilizing the characteristic that the vertical vibration velocity detector has opposite polar response to the uplink wave field and the downlink wave field from the deep sea sound field propagation angle. The method for performing the wave field separation by using the single vector detector is beneficial to the conventional method for performing the wave field separation based on the vertical array, and solves the problems of large calculation amount of the vertical array scheme and difficulty in offshore construction layout.
The method provided by the invention is an experimental configuration for transmitting the underwater sound source and receiving the underwater vector detector, and carries out multi-path separation on the received wave field on the premise that the depth and the distance of the sound source are known and the depth of the sound source is greater than the depth of the detector. Because the attenuation of multiple reflection energy in deep sea is large, the wave field received by the deep sea direct wave area mainly comprises direct waves and sea surface reflection, and only the two main paths are considered in the invention. Under the deep sea environment, the single-vector hydrophone submerged buoy is arranged near the sea surface and receives broadband signals radiated by an underwater sound source. The method comprises the following steps of estimating a sound source azimuth angle through sound energy flow, estimating the amplitude ratio of a horizontal component and a vertical component through sound source depth and distance, correcting vibration velocity component energy, merging and the like, and realizing sound ray separation of different paths in deep sea.
Step 1: FIG. 1 is a schematic diagram of the simulation environment in this embodiment, the depth of seawater is 1600m, the level of the seawater depth is unchanged, and the sound velocity is shown in FIG. 2; in simulation, the distribution depth of a vector hydrophone submerged buoy system is 300m, and the distance from the seabed is 1300m; the underwater sound source is a pulse sound source, the depth is 600m, the sound source wavelet is a Rake wavelet, the frequency of the sound source radiation signal is 1-50Hz, and the dominant frequency is 20Hz; the horizontal distance range between the underwater sound source and the vector hydrophone is the arrival range of the direct wave, and the distance is related to the sea depth, the sound velocity profile, the sound source depth and the receiving point depth, in the embodiment, the sea depth is 1600 meters, the sound source is positioned at the depth of 600 meters, the receiving point is positioned at the depth of 300 meters, and the arrival range of the direct wave is 0-6.7km.
In simulation, the horizontal distance of an underwater sound source is increased to 10km from 0km, the distance interval is 0.05km, the azimuth angle of a water surface sound source is set to 60 degrees, and software adopted by time domain wave field simulation is SPECFAM 2D;
only two paths of direct waves (sound rays 1) and sea surface reflected waves (sound rays 2) which play a main role in a sound field are considered, and other paths which pass through the sea surface and sea bottom and are reflected for multiple times are not considered. The purpose of this embodiment is to separate the two sound rays by wave field combination using a single vector hydrophone.
Step 2: 3 (a) -3 (d) show the recorded time-domain sound pressure signal p (t) and three-component particle velocity signal V of the simulated vector hydrophone x (t)、V y (t) and V z (t); obtaining the frequency spectrum of the sound pressure and the horizontal vibration velocity component by using fast Fourier transform, and then calculating the sound energy flow I in the x direction by using the following formula x And the flow of acoustic energy I in the y direction y ：
Wherein the superscript denotes a complex conjugate operator, the symbol
Representing the real part of the data; p (f) i )、V x (f i )、V y (f i ) Are respectively frequency points f i A sound pressure signal frequency spectrum, an x-direction mass point vibration velocity signal frequency spectrum and a y-direction mass point vibration velocity signal frequency spectrum, wherein i =1,2,3 1 And f L To select the upper and lower boundaries of the frequency range, in this embodiment, f 1 And f L The frequency range in between is chosen to be 1-50Hz.
And step 3: calculating the azimuth angle of the underwater sound source by using the following formula according to the sound energy flow in the x direction and the y direction:
calculating to obtain 60 degrees of azimuth angle of the sound source;
the vibration velocity component in the radial direction is synthesized using the following formula:
and 4, step 4: simulating by using BELLHOP according to the known sound source depth Ds, the known wave detection point depth Dr, the known transceiving distance r, the known sea depth H and the known sea depth profile; calculating the variation curve of energy along with the transmitting-receiving distance of the direct wave (sound ray 1) of the underwater sound source on the vertical component and the radial component by using the known sound source depth and the known sound source distance, and obtaining the amplitude of the direct wave on the horizontal component and the vertical component as shown in figure 4; calculating the change curve of the energy of the sea surface reflected wave (sound ray 2) of the underwater sound source along with the receiving and transmitting distances on the vertical component and the radial component by using the known sound source depth and the known sound source distance, and obtaining the amplitude of the sea surface reflected wave on the horizontal component and the vertical component as shown in figure 4;
and 5: calculating a first amplitude ratio alpha of the direct sound wave line between the vertical component and the horizontal component by the following formula 1 ：
α 1 ＝Amp z1 /Amp x1 ；
Calculating a second amplitude ratio alpha of the sea surface reflected wave in the vertical component and the horizontal component by the following formula 2 ：
α 2 ＝Amp z2 /Amp x2 ；
Step 6: correcting radial component signals of a time domain respectively aiming at direct waves and sea surface reflected waves according to a first amplitude ratio and a second amplitude ratio of a vertical component and a radial component which are obtained by calculating the depth and the distance of a sound source, and performing wave field combination on the radial component signals and the vertical component; wherein
And (3) performing addition operation on the radial component and the vertical component to obtain an upgoing wave field:
V up ＝α 2 *V r +V z (t)；
wherein v is r Is the radial vibration velocity component of the time domain, v z For the vertical velocity component in the time domain, take the wavefield merging process with the transceiving distances of 2.2693km and 3.5161km as an example, see fig. 5 (a) and fig. 5 (b), respectively. It can be seen that in fig. 5 (a), after the wavefield is combined, the direct arrival near the arrival time of 1.6s is strengthened, and the sea surface reflection near 1.65s is effectively suppressed, which corresponds to only the up-going wavefield left after the wavefield is combined,while the down going wavefield is removed. In fig. 5 (b), the direct wave and the sea surface reflection wave are aliased together in time, the arrival time is about 2.4s, after the wave fields are combined, the direct wave about 2.4s is enhanced, and the sea surface reflection wave immediately after the direct wave and the sea surface reflection wave are suppressed. It should be noted here that, due to the existence of the sound velocity profile, the arrival angle of the direct wave is inverted in the process of increasing the transceiving distance from small to large, the arrival angle of the direct wave in this embodiment is changed from negative to positive at 3.96km, and correspondingly, after the wave field merging is performed, the sea surface reflection is still effectively suppressed, and the direct wave is not reinforced as within 3.96km but is also weakened, but this embodiment is only used to illustrate the technical solution of the present invention and is not limited, and the method provided by the present invention is still effective for the transceiving distance beyond 3.96km in the wave field separation.
And (3) carrying out subtraction operation on the radial component and the vertical component to obtain a downstream wave field:
V down ＝α 1 *V r -V z (t)；
also take the wave field merging process with the transceiving distance of 2.2693km and 3.5161km as an example, see fig. 6 (a) and fig. 6 (b), respectively. It can be seen that in fig. 6 (a), after the wavefield is combined, the direct arrival at the arrival time of around 1.6s is suppressed, and the sea surface reflection around 1.65s is enhanced, which corresponds to the fact that only the down-going wavefield is left after the wavefield is combined, and the up-going wavefield is removed. In fig. 6 (b), the direct wave and the sea surface reflection wave are aliased together in time, the arrival time is between 2.4s, after the wavefield is combined, the direct wave near 2.4s is effectively suppressed, and the sea surface reflection wave immediately after the direct wave is enhanced. It should be noted here that, due to the existence of the sound velocity profile, the arrival angle of the direct wave is inverted in the process of increasing the transceiving distance from small to large, the arrival angle of the direct wave in this embodiment is changed from negative to positive at 3.96km, and correspondingly, beyond 3.96km, after implementing the wave field merging, the direct wave is effectively suppressed, and the sea surface reflected wave is not strengthened as within 3.96km, but is weakened as well, but this embodiment is only used to illustrate the technical solution of the present invention and is not limited, and the method provided by the present invention is still effective for the transceiving distance beyond 3.96km, in the wave field separating.
In order to examine the effect of combining and separating the uplink wave field and the downlink wave field of the wave field, the correlation coefficient of the direct wave obtained after the wave field separation and the correlation coefficient of the pure direct wave obtained by theoretical calculation are calculated, and the correlation coefficients at different distances are shown in a figure 7. It can be seen that the correlation coefficient between the upgoing waves obtained after the wave field separation and the correlation coefficient obtained by theoretical calculation is greater than 0.8 in the range of 6.5 km. At this sea depth and sound velocity profile, the direct wave is almost absent beyond a distance of 6.5 km. Simulation data processing results show that the method for separating the uplink wave field and the downlink wave field is effective, and two paths of direct waves radiated by an underwater sound source and sea surface reflection can be effectively separated under the condition that the depth of a detector is smaller than the depth of the sound source.
It can be seen from the above description that the present invention can obtain an up-going wavefield or a down-going wavefield by simple operation of the radial vibration velocity component and the vertical vibration velocity component recorded by the subsurface buoy of the single vector hydrophone placed near the sea surface, and achieve the purpose of wavefield separation.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention and are not limited. Although the present invention has been described in detail with reference to the embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (7)
1. A separation method of a deep-sea sound field multi-path arrival structure based on a single-vector hydrophone is used for separating an uplink wave field and a downlink wave field of a deep-sea sound field by processing a received broadband signal radiated by an underwater sound source, and is characterized by comprising the following steps of:
step 1): collecting broadband signals radiated by an underwater sound source with the depth larger than the deployment position through a single vector hydrophone deployed near the sea surface in the deep sea and recording scalar sound in a time domainThe component V of the pressure signal p (t) and the particle vibration velocity signal in the vertical z direction z (t) component V of the particle velocity signal in the horizontal x-direction x (t) and the component V of the particle velocity signal in the horizontal y-direction y (t) wherein the x, y and z directions are perpendicular to each other;
step 2): obtaining the sound pressure signal p (t) and particle vibration velocity components V in the directions of x, y and z by using fast Fourier transform x (t)、V y (t) and V z (t) spectra P (f), V corresponding to x (f)、V y (f) And V z (f) And calculating the acoustic energy flow I in the x, y and z directions x 、I y And I z ；
And step 3): using the acoustic energy flow I in the x and y directions x And I y To obtain the azimuth angle of the underwater sound source
Using the component V of the time domain particle velocity signal in the horizontal x direction x (t) a component V of the particle velocity signal in the horizontal y-direction y (t) and azimuth of underwater sound source
Radial vibration velocity component V of synthetic time domain r ；
Step 4): simulating by using a ray method BELLHOP to respectively obtain the amplitudes of the direct wave and the sea surface reflected wave on a horizontal component and a vertical component;
and step 5): calculating a first amplitude ratio alpha of a vertical component and a horizontal component of a direct wave 1 And a second amplitude ratio alpha of the vertical component to the horizontal component of the sea surface reflected wave 2 ：
Wherein, amp z1 Representing the amplitude of the direct wave in the vertical component, amp x1 Representing the amplitude of the direct wave in the horizontal component，Amp z2 Representing the amplitude of the sea surface reflection wave in the vertical component, amp x2 Representing the amplitude of said sea surface reflection waves in a horizontal component;
step 6): respectively acquiring a downlink wave field and an uplink wave field of the deep sea sound field to achieve the purpose of separating the uplink wave field and the downlink wave field of the deep sea sound field; wherein the first amplitude ratio alpha is obtained by passing through 1 Corrected radial vibration velocity component V of time domain r Vertical component V of particle velocity in time domain z (t) subtracting to obtain the downgoing wave field V of the sound field down ：
V down ＝α 1 *V r -V z (t)；
By passing through said second amplitude ratio alpha 2 Corrected radial vibration velocity component V of time domain r Perpendicular component V to particle velocity z (t) adding to obtain the up-going wavefield V of the sound field up ：
V up ＝α 2 *V r +V z (t)；
Wherein, V r Is the radial vibration velocity component of the time domain.
2. The separation method of the multi-path arrival structure of the deep-sea sound field based on the single-vector hydrophone as claimed in claim 1, wherein the acoustic energy flow I in x, y and z directions is calculated in the step 2) x 、I y And I z The method comprises the following specific steps:
acoustic power flow in x direction:
acoustic power flow in the y direction:
acoustic power flow in z direction:
wherein the superscripts denote complex conjugate operators, symbols
Representing the real part of the fetched data, P (f) i ) Is a frequency point f i Acoustic pressure signal spectrum of (V) x (f i ) Is a particle vibration velocity signal in the x direction at a frequency point f i Frequency spectrum of (V) y (f i ) Is a mass point vibration velocity signal in the y direction at a frequency point f i Spectrum of (a), i =1,2,3 1 And f L Respectively, the minimum and maximum frequencies involved in the operation.
3. The method for separating the multi-path arrival structure of the deep-sea sound field based on the single-vector hydrophone according to claim 1, wherein the step 3) specifically comprises the following steps:
using acoustic energy flow I in x and y directions x And I y Calculating azimuth of underwater sound source
Using the component V of the time domain particle velocity signal in the horizontal x direction x (t) component V of the particle velocity signal in the horizontal y-direction y (t) and underwater sound source azimuth
Radial vibration velocity component V of synthetic time domain r ：
4. The method for separating the multi-path arrival structure of the deep sea sound field based on the single vector hydrophone as claimed in the claim 1, wherein the deep sea has a sea depth H ranging from 1000 m to 6000m.
5. The method for separating the deep-sea sound field multi-path arrival structure based on the single-vector hydrophone as claimed in claim 1, wherein the vector hydrophone is arranged in a depth range of 0-300m.
6. The method for separating the multi-path arrival structure of the deep-sea sound field based on the single-vector hydrophone as claimed in claim 1, wherein a horizontal distance range between the underwater sound source and the vector hydrophone is an arrival range of the direct wave.
7. The separation method of the multi-path arrival structure of the deep-sea sound field based on the single-vector hydrophone as claimed in claim 1, wherein the depth of the underwater sound source is 300-600m.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210168409.9A CN114563075B (en) | 2022-02-23 | 2022-02-23 | Separation method of deep sea sound field multi-path arrival structure based on single-vector hydrophone |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210168409.9A CN114563075B (en) | 2022-02-23 | 2022-02-23 | Separation method of deep sea sound field multi-path arrival structure based on single-vector hydrophone |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114563075A CN114563075A (en) | 2022-05-31 |
| CN114563075B true CN114563075B (en) | 2022-11-22 |
Family
ID=81713450
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210168409.9A Active CN114563075B (en) | 2022-02-23 | 2022-02-23 | Separation method of deep sea sound field multi-path arrival structure based on single-vector hydrophone |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114563075B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO20025328D0 (en) * | 2001-11-07 | 2002-11-06 | Inst Francais Du Petrole | Procedure for absolute preserved amplitude processing of seismic well data |
| CN110488253A (en) * | 2019-09-05 | 2019-11-22 | 西北工业大学 | More way time delay estimation methods based on cross-correlation function delay summation |
| CN113063490A (en) * | 2021-03-12 | 2021-07-02 | 北京科技大学 | Sound field separation method based on sound pressure and particle vibration velocity double-sided measurement |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020118602A1 (en) * | 2001-02-27 | 2002-08-29 | Sen Mrinal K. | Angle dependent surface multiple attenuation for two-component marine bottom sensor data |
-
2022
- 2022-02-23 CN CN202210168409.9A patent/CN114563075B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO20025328D0 (en) * | 2001-11-07 | 2002-11-06 | Inst Francais Du Petrole | Procedure for absolute preserved amplitude processing of seismic well data |
| CN110488253A (en) * | 2019-09-05 | 2019-11-22 | 西北工业大学 | More way time delay estimation methods based on cross-correlation function delay summation |
| CN113063490A (en) * | 2021-03-12 | 2021-07-02 | 北京科技大学 | Sound field separation method based on sound pressure and particle vibration velocity double-sided measurement |
Non-Patent Citations (1)
| Title |
|---|
| 矢量水听器宽带多途信号仿真及俯仰角估计性能研究;吴艳群 等;《声学技术》;20091031;第28卷(第05期);第577-581页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114563075A (en) | 2022-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112083404B (en) | Single-vector hydrophone sound source depth estimation method based on multi-path feature matching | |
| CN103048642B (en) | Method for positioning water sound pulse signal matching field based on frequency domain least squares method | |
| CN111580048B (en) | Broadband sound source depth estimation method using single-vector hydrophone | |
| CN111308474B (en) | Towed deep sea seabed shallow structure acoustic detection system and method | |
| CN113109817B (en) | Vector hydrophone deployment depth estimation method | |
| CN107561486A (en) | It is a kind of based on active when anti-Shallow Water Target Wave arrival direction estimating method | |
| CN111458008B (en) | Ocean earth sound parameter inversion method based on single fixed hydrophone and single mobile sound source | |
| CN112269164A (en) | Weak target positioning method based on interference structure matching processing under deep sea reliable acoustic path | |
| Wu et al. | Side-scan sonar and sub-bottom profiler surveying | |
| Shimura et al. | Time-reversal communication with moving source-receiver in shallow water | |
| CN112987003B (en) | HFM signal separation method and system in active sonar | |
| CN102353942B (en) | Adaptive cancellation method for interference of underwater reverberation | |
| CN114563075B (en) | Separation method of deep sea sound field multi-path arrival structure based on single-vector hydrophone | |
| Josso et al. | On the consideration of motion effects in the computation of impulse response for underwater acoustics inversion | |
| CN115561764A (en) | Moving target depth estimation method based on single-vector hydrophone | |
| CN112924926B (en) | DOA estimation method under multi-path channel | |
| CN115220026A (en) | Matching field ice-crossing positioning method based on ice surface horizontal array | |
| CN113704685A (en) | Deep sea blind deconvolution method based on vertical line array | |
| CN113126029A (en) | Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment | |
| CN112684437A (en) | Passive distance measurement method based on time domain warping transformation | |
| Yayu et al. | Research on Location of Underwater Sound Source Target in Deep Sea Sound Field Based on Bellhop Model | |
| CN110133603A (en) | High-frequency ground wave radar ocean clutter cancellation method based on rooting Euclidean geometry center of gravity | |
| CN112415495B (en) | Submarine reverberation signal simulation method | |
| Xiang-jian et al. | An Imaging Algorithm for High-speed Side-scan Sonar Based on Multi-beam Forming Technology | |
| CN115980661A (en) | Broadband sound source DOA estimation method based on shallow sea channel multipath characteristics |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |