云南壳幔S波速度结构与强震的构造背景

本文选取云南及周边65个台站记录到的47个地震事件,利用相匹配滤波技术分离出了基阶Rayleigh面波信号.选取与震中处于同一大圆弧上的两个台站,利用双台格林函数法获取了台间相速度频散,频散的周期范围在10～80 s之间.从2000个波形记录中提取了152个台站对之间的相速度频散,最后,利用台间的相速度频散反演得到云南...     

四川盆地中部浅层地壳一维剪切波速度结构初步研究

2010年1月31日四川盆地中部的遂宁发生了M5.0地震,四川盆地内部的台站观测到了很强的短周期面波波形数据。本文对地震数据的Rayleigh波和Love波提取基阶群速度频散曲线并反演得到四川盆地中部浅层地壳一维剪切波波速度结构。研究结果表明,四川盆地中部近地表剪切波平均速度约为2km/s,并且随深度逐渐增加,地壳深度在10km左右时,剪切波速度达到了3.5km/s。此结果适用于四川盆地中部以西,从遂宁到龙泉山附近,而四川盆地东部,从遂宁到华蓥山断裂不适用。该结果可为龙门山断裂附近的三维结构研究提供参考,并对强地面运动计算和区域内地震定位的研究有一定参考价值。     

基于改进Frequency-Bessel变换法的主动源Love波频散分析

近年来,Love波多道分析技术在浅地表结构探测中受到越来越广泛的应用.与Rayleigh波相比,Love波相速度不受纵波速度影响,其反演参数较少,可使反演过程更加稳定、求解的横波速度模型更加可靠.高分辨率和多模式面波频散分析是面波多道分析技术中至关重要的一环.本文利用改进Frequency-Bessel变换法对主动源Love波进行频散分析.通过改进的0阶Frequency-Bessel变换,将时间-空间域多道Love波记录变换到频率-波数域,获得其频散能量谱.通过公式推导和数值模拟,验证了该方法的有效性.实例测试表明,该方法具有较高的分辨率和多模式分辨能力,为多分量面波勘探提供了一种有效的频散分析手段.     

利用Rayleigh波反演浅土层的剪切波速度结构

根据均匀半空间Rayleigh波相速度与介质剪切波速度之间的关系,对某一频率面波的影响深度内各层介质进行"均匀"化,使其等效于均匀半空间．并利用模型进行正演,以确定在这种均匀化的前提下,面波勘探深度与波长的关系．应用相匹配滤波技术从实测面波信号中分离出基阶Rayleigh波信号,对它进行多重滤波和叠加处理,精确地计算出5．0-30．0Hz之间的基阶面波相、群速度频散．使用"均匀"化的方法,从相速度频散曲线中获得反演的初始模型,利用群速度频散反演得到35m以上各土层的剪切波速度结构．结果表明,反演结果与钻孔资料较为吻合     

基于贝叶斯理论的接收函数与环境噪声联合反演

基于Bayes反演理论(Tarantola,1987,2005),在接收函数非线性复谱比反演方法基础上(刘启元等,1996),本文讨论了接收函数与地震环境噪声Rayleigh波相速度频散的联合反演.本文采用修正后的快速广义反射/透射系数方法(Pei et al., 2008,2009) 计算Rayleigh波相速度频散, 并引入地壳泊松比的全局性搜索.数值检验表明:(1)接收函数与环境噪声的联合反演能够有效地解决反演结果对初始模型依赖的问题,即使对地壳速度结构仅有非常粗略的初始估计(例如,垂向均匀模型),本文方法仍能给出模型参数的可靠估计;(2)由于环境噪声与接收函数在频带上的适配性明显优于地震面波,接收函数与环境噪声的非线性联合反演能更好地约束台站下方近地表的速度结构;对于周期范围为2~40s的环境噪声相速度频散,利用本文方法能够可靠推测台站下方0~80 km深度范围的S波速度结构, 其浅表速度结构的分辨率可达到1 km; (3)本文方法能够可靠地估计地壳泊松比,泊松比的全局性搜索有助于合理解释接收函数和环境噪声的面波频散数据.利用本文方法对川西台阵KWC05台站观测的接收函数与环境噪声的联合反演表明,该台站下方地壳厚度为44 km,上地壳具有明显的高速结构,24~42 km范围的中下地壳具有低速结构.该台站下方地壳的平均泊松比为0.262,壳内低速带的泊松比为0.27.     

基于矢量波数变换法(VWTM)的多道Rayleigh波分析方法

在近二十年来,多道面波分析法(MASW)由于其便捷、高效等特性在浅层地震勘探领域得到了广泛的应用.本文基于多道面波勘探的采集方式,提出了一种新的面波多道分析方法——矢量波数变换法(VWTM).该方法通过对震源的近似,基于水平层状模型得到台站与震源间近似格林函数,然后进行矢量波数变换得到含有高阶模态Rayleigh波(频率-相速度)频散能量图.本研究首先利用合成地震数据到频散能量图与理论频散曲线进行叠加分析该方法的有效性和正确性;然后与相移法进行对比分析,我们发现在频散能量图中VWTM法对基阶、高阶模态成像均具有更高的分辨率和成像质量;最后我们将其应用于实际多道瞬态面波探测中,通过与相移法进行对比分析,发现VWTM法是一种方便、实用、有效的Rayleigh波频散提取方法.VWTM法提取多模态的Rayleigh波频散特征具有巨大潜力,可为基阶、高阶面波频散联合反演提供丰富的高阶模态频散信息.     

HTI煤层介质槽波波场与频散特征初步研究

煤层内裂隙较为发育,具有明显的各向异性.目前槽波理论研究以各向同性介质为主,对HTI介质中槽波及其频散性质研究很少.本文以弱各向异性、含垂直裂隙HTI煤层介质为研究对象,研究了HTI煤层介质中的三维槽波波场,采用交错网格高阶有限差分法模拟槽波,推导了三层水平层状HTI煤层介质的Love型槽波理论频散公式和振幅深度分布,分析了HTI各弹性参数对频散曲线的影响.HTI介质和各向同性介质基阶Love槽波频散曲线差异较小,高阶较大;煤厚主要影响Airy相频率,而Airy相速度不变;煤层vs对Airy相速度影响很大;煤层γ对基阶Love槽波影响很小,高阶稍大.各波偏振方向不再与波的传播方向平行或垂直,而是呈一定夹角.利用基阶Love槽波频散曲线推测裂隙发育较为困难,可利用高阶频散曲线.     

中国大陆及海域Love波层析成像

收集了研究区域(68°-150°E,5°-55°N)内33个数字地震台站记录的面波资料,利用多重滤波技术提取了4000余条路径上的Love波群速度频散曲线. 将研究区域划分成1°×1°网格,采取Occam反演方法得到了7.3-184s共43个周期的Love波群速度分布图;然后对网格结点进行S波速度结构反演,得出研究区域内420km深度内的地壳上地幔三维速度结构. 并采用Checkerboard方法对分辨率进行检验,得到横向的分辨率约为3°-5°. 研究结果表明:中国大陆地壳上地幔结构的横向不均匀性非常明显,内部结构与地表特征的相关性可以达到0-150km深度. 大陆地区东西分带、南北分块,块体的边界反映比较清晰.     

中国西部地区地壳上地幔S波速度结构

利用10～184 s基阶瑞利波频散曲线反演得到了中国西部及邻近区域地壳上地幔(0～300 km)的S波速度结构.反演采用了传统的两步法,即通过时频分析获得频散曲线后,首先利用Occam方法反演出各周期面波速度在二维网格结点上的分布,然后反演每个结点下方的S波速度结构,从而给出研究区域的三维速度结构.     

基于频率-贝塞尔变换法的关东盆地S波速度成像

将一种新的方法——频率-贝塞尔变换法（F-Jmethod）应用于日本NIED在关东盆地布设的MeSO-net台网的背景噪声数据中,证明频率-贝塞尔变换法可以有效地从背景噪声中提取高阶频散曲线.利用提取的基阶和高阶频散曲线反演关东盆地区域的沉积层和地震基岩层的S波速度结构,并将我们反演得到的S波速度结构与Koketsu等提出的日本综合速度结构模型进行对比讨论.我们的例子证明,在基阶面波的基础上,高阶面波能减少在反演中的非唯一性,得到更为准确的S波速度结构.     

唐山北部地区浅层一维剪切波速度结构初步研究

2019年8月2日河北省唐山市路南区某矿井区发生M_L2.4巷道塌陷,基于塌陷周边台站观测到的短周期瑞雷波,提取面波基阶群速度频散曲线,并利用迭代反演方法得到研究区域地下10 km深度范围内的一维剪切波速度结构,用于精定位分析。速度分析结果表明,研究区域浅表剪切波速度约为2.46 km/s;深度为2 km时,塌陷周边存在小范围的低速区,速度约为2.57 km/s;深度约为4 km时,剪切波速度达3.47 km/s;深度为5～9 km时,唐山东部沉积盆地内存在1个剪切波低速层。精定位分析结果表明,增加浅层速度模型有助于提高深度较小的地震定位精度;塌陷周边的低速区向下延伸近20 km,为地震多发区。     

高频面波方法的若干新进展

面波多道分析方法(MASW)通过分析高频瑞雷波确定浅地表剪切波速度.在过去的20年中,由于该方法具有非侵入性、无损、高效及价格低的特点,越来越受到浅地表地球物理和地质工程学界的重视,视为未来最有希望的技术之一.这篇综述论文将介绍中国地质大学(武汉)浅地表地球物理团队近年来在研究高频面波的传播理论和应用中取得的部分成果.非几何波是一种仅存在于浅地表介质,尤其是未固结的沉积物中的独特的地震波.它的存在对快速而准确地获得表层S波速度有一定价值.我们的研究表明非几何波是一种具有频散特性的泄漏波.泄漏波的存在可能导致将其误认为瑞雷波的基阶或高阶能量,从而造成模式误判.这种模式误判会导致错误的反演结果.我们通过求取高基阶分离后的瑞雷波格林函数证明虚震源法瑞雷波勘探的可行性.这个结果将极大地降低野外瑞雷波勘探成本.勒夫波多道分析方法(MALW)中未知参数比瑞雷波的少,这使得勒夫波的频散曲线比瑞雷波的简单.因此,勒夫波反演更稳定,非唯一性更低.勒夫波数据生成的能量图像通常比瑞雷波的清晰,并具有更高的分辨率,从而可以更容易地拾取精确的勒夫波的相速度.利用雅克比矩阵分析波长与探测深度的关系表明对相同波长的基阶模式而言,瑞雷波的探测深度是勒夫波的1.3~1.4倍;而两种波的相同波长的高阶模式波的探测深度相同.我们也尝试了时间域勒夫波反演.按照勒夫波分辨率将地球模型剖分成了不同尺寸的块体,利用反卷积消除了地震子波对勒夫波波形的影响,通过更新每个块体的S波速度来拟合勒夫波波形,从而获得地下S波速度模型.该方法不基于水平层状模型假设,适用于任意二维介质模型.     

Feasibility of waveform inversion of Rayleigh waves for shallow shear-wave velocity using a genetic algorithm

Conventional surface wave inversion for shallow shear (S)-wave velocity relies on the generation of dispersion curves of Rayleigh waves. This constrains the method to only laterally homogeneous (or very smooth laterally heterogeneous) earth models. Waveform inversion directly fits waveforms on seismograms, hence, does not have such a limitation. Waveforms of Rayleigh waves are highly related to S-wave velocities. By inverting the waveforms of Rayleigh waves on a near-surface seismogram, shallow S-wave velocities can be estimated for earth models with strong lateral heterogeneity. We employ genetic algorithm (GA) to perform waveform inversion of Rayleigh waves for S-wave velocities. The forward problem is solved by finite-difference modeling in the time domain. The model space is updated by generating offspring models using GA. Final solutions can be found through an iterative waveform-fitting scheme. Inversions based on synthetic records show that the S-wave velocities can be recovered successfully with errors no more than 10% for several typical near-surface earth models. For layered earth models, the proposed method can generate one-dimensional S-wave velocity profiles without the knowledge of initial models. For earth models containing lateral heterogeneity in which case conventional dispersion-curve-based inversion methods are challenging, it is feasible to produce high-resolution S-wave velocity sections by GA waveform inversion with appropriate priori information. The synthetic tests indicate that the GA waveform inversion of Rayleigh waves has the great potential for shallow S-wave velocity imaging with the existence of strong lateral heterogeneity.     

基阶与高阶瑞利波联合反演研究

研究了六层层状介质模型瑞利波基阶和高模式波相速度对横波速度、深度的敏感性,结果表明：基阶波较高模式波对7 m以内浅部地层的横波速度更敏感,敏感性频带在10～25 Hz范围内,峰值频带集中在18 Hz左右;高模式波较基阶波对深部地层的横波速度更敏感,敏感性频带宽,峰值分散.基阶波对浅层的敏感性和高模式波穿透深度更深的特点为近地表岩土层二维横波速度结构的联合反演提供了前提条件.利用阻尼最小二乘SVD（Singular Value Decomposition）算法联合基阶与高模式波对理论模型和实例数据进行横波速度反演,反演结果表明联合反演增强了反演的稳定性,提高了反演的精度.     

Seismic crust structure beneath the Aegean region in southwest Turkey from radial anisotropic inversion of Rayleigh and Love surface waves

The Turkish plate is covered by hundreds of accelerometer and broadband seismic stations with less than 50 km inter-station distance providing high-quality earthquake recordings within the last decade. We utilize part of these stations to extract the fundamental mode Rayleigh and Love surface wave phase and group velocity data in the period range 5–20 s to determine the crust structure beneath the Aegean region in southwest Turkey. The observed surface wave signals are interpreted using both single-station and two-station techniques. A tomographic inversion technique is employed to obtain the two-dimensional group velocity maps from the single-station group velocities. One-dimensional velocity–depth profiles under each two-dimensional mesh point, which are jointly interpreted to acquire the three-dimensional image of the shear-wave velocities underneath the study area, are attained by utilizing the least-squares inversion technique, which is repeated for both Rayleigh and Love surface waves. The isotropic crust structure cannot jointly invert the observed Rayleigh and Love surface waves where the radial anisotropic crust better describes the observed surface wave data. The intrusive magmatic activity related to the northward subducting African plate under the Turkish plate results the crust structure deformations, which we think, causing the observed radial anisotropy throughout complex pattern of dykes and sills. The magma flow resulting in the mineral alignment within dykes and sills contributes to the observed anisotropy. Due to the existence of dykes, the radial anisotropy in the upper crust is generally negative, i.e., vertically polarized S-waves (Vsv) are faster than horizontally polarized S-waves (Vsh). Due to the existence of sills, the radial anisotropy in the middle-to-lower crust is generally positive, i.e., horizontally polarized S-waves (Vsh) are faster than vertically polarized S-waves (Vsv). Similar radial anisotropic results to those of the single-station analyses are obtained by the two-station analyses utilizing the cross-correlograms. The widespread volcanic and plutonic rocks in the region are consistent with the current seismic interpretations of the crustal deformations.     

Shallow geology characterization using Rayleigh and Love wave dispersion curves derived from seismic noise array measurements

The local geology and shallow S-wave velocity structure of a site are recognized to be key factors for the increase in the damaging potential of seismic waves. Indeed, seismic amplitudes may be amplified in frequency ranges unfavorable for building stock by the presence of soft sedimentary covers over lying hard bedrock. Hence, microzonation activities, which aim at assessing the site response as accurately as possible, have become a fundamental task for the seismic risk reduction of urbanized areas. Methods based on the measurement of seismic noise, which typically are fast, non-invasive, and low cost, have become a very attractive option in microzonation studies.Using observations derived from seismic noise recordings collected by two-dimensional arrays of seismic stations, we present a novel joint inversion scheme for surface wave curves. In particular, the Love wave, the Rayleigh wave dispersion and the HVSR curves are innovatively combined in a joint inversion procedure carried out following a global search approach (i.e., the Genetic Algorithm).The procedure is tested using a data set of seismic noise recordings collected at the Bevagna (Italy) test-site. The results of the novel inversion scheme are compared with the inversion scheme proposed by Parolai et al. (2005), where only Rayleigh wave dispersion and HVSR curves are used, and with a cross-hole survey.     

Dispersion Energy Analysis of Rayleigh and Love Waves in the Presence of Low-Velocity Layers in Near-Surface Seismic Surveys

High-frequency surface-wave analysis methods have been effectively and widely used to determine near-surface shear (S) wave velocity. To image the dispersion energy and identify different dispersive modes of surface waves accurately is one of key steps of using surface-wave methods. We analyzed the dispersion energy characteristics of Rayleigh and Love waves in near-surface layered models based on numerical simulations. It has been found that if there is a low-velocity layer (LVL) in the half-space, the dispersion energy of Rayleigh or Love waves is discontinuous and ‘‘jumping’’ appears from the fundamental mode to higher modes on dispersive images. We introduce the guided waves generated in an LVL (LVL-guided waves, a trapped wave mode) to clarify the complexity of the dispersion energy. We confirm the LVL-guided waves by analyzing the snapshots of SH and P–SV wavefield and comparing the dispersive energy with theoretical values of phase velocities. Results demonstrate that LVL-guided waves possess energy on dispersive images, which can interfere with the normal dispersion energy of Rayleigh or Love waves. Each mode of LVL-guided waves having lack of energy at the free surface in some high frequency range causes the discontinuity of dispersive energy on dispersive images, which is because shorter wavelengths (generally with lower phase velocities and higher frequencies) of LVL-guided waves cannot penetrate to the free surface. If the S wave velocity of the LVL is higher than that of the surface layer, the energy of LVL-guided waves only contaminates higher mode energy of surface waves and there is no interlacement with the fundamental mode of surface waves, while if the S wave velocity of the LVL is lower than that of the surface layer, the energy of LVL-guided waves may interlace with the fundamental mode of surface waves. Both of the interlacements with the fundamental mode or higher mode energy may cause misidentification for the dispersion curves of surface waves.     

Advantages of Using Multichannel Analysis of Love Waves (MALW) to Estimate Near-Surface Shear-Wave Velocity

As theory dictates, for a series of horizontal layers, a pure, plane, horizontally polarized shear (SH) wave refracts and reflects only SH waves and does not undergo wave-type conversion as do incident P or Sv waves. This is one reason the shallow SH-wave refraction method is popular. SH-wave refraction method usually works well defining near-surface shear-wave velocities. Only first arrival information is used in the SH-wave refraction method. Most SH-wave data contain a strong component of Love-wave energy. Love waves are surface waves that are formed from the constructive interference of multiple reflections of SH waves in the shallow subsurface. Unlike Rayleigh waves, the dispersive nature of Love waves is independent of P-wave velocity. Love-wave phase velocities of a layered earth model are a function of frequency and three groups of earth properties: SH-wave velocity, density, and thickness of layers. In theory, a fewer parameters make the inversion of Love waves more stable and reduce the degree of nonuniqueness. Approximating SH-wave velocity using Love-wave inversion for near-surface applications may become more appealing than Rayleigh-wave inversion because it possesses the following three advantages. (1) Numerical modeling results suggest the independence of P-wave velocity makes Love-wave dispersion curves simpler than Rayleigh waves. A complication of “Mode kissing” is an undesired and frequently occurring phenomenon in Rayleigh-wave analysis that causes mode misidentification. This phenomenon is less common in dispersion images of Love-wave energy. (2) Real-world examples demonstrated that dispersion images of Love-wave energy have a higher signal-to-noise ratio and more focus than those generated from Rayleigh waves. This advantage is related to the long geophone spreads commonly used for SH-wave refraction surveys, images of Love-wave energy from longer offsets are much cleaner and sharper than for closer offsets, which makes picking phase velocities of Love waves easier and more accurate. (3) Real-world examples demonstrated that inversion of Love-wave dispersion curves is less dependent on initial models and more stable than Rayleigh waves. This is due to Love-wave’s independence of P-wave velocity, which results in fewer unknowns in the MALW method compared to inversion methods of Rayleigh waves. This characteristic not only makes Love-wave dispersion curves simpler but also reduces the degree of nonuniqueness leading to more stable inversion of Love-wave dispersion curves.     

3-D Velocity Structure of the Middle and Upper Crust in the Yunnan Region,China

Using short-period (1–18 s) surface wave data recorded by 23 stations of the Yunnan Digital Seismic Network of China we determined phase velocities of the fundamental Rayleigh wave along 209 paths by the two-station narrowband filtering and cross-correlation method, followed by an inversion for phase velocity distributions at various periods using the Ditmar-Yanovskaya method. We then obtained a 3-D S-wave velocity structure of the middle and upper crust in the Yunnan region using the genetic algorithm. The results show strong lateral variation of phase velocity in the region. The short-period phase velocity variation is closely related to thickness variation of sedimentary layer in the shallow crust. Within the depth range of 26–30 km, the S-wave velocity in the Sichuan-Yunnan rhombic block is lower than in the surrounding areas. Most large earthquakes of M > 6.0 in Yunnan occurred in the transition zones between low and high velocities.     

3-D surface wave group velocity distribution in Central-Southern Africa

Group velocities estimated from fundamental mode Love and Rayleigh waves are used in a tomography process in central-southern Africa. The waves were generated by eighteen earthquakes, which occurred along the East African Rift and recorded at BOSA, LBTB and SLR seismic stations in southern Africa. The group velocities from Love and Rayleigh waves were isolated using the Multiple Filter Technique (MFT) at the period range of 10 to 50 seconds. The tomography method developed by Ditmar and Yanovskaya (1987) and Yanovskaya and Ditmar (1990), was applied to calculate the lateral distribution of surface wave group velocities in central-southern Africa. The results of the tomographic inversion were plotted as distribution maps. In addition to the maps, I also produced two velocity cross-sections across the area of study. The velocity distribution maps show the regional tectonic units, though with poor resolution. The azimuthal bias of the surface wave paths is reflected in the distribution of the group velocities. The Moho depth appears to correlate with velocities at a period of about 30 s. A low velocity feature observed beneath the Zimbabwe craton implies a thickening upper asthenosphere and lithospheric thinning beneath the Zimbabwe craton. Also estimated was a shear wave velocity model beneath the Zimbabwe craton.