Next Article in Journal
Specific Loss Power of Co/Li/Zn-Mixed Ferrite Powders for Magnetic Hyperthermia
Previous Article in Journal
Monitoring of Interfacial Debonding of Concrete Filled Pultrusion-GFRP Tubular Column Based on Piezoelectric Smart Aggregate and Wavelet Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Inversion for Geoacoustic Parameters in Shallow Sea †

1
College of Underwater Acoustic Engineering, Harbin Engineering University, Harbin 150001, China
2
Institute of Marine Science and Technology, Zhejiang Ocean University, Zhoushan 316022, China
3
State Key Laboratory of Acoustics, Chinese Academy of Sciences, Beijing 100190, China
4
Acoustic Science and Technology Laboratory, Harbin Engineering University, Harbin 150001, China
5
Key Laboratory of Marine Information Acquisition and Security Industry and Information Technology, Harbin Engineering University, Harbin 150001, China
6
Institute of Naval Architecture and Mechanical-electrical Engineering, Zhejiang Ocean University, Zhoushan 316022, China
*
Author to whom correspondence should be addressed.
This manuscript is extension version of the conference paper: Zheng, G.; Zhu, H.; Zhu, J. A method of geo-acoustic parameter inversion in shallow sea by the Bayesian theory and the acoustic pressure field. In the proceedings of the 2nd Franco-Chinese Acoustic Conference (FCAC2018), Le Mans, French, 29–31 October 2018.
These Authors contributed equally to this work.
Sensors 2020, 20(7), 2150; https://doi.org/10.3390/s20072150
Submission received: 11 March 2020 / Revised: 4 April 2020 / Accepted: 7 April 2020 / Published: 10 April 2020
(This article belongs to the Section Physical Sensors)

Abstract

:
Geoacoustic parameter inversion is a crucial issue in underwater acoustic research for shallow sea environments and has increasingly become popular in the recent past. This paper investigates the geoacoustic parameters in a shallow sea environment using a single-receiver geoacoustic inversion method based on Bayesian theory. In this context, the seabed is regarded as an elastic medium, the acoustic pressure at different positions under low-frequency is chosen as the study object, and the theoretical prediction value of the acoustic pressure is described by the Fast Field Method (FFM). The cost function between the measured and modeled acoustic fields is established under the assumption of Gaussian data errors using Bayesian methodology. The Bayesian inversion method enables the inference of the seabed geoacoustic parameters from the experimental data, including the optimal estimates of these parameters, such as density, sound speed and sound speed attenuation, and quantitative uncertainty estimates. The optimization is carried out by simulated annealing (SA), and the Posterior Probability Density (PPD) is given as the inversion result based on the Gibbs Sampler (GS) algorithm. Inversion results of the experimental data are in good agreement with both measured values and estimates from Genetic Algorithm (GA) inversion result in the same environment. Furthermore, the results also indicate that the sound speed and density in the seabed have fewer uncertainties and are more sensitive to acoustic pressure than the sound speed attenuation. The sea noise could increase the variance of PPD, which has less influence on the sensitive parameters. The mean value of PPD could still reflect the true values of geoacoustic parameters in simulation.

1. Introduction

In shallow sea environments, sound propagation is usually constrained and significantly affected by the geoacoustic parameters such as density, sound speed, and sound speed attenuation in the sea bottom [1]. The recent importance of study on acoustic problems in shallow sea has magnified the requirement for correct geoacoustic parameter estimates [2]. To obtain the sea bottom geoacoustic parameters, there are two methods. One uses the kinds of samplers to get the sea bottom properties directly, the other uses the acoustic signals to invert for these properties based on the sea bottom information implicitly carried in the signals, known as geoacousitc inversion. Furthermore, compared to the method of direct measurement, using acoustic signals to invert geoacoustic parameters are more efficient, and have attracted wide research interest over the past decade [3,4,5,6,7,8].
Nowadays, there are two ways to obtain marine environmental parameters. One is direct measurement, including on-site measurement and post-sampling analysis. The other method is acoustic inversion. Figure 1a,b show the on-site measurement and geoacoustic parameters inversion respectively. Compared with direct measurement, the acoustic inversion method can obtain the acoustic parameters quickly and efficiently, and avoid the waste of manpower and material resources.
Past studies revealed that geoacoustic inversion is mainly based on data measured by Vertical Linear Array (VLA) or Horizontal Linear Array (HLA). The sea bottom was usually regarded as a liquid medium, which only considers three geoacoustic parameters, namely, the compression wave (P-wave) speed, the density of the sea bottom, and P-wave speed attenuation. Most of these studies used the optimize methods such as Genetic Algorithm (GA) and Simulated Annealing (SA), etc., to solve the cost functions and thus get the optimal solution. But in practical application, a VLA or HLA carries a high computational cost, and the accuracy of the measured data can be affected greatly by the array shape, which is influenced by ocean currents and storms [9]. Therefore, more and more geoacoustic inversion based on low-frequency sound propagation characteristics have been proposed, and these methods greatly reduce the cost and complexity for recording systems in measurement [10,11,12,13]. However, the sea bottom was still set like a liquid medium. It has been proved that, when the low-frequency acoustic signal propagating in a shallow sea, the influence of shear wave (S-wave) in the bottom cannot be ignored [14], the sea bottom should be set as an elastic medium.
In addition, the geoacousitc inversion is typically a highly nonlinear, multi-parameter, and multi-optima inversion problem. There is usually not a single optimal solution due to the measurement and calculation errors. The inversion results of traditional optimize methods are all prone to being trapped by local minima, can only determine the best-fit model and not provide quantitatively nonlinear uncertainty estimation of the awaiting inversion parameters. In recent years, the Bayesian inversion method underwent considerable development and significantly contributed to estimating bottom properties and their uncertainties based on a Bayesian formulation [15,16,17,18,19]. The Bayesian inversion method is a global optimization algorithm based on the probability theory, applied mathematics, and optimization theory. It can efficaciously estimate Maximum Posterior Probability (MAP) model parameters, and then qualitatively and quantitatively analyze the uncertainty of parameter inversion results from a statistical perspective.
To close the gap, an effort has been made in this study to develop an inversion method for geoacoustic parameters by using a single acoustic sensor based on Bayesian theory. The significant impact from the elastic sea bottom on the low-frequency sound propagation in the shallow sea can be used to inference the geoacoustic parameters. According to Gibbs Sampler (GS) algorithm, the Posterior Probability Densities (PPDs) of five geoacoustic parameters, including the P-wave speed, S-wave speed, the density of the sea bottom, P-wave speed attenuation, and S-wave speed attenuation are given as the inversion results. Further, the effectiveness of this method is verified by the simulation analysis and an anechoic tank experiment. Before verification, the sound propagation characteristics were analyzed, which corresponded to the inversion result. For further confirmation, the Transmission Loss (TL) and the time domain diagram were compared, respectively, which include the measured results and the simulation results. The comparison shows that measured results are consistent with simulation results, which proves the reliability of the inversion method.
This paper is organized as follows: the prediction model of the underwater acoustic-field is described in Section 2 based on Pekeris waveguide. In Section 3.1 the inversion method based on Bayesian theory is introduced. The cost function is established by measuring sound pressure and the prediction sound pressure and presented in Section 3.2. The analysis of sound propagation characteristics is given in Section 4.1 followed by the inversion parameter’s sensitivity and inversion results of simulation data (simulated in the noiseless and noisy sea environment) in Section 4.2 and Section 4.3 respectively. The scaling experiment carried out in the laboratory is introduced in Section 5.1 followed by the procedure of the scaling experiment in Section 5.2. The inversion results of the scaling experiment and verification of the results are presented in Section 5.3 and Section 5.4, respectively. Finally, conclusions are given in Section 6.

2. Prediction Modeling of a Shallow Sea Acoustic Field

The S-wave of the seabed has a significant impact on the propagation of low-frequency acoustic signals in the shallow sea environments. The environment model in this study is developed based on Pekeris waveguide [20], and the seabed is regarded as an N-layered elastic medium. A schematic diagram of the waveguide model is shown in Figure 2. The z-axis indicates the depth of the sea, and the r-axis represents the propagation direction of acoustic signals. As shown in the figure, φ1, φpn, and ψsn are the displacement potential function in the water column, P-wave displacement potential function, and S-wave displacement potential function in the layered seabed, respectively. ρ1 and c1 are the density and sound speed in the water column. ρbn, cpn, and csn are the density, compression wave (P-wave) speed, and shear wave (S-wave) speed, respectively, of the n layer of seabed, αpn and αsn represent the attenuation of P-wave and S-wave, respectively; these five parameters in each layer of seabed are the expected inversion objects in this study. Moreover, zs is the depth of the point acoustic source, z = 0 and z = H are set as the sea surface and seabed in the model.
According to the wave theory, in the frequency domain, the displacement potential functions of each layer in Figure 2 satisfy the following equations:
1 r r ( r φ 1 r ) + 2 φ 1 z 2 + k 1 2 φ 1 = 4 π δ ( r , z z s ) ,   0 z < H 1
{ 1 r r ( r φ p n r ) + 2 φ p n z 2 + k p n 2 φ p n = 0 × × ψ s n k s n 2 ψ s n = 0 ,   H n z < H n + 1   ,   n 1
{ 1 r r ( r φ p N r ) + 2 φ p N z 2 + k p N 2 φ p N = 0 × × ψ s N k s N 2 ψ s N = 0 ,   z H N
where δ(r, z) is the source function, kmn=ω/cm (m=1,p,s and n=1…N) is the wave number of each layer and ω=2πf0 is the angular frequency of the point acoustic source at f0.
The general solutions of Equations (1)–(3) can be expressed as follows:
φ 1 ( r , z , ω ) = 0 Z 1 ( z , ξ , ω ) J 0 ( ξ r ) ξ d ξ ,   Z 1 ( z , ξ , ω ) = { A sin β z 0 z < z s B sin β 1 z + C cos β 1 z z s < z < H 1
{ φ p n ( r , z , ω ) = 0 Z p n ( z , ξ , ω ) J 0 ( ξ r ) ξ d ξ ψ s n ( r , z , ω ) = 0 Z s n ( z , ξ , ω ) J 1 ( ξ r ) ξ d ξ ,   { Z p n ( z , ξ , ω ) = P n u p e i β p n z + P n d o w n e i β p n z Z s n ( z , ξ , ω ) = S n u p e i β s n z + S n d o w n e i β s n z H n z < H n + 1   ,   1 n < N
{ φ p N ( r , z , ω ) = 0 Z p N ( z , ξ , ω ) J 0 ( ξ r ) ξ d ξ ψ s N ( r , z , ω ) = 0 Z s N ( z , ξ , ω ) J 1 ( ξ r ) ξ d ξ ,   { Z p N ( z , ξ , ω ) = P N e i β p N z Z s N ( z , ξ , ω ) = S N e i β sN z z H N
where β m n = k m n 2 ξ 2 (m = 1,p,s and n = 1…N), ξ is the horizontal wavenumber, J0 and J1 are the first order and second-order Bessel functions, respectively, and A, B, C, Pn and Sn are the indeterminate coefficients of the potential functions in each layer, the high mark ‘up’ and ‘down’ are used to representative the upgoing waves and downgoing waves in these potential functions. In the water column, the relationship between sound pressure p and potential function φ1 is p = ρ 1 ω 2 φ 1 .
Based on the point source condition, the boundary conditions on the fluid/solid interface for H1 (continuous normal displacement, continuous normal stress, and zero tangential stress) and Hn+1(continuous normal displacement, continuous normal stress, continuous tangential displacement, and continuous tangential stress), the relationship between these coefficients can be written as a (4N+1)th order matrix equations, which is shown in Equation (7).
( a i j ) ( 4 N + 1 ) × ( 4 N + 1 ) ( b i j ) ( 4 N + 1 ) × 1 = ( c i j ) ( 4 N + 1 ) × 1   N 1
The various indeterminate coefficients, such as A, B, C, Pn and Sn, can be solved by (bij) (4N+1) ×1 = [(aij)(4N+1) ×(4N+1)]−1. (cij) (4N+1) ×1. Substituting bij into Equations (4)–(6) further result in the potential functions in each layer. When replacing the coefficients in Equations (4)–(6), the expressions of displacement potential function in Figure 2 can be obtained. The expressions of sound pressure p in the water column, as the study object in inversion, can also be obtained with Equation (8).
p ( r , z , ω ) = ρ 1 ω 2 0 Z 1 ( z , ξ , ω ) J 0 ( ξ r ) ξ d ξ
p ( r j , z , ω ) = S ( ω ) ρ 1 ω 2 Δ ξ 2 π r e i ( ξ min r j π 4 ) l = 0 N s 1 [ Z 1 ( z , ξ l ) e i r min l Δ ξ ξ l ] e i 2 π l j N s
ξl = ξmin + lΔξ, l = 0, 1, 2, …, (NS−1), rj = rmin + jΔr, j = 0, 1, 2, …, (NS−1), Δr · Δξ = 2π/NS.
For the integral term of sound pressure in Equation (8), the Fast Field Method (FFM) considers the impact of lateral waves in the calculation and can measure the sound pressure easily and reliably [21]. It is chosen to describe the p for inversion in this study. It first discretizes the horizontal wavenumber ξ and the propagation distance r of Equation (9), and then the integral term is solved directly by using Fast Fourier Transform (FFT) [22].

3. Inversion Method

3.1. Bayesian Inversion Theory

Set the experiment measured data vector as d with elements di, and let m represent the model vector composed of the awaiting inversion geoacoustic parameters mi. Both di and mi are considered random variables that are related via Bayes’ rule [23]
P ( m | d ) = P ( d | m ) P ( m ) / P ( d )
The P(m|d) is the posterior probability density (PPD). P(d|m) is the conditional probability density function (PDF) of m under given d, P(m) is the prior PDF of m and representing the available parameter information independent of the data d, P(d) is the PDF of d. As the P(d) is independent of m, and the P(d|m) can be regarded as the likelihood function L(d|m) for the measured data [23], the Equation (10) can be written as:
P ( m | d ) L ( m ) P ( m )
The L(m) is determined by the form of the data and the statistical distribution of the data errors, including both measurement and model errors. Considering that it is difficult to obtain an independent estimate of the error statistics in practice, the assumption of unbiased Gaussian errors is used in processing, the form of the likelihood function is
L ( m ) = P ( d | m ) exp [ E ( m ) ]
where E(m) is the cost function. After normalizing
P ( m | d ) = exp [ E ( m ) ] P ( m ) exp [ E ( m ) ] P ( m ) d m
The integration domain spans an M-dimensional parameter space, M is the number of awaiting inversion parameters.
In Bayesian theory, the posterior probability density (PPD) of m can be regarded as the inversion results. To interpret the M-dimensional PPD requires us to estimate the properties of the parameter value, uncertainties, and inter-relationships, such as the MAP model m ^ , mean model m ¯ and marginal probability distributions P(mi|d). These properties are respectively defined as:
m ^ = Arg max { P ( m | d ) }
m ¯ = m P ( m | d ) d m
P ( m i | d ) = δ ( m i m i ) P ( m | d ) d m

3.2. Cost Function

A sufficient cost function is necessary to evaluate the geoacoustic parameters using the Bayesian inversion method. Via Bayes’ rule, the cost function is developed based on the likelihood function L(m)under the assumption of the Gaussian data errors [24].
L ( m ) = f = 1 F 1 π K | C m f | exp { [ p mea f p pre f ( m ) ] T ( C m f ) 1 [ p mea f p pre f ( m ) ] }
where p mea f ( m ) presents the measured sound pressure data at Kreceive positions by a single sensor under the fth frequency, p pre f ( m ) and C m f respectively present the model prediction sound pressure data and the covariance matrix under the same condition.
Consider the model prediction sound pressure p pre f ( m ) can be expressed as:
p pre f ( m ) = A f e i θ f p FFM f ( m )
where p FFM f ( m ) is the sound pressure computed via the FFM, Af and θf is the magnitude and phase of the unknown complex source at each frequency. To remove the dependence on Af and θf by setting the L ( m ) / A f = L ( m ) / θ f = 0 [25], the maximum likelihood solution of the source is
A f e i θ f = [ p FFM f ( m ) ] p mea f | p FFM f ( m ) | 2
where * denotes the conjugate transpose.
Ignoring the spatial correlation of the data, the common approximation of diagonal covariance C m f = v f I , where vf is the unknown variance at the fth frequency and I is the identity matrix, the L(m)can be simplified written as:
L ( m ) = f = 1 F 1 ( π ν f ) K exp [ B f ( m ) | p mea f | 2 ν f ]
where Equation (20) represents the normalized Bartlett disqualification:
B f ( m ) = 1 | [ p FFM f ( m ) ] p mea f | 2 | p mea f | 2 | p FFM f ( m ) | 2
To obtain a maximum-likelihood estimate of the data variance by setting the L ( m ) / v f = 0 , the maximum likelihood solution of the variance estimate is:
ν ^ f = B f ( m ) | p mea f | 2 K
Substituting Equation (22) into Equation (20) and Equation (12), the cost function E(m) becomes
E ( m ) = K f = 1 F In [ B f ( m ) | p mea f | 2 ]
where the Maximum Posterior Probability model m ^ for measured data with unknown variances can be found by minimizing the E(m), using a global optimization scheme.
As can be seen from Figure 3, in the research, the simulated annealing (SA) is used to optimize the m ^ of expecting geoacoustic parameters from Equation (23). The Gibbs sampler (GS) algorithm is applied to compute the PPD moments as defined in Equations (10)–(13), and the sample model is accumulated by SA at temperature T = 1 [26]. In Section 5, the genetic algorithm (GA) is also used to optimize the m ^ , which is regarded as a compression to verify the inversion method proposed in this paper. For SA settings, the starting T0 = 100, cooling rate ξ = 0.99, the number of samples n = 10000, and the type of perturbation is non-uniform. For GA settings, the population of each generation is pop = 1000, the type of coding is binary coding, the maximum generation is 100, the type of selection is roulette wheel selection, the cross-over rate is 0.95, and the mutation rate is 0.05. GS settings corresponded to the SA settings. A block diagram of the research is shown in Figure 3.

4. Simulation Results

Considering the shallow sea environment, where the seabed is regarded as a semi-infinite elastic medium, the following research is conducted based on the model with a 1-layered elastic bottom. In this model, the inversion objects are five geoacoustic parameters, namely, the P-wave speed cp, S-wave speed cs, the attenuation coefficient for the two speeds αp and αs, and seabed density ρb. The true values of waveguide parameters in the simulation are presented in Table 1.

4.1. The Analysis of the Acoustic Propagation Characteristics

In a simulation, the transmission loss (TL) was chosen to reveal the change of acoustic propagation characteristics, and the TL is defined as Equation (24).
TL = 20 lg | p ( r j , z , ω ) p r e f | r = 1 m | ,   p r e f = e i k r r
where the p(rj, z, ω) was calculated by Equation (22).
Figure 4 reveals the comparison results of the TL under different geo-acoustic parameters in the simulation. Figure 4a to Figure 3e show the influence of cp, cs, ρb, αp, and αs on acoustic propagation separately. The changing values of each parameter were set to the simulation true value deviation ±10%; when one of the above parameters changed, the other parameters remained fixed. In Figure 4a–e, the solid blue line represents the calculation results under the simulated true values, where the dotted black line and the dashed red line indicate the calculation results under the changing values. Figure 4f revealed the comparison of TLs’ anomalies when the value of each parameter was changed, and the change bounds were set between −10% to +10% of the true value [27].
It can be deduced from Figure 4a–e that the change of five parameters has a different impact on the acoustic propagation characteristics. Keeping the other simulation conditions fixed, the Transmission Loss (TL) changes the most when the cp and cs are varied. When ρb changes, the variation in the Transmission Loss (TL) is relatively prominent. When the αp or αs is changed, the variation in the Transmission Loss (TL) is the least obvious. The anomaly values of the five parameters in Figure 4f reveal that, in the discus bounds, these parameters are in the descending order of degree of influence on TL are cp, cs, ρb, αp, and αs. In this situation, the influence order of degree of five parameters on acoustic propagation characteristics can be preliminarily summarized as: cp>cs>ρb>αp>αs.

4.2. The Analysis of the Inversion Parameters’ Sensitivity

Using the simulation environmental conditions, as mentioned above, Figure 5a–e respectively represent the numerical variations of the cost function E(m) with the change of a single parameter. In each search bound of the parameter, E(m) touches the minimum value only at the true simulated value of the parameter, which can avoid the impact of the local optimum solution on the optimization of the cost function in the subsequent algorithm. Nonetheless, with the change of those five parameters, the range of E(m) variation is different. It can be seen from Figure 5f, that within the search bounds of five parameters, the parameters are in the descending order of degree of influence on E(m) are cp, cs, ρb, αp and αs. Therefore, the influence order of the degree of five geoacoustic parameters on p(r, z) can be defined as: cp > cs > ρb > αp > αs. This further verifies our results discussed in Section 4.1.

4.3. The Inversion Results

In this section, based on the research basis in Section 4.1 and Section 4.2, the feasibility and reliability of the inversion method are discussed with simulation. The focus of discussion is the application effect of the inversion method in noisy shallow sea environment. It can be seen from the literature [28], that the energy spectrum level of the noise in shallow sea could be −40 dB (0dB@1Pa), the energy spectrum level in simulation is set as −40 dB.
Figure 6 and Figure 7 show the PPDs of the five geoacoustic parameters in the noiseless and noisy sea environment respectively. In each figure, the vertical axis represents the probability, and true value means the simulated value indicated by the red line. The mean value and variance of inversion results are represented by the green segment in which the mean value is shown at the middle of the segment. The length of the green segment reflects the variance of each parameter.
It can be inferred that, in simulation, the maximum value of the probability density of each parameter is near to its true value, and the sharpness of the probability density curve shows the sensitivity of the cost function to each parameter: cp, cs > ρb > αp and αs (Figure 6 and Figure 7). Further, the results are lining with the findings presented in Section 4.1 and Section 4.2. Comparing Figure 6 and Figure 7, we can find that for the same parameter the PPD in Figure 7 is wider.
The marginal probability distribution between two parameters is presented in Figure 8 and Figure 9, which shows the uncertainty of different parameters. The true value of each parameter in the simulation is indicated with a white dashed line. It can be seen clearly from the figures that all the true values are close to the highest probability of PPD.
For the parameters cp, cs and ρb, the inversion results are similar in the two different marine environments, such as Figure 6a–c, Figure 7a–c, Figure 8a–c,e,f and Figure 9a–c,e,f. It is summarized that noise has less effect on the inversion results of sensitive parameters.
To clearly figure out the effect of noise on inversion results, Table 2 gives various inversion results between above two conditions. The value of each parameter is within the mean ± variance in the table. We can easily see that in the noiseless environment, the mean values of the inversion results are closer to the true value, and the variances are smaller. The sea noise could increase the variance of inversion results, but the mean of inversion results still could reflect the true values in simulation.
Figure 10 shows the comparison of TLs in the different sea environment. As can be seen from Figure 10, the inversion results are both consistent with the simulation results. The results in Figure 10 and Table 2 proved that no matter whether it is applied to a noiseless or noisy sea environment, the developed inversion method is reliable.

5. Results of Measured Data

5.1. Introduction to the Scaling Experiment

Suppose in a shallow sea waveguide, the various acoustic parameters (sound speed, sound speed attenuation and density) remain unchanged, the geometric parameters for the waveguide are scaled down by N times whereas the sound frequency is scaled up by N times, as H′ = H/N, zs = z/N, z′ = z/N, r′ = r/N, and f′ = Nf. In this case, the relationship between the original pressure p and the scaled pressure p′ can be obtained: p′(r′, z′, ω′) = Np(r, z, ω), and both the fluctuation and distribution of the original acoustic field remain unchanged in the scaled acoustic field [29].
Figure 11 displays the transmission losses under four scaling conditions, namely N = 0.1, 1, 10, 100. The corresponding sound frequency and water depth are 150 × 0.1 Hz, 150 × 1 Hz, 150 × 10 Hz, 150 × 100 Hz, and 100/0.1 m, 100/1 m, 100/10 m and 100/100 m, respectively. The acoustic source depth zs and receive depth zr are 20/0.1 m, 20/1 m, 20/10 m, 20/100 m and 10/0.1 m, 10/1 m, 10/10 m, 10/100 m respectively. The values of various acoustic parameters are set to be the same as the true value in Section 4.
The results of the theoretical analysis and the simulation indicate that the propagation of a high-frequency underwater acoustic signal in small-sized acoustic fields such as in a tank can be used to simulate the propagation of low-frequency underwater acoustic signals in large-sized environments such as the ocean. Thus, we carried out a scaling experiment, keeping in mind that the experiments performed in a tank in the laboratory have the advantages of stable environmental conditions, simple equipment distribution, and low cost.

5.2. Procedure and Data Processing of the Scaling Experiments

The experimental data were eventually utilized to validate the above-mentioned inversion method. The automatic measurement system is shown in Figure 12a,b. The experiment used a Poly Vinyl Chloride (PVC) slab (size: 1.53 m × 1.1 m × 0.105 m) to simulate the semi-infinite elastic bottom. A high-frequency underwater sound wave is transmitted by a source at a fixed position, and the wave is received by a single acoustic sensor at different positions with equal lengths. To increase the accuracy of the measurement, the acoustic sensor was fixed on a precision micro-worktable. With the micro-worktable, the acoustic sensor can move 2 mm every time with an error of less than 20 um, and a computer was used to control the micro-worktable and acquire data. When the measurement finishes at one position, the micro-worktable automatically takes the acoustic sensor to the next position according to the desired direction.
Table 3 shows the measurement parameters used in the experiment. The depth of source zs = 87 mm, the depth of acoustic sensor zr = 10 mm, the depth of water layer H = 182 mm. The sound speed c1 = 1450.212 m·s−1(obtained by the empirical formula with the water temperature in tank, the temperature is measured by the temperature recorder of Star-Oddi company [30]).
During the measurement experiment process, the source was fixed and transmitted the CW signal with the center frequency 135 kHz. The acoustic sensor was fixed on the movable walkway with a sampling rate at fs = 20 MHz. The starting point was 60 m away from the source position. The sampling rate of each measurement point was repeated ten times to reduce the influences of random error. After completing the one-point measurement, the walkway drives the acoustic sensor to move away from the transmitter by 2 mm, and repeated the same for 719 points. Figure 13a shows the TL measured in the experiment. Figure 13b shows the time-domain diagram of the 50th to 150th receiving signal received during the experiment. The red lines in Figure 13b represent the arrival time of the direct signal, the surface reflection signal, and the bottom reflection signal at the reception point from the 50th to the 150th. The measurement parameters are shown in Table 3. It can be inferred from Figure 13b that the arrival time of each path signal achieved by the simulation is consistent with the actual arrival time, which verifies the reliability of the measurement experiment.

5.3. Analysis of Inversion Results

Table 4 shows the inversion values of the seabed and search bounds of five parameters. In the processing of experimental data, the SA Settings are identical to the ones in Section 4.3.
Figure 14 shows the PPDs of the five geoacoustic parameters. The vertical axis represents the probability, and GA means the inversion result of GA indicated by the red line. The mean value and variance of inversion results are indicated by the green segment in which the mean value is shown at the middle of the segment. The length of the green segment reflects the variance of each parameter. We can easily see from Figure 14 that different parameters have different sensitivities.
The PPDs of cp, cs and ρb are relatively narrow in their prior search bounds, which proves cp, cs, ρb are more sensitive to the cost function, and have fewer uncertainties. The PPDs of αp, αs are flat in their prior search bounds, which means that αp and αs are not sensitive to the cost function. The sharpness of the probability density curve shows the sensitivity of the cost function to every parameter: cp, cs > ρb > αp, and αs. The inversion results obtained via the measurement experiment data are consistent with the simulation results discussed in Section 4.
Figure 15 represents the 2D marginal of the PPDs, where the dashed white line indicates the GA inversion result of each parameter, and the values of the variance represent the magnitude of the uncertainty. As shown in Figure 15a,b,e, the result of GA is near or at the highest probability. While in the Figure 15c,d,f–j, we can see that there is a little deviation between those two kinds of inversion methods, but for each parameter, the result of GA is within the margin of error. The density of the plastic plates was known in the measurement experiment, which was 1.2 g·cm−3. It is very near to the inversion value, which proves the viability of the inversion results. Thus, it is concluded that: cp, cs, and ρb are more sensitive and have fewer uncertainties to acoustic pressure than αp and αs.

5.4. The Verification

To verify the effectiveness of the inversion method, two different approaches are used as shown in Figure 16a,b. Figure 16a shows the comparison of the TLs obtained through measurement and simulation. As can be seen from Figure 16a, the measurement results are consistent with the simulation results. Figure 16b represents the signal amplitude of the measurement and simulation results in the time domain. Here, the two kinds of results showed a very high correlation, which reaches 0.89 when the effect of the transmitter, after it is off, is not considered. The after off effects of the transmitter are shown in Figure 16b. Thus, the approaches used here prove that the inversion results are credible.

6. Conclusions

A method of geoacoustic parameter inversion is proposed based on Bayesian Theory. Considering the S-wave of the seabed has a significant impact on the propagation of low-frequency acoustic signals in the shallow sea environments. Here, the five geoacoustic parameters, namely P-wave speed (cp), S-wave speed (cs), seabed density (ρb), p-wave attenuation (αp), and s-wave attenuation (αs) are inverted. In this paper, we used simulated data and the experimental data to verify the feasibility of the inversion method. The Posterior Probability Densities (PPDs) of the five parameters are obtained using the cost function based on the Bayesian inversion theory. The PPDs and 2D marginal PPDs provide parameter estimates and uncertainties. It is summarized that cp, cs, and ρb in the seabed have fewer uncertainties and are more sensitive to acoustic pressure than αp and αs.
Before verification, the acoustic propagation characteristics are analyzed, which corresponded to the inversion result. Our analysis results showed that in shallow sea environments, the influence order of the degree of five parameters on acoustic propagation characteristics could be summarized as cp > cs > ρb > αp > αs.
In the simulation, we find that the sea noise could affect the inversion results. In the noisy environment, the PPDs of every parameter will be wider than in the noiseless environment, the mean of PPD still could reflect the true values of geoacoustic parameters in simulation, and noise has less influence on the inversion results of sensitive parameters.
For further confirmation, the Transmission Loss (TL) and the normalized signal amplitude were compared, respectively, which include the measured results and the simulation results. The comparison shows that measured results are consistent with simulation results. Additionally, the density of the Poly Vinyl Chloride (PVC) plates was obtained in the measurement experiment; it was 1.2 g/m3, which is very close to the inversion value (1.2451 ± 0.0758g/m3). All these prove the credibility of the inversion method.
These results obtained from the geoacoustic inversion method will provide in-depth knowledge ‘from the basic principle to the advanced application of underwater acoustic sensors in the shallow sea acoustic propagation’ and will meet the growing interest in the need for correct geoacoustic parameter estimates and the uncertainties of the parameters.

Author Contributions

H.Z. consulted and arranged the theories in a considerable amount of literature. H.Z., X.W. and N.L. designed and took part in the experiments. H.Z. and G.Z. proposed the inversion method described in this paper. H.Z., G.Z. and Y.X. processed the tank experiment data. H.Z., G.Z. and S.K. wrote the whole paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was found by the National Natural Science Foundation of China (Grant No.: 11704337), State Key Laboratory of Acoustics, Chinese Academy of Sciences (SKLA201901) and the Open Foundation from Marine Sciences in the First-Class Subjects of Zhejiang (20190001).

Acknowledgments

We acknowledge the National Key Laboratory, College of Underwater Acoustic Engineering, Harbin Engineering University, China for supporting in this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jensen, F.B.; Kuperman, W.A.; Porter, M.B. ComputationalOcean Acoustic, 2nd ed.; Am Inst. Physics: New York, NY, USA, 1994. [Google Scholar]
  2. Wan, L.; Badiey, M.; Knobles, D.P. Geoacoustic inversion using low frequency broadband acoustic measurements from L-shaped arrays in the Shallow sea 2006 Experiment. J. Acoust. Soc. Am. 2016, 140, 2358–2373. [Google Scholar] [CrossRef] [PubMed]
  3. Gao, F.; Pan, C.M.; Sun, L. Geoacoustic parameters inversion of bayes matched-field: A multi-annealing gibbs sampling algorithm. Acta Armamentarii 2017, 38, 1385–1394. [Google Scholar] [CrossRef]
  4. Bevans, D.A.; Buckingham, M.J.; Hursky, P. A geoacoustic inversion technique using the low-frequency sound from the main rotor of a Robinson R44 helicopter. J. Acoust. Soc. Am. 2016, 140, 3169. [Google Scholar] [CrossRef]
  5. Sheng, X.; Dong, C.; Guo, L.; Li, L. A bioinspired twin inverted multi scale matched filtering method for detecting an underwater moving target in a reverberant environment. Sensors 2019, 19, 5305. [Google Scholar] [CrossRef] [Green Version]
  6. Bonnel, J.; Lin, Y.T.; Eleftherakis, D.; Goff, J.A.; Dosso, S.; Chapman, R.; Miller, J.H.; Potty, G.R. Geoacoustic inversion on the New England Mud Patch using warping and dispersion curves of high-order modes. J. Acoust. Soc. Am. 2018, 143, EL405–EL411. [Google Scholar] [CrossRef] [Green Version]
  7. Li, C.Z.; Yang, Y.; Wang, R.; Yan, X.F. Acoustic parameters inversion and sediment properties in the Yellow River reservoir. Appl. Geophys. 2018, 15, 78–90. [Google Scholar] [CrossRef]
  8. Michalopoulou, Z.; Gerstoft, P. Multipath Broadband Localization, Bathymetry, and Sediment Inversion. IEEE J. Ocean. Eng. 2019, 1–11. [Google Scholar] [CrossRef]
  9. Li, X.; Piao, S.; Zhang, M.; Liu, Y. A passive source location method in a shallow sea waveguide with a single sensor based on bayesian theory. Sensors 2019, 19, 1452. [Google Scholar] [CrossRef] [Green Version]
  10. Wan, L.; Badiey, M.; Knobles, D.P. Study on parameter correlations in the modal dispersion based geoacoustic inversion. J. Acoust. Soc. Am. 2017, 141, 3487. [Google Scholar] [CrossRef]
  11. Li, Z.; Zhang, R. Hybrid geoacoustic inversion method and its application to different sediments. J. Acoust. Soc. Am. 2017, 142, 2558. [Google Scholar] [CrossRef]
  12. Zheng, Z.; Yang, T.C.; Pan, X. Geoacoustic inversion using an autonomous underwater vehicle in conjunction with distributed sensors. IEEE J. Ocean. Eng. 2018, 1–23. [Google Scholar] [CrossRef]
  13. Liu, Y.; Zhang, H.; Li, Z.; Wang, X.; Ma, J. Particle Filtering for Localization of Broadband Sound Source Using an Ocean-Bottom Seismometer Sensor. Sensors 2019, 19, 2236. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, Y.; Yang, S.; Zhang, H.; Wang, X. Compressional-shear wave coupling induced by velocity gradient in elastic medium. Acta Phys. Sin. 2018, 67. [Google Scholar] [CrossRef]
  15. Bonnel, J.; Dosso, S.E.; Eleftherakis, D.; Chapman, N.R. Trans-Dimensional Inversion of Modal Dispersion Data on the New England Mud Patch. IEEE J. Ocean. Eng. 2019, 1–15. [Google Scholar] [CrossRef] [Green Version]
  16. Guo, X.; Yang, K.; Ma, Y. Geoacoustic inversion for bottom parameters via Bayesian theory in deep ocean. Chin. Phys. Lett. 2017, 34, 68–72. [Google Scholar] [CrossRef]
  17. Li, Q.; Yang, F.; Zhang, K. Moving source parameter estimation in an uncertain environment. Acta Phys. Sin. 2016, 65. [Google Scholar] [CrossRef]
  18. Ravenna, M.; Lebedev, S. Bayesian inversion of surface-wave data for radial and azimuthal shear-wave anisotropy, with applications to central Mongolia and west-central Italy. Geophys. J. Int. 2017, 213, 278–300. [Google Scholar] [CrossRef] [Green Version]
  19. Carlson, G.; Johnson, K.; Chuang, R.; Rupp, J. Spatially Varying Stress State in the Central U.S. From Bayesian Inversion of Focal Mechanism and In Situ Maximum Horizontal Stress Orientation Data. J. Geophys. Res. 2018, 123, 3871–3890. [Google Scholar] [CrossRef]
  20. Pekeris, C.L. Theory of propagation of explosive sound in shallow sea. Geol. Soc. Am. Mem. 1948, 27, 1–116. [Google Scholar] [CrossRef]
  21. Zhu, H.-H. Geoacoustic Parameters Inversion Based on Waveguide Impedance in Acoustic Vector Field. Ph.D. Dissertation, Harbin Engineering University, Harbin, China, 2014. [Google Scholar]
  22. Dinapoli, F.R. Fast Field Program for Multi-Layered Media; Naval Underwater Systems Center Newport Ri: Middletown, Rhode Island, 1971; pp. 2–10. [Google Scholar]
  23. Dosso, S.E.; Dettmer, J. Bayesian matched-field geoacoustic inversion. Inverse Probl. 2011, 27, 055009. [Google Scholar] [CrossRef] [Green Version]
  24. Tollefsen, D.; Dosso, S.E.; Knobles, D.P. Ship-of-opportunity noise inversions for geoacoustic profiles of a layered mud-sand seabed. IEEE J. Ocean. Eng. 2019, 1–12. [Google Scholar] [CrossRef]
  25. Yang, K.; Xiao, P.; Duan, R.; Ma, Y. Bayesian inversion for geoacoustic parameters from ocean bottom reflection loss. J. Comput. Acoust. 2017, 25, 1750019–1–1750019–17. [Google Scholar] [CrossRef]
  26. Metropolis, N. Equation of state calculations by fast computing machines. J. Chem. Phys. 1953, 21, 1087–1091. [Google Scholar] [CrossRef] [Green Version]
  27. Zheng, G.; Zhu, H.; Zhu, J. A method of geo-acoustic parameter inversion in shallow sea by the Bayesian theory and the acoustic pressure field. In Proceedings of the 2nd Franco-Chinese Acoustic Conference (FCAC2018), Le Mans, French, 29–31 October 2018. [Google Scholar]
  28. Zhao, W.; Chang, S. Underwater Acoustic, 2nd ed.; Science Press: Beijing, China, 1981. [Google Scholar]
  29. Zhu, H.; Zhu, J.; Zheng, G. A separation method for normal modes in shallow sea under near field. Acta Acust. 2019, 44, 41–50. [Google Scholar] [CrossRef]
  30. Liu, B.S.; Lei, J.Y. Principles of Underwater Sound, 2nd ed.; Harbin Engineering University Publications: Harbin, China, 2006. [Google Scholar]
Figure 1. (a) On-site measurement, (b) Geoacoustic parameters inversion.
Figure 1. (a) On-site measurement, (b) Geoacoustic parameters inversion.
Sensors 20 02150 g001
Figure 2. A layered shallow sea waveguide model.
Figure 2. A layered shallow sea waveguide model.
Sensors 20 02150 g002
Figure 3. The diagram of the research.
Figure 3. The diagram of the research.
Sensors 20 02150 g003
Figure 4. The influence of the five geoacoustic parameters on acoustic propagation, (a)-(e) corresponds to the cp, cs, ρb, αp and αs respectively, (f) reveals the comparison of TLs’ anomalies when the value of each parameter was changed.
Figure 4. The influence of the five geoacoustic parameters on acoustic propagation, (a)-(e) corresponds to the cp, cs, ρb, αp and αs respectively, (f) reveals the comparison of TLs’ anomalies when the value of each parameter was changed.
Sensors 20 02150 g004
Figure 5. The sensitivity of the cost function (E(m)) for the five geoacoustic parameters. (a) to (e)corresponds to the cp, cs, ρb, αp, and αs, respectively. (f) corresponds to the comparison of five parameter’s influence on E(m).
Figure 5. The sensitivity of the cost function (E(m)) for the five geoacoustic parameters. (a) to (e)corresponds to the cp, cs, ρb, αp, and αs, respectively. (f) corresponds to the comparison of five parameter’s influence on E(m).
Sensors 20 02150 g005
Figure 6. Five geo-acoustic parameters’ Posterior Probability Density (PPD) in the noiseless sea environment. The red lines mean the true value of each parameter in simulation. The green segments represent the mean value of the inversion results and their variance.
Figure 6. Five geo-acoustic parameters’ Posterior Probability Density (PPD) in the noiseless sea environment. The red lines mean the true value of each parameter in simulation. The green segments represent the mean value of the inversion results and their variance.
Sensors 20 02150 g006aSensors 20 02150 g006b
Figure 7. Five geo-acoustic parameters’ PPD in the noisy sea environment. The red lines mean the true value of each parameter in simulation. The green segments represent the mean value of the inversion results and their variance.
Figure 7. Five geo-acoustic parameters’ PPD in the noisy sea environment. The red lines mean the true value of each parameter in simulation. The green segments represent the mean value of the inversion results and their variance.
Sensors 20 02150 g007
Figure 8. 2D marginal PPDs between different parameters in the noiseless sea environment. White dashed lines mark the true values of each parameter in the simulation.
Figure 8. 2D marginal PPDs between different parameters in the noiseless sea environment. White dashed lines mark the true values of each parameter in the simulation.
Sensors 20 02150 g008
Figure 9. 2D marginal PPDs between different parameters in the noisy sea environment. White dashed lines mark the true values of each parameter in the simulation.
Figure 9. 2D marginal PPDs between different parameters in the noisy sea environment. White dashed lines mark the true values of each parameter in the simulation.
Sensors 20 02150 g009
Figure 10. The verification of simulation results. The solid blue line means the Transmission Loss (TL) measured in the simulation and the dashed red line represents the TL simulated by inversion result. (a) the comparison of TL in the noiseless sea environment; (b) the comparison of TL in the noisy sea environment.
Figure 10. The verification of simulation results. The solid blue line means the Transmission Loss (TL) measured in the simulation and the dashed red line represents the TL simulated by inversion result. (a) the comparison of TL in the noiseless sea environment; (b) the comparison of TL in the noisy sea environment.
Sensors 20 02150 g010
Figure 11. Transmission losses under the four scaling conditions. The solid blue line represents the frequency 150 × 0.1 Hz, the dashed red line represents the frequency 150 × 1 Hz, the dotted black line represents the frequency 150 × 10 Hz, and the dashed dotted green line represents the frequency 150 × 100 Hz.
Figure 11. Transmission losses under the four scaling conditions. The solid blue line represents the frequency 150 × 0.1 Hz, the dashed red line represents the frequency 150 × 1 Hz, the dotted black line represents the frequency 150 × 10 Hz, and the dashed dotted green line represents the frequency 150 × 100 Hz.
Sensors 20 02150 g011
Figure 12. (a) Experimental measurement system, (b) The movable micro-worktable for measurement equipment.
Figure 12. (a) Experimental measurement system, (b) The movable micro-worktable for measurement equipment.
Sensors 20 02150 g012
Figure 13. (a) The TL measured in the experiment; (b) The arrival time of the signal. The red line represents the arrival time of the direct signal, the dashed red line means the arrival time of surface reflection signal, and the dotted red line means the arrival time of the bottom reflection signal.
Figure 13. (a) The TL measured in the experiment; (b) The arrival time of the signal. The red line represents the arrival time of the direct signal, the dashed red line means the arrival time of surface reflection signal, and the dotted red line means the arrival time of the bottom reflection signal.
Sensors 20 02150 g013
Figure 14. Five geo-acoustic parameters’ PPD. Red lines mean the inversion result of each parameter by GA. Green segments represent the mean value of the inversion results and their variance.
Figure 14. Five geo-acoustic parameters’ PPD. Red lines mean the inversion result of each parameter by GA. Green segments represent the mean value of the inversion results and their variance.
Sensors 20 02150 g014
Figure 15. 2D marginal PPDs between different parameters. White dashed lines mark the Genetic Algorithm (GA) results of each parameter.
Figure 15. 2D marginal PPDs between different parameters. White dashed lines mark the Genetic Algorithm (GA) results of each parameter.
Sensors 20 02150 g015
Figure 16. The verification approaches, (a) solid blue line means the TL measured in experiment and the dashed red line represents the TL simulated by inversion result; (b) solid blue line means the normalized signal amplitude in the time domain at the 470th reception point (one meter from the starting point), and the dashed red line represents the normalized signal amplitude in the time domain simulated by inversion result. The after off effects of the transmitter are highlighted by arrows.
Figure 16. The verification approaches, (a) solid blue line means the TL measured in experiment and the dashed red line represents the TL simulated by inversion result; (b) solid blue line means the normalized signal amplitude in the time domain at the 470th reception point (one meter from the starting point), and the dashed red line represents the normalized signal amplitude in the time domain simulated by inversion result. The after off effects of the transmitter are highlighted by arrows.
Sensors 20 02150 g016
Table 1. The simulation parameters of ocean environment.
Table 1. The simulation parameters of ocean environment.
ParametersSimulated Value
Depth H/m100
Sound speed c1/m·s−11500
Sea-water density ρ1/ g·cm−31.025
P-wave speed cp/m·s−12000
S-wave speed cs/ m·s−11000
Seabed density ρb/g·cm−31.5
P-wave attenuation αp/dB·λ−10.5
S-wave attenuation αs/ dB·λ−10.5
Table 2. Inversion parameter list. The true values, search bounds and inversion values are shown.
Table 2. Inversion parameter list. The true values, search bounds and inversion values are shown.
ParametersTrue ValuesSearch BoundsInversion Values (Noiseless)Inversion Values (Noisy)
cp/m·s−120001800, 22002000.9367 ± 14.50022001.7484 ± 25.9427
cs/m·s−11000900, 1100996.5609 ± 9.19821000.9983 ± 11.7670
ρb/g·cm−31.51.35, 1.651.5040 ± 0.01211.5063 ± 0.0152
αp/dB·λ−10.50.45, 0.550.4974 ± 0.00460.5023 ± 0.0053
αs/dB·λ−10.50.45, 0.550.4989 ± 0.00500.5015 ± 0.0.0067
Table 3. The measurement parameters of the anechoic tank experiment.
Table 3. The measurement parameters of the anechoic tank experiment.
zs /mmzr /mmH /mmc1/m·s−1
87841821450.212
Table 4. Inversion parameters: the search bounds and inversion values are shown.
Table 4. Inversion parameters: the search bounds and inversion values are shown.
ParametersSearch BoundsInversion Values
cp/m·s−12200–25002397.3563 ± 31.9997
cs/m·s−11100–13001187.9400 ± 6.8722
ρb/g·cm−31.0–1.81.2451 ± 0.0758
αp/dB·λ−10.1–1.10.6616 ± 0.2489
αs/dB·λ−10.1-1.10.8705 ± 0.1468

Share and Cite

MDPI and ACS Style

Zheng, G.; Zhu, H.; Wang, X.; Khan, S.; Li, N.; Xue, Y. Bayesian Inversion for Geoacoustic Parameters in Shallow Sea. Sensors 2020, 20, 2150. https://doi.org/10.3390/s20072150

AMA Style

Zheng G, Zhu H, Wang X, Khan S, Li N, Xue Y. Bayesian Inversion for Geoacoustic Parameters in Shallow Sea. Sensors. 2020; 20(7):2150. https://doi.org/10.3390/s20072150

Chicago/Turabian Style

Zheng, Guangxue, Hanhao Zhu, Xiaohan Wang, Sartaj Khan, Nansong Li, and Yangyang Xue. 2020. "Bayesian Inversion for Geoacoustic Parameters in Shallow Sea" Sensors 20, no. 7: 2150. https://doi.org/10.3390/s20072150

APA Style

Zheng, G., Zhu, H., Wang, X., Khan, S., Li, N., & Xue, Y. (2020). Bayesian Inversion for Geoacoustic Parameters in Shallow Sea. Sensors, 20(7), 2150. https://doi.org/10.3390/s20072150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop