地球物理学报  2018, Vol. 61 Issue (10): 3951-3963   PDF    
利用接收函数方法确定青藏高原东北缘近地表S波速度
钱银苹1,2, 沈旭章1,3, 李翠芹4, 梅秀苹1, 焦煜媛1     
1. 中国地震局兰州地震研究所, 兰州 730000;
2. 广东省地震局, 广州 510070;
3. 中山大学地球科学与工程学院, 广州 510275;
4. 山东省地震局, 济南 250014
摘要:公式推导和数值实验都表明远震接收函数直达波幅度对近地表S波速度比较敏感,可以用于确定近地表(如沉积层)的S波速度值.本文利用该方法分析了青藏高原东北缘区域数字地震台网36个台站的远震资料,得到了研究区近地表S波速度分布特征.研究结果表明:在深度5±2.5 km范围内,祁连山—河西走廊区域近地表S波速度均值较小;位于地形梯度带的甘东南区域近地表S波速度呈现较为复杂的分布.通过与该区域沉积层厚度结果对比,发现近地表S波速度与沉积层厚度呈现较好的负相关.本研究是对目前接收函数方法的有效补充,在探测研究近地表速度结构方面,相比人工地震方法其成本几乎可忽略且环保.在目前地震波观测数据海量增长的背景下,该方法对于确定不同区域近地表结构具有较大潜力.
关键词: 接收函数      青藏高原东北缘      近地表S波速度     
Constraining the sub-surface S-wave velocity of the northeastern margin of Tibetan Plateau with receiver functions
QIAN YinPing1,2, SHEN XuZhang1,3, LI CuiQin4, MEI XiuPing1, JIAO YuYuan1     
1. Lanzhou Institute of Seismology, China Earthquake Administration, Lanzhou 730000, China;
2. Earthquake Administration of Guangdong Province, Guangzhou 510070, China;
3. School of Earth Science and Engineering, SunYat-sen University, Guangzhou 510275, China;
4. Earthquake Administration of Shandong Province, Jinan 250014, China
Abstract: We systematically investigated the effect of sub-surface S-wave velocity on the direct wave (Pp) of receiver function by the formula derivation and numerical tests. The results indicated that the sub-surface S-wave velocity can be determined with the receiver function amplitude. This method has been applied to the teleseismic waveforms of 36 stations in the northeastern margin of Tibetan Plateau, and obtained the sub-surface S-wave velocity distribution within the depth range of 5±2.5 km. The results show that the sub-surface S-wave velocity is smaller in the Qilianshan-Hexi corridor region, and is more complicated in the southeastern Gansu region of terrain gradient zone. The comparison of the sedimentary thickness and the sub-surface S-wave velocity implies that the sub-surface S-wave velocity is negatively correlated with the sedimentary thickness. This method is an effective supplement to the application of receiver function. Comparing with the traditional active-source seismic exploration, the cost of our method is nearly ignorable, as well environmentally friendly. So this method also has great potential with the exponential increasing of seismological observations in the future.
Keywords: Receiver function    The northeastern margin of Tibetan Plateau    Sub-surface S-wave velocity    
0 引言

近地表速度结构特征的研究在石油天然气勘探、有色金属勘查和防震减灾中都具有重要作用(法林等,1996滕吉文,2006王伟君等,2012张维等, 2012, 2013).目前,对近地表结构的探测主要通过人工震源(如爆炸源,震源车等)或有限深度的钻探获取(王有学等,2002何正勤等,2013邓小娟等,2017).由于相关研究方法成熟,且钻探方法更是可以直接取得深部岩石样本的信息,因此得到的结果精度也较高.但是,这些手段的成本都比较昂贵、对环境有一定的破坏性,而且在人口密集地区开展此类探测存在局限.在一些对结果要求精度不高的工作中,如果能利用天然地震波资料实现近地表结构的研究将会节约大量的成本.

天然地震波在传播过程中穿透了地球内部,因此通过分析天然地震波形资料可以对地球深部结构进行研究.接收函数是地震学中比较常用的利用天然地震波资料确定接收台站下方速度结构的方法(Yuan et al., 1997Niu and James, 2002Shen et al., 2011Pan et al., 2011韩明等,2017),该方法被广泛应用于地壳、岩石圈和上地幔间断面等的探测.如Zhu和Kanamori(2000)提出的确定地壳厚度和波速比的H-k搜索叠加方法;Yuan等(2006)利用理论地震图论证了S波接收函数研究岩石圈软流圈分界面的可靠性和局限性;钱银苹和沈旭章(2017)发展了利用接收函数幅度来约束间断面速度和密度跃变的δβ-δρ扫描叠加方法.这些研究对了解地球深部结构提供了重要佐证.

此外,前人研究结果表明,在近地表(10km以内深度)结构的探测中,接收函数也是一种比较有效的方法(Aki and Richards, 1980Julià et al., 2000;Julià,2007;彭恒初等,2012).Ammon(1991)对接收函数的震相幅度和近地表速度结构的关系进行了研究,表明直达P波垂向与径向的振幅比是关于近地表S波速度的函数;Julià(2007)对此开展了进一步的研究,利用接收函数的震相幅度计算了西班牙中部和印度地盾地区PAB和HYB两个台近地表S波速度,但是对S波速度对应的深度范围没有具体分析.

青藏高原东北缘是地学界研究前沿性重大科学问题的热点地区之一,该区域覆盖了很厚的沉积层能够放大地震的破坏,深入研究其近地表结构在防震减灾和工程建设具有重要意义(邓起东等,2002李永华等2006周民都等,2006Wang et al., 2014).同时,获取可靠稳定的近地表速度结构对地震定位、地壳上地幔深部结构成像结果分析、地震场地效应和改善深部结构研究的近地表约束也可提供有效支撑(何正勤等, 2007Zhang et al., 2013Ni et al., 2014张宝龙等,2016Li et al., 2018).此外,利用接收函数方法研究青藏高原东北缘近地表沉积层相关的结构,对于理解青藏高原隆升过程以及与之相伴的浅表风化剥蚀与物质运移过程非常重要.

本文通过数值实验,对Julià(2007)提出的利用接收函数幅度来约束近地表速度结构的方法进行了改进,并对该方法的稳定性和可靠性进行了测试.本文将该方法应用到青藏高原东北缘36个固定台站的远震资料,确定了该区域近地表S波速度,并结合研究区沉积层厚度对结果进行了讨论,特别是基于不同模型对近地表S波速度对应的深度进行了分析.

1 方法

由于远震事件到台站的距离足够远,地震波传播到台站下方时,以很陡的角度入射到台站(图 1a),可近似为平面波(图 1b).地震波到达接收台站会形成直达P波(Pp),也会在界面转换成S波(Ps),还会在地表反射遇到下方的间断面后再次反射到达地表就会形成多次波,在多次反射中也会产生P波和S波之间的转换以及PpPs和PpSs+ PsPs震相.图 1c为以爆炸源为震源的理论地震图,通过水平和垂直分量的反褶积可得到反映接收台站下方速度结构的接收函数时间序列(图 1d).下面结合接收函数计算过程,对震相幅度进行深入分析.

图 1 计算接收函数的流程图 (a)远震波形在地球介质中传播的示意图;(b)波形在介质中传播的射线路径;(c)以爆炸源为模型的理论地震图;(d)简单速度模型计算的接收函数. Fig. 1 The flow chart of calculating receiver function (a) The diagram of the propagation of the wave in the earth; (b) The ray path that the waveform propagates in the medium; (c) The theoretical seismogram with model of explosion source; (d) The receiver function for simple velocity model calculation.
1.1 原理

P波从水平分层底部入射,在穿过一维速度模型时会产生地表位移响应(Langston,1979图 1b),介质结构位移响应在时间域可表示为

(1)

(2)

其中Z(t)和R(t)分别表示介质结构响应垂直分量和径向分量的位移,s(t)表示入射平面波的有效震源时间函数,tk为震相到时(当k=0时为直达P波),zkrk分别表示震相垂直分量和径向分量的幅度,zkrk是与反射、透射系数乘积相关的数值.用地表位移的垂直分量对其水平分量作反褶积处理就可得到介质结构响应,在频率域中表示为

(3)

其中S(ω)为震源谱,H(ω)为介质结构响应,也就是径向接收函数(图 1d),ω为射线频率.对接收函数作傅里叶变换:

(4)

(5)

其中表示在垂直分量上通过直达P波脉冲幅度归一后的k震相幅度;.

接收函数能够反映地震波的属性,通过计算径向接收函数可以反推间断面结构.考虑水平单层的简单模型,仅存在直达P波、P-S转换波和多次波,(4)式和(5)式可简化为

(6)

(7)

把(6)式和(7)式代入(3)式,忽略高阶项后可得到

(8)

上式(8)中,当在自由表面时(t=0)则表示直达P波的绝对振幅.r0/z0的振幅比表示近地表的速度,在接收函数的建模中加入该幅度比信息非常重要.假设存在一个近地表的速度常数,我们可以得出入射P波引起的自由表面位移方程(Aki and Richards, 1980,p.190):

(9)

其中ηVS0=(1/VS02-2p2)1/2为S波的垂直慢度,p是水平慢度(即射线参数),VS0是近地表的S波速度.为了更好地解释近地表速度结构,(9)式表示在低频时震相振幅与平均速度之间的关系.对于更高频,振幅必须使用合成地震记录计算,但直达P波幅度比仍依赖于近地表速度(Ammon,1991).

直达P波振幅仅与近地表S波速度(VS0)和射线参数(p)有关,其中射线参数可由震中距和震源深度计算得到(Crotwell et al., 1999).当给定震中距时,直达P波幅度主要受近地表S波速度的影响,根据P波接收函数的幅度值即可反推近地表S波速度.由于接收函数只和台站下方速度结构有关,对于地壳横向结构较均匀的台站,不同震中距接收函数的波形资料是相近的、存在一定规律,因此可利用同一个台站不同震中距接收函数的直达P波幅度来约束台站下方近地表S波速度.

1.2 数值试验

上述过程基于公式推导给出了Pp震相幅度对VS0的影响,在实际工作中,我们通过对地震波形的反褶积和滤波来得到观测接收函数.为了与实际资料的处理过程具有可比性,我们利用理论地震图(Ma,2013)和处理实际资料相同的手段对上述的方法进行了测试.

首先建立单层地壳速度模型(图 2a),单层地壳S波速度值为3.75 km·s-1,波速比为1.732,地壳之下速度结构参考IASP91模型(Kennett et al., 1995).选取震中距为70°,利用反射透射率法计算理论地震图(Ma,2013).然后选取1.5的高斯滤波因子(alpha)通过时间域迭代反褶积方法计算理论接收函数(Ligorria and Ammon, 1999),在理论接收函数中测量Pp波幅度.我们以P波初动为起始时间,在起始时间附近(±2 s)提取其波峰的最大值,即为接收函数Pp波幅度的观测值.

图 2 单个震中距的接收函数(a)及对应的近地表S波速度(b) APp表示直达P波幅度,横坐标中“+”符号表示搜索到的最优的VS0. Fig. 2 The receiver function of single epicentral distance (a) and the sub-surface S-wave velocity (b) APp is the direct wave (Pp) amplitude of receiver function, and the "+" symbol is the optimal VS0 to be searched in the abscissa.

为了利用Pp波幅度确定近地表速度值,我们采取简单、直接的网格搜索法.选取VS0的变化范围为2.5~3.5 km·s-1,搜索步长为0.01 km·s-1,每搜索一次,可根据公式(9)计算Pp波绝对幅度的理论值.在接收函数时间轴(横轴)零时刻附近±2 s提取其幅度的最大值,即为接收函数Pp波幅度的观测值.

对于每一步搜索过程中的VS0都可算出不同震中距Pp幅度的理论值,其观测值和理论值之间的差值用适配函数misfit(VS0)表示,可作为衡量VS0是否接近实际值的标准.misfit(VS0)定义为

(10)

其中Aobs为Pp幅度的观测值,Asyn为Pp幅度理论值,Σ表示对不同射线参数的网格扫描叠加.根据公式(10)计算misfit(VS0),其最小值所对应横坐标为VS0的最优值(图 2b).如图 2a所示单个震中距的R向接收函数(RRF)测量其Pp波幅度,从而得到VS0.

对于理论接收函数,由于无任何干扰,利用单条波形就可以确定第一层的S波速度.实际观测中,会利用到不同震中距的资料,因此我们以不同震中距的理论接收函数为例,对数据处理过程进行表述.基于IASP91模型,选取震中距范围为45°~95°,步长为1°,用反射率法计算理论地震图(Ma,2013),并利用时间域迭代反褶积方法计算接收函数(Ligorria and Ammon, 1999).由于不同alpha对应的接收函数幅度存在差异,进一步验证不同alpha的影响,选取alpha为1.5和5时计算两组接收函数,以讨论和检测不同滤波频率窗对接收函数幅度的影响.对于每个震中距的每条接收函数,都可以确定一个VS0,利用所有接收函数确定的VS0,其均方误差(ε)的计算为

(11)

其中N为每个台站不同震中距的波形条数,VS0是所有VS0(i)的平均值,为剩余均方差,U为不同VS0的回归平方和,误差百分比δ=ε/VS0·100%.

图 3是不同震中距的理论接收函数及根据上述步骤得到的VS0.图 3a图 3c分别是alpha为1.5和5时的接收函数,图 3b图 3d分别为对应的VS0.当alpha=1.5时,VS0=3.79±0.1 km·s-1;当alpha=5时,VS0=3.78±0.1 km·s-1,初始模型中VS0=3.75 km·s-1,计算结果与初始模型基本一致,说明该方法对于研究近地表S波速度是有效的.

图 3 基于IASP91模型的理论接收函数及近地表S波速度 (a—b)分别是alpha=1.5时不同震中距的接收函数和VS0;(c—d)分别是alpha=5时不同震中距的接收函数和VS0,细线表示不同震中距的接收函数对应的VS0搜索过程,横坐标中“+”符号表示搜索到的最优的VS0. Fig. 3 Synthetic receiver functions based on IASP91 model (left) and the sub-surface S-wave velocity (right) The receiver functions (a) and VS0 (b) of different epicentral distances when alpha is 1.5;(c—d) The receiver functions (c) and VS0 (d) of different epicentral distances when alpha is 5, and the thin line indicates the VS0 search process, and the "+" symbol is the optimal VS0 to be searched in the abscissa.
2 观测资料处理

本研究收集了青藏高原东北缘36个宽频带三分量固定台站记录的远震资料,选取了时间为2009—2011年、震中距范围在40°~95°、MS>5.0的地震记录.由于P波接收函数上的信息是入射P波发生折射转换而形成,只有P波初动清晰的波形才能在穿过间断面时转化成可识别的震相,因此对远震资料进行预筛选.具体数据处理流程如下:

(1) 选取Z分量上P波初动清晰的资料,截取时间窗P波初动前20 s和后150 s的资料.

(2) 选取alpha为1.5对信号进行滤波,使用时间域反褶积方法(Ligorria and Ammon, 1999)得到11711条径向观测接收函数.

(3) 用观测资料的垂向分量对本身作反褶积,然后用其最大幅度值对径向接收函数作归一后得到绝对幅度(Ammon et al., 1990).

(4) 测量每条接收函数最大波峰(Pp)的幅度,采用和数值实验中相同的方法和步骤,计算得到研究区台站下方的近地表S波速度.

为了提高信噪比,将每个台站所有接收函数按照震中距在每1°范围内进行叠加(图 4).图 4a为GTA台最终所选的462条观测接收函数在不同震中距的叠加结果,其波形信噪比较高,各震相清晰;图 4b显示了GTA台根据上述方法得到的近地表S波速度VS0=3.21±0.27 km·s-1,百分误差为δ=8.6%;图 4c图 4d为WXT台不同震中距叠加后的401条观测接收函数和该台下方VS0=3.08±0.23 km·s-1δ=7.1%.

图 4 GTA和WXT台观测接收函数(左)及近地表S波速度(右) (a)和(c)分别GTA和WXT台不同震中距的接收函数;(b)和(d)为对应的VS0,细线表示不同震中距的接收函数对应的VS0结果,横坐标中“+”符号表示搜索到的最优的VS0. Fig. 4 The receiver functions of GTA and WXT stations (left) as well as their sub-surface S-wave velocity (rigth) The receiver functions of different epicentral distances to the GTA (a) and WXT stations (c); The sub_surface S-wave velocity of the GTA (b) and WXT stations (d), and the thin line indicates the VS0 search process, and the "+" symbol is the optimal VS0 to be searched in the abscissa.

相同的步骤应用于各个台站,所有台站的VS0结果及相关信息列在表 1.从表 1可看出,研究区域内VS0变化范围为2.59 km·s-1到3.50 km·s-1,均值为3.11 km·s-1,显示出0.92km·s-1的差异.此外,各台站VS0的误差小于±0.29 km·s-1,误差百分比都小于±10%,说明该方法对观测资料处理所得接收函数收敛性很好,近地表S波速度是可靠且稳定的.

表 1 各台站所得近地表S波速度(VS0)结果 Table 1 The results of the sub-surface S-wave velocity(VS0)

根据表 1中每个台站的结果,我们拟绘制了研究区域内VS0分布图(图 5),根据VS0分布特征并结合前人研究结果把研究区分为如图 5所示三部分.

图 5 台站分布和近地表S波速度 Ⅰ区域为99°E以西的塔里木盆地和祁连地块的交汇区域,Ⅱ区域为99°E—103°E的祁连—河西走廊区域,Ⅲ区域为103°E以东的甘东南区域. Fig. 5 The distribution of stations and the sub-surface S-wave velocity in the region The area Ⅰ is intersection of the Tarim Basin and the Qilian block west of 99°E; The Ⅱ area is the Qilian-Hexi corridor region of 99°E—103°E, and the Ⅲ area is the southeastern Gansu region of terrain gradient zone east of 103°E.

(1)在99°E以西的塔里木盆地和祁连地块的交汇区域VS0出现小幅度差异(图 5中Ⅰ区域),VS0从2.89 km·s-1变化到3.39 km·s-1,均值为3.1±0.25 km·s-1,相对较小.

(2) 祁连—河西走廊区域(图 5中Ⅱ区域,99°E—103°E)的VS0较稳定,均维持在3.18±0.26 km·s-1,与山脉走向也一致.

(3) 与祁连—河西走廊区域相比,在103°E以东的甘东南区域的VS0较为复杂(图 5中Ⅲ区域),最小为2.68 km·s-1,最大达到3.51 km·s-1,差异十分明显.导致VS0变化较大的可能原因在于该区域构造活动特殊,处于青藏块体、鄂尔多斯块体和扬子块体这3个块体碰撞作用的交接区域,成为青藏高原东部边缘最大的地形梯度带(Pan and Niu, 2011沈旭章等,2013Shen et al., 2017),因而相关的介质结构稳定性差,会在局部区域出现S波速度变化较大的结果(Ⅲ区域).

此外,Ⅱ区域VS0在祁连—河西走廊区域具有很好的稳定性和一致性,与Ⅲ区域相比其均值较小,而Ⅲ区域不仅均值较大且相邻台站间存在显著差异.

3 讨论

Wang等(2017)利用大量流动台阵资料的远震接收函数确定了青藏高原东北缘地区沉积层厚度,我们将其与本文地表S波速度结果进行了对比,如图 6所示.结果显示,祁连—河西走廊区域沉积层厚度相对周边较厚,而甘东南区域沉积层厚度变化较为复杂,局部呈现较大的起伏.该沉积层厚度结果和我们所得到的VS0结果成负相关:当沉积层较厚时,VS0较小;反之当沉积层较薄时,VS0较大,这也与人工地震得到的沉积层厚度结果相一致(滕吉文等,2010Wang et al., 2014).肖卓和高原(2017)用双差层析成像方法对青藏高原东北缘地壳速度结构进行反演发现近地表的地震波速度与地形和沉积层厚度对应较好.李敏娟等(2018)利用高密度台阵资料和双差层析成像技术反演了甘肃省东南部中上地壳速度结构,其深度在5 km处的S波速度约为3.5 km·s-1,本文对应区域S波速度为3.4 km·s-1,两者结果较为接近,因此本文结果可能反映了5 km深度左右的速度信息.对此我们进行进一步讨论和分析.

图 6 研究区域内沉积层厚度(Wang et al., 2017)和近地表S波速度(VS0)比较 背景图为Wang等(2017)沉积层厚度结果,图下方为沉积层厚度色标,图右方为S波速度色标,正方形为台站下方VS0,Ⅱ域、Ⅱ区域、Ⅲ区域同图 5. Fig. 6 Comparison of the sedimentary thickness (Wang et al., 2017) and VS0 in the study region The background map is Wang et al. (2017) sedimentary thickness results; At the bottom of the graph is a deep color mark; The square is VS0, and the right side is the S-wave velocity color mark. The Ⅰ area, the Ⅱ region and the Ⅲ region are the same as Fig. 5.

在计算VS0的方法中,公式(9)不含深度信息,因此我们无法直接确定其具体深度,为了测试S波速度的深度分布范围,我们利用两个模型的数值实验进行了分析和讨论(图 7图 8).图 7a模型1中在35 km间断面以上到地表速度呈梯度递增,图 8a模型2间断面深度50 km,地表以下0~25 km速度呈线性增加,在深度为25~50 km范围速度不变.图 7b图 8b是根据模型计算的两个高斯滤波因子的理论接收函数;图 7c图 8c为本文方法得到的两个高斯滤波因子对应的近地表S波速度,当高斯滤波因子变大(波长变短)时,对应的S波速度变小,所得到的近地表S波速度的深度逐渐变浅.把计算得到的S波速度与模型中平均S波速度比较,当计算的S波速度和模型某深度的平均速度接近,则可得到不同滤波因子的近地表S波速度所反映的深度(如图 7d图 8d).本文计算接收函数时选取alpha=1.5,对应深度约为7 km,据此估算本文得到的近地表S波速度对应深度约为5±2.5 km.

图 7 数值模型1的接收函数及alpha和反映深度的关系 (a)地表 0~35 km梯度变化的速度模型;(b) alpha为1.5和6.0的接收函数;(c) alpha为1.5和6.0时对应的S波速度;(d) alpha和反映深度的关系. Fig. 7 Relationship between the Gaussian filter with half-width (alpha) of calculate receiver function and the S-wave velocity corresponding depth for velocity model_1 (a) The velocity model of 0~35 km gradient change on the surface; (b) The receiver function of alpha is 1.5 and 6.0, and the S-wave velocity (c); (d) The relation between the alpha and the depth corresponding to the S-velocity structure.
图 8 数值模型2的接收函数及滤波因子和深度的关系 (a)地表 0~25 km梯度变化的速度模型; (b) alpha为1.5和6.0的接收函数; (c) alpha为1.5和6.0时对应的VS; (d)滤波因子和反映深度的关系. Fig. 8 Relationship between thealpha of calculate receiver function and the S-velocity corresponding depth for velocity model_2 Thevelocity model of 0~25 km gradient change on the surface (a); other illustrations are the same as Fig. 7.

彭恒初等(2012)提出通过对径向接收函数和垂向接收函数作低通滤波得到了S波视速度和低通滤波参数的变化曲线,然后转换成台站下方不同深度的S波速度结构.该方法中显示对于高频接收函数只有直达P波对径向接收函数(t=0)有贡献,可得到近地表的绝对S波速度.这与本文的分析一致,两种方法在一定程度上是互补的.

地震波速度的深度信息本质上是一种不同周期的体波对深部速度的反映,该方法可能具有和面波频散相同的特性,在未来深部结构的研究中具有较大的潜力.由于本文的近地表S速度结果和沉积层厚度(Wang et al., 2017)呈反相关,据此推断该结果可能主要反映了沉积层的S波速度.本文中,基于不同模型数值实验,对所得S波速度的深度范围进行了分析,但是尚未得出准确的深度信息.未来工作中我们将进一步对S波速度的深度信息进行深入研究,期望通过不同频率和大量模型结合能达到和面波频散类似的效果.

4 结论

本文对接收函数直达波幅度确定近地表S波速度的方法进行了细致分析,并通过数值实验对该方法进行了发展.基于青藏高原东北缘区域36个台站的远震接收函数,得到了该区域在5±2.5 km深度处收敛性较好的S波速度.

本文得到如下具体结论:

(1) 通过理论推导、数值实验及与前人研究结果的对比,本文提出的利用接收函数直达波绝对幅度研究近地表S波速度的方法具有较好的可靠性和稳定性.

(2) 本文得到的青藏高原东北缘近地表S波速度和沉积层厚度成负相关,说明本文结果可能主要反映了沉积层的S波速度.

(3) 研究区域的近地表S波速度结果显示,在研究区99°E以西的塔里木盆地和祁连地块的交汇区域VS0出现小幅度差异,VS0从2.89 km·s-1变化到3.39 km·s-1;在祁连—河西走廊区域VS0维持在3.18±0.26 km·s-1,与山脉走向也一致;在103°E以东的甘东南区域VS0最小为2.68 km·s-1,最大达到3.51 km·s-1,相对其他两个区域其结果横向差异较大.

(4) 在近地表结构的研究中,利用天然地震的接收函数方法相对人工震源探测成本可以忽略不计,而且绝对环保.随着天然地震波观测的大规模展开,本文提出的方法在对分辨率要求不太高区域近地表速度的研究中将发挥巨大作用.

致谢  感谢审稿专家宝贵的修改意见和建议.地震波形观测数据由甘肃省地震台网中心提供,中国地震台网中心的马延路博士提供了计算理论地震图和理论接收函数的程序,资料分析中使用了Lawrence Livermore实验室的SAC软件,在此一并表示感谢.
References
Aki K, Richards P G. 1980. Quantitative Seismology:Theory and Methods, Vol. 1 and 2. San Francisco:W. H. Freeman.
Ammon C J, Randall G E, Zandt G. 1990. On the nonuniqueness of receiver function inversions. Journal of Geophysical Research, 95(B10): 15303-15318. DOI:10.1029/JB095iB10p15303
Ammon C J. 1991. The isolation of receiver effects from teleseismic P waveforms. Bulletin of the Seismological Society of America, 81(6): 2504-2510.
Crotwell H P, Owens T J, Ritsema J. 1999. The TauP Toolkit:flexible seismic travel-time and ray-path utilities. Seismological Research Letters, 70(2): 154-160. DOI:10.1785/gssrl.70.2.154
Deng Q D, Zhang P Z, Ran Y K, et al. 2003. Basic characteristics of active tectonics of China. Science in China Series D:Earth Sciences, 46(4): 356-372.
Deng X J, Feng S Y, Zuo Y, et al. 2017. Velocity structure of the shallow part of the profile using surface wave of reflection seismic data. Journal of Geodesy and Geodynamics , 37(3): 256-260.
Fa L, Lin F, Chen W H. 1996. Source function derivation and property analysis of ceramic ring transducer for petroleum exploration. Acta Geophysica Sinica , 39(S1): 387-399.
Han M, Li J Y, Xu X Y, et al. 2017. Analysis for crustal anisotropy beneath the southeastern margin of Tibet by stacking azimuthal receiver functions. Chinese Journal of Geophysics , 60(12): 4537-4556. DOI:10.6038/cjg20171202
He Z Q, Ding Z F, Jia H, et al. 2007. To determine the velocity structure of shallow crust with surface wave information in microtremors. Chinese Journal of Geophysics , 50(2): 492-498.
He Z Q, Hu G, Lu L Y, et al. 2013. The shallow velocity structure for the Tonghai basin in Yunnan. Chinese Journal of Geophysics , 56(11): 3819-3827. DOI:10.6038/cjg20131123
Julià J, Ammon C J, Herrmann R B, et al. 2000. Joint inversion of receiver function and surface wave dispersion observations. Geophysical Journal International, 143(1): 99-112. DOI:10.1046/j.1365-246x.2000.00217.x
Julià J. 2007. Constraining velocity and density contrasts across the crust-mantle boundary with receiver function amplitudes. Geophysical Journal International, 171(1): 286-301. DOI:10.1111/gji.2007.171.issue-1
Kennett B L N, Engdahl E R, Buland R. 1995. Constraints on seismic velocities in the Earth from traveltimes. Geophysical Journal International, 122(1): 108-124. DOI:10.1111/gji.1995.122.issue-1
Langston C A. 1979. Structure under Mount Rainier, Washington, inferred form teleseismic body waves. Journal of Geophysical Research, 84(B9): 4749-4762. DOI:10.1029/JB084iB09p04749
Li M J, Shen X Z, Zhang Y S, et al. 2018. Fine crustal structures of northeast margin of the Tibetan Plateau and structural features of Jiuzhaigou earthquake focal area constrained by the data from a high-density seismic array. Chinese Journal of Geophysics , 61(5): 2075-2087. DOI:10.6038/cjg2018L0720
Li Y H, Wu Q J, An Z H, et al. 2006. The Poisson ratio and crustal structure across the NE Tibetan Plateau determined from receiver functions. Chinese Journal of Geophysics , 49(5): 1359-1368.
Li Y L, Wang B S, He R Z, et al. 2018. Fine relocation, mechanism, and tectonic indications of middle-small earthquakes in the Central Tibetan Plateau. Earth Planet. Phys., 2(5): 406-419. DOI:10.26464/epp2018038
Ligorria J P, Ammon C J. 1999. Iterative deconvolution and receiver-function estimation. Bulletin of the Seismological Society of America, 89(5): 1395-1400.
Ma Y L. 2013. YASEIS:Yet Another computer program to calculate synthetic SEIS mograms for a spherically multi-layered Earth model. //EGU General Assembly Conference Abstracts. Vienna, 15: 5596.
Ni S D, Li Z W, Somerville P. 2014. Estimating subsurface shear velocity with radial to vertical ratio of local P waves. Seismological Research Letters, 85(1): 82-90. DOI:10.1785/0220130128
Niu F L, James D E. 2002. Fine structure of the lowermost crust beneath the Kaapvaal craton and its implications for crustal formation and evolution. Earth and Planetary Science Letters, 200(1-2): 121-130. DOI:10.1016/S0012-821X(02)00584-8
Pan S Z, Niu F L. 2011. Large contrasts in crustal structure and composition between the Ordos Plateau and the NE Tibetan Plateau from receiver function analysis. Earth and Planetary Science Letters, 303(3-4): 291-298. DOI:10.1016/j.epsl.2011.01.007
Peng H C, Hu J F, Yang H Y, et al. 2012. An effective technique to constrain the non-uniqueness of receiver function inversion. Chinese Journal of Geophysics , 55(4): 1168-1178. DOI:10.6038/j.issn.0001-5733.2012.04.013
Qian Y P, Shen X Z. 2017. The approach and application of constraining velocity and density contrasts across the Moho using receiver functions. Chinese Journal of Geophysics , 60(8): 2980-2992. DOI:10.6038/cjg20170807
Shen X Z, Mei X P, Zhang Y S. 2011. The crustal and upper-mantle structures beneath the Northeastern Margin of Tibet. Bulletin of the Seismological Society of America, 101(6): 2782-2795. DOI:10.1785/0120100112
Shen X Z, Zhou Y Z, Zhang Y S, et al. 2013. Geodynamic significance of the crust structure beneath the northeastern margin of Tibet. Progress in Geophys. , 28(5): 2273-2282. DOI:10.6038/pg20130509
Shen X Z, Liu M, Gao Y, et al. 2017. Lithospheric structure across the northeastern margin of the Tibetan Plateau:Implications for the plateau's lateral growth. Earth and Planetary Science Letters, 459: 80-92. DOI:10.1016/j.epsl.2017.11.027
Teng J W. 2006. Strengthening geophysical exploration and exploitation of metallic minerals in the second deep space of the crustal interior. Geological Bulletin of China , 25(7): 767-771.
Teng J W, Wang F Y, Zhao W Z, et al. 2010. Velocity structure of layered block and deep dynamic process in the lithosphere beneath the Yinshan orogenic belt and Ordos Basin. Chinese Journal of Geophysics , 53(1): 67-85. DOI:10.3969/j.issn.0001-5733.2010.01.008
Wang S J, Wang F Y, Zhang J S, et al. 2014. The P-wave velocity structure of the lithosphere of the North China craton-results from the Wendeng-Alxa Left Banner deep seismic sounding profile. Science China:Earth Sciences, 57(9): 2053-2063. DOI:10.1007/s11430-014-4903-7
Wang W J, Chen Q F, Chen L. 2012. Review on near-surface soft structure investigation with seismic noise. Earthquake Research in China , 28(4): 335-350.
Wang W L, Wu J P, Fang L H, et al. 2017. Sedimentary and crustal thicknesses and Poisson's ratios for the NE Tibetan Plateau and its adjacent regions based on dense seismic arrays. Earth and Planetary Science Letters, 462: 76-85. DOI:10.1016/j.epsl.2016.12.040
Wang Y X, Yu Q F, Han G H. 2002. The fine velocity structure in the shallow crust along the HQ-E seismic profile in Southeastern China. Geophysical & Geochemical Exploration , 26(2): 91-96.
Xiao Z, Gao Y. 2017. Crustal velocity structure beneath the northeastern Tibetan plateau and adjacent regions derived from double difference tomography. Chinese Journal of Geophysics , 60(6): 2213-2225. DOI:10.6038/cjg20170615
Yuan X H, Ni J, Kind R, et al. 1997. Lithospheric and upper mantle structure of southern Tibet from a seismological passive source experiment. Journal of Geophysical Research:Solid Earth, 102(B12): 27491-27500. DOI:10.1029/97JB02379
Yuan X H, Kind R, Li X Q, et al. 2006. The S receiver functions:synthetics and data example. Geophysical Journal International, 165(2): 555-564. DOI:10.1111/gji.2006.165.issue-2
Zhang B L, Li Z W, Bao F, et al. 2016. Shallow shear-wave velocity structures under the Weishan volcanic cone in Wudalianchi volcano field by microtremor survey. Chinese Journal of Geophysics , 59(10): 3662-3673. DOI:10.6038/cjg20161013
Zhang W, He Z Q, Hu G, et al. 2012. Detect the velocity structure of shallow crust with artificial and nature source Rayleigh wave technology. Technology for Earthquake Disaster Prevention , 7(1): 26-36.
Zhang W, He Z Q, Hu G, et al. 2013. Detection of the shallow velocity structure with surface wave prospection method. Progress in Geophys. , 28(4): 2199-2206. DOI:10.6038/pg20130467
Zhang Y Q, Guo Y R, Hou W, et al. 2013. Geochemical characteristics and exploration potential of the middle-upper Ordovician source rocks on the Western and Southern Margin of Ordos Basin. Natural Gas Geoscience, 24(5): 894-904.
Zhou M D, Zhang Y S, Shi Y L, et al. 2006. Three-dimensional crustal velocity structure in the northeastern margin of the Qinghai-Tibetan plateau. Progress in Geophys. , 21(1): 127-134.
Zhu L P, Kanamori H. 2000. Moho depth variation in southern California from teleseismic receiver functions. Journal of Geophysical Research:Solid Earth, 105(B2): 2969-2980. DOI:10.1029/1999JB900322
邓起东, 张培震, 冉勇康, 等. 2002. 中国活动构造基本特征. 中国科学:地球科学, 32(12): 1020-1030.
邓小娟, 酆少英, 左莹, 等. 2017. 利用反射地震资料中的面波研究剖面浅部结构. 大地测量与地球动力学, 37(3): 256-260.
法林, 林峰, 陈文辉. 1996. 石油勘探用薄圆环压电换能器的瞬态特性分析. 地球物理学报, 39(S1): 387-399.
韩明, 李建有, 徐晓雅, 等. 2017. 按方位叠加接收函数分析青藏高原东南缘的地壳各向异性. 地球物理学报, 60(12): 4537-4556. DOI:10.6038/cjg20171202
何正勤, 丁志峰, 贾辉, 等. 2007. 用微动中的面波信息探测地壳浅部的速度结构. 地球物理学报, 50(2): 492-498. DOI:10.3321/j.issn:0001-5733.2007.02.021
何正勤, 胡刚, 鲁来玉, 等. 2013. 云南通海盆地的浅层速度结构. 地球物理学报, 56(11): 3819-3827. DOI:10.6038/cjg20131123
李敏娟, 沈旭章, 张元生, 等. 2018. 基于密集台阵的青藏高原东北缘地壳精细结构及九寨沟地震震源区结构特征分析. 地球物理学报, 61(5): 2075-2087. DOI:10.6038/cjg2018L0720
李永华, 吴庆举, 安张辉, 等. 2006. 青藏高原东北缘地壳S波速度结构与泊松比及其意义. 地球物理学报, 49(5): 1359-1368. DOI:10.3321/j.issn:0001-5733.2006.05.015
彭恒初, 胡家富, 杨海燕, 等. 2012. 接收函数反演地壳S波速度结构的有效约束方法. 地球物理学报, 55(4): 1168-1178. DOI:10.6038/j.issn.0001-5733.2012.04.013
钱银苹, 沈旭章. 2017. 接收函数确定Moho面速度和密度跃变的方法研究及应用. 地球物理学报, 60(8): 2980-2992. DOI:10.6038/cjg20170807
沈旭章, 周元泽, 张元生, 等. 2013. 青藏高原东北缘地壳结构变化的地球动力学意义. 地球物理学进展, 28(5): 2273-2282. DOI:10.6038/pg20130509
滕吉文. 2006. 强化开展地壳内部第二深度空间金属矿产资源地球物理找矿、勘探和开发. 地质通报, 25(7): 767-771. DOI:10.3969/j.issn.1671-2552.2006.07.001
滕吉文, 王夫运, 赵文智, 等. 2010. 阴山造山带-鄂尔多斯盆地岩石圈层、块速度结构与深层动力过程. 地球物理学报, 53(1): 67-85. DOI:10.3969/j.issn.0001-5733.2010.01.008
王伟君, 陈棋福, 陈凌. 2012. 松散浅层结构的地震噪声探测综述. 中国地震, 28(4): 335-350. DOI:10.3969/j.issn.1001-4683.2012.04.001
王有学, 余钦范, 韩国华. 2002. HQ-E爆破地震测深剖面的地壳浅部精细结构及其地质构造研究. 物探与化探, 26(2): 91-96. DOI:10.3969/j.issn.1000-8918.2002.02.002
肖卓, 高原. 2017. 利用双差成像方法反演青藏高原东北缘及其邻区地壳速度结构. 地球物理学报, 60(6): 2213-2225. DOI:10.6038/cjg20170615
张宝龙, 李志伟, 包丰, 等. 2016. 基于微动方法研究五大连池火山区尾山火山锥浅层剪切波速度结构. 地球物理学报, 59(10): 3662-3673. DOI:10.6038/cjg20161013
张维, 何正勤, 胡刚, 等. 2012. 用人工源和天然源面波联合探测浅层速度结构. 震灾防御技术, 7(1): 26-36. DOI:10.3969/j.issn.1673-5722.2012.01.003
张维, 何正勤, 胡刚, 等. 2013. 用面波联合勘探技术探测浅部速度结构. 地球物理学进展, 28(4): 2199-2206. DOI:10.6038/pg20130467
周民都, 张元生, 石雅鏐, 等. 2006. 青藏高原东北缘地壳三维速度结构. 地球物理学进展, 21(1): 127-134. DOI:10.3969/j.issn.1004-2903.2006.01.019