基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.     

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.

     

S-velocity structure beneath the Bohemian Massif from Monte Carlo inversion of seismic receiver function

The paper contains an analysis of S-velocity distribution in the crust and upper mantle beneath the Bohemian Massif, which is the second biggest Variscan outcrop in Europe. It occupies mainly the west part of Czech Republic and also part of south-west Poland and south-east Germany. We use data from 10 permanent stations set in the region. Some previous papers relate to the same scope but use linear methods to inverse receiver function. Our new approach involves Monte Carlo techniques for inversion procedure, which is more convenient and robust for such a non-linear task. The result of Monte Carlo inversion is compared with the previously achieved one. The obtained Moho depths vary from 29 km in the north-west part of the Bohemian Massif to 38 km in the south and south-east and are consistent with other papers. Some discrepancies occur in the middle and upper crust.     

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

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

2010年高雄地震震源参数的近远震波形联合反演

本文改进了传统基于近震波形数据的点源震源参数反演的Cut And Paste(CAP)方法,实现了近震Pnl波、面波和远震P波、SH波的联合反演的CAPjoint算法.对2010年3月高雄地震,分别进行单独反演以及联合反演,获得各自的震源机制解及深度,其中联合反演所得的最佳双力偶机制解参数为,节面1:走向317°,倾角36°,滑移角52°,节面2:走向181°,倾角62°,滑移角114°,深度为21 km.并对不同震中距波形对本次地震以及几种典型机制解断层几何参数的敏感性进行测试.为验证联合反演方法的可靠性,本文采用重抽样思想发展而来的Bootstrap方法,对近震数据的子集及其与远震数据的联合反演所得的参数进行统计,验证了在稀疏近台条件下联合反演中添加远震数据对地震震源参数约束的作用.     

2010年高雄地震震源参数的近远震波形联合反演

本文改进了传统基于近震波形数据的点源震源参数反演的Cut And Paste(CAP)方法,实现了近震Pnl波、面波和远震P波、SH波的联合反演的CAPjoint算法.对2010年3月高雄地震,分别进行单独反演以及联合反演,获得各自的震源机制解及深度,其中联合反演所得的最佳双力偶机制解参数为,节面1:走向317°,倾角36°,滑移角52°,节面2:走向181°,倾角62°,滑移角114°,深度为21 km.并对不同震中距波形对本次地震以及几种典型机制解断层几何参数的敏感性进行测试.为验证联合反演方法的可靠性,本文采用重抽样思想发展而来的Bootstrap方法,对近震数据的子集及其与远震数据的联合反演所得的参数进行统计,验证了在稀疏近台条件下联合反演中添加远震数据对地震震源参数约束的作用.     

华北地区基于噪声的瑞利面波群速度层析成像

本文利用华北地震科学台阵2007年1~4月份、190个宽频带和10个甚宽带垂直分量的地震噪声数据,通过互相关方法提取了瑞利面波的经验格林函数,用多重滤波方法测量了瑞利面波的群速度频散曲线.我们将研究区域划分为0.5°×0.5°的网格,利用噪声层析成像方法得到了研究区域7 s、12 s、16 s、23 s的瑞利面波群速度分布图像,所得结果较好地揭示了地壳内部、尤其是浅部地壳的横向速度变化.研究表明,短周期的群速度分布同地表地质结构、地形密切相关;华北地区的地壳结构具有明显的横向不均匀性,华北盆地及山间的沉积盆地显示出低速异常,而基岩广泛出露的太行山和燕山隆起区,呈现高速异常;多数强震(M≥6.0级)都发生在高群速度与低群速度的过渡地带.     

Surface wave tomography on a large-scale seismic array combining ambient noise and teleseismic earthquake data

We discuss two array-based tomography methods,ambient noise tomography (ANT) and two-planewave earthquake tomography (TPWT),which are capable of taking advantage of emerging large-scale broadband seismic arrays to generate high resolution phase velocity maps,but in complementary period band:ANT at 8-40 s and TPWT at 25-100 s period.Combining these two methods generates surface wave dispersion maps from 8 to 100 s periods,which can be used to construct a 3D v S model from the surface to ～200 km depth.As an i...     

基于接收函数与背景噪声联合反演的研究与应用

基于唐海—商都宽频带地震台阵2006—2009年连续三年的波形记录,利用环境噪声相位自相关函数对台阵下方的莫霍面反射P波进行分析。通过对同一个台站多个时间段的自相关结果进行分组、采用两步叠加处理增加信号强度：① 在组内进行线性叠加,对组间的叠加结果进行相位加权叠加;② 基于华北地区的背景速度结构信息,在地壳平均速度5%不确定性的时窗内,根据自相关函数包络线的二阶导数最大值确定P波的莫霍面反射时间,经时间−深度转化,获得台阵下方的莫霍面深度。结果显示,莫霍面从东南向西北总体由浅变深,中间有小幅度的起伏。噪声自相关方法确定的莫霍面平均深度相较于参考的接收函数结果的偏差为0.8 km,相应的双程走时偏差约为0.3 s。以月份叠加的自相关函数结果显示,PmP信号的噪声源具有显著的季节性变化。自相关函数的波形特征显示华北地区的地壳−地幔转换带的速度梯度模式不同。

     

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

基于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.     

Crust and upper mantle structure of the Bohemian Massif from the dispersion of seismic surface waves

Summary The phase velocity dispersion of Rayleigh waves for the Moxa-Vienna (MOX-VIE) and Moxa-Kaperské Hory (MOX-KHC) profiles, and of both Rayleigh and Love waves for the Kaperské Hory-Ksi (KHC-KSP) profile have been measured and inverted into models of shearwave velocity vs. depth. The three paths cross, respectively, the central part of the Bohemian Massif, its western margin, and the Bohemian Pluton and Cretaceous. For the MOX-VIE profile mean and lower crustal shear wave velocities of 3.7 and 3.9 km/s, respectively, a mean Moho depth of 34 km, and no existence of a low-velocity layer in the lower crust were found. The model obtained for the MOX-KHC profile is characterized by a slightly lower velocity in the lower crust (3.8 km/s), by a slightly lower Moho depth (32 km), and by the appearance of a weak low-velocity channel between 55 and 140 km. The crustal section of the final model for the KHC-KSP profile agrees well with the KHKS82 model derived by Novotný from results of DSS along international profile VII. Our final Rayleigh-wave model has significantly lower shear-wave velocities down to 215 km in the mantle. A systematic difference of 0.18 km/s between the average velocities of Rayleigh and Love waves has been revealed for the depth range from 30 to 215 km. Since almost no contamination of the fundamental Love mode with higher modes has been observed, and since the investigated structure hardly contains an unresolved system of thin, alternately low- and high-velocity layers, the cause of the difference is evidently polarization anisotropy of the upper mantle beneath the Bohemian Massif. It is recommended that the discussed investigations should be supplemented with data from the fan of KSP-GRF (Gräfenberg Array, Germany) paths and from the KHC-BRG (Berggiesshübel, Germany) profile.     

大别-苏鲁及其邻近地区基于背景 噪声的勒夫波群速度成像

大别-苏鲁造山带是中国大陆东部地区最重要的构造之一. 为了研究该地区的地壳上地幔速度结构,本文收集了国家数字地震台网和中国区域地震台网的山东、河南、安徽、江苏和湖北等省的144个宽频带地震台连续两年(2009年5月—2011年5月)的水平向地震记录(E分量和N分量)数据,首先对台站对之间E-E,E-N,N-N和N-E分量进行互相关,然后分别对这4个互相关分量采用相位权重叠加法进行叠加,最后旋转到横向分量(T-T)获得勒夫波经验格林函数(EGF);用频时分析(FTAN)方法获得4000余条勒夫波群速度频散曲线,并进一步反演得到了周期为6—40s的勒夫波群速度分布图.结果表明,周期为6—10s的勒夫波群速度分布与地表构造特征相吻合.大别造山带、苏鲁造山带、湖北西部隆起均表现为高速;华北盆地发育,表现为大面积的低速;江汉盆地、南襄盆地、合肥盆地等因其规模不同而显示不同程度的低速.在周期为6—30s的勒夫波群速度分布图上,大别和苏鲁地区均显示高速,已有的研究结果中地壳的低速并没有得到反映.其原因一方面可能与勒夫波群速度纵向分辨不高有关,另一方面高压变质岩深度分布可能比已有研究结果给出的要深. 6—30 s的分布图上郯庐断裂带及其邻近地区表现为不同程度的高速,可能与该地区白垩纪以来处于拉张构造体制,地幔物质受到扰动,造成物质上涌有关.     

广东省东部地区背景噪声面波群速度层析成像

选取广东省东部及福建和江西测震台网的26个固定台站2011年3～9月的地震背景噪声的垂直分量数据,用互相关方法获得瑞利面波的格林函数,再用频时分析方法提取瑞利面波的群速度频散曲线,最后用面波层析成像方法反演得到研究区周期为5～15s的瑞利面波群速度图像.反演结果揭示了广东省东部地区浅部地壳群速度结构存在明显横向不均匀性.短周期群速度异常与研究区内山脉、盆地的分布有较好的对应关系,沉积厚度较大的盆地地区表现为低速异常,而基底埋深较浅的山脉则显示了高速异常;地热分布对面波群速度产生重要影响,地热值较高的区域面波群速度一般表现为低速异常;粤东地区特别是粤东沿海地区地壳中部可能存在低速层.     

基于背景噪声研究大别-苏鲁及其邻区的瑞雷波群速度结构

大别-苏鲁造山带是中国大陆东部地区最重要的构造带之一. 为了研究该地区上中下地壳的速度结构, 选用国家数字地震台网和中国区域地震数据台网5省(山东、 安徽、 江苏、 河南和湖北)连续两年(2009年5月—2011年5月)的垂直向地震记录, 进行背景噪声互相关处理, 叠加得到了台站对间的面波经验格林函数. 采用多重滤波法提取了近4000条频散曲线, 并反演得到了研究区10—25 s的瑞雷波群速度分布结果. 通过分析大别-苏鲁及其邻区的瑞雷波群速度结构图像, 发现不同构造块体具有不同的瑞雷波群速度结构: ① 研究区内的郯庐断裂带及其周边地区包括鲁西地块和胶北地块上中下地壳均表现出明显的高速异常, 可能是在拉张环境下岩浆岩上涌, 导致高速、 高密度的变质岩在地壳富集而形成; ② 苏鲁高压变质带的瑞雷波速度在10—25 s周期内明显高于其它地区, 其上中下地壳均表现出较高的群速度结构特征, 认为苏鲁高压变质带至少延伸到下地壳, 而大别造山带在10 s时表现出高速特征, 但在15—25 s没有明显的高速特征, 故无法从其结果中判断大别高压变质带的垂向延伸范围; ③ 华北板块上中下地壳均表现为低速特征, 体现了研究区内华北板块的大陆地壳减薄特征.      

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.     

Non-uniqueness in surface-wave inversion and consequences on seismic site response analyses

Surface-wave tests are based on the solution of an inverse problem for shear-wave velocity profile identification from the experimentally measured dispersion curve. The main criticisms for these testing methodologies are related to the inverse problem solution and arise from the possible equivalence of different shear-wave velocity profiles. In this paper, some implications of solution non-uniqueness for seismic response studies are investigated using both numerical simulations and experimental data. A Monte Carlo approach for the inversion problem has been used to obtain a set of equivalent shear-wave velocity models. This selection is based on a statistical test which takes into account both data uncertainty and model parameterization. This set of solutions (i.e., soil profiles) is then used to evaluate the seismic response with a conventional one-dimensional analysis. It is shown that equivalent profiles with respect to surface-wave testing are equivalent also with respect to site amplification, thus countering the criticism related to inversion uncertainty for the engineering use of surface-wave tests.     

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

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

Source parameters inversion of the 2013 Lushan earthquake by combining teleseismic waveforms and local seismograms

On April 20, 2013, a magnitude M s 7.0 earthquake occurred in Lushan, Sichuan Province, China, and caused heavy casualties and economic losses. Based on the local broadband waveforms in Sichuan and adjacent provinces regional networks and teleseismic broadband records from IRIS stations, the focal mechanism and the focal depth are determined by the CAP (Cut And Paste) and its upgraded methods, CAPtele and CAPjoint, respectively. The results show that the focal mechanisms and depth from different methods are steady, and the best double couple solution derived from the joint inversion is 210°, 44°, and 91° for strike, dip, and rake angles respectively for one nodal plane and 29°, 46°, and 89° for another with 16 km focal depth and M w 6.66 moment magnitude. In order to verify the reliability of the results, a number of tests are performed based on local seismograms with different velocity models. They indicate that there is about 10 degree’s fluctuation in focal mechanisms and about 2 km variation in focal depth with a complex velocity structure. Furthermore, inverted by re-sampling the teleseismic waveforms on the basis of epicentral distance, the solutions are consistent with each other, which manifests that the teleseismic records are effective for constraining source parameters of the Lushan earthquake.     

Theoretical background for the inversion of seismic waveforms including elasticity and attenuation

To account for elastic and attenuating effects in the elastic wave equation, the stress-strain relationship can be defined through a general, anisotropic, causal relaxation function ^ijkl (x, ). Then, the wave equation operator is not necessarily symmetric (self-adjoint), but the reciprocity property is still satisfied. The representation theorem contains a term proportional to the history of strain. The dual problem consists of solving the wave equation withfinal time conditions and an anti-causal relaxation function. The problem of interpretation of seismic waveforms can be set as the nonlinear inverse problem of estimating the matter density (x) and all the functions ^ijkl (x, ). This inverse problem can be solved using iterative gradient methods, each iteration consisting of the propagation of the actual source in the current medium, with causal attenuation, the propagation of the residuals—acting as if they were sources—backwards in time, with anti-causal attenuation, and the correlation of the two wavefields thus obtained.