CN116299156A - Hydrophone vertical array element position estimation method and optimization strategy thereof - Google Patents
Hydrophone vertical array element position estimation method and optimization strategy thereof Download PDFInfo
- Publication number
- CN116299156A CN116299156A CN202310601868.6A CN202310601868A CN116299156A CN 116299156 A CN116299156 A CN 116299156A CN 202310601868 A CN202310601868 A CN 202310601868A CN 116299156 A CN116299156 A CN 116299156A
- Authority
- CN
- China
- Prior art keywords
- far
- hydrophone
- vertical array
- array
- field
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005457 optimization Methods 0.000 title claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000005315 distribution function Methods 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 230000005236 sound signal Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000009304 pastoral farming Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- 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
- G01S3/802—Systems for determining direction or deviation from predetermined direction
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The application belongs to the technical field of underwater acoustic signal processing, and provides a hydrophone vertical array element position estimation method and an optimization strategy thereof, wherein the estimation method comprises the following steps: the method comprises the steps of respectively transmitting calibration signals to a hydrophone vertical array through at least two far-field sound sources, respectively receiving the calibration signals transmitted by each far-field sound source through each array element, carrying out beam forming processing on the received signals, and estimating the pitch angle of the actual layout track of the hydrophone vertical array based on the result of the beam forming processing
The corresponding projection pitch angle is formed on a far-field plane where the far-field sound source is located; acquisition based on at least two of said projected pitch angles
Estimate of (2)
Vertical array of hydrophonesAzimuth angle of actual layout track
Estimate of (2)
The method comprises the steps of carrying out a first treatment on the surface of the Based on the following
、
The actual position of each array element is estimated. The method and the device utilize the direction characteristics of the simple wave to form the wave beam through the calibration signals emitted by the two non-coplanar far-field sound sources, and can obtain the complete array type information of the hydrophone vertical array at the same time.
Description
Technical Field
The application belongs to the technical field of underwater acoustic signal processing, and particularly provides a hydrophone vertical array element position estimation method and an optimization strategy thereof.
Background
With the continuous development of sonar equipment and the maturation of beam forming technology, multi-element hydrophone arrays are increasingly used in the field of target detection and positioning, and hydrophone vertical arrays play an important role as a common distribution mode of sensors in water. For a hydrophone vertical array, the accurate acquisition of the positions of each array element has obvious influence on the gain processing and the detection positioning performance of the vertical array, particularly for high-frequency signals, even if small deviation exists in the estimation of the positions of the array and the array elements, the performance of an array signal processing algorithm is greatly reduced, particularly when the array size of the vertical array is large, obvious phenomena such as inclination and suspension can be generated under the action of ocean currents and gravity, and therefore, the actual positions of the array elements must be accurately estimated for signal calibration before the hydrophone vertical array is used for underwater acoustic signal processing.
The related research of the prior estimation method of the vertical array type is less, and the method is mainly divided into two types, namely, the vertical array is used for arranging a pressure sensor, and the inclination degree of the vertical array is estimated by measuring the depths of the head and the tail of the vertical array and the middle position, but the method can only obtain the pitch angle of the vertical array and cannot estimate the inclination azimuth angle of the vertical array; the other is to arrange sound source emission signals right above the vertical array or to estimate the relative arrival time among the array elements by acoustic means such as ambient noise cross correlation processing and the like, and to estimate the inclination degree of the vertical array by combining the sound velocity distribution in water.
Therefore, a method which is simple in algorithm, small in calculation amount and capable of simultaneously estimating the pitch angle and the azimuth angle of the vertical array of the hydrophone is needed. The algorithm provided by the patent is based on the characteristic of the simple wave, the calibration signal is transmitted by utilizing far-field sound sources in two orthogonal planes or two non-orthogonal planes with known included angles, and the tilt pitch angle and the azimuth angle of the vertical array can be estimated by performing conventional wave beam forming processing on the received signals of the vertical array of the hydrophone, so that the algorithm is simple, the calculated amount is small, and the method has obvious advantages compared with other methods.
Disclosure of Invention
The method and the strategy for optimizing the method utilize the direction characteristics of the simple waves to obtain the complete array type information of the hydrophone vertical array simultaneously through the beam forming results of far-field sound source radiation signals in two orthogonal planes or two non-orthogonal planes with known included angles.
The first aspect of the present application provides a hydrophone vertical array element position estimation method, which is used for estimating the true positions of each array element forming a hydrophone vertical array, and includes the following steps:
s1, respectively transmitting calibration signals to the hydrophone vertical array through at least two far-field sound sources, wherein the positions of any two far-field sound sources are not coplanar with a plane formed by ideal layout tracks of the hydrophone vertical array;
s2, for the calibration signals transmitted by each far-field sound source, respectively receiving through each array element and carrying out beam forming processing on the received signals, and estimating the pitch angle of the actual layout track of the hydrophone vertical array based on the result of the beam forming processing
The method comprises the steps of projecting pitch angles corresponding to far-field planes of far-field sound sources, wherein the far-field plane of each far-field sound source is a plane formed by ideal layout tracks of the far-field sound sources and a hydrophone vertical array;
s3, acquiring based on at least two projection pitch angles
Estimate of +.>
And azimuth angle +.of actual layout track of hydrophone vertical array>
Estimate of +.>
；
S4, based on the
、/>
The actual position of each array element is estimated.
Preferably, the calibration signal is a broadband short pulse acoustic signal; and the distance between the far-field sound source and the hydrophone vertical array is more than or equal to 10 times of the water depth of the position where the hydrophone vertical array is positioned.
Preferably, the far-field planes corresponding to the at least two far-field sound sources are perpendicular to each other.
Preferably, the number of array elements of the hydrophone vertical array is greater than or equal to 8, and each array element is arranged at equal intervals.
Further, the beamforming process is performed in step S2 for each calibration signal emitted by a far-field sound source based on the following steps:
s21, receiving the calibration signals through each array element and performing Fourier transform of the following formula:
wherein ,
for the number of array elements of the hydrophone vertical array, +.>
For the serial number of array element, < > for>
Is->
Depth of individual array elements->
For integration time +.>
To the +.>
The time domain received signals generated by the array elements,
is->
Is a frequency spectrum of (2);
s22, will
Transform Jian Zhengbo represents:
wherein 
For water density->
In order to calibrate the frequency spectrum of the signal,/>
depth of far-field sound source +.>
Horizontal spacing for vertical array of far-field sound source and hydrophone, < >>
Serial number of simple wave, +.>
Maximum sequence number for effective reduced wave, < ->
、/>
Respectively +.>
Number Jian Zhengbo, eigenvalue number and eigenvalue function;
s23, further to
The conversion is performed to the following formula:
wherein ,
is->
Glancing angle of the mode ray no Jian Zhengbo, < >>
，/>
Is the average sound velocity of the water body;
s24, for each array element based on the following formula
Beamforming is carried out to obtain a frequency-glancing angle two-dimensional distribution function of the received signal:
wherein ,
for receiving signals about frequency->
And glancing angle->
Is a two-dimensional distribution function of>
Is the interval of array elements;
s25, based on the following formula
Performing inverse Fourier transform to obtain a time-glancing angle two-dimensional distribution function of the received signal:
wherein ,
for receiving signals about the time of reception->
And glancing angle->
Is a two-dimensional distribution function of>
、/>
The upper and lower limits of the frequency integral, respectively;
s26, estimating the following formula
And the corresponding projection pitch angle on the far-field plane where each far-field sound source is located:
wherein ,
serial number of far-field sound source, +.>
For said->
In->
Projection pitch angle corresponding to far-field plane where each far-field sound source is located, < >>
Respectively +.>
Glancing angle of up-going wave and glancing angle of down-going wave of reduced number wave.
Further, the determination is based on the following formula
、/>
：
wherein ,
the sequence numbers of any two different far-field sound sources are respectively given.
The second aspect of the present application further provides an optimization strategy of a hydrophone vertical array element position estimation method, which is used for optimizing the hydrophone vertical array element position estimation method, and includes the following steps:
a1, estimating the actual positions of all array elements by using the hydrophone vertical array element position estimation method;
a2 front based on hydrophone vertical array
The first vertical array subarray is constructed by array elements, and the first vertical array subarray is constructed according to the first vertical array subarray of the hydrophone>
Constructing a second vertical subarray from the array element to the last array element, wherein the number of the array elements of the first vertical subarray and the second vertical subarray is more than or equal to 8;
a3, executing steps S1 to S3 on the first vertical array and the second vertical array respectively, and obtaining the pitch angle of the actual layout track of the first vertical array
Azimuth angle->
Estimate of +.>
、/>
And acquiring pitch angle of actual layout track of second vertical subarray +.>
Azimuth angle->
Estimate of +.>
、/>
；
A4, if
And->
Difference or +.>
And->
If the difference of (2) is greater than the preset threshold, executing step A5, otherwise changing +.>
And returns to step A1 until +.>
Is a value range of (a);
a5 based on the following
、/>
、/>
、/>
The actual positions of the individual array elements are re-estimated.
Preferably, the method comprises the steps of,
the value range of the hydrophone is 1/4 to 3/4 of the number of the array elements of the vertical array of the hydrophone.
According to the hydrophone vertical array element position estimation method and the hydrophone vertical array element position optimization strategy, the hydrophone vertical array in an inclined state is utilized to carry out wave number formation on calibration signals emitted by a plurality of non-coplanar far-field sound sources, and projection information of actual layout tracks of the hydrophone vertical array on each far-field plane is conveniently determined by utilizing symmetry of the forward wave upward traveling wave and the downward traveling wave relative to the incident direction, so that estimation results of pitch angles and azimuth angles of the hydrophone vertical array layout tracks can be simultaneously obtained.
Drawings
FIG. 1 is a schematic diagram of a specific vertical array of hydrophones deployed below the sea surface;
FIG. 2 is a flow chart of a hydrophone vertical array element position estimation method provided in accordance with an embodiment of the present application;
FIG. 3 is a schematic illustration of two orthogonal planes as far field planes in which far field sound sources are disposed, according to a preferred embodiment of the present application;
FIG. 4 is a schematic diagram of a far-field sound source arrangement with an XZ plane as a far-field plane in accordance with a preferred embodiment of the present application;
FIG. 5 is a schematic illustration of far-field acoustic signal incidence when the projection of the hydrophone vertical array actual layout trace on the XZ plane has a counterclockwise offset relative to the OZ axis;
FIG. 6 is a schematic illustration of far-field acoustic signal incidence when the projection of the hydrophone vertical array actual deployment trajectory on the XZ plane has a clockwise offset relative to the OZ axis;
fig. 7 is a schematic diagram of a Jian Zhengbo exploded principle in some specific embodiments;
FIG. 8 is a two-dimensional distribution function of time-glancing angle for a projected pitch angle of a vertical array of hydrophones in the far-field plane of 0℃in some embodiments
Is a picture of (1);
FIG. 9 is a graph showing a two-dimensional distribution function of time-glancing angle for a projected pitch angle of a vertical array of hydrophones in a far-field plane of 5℃in some embodiments
Is a picture of (1);
fig. 10 is a schematic diagram of a three-dimensional trajectory formed by the elements of a vertical array of hydrophones under the combined action of ocean current impact and cable drag in some embodiments.
Detailed Description
The present application will be further described below based on preferred embodiments with reference to the accompanying drawings.
In the description of the embodiments of the present application, it should be noted that, if the terms "upper," "lower," "inner," "outer," and the like indicate an azimuth or a positional relationship based on the azimuth or the positional relationship shown in the drawings, or an azimuth or a positional relationship that a product of the embodiments of the present application conventionally puts in use, it is merely for convenience of describing the present application and simplifying the description, and does not indicate or imply that the device or element to be referred to must have a specific azimuth, be configured and operated in a specific azimuth, and therefore should not be construed as limiting the present application. Furthermore, in the description of the present application, the terms first, second, etc. are used herein for distinguishing between different elements, but not necessarily for describing a sequential or chronological order of manufacture, and may not be construed to indicate or imply a relative importance, and their names may be different in the detailed description of the present application and the claims. In addition, various components on the drawings are enlarged or reduced for ease of understanding, but this is not intended to limit the scope of the present application.
Fig. 1 shows a schematic track of a hydrophone vertical array arranged below the sea surface, in which in order to facilitate the representation of the position and direction in the process of using the hydrophone vertical array to perform underwater acoustic signal processing, an orthogonal coordinate system is generally established by taking the array element at the lowest end as an origin O, and a Z-axis is taken as an ideal layout track of the hydrophone vertical array.
In the actual layout and use process of the hydrophone vertical array, the actual layout track is generally inclined under the influence of ocean currents, so that the actual layout track is not coincident with the ideal layout track. As shown in FIG. 1, the connection line between the origin O and the uppermost array element A is the actual layout track of the hydrophone vertical array, and the attitude deviation from the ideal layout track (namely the OZ axis) can be realized through a pitch angle
Azimuth angle->
Representation, wherein->
Is the included angle between the positive direction of the Z axis and the OA connection line, < >>
Is the angle between the projection of the OA line on the horizontal plane (i.e. XY plane) and the positive direction of the X axis. The deviation of the actual layout track and the ideal layout track of the vertical array of the hydrophone greatly influences the processing of the follow-up vertical array data, and particularly, the high-frequency acoustic signal processing is more sensitive to the array type of the vertical array, so that the accurate data of the actual layout track of the vertical array of the hydrophone is required to be acquired to ensure that the accurate correction is provided for the follow-up acoustic signal processing.
Due to the fact that the multi-path sound propagation characteristic exists in the ocean channel, the ocean channel forms a natural angle filter through multi-path interference, so that the strength of sound signals received by all array elements of the vertical array of the hydrophone is inconsistent at different arrival moments, the strength distribution of the sound signals is related to the depth and the distance of a sound source, and the estimation of the inclination angle of the vertical array by utilizing the relative arrival time of pulses of all array elements is difficult; in addition, there is a method for estimating the vertical array layout track by using the pressure sensor sea test acoustic signal, however, the method can only estimate the pitch angle of the vertical array, but can not effectively estimate the azimuth angle of the vertical array at the same time.
For this purpose, the present application provides a hydrophone vertical array element position estimation method, which can simultaneously obtain estimation results of a pitch angle and an azimuth angle of a hydrophone vertical array layout track, and fig. 2 shows a schematic implementation flow chart of the estimation method in some preferred embodiments, and as shown in fig. 2, the estimation method includes the following steps:
s1, respectively transmitting calibration signals to the hydrophone vertical array through at least two far-field sound sources, wherein the positions of any two far-field sound sources are not coplanar with a plane formed by ideal layout tracks of the hydrophone vertical array;
s2, for the calibration signals transmitted by each far-field sound source, respectively receiving through each array element and carrying out beam forming processing on the received signals, and estimating the pitch angle of the actual layout track of the hydrophone vertical array based on the result of the beam forming processing
The method comprises the steps of projecting pitch angles corresponding to far-field planes of far-field sound sources, wherein the far-field plane of each far-field sound source is a plane formed by ideal layout tracks of the far-field sound sources and a hydrophone vertical array;
s3, acquiring based on at least two projection pitch angles
Estimate of +.>
And azimuth angle +.of actual layout track of hydrophone vertical array>
Estimate of +.>
；
S4, based on the
、/>
The actual position of each array element is estimated.
According to the hydrophone vertical array element position estimation method, firstly, calibration signals are transmitted in at least two far-field planes which are not coplanar with ideal layout tracks of a vertical array through a step S1, then projections of actual pitch angles of the vertical array in all the far-field planes are estimated through a step S2, finally estimation results of at least two far-field planes are utilized to obtain estimation of pitch angles and azimuth angles of the vertical array through steps S3 and S4, and detailed description is given below for specific implementation modes of the steps S1 to S4 by combining drawings and specific embodiments.
FIG. 3 shows a schematic diagram of two orthogonal planes, namely an XZ plane and a YZ plane, respectively, as far-field planes for disposing far-field sound sources in a preferred embodiment, projections of actual layout tracks OA of a hydrophone vertical array on the XZ plane and the YZ plane respectively forming angles with OZ
And->
In the present application, < >>
And->
Pitch angle>
Projected pitch angle in far field plane, when +.>
、/>
Estimate of +.>
、/>
Then, the pitch angle can be obtained by utilizing the geometric relationship>
Azimuth angle
Estimate of +.>
、/>
And then the positions of the array elements are estimated by utilizing the interval information among the array elements.
It should be noted that the embodiment shown in fig. 3, in which two orthogonal planes are taken as far-field planes, is only a preferred embodiment of the present application, in other embodiments of the present application, the number of far-field planes may also be greater than 2, and the angle between any two far-field planes may be an included angle different from 90 degrees, so long as the number of far-field planes is not less than 2, and any two far-field planes are not coplanar with respect to an ideal layout track of a vertical array of hydrophones, without departing from the inventive concept of the present application.
Fig. 4 illustrates schematically an embodiment of the arrangement of the far-field sound source in the far-field plane, taking the XZ plane as an example, and a person skilled in the art can imagine by means of fig. 4 a specific arrangement of the far-field sound source in other far-field planes (such as YZ plane or other far-field planes) formed by an ideal layout trajectory of the far-field sound source and the vertical array of hydrophones.
As shown in fig. 4, the distance between the sound source 1 and the hydrophone vertical array is far greater than the water depth of the position where the hydrophone vertical array is located, for example, in some preferred embodiments, the ratio of the two is greater than or equal to 10, so that the sound signal emitted from the sound source 1 can be regarded as a far-field sound source in the XZ plane, and the sound signal emitted from the sound source can be regarded as a far-field sound signal incident from the horizontal direction; further, to ensure accuracy of estimation, in some preferred embodiments of the present application, the number of array elements of the vertical array of hydrophones is 8 or more (e.g., the vertical array of hydrophones may be an 8-element array, a 16-element array, a 32-element array, etc.), and the individual array elements are disposed at equal intervals.
Fig. 5 and 6 further show the scenario of far-field acoustic signal incidence when the projection of the actual layout track of the hydrophone vertical array on the XZ plane has a counterclockwise and clockwise offset relative to the OZ axis, as shown in fig. 5 and 6, the projection of the hydrophone vertical array on the XZ plane has a projection pitch angle
And the inclination is such that when the vertical array performs underwater sound positioning, it is used for calibrating the horizontal directionThe X' axis forms an angle +.>
And the included angle->
=/>
And in this case the far field acoustic signal incident horizontally on the XZ plane (i.e. incident at 0 grazing angle to the X-axis) will be "considered" by the tilted hydrophone vertical array as being along the grazing angle +.>
Incidence; it follows that a horizontal incoming far field acoustic signal can be used as the calibration signal, and a tilted hydrophone vertical array will have a glancing angle to the calibration signal relative to the X' axis>
Estimation is performed and the ∈ ->
Estimated value of (i.e.)>
Estimate of +.>
。
After the calibration signals are transmitted by the plurality of far-field sound sources through the step S1, the step S2 processes the calibration signals transmitted by the far-field sound sources respectively to obtain estimated values of projection pitch angles corresponding to the hydrophone vertical arrays in the far-field planes, and specifically, in the embodiment of the application, the step S2 comprises the following steps:
s21, receiving the calibration signals through each array element and performing Fourier transform of the following formula:
wherein ,
for the number of array elements of the hydrophone vertical array, +.>
For the serial number of array element, < > for>
Is->
The depth of the individual array elements is determined,
for integration time +.>
To the +.>
The time domain received signals generated by the array elements,
is->
Is a frequency spectrum of (c).
S22, will
Transform Jian Zhengbo represents:
wherein 
For water density->
For calibrating the frequency spectrum of the signal +.>
Depth of far-field sound source +.>
Horizontal spacing for vertical array of far-field sound source and hydrophone, < >>
Serial number of simple wave, +.>
Maximum sequence number for effective reduced wave, < ->
、/>
Respectively +.>
Number Jian Zhengbo, eigenvalue wavenumber and eigenvalue function.
S23, further to
The conversion is performed to the following formula:
wherein ,
is->
Glancing angle of the mode ray no Jian Zhengbo, < >>
，/>
Is the average sound velocity of the water body;
s24, for each array element based on the following formula
Beamforming is carried out to obtain a frequency-glancing angle two-dimensional distribution function of the received signal:
wherein ,
for receiving signals about frequency->
And glancing angle->
Is a two-dimensional distribution function of>
Is the spacing of the array elements.
S25, based on the following formula
Performing inverse Fourier transform to obtain a time-glancing angle two-dimensional distribution function of the received signal:
wherein ,
for receiving signals about the time of reception->
And glancing angle->
Is a two-dimensional distribution function of>
、/>
Respectively isUpper and lower limits of frequency integration;
steps S22 to S25 obtain glancing angle information of forward wave up-going wave and down-going wave by means of forward wave beam forming, fig. 7 illustrates the principle of forward wave decomposition by taking OX plane as an example, as shown in fig. 7, for the case of vertical array receiving far-field acoustic signals, ocean environment can be simplified into a layered dielectric waveguide with unchanged level, under the model, frequency domain signals received by each array element are decomposed into a plurality of forward waves, and modal rays can be obtained by synthesizing the forward wave and the down-going wave.
Based on the principle, the frequency-glancing angle two-dimensional distribution function of the received signal can be obtained by carrying out wave beam formation on the frequency domain signals of each array element of the vertical array
When the number of array elements is far greater than 1, the glancing angles of the forward wave upward traveling wave and the downward traveling wave are approximately in a symmetrical relation relative to the glancing angles of the modal rays, so that glancing angles of the forward wave upward traveling wave and the downward traveling wave can be obtained through a frequency-glancing angle two-dimensional distribution function, and glancing angle information of the modal rays of each number Jian Zhengbo is obtained.
For the broadband short pulse signal, when the signal frequency is far greater than the Jian Zhengbo cut-off frequency, the eigenfunction, group velocity and modal ray pitch angle of the simple wave are almost unchanged with the frequency, so that the method can be further applied to
Performing an inverse Fourier transform process to obtain a time-glancing angle two-dimensional distribution function of the received signal>
FIGS. 8 and 9 show the time-glancing angle two-dimensional distribution function +.A projection pitch angle of the hydrophone vertical array in a far-field plane is 0 degree and 5 degrees, respectively, obtained by the above steps>
The areas with higher sound intensity in the figure represent the arrival of the upward wave and the downward wave of the multiple number of the normal waves respectivelyReaching the moment and glancing angle.
As can be seen from fig. 8 and 9, due to the different group velocities of the effective Jian Zhengbo, the arrival time of each simple wave is different, and each simple wave uplink wave and downlink wave are separable not only in the glancing angle, but also in the relative arrival time, so that the estimation accuracy of the glancing angle of the simple wave uplink wave and the glancing angle of the downlink wave is improved.
Further, the previous image can be extracted from fig. 8 and 9 by manual recognition or various existing image recognition algorithms
Glancing angles of the up-going and down-going waves and estimating said +.>
Projection pitch angles corresponding to far-field planes where all far-field sound sources are located:
wherein ,
serial number of far-field sound source, +.>
For said->
In->
Projection pitch angle corresponding to far-field plane where each far-field sound source is located, < >>
Respectively +.>
Glancing angle of up-going wave and glancing angle of down-going wave of reduced number wave.
Further, if the estimated values of the projected pitch angles of at least two far-field planes are obtained, the pitch angle and the azimuth angle of the actual layout track of the hydrophone vertical array can be estimated through the step S3.
Specifically, in some preferred embodiments, for two mutually perpendicular far field planes, the estimation can be made by:
wherein ,
the serial numbers of any two different far-field sound sources are respectively shown in the figure 3 as an example, and the projection pitch angles of the actual layout tracks OA of the hydrophone vertical array in the figure on the XZ plane and the YZ plane are respectively +.>
、/>
The estimated values are +.>
、
The following formula is introduced:
in addition, for two far-field planes with angles other than 90 °, the relationship between the angles of the two far-field planes needs to be considered when estimating the pitch angle and the azimuth angle, and the above processes are known to those skilled in the art and are not described herein.
Although the estimation of the pitch angle and the azimuth angle of the overall array can be obtained by the hydrophone vertical array element position optimization method, in some areas where complex sea currents exist, as shown in fig. 10, the array types of the vertical array are difficult to keep consistent because the degree of the sea water flow velocity changes with depth, and each array element forms an "inflection point" shown as a point A' in the figure on the three-dimensional track under the combined action of sea current impact and cable dragging at a part of the depth with a larger flow velocity, the inflection point can enable the upper and lower array type tracks to be obviously different, even a single pitch angle and a single azimuth angle cannot be used for representing the same, and therefore, some preferred embodiments of the application also provide an optimization strategy for the hydrophone vertical array element position estimation method, which comprises the following steps:
a1, estimating the actual positions of all array elements by using the hydrophone vertical array element position estimation method;
a2 front based on hydrophone vertical array
The first vertical array subarray is constructed by array elements, and the first vertical array subarray is constructed according to the first vertical array subarray of the hydrophone>
Constructing a second vertical subarray from the array element to the last array element, wherein the number of the array elements of the first vertical subarray and the second vertical subarray is more than or equal to 8;
a3, executing steps S1 to S3 on the first vertical array and the second vertical array respectively, and obtaining the pitch angle of the actual layout track of the first vertical array
Azimuth angle->
Estimate of +.>
、/>
And acquiring pitch angle of actual layout track of second vertical subarray +.>
Azimuth angle->
Estimate of +.>
、/>
；
A4, if
And->
Difference or +.>
And->
If the difference of (2) is greater than the preset threshold, executing step A5, otherwise changing +.>
And returns to step A1 until +.>
Is a value range of (a);
a5 based on the following
、/>
、/>
、/>
The actual positions of the individual array elements are re-estimated.
In the steps, firstly, the whole array type information of the vertical array of the hydrophone is obtained, then the vertical array is divided into two subarrays in a certain traversing range, the array type information of the vertical array is obtained respectively, whether inflection points exist or not is judged, if yes, the inflection points are used for estimating the array element positions of the upper subarray and the lower subarray respectively, and otherwise, the array element positions are estimated according to the array type information of the whole vertical array.
Analysis of the matrix type change of the vertical matrix under ocean current impact has shown that such inflection points tend to occur in the middle region of the overall vertical matrix, which, in some preferred embodiments,
the value range of the hydrophone is 1/4 to 3/4 of the number of the array elements of the vertical array of the hydrophone.
While the foregoing is directed to embodiments of the present application, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
Claims (8)
1. The hydrophone vertical array element position estimation method is used for estimating the real positions of all array elements forming the hydrophone vertical array and is characterized by comprising the following steps:
s1, respectively transmitting calibration signals to the hydrophone vertical array through at least two far-field sound sources, wherein the positions of any two far-field sound sources are not coplanar with a plane formed by ideal layout tracks of the hydrophone vertical array;
s2, for the calibration signals transmitted by each far-field sound source, respectively receiving through each array element and carrying out beam forming processing on the received signals, and estimating the pitch angle of the actual layout track of the hydrophone vertical array based on the result of the beam forming processing
The method comprises the steps of projecting pitch angles corresponding to far-field planes of far-field sound sources, wherein the far-field plane of each far-field sound source is a plane formed by ideal layout tracks of the far-field sound sources and a hydrophone vertical array;
s3, based on at least two projection noddingElevation acquisition
Estimate of +.>
And azimuth angle +.of actual layout track of hydrophone vertical array>
Estimate of +.>
；
S4, based on the
、/>
The actual position of each array element is estimated.
2. The hydrophone vertical array element position estimation method of claim 1, wherein:
the calibration signal is a broadband short pulse sound signal; the method comprises the steps of,
the distance between the far-field sound source and the hydrophone vertical array is more than or equal to 10 times of the water depth of the position where the hydrophone vertical array is located.
3. The hydrophone vertical array element position estimation method of claim 1, wherein:
the far-field planes corresponding to at least two far-field sound sources are perpendicular to each other.
4. The hydrophone vertical array element position estimation method of claim 1, wherein:
the number of array elements of the hydrophone vertical array is more than or equal to 8, and all the array elements are arranged at equal intervals.
5. The hydrophone vertical array element position estimation method according to claim 4, wherein the beamforming processing is performed in step S2 for each calibration signal emitted by a far-field sound source based on the steps of:
s21, receiving the calibration signals through each array element and performing Fourier transform of the following formula:
wherein ,
for the number of array elements of the hydrophone vertical array, +.>
For the serial number of array element, < > for>
Is->
Depth of individual array elements->
For integration time +.>
To the +.>
Time domain received signal generated by each array element, < >>
Is that
Is a frequency spectrum of (2);
s22, will
Transform Jian Zhengbo represents:
wherein 
For water density->
For calibrating the frequency spectrum of the signal +.>
Depth of far-field sound source +.>
Horizontal spacing for vertical array of far-field sound source and hydrophone, < >>
Serial number of simple wave, +.>
Maximum sequence number for effective reduced wave, < ->
、/>
Respectively +.>
Number Jian Zhengbo, eigenvalue number and eigenvalue function;
s23, further to
The conversion is performed to the following formula:
wherein ,
is->
Glancing angle of the mode ray no Jian Zhengbo, < >>
，/>
Is the average sound velocity of the water body;
s24, for each array element based on the following formula
Beamforming is carried out to obtain a frequency-glancing angle two-dimensional distribution function of the received signal:
wherein ,
for receiving signals about frequency->
And glancing angle->
Is a two-dimensional distribution function of>
Is the interval of array elements;
s25, based on the following formula
Performing inverse Fourier transform to obtain a time-glancing angle two-dimensional distribution function of the received signal:
wherein ,
for receiving signals about the time of reception->
And glancing angle->
Is a two-dimensional distribution function of>
、/>
The upper and lower limits of the frequency integral, respectively;
s26, estimating the following formula
And the corresponding projection pitch angle on the far-field plane where each far-field sound source is located:
wherein ,
serial number of far-field sound source, +.>
For said->
In->
Projection pitch angle corresponding to far-field plane where each far-field sound source is located, < >>
Respectively +.>
Glancing angle of up-going wave and glancing angle of down-going wave of reduced number wave.
6. The hydrophone vertical array element position estimation method of claim 5, wherein the determination is based on the following formula
、/>
：
wherein ,
the sequence numbers of any two different far-field sound sources are respectively given.
7. An optimization strategy for a hydrophone vertical array element position estimation method, which is used for optimizing the hydrophone vertical array element position estimation method according to claim 1, and is characterized by comprising the following steps:
a1, estimating the actual positions of the array elements by using the hydrophone vertical array element position estimation method according to claim 1;
a2 front based on hydrophone vertical array
The first vertical array subarray is constructed by array elements, and the first vertical array subarray is constructed according to the first vertical array subarray of the hydrophone>
Constructing a second vertical subarray from the array element to the last array element, wherein the number of the array elements of the first vertical subarray and the second vertical subarray is more than or equal to 8;
a3, executing steps S1 to S3 on the first vertical array and the second vertical array respectively, and obtaining the pitch angle of the actual layout track of the first vertical array
Azimuth angle->
Estimate of +.>
、/>
And acquiring pitch angle of actual layout track of second vertical subarray +.>
Azimuth angle->
Estimate of +.>
、/>
；
A4, if
And->
Difference or +.>
And->
If the difference of (2) is greater than the preset threshold, executing step A5, otherwise changing +.>
And returns to step A1 until +.>
Is a value range of (a);
a5 based on the following
、/>
、/>
、/>
The actual positions of the individual array elements are re-estimated.
8. The optimization strategy of the hydrophone vertical array element position estimation method according to claim 7, wherein:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310601868.6A CN116299156B (en) | 2023-05-26 | 2023-05-26 | Hydrophone vertical array element position estimation method and optimization strategy thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310601868.6A CN116299156B (en) | 2023-05-26 | 2023-05-26 | Hydrophone vertical array element position estimation method and optimization strategy thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116299156A true CN116299156A (en) | 2023-06-23 |
| CN116299156B CN116299156B (en) | 2023-08-15 |
Family
ID=86824342
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310601868.6A Active CN116299156B (en) | 2023-05-26 | 2023-05-26 | Hydrophone vertical array element position estimation method and optimization strategy thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116299156B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116593965A (en) * | 2023-07-19 | 2023-08-15 | 中国海洋大学 | Horizontal array type estimation method based on far-field opportunity sound source |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160091361A1 (en) * | 2014-09-25 | 2016-03-31 | Steven E. Crocker | System and method for the calibration of a hydrophone line array |
| CN109270518A (en) * | 2018-10-29 | 2019-01-25 | 中国电子科技集团公司第三研究所 | A kind of the multi-target simulation method, apparatus and sonar of underwater sound vertical vector array sonar |
| CN109283492A (en) * | 2018-10-29 | 2019-01-29 | 中国电子科技集团公司第三研究所 | Multi-target DOA estimation method and underwater sound vertical vector array system |
| CN111398902A (en) * | 2020-03-11 | 2020-07-10 | 哈尔滨工程大学 | Passive distance and direction measuring method for right-angled triangle ternary hydrophone array |
| CN112098938A (en) * | 2020-08-31 | 2020-12-18 | 黑龙江工程学院 | Six-element cone vector array-based underwater acoustic target dimension reduction matching sound field positioning method |
| CN113495275A (en) * | 2021-05-24 | 2021-10-12 | 中国海洋大学 | Passive positioning method, system and application for vertical synthetic aperture of single hydrophone |
-
2023
- 2023-05-26 CN CN202310601868.6A patent/CN116299156B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160091361A1 (en) * | 2014-09-25 | 2016-03-31 | Steven E. Crocker | System and method for the calibration of a hydrophone line array |
| CN109270518A (en) * | 2018-10-29 | 2019-01-25 | 中国电子科技集团公司第三研究所 | A kind of the multi-target simulation method, apparatus and sonar of underwater sound vertical vector array sonar |
| CN109283492A (en) * | 2018-10-29 | 2019-01-29 | 中国电子科技集团公司第三研究所 | Multi-target DOA estimation method and underwater sound vertical vector array system |
| CN111398902A (en) * | 2020-03-11 | 2020-07-10 | 哈尔滨工程大学 | Passive distance and direction measuring method for right-angled triangle ternary hydrophone array |
| CN112098938A (en) * | 2020-08-31 | 2020-12-18 | 黑龙江工程学院 | Six-element cone vector array-based underwater acoustic target dimension reduction matching sound field positioning method |
| CN113495275A (en) * | 2021-05-24 | 2021-10-12 | 中国海洋大学 | Passive positioning method, system and application for vertical synthetic aperture of single hydrophone |
Non-Patent Citations (2)
| Title |
|---|
| YAN LIANG等: "A DOA Estimation algorithm for the vertical line array of vector hydrophone based on data fusion method", 《2020 3RD IEEE INTERNATIONAL CONFERENCE ON INFORMATION COMMUNICATION AND SIGNAL PROCESSING》 * |
| 史文佳: "基于垂直阵的差频匹配场声源距离估计方法", 《声学技术》, vol. 38, no. 2 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116593965A (en) * | 2023-07-19 | 2023-08-15 | 中国海洋大学 | Horizontal array type estimation method based on far-field opportunity sound source |
| CN116593965B (en) * | 2023-07-19 | 2023-09-29 | 中国海洋大学 | Horizontal array type estimation method based on far-field opportunity sound source |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116299156B (en) | 2023-08-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7446262B2 (en) | Systems and methods for synthetic aperture sonar | |
| CN109900256B (en) | Self-adaptive ocean mobile acoustic tomography system and method | |
| CN113702909B (en) | Sound source localization analysis solution calculation method and device based on arrival time difference of sound signals | |
| CN116299156B (en) | Hydrophone vertical array element position estimation method and optimization strategy thereof | |
| EP2205995A2 (en) | Controlling seismic source elements based on determining a three-dimensional geometry of the seismic source elements | |
| CN108845325A (en) | Towed linear-array sonar submatrix error misfits estimation method | |
| Zhao et al. | Open-lake experimental investigation of azimuth angle estimation using a single acoustic vector sensor | |
| CN109725285B (en) | DOA estimation method based on MVDR covariance matrix element self-adaptive phase angle conversion | |
| CN113281706A (en) | Target positioning method and device and computer readable storage medium | |
| CN114325584B (en) | Synthetic aperture-based multi-array-element ultrasonic sound source three-dimensional imaging method and system | |
| US6018494A (en) | Signature estimation of a seismic source | |
| CN111551942B (en) | Underwater autonomous vehicle docking method based on deconvolution algorithm | |
| CN116068493A (en) | Passive sound source positioning method for deep sea large-depth vertical distributed hydrophone | |
| US9578433B2 (en) | Method for self-calibration of a set of sensors, in particular microphones, and corresponding system | |
| CN112098938B (en) | Six-element cone vector array-based underwater acoustic target dimension reduction matching sound field positioning method | |
| CN109669172B (en) | Weak target direction estimation method based on strong interference suppression in main lobe | |
| Santos et al. | Dual accelerometer vector sensor mounted on an autonomous underwater vehicle (AUV)-Experimental results | |
| CN107202975A (en) | A kind of a burst of first attitude error rectification method of two-dimensional vector | |
| Kerstens et al. | An optimized planar MIMO array approach to in-air synthetic aperture sonar | |
| CN115390052A (en) | Imaging sonar object simulation device and method | |
| Sinha et al. | Study of acoustic vector sensor based direction of arrival estimation of in-air maneuvering tonal source | |
| CN116500625B (en) | Recovery imaging method, device, system, electronic equipment and readable storage medium | |
| Li et al. | Underwater data fusion of multiple arrays in a mono-platform for low SNR | |
| Liang et al. | An Improved Broadband Target Orientation Method Using a Vector Hydrophone Vertical Array | |
| CN111505578B (en) | ULA (ultra-low-resolution) target multi-source positioning method and device based on time reversal focusing |
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 |