CN106154276B - Deep seafloor parameter inversion method based on bottom reverberation and propagation loss - Google Patents
Deep seafloor parameter inversion method based on bottom reverberation and propagation loss Download PDFInfo
- Publication number
- CN106154276B CN106154276B CN201610478946.8A CN201610478946A CN106154276B CN 106154276 B CN106154276 B CN 106154276B CN 201610478946 A CN201610478946 A CN 201610478946A CN 106154276 B CN106154276 B CN 106154276B
- Authority
- CN
- China
- Prior art keywords
- reverberation
- propagation loss
- inverting
- seafloor
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H7/00—Measuring reverberation time ; room acoustic measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
- G01N9/24—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity by observing the transmission of wave or particle radiation through the material
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Chemical & Material Sciences (AREA)
- Acoustics & Sound (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
In order to effectively obtain the bottom parameters under abyssal environment; the present invention proposes a kind of deep seafloor parameter inversion method based on bottom reverberation and propagation loss; single hydrophone cloth is put into water first; ship navigates in investigation marine site into walking; emission sound source signal simultaneously; deep-sea experimental data is acquired, bottom reverberation and propagation loss experimental data are extracted.Modeling and simulating is carried out to bottom reverberation model, experimental data is matched with the bottom reverberation time series of model emulation, inverting obtains Bottom sound speed and density.Propagation loss is calculated using ray model to be matched with experiment value, inverting obtains seabed attenuation coefficient.The purpose of the present invention is provide effective technical method to obtain deep seafloor parameter.
Description
Technical field
The present invention relates to a kind of methods carrying out inverting to deep seafloor parameter using bottom reverberation and propagation loss, are applicable in
In abyssal environment, belong to the fields such as Underwater Acoustics Engineering, ocean engineering and sonar technique.
Background technology
In abyssal environment, the acoustics parameters (such as Bottom sound speed, density and attenuation coefficient) in seabed play acoustic propagation
Great influence.Deep sea in-situ measures and the difficulty and cost of sea floor sampling all increase, and can only often obtain local seabed
Parameter, and sea floor sampling is hardly possible in abyssal environment on a large scale.Therefore the method for obtaining deep seafloor parameter is limited,
Seabed geoacoustic inversion is one of feasible means.Inverting there are many method for ground sound parameter, such as Matched Field inverting,
Propagation loss and Waveform Matching inverting, normal mode mode Dispersion inverting, Bottom-Reflection-Loss inverting etc., these methods are all
It is applied in neritic environment, but in abyssal environment and non-verified.Each method has the advantage and disadvantage of oneself.As
It is usually used in remote experimental data with field inverting, reflects the average effect of water body and seabed spatial variations environment, but needs big
Array aperture, and it is computationally intensive, there are nonuniqueness.Propagation loss has good sensitivity with Waveform Matching inverting to seabed decaying
Property, but accurate acoustic environment information is needed to obtain actual measurement propagation loss, such as receive the geometric position of battle array, Sound speed profile.
Normal mode mode Dispersion inverting can obtain high-resolution bottom parameters, but characteristic unobvious in abyssal environment, and
Normal wave pattern near field when the depth of water is larger calculates inaccurate.Inverting based on Bottom-Reflection-Loss, can be with high-resolution anti-
The lift height in seabed, density and the velocity of sound are performed, but insensitive to seabed attenuation coefficient.Due to the difference of acoustic propagation characteristic, perhaps
The ground sound inversion method in more shallow seas deep-sea no longer be applicable in, and without a kind of inversion method can Simultaneous Inversion go out Bottom sound speed,
Density and attenuation coefficient.
Reverb signal is that heterogeneity, sea/seabed interface of medium in ocean and sedimentary during acoustic propagation rise and fall
The set of caused all scattered signals.Deep seafloor reverberation has close relationship with bottom parameters.Deep-sea multi-path signals it
Between delay inequality be more than shallow sea, therefore the bottom reverberation of opposite " pure " can be obtained from underwater sound signal for the anti-of bottom parameters
It drills.When transmitting-receiving, which is closed, sets, sound source is co-located with receiving hydrophone, therefore single base reverberation is mainly main by back scattering
Effect.When bistatic, sound source and receiving hydrophone can be located at the different depth of same vertical line, or at a distance of it is certain away from
From.The present invention uses sound source and receiving hydrophone distance away, and the reverb signal received is mainly by lateral scattering
Caused by forward scattering.It has been investigated that there is centainly the velocity of sound and density in the attenuation law of bottom reverberation intensity at any time
Sensibility, with the propagation loss of distance change, to acoustic attenuation coefficient, there are certain sensibility.
Invention content
Technical problems to be solved
In order to avoid the shortcomings of the prior art, the present invention proposes a kind of deep-sea based on bottom reverberation and propagation loss
Bottom parameters inversion method provides effective technical method to obtain deep seafloor parameter.
Technical solution
A kind of deep seafloor parameter inversion method based on bottom reverberation and propagation loss, it is characterised in that steps are as follows：
Step 1：Single hydrophone cloth is put into water, ship walks boat, with Fixed Time Interval to throwing explosive sound source in the sea, time
Return after investigating marine site is gone through, single hydrophone is fished for and obtains bottom reverberation and propagation loss experimental data；
Step 2：Establish seafloor model：Three-dimensional system of coordinate is established in ocean space, wherein sound source S is located at (0,0, zs), it is single
Hydrophone R is located at (rr,0,zr), the position of arbitrary scatterer is located at (xb,yb, 0), θijIt is the incident graze of j-th of scatterer
Angle, θsjIt is scattering glancing angle,It is the deflection of scattered wave；7 bottom parameters：Seafloor density ρ, Bottom sound speed cp, seabed is declined
Subtract factor alphaλ, seabed surface roughness spectral intensity w2With spectrum index γ2, seabed inhomogeneities spectral intensity w3With spectrum index γ3；
Step 3：Seafloor model based on step 2 carries out theoretical modeling to deep seafloor reverberation：
RLmodel(f,zs,zr,tj)=Intensity (f, zs,zr,tj)+NL(f,zr) (1)
tj=tsj+trj (3)
Wherein, RLmodel(f,zs,zr,tj) it is theoretical reverberation level, Intensity (f, zs,zr,tj) it is reverrberation intensity；F is
The centre frequency of signal, NL are to receive depth as zr, frequency be f when seanoise intensity, I0It is sound source level, cwIt is in seawater
The velocity of sound, αwIt is the absorption of seawater, tsjAnd trjIt is sound source respectively to scattering point and scattering point to the time of receiving hydrophone, σjIt is
Scattering section, δ AjIt is the insonify area of scatterer, tjMeet formula (4) from the relationship in different scatterer propagation times；
Step 4：Bandwidth is carried out to deep seafloor reverberation experimental data to be filtered for 100Hz：
RLexp(f,zs,zr,tj)=20log10(Veff(tj))-mv-10log10(Δf) (4)
Wherein, RLexp(f,zs,zr,tj) it is experiment reverberation level, VeffIt is effective voltage value, mvIt is hydrophone sensitivity, Δ f
It is the frequency bandwidth of filter；
Step 5：Using multifrequency joint inversion, the cost function for Inversion for bottom density is established：
Wherein, m is the number of the frequency for inverting, and n is the number of scatterer in reverberation modeling；Using Parallel Particle Swarm Optimization
Genetic Optimization Algorithm calculates accumulation least square and error under m frequency, obtains the minimum value of object function, correspond at this time
Optimizing parameter optimal value be inverting optimal value；Bottom sound speed is obtained by Hamilton empirical equations：
cp+ 487.7 ρ of=2330.4-1257.0 ρ2 (6)
Step 6：As Given information, foundation declines the seafloor density and Bottom sound speed that step 5 inverting is obtained for inverting
Subtract the cost function of coefficient：
Wherein, l is the number of explosive sound source, TLexp(f,zs,zr,rrk) it is propagation loss measured value, TLmodel(f,zs,zr,
rrk,αλ) be propagation loss model calculation value, when cost function is minimized be seabed attenuation coefficient inverting value.
Single hydrophone can be replaced the array formed with multiple hydrophones.
Hydrophone quantitative range is 1-100, forms line array or ring array.
Advantageous effect
A kind of deep seafloor parameter inversion method based on bottom reverberation and propagation loss proposed by the present invention, first in step
Deep-sea Inversion System is laid in rapid one, acquires experimental data.Then ocean space is established in step 2 and step 3 three-dimensional
Coordinate system models deep seafloor reverberation, obtains the bottom reverberation strength time sequence of model emulation；In step 4 to reality
It tests data processing and obtains bottom reverberation strength time sequence；In step 5, propose based on the more of bottom reverberation time attenuation law
Frequency inverting cost function, inverting obtain seafloor density, and Bottom sound speed is obtained by Hamilton empirical equations.In step 6, utilize
Propagation loss inverting obtains attenuation coefficient.This method has obtained sea in step 1 to step 6 using deep-sea experimental data inverting
Bottom parameter, to provide effective technical method to obtain deep seafloor parameter.It has the beneficial effect that：
1) single hydrophone is utilized to obtain reverberation and propagation loss, you can carry out the inverting of bottom parameters, subsurface buoy scale very little,
It is possible to prevente effectively from the use of array.
2) Bottom sound speed, density and attenuation coefficient can be obtained simultaneously, improves the forecast precision of sound-field model, had important
Acoustic applications are worth.
Description of the drawings
Fig. 1 is Sound speed profile of the method for the present invention using CTD samplings and the experiment marine site of WOA09 database interfusions
Fig. 2 is the method for the present invention ocean space three-dimensional system of coordinate schematic diagram
Fig. 3 is that the method for the present invention carries out sensitivity analysis using 7 bottom parameters of bottom reverberation model pair
Fig. 4 is that the method for the present invention hydrophone depth is 820m, and sound source is 2.8583km with hydrophone distance, and centre frequency is
(a) reverb signal for being received when 1000Hz and (b) bottom reverberation signal
Fig. 5 is the comparison diagram of bottom reverberation intensity and experimental data that inversion result is brought into model calculating by the method for the present invention
(a) other frequencies of inverting frequency (b)
Fig. 6 is the scatter plot of refutation process of the method for the present invention based on bottom reverberation
Fig. 7 is the method for the present invention inverting value and experiment sampled value comparison diagram (a) density (b) velocity of sound
Fig. 8 is that the method for the present invention (a) includes that the propagation loss figure of landform and (b) utilize propagation loss to the quick of attenuation coefficient
Perceptual analysis
Fig. 9 is the inverting cost function value of different seabed attenuation coefficients of the method for the present invention based on propagation loss
Figure 10 is the method for the present invention propagation loss model different frequency calculated value and experiment value comparison diagram
Specific implementation mode
In conjunction with embodiment, attached drawing, the invention will be further described：
Single hydrophone cloth is put into water, receives the sound-source signal of transmitting, extracts bottom reverberation experimental data and modeling and simulating
Bottom reverberation time series is matched, and inverting obtains Bottom sound speed and density.By propagation loss experiment value and sound-field model
Calculated value is matched, and inverting obtains seabed attenuation coefficient.Its process is divided into following steps：
Step 1 single hydrophone cloth is put into water, and ship walks boat, during which with Fixed Time Interval to throwing explosive sound source in the sea,
Return behind traversal investigation marine site, fishes for hydrophone and obtains bottom reverberation and propagation loss experimental data.
Step 2 establishes three-dimensional system of coordinate in ocean space, and sound source S is located at (0,0, zs), it receives single hydrophone R and is located at (rr,
0,zr), the position of arbitrary scatterer is located at (xb,yb,0)。θijIt is the incident glancing angle of j-th of scatterer, θsjIt is scattering graze
Angle,It is the deflection of scattered wave.Seabed is assumed to be the semi-infinite half-space, and the advantages of being assumed using this seafloor model is that can have
Effect improves the efficiency that model calculates, and can reflect the average structure in seabed.Unknown bottom parameters share 7 under model
It is a：Seafloor density ρ, Bottom sound speed cp, seabed attenuation coefficient αλ, seabed surface roughness spectral intensity w2With spectrum index γ2, seabed
Inhomogeneities spectral intensity w3With spectrum index γ3。
Step 3 models deep seafloor reverberation.Bottom reverberation is the function decayed about the time, ignores seawater
Inhomogeneities, expression formula can be written as
RLmodel(f,zs,zr,tj)=Intensity (f, zs,zr,tj)+NL(f,zr) (8)
tj=tsj+trj (10)
Wherein, RLmodel(f,zs,zr,tj) it is theoretical reverberation level, Intensity (f, zs,zr,tj) it is reverrberation intensity；F is
The centre frequency of signal, NL are to receive depth as zr, frequency be f when seanoise intensity, can be obtained by experimental data.I0It is
Sound source level, cwIt is the velocity of sound in seawater, αwIt is the absorption of seawater, tsjAnd trjIt is sound source respectively to scattering point and scattering point to receiving
The time of hydrophone, σjIt is scattering section, δ AjIt is the insonify area of scatterer, tjIt is full from the relationship in different scatterer propagation times
Sufficient formula (4).
Step 4 carries out bandwidth to deep-sea reverberation experimental data and is filtered for 100Hz, and reverberation level is tested by following formula in seabed
It calculates：
RLexp(f,zs,zr,tj)=20log10(Veff(tj))-mv-10log10(Δf) (11)
Wherein, VeffIt is effective voltage value, mvIt is hydrophone sensitivity, Δ f is the frequency bandwidth of filter.
Step 5 is the nonuniqueness that cost function convergence speed is conciliate slowly when solving the problems, such as single-frequency independence inverting, using more
Frequency joint inversion, the cost function for Inversion for bottom density be,
Wherein, m is the number of the frequency for inverting, and n is the number of scatterer in reverberation modeling.To improve inverting effect
Rate is calculated accumulation least square and error under m frequency, is obtained object function using Parallel Particle Swarm Optimization genetic Optimization Algorithm
Minimum value, corresponding optimizing parameter optimal value is inverting optimal value at this time.Bottom sound speed is obtained by Hamilton empirical equations
It arrives：
cp+ 487.7 ρ of=2330.4-1257.0 ρ2 (13)
The seafloor density and the velocity of sound that step 6 obtains front inverting are used for the generation of inverting attenuation coefficient as Given information
Valence function is
Wherein, l is the number of explosive sound source, TLexp(f,zs,zr,rrk) it is propagation loss measured value, TLmodel(f,zs,zr,
rrk,αλ) be propagation loss model calculation value, when cost function is minimized be seabed attenuation coefficient inverting value.
The system configuration of wherein step 1 includes single hydrophone and loads explosive sound source or pull the ship of sound source.First will
Single hydrophone cloth is put into water, and ship navigates in investigation marine site into walking, while emission sound source signal, acquires deep-sea experimental data, carries
Take bottom reverberation and propagation loss experimental data.By bottom reverberation experimental data and the bottom reverberation time series of modeling and simulating into
Row matching, inverting obtain Bottom sound speed and density.Propagation loss experiment value is matched with the calculated value of sound field ray model,
Inverting obtains seabed attenuation coefficient.
The array of the arbitrary formation such as line array, ring array with multiple hydrophones, frequency range 100- can also be used
10kHz, hydrophone quantitative range are 1-100, and reception depth is 10-5000m, obtains deep-sea experimental data, utilizes each water
The inversion result of device is listened averagely to obtain bottom parameters.
The present embodiment Fig. 1 gives the Sound speed profile in the experiment marine site using CTD samplings and WOA09 database interfusions, water
Depth is 4008m.Selection and sound source warp, the immediate Climatological section of Position Latitude in WOA09 databases, and according to
Mackenzie experiential sound speed formulas calculate Sound speed profile (dotted line).The Sound speed profile sampled using CTD on sound channel axis
(solid line), using the Sound speed profile being calculated by WOA09 under CTD maximum sampling depths, sound channel axis is with CTD maximums using deep
The Sound speed profile (dotted line) merged using two sections between degree.Refutation process is completed by following six step：
(1) hydrophone cloth is put into water, and ship walks boat, is during which traversed with Fixed Time Interval to explosive sound source is thrown in the sea
Return behind investigation marine site fishes for hydrophone and obtains experimental data.
(2) ocean space three-dimensional system of coordinate establishes process：Fig. 2 is to establish three-dimensional system of coordinate in ocean space, and sound source S is located at
(0,0,zs), receiving hydrophone R is located at (rr,0,zr), the position of arbitrary scatterer is located at (xb,yb,0)。θijIt is j-th of scattering
The incident glancing angle of body, θsjIt is scattering glancing angle,It is the deflection of scattered wave.Seabed is assumed to be the semi-infinite half-space, using this
The advantages of kind seafloor model is assumed can effectively improve the efficiency of model calculating, and can reflect the average structure in seabed.
Unknown bottom parameters share 7 under model：Seafloor density ρ, Bottom sound speed cp, seabed attenuation coefficient αλ, seabed rough surface
Spend spectral intensity w2With spectrum index γ2, seabed inhomogeneities spectral intensity w3With spectrum index γ3。
(3) modeling process of deep seafloor reverberation：Bottom reverberation intensity is the function decayed about the time, ignores seawater
Inhomogeneities, expression formula can be written as
RLmodel(f,zs,zr,tj)=Intensity (f, zs,zr,tj)+NL(f,zr) (15)
tj=tsj+trj (17)
Wherein, f is the centre frequency of signal, and NL is to receive depth as zr, seanoise intensity of frequency when being f, can be by
Experimental data obtains.I0It is sound source level, cwIt is the velocity of sound in seawater, αwIt is the absorption of seawater, tsjAnd trjIt is sound source respectively to dissipating
Exit point and scattering point are to the time of receiving hydrophone, tjMeet formula (10), δ A from the relationship in different scatterer propagation timesjIt is
The insonify area of j-th of scatterer, can be found out by following formula：
Wherein, φ is the angle of sound source, scatterer and receiver, and L is the distance between sound source and receiver, and δ t are sound sources
Pulse width, δ φ are beam angles.σjIt is scattering section, can be found out by following formula：
Wherein,It is rough interface scattering section,It is volume scattering section.Its computational methods
Referring to " High-frequency bottom backscattering:Roughness versus sediment volume
Scattering ", this article in August, 1992 are published in《The Journal of the Acoustical Society of
America》92nd phase, first page number 962.Fig. 3 is to carry out sensitivity analysis to bottom parameters using bottom reverberation model.By
For sensitivity curves analysis result it is found that bottom reverberation is insensitive to attenuation coefficient, the sensibility of seafloor density and the velocity of sound is all relatively strong.
(4) bottom reverberation Data Processing in Experiment：To deep-sea reverberation experimental data effective voltage value VeffCarry out bandwidth 100Hz
To be filtered, bottom reverberation signal is calculated by following formula：
RLexp(f,zs,zr,tj)=20log10(Veff(tj))-mv-10log10(Δf) (20)
Wherein, mvIt is hydrophone sensitivity, Δ f is the frequency bandwidth of filter.Fig. 4 is that hydrophone depth is 820m, sound
Source is 2.8583km with hydrophone distance, and obtained reverb signal and sea are handled according to formula (8) when centre frequency is 1000Hz
Bottom reverb signal, the as can be seen from FIG. much more apparent way of deep-sea underwater sound signal and reverberation trailing phenomenon.By model and experimental data
The comparison of more way delay inequalitys it is available, direct wave (D) and the delay inequality of sea surface reflection wave (SR) are closer to, and what is respectively followed mixes
Sound overlaps, and what is reached successively later is bottom echo (BR), sea bottom echo (SBR), bottom and surface of sea back wave
(BSR) and sea bottom and surface of sea back wave (SBSR) and the reverb signal that respectively follows.
(5) seafloor density and Bottom sound speed refutation process：Cost function convergence speed is slow when to solve single-frequency independence inverting
The nonuniqueness problem of reconciliation, using multifrequency joint inversion, the cost function for Inversion for bottom density is,
Wherein, m is the number of the frequency for inverting, and n is the number of scatterer in reverberation modeling.To improve inverting effect
Rate is calculated accumulation least square and error under m frequency, is obtained object function using Parallel Particle Swarm Optimization genetic Optimization Algorithm
Minimum value, corresponding optimizing parameter optimal value is inverting optimal value at this time.Bottom sound speed is obtained by Hamilton empirical equations
It arrives：
cp+ 487.7 ρ of=2330.4-1257.0 ρ2 (22)
Fig. 5 is to bring inversion result bottom reverberation intensity and the experimental data of model calculating into inverting frequency and other frequencies
The comparison diagram of rate, the two are consistent substantially.Fig. 6 is the scatter plot of the refutation process based on bottom reverberation, it can be seen that inverted parameters to
The convergent process of inverting optimal value.Fig. 7 is inverting value and experiment sampled value comparison diagram.As a result illustrate, Bottom sound speed and density with
Experiment sampled value is consistent substantially, has higher confidence level.
(6) refutation process of seabed attenuation coefficient：Fig. 8 be comprising landform propagation loss figure and using propagation loss to declining
The sensitivity analysis for subtracting coefficient, it is insensitive because there is the presence of seamount to decay seabed in 20km, decay seabed in 20-
Decay to seabed more sensitive in the first and second shadow zones in 100km distances.Cost function for inverting attenuation coefficient is
Wherein, l is the number of explosive sound source, TLexp(f,zs,zr,rrk) it is propagation loss measured value, TLmodel(f,zs,zr,
rrk,αλ) be propagation loss model calculation value.Fig. 9 is the inverting cost letter of the different seabed attenuation coefficients based on propagation loss
Numerical value, attenuation coefficient when cost function is minimized are attenuation coefficient inverting value.Figure 10 is to bring inverting value into sound field mould
Type obtains the comparison of model calculation value (dot) and experimental measurements (solid line), due to the presence of seamount, the sea within the scope of this
Bottom parameter is to apart from related, and calculated value and measured value mismatch within 10-20km, landform is relatively flat other than 20km, calculated value with
Measured value coincide substantially.
The present invention achieves apparent implementation result in an exemplary embodiment, in situ measurement compared with sea floor sampling, behaviour
Make simple, easy to implement, compared with existing traditional Matched Field inversion method, complicated array, superiority need not be used
It is that inverting can be realized using single hydrophone, and obtains Bottom sound speed, density and attenuation coefficient simultaneously, to obtains deep seafloor
Parameter provides effective technical method, and there are important acoustic applications to be worth.
Claims (3)
1. a kind of deep seafloor parameter inversion method based on bottom reverberation and propagation loss, it is characterised in that steps are as follows：
Step 1：Single hydrophone cloth is put into water, ship walks boat, is adjusted with Fixed Time Interval to explosive sound source, traversal is thrown in the sea
Return behind marine site is ground, single hydrophone is fished for and obtains bottom reverberation and propagation loss experimental data；
Step 2：Establish seafloor model：Three-dimensional system of coordinate is established in ocean space, wherein sound source S is located at (0,0, zs), single hydrophone
R is located at (rr,0,zr), the position of arbitrary scatterer is located at (xb,yb, 0), θijIt is the incident glancing angle of j-th of scatterer, θsjIt is
Glancing angle is scattered,It is the deflection of scattered wave；7 bottom parameters：Seafloor density ρ, Bottom sound speed cp, seabed attenuation coefficient
αλ, seabed surface roughness spectral intensity w2With spectrum index γ2, seabed inhomogeneities spectral intensity w3With spectrum index γ3；
Step 3：Seafloor model based on step 2 carries out theoretical modeling to deep seafloor reverberation：
RLmodel(f,zs,zr,tj)=Intensity (f, zs,zr,tj)+NL(f,zr) (1)
tj=tsj+trj (3)
Wherein, RLmodel(f,zs,zr,tj) it is theoretical reverberation level, Intensity (f, zs,zr,tj) it is reverrberation intensity；F is signal
Centre frequency, NL be receive depth be zr, frequency be f when seanoise intensity, I0It is sound source level, cwIt is the sound in seawater
Speed, αwIt is the absorption of seawater, tsjAnd trjIt is sound source respectively to scattering point and scattering point to the time of receiving hydrophone, σjIt is scattering
Section, δ AjIt is the insonify area of scatterer, tjMeet formula (3) from the relationship in different scatterer propagation times；
Step 4：Bandwidth is carried out to deep seafloor reverberation experimental data to be filtered for 100Hz：
RLexp(f,zs,zr,tj)=20log10(Veff(tj))-mv-10log10(Δf) (4)
Wherein, RLexp(f,zs,zr,tj) it is experiment reverberation level, VeffIt is effective voltage value, mvIt is hydrophone sensitivity, Δ f is filter
The frequency bandwidth of wave device；
Step 5：Using multifrequency joint inversion, the cost function for Inversion for bottom density is established：
Wherein, m is the number of the frequency for inverting, and n is the number of scatterer in reverberation modeling；Using Parallel Particle Swarm Optimization heredity
Optimization algorithm calculates accumulation least square and error under m frequency, obtains the minimum value of object function, corresponding at this time to seek
Excellent parameter optimal value is inverting optimal value；Bottom sound speed is obtained by Hamilton empirical equations：
cp+ 487.7 ρ of=2330.4-1257.0 ρ2 (6)
Step 6：The seafloor density and Bottom sound speed that step 5 inverting is obtained are established as Given information for inverting decaying system
Several cost functions：
Wherein, l is the number of explosive sound source, TLexp(f,zs,zr,rrk) it is propagation loss measured value, TLmodel(f,zs,zr,rrk,
αλ) be propagation loss model calculation value, when cost function is minimized be seabed attenuation coefficient inverting value.
2. a kind of deep seafloor parameter inversion method based on bottom reverberation and propagation loss according to claim 1,
It is characterized in that single hydrophone can be replaced the array being made of multiple hydrophones.
3. a kind of deep seafloor parameter inversion method based on bottom reverberation and propagation loss according to claim 2,
It is characterized in that hydrophone quantitative range is 1-100, forms line array or ring array.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610478946.8A CN106154276B (en) | 2016-06-27 | 2016-06-27 | Deep seafloor parameter inversion method based on bottom reverberation and propagation loss |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610478946.8A CN106154276B (en) | 2016-06-27 | 2016-06-27 | Deep seafloor parameter inversion method based on bottom reverberation and propagation loss |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106154276A CN106154276A (en) | 2016-11-23 |
| CN106154276B true CN106154276B (en) | 2018-10-19 |
Family
ID=57350141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610478946.8A Active CN106154276B (en) | 2016-06-27 | 2016-06-27 | Deep seafloor parameter inversion method based on bottom reverberation and propagation loss |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106154276B (en) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106441543A (en) * | 2016-12-09 | 2017-02-22 | 华南理工大学 | Underwater detection path acoustic velocity measurement method based on three-dimensional orthogonal array and underwater detection path acoustic velocity measurement device thereof |
| CN107656243A (en) * | 2017-08-25 | 2018-02-02 | 天津大学 | Combine DOA/TOA oceans multi-path environment localization method in inhomogeneous medium |
| CN108398690B (en) * | 2017-12-08 | 2020-08-11 | 国家深海基地管理中心 | Submarine backscattering intensity measuring method |
| CN109975812B (en) * | 2017-12-27 | 2020-10-23 | 中国科学院声学研究所 | Method for calculating sound source depth and sound source level of explosive sound |
| CN109556701A (en) * | 2018-11-01 | 2019-04-02 | 浙江海洋大学 | A kind of shallow sea geoacoustic inversion method based on broadband vertical wave impedance |
| CN109724684B (en) * | 2019-01-03 | 2020-10-13 | 武汉大学 | Direct signal propagation time measuring method based on underwater autonomous vehicle |
| CN109725053B (en) * | 2019-01-22 | 2021-03-12 | 中国人民解放军国防科技大学 | Method for acquiring water sound field characteristic data based on high-resolution ocean re-analysis product |
| CN110542919A (en) * | 2019-07-10 | 2019-12-06 | 自然资源部第三海洋研究所 | Sea sand resource detection method based on electric spark source seismic data inversion |
| CN110554394B (en) * | 2019-09-05 | 2021-04-20 | 杭州瑞利海洋装备有限公司 | ADCP flow measurement signal-based TVC curve accurate setting method |
| CN110568427B (en) * | 2019-09-05 | 2021-04-16 | 杭州瑞利海洋装备有限公司 | Reverberation parameter inversion method based on ADCP flow measurement signal |
| CN110824485A (en) * | 2019-09-25 | 2020-02-21 | 中国人民解放军军事科学院国防科技创新研究院 | Bistatic sonar underwater target feature extraction and classification method based on AUV |
| CN111458008B (en) * | 2020-04-16 | 2021-03-26 | 浙江大学 | Ocean earth sound parameter inversion method based on single fixed hydrophone and single mobile sound source |
| CN111650159B (en) * | 2020-06-17 | 2023-07-18 | 自然资源部第一海洋研究所 | Sea surface backscattering strength measuring method |
| CN113984179B (en) * | 2021-10-25 | 2023-09-19 | 中国人民解放军海军工程大学 | Static ocean reverberation measurement acquisition method and device |
| CN115616073B (en) * | 2022-12-19 | 2023-05-16 | 广州海洋地质调查局三亚南海地质研究所 | Method and device for estimating sea grass bed oxygen release amount by using hydrophone |
| CN116859460B (en) * | 2023-09-04 | 2023-11-21 | 自然资源部第一海洋研究所 | Submarine physical property parameter inversion method suitable for elastic wave near-field reflection |
| CN117629382B (en) * | 2024-01-25 | 2024-04-19 | 山东科技大学 | Underwater sound velocity and topography distance measuring method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183435A (en) * | 2011-01-25 | 2011-09-14 | 中国船舶重工集团公司第七一五研究所 | Method for measuring submarine density and sound velocity based on multi-path reflection theory |
| CN103487793A (en) * | 2013-09-22 | 2014-01-01 | 中国人民解放军海军工程大学 | Broadband reverberation waveform simulation method based on normal mode theory |
-
2016
- 2016-06-27 CN CN201610478946.8A patent/CN106154276B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183435A (en) * | 2011-01-25 | 2011-09-14 | 中国船舶重工集团公司第七一五研究所 | Method for measuring submarine density and sound velocity based on multi-path reflection theory |
| CN103487793A (en) * | 2013-09-22 | 2014-01-01 | 中国人民解放军海军工程大学 | Broadband reverberation waveform simulation method based on normal mode theory |
Non-Patent Citations (4)
| Title |
|---|
| Bistatic bottom reverberation in deep ocean: modeling and data comparison;Liya Xu 等;《IEEE》;20160413;1-5 * |
| 利用海底反射信号进行地声参数反演的方法;杨坤德 等;《物理学报》;20090331;第58卷(第3期);1798-1805 * |
| 深海海底混响模型初步研究;翁晋宝 等;《声学技术》;20141231;第33卷(第S2期);67-69 * |
| 深海海底混响的概率密度分布特性研究;徐丽亚 等;《声学技术》;20151231;第34卷(第6期);100-103 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106154276A (en) | 2016-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106154276B (en) | Deep seafloor parameter inversion method based on bottom reverberation and propagation loss | |
| CN108181651B (en) | A kind of deep-sea geoacoustic inversion method based on reliable acoustic path | |
| CN107247291B (en) | The shallow stratum sound energy attenuation model construction in seabed and two important sound energy attenuation characteristic parameter extraction method | |
| RU2603724C2 (en) | Method and device to control acoustic characteristics of network of acoustic nodes located along towed acoustic linear antennae | |
| Wang et al. | Review of underwater acoustic propagation models. | |
| Knobles et al. | Seabed acoustics of a sand ridge on the New Jersey continental shelf | |
| Seri et al. | Long range passive ocean acoustic waveguide remote sensing (poawrs) of seismo-acoustic airgun signals received on a coherent hydrophone array | |
| Rutenko et al. | A method for estimating the characteristics of acoustic pulses recorded on the sakhalin shelf for multivariate analysis of their effect on the behavior of gray whales | |
| WO2013055236A1 (en) | System and method for modelling multibeam backscatter from a seafloor | |
| Yu et al. | Measurements of midfrequency acoustic backscattering from a sandy bottom in the South Yellow Sea of China | |
| Kozaczka et al. | Processing data on sea bottom structure obtained by means of the parametric sounding | |
| Rui et al. | Application of Sub-Bottom Profiler to Study Riverbed Structure and Sediment Density | |
| Moon et al. | Imaging the yellow sea bottom cold water from multichannel seismic data | |
| Holland | Mapping seabed variability: Rapid surveying of coastal regions | |
| Ohkawa et al. | Acoustic backscattering from a sandy seabed | |
| Barlee et al. | Geoacoustic model parameter estimation using a bottom-moored hydrophone array | |
| Leighton et al. | Acoustic propagation in gassy intertidal marine sediments: An experimental study | |
| Crawford et al. | Development of a sound propagation model to monitor seismic survey noise in Irish waters | |
| Chen et al. | Effects of nonlinear internal gravity waves on normal-incident reflection measurements of seafloor sediments | |
| RU2814150C1 (en) | Parametric profiler | |
| Xing et al. | Construction of the Central Arctic Sea Ice Structure and Acoustic Velocity Model at the Short-Term Ice Station During N11 CHINARE | |
| Li et al. | Modeling of bottom backscattering from three-dimensional volume inhomogeneities and comparisons with experimental data | |
| RU126146U1 (en) | MULTI-BEAM Echo Sounder | |
| Belov et al. | Acoustic calibration of the signal propagation path in the shallow water | |
| Uzhansky et al. | Using a kilometers-long horizontal multichannel array for estimating seafloor sound speed in a passive margin setting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |