摘要利用鄂尔多斯及邻区的固定台站和流动台站记录到的面波资料，使用基于程函方程的面波层析成像的新方法，得到了鄂尔多斯及邻区12～150s周期的瑞利波相速度的分布图.短周期的相速度很好的揭示出鄂尔多斯及邻区的山脉隆起和拉张盆地等引起的高低速异常的分布特征.中等周期的相速度显示鄂尔多斯盆地内部为显著的高速异常，而它的四周逐渐地被低速异常所包围.大同火山区的低速异常在垂向上进行不断的延伸和扩展，推测是由于其下方的岩浆上涌引起地壳和上地幔顶部升温造成的.长周期的相速度显示鄂尔多斯盆地下方在岩石圈深度范围内，以北纬38°为界南北部岩石圈的差异比较明显，呈现南部厚北部薄的特点.研究发现研究区从西到东岩石圈逐渐减薄，这可能是由于受到太平洋板块俯冲引起的地幔热物质上涌和运移的影响.
关键词鄂尔多斯地区;面波层析成像;程函方程;瑞利波;相速度
0引言
鄂尔多斯及其邻区位于华北克拉通的西部，其中鄂尔多斯盆地具有稳定的克拉通基底的特性.鄂尔多斯地块受到来自青藏高原北东向的挤压作用，块体的四周被一系列的断陷盆地带所围，其地块内部比较稳定，而周围盆地带的运动比较复杂(张培震等，2002).鄂尔多斯块体的周围分布着银川—吉兰泰断陷带、河套断陷带、山西断陷带以及渭河断陷带，断陷带一般由若干条断裂带组成，这些断陷带的形成时期和运动方式不完全相同(国家地震局，1988;邓起东等，1999).鄂尔多斯盆地的周围分布着强震活动带，其盆地周缘发生过多次大地震，例如1303年的洪洞大地震、1556年的华县大地震、1695年的临汾大地震、1739年的平罗大地震、1920年的海原大地震等等(顾功叙，1983).
鄂尔多斯地块因其独特的地质条件，引起较多的地质学家和地球物理学家对该区域的研究.大尺度体波走时及远震体波走时层析成像揭示，鄂尔多斯下方的高速异常可以延续至200～300km深度，大同和山西地堑之下的低速异常可一直延伸至400～500km(HuangandZhao，2006;Sunetal.，2008;Zhaoetal.，2009;Huangetal.，2009;Lei，2012;TianandZhao，2013;Weietal.，2016;郭慧丽等，2017;陈兆辉等，2018).区域地震波走时层析成像揭示，青藏高原东北缘和鄂尔多斯块体西缘地区从地壳到上地幔顶部范围内的速度结构存在显著的横向不均匀性(杨婷等，2012;王椿镛等，2015;Chengetal.，2016;SunandKennett，2016;汪晟等，2017;郭慧丽和丁志峰，2018).地震面波和噪声层析成像(AnandShi，2006;Zhengetal.，2008，2011;何正勤等，2009;房立华等，2009;黄忠贤，2011;Wangetal.，2012;郑现等，2012;Tangetal.，2013;Legendreetal.，2014;Baoetal.，2015;范文渊等，2015;Shenetal.，2016)揭示了鄂尔多斯中下地壳及上地幔岩石圈具有较快的S波速度.由于受地震台站分布不均匀和密集程度的影响，有关鄂尔多斯块体内部及其周边拉张盆地的速度横向变化仍然只能获得粗略的图像，很难得到更加精细的深部结构的变化特征.
相关知识推荐：地球物理学期刊有哪些
中国地震局利用“地震科学探测台阵系统”的观测设备在鄂尔多斯及邻区进行了地震台阵观测.目前已有较多科研人员利用这些地震台阵观测的数据，运用不同的方法对鄂尔多斯及邻区的深部结构进行研究(郭慧丽等，2017;潘佳铁等，2017;Wangetal.，2017;王兴臣等，2017;钟世军等，2017;高翔等，2018).近年来，与鄂尔多斯相邻的华北地区增设了大量的流动台站，台站的覆盖性更好，分布更广泛，数据更加丰富.本文充分利用固定台站和流动台站的观测资料，使用基于程函方程面波成像的方法(Linetal.，2009;JinandGaherty，2015;钟世军等，2017)来研究鄂尔多斯及邻区的深部结构特征，进一步探讨鄂尔多斯盆地南北部岩石圈的差异以及大同火山活动相关的动力学问题.
1数据
本文使用的数据是“中国地震台阵探测”二期和三期以及其他布设的流动台站再加上固定台站记录的地震数据.本文总共使用了848个地震台站，其中固定台站有131个，流动台站则包括“中国地震台阵探测”二期246个和“中国地震台阵探测”三期309个及其他流动观测台站162个.研究区内的流动地震仪器主要包括两类：宽频带地震仪(地震计CMG-3ESPC和数采Reftek-130，频带范围60s～50Hz)和甚宽频带地震仪(地震计CMG-3T和数采Reftek-130S，频带范围120s～50Hz).流动台站和固定台站的分布位置较好的覆盖整个鄂尔多斯地块及其邻边地区(图1).迄今为止研究区的台站分布最完整、最密集，台站间距约35km.首先得到2006年10月至2018年5月之间的地震事件的波形数据，然后将满足震中距在10°～160°的范围内、震源深度小于50km，震级大于MS5.5条件下的事件波形数据筛选出来.图2给出了经过筛选得到的657个地震事件相对于台阵中心的位置分布.将筛选得到的地震波形数据进行重采样(1Hz)、去仪器响应、去均值、去倾斜分量、带通滤波等处理.在处理中，采用了3个频段频率范围分别是：12～20s、25～55s、60～150s，带通滤波参数分别为8～25s、20～60s、60～200s，分别参考相速度为3.4km·s-1、4km·s-1、4.3km·s-1，最大台站对间距分别为75km、150km和150km，光滑长度分别为0.25°.数据的初步处理可以减小一系列的误差，有利于得到分辨率较高的成像结果.
2方法
Lin等(2009)提出了一种高频近似条件下程函方程的面波层析成像方法，并且使用美国西部地区的噪声数据得到了分辨率较高的相速度分布图.基于程函方程面波层析成像的新方法不经过反演，直接测量不同频率下面波的走时场，通过计算走时场的梯度，然后确定不同周期空间慢度场的分布，最后得到某一频率下面波相速度的二维分布(冯吉坤等，2014).Jin和Gaherty(2015)提出了利用密集台阵阵列反演面波的相速度，在程函方程面波层析成像的基础上进行改进并且开发了基于多通道互相关的面波相速度自动测量系统(ASWMS)，通过该系统使用2006年到2014年美国台阵记录到的800多个事件的宽频带地震数据，得到美国大陆20～100s之间的瑞雷波的相速度分布图.其基本流程可分为两大部分：台站间的相延迟测量;通过程函方程进行相速度的计算.
3结果
本文利用基于程函方程的面波层析成像的新方法，得到鄂尔多斯及邻区12～150s周期的相速度，图5中展示了12s、16s、20s、30s、45s、60s、90s、120s、150s周期的相速度图像，在一定程度上很好地揭示了深部结构的横向变化特征.基阶瑞利波的相速度大约对1/3波长深度的S波比较敏感，我们计算了不同周期、不同深度的频散灵敏度系数曲线(图6)，其中一维初始速度模型(图6a)参考了Crust1.0.根据频散灵敏度系数曲线易知，不同周期的相速度图像很好的与地表沉积层、地壳、上地幔等不同深部结构的特征对应起来.
为了更直观的显示相速度结果的可靠性，我们画了相速度的误差统计图(图7).在12s的图像中，河套盆地附近存在较大的误差.短周期面波资料不够发育信噪比相对较差;沉积层区域的台站噪声相对较高;河套盆地位于研究区边缘地区，台站对覆盖不好，因此，该地区更容易受到影响.整体来讲，12～60s相速度的误差比较小能够很好的反映地壳和上地幔深部精细结构的特征，随着周期的增加，由于受到波长的影响导致相速度的误差会有所增加，但是在同一周期研究区域内的相速度误差呈均匀分布，这说明在一定程度上仍能很好的反映地下深部结构的特征.
根据频散灵敏度系数曲线可知，周期为12～20s的图像主要对深度在12～34km范围的S波速度敏感，能够反映地表沉积层和地壳的深部结构特征.河套盆地和银川盆地呈现出明显的低速异常，和钟世军等(2017)得到的结果比较吻合.鄂尔多斯盆地东南部的相速度高于西北部，这和郑现等(2012)利用背景噪声得到的结果是一致的.随着周期的增加，由于各盆地沉积层厚度的不同导致低速异常的变化幅度不同.山西断陷盆地由于具有大量的沉积层而显示低速;两侧的太行山和吕梁山具有较浅或者出露的基岩而显示较高的速度分布(周铭等，2016;王霞等，2019).短周期相速度的高低速异常的分布特征很好的将地形地貌的特征刻画出来.
周期为30～60s的图像主要对深度在28～136km范围内的S波速度敏感，能够反映下地壳和上地幔的深部结构特征.大同区域的中下地壳存在大面积低速异常的现象，这和郭震等(2015)利用噪声面波和布格重力异常联合反演得到的结果是一致的.大同火山区的低速异常从中下地壳到上地幔的顶部进行不断的延伸和扩展，低速异常变得越来越明显且范围不断扩大，与已有的研究结果一致(Wangetal.，2012;杨婷等，2012;宋美琴等，2013;王霞等，2019).利用接收函数得到鄂尔多斯块体的地壳厚度变化比较平缓约为44km，华北克拉通中部地壳厚度相对鄂尔多斯地块较薄，在鄂尔多斯块体及其周围整体上表现出西部厚，东部薄，北部厚，南部薄的分布特点(Wangetal.，2017;娄辛辉等，2017).根据成像结果可知，相速度的图像在一定程度上大致反映出研究区地壳厚度差异的分布特征.山西断陷带以北纬38°为分界线在下地壳和上地幔的南北部存在明显的速度差异的分布特征，与现有结果具有较好的一致性(宋美琴等，2013;周铭等，2016;王霞等，2019).随着周期的增加，鄂尔多斯盆地的四周逐渐连成一片呈低速异常的分布特征，内部整体上呈现明显的高速异常，这与盆地内部相对稳定、地震活动比较弱、断裂比较少，而周围分布一系列断陷盆地的特征相吻合.
周期为90～150s的图像对深度在82～330km范围的S波速度比较敏感，主要反映了深部岩石圈的结构特征.鄂尔多斯盆地下方的高速异常能够延续至200km，说明盆地下方具有较厚的岩石圈结构，这和已有的研究成果吻合较好(AnandShi，2006;Chenetal.，2009;Chen，2010;李多等，2012;TianandZhao，2013;于勇等，2016;陈兆辉等，2018;高翔等，2018).鄂尔多斯盆地在岩石圈范围内以北纬38°为分界线，北部是低速异常，南部是高速异常.随着周期的增加，鄂尔多斯盆地在岩石圈范围分界线两侧的速度差异变得越来越明显，这和已有研究具有较好的一致性(Tianetal.，2009;高翔等，2018).
4讨论
鄂尔多斯块体周边断陷系的动力学成因问题一直受到地学界的高度关注，学者们提出的不同模式曾引起强烈争论.其中有两种端元模式，一种强调大型走滑断裂拉分作用的重要性，认为沿青藏高原东北缘发育的断裂带及其向东扩展主导了华北地区晚新生代盆地的发育和构造变形特征(Tapponnieretal.，1982;Peltzeretal.，1985;Gaudemeretal.，1995;张岳桥等，1999)，印度—欧亚陆陆碰撞作用是其主要成因机制.另一种观点认为，鄂尔多斯盆地周缘新生代断陷盆地与太平洋板块向亚洲大陆俯冲产生的弧后扩张有关，属于弧后盆地(UyedaandKanamori，1979;Tianetal.，1992;Northrupetal.，1995).也有学者认为，断陷盆地的发生是上述两种板块动力系统相互作用的结果(张岳桥等，2006;Chen，2010;陈兆辉等，2018).
相速度的成像结果表明，研究区的主要特征是鄂尔多斯盆地逐渐的被低速异常包围，整个华北克拉通中部为显著的低速异常.其次以北纬38°为分界线，山西断陷带的南北部速度差异在下地壳和上地幔比较明显，而鄂尔多斯盆地的南北部速度差异在深度200km附近比较明显.
研究发现整个研究区以北纬38°为分界线的南北部从浅部沉积层到深部岩石圈范围内均存在明显的速度差异的分布特征.分界线的北部有银川盆地、河套盆地、大同火山区等一系列的断陷盆地，上地壳相速度的分布特征主要受到盆地沉积层的影响.方盛明等(2009)利用深地震反射探测剖面获得银川地堑内新生代沉积最厚约9km.Wang等(2017)利用接收函数研究表明，分界线北部的沉积层比南部的厚，其中河套盆地沉积层的厚度超过5km.对于山西断陷带在下地壳和上地幔范围内南北部差异的成因，周铭等(2016)推测山西断陷带两侧的构造应力和构造运动方式不同导致地幔物质上涌.鄂尔多斯盆地在岩石圈范围内的南北部差异，Liu等(2004)认为鄂尔多斯盆地北部上地幔存在的低速体是由于青藏高原上地幔物质运移造成的，已有研究(Tianetal.，2011;Dongetal.，2014)认为是太平板块俯冲引起地幔热物质上涌造成的.华北克拉通下方的岩石圈整体上从西到东逐渐减薄(Chen，2010;朱日祥等，2011;黄方等，2015;于勇等，2016)，这说明华北克拉通的岩石圈正在遭受不同程度的破坏.对于鄂尔多斯地区的构造应力场，目前认为鄂尔多斯块体主要受到青藏高原的挤压及太平洋板块俯冲的共同作用.前人研究(Huangetal.，2008;常利军等，2011;于勇等，2016;GuoandChen，2017)在鄂尔多斯盆地南部存在一个狭隘的通道，来自青藏高原的软流圈通过绕流进入鄂尔多斯盆地的东南部，对鄂尔多斯盆地东南部的岩石圈进行侵蚀和改造.目前我们对这一说法持有怀疑的态度，假设来自青藏高原软流圈绕流观点成立，软流圈经过鄂尔多斯盆地南部时应该对其破坏更加严重，然而成像结果表明鄂尔多斯盆地南部的岩石圈比北部的厚.已有研究表明大同火山群具有大范围的新生代幔源岩浆岩出露(Liuetal.，2004;Xuetal.，2005)，大同火山区从下地壳开始不断的向周围延伸，随后整个华北克拉通中部为显著的低速异常，相比于来自青藏高原的软流圈绕流到鄂尔多斯盆地东南部的观点，我们更倾向于太平洋板块俯冲引起深部热物质上涌和运移这一观点.Lei(2012)成像结果表明大同火山和渤海海底的低速异常在200km以下相连呈一个Y形结构且延伸至下地幔，认为大同火山尽管与太平洋板块深俯冲形成的“大地幔楔结构与动力学”(LeiandZhao，2005，2006;HuangandZhao，2006;Leietal.，2013;Zhangetal.，2020)有关，但更可能的是起源于下地幔的地幔柱(Lei，2012).本文结果也揭示出大同火山下方从地壳到上地幔存在一个显著的柱状低速异常，因此，我们认同大同火山更有可能起源于下地幔的地幔柱这一观点.
5结论
本文使用基于程函方程的面波层析成像的方法，利用流动台站和固定台记录到的地震面波数据，得到了鄂尔多斯及邻区瑞利波相速度的分布图.相速度分布图很清晰的反映出研究区范围内的地壳和上地幔的速度结构的精细特征，与地质特征吻合较好.
短周期时，瑞利波相速度的高低速异常的分布特征与地形地貌很好的对应.中等周期时，鄂尔多斯盆地的四周逐渐的被低速异常包围，这和鄂尔多斯盆地周围分布一系列的断陷盆地相吻合.长周期时，鄂尔多斯盆地和山西断陷带以北纬38°为分界线南北部的岩石圈厚度差异比较明显.
大同火山区从下地壳到上地幔不断地延伸和扩展随后整个华北克拉通中部为显著的低速异常，其次整个研究区以北纬38°为界从浅部沉积层到深部岩石圈的范围内南北部的结构特征存在明显的差异，我们认为主要由于太平洋板块俯冲的作用引起热物质上涌和运移从而导致鄂尔多斯及邻区的岩石圈遭到不同程度的破坏.致谢感谢中国地震局地球物理研究所“中国地震科学探测台阵数据中心”为本文提供的地震数据资料;感谢两名审稿专家的评阅.——论文作者：蔡光耀1，王未来1*，吴建平1，2，房立华1，2
