Scattering of surface waves modelled by the integral equation method

The integral equation method is used to model the propagation of surface waves in 3-D structures. The wavefield is represented by the Fredholm integral equation, and the scattered surface waves are calculated by solving the integral equation numerically. The integration of the Green's function elements is given analytically by treating the singularity of the Hankel function at   R = 0  , based on the proper expression of the Green's function and the addition theorem of the Hankel function. No far-field and Born approximation is made. We investigate the scattering of surface waves propagating in layered reference models imbedding a heterogeneity with different density, as well as Lamé constant contrasts, both in frequency and time domains, for incident plane waves and point sources.     

An analytical approach for the scattering of SH waves by a symmetrical V-shaped canyon: shallow case

Based on the application of the region-matching technique, an analytical approach is presented for the scattering of plane SH waves from a shallow symmetrical V-shaped canyon, and then a series solution is derived. The analysed region is divided into an enclosed and an open region by introducing a semi-circular auxiliary boundary. In each region, the displacement field can be expressed as infinite sum of appropriate wavefunctions satisfying partial boundary conditions, respectively. The unknown coefficients can be determined by enforcing the continuity conditions in connection with the Graf's addition formula. The frequency- and time-domain responses are both evaluated and displayed for several physical parameters. From graphical results, the effects of the canyon depth on surface ground motion are conspicuous. The proposed series solutions can serve as benchmark for numerical methods, in particular for those at much higher frequencies.     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.