Reflection response for normal incidence probing of vertically nonhomogeneous media

A time domain synthetic reflection seismogram is detailed and, as a limiting condition on this model, the analytic reflection impulse response for a one-dimensional lossless acoustic medium with piecewise continuous acoustic impedance is obtained. This analytic impulse response solution, in the structure of a sum of terms by order of reflection, provides insight to some of the poorly understood aspects of acoustic reflections from stratified and smoothly varying media as may be encountered in shallow marine sediments and elsewhere. It offers as well an approach for the inverse problem of extracting acoustic impedance profiles from reflection response data, though other effects (such as wavefront spreading, dispersive and absorptive attenuation, and wavelet broadening attendant with pulse propagation through a medium) need to be accommodated.     

The effect of a layer of varying density on high-frequencyreflection, forward loss, and backscatter [seafloor acoustics]

This paper examines the impact that a thin layer of varying density would have on high-frequency reflection, forward loss, and backscattering of acoustic plane waves from the seafloor. A functional form for density stratification was found by examination of several high-resolution density profiles obtained from X-ray computed tomography (CT) scans of seafloor cores. A solution based on these general profiles was used to estimate the reflection coefficient. The influence of the density profile on reflection loss and backscatter was then calculated using the estimated reflection coefficient. Parameter values used in simulations were obtained from the literature or from the CT scans of cores. It was found that inclusion of a density profile adds a strong frequency dependence to estimates of the reflection coefficient and forward loss. The largest effect on total scattering strength is near normal incidence where returns are dominated by interface scattering. The effect of the density profile on the strength of acoustic returns suggests that care should be taken when using high-frequency systems for measuring sediment properties, especially near normal incidence     

Effects of seabed properties on acoustic wave fields in a seismo-acoustic ocean waveguide

Acoustic wave fields in an ocean waveguide with a sediment layer having continuously varying density and sound speed overlying an elastic subbottom are considered in this analysis. The objective of this study is to investigate the effects of seabed acoustic properties, including the density and sound speed of the sediment layer and subbottom, on the characteristics of the wave fields. Examination of the reflection coefficient, wavenumber spectrum, and noise intensity of the sound field through numerical analysis has shown that the variation in the acoustic properties in the sediment layer is an important factor in determining the reflected or noise sound fields. In particular, the sediment thickness-to-wavelength ratio and the types of variation of acoustic properties inside the layer give rise to many characteristics that potentially allow for acoustic inversion of the seabed properties. With regard to the wave-field components in a shallow-water environment, the various types of waves existing in a seismo-acoustic waveguide have been illustrated. The results indicate that the effects of the sediment properties on the wavenumber spectrum are primarily on the continuous and evanescent regimes of the wave field. The noise intensity generated by distributive random monopoles at various depths, together with the effect of refractive sound-speed distribution in the water column, has been obtained and analyzed.     

分层介质声特征阻抗反演 -预处理方法实验验证与比较

数值不稳定性长期困扰着声反演方法的应用.避免该问题的方法之一是进行冲激响应提取的预处理,由冲激响应反演实验材料的各层反射系数序列,重建特性阻抗剖面.该文采用不同反卷积的方法提取冲激响应,反演阻抗剖面,并就精度及运算时间进行详细比较,反演结果与实测及手册数值相当吻合.     

Acoustic density measurements of consolidating cohesive sediment beds by means of a non-intrusive “Micro-Chirp” acoustic system

A non-intrusive “Micro-Chirp” acoustic system and a signal-processing protocol have been developed to estimate the bulk density of consolidating cohesive sediment beds. Using high-frequency (300–700 kHz) Chirp acoustic waves, laboratory measurements were conducted with clay–water mixtures. Because acoustic echo strength is proportional to variations in acoustic impedance, and the speed of sound in the clay bed hardly changed during consolidation, the bulk density could be successfully estimated without disturbing the sediment bed. Based on acoustic signal analysis, this study demonstrates that the reflection coefficient and bulk density at the water–sediment interface increase with consolidation time, and that a single speed of sound value can be used for practical bulk density estimation in muddy environments.     

Ocean acoustic tomography as a data assimilation problem

The ocean acoustic tomographic (OAT) approach to sound speed field estimation is generalized to include a variety of sources of information of interest such as an oceanographic model of the sound speed field, direct local sound speed measurements, and a full field acoustic propagation model as well as measurements. The inverse problem is presented as a four-dimensional field estimation problem using a variational approach commonly used in oceanographic data assimilation. The current OAT approach is shown to be a special case of the general framework. The matched-field tomography (MFT) approach is also discussed within this context. A simple implementation of this novel approach is then investigated in the absence of a suitable oceanographic model, and acoustic propagation is accounted for using a standard parabolic equation model. The inverse equations derived are validated numerically through a simple inversion example, and some issues on environmental mismatch and computations are discussed. The developments then provide a basic framework for ongoing data-model melding in acoustically focused oceanographic sampling (AFOS) network     

Relationships between the sound speed ratio and physical properties of surface sediments in the South Yellow Sea

Building empirical equations is an effective way to link the acoustic and physical properties of sediments. These equations play an important role in the prediction of sediments sound speeds required in underwater acoustics.Although many empirical equations coupling acoustic and physical properties have been developed over the past few decades, further confirmation of their applicability by obtaining large amounts of data, especially for equations based on in situ acoustic measurement techniques, is required. A sediment acoustic survey in the South Yellow Sea from 2009 to 2010 revealed statistical relationships between the in situ sound speed and sediment physical properties. To improve the comparability of these relationships with existing empirical equations, the present study calculated the ratio of the in situ sediment sound speed to the bottom seawater sound speed, and established the relationships between the sound speed ratio and the mean grain size, density and porosity of the sediment. The sound speed of seawater at in situ measurement stations was calculated using a perennially averaged seawater sound speed map by an interpolation method. Moreover, empirical relations between the index of impedance and the sound speed and the physical properties were established. The results confirmed that the existing empirical equations between the in situ sound speed ratio and the density and porosity have general suitability for application. This study also considered that a multiple-parameter equation coupling the sound speed ratio to both the porosity and the mean grain size may be more useful for predicting the sound speed than an equation coupling the sound speed ratio to the mean grain size.     

A simplified pseudo-fluid model derived from biot theory through low grazing angle approximation

We propose an approximation technique with which the porous Biot model can be converted into a pseudo-fluid medium model, i.e., a medium represented by its sound speed and density. This technique begins from an analytic derivation of the reflection coefficient on a planar interface of fluid and porous ocean sediment. Invoking the low grazing angle approximation, useful for solving long range propagation problems, the pseudo-fluid medium is obtained. Alternate pseudo-fluid medium can be obtained through the weak frame approximation (Williams, J. Acoust. Soc. Am. 110, 2276-2281, 2001). In this paper, we discuss the accuracy and limitations of the low grazing angle approximation by numerically comparing the reflection coefficient to the full Biot model cases as functions of various Biot parameters, frequency, and water sound speed. The usefulness of the present low grazing angle approximated pseudo-fluid medium in ocean wave propagation modeling is demonstrated by comparing the transmission loss results with that of the full Biot model.     

Remote sensing of ocean sound speed profiles by a perceptron neuralnetwork

A method is developed to estimate ocean sound speed profiles through synthesis of remotely measured environmental data and historical statistics of sound speed obtained at a remotely sensed location. Sound speed profiles are represented by an expansion of empirical orthogonal functions (EOF) of the historical sound speed variation, while the remotely sensed environmental data provide real-time information to determine the expansion coefficients. Environmental inputs are limited to sea surface temperature available from satellite infrared sensors, acoustic time-of-flight and ocean bottom temperature measurable from bottom mounted acoustic and thermal transducers. A multilayer perceptron neural network is implemented to learn the functional transformation from the measured environmental input to the desired EOF coefficient output on a set of representative sound speed profiles. Sea surface temperature, time-of-year, and time-of-flight from the acoustic multipath that maximally samples the vertical sound speed are found to be the dominant inputs. The trained network is computationally efficient and produces estimates for untrained environmental inputs with a mean error of 1.1-4.4 m/s     

Laboratory experiment to measure 5-MHz volume backscattering strengths from Red-tide causing microalgae Chattonella antiqua

An acoustic laboratory experiment using 5-MHz signals was conducted to measure the volume backscattering strengths of red-tide causing microalgae, Chattonella antiqua, which is one of the species of harmful algal blooms in the coastal waters of Korea and Japan. The measured backscattering strengths increased with cell abundance, with a slope of approximately 10 dB per decade increase in cell numbers. The density and sound speed ratios of the Chattonella cell to the water medium were estimated via the density gradient centrifugation method and the time-travel difference method, respectively. Finally, the measured backscattering strengths were compared to those predicted by a fluid-sphere scattering model, in which the estimated sound speed and density contrasts were used as input parameters.     

Analysis of optimal conditions for recording signals when studying the atmospheric boundary layer with the acoustic method of partial reflection

Equations for the coefficient of partial reflection K from stratified inhomogeneities in the atmospheric boundary layer have been derived on the basis of the Epstein transition and symmetrical layer models as functions of three dimensionless parameters, i.e., the relative layer altitude, its relative thickness, and the relative variations in the effective sound speed in a layer. The equations have been obtained for the relative layer altitude at which the total internal reflection appears; the behavior of the function K is studied at close altitudes. Significant weakening of the dependence of coefficient K on the relative layer thickness in these conditions is shown, which makes it possible to record partially reflected signals in a wide range of wave-lengths or frequencies of the sounding signal. In other cases, the coefficient of partial reflection K strongly depends on the layer thickness. According to experimental data on variations in the amplitude of received acoustic signals with an increase in the source-detector distance, a technique for the parameterization of the additional impedance attenuation of sound that propagates over the earth’s surface has been developed, and these parameters have been experimentally estimated for different stratification conditions and sounding signal frequencies. Many records of background acoustic noises typical for one or another measurement sites have been distinguished and classified, a technique for estimating the minimum signal amplitude distinguishable against noises has been developed, and the corresponding estimates have been made. Based on these data and the specifications of three different industrial acoustic sources, the parameter limits provided by these sources have been estimated for the method of partial reflection.     

Oceanic sound speed profile inversion

A modal (full-wave) method has been developed to predict ocean sound speed profiles from propagated acoustic field data. The method assumes a point source of sound in the ocean and uses as data the values of the transmitted acoustic field at an array. The formalism for depth-dependent sound speeds consists of the standard Hankel integral transform of the depth solution. In the travel length coordinate, the latter is written exactly, using the Green's function, in terms of an integral equation whose kernel includes the sound speed profile correction. A Born approximation to this equation is used. This is just the WKB solution, and permits the use of a nontrivial input (or guess) profile, here chosen as bilinear. The use of asymptotic methods enables us to write the data as an integral transform over the profile correction. The transform can be inverted. An example is presented for full-bandwidth inversion.     

Optimal localization of a seafloor transponder in shallow water using acoustic ranging and GPS observations

This study utilized circular and straight-line survey patterns for acoustic ranging to determine the position of a seafloor transponder and mean sound speed of the water column. To reduce the considerable computational burden and eliminate the risk of arriving at a local minimum on least-squares inversion, the position of a seafloor transponder was estimated by utilizing optimization approaches. Based on the implicit function theorem, the Jacobian for this inverse problem was derived to investigate the constraints of employing circular and straight-line survey patterns to estimate the position of a transponder. Both cases, with and without knowledge of the vertical sound speed profile, were considered. A transponder positioning experiment was conducted at sea to collect acoustic and GPS observations. With significant uncertainties inherent in GPS measurements and the use of a commercial acoustic transponder not designed for precise ranging, experimental results indicate that the transponder position can be estimated accurately on the order of decimeters. Moreover, the mean sound speed of the water column estimated by the proposed optimization scheme is in agreement with that derived from conductivity, temperature, and density (CTD) measurements.     

Acoustic bubble counting technique using sound speed extracted fromsound attenuation

Sound attenuation has been solely used to estimate bubble size distributions of bubbly water in the conventional acoustic bubble sizing methods. These conventional methods are useful for the void fraction around 10^-6 or lower. However, the change of compressibility in the bubbly water also should be considered in bubble sizing for the void fraction around 10^-5 or higher. Recently the sound speed as well as sound attenuation was considered for acoustic bubble size estimation in bubbly water. In this paper, the sound speed estimated from sound attenuation in bubbly water by an iterative method is used for a bubble counting. This new iterative inverse bubble sizing technique is numerically tested for bubble distributions of single-size Gaussian, and power-law functions. The numerical simulation results are in agreement with the given bubble distributions even for the high void fractions of 10^-4-10^-3. It suggests that the iterative inverse technique can be a very powerful tool for practical use in acoustic bubble counting in the ocean     

Acoustic backscattering from a sandy seabed

Acoustic backscattering from a sandy seabed was measured at a frequency of 5.5 kHz at a wide range of grazing angles. The measurement system used was the University of Miami's sonar tower, consisting of an omni-directional broadband source and two 16-channel hydrophone receiver arrays. A volume scattering model, which combines a fluid model with reflection/transmission coefficients derived from the Biot theory, is used. This model allows energy penetration into the bottom, calculations of the volume scattering at all grazing angles, and the frequency dependence of the sound speed in the water-saturated sediment. In the model, rather than assume sound-speed correlation length in sedimentary volume, core data were used to assimilate a 3-D fluctuation spectrum of the density. The numerical results showed excellent agreement with the measurement at lower grazing angles. We concluded that the interface roughness scattering was dominant at lower grazing angles, while the volume scattering is dominant at higher grazing angles at the sandy site. The border of the dominance of the interface and volume scattering was the so-called critical angle at this frequency. The frequency dependence of sound speeds is also discussed.     

Kraken海洋声学模型及其声传播与衰减的数值试验

针对射线、简正波、PE、FFP等传播模型的算法原理及其适用的海洋环境,建立了以Kraken声学模型计算软件为基础的海洋声场数值预报系统.应用该预报系统对4组典型的海洋声场进行了数值试验,结果表明:在相同的海面和海底边界条件下,声场分布是由声速剖面和声源位置决定的.在负梯度声场中,所有声线都折向海底,在极限声线外产生阴影区.声源位于声道轴附近的温跃层中会产生波导传播.用射线理论解释了上述现象的成因,指出了其实际应用价值.     

Seismo-acoustic field statistics in shallow water

The spatial statistics of the acoustic field in shallow water are strongly affected by interfacial roughness and volume fluctuations in the water column or the seabed. These features scatter energy, reducing the coherence of the acoustic field. This paper introduces a consistent, mode-based modeling framework for ocean scattering. First, the rough surface scattering theory of Kuperman and Schmidt is reformulated in terms of normal modes, resulting in computation times which are reduced by several orders of magnitude. Next, a perturbation theory describing scattering from sound speed and density fluctuations in acoustic media is developed. The scattering theories are combined with KRAKEN, creating a unified normal mode code for wave theory modeling of shallow-water spatial statistics. The scattered field statistics are found to be a complicated function of scattering mechanism, scatterer statistics, and acoustic environment. Bottom properties, including elasticity, strongly influence the scattered field     

Acoustic plane-wave reflection from a rough surface over a random fluid half-space

Plane-wave reflection from a rough surface overlying a fluid half-space, with a sound speed distribution subject to a small and random perturbation, is considered. A theory based upon a boundary perturbation method in conjunction with a formulation derived from Green's function for the coherent field in the random medium have been applied to a typical oceanic environment to study their effects on the plane-wave reflection. By considering the coherent field itself, the plane-wave reflection may be obtained straightforwardly through a procedure consistent with the formalisms currently employed in rough surface scattering. The results show that both the rough surface and medium inhomogenieties may reduce the plane-wave reflection, however, the characteristics of the curves representing their effects are different, enabling us to identify the dominant scattering mechanism. The results for the coherent reflection due to the individual scattering mechanism are compatible with those found in the existing literature.     

振动系统的逆与优化问题(Part-Ⅱ)部分边界值再构内部场(Ⅰ)理论

研究由一振动体的部分边界测量数据再构振动体内部场的逆问题.提出基于偏微分方程控制方法的一种算法,并证明原用于处理散射问题的Null-field method 可用来求解优化问题的随伴问题.     

Assessment of acoustic iterative inverse method for bubble sizing to experimental data

Comparative study was carried out for an acoustic iterative inverse method to estimate bubble size distributions in water. Conventional bubble sizing methods consider only sound attenuation for sizing. Choi and Yoon [IEEE, 26(1), 125–130 (2001)] reported an acoustic iterative inverse method, which extracts the sound speed component from the measured sound attenuation. It can more accurately estimate the bubble size distributions in water than do the conventional methods. The estimation results of acoustic iterative inverse method were compared with other experimental data. The experimental data show good agreement with the estimation from the acoustic iterative inverse method. This iterative technique can be utilized for bubble sizing in the ocean.