引言

地震各向异性是地震学和地球动力学重要研究内容，可反映地壳介质变形和构造运动（Silver等，1988, 1991；Kaneshima，1990）。剪切波分裂方法可用于研究地震各向异性，Crampin提出EDA（extensive-dilatancy anisotropy）理论后（Crampin, 1978, 1981），部分学者进行了相关研究。已有研究表明，在应力作用下，EDA裂隙面平行于最大主压应力方向，使地震波呈各向异性（高原等，1999），剪切波进入各向异性介质裂隙中传播时，会分裂成快剪切波和慢剪切波，快剪切波优势偏振方向反映区域主压应力方向，慢剪切波时间延迟反映各向异性程度，时间延迟变化反映应力变化（Crampin, 1978, 1981；高原等，1999；Crampin等，2006）。因此，通过研究剪切波分裂可揭示研究区复杂的构造背景。

乌海市地处华北地块西北部、青藏高原块体东北部、鄂尔多斯块体西北部、银川地堑北端。银川地堑内存在第四纪以来活动性最强的NNE、NE向断裂带，这些断裂带控制了地堑的形成，且在青藏高原块体NE向挤压应力作用下，使地堑做右旋剪切运动（徐锡伟等，1994），使盆地内部和边缘发生走滑和正断层型地震（郭祥云等，2017）。

近年来，部分学者使用甘肃省、宁夏回族自治区、内蒙古自治区等地区地震台网连续多年近场地震数字观测资料，对鄂尔多斯块体西缘及周边地区开展地壳各向异性研究，分析鄂尔多斯块体西缘南北部剪切波分裂特征和区域主压应力场特征（张辉等，2012；太龄雪等，2017；许英才等，2019；张晖等，2020），但对鄂尔多斯块体西缘银川地堑北部剪切波分裂特征和地壳各向异性特征的研究仍不够深入。多位学者采用单台连续多年数字观测资料对特定地区地壳各向异性进行研究，取得了一定成果（邓嘉美等，2015；吴鹏等，2017）。本文采用SAM方法（高原等，2008a），利用乌海地震台（台站代码WUH，39.66°N，106.83°E）记录到的地震波形数据，对乌海地区地壳各向异性进行研究，得到乌海地区地壳各向异性特征和应力场特征，并对该地区剪切波分裂参数空间分布特征进行讨论，分析乌海地区地壳各向异性与当地应力场和周边断裂带的关系。

1.   地质构造背景

研究区位于银川地堑北端（图 1），银川地堑被阿拉善块体和鄂尔多斯块体所夹，西边界为NE走向贺兰山东麓断裂带，北边界为EW走向宗别立-正谊关断裂带（黄雄南等，2012）。银川地堑以东地区与鄂尔多斯块体相连，与吉兰泰盆地、贺兰山隆起共同构成鄂尔多斯块体断陷盆地的一部分。银川地堑受青藏高原块体NE向推挤作用，使银川地堑外侧NNE向断层发生了引张，形成一系列正断层，从而形成银川地堑。第四纪以来，由于青藏高原向东滑移，造成银川地堑南部NNW走向断层产生左旋走滑，形成了新褶皱，使银川地堑进一步扩张（马杏垣，1989）。银川地堑被三关口-牛首山断裂带划分为南北两个部分，两部分地质构造条件、地震活动性及应力场明显不同，南部属于青藏高原块体东北缘，北部属于华北块体西缘，北部主压应力方向为NNE向，南部主压应力方向为NEE近EW向（马禾青等，2011；太龄雪等，2017；许英才等，2019）。

图 1  研究区域构造背景

Figure 1.  Tectonics background and seismic stations in the study area

注：黑色细曲线表示市级行政边界，F1为贺兰东麓断裂带，F2为崇岗-芦花台断裂带，F3为银川-平罗隐伏断裂带，F4为黄河断裂带，F5为正谊关断裂带，F6为和屯-本井断裂带，F7为卓子山断裂带，F8为磴口-本井断裂带，F9为磴口-徐力斯特乌拉断裂带。

下载: 全尺寸图片 幻灯片

乌海地震台地处多个断裂带交汇地区，位于NNE向卓子山断裂带南端，西侧为NE向磴口-本井断裂带，地质构造环境较复杂。目前，对乌海地区地壳各向异性的研究较少，因此，对乌海地区剪切波分裂进行研究，以深入了解该地区地壳各向异性、地质构造和应力场特征。

2.   研究数据与方法

收集整理乌海地震台2014年1月至2020年6月地震观测资料，为消除地面剪切波全反射的影响，并考虑研究区速度结构（周民都等，2006；张永谦等，2011），选择地震射线临界角即剪切波窗≤50°的观测资料进行分析（Crampin等，2006；张晖等，2020），乌海地震台于2020年2月19日17:47:58记录到的震源深度21km、震中距16km、M_L2.1地震波形带通滤波后的波形如图 2所示。

图 2  带通滤波后的地震波形

Figure 2.  Bandpass filtered seismic waveform

下载: 全尺寸图片 幻灯片

研究方法为剪切波分裂系统分析方法，即SAM方法（高原等，2008a），研究形态地壳介质各向异性特征。SAM方法主要包括相关函数计算、时间延迟校正和偏振分析功能。该方法首先旋转2个水平分量，并进行相关函数计算，得到快剪切波优势偏振方向和慢剪切波时间延迟初值，然后依次分离快、慢剪切波，并对快、慢剪切波进行时间延迟校正，校正后的偏振图呈线性，乌海地震台于2020年2月19日17:47:58记录的震源深度21km、震中距16km、M_L2.1地震波形剪切波分裂过程如图 3所示，图中S1和S2分别表示快、慢波位置。对乌海地区地震资料进行剪切波窗挑选，剔除断记、干扰和无显著剪切波分裂特征的事件，最后得到65条地震波有效剪切波分裂计算结果，得到有效计算结果的地震事件震源深度主要分布在15—32km，震级范围为M_L0.1—2.9。

图 3  地震波形剪切波分裂过程

Figure 3.  Shear-waves splitting analysis of seismic wave recorded at WUH station

下载: 全尺寸图片 幻灯片

3.   研究结果与分析

3.1   剪切波分裂结果

研究得到乌海地区剪切波分裂参数，其中快剪切波平均优势偏振方向为NE63.1°±46.4°，慢剪切波平均时间延迟为（1.13±0.66）ms/km，如图 4所示。快剪切波平均优势偏振方向与鄂尔多斯块体西缘构造带走向一致（马禾青等，2011），与银川地堑北部主压应力方向一致，说明银川地堑北部应属于华北块体，与鄂尔多斯块体共同构成华北块体西缘。银川地堑主压应力方向主要与青藏高原块体NE向挤压作用有关（徐纪人等，2006），受始新世以来青藏高原NE向挤压作用，使鄂尔多斯块体西缘发生裂陷，燕山运动又导致古贺兰山解体（王美芳等，2008；刘建辉等，2010），从而使陷落区域形成银川地堑。

图 4  快剪切波优势偏振方向等面积投影玫瑰图及地震分布图

Figure 4.  WUH station fast wave polarization direction equal area projection rose diagram and earthquake distribution diagram

下载: 全尺寸图片 幻灯片

由图 4可知，乌海地震台快剪切波优势偏振方向表现出复杂的特征，基本存在4个优势偏振方向，分别为NE、EW、NNE、NNW向。为进一步分析快剪切波优势偏振方向与断裂带和地质构造的关系，分别对4个优势偏振方向进行分析，结果见表 1。由表 1可知，乌海地震台快剪切波平均优势偏振方向标准差较大，这主要因为存在多个优势偏振方向。

表 1  快剪切波优势偏振方向分析结果

Table 1.  Statistics of different fast wave polarization directions at WUH station

项目	地震波条数	快剪切波平均优势偏振方向/°	标准差	慢剪切波平均时间延迟/ms·km-1	标准差
NE向	24	51.50	10.88	1.23	0.71
EW向	13	88.08	13.23	1.04	0.53
NNE向	19	16.32	8.15	0.92	0.54
NNW向	9	156.56	15.51	1.40	0.75

下载: 导出CSV 

| 显示表格

3.2   快剪切波优势偏振方向与断裂带

乌海地震台共有19条地震波快剪切波优势偏振方向为NNE向，对应的19个地震事件震中分布如图 5所示。由图 5可知，地震主要集中在NNE向卓子山断裂带周边和NE向磴口-本井断裂带东侧，乌海地震台NNE向快剪切波优势偏振方向主要受卓子山断裂带和磴口-本井断裂带的影响。

图 5  NNE向快剪切波优势偏振方向等面积投影玫瑰图及地震分布图

Figure 5.  WUH station NNE direction fast wave polarization direction equal area projection rose diagram and earthquake distribution diagram

下载: 全尺寸图片 幻灯片

乌海地震台NE向快剪切波优势偏振方向等面积投影玫瑰图及地震分布图如图 6所示，由图 6可知，地震主要分布在卓子山断裂带，但磴口-本井断裂带周边地震显著增多，NE向作为优势偏振方向，同时受卓子山断裂带和磴口-本井断裂带的影响。

图 6  NE向快剪切波优势偏振方向等面积投影玫瑰图及地震分布图

Figure 6.  WUH station NE direction fast wave polarization direction equal area projection rose diagram and earthquake distribution diagram

下载: 全尺寸图片 幻灯片

磴口-本井断裂带和卓子山断裂带均位于银川地堑北部，为乌海地区主要断裂带，其中磴口-本井断裂带是晚更新世隐伏正断裂带，走向NE，倾向NW（邓起东等, 1994, 1999）；卓子山断裂带是燕山期断裂带，位于卓子山东麓地区，走向近NS，长约40km，断面倾向西，倾角约45°—65°，为高角度逆冲断层（卓鱼周，2015），被正谊关断裂带右行错开，分为南北两段（赵红格，2003）。由图 5、图 6可知，NNE向和NE向快剪切波优势偏振方向与断裂带分布及特征关系密切，这主要是因为局部构造和断裂带特性引起区域应力发生变化（高原等, 1999, 2010, 2008b；吴晶等，2007）。

乌海地震台EW向快剪切波优势偏振方向等面积投影玫瑰图及地震分布图如图 7所示，由图 7可知，地震主要分布在卓子山断裂带和磴口-本井断裂带南部、和屯-本井断裂带和正谊关断裂带北部。和屯-本井断裂带和正谊关断裂带均呈近EW走向（邓起东等, 1994, 1999），且正谊关断裂带横切了卓子山断裂带，将卓子山断裂带分为南北两部分，同时，已有研究认为卓子山断裂带逐渐由NS向转变为EW向（赵红格，2003），与乌海地震台快剪切波EW向优势偏振方向吻合，表明乌海地震台EW向优势偏振方向应受断裂带走向和局部构造共同影响。

图 7  EW向快剪切波优势偏振方向等面积投影玫瑰图及地震分布图

Figure 7.  WUH station EW direction fast wave polarization direction equal area projection rose diagram and earthquake distribution diagram

下载: 全尺寸图片 幻灯片

乌海地震台NNW向快剪切波优势偏振方向等面积投影玫瑰图及地震分布图如图 8所示，由图 8可知，地震主要分布在卓子山断裂带西侧、正谊关断裂带和贺兰山东麓断裂带东北侧。此优势偏振方向与周边断裂带走向不一致，也与主压应力方向不一致，部分已有研究得出类似结果（许英才等，2019）。贺兰山东麓断裂带东北方向延长线经乌海地震台附近，因此乌海地震台快剪切波优势偏振方向可能受贺兰山东麓断裂带构造局部应力的影响（张晖等，2020），已有研究表明在垂直运动大于水平运动的高倾角正断层所夹区域及逆冲断裂带均对剪切波分裂结果产生影响，造成剪切波分裂结果变化复杂（石玉涛等，2009；张晖等，2020），进而反映出研究区原地主压应力复杂的变化特征。

图 8  NNW向快剪切波优势偏振方向等面积投影玫瑰图及地震分布图

Figure 8.  WUH station NNW direction fast wave polarization direction equal area projection rose diagram and earthquake distribution diagram

下载: 全尺寸图片 幻灯片

4.   结论与讨论

基于内蒙古自治区乌海地震台2014年1月至2020年6月地震观测资料，采用SAM方法进行剪切波分裂研究。研究结果表明，乌海地区快剪切波平均优势偏振方向为NE63.1°±46.4°，慢剪切波平均时间延迟为（1.13±0.66）ms/km。乌海地震台快剪切波偏振存在4个优势偏振方向，分别为NE、EW、NNE、NNW向。由快剪切波优势偏振方向等面积投影玫瑰图及地震分布图可知，优势偏振方向NE和NNE向主要受主压应力场和磴口-本井断裂带及卓子山断裂带的影响，与NE向磴口-本井断裂带和NNE向卓子山断裂带走向一致；优势偏振方向EW向与卓子山断裂带南端被正谊关断裂带右旋所切有关，可能受卓子山和近EW向正谊关断裂带及和屯-本井断裂带共同影响；优势偏振方向NNW向与贺兰山东麓断裂带拉张方向一致，受贺兰山东麓断裂带第四纪以来抬升作用的影响，造成乌海地区复杂的原地主压应力特征，使乌海地震台剪切波分裂结果变化复杂。

乌海地区应力场主要受青藏高原块体NE向推挤作用的影响，同时还受磴口-本井断裂带、卓子山断裂带、和屯-本井断裂带、正谊关断裂带及贺兰山东麓断裂带的影响，因此，乌海地震台快剪切波优势偏振方向结果是主压应力场、断裂带、地质构造、原地主压应力共同作用的结果，显示出乌海地区复杂的地质构造和应力场特征。

由于本次研究仅使用了1个台站数据，数据量较少，因此，本文研究结果还需通过更多有效资料进行验证。