Benchmarking geoacoustic inversion methods using range-dependent field data

Over the past decade, inversion methods have been developed and applied to acoustic field data to provide information about unknown ocean-bottom environments. An effective inversion must provide both an estimate of the bottom parameters and a measure of the uncertainty of the estimated values. This paper summarizes results from the Office of Naval Research (ONR)/Space and Naval Warfare Systems Command (SPAWAR) Geoacoustic Inversion Techniques Workshop, test cases 4 and 5. The workshop was held to benchmark present-day inversion methods for estimating geoacoustic profiles in shallow water. The format of the workshop was a blind test to estimate unknown geoacoustic profiles by inversion of measured acoustic transmission loss data in octave bands and reverberation envelopes. The data sets for test cases 4 and 5 were taken at two locations in shallow water, one in the East China Sea and the other along the southwest coast of Florida. The limitations of the data and the limits to the knowledge of the sites are discussed. In both cases, impulsive sources were used in conjunction with air-deployed sonobuoys. Since the measured data was incoherent, only methods consistent with total energy matching were applicable. Comparisons between the different inversion techniques presented at the workshop are discussed. For test cases 4 and 5, a precise metric was unavailable for comparison.     

Combined Geoacoustic Inversion of Propagation and Reverberation Data

Sonar performance predictions in shallow water are strongly dependent on good knowledge of the geoacoustic and scattering properties of the seabed. One technique to extract information about the bottom is to use a towed source and a towed horizontal array. This towed system has been shown to be applicable for characterizing the bottom properties locally by inversion of the acoustic signals received directly on the towed array at short ranges. The same towed system has also been applied to extract bottom properties from long-range reverberation data providing effective bottom properties over a large area. However, independent geoacoustic inversion of the short-range propagation and long-range reverberation data can introduce low sensitivity and uncertainty in the extracted bottom properties. An attempt to resolve this low sensitivity and ambiguity is made by a simultaneous geoacoustic inversion of short-range propagation and long-range reverberation data with the intention of constraining the possible solutions of the bottom properties.      

多频海底声学原位测试系统研制和试用

海底沉积物的声学特性(最重要的是声速和声衰减)以及它们与物理(包括土力学)特性之间的关系是沉积物声学中两个重要的研究项目.介绍了新研制的实时监控多频海底声学原位测试系统.该系统可测量浅表层沉积物的声速.探测频率为8,10,12,15 kHz,可根据实际情况选择发射波形、接收增益和采样长度,采样率为0.5~2.0 MHz,工作水深为300 m.系统具有倾斜传感器、8通道扩充等功能.用该系统在杭州湾测得了四种频率的沉积物原位声速.     

Geoacoustic seafloor exploration with a towed array in a shallowwater area of the Strait of Sicily

Acoustic propagation in shallow water is greatly dependent on the geoacoustic properties of the seabottom. This paper exploits this dependence for estimating geoacoustic sediment properties from the bottom acoustic returns of known signals received on a hydrophone line array. There are two major issues in this approach: one is the feasibility of acoustic inversion with a limited aperture line array, the other is related to the knowledge of the geometry of the experimental configuration. To test the feasibility of this approach, a 40-hydrophone-4-m spaced towed array together with a low-frequency acoustic source, was operated at a shallow water site in the Strait of Sicily. In order to estimate the array deformation in real time, it has been equipped with a set of nonacoustic positioning sensors (compasses, tilt-meters, pressure gauges). The acoustic data were inverted using two complementary approaches: a genetic algorithm (GA) like approach and a radial basis functions (RBF) inversion scheme. More traditional methods, based on core sampling, seismic survey and geophone data, together with Hamilton's regression curves, have also been employed on the same tracks, in order to provide a ground truth reference environment. The results of the experiment, can be summarized as follows: 1) the towed array movement is not negligible for the application considered and the use of positioning sensors are essential for a proper acoustic inversion, 2) the inversion with GA and RBF are in good qualitative agreement with the ground truth model, and 3) the GA scheme tends to have better stability properties. On the other hand, repeated in version of successive field measurements requires much less computational effort with RBF     

Perturbative inversion of geoacoustic parameters in a shallow waterenvironment

In many strategic shallow water areas, the geoacoustic properties of the sub-bottom are largely unknown. This paper demonstrates that inverse theory and measured data from a single hydrophone can be used to accurately deduce the geoacoustic properties of the sub-bottom, even when the initial background geoacoustic model is a highly inaccurate estimate. Since propagation in shallow water is very sensitive to the geoacoustic properties of the sub-bottom the inverse technique is a vitally important, practical, and inexpensive means to improve sonar performance prediction in a potentially hostile environment. To provide ground truth for the inverse technique, measured data collected during Project GEMINI were compared to the inverse solutions. Detailed, site-specific geoacoustic models were developed for two array locations and the finite-element parabolic equation (FEPE) model was used to estimate transmission loss (TL). The model estimates from FEPE compared well with the measured data and the detailed geoacoustic models were considered as “ground truth.” To test the efficacy of the technique, initial background geoacoustic models were constructed assuming no a priori information of the bottom. The resultant inverse solution was used to predict the geoacoustic properties at each of the sites. The final results were in excellent agreement with the measured data and the resulting TL estimates derived from the inverse technique were as good or better than the TL estimates obtained from the detailed, site-specific geoacoustic models     

Reverberation vertical coherence and sea-bottom geoacoustic inversion in shallow water

Optimal array-processing techniques in the ocean often require knowledge of the spatial coherence of the reverberation. A mathematical model is derived for the reverberation vertical coherence (RVC) in shallow water (SW). A method for analysis of RVC data is introduced. Measured reverberation cross-correlation coefficients as a function of time and frequency, obtained during the Asian Seas International Acoustic Experiment (ASIAEX) in the East China Sea, are reported. SW reverberation from a single shot provides a continuous spatial sampling of the surrounding sound field up to several tens of kilometers and holds valuable information on the geoacoustic properties of the sea floor over this distance. SW reverberation data can, therefore, be used as the basis for a quick and inexpensive method for geoacoustic inversion and has the obvious advantage that acquiring the data in situ requires only a single platform. This paper considers the use of the vertical coherence of the reverberation as the starting point for such an inversion. Sound speed and attenuation in the sea bottom at the ASIAEX site are obtained over a frequency range of 100-1500 Hz by finding values that provide the best match between the measured and predicted RVC.     

The Effects of Microstructure on Shear Properties of Shallow Marine Sediments

This study was undertaken to investigate the implication of geoacoustic behaviors in the shallow marine sediments associated with the changes in geotechnical index properties. Two piston cores (270 cm and 400 cm in core length) used in this study were recovered from stations 1 and 2, the western continental margin, the East Sea. Scanning electron microscopy (SEM) was employed to illustrate the effects of microstructure on shear properties. The direct SEM observation of sediment fabrics is inevitable to understand the correlation of the changes in geoacoustic properties to the sediment structure. The consolidation of sediments by overburden stress resulting in the clay fabric alteration appears to play an important role in changing shear properties. Water contents and porosity of sediments gradually decreases with increasing depth, whereas wet bulk density shows a reverse trend. It is interesting to note that shear wave velocities increase rapidly from 8 to 20 m/s while compressional wave velocities significantly fluctuate, ranging from 1450 to 1550 m/s with depth. The fabric changes in sediment with increasing depth for example, uniform grain size and well oriented clay fabrics may cause the shear strength increase from 1 to 12 kPa. Shear wave velocity is, therefore, shown to be very sensitive to the changes in undrained strength for unconsolidated marine sediments. This correlation allows an in-situ estimation of shear stress in the subsurface from shear wave velocity data.     

海水深度对两种可控源电磁勘探方法的影响

在过去的十几年中,海洋可控源电磁法的迅猛发展在油气勘探中起到了重要的作用。在浅海环境中,常见的两种可控源电磁法的数据采集系统采用海底基站式（Sea Bed Logging（SBL））,以及拖缆式（Towed Streamer electromagnetics（TSEM））两种模式。本文的研究目的是研究在不同水深条件下,两种数据采集系统的工作模式以及海水对其的影响。为了更好的比较两种数据采集方式,本文采用了三种不同的分析方法：一维敏感度模型,二维油藏敏感指数,Occam反演。一维敏感度模型能够迅速的给出结果,但是由于结果的信息量有限,不足以给定准确的结论。二维油藏敏感指数可以依据正演的结果快速的给定一个水深和油藏敏感度的关系。Occam反演方法用于对合成数据进行反演计算,估计地下介质电阻率参数。为了更好的评价反演结果,本文引入了横向电阻异常率（ATR）。通过对合成数据的反演结果的分析,本文得出以下结论：在浅水环境中,SBL和TSEM数据采集系统都可以提供较高的分辨率,对高电阻层有很好的敏感度。在深水环境中,SBL数据采集系统具有明显的优势。ATR对于反演的电阻率参数结果具有很好的评价作用。     

In situ and laboratory geoacoustic measurements in soft mud and hard-packed sand sediments: Implications for high-frequency acoustic propagation and scattering

Near-surface sediment geoacoustic and physical properties were measured in gas-rich, muddy sediments of Eckernförde Bay, Baltic Sea, and in hard-packed, sandy sediments of the northeastern Gulf of Mexico. Values of compressional and shear wave velocity are much lower in muddy compared to sandy sediments. The spatial and temporal variability of sediment physical and geoacoustic properties and, as a consequence, the scattering and propagation of high-frequency acoustic waves are primarily related to the presence and absence of free methane gas bubbles at the muddy site and to the abundance and distribution of shell material on sandy sediments.     

Generating interface waves using a freely falling, instrumentedsource

In recent years, interface waves such as the Scholte wave have become important tools in the study of the geoacoustic properties of near-bottom seafloor sediments. Traditionally, these waves have been generated by explosive or pneumatic sources deployed at or near the seafloor and monitored by ocean-bottom seismographs or geophone arrays. While these sources generate the requisite interface waves, they also produce higher frequency compressional waves in the water and sediment that tend to contaminate the surface wave and make inversion of the data difficult in the near field. In this paper, a new source consisting of a freely falling projectile instrumented with an accelerometer is described. When the projectile impacts the bottom, the exact time history of the vertical force applied to the sediment is known and therefore may be convolved with the transfer function of a sediment geoacoustic model to produce accurate synthetic seismograms. Moreover, the vertical force applied to the seafloor is very efficient in generating surface wave motion while producing very little compressional wave energy so that the near-field signals are much more easily analyzed. An example of the use of the new source is presented including inversion of the received signals to obtain shear-wave velocity and attenuation as a function of depth in the near bottom sediments at a shallow-water site     

Matched-field inversion for geoacoustic model parameters in shallowwater

Matched-field inversion is used to, estimate geoacoustic properties from data obtained in an experiment with a vertical line array (VLA). The experiment was carried out using broad-band sources (shots) in water depths of about 200 m on the continental shelf off Vancouver Island. The data were processed to obtain spectral components of the field for frequencies near the bubble frequency for the shot. The ocean bottom in this region consists of a layer of mainly sandy sediments (about 100 m thick) overlying older consolidated material. Consequently, the inversion was designed to estimate the parameters of a two-layer elastic sediment model. In the inversion, an adaptive global search algorithm was used to investigate the multidimensional space of geoacoustic models in order to determine the set of values corresponding to the best replica field. Convergence is driven by adaptively guiding the search to regions of the parameter space associated with above-average values of the matched field correlation between the measured and replica fields. The geoacoustic profile estimated by the inversion consisted of a 125-m layer with compressional speed ~1700 m/s and shear speed ~400 m/s, overlying a layer with compressional speed ~1900 m/s. This model is consistent with the results from conventional seismic experiments carried out in the same region     

基于浅地层剖面的海底浅表层沉积物物理性质参数反演技术研究

海底浅表层（小于1 m）沉积物的物理性质,如粒度、孔隙度、密度等是海洋沉积学研究和海洋工程地质分析的重要内容,目前主要基于有限的海底取样或原位测试获取这些沉积物的物理性质。浅地层剖面是基于声学信号（频率几千赫兹）在沉积物中的传播得到可反映沉积地层结构的数据,其中的一些声学参数,如海底反射系数、波阻抗等与沉积物物理性质密切相关。如何充分而有效地利用浅地层剖面资料,反演得到剖面覆盖区海底浅表层沉积物的物理性质参数,极具科学意义和应用价值,且基于声学属性反演沉积物物理性质是当前研究的热点。为此,本文基于渤海LD16-3CEPA至LD10-1PAPD路由段的浅地层剖面数据和海底表层沉积物的实测物理参数,利用Biot-Stoll模型建立研究区海底反射系数和沉积物物理性质之间的关系,并基于浅地层剖面数据计算得到的海底反射系数,反演了研究区海底浅表层沉积物的孔隙度、密度、平均粒径等物理性质参数。其中反演的孔隙度、密度、平均粒径与实测孔隙度、密度、平均粒径基本相符,偏差度基本都在20%的偏差范围内,表明该反演方法在该区的应用是可行的。     

海底沉积物声学特征定量分析及其智能分类研究

海底底质特性描述及分类是当今浅海声学的研究热点,海底沉积物的物理结构特性与其声学响应特征密切相关。在分析海底沉积物声传播特性的基础上,应用现代计算机信号分析技术手段,对海底沉积物声学响应波形提取了4个特征参数:声速、波幅指数、波形关联维分形指数和声波频谱的频率矩。以这4个特征参数作为输入向量,海底沉积物的结构类型作为输出向量,建立径向基概率神经网络模型。研究表明建立的神经网络模型具有较强的海底沉积物分类预报能力。     

Effects of sea states on seafloor compliance studies

Gas hydrates affect the bulk physical properties of marine sediments, in particular, elastic parameters. Shear modulus is an important parameter for estimating the distribution of hydrates in the marine sediments. However, S-wave information is difficult to recover without proper datasets. Seafloor compliance, the transfer function between pressure induced by surface gravity waves and the associated seafloor deformation, is one of few techniques to study shear modulus in the marine sediments. The coherence between recorded time series of displacement and pressure provides a measure of the quality of the calculated transfer function, the seafloor compliance. Thus, it is important to understand how to collect high coherence datasets. Here we conducted a 10-month pilot experiment using broadband seismic sensors and differential pressure gauges. We found that data collected in shallow water depth and during rough seas gave high coherence. This study is the first time long-term data sets have been employed to investigate seafloor compliance data quality and its dependence on sea state. These results will help designing future large-scale compliance experiments to study anomalously high shear moduli associated with the presence of gas hydrate or cold vents, or alternatively anomalously low shear moduli, associated with partial melt and magma chamber.     

Estimates of geoacoustic model parameters from inversions of horizontal and vertical line array data

This paper describes results from geoacoustic inversion of low-frequency acoustic data recorded at a receiving array divided into two sections, a sparse bottom laid horizontal array (HLA) and a vertical array (VLA) deployed in shallow water. The data are from an experiment conducted by the Norwegian Defence Research Establishment (FFI) in the Barents Sea, using broadband explosives (shot) sources. A two-layer range-independent geoacoustic model, consistent with seismic profiles from the area, described the environment. Inversion for geoacoustic model parameters was carried out using a fast implementation of the hybrid adaptive simplex simulated annealing (ASSA) inversion algorithm, with replica fields computed by the ORCA normal mode code. Low-frequency (40-128 Hz) data from six shot sources at ranges 3-9 km from the array were considered. Estimates of sediment and substrate p-wave velocities and sediment thickness were found to be consistent between independent inversions of data from the two sections of the array.     

Geoacoustic Uncertainties From Viscoelastic Inversion of Seabed Reflection Data

This paper applies nonlinear Bayesian inversion to seabed reflection data to estimate viscoelastic parameters of the upper sediments. The inversion provides maximum a posteriori probability (MAP) parameter estimates with uncertainties quantified in terms of marginal probability distributions, variances, and credibility intervals; interparameter relationships are quantified by correlations and joint marginal distributions. The inversion is applied to high-resolution reflectivity data from two sites in the Strait of Sicily. One site is characterized by low-speed sediments, resulting in data with a well-defined angle of intromission; the second is characterized by high-speed sediments, resulting in a critical angle. Data uncertainties are quantified using several approaches, including maximum-likelihood (ML) estimation, treating uncertainties as nuisance parameters in the inversion, and analysis of experimental errors. Statistical tests are applied to the data residuals to validate the assumed uncertainty distributions. Excellent results (i.e., small uncertainties) are obtained for sediment compressional-wave speed, compressional attenuation, and density; shear parameters are less well determined although low shear-wave speeds are indicated. The Bayesian analysis provides a quantitative comparison of inversion results for the two sites in terms of the resolution of specific geoacoustic parameters, and indicates that the geoacoustic information content is significantly higher for intromission data     

含砂量变化影响海底沉积物压缩波速度的分析研究

本文通过对南海海底沉积物样品的声学物理参数和沉积粒度特征统计分析,发现了高、低含砂量沉积物的声学物理特征存在明显差异,建立了海底沉积物的含砂量与压缩波速度、孔隙度、含水量和密度等经验公式,分析了含砂量变化与沉积物的体积压缩模量和密度变化的关系,从声速理论基础上阐明了含砂量变化引起沉积物压缩波速度变化的内在原因是含砂量变化引起了体积压缩模量和密度发生了变化,说明了含砂量增大引起沉积物压缩波速度增大的内在原因是含砂量增大引起了体积压缩模量变化量大于密度变化量,从而在数据统计和理论分析结合基础上,论证了含砂量是影响海底沉积物压缩波速度的重要因素之一。这一研究对声学方法反演海底沉积物粒度参数和沉积物类型、地声参数转换模型的建立,以及对水声反演海底和海底资源勘探等方面都具有重要理论意义和应用价值。     

Gassy sediment occurrence and properties: Northern Gulf of Mexico

Free gas is ubiquitous at shallow sediment depths of the northern margin of the Gulf of Mexico. Gassy sediment patches are between 250 and 500 m in horizontal size. Often the gassy layers are within 100 m from the sea floor and are only a few meters thick. Both biogenic and thermogenic gas hydrates have been recovered. Stability values of temperature and pressure indicate that hydrates can exist in water depths less than 500 m. Gassy sediment geoacoustic parameter values are not well constrained because of a lack of concurrent measurements of acoustic properties and sediment gas content. For Gulf of Mexico gassy sediment, some reportedin situ values of sound speed are reduced by an order of magnitude below values for water saturated sediments. More commonly, sound speed is reduced from water saturated sediment values by only 15 to 50 percent.     

The dynamic balance of manganese transport and diagenesis in the Otago Harbour,New Zealand

The input and internal cycling of manganese in the Otago Harbour, a shallow marine embayment, have been studied. Freshwater transport of dissolved and particulate Mn was measured during sustained steady flows and during a major flood event during which particulate mass transport increased by orders of magnitude. Flux calculations derived from these data revealed that input to the harbour system is dominated by particulate transport during flood events occupying a minor part of the flow record. In contrast, under average flow conditions the major input was in the form of dissolved Mn. Within the harbour, Mn transport in the water column was dominated by essentially conservative mixing of dissolved Mn formed through diagenetic remobilisation of riverine particulate Mn deposited in sediments in the upper basin of the harbour. The distribution of Mn in suspended matter showed a marked enrichment compared to bottom sediments indicative of oxidative precipitation of Mn(II) onto particle surfaces in the water column and/or selective resuspension processes. Within the margins of uncertainty in various flux estimates, the loss rate of Mn from the harbour balanced the terrestrial source function, illustrating the dynamics of Mn transport and diagenesis important to this element's oceanic geochemistry.     

Contaminant transport modeling in marine clays

The need to develop various offshore containment systems encourages the contaminant transport study in marine clays. In the present investigation, an attempt has been made to examine six cases of contaminant transport in marine clays. Adequate literature background on the chemical diffusion of soils, and the earlier reported modeling techniques relevant to the present work have been highlighted. These methods enable a rapid examination of the impact of contaminant physical properties on the environment over a period of time. The numerical results of this study help to understand the chemicals transport phenomena in the marine environment especially at sea bed.