CN110398744A - Ocean thermocline characteristic parameter optimizing and inverting method based on acoustic signals - Google Patents
Ocean thermocline characteristic parameter optimizing and inverting method based on acoustic signals Download PDFInfo
- Publication number
- CN110398744A CN110398744A CN201910684576.7A CN201910684576A CN110398744A CN 110398744 A CN110398744 A CN 110398744A CN 201910684576 A CN201910684576 A CN 201910684576A CN 110398744 A CN110398744 A CN 110398744A
- Authority
- CN
- China
- Prior art keywords
- signal
- thermocline
- characteristic parameter
- temperature
- formula
- 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.)
- Pending
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
- 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
- G01S15/885—Meteorological systems
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention provides an ocean thermocline characteristic parameter optimizing and inverting method based on acoustic signals, which comprises the following steps: acquiring the vertical profile distribution of the temperature and salinity of the marine environment; constructing an acoustic echo prediction model based on an acoustic reflection principle, and obtaining an acoustic echo prediction signal according to the acoustic echo prediction model; constructing a thermocline characteristic parameter inversion model by using the temperature vertical profile distribution, the salinity vertical profile distribution and the acoustic echo prediction signal of the marine environment; and solving the constructed thermocline characteristic parameter inversion model to obtain the characteristic parameters of the ocean thermocline. The ocean thermocline characteristic parameter optimizing and inverting method effectively solves the problem of low inversion efficiency of traditional acoustic methods, and provides an effective method for efficient inversion of the ocean thermocline.
Description
Technical field
The present invention relates to a kind of thermocline layer characteristic parameter optimizing inversion method, especially a kind of seas based on acoustical signal
Foreign thermocline characteristic parameter optimizing inversion method.
Background technique
Thermocline layer as a kind of important marine environmental phenomenon, for climate change, sea fishery, underwater communication and
Submarine activity has great influence.Obtaining thermocline information using conventional conventional method has the deficiency that time-consuming, at high cost,
Then a kind of effective method will be become by obtaining thermocline layer characteristic parameter using acoustic method.Therefore, it is necessary to design
A kind of thermocline layer characteristic parameter optimizing inversion method based on acoustical signal out, can be by constructing the ocean based on acoustical signal
Thermocline parameter optimization inverting Optimized model, and solved using population-genetic fusion algorithm, obtain thermocline layer spy
Levy parameter.
Summary of the invention
It is an object of the invention to: a kind of thermocline layer characteristic parameter optimizing inversion method based on acoustical signal is provided,
It can be by constructing the thermocline layer parameter optimization inverting Optimized model based on acoustical signal, and use population-genetic fusion
Algorithm is solved, and thermocline layer characteristic parameter is obtained.
In order to achieve the above-mentioned object of the invention, the present invention provides a kind of, and the thermocline layer characteristic parameter based on acoustical signal is sought
Excellent inversion method, includes the following steps:
Step 1, the vertical section distribution of the temperature and salinity of marine environment is obtained;
Step 2, it is based on sound reflecting principle building sound echo prediction model, it is pre- according to sound echo prediction model acquisition sound echo
Survey signal;
Step 3, believed using the temperature vertical section distribution of marine environment, the distribution of salinity vertical section and sound echo prediction
Number building thermocline characteristic parameter inverse model;
Step 4, the thermocline characteristic parameter inverse model of building is solved, obtains the feature ginseng of thermocline layer
Number.
Further, in step 1, the temperature vertical section of marine environment is distributed expression are as follows:
In formula (1), z is the depth of water, αiFor temperature EOF coefficient, T0It (z) is mean temperature section, fiIt (z) is feature vector.
Further, mean temperature section T0(z) calculation formula are as follows:
In formula (2), N is the number of temperature profile, tiIt (z) is the temperature value of different sections.
Further, feature vector fi(z) calculating process are as follows:
First by obtaining the anomaly value of temperature profile:
ΔTi(z)=ti(z)-T0(z) (3)
Then anomaly value is formed into matrix X:
The covariance matrix S of calculating matrix X again:
S=XXT (5)
Eigenvalues Decomposition is done to covariance matrix S again and obtains feature value vector σ, then eigenvectors matrix F has:
SF=F σ (6)
F=[f1(z), f2(z) ..., fN(z)] (7)
σ=[σ1, σ2..., σN] (8)
Significance test finally is carried out to the feature vector σ decomposited, thus each characteristic value in judging characteristic vector σ
There are physical significance signal or noise signal, examine formula are as follows:
When formula (9) are set up, then characteristic value σkAnd σk+1It is separable, then after testing to Eigenvalues Decomposition, by feature
Value is ranked up by sequence from big to small, and selected characteristic is worth the corresponding feature vector f of biggish preceding m rank characteristic valuei(z)。
Further, the salinity vertical section distribution of marine environment is to carry out fitting of a polynomial to salinity according to T-S relationship,
Fitting formula are as follows:
In formula (10), z is the depth of water, PkFor salinity EOF coefficient, T (z) is the vertical section temperature value of temperature.
Further, the specific steps in step 2, when according to sound echo prediction model acquisition sound echo prediction signal
Are as follows:
Building sound echo prediction model first are as follows:
Y (n)=h (n) * x (n) (11)
In formula (11), h (n) is the impulse response of briny environment medium, i.e., received when incoming signal is δ pulse
Echo record, impulse response h (n) includes the characteristic information of thermocline layer structure, and y (n) is the sound echo-signal received,
X (n) is incident acoustical signal, and * indicates convolution;
Then the layered medium model of the times thickness such as foundation is respectively layered when input signal is impulse signal δ (t) again
Impulse response signals indicate are as follows:
In formula (12), f (0,0)=1, f (n, 0)=1, f (n, n)=r0rn, n >=1, riFor the acoustical reflection factor of each layering,
D (0)=δ (t) is impulse signal, and the h (n) being calculated, which is updated to formula (11) just, can obtain sound echo prediction signal.
Further, in step 3, the thermocline characteristic parameter inverse model of building are as follows:
In formula (14), T (z) and S (z) are respectively two elements of temperature and salinity of inverting, and T (z) and S (z) are water
The function of deep z, H are the number of samples of sound echo-signal, and i is the sequence number of sound echo signal sample, yiIndicate what prediction obtained
Sound echo-signal, for the convolution of incident acoustical signal and briny environment medium impulse response, yi' indicate that the sound that actual measurement obtains returns
Wave signal.
Further, in step 4, when solving to the thermocline characteristic parameter inverse model of building, population-is utilized
Genetic fusion algorithm carries out optimizing iterative solution to thermocline characteristic parameter inverse model, obtains temperature profile parameter and salinity is cutd open
Face parameter calculates Sound speed profile parameter further according to temperature profile parameter and salt profile parameter are as follows:
C=1449.14+ Δ cT+ΔcS+ΔcP+ΔcSTP (17)
In formula:
ΔcT=4.5721T-4.4532 × 10-2T2-2.6045×10-4T3+7.985×10-6T4
ΔcS=1.3980 (S-35)+1.692 × 10-3(S-35)2
ΔcP=1.60272 × 10-1P+1.0268×10-5P2+3.5216×10-9P3-3.3603×10-12P4
ΔcSTP=(S-35) (- 1.1244 × 10-2T+7.7711×10-7T2+7.7016×10-5P-1.2943×10-7P2
+3.1580×10-8PT+1.5790×10-9PT2)+P(-1.8607×10-4T+7.4812×10-6T2+4.5283×10-8T3)+
P2(-2.5294×10-7T+1.8563×10-9T2)+P3(-1.9646×10-10T)
P=1.033+1.028126 × 10-1Z+2.38×10-7Z2-6.8×10-17Z4
T is temperature, and unit is degree Celsius;S is salinity, unit ‰;P is pressure, unit dimension gram/cm;Z is deep
Degree, unit is rice.
The beneficial effects of the present invention are: it is anti-that thermocline layer characteristic parameter high efficiency is realized based on sound echo-signal
It drills, solves the problems, such as that conventional acoustic method extracts inverting low efficiency, provide one for the efficient inverting of thermocline layer
The effective method of kind.
Detailed description of the invention
Fig. 1 is inversion method flow chart of the invention;
Fig. 2 is reflection echo schematic illustration of the invention;
Fig. 3 is particle of the invention-genetic dynamics parallel-melt method flow chart.
Specific embodiment
Technical solution of the present invention is described in detail with reference to the accompanying drawing, but protection scope of the present invention is not limited to
The embodiment.
Embodiment 1:
As shown in Figure 1, the thermocline layer characteristic parameter optimizing inversion method disclosed by the invention based on acoustical signal, including
Following steps:
Step 1, the vertical section distribution of the temperature and salinity of marine environment is obtained;
Step 2, it is based on sound reflecting principle building sound echo prediction model, it is pre- according to sound echo prediction model acquisition sound echo
Survey signal;
Step 3, believed using the temperature vertical section distribution of marine environment, the distribution of salinity vertical section and sound echo prediction
Number building thermocline characteristic parameter inverse model;
Step 4, the thermocline characteristic parameter inverse model of building is solved, obtains the feature ginseng of thermocline layer
Number.
Sound wave is influenced maximum to be temperature in ocean, followed by salinity.In common inversion, usually with the velocity of sound
For inverting object, and the basic parameter of Seawater is temperature and salinity etc., can conveniently be sent using temperature and salinity
Livelihood calculates the physical quantity of other care, therefore inversion result can be summed up in the point that the basic parameter of Seawater itself as much as possible
Form;Meanwhile the unknown number of inverse model should lack as best one can, to improve the efficiency of inverting.Based on these considerations, to ocean
On the basis of temperature profile and salt profile property research, temperature profile can be characterized using Empirical Orthogonal Function method and salinity is cutd open
EDS maps, therefore, in step 1, the temperature vertical section of marine environment is distributed expression are as follows:
In formula (1), z is the depth of water, αiFor temperature EOF coefficient, T0It (z) is mean temperature section, fiIt (z) is feature vector.
Further, mean temperature section T0(z) calculation formula are as follows:
In formula (2), N is the number of temperature profile, tiIt (z) is the temperature value of different sections.
Further, feature vector fi(z) calculating process are as follows:
First by obtaining the anomaly value of temperature profile:
ΔTi(z)=ti(z)-T0(z) (3)
Then anomaly value is formed into matrix X:
The covariance matrix S of calculating matrix X again:
S=XXT (5)
Eigenvalues Decomposition is done to covariance matrix S again and obtains feature value vector σ, then eigenvectors matrix F has:
SF=F σ (6)
F=[f1(z), f2(z) ..., fN(z)] (7)
σ=[σ1, σ2..., σN] (8)
Significance test finally is carried out to the feature vector σ decomposited, thus each characteristic value in judging characteristic vector σ
There are physical significance signal or noise signal, examine formula are as follows:
When formula (9) are set up, then characteristic value σkAnd σk+1Be it is separable, characteristic value has physical significance, to characteristic value
After decomposition is tested, characteristic value is ranked up by sequence from big to small, selected characteristic is worth biggish preceding m rank characteristic value pair
The feature vector f answeredi(z)。
Further, the salinity vertical section distribution of marine environment is to carry out fitting of a polynomial to salinity according to T-S relationship,
Fitting formula are as follows:
In formula (10), z is the depth of water, PkFor salinity EOF coefficient, T (z) is the vertical section temperature value of temperature.It establishes with temperature
Degree and salinity are the marine environment dielectric model of major parameter, pass through experimental analysis, temperature, seven rank experiences of saline environment parameter
Orthogonal function more can accurately describe temperature, the vertical distribution of salinity.
Have the characteristics that heterogeneity for seawater, then use for reference layered medium model, a small range can be by seawater
It is approximately layered medium model, as shown in Fig. 2, thinking is uniform in Seawater horizontal direction, and in vertical direction
Using etc. whens thickness model, seawater is considered as to the fine and closely woven plane layer of U equal times thickness, the inner parameter in each layer is identical
, i.e. conforming layer;And between every layer it is mutation, in each layering transmission time Δ t having the same, rather than identical sky
Between thickness deltat x.However, this and real medium are that continuous situation has a certain difference, if by fine and closely woven plane layer thickness control
System is in incidence wave wavelengthThe reflection characteristics of so layered medium are consistent with continuous media, if layering number U obtains foot
Enough big, then thickness layering and actual seawater section are close whens waiting.
For primary detection, variation of the Seawater in time scale is much smaller than the variation on space scale, thus,
Can be by Seawater as a linear time invariant system, the reflection echo of thermocline layer is considered as incident acoustic wave and sea
The convolution of the impulse response of foreign environmental parameter.Then in step 2, according to sound echo prediction model acquisition sound echo prediction signal
When specific steps are as follows:
Building sound echo prediction model first are as follows:
Y (n)=h (n) * x (n) (11)
In formula (11), h (n) is the impulse response of briny environment medium, i.e., received when incoming signal is δ pulse
Echo record, impulse response h (n) includes the characteristic information of thermocline layer structure, and y (n) is the sound echo-signal received,
X (n) is incident acoustical signal, and * indicates convolution;
Then the layered medium model of the times thickness such as foundation is respectively layered when input signal is impulse signal δ (t) again
Impulse response signals indicate are as follows:
In formula (12), f (0,0)=1, f (n, 0)=1, f (n, n)=r0rn, n >=1, riFor the acoustical reflection factor of each layering,
D (0)=δ (t) is impulse signal, and the h (n) being calculated, which is updated to formula (11) just, can obtain sound echo prediction signal.
Since sound wave propagates the mainly influence by ocean environment parameters such as temperature, salinity, present invention building in the seawater
Thermocline layer parametric inversion mathematical model, using temperature, salinity as research object, essence is namely based on to be worked as inverting object
Simulate inverting sound analogue echoes signal and measurement field sound echo-signal closest to when, it is believed that the thermocline parameter of simulation with
Actual measurement field thermocline layer parameter optimal approximation, therefore, in step 3, the thermocline characteristic parameter inverse model of building are as follows:
In formula (14), T (z) and S (z) are respectively two elements of temperature and salinity of inverting, and T (z) and S (z) are water
The function of deep z, H are the number of samples of sound echo-signal, and i is the sequence number of sound echo signal sample, yiIndicate what prediction obtained
Sound echo-signal, for the convolution of incident acoustical signal and briny environment medium impulse response, yi' indicate that the sound that actual measurement obtains returns
Wave signal.
Further, in step 4, when solving to the thermocline characteristic parameter inverse model of building, population-is utilized
Genetic fusion algorithm carries out optimizing iterative solution to thermocline characteristic parameter inverse model, obtains temperature profile parameter and salinity is cutd open
Face parameter calculates Sound speed profile parameter further according to temperature profile parameter and salt profile parameter are as follows:
C=1449.14+ Δ cT+ΔcS+ΔcP+ΔcSTP (17)
In formula:
ΔcT=4.5721T-4.4532 × 10-2T2-2.6045×10-4T3+7.985×10-6T4
ΔcS=1.3980 (S-35)+1.692 × 10-3(S-35)2
ΔcP=1.60272 × 10-1P+1.0268×10-5P2+3.5216×10-9P3-3.3603×10-12P4
ΔcSTP=(S-35) (- 1.1244 × 10-2T+7.7711×10-7T2+7.7016×10-5P-1.2943×10-7P2
+3.1580×10-8PT+1.5790×10-9PT2)+P(-1.8607×10-4T+7.4812×10-6T2+4.5283×10-8T3)+
P2(-2.5294×10-7T+1.8563×10-9T2)+P3(-1.9646×10-10T)
P=1.033+1.028126 × 10-1Z+2.38×10-7Z2-6.8×10-17Z4
C is the velocity of sound, and unit is meter per second;T is temperature, and unit is degree Celsius;S is salinity, unit ‰;P is pressure, single
Position dimension gram/cm;Z is depth, and unit is rice.
As shown in figure 3, further, population-genetic fusion algorithm basic principle are as follows: the individual in group is each
In secondary optimizing iteration journey, firstly, selecting the big individual of certain amount fitness value using sequencing selection method carries out crossover operation
The next generation is generated, then, mutation operation is carried out according to certain mutation probability to the filial generation that crossover operation generates, obtains filial generation
Position and speed, meanwhile, for being not carried out father's individual of crossover operation still according to particle swarm algorithm speed and location updating
Formula generates the next generation, and the optimal location of optimal location and group individual in two kinds of modes of operation is shared.Population-something lost
Pass the execution step of blending algorithm are as follows:
Step (I), to position and speed individual in population size, initial population, weight factor, maximum speed, intersection
Probability, mutation probability and maximum number of iterations are initialized, and individual position and speed is to be randomly generated in initial population
Temperature EOF coefficient and salinity EOF coefficient, initialization population size is set as 50, weight factor 1.0, maximum speed 5.0,
Crossover probability is 0.8, mutation probability 0.05;
Step (II) calculates the functional value of each individual according to formula (14), chooses the current function value of an individual,
And be compared with history optimal function value, if current function value is higher than history optimal function value, updated with current function value
History optimal function value, conversely, then retaining history optimal function value;
Step (III) one by one compares the functional value of the current function value of each individual and group's overall situation optimum position
Compared with if current function value is higher, global optimum position being updated with the position of current function value, conversely, then retaining original
Global optimum position;
Step (IV) is ranked up all individuals according to the size of functional value, 2M optimal individual is selected to be handed over
Fork and mutation operator;
Step (V), the individual in group in addition to 2M individual carry out more according to population position and speed more new formula
Newly;
Step (VI), if meeting the termination condition of maximum number of iterations, optimizing iteration terminates, otherwise, return step
(Ⅱ)。
Further, in step (V), when carrying out crossing operation, individual adaptation degree is ranked up, obtains coming front
2M individual, by 2M it is individual match two-by-two, the formula of crossing operation are as follows:
In formula (15),Speed and the position of filial generation are represented,Represent the speed of parent
Degree and position, k are the number of iteration, and random number of the r between [0,1] is indicated with the random number for meeting Gaussian Profile;
In step (V), the formula of mutation operator are as follows:
In formula (16),Speed and the position of filial generation are represented,Speed and the position of parent are represented, k is repeatedly
The number in generation, CiFor the random number between [0,1].
As described above, must not be explained although the present invention has been indicated and described referring to specific preferred embodiment
For the limitation to invention itself.It without prejudice to the spirit and scope of the invention as defined in the appended claims, can be right
Various changes can be made in the form and details for it.
Claims (8)
1. a kind of thermocline layer characteristic parameter optimizing inversion method based on acoustical signal, which comprises the steps of:
Step 1, the vertical section distribution of the temperature and salinity of marine environment is obtained;
Step 2, it is based on sound reflecting principle building sound echo prediction model, is believed according to sound echo prediction model acquisition sound echo prediction
Number;
Step 3, the temperature vertical section distribution of marine environment, the distribution of salinity vertical section and sound echo prediction signal structure are utilized
Build thermocline characteristic parameter inverse model;
Step 4, the thermocline characteristic parameter inverse model of building is solved, obtains the characteristic parameter of thermocline layer.
2. the thermocline layer characteristic parameter optimizing inversion method according to claim 1 based on acoustical signal, feature exist
In in step 1, the temperature vertical section of marine environment is distributed expression are as follows:
In formula (1), z is the depth of water, αiFor temperature EOF coefficient, T0It (z) is mean temperature section, fiIt (z) is feature vector.
3. the thermocline layer characteristic parameter optimizing inversion method according to claim 2 based on acoustical signal, feature exist
In mean temperature section T0(z) calculation formula are as follows:
In formula (2), N is the number of temperature profile, tiIt (z) is the temperature value of different sections.
4. the thermocline layer characteristic parameter optimizing inversion method according to claim 2 based on acoustical signal, feature exist
In feature vector fi(z) calculating process are as follows:
First by obtaining the anomaly value of temperature profile:
ΔTi(z)=ti(z)-T0(z) (3)
Then anomaly value is formed into matrix X:
The covariance matrix S of calculating matrix X again:
S=XXT (5)
Eigenvalues Decomposition is done to covariance matrix S again and obtains feature value vector σ, then eigenvectors matrix F has:
SF=F σ (6)
F=[f1(z), f2(z) ..., fN(z)] (7)
σ=[σ1,σ2..., σN] (8)
Significance test finally is carried out to the feature vector σ decomposited, so that each characteristic value in judging characteristic vector σ is that have
Physical significance signal or noise signal examine formula are as follows:
When formula (9) are set up, then characteristic value σkAnd σk+1It is separable, then after testing to Eigenvalues Decomposition, characteristic value is pressed
Sequence from big to small is ranked up, and selected characteristic is worth the corresponding feature vector f of biggish preceding m rank characteristic valuei(z)。
5. the thermocline layer characteristic parameter optimizing inversion method according to claim 2 based on acoustical signal, feature exist
In the salinity vertical section distribution of marine environment is to carry out fitting of a polynomial, fitting formula to salinity according to T-S relationship are as follows:
In formula (10), z is the depth of water, PkFor salinity EOF coefficient, T (z) is the vertical section temperature value of temperature.
6. the thermocline layer characteristic parameter optimizing inversion method according to claim 5 based on acoustical signal, feature exist
In specific steps in step 2, when according to sound echo prediction model acquisition sound echo prediction signal are as follows:
Building sound echo prediction model first are as follows:
Y (n)=h (n) * x (n) (11)
In formula (11), h (n) is the impulse response of briny environment medium, i.e. the echo received when incoming signal is δ pulse
Record, impulse response h (n) includes the characteristic information of thermocline layer structure, and y (n) is the sound echo-signal received, x (n)
For incident acoustical signal, * indicates convolution;
Then again according to etc. times thickness layered medium model, when input signal be impulse signal δ (t) when, what is be respectively layered rushes
Swashing response signal indicates are as follows:
In formula (12), f (0,0)=1, f (n, 0)=1, f (n, n)=r0rn, n >=1, riFor the acoustical reflection factor of each layering, D (0)
=δ (t) is impulse signal, and the h (n) being calculated, which is updated to formula (11) just, can obtain sound echo prediction signal.
7. the thermocline layer characteristic parameter optimizing inversion method according to claim 1 based on acoustical signal, feature exist
In, in step 3, the thermocline characteristic parameter inverse model of building are as follows:
In formula (14), T (z) and S (z) are respectively two elements of temperature and salinity of inverting, and T (z) and S (z) are depth of water z
Function, H are the number of samples of sound echo-signal, and i is the sequence number of sound echo signal sample, yiIndicate the sound echo that prediction obtains
Signal, for the convolution of incident acoustical signal and briny environment medium impulse response, yi' indicate the sound echo letter that actual measurement obtains
Number.
8. the thermocline layer characteristic parameter optimizing inversion method according to claim 7 based on acoustical signal, feature exist
In when solving to the thermocline characteristic parameter inverse model of building, utilizing population-genetic fusion algorithm pair in step 4
Thermocline characteristic parameter inverse model carries out optimizing iterative solution, obtains temperature profile parameter and salt profile parameter, further according to
Temperature profile parameter and salt profile parameter calculate Sound speed profile parameter are as follows:
C=1449.14+ Δ cT+ΔcS+ΔcP+ΔcSTP (17)
In formula:
ΔcT=4.5721T-4.4532 × 10-2T2-2.6045×10-4T3+7.985×10-6T4
ΔcS=1.3980 (S-35)+1.692 × 10-3(S-35)2
ΔcP=1.60272 × 10-1P+1.0268×10-5P2+3.5216×10-9P3-3.3603×10-12P4
ΔcSTP=(S-35) (- 1.1244 × 10-2T+7.7711×10-7T2+7.7016×10-5P-1.2943×10-7P2+
3.1580×10-8PT+1.5790×10-9PT2)+P(-1.8607×10-4T+7.4812×10-6T2+4.5283×10-8T3)+P2
(-2.5294×10-7T+1.8563×10-9T2)+P3(-1.9646×10-10T)
P=1.033+1.028126 × 10-1Z+2.38×10-7Z2-6.8×10-17Z4
T is temperature, and unit is degree Celsius;S is salinity, unit ‰;P is pressure, unit dimension gram/cm;Z is depth,
Unit is rice.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910684576.7A CN110398744A (en) | 2019-07-26 | 2019-07-26 | Ocean thermocline characteristic parameter optimizing and inverting method based on acoustic signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910684576.7A CN110398744A (en) | 2019-07-26 | 2019-07-26 | Ocean thermocline characteristic parameter optimizing and inverting method based on acoustic signals |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110398744A true CN110398744A (en) | 2019-11-01 |
Family
ID=68325213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910684576.7A Pending CN110398744A (en) | 2019-07-26 | 2019-07-26 | Ocean thermocline characteristic parameter optimizing and inverting method based on acoustic signals |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110398744A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111144666A (en) * | 2020-01-02 | 2020-05-12 | 吉林大学 | Ocean thermocline prediction method based on deep space-time residual error network |
CN111259943A (en) * | 2020-01-10 | 2020-06-09 | 天津大学 | Thermocline prediction method based on machine learning |
CN112115406A (en) * | 2020-09-28 | 2020-12-22 | 自然资源部第二海洋研究所 | Ocean internal mesoscale vortex inversion method and system based on remote sensing sea surface data |
CN113051795A (en) * | 2021-03-15 | 2021-06-29 | 哈尔滨工程大学 | Three-dimensional temperature-salinity field analysis and prediction method for offshore platform guarantee |
CN113218493A (en) * | 2021-04-08 | 2021-08-06 | 中国人民解放军国防科技大学 | Sound velocity profile inversion method based on empirical orthogonal function method |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1078349A (en) * | 1996-09-03 | 1998-03-24 | Oki Electric Ind Co Ltd | Marine acoustic tomographic field information reading device |
CN104063563A (en) * | 2014-07-16 | 2014-09-24 | 国家海洋局第一海洋研究所 | Method for calculating ocean spring layer characteristic values through multi-line-segment least square fitting |
KR101621143B1 (en) * | 2014-06-19 | 2016-05-16 | 포항공과대학교 산학협력단 | Underwater docking system based on underwater agent and the method of docking using thereof |
CN105911551A (en) * | 2016-05-09 | 2016-08-31 | 浙江大学 | Acoustic velocity profile inversion method based on weighted-EnKF algorithm |
CN108429591A (en) * | 2017-11-13 | 2018-08-21 | 西北工业大学 | A kind of multicarrier underwater acoustic communication method suitable for deep-sea channel |
CN108630209A (en) * | 2018-04-24 | 2018-10-09 | 中国科学院深海科学与工程研究所 | A kind of marine organisms recognition methods of feature based fusion and depth confidence network |
CN109145486A (en) * | 2018-09-05 | 2019-01-04 | 南通晟霖格尔电子科技有限公司 | The method of multi-line section least square fitting calculating ocean spring layer characteristic value |
CN109725053A (en) * | 2019-01-22 | 2019-05-07 | 中国人民解放军国防科技大学 | Method for acquiring water sound field characteristic data based on high-resolution ocean re-analysis product |
-
2019
- 2019-07-26 CN CN201910684576.7A patent/CN110398744A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1078349A (en) * | 1996-09-03 | 1998-03-24 | Oki Electric Ind Co Ltd | Marine acoustic tomographic field information reading device |
KR101621143B1 (en) * | 2014-06-19 | 2016-05-16 | 포항공과대학교 산학협력단 | Underwater docking system based on underwater agent and the method of docking using thereof |
CN104063563A (en) * | 2014-07-16 | 2014-09-24 | 国家海洋局第一海洋研究所 | Method for calculating ocean spring layer characteristic values through multi-line-segment least square fitting |
CN105911551A (en) * | 2016-05-09 | 2016-08-31 | 浙江大学 | Acoustic velocity profile inversion method based on weighted-EnKF algorithm |
CN108429591A (en) * | 2017-11-13 | 2018-08-21 | 西北工业大学 | A kind of multicarrier underwater acoustic communication method suitable for deep-sea channel |
CN108630209A (en) * | 2018-04-24 | 2018-10-09 | 中国科学院深海科学与工程研究所 | A kind of marine organisms recognition methods of feature based fusion and depth confidence network |
CN109145486A (en) * | 2018-09-05 | 2019-01-04 | 南通晟霖格尔电子科技有限公司 | The method of multi-line section least square fitting calculating ocean spring layer characteristic value |
CN109725053A (en) * | 2019-01-22 | 2019-05-07 | 中国人民解放军国防科技大学 | Method for acquiring water sound field characteristic data based on high-resolution ocean re-analysis product |
Non-Patent Citations (5)
Title |
---|
S. ZHENG AND X. YANGUANG: "Estimation of subsurface vertical thermal structure from sea surface temperature", 《2012 IEEE SYMPOSIUM ON ROBOTICS AND APPLICATIONS (ISRA), KUALA LUMPUR》 * |
侍茂崇: "《海洋调查方法导论》", 30 April 2008, 中国海洋大学出版社 * |
毛科峰,陈希,李妍,李振锋,梁毅: "基于声信号海温垂直结构反演试验", 《解放军理工大学学报(自然科学版)》 * |
毛科峰,陈希,李振锋,丁亚梅: "海水温盐垂直结构的声反射回波模型及应用", 《海洋学报(中文版)》 * |
邢传玺: "《浅海环境参数反演及声信号处理技术》", 30 June 2018, 北京理工大学出版社 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111144666A (en) * | 2020-01-02 | 2020-05-12 | 吉林大学 | Ocean thermocline prediction method based on deep space-time residual error network |
CN111259943A (en) * | 2020-01-10 | 2020-06-09 | 天津大学 | Thermocline prediction method based on machine learning |
CN112115406A (en) * | 2020-09-28 | 2020-12-22 | 自然资源部第二海洋研究所 | Ocean internal mesoscale vortex inversion method and system based on remote sensing sea surface data |
CN112115406B (en) * | 2020-09-28 | 2024-01-12 | 自然资源部第二海洋研究所 | Ocean internal mesoscale vortex inversion method and system based on remote sensing sea surface data |
CN113051795A (en) * | 2021-03-15 | 2021-06-29 | 哈尔滨工程大学 | Three-dimensional temperature-salinity field analysis and prediction method for offshore platform guarantee |
CN113051795B (en) * | 2021-03-15 | 2023-04-28 | 哈尔滨工程大学 | Three-dimensional Wen Yanchang analysis and prediction method for offshore platform guarantee |
CN113218493A (en) * | 2021-04-08 | 2021-08-06 | 中国人民解放军国防科技大学 | Sound velocity profile inversion method based on empirical orthogonal function method |
CN113218493B (en) * | 2021-04-08 | 2023-12-08 | 中国人民解放军国防科技大学 | Acoustic velocity profile inversion method based on empirical orthogonal function method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110398744A (en) | Ocean thermocline characteristic parameter optimizing and inverting method based on acoustic signals | |
CN109993280B (en) | Underwater sound source positioning method based on deep learning | |
Buchanan et al. | Marine acoustics: direct and inverse problems | |
CN111639747B (en) | GNSS-R sea surface wind speed inversion method and system based on BP neural network | |
CN113568055B (en) | Aviation transient electromagnetic data inversion method based on LSTM network | |
CN106932493B (en) | A kind of analysis method of the ultrasonic sensor array parameter based on genetic algorithm | |
CN109443516A (en) | A kind of passive acquisition methods of Bottom sound speed based on the vertical vibration velocity signal of noise field | |
CN113568056B (en) | Aviation transient electromagnetic data inversion method based on convolutional neural network | |
CN117198330B (en) | Sound source identification method and system and electronic equipment | |
CN109521469B (en) | Regularization inversion method for elastic parameters of submarine sediments | |
CN115841076A (en) | Shallow sea layered seabed ground sound parameter inversion method based on BP neural network model | |
Lermusiaux et al. | Multi-resolution probabilistic ocean physics-acoustics modeling: Validation in the New Jersey continental shelf | |
CN114970378A (en) | Sea clutter sample library construction method based on GAN network | |
CN117010132B (en) | Space array position optimization method and system of underwater multi-base sound system | |
CN116859460B (en) | Submarine physical property parameter inversion method suitable for elastic wave near-field reflection | |
CN116660996B (en) | Drifting type shallow sea local earth sound parameter prediction method based on deep learning | |
CN116956737A (en) | Method for inverting electron density of low ionosphere based on neural network and genetic algorithm | |
CN114974298A (en) | Ocean sound source positioning and ocean parameter inversion method based on deep learning | |
TANG et al. | Application of LVQ neural network combined with the genetic algorithm in acoustic seafloor classification | |
Kalikhman | Patchy distribution fields: sampling distance unit and reconstruction adequacy | |
Dasgupta et al. | Time-reversal imaging for classification of submerged elastic targets via Gibbs sampling and the relevance vector machine | |
Simons et al. | The potential of inverting geo-technical and geo-acoustic sediment parameters from single-beam echo sounder returns | |
Rixen et al. | Along or across front survey strategy? An operational example at an unstable front | |
CN109342568B (en) | Full-period voiceprint recognition method for instability and collapse of muddy submarine channel | |
Wichern et al. | Environmentally adaptive acoustic transmission loss prediction in turbulent and nonturbulent atmospheres |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191101 |
|
RJ01 | Rejection of invention patent application after publication |