Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

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.     

Comparison of Seismic Dispersion and Attenuation Models

The frequency-dependent attenuation of seismic waves causes decreased resolution of seismic images with depth, and the difference in transmission losses induces amplitude variations with offset. Transmission losses may occur due to friction or fluid movement, or may result from scattering in thin-layer. Whatever the physical mechanism, they can often be conveniently described using an empirical formulation wherein the elastic moduli and propagation velocity are complex functions of frequency.We have compiled and compared algebraically and numerically eight different models involving complex velocity: the Kolsky-Futterman model, the power-law model, Kjartansson's model, Müller's model, Azimi's second and third model, the Cole-Cole model, and the standard linear-solid model.For two different parameter sets, the attenuation and phase velocity are computed in the seismic frequency band, and the plane-wave propagation of a Ricker wavelet for the other models is compared with that for the Kolsky-Futterman model. The first parameter set consists of parameters for each of the models calculated from expressions given in the appendix. These expressions make the different models behave similarly to the KF model. The second parameter set consists of model parameters that are numerically adapted to the KF model.By selecting proper parameters, all models, except the standard linear-solid model, show behavior similar to that of the Kolsky-Futterman model. The SLS model behaves differently from the other models as the frequency goes to zero or infinity. Broadband measurement data is needed to select a specific model for a given seismic experiment.     

Erratum: Comparison of Seismic Dispersion and Attenuation Models: Stud. Geophys. Geod., 46(2002), 293-320

Studia Geophysica et Geodaetica -     

Countergradient Fluxes of Conserved Variables in the Clear convective and Stratocumulus-topped Boundary Layer: The role of the Entrainment Flux

Large-eddy simulations of a clear convective boundary layer (CBL)and a stratocumulus-topped boundary layer are studied. Bottom-upand a top-down scalars were included in the simulations, and theprinciple of linear superposition of variables was applied toreconstruct the fields of any arbitrary conserved variable.This approach allows a systematic analysis of countergradient fluxesas a function of the flux ratio, which is defined as the ratio betweenthe entrainment flux and the surface flux of the conserved quantity.In general, the turbulent flux of an arbitrary conserved quantityis counter to the mean vertical gradient if the heights where thevertical flux and the mean vertical gradient change sign do notcoincide. The regime where the flux is countergradient is thereforebounded by the so-called zero-flux and zero-gradient heights. Becausethe vertical flux changes sign only if the entrainment flux has anopposite sign to the surface flux, countergradient fluxes arepredominantly found for negative flux ratios. In the CBL the fluxratio for the virtual potential temperature is, to a good approximation,constant, and equal to -0.2. Only if the moisture contribution to thevirtual potential temperature is negligibly small will the flux ratio forthe potential temperature be equal to this value. Otherwise, theflux ratio for the potential temperature can have any arbitrary(negative) value, and, as a consequence, the fluxes for thepotential temperature and the virtual potential temperature willbe countergradient at different heights. As a practical application ofthe results, vertical profiles of the countergradient correction termfor different entrainment-to-surface-flux ratios are discussed.