地球物理学报  2020, Vol. 63 Issue (3): 1104-1116   PDF    
青藏高原东南缘三江地区上地壳各向异性初步研究
吴鹏1,3,4, 高原2, 陈安国2, 石玉涛2     
1. 中国地震局地球物理研究所, 北京 100081;
2. 中国地震局地震预测研究所(地震预测重点实验室), 北京 100036;
3. 河北省地震局邯郸中心台, 河北邯郸 056001;
4. 河北省地震局测震专业仪器系统检测评估与观测场地遴选创新团队, 石家庄 050021
摘要:三江地区位于青藏高原东南缘,川滇菱形块体的西侧.受陆-陆板块碰撞作用的影响,构造活动强烈,地震活动频繁.为研究该区的应力环境、构造特征及二者间的关系,本文使用三江流动线性地震台阵SL-Array(2016-12—2017-05)和国家固定地震台网(2015-01—2017-12)记录到的地震波形资料,运用剪切波分裂分析技术研究三江地区上地壳各向异性研究.计算得到该区域快剪切波优势偏振方向为NNW向,与区域主压应力方向一致.结果显示各向异性特征有分区性,以维西—乔后断裂和小金河—丽江断裂为界,将线性台阵划分为A、B和C三个区域.A区快波优势偏振方向表现出与区域主压应力方向的一致性.B区局部构造复杂,快波优势偏振方向表现为近NS向.C区结果比较离散,无明显快波优势方向.自西向东,研究区域快波优势偏振方向表现为NNW至近NS向的变化趋势.计算得到研究区域慢剪切波时间延迟为2.8±1.7 ms·km-1,其中B区最大,A区最小,反映了该区地壳介质各向异性强度的不均匀分布,也揭示了区域构造复杂程度与地震各向异性强度的关系.
关键词: 青藏高原东南缘      三江地区      地壳各向异性      剪切波分裂     
Preliminary study on the anisotropy of the upper crust in the Sanjiang area, southeastern margin of the Tibetan Plateau
WU Peng1,3,4, GAO Yuan2, CHEN AnGuo2, SHI YuTao2     
1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Key Laboratory of Earthquake, Institute of Earthquake Science, China Earthquake Administration, Beijing 100036, China;
3. Handan Central Seismic Station, Hebei Earthquake Agency, Handan Hebei 056001, China;
4. Innocative Research Team of Hebei Earthquake Agency for Seismograghs Calibration, Valuation and Selection of Observation Sites, Shijiazhuang 050021, China
Abstract: The Sanjiang region, located on the southeastern margin of the Tibetan Plateau and the southwestern Sichuan-Yunnan rhombic block, has strong tectonic activity and frequent seismic activity due to the effect of continental collision. In order to study the stress environment, structural characteristics and their relationship in this region, this paper uses the seismic waveform data recorded by the Sanjiang Temporary Linear Seismic Array-SL Array (2016.12-2017.05) and the National Seismic Network (2015-01-2017-12) to study the anisotropy of the upper crust by applying the systematic analysis method (SAM). The results show that the dominant polarization direction of the fast shear waves in this region is NNW, which is consistent with the direction of the regional principal compressive stress. Bounded by the Weixi-Qiaohou fault and the Xiaojinhe-Lijiang fault, the study region can be divided into three subregions:A, B and C. The polarization direction of the fast waves in subregion A is NNW, which is consistent with the direction of the regional principal compressive stress. The local structure of subregion B is complex, and the polarization direction of the fast waves in this subregion is NS. The results of the subregion C are discrete, with no obvious fast wave dominant direction. From west to east, the dominant polarization direction of the fast waves in the SL Array area varies from NNW to NS. The time delay of the slow waves is 2.8±1.7 ms·km-1. Time delay in subregion B is the largest, however that in subregion A is the smallest. This phenomenon indicates that the anisotropy intensity in the crust is distributed unevenly and may be closely related to local tectonic complexity.
Keywords: Southeastern margin of the Tibetan Plateau    Sanjiang region    Seismic anisotropy    Shear-wave splitting    
0 引言

三江地区是由怒江、澜沧江和金沙江西向东展布、并流形成的区域,位于青藏高原东南缘,川滇菱形块体西南,是印度板块和欧亚板块碰撞的侧向地带,也是青藏高原物质向东南逃逸的重要区域(Tapponnier et al, 1982曾融生等,1992Royden et al., 2008Zhang et al., 2004Hou et al., 2004).由于晚古生代-中生代的构造演化和新生代大陆碰撞,这里发生了大规模成矿作用,形成了多金属富集区(邓军等,2014).三江及周边地区变形强烈,地震活跃,监控该区现今的应力状态,分析其与区域构造的关系,对研究青藏高原侧向生长模式和岩石圈变形机制具有重要意义.

地球内部广泛地存在各向异性现象,获取介质各向异性特征的途径之一是分析S波分裂参数(Crampin, 1978, 1981Vinnik et al., 1989Silver and Chan, 1991),包括快剪切波(快波)的偏振方向和慢剪切波(慢波)的时间延迟两个指标.研究表明,位于活动断裂上台站记录到的快波偏振方向与断裂走向一致,而远离断裂的台站记录到的快波偏振方向与区域主压应力方向平行(Crampin et al., 2003; 石玉涛等,2006Boness and Zoback, 2006Gao et al., 2011);慢波时间延迟反映了介质各向异性程度和应力场的变化(Crampin and Peacock, 2005Gao and Crampin, 20062008Gao et al., 2011, 2012Shi et al., 2012).因此,S波分裂结果可以反映区域应力状态与局部构造特征.通过S波分裂方法不但能够研究大尺度的地壳介质特征和应力状态,研究小尺度的应力环境特征也同样有效(石玉涛等,2009郑拓等,2017).利用介质的各向异性特征可以探讨区域构造的变形、区域应力分布以及断裂性质等问题(高原等, 2010, 2013).

地壳各向异性结果的可靠性,很大程度上取决于地震事件波形的数量和质量.青藏高原东南缘的固定台站数量较少,台间距较大,地震台站覆盖密度难以满足小尺度区域地壳各向异性研究的需要.因此,在重要的构造区需要额外架设流动地震台站,以获得密集的高质量波形数据.目前,由于缺乏足够多的高质量地震波形数据,该区域的地壳S波分裂研究成果较少,只是针对川滇地区强震震源区或大区域的各向异性特征开展了研究(雷军等,1997高原等,2004太龄雪等,2015李金等,2015).基于密集台阵观测研究三江地区各向异性特征的研究还是空白,需要做进一步的研究.

1 构造背景

三江地区位于欧亚板块和印度板块的侧向碰撞地带.在特提斯构造背景下,该区域内的不同陆块在不同时期经历了强烈的变形作用(管烨等,2004).受印度板块的推挤作用,第四纪以来,新构造活动强烈,深大断裂体系发育,地震活动频繁(Wang and Gao, 2014Li et al., 2018).

受大陆板块的侧向碰撞影响,三江地区形成了大规模逆冲断层(刘俊来等,2006),其东部为川滇菱形块体,中部为兰坪—思茅盆地,西部是腾冲—保山地块.金沙江缝合带和怒江缝合带从研究区内穿过.区域内有多个活动断裂展布,其中有近NS向的怒江西断裂、怒江断裂、澜沧江断裂、兰坪—云龙断裂等,NNW向的永平断裂、漕涧断裂、维西—乔后断裂、南涧—巍山断裂、红河断裂等,NE向的中甸—龙蟠—乔后断裂、小金河—丽江断裂、程海断裂等(图 1).其中红河断裂北起洱源,向东南方向延伸,全长970 km,是川滇菱形块体的西南边界,红河断裂向北延伸是维西—乔后断裂.丽江—小金河断裂把川滇块体分成南北两区,北区为川西北块体,南区是滇中块体.

图 1 三江地区构造背景与台站、震中分布图 红色方框代表研究区域;橙色粗线段代表一级块体边界,橙色细线段代表二级块体边界.F1怒江西支断裂;F2怒江断裂;F3漕涧断裂;F4澜沧江断裂;F5永平断裂;F6兰坪—云龙断裂;F7南涧—巍山断裂;F8维西—乔后断裂;F9中甸—龙蟠—乔后断裂;F10红河断裂北段;F11小金河—丽江断裂;F12程海断裂.图中地震来自中国地震台网中心2015年01月—2017年12月记录的ML4.0以下的地震目录. Fig. 1 Faults, seismic station and earthquakes distribution in Sanjiang area The Red box represents research area, the orange thick line indicates 1st-order boundary of tectonic blocks, and the orange light lines indicate 2nd-order boundary of tectonic blocks. The inset is tectonic background of this study area. Notes: F1 West Nujiang Fault; F2 Nujiang Fault; F3 Caojian Fault; F4 Lancang River Fault; F5 Yongping Fault; F6 Lanping-Yunlong Fault; F7 Nanjian-Weishan Fault; F8 Weixi-Qiaohou Fault; F9 Zhongdian-Longpan-Qiaohou Fault; F10 Red River Fault; F11 Xiaojinhe-Lijiang Fault; F12 Chenghai Fault. The earthquakes are from seismic catalogue recorded by China Earthquake Network Centre (CENC) from Jan., 2015 to Dec., 2017.

三江地区复杂的应力环境造成该区断裂构造极其发育,其走向大致分为NW、近NS和NE向.震源机制和GPS观测研究结果表明,三江地区的主压应力方向为近NS向(阚荣举等,1983杨国华等,2003).吴建平等(2004)利用中小地震求解震源机制解对云南地区构造应力场开展研究却显示区域主压应力场有顺时针旋转的趋势.此外,由于研究区位于中国大陆一、二级构造块体的过渡地带,地壳厚度变化剧烈,莫霍面从西北到东南逐渐变浅.研究结果表明云南西北地区,地壳厚度达到62 km,而南部地区,地壳厚度仅为33 km,S波速度结构具有很强的横向不均匀性(吴建平等,2001).Zhang和Gao(2019)利用中国科学台阵一期的流动地震台站数据,得到研究区域地壳厚度范围是46~55 km.由于研究区域构造复杂,对于该地区应力状态的认识尚未达成共识.因此,本文利用三江地区的流动密集地震台阵(SL-Array)和固定地震台站的波形数据,采用剪切波系统分析方法(SAM),对该地区开展地壳各向异性研究(高原等,2004b),试图为三江地区的地壳变形、应力和地壳结构特征提供新的约束.

2 资料与研究方法 2.1 数据资料

2016年11月,中国地震局地震预测研究所在三江地区架设流动线性地震台阵(SL-Array),共计28个流动台站,采用Nanometics公司的仪器,其中地震计型号Trillium 120P,周期为120 s;数据采集器型号Taurus,采样率为100 Hz.本文主要利用SL-Array线性台阵(2016年12月—2017年5月)和国家地震台网(2015年01月—2017年12月)记录到的小震波形资料进行地壳剪切波分裂研究.在此期间,根据中国地震台网中心提供的地震目录,在研究区内共获得ML4.0以下小震波形记录7548条(图 1).在计算过程中,本文使用40 km厚的单层速度模型,剪切波窗口设为50°.

2.2 研究方法

在本研究中,采用剪切波分裂系统分析方法-SAM方法开展地壳各向异性研究.SAM方法包含相关函数计算、时间延迟校正和偏振分析检验.该方法具有自我检验的功能,得到的结果具有很高的可靠性.图 2是进行滤波处理后台站SL06和EYA记录到的地震波形:(a)为台站SL06记录到的2017年01月19日ML1.6地震事件,震源深度13.0 km,震中距7.47 km;(b)图为台站EYA记录到的2016年11月15日ML2.0地震事件,震源深度9.0 km,震中距9.13 km.图 3图 4是利用SAM方法对台站SL06和EYA记录的地震波形进行S波分裂参数计算分析的实例.S波在通过各向异性介质时,会分裂成两列新的波列(即快波和慢波),对这两列波重新投影,分别计算它们在各个方位、不同时间延迟的相关函数,相关系数最大所对应的旋转角度和时间延迟即为快波的偏振方向和慢波的时间延迟.在实际计算过程中,由于观测波形存在未知的干扰分量,相关函数计算结果经常偏离较大,特别是时间延迟的误差较大,因此仅仅把相关函数计算结果作为参考,最后用可视化的偏振图进行检测得到最终结果.在偏振分析图中可以看到,在快波到达后、慢波未到时,S波质点运动轨迹的方向与正北方向成一夹角α,即为快波的偏振方向.慢波到达后,叠加到快波上,质点运动轨迹会发生变化(图 3).将两个分量的剪切波同时旋转α角度可以分离快、慢剪切波.理论上,快、慢剪切波来自同一震源的同一列剪切波,消除时间延迟后,质点的运动轨迹具有线性相关(图 4).经过时间延迟校正后,两个记录S波质点的运动轨迹线性度很高,证明了分析结果的可靠性.

图 2 波形记录图 (a)台站SL06记录到的2017年01月19日ML1.6地震事件波形,震源深度13.0 km,震中距为7.47 km;(b)台站EYA记录到的2016年11月15日ML2.0地震事件波形,震源深度9.0 km,震中距9.13 km;图中横坐标表示时间窗长度,纵坐标表示振幅count值. Fig. 2 Seismic waveforms at two seismic stations (a) Seismic waveforms of a small earthquake at station SL06. This earthquake occurred on Jan. 19, 2017, magnitude ML1.6, focal depth 13.0 km, epicenter distance 7.47 km; (b) Seismic waveforms of a small earthquake at station EYA. This earthquake occurred on NOV. 15, 2016, magnitude ML2.0, focal depth 9.0 km, epicenter distance 9.13 km. Abscissa is time window. Ordinate is amplitude in counts. Others are same as Fig. 2.
图 3 偏振分析图 (a)、(b)上图分别为台站SL06和台站EYA记录的快、慢剪切波质点运动轨迹图.图中S1、S2分别为快、慢剪切波初至点;下图分别为对应台站南北、东西向分量的剪切波波形,横坐标表示采样点的个数,纵坐标表示振幅count值,阴影矩形为所用的波形段. Fig. 3 Polarization analysis of shear-wave splitting (SWS) (a) The left is records at station SL06. (b)The right is records at station EYA. The upper is polarization diagrams of shear-waves, which plotting data are same as the shaded data in the lower. S1 and S2 point to arrivals of fast shear-wave and slow shear-wave, respectively. The lower is shear-waves of two horizontal components. Abscissa is number of sampling points. Ordinate is amplitude in counts.
图 4 经过时间延迟校正后的偏振分析检验图 (a)、(b)上图分别为两台站记录的波形经时间延迟校正后得到的快慢剪切波质点运动轨迹图.其中,台站SL06和台站EYA的时间延迟分别为0.09 s和0.06 s;下图分别为对应台站的快(F)、慢(S)剪切波的波形,其余同图 3. Fig. 4 Polarization analysis of fast and slow shear-waves after eliminating time delay All are the same as in Figure 3 only after eliminating time difference between fast shear-wave and slow one (i.e. time delay). (a) Station SL06, time delay is 0.09 s. (b) Station EYA, time delay is 0.06 s.
3 研究结果

通过计算分析,本文得到34个台站的剪切波分裂参数(表 1),其中有效记录在3条以上(含3条)的台站有25个.统计各台站的快剪切波偏振结果,得到了研究区域快剪切波偏振方向的下半球等面积投影玫瑰图(图 5).结合研究区各台站下方的剪切波分裂结果,计算得到区域快剪切波平均偏振方向为165°±30°.通过各台站下方快剪切波偏振方向的下半球等面积投影玫瑰图(图 6)可以看出,大部分台站的快剪切波优势偏振方向一致性较好,例如SL04-SL09、SL11、SL17和TUS、YOS等;一部分台站有两个快剪切波的优势偏振方向,例如SL01、SL13、SL14、SL19、SL21、EYA、HEQ、LUS等;个别台站快剪切波结果比较复杂,如SL15、SL18、SL21等,还有的台站有效结果较少(1~2条有效记录),例如SL03、SL20、SL23-SL26、SJ01、LIJ、YUL等.

表 1 三江地区台站下方地壳剪切波分裂结果 Table 1 SWS parameters at stations in Sanjiang zone
图 5 研究区快剪切波偏振方向的下半球等面积投影玫瑰图 Fig. 5 Lower hemispherical project and equal-area rose diagram of fast shear-wave polarization beneath seismic stations
图 6 台站下方快剪切波偏振方向的下半球等面积投影玫瑰图 图中蓝色和绿色玫瑰图分别表示流动地震台站和固定地震台站的快剪切波偏振方向. Fig. 6 Lower hemispherical project and equal-area rose diagram of fast shear-wave polarization beneath seismic stations Blue diagrams are fast shear-wave polarization from the temporary stations, and green diagrams are fast shear-wave polarization from the permanent stations.
4 讨论 4.1 快剪切波偏振方向

研究区内主要发育NW、近NS和NE向的断裂构造,SL测线横跨研究区的多条断裂(图 7).由计算结果可知,研究区台站下方的快剪切波偏振方向主要为NNW、近NS和NE向,与该区域内断裂走向具有较好的一致性,与太龄雪等(2015)利用中国地震科学台阵一期的观测数据在本区域得到的结果一致.结果显示,研究区的各向异性特征具有分区性.根据快剪切波偏振方向的空间分布特征和三江地区的构造背景,大致以维西—乔后断裂和小金河—丽江断裂为界,将SL测线上的台站(图 7中的阴影部分)分为A、B和C三个子区域,分别讨论其剪切波分裂特征.

图 7 快剪切波偏振方向分布 图中,黄色实线表示只有1~2条有效记录的结果;浅蓝色阴影部分为SL地震台站测线,蓝色虚线为分区界限;白色箭头为区域最大主压应力方向(钱晓东等,2011);黑色线段为背景噪声(8s)得到的方位各向异性结果(鲁来玉等,2014);地图上方的三个等面积投影玫瑰图为快波偏振方向分区玫瑰图,其他含义同图 1. Fig. 7 Spatial patterns of fast shear-wave dominant polarization in research area The yellow lines means number of SWS data is less than 3. The light blue shaded part is seismic stations from SL line, and the thick light-blue lines are boundary of subzones. A pair of large white arrows represents regional principal compressive stress. The black bars denote the azimuthal anisotropy in research area studied by surface wave tomography (8 s). Three rose diagrams are fast polarizations in three subzones. Others are same as in Fig. 1.

表 2给出了A、B、C三个子区域的剪切波分裂结果.可以看出(图 7),A区的快剪切波偏振方向为NNW向,与区域主压应力场方向一致;B区快波优势偏振方向为近NS向,有些台站具有两个优势偏振方向,有些台站结果比较离散,这一现象反映了该区域局部构造的复杂性;C区的快波偏振方向比较离散,没有明显的优势偏振方向,原因是台站获得的有效记录条数较少,近半数台站的结果只有1~2条,需要后期补充更多的数据来分析.研究区域内,自西向东,快波优势偏振方向表现为从NNW至近NS向的变化趋势,这与区域主压应力的结果(钱晓东等,2011)、GPS得到的地表应变结果一致(Gupta et al., 2015Jin et al., 2019).鲁来玉等(2014)利用背景噪声得到的方位各向异性结果的变化趋势与本文的结果也基本一致.

表 2 三江地区三个子区(图 7的阴影部分)剪切波分裂结果 Table 2 SWS parameters in three subzones in Sanjiang zone(The shaded part in Fig. 7)

A区台站位于维西—乔后断裂以西,区域内的构造单元从西往东依次是保山地块、澜沧江断裂带和兰坪—思茅盆地,断裂走向主要是NNW和近NS向,断裂性质多为逆冲型(刘俊来等,2006).该区的各台站下方快波偏振方向一致性较好,快波优势偏振方向为NNW向,与云南西部区域主压应力方向一致(图 7).

台站SL01位于怒江断裂东侧,其快剪切波有两个优势方向,第一优势方向为近EW向,与怒江断裂走向几乎垂直,第二优势方向为NNW向,与怒江断裂走向一致,但第二优势方向数据只有一条,其结果不具有代表性.因此,我们推测SL01主要受附近逆冲断裂的影响.台站SL02位于漕涧断裂西侧,快剪切波偏振方向为近NS向,与断裂走向、区域主压应力方向基本一致.台站SL05位于兰坪—云龙断裂上,其快波NNW向,与兰坪—云龙断裂走向平行.台站SL04和SL06分别位于兰坪—云龙断裂的两侧,距离断裂有一定距离,这两个台的快波偏振方向一致性好,均为NNW向,与区域主压应力方向平行.台站SL07-09远离断裂,他们的快波偏振方向为NNW向或NS向,这与区域主压应力方向一致或近似平行.台站SL03和SL26的有效记录条数太少,暂不做分析.研究表明,快剪切波偏振方向受区域应力场控制,与台站原地的主压应力方向一致;活动断裂上的台站的快波偏振方向与断裂走向一致(Crampin et al., 2003石玉涛等,2006吴晶等,2007).台站靠近断裂,其快波偏振方向会受到断裂的影响.A区台站基本上距离附近的断裂都有一定距离,个别台站(如SL05)离断裂较近,但该断裂走向也与其快波偏振方向基本一致,故A区的快波偏振优势方向主要反映了区域主压应力方向,这里显示的上地壳各向异性主要是应力作用下的大量充液微裂隙定向排列所致.

B区位于川滇菱形块体的边界,是多条深大断裂的交汇区域,构造环境非常复杂.区域内有比较发育的红河断裂及其向北延伸的维西—乔后断裂,NW走向的南涧—巍山断裂,NE走向的中甸—龙蟠—乔后断裂和小金河—丽江断裂等.观测台站分别位于红河断裂的两侧.受复杂构造影响,各台站的剪切波分裂结果相差较大.获得的B区快波优势偏振方向为NNE向(图 7).

台站SL11-向,直接反映了在块体边界、断裂交汇区域复杂的应力环境(图 6).台站SL11位于维西—乔后断裂的西侧,快波偏振方向与维西—乔后断裂走向一致;台站SL12虽然在维西—乔后断裂上,但是其快波偏振方向为NNE向,与断裂走向不一致,反而与东侧的中甸—龙蟠—乔后断裂走向平行,结合台站附近的地形特征可知,SL12台仪器布设于NNE走向的山脉脚下,其快波偏振方向与山脉走向平行,反映了地形地貌特征对快波偏振方向的影响.台站SL13位于交汇带内部,其有两个偏振优势方向,第一优势方向为NNW向,与台站西侧的维西—乔后断裂走向一致,同时与区域主要应力方向一致;第二优势方向为NEE向.SL14位于交汇地带的东侧,同样具有两个优势偏振方向,第一偏振方向为NNE向,与附近的中甸—龙蟠—乔后断裂走向平行,第二优势方向有效记录较少,暂不做分析.台站SL15-18处于中甸—龙蟠—乔后断裂与小金河—丽江断裂之间,属于川滇菱形块体内部,地形明显起伏.SL15快波第一优势方向为近NS向,与SL17的偏振方向一致,与小金河—丽江断裂走向几乎平行;第二优势方向为NWW向,与中甸—龙蟠—乔后断裂垂直,反映了台站原地的应力特征.台站SL16的快波偏振方向NE向,与附近的中甸—龙蟠—乔后断裂走向呈锐角相交.SL18结果较离散,但其快波优势偏振方向为近EW向.位于两条断裂之间的SL15-SL18,快波偏振受到区域应力场、断裂和地形的多重影响,故偏振方向比较离散.但B区全部有效数据放在一起,依然显示出明显的优势方向,且与区域主压应力方向一致.

C区位于丽江—小金河断裂以西,川滇菱形块体内部,区域内的主要断裂走向为NNE或NE向,该区域剪切波分裂结果比较离散且有效结果较少,无显著的优势偏振方向.台站SL19位于丽江—小金河断裂的南侧,有两个优势偏振方向,第一优势偏振方向均为近NS向,与区域主压应力方向基本一致.台站SL21同样具有两个优势偏振方向,其第一优势方向与SL22的偏振方向具有很好的一致性,都为NE向,与丽江—小金河断裂走向一致.该区的其他台站有效记录均在1~2条,需要后期更多的数据来分析.C区的快波偏振方向分布特征反映了该区域上地壳各向异性更多受到断裂的影响.当然,进一步的研究还需更多数据的支持.

除SL测线上的流动台站以外,我们还得到了7个固定地震台站和测线以外的2个流动台站SJ01、SJ02的剪切波分裂结果.台站LUS位于怒江西断裂西侧,具有NNW和NEE向两个快波优势偏振方向,第一优势方向与怒江西断裂的走向一致,也与钱晓东等(2011)得到的该地区主压应力方向一致,而NEE向的第二优势方向与怒江西断裂走向垂直,可能是逆冲断裂的影响.台站TUS位于红河断裂带上,其快波偏振方向为NNW向且一致性较好,与红河断裂走向一致.EYA台位于川滇菱形块体的西南边缘,同时也是红河断裂、南涧—巍山断裂、维西—乔后断裂和丽江—小金河断裂等多条断裂的交汇处,两个优势偏振方向为近NS向和NEE向,反映了台站下方复杂的应力环境.HEQ台位于川滇菱形块体的内部,其第一优势方向为近NS向,与区域主压应力方向一致;第二优势偏振方向为NE向,与丽江—小金河断裂走向一致.台站YOS位于程海断裂附近,快波偏振方向NNE向,与程海断裂走向平行.台站SJ02的偏振方向为NE向,与北侧丽江—小金河断裂走向平行,体现了局部构造对台站下方各向异性的影响.SJ01、LIJ、YUL等台站有效结果数量很少,需要后期更多的资料来分析.

研究区域全部有效数据的快剪切波平均偏振方向为165°±30°(参见图 5),与区域主压应力方向有很好的一致性(钱晓东等,2011).SL测线上得到的A、B、C三个分区中,A区各台站的快波偏振方向一致性较好,总体呈NNW向,与主压应力方向一致;B区快波偏振方向虽然较离散,但快波偏振优势方向明显为近NS向;表明在A、B两区本研究获得的上地壳各向异性主要由应力引起.C区无显著的优势偏振方向,反映了断裂等构造因素的明显作用,该区每个台站下方的各向异性反映了应力与原地构造的多重影响.

4.2 慢剪切波时间延迟

时间延迟反映了地下介质各向异性的程度,它是震源到台站路径上的剪切波传播路径的整体效应.研究认为,慢剪切波延迟时间一般在0.04~0.30 s之间(Crampin et al., 1984Silver and Chan, 1991),而位于复杂构造上的台站,时间延迟可能会增大至0.5 s(Savage et al., 2012).由于每个地震具有不同的震中距,为了统一描述介质各向异性程度,同时也方便比较,我们把时间延迟归一化,即计算单位射线长度上的时间延迟(单位为ms·km-1).

表 1可知,研究区域大部分台站的时间延迟在1.4~5.7 ms·km-1范围内,计算得到三江地区慢剪切波的平均延迟时间为2.8±1.7 ms·km-1.表 2给出了SL测线上的A、B、C三个区域的慢剪切波时间延迟,其中时间延迟最小的是A区,最大的是B区.在川滇块体的西部,兰坪—思茅地块和保山地块,主要受区域构造应力场的控制,断裂以及山脉走向与区域构造应力场方向基本一致,表现出的各向异性强度较弱.而在川滇块体的西边界区域,既是块体之间的过渡区域,又是多条断裂汇聚的地带,断裂之间的性质特征各不相同,表现出较强的各向异性特征.反映了研究区不同构造块体间地震各向异性强度的不均匀分布,同时揭示了区域内构造复杂程度与地壳各向异性强度的关系.研究表明,断裂构造及因断裂作用而在断裂附近产生的裂缝对地壳各向异性有很大影响(Savage et al., 2012),断裂附近的台站记录到的地震射线传播路径在通过断裂区时对时间延迟会产生影响,如怒江西断裂附近的LUS台,维西—乔后断裂与丽江—小金河断裂之间SL11-SL19台,位于多条断裂交汇处的EYA台等,其时间延迟均明显超过其他台的时间延迟.

Sun等(2012)利用中国地震局在青藏高原东南缘地区的固定台站,采用Pms转换波接收函数方法,获得12个台站的地壳各向异性参数,发现在青藏高原东南缘地区下方存在显著的地震各向异性,时间延迟范围在0.5~0.9 s.孙长青等(2013)利用云南地区布设的21个流动宽频带地震台,利用切向能量最小化方法获得Pms分裂参数,结果显示慢波时间延迟平均为0.25 s.接收函数得到的慢波时间延迟明显大于本文得到的结果.接收函数方法使用远震波形,损失了高频信号.因此,我们推测本文的各向异性结果主要来自于上地壳定向排列的微裂隙.

5 结论

本文采用SAM方法分析了云南流动线性地震台阵(2016-12—2017-05)和国家固定地震台网(2015-01—2017-12)记录的小震数据,获得了34个台站的314条有效记录,得到了研究区的剪切波分裂参数和上地壳各向异性的空间分布特征.

三江地区快剪切波平均偏振方向为165°±30°,与区域主压应力场方向一致,均为NNW向.研究结果显示,研究区的各向异性特征具有分区性.本文以维西—乔后断裂和小金河—丽江断裂为界,将SL测线上的流动地震台站划分为A、B和C三个区域.A区快波优势偏振方向为NNW向,且各台站结果的一致性较好,与区域主压应力方向平行.B区既是构造块体的交界,也是多条深大断裂的交汇区域,快剪切波偏振方向结果分布较凌乱,但优势方向表现为近NS向.C区快剪切波偏振方向比较离散,无明显优势方向.三江地区慢剪切波时间延迟为2.8±1.7 ms·km-1,其中时间延迟最大的是B区,最小的是A区,反映了研究区的地壳各向异性强度存在东西方向上的差异,揭示了活动断裂、地表构造复杂程度与地壳各向异性强度的联系.

受大陆侧向碰撞影响,研究区构造活动强烈,断裂发育.研究区复杂的构造环境决定了区域内复杂的应力状态.本文得到的剪切波分裂参数反映了研究区内上地壳各向异性的空间分布特征.自西向东,快波优势偏振方向表现出NNW向至NS向的变化,C区由于有效结果较少,无优势偏振方向.本文的结果使用的流动观测数据较少,只是初步结果,后期将利用更多的流动观测数据对该地区开展进一步的研究工作.

致谢  云南省地震局为本研究提供了固定地震台站的地震波形数据,两位审稿人为本文提出了宝贵意见,在此一并感谢.
References
Boness N L, Zoback M D. 2006. Mapping stress and structurally controlled crustal shear velocity anisotropy in California. Geology, 34(10): 825-828. DOI:10.1130/g22309.1
Crampin S. 1978. Seismic-wave propagation through a cracked solid:Polarization as a possible dilatancy diagnostic. Geophysical Journal International, 53(3): 467-496. DOI:10.1111/j.1365-246X.1978.tb03754.x
Crampin S. 1981. A review of wave motion in anisotropic and cracked elastic-media. Wave Motion, 3(4): 343-391. DOI:10.1016/0165-2125(81)90026-3
Crampin S, Chesnokov E M, Hipkin R G. 1984. Seismic anisotropy-the state of the art:Ⅱ. Geophysical Journal International, 76(1): 1-16. DOI:10.1111/j.1365-246X.1984.tb05017.x
Crampin S, Chastin S, Gao Y. 2003. Shear-wave splitting in a critical crust:Ⅲ. Preliminary report of multi-variable measurements in active tectonics. Journal of Applied Geophysics, 54(3-4): 265-277. DOI:10.1016/j.jappgeo.2003.01.001
Crampin S, Peacock S. 2005. A review of shear-wave splitting in the compliant crack-critical anisotropic Earth. Wave Motion, 41(1): 59-77. DOI:10.1016/j.wavemoti.2004.05.006
Deng J, Wang C M, Li W C, et al. 2014. The situation and enlightenment of the research of the tectonic evolution and metallogenesis in the Sanjiang Tethys. Earth Science Frontiers (in Chinese), 21(1): 52-64.
Gao Y, Liang W, Ding X, et al. 2004a. Variational characteristics of shear-wave splitting on the 2001 Shidian earthquakes in Yunnan, China. Acta Seismologica Sinica (in Chinese), 26(6): 576-582.
Gao Y, Liu X Q, Liang W, et al. 2004b. Systematic analysis method of shear-wave splitting:SAM software system. Earthquake Research in China (in Chinese), 20(1): 101-107.
Gao Y, Crampin S. 2006. A stress-forecast earthquake (with hindsight), where migration of source earthquakes causes anomalies in shear-wave polaridations. Tectonophysics, 426(3-4): 253-262. DOI:10.1016/j.tecto.2006.07.013
Gao Y, Crampin S. 2008. Shear-wave splitting and Earthquake forecasting. Terra Nova, 20(6): 440-448. DOI:10.1111/j.1365-3121.2008.00836.x
Gao Y, Wu J, Yi G X, et al. 2010. Crust-mantle coupling in North China:Preliminary analysis from seismic anisotropy. Chinese Science Bulletin (in Chinese), 55(31): 3599-3605. DOI:10.1007/s11434-010-4135-y
Gao Y, Wu J, Fukao Y, et al. 2011. Shear wave splitting in the crust in North China:stress, faults and tectonic implications. Geophysical Journal International, 187(2): 642-654. DOI:10.1111/j.1365-246X.2011.05200.x
Gao Y, Wu J T, Tai L X. 2012. Shear-wave splitting in the crust:regional compressive stress from polarizations of fast shear-waves. Earthquake Science, 25(1): 35-45. DOI:10.1007/s11589-012-0829-3
Gao Y, Wang Q, Zhao B, et al. 2014. A rupture blank zone in middle south part of Longmenshan Faults:Effect after Lushan MS7.0 earthquake of 20 April 2013 in Sichuan, China. Science China Earth Sciences, 57(9): 2036-2044. DOI:10.1007/s11430-014-4827-2
Guan Y, Wang A J, Cao D H, et al. 2004. East-West Structures in the Sanjiang Region, Western Yunnan and Their Implications. Acta Geologica Sinica (in Chinese), 78(4): 494-499.
Gupta T D, Riguzzi F, Dasgupta S, et al. 2015. Kinematics and strain rates of the Eastern Himalayan Syntaxis from new GPS campaigns in Northeast India. Tectonophysics, 655: 15-26. DOI:10.1016/j.tecto.2015.04.017
Hou Z Q, Gao Y F, Qu X M, et al. 2004. Origin of adakitic intrusives generated during mid-Miocene east-west extension in southern Tibet. Earth and Planetary Science Letters, 220(1-2): 139-155. DOI:10.1016/S0012-821X(04)00007-X
Jin H L, Gao Y, Su X N, et al. 2019. Contemporary crustal tectonic movement in the southern Sichuan-Yunnan block based on dense GPS observation data. Earth and Planetary Physics, 3(1): 53-61. DOI:10.26464/epp2019006
Kan R J, Wang S J, Huang K, et al. 1983. Modern tectonic stress field and relative motion of intraplate block in southwestern China. Seismology and Geology (in Chinese), 5(2): 79-90.
Lei J, Wang P D, Yao C, et al. 1997. The near-field shear wave splitting and its relation with structure in Jianchuan, Yunnan province. Acta Geophysica Sinica (in Chinese), 40(6): 791-801.
Li J, Gao Y, Xu F K, et al. 2015. The shear-wave splitting study on the Yingjiang MS6.1 earthquake sequence on May 30, 2014.. Earthquake Research in China (in Chinese), 31(2): 245-252.
Li Q S, Feng S Y, Bai Z M, et al. 2018. New research progress on crust-upper mantle structure in the southeastern margin of Qinghai-Tibet Plateau, China. Journal of Earth Sciences and Environment (in Chinese), 40(6): 757-778.
Liu J L, Song Z J, Cao S Y, et al. 2006. The dynamic setting and processes of tectonic and magmatic evolution of the oblique collision zone between Indian and Eurasian plates:Exemplified by the tectonic evolution of the Three River region, eastern Tibet. Acta Petrologica Sinica (in Chinese), 22(4): 775-786.
Lu L Y, He Z Q, Ding Z F, et al. 2014. Azimuth anisotropy and velocity heterogeneity of Yunnan area based on seismic ambient noise. Chinese Journal of Geophysics (in Chinese), 57(3): 822-836. DOI:10.6038/cjg20140312
Peng Z G, Ben-Zion Y. 2004. Systematic analysis of crustal anisotropy along the Karadere-Düzce branch of the North Anatolian fault. Geophysical Journal International, 159(1): 253-274. DOI:10.1111/j.1365-246X.2004.02379.x
Qian X D, Qin J Z, Liu L F. 2011. Study on recent tectonic stress field in Yunnan region. Seismology and Geology (in Chinese), 33(1): 91-106.
Royden L H, Burchfiel B C, Van Der Hilst R D. 2008. The geological evolution of the Tibetan Plateau. Science, 321(5892): 1054-1058. DOI:10.1126/science.1155371
Savage M K, Peppin W A, Vetter U R. 2012. Shear wave anisotropy and stress direction in and near Long Valley Caldera, California, 1979-1988. Journal of Geophysical Research:Solid Earth, 951(B7): 11165-11177.
Shi Y T, Gao Y, Wu J, et al. 2006. Seismic anisotropy of the crust in Yunnan, China:Polarizations of fast shear-waves. Acta Seismologica Sinica (in Chinese), 28(6): 574-585.
Shi Y T, Gao Y, Zhao C P, et al. 2009. A study of seismic anisotropy of Wenchuan earthquake sequence. Chinese Journal of Geophysics (in Chinese), 52(2): 389-407.
Shi Y T, Gao Y, Su Y J, et al. 2012. Shear-wave splitting beneath Yunnan area of southwest China. Earthquake Science, 25(1): 25-34. DOI:10.1007/s11589-012-0828-4
Silver P G, Chan W W. 1991. Shear wave splitting and subcontinental mantle deformation. Journal of Geophysical Research:Solid Earth, 96(B10): 16429-16454. DOI:10.1029/91JB00899
Sun C Q, Lei J S, Li C, et al. 2013. Crustal anisotropy beneath the Yunnan region and dynamic implications. Chinese Journal of Geophysics (in Chinese), 56(12): 4095-4105. DOI:10.6038/cjg20131214
Sun Y, Niu F L, Liu H F, et al. 2012. Crustal structure and deformation of the SE Tibetan plateau revealed by receiver function data. Earth and Planetary Science Letters, 349-350: 186-197. DOI:10.1016/j.epsl.2012.07.007
Tai L X, Gao Y, Liu G, et al. 2015. Crustal seismic anisotropy in the southeastern margin of Tibetan Plateau by ChinArray data:shear-wave splitting from temporary observations of the first phase. Chinese Journal of Geophysics (in Chinese), 58(11): 4079-4091. DOI:10.6038/cjg20151116
Tapponnier P, Peltzer G, Le Dain A, et al. 1982. Propagating extrusion tectonics in Asia:new insights from simple experiments with plasticine. Geology, 10(12): 611-616. DOI:10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2
Vinnik L P, Kind R, Kosarev G L, et al. 1989. Azimuthal anisotropy in the lithosphere from observations of long-period S-waves. Geophysical Journal International, 99(3): 549-559. DOI:10.1111/j.1365-246X.1989.tb02039.x
Wang Q, Gao Y. 2014. Rayleigh wave phase velocity tomography and strong earthquake activity on the southeastern front of the Tibetan Plateau. Science China Earth Sciences, 57(10): 2532-2542. DOI:10.1007/s11430-014-4908-2
Wu J, Gao Y, Chen Y T. 2007. Seismic anisotropy in the crust in northwestern capital area of China. Chinese Journal of Geophysics (in Chinese), 50(1): 209-220.
Wu J P, Ming Y H, Wang C Y. 2001. The S wave velocity structure beneath digital seismic stations of Yunnan province inferred from teleseismic receiver function modelling. Chinese Journal of Geophysics (in Chinese), 44(2): 228-237.
Wu J P, Ming Y H, Wang C Y. 2004. Source mechanism of small-to-moderate earthquakes and tectonic stress field in Yunnan Province. Acta Seismologica Sinica (in Chinese), 26(5): 457-465.
Yang G H, Wang Q, Wang M, et al. 2003. Current horizontal crustal movement and deformation characteristics in Yunnan area. Journal of Geodesy and Geodynamics (in Chinese), 23(2): 7-14.
Zeng R S, Zhu J S, Zhou B, et al. 1992. Seismicity and focal mechanism in Tibetan Plateau and its implications to lithospheric flow. Acta Seismologica Sinica (in Chinese), 14(S1): 523-533.
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 Z Q, Gao Y. 2019. Crustal thicknesses and Poisson's ratios beneath the Chuxiong-Simao Basin in the Southeast Margin of the Tibetan Plateau. Earth and Planetary Physics., 3(1): 69-84. DOI:10.26464/epp2019008
Zheng T, Ding Z F, Chang L J, et al. 2017. S-wave splitting in upper crust near the scientific drilling WFSD-3 at the Wenchuan earthquake fault zone. Chinese Journal of Geophysics (in Chinese), 60(5): 1690-1702. DOI:10.6038/cjg20170507
邓军, 王长明, 李文昌, 等. 2014. 三江特提斯复合造山与成矿作用研究态势及启示. 地学前缘, 21(1): 52-64.
高原, 梁维, 丁香, 等. 2004a. 云南2001年施甸地震的剪切波分裂参数变化特征. 地震学报, 26(6): 576-582.
高原, 刘希强, 梁维, 等. 2004b. 剪切波分裂系统分析方法(SAM)软件系统. 中国地震, 20(1): 101-107.
高原, 吴晶, 易桂喜, 等. 2010. 从壳幔地震各向异性初探华北地区壳幔耦合关系. 科学通报, 55(29): 2837-2843.
高原, 王琼, 赵博, 等. 2013. 龙门山断裂带中南段的一个破裂空段-芦山地震的震后效应. 中国科学:地球科学, 43(6): 1038-1046.
管烨, 王安建, 曹殿华, 等. 2004. 云南三江造山带近东西向构造特征及其研究意义. 地质学报, 78(4): 494-499. DOI:10.3321/j.issn:0001-5717.2004.04.008
阚荣举, 王邵晋, 黄琨, 等. 1983. 中国西南地区现代构造应力场与板内断块相对运动. 地震地质, 5(2): 79-90.
雷军, 王培德, 姚陈, 等. 1997. 云南剑川近场横波特征及其与构造的关系. 地球物理学报, 40(6): 791-801. DOI:10.3321/j.issn:0001-5733.1997.06.008
李金, 高原, 徐甫坤, 等. 2015. 2014年5月30日盈江6. 1级地震序列剪切波分裂研究.中国地震, 31(2): 245-252.
李秋生, 酆少英, 白志明, 等. 2018. 青藏高原东南缘壳幔结构探测研究新进展. 地球科学与环境学报, 40(6): 757-778. DOI:10.3969/j.issn.1672-6561.2018.06.006
刘俊来, 宋志杰, 曹淑云, 等. 2006. 印度-欧亚侧向碰撞带构造-岩浆演化的动力学背景与过程-以藏东三江地区构造演化为例. 岩石学报, 22(6): 775-786.
鲁来玉, 何正勤, 丁志峰, 等. 2014. 基于背景噪声研究云南地区面波速度非均匀性和方位各向异性. 地球物理学报, 57(3): 822-836. DOI:10.6038/cjg20140312
钱晓东, 秦嘉政, 刘丽芳. 2011. 云南地区现代构造应力场研究. 地震地质, 33(1): 91-106. DOI:10.3969/j.issn.0253-4967.2011.01.009
石玉涛, 高原, 吴晶, 等. 2006. 云南地区地壳介质各向异性-快剪切波偏振特性. 地震学报, 28(6): 574-585. DOI:10.3321/j.issn:0253-3782.2006.06.002
石玉涛, 高原, 赵翠萍, 等. 2009. 汶川地震余震序列的地震各向异性. 地球物理学报, 52(2): 398-407.
孙长青, 雷建设, 李聪, 等. 2013. 云南地区地壳各向异性及其动力学意义. 地球物理学报, 56(12): 4095-4105. DOI:10.6038/cjg20131214
太龄雪, 高原, 刘庚, 等. 2015. 利用中国地震科学台阵研究青藏高原东南缘地壳各向异性:第一期观测资料的剪切波分裂特征. 地球物理学报, 58(11): 4079-4091. DOI:10.6038/cjg20151116
吴晶, 高原, 陈运泰. 2007. 首都圈西北部地区地壳介质地震各向异性特征初步研究. 地球物理学报, 50(1): 209-220. DOI:10.3321/j.issn:0001-5733.2007.01.027
吴建平, 明跃红, 王椿镛. 2001. 云南数字地震台站下方的S波速度结构研究. 地球物理学报, 44(2): 228-237. DOI:10.3321/j.issn:0001-5733.2001.02.010
吴建平, 明跃红, 王椿镛. 2004. 云南地区中小地震震源机制及构造应力场研究. 地震学报, 26(5): 457-465. DOI:10.3321/j.issn:0253-3782.2004.05.001
杨国华, 王琪, 王敏, 等. 2003. 云南地区现今地壳水平运动与变形特征. 大地测量与地球动力学, 23(2): 7-14.
曾融生, 朱介寿, 周兵, 等. 1992. 青藏高原及其东部邻区的三维地震波速度结构与大陆碰撞模型. 地震学报, 14(S1): 523-533.
郑拓, 丁志峰, 常利军, 等. 2017. 汶川地震断裂带科学深钻WFSD-3附近上地壳S波分裂特征. 地球物理学报, 60(5): 1690-1702. DOI:10.6038/cjg20170507