地球物理学报  2019, Vol. 62 Issue (12): 4571-4587   PDF    
川滇交界东段昭通、莲峰断裂带的深部结构特征与2014年鲁甸MS6.5地震
李大虎1,2, 丁志峰2, 吴萍萍3, 梁明剑1, 吴朋1, 顾勤平4, 康清清4     
1. 川滇国家地震预报实验场四川分中心, 成都 610041;
2. 中国地震局地球物理研究所, 北京 100081;
3. 防灾科技学院地震科学系, 河北燕郊 056201;
4. 江苏省地震局, 南京 210014
摘要:2014年鲁甸MS6.5地震位于川滇菱形块体向东突出的过渡变形区大凉山次级块体南东缘的昭通、莲峰断裂带内部,属于青藏高原东南缘南北地震带的中南段,近十多年来,该断裂带及其周边中强地震的发生频次明显增多,昭通、莲峰断裂带是否具备孕育和发生强震的深部构造背景成为一个亟待研究的问题.为了研究昭通、莲峰断裂带的深部结构特征及孕震背景,探求2014年鲁甸MS6.5地震的成因的深部动力机制,本文充分收集了四川、云南等区域数字地震台网和"中国地震科学台阵探测-南北地震带南段"("喜马拉雅"项目Ⅰ期)流动地震台阵的观测数据,应用区域震和远震联合反演的方法得到川滇地区三维速度结构图像,在此基础上重点剖析和研究了昭通、莲峰断裂带P波速度结构;再对昭通、莲峰断裂带及周边区域的重力、航磁数据进行三维视密度和视磁化强度反演,得到了壳内不同深度层视密度的横向变化特征和反映壳内磁性物质的分布范围以及结晶基底的视磁化强度异常分布情况,综合分析研究昭通、莲峰断裂带的深部结构特征及孕震动力环境.研究结果表明:川滇交界东部昭通、莲峰断裂带及其周边地区上地壳物质存在显著的横向介质差异,中下地壳深度范围大凉山次级块体西南缘存在低速异常分布,并呈现出近SN向的展布特征,2014年鲁甸MS6.5地震位于该高低速异常的分界线附近略偏向高速体一侧.P波速度结构还揭示了鲁甸MS6.5主震震源体下方中下地壳存在大范围低速异常分布,P波速度异常扰动与重磁异常的展布特征、梯度变化在深度和分区特征上均具有较好的联系和可比性,结合昭通、莲峰断裂带中下地壳范围内存在大范围的低密度弱磁性异常分布,综合表明了该区中下地壳物质相对较为软弱,这种特有的深部物性结构特征有利于应力在脆性的上地壳内积累和集中.研究结果还揭示了共轭断裂的深部构造形态,高低航磁异常边界与NW向的苞谷脑—小河断裂的深部展布形态相一致,苞谷脑—小河断裂处于航磁异常突变带附近,昭通断裂北段(昭通—鲁甸段)位于上地壳强磁性、高波速异常区内且具有深大断裂的深部地球物理场响应特征,因此该断裂段(昭通—鲁甸段)具备发生7级及以上强震的深部构造背景.当大凉山次级块体内部的中下地壳低速管流层自NW向SE方向运动到昭通、莲峰断裂带附近时,受到华南块体的强烈阻挡,应力在昭通、莲峰断裂附近基底性质存在差异处集中,脆性上地壳中低强度区域在横向挤压的构造应力场作用下易于破裂从而引发强震,这也正是昭通、莲峰断裂带内部鲁甸MS6.5地震孕育和发生的深部构造环境.
关键词: 昭通、莲峰断裂带      鲁甸MS6.5地震      速度结构      视密度、视磁化强度反演     
Deep structure of the Zhaotong and Lianfeng fault zones in the eastern segment of the Sichuan-Yunnan border and the 2014 Ludian MS6.5 earthquake
LI DaHu1,2, DING ZhiFeng2, WU PingPing3, LIANG MingJian1, WU Peng1, GU QinPing4, KANG QingQing4     
1. Sichuan-Yunnan Earthquake Prediction Experimental Site, Sichuan Branch, Chengdu 610041, China;
2. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
3. Institute of Disaster Prevention, Langfang Hebei 056201, China;
4. Earthquake Administration of Jiangsu Province, Nanjing 210014, China
Abstract: The 2014 Ludian MS6.5 earthquake occurred in the Zhaotong and Lianfeng fault zones on the southeastern margin of the Daliangshan secondary block in the eastward-protruding transitional deformation zone of the Chuandian diamond block, part of the north-south seismic belt in the southeastern margin of the Tibetan plateau. Over the past ten years, moderate-sized earthquakes on this fault zone and its surrounding areas have increased notably. Whether the deep structure of this fault zone can generate major earthquakes is an urgent issue to be addressed. This work attempted to answer this question and examines the deep dynamic mechanism of the 2014 Ludian MS6.5 earthquake. We collected data from seismic stations of the regional digital seismic networks of Yunnan and Sichuan and mobile seismic arrays of Chinese seismic array detection-the southern section of North-South seismic belt ("Himalaya" project phase 1). Then we imaged the three-dimensional velocity structure beneath the Sichuan-Yunnan area using joint inversion of regional and teleseismic data. On this basis, the P-wave velocity structure of the Lianfeng and Zhaotong fault zones was studied exhaustively. And then we performed inversion gravity and aeromagnetic data to obtain apparent density, apparent magnetization intensity in the study area, which permits to reveal the deep dynamic environment of earthquakes. The results show that there are significant lateral medium differences in the upper crust of beneath the Zhaotong and Lianfeng fault zones and adjacent areas. The low velocity anomalies are present in the middle-lower crust below the southwestern edge of the Daliang Shan secondary block, showing a nearly NS trend. The 2014 MS6.5 earthquake is located near the boundary between the high and low velocity anomalies and slightly tend toward the side of the high-velocity body. The P-wave velocity structure also reveals a large-scale low velocity anomaly distribution in the middle-lower crust below the MS6.5 mainshock source. The distribution of P-wave velocity perturbation and gravitational and magnetic anomalies have better connection and comparability in both the depth and the partition. In combination with the distribution of low-density and weak magnetic anomalies in the middle-lower crust of the Zhaotong-Lianfeng fault zone, it is inferred that the middle-lower crust in this area is relatively weak. Such a deep structure is easy to accumulate stress in the brittle upper crust. The results also reveal deep structure of conjugate faults, and the boundary between high and low aeromagnetic magnetic anomalies is consistent with the NW-trending Baogunao-Xiaohe fault. The Zhaotong-Ludian section of the Zhaotong fault zone is located in the strong magnetic and high-velocity anomaly zone of the upper crust and has response to deep geophysical fields, so the fault segment (Zhaotong-Ludian section) has deep structural background of major earthquakes. Research suggests that the low-velocity channel flow in the middle-lower crust of the Daliang Shan block moves from NW to SE, and is strongly blocked by the South China block in the vicinity of the Lianfeng and Zhaotong fault zones, thus resulting in stress accumulation near the basement of the Zhaotong and Lianfeng faults where there are significant differences of properties. Consequently the low-strength regions of the brittle upper crust tend to rupture under the lateral compressive tectonic stress field leading to strong earthquakes. This is just the deep tectonic setting of the 2014 Ludian MS6.5 earthquake in the Zhaotong and Lianfeng fault zones.
Keywords: Zhaotong and Lianfeng fault zones    Ludian MS6.5 earthquake    Velocity structure    Apparent density-Apparent magnetization inversion    
0 引言

据中国国家地震台网测定,北京时间2014年8月3日16时30分,云南省昭通市鲁甸县境内发生了MS6.5地震,震中位于鲁甸县龙头山镇(103.3°E,27.1°N),震源深度约12 km(CENC, http://www.cenc.ac.cn/).此次地震造成昭通市鲁甸县、巧家县、昭阳区、永善县和曲靖市会泽县等地大量的人员伤亡和财产损失,是2008年以来伤亡人数仅次于汶川8.0级和玉树7.1级的我国大陆破坏性地震,据野外考察可知鲁甸MS6.5地震烈度最高达Ⅸ,等震线长轴呈NW-NNW向(徐锡伟等,2014).

从区域构造位置来看,昭通、莲峰断裂带由昭通、莲峰两条NE向断裂组成,平行展布于川滇交界东段、四川大凉山南部至云南昭通之间的地区,局部进入贵州威宁境内(图 1).该断裂带是在古生代时与四川盆地的华蓥山断裂带同期发育、中生代时进一步发展的以挤压逆冲为主的区域性大断裂,属于古华蓥山断裂带的一部分(闻学泽等,2013).其中,莲峰断裂东北起自盐津以西,向西南经莲峰,止于巧家以北,全长150 km,总体走向35°~45°,倾向NW.莲峰断裂作为一条区域性的大断裂,对地层发育与区域构造变形起控制作用.昭通断裂东北起自盐津东南,向西南经彝良、昭通、鲁甸、会泽,止于巧家以南小江断裂带东侧,长150 km,昭通断裂与莲峰断裂平行,对地层发育与区域构造变形也起到明显控制作用.川滇活动地块东边界的中-长期大地震潜势,无论在地震预测研究还是在防震减灾的需求上都非常值得探索与研究(闻学泽等,2008; 易桂喜等, 2008).昭通、莲峰断裂带所属的川滇块体东边界及其以东地区是全国地震重点危险区之一,据2012年中国地震台网中心中国大陆及川滇地区地震大形势分析,与去年相比该危险区北边界向南收缩,预测的地震水平为6~7级,强震危险性不容忽视.尤其是近十多年来,该断裂带及其附近中强地震的发生频次明显增多,如2003年鲁甸MS5.0、5.1地震;2004年鲁甸MS5.6地震、2006年盐津MS5.1地震;2010年盐津MS4.6地震;2012年彝良MS5.7、5.6地震;2013年盐津MS4.2地震;2013年宁南—云南巧家交界MS4.9、3.1地震;2014年4月5日永善县MS5.3地震;2014年8月3日鲁甸MS6.5地震、8月17日永善MS5.0地震,中强地震的频发使得昭通、莲峰断裂带是否具备孕育和发生强震/大地震的深部构造背景成为一个亟待研究的问题.

易桂喜等(2010a)利用地震活动性参数b值的空间分布,结合历史强震与现今地震活动背景,分析了莲峰断裂带的现今活动习性和潜在强震危险性,研究结果表明存在于主干断裂附近巧家与东川间以及嵩明北西的凹凸体将是川滇活动地块东边界未来大震或强震的震源区.易桂喜等(2016)基于波形资料计算获得的2014年10月越西M5.0地震的视应力值(0.99 MPa)高于汶川8.0级地震后其附近区域的地震视应力,表明了越西地震的视应力高值可能反映了川滇块体整个东边界区域当前的应力水平在逐渐升高,发生强震的危险性也在升高,其潜在的地震危险背景值得关注.闻学泽等(2013)的研究结果表明,莲峰断裂带和昭通断裂带为晚第四纪-全新世活动断裂,具逆冲-右旋走滑运动性质,且GPS观测结果也表明该地区水平缩短速率南段5~7 mm·a-1,中-北段3~4 mm·a-1,有可能是被莲峰断裂带和昭通断裂带所吸收,这些证据均表明昭通、莲峰断裂带南、北两个段落已不同程度闭锁,其中,昭通断裂带具有更高的闭锁与应变积累程度,具有发生强震-大地震的构造能力.鲁甸附近的高应力段很可能与鲁甸—彝良之间的小震空白段组成一个强闭锁的断裂段,其潜在地震的最大矩震级估值为MW7.4左右.陈石等(2014)通过对鲁甸地震震源区及周边的三维密度结构特征研究结果表明震源区南北两侧块体运动受到昭通断裂不同程度阻隔,在川滇交界复杂构造体系下的地壳物性结构差异性非常显著,并指出需要结合区域布格异常等其他地球物理场资料来对地壳深部结构与差异性特点进行评估,发现可能存在的深部隐伏断裂.Xu等(2015)利用丽江—攀枝花—清镇剖面的人工源宽角地震数据重建峨眉山大火成岩省的地壳速度结构,研究结果对理解峨眉山大火成岩省的深部结构与属性特征以及2014年鲁甸MS6.5地震孕震环境具有重要的意义.

图 1 昭通、莲峰断裂带区域地震构造环境图 Fig. 1 Regional seismotectonic environment of the Zhaotong、Lianfeng fault zone

已有的人工地震测深资料由于受到了剖面数量不足、分布不均匀且剖面位置并未跨越NE向的昭通、莲峰断裂带等因素的限制(熊绍柏等, 1986, 1993王夫运等,2008),未能全面了解和深入认识川滇交界东部地区速度结构的横向变化特征.早期一些地震学者先后利用天然地震层析成像方法在四川、云南地区已取得一系列的研究成果(孙若昧等, 1991; 王椿镛等, 1994, 1999, 2002; Huang et al., 2002; Yao et al., 2006, 2008; Lei and Zhao, 2009, 2013, 2016; Wang et al., 2010; 易桂喜等,2010b王志等, 2017; 朱介寿等,2017),这些研究成果对揭示川滇地区地壳上地幔的构造环境及其动力学机制都起到了非常积极的推动作用,但由于所使用到的川滇交界东部地震台站分布较为稀疏(如昭通、莲峰断裂带及邻区仅有几个地震台站),且地震射线覆盖程度不高,成像结果未能充分剖析和反映昭通、莲峰断裂带及其邻区的深部结构特征和孕震环境.需要提及的是,2006年中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°N—32°N,100°E—105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵——川西台阵就覆盖了川滇交界东部地区的昭通、莲峰断裂带地区,研究人员先后利用川西台阵观测到的大量波形数据开展了接收函数反演、环境噪声数据高分辨率地壳速度结构成像等工作(刘启元等, 2008, 2009李昱等,2010Liu et al., 2014),川西台阵主要关注了松潘甘孜地块、四川盆地以及川滇块体的深部结构特征及各向异性,并首次证实青藏高原通过下部地壳物质流动和上部地壳沿断层块体滑移两种方式向东扩张,尚未对川滇交界东部地区的昭通、莲峰断裂带及其周边地区深部构造背景及其构造涵义进行过详细的讨论.由于昭通、莲峰断裂带及周边长期缺失详细的深部地球物理场响应资料,这也就给准确探讨该区的地震发震构造背景、运动学特征及其孕震机理带来较大困难.

鉴于此,本文充分收集了四川、云南等区域数字地震台网和“中国地震科学台阵探测—南北地震带南段”(“喜马拉雅”项目Ⅰ期)流动地震台阵等共计634个台站(阵)的观测数据,应用区域震和远震联合反演的方法得到了川滇地区速度结构特征,在此基础上,重点剖析和研究了昭通、莲峰断裂带及邻区的三维P波速度结构;再对昭通、莲峰断裂带及周边区域的重力、航磁数据进行三维视密度、视磁化强度反演,得到了壳内不同深度层视密度的横向变化特征和反映壳内磁性物质的分布范围以及结晶基底的视磁化强度异常展布情况.最后,我们综合分析研究了昭通、莲峰断裂带的深部构造环境及其动力学背景、探求鲁甸MS6.5地震的成因机制.该研究成果对于理解昭通、莲峰断裂带地震孕育的深部地球动力学机制、讨论鲁甸MS6.5的孕震构造背景以及评估未来地震活动趋势都具有重要的科学意义.

1 地震观测资料与成像方法

近年来,随着四川和云南区域数字地震观测台网的建设以及南北地震带南段布设的密集宽频带大型流动地震台阵,为研究川滇地区三维速度结构、介质物性特征和强震孕育的深部动力学模型等问题奠定了丰富的波形数据基础.为了获取高质量的地震观测数据以及避免仪器极性和方位角偏差等问题,我们主要选取了在区域测震地震计方位角普查工作后的地震观测资料,即对于四川和云南区域数字波形台网的数据观测时段均从2009年1月开始.本文使用的资料来自于四川省和云南省及其周边地区所布设的区域数字地震台网和中国地震科学台阵共计634个台(台站的分布如图 2a所示)所记录到的18530个近震事件和754个远震事件.首先,对近震的波形数据进行了到时提取,对于来自于区域数字地震台站的近震数据,我们选取震相观测报告中震级大于ML2.0的地震事件17734个,对来自于中国地震科学台阵记录的近震数据,选取震级大于ML3.0的地震事件796个,每个地震的P波到时观测数据均不少于15个.从图 2还可以看出,龙门山断裂带、木里—盐源弧形构造带以及本文所重点关注的川滇交界东段的昭通、莲峰断裂带等区域均有较好的地震事件分布.其次,在选取远震事件时,为保证地震资料的信噪比和可靠性,我们所选地震事件的震级均≥5.0,以确保地震波到达台站时具有较强的能量和清晰的初动;同时要求远震事件的震中距30°≤Δ≤90°,以确保地震波射线路径主要集中在地壳和上地幔;除此之外,要求每个远震事件被20个以上的地震台站(阵)记录到.根据以上原则,我们对远震事件进行分析和筛选,最终筛选得到的远震事件共计754个,可以看出成像所用的远震事件具有较好的方位角分布(图 2b).

图 2 本研究所用的地震台站(a)和地震事件分布图(b) Fig. 2 Distribution of seismic stations (a) and earthquake events (b) used in this study

本文采用了Zhao等(1994)提出的体波走时层析成像方法来反演川滇地区三维P波速度结构,在此基础上重点剖析和研究川滇交界东部地区的昭通、莲峰断裂带及其周边地区的速度结构特征.该方法允许速度在三维空间内任意变化,并通过在模型空间中设置一系列的三维网格节点,节点处的速度扰动作为反演中的未知数被求解,而模型中其他任意点的速度扰动可由与之相邻的8个节点的速度扰动线性插值得到.在射线追踪过程中对Um和Thurber (Um and Thurber, 1987)提出的近似弯曲算法进行了改进,迭代地应用伪弯曲技术和斯奈尔定律进行三维射线跟踪,使之适用于复杂的速度间断面存在的情况,在反演过程中,采用带阻尼因子的LSQR方法(Paige and Saunders, 1982)求解大型稀疏的观测方程组,且阻尼满足了模型和数据方差均为最小.我们对反演结果分别采用检测板测试(Zhao et al., 1994丁志峰,1999)和地震射线密度等方式进行评估,测试步骤是先建立三维空间网格点,在一维速度模型基础上加上正负相间的速度扰动(本文采用扰动值为±3%),然后计算与实际观测资料射线路径分布相同的理论走时,将生成的理论走时进行反演,通过对反演结果中每一个节点上的速度值与理论值的对比,可以获得观测数据成像的分辨能力.检测板结果显示研究区不同深度层的扰动幅度恢复较好,尤其是本研究区内部活动块体边界断裂及其周边等区域的测试结果较为理想(Li et al., 2018),这主要是由于大型边界断裂构造附近架设了较为密集的台站(阵)以及大量地震事件集中分布造成的.结合地震射线密度结果,同样发现本文所重点关注的川滇交界东段的昭通、莲峰断裂带与龙门山断裂带和怒江断裂南段均位于射线密度较高的区域.

2 P波速度结构特征

由于2014年鲁甸MS6.5主震发生在川滇块体东侧的大凉山次级块体南东缘昭通、莲峰断裂带内NW向的苞谷脑—小河断裂附近,因此截取了川滇交界东部地区的昭通、莲峰断裂带及其周边地区的三维P波速度结构进行重点剖析和研究.由于反演结果空间分辨率受到地震台站密度、事件数量以及射线覆盖密度的影响,在模型网格划分的时候,水平方向采用0.5°×0.5°均匀划分网格,在阻尼系数的选取时,通过模型变化和数据走时方差的折衷关系曲线来确定.分别用不同的阻尼值进行了多次反演试算,最后确定了反演阻尼系数为15.分辨率测试结果也表明研究区P波速度结构的空间分辨能力可以达到0.5°,特别是昭通断裂带鲁甸MS6.5地震震源区附近.图 3给出了研究区不同深度范围三维P波速度异常分布图,在浅部上地壳深度范围内,P波速度异常分布特征与地表地质构造、地形地貌密切相关,由1 km速度结构分布图可以看出马边构造带和金阳—峨边断裂所夹持的地区表现出高速异常特征,这主要是因为该处新构造位置处于凉山中升区,属于青藏高原东缘向华南褶皱系的过渡地带,故表现出了与第四纪以来强烈隆升相关的高速异常.位于研究区北东侧的华蓥山断裂南段和马边断裂东南段交汇处的宜宾、盐津等地也表现出与区内间歇性大面积整体抬升相关的高速异常.莲峰断裂带作为一条区域性的大断裂,对地层发育与区域构造变形起控制作用,莲峰断裂西侧的布拖盆地呈现出圈闭状的低速异常;昭通断裂控制了昭通、鲁甸新生代盆地的发育,因此,昭通、鲁甸地区的速度结构均表现出与地形相关的低速异常特征,该低速异常的边界受到昭通断裂构造所控制,其延伸方向与断裂走向基本一致.在10 km深度层上,昭通、莲峰断裂带及其周边上地壳速度结构呈现出明显的横向不均匀分布,并形成了尺度不同、高低速相间的分块结构,其中低速异常主要分布在大凉山断裂南段以及小江断裂的北段区域,这与吴建平(2013)的三维P波速度结构成果和徐涛等(2014, 2015)利用丽江—清镇650 km长深地震测深剖面,通过有限差分反演揭示该地区地壳速度结构研究成果相一致.昭通、莲峰断裂带及其周边上地壳物质存在显著的横向介质差异,昭通、莲峰断裂带处于高波速异常区,其间存在局部的相对低速异常体,2014年鲁甸MS6.5地震位于该高低速异常的分界线附近略偏向高速体一侧,反L型的余震沿着高低速异常的分界线NW向聚集分布,这种特有的速度结构特征有利于应力在脆性的上地壳内积累和集中.

图 3 三维P波速度异常分布 图中黑线为区内主要断裂,白色星号代表鲁甸MS6.5主震. Fig. 3 The images of 3D P wave velocity perturbation Black lines denote main faults in the study region, White star denotes the Ludian MS6.5 mainshock.

随着反演深度的增加,研究区中下地壳速度异常分布特征呈现一定的趋势性变化,其中15 km深度层上大凉山次级块体西南缘的低速异常分布特征明显,赵国泽等(2008)通过对青藏高原东缘及附近地区石棉—乐山剖面大地电磁资料的研究揭示了大凉山次级块体在约15 km深度存在“管道流”(低阻层),这与我们所揭示出的P波低速层展布范围相符.通过三维P波速度结构还可以看出,小江断裂北段和昭通断裂南段则处于低速异常区内,鲁甸MS6.5地震震源区附近出现大范围低速异常分布,这与李永华等(2014)的Rayleigh波群速度频散和远震接收函数联合反演所得到的S波速度结构结果相一致.30 km深度处的昭通、莲峰断裂带部分区域的速度异常分布形态和范围均有所改变,大凉山次级块体西南缘的低速异常仍然存在,并呈现出近SN向的展布特征.莲峰断裂北段的高速异常特征尤为明显,且正异常前缘向南扩展至苞谷脑—小河断裂附近,向东延伸至鲁甸、昭通和彝良一带,并被昭通断裂北段限制住了其继续东扩的范围.40 km和50 km深度处的三维P波速度结构反映了莲峰、昭通断裂带下地壳的速度结构特征,昭通断裂北段作为高低速异常分界线特征愈发明显,滇东北块体分布大范围的低速异常,表明该区下地壳物质相对较为软弱.Huang等(2015, EPSL)利用ChinArray地震台阵数据研究青藏高原东南缘不同深度P波各项异性结果也表明了该区40 km处低速异常的存在,并提出了稳定的扬子克拉通块体深部受到青藏高原东南缘活动块体的作用而遭到破坏的观点.丽江—新市镇的人工地震测深剖面的结果同样揭示了该区域地壳为高、低速相间的多层层状结构,上地壳和上地壳底部低速层微向西倾,显塑性的低速层埋于壳下20 km附近(崔作舟等, 1987).结合已有的冕宁—宜宾大地电磁测深剖面资料(万战生等, 2010)综合分析表明该低速高导异常与壳内物质的部分熔融或者壳内流体的存在有关.我们的P波成像结果还揭示了鲁甸MS6.5地震震源体下方中下地壳存在大范围低速异常分布,相对于中下地壳分布的低速异常而言,鲁甸MS6.5地震的震源体则处在坚硬脆性的上地壳介质范围内.

3 视密度、视磁化强度反演

基于重磁位场的分解与反演成像的异常特征研究可以为发震构造的深部物性特征与孕震环境研究提供重磁场响应依据,因此,针对特定的发震构造或潜在强震危险区,开展多方法、多层次的位场数据处理和异常特征提取研究,实现分解出的重磁特征能够反映发震断裂的深部结构和孕震背景显得尤为重要.然而,迄今为止,仍未见对昭通、莲峰断裂带及其邻区的重磁场分布特征进行过细致地分析研究,从重磁场响应特征来揭示昭通、莲峰断裂带深部构造特征、动力背景及其孕震环境仍需要深入探究.因此,为了揭示昭通、莲峰断裂带壳内不同深度范围内密度的横向变化情况和磁化强度分布差异特征,分析与地震孕育和发生过程密切相关的地壳结构和孕震构造特征以及块体内部的介质特性,同时也为了进一步验证三维P波成像结果的可靠性与合理性,需要我们提取出不同深度层次的重磁位场异常所含有的深部重磁异常响应信息,进而系统地分析昭通、莲峰断裂带的深部孕震环境和介质物性分布特征与2014年鲁甸MS6.5地震之间存在的关系.

本文所采用的视密度反演是一种基于位场分离和延拓的三维视密度反演方法,该方法已在不同地区的重力资料处理和研究工作中均取得了很好的效果(徐世浙,2007徐世浙等,2009杨金玉等,2008).首先基于场分离的切割法对观测到的地面重力场进行不同深度层源的切割分离,其步骤主要是实现对位场曲面上的局部凸起进行多次切割,达到对区域场和局部场进行分离的目的;由于位场数据处理中常用的向下延拓深度有限,通常只能达到3~5倍的点距(Pawlowski, 1995Fedi and Florio, 2002Cooper, 2004),为此,我们又采用了由徐世浙(2007)提出的基于位场大深度向下延拓的迭代法,将各层的场延拓至相应的深度后再反演得到各深度层的视密度分布情况.

图 4显示了不同深度层的视密度反演结果,从5 km深度图可知攀西构造带轴部的西昌—米易—攀枝花地带表现出明显的相对视密度正值异常,并随着深度的增加,其轴部高值正异常分布的特征依旧存在,这反映了轴部地带的地壳内部存在着来自深部的基性和超基性岩的高密度侵入体,而大凉山次级块体内部的美姑—金阳陷褶束出现了明显的条带状视密度低异常区.从10 km深度图上可以看出攀西构造带正值异常依然存在,而布拖盆地表现出串珠状圈闭低值异常区,这与前面三维P波速度结构所揭示出莲峰断裂西侧的布拖盆地呈现出圈闭状的低速异常分布相一致;随着反演深度的增加,研究区中下地壳视密度异常分布特征呈现一定的趋势性变化,其中15 km深度层的视密度反演结果显示了昭通、莲峰断裂带及其周边区域视密度变化异常等值线在昭通断裂附近形成了宽缓的梯级带.从15 km深度图中还可以看出凉山次级块体内部的美姑—金阳陷褶束和滇东块体内部会泽—曲靖台褶束均显示出明显的串珠状视密度低的异常区,西昌—米易—攀枝花等地区则显示出条带状相连的正视密度异常,昭通断裂带正好位于正负异常的分界线上.从新构造分区来看,由于位于大凉山断裂以东、荥经—马边—盐津断裂带以西的凉山中升区第四纪以来表现为整体性的大面积抬升,抬升幅度在2000 m左右,滇东—黔西掀升区整体性好,抬升幅度在1500 m左右,区内地势总体仍是北高南低,西高东低,并向东南掀斜,昭通断裂正处于这两个新构造分区的边界带附近,其特有的深部构造形态与该区的构造变形和地震的孕育及发生密切相关.随着反演深度的增加(30~50 km),昭通、莲峰断裂带及其周边中下地壳视密度变化区域愈发明显,川东南宜宾—自贡地区以及攀西地区的视密度正异常范围进一步扩大,两区中间夹持着负视密度异常区.30 km深度处滇东—黔西中等掀升区在整体上表现出了大范围的低值异常分布特征,并形成了明显的宣威低视密度中心带,这一趋势性变化在40 km、50 km深度层视密度反演图上则更为明显,这与之前三维P波速度结构所揭示出速度异常分布范围相对应.我们的反演结果也与石磊等(2015)研究结果相一致,两者都是以则木河断裂带为界,以东表现为低密度特征.Zhang等(2013)和Deng等(2016)的丽江—贵阳重力剖面的研究结果以人工地震资料为约束进行重力反演得到剖面实际密度也同样表明从地表平均到30 km、40 km和50 km深度层同样显示为低密度特征,这与我们本文研究结果认识也是相一致的.因此,基于P波速度结构和三维视密度反演结果均表明了昭通、莲峰断裂带中下地壳深度范围内存在低速、低密度的异常分布,进而说明了鲁甸MS6.5地震震源体处在坚硬脆性的上地壳介质内部.

图 4 不同深度的三维视密度反演图 Fig. 4 The figures of 3D apparent density inversion at different depth

接着,我们对研究区1:5万高精度航磁数据进行分析研究,从航磁ΔT异常等值线平面图可以看出昭通、莲峰断裂及周边地区的异常特征以串珠圈闭状为主(图 5),由于航磁ΔT数据反映的是不同深度磁性场源体(构造)综合叠加效应,故与孕震背景密切相关的区域性深部场源体的异常分布形态并不明显.为了压制浅部磁性场源体产生的局部异常或高频干扰信息、消除地磁场倾斜磁化对航磁异常所造成的干扰,提取出不同深度的航磁异常信息,为系统地分析强震危险区内的深部孕震环境、介质磁性分布特征与强震活动关系之间的关系提供了可靠的深部地球物理场依据,我们又进行了三维磁化强度的反演计算工作,具体步骤是先对原始磁异常数据进行相关预处理和化极处理,并结合本研究区的地质情况对化极磁力异常数据进行不同高度向上延拓处理的对比分析,选取合理延拓高度的磁异常作为磁化强度反演计算的基础数据,最后进行三维反演得到不同深度处的视磁化强度反演图(图 6).据Chen等(2015)Xu等(2015)对峨眉山大成火成岩研究表明其内带地幔大程度的部分熔融会形成大量的岩浆,除一部分喷出地表外,另一部分以侵入岩和底侵岩浆的形式停留在地壳的不同深度,故这些基性岩浆的在不同深度的航磁异常反演分布图中均会有所体现.

图 5 航磁ΔT异常平面等值线图 Fig. 5 The aeromagnetic ΔT anomalous contour map
图 6 三维视磁化强度反演结果 Fig. 6 The results of apparent magnetization inversion

图 6a为6 km反演深度图,从图中可以看出西昌—米易—攀枝花地带存在一系列梯度变化较大的短轴状圈闭正异常,总体呈NS向串珠状连续排列,由于该地区从北向南出露了一系列以早元古界-晚太古界的穹隆体为主体的各类岩浆岩体,因此强磁性异常多与这些基性、超基性岩浆岩的侵入有关,并且在不同深度的视磁化强度图中均有很好的航磁场响应.除此之外,我们还发现了一条展布于凉山块体内部及周缘的磁性隆起带,该磁性带主要表现出正异常分布特征,该带主要分布于昭觉、布拖、巧家、会东、托古等地,地面出露有一系列由震旦纪结晶杂岩、震旦-寒武系组成的断块和晚古生代辉绿岩以及玄武岩.昭通、莲峰断裂带南、北两段视磁化强度特征差异明显,这一趋势在9 km深度处的视磁化强度变化图中得到较好地体现,串珠状异常展布范围进一步增大,深层视磁化强度反演结果主要反映上地壳深部磁性基底特征,同时也揭示出共轭断裂的深部构造形态(图 6b).NW向的苞谷脑—小河断裂是与NE向昭通—鲁甸断裂相配套的共轭断裂,由数条断续展布的断层组成,反演结果表明该高低磁异常的分界线与苞谷脑—小河断裂的深部展布形态相一致,是其深部构造形态的反映,2014年鲁甸MS6.5地震位于该高低磁异常的分界线附近略偏向强磁异常体一侧,反L型的余震沿着该磁异常的分界线NW向聚集分布.除此之外,昭通断裂北段(昭通—鲁甸段)位于上地壳强磁性异常区内部,而断裂南段则处于负磁异常区内.鲁甸地震重新定位后的余震序列也是沿着航磁异常突变带分布,说明了该处共轭断裂的深部构造形态对地震发生和空间展布具有明显的控制作用.随着反演深度的增加,研究区中下地壳深度图中降低的磁场强度背景反映了昭通、莲峰断裂带及其以东地区深部物质磁性较弱的特性,其中图 6c 15 km深度图上揭示了昭通、莲峰断裂带南段地区较为平滑的弱磁性基底背景场且区域视磁化强度由北向南逐渐降低,昭通、莲峰断裂带及其周边视磁化强度变化图的等值线形态比浅层磁异常分布更为宽缓(图 6(d—f)),昭通、莲峰断裂带及其周边地区上地壳物质存在的横向介质差异,随着大凉山次级块体的南东向运动,受到华南块体边界部分地段上地壳内部刚性磁性基底的强烈阻挡,从而加剧了昭通—鲁甸地区基底岩层的褶皱变形,并产生了强烈的应力积累,也正是由于块体内部或者块体之间基底性质存在的局部差异,强磁性、高波速坚硬介质发育的昭通断裂段鲁甸地区往往有利于上地壳内部出现应力相对集中,脆性上地壳内部中低强度区域在横向挤压的构造应力场作用下易于发生破裂,从而对鲁甸MS6.5地震的孕育和发生起到重要的作用.

4 讨论

(1) 川滇交界东部昭通、莲峰断裂带的上地壳速度结构呈现出不均匀性分布特征,其中低速异常主要分布在大凉山断裂南段区域,昭通、莲峰断裂带处于高波速异常区,其间存在局部的相对低速异常体,表明了昭通、莲峰断裂带及其周边地区上地壳物质存在显著的横向介质差异,2014年鲁甸MS6.5地震位于该高低速异常的分界线附近略偏向高速体一侧,这与视磁化强度反演结果所揭示出的鲁甸MS6.5地震的震中位于正负磁异常的分界线附近略偏向强磁异常一侧相互吻合.为了进一步揭示昭通、莲峰断裂带内部2014年鲁甸MS6.5地震发生的深部介质特征和地震分布之间的关系,利用四川、云南省固定地震台网观测数据和震后流动台站余震监测资料对2014年8月3日至2014年8月11日期间发生的、震级ML≥1.5的506个地震序列进行了重新定位(图 7a),获得了359个事件的重定位结果,其震中分布图显示鲁甸地震序列平面上呈现共轭型分布特征,即沿着近EW向和NW向呈条带状聚集展布.其中,优势展布方向为NW向的余震序列主要集中分布在昭通—鲁甸断裂的西侧,呈密集条带状分布特征(图 7b),共轭断裂特定的深部构造形态对地震发生和空间展布具有明显的控制作用,鲁甸MS6.5主震位于高低异常分界线附近,苞谷脑—小河断裂处于正负航磁异常突变带附近.除此之外,还分析统计了地震序列震源深度分布特征发现鲁甸地震序列的震源深度优势分布层位主要集中在3~15 km之间(图 7cd),80.3%的地震震源深度较浅,平均深度为仅9.32 km,说明绝大多数地震主要发生在脆性的上地壳范围内.同时,昭通、莲峰断裂带及周边的P波速度结构还揭示了鲁甸MS6.5主震震源体下方中下地壳深度存在大范围低速异常分布,这种特有的深部物性结构特征有利于应力在脆性的上地壳内积累和集中(易桂喜等, 2010b, 2011).

图 7 精定位前后的地震震中分布和沿A-A′、B-B′剖面震源深度分布图 Fig. 7 Distributions of epicenter before and after precise relocation and the A-A′、B-B′profile

(2) 基于重磁异常反演结果也可用于分析昭通、莲峰断裂带的深部孕震环境和介质物性分布特征与2014年鲁甸MS6.5地震之间存在的关系.三维视密度反演结果表明了昭通、莲峰断裂带中下地壳深度范围内存在低密度异常分布,进而说明了鲁甸MS6.5地震震源体处在坚硬脆性的上地壳介质内部.同时,研究区中下地壳深度图中降低的磁场强度背景反映了昭通、莲峰断裂带及其以东地区深部物质磁性较弱的特性,综合说明了该区中下地壳物质相对较为软弱,强磁性坚硬介质发育的昭通断裂段鲁甸地区往往有利于上地壳内部出现应力相对集中,脆性上地壳内部中低强度区域在横向挤压的构造应力场作用下易于发生破裂,从而对鲁甸MS6.5地震的孕育和发生起到重要的作用.

(3) 2014年鲁甸MS6.5地震所在的昭通、莲峰断裂带及与其交汇、切割的马边—盐津断裂均是6.0级左右中等强度地震的频发区域,据史料记载到的最大地震为发生在马边—盐津断裂与莲峰断裂交汇部位的1974年大关北7.1级地震和1216年3月马湖7.0级地震,主要发生在马边—盐津断裂上,由于马边断裂带连续性差,由一系列相互切割的NE向和NW向断层组合而成,有研究者称之为马边新生地震构造带(韩竹军等,2009).有史料记载以来,沿大凉山断裂带曾发生1480年9月越西5.5级、1881年6月越西5.0级和1944年6月昭觉51/4级地震、2014年10月越西M5.0地震(易桂喜等,2016),但没有记录到6.5级以上破坏型地震发生,且该断裂带地表迹线连续性差,也是一条形成时代较安宁河断裂要晚的新生断裂带(Xu et al., 2003),探槽开挖结果揭示出大凉山断裂带存在多次古地震事件(宋方敏等,2002何宏林等,2008).因此,驱动大凉山次级块体内部断裂及其南东缘的莲峰、昭通断裂构造带的构造变形和地震活动的深部动力来源是一个值得探究的问题.

赵国泽、万战生等(赵国泽等,2008万战生等,2010)在川滇块体东边界的大凉山次级块体内部也发现了高导层,并认为是青藏高原东边缘带向东南方向挤出作用下形成的“管流”层.我们的前期研究成果也表明了大凉山次级块体中下地壳深度范围内低速异常的优势展布方向为近SN向,与大凉山断裂的走向基本一致,在沿着EW方向上,川西北次级块体和大凉山块体在中下地壳深度范围内均存在连续的低速层分布,我们推测该低速层自西向东越过大凉山断裂,最终止于荥经—马边断裂构造带及其附近(Li et al., 2018).结合前人的震源机制研究结果可知,昭通、莲峰断裂带活动与变形的动力源应是直接来自大凉山次级块体的南东向运动,而间接来自川滇块体的南南东向运动(阚荣举等, 1977, 1983成尔林,1981张培震,2008闻学泽等,2013).以上研究结果与我们的P波成像结果所显示的低速管流层进入大凉山次级块体内部转向南东方向运动相一致.本文的三维视密度分布特征和视磁化强度反演结果共同揭示了昭通、莲峰断裂带中下地壳范围内存在大范围低密度、弱磁性异常分布,综合说明了该区中下地壳物质相对较为软弱.因此,我们认为发生在青藏高原东缘中下地壳物质的塑性流展为大凉山次级块体内部及其南东缘的莲峰、昭通断裂构造带的构造变形和地震活动提供了深部动力来源,当大凉山块体内部的中下地壳低速层自NW向SE方向运动到昭通、莲峰断裂带附近时,当大凉山块体内部的中下地壳低速管流层自NW向SE方向运动到昭通、莲峰断裂带附近时,受到华南块体的强烈阻挡,应力在莲峰、昭通断裂附近基底性质存在明显差异处集中,脆性上地壳中低强度的区域在横向挤压的构造应力场作用下易于破裂从而引发强震,这也正是昭通、莲峰断裂带地震孕育和发生的深部构造环境.

(4) 从区域构造部位和地理位置分析,龙门山断裂带是位于巴颜喀拉块体东缘松潘—甘孜块体的边界断裂带,2008年汶川MS8.0地震、2013年芦山MS7.0地震发生在龙门山断裂带上.而昭通、莲峰断裂带是位于大凉山次级块体南东缘的边界断裂带,2014年8月3日鲁甸MS6.5地震发生在昭通断裂带附近,除此之外,昭通、莲峰断裂带及其周边地区在历史上曾发生过多次强震,如1917年大关7¼级地震和1974年大关北7.1级地震.尤其是近些年来,昭通、莲峰断裂带不但中强地震活动频次增多、频度增大,而且该区中小地震活动也十分密集,空间分布呈现一定的不均匀性特征,绝大多数3.0级以上地震发生在昭通、莲峰断裂附近及其所夹持的地区,因此,该区地震危险背景和潜在发震能力已引起广泛关注.

据我们的野外地质调查工作发现,莲峰断裂和昭通断裂切过一系列山脊和高台地,并形成断层垭口、断层槽谷等地貌现象,追索莲峰断裂于永善县水竹乡附近形成一系列断层垭口和断层槽谷地貌,昭通断裂北段于邀集块、赵家垭口一带形成断层槽谷地貌,断层槽谷延伸数百米,甚至1~2 km长,偶见有泉眼出露.已有研究结果还表明了昭通断裂的鲁甸附近段存在异常低b值或高应力的断裂段,昭通、莲峰断裂带南、北两个段落业已不同程度闭锁,并且昭通断裂带具有更高的闭锁和应变积累程度(闻学泽等,2013).根据本文对川滇交界东部昭通、莲峰断裂带深部结构特征和孕震环境的综合研究成果可知,昭通、莲峰断裂带上地壳介质具有明显的横向不均匀性,中下地壳深度范围内存在低波速、低密度和弱磁性的异常分布.需要关注的是,昭通断裂北段(昭通—鲁甸段)位于上地壳强磁性、高波速异常区内且具有深大断裂的深部地球物理场响应特征,因此昭通断裂北段(昭通—鲁甸段)具备发生7.0级及以上强震的深部构造(孕震)背景.本文研究成果为2014年鲁甸MS6.5地震的深部介质条件和孕震机制提供了可靠的深部地球物理场依据,也将有助于研判该区地震活动趋势和深入理解驱动昭通、莲峰断裂带构造变形及演化的深部动力机制等科学问题.

5 结论

本文充分收集了四川、云南等区域数字地震台网和“中国地震科学台阵探测-南北地震带南段”(“喜马拉雅”项目Ⅰ期)流动地震台阵的观测数据,应用区域震和远震联合反演的方法重点剖析和研究了昭通、莲峰断裂带三维P波速度结构;同时,我们对昭通、莲峰断裂带及周边区域的重力、航磁数据进行了三维视密度、视磁化强度反演,得到了壳内不同深度层视密度的横向变化特征和反映壳内磁性物质的分布范围以及结晶基底的视磁化强度异常分布情况,在此基础上综合分析研究昭通、莲峰断裂带的深部介质性质、构造形态及深部孕震环境,探求2014年鲁甸MS6.5地震成因的深部动力机制,具体结论如下:

(1) 川滇交界东部昭通、莲峰断裂带及其周边地区上地壳物质存在显著的横向介质差异,中下地壳深度范围大凉山次级块体西南缘存在低速异常分布,并呈现出近SN向的展布特征,研究区正异常前缘向南扩展至苞谷脑—小河断裂附近,向东延伸至鲁甸、昭通和彝良一带,并被昭通断裂北段限制住了其继续东扩的范围,2014年鲁甸MS6.5地震位于该高低速异常的分界线附近略偏向高速体一侧,昭通、莲峰断裂带及周边绝大多数地震主要发生在脆性的上地壳范围内.

(2) 成像结果还揭示了鲁甸MS6.5主震震源体下方中下地壳存在大范围低速异常分布,三维P波速度结构与重磁异常的展布特征、梯度变化在深度和分区特征上均具有较好的联系和可比性,结合昭通、莲峰断裂带中下地壳范围内存在大范围的低密度、弱磁性异常分布,综合表明了该区中下地壳物质相对软弱,这种特有的深部物性结构特征有利于应力在脆性上地壳内积累和集中.

(3) 我们的研究结果还揭示了与NE向昭通—鲁甸断裂相配套的共轭断裂—NW向苞谷脑—小河断裂的深部构造形态,三维视磁化强度反演结果表明高低航磁异常的分界线与NW向苞谷脑—小河断裂的深部展布形态相一致,苞谷脑—小河断裂处于该正负航磁异常突变带附近,2014年鲁甸MS6.5地震位于该高低磁异常的分界线附近略偏向强磁异常一侧,昭通断裂北段(昭通—鲁甸段)位于上地壳强磁性、高波速异常区内且具有深大断裂的深部地球物理场响应特征,因此该断裂段(昭通—鲁甸段)具备发生7级及以上强震的深部构造(孕震)背景.

(4) 我们认为发生在青藏高原东缘中下地壳低速物质的塑性流展为大凉山次级块体内部及其南东缘的昭通、莲峰断裂带的构造变形和地震活动提供了深部动力来源,当大凉山次级块体内部的中下地壳低速管流层自NW向SE方向运动到昭通、莲峰断裂带附近时,受到华南块体的强烈阻挡,应力在莲峰—昭通断裂附近基底性质存在差异处集中,脆性上地壳中低强度的区域在横向挤压的构造应力场作用下易于破裂从而引发强震,这也正是昭通、莲峰断裂带地震孕育和发生的深部构造环境.

致谢  感谢日本东北大学赵大鹏教授提供的地震体波层析成像程序,感谢美国地质调查局(U.S. Geological Survey)地震科学中心Walter D Mooney教授和中国地震局地质研究所詹艳研究员对文章给予的指导和帮助,文章的撰写过程先后与杨歧焱博士、吴微微高级工程师进行了多次有益的交流和探讨,作者一并表示最衷心的感谢.
References
Chen S, Wang Q H, Wang Q S, et al. 2014. The 3D density structure and gravity change of Ludian MS6.5 Yunnan epicenter and surrounding regions. Chinese Journal of Geophysics (in Chinese), 57(9): 3080-3090. DOI:10.6038/cjg20140933
Chen Y, Xu YG, Xu T, et al. 2015. Magmatic underplating and crustal growth in the Emeishan Large Igneous Province, SW China, revealed by a passive seismic experiment. Earth Planet. Sci. Lett, 432: 103-114. DOI:10.1016/j.epsl.2015.09.048
Cheng E L. 1981. Recent tectonic stress field and tectonic movement of the Sichuan province and its vicinity. Acta Seimologica Sinica (in Chinese), 3(3): 231-241.
Cooper G. 2004. The stable downward continuation of potential field data. Exploration Geophysics, 35(4): 260-265. DOI:10.1071/EG04260
Cui Z Z, Lu D Y, Chen J P, et al. 1987. The deep structure and tectonic features of the crust in Panxi area. Acta Geophysica Sinica (in Chinese), 30(6): 566-580.
Deng Y, Chen Y, Wang P, et al. 2016. Magmatic underplating beneath the Emeishan large igneous province (South China) revealed by the COMGRAELIP experiment. Tectonophysics, 672-673: 16-23.
Ding Z F. 1999. The theory and application of local earthquake tomography [Ph. D. thesis]. Beijing: Institute of Geophysics, China Earthquake Administration.
Fedi M, Florio G. 2002. A stable downward continuation by using the ISVD method. Geophysical Journal International, 151(1): 146-156. DOI:10.1046/j.1365-246X.2002.01767.x
Han Z J, He Y L, An Y F, et al. 2009. A new seismotectonic belt:features of the latest structural deformation style in the Mabian seismotectonic zone. Acta Geologica Sinica (in Chinese), 83(2): 218-229.
He H L, Ikeda Y, He Y L, et al. 2008. Newly-generated Daliangshan fault zone-Shortcutting on the central section of Xianshuihe-Xiaojiang fault system. Science in China Series D:Earth Sciences, 51(9): 1248-1258. DOI:10.1007/s11430-008-0094-4
Huang J L, Zhao D P, Zheng S H. 2002. Lithospheric structure and its relationship to seismic and volcanic activity in southwest China. Journal of Geophysical Research:Solid Earth, 107(B10): ESE 13-1-ESE 13-14. DOI:10.1029/2000JB000137
Huang Z C, Wang P, Xu M J, et al. 2015. Mantle structure and dynamics beneath SE Tibet revealed by new seismic images. Earth and Planetary Science Letters, 411: 100-111. DOI:10.1016/j.epsl.2014.11.040
Kan R J, Zhang S C, Yan F T, et al. 1977. Present tectonic stress field and its relation to the characteristics of recent tectonic activity in southwestern China. Acta Geophysica Sinica (in Chinese), 20(2): 96-109.
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 S, Zhao D P. 2009. Structural heterogeneity of the Longmenshan fault zone and the mechanism of the 2008 Wenchuan earthquake (MS8.0). Geochemistry, Geophysics, Geosystems, 10(10): Q10010. DOI:10.1029/2009GC002590
Lei J S, Zhao D P, Su J R, et al. 2009. Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake. Chinese Journal of Geophysics (in Chinese), 52(2): 339-345.
Lei J S, Zhao D P. 2013. Structural heterogeneity of the Longmenshan fault zone and the mechanism of the 2008 Wenchuan earthquake (MS8.0). Geochemistry, Geophysics, Geosystems, 10(10): Q10010. DOI:10.1029/2009GC002590
Lei J S, Zhao D P. 2016. Teleseismic P-wave tomography and mantle dynamics beneath Eastern Tibet. Geochemistry, Geophysics, Geosystems, 17(5): 1861-1884. DOI:10.1002/2016GC006262
Li D H, Liao H, Ding Z F, et al. 2018. Joint Inversion of the 3D P wave velocity structure of the crust and upper mantle under the southeastern margin of the Tibetan plateau using regional earthquake and teleseismic data. Acta Geologica Sinica (English Edition), 92(1): 16-33. DOI:10.1111/1755-6724.13492
Li Y, Yao H J, Liu Q Y, et al. 2010. Phase velocity array tomography of Rayleigh waves in western Sichuan from ambient seismic noise. Chinese Journal of Geophysics (in Chinese), 53(4): 842-852. DOI:10.3969/j.issn.0001-5733.2010.04.009
Li Y H, Xu X M, Zhang E H, et fal. 2014. Three-dimensional crust structure beneath SE Tibetan plateau and its seismotectonic implications for the Ludian and Jinggu earthquake. Seismology and Geology (in Chinese), 36(4): 1204-1216.
Liu Q Y, Chen J H, Li S C, et al. 2008. The MS8.0 Wenchuan earthquake:Preliminary results from the western Sichuan mobile seismic array observations. Seismology and Geology (in Chinese), 30(3): 584-596.
Liu Q Y, Li Y, Chen J H, et al. 2009. Wenchuan MS8.0 earthquake:preliminary study of the S-wave velocity structure of the crust and upper mantle. Chinese Journal of Geophysics (in Chinese), 52(2): 309-319.
Liu Q Y, Van Der Hilst R D, Li Y, et al. 2014. Eastward expansion of the Tibetan Plateau by crustal flow and strain partitioning across faults. Nature Geoscience, 7(5): 361-365. DOI:10.1038/ngeo2130
Paige C C, Saunders M A. 1982. LSQR:An algorithm for sparse linear equations and sparse least squares. ACM Transactions on Mathematical Software, 8(1): 43-71. DOI:10.1145/355984.355989
Pawlowski R S. 1995. Preferential continuation for potential-field anomaly enhancement. Geophysics, 60(2): 390-398. DOI:10.1190/1.1443775
Song F M, Li R C, Xu X W. 2002. Preliminary results of the investigation of paleo-earthquakes along the Daliangshan fault, zone, Sichuan province, China. Seismology and Geology (in Chinese), 24(1): 27-34.
Sun R M, Liu F T, Liu J H. 1991. Seismic tomography of Sichuan. Acta Geophysica Sinica (in Chinese), 34(6): 708-716.
Um J. Thurber C. 1987. A fast algorithm for two-point seismic ray tracing. Bulletin of the Seismological Society of America, 77(3): 972-986.
Wan Z S, Zhao G Z, Tang J, et al. 2010. The electrical structure of the crust along Mianning-Yibin profile in the eastern edge of Tibetan plateau and its tectonic implications. Chinese Journal of Geophysics (in Chinese), 53(3): 585-594. DOI:10.3969/j.issn.0001-5733.2010.03.012
Wang C Y, Wang X L, Yan Q Z. 1994. Three-dimensional velocity structure below the Kunming seismic network. Acta Seismologica Sinica (in Chinese), 16(2): 167-175.
Wang C Y, Wu J P, Lou H. 1999. Research on relation of deep structure and strong earthquake activity in Sichuan-Yunnan region. Seismological and Geomagnetic Observation and Research (in Chinese), 20(5): 80-87.
Wang C Y, Mooney W D, Wang X L, et al. 2002. Study on 3-D velocity structure of crust and upper mantle in Sichuan-Yunnan region, China. Acta Seismologica Sinica (in Chinese), 24(1): 1-16.
Wang F Y, Duan Y H, Yang Z X, et al. 2008. Velocity structure and active fault of Yanyuan-Mabian seismic zone-The result of high-resolution seismic refraction experiment. Science in China Series D:Earth Sciences, 51(9): 1284-1296. DOI:10.1007/s11430-008-0098-0
Wang Z, Zhao D P, Wang J. 2010. Deep structure and seismogenesis of the north-south seismic zone in southwest China. Journal of Geophysical Research:Solid Earth, 115(B12): B12334. DOI:10.1029/2010JB007797
Wang Z, Wang X B, Huang R Q, et al. 2017. Deep structure imaging of multi-geophysical parameters and seismogenesis in the Longmenshan fault zone. Chinese Journal of Geophysics (in Chinese), 60(6): 2068-2079. DOI:10.6038/cjg20170604
Wen X Z, Fan J, Yi G X, et al. 2008. A seismic gap on the Anninghe fault in western Sichuan, China. Science in China Series D:Earth Sciences, 51(10): 1375-1387. DOI:10.1007/s11430-008-0114-4
Wen X Z, Du F, Yi G X, et al. 2013. Earthquake potential of the Zhaotong and Lianfeng fault zones of the eastern Sichuan-Yunnan border region. Chinese Journal of Geophysics (in Chinese), 56(10): 3361-3372. DOI:10.6038/cjg20131012
Wu J P, Yang T, Wang W L, et al. 2013. Three dimensional P-wave velocity structure around Xiaojiang fault system and its tectonic implications. Chinese Journal of Geophysics (in Chinese), 56(7): 2257-2267. DOI:10.6038/cjg20130713
Xiong S B, Teng J W, Yin Z W, et al. 1986. Explosion seismological study of the structure of the crust and upper mantle at southern part of the Panxi tectonic belt. Acta Geophysica Sinica (in Chinese), 29(3): 233-244.
Xiong S B, Zheng Y, Yin Z X, et al. 1993. The 2-D structure and it's tectonic implications of the crust in the Lijiang-Pan Zhihua-Zhehai region. Acta Geophysica Sinica (in Chinese), 36(4): 434-444.
Xu S Z. 2007. A comparison of effects between the iteration method and FFT for downward continuation of potential fields. Chinese Journal of Geophysics (in Chinese), 50(1): 285-289.
Xu S Z, Yu H L, Li H X, et al. 2009. The inversion of apparent density based on the separation and continuation of potential field. Chinese Journal of Geophysics (in Chinese), 52(6): 1592-1598. DOI:10.3969/j.issn.0001-5733.2009.06.021
Xu T, Zhang Z J, Liu B F, et al. 2015. Crustal velocity structure in the Emeishan Large Igneous Province and evidence of the Permian mantle plume activity. Science China:Earth Sciences, 58(7): 1133-1147. DOI:10.1007/s11430-015-5094-6
Xu T, Zhang M H, Tian X B, et al. 2014. Upper crustal velocity of Lijiang-Qingzhen profile and its relationship with the seismogenic environment of the MS6.5 Ludian earthquake. Chinese Journal of Geophysics, Science China:Earth Sciences, 57(9): 3069-3079.
Xu T, Zhang Z J, Liu B F, et al. 2015. Crustal velocity structure in the Emeishan large igneous province and evidence of the Permian mantle plume activity. Constraints of wide-angle seismic data from Lijiang-Qingzhen, 45(5): 561-576.
Xu X W, Wen X Z, Zheng R Z, et al. 2003. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region, China. Science in China Series D:Earth Sciences, 46(S2): 210-226.
Xu X W, Jiang G Y, Yu G H, et al. 2014. Discussion on seismogenic fault of the Ludian MS6.5 earthquake and its tectonic attribution. Chinese Journal of Geophysics (in Chinese), 57(9): 3060-3068. DOI:10.6038/cjg20140931
Yang J Y, Xu S Z, Yu H L, et al. 2008. Application of apparent density inversion method in the East China Sea and its adjacent area. Chinese Journal of Geophysics (in Chinese), 51(6): 1909-1916.
Yao H, Van Der Hilst R D, De Hoop M V. 2006. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-Ⅰ Phase velocity maps. Geophysical Journal International, 166(2): 732-744. DOI:10.1111/j.1365-246X.2006.03028.x
Yao H J, Beghein C, Van Der Hilst R D. 2008. Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-Ⅱ.Crustal and upper-mantle structure. Geophysical Journal International, 173(1): 205-219. DOI:10.1111/j.1365-246X.2007.03696.x
Yi G X, Wen X Z, Su Y J. 2008. Study on the potential strong-earthquake risk for the eastern boundary of the Sichuan-Yunnan active faulted-block, China. Chinese Journal of Geophysics (in Chinese), 51(6): 1719-1725.
Yi G X, Wen X Z, Zhang Z W, et al. 2010a. Study on potential strong earthquake risk in Mabian area, southern Sichuan. Seismology and Geology (in Chinese), 32(2): 282-293.
Yi G X, Yao H J, Zhu J S, et al. 2010b. Lithospheric deformation of continental China from Rayleigh wave azimuthal anisotropy. Chinese Journal of Geophysics (in Chinese), 53(2): 256-268. DOI:10.3969/j.issn.0001-5733.2010.02.004
Yi G X, Wen X Z, Xin H, et al. 2011. Distributions of seismicity parameters and seismic apparent stresses on the Longmenshan-Minshan tectonic zone before the 2008 MS8.0 Wenchuan earthquake. Chinese Journal of Geophysics (in Chinese), 54(6): 1490-1500. DOI:10.3969/j.issn.0001-5733.2011.06.008
Yi G X, Long F, Zhao M, et al. 2016. Focal mechanism and seismogenic structure of the M5.0 Yuexi earthquake on 1 Oct. 2014, southwestern China. Seismology and Geology (in Chinese), 38(4): 1124-1136.
Zhang P Z. 2008. The tectonic deformation, strain distribution and deep dynamic processes in the eastern margin of the Qinghai-Tibetan plateau. Science in China (Series D) (in Chinese), 38(9): 1041-1056.
Zhang X, Wang P, Xu T, et al. 2018. Density structure of the crust in the Emeishan large igneous province revealed by the Lijiang-Guiyang gravity profile. Earth and Planetary Physics, 2: 74-81.
Zhao D P, Hasegawa A, Kanamori H. 1994. Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events. Journal of Geophysical Research:Solid Earth, 99(B11): 22313-22329. DOI:10.1029/94JB01149
Zhao G Z, Chen X B, Wang L F, et al. 2008. Evidence of crustal 'channel flow' in the eastern margin of Tibetan Plateau from MT measurements. Chinese Science Bulletin, 53(12): 1887-1893.
Zhu J S, Wang X B, Yang Y H, et al. 2017. The crustal flow beneath the eastern margin of the Tibetan Plateau and its process of dynamics. Chinese Journal of Geophysics (in Chinese), 60(6): 2038-2057. DOI:10.6038/cjg20170602
陈石, 王青华, 王谦身, 等. 2014. 云南鲁甸MS6.5地震震源区和周边三维密度结构及重力场变化. 地球物理学报, 57(9): 3080-3090. DOI:10.6038/cjg20140933
成尔林. 1981. 四川及其邻区现代构造应力场和现代构造运动特征. 地震学报, 3(3): 231-241.
崔作舟, 卢德源, 陈纪平, 等. 1987. 攀西地区的深部地壳结构与构造. 地球物理学报, 30(2): 566-580.
丁志峰. 1999.近震层析成像的理论与应用[博士论文].北京: 中国地震局地球物理研究所.
韩竹军, 何玉林, 安艳芬, 等. 2009. 新生地震构造带-马边地震构造带最新构造变形样式的初步研究. 地质学报, 83(2): 218-229.
何宏林, 池田安隆, 何玉林, 等. 2008. 新生的大凉山断裂带-鲜水河-小江断裂系中段的裁弯取直. 中国科学D辑:地球科学, 38(5): 564-574.
阚荣举, 张四昌, 晏凤桐, 等. 1977. 我国西南地区现代构造应力场与现代构造活动特征的探讨. 地球物理学报, 20(2): 96-109.
阚荣举, 王绍晋, 黄崐, 等. 1983. 中国西南地区现代构造应力场与板内断块相对运动. 地震地质, 5(2): 79-90.
李昱, 姚华建, 刘启元, 等. 2010. 川西地区台阵环境噪声瑞利波相速度层析成像. 地球物理学报, 53(4): 842-852. DOI:10.3969/j.issn.0001-5733.2010.04.009
李永华, 徐小明, 张恩会, 等. 2014. 青藏高原东南缘地壳结构及云南鲁甸、景谷地震深部孕震环境. 地震地质, 36(4): 1204-1216. DOI:10.3969/j.issn.0253-4967.2014.04.021
刘启元, 陈九辉, 李顺成, 等. 2008. 汶川MS8.0地震:川西流动地震台阵观测数据的初步分析. 地震地质, 30(3): 584-596. DOI:10.3969/j.issn.0253-4967.2008.03.002
刘启元, 李昱, 陈九辉, 等. 2009. 汶川MS8.0地震:地壳上地幔S波速度结构的初步研究. 地球物理学报, 52(2): 309-319.
石磊, 楼海, 王谦身, 等. 2015. 攀西地区重力场特征及地壳密度结构. 地球物理学报, 58(7): 2402-2412.
宋方敏, 李如成, 徐锡伟. 2002. 四川大凉山断裂带古地震研究初步结果. 地震地质, 24(1): 27-34. DOI:10.3969/j.issn.0253-4967.2002.01.003
孙若昧, 刘福田, 刘建华. 1991. 四川地区的地震层析成像. 地球物理学报, 34(6): 708-716. DOI:10.3321/j.issn:0001-5733.1991.06.005
万战生, 赵国泽, 汤吉, 等. 2010. 青藏高原东边缘冕宁-宜宾剖面电性结构及其构造意义. 地球物理学报, 53(3): 585-594. DOI:10.3969/j.issn.0001-5733.2010.03.012
王椿铺, Mooney W D, 王溪莉, 等. 2002. 川滇地区地壳上地幔三维速度结构研究. 地震学报, 24(1): 1-16. DOI:10.3321/j.issn:0253-3782.2002.01.001
王椿镛, 王溪莉, 颜其中. 1994. 昆明地震台网下方的三维速度结构. 地震学报, 16(2): 167-175.
王椿镛, 吴建平, 楼海. 1999. 关于川滇地区深部结构与强震活动关系的研究. 地震地磁观测与研究, 20(5): 80-87. DOI:10.3969/j.issn.1003-3246.1999.05.013
王夫运, 段永红, 杨卓欣, 等. 2008. 川西盐源-马边地震带上地壳速度结构和活动断裂研究高分辨率地震折射实验结果. 中国科学D辑:地球科学, 38(5): 611-621.
王志, 王绪本, 黄润秋, 等. 2017. 龙门山断裂带多参数深部结构成像与地震成因研究. 地球物理学报, 60(6): 2068-2079. DOI:10.6038/cjg20170604
闻学泽, 范军, 易桂喜, 等. 2008. 川西安宁河断裂上的地震空区. 中国科学D辑:地球科学, 38(7): 797-807.
闻学泽, 杜方, 易桂喜, 等. 2013. 川滇交界东段昭通、莲峰断裂带的地震危险背景. 地球物理学报, 56(10): 3361-3372. DOI:10.6038/cjg20131012
吴建平, 杨婷, 王未来, 等. 2013. 小江断裂带周边地区三维P波速度结构及其构造意义. 地球物理学报, 56(7): 2257-2267. DOI:10.6038/cjg20130713
熊绍柏, 滕吉文, 尹周勋, 等. 1986. 攀西构造带南部地壳与上地幔结构的爆炸地震研究. 地球物理学报, 29(3): 233-244.
熊绍柏, 郑晔, 尹周勋, 等. 1993. 丽江-攀枝花-者海地带二维地壳结构及其构造意义. 地球物理学报, 36(4): 434-444. DOI:10.3321/j.issn:0001-5733.1993.04.004
徐世浙. 2007. 迭代法与FFT法位场向下延拓效果的比较. 地球物理学报, 50(1): 285-289. DOI:10.3321/j.issn:0001-5733.2007.01.035
徐世浙, 余海龙, 李海侠, 等. 2009. 基于位场分离与延拓的视密度反演. 地球物理学报, 52(6): 1592-1598. DOI:10.3969/j.issn.0001-5733.2009.06.021
徐涛, 张明辉, 田小波, 等. 2014. 丽江-清镇剖面上地壳速度结构及其与鲁甸MS6.5级地震孕震环境的关系. 地球物理学报, 57(9): 3069-3079.
徐涛, 张忠杰, 刘宝峰, 等. 2015. 峨眉山大火成岩省地壳速度结构与古地幔柱活动遗迹:来自丽江-清镇宽角地震资料的约束. 中国科学:地球科学, 45(5): 561-576.
徐锡伟, 江国焰, 于贵华, 等. 2014. 鲁甸6.5级地震发震断层判定及其构造属性讨论. 地球物理学报, 57(9): 3060-3068. DOI:10.6038/cjg20140931
杨金玉, 徐世浙, 余海龙, 等. 2008. 视密度反演在东海及邻区重力异常解释中的应用. 地球物理学报, 51(6): 1909-1916. DOI:10.3321/j.issn:0001-5733.2008.06.034
易桂喜, 闻学泽, 苏有锦. 2008. 川滇活动地块东边界强震危险性研究. 地球物理学报, 51(6): 1719-1725. DOI:10.3321/j.issn:0001-5733.2008.06.012
易桂喜, 闻学泽, 张致伟, 等. 2010a. 川南马边地区强震危险性分析. 地震地质, 32(2): 282-293.
易桂喜, 姚华建, 朱介寿, 等. 2010b. 用Rayleigh面波方位各向异性研究中国大陆岩石圈形变特征. 地球物理学报, 53(2): 256-268. DOI:10.3969/j.issn.0001-5733.2010.02.004
易桂喜, 闻学泽, 辛华, 等. 2011. 2008年汶川MS8.0地震前龙门山-岷山构造带的地震活动性参数与地震视应力分布. 地球物理学报, 54(6): 1490-1500. DOI:10.3969/j.issn.0001-5733.2011.06.008
易桂喜, 龙锋, 赵敏, 等. 2016. 2014年10月1日越西M5.0地震震源机制与发震构造分析. 地震地质, 38(4): 1124-1136. DOI:10.3969/j.issn.0253-4967.2016.04.025
张培震. 2008. 青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程. 中国科学D辑:地球科学, 38(9): 1041-1056.
赵国泽, 陈小斌, 王立风, 等. 2008. 青藏高原东边缘地壳"管流"层的电磁探测证据. 科学通报, 53(3): 345-350. DOI:10.3321/j.issn:0023-074X.2008.03.011
朱介寿, 王绪本, 杨宜海, 等. 2017. 青藏高原东缘的地壳流及动力过程. 地球物理学报, 60(6): 2038-2057. DOI:10.6038/cjg20170602