CN115687901A - Water surface and underwater target distinguishing method and equipment based on shallow water sound field correlation - Google Patents
Water surface and underwater target distinguishing method and equipment based on shallow water sound field correlation Download PDFInfo
- Publication number
- CN115687901A CN115687901A CN202211273705.1A CN202211273705A CN115687901A CN 115687901 A CN115687901 A CN 115687901A CN 202211273705 A CN202211273705 A CN 202211273705A CN 115687901 A CN115687901 A CN 115687901A
- Authority
- CN
- China
- Prior art keywords
- target
- shallow
- underwater
- depth
- water
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000001228 spectrum Methods 0.000 claims abstract description 14
- 230000006870 function Effects 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 10
- 238000009304 pastoral farming Methods 0.000 claims description 10
- 230000001419 dependent effect Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 3
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000009191 jumping Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- IBBLRJGOOANPTQ-JKVLGAQCSA-N quinapril hydrochloride Chemical compound Cl.C([C@@H](C(=O)OCC)N[C@@H](C)C(=O)N1[C@@H](CC2=CC=CC=C2C1)C(O)=O)CC1=CC=CC=C1 IBBLRJGOOANPTQ-JKVLGAQCSA-N 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
A method and a device for distinguishing underwater and water surface targets based on shallow water sound field correlation belong to the technical field of underwater and water surface target distinguishing. The method aims to solve the problems that the resolution effect is influenced due to the fact that the posture of a vertical array is difficult to maintain in the existing method for resolving the underwater targets on the water surface based on the shallow sea waveguide theory, and the energy attenuation of a high-order normal wave mode along with the distance is not beneficial to the use of a mode characteristic method. The invention respectively uses the front section and the rear section of the horizontal array to respectively form beams for a target, makes the output signals formed by the two sections of beams into a cross spectrum, then carries out theoretical analysis on the change of the cross spectrum phase along with the frequency, and finally utilizes the phase statistical characteristics which are related in depth and are not influenced by the fluctuation of the frequency spectrum amplitude for depth resolution. The invention is used for distinguishing the underwater target on the water surface.
Description
Technical Field
The invention belongs to the technical field of underwater sound surface and underwater target resolution, and relates to a surface and underwater target resolution method and device based on shallow water sound field correlation.
Background
The passive sonar target depth resolution problem is a difficult problem of an underwater acoustic boundary all the time, the main content of the passive sonar target depth resolution problem is to judge whether a target is on the water surface or underwater in noise and interference, and the passive sonar target depth resolution problem has important practical significance for anti-submarine battles. At present, the methods for distinguishing underwater targets on water surface are divided into two categories: one is a matching field based target depth resolution method; another is a depth-resolved method based on normal wave modal characteristics.
The matching field processing method matches the sound field model with the received signals through the relation between the underwater sound channel and the sound source position, further inverts the position of the sound source, and distinguishes underwater targets on the water surface through the depth estimation result. Such methods, which are currently more effective, are matched field techniques (MFPs) using a modeled acoustic field for matching and matched mode techniques (MMPs) based on the shallow sea waveguide theory. Although the MFP method can realize the estimation of the underwater sound source distance and depth, the estimation result is affected by the problems of lack of robustness, environmental mismatch, target motion and the like. In order to improve the robustness of the MFP under adverse conditions, subsequent scholars propose improvements such as spatial filtering and adaptive methods, but in actual use, the MFP still cannot achieve ideal robustness. Until the modal filtering method is applied to the matching field and a good effect is obtained, the MMP technique is developed accordingly. The researchers began to devote attention to the study of a target depth-resolved technique based on the MMP method.
The normal wave in the shallow sea waveguide has the influence of depth on each order of modal function, and a shallow source is difficult to excite the low-order normal wave, and particularly under the typical negative jump layer hydrological condition, the low-order modal energy occupation ratio of the shallow source excitation is low, and the high-order modal energy occupation ratio is high. According to the characteristic, V.E Premus and Conan simplify the sound source positioning problem of the matching field method into two classification problems, provide modal filtering methods with different structures, and distinguish underwater targets on the water surface by taking the energy ratio of modal subspace as statistic. However, in an underwater acoustic environment, array aperture limitations result in overlap between modal subspaces. For this problem, premus introduced Scharf-Friedlander matched subspace Detector, however, this method was not ideal in horizontal matrix effect. Y.C Yang provides a method for positioning a moving sound source based on the formation of synthetic aperture modal beams, does not need environmental information to avoid the problem of environmental adaptation, can use a single array element for estimation, but needs a certain horizontal relative movement distance between a target and a receiving array to avoid the interference between the modalities.
The method for resolving the targets on the water surface and under the water based on the shallow sea waveguide theory still needs to rely on a vertical array or an array with effective vertical apertures (a horizontal array with targets in the array end-fire direction) to obtain the ideal effect. However, the attitude of the vertical array is difficult to maintain, and the attenuation of the energy of the higher-order normal wave mode with distance is not favorable for the use of the mode characteristic method. Therefore, the method for distinguishing the underwater targets on the water surface needs more researches on a method for extracting modal characteristics by using a horizontal array.
Disclosure of Invention
The invention aims to solve the problems that the resolution effect is influenced due to the fact that the posture of a vertical array is difficult to maintain in the existing method for resolving the underwater targets on the water surface based on the shallow sea waveguide theory, and the energy attenuation of a high-order normal wave mode is not beneficial to the use of a mode characteristic method along with the distance.
A method for distinguishing water surface and underwater targets based on shallow water sound field correlation comprises the following steps:
in a shallow water waveguide, the distance r, depth z will be r Where the sound field is represented as a superposition of a series of simple normal waves
wherein ,zs X (f) is the frequency spectrum of the sound source radiation noise, j represents an imaginary number, ρ is the water density,is a depth-dependent modal function, k m and βm Respectively the horizontal wave number and the attenuation coefficient of the mth order normal wave, wherein M is the order of the normal wave;
aiming at an N-element horizontal array near the water bottom, the first N/2 array elements of the horizontal array are used for carrying out beam forming on a target to obtain a complex sound pressure signal p 1 Then N/2 array elements carry out beam forming on the target to obtain a complex sound pressure signal p 2 ;
Based on complex sound pressure p 1 and p2 To obtainRepresenting complex sound pressure p 2 Is further conjugated to obtain
wherein ,r a 、r b a and b represent array element serial numbers which are horizontal distances between the array elements and the target; z is a radical of formula r 、z s In order to be the depth of the film,is a modal function with respect to depth; k is the wave number, k m and kn Respectively of m and n ordersHorizontal wave number of normal wave; beta is a beta m and βn Attenuation coefficients of the mth order and nth order normal waves respectively; d is the array element spacing, and theta is the included angle between the target and the horizontal direction;
when the horizontal distance r between the reference array element and the target is more than ten times larger than the aperture of the array, the reference array element and the target have the same horizontal distance
when m is not equal to n, the distance and the depth are smoothed to obtain
depth of inversion [ eta ] in mth order mode m ,ξ m ]The square average of the internal mode function is
wherein ,SRm The horizontal distance traveled by the intrinsic sound ray in a circulation region;c (z) is the speed of sound at depth z,omega isAngular frequency, f frequency, E constant;
approximating the accumulation of equation (9) as an integral
wherein ,α0 For grazing angles expressed as integral objects, the lower limit of integration α 0I =arccos(max(k(z s ),k(z r ))/k 0 ),k 0 The wave number at the depth where the sound velocity is minimum;
finally obtain p 1 and p2 Correlation coefficient of
wherein ,W(α0 ) Is the angular energy density;
a wideband signal having L frequency domain samples, the horizontal correlation coefficient being expressed as
Γ(f l ) The l term of the broadband signal cross spectrum Γ; w is a group of l (α 0 ) Energy density with respect to frequency-dependent grazing angle;
to determine the phase of gamma
φ(z s ,ω)=angle(Γ) (17)
Finally based on the phase phi (z) s Omega), the binary signal detection model is adopted to realize the identification of shallow source or deep source.
Further, based on the phase phi (z) s Omega) the process of identifying shallow or deep sources using a binary signal detection model comprises the following steps:
based on the phase phi (z) s ω) is made as the frequency domain variance of the phaseMake it intoStatistics for sound source depth discrimination, denoted as R var ;
When R is var Greater than or equal to the frequency domain variance threshold eta var When, the target is a shallow source; otherwise, the target is a deep source.
or ,
based on the phase phi (z) s Omega) the process of identifying shallow or deep sources by using a binary signal detection model comprises the following steps:
based on the phase phi (z) s ω), make the frequency domain first order difference diff (φ (z) of the phase s ω)) as a statistic for sound source depth discrimination, and is denoted as R diff ;
When R is diff Greater than or equal to frequency domain first order difference threshold eta diff When, the target is a shallow source; otherwise, the target is a deep source.
or ,
based on the phase phi (z) s Omega) the process of identifying shallow or deep sources using a binary signal detection model comprises the following steps:
based on the phase phi (z) s ω), the variance var (diff (Φ (zs, ω))) of the first order difference of the phase frequency domain is taken as the statistic of sound source depth judgment and is recorded as R vdiff ;
When R is vdiff Variance threshold eta greater than or equal to first-order difference of phase frequency domain vdiff When, the target is a shallow source; otherwise, the target is a deep source.
Further, the angular energy density W (α) 0 ) The following were used:
wherein, SR (alpha) 0 ) Is an angle alpha 0 Corresponding to the horizontal distance traveled by the intrinsic ray in one of the loop regions,andrespectively a depth z s and zr Grazing angle of (d), realβ(α 0 ) Is SR (alpha) 0 ) Corresponding energy loss.
Further, the complex sound pressure p 1 and p2 The following were used:
wherein k is the wave number.
Further, based on complex sound pressure p 1 and p2 To obtainBased on complex sound pressure p 1 and p2 Determined by the correlation coefficient of (2), complex sound pressure p 1 and p2 The correlation coefficient of (a) is as follows:
Further, the horizontal distance between the array element and the target is as follows:
r i =r+(i-1)dcosθ
wherein ,ri Indicating the horizontal distance of the ith array element from the target.
Further, the wave numberω =2 π f is the angular frequency, f is the frequency, c is the speed of sound in water.
A device for resolving targets on water surface and underwater based on shallow water sound field correlation comprises a processor and a memory, wherein at least one instruction is stored in the memory, and the at least one instruction is loaded and executed by the processor to realize the method for resolving the targets on water surface and underwater based on the shallow water sound field correlation.
Has the advantages that:
the invention respectively uses the front section and the rear section of the horizontal array to respectively form beams for a target, and makes the output signals of the two sections of beam forming into a cross spectrum; through the theoretical analysis of the change of the cross-spectrum phase along with the frequency, the depth resolution is carried out by utilizing a plurality of phase statistical characteristics which are related in depth and hardly influenced by the fluctuation of the frequency spectrum amplitude, the accuracy of the target resolution can be improved, and the depth resolution depends on the phase change of the cross-power spectrum without mode space information or earth-sound parameters. Therefore, the method can well solve the problem that the energy attenuation of the high-order normal wave mode along with the distance is not beneficial to the use of the mode characteristic method, and can well solve the problem that the resolving effect is influenced due to the fact that the posture of the vertical array is difficult to maintain in the water surface underwater target resolving method based on the shallow sea waveguide theory.
Drawings
Fig. 1 is a sound velocity cross-sectional view.
Fig. 2 is a horizontal correlation diagram for different sound source depths.
Fig. 3 (a) to 3 (d) are graphs of phase with frequency for different depths (3 m, 15m, 49m, and 139 m).
Detailed Description
The first specific implementation way is as follows:
the embodiment is a method for distinguishing underwater targets on water surface based on shallow water sound field correlation, which comprises the following steps:
in a shallow water waveguide, the distance r, depth z will be r Where the sound field is represented as a superposition of a series of simple normal waves
wherein ,zs X (f) is the frequency spectrum of the sound source radiation noise, j represents an imaginary number, ρ is the water density,is a depth-dependent modal function, k m and βm The horizontal wave number and the attenuation coefficient of the mth order normal wave are respectively, and M is the order of the normal wave.
Consider an N-element horizontal array close to the water bottom, with the array element spacing d and the target angle theta to the horizontal. Let the horizontal distance between the reference array element and the target be r, the horizontal distance between the ith array element and the target can be expressed as
r i =r+(i-1)dcosθ (2)
Using the first N/2 array elements of the horizontal array to perform wave beam forming on the target to obtain a complex sound pressure signal p 1 Then N/2 array elements form wave beams to the target to obtain a complex sound pressure signal p 2 Is provided with
wherein ,in wave number, ω =2 π f is the angular frequency, f is the frequency, and c is the speed of sound in water. a. b represents array element serial number;is a modal function with respect to depth; k is a radical of formula n and βn Respectively the horizontal wave number and the attenuation coefficient of the nth order normal wave;
wherein ,representing complex sound pressure p 2 The conjugate, | - |, represents the absolute value,representing a smooth average over distance.
when the horizontal distance r between the target and the array is much larger than the aperture of the array (the horizontal distance between the target and the array is more than ten times the aperture of the array), the following approximation can be made
Substituting the formula (6) and the formula (7) into the formula (5), and finishing to obtain
When m = n, the first term in parentheses of equation (8) is a slow varying function of distance and depth; when m ≠ n, the second term in the parenthesis of equation (8) varies strongly with distance and depth due to weak correlation between different modalities, and disappears after smoothing of distance and depth
depth of inversion [ eta ] in mth order mode m ,ξ m ]The square average of the internal mode function is
In the above formula, SR m The horizontal distance traveled by the intrinsic sound ray in a circulation region;c (z) is the speed of sound at depth z;constant E =0.875.
When the frequency of the sound source is higher, multi-order normal waves exist in shallow water, and the grazing angle approaches to be continuous. The accumulation of equation (9) can be approximated as an integral
k m =k(z)cosα z (11)
wherein ,αz Is the glancing angle at depth z. Substituting the formulas (10), (11) and (12) into the formula (9), and finishing to obtain
wherein ,α0 For grazing angles expressed as integration objects, lower integration limit α 0I =arccos(max(k(z s ),k(z r ))/k 0 ),k 0 The wave number at the depth where the speed of sound is the minimum.
Thus, p can be obtained 1 and p2 Correlation coefficient of
Wherein the angular energy density
SR(α 0 ) Is an angle alpha 0 Corresponding to the horizontal distance, alpha, traveled by the intrinsic ray in a region of the cycle zs and αzr Respectively a depth z s and zr Grazing angle of (d), realβ(α 0 ) For a distance SR (alpha) per cycle 0 ) Energy loss of (2).
For wideband signals having L frequency domain samples, the horizontal correlation coefficient can be expressed as
Γ(f l ) Is the l-th term of the wideband signal cross-spectrum Γ. W is a group of l (α 0 ) For the energy density at frequency-dependent grazing angles (the expression over a wide band is denoted as W l (·));
Recording the phase of gamma as
φ(z s ,ω)=angle(Γ) (17)
φ(z s ω) and the correlation γ (z) s ,z r And f) related.
Under the conditions of negative sound velocity gradient (thermocline) and receiving depth close to the water bottom, because a shallow source (surface source) is difficult to excite a low-order mode and a deeper source (submerged source), the high-order mode of the shallow source has higher proportion, the horizontal correlation is low (the correlation coefficient is small), and the phase phi (z) of the cross spectrum is small s ω) sharp jump with frequency; after long-distance transmission, due to boundary reflection and other reasons, the loss of the high-order mode is large, in the deep source excited modes, the low-order mode occupies the dominant position, the horizontal correlation is higher (the correlation coefficient is large), and the phase phi (z) is large s ω) varies relatively smoothly with frequency.
To illustrate this characteristic, simulations were performed in an environment having a sound velocity gradient as shown in fig. 1. The sound source frequency is 200Hz to 400Hz, and the sound source signal is incident from 0 degrees. A uniform horizontal array with array element spacing of 5m, array element number N of 80, horizontal distance r between sound source and reference array element of 10km, and receiving depth z r Is 160m. Fig. 2 is a horizontal correlation coefficient at a receiving position for different sound source depths. The phase variation with frequency at different sound source depths is calculated using equation (17) as shown in fig. 3 (a) to 3 (d).
FIGS. 3 (a) to 3 (d) show phi (z) at sound source depths of 3m, 15m, 49m and 139m s ω) as a function of frequency. As can be seen from the comparison, the shallow source (3 m, 15 m) has a phase which varies strongly with frequency, at [ - π, π]Jumping rapidly; the phase of the deep sources (49 m and 139 m) is flat with frequency, the phase changes almost slowly, and jumps occur only at certain frequencies. The law of the phase change of the shallow source and the deep source is in accordance with the horizontal correlation of the receiving depth. To describe this change rule more intuitively:
frequency domain variance of A-phaseDescribing the discrete degree of the phase in the frequency band, wherein the variance of the phase which changes rapidly along with the frequency is large, and the variance which changes slowly is small;
b is the frequency domain first order difference diff (phi (z) of the phase s ω)), the change speed of the phase within the band is described. Along with the phase of the rapid change of the frequency, more spikes representing the speed jump appear in the first-order difference, the slowly-changed spikes are fewer and are relatively flat, and the underwater target on the water surface can be judged by counting the number of the jumping spikes, diff (phi (z) s ,ω)) l =φ(z s ,ω l+1 )-φ(z s ,ω l ) L =1,2, \8230, L-1 is the L-th term of diff (Φ (zs, ω));
c is the variance var (diff (phi (z) of the first order difference of the phase frequency domain s Omega)) of the phase, and the dispersion degree of the first-order difference of the phase is described, along with the phase which changes rapidly with frequency, the first-order difference has more spikes and large dispersion, and the spikes which change slowly with large variance are fewer, the dispersion is small, and the variance is small.
And judging whether the sound source is on the water surface or under the water by using a binary signal detection model in a detection theory and a noisy signal received by a horizontal array. The problem can be described as follows
wherein ,zlim To resolve the depth, typically set to 15m to 30m below the water, the positive direction of the depth axis z points in the direction of increasing water depth. Depth z of sound source s Less than z lim With the target being a surface source (shallow source), say H 0 (ii) a Depth z of sound source s Greater than z lim With the target being an underwater source (deep source), i.e. assuming H 1 。
The phase change is described using three kinds of feature quantities, and each of them is regarded as a statistic for sound source depth discrimination and is denoted as R var ,R diff and Rvdiff According to the description of the phase law in the previous section, there are
R mod For one of the above statistics, when R mod When the target is greater than or equal to the threshold eta, the target is a shallow source (a water surface target); when the statistic is less than the threshold eta, the target is a deep source (underwater target).
The second embodiment is as follows:
the embodiment is a surface and underwater target distinguishing device based on shallow water sound field correlation, the device comprises a processor and a memory, and it should be understood that the device comprises any device described in the present invention, including the processor and the memory, and the device may also comprise other units and modules which perform display, interaction, processing, control and other functions through signals or instructions;
it should also be appreciated that the memory may include a non-transitory machine-readable medium having stored thereon instructions, which may be used to program a computer system, or other electronic device. Storage media may include, but is not limited to, magnetic storage media, optical storage media; a magneto-optical storage medium comprising: read-only memory ROM, random access memory RAM, erasable programmable memory (e.g., EPROM and EEPROM), and flash memory layers; or other type of media suitable for storing electronic instructions.
The memory stores at least one instruction, and the at least one instruction is loaded and executed by the processor to realize the method for resolving the underwater and water surface targets based on the shallow water sound field correlation.
The above-described calculation examples of the present invention are merely to describe the calculation model and the calculation flow of the present invention in detail, and are not intended to limit the embodiments of the present invention. It will be apparent to those skilled in the art that other variations and modifications can be made on the basis of the foregoing description, and it is not intended to exhaust all of the embodiments, and all obvious variations and modifications which fall within the scope of the invention are intended to be included within the scope of the invention.
Claims (10)
1. A method for distinguishing water surface and underwater targets based on shallow water sound field correlation is characterized by comprising the following steps:
in a shallow water waveguide, the distance r, depth z will be r Where the sound field is represented as a superposition of a series of normal waves
wherein ,zs X (f) is the frequency spectrum of the sound source radiation noise, j represents an imaginary number, rho is the water density,is a depth-dependent modal function, k m and βm Respectively the horizontal wave number and the attenuation coefficient of the mth order normal wave, wherein M is the order of the normal wave;
aiming at an N-element horizontal array near the water bottom, the first N/2 array elements of the horizontal array are used for carrying out beam forming on a target to obtain a complex sound pressure signal p 1 Then N/2 array elements form wave beams to the target to obtain a complex sound pressure signal p 2 ;
Based on complex sound pressure p 1 and p2 To obtain Representing complex sound pressure p 2 Is further conjugated to obtain
wherein ,r a 、r b a and b represent array element serial numbers which are horizontal distances between the array elements and the target; z is a radical of r 、z s In order to be the depth of the film,is a modal function with respect to depth; k is the wave number, k m and kn The horizontal wave numbers of the mth order normal wave and the nth order normal wave respectively; beta is a m and βn Attenuation coefficients of the mth order and nth order normal waves respectively; d is the array element spacing, and theta is the included angle between the target and the horizontal direction;
when the horizontal distance r between the reference array element and the target is more than ten times larger than the aperture of the array, the reference array element and the target have the same structure
wherein ,G mn the real number part is the real number part which contains the array element ordinal number term after being added,
when m is not equal to n, the distance and the depth are smoothed to obtain
depth of inversion [ eta ] in mth order mode m ,ξ m ]The square average of the internal mode function is
wherein ,SRm Water passing through a circulation zone for intrinsic sound raysA flat distance;c (z) is the speed of sound at depth z,omega is angular frequency, f is frequency, E is constant;
approximating the accumulation of equation (9) as an integral
wherein ,α0 For grazing angles expressed as integral objects, the lower limit of integration α 0I =arccos(max(k(z s ),k(z r ))/k 0 ),k 0 Is the wave number at the depth where the sound velocity is the minimum;
finally obtain p 1 and p2 Correlation coefficient of (2)
wherein ,W(α0 ) Is the angular energy density;
a wideband signal having L frequency domain samples, the horizontal correlation coefficient being expressed as
Γ(f l ) The l-th term of the broadband signal cross spectrum Γ; w l (α 0 ) Energy density with respect to frequency-dependent grazing angle;
to determine the phase of gamma
Last radicalIn phase phi (z) s ω), a binary signal detection model is adopted to realize the identification of shallow sources or deep sources.
2. The method for resolving the underwater and water surface targets based on the shallow water sound field correlation as claimed in claim 1, characterized by being based on the phase phi (z) s Omega) the process of identifying shallow or deep sources using a binary signal detection model comprises the following steps:
based on the phase phi (z) s ω) as a frequency domain variance of the phaseThis is taken as the statistic for sound source depth discrimination and is denoted as R var ;
When R is var Greater than or equal to the frequency domain variance threshold eta var When the target is a shallow source; otherwise, the target is a deep source.
3. The method for resolving the underwater and water surface targets based on the shallow water sound field correlation as claimed in claim 1, wherein the method is based on phase phi (z) s Omega) the process of identifying shallow or deep sources using a binary signal detection model comprises the following steps:
based on the phase phi (z) s ω), make the frequency domain first order difference diff (φ (z) of the phase s ω)) as a statistic for sound source depth discrimination, and is denoted as R diff ;
When R is diff Greater than or equal to frequency domain first-order difference threshold eta diff When, the target is a shallow source; otherwise, the target is a deep source.
4. The method for resolving the underwater and water surface targets based on the shallow water sound field correlation as claimed in claim 1, characterized by being based on the phase phi (z) s Omega) the process of identifying shallow or deep sources using a binary signal detection model comprises the following steps:
based on the phase phi (z) s ω), variance var (diff (φ (z)) of the first order difference of the phase frequency domain s ω))) as the sound source depth criterionOther statistics, denoted as R vdiff ;
When R is vdiff Variance threshold eta greater than or equal to first-order difference of phase frequency domain vdiff When the target is a shallow source; otherwise, the target is a deep source.
5. Method for resolving underwater and surface targets based on shallow water sound field correlation according to one of claims 1 to 4, wherein the angular energy density W (α) 0 ) The following were used:
7. The method for resolving the underwater and water surface targets based on the shallow water sound field correlation as claimed in claim 6, wherein the method is based on complex sound pressure p 1 and p2 To obtainBased on complex sound pressure p 1 and p2 Determined by the correlation coefficient of (2), complex sound pressure p 1 and p2 The correlation coefficient of (a) is as follows:
8. The method for resolving the underwater and water surface targets based on the shallow water sound field correlation as claimed in claim 7, wherein the horizontal distance between the array element and the target is as follows:
r i =r+(i-1)d cosθ
wherein ,ri Indicating the horizontal distance of the ith array element from the target.
10. A shallow water acoustic field correlation based water surface and underwater target resolution device, which comprises a processor and a memory, wherein the memory stores at least one instruction, and the at least one instruction is loaded and executed by the processor to implement a shallow water acoustic field correlation based water surface and underwater target resolution method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211273705.1A CN115687901B (en) | 2022-10-18 | 2022-10-18 | Water surface and underwater target distinguishing method and device based on shallow water sound field correlation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211273705.1A CN115687901B (en) | 2022-10-18 | 2022-10-18 | Water surface and underwater target distinguishing method and device based on shallow water sound field correlation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115687901A true CN115687901A (en) | 2023-02-03 |
CN115687901B CN115687901B (en) | 2023-06-16 |
Family
ID=85067117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211273705.1A Active CN115687901B (en) | 2022-10-18 | 2022-10-18 | Water surface and underwater target distinguishing method and device based on shallow water sound field correlation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115687901B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117991241A (en) * | 2024-01-19 | 2024-05-07 | 中国科学院声学研究所 | Target depth information extraction method and system based on modal wave number spectrum transformation |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170031020A1 (en) * | 2015-07-31 | 2017-02-02 | Teledyne Instruments, Inc. | Small aperture acoustic velocity sensor |
CN107765310A (en) * | 2017-09-26 | 2018-03-06 | 浙江海洋大学 | A kind of shallow sound field computational methods for shallow sea layering elastic bottom |
US20180128909A1 (en) * | 2016-11-10 | 2018-05-10 | Northwestern Polytechnical University | Method for obtaining horizontal longitudinal correlation of deep-sea great-depth sound field |
CN109975815A (en) * | 2019-03-22 | 2019-07-05 | 武汉源海博创科技有限公司 | A kind of submarine target multi-beam sonar detection system and method |
CN111708007A (en) * | 2020-06-30 | 2020-09-25 | 中国科学院声学研究所 | Target depth identification method and system based on modal flicker index matching analysis |
CN112098983A (en) * | 2020-09-14 | 2020-12-18 | 中国科学院声学研究所 | Shallow sea horizontal array passive positioning method and system based on spatial domain deconvolution processing |
CN112987004A (en) * | 2021-02-05 | 2021-06-18 | 中国人民解放军国防科技大学 | Water surface and underwater target classification method based on horizontal array in shallow sea environment |
US11209571B1 (en) * | 2018-09-25 | 2021-12-28 | The United States Of America As Represented By The Secretary Of The Navy | System and method for the calibration of a hydrophone line array in a quasi-diffuse ambient sound field |
CN114280615A (en) * | 2021-11-11 | 2022-04-05 | 中国船舶重工集团公司第七一五研究所 | Shallow sea surface underwater target distinguishing method based on correlation matrix characteristic value attenuation coefficient |
CN114280584A (en) * | 2021-11-11 | 2022-04-05 | 中国船舶重工集团公司第七一五研究所 | Convergent region and shadow region distinguishing method based on correlation matrix eigenvalue fading coefficient |
CN114779170A (en) * | 2022-06-05 | 2022-07-22 | 哈尔滨理工大学 | Shallow sea near-field sound source positioning method |
CN115130362A (en) * | 2022-07-08 | 2022-09-30 | 中国科学院声学研究所 | Array element position calibration method and system based on channel matching |
-
2022
- 2022-10-18 CN CN202211273705.1A patent/CN115687901B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170031020A1 (en) * | 2015-07-31 | 2017-02-02 | Teledyne Instruments, Inc. | Small aperture acoustic velocity sensor |
US20180128909A1 (en) * | 2016-11-10 | 2018-05-10 | Northwestern Polytechnical University | Method for obtaining horizontal longitudinal correlation of deep-sea great-depth sound field |
CN107765310A (en) * | 2017-09-26 | 2018-03-06 | 浙江海洋大学 | A kind of shallow sound field computational methods for shallow sea layering elastic bottom |
US11209571B1 (en) * | 2018-09-25 | 2021-12-28 | The United States Of America As Represented By The Secretary Of The Navy | System and method for the calibration of a hydrophone line array in a quasi-diffuse ambient sound field |
CN109975815A (en) * | 2019-03-22 | 2019-07-05 | 武汉源海博创科技有限公司 | A kind of submarine target multi-beam sonar detection system and method |
CN111708007A (en) * | 2020-06-30 | 2020-09-25 | 中国科学院声学研究所 | Target depth identification method and system based on modal flicker index matching analysis |
CN112098983A (en) * | 2020-09-14 | 2020-12-18 | 中国科学院声学研究所 | Shallow sea horizontal array passive positioning method and system based on spatial domain deconvolution processing |
CN112987004A (en) * | 2021-02-05 | 2021-06-18 | 中国人民解放军国防科技大学 | Water surface and underwater target classification method based on horizontal array in shallow sea environment |
CN114280615A (en) * | 2021-11-11 | 2022-04-05 | 中国船舶重工集团公司第七一五研究所 | Shallow sea surface underwater target distinguishing method based on correlation matrix characteristic value attenuation coefficient |
CN114280584A (en) * | 2021-11-11 | 2022-04-05 | 中国船舶重工集团公司第七一五研究所 | Convergent region and shadow region distinguishing method based on correlation matrix eigenvalue fading coefficient |
CN114779170A (en) * | 2022-06-05 | 2022-07-22 | 哈尔滨理工大学 | Shallow sea near-field sound source positioning method |
CN115130362A (en) * | 2022-07-08 | 2022-09-30 | 中国科学院声学研究所 | Array element position calibration method and system based on channel matching |
Non-Patent Citations (4)
Title |
---|
GANG YANG 等: "A Study of the Sound Vertical Correlation in Shallow Water", 《THE 2017 4TH INTERNATIONAL CONFERENCE ON SYSTEMS AND INFORMATICS》 * |
周思同;何琳;帅长庚;杨家轩;: "基于简正波和波叠加法的水下非自由声场重建技术仿真研究", 兵工学报, no. 02 * |
张灵珊;彭朝晖;王光旭;: "空气中声源激发水下声场的水平纵向相关性", 声学学报, no. 05 * |
邱龙皓 等: "浅海宽带声源深度判决方法", 《船舶力学》, vol. 24, no. 2 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117991241A (en) * | 2024-01-19 | 2024-05-07 | 中国科学院声学研究所 | Target depth information extraction method and system based on modal wave number spectrum transformation |
Also Published As
Publication number | Publication date |
---|---|
CN115687901B (en) | 2023-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pailhas et al. | Analysis and classification of broadband echoes using bio-inspired dolphin pulses | |
CN109459744B (en) | Steady self-adaptive beam forming method for realizing multi-interference suppression | |
CN112987004B (en) | Water surface and underwater target classification method based on horizontal array in shallow sea environment | |
CN111025273B (en) | Distortion drag array line spectrum feature enhancement method and system | |
CN111693971B (en) | Wide-beam interference suppression method for weak target detection | |
Kim et al. | Adaptive surface interference suppression for matched-mode source localization | |
CN112255607B (en) | Sea clutter suppression method | |
CN115687901A (en) | Water surface and underwater target distinguishing method and equipment based on shallow water sound field correlation | |
Goldhahn et al. | Waveguide invariant broadband target detection and reverberation estimation | |
US6424596B1 (en) | Method and apparatus for reducing noise from near ocean surface sources | |
CN115169216A (en) | High-robustness confrontation sample generation method and system for underwater sound intelligent camouflage | |
CN106019290B (en) | Weighted broadband time reversal operator decomposition multi-target acoustic imaging method | |
CN111650574B (en) | Underwater space-time self-adaptive processing method and system based on sparse recovery | |
CN117890853A (en) | Single-vector hydrophone shallow sea multi-target direction finding method based on particle swarm optimization compensation | |
US7319640B1 (en) | Noise suppression system | |
CN110850421A (en) | Underwater target detection method based on space-time adaptive processing of reverberation symmetric spectrum | |
CN115201821B (en) | Small target detection method based on strong target imaging cancellation | |
CN113075645B (en) | Distorted matrix line spectrum enhancement method based on principal component analysis-density clustering | |
CN113126029B (en) | Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment | |
Malyshkin | Experimental testing of the efficiency of fast projective adaptive algorithms | |
CN113671507A (en) | Waveguide invariant estimation method based on deep sea vertical array | |
CN113221651A (en) | Seafloor sediment classification method using acoustic propagation data and unsupervised machine learning | |
Jenkins et al. | Bayesian optimization with Gaussian process surrogate model for source localization | |
CN117907998B (en) | Shallow sea broadband sound source ranging method, medium and system | |
Zala et al. | Prewhitening for improved detection by matched‐field processing in ice‐ridging correlated noise |
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 |