Estimation of phase and group velocities for multi-modal ground roll using the ‘ phase shift ’ and ‘ slant stack generalized S transform based ’ methods
Estimation of phase and group velocities for multi-modal ground roll using the ‘ phase shift ’ and ‘ slant stack generalized S transform based ’ methods
复制标题
DOI:
--
复制
发表时间:
2011
期刊:
影响因子:
--
通讯作者:
R. Askari;R. Ferguson;K. DeMeersman
中科院分区:
文献类型:
--
作者:
R. Askari;R. Ferguson;K. DeMeersman
Phase and group velocities are two important factors that determine shear wave velocity in Multi-channel Analysis of Surface Wave (MASW) surveys. In this study, we present two different methods. The phase shift method and the slant stack generalized S transform based method for the estimation of the phase and group velocities respectively. The phase shift method uses the idea of the slant stack in the Fourier domain to estimate the phase velocity. The slant stack generalized S transform based method uses the slant stack idea in the time-frequency domain based on the generalized S transform. These methods are robust to estimate phase and group velocities where ground roll is multimodal. We anticipate that, through inversion of the phase and group velocities of multimodal ground rolls, a better estimation of near surface shear wave velocity will be obtainable. INTRODUCTION One of the most important goals in engineering studies is the estimation of soil rigidity for the near surface (Xia et al., 2002a). Due to the dependency of soil rigidity and shear wave velocity, accurate estimation of shear wave velocity is a key factor in engineering studies. Although, seismic refraction methods are one of the conventional methods used for shear wave velocity estimation (Palmer, 1980), they fail to estimate S-velocity where geological structure is complex (Xia et al., 2002b) or where hidden layers (a layer whose velocity is less than its upper layer) are present (Sheriff and Geldart, 1986). As an alternative to the refraction methods, the analysis of dispersed surface waves (Rayleigh and Love waves) is a well-known procedure to estimate shear wave velocity (e.g. Evison et al., 1959; Stokoe et al., 1988; KeilisBorok, 1989; Lay and Wallace, 1995; Xia et al., 1999). The method is a based on the inversion of a surface wave dispersion curve (either phase or group velocity curves) calculated from a real data to a vertical shear wave velocity profile. Some studies indicate that the method overcomes the pitfalls of the refraction methods such as geological complexity (Xia et al., 2002b) or the presence of hidden layers (Feng et al., 2005). Since surface wave methods are based on the inversion of phase and group velocities to a shear wave velocity profile, the processing of data in order to truly estimate a dispersion curve is very crucial. Many algorithms have been developed to address this issue -p transform (McMechan and Yedlin, 1981), the wavelet transform (Kulesh et al., 2005; Holschneider et al., 2005; Kulesh et all, 2008) and the generalized S transform (Askari and Ferguson; 2011) for phase velocity, and narrow band-pass filtering (Herrmann, 1973), the wavelet transform (Kulesh et al., 2005; Holschneider et al., 2005) and the generalized S transform (Askari and Ferguson, 2011) for group velocity. Askari, Ferguson, and DeMeersman 2 CREWES Research Report — Volume 23 (2011) Multi-channel Analysis of Surface Waves (MASW) is a surface wave method based on multi-channel seismic recorded ground roll for near surface studies (Park, 1999). In this study, we present two methods for the estimation of the phase and group velocities for the MASW survey respectively. The phase shift (Park et al., 1998) is a method which is based on the estimation of the Fourier phase spectrum difference from one trace to another for the estimation of the phase velocity. Then, we introduce a slant stack model based on the generalized S transform (Pinnegar and Mansinha, 2003) to estimate the group velocity of multi-channel ground roll. PHASE SHIFT METHOD Considering a seismic record in the time-offset domain u(t,x) containing the ground roll, the Fourier transform for each trace is expressed U(x, ) = u(x, t)e dt. (1) The equation above could be rewritten as