Structural characteristics of Yanyuan basin deduced from seismic survey and its dynamic implication
-
摘要: 基于覆盖盐源盆地的短周期天然地震台阵和布设的一条人工地震测线所获得的地震数据,从中提取地震能量属性,并通过地震层析成像获得该地区的浅部地震速度结构,继而对短周期地震台阵一个月的噪声数据进行互相关得到经验格林函数,再通过时频分析获得相速度频散曲线,反演获得不同深度的S波速度分布。研究结果显示,盐源盆地地震特征主要分为三层:浅部为低速新生代沉积地层,P波速度为1.0—2.1 km/s,反射轴明显、连续性好,反射能量强,S波速度为中高南低,南部低速体与盐源断裂走向一致,新生代盆地整体呈南深北浅,沉积发育和构造形态受盐源断裂控制;中部为中低速三叠系地层,P波速度为3.5—4.5 km/s,反射轴不连续,反射能量较弱,S波速度逐渐变大,整体变化变小;深部为高速古生代地层,地震反射较为杂乱,反射能量弱;上地壳存在滑脱面,该界面为沉积盖层与结晶基底的分界面,且向浅部发育一系列断层,其中金河—箐河断裂为盐源盆地与康滇地块的分界线,这些断裂带也是盐源盆地地震频发的部位。Abstract: In this paper, an artificial seismic line and a short period natural seismic array covering Yanyuan basin are set up. The seismic energy attribute is extracted from the artificial seismic data, and the seismic velocity structure is obtained by seismic tomography. Meanwhile the S-wave velocity distribution at different depths is obtained by inverting the empirical Green’s functions retrieved from cross-correlation of the ambient noise data recorded by the short period seismic array. The results show that Yanyuan basin is mainly divided into three layers by seismic characteristics. The shallow layer is the Cenozoic sedimentary stratum with low P wave velocity ranging from 1.0 km/s to 2.1 km/s, and it is characterized by good continuity and strong reflection. The S wave velocity is low in the south and high in the middle of Yanyuan basin, and the trend of the low-velocity anomaly in the south is in agreement with the strike of Yanyuan fault. Cenozoic basin is thick in the south and thin in the north, where the sedimentary and structure are controlled by Yanyuan fault. The middle layer is medium-low P wave velocity ranging from 3.5 km/s to 4.5 km/s, and it is characterized by discontinuous reflection axis and weak reflection, which can be interpreted as Triassic strata. The S-velocity becomes higher gradually and the change becomes smaller in general. The deep layer is high velocity, characterized by chaotic seismic reflection and weak reflection energy, interpreted as Paleozoic strata. Furthermone, there is a detachment surface in the upper crust, which is the interface between sedimentary and basement. A series of faults are developed from the detachment to the surface, and these faults are also an earthquake-prone area in Yanyuan basin. One of the faults is Jinhe-Qinghe fault, which is the boundary fault between Yanyuan basin and Kangdian block.
-
引言
四川地处青藏高原东南缘,位于我国南北地震带中南段,由于受青藏高原东向的挤压作用,历史上曾发生多次破坏性地震。该地区地质构造复杂多样,地震活动性强,宜于我国大陆内部地震成因及其活动规律的探索。威远位于四川盆地的中西部,盆地结构相对稳定,地表构造形态主要受到周边大型构造带活动的影响(Zheng et al,2013)。威远作为我国长宁−威远国家级页岩气示范区,历史上地震活动较弱(董大忠等,2014;谢军,2018),但近年来,特别是2015年以来(图1),威远地区地震活动性逐渐增强(易桂喜等,2019,2020;Yang et al,2020)。较强的地震活动性对页岩气开采和平台设施的影响较大,特别是反复发生的中强地震对其可能会造成严重的破坏。地壳应力场反映了浅部岩石的受力情况,是地震学研究的重要内容,也是探索地震、滑坡等地质灾害的重要手段。地壳应力场的空间分布和变化特征可以为页岩气开采提供参考,有效规避由于水力压裂施工注水而诱发强震的潜在风险。研究区内近年来频繁发生的中小地震为初步探讨地震各向异性和地壳应力场特征提供了丰富的数据资源。
剪切波分裂是获得地球内部介质各向异性和地壳应力场特征的重要手段之一(Crampin,1981;Vinnik et al,1989;Gao et al,2011)。理论分析与高压高温实验表明,低压条件下的地壳浅层的各向异性主要受裂隙定向排列的影响(张国苓等,2010)。当剪切波在各向异性介质中传播时,会分裂成快、慢剪切波,这两列波传播速度不同且偏振方向近似垂直。进行剪切波分裂计算可获得快波偏振方向和慢波延迟时间两个重要参数,其中:快波偏振方向反映了原地主压应力方向;而慢波延迟时间反映了介质的各向异性程度,对地壳应力场的变化非常敏感(Gao,Crampin,2004,2008)。因此,借由剪切波分裂方法可以获得地壳介质各向异性分布,探讨地壳应力场的特征和变化(吴晶等,2007),从而为页岩气开采提供参考,有效规避地震灾害风险。目前,针对四川地区的各向异性研究,无论在活动断裂带附近还是水库库区,已取得了一定的成果,例如石玉涛等(2009,2013)研究了龙门山断裂带附近剪切波分裂参数的区域特征以及与活动断裂之间的关系;张永久等(2008)利用四川区域地震台网数据获得了四川地区的地壳各向异性。根据剪切波分裂参数的变化还对四川紫坪铺水库库区地壳应力场的特征性变化进行了探讨(张永久等,2010;刘莎,吴朋,2015)。剪切波分裂方法用于地壳应力场的探索已经发展得较为成熟,现拟将该方法应用于工业开采区,从而获得页岩气开采区的地壳各向异性空间分布特征。在威远页岩气开采区小尺度范围内密集台阵的地震观测积累了丰富的数据资料(图2)。本文拟利用威远地区40个密集流动地震台站在2019年11月—2020年5月间记录的地震波形数据,进行剪切波分裂计算并分析剪切波分裂参数的特征,探讨威远地区地壳各向异性的空间分布。
图 2 威远地区构造背景、台站和2019年11月—2020年5月地震分布F1:牛颈垇断层;F2:大垭口断层;F3:东兴场断层;F4:高桥断层;F5:资中断层;F6:墨林场断层;F7:重滩断层;F8:黄桷坡断层;F9:长岩断层,下同Figure 2. Tectonic settings,the distribution of seismic stations and the epicenters between November 2019 and May 2020 in Weiyuan regionF1:Niujing’ao fault;F2:Dayakou fault;F3:Dongxingchang fault;F4:Gaoqiao fault;F5:Zizhong fault;F6:Molinchang fault;F7:Chongtan fault; F8:Huangjuepo fault;F9:Changyan fault,the same below1. 地震数据与剪切波分裂
1.1 构造背景与地震数据
四川盆地地区的断裂构造方向以北北东或北东方向为主,在北西向的区域主压应力作用下,活动断裂一直处于正压状态,差异活动不明显(唐荣昌等,1995)。威远地区位于华蓥山断裂西侧,四川盆地低缓构造带的中南部,该区域地势西北高、东南低(何超红等,2015)。横贯威远县西北部,在地质构造上隶属于乐山—龙女寺古隆起的威远背斜为一大型穹隆构造,南陡北缓、西窄紧、东开阔,长轴呈北东向,大体与龙门山一致(朱瑜等,2010)。威远背斜主要形成于新生代,在经历了加里东运动、海西运动、印支运动、燕山运动和喜马拉雅运动的叠加作用后,形成了现今的大型穹窿构造(刘顺等,2000;喻颐等,2013,曾求等,2020)。威远地区的断裂性质主要以逆冲为主,如大垭口断层、墨林场断层、资中断层、黄桷坡断层、长岩断层、重滩断层等,断层的走向主要呈北西向、北东向和近似东西向(图2)。
以往的研究结果显示威远地区快剪切波偏振方向为北西向,并推测该区内的地壳应力场在长时间内保持稳定状态(陈安国等,2019);基于地震背景噪声得到的威远地区浅层速度结构成像结果表明剪切波速度具有与威远背斜相对应的横向不均匀性分布特征(曾求等,2020)。地震震源机制解是确定构造变形特征和震源区应力场特征的基础性资料。利用震源机制解结果推断威远地区构造应力环境相对简单,区域最大主压应力方向为北西向(易桂喜等,2020);威远地区地震重定位结果显示地震多集中发生在P波和S波的高低速异常转换带(杜广宝等,2021);基于四川地区GPS水平运动速度场数据对四川地壳应变场特征的分析结果表明,威远一带地壳运动速度场方向为南东方向(白卓立等,2020)。
四川盆地页岩气资源丰富,特别是威远一带有良好的页岩气开发前景(石强等,2017)。近五年来,威远一带频繁发生中小地震,地震活动性逐步增强,特别是2019年9月8日连续发生了MS5.4和MS4.3中强地震。该地区的地震定位和震源机制解等研究指出威远地区中强地震可能与页岩气水力压裂施工有关,是由孔隙压力变化诱发所致(Yang et al,2020)。易桂喜等(2020)的震源机制解研究结果认为威远中强地震发震构造可能是威远背斜南翼的一系列缓倾角盲冲断层,威远地区的构造活动主要受区域应力场控制。为了进一步探讨威远地区的结构特征,2019年中国地震局地球物理研究所在四川威远一带布设了40个流动地震台站进行地震活动性监测,这些台站均匀地分布在威远震区,平均台间距为10 km (图2)。各地震台站均使用三分向地震计,采样率为100 Hz。由于上地壳的近震剪切波主要频段为3—10 Hz,因此在进行剪切波分裂数据处理过程中,使用2—20 Hz的带通滤波。
1.2 数据的选择与处理
本文利用威远地区流动密集台阵在2019年11月至2020年6月记录的近震波形资料进行剪切波分裂研究。剪切波分裂计算需要对地震波形数据进行严格的筛选。首先,为了保证地震台站记录到的地震波形与入射波相类似,需要选取位于剪切波窗口内的地震波形数据。理论上,对于泊松比为0.25的介质,其剪切波窗口约为35°。实际上,由于地表低速层的影响,剪切波窗口可有效地扩大至45°—50° (Crampin,Peacock,2005)。根据荣县—威远—资中地区一维速度模型(易桂喜等,2020),本文选取的剪切波窗口为50°。在此基础上,挑选出剪切波波形清晰且信噪比较高的地震数据。本文最终选用了研究区内23个流动地震台站的波形数据。之后采用偏振分析法获得剪切波分裂参数,截取剪切波到时前20个采样点和后40个采样点的数据进行偏振分析(图3)。根据剪切波质点偏振图的变化来确定快波偏振方向和慢波延迟时间。由于地壳各向异性介质的存在,剪切波质点偏振图的方向会突然呈近垂直向的偏转,据此来获得快剪切波和慢剪切波的到时(刘莎等,2015)。
图 3 地震波形及偏振分析示例(a) 台站SC38记录的2020-12-06 17:05:08发生的ML2.1地震波形;(b) 台站SC12记录的2021-05-11 21:40:46发生的ML1.5地震波形;(c) 台站SC06记录的2021-04-07 11:26:42发生的ML1.9地震波形。左图为三分量地震波形,t1,t2,t3为用于偏振分析的剪切波时间段;右图为三个时段剪切波质点偏振图,S1和S2分别为快剪切波和慢剪切波的到时Figure 3. Seismic waveforms and diagrams of the polarization analysis(a) Seismic waves of the ML2.1 earthquake at 17:05:08 on December 6,2020 recorded by the station SC38;(b) Seismic waves of the ML1.5 earthquake at 21:40:46 on May 11,2021 recorded by the station SC12;;(c) Seismic waves of the ML1.9 earthquake at 11:26:42 on April 7,2021 recorded by the station SC06 ;and the left panels are the three-component records of seismic waveform. t1,t2 and t3 are the time periods of shear wave for polarization analysis,respectively. The right panels are the particle trails of the shear waves in three time periods. S1 and S2 are the arrival times of fast and slow shear waves,respectively.2. 计算结果
通过对威远地区流动地震台站剪切波窗口数据的筛选并进行剪切波分裂计算,得到了23个地震台站的快波偏振方向和慢波延迟时间,剪切波分裂参数共813对(表1)。多数台站获得了10对以上的剪切波分裂参数,特别是SC08,SC15,SC16,SC27和SC38台站均获得了50对以上的剪切波分裂参数结果,保证了数据结果的可靠性。SC03,SC20,SC21和SC30台站虽然有较少的剪切波分裂参数结果,但是这些台站剪切波分裂参数结果与其它台站具有较好的一致性。
表 1 威远地区各台站剪切波分裂参数结果Table 1. The results of shear wave splitting parameters at the seismic stations in Weiyuan area台站
代码有效事件
个数优势方向
个数快波偏振
方向/°慢波延迟时间
/(ms·km−1)台站
代码有效事件
个数优势方向
个数快波偏振
方向/°慢波延迟时间
/(ms·km−1)SC01 34 1 −39±17.8 4.85±1.21 SC15 78 1 −42±18.3 5.23±1.72 SC02 32 1 −34±30.1 5.12±1.04 SC16 43 1 −46±15.7 4.66±1.38 SC03 5 1 −55±4.0 3.74±0.85 46 2 51±13.2 4.73±1.19 SC05 17 1 −46±20.8 3.98±1.01 SC17 35 1 −41±16.8 4.49±1.54 SC06 38 1 −49±10.9 4.13±1.28 SC19 15 1 −17±23.7 3.31±0.36 11 2 37±18.3 3.99±1.06 SC20 7 1 −39±11.7 3.41±0.71 SC07 36 1 20±14.6 4.35±1.13 SC21 6 1 −36±12.9 3.48±0.76 SC08 41 1 −46±15.3 4.76±1.42 SC22 31 1 −33±25.6 3.65±1.03 12 2 15±9.9 4.21±1.16 SC27 67 1 −48±14.1 4.87±1.20 SC09 16 1 −39±11.6 4.40±1.02 SC28 48 1 −44±18.2 4.34±0.68 7 2 25±17.5 3.89±0.88 SC30 5 1 −57±13.5 4.92±1.38 SC10 13 1 −43±16.1 3.37±1.18 SC38 53 1 −45±14.7 4.52±1.15 SC11 15 1 −31±20.5 4.24±1.61 14 2 56±11.8 4.33±1.07 SC12 24 1 −42±15.8 5.06±1.80 SC39 47 1 −36±22.3 4.37±1.06 17 2 44±13.2 4.62±1.21 研究区域内23个台站快波偏振方向等面积投影玫瑰图如图4所示,可以看出多数台站快波偏振优势方向较为明显,有部分台站的剪切波分裂参数结果显示出了两个快波偏振优势方向。而慢波延迟时间值代表了震源到台站整个射线路径上慢剪切波相对于快剪切波的时间差,反映了整个射线路径上的各向异性强度。考虑到不同地震事件的震源距差别较大,为了便于分析比较,采用归一化慢波延迟时间来反映地壳介质的各向异性程度,即单位射线路径长度的慢波延迟时间,归一化的慢波延迟时间单位为ms/km。本文获得的威远地区慢波延迟时间均值为(4.43±1.15) ms/km。SC19台站慢波延迟时间最小,为3.31 ms/km;SC15台站最大,为5.23 ms/km。
各台站快波偏振方向等面积投影玫瑰图显示,威远地区台站的快波偏振方向一致性较好,其中:SC01,SC02,SC03,SC05,SC10,SC11,SC15,SC17,SC20,SC21,SC22,SC27,SC28,SC30和SC39等15个台站具有单一的快波偏振优势方向且显示为北西向;SC19台站的快波偏振方向结果较为零散,误差较大,其单一的快波偏振优势方向为北北西向;SC07台站快波偏振优势方向只有一个,为北北东向;SC06,SC12,SC16,SC38台站快波偏振优势方向有两个,为北西向和北东向;SC08和SC09台站也有两个快波偏振优势方向,为北西向和北北东向。
3. 地壳各向异性特征
3.1 快波偏振方向的空间分布
快波偏振方向反映了地壳中微裂隙的走向,而区域应力场控制着地壳微裂隙的排列方向,使其平行于最大主压应力的方向,最终使得地壳各向异性介质的快波偏振方向平行于最大主压应力方向(王雪鹤等,2021)。快波偏振方向是剪切波在地下传播路径上对各向异性介质的综合表现,定性地反映了地壳应力场的方向(陈安国等,2019)。但是在许多构造复杂的区域,地壳剪切波分裂观测研究结果显示:部分台站的快波偏振方向出现了两个或多个优势方向,这通常是因为除了受到区域应力的控制之外,还受到附近活动断裂、不规则地表地形等因素的影响(陈安国等,2019;Gao et al,2011)。
为了更好地分析威远地区的地壳各向异性和区域应力分布情况,将每个台站的快波偏振优势方向结果展示在空间分布图上,得到了威远地区快波偏振方向空间分布特征,如图5所示。
位于牛颈垇断层附近的台站SC19快波偏振方向呈北北西向,与断裂走向近似平行。位于该断层西侧的台站SC21和东侧的台站SC03快波偏振方向均为北西向,与断层走向斜交,但是与区域主压应力方向平行。台站SC20快波偏振方向为北西向,其东侧的固定台站HMS在以往的剪切波分裂研究中的快波偏振方向显示为北西向(陈安国等,2019),这两个台站快波偏振方向一致。但是台站HMS附近的大垭口断层走向为北东东向,与附近台站的快波偏振方向斜交。东兴场断层走向近似为东西向,其东侧的台站SC02快波偏振方向呈北西向,表明该台站快波偏振方向受断层的影响不大。台站SC07位于高桥断层南段,快波偏振方向与断层的走向一致,均为北北东向;而位于SC07附近的台站SC01快波偏振方向为北西向,与台站SC07结果不一致,这表明台站SC07快波偏振方向受断层的影响较大,而其西侧的台站SC01则主要受区域应力场的控制。台站SC17快波偏振方向呈北西向,与区域主压应力的方向一致。位于威远背斜附近的台站SC16,SC06,SC38,SC08和SC09均表现出两个快波偏振优势方向,北西向的快波偏振方向与区域主压应力的方向一致,北东向的快波偏振方向与主压应力方向近似垂直,这可能与局部断裂构造有关。位于墨林场断层以西的台站SC22和SC39,其快波偏振优势方向都显示为北西向,与断层走向和主压应力方向均一致。重滩断层西南侧的台站SC30快波偏振方向与断层走向一致,为北西向。北西向的黄桷坡断层和北东向的长岩断层斜交,位于黄桷坡断层北侧的台站SC27和SC28快波偏振优势方向显示为北西方向,快波偏振方向受断层影响较小,与区域主压应力方向平行。台站SC15快波偏振方向表现为北西向,与西侧的墨林场断层方向一致。台站SC12快波偏振优势方向显示为北西向和北东向。其东侧的台站SC11和SC10快波偏振方向为单一的北西向,与区域主压应力方向平行。由此可见,威远地区台站快波偏振优势方向多表现为北西向,也有个别结果为北东向。该地区主压应力方向为北西向,区域应力场对快波偏振方向的作用较为明显。部分台站表现出的北东向的快波偏振优势方向多是受到了局部断裂构造的影响。当地震台站位于断裂上或距离断裂足够近,地震波射线路径穿过断层,那么快波偏振方向多平行于断裂;随着地震台站远离断裂位置,快波偏振方向则可能与断裂走向不同(郭铁龙,高原,2020)。在构造环境复杂的地区,快波偏振方向往往较为复杂,通常是应力场与局部构造双重作用的结果。威远地区地壳各向异性的空间分布除了受到区域应力场的影响,还受到复杂的局部构造作用的影响。
威远位于四川盆地中西部,而四川盆地属于扬子地块,位于青藏高原东南缘。青藏高原东南缘地处青藏高原侧向挤出变形的前缘,在高原的边缘过渡带地壳应力场以挤压为主,主压应力方向为北西—南东(李长军等,2019;Han et al,2019;黄周传等,2021)。威远地区的快波偏振方向多呈北西方向,与青藏高原东南缘向东推挤的构造特征相一致(黄周传等,2021)。对比四川地区中下地壳各向异性研究结果中得到的Pms波分裂的快波偏振方向为近似北西向(杨妍,2019),本文获得的上地壳各向异性研究结果与其近似。这说明在研究区域内上下地壳的各向异性快波偏振方向相近,威远地区的地壳内部变形特征一致。上地幔远震SKS波分裂结果显示威远地区上地幔各向异性方向为北西向(常利军等,2008),这与本文得到的上地壳各向异性方向相似。综上分析,威远地区岩石圈变形可能是遵循壳幔耦合的垂直连贯变形模式。
从快波偏振优势方向空间分布可以看出围绕威远背斜构造单元的地震台站快波偏振方向较为复杂,多数台站显示有两个快波偏振优势方向。威远背斜是研究区域的主要构造单元,主要的构造运动是较为强烈的挤压褶皱,使古老沉积地层弯曲(徐炳高,1997),从而致使威远背斜区域的沉积盖层起伏变化较大(刘顺,2001)。背景噪声三维成像结果显示以威远背斜为分界线,其两侧的横波波速变化趋势具有很明显的差异,西北部高于东南(曾求等,2020)。本文中获得具有两个快波偏振优势方向的台站集中于威远背斜的两侧,可见在波速变化分界的区域快波偏振方向较为零散,北东向的快波偏振方向受到局部构造环境的影响,表现出与区域主压应力不一致的各向异性方向。
3.2 慢波延迟时间分布
威远地区慢波延迟时间的空间分布有一定的区域特征(图6)。本文中所获得的有两个快波偏振优势方向的台站有6个,其中SC06,SC16和SC38台站在两个快波偏振优势方向上的慢波延迟时间值基本一致,相差很小;而SC08,SC09和SC12台站在两个快波偏振优势方向上的慢波延迟时间有一定的差距,其中SC08和SC09台站位于威远背斜褶皱拱起的位置。慢波延迟时间对台站下方地壳中微裂隙的几何形态变化较为敏感,反映了该处的各向异性强度(刘莎等,2015)。在两个不同的快波偏振优势方向上慢波延迟时间的差异表明了各向异性强弱不同,推测该台站下方地壳各向异性的复杂性。整个研究区域内,慢波延迟时间最大值出现在台站SC15,该台站位于威远背斜南段南侧。在空间分布上,威远背斜周围及其以南区域的各台站慢波延迟时间普遍大于研究区域北部地区的各台站,如SC19,SC20以及SC21台站的慢波延迟时间值相对较小。慢波延迟时间的空间分布表明研究区域内地壳各向异性强度有所差异,威远背斜以南区域的地壳各向异性强于北部区域。
4. 讨论与结论
四川威远地区是我国页岩气开采的主要区域之一。近年来,特别是2015年以来,威远地区地震活动性显著增强。较高的地震活动性以及不断发生的中小地震对页岩气开采和平台设施存在较大的安全隐患。本文对威远地区流动密集台阵在2019年11月至2020年5月期间记录的近震波形数据经过剪切波窗口数据的筛选和剪切波分裂计算,共获得了威远地区23个台站的剪切波分裂参数。根据快波偏振方向和慢波延迟时间的空间分布,探讨了该地区地壳各向异性的区域性特征。
研究区域内多数台站快波偏振优势方向明显,其中6个台站有两个快波偏振优势方向。整体上看,威远地区地震台站快波偏振方向表现得较为一致,多数台站快波偏振优势方向呈北西向,与区域主压应力方向一致,只有个别台站快波偏振优势方向显示为北东向。快波偏振方向的空间分布结果表明,有两个快波偏振优势方向的台站集中分布在威远背斜周围。在构造复杂的区域,个别台站的快波偏振方向出现两个或多个优势方向,通常是因为快波偏振方向除了受到区域应力场的控制之外,还受到附近活动断裂、不规则地表地形等因素的影响。威远背斜是研究区域的主要构造单元,其主要的构造运动是较为强烈的挤压褶皱(曾求等,2020)。而且在威远背斜两侧剪切波速度有明显差异,所以在威远背斜这一复杂的构造单元周围地壳介质的各向异性程度较为复杂。
威远地区慢波延迟时间主要分布在3.31—5.23 ms/km范围内,均值为4.43 ms/km。慢波延迟时间最小值出现在SC19台站,位于研究区域北部的牛颈垇断层南端。慢波延迟时间最大值出现在SC15台站,该台站位于威远背斜南段南侧。从慢波延迟时间的空间分布可以看出,威远背斜南侧区域的慢波延迟时间普遍大于北侧区域,反映出研究区域南侧地壳各向异性强度要高于北侧地区。在威远背斜褶皱拱起的位置上,具有两个快波偏振优势方向的台站在两个偏振优势方向上的慢波延迟时间值差异较大。该结果进一步说明了威远背斜周围地壳介质的各向异性特征较为复杂。
在威远页岩气开采区,地震分布较为密集。已有研究表明,页岩气水力压裂的工业流体注入能显著改变当地的构造应力环境。在一定条件下,流体的注入可能使过去不活动的老断层重新活化或者不成熟的断层发生破裂,从而诱发大量地震(Lei et al,2019a,b;易桂喜等,2020)。地震重定位的结果显示在威远地区多数地震事件的空间位置与地表已知断层无明显对应关系,但是中小地震活动与地壳速度结构变化具有一定的相关性(杜广宝等,2021)。本文中获得的快波偏振方向在地壳速度结构显著变化的威远背斜区域表现出复杂的特征,而且慢波延迟时间基本以威远背斜为界,呈现出南高北低的各向异性强度。剪切波分裂参数复杂的区域对应于中小地震发生的聚集区。而该区域也是水力压裂施工活动区域,因此,威远地区快波偏振方向的复杂性与水力压裂施工可能有一定相关性。但是,由于本文数据有限,此方面的研究还需要后期结合其它资料进行进一步论证。
中国地震局地球物理研究所地震科学数据中心为本研究提供波形数据服务,三位审稿专家为本文提出了非常有益的意见,作者在此一并表示感谢。
-
图 1 盐源盆地所处位置及其地质构造 [ 修改自卢海建等(2015) ]
Figure 1. Location of Yanyuan basin and its geological structure modified from Lu et al (2015)
图 3 (a) 部分台站对的互相关波形(0.2—2 Hz带通滤波);(b) 不同周期T的瑞雷波群速度对S波速度的敏感核;(c) 瑞雷波群速度频散曲线测量图,红色圆点表示所提取的频散点
Figure 3. (a) The cross-correlation seismograms between some stations (0.2—2 Hz);(b) Depth sensitivity kernels of Rayleigh wave group velocity to the S-wave velocity at different periods T;(c) Rayleigh wave group velocity measurements from empirical Green’s functions with red dots as the extracted dispersion points
图 6 盐源盆地地震剖面处理与解释图
(a) 地震能量剖面图;(b) 地震反演所得P波速度结构图;(c) 地震解释剖面,图中震源机制解引自GCMT (2018)
Figure 6. Seismic profile processing and interpretation for Yanyuan basin
(a) Seismic energy profile;(b) The seismic P-wave velocity structure profile by inversion;(c) Interpretation profile where the focal mechanism solution after GCMT (2018)
-
程万正,杨永林. 2002. 川滇地块边界构造带形变速率变化与成组强震[J]. 大地测量与地球动力学,22(4):21–25. Cheng W Z,Yang Y L. 2002. Deformation rate changes of tectonic belts along boundaries of Yunnan-Sichuan block and their relation to grouped strong earthquakes[J]. Journal of Geodesy and Geodynamics,22(4):21–25 (in Chinese).
董蕾,沈旭章,钱银苹. 2020. 青藏高原东南缘Moho面速度密度跃变研究[J]. 地球物理学报,63(3):915–927. doi: 10.6038/cjg2020N0168 Dong L,Shen X Z,Qian Y P. 2020. Study on velocity and density contrasts across the Moho in the southeastern margin of the Tibetan Plateau[J]. Chinese Journal of Geophysics,63(3):915–927 (in Chinese).
葛肖虹. 1984. 川西盐源推覆构造的探讨[J]. 长春地质学院学报,(1):36–43. Ge X H. 1984. A discussion on nappe structure in Yanyuan,west Sichuan[J]. Journal of Changchun College of Geology,(1):36–43 (in Chinese).
金胜,魏文博,汪硕,叶高峰,邓明,谭捍东. 2010. 青藏高原地壳高导层的成因及动力学意义探讨:大地电磁探测提供的证据[J]. 地球物理学报,53(10):2376–2385. Jin S,Wei W B,Wang S,Ye G F,Deng M,Tan H D. 2010. Discussion of the formation and dynamic signification of the high conductive layer in Tibetan crust[J]. Chinese Journal of Geophysics,53(10):2376–2385 (in Chinese).
李立,金国元. 1987. 攀西裂谷带及龙门山断裂带地壳上地幔的大地电磁测深研究[J]. 物探与化探,11(3):161–169. Li L,Jin G Y. 1987. Telluric electromagnetic sounding study of crust and upper mantle in the Panxi “rift zone” and the Longmenshan faulted zone[J]. Geophysical &Geochemical Exploration,11(3):161–169 (in Chinese).
李生. 2004. 四川锦屏山地区推覆构造带特征及其研究意义[J]. 沉积与特提斯地质,24(1):70–77. doi: 10.3969/j.issn.1009-3850.2004.01.011 Li S. 2004. Nappe structures in the Jinpingshan area,Sichuan and their geological implications[J]. Sedimentary Geology and Tethyan Geology,24(1):70–77 (in Chinese).
李勇,侯中健,司光影,周荣军,梁兴中. 2001. 青藏高原东南缘晚第三纪盐源构造逸出盆地的沉积特征与构造控制[J]. 矿物岩石,21(3):34–43. doi: 10.3969/j.issn.1001-6872.2001.03.006 Li Y,Hou Z J,Si G Y,Zhou R J,Liang X Z. 2001. Sedimentary characteristics and tectonic controls in Neogene Yanyuan tectonic escape basin in southeastern Qinghai-Tibet Plateau[J]. Journal of Mineralogy and Petrology,21(3):34–43 (in Chinese).
廖忠礼,邓永福,廖光宇. 2003. 四川锦屏地区新生代冲断作用[J]. 大地构造与成矿学,27(2):152–159. doi: 10.3969/j.issn.1001-1552.2003.02.008 Liao Z L,Deng Y F,Liao G Y. 2003. Formation and evolution of Mesozoic thrust fold belt in Jinping area,Sichuan[J]. Geotectonica et Metallogenia,27(2):152–159 (in Chinese).
刘家铎,刘文周. 1995. 盐源盆地演化及早三叠世青天堡组沉积环境[J]. 成都理工学院学报,22(3):70–74. Liu J D,Liu W Z. 1995. The evolution of Yanyuan basin and the sedimentary environment of the Early Triassic Qingtianbao group[J]. Journal of Chengdu Institute of Technology,22(3):70–74 (in Chinese).
刘薇,邵昌盛,李大虎,顾勤平. 2018. 利用重力、航磁资料研究木里—盐源弧形构造带及邻区的深部孕震环境[J]. 成都理工大学学报(自然科学版),45(2):250–256. doi: 10.3969/j.issn.1671-9727.2018.02.12 Liu W,Shao C S,Li D H,Gu Q P. 2018. Application of gravity and aeromagnetic data to the study of deep seismogenic environment of the Muli-Yanyuan arc structure belt and its adjacent area in Sichuan,China[J]. Journal of Chengdu University of Technology (Science &Technology Edition)
,45(2):250–256 (in Chinese). 卢海建,王二七,李仕虎,李海兵. 2015. 青藏高原东南缘构造旋转变形分析:以四川盐源盆地古地磁研究为例[J]. 中国地质,42(5):1188–1201. doi: 10.3969/j.issn.1000-3657.2015.05.002 Lu H J,Wang E Q,Li S H,Li H B. 2015. Rotational deformation of the southeastern margin of Tibet:A paleomagnetic study of the Yanyuan basin,Sichuan Province[J]. Geology in China,42(5):1188–1201 (in Chinese).
骆耀南,俞如龙. 2002. 西南三江地区造山演化过程及成矿时空分布[J]. 地球学报,23(5):417–422. doi: 10.3321/j.issn:1006-3021.2002.05.005 Luo Y N,Yu R L. 2002. Orogenic evolution and Metallogenic time-space distribution in Jinshajiang-Lancangjiang-Nujiang region,Southwest China[J]. Acta Geoscientia Sinica,23(5):417–422 (in Chinese).
牟雅元,代宽宏,张振勋,邱发青,吴今生,胡德军,王松,代放,刘德安,胡方良. 2004. 2003年8月21日盐源5.0级地震的预测与分析[J]. 国际地震动态,(2):9–13. doi: 10.3969/j.issn.0253-4975.2004.02.003 Mou Y Y,Dai K H,Zhang Z X,Qiu F Q,Wu J S,Hu D J,Wang S,Dai F,Liu D A,Hu F L. 2004. Prediction of the Yanyuan MS5.0 earthquake on Aug. 21,2003 and the relevant analysis[J]. Recent Developments in World Seismology,(3):9–13 (in Chinese).
四川省地质矿产研究所专题研究组. 1987. 盐源—丽江地区三叠纪地层及沉积相[M]. 北京: 地质出版社: 107−116. Special Research Group of Sichuan Institute of Geology and Mineral Resources. 1987. The Triassic Strata and Sedimentary Facies in Yanyuan-Lijiang District[M]. Beijing: Geological Publishing House: 107−116 (in Chinese).
孙洁,晋光文,白登海,王立凤. 2003. 青藏高原东缘地壳、上地幔电性结构探测及其构造意义[J]. 中国科学:D辑,33(增刊1):173–180. Sun J,Jin G W,Bai D H,Wang L F. 2003. Sounding of electrical structure of the crust and upper mantle along the eastern border of Qinghai-Tibet Plateau and its tectonic significance[J]. Science in China:Series D,46(2):243–253.
谭捍东,姜枚,吴良士,魏文博. 2006. 青藏高原电性结构及其对岩石圈研究的意义[J]. 中国地质,33(4):906–911. doi: 10.3969/j.issn.1000-3657.2006.04.020 Tan H D,Jiang M,Wu L S,Wei W B. 2006. Electrical conductivity structure of the Qinghai-Tibet Plateau and its significance for the study of the lithosphere[J]. Geology in China,33(4):906–911 (in Chinese).
唐若龙. 1987. 木里—盐源推覆构造特征与金铜、铅锌的成矿关系[J]. 四川地质学报,7(2):3–9. Tang R L. 1987. Relationship between the characteristics of Muli-Yanyuan nappe structure and the mineralization of gold,copper,lead and zinc[J]. Acta Geologica Sichuan,7(2):3–9 (in Chinese).
滕吉文,尹周勋,刘宏兵,张中杰,胡家富,孙克忠,魏计春. 1994. 青藏高原岩石层三维和二维结构与大陆动力学[J]. 地球物理学报,37(增刊2):117–130. Teng J W,Yin Z X,Liu H B,Zhang Z J,Hu J F,Sun K Z,Wei J C. 1994. The 2D and 3D lithosphere structure and continental dynamics of Qinghai-Xizang Plateau[J]. Acta Geophysica Sinica,37(S2):117–130 (in Chinese).
王椿镛,吴建平,楼海,常利军,苏伟. 2006. 青藏高原东部壳幔速度结构和地幔变形场的研究[J]. 地学前缘,13(5):349–359. doi: 10.3321/j.issn:1005-2321.2006.05.010 Wang C Y,Wu J P,Lou H,Chang L J,Su W. 2006. Study of crustal and upper mantle’s structure and mantle deformation field beneath the eastern Tibetan Plateau[J]. Earth Science Frontiers,13(5):349–359 (in Chinese).
王夫运,段永红,杨卓欣,张成科,赵金仁,张建狮,张先康,刘启元,朱艾斓,徐锡伟,刘宝峰. 2008. 川西盐源—马边地震带上地壳速度结构和活动断裂研究:高分辨率地震折射实验结果[J]. 中国科学:D辑,38(5):611–621. Wang F Y,Duan Y H,Yang Z Y,Zhang C K,Zhao J R,Zhang J S,Zhang X K,Liu Q Y,Zhu A L,Xu X W,Liu B F. 2008. Velocity structure and active fault of Yanyuan-Mabian seismic zone:The result of high-resolution seismic refraction experiment[J]. Science in China:Series D,51(9):1284–1296. doi: 10.1007/s11430-008-0098-0
王琼,高原. 2014. 青藏东南缘背景噪声的瑞利波相速度层析成像及强震活动[J]. 中国科学:地球科学,44(11):2440–2450. Wang Q,Gao Y. 2014. Rayleigh wave phase velocity tomography and strong earthquake activity on the southeastern front of the Tibetan Plateau[J]. Science China Earth Sciences,57(10):2532–2542. doi: 10.1007/s11430-014-4908-2
汪一鹏,沈军,王琪,熊熊. 2003. 川滇块体的侧向挤出问题[J]. 地学前缘,10(增刊1):188–192. Wang Y P,Shen J,Wang Q,Xiong X. 2003. On the lateral extrusion of Sichuan-Yunnan block (Chuan-Dian block)[J]. Earth Science Frontiers,10(S1):188–192 (in Chinese).
王正和,邓敏,程锦翔,张海全. 2018. 康滇古陆西侧断裂及岩浆活动对油气保存条件的影响:以盐源盆地为例[J]. 地球科学,43(10):3616–3624. Wang Z H,Deng M,Cheng J X,Zhang H Q. 2018. Influence of fault and magmatism on oil and gas preservation condition,to the west of Kangdian ancient continent:Taking Yanyuan basin as an example[J]. Earth Science,43(10):3616–3624 (in Chinese).
韦伟,孙若昧,石耀霖. 2010. 青藏高原东南缘地震层析成像及汶川地震成因探讨[J]. 中国科学:地球科学,40(7):831–839. Wei W,Sun R M,Shi Y L. 2010. P-wave tomographic images beneath southeastern Tibet:Investigating the mechanism of the 2008 Wenchuan earthquake[J]. Science China Earth Sciences,53(9):1252–1259. doi: 10.1007/s11430-010-4037-5
魏文博,金胜,叶高峰,邓明,景建恩,Martyn U,Alan J G. 2009. 藏南岩石圈导电性结构与流变性:超宽频带大地电磁测深研究结果[J]. 中国科学:D辑,39(11):1591–1606. Wei W B,Jin S,Ye G F,Deng M,Jing J E,Martyn U,Alan J G. 2010. Conductivity structure and rheological property of lithosphere in southern Tibet inferred from super-broadband magnetotelluric sounding[J]. Science in China:Series D,53(2):189–202. doi: 10.1007/s11430-010-0001-7
韦一,张宗言,何卫红,吴年文,杨兵. 2014. 上扬子地区中生代沉积盆地演化[J]. 地球科学:中国地质大学学报,39(8):1065–1078. Wei Y,Zhang Z Y,He W H,Wu N W,Yang B. 2014. Evolution of sedimentary basins in the Upper Yangtze during Mesozoic[J]. Earth Science:Journal of China University of Geosciences,39(8):1065–1078 (in Chinese). doi: 10.3799/dqkx.2014.095
向宏发,徐锡伟,虢顺民,张晚霞,李洪武,于贵华. 2002. 丽江—小金河断裂第四纪以来的左旋逆推运动及其构造地质意义:陆内活动地块横向构造的屏蔽作用[J]. 地震地质,24(2):188–198. doi: 10.3969/j.issn.0253-4967.2002.02.006 Xiang H F,Xu X W,Guo S M,Zhang W X,Li H W,Yu G H. 2002. Sinistral thrusting along the Lijiang-Xiaojinhe fault since Quaternary and its geologic-tectonic significance:Shielding effect of transverse structure of intracontinental active block[J]. Seismology and Geology,24(2):188–198 (in Chinese).
徐锡伟,闻学泽,郑荣章,马文涛,宋方敏,于贵华. 2003. 川滇地区活动块体最新构造变动样式及其动力来源[J]. 中国科学:D辑,33(增刊1):151–162. Xu X W,We X Z,Zheng R Z,Ma W T,Song F M,Yu G H. 2003. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region,China[J]. Science in China:Series D,46(2):210–226.
杨文采,侯遵泽,于常青. 2015. 青藏高原地壳的三维密度结构和物质运动[J]. 地球物理学报,58(11):4223–4234. Yang W C,Hou Z Z,Yu C Q. 2015. Three-dimensional density structure of the Tibetan Plateau and crustal mass movement[J]. Chinese Journal of Geophysics,58(11):4223–4234 (in Chinese).
杨文采,侯遵泽,徐义贤,颜萍. 2017. 青藏高原下地壳热变形和管道流研究[J]. 地质论评,63(5):1141–1152. Yang W C,Hou Z Z,Xu Y X,Yan P. 2017. A study on thermal deformation and lower crust channel flows in Qinghai-Xizang (Tibet) Plateau[J]. Geological Review,63(5):1141–1152 (in Chinese).
杨文采,瞿辰,任浩然,黄联捷,胥颐,于常青,刘晓宇. 2019a. 青藏高原软流圈与特提斯洋板块俯冲[J]. 地质论评,65(3):521–532. Yang W C,Qu C,Ren H R,Huang L J,Xu Y,Yu C Q,Liu X Y. 2019a. The asthenosphere of the Qinghai-Xizang (Tibetan) Plateau and subduction of the Tethys Ocean[J]. Geological Review,65(3):521–532 (in Chinese).
杨文采,瞿辰,任浩然,黄连捷,胥颐,于常青. 2019b. 青藏高原地壳地震纵波速度的层析成像[J]. 地质论评,65(1):2–14. Yang W C,Qu C,Ren H R,Huang L J,Xu Y,Yu C Q. 2019b. Crustal P-wave seismic tomography of the Qinghai-Xizang (Tibetan) Plateau[J]. Geological Review,65(1):2–14 (in Chinese).
杨文采,金胜,张罗磊,瞿辰,胡祥云,魏文博,于常青,于鹏. 2020. 青藏高原岩石圈三维电性结构[J]. 地球物理学报,63(3):817–827. doi: 10.6038/cjg2020N0197 Yang W C,Jin S,Zhang L L,Qu C,Hu X Y,Wei W B,Yu C Q,Yu P. 2020. The three-dimensional resistivity structures of the lithosphere beneath the Qinghai-Tibet Plateau[J]. Chinese Journal of Geophysics,63(3):817–827 (in Chinese).
杨妍. 2019. 川滇地区地壳速度结构及方位各向异性的联合反演研究[D]. 合肥: 中国科学技术大学: 21−24. Yang Y. 2019. Joint Inversion for Crustal Velocity Structure and Azimuthal Anisotropy: Application to the Sichuan-Yunnan Region, SW China[D]. Hefei: University of Science and Technology of China: 21−24 (in Chinese).
杨卓欣,王夫运,段永红,张成科,赵金仁,张建狮,刘宝峰. 2011. 川滇活动地块东南边界基底结构:盐源—西昌—昭觉—马湖深地震测深剖面结果[J]. 地震学报,33(4):431–442. doi: 10.3969/j.issn.0253-3782.2011.04.003 Yang Z X,Wang F Y,Duan Y H,Zhang C K,Zhao J R,Zhang J S,Liu B F. 2011. Basement structure of southeastern boundary region of Sichuan-Yunnan active block:Analysis result of Yanyuan-Xichang-Zhaojue-Mahu deep seismic sounding profile[J]. Acta Seismologica Sinica,33(4):431–442 (in Chinese).
曾融生,朱介寿,周兵,丁志峰,何正勤,朱露培,骆循,孙为国. 1992. 青藏高原及其东部邻区的三维地震波速度结构与大陆碰撞模型[J]. 地震学报,11(增刊1):523–533. Zeng R S,Zhu J S,Zhou B,Ding Z F,He Z Q,Zhu L P,Luo X,Sun W G. 1992. Three dimensional seismic wave velocity structure and continental collision model in Qinghai-Tibet Plateau and its eastern neighbor area[J]. Acta Seismologica Sinica,11(S1):523–533 (in Chinese).
张丽敏,张志斌,崔文玲. 2014. 川西盐源两个碳酸岩杂岩体的厘定及其地质意义[J]. 大地构造与成矿学,38(1):131–139. Zhang L M,Zhang Z B,Cui W L. 2014. Discovery of two carbonatite intrusive complexes in Yanyuan area of western Sichuan and its geological significance[J]. Geotectonica et Metallogenia,38(1):131–139 (in Chinese).
张岳桥,李海龙. 2016. 青藏高原东部晚新生代重大构造事件与挤出造山构造体系[J]. 中国地质,43(6):1829–1852. Zhang Y Q,Li H L. 2016. Late Cenozoic tectonic events in east Tibetan Plateau and extrusion-related orogenic system[J]. Geology in China,43(6):1829–1852 (in Chinese).
钟康惠,刘肇昌,施央申,李凡友,舒良树. 2004. 盐源—丽江构造带是新生代陆内造山带[J]. 地质学报,78(1):36–43. doi: 10.3321/j.issn:0001-5717.2004.01.005 Zhong K H,Liu Z C,Shi Y S,Li F Y,Shu L S. 2004. Yanyuan-Lijiang tectonic zone:A Cenozoic intracontinental orogenic belt[J]. Acta Geologica Sinica,78(1):36–43 (in Chinese).
中国地震台网中心. 2020. 历史查询[EB/OL]. [2020-12-12]. http://www.ceic.ac.cn/history. China Earthquake Networks Center. 2020. History query[EB/OL]. [2020-12-12]. http://www.ceic.ac.cn/history (in Chinese).
朱民,陈汉林,周静,杨树锋. 2016. 上扬子西南盐源盆地早三叠世物源体系及构造意义[J]. 地球科学,41(8):1309–1321. Zhu M,Chen H L,Zhou J,Yang S F. 2016. Provenance of Early Triassic in Yanyuan basin,upper Yangtze and its implication for the tectonic evolution[J]. Earth Science,41(8):1309–1321 (in Chinese).
周友华. 1986. 我国川滇菱形地块运动的探讨[J]. 地震研究,9(1):89–99. Zhou Y H. 1986. Discussion on the movement of Sichuan-Yunnan rhombus block in China[J]. Journal of Seismological Research,9(1):89–99 (in Chinese).
Bai D H,Unsworth M J,Meju M A,Ma X B,Teng J W,Kong X R,Sun Y,Sun J,Wang L F,Jiang C S,Zhao C P,Xiao P F,Liu M. 2010. Crustal deformation of the eastern Tibetan Plateau revealed by magnetotelluric imaging[J]. Nat Geosci,3(5):358–362. doi: 10.1038/ngeo830
Bao X W,Sun X X,Xu M J,Eaton D W,Song X D,Wang L S,Ding Z F,Mi N,Li H,Yu D Y,Huang Z C,Wang P. 2015. Two crustal low-velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave dispersion and receiver functions[J]. Earth Planet Sci Lett,415:16–24. doi: 10.1016/j.jpgl.2015.01.020
Beaumont C,Jamieson R A,Nguyen M H,Medvedev S. 2004. Crustal channel flows:1. Numerical models with applications to the tectonics of the Himalayan-Tibetan orogen[J]. J Geophys Res:Solid Earth,109(B6):B06406.
Burg J P,Chen G M. 1984. Tectonics and structural zonation of southern Tibet,China[J]. Nature,311(5983):219–223. doi: 10.1038/311219a0
GCMT. 2018. Global CMT catalog search[EB/OL]. [2018-12-04]. https://www.globalcmt.org/CMTsearch.html.
Hodges K V,Parrish R R,Housh T B,Lux D R,Burchfiel B C,Royden L H,Chen Z. 1992. Simultaneous Miocene extension and shortening in the Himalayan orogen[J]. Science,258(5087):1466–1470. doi: 10.1126/science.258.5087.1466
Molnar P,Tapponnier P. 1977. Relation of the tectonics of eastern China to the India-Eurasia collision:Application of slip-line field theory to large-scale continental tectonics[J]. Geology,5(4):212–216. doi: 10.1130/0091-7613(1977)5<212:ROTTOE>2.0.CO;2
Molnar P,Tapponnier P. 1978. Active tectonics of Tibet[J]. J Geophys Res:Solid Earth,83(B11):5361–5375. doi: 10.1029/JB083iB11p05361
Replumaz A,Tapponnier P. 2003. Reconstruction of the deformed collision zone between India and Asia by backward motion of lithospheric blocks[J]. J Geophys Res:Solid Earth,108(B6):2285.
Royden L H,Burchfiel B C,King R W,Wang E,Chen Z L,Shen F,Liu Y P. 1997. Surface deformation and lower crustal flow in eastern Tibet[J]. Science,276(5313):788–790. doi: 10.1126/science.276.5313.788
Schoenbohm L M,Burchfiel B C,Chen L Z. 2006. Propagation of surface uplift,lower crustal flow,and Cenozoic tectonics of the southeast margin of the Tibetan Plateau[J]. Geology,34(10):813–816. doi: 10.1130/G22679.1
Shapiro N M,Ritzwoller M H,Molnar P,Levin V. 2004. Thinning and flow of Tibetan crust constrained by seismic anisotropy[J]. Science,305(5681):233–236. doi: 10.1126/science.1098276
Yang L,Liu H B,Zhao J M. 2009. Characteristic of crustal structure beneath the rifts in southern Tibetan Plateau[J]. Earthquake Science,22(4):373–377. doi: 10.1007/s11589-009-0373-y
Zhang G,Wang X B,Fang H,Guo Z M,Zhang Z B,Luo W,Cai X L,Li J,Li Z,Wu X. 2015. Crust and upper mantle electrical resistivity structure in the Panxi region of the eastern Tibetan Plateau and its significance[J]. Acta Geologica Sinica,89(2):531–541. doi: 10.1111/1755-6724.12445
Zhao W L,Yuen D A. 1987. Injection of Indian crust into Tibetan lower crust:A temperature-dependent viscous model[J]. Tectonics,6(4):505–514. doi: 10.1029/TC006i004p00505