The study solves 2D elastodynamic model for seismic in-plane wave propagation in laterally inhomogeneous geological profiles situated in a vertically inhomogeneous half-space with a seismic source. For the aim, an efficient hybrid Modal SummationBoundary Integral Equation Method (MSM-BIEM) is applied. The MSM is used as a tool for the simulation of wave propagation from the source position to the local multilayered laterally inhomogeneous geological profile, where the BIEM is applied. The BIEM is based on the frequency-dependent fundamental solution of the governing equation in elastodynamics and the hybrid method works in the frequency domain. The inverse FFT solution is applied to obtain time histories. The hybrid tool is applied to several models for the investigation of local site effects due to: (a) the impedance contrasts between soil layers, (b) surface topography and lateral inhomogeneity, (c) the seismic source properties and (d) the existence of water saturation in soils. The application of the modeling tool is a contribution to the seismic risk analysis of Sofia city.
Abbreviations
RMRay methods FDM Finite difference method MSM Modal summation method BIEM Boundary integral equation method DBIEM Direct boundary integral equation method MS-BIEM Modal summation-boundary integral equation method FFT Fast Fourier transform CPV integrals Cauchy Principal Value integrals List of symbolsVertically inhomogeneous half-space modeled with a series of M homogeneous flat layers C i , parallel to the free surfaceCommon boundary between X and X 0 q k , k k , l k Density and lame constants of each kth layer in X 0 a k , b k Complex longitudinal and shear wave velocities of each anelastic kth layer in X 0 V Pk , V Sk Real part of complex longitudinal and shear wave velocities of each anelastic kth layer in X 0