Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.     

On the Acoustic-Electromagnetic Analogy for the Reflection-Refraction Problem

The same mathematical theory can be used to describe physical phenomena of different nature. For instance, the wave equation and the related mathematical developments can be used to describe elastic and electromagnetic wave propagation, and it is also extensively used in quantum mechanics. Fresnel's equations are a classical example of the analogy between shear waves and light waves. George Green in the nineteenth century, used analogies to obtain the reflection coefficients for sound waves and light waves, before the advent of the electromagnetic theory of light.In this work, we investigate the mathematical analogy between elastic waves and electromagentic waves. We obtain a complete parallelism for the reflection and refraction problem, considering the most general situation, that is, the presence of anisotropy and attenuation—viscosity in the elastic case and conductivity in the electromagnetic case. The analogy is illustrated with Fresnel's equations, the Brewster and critical angles, the concept of reflectivity and transmissivity, and the corresponding duals fields. The analysis of the elastic-solid theory of reflection applied by Green to light waves, and a brief historical review of wave propagation through the ether, further illustrate the analogy.     

Acoustic and electromagnetic properties of soils saturated with salt water and NAPL

Electrical, seismic, and electromagnetic methods can be used for noninvasive determination of subsurface physical and chemical properties. In particular, we consider the evaluation of water salinity and the detection of surface contaminants. Most of the relevant properties are represented by electric conductivity, P-wave velocity, and dielectric permittivity. Hence, it is important to obtain relationships between these measurable physical quantities and soil composition, saturation, and frequency. Conductivity in the geoelectric frequency range is obtained with Pride's model for a porous rock. (The model considers salinity and permeability.) White's model of patchy saturation is used to calculate the P-wave velocity and attenuation. Four cases are considered: light nonaqueous phase liquid (LNAPL) pockets in water, dense nonaqueous phase liquid (DNAPL) pockets in water, LNAPL pockets in air, and DNAPL pockets in air. The size of the pockets (or pools), with respect to the signal wavelength, is modeled by the theory. The electromagnetic properties in the GPR frequency range are obtained by using the Hanai–Bruggeman equation for two solids (sand and clay grains) and two fluids (LNAPL or DNAPL in water or air). The Hanai–Bruggeman exponent (1/3 for spherical particles) is used as a fitting parameter and evaluated for a sand/clay mixture saturated with water.Pride's model predicts increasing conductivity for increasing salinity and decreasing permeability. The best-fit exponent of the Hanai–Bruggeman equation for a sand/clay mixture saturated with water is 0.61, indicating that the shape of the grains has a significant influence on the electromagnetic properties. At radar frequencies, it is possible to distinguish between a water-saturated medium and a NAPL-saturated medium, but LNAPL- and DNAPL-saturated media have very similar electromagnetic properties. The type of contaminant can be better distinguished from the acoustic properties. P-wave velocity increases with frequency, and has dissimilar behaviour for wet and dry soils.     

Estimation of gas-hydrate concentration and free-gas saturation at the Norwegian-Svalbard continental margin

We estimate the concentration of gas hydrate and free gas at an area located to the north of the Knipovich Ridge (western Svalbard margin). The method is based on P-wave velocities computed by reflection tomography applied to multicomponent ocean-bottom seismometer data. The tomographic velocity field is fitted to theoretical velocities obtained from a poro-elastic model based on a Biot-type approach (the interaction between the rock frame, gas hydrate and fluid is modelled from first physical principles). We obtain average hydrate concentrations of 7% and maximum free-gas saturations of 0.4% and 9%, depending on the saturation model.     

Comparison of P‐wave attenuation models of wave‐induced flow

Wave‐induced oscillatory fluid flow in the vicinity of inclusions embedded in porous rocks is one of the main causes for P‐wave dispersion and attenuation at seismic frequencies. Hence, the P‐wave velocity depends on wave frequency, porosity, saturation, and other rock parameters. Several analytical models quantify this wave‐induced flow attenuation and result in characteristic velocity–saturation relations. Here, we compare some of these models by analyzing their low‐ and high‐frequency asymptotic behaviours and by applying them to measured velocity–saturation relations. Specifically, the Biot–Rayleigh model considering spherical inclusions embedded in an isotropic rock matrix is compared with White's and Johnson's models of patchy saturation. The modeling of laboratory data for tight sandstone and limestone indicates that, by selecting appropriate inclusion size, the Biot‐Rayleigh predictions are close to the measured values, particularly for intermediate and high water saturations.     

Approximating constant‐Q seismic propagation in the time domain

In this study, we investigate the accuracy of approximating constant‐Q wave propagation by series of Zener or standard linear solid (SLS) mechanisms. Modelling in viscoacoustic and viscoelastic media is implemented in the time domain using the finite‐difference (FD) method. The accuracy of numerical solutions is evaluated by comparison with the analytical solution in homogeneous media. We found that the FD solutions using three SLS relaxation mechanisms as well as a single SLS mechanism, with properly chosen relaxation times, are quite accurate for both weak and strong attenuation. Although the RMS errors of FD simulations using a single relaxation mechanism increase with increasing offset, especially for strong attenuation (Q = 20), the results are still acceptable for practical applications. The synthetic data of the Marmousi‐II model further illustrate that the single SLS mechanism, to model constant Q, is efficient and sufficiently accurate. Moreover, it benefits from less computational costs in computer time and memory.