接收函数方法的研究综述

从远震体波波形数据中提取的P波接收函数和S波接收函数已经成为研究台站下方地壳上地幔速度间断面的最有效的方法之一.P波接收函数已经被广泛应用于获取地壳内部S波速度结构、地壳厚度及物质成分组成、地幔过渡带的厚度变化以及岩石圈地幔的间断面等,而S波接收函数是P波接收函数的一个很好的补充,因为在Moho和地幔过渡带之间的深度范围内,Sp转换波比来自浅部间断面的多次波到达早,而在P波接收函数中,相同深度范围的间断面的Ps转换波往往被来自浅部间断面的多次波干扰或者淹没,因此S波接收函数是目前获取岩石圈地幔深度范围内速度间断面结构(如Moho和LAB)的比较有效的方法,比面波观测具有更高的分辨率.本文详细阐述了P波接收函数和S波接收函数的方法原理以及在地壳上地幔速度间断面的研究中所采用的研究思路.     

接收函数方法与华北克拉通北部岩石圈结构的研究

正自从Langston提出接收函数以来,随着地震观测台网越来越密集,利用地震台阵观测研究构造带下方速度结构取得了重要进展,很多学者利用接收函数研究地壳上地幔速度结构。虽然P波接收函数可以获得地壳上地幔结构,但由于莫霍面和壳内间断面多次反射震相的干扰,单纯考虑Ps转换震相难以精确地确定岩石圈边界。Farra和Vinnik利用类似提取P波接收函数的方法得到Sp转换震相的S波接收函数。相比P波接收函数,由于S波接收函数不受间断面多次反射震相的干扰,因而在岩石圈—软     

长白山火山区壳幔S波速度结构研究

利用面波层析成像和远震接收函数方法对长白山地区的地壳上地幔速度结构进行了研究。结果表明:长白山火山区附近存在岩石圈减薄、上地幔软流圈增厚以及上地幔S波速度降低等与上地幔高温物质有关的现象,它表明长白山的岩浆系统一直延伸到上地幔软流圈范围。天池火山区地壳内部存在明显的S波低速层,在离天池火山口较近的WQD台附近,低速层顶部埋深约8km,厚度近20km,S波最小速度约2.2km/s。在距离天池火山北部50km的EDO台地壳中没有明显的低速层。火山区S波速度结构总体表现出距离天池越近,地壳的V_P/V_S越大,低速层的厚度和幅度增加的特征,表明天池火山口附近地壳内部存在高温物质或岩浆囊。CBS台站不同方位的接收函数及反演结果表明,地表低速层厚度以及莫霍面深度存在随方位的变化。地表低速层在南部方向明显较厚,莫霍面深度在南部天池火山口方向存在小幅度抬升。CBS台站附近特殊的近地表速度结构可能是该台站记录的火山地震波形主频较低的主要因素。天池火山口附近莫霍面的小幅度抬升意味着存在与火山作用有关的壳幔物质交换通道     

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

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

大陆深俯冲的深浅动力学响应:来自帕米尔高原地壳精细结构的约束

帕米尔高原位于青藏高原西构造结,强烈的陆内中源地震活动指示帕米尔高原下方正发生着大陆深俯冲过程.文章基于帕米尔高原现有流动地震台阵观测资料,利用接收函数谐波分析和环境噪声方法,在有效协调接收函数与面波频散分辨尺度的情况下,通过接收函数和面波频散联合反演,重建了横跨帕米尔高原到费尔干纳盆地等主要构造单元的壳幔二维S波速度结构.通过与接收函数CCP叠加成像剖面之间的对比分析,可靠地揭示了帕米尔高原中源地震分布、Moho面形态、上地幔和壳内低速异常分区/分层之间的空间配置关系,不仅进一步证实了亚洲大陆下地壳在帕米尔高原下方与岩石圈地幔耦合在一起发生了深俯冲,而且指示了深俯冲地壳物质的变质脱水反应在该区中源地震活动和壳内变形中起到的重要作用,从而为认识和理解大陆深俯冲过程及其动力学响应提供了新的地震学制约.     

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

本文回顾了接收函数线性反演及接收函数与频散曲线联合反演的理论.为了检验两种方法在实测数据中的应用效果,我们选择云南区域台网位于四个不同次级块体的四个台站为例:首先,计算了四个台站的远震P波接收函数;然后,分别运用接收函数线性反演和接收函数与背景噪声数据联合反演两种方法分别得到台站下方的S波速度结构;另外,运用H-k扫描方法验证两种反演方法获取的莫霍面深度的可靠性.在反演过程中,当我们提供具有深一些莫霍面的初始模型时,联合反演获取的莫霍面深度与H-k扫描结果具有很好的一致性,同时联合反演也能有效地捕捉低速层,对反演的S波速度也有较好的约束.由此推知,加入背景噪声相速度数据后的联合反演对初始模型的依赖性低于单独用接收函数反演的结果.     

利用接收函数计算地壳各向异性的可靠性分析及倾斜界面的影响*

使用接收函数研究壳幔速度间断面和速度结构已是常用的技术,但介质各向异性或倾斜的莫霍(Moho)界面,都会造成接收函数波形的复杂性。本文利用远震P波接收函数计算台站下方的地壳各向异性,通过信噪比测试和谐波分析两种方法来验证各向异性结果的可靠性。通过甘肃地震台网的两个台站记录,讨论各向异性和倾斜界面对接收函数的影响,结果显示,台站BYT(白银)下方具有各向异性,而台站WYT(渭源)下方由于可能存在倾斜界面,得到的各向异性结果则有待进一步分析。为了更好地认识倾斜界面对各向异性计算结果的影响,采用合成理论地震图,计算接收函数,然后利用合成接收函数进行各向异性分析。结果表明,计算得到的快波方向不会受到倾斜界面的影响,但是时间延迟会受到影响。     

接收函数曲波变换去噪与偏移成像

增强接收函数偏移图像的垂向分辨率意味着提高参与叠加的接收函数的频率,但是采用高频接收函数通常伴随着对接收函数质量和参考速度模型的更高要求.通过叠加处理可去除部分接收函数中的随机噪声干扰,但同一台站的接收函数之间经常存在难以通过简单叠加消除的噪声信号.压制接收函数随机噪声的干扰可加强成像效果和提高图像分辨率,对推进叠加偏移成像质量的提高有重要的实际意义.本文利用在川西地区布设的31个流动台站所记录的远震波形数据,使用曲波变换去噪后信噪比增强的接收函数进行共转换点叠加(CCP),获得沿北纬31°线下方800km深度范围内速度间断面图像.研究结果表明:(1)对接收函数进行曲波变换去噪,可压制随机噪声,增强转换震相的追踪性,提高数据信噪比;(2)通过去噪处理,大幅提高接收函数用于偏移成像的主频率;(3)偏移结果确认了接收函数反演得到的松潘和川滇块体下方具有厚度约10～20km的过渡性Moho的认识;(4)上地幔过渡带的结果预示在龙门山断裂带以西的小范围内有可能存在下地壳或上地幔物质的拆沉.     

利用接收函数计算地壳各向异性的可靠性分析及倾斜界面的影响

使用接收函数研究壳幔速度间断面和速度结构已是常用的技术,但介质各向异性或倾斜的莫霍(Moho)界面,都会造成接收函数波形的复杂性。本文利用远震P波接收函数计算台站下方的地壳各向异性,通过信噪比测试和谐波分析两种方法来验证各向异性结果的可靠性。通过甘肃地震台网的两个台站记录,讨论各向异性和倾斜界面对接收函数的影响,结果显示,台站BYT(白银)下方具有各向异性,而台站WYT(渭源)下方由于可能存在倾斜界面,得到的各向异性结果则有待进一步分析。为了更好地认识倾斜界面对各向异性计算结果的影响,采用合成理论地震图,计算接收函数,然后利用合成接收函数进行各向异性分析。结果表明,计算得到的快波方向不会受到倾斜界面的影响,但是时间延迟会受到影响。     

远震接收函数方法及其在海域岛礁区地壳结构研究中的应用

远震接收函数方法从远震P波波形中提取出台站下方主要速度间断面的信息,是研究地壳、上地幔结构的有效方法。该方法通过拟合转换波波形来约束间断面深度和横波速度。传统的接收函数线性反演方法强烈依赖于初始模型的选取,其反演结果存在较大的不确定性,为此新的基于非线性的反演方法逐渐发展起来,另外,接收函数与面波联合反演以及S波接收函数方法的提出也在一定程度上降低了反演结果的非唯一性。接收函数方法在大陆区地壳结构的研究中得到了广泛运用,但对于海区和岛礁区的研究却非常缺乏,本文回顾了国内外一些岛礁区的接收函数研究实例,介绍了西沙群岛琛航岛的天然地震观测及其结果,探讨了接收函数方法在岛礁区地壳结构研究中的应用前景。     

Crustal Thickness and Velocity Ratio beneath National Seismic Stations in the Sichuan-Yunnan Area Based on Receiver Function Analysis

By using the teleseismic receiver function method, this paper analyzes the crustal thickness and v_P/v_S ratios beneath the 4 National seismic stations (KMI, TNC, CD2 and PZH) in the Sichuan-Yunnan area. This study gives the variance of Moho depths and velocity ratios of the 4 stations in different directions. The results show that the Moho depth beneath the Kunming station is around 50km, and the velocity ratio varies between 1.62 and 1.69. The thickness of crust and the velocity ratio do not change much with the direction. The crust beneath Tengchong station shows clear directivity, being 40.7km thick in the northeast and 49.7km thick in the southeast. The difference of the v_P/v_S values is remarkable between the two directions, reaching 0.2. The Chengdu station also has shallow Moho, about 40km, but is 8km deeper in the northeast and southwest and the velocity ratio has a change of 0.13 between the two directions. The crust beneath the Panzhihua station is stable. In all directions, the Moho depth is around 60km and the v_P/v_S ratio doesn't change significantly.     

利用远震接收函数研究黄海周边地区的地壳上地幔结构

利用宽频带流动台站(YSBSN)记录的远震波形数据和远震接收函数方法,反演了黄海东、西两侧地壳上地幔的S波速度结构.结果表明,莫霍面深度在30~38 km之间变化,位于中方一侧的JNN台下方地壳厚度最大,可以归因于华北板块和扬子板块的碰撞;韩方一侧的地壳厚度自北向南逐渐变厚,但仍然难以厘定朝鲜半岛南部潜在碰撞带的位置,这些问题的解决需要更大范围的流动台站观测.由于部分台站位于巨厚的沉积层和多孔的火山岩之上,与浅部构造的相关性使得接收函数表现出较大振幅的混响,从而影响了来自深部结构的转换震相.     

接收函数确定Moho面速度和密度跃变的方法研究及应用

通过对单层模型反射和透射系数的推导,提出了利用接收函数一次转换波和多次波确定Moho面速度和密度跃变的速度-密度跃变(δβ-δρ)扫描叠加方法.利用反射率法计算了不同模型的远震理论地震图,按照与处理实际观测波形一致的方法和流程计算了理论接收函数;根据不同模型数值试验结果,深入分析了界面速度和密度跃变对接收函数震相幅度的影响.利用(δβ-δρ)扫描叠加方法,对理论接收函数进行了数值试验,结果证明了该方法的可行性.最后将该方法应用于位于青藏高原东北缘的高台(GTA)台和兰州(LZH)台,确定了两个台站下方Moho面的速度跃变分别约为(19±1)%和(20±1)%,密度跃变最小值为(4±2)%和(6±2)%.     

Crustal structure of Northwest Zagros (Kermanshah) and Central Iran (Yazd and Isfahan) using teleseismic Ps converted phases

Receiver functions are widely employed to detect P-to-S converted waves and are especially useful to image seismic discontinuities in the crust. In this study we used the P receiver function technique to investigate the velocity structure of the crust beneath the Northwest Zagros and Central Iran and map out the lateral variation of the Moho boundary within this area. Our dataset includes teleseismic data (M _b ≥ 5.5, epicentral distance from 30° to 95°) recorded at 12 three-component short-period stations of Kermanshah, Isfahan and Yazd telemetry seismic networks. Our results obtained from P receiver functions indicate clear Ps conversions at the Moho boundary. The Moho depths were firstly estimated from the delay time of the Moho converted phase relative to the direct P wave beneath each network. Then, we used the P receiver function inversion to find the properties of the Moho discontinuity such as depth and velocity contrast. Our results obtained from PRF are in good agreement with those obtained from the P receiver function modeling. We found an average Moho depth of about 42 km beneath the Northwest Zagros increasing toward the Sanandaj-Sirjan Metamorphic Zone and reaches 51 km, where two crusts (Zagros and Central Iran) are assumed to be superposed. The Moho depth decreases toward the Urmieh-Dokhtar Cenozoic volcanic belt and reaches 43 km beneath this area. We found a relatively flat Moho beneath the Central Iran where, the average crustal thickness is about 42 km. Our P receiver function modeling revealed a shear wave velocity of 3.6 km/s in the crust of Northwest Zagros and Central Iran increasing to 4.5 km/s beneath the Moho boundary. The average shear wave velocity in the crust of UDMA as SSZ is 3.6 km/s, which reaches to 4.0 km/s while in SSZ increases to 4.3 km/s beneath the Moho.     

长白山-镜泊湖火山区地壳结构接收函数研究

利用71个远震的波形资料,用接收函数方法提取了布设在长白山—镜泊湖火山区的34个宽频带流动数字地震台站的接收函数,通过对接收函数反演,获得了台站下方的S波速度结构.研究结果表明,沈阳—敦化一线莫霍面深度32～33km,向西地壳厚度加厚,到长春附近地壳厚度约为36km.在天池火山口莫霍面深度为达38km,而镜泊湖火山口森林的莫霍面深度约为39km.总体看研究区的地壳厚度是南浅北深.长白山天池火山口附近地下10km左右有一明显的低速层存在;镜泊湖火山口森林附近30km也可能有低速体存在;研究发现莫霍面上S波速度梯度在火山口附近和远离火山口有明显区别.在火山口附近其莫霍面的S波速度梯度比非火山口地区的S波速度梯度明显小,说明火山口下与一般的地壳莫霍面结构有差别.研究发现沈阳—敦化一线两侧的莫霍面深度有较大变化,其位置与地表的敦化—密山断裂基本一致,说明敦化—密山断裂是研究区的一条非常重要的地质构造带.     

Lithospheric structure of NW Iran from P and S receiver functions

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b?>?5.5, epicentral distance between 30° and 95° for P receiver functions and M_b?>?5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.     

Erratum to: Lithospheric structure of NW Iran from P and S receiver functions

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b > 5.5, epicentral distance between 30° and 95° for P receiver functions and M_b > 5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.     

沿107.6°E南北向剖面鄂尔多斯地块及周缘地区地壳结构

鄂尔多斯地块紧邻青藏高原东北缘,位于华北克拉通的西部,在我国中生代、新生代以来东部地区的构造活动中起到了重要作用.对鄂尔多斯及其周缘地区的研究可以提供有关华北克拉通的形成、演化和破坏过程的重要信息.本文选取了纵贯鄂尔多斯的107.6°E附近南北剖面上的44个流动地震台站进行分析,采用接收函数方法,进行Kirchhoff偏移成像,并且结合在该区域内前人的地震面波频散进行联合反演,获得剖面下方的地壳内部精细结构.研究结果显示:(1)莫霍面在鄂尔多斯北部较平缓,约45km深;在鄂尔多斯南部有所加深,达到50km;其北边的河套盆地的地壳厚度约为50km;南边的渭河盆地到秦岭地区及四川盆地的地壳厚度从约为40km增厚到47～50km.(2)河套盆地下方存在大规模的低速异常,最深可达25km,反映了其显著的拉张构造和沉积历史.(3)秦岭造山带下方的低速异常对应于其主要为长英质的地壳组分,可能是由于中生代的拆沉作用导致的地壳下部基性岩石层的缺失.(4)以38°N为界的鄂尔多斯地块,南北部地壳速度结构存在差异,可能表明了这两部分经历的构造历史不同.     

南海西沙石岛地震台下的地壳结构研究

为加强海上天然地震观测和岩石圈结构研究,我们在南海西沙群岛建立了石岛流动地震台,进行了为期一年多的观测试验.观测结果表明,虽然岛礁区的地震数据受热带气旋影响期间背景噪声较大,但仍能记录到MW 6级以上的地震.本文对具有清晰P波波形的远震记录进行了接收函数处理和计算模拟,得到台站下简单的地壳结构模型,其莫霍面深度为28 km,上地壳顶部有一层2 km厚的低速层,横波速度只有23 km/s,向下逐渐过渡到横波速度为38 km/s的下地壳.与已有的研究结果比较,发现石岛台与琼中台的地壳结构模型是西沙海槽剖面的自然延伸,西沙地块的地壳结构属减薄型陆壳,可与华南地块的正常陆壳对比.     

Receiver functions of CCDSN and crustal structure of Chinese mainland

The teleseismic receiver functions of 48 stations belonging to the CCDSN are used to invert the crustal structure beneath each station with the neighborhood algorithm. Thin layers with low velocity have been found beneath eight stations with "abnormal" observed receiver functions. Unreasonable results of few stations have been adjusted lightly with the trial-and-error method. The final result indicates that the crust in the western China is relatively thicker than the eastern China. The crust thickness beneath the Tibetan plateau is very large, which reaches 84 km at the station LSA. Double-crust structure exists below the stations LSA and CAD in Tibet, which might imply the collision between the Indian and Eurasian plates. A pronounced low velocity zone in the lower crust beneath the station TNC of Yunnan province might relate to the high temperature or emergence of partially molten material caused by Quaternary volcano, magma and geothermal activities in this area. The Moho is a transitional zone made up of thin layers instead of simple sharp discontinuity beneath several stations. The Conrad discontinuity is clearly identified beneath 20 stations mainly in the southeastern China, whereas it is blurry beneath 14 stations and uncertain beneath remaining stations.