Cross-Ice Acoustic Communication:Cascade Acoustic Channel Model and Experimental Results

Cross-ice acoustic information transmission is an effective means of communication in polar sea areas covered by ice.However,the channel is extremely complicated because of the combined influence of water,ice,and air.Based on the normalmode theory,this paper establishes a cascade acoustic channel(CAC)model for the transmission of underwater acoustic waves across ice layer.The model can calculate the displacement response of the ice layer’s upper surface by separating the upward waves from normal modes in the water and multiplying it by a transmission coefficient matrix.The relationship between the displacement response of the upper surface of ice layer and the acoustic frequency is calculated by the finite-element method,and the calculation result was consistent with that of the CAC model.To verify the applicability of the model,a cross-ice acoustic communication experiment was conducted in Songhua River in January 2019.Experimental results show the energy of the acoustic signals received by geophones is closely related to sound frequency and crossice acoustic communication is feasible.The result of present research is important for understanding crossice acoustic channel characteristics and developing future cross-ice acoustic communication in polar sea areas.

Three-dimensional parabolic equation model for seismo-acoustic propagation: Theoretical development and preliminary numerical implementation

A three-dimensional（3D） parabolic equation（PE） model for sound propagation in a seismo-acoustic waveguide is developed in Cartesian coordinates, with x, y, and z representing the marching direction, the longitudinal direction, and the depth direction, respectively. Two sets of 3D PEs for horizontally homogenous media are derived by rewriting the 3D elastic motion equations and simultaneously choosing proper dependent variables. The numerical scheme is for now restricted to the y-independent bathymetry. Accuracy of the numerical scheme is validated, and its azimuthal limitation is analyzed. In addition, effects of horizontal refraction in a wedge-shaped waveguide and another waveguide with a polyline bottom are illustrated. Great efforts should be made in future to provide this model with the ability to handle arbitrarily irregular fluid-elastic interfaces.

Perfectly Matched Layer for an Elastic Parabolic Equation Model in Ocean Acoustics

The perfectly matched layer(PML) is an effective technique for truncating unbounded domains with minimal spurious reflections. A fluid parabolic equation(PE) model applying PML technique was previously used to analyze the sound propagation problem in a range-dependent waveguide(Lu and Zhu, 2007). However, Lu and Zhu only considered a standard fluid PE to demonstrate the capability of the PML and did not take improved one-way models into consideration. They applied a [1/1] Padé approximant to the parabolic equation. The higher-order PEs are more accurate than standard ones when a very large angle propagation is considered. As for range-dependent problems, the techniques to handle the vertical interface between adjacent regions are mainly energy conserving and single-scattering. In this paper, the PML technique is generalized to the higher order elastic PE, as is to the higher order fluid PE. The correction of energy conserving is used in range-dependent waveguides. Simulation is made in both acoustic cases and seismo-acoustic cases. Range-independent and range-dependent waveguides are both adopted to test the accuracy and efficiency of this method. The numerical results illustrate that a PML is much more effective than an artificial absorbing layer(ABL) both in acoustic and seismo-acoustic sound propagation modeling.

Depth Classification of Underwater Targets Based on Complex Acoustic Intensity of Normal Modes

In order to solve the problem of depth classification of the underwater target in a very low frequency acoustic field, the active component of cross spectra of particle pressure and horizontal velocity （ACCSPPHV） is adopted to distinguish the surface vessel and the underwater target. According to the effective depth of a Pekeris waveguide, the placing depth forecasting equations of passive vertical double vector hydrophones are proposed. Numerical examples show that when the sum of depths of two hydro- phones is the effective depth, the sign distribution of ACCSPPHV has nothing to do with horizontal distance; in addition, the sum of the first critical surface and the second critical surface is equal to the effective depth. By setting the first critical surface less than the difference between the effective water depth and the actual water depth, that is, the second critical surface is greater than the actual depth, the three positive and negative regions of the whole ocean volume are equivalent to two positive and negative regions and therefore the depth classification of the underwater target is obtained. Besides, when the 20 m water depth is taken as the first critical surface in the simulation of underwater targets （40 Hz, 50 Hz, and 60 Hz respectively）, the effectiveness of the algorithm and the cor- reemess of relevant conclusions are verified, and the analysis of the corresponding forecasting performance is conducted.

The Sensitive Properties of Hydrate Reservoirs Based on Seismic Stereoscopic Detection Technology

Higher-precision determinations of hydrate reservoirs, hydrate saturation levels and storage estimations are important for guaranteeing the ability to continuously research, develop and utilize natural gas hydrate resources in China. With seismic stereoscopic detection technology, which fully combines the advantages of different seismic detection models, hydrate formation layers can be observed with multiangle, wide-azimuth, wide-band data with a high precision. This technique provides more reliable data for analyzing the distribution characteristics of gas hydrate reservoirs, establishing velocity models, and studying the hydrate-sensitive properties of petrophysical parameters;these data are of great significance for the exploration and development of natural gas hydrate resources. Based on a velocity model obtained from the analysis of horizontal streamer velocity data in the hydrate-bearing area of the Shenhu Sea, this paper uses three VCs(longitudinal spacing of 25 m) and four OBSs(transverse spacing of 200 m) to jointly detect seismic datasets consisting of wave points based on an inversion of traveltime imaging sections. Accordingly, by comparing the differences between the seismic phases in the original data and the forward-modeled seismic phases, multiangle coverage constraint corrections are applied to the initial velocity model, and the initial model is further optimized, thereby improving the imaging quality of the streamer data. Petrophysical elastic parameters are the physical parameters that are most directly and closely related to rock formations and reservoir physical properties. Based on the optimized velocity model, the rock elastic hydrate-sensitive parameters of the hydrate reservoirs in the study area are inverted, and the sensitivities of the petrophysical parameters to natural gas hydrates are investigated. According to an analysis of the inversion results obtained from these sensitive parameters, λρ, Vp and λμ are simultaneously controlled by the bulk modulus and shear modulus, while Vs and μρ are controlled only by the shear modulus, and the latter two parameters are less sensitive to hydrate-bearing layers. The bulk modulus is speculated to be more sensitive than the shear modulus to hydrates. In other words, estimating the specific gravity of the shear modulus among the combined parameters can affect the results from the combined elastic parameters regarding hydrate reservoirs.

Fast Prediction of Acoustic Radiation from a Hemi-capped Cylindrical Shell in Waveguide

In order to predict acoustic radiation from a structure in waveguide, a method based on wave superposition is proposed, in which the free-space Green＇s function is used to match the strength of equivalent sources. In addition, in order to neglect the effect of sound reflection from boundaries, necessary treatment is conducted, which makes the method more efficient. Moreover, this method is combined with the sound propagation algorithms to predict the sound radiated from a cylindrical shell in waveguide. Numerical simulations show the effect of how reflections can be neglected if the distance between the structure and the boundary exceeds the maximum linear dimension of the structure. It also shows that the reflection from the bottom of the waveguide can be approximated by plane wave conditionally. The proposed method is more robust and efficient in computation, which can be used to predict the acoustic radiation in waveguide.

Feature Extraction for Acoustic Scattering from a Buried Target

Elastic acoustic scattering is important for buried target detection and identification. For elastic spherical objects, studies have shown that a series of narrowband energetic arrivals follow the first specular one. However, in practice, the elastic echo is rather weak because of the acoustic absorption, propagation loss, and reverberation, which makes it difficult to extract elastic scattering features, especially for buried targets. To remove the interference and enhance the elastic scattering, the de-chirping method was adopted here to address the target scattering echo when a linear frequency modulation (LFM) signal is transmitted. The parameters of the incident signal were known. With the de-chirping operation, a target echo was transformed into a cluster of narrowband signals, and the elastic components could be extracted with a band-pass filter and then recovered by remodulation. The simulation results indicate the feasibility of the elastic scattering extraction and recovery. The experimental result demonstrates that the interference was removed and the elastic scattering was visibly enhanced after de-chirping, which facilitates the subsequent resonance feature extraction for target classification and recognition.

Underwater square-root cubature attitude estimator by use of quaternion-vector switching and geomagnetic field tensor

This paper presents a kind of attitude estimation algorithm based on quaternion-vector switching and square-root cubature Kalman filter for autonomous underwater vehicle(AUV).The filter formulation is based on geomagnetic field tensor measurement dependent on the attitude and a gyro-based model for attitude propagation. In this algorithm, switching between the quaternion and the three-component vector is done by a couple of the mathematical transformations. Quaternion is chosen as the state variable of attitude in the kinematics equation to time update, while the mean value and covariance of the quaternion are computed by the three-component vector to avoid the normalization constraint of quaternion. The square-root forms enjoy a continuous and improved numerical stability because all the resulting covariance matrices are guaranteed to stay positively semidefinite. The entire square-root cubature attitude estimation algorithm with quaternion-vector switching for the nonlinear equality constraint of quaternion is given. The numerical simulation of simultaneous swing motions in the three directions is performed to compare with the three kinds of filters and the results indicate that the proposed filter provides lower attitude estimation errors than the other two kinds of filters and a good convergence rate.

Simulation and verifi cation of an air-gun array wavelet in time-frequency domain based on van der waals gas equation

An air gun generates acoustic signals for seismic exploration by releasing a high-pressure gas.A large error is always gradually introduced into the ideal-gas model when the pressure in the air-gun chamber exceeds 100 atm.In the van der Waals non-ideal-gas theory,the gas in the air gun can be regarded as an actual gas,and the error is less than 2%.The van der Waals model is established in combination with the quasi-static open thermodynamic system and bubble-motion equation by considering the bubble rise,bubble interaction,and throttling eff ect.The mismatch between the van der Waals and ideal-gas models is related to the pressure.Theoretically,under high-pressure conditions,the van der Waals air-gun model yields results that are closer to the measured results.Marine vertical cables are extended to the seafl oor using steel cables that connect the cement blocks,but the corresponding hydrophones are suspended in the seawater.Thus,noise associated with ships,ocean surges,and coupling problems is avoided,and the signal-to-noise ratio and resolution of marine seismic data are improved.This acquisition method satisfies the conditions of recording air-gun far-fi eld wavelets.According to an actual vertical-cable observation system,the van der Waals air-gun model is used to model the wavelet of different azimuth and take-off angles.The characteristics of the experimental and simulated data demonstrate good agreement,which indicates that the van der Waals method is accurate and reliable.The accuracy of the model is directly related to the resolution,thus aff ecting the resolution ability of the stratum.

MIMO Underwater Acoustic Communication in Shallow Water with Ice Cover

Although multiple-input multiple-output(MIMO) underwater acoustic(UWA) communication has been intensively investigated in the past years, existing works mainly focus on open-water environment. There is no work reporting MIMO acoustic communication in under-ice environment. This paper presents results from a recent MIMO acoustic communication experiment which was conducted in Bohai Gulf during winter. In this experiment, high frequency MIMO signals centered at 10 kHz were transmitted from a two-element source array to a four-element vertical receiving array at 1 km range. According to the received signal of different array elements, MIMO acoustic communication in under-ice environment suffers less effect from co-channel interference compared with that in open-water environment. In this paper, time reversal followed by a single channel decision feedback equalizer is used to process the experimental data. It is demonstrated that this simple receiver is capable of realizing robust performance using fewer hydrophones(i.e. 2) without the explicit use of complex co-channel interference cancelation algorithms, such as parallel interference cancelation or serial interference cancelation.

A Block Parallel l_0-Norm Penalized Shrinkage and Widely Linear Affine Projection Algorithm for Adaptive Filter

To improve the identification capability of AP algorithm in time-varying sparse system, we propose a block parallel l_0-SWL-DCD-AP algorithm in this paper. In the proposed algorithm, we first introduce the l_0-norm constraint to promote its application for sparse system. Second, we use the shrinkage denoising method to improve its track ability. Third, we adopt the widely linear processing to take advantage of the non-circular properties of communication signals. Last, to reduce the high computational complexity and make it easy to implemented, we utilize the dichotomous coordinate descent(DCD) iterations and the parallel processing to deal with the tapweight update in the proposed algorithm. To verify the convergence condition of the proposed algorithm, we also analyze its steadystate behavior. Several simulation are done and results show that the proposed algorithm can achieve a faster convergence speed and a lower steady-state misalignment than similar APA-type algorithm. When apply the proposed algorithm in the decision feedback equalizer(DFE), the bite error rate(BER) decreases obviously.

Blind Source Separation based on Time-Frequency Morphological Characteristics for Rigid Acoustic Scattering by Underwater Objects

Separation of the components of rigid acoustic scattering by underwater objects is essential in obtaining the structural characteristics of such objects. To overcome the problem of rigid structures appearing to have the same spectral structure in the time domain, time-frequency Blind Source Separation （BSS） can be used in combination with image morphology to separate the rigid scattering components of different objects. Based on a highlight model, the separation of the rigid scattering structure of objects with time-frequency distribution is deduced. Using a morphological filter, different characteristics in a Wigner-Ville Distribution （WVD） observed for single auto term and cross terms can be simplified to remove any cross-term interference. By selecting time and frequency points of the auto terms signal, the accuracy of BSS can be improved. A simulation experimental has been used to analyze the feasibility of the new method, with changing the pulse width of the transmitted signal, the relative amplitude and the time delay parameter. And simulation results show that the new method can not only separate rigid scattering components, but can also separate the components when elastic scattering and rigid scattering exist at the same time. Experimental results confirm that the new method can be used in separating the rigid scattering structure of underwater objects.

Bubble acoustical scattering cross section under multi-frequency acoustic excitation

The acoustical scattering cross section is usually employed to evaluate the scattering ability of the bubbles when they are excited by the incident acoustic waves. This parameter is strongly related to many important applications of performance prediction for search sonar or underwater telemetry, acoustical oceanography, acoustic cavitation, volcanology, and medical and industrial ultrasound. In the present paper, both the analytical and numerical analysis results of the acoustical scattering cross section of a single bubble under multi-frequency excitation are obtained. The nonlinear characteristics（e.g.,harmonics, subharmonics, and ultraharmonics） of the scattering cross section curve under multi-frequency excitation are investigated compared with single-frequency excitation. The influence of several paramount parameters（e.g., bubble equilibrium radius, acoustic pressure amplitude, and acoustic frequencies） in the multi-frequency system on the predictions of scattering cross section is discussed. It is shown that the combination resonances become significant in the multi-frequency system when the acoustic power is big enough, and the acoustical scattering cross section is promoted significantly within a much broader range of bubble sizes and acoustic frequencies due to the generation of more resonances.

Characteristic Analysis of Underwater Acoustic Scattering Echoes in the Wavelet Transform Domain

Underwater acoustic scattering echoes have time–space structures and are aliasing in time and frequency domains. Different series of echoes properties are not identified when incident angle is unknown. This article investigates variations in target echoes of monostatic sonar to address this problem. The mother wavelet with similar structures has been proposed on the basis of preprocessing signal waveform using matched filter, and the theoretical expressions between delay factor and incident angle are derived in the wavelet domain. Analysis of simulation data and experimental results in free-field pool show that this method can effectively separate geometrical scattering components of target echoes. The time delay estimation obtained from geometrical echoes at a single angle is consistent with target geometrical features, which provides a basis for object recognition without angle information. The findings provide valuable insights for analyzing elastic scattering echoes in actual ocean environment.

An improved least mean square/fourth direct adaptive equalizer for under-water acoustic communications in the Arctic

An improved least mean square/fourth direct adaptive equalizer(LMS/F-DAE)is proposed in this paper for underwater acoustic communication in the Arctic.It is able to process complex-valued baseband signals and has better equalization performance than LMS.Considering the sparsity feature of equalizer tap coefficients,an adaptive norm(AN)is incorporated into the cost function which is utilized as a sparse regularization.The norm constraint changes adaptively according to the amplitude of each coefficient.For small-scale coefficients,the sparse constraint exists to accelerate the convergence speed.For large-scale coefficients,it disappears to ensure smaller equalization error.The performance of the proposed AN-LMS/F-DAE is verified by the experimental data from the 9th Chinese National Arctic Research Expedition.The results show that compared with the standard LMS/F-DAE,AN-LMS/F-DAE can promote the sparse level of the equalizer and achieve better performance.

Space-division multiple access for CDMA multiuser underwater acoustic communications

Time reversal mirror （TRM） can use the physical characteristics of the underwater acoustic （UWA） channel to focus on the desired user in multi-user UWA communication. The active average sound intensity （AASI） detector can estimate all azimuths of users with the same frequency band at the same time in order to achieve directional communication by vector combination. Space-division multiple access （SDMA） based on TRM combined with the AASI detector is proposed in this paper, which can make the capacity of the code division multiple access （CDMA） UWA system significantly increase. The simulation and lake test results show that the 7-user UWA multi-user system can achieve low bit error communication.

Bubble nonlinear dynamics and stimulated scattering process

A complete understanding of the bubble dynamics is deemed necessary in order to achieve their full potential applications in industry and medicine. For this purpose it is first needed to expand our knowledge of a single bubble behavior under different possible conditions including the frequency and pressure variations of the sound field. In addition, stimulated scattering of sound on a bubble is a special effect in sound field, and its characteristics are associated with bubble oscillation mode. A bubble in liquid can be considered as a representative example of nonlinear dynamical system theory with its resonance, and its dynamics characteristics can be described by the Keller–Miksis equation. The nonlinear dynamics of an acoustically excited gas bubble in water is investigated by using theoretical and numerical analysis methods. Our results show its strongly nonlinear behavior with respect to the pressure amplitude and excitation frequency as the control parameters, and give an intuitive insight into stimulated sound scattering on a bubble. It is seen that the stimulated sound scattering is different from common dynamical behaviors, such as bifurcation and chaos, which is the result of the nonlinear resonance of a bubble under the excitation of a high amplitude acoustic sound wave essentially. The numerical analysis results show that the threshold of stimulated sound scattering is smaller than those of bifurcation and chaos in the common condition.

Parallel optimization of underwater acoustic models:A survey

Underwater acoustic models are effective tools for simulating underwater sound propagation.More than 50 years of research have been conducted on the theory and computational models of sound propagation in the ocean.Unfortunately,underwater sound propagation models were unable to solve practical large-scale three-dimensional problems for many years due to limited computing power and hardware conditions.Since the mid-1980s,research on high performance computing for acoustic propagation models in the field of underwater acoustics has flourished with the emergence of high-performance computing platforms,enabling underwater acoustic propagation models to solve many practical application problems that could not be solved before.In this paper,the contributions of research on high-performance computing for underwater acoustic propagation models since the 1980s are thoroughly reviewed and the possible development directions for the future are outlined.

Improving sound diffusion in a reverberation tank using a randomly fluctuating surface

Underwater reverberation environments that satisfy the conditions of uniformity and isotropy of the diffuse field can be used to measure the acoustic characteristics of underwater targets.This study combines two practical indicators—the standard deviation of the absolute sound pressure field(to indicate uniformity)and the analysis of the wavenumber spectrum in the spherical harmonics domain(to indicate isotropy)—for an accurate evaluation of the diffusion of the sound field in a reverberation tank.A method is proposed that can improve the narrow-band diffusion of the sound field by employing a randomly fluctuating surface.An acoustic experiment was performed in a reverberation water tank(1.2 m×1 m×0.8 m),where a randomly fluctuating surface was generated by making waves.The experimental results show that as the wave motion contributes effectively to the random reflection of sound rays in all directions,the uniformity and isotropy are improved significantly when the surface is fluctuating randomly.This work helps to ensure accurate measurements of the characteristics of underwater targets in reverberation tanks.

Radiation from finite cylindrical shell with irregular-shaped acoustic enclosure

In practical situations, large machinery is usually placed in an underwater vessel and changes the acoustic enclosure shape into an irregular one. The existence of machinery causes the difficulties in expressing sound transmission and radiation analytically. In this study, the sound radiation of a cylindrical shell excited by an internal acoustic source is modeled and analyzed. The cylindrical shell contains a machine modeled as a rectangular object, which is attached to a shell with a spring-mass system. The acoustic field of the cavity is computed by the integro-modal approach. The effect of object size on the coupling between acoustic mode and structural mode is investigated. The relationship between object volume and sound radiation is also studied. Numerical results show that the existence of objects inside vessels leads to a more effective coupling between the structure and acoustic enclosure than the existence of no objects in a regular-shaped cavity（i.e. empty vessel）.