凹陷地形对Rayleigh面波传播影响的研究

利用局域离散波数法,研究了地形对近地表爆炸源产生的Rayleigh面波的影响,展示了Rayleigh面波穿越地形过程中波的传播进程,研究了穿越地形前后Rayleigh面波能量、频率的变化,并讨论了该变化与地形起伏程度的关系,从中可见,Rayleigh面波穿越地形时,激发出新的体波和面波,并且越陡的地形对面波传播的阻碍性越强,且面波高频部分能量损失的越多.     

The velocity structure of the Carpathian zone from the ambient noise surface wave tomography

The dispersion curves of the Rayleigh wave group velocities are constructed along 60 interstation seismic paths in Central Europe based on the cross-correlation function of seismic noise. Together with the previous data (Yanovskaya and Lyskova, 2013), this information was used for reconstructing the three-dimensional distribution of S-wave velocities in the upper mantle of the Carpathian region. In the present work, the previous results are refined by expanding the data set by the additional seismic paths that intersect the Carpathian region and by modifying the procedure for constructing the locally averaged dispersion curves so as to obtain a more compact resolution. The results of the study suggest the complex, multidirectional character of the plate motion in the region.     

基于噪声面波直接反演法研究天山中段地壳S波速度结构

天山造山带是新生代以来复活隆升的陆内造山带,强烈的地震活动性使得理解和认识天山造山带深部结构及盆山耦合关系尤为重要。文章中使用天山造山带及邻区(40°～49°N,79°～93°E)85个台站2017—2019年的背景噪声资料,结合背景噪声互相关方法获得了6～52 s瑞利波相速度频散曲线,利用基于射线追踪的面波直接反演法对天山中段地壳三维S波速度结构及盆山耦合关系进行研究。结果显示：地壳浅层S波速度分布与构造单元中沉积层厚度相关,塔里木盆地北缘、准噶尔盆地南缘表现为低速,天山造山带表现为高速;到了中下地壳,天山造山带下方存在被高速异常包裹的低速体;莫霍面附近,天山造山带表现出相对低速;准噶尔盆地南缘和天山造山带的地壳厚度分别在45～50 km、50～62 km之间,沿南北向,天山造山带莫霍面呈现较为宽缓的形态;在82°～86.5°E之间,塔里木盆地和准噶尔盆地向天山下方双向俯冲,86.5°～88°E之间,准噶尔盆地向天山南向俯冲,由西向东,不同盆山耦合关系揭示了新生代以来天山中段不同区域构造运动差异,为进一步探讨造山动力过程提供参考。     

The velocity structure of the upper mantle of europe from the ambient noise surface wave tomography

The method for surface wave tomography based on the records of ambient seismic noise (Ambient Noise Tomography, ANT) is applied to the data from the East European and West European stations. In order to reduce the effects of the earthquakes at long periods, the cross correlation functions were calculated for the time interval of 2001–2003, when distinct clusters of the earthquakes were absent. Using the local dispersion curves in the range of 10–100 s, we reconstructed the vertical velocity sections at the nodes of the 3° × 3° grid. On the basis on these curves, we calculated the horizontal distributions of S-velocity variations in the upper mantle in the depth interval of 75–275 km and the vertical velocity sections along the profiles across the Vrancea zone and the region of the Baltic and Ukrainian shields. The velocity distribution in the Vrancea zone confirms the subduction of the ancient oceanic plate from the east westwards and the detachment of its bottom part, as hypothesized by some authors. Beneath the Baltic Shield lithosphere, there is a low-velocity zone, which can be interpreted as the asthenospheric layer. It is noted that the velocity distributions beneath the Baltic and Ukrainian shields are similar, which probably points to the genetic relationship between these two structures.     

Love wave tomography in Italy from seismic ambient noise

We estimate Love wave empirical Green's functions from cross-correlations of ambient seismic noise to study the crust and uppermost mantle structure in Italy.Transverse-component ambient noise data from October 2005 through March 2007 recorded at 114 seismic stations from the Istituto Nazionale di Geofisica e Vulcanologia (INGV) national broadband network,the Mediterranean Very Broadband Seismographic Network (MedNet) and the Austrian Central Institute for Meteorology and Geodynamics (ZAMG) yield more than 2 000 Love wave group velocity measurements using the multiple-filter analysis technique.In the short period band (5-20 s),the cross-correlations show clearly one-sided asymmetric feature due to non-uniform noise distribution and high local activities,and in the long period band (20 s) this feature becomes weak owing to more diffusive noise distribution.Based on these measurements,Love wave group velocity dispersion maps in the 8-34 s period band are constructed,then the SH wave velocity structures from the Love wave dispersions are inverted.The final results obtained from Love wave data are overall in good agreement with those from Rayleigh waves.Both Love and Rayleigh wave inversions all reveal that the Po plain basin is resolved with low velocity at shallow depth,and the Tyrrhenian sea is characterized with higher velocity below 8 km due to its thin oceanic crust.     

基于背景噪声的川滇地区勒夫波层析成像

本文选取了川滇地区98个固定台站记录到的三分量地震背景噪声数据,采用互相关方法提取了勒夫波互相关函数,并利用自适应时频分析方法获取了勒夫波群速度频散曲线,经反演得到周期为8—30 s的勒夫波群速度分布图像.层析成像结果显示:短周期的勒夫波群速度分布图像呈明显的横向不均匀性,且与地表地质和构造特征基本一致,其中四川盆地呈不均匀的低速异常,盆地内成都平原地区的群速度相对低于盆地中部的丘陵地区,速度分界线为遂宁与峨眉山之间的连线,四川盆地内的群速度变化反映出沉积层厚度的变化情况;攀枝花地区呈高速异常,可能与古地幔活动有关,幔源物质以侵入岩和底侵岩浆的形式停留在地壳的不同深度,从而形成高速异常的特征.本文结果为了解川滇地区的构造运动提供了地震学线索,并为下一步研究地壳径向各向异性奠定了基础.      

华北地区勒夫波噪声层析成像研究

利用华北地震科学台阵200个流动地震台站、14个月的连续波形数据,通过互相关方法提取了勒夫波的经验格林函数,使用多重滤波方法测量了5229条勒夫波的群速度频散曲线,采用噪声层析成像方法得到了研究区域4~30 s的勒夫波群速度分布图像,横向分辨率在多数区域可以达到0.25°×0.25°. 层析成像结果显示,短周期的群速度分布特征与地表地质和构造特征基本一致,华北盆地和山西断陷带内的盆地呈现低速异常,燕山隆起和太行山隆起表现为高速异常;中周期的群速度分布图揭示了华北盆地内部隆起和坳陷的空间分布范围及沉积层的厚度差异. 勒夫波频散曲线具有明显的分区特征,太行山隆起、燕山隆起和鄂尔多斯块体东北缘的频散曲线形态基本一致,和其它典型克拉通相似;张渤地震带和山西裂谷盆地的频散曲线基本一致,接近于埃塞俄比亚裂谷. 中西部块体的平均频散曲线和其它典型克拉通相似,而东部块体的频散曲线和中西部块体存在较大差异,且低于其它典型克拉通,表明东部块体的地壳受到了强烈的破坏和改造,而中西部块体受改造的程度较低,仍具有稳定克拉通的物理性质. 太行山重力梯度带的东西两侧在地形地貌、速度结构、频散特征、地壳厚度、岩石圈厚度以及地幔过渡带厚度均存在显著差异,是划分华北克拉通破坏空间范围的一条重要界线,推测华北克拉通的破坏范围主要在太行山重力梯度带以东.     

Three dimensional shear wave velocity structure of the crust and upper mantle beneath China from ambient noise surface wave tomography

We determine the three-dimensional shear wave velocity structure of the crust and upper mantle in China using Green's functions obtained from seismic ambient noise cross-correlation.The data we use are from the China National Seismic Network,global and regional networks and PASSCAL stations in the region.We first acquire cross-correlation seismograms between all possible station pairs.We then measure the Rayleigh wave group and phase dispersion curves using a frequency-time analysis method from 8 s to 60 s.After that,Rayleigh wave group and phase velocity dispersion maps on 1° by 1° spatial grids are obtained at different periods.Finally,we invert these maps for the 3-D shear wave velocity structure of the crust and upper mantle beneath China at each grid node.The inversion results show large-scale structures that correlate well with surface geology.Near the surface,velocities in major basins are anomalously slow,consistent with the thick sediments.East-west contrasts are striking in Moho depth.There is also a fast mid-to-lower crust and mantle lithosphere beneath the major basins surrounding the Tibetan plateau (TP) and Tianshan (Junggar,Tarim,Ordos,and Sichuan).These strong blocks,therefore,appear to play an important role in confining the deformation of the TP and constraining its geometry to form its current triangular shape.In northwest TP in Qiangtang,slow anomalies extend from the crust to the mantle lithosphere.Meanwhile,widespread,a prominent low-velocity zone is observed in the middle crust beneath most of the central,eastern and southeastern Tibetan plateau,consistent with a weak (and perhaps mobile) middle crust.     

Crust and upper mantle structure of East Asia from ambient noise and earthquake surface wave tomography

The complex tectonic background of East Asia makes it an ideal region for investigating the evolution of the continental lithosphere, for which high-resolution lithospheric structural models are essential. In this study, we measured Rayleigh-wave phase-velocity dispersion curves at periods of 10–120 s and group velocity dispersion curves at periods of 10–140 s using event records from more than 1,000 seismic stations in and around China. By jointly inverting new and previously published dispersion data from ambient noise and earthquakes, we developed a high-resolution shear-wave velocity model down to a depth of ～300 km beneath East Asia. Our model revealed heterogeneous lithospheric structures beneath East Asia, and allowed us to investigate the velocity structure of the entire lithosphere. We also derived crustal and lithospheric thickness models from the three-dimensional (3D) shear-wave model, revealing strong spatial heterogeneity and a general thinning trend of lithospheric thickness from west to east across the study region. Overall, our models reveal important lithospheric features beneath East Asia and provide a valuable baseline dataset for understanding continental-scale dynamics and evolution.     

基于背景噪声和地震面波联合反演华北克拉通中部岩石圈结构

基于ChinArray三期项目布设于华北克拉通中部的流动台阵观测数据,利用背景噪声互相关和地震面波层析成像获取了研究区内6—140 s周期的瑞雷面波频散,使用蒙特卡罗非线性反演方法获得了华北克拉通中部岩石圈的高分辨率三维S波速度结构。结果显示华北克拉通不同地块的岩石圈速度结构存在显著的横向差异：其中鄂尔多斯盆地腹地整体表现为高速特征,延伸至200 km以下,但其东南缘存在小范围的低速异常;东部的华北盆地整体表现为低速特征,具有较薄的地壳和岩石圈厚度;中部造山带南北两端以及南北重力梯度线下方存在相连接的低速区域,在深处延伸至华北盆地下方;在下地壳和上地幔顶部,大同火山群区域的低速体逐渐向西偏移至鄂尔多斯盆地东北角下方;而在上地幔中,该区域的低速异常随深度增加而逐渐减弱,低速体延伸至东南方向的华北盆地下方。基于本研究获得的S波速度模型,我们认为：鄂尔多斯盆地腹地保持了克拉通特性,但其东南缘存在局部的岩石圈改造作用;华北盆地发生了强烈的岩石圈破坏减薄和地壳伸展变形;中部造山带南北端以及南北重力梯度线下方的岩石圈发生了局部的改造减薄,其机制可能都来源于华北盆地下方地幔热物质的上涌;大同火山群下...     

新疆地区环境噪声层析成像研究

本文利用北京大学宽频带流动地震台阵和新疆地震监测台网12个月连续地震记录数据,采用环境噪声层析成像方法获得了新疆地区的周期从10~35 s范围内的瑞利面波相速度异常分布图像.研究结果与地表地质构造相一致,新疆地区的天山、两盆大地构造与瑞利面波相速度异常有较好的对应关系.短周期10~20 s的图像显示两个盆地内都呈现低速异常,对应盆地较厚的沉积层.中长周期25~35 s的层析成像结果显示东西天山地区存在明显异常差异.本文获得的新疆地区地壳和上地幔相速度结果,也为进一步结合传统天然地震面波成像方法,最终获得新疆地区的三维剪切波速度结构奠定基础.     

福建地区环境噪声瑞利波群速度层析成像

本文利用福建、广东、江西和浙江等4省数字地震监测台网中69个宽频带台站2010年4月13日～10月13日的速度型噪声记录和互相关技术,提取两两台站间的瑞利面波格林函数,并利用层析成像技术最终反演获得了福建地区(115.7～120.5°E、23.6～28.35°N)瑞利面波群速度的分布图像。结果表明,利用该方法得到的福建地区周期为3～5s的瑞利面波群速度在2.73～3.17km/s之间,平均速度为2.96km/s。瑞利波群速度横向变化与福建地区地表构造特征和地热异常区分布情况是基本吻合的;该群速度分布图在漳州盆地和福州盆地区域表现出非常明显的低速异常,这主要是因为受到盆地沉积层和高地热的影响;杉岭地块属于山地地形,但由于受到高地热的影响,因此表现为低速异常;永安-玳瑁山地块为山地区,而受到中地热的影响,速度异常表现较为复杂。     

地震波模拟中多轴复频移近似完全匹配层的吸收效果及稳定性

在进行地震波模拟计算的过程中用有限的计算区域模拟地下无限空间,需要进行边界截断.为了在边界处不产生虚假反射影响模拟结果,需要引入吸收边界条件.本文采用的近似完全匹配层是一种新型非分裂完全匹配层,计算效率较高.同时相比于其他非分裂完全匹配层,其还具有不改变方程的形式、易于实现等优势.但是当入射波角度较大,边界吸收效果变弱,且残留在边界中的能量使近似完全匹配层变得极其不稳定.多轴复频移近似完全匹配层的提出就是为了改善对大角度入射波的吸收并且提高边界的稳定性.通过实验模拟和矩阵特征值灵敏度来研究多轴复频移近似完全匹配层的吸收效果及稳定性.结果表明该方法不仅能够吸收掠入波,而且对常规入射波的吸收也得到提升,同时拥有更好的稳定性.

     

Fault detection by reflected surface waves based on ambient noise interferometry

Detecting subsurface fault structure is important for evaluating potential earthquake risks associated with active faults. In this study, we propose a new method to detect faults using reflected surface waves observed in ambient noise cross correlation functions. Ambient noise tomography using direct surface waves obtained from ambient noise interferometry has been widely used to characterize active fault zones. In cases where a strong velocity contrast exists across the fault interface, fault-reflected surface waves are expected. We test this idea using a linear array deployed in the Suqian segment of Tanlu fault zone in Eastern China. The fault-reflected surface waves can be clearly seen in the cross-correlation functions of the ambient noise data, and the spatial position of the fault on the surface is close to the stations where the reflected signals first appear. Potentially reflected surface waves could also be used to infer the dip angle, fault zone thickness and the degree of velocity contrast across the fault by comparing synthetic and observed waveforms.     

2D surface topography boundary conditions in seismic wave modelling

New formulations of boundary conditions at an arbitrary two-dimensional (2D) free-surface topography are derived. The top of a curved grid represents the free-surface topography while the grid's interior represents the physical medium. The velocity–stress version of the viscoelastic wave equations is assumed to be valid in this grid. However, the rectangular grid version attained by grid transformation is used to model wave propagation in this work in order to achieve the numerical discretization. We show the detailed solution of the particle velocities at the free surface resulting from discretizing the boundary conditions by second-order finite-differences (FDs). The resulting system of equations is spatially unconditionally stable. The FD order is gradually increased with depth up to eighth order inside the medium. Staggered grids are used in both space and time, and the second-order leap-frog and Crank–Nicholson methods are used for time-stepping. We simulate point sources at the surface of a homogeneous medium with a plane free surface containing a hill and a trench. Applying parameters representing exploration surveys, we present examples with a randomly realized surface topography generated by a 1D von Kármán function of order 1. Viscoelastic simulations are presented using this surface with a homogeneous medium and with a layered, randomized medium realization, all generating significant scattering.     

On the wave propagation in inhomogeneous non-isotropic media

Summary FollowingEason, we have discussed here the propagation of elastic waves in non-homogeneous spheres and cylinders when the curved surface is given a uniform normal loading. The material is assumed to be transversely isotropic with respect to a direction of symmetry, the stress and displacement components within the body may be assumed to depend on one space co-ordinate and time alone. The particular case in which the elastic parameters are proportional to (radius) ⁿ has been considered.     

中国大陆中东部地区基于背景噪声的瑞利波层析成像

本研究使用了中国大陆中东部地区494个分布基本均匀的宽频带地震台站和7个中国大陆周边地区IRIS台站资料,反演得到了中东部地区高分辨率的瑞利面波层析成像结果.本文使用这些台站记录到的从2009年1月到2010年9月的垂直分量连续波形数据,首先通过对台站对间进行波形互相关和叠加运算,计算得到各台站对间的经验格林函数.然后用时频分析法提取了约125000条台站对间的频散曲线,并剔除了经验格林函数信噪比小于10的频散曲线.最后反演得到了研究区周期8~40 s、分辨率达0.5°的瑞利波群速度分布图像.不同周期的速度分布图像显示,研究区瑞利波群速度分布与地质构造特征具有较好的相关性.8~20 s的瑞利波群速度在研究区内主要盆地表现为低速分布,而在造山带呈现高速分布;25~40 s的瑞利波群速度图中,存在一条北北东—南南西向的分界线,该分界线与中国大陆东部的地壳厚度突变带基本吻合.25 s以下周期,华北平原的显著低速区形态与该地区早第三纪以来的断块分布构造一致.揭示了盆地下方介质结构强烈的非均匀性,也与较厚的沉积层分布有关.低速的四川盆地中部,显示出显著的高速特征,揭示了四川盆地下方基底的上隆特征;20 s以下周期的群速度图像中,鄂尔多斯盆地西北部速度低于东南部,揭示出其地壳中上部介质结构的横向不均匀性.