地球物理学报  2017, Vol. 60 Issue (6): 2304-2314   PDF    
青藏高原东北缘及周边地区基于程函方程的面波层析成像
钟世军1,2, 吴建平1 , 房立华1, 王未来1, 范莉苹1, 王怀富1     
1. 中国地震局地球物理研究所, 北京 100081;
2. 北京市地震局, 北京 100080
摘要: 利用"中国地震科学台阵探测"在南北地震带北段布设的670套宽频带地震台站记录到的面波资料,使用新近发展的程函方程面波层析成像方法,获得了青藏高原东北缘及周边地区12~60 s周期范围比以往成像结果具有更高分辨率的瑞利面波相速度分布图像.青藏高原东北缘的祁连褶皱系西段、秦岭褶皱系西段和松潘—甘孜褶皱系,在16~60 s周期范围内均显示出明显的低速异常分布,表明该地区的地壳力学强度较低,在强烈的构造应力作用下易发生形变.与西段不同,祁连褶皱系东段和秦岭褶皱系中段的相速度分布特征揭示,其中下地壳的速度明显高于高原内部区域.鄂尔多斯块体整体上表现为稳定块体具有的高速特征,但其西部边缘在中上地壳的速度比块体中部地区偏低,且存在一定的不均匀性.鄂尔多斯块体西北缘的临河断陷盆地和西缘的银川断陷盆地,在较短的周期范围内(12~20 s)表现为局部低速特征,但与银川断陷盆地不同,临河断陷盆地的低速特征可一直延续至60 s周期以上,表明该盆地下方地壳及上地幔速度明显偏低,可能与深部热作用有关.阿拉善块体与其北部地区的速度差异主要表现在中上地壳,这一现象值得今后进一步探讨.基于程函方程面波层析成像方法给出了青藏高原东北缘及周边地区高分辨率的成像结果,揭示了以往面波层析成像难以获得的深部细节特征,为该地区的深部构造研究提供了新的信息.
关键词: 青藏高原东北缘      面波层析成像      程函方程      瑞利波      相速度     
Surface wave Eikonal tomography in and around the northeastern margin of the Tibetan plateau
ZHONG Shi-Jun1,2, WU Jian-Ping1, FANG Li-Hua1, WANG Wei-Lai1, FAN Li-Ping1, WANG Huai-Fu1     
1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Beijing Earthquake Administration, Beijing 100080, China
Abstract: Based on surface wave data recorded by 670 dense deployed seismic stations of ChinArray, high-resolution Rayleigh wave phase velocity images between 12 s and 60 s are obtained using the recently developed Eikonal tomography method.The lateral resolution is higher than ever before.The results show that from periods of 16 s to 60 s, distinct low-velocity anomalies exist in the western and middle Qilian fold system and Songpan-Garzê fold system.It suggests low mechanical strength of the crust in these areas, which makes the crust easy to deform under strong tectonic stress.Different from western parts, phase velocities of middle-lower crust in the eastern Qilian fold system and eastern Qinling fold system are higher than that in the interior of the Tibetan plateau.The Ordos block behaves as a stable massif, featured by high-velocity anomalies.However, in upper and middle crust, local low-velocity anomalies at periods 12 s to 20 s are present in the Linhe rift basin and Yinchuan rift basin.Different from the Yinchuan rift basin, low-velocity anomalies in the Linhe rift basin extend to 60 s or more.It indicates low velocity in crust and upper mantle beneath the Linhe rift basin, which may be associated with the thermal process from deep mantle.The velocity difference between the Alxa block and its adjacent northern area mainly exists in upper-middle crust, which needs to be further discussed.The high-resolution images determined with Eikonal tomography in and around the northeastern margin of the Tibetan plateau reveal more detailed deep structure than previous surface wave tomography results and provide new information for research of deep structure in this region.
Key words: Northeastern margin of the Tibetan plateau      Surface wave tomography      Eikonal equation      Rayleigh wave      Phase velocity     
1 引言

青藏高原东北缘是青藏高原向内陆扩展的前缘部分, 中、晚新世以来该地区经历了强烈的造山运动(Tapponnier et al., 2001), 内部发育了多条规模较大的断裂, 如祁连山北缘断裂、海原断裂、西秦岭北缘断裂等(图 1).青藏高原东北部和鄂尔多斯地块之间有宽度为200 km的巨大边界带, 是中国最活跃的地震带之一(张培震等, 2003).GPS观测研究表明, 青藏高原东北缘地壳缩短强烈, 现今运动方向存在明显的横向变化(Zhang et al., 2004Gan et al., 2007王双绪等, 2013).人工地震剖面揭示青藏高原东北缘中下地壳存在高低速相间的低速结构(李松林等, 2002刘宝峰等, 2003张先康等, 2003Liu et al., 2006嘉世旭和张先康, 2008);接收函数研究结果(陈九辉等, 2005李永华等, 2006段永红等, 2007刘启民等, 2014姚志祥等, 2014张洪双等, 2015)揭示该地区地壳厚度变化范围超过20 km, 地壳内部速度结构与一些重要断裂带和地质块体之间存在明显的相关性.地震体波、面波和噪声层析成像从更宽的范围(郭飚等, 2004徐果明等, 2007Zheng et al., 2008, 2010Yang et al., 2010Sun et al., 2010周民都等, 2012黄忠贤等, 2013Li et al., 2013Bao et al., 2013Li et al., 2014徐小明等, 2015)揭示了青藏高原东北缘及周边地区速度结构强烈的横向变化, 但不同的结果之间存在较大的差异(徐果明等, 2007Yang et al., 2010Sun et al., 2010Zhang et al., 2011黄忠贤等, 2013Li et al., 2013, 2014Bao et al., 2013徐小明等, 2015).

图 1 青藏高原东北缘及周边地区构造背景分布图(图 1图 2图 4的断裂分布邓起东等(2002)) F1:狼山山前断裂;F2:噔口—本井断裂;F3:贺兰山东麓断裂;F4:黄河—灵武断裂;F5:祁连山北缘断裂;F6:海原断裂;F7:云雾山断裂;F8:六盘山断裂;F9:西秦岭北缘断裂;F10:迭部—白龙江断裂;F11:龙门山断裂带. Fig. 1 Tectonic setting in and around the northeastern margin of the Tibetan plateau F1:Langshan range-front fault; F2:Dengkou-Benjing fault; F3:Helanshan east piedmont fault; F4:Yellow River-Lingwu fault; F5:Qilianshan north edge fault; F6:Haiyuan fault; F7:Yunwushan fault; F8:Liupanshan fault; F9:West Qinling north edge fault; F10:Diebu-Bailongjiang fault; F11:Longmenshan fault zone.
图 2 青藏高原东北缘及周边地区流动台站分布图 蓝色三角形代表流动台站, 黑色圆圈代表右下角地名所在位置. Fig. 2 Map showing distribution of temporal seismic stations in and around the northeastern margin of the Tibetan plateau Blue triangles indicate temporary stations. Black circles represent the locations of places on the lower right.
图 4 青藏高原东北缘不同周期的瑞利面波相速度分布图 Fig. 4 Rayleigh wave phase velocities at different periods in the northeastern margin of the Tibetan plateau

尽管在青藏高原东北缘及周边地区已开展了一系列深部地球物理研究, 但现有的结果或受剖面位置的限制难以给出该地区整体的深部结构特征, 或者受台站分布的限制成像分辨率较低, 很难获得块体内部以及不同块体之间深部结构的变化细节特征, 制约了对该地区深部动力学过程的理解与认识.最近“中国地震科学台阵探测”项目在南北地震带北段地区布设的高密度流动台阵, 为青藏高原东北缘及周边地区高分辨率地震波成像提供了难得的机遇.本文利用该密集台阵记录的面波资料, 使用面波层析成像方法揭示该地区高分辨率的瑞利波相速度分布特征.

传统的面波成像方法均假定面波沿大圆路径传播, 而忽略了复杂介质引起传播路径弯曲对成像结果造成的影响.Lin等(2009)提出了基于程函方程的面波成像方法, 该方法通过重建面波不同周期的空间走时场, 应用程函方程计算各周期走时面的空间梯度, 从而得到对应周期空间慢度场的分布.它不但考虑了偏离大圆路径的传播效应, 而且避免了常规成像方法的正演和反演过程.基于程函方程的面波成像方法在台站比较均匀密集的区域可以获得与台站间距相当的分辨率(Lin et al., 2009Lin and Ritzwoller, 2011Gouédard et al., 2012Zhou et al., 2012), 取得很好的成像效果.Jin和Gaherty (2015)应用基于程函方程的面波成像方法, 获得了整个美国大陆高分辨率的相速度分布图像.本文利用青藏高原东北缘及周边地区布置的670套宽频带流动地震台站(图 2)资料(喜马拉雅地震科学台阵, 2011), 应用基于程函方程的面波成像方法, 得到了研究区域内12~60 s周期瑞利面波相速度分布图像, 分辨率可达30~50 km.高分辨率的相速度分布图像将会为研究该地区的壳幔结构、构造变形及地震活动性等研究提供重要约束.

2 数据和方法

由中国地震局地球物理研究所牵头, 共15家单位参加的地震行业科研专项“中国地震科学台阵探测——南北地震带北段”项目, 在青藏高原东北缘及周边地区布设了670套宽频带流动地震台站, 这是迄今为止在该地区布设台站数目最多、密度最高的地震观测台阵, 平均台间距为35 km.我们收集了2013年9月至2015年8月期间记录到的849个地震事件(图 3)的垂直向波形记录.在对数据正式处理前, 我们对地震事件的波形记录进行了预挑选, 震中距限制在10°~160°范围内、震源深度小于50 km、震级大于MS5.0.然后, 对挑选出的垂直向记录进行预处理, 包括去均值、去倾斜、去仪器响应和数据重采样(1 Hz).最终选择了442个地震事件(图 3)用于面波层析成像研究.

图 3 地震事件相对于台阵中心的方位分布图 黑色三角形代表台阵中心位置, 红色实心圆圈代表挑选出的地震事件的震中位置. Fig. 3 Distributions of azimuths of seismic events with respect to array center circles represent epicenters of selected seismic events. Black triangle represents the center of the array. Red solid

Gee和Jordan (1992)提出了基于广义地震数据泛函分析的方法, 该方法被广泛应用于分析构建地壳及上地幔模型.这些分析的共同点是将观测地震记录与由参考地球模型合成的理论地震记录做互相关计算, 然后应用五参数高斯小波拟合经过窄通带滤波后的相关波形, 提取出一系列频散参数.然后应用这一系列频散参数反演得到改善的参考地球模型(Gaherty and Jordan, 1995Gaherty et al., 1996Gaherty, 2001, 2004Gaherty and Dunn, 2007Chen et al., 2007).Jin和Gaherty (2015)将广义地震数据泛函分析方法应用到多道上, 不同的是应用邻近台站对实际观测波形的互相关函数提取频散参数.五参数高斯小波拟合表达式为:

(1)

其中Fi(ωi)为中心频率为ωi的零相位高斯窄带滤波器, 带宽随中心频率线性变化;C (t)为互相关函数, 应用汉宁窗WC隔离出互相关函数主要能量部分;Ga(x)为高斯函数e-x2/2;拟合波形的五个参数为Aσωtgtp, 分别代表幅度因子、窄带波形的半带宽、窄带波形的中心频率、窄带波形的群延迟时间及ωi频率对应的相延迟时间.

由程函方程的表达式可知(Wielandt, 1993Shearer, 1999), 单频面波空间走时面梯度场的倒数即为相速度的分布.Lin等(2009, 2011)应用最小曲率面拟合技术(Smith and Wessel, 1990)来构造走时场的拟合曲面, 而Jin和Gaherty (2015)并没有直接应用该方法得到相速度, 而是采用了慢度向量反演的方法.因为应用慢度向量反演的方法能够压制噪声引起的高波数、低幅度的相速度变化, 而且相对于前者得到的相速度分辨率更好.本文应用Jin和Gaherty (2015)改进后的基于程函方程的面波层析成像方法来研究青藏高原东北缘及周边地区的瑞利波相速度分布.

反演慢度向量的目标函数由两项组成:第一项反映观测得到的相延迟与预测值之间的失配值, 第二项反映对慢度向量的光滑约束.其表达式如下:

(2)

其中S(r)为慢度向量, ri是连接两个邻近台站的球面路径, δτpi为观测得到的两个邻近台站的相延迟时间, SR为慢度向量沿大圆路径的分量, ST为垂直于大圆路径的分量, λ为光滑因子.

与传统面波成像中采用双台法测量相速度不同, 程函方程面波成像方法利用台站周围所有可能的台站对进行面波传播速度测量与反演, 它不需要对面波传播方向进行假设, 不必做射线追踪, 可以自然地考虑射线可能存在的弯曲现象, 获得更可靠的相速度测量结果(Lin et al., 2009);程函方程面波成像方法对地震资料的要求与传统面波方法相似, 要求地震记录中具有明显面波信号, 震中距不宜过近(与震源深度有关).该方法的不足之处是对台站分布要求较高, 只有台站分布相对密集均匀的情况下才能获得较好的成像结果.

我们选择了Jin和Gaherty (2015)开发的面波相速度自动测量系统来计算瑞利面波相速度.整个处理流程分为三大步:基于(1) 式, 利用单个地震的台阵观测资料, 测量邻近台站对的相延迟时间, 得到观测区域内的相延迟时间分布;基于(2) 式, 利用单个地震事件测量得到的相延迟时间分布, 计算观测区域内的相速度分布;对区域内各网格节点, 由不同地震事件计算得到的相速度加权平均, 权重为慢度反演的射线密度, 最终得到综合各选定地震事件的相速度分布结果.

3 结果

应用基于程函方程的面波层析成像方法, 我们得到了12~60 s周期范围内的瑞利面波相速度分布(图 4).可以看出, 各周期不同构造单元相速度差异明显, 具有显著的横向非均匀性.

面波相速度频散对于纵波速度、横波速度及密度的敏感程度不同, 理论研究表明瑞利面波相速度对横波速度最为敏感, 而纵波和密度对瑞利面波的相速度影响较弱.某一周期的相速度大致对其波长三分之一深度处的S波速度最敏感, 因此不同周期的相速度分布图反映了不同深度范围内的S波速度结构变化.我们计算了基阶瑞利面波对横波速度的敏感核(图 5), 地球模型参考了CRUST2.0和AK135, 地壳分为三层, 厚度为50 km.由图 5b所示, 短周期面波相速度的敏感深度较浅、范围较窄, 中长周期的敏感深度较深、范围较宽.因此短周期面波探测深度较浅、在深度方向上的分辨率较高, 而中长周期面波探测深度较深、但在深度方向上的分辨率会逐渐降低.

图 5 (a)计算敏感核所用的参考模型;(b)基阶Rayleigh面波不同周期相速度在不同深度对S波速度的敏感程度分布 Fig. 5 (a)Reference model used in sensitivity kernel calculation; (b) Sensitivity kernels of fundamental Rayleigh wave phase velocity with respect to shear wave velocity at different periods

周期为12 s的瑞利波相速度对8~18 km深度的S波速度比较敏感, 主要反映了研究区域内上地壳的速度变化.临河断陷盆地和银川断陷盆地表现为显著的低速异常, 祁连褶皱系西段、秦岭褶皱系西段和松潘—甘孜褶皱系均呈现出微弱的低速异常分布.鄂尔多斯块体西部靠近边界附近存在中低速异常分布, 其南部靠近边界处存在小范围高速异常.秦岭褶皱系中段和扬子地块北部存在大范围的高速异常分布.阿拉善地块的中低速分布特征与其北部地区的高速异常形成了鲜明对比.

周期为16 s的瑞利波相速度对15~25 km深度的S波速度比较敏感, 主要反映了研究区域内上地壳下部和中地壳上部的速度变化.祁连褶皱系西段、秦岭褶皱系西段和松潘—甘孜褶皱系的低速异常比较明显, 并且范围也有所扩大.鄂尔多斯块体内部呈现大面积高速分布特征, 其西缘银川断陷盆地的低速异常明显减弱.

周期为20~25 s的瑞利波相速度对20~40 km深度的S波速度比较敏感, 主要反映了研究区域内中地壳的速度变化.祁连褶皱系西段、秦岭褶皱系西段和松潘—甘孜褶皱系的低速异常变得十分显著, 并且祁连褶皱系的低速异常呈现出沿该区断裂带走向延伸的分布特征.鄂尔多斯块体内部、阿拉善地块以北地区、秦岭褶皱系中段和扬子块体北部均显示出大范围的高速分布特征, 并且鄂尔多斯块体南部和秦岭褶皱系中段的高速异常连成一体.

周期为32~40 s的瑞利波相速度对30~70 km深度的S波速度比较敏感, 主要反映了研究区域内下地壳及上地幔顶部的速度变化.在祁连褶皱系西段、秦岭褶皱系西段和松潘—甘孜褶皱系, 低速异常范围均有所增大.阿拉善块体及其以北地区, 相速度横向变化明显减弱, 整体上呈现出比较均匀的中高速分布特征.鄂尔多斯块体西部边界附近相速度明显低于块体内侧, 其南部边界至秦岭褶皱系中段的高速异常逐渐减弱.

周期为50~60 s的瑞利波相速度对60~140 km深度的S波速度比较敏感, 主要反映了研究区域内上地幔浅部的速度变化.祁连褶皱系西段、秦岭褶皱系西段和松潘—甘孜褶皱系的低速异常几乎连成一体, 而且该低速异常向北东方向扩展到了海原断裂附近.鄂尔多斯块体南部边界至秦岭褶皱系中段已表现出微弱的低速异常分布.扬子块体内部的高速异常大致以其边缘的断裂带为界与松潘—甘孜褶皱系和秦岭褶皱系中段分隔开来.

4 讨论

作为青藏高原的重要边界带及强震活动带, 许多学者利用固定台站和流动地震台站的观测资料对青藏高原东北缘地区的深部结构进行了研究, 获得了一些很有意义的认识(何正勤等, 2002苏伟等, 2008Zheng et al., 2008, 2010Yang et al., 2010Sun et al., 2010Li et al., 2014徐小明等, 2015).但受台站分布的限制, 现有结果在揭示块体边界以及块体内部深部结构的非均匀变化特征方面仍然存在不足.本文给出了横跨整个青藏高原东北缘地区, 比以往结果具有更高横向分辨率的面波相速度成像结果, 展示了一些低分辨率成像结果难以获得的细节特征.成像结果揭示出地表地质构造与深部速度结构具有很好的相关性.在青藏高原东北缘, 整体表现出高原内部速度低, 边缘及外部速度高的变化趋势;在祁连褶皱系西段与东段、秦岭褶皱系西段与中段、阿拉善块体与其北部块体、鄂尔多斯与周边拉张盆地等均表现出明显的深部结构差异.

以祁连褶皱系、秦岭褶皱系西段为界, 其南部不同周期均表现出低速特征, 表明青藏高原内部地壳与上地幔顶部S波速度较低.基于固定地震台站的一些研究结果也揭示该地区地壳速度较低, 且存在中下地壳低速层(李永华等, 2006吴建平等, 2009Bao et al., 2013Li et al., 2014王椿镛等, 2015).地壳内部速度的高低, 与地壳介质力学强度密切相关, 低速异常通常代表介质的力学强度较低, 高速异常往往代表介质具有较大的力学强度.力学强度较低的地壳介质, 在应力作用下更容易发生变形.青藏高原东北部的快速隆升与力学强度较低的地壳上地幔介质在强烈的区域构造应力作用下易发生变形有关.我们的成像结果虽然与Li等(2013)的结果在东昆仑断裂以南地区具有相似性, 但在东昆仑和海原断裂之间40~60 s周期的面波相速度并不存在明显的局部高速异常, 12~32 s周期的面波相速度在该地区表现为低速区内小幅度增高的特征, 表明该地区地壳平均速度比青藏高原东北缘其他区域略高, 我们推测东昆仑与祁连山之间平缓而略低的地形特征可能与该地区地壳介质强度略大, 地壳变形相对缓慢有关.

在秦岭褶皱系中段, 12~40 s周期的相速度表现为明显的高相速度分布特征, 只是在50~60 s周期范围内有所减弱, 但速度值仍然明显高于秦岭褶皱系西段和松潘—甘孜褶皱系, 32 s以上周期的相速度比鄂尔多斯内部略低.Zhang等(2011)获得的成像结果也揭示该区域速度明显高于青藏高原内部地区, 但主要受台站分布和成像分辨率的限制, 其结果未能揭示出鄂尔多斯块体的高速特征, 以及秦岭褶皱系中段与鄂尔多斯块体之间的速度过渡变化特征.Clark和Royden (2000)通过拟合青藏高原东部边缘地形变化, 对青藏高原下地壳的柔性进行了估计, 并提出在青藏高原内部中下地壳存在黏滞系数很低的通道流, 围绕四川盆地向外流动.本文的结果显示, 四川盆地与鄂尔多斯之间的秦岭褶皱系及其以西的近200 km范围内, 地壳结构应以高速特征为主, 中下地壳不存在明显的低速层结构.因此, 如果青藏高原内部存在地壳通道流, 也很难从四川盆地与鄂尔多斯块体之间流过, 或者说通道流在这一地区并不明显.

鄂尔多斯块体西北缘的临河断陷盆地在12~60 s周期范围内均表现为显著的低速异常, 低速异常形态与断陷盆地的展布相一致.地质研究表明, 临河凹陷于渐新世开始形成, 在新生代沉积厚度达14.8 km, 沉积层厚度呈现北深南浅、西深东浅的分布特征(国家地震局《鄂尔多斯周缘活动断裂系》课题组, 1988).我们推测, 该区域较短周期(12~16 s)的低速异常是由巨厚沉积层引起的.临河断陷盆地存在长周期的显著低速异常, 则表明下地壳及上地幔顶部S波速度仍然较低, 可能与上地幔热作用有关.在12~25 s周期, 鄂尔多斯块体西缘的银川断陷盆地也显示出比较明显的低速异常分布, 但在其西侧小范围的弱低速异常在40 s周期仍有明显显示.银川断陷是新生代时期贺兰山解体而形成的次级构造, 为陆相断陷盆地.而断陷盆地内部结构比较复杂, 受到多条正断层的控制, 形成中央凹陷.银川断陷盆地较短周期的低速异常同样也与较厚的低速沉积层有关, 而中长周期的低速异常揭示该盆地下方中下地壳及上地幔顶部存在弱低速结构.

鄂尔多斯块体内部相速度显著高于青藏高原东北部, 整体上呈现出高速异常分布特征, 但其西部边界附近相速度低于块体内部, 且存在明显的横向不均匀特征, 表明稳定的鄂尔多斯块体在西部边缘地区可能已受到构造作用的影响, 产生了一定程度的构造变形.在12~40 s周期, 祁连山以北和阿拉善块体以北地区的相速度介于鄂尔多斯与青藏高原地区之间, 但更接近稳定的鄂尔多斯块体.并且在20 s周期后, 相速度横向变化较小, 表明中下地壳具有高速、横向均匀的特征, 是相对稳定的块体.而阿拉善地块, 在12~25 s周期, 相速度介于青藏高原与其北部块体之间, 内部存在一些局部低速异常, 表明块体可能正处在经历早期破坏阶段.地质研究表明, 阿拉善地块为中朝古板块的西延部分, 其南侧为正在隆升的青藏高原, 东侧为稳定的鄂尔多斯块体.喜马拉雅运动时期, 在青藏高原NNE方向推挤作用下, 阿拉善地块周边变形比较明显(张进等, 2007), 其西缘的龙首山地区、东缘的贺兰山地区、以及南缘的香山—天景山地区构造活动强烈, 在这些地区形成了新生代活动断裂带.阿拉善地块内部地震记载历史短暂, 可中强地震屡次发生, 最大地震为1954年的民勤7.0级地震.块体内部地震呈NEE向密集条带, 总体上形成了一个三角状的地震“填充”区, 这与稳定的鄂尔多斯块体内部缺震或少震现象形成了鲜明对比(王萍和王增光, 1997), 表明阿拉善地块内部正在经历复杂的构造变形.

5 结论

本文使用近年来新发展起来的程函方程面波层析成像方法, 采用密集宽频带地震台阵观测资料, 获得了青藏高原东北缘及周边地区高分辨率瑞利面波相速度分布图像, 揭示了以往面波层析成像难以获得的深部细节特征, 为进一步研究该地区的深部构造提供了新的重要信息.

青藏高原东北缘祁连褶皱系西段、秦岭褶皱系西段和松潘—甘孜褶皱系在16~60 s周期范围内均显示出明显的低速异常分布, 表明该地区地壳及上地幔顶部附近的S波速度较低, 介质的力学强度较弱, 与周边块体的速度结构特征形成鲜明的对比.力学强度较低的青藏高原壳幔介质在板块碰撞导致的强烈构造挤压应力作用下, 易发生变形、褶皱、增厚和断层活动.

鄂尔多斯块体内部相速度显著高于青藏高原东北缘, 整体上呈现出高速异常分布特征, 但其西部边界附近相速度低于块体内部, 且存在明显的横向不均匀特征, 表明稳定的鄂尔多斯块体在西部边缘地区可能已受到构造作用的影响, 产生了一定程度的构造变形.祁连山以北和阿拉善块体以北地区的相速度介于鄂尔多斯与青藏高原地区之间, 但更接近稳定的鄂尔多斯块体, 可以确认其中下地壳具有高速、横向变化相对较小的特征, 是相对稳定的块体.阿拉善地块中上地壳的相速度介于青藏高原与其北部块体之间, 可能是正在经历早期破坏的块体.

鄂尔多斯块体西北缘的临河断陷盆地和西缘的银川断陷盆地, 在较短的周期范围内(12~20 s)表现为局部低速特征, 但与银川断陷盆地不同, 临河断陷盆地的低速特征可一直延续至60 s周期以上, 表明该盆地下方地壳及上地幔速度明显偏低, 可能与深部热作用有关.

致谢

感谢“中国地震科学探测台阵数据中心”为本文提供观测资料.感谢两位匿名审稿专家及编辑部老师提出的宝贵意见.

参考文献
Bao X W, Song X D, Xu M J, et al. 2013. Crust and upper mantle structure of the North China Craton and the NE Tibetan Plateau and its tectonic implications. Earth Planet. Sci. Lett., 369-370: 129-137. DOI:10.1016/j.epsl.2013.03.015
Chen J H, Liu Q Y, Li S C, et al. 2005. Crust and upper mantle S-wave velocity structure across Northeastern Tibetan Plateau and Ordos block. Chinese J. Geophys., 48(2): 333-342. DOI:10.3321/j.issn:0001-5733.2005.02.015
Chen P, Zhao L, Jordan T H. 2007. Full 3D tomography for the crustal structure of the Los Angeles region. Bull. Seismol. Soc. Am., 97(4): 1094-1120. DOI:10.1785/0120060222
ChinArray-Himalaya. 2011. China Seismic Array waveform data of Himalaya Project (in Chinese). Institute of Geophysics, China Earthquake Administration.
Clark M K, Royden L H. 2000. Topographic ooze: Building the eastern margin of Tibet by lower crustal flow. Geology, 28(8): 703-706. DOI:10.1130/0091-7613(2000)28〈703:TOBTEM〉2.0.CO;2
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.
Duan Y H, Zhang X K, Liu Z, et al. 2007. Crustal structure using receiver function in the East part of A'Nyêmaqên suture belt. Acta Seismologica Sinica, 29(5): 483-491.
Gaherty J B, Jordan T H. 1995. Lehmann discontinuity as the base of an anisotropic layer beneath continents. Science, 268(5216): 1468-1471. DOI:10.1126/science.268.5216.1468
Gaherty J B, Jordan T H, Gee L S. 1996. Seismic structure of the upper mantle in a central Pacific corridor. J. Geophys. Res., 101(B10): 22291-22309. DOI:10.1029/96JB01882
Gaherty J B. 2001. Seismic evidence for hotspot-induced buoyant flow beneath the Reykjanes Ridge. Science, 293(5535): 1645-1647. DOI:10.1126/science.1061565
Gaherty J B. 2004. A surface wave analysis of seismic anisotropy beneath eastern North America. Geophys. J. Int., 158(3): 1053-1066. DOI:10.1111/j.1365-246X.2004.02371.x
Gaherty J B, Dunn R A. 2007. Evaluating hot spot-ridge interaction in the Atlantic from regional-scale seismic observations. Geochem. Geophys. Geosyst., 8(5): Q05006. DOI:10.1029/2006GC001533
Gan W J, Zhang P Z, Shen Z K, et al. 2007. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements. J. Geophys. Res., 112(B8): B08416. DOI:10.1029/2005JB004120
Gee L S, Jordan T H. 1992. Generalized seismological data functionals. Geophys. J. Int., 111(2): 363-390. DOI:10.1111/j.1365-246X.1992.tb00584.x
Gouédard P, Yao H J, Ernst F, et al. 2012. Surface wave eikonal tomography in heterogeneous media using exploration data. Geophys. J. Int., 191(2): 781-788. DOI:10.1111/j.1365-246X.2012.05652.x
Guo B, Liu Q Y, Chen J H, et al. 2004. Seismic tomographic imaging of the crust and upper mantle beneath the Northeastern edge of the Qinghai-Xizang plateau and the Ordos area. Chinese J. Geophys., 47(5): 790-797. DOI:10.3321/j.issn:0001-5733.2004.05.009
He Z Q, Ding Z F, Ye T L, et al. 2002. The group velocity distribution of Rayleigh waves and crustal and upper mantle velocity structure of the Chinese mainland and its vicinity. Acta Seismologica Sinica, 24(3): 252-259.
Huang Z X, Li H Y, Xu Y. 2013. Lithospheric S-wave velocity structure of the North-South seismic belt of China from surface wave tomography. Chinese J. Geophys., 56(4): 1121-1131. DOI:10.6038/cjg20130408
Jia S X, Zhang X K. 2008. Study on the crust phases of deep seismic sounding experiments and fine crust structures in the northeast margin of Tibetan plateau. Chinese J. Geophys., 51(5): 1431-1443. DOI:10.3321/j.issn:0001-5733.2008.05.016
Jin G, Gaherty J B. 2015. Surface wave phase-velocity tomography based on multichannel cross-correlation. Geophys. J. Int., 201(3): 1383-1398. DOI:10.1093/gji/ggv079
Li L, Li A B, Shen Y, et al. 2013. Shear wave structure in the northeastern Tibetan Plateau from Rayleigh wave tomography. J. Geophys. Res., 118(8): 4170-4183. DOI:10.1002/jgrb.50292
Li S L, Zhang X K, Zhang C K, et al. 2002. A preliminary study on the crustal velocity structure of Maqin-Lanzhou-Jingbian by means of deep seismic sounding profile. Chinese J. Geophys., 45(2): 210-217. DOI:10.3321/j.issn:0001-5733.2002.02.007
Li X F, Li H Y, Shen Y, et al. 2014. Crustal velocity structure of the northeastern Tibetan Plateau from ambient noise surface-wave tomography and its tectonic implications. Bull. Seismol. Soc. Am., 104(3): 1045-1055. DOI:10.1785/0120130019
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 J. Geophys., 49(5): 1359-1368. DOI:10.3321/j.issn:0001-5733.2006.05.015
Lin F C, Ritzwoller M H, Snieder R. 2009. Eikonal tomography: surface wave tomography by phase front tracking across a regional broad-band seismic array. Geophys. J. Int., 177(3): 1091-1110. DOI:10.1111/j.1365-246X.2009.04105.x
Lin F C, Ritzwoller M H. 2011. Helmholtz surface wave tomography for isotropic and azimuthally anisotropic structure. Geophys. J. Int., 186(3): 1104-1120. DOI:10.1111/j.1365-246X.2011.05070.x
Liu B F, Li S L, Zhang X K, et al. 2003. Study of crustal structure with S-wave data from Maqin-Jingbian profile. Acta Seismologica Sinica, 25(1): 82-88.
Liu M J, Mooney W D, Li S L, et al. 2006. Crustal structure of the northeastern margin of the Tibetan Plateau from the Songpan-Garzê terrane to the Ordos basin. Tectonophysics, 420(1-2): 253-266. DOI:10.1016/j.tecto.2006.01.025
Liu Q M, Zhao J M, Lu F, et al. 2014. Crustal structure of northeastern margin of the Tibetan Plateau by receiver function inversion. Science China Earth Sciences, 57(4): 741-750. DOI:10.1007/s11430-013-4772-5
Shearer P M. 1999. Introduction to Seismology. Cambridge: Cambridge University Press: 361-365.
Smith W H F, Wessel P. 1990. Gridding with continuous curvature splines in tension. Geophysics, 55(3): 293-305. DOI:10.1190/1.1442837
Su W, Wang C Y, Huang Z X. 2008. Azimuthal anisotropy of Rayleigh waves beneath the Tibetan Plateau and adjacent areas. Science in China Series D: Earth Sciences, 51(12): 1717-1725. DOI:10.1007/s11430-008-0137-x
Sun X L, Song X D, Zheng S H, et al. 2010. Three dimensional shear wave velocity structure of the crust and upper mantle beneath China from ambient noise surface wave tomography. Earthq. Sci., 23(5): 449-463. DOI:10.1007/s11589-010-0744-4
Tapponnier P, Xu Z Q, Roger F, et al. 2001. Oblique stepwise rise and growth of the Tibet plateau. Science, 294(5547): 1671-1677. DOI:10.1126/science.105978
The Research Group on Active Fault System around Ordos Massif, State Seismological Bureau. 1988. Active Fault System around Ordos Massif (in Chinese). Beijing: Seismological Press.
Wang C Y, Yang W C, Wu J P, et al. 2015. Study on the lithospheric structure and earthquakes in North-South Tectonic Belt. Chinese J. Geophys., 58(11): 3867-3901. DOI:10.6038/cjg20151101
Wang P, Wang Z G. 1997. Division of the Alxa block and its attribution. Earthquake, 17(1): 103-112.
Wang S X, Jiang F Y, Hao M, et al. 2013. Investigation of features of present 3D crustal movement in eastern edge of Tibet plateau. Chinese J. Geophys., 56(10): 3334-3345. DOI:10.6038/cjg20131010
Wielandt E. 1993. Propagation and structural interpretation of non-plane waves. Geophys. J. Int., 113(1): 45-53. DOI:10.1111/j.1365-246X.1993.tb02527.x
Wu J P, Huang Y, Zhang T Z, et al. 2009. Aftershock distribution of the Ms8.0 Wenchuan earthquake and three dimensional P-wave velocity structure in and around source region. Chinese J. Geophys., 52(2): 320-328.
Xu G M, Yao H J, Zhu L B, et al. 2007. Shear wave velocity structure of the crust and upper mantle in western China and its adjacent area. Chinese J. Geophys., 50(1): 193-208. DOI:10.3321/j.issn:0001-5733.2007.01.026
Xu X M, Ding Z F, Ye Q D, et al. 2015. The crustal and upper mantle structure beneath the South-North seismic zone from the inversion of Love wave phase velocity. Chinese J. Geophys., 58(11): 3928-3940. DOI:10.6038/cjg20151104
Yang Y J, Zheng Y, Chen J, et al. 2010. Rayleigh wave phase velocity maps of Tibet and the surrounding regions from ambient seismic noise tomography. Geochem. Geophys. Geosyst., 11(8): Q08010. DOI:10.1029/2010GC003119
Yao Z X, Wang C Y, Zeng R S, et al. 2014. Crustal structure in western Qinling tectonic belt and its adjacent regions deduced from receiver functions. Acta Seismologica Sinica, 36(1): 1-19. DOI:10.3969/j.issn.0253-3782.2014.01.001
Zhang H S, Gao R, Tian X B, et al. 2015. Crustal S-wave velocity beneath the northeastern Tibetan plateau inferred from teleseismic P-wave receiver functions. Chinese J. Geophys., 58(11): 3982-3992. DOI:10.6038/cjg20151108
Zhang J, Li J Y, Li Y F, et al. 2007. The Cenozoic deformation of the Alxa block in Central Asia-question on the northeastern extension of the Altyn Tagh Fault in Cenozoic Time. Acta Geologica Sinica, 81(11): 1481-1497.
Zhang P Z, Deng Q D, Zhang G M, et al. 2003. Active tectonic blocks and strong earthquakes in the continent of China. Science in China Series D: Earth Sciences, 46(S2): 13-24.
Zhang P Z, Shen Z K, Wang M, et al. 2004. Continuous deformation of the Tibetan Plateau from global positioning system data. Geology, 32(9): 809-812. DOI:10.1130/G20554.1
Zhang Q, Sandvol E, Ni J, et al. 2011. Rayleigh wave tomography of the northeastern margin of the Tibetan Plateau. Earth Planet. Sci. Lett., 304(1-2): 103-112. DOI:10.1016/j.epsl.2011.01.021
Zhang X K, Li S L, Wang F Y, et al. 2003. Differences of crustal structures in northeastern edge of Tibet Plateau, Ordos and Tangshan earthquake region in North China-results of deep seismic sounding. Seismology and Geology, 25(1): 52-60.
Zheng S H, Sun X L, Song X D, et al. 2008. Surface wave tomography of china from ambient seismic noise correlation. Geochem. Geophys. Geosyst., 9(5): Q05020. DOI:10.1029/2008GC001981
Zheng Y, Yang Y J, Ritzwoller M H, et al. 2010. Crustal structure of the northeastern Tibetan plateau, the Ordos block and the Sichuan basin from ambient noise tomography. Earthq. Sci., 23(5): 465-476. DOI:10.1007/s11589-010-0745-3
Zhou L Q, Xi J Y, Shen W S, et al. 2012. The structure of the crust and uppermost mantle beneath South China from ambient noise and earthquake tomography. Geophys. J. Int., 189(3): 1565-1583. DOI:10.1111/j.1365-246X.2012.05423.x
Zhou M D, Wang C Y, Zeng R S. 2012. Seismic tomography of the velocity structure of the crust and upper mantle in northeastern margin of the Qinghai-Tibet Plateau. Northwestern Seismological Journal, 34(3): 224-233. DOI:10.3969/j.issn.1000-0844.2012.03.0224
陈九辉, 刘启元, 李顺成, 等. 2005. 青藏高原东北缘—鄂尔多斯地块地壳上地幔S波速度结构. 地球物理学报, 48(2): 333–342. DOI:10.3321/j.issn:0001-5733.2005.02.015
邓起东, 张培震, 冉勇康, 等. 2002. 中国活动构造基本特征. 中国科学(D辑), 32(12): 1020–1030. DOI:10.3321/j.issn:1006-9267.2002.12.007
段永红, 张先康, 刘志, 等. 2007. 阿尼玛卿缝合带东段地壳结构的接收函数研究. 地震学报, 29(5): 483–491.
国家地震局《鄂尔多斯周缘活动断裂系》课题组. 1988. 鄂尔多斯周缘活动断裂系. 北京: 地震出版社.
郭飚, 刘启元, 陈九辉, 等. 2004. 青藏高原东北缘—鄂尔多斯地壳上地幔地震层析成像研究. 地球物理学报, 47(5): 790–797. DOI:10.3321/j.issn:0001-5733.2004.05.009
何正勤, 丁志峰, 叶太兰, 等. 2002. 中国大陆及其邻域的瑞利波群速度分布图象与地壳上地幔速度结构. 地震学报, 24(3): 252–259.
黄忠贤, 李红谊, 胥颐. 2013. 南北地震带岩石圈S波速度结构面波层析成像. 地球物理学报, 56(4): 1121–1131. DOI:10.6038/cjg20130408
嘉世旭, 张先康. 2008. 青藏高原东北缘深地震测深震相研究与地壳细结构. 地球物理学报, 51(5): 1431–1443. DOI:10.3321/j.issn:0001-5733.2008.05.016
李松林, 张先康, 张成科, 等. 2002. 玛沁—兰州—靖边地震测深剖面地壳速度结构的初步研究. 地球物理学报, 45(2): 210–217. DOI:10.3321/j.issn:0001-5733.2002.02.007
李永华, 吴庆举, 安张辉, 等. 2006. 青藏高原东北缘地壳S波速度结构与泊松比及其意义. 地球物理学报, 49(5): 1359–1368. DOI:10.3321/j.issn:0001-5733.2006.05.015
刘宝峰, 李松林, 张先康, 等. 2003. 玛沁—靖边剖面S波资料研究与探讨. 地震学报, 25(1): 82–88.
刘启民, 赵俊猛, 卢芳, 等. 2014. 用接收函数方法反演青藏高原东北缘地壳结构. 中国科学:地球科学, 44(4): 668–679.
苏伟, 王椿镛, 黄忠贤. 2008. 青藏高原及邻区的Rayleigh面波的方位各向异性. 中国科学D辑:地球科学, 38(6): 674–682.
王椿镛, 杨文采, 吴建平, 等. 2015. 南北构造带岩石圈结构与地震的研究. 地球物理学报, 58(11): 3867–3901. DOI:10.6038/cjg20151101
王萍, 王增光. 1997. 阿拉善活动块体的划分及归宿. 地震, 17(1): 103–112.
王双绪, 蒋锋云, 郝明, 等. 2013. 青藏高原东缘现今三维地壳运动特征研究. 地球物理学报, 56(10): 3334–3345. DOI:10.6038/cjg20131010
吴建平, 黄媛, 张天中, 等. 2009. 汶川Ms8.0级地震余震分布及周边区域P波三维速度结构研究. 地球物理学报, 52(2): 320–328.
喜马拉雅地震科学台阵. 2011. 中国地震科学探测台阵波形数据——喜马拉雅计划. 中国地震局地球物理研究所.
徐果明, 姚华建, 朱良保, 等. 2007. 中国西部及其邻域地壳上地幔横波速度结构. 地球物理学报, 50(1): 193–208.
徐小明, 丁志峰, 叶庆东, 等. 2015. 基于Love波相速度反演南北地震带地壳上地幔结构. 地球物理学报, 58(11): 3928–3940. DOI:10.6038/cjg20151104
姚志祥, 王椿镛, 曾融生, 等. 2014. 利用接收函数方法研究西秦岭构造带及其邻区地壳结构. 地震学报, 36(1): 1–19. DOI:10.3969/j.issn.0253-3782.2014.01.001
张洪双, 高锐, 田小波, 等. 2015. 青藏高原东北缘地壳S波速度结构及其动力学含义——远震接收函数提供的证据. 地球物理学报, 58(11): 3982–3992. DOI:10.6038/cjg20151108
张进, 李锦轶, 李彦峰, 等. 2007. 阿拉善地块新生代构造作用——兼论阿尔金断裂新生代东向延伸问题. 地质学报, 81(11): 1481–1497. DOI:10.3321/j.issn:0001-5717.2007.11.003
张培震, 邓起东, 张国民, 等. 2003. 中国大陆的强震活动与活动地块. 中国科学(D辑), 33(增): 12–20.
张先康, 李松林, 王夫运, 等. 2003. 青藏高原东北缘、鄂尔多斯和华北唐山震区的地壳结构差异——深地震测深的结果. 地震地质, 25(1): 52–60.
周民都, 王椿镛, 曾融生. 2012. 青藏高原东北缘地壳上地幔速度结构的地震层析成像研究. 西北地震学报, 34(3): 224–233. DOI:10.3969/j.issn.1000-0844.2012.03.0224