Erik H. Saenger

Bio: Erik H. Saenger is an academic researcher from Ruhr University Bochum. The author has contributed to research in topics: Wave propagation & Attenuation. The author has an hindex of 32, co-authored 189 publications receiving 4243 citations. Previous affiliations of Erik H. Saenger include Karlsruhe Institute of Technology & ETH Zurich.

Finite‐difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid

Abstract: We describe the application of the rotated staggered-grid (RSG) finite-difference technique to the wave equations for anisotropic and viscoelastic media. The RSG uses rotated finite-difference operators, leading to a distribution of modeling parameters in an elementary cell where all components of one physical property are located only at one single position. This can be advantageous for modeling wave propagation in anisotropic media or complex media, including high-contrast discontinuities, because no averaging of elastic moduli is needed. The RSG can be applied both to displacement-stress and to velocity-stress finite-difference (FD) schemes, whereby the latter are commonly used to model viscoelastic wave propagation. With a von Neumann-style anlysis, we estimate the dispersion error of the RSG scheme in general anisotropic media. In three different simulation examples, all based on previously published problems, we demonstrate the application and the accuracy of the proposed numerical approach.

Accuracy of heterogeneous staggered-grid finite-difference modeling of Rayleigh waves

Abstract: Heterogeneous finite-difference (FD) modeling assumes that the boundary conditions of the elastic wavefield between material discontinuities are implicitly fulfilled by the distribution of the elastic parameters on the numerical grid. It is widely applied to weak elastic contrasts between geologic formations inside the earth. We test the accuracy at the free surface of the earth. The accuracy for modeling Rayleigh waves using the conventional standard staggered-grid (SSG) and the rotated staggered grid (RSG) is investigated. The accuracy tests reveal that one cannot rely on conventional numerical dispersion discretization criteria. A higher sampling is necessary to obtain acceptable accuracy. In the case of planar free surfaces aligned with the grid, 15 to 30 grid points per minimum wavelength of the Rayleigh wave are required. The widely used explicit boundary condition, the so-called image method, produces similar accuracy and requires approximately half the sampling of the wavefield compared to heterogeneous free-surface modeling. For a free-surface not aligned with the grid (surface topography), the error increases significantly and varies with the dip angle of the interface. For an irregular interface, the RSG scheme is more accurate than the SSG scheme. The RSG scheme, however, requires 60 grid points per minimum wavelength to achieve good accuracy for all dip angles. The high computation requirements for 3D simulations on such fine grids limit the application of heterogenous modeling in the presence of complex surface topography.

This is How the Discussion Started

Abstract: A copy of Guangbo jiemu bao [Broadcast Program Report] was being passed from hand to hand among a group of young people eager to be the first to read the article introducing the program "What Is Revolutionary Love?" It said: "… Young friends, you are certainly very concerned about this problem'. So, we would like you to meet the young women workers Meng Xiaoyu and Meng Yamei and the older cadre Miss Feng. They are the three leading characters in the short story ‘The Place of Love.’ Through the description of the love lives of these three, the story induces us to think deeply about two questions that merit further examination.

Ground Motion Characteristics Estimated from Spectral Ratio between Horizontal and Verticcl Components of Mietremors.

Abstract: The spectral ratio between horizontal and vertical components (H/V ratio) of microtremors measured at the ground surface has been used to estimate fundamental periods and amplification factors of a site, although this technique lacks theoretical background. The aim of this article is to formulate the H/V technique in terms of the characteristics of Rayleigh and Love waves, and to contribute to improve the technique. The improvement includes use of not only peaks but also troughs in the H/V ratio for reliable estimation of the period and use of a newly proposed smoothing function for better estimation of the amplification factor. The formulation leads to a simple formula for the amplification factor expressed with the H/V ratio. With microtremor data measured at 546 junior high schools in 23 wards of Tokyo, the improved technique is applied to mapping site periods and amplification factors in the area.

Seismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks — A review

Abstract: One major cause of elastic wave attenuation in heterogeneous porous media is wave-induced flow of the pore fluid between heterogeneities of various scales. It is believed that for frequencies below 1 kHz, the most important cause is the wave-induced flow between mesoscopic inhomogeneities, which are large compared with the typical individual pore size but small compared to the wavelength. Various laboratory experiments in some natural porous materials provide evidence for the presence of centimeter-scale mesoscopic heterogeneities. Laboratory and field measurements of seismic attenuation in fluid-saturated rocks provide indications of the role of the wave-induced flow. Signatures of wave-induced flow include the frequency and saturation dependence of P-wave attenuation and its associated velocity dispersion, frequency-dependent shear-wave splitting, and attenuation anisotropy. During the last four decades, numerous models for attenuation and velocity dispersion from wave-induced flow have been developed with varying degrees of rigor and complexity. These models can be categorized roughly into three groups according to their underlying theoretical framework. The first group of models is based on Biot’s theory of poroelasticity. The second group is based on elastodynamic theory where local fluid flow is incorporated through an additional hydrodynamic equation. Another group of models is derived using the theory of viscoelasticity. Though all models predict attenuation and velocity dispersion typical for a relaxation process, there exist differences that can be related to the type of disorder periodic, random, space dimension and to the way the local flow is incorporated. The differences manifest themselves in different asymptotic scaling laws for attenuation and in different expressions for characteristic frequencies. In recent years, some theoretical models of wave-induced fluid flow have been validated numerically, using finite-difference, finite-element, and reflectivity algorithms applied to Biot’s equations of poroelasticity. Application of theoretical models to real seismic data requires further studies using broadband laboratory and field measurements of attenuation and dispersion for different rocks as well as development of more robust methods for estimating dissipation attributes from field data.