滇西北地区云县—宁蒗宽角反射/折射剖面结果

陈思文　王宝善*　田晓峰　王夫运　刘宝峰　李　璐

1)中国地震局地震观测与地球物理成像重点实验室(中国地震局地球物理研究所)、北京　100081　2)中国地震局地球物理勘探中心、郑州　450002



滇西北地区云县—宁蒗宽角反射/折射剖面结果

陈思文1)王宝善1)*田晓峰2)王夫运2)刘宝峰2)李璐1)

1)中国地震局地震观测与地球物理成像重点实验室(中国地震局地球物理研究所)、北京1000812)中国地震局地球物理勘探中心、郑州450002

位于滇西地区的红河断裂是中国最长的走滑断裂之一、具有很高的地震潜势。为了调查红河断裂的复杂结构、在滇西地区穿过红河深大断裂带完成了1条由云县至宁蒗近SN向长300km的宽角地震反射/折射剖面。结合初至波走时成像及正演建模方法、对该测线观测数据进行了一维、二维分析拟合、获得了该地区沿测线的二维地壳速度结构模型。结果显示：地壳P波平均速度为6.2～6.3km/s、基本呈现为1个均匀的正速度梯度结构、但在中地壳和下地壳不同区域有部分低速异常；沿测线Moho界面埋深存在较大的横向变化、自南向北明显变深、南侧Moho埋深约为45km、北侧Moho埋深可达54km、较之典型的大陆地壳、存在明显的增厚现象；而沿测线中上地壳厚度变化不大、地壳增厚主要缘于下地壳厚度的增加；以红河断裂为界、地壳结构存在明显的横向差异、暗示了红河断裂作为扬子准地台和三江地槽系构造边界的作用；测线穿过区域红河断裂两侧没有明显的Moho埋深变化、结合临近区域的其他研究结果、表明红河断裂在不同区段存在较大的差异。

二维地壳速度结构宽角反射/折射红河断裂滇西地区

0　引言

滇西地区位于青藏高原东南缘、东邻扬子地块；受印度-欧亚大陆碰撞影响、区内地壳强烈变形、深大断裂纵横交错。其中红河断裂作为青藏东南缘1条醒目的大地构造边界、具有漫长的发育历史和十分复杂的演化过程(Peltzeretal.、1988；Schäreretal.、1990)。红河断裂带东北侧为扬子准地台、西南侧为三江地槽系。在这2个一级构造单元的构造演化和地层发育以及地壳剧烈运动和变形的整个过程中、红河断裂带始终起到了活动块体间相对运动与变形的转换或调节作用(Tapponnieretal.、1990)。根据断裂带两侧地震和震源机制解分析、红河断裂的活动方式主要为走滑型、同时伴有正断型和少许逆冲型、且区域应力反演结果表明红河断裂周围最大主应力轴呈现顺时针旋转的空间变化规律(Zhaoetal.、2013)。而滇西地区的澜沧江断裂带属于深大断裂带、总体走向近SN、在云县北部形成明显的大拐弯(王绍晋等、2007)、是控制滇西地区地质发展的主干断裂。滇西北地区的程海断裂北起宁蒗附近、向SW延伸经程海、宾川、再向S交会于红河断裂带上、是1条新生代强烈活动的左旋走滑兼逆冲断裂(谭筱红等、1999)。滇西地区特殊的构造环境、复杂的地质背景和活跃的地震活动、使得该区成为大陆地球动力学、强震孕育及监视、地壳演化等地学问题的天然试验场。

为了解云南地区深部构造背景、中国在20世纪80年代开展了 “滇深82”和 “滇深86”(Huetal.、1986；Kanetal.、1986)2项研究工程、布设了4条深地震探测测线、获得了地壳上地幔速度结构的一些重要特征。除了地震方法外、科学家们运用各种地质和地球物理方法在此处开展了大量的研究、如利用天然地震波资料开展的体波与面波层析成像(Yaoetal.、2006；Huangetal.、2012)、利用接收函数方法得到地壳S波波速和泊松比(吴建平等、2001；Huetal.、2005)、横波分裂(常利军等、2006；Huangetal.、2007)、噪声成像(Yangetal.、2012)等、这些研究成果在云南地区Moho界面埋深南浅北深、红河断裂是重要的构造边界断裂等方面具有较好的共识。然而、关于该区壳内低速层的分布及其形态、空间展布、红河断裂带不同区段的结构差异等方面的研究仍存在较大争议(吴建平等、2001；Huetal.、2005；Zhangetal.、2009；Huangetal.、2012)。

为进一步揭示研究区域的地壳深部结构、同时为宾川固定气枪发射台开展的人工震源地下介质监测研究提供参考(王宝善等、2011)、我们于2012年6月完成了自云县、经宾川到宁蒗的长约300km的人工地震宽角反射/折射剖面、通过对该深地震测深数据进行分析处理、建立了研究区域沿剖面的二维精细地壳速度模型、并对该区的地壳速度结构进行了研究、对其可能的地质环境和构造背景含义进行了解释、以期深化对本区深部结构及其动力学机制的理解。

1　人工地震测深剖面位置与观测系统

本次深地震宽角反射/折射测线南端始于云南云县(100.14°E、24.45°N)、经过弥渡、宾川、金沙江、永胜、止于宁蒗(101.00°E、27.03°N)、全长约300km(图1)。测线先后穿过三江地槽系和扬子准地台2个一级构造单元、地形高差悬殊、海拔1020～2986m。该测线南端与南定河断裂北段斜交、南段近垂直于澜沧江断裂、无量山断裂和红河深断裂、北段与复杂的程海断裂走向基本一致。因而测线资料的解译适用于研究区不同块体速度结构特征及其演化的研究与理解、也适用于宾川、云县一带深部构造特征的探测。

图1　研究区域地形构造(邓起东、2007)和宽角折/反射剖面位置示意图Fig. 1　Major active tectonics(after DENG Qi-dong、2007)and schematic locations of shot points (yellow pentagrams)and wide-angle reflection/refraction seismic profile in the study area.黑色实线为活动断层：F1南定河断裂、F2澜沧江断裂、F3无量山断裂、F4红河断裂、F5程海断裂；红色和白色圆分别代表1515—2007年5级以上地震的震中位置和测线检波器位置、黄色星号代表炮点位置；A 三江地槽系、B 扬子准地台、C 松潘-甘孜褶皱系、D 左贡-耿马褶皱系；蓝色虚线代表 “滇深82”、“滇深86”中与本文剖面有交会的测线位置(L1 宾川—遮放测线，　　　　　　　 L2 洱源—江川测线、L3 中甸—思茅测线；本文剖面与L1、L2相交的位置为宾川地区、与L3交于景云桥和大仓地区)

表1　炮点参数一览表

Table1　Shot parameters list

炮点编号桩号/km2012年爆炸时间(月-日时:分:秒)炮点坐标药量/kg井数(口)平均井深/m高程/m岩性炮点地名经度纬度SP1107.22606-1801:20:15.423100°11.363'24°30.464'2793660.931220耕土云南省临江市云县爱华镇永胜村SP2253.52306-1901:20:15.540100°30.041'25°48.330'987364.501569风化岩云南省大理市宾川县大银镇甸尾村SP5375.59006-1801:00:15.939100°59.990'26°49.303'3003763.712986花岗岩云南省丽江市宁蒗县战河乡万河村

本次探测沿剖面分别在云县、宾川、宁蒗实施3次吨级井下爆破激发(表1)。这几次爆破的炸药量不等(987～3303kg)、炮点激发单井深度为60～65m、共布设150个PDS-2轻便数字地震仪、观测点距约2km、构成了较为简单的追逐相遇观测系统(图2)。为了消除和减轻人为活动造成的干扰、本次探测中所有地震波激发和接收工作都在凌晨时分进行。

图2　滇西地区云县—宁蒗深地震测深剖面观测系统图Fig. 2　Observation layout of wide-angle reflection/refraction profile of Yunxian-Ninglang in West Yunnan.

2　地震测深剖面震相特征

从炮点记录截面图中标定震相、并对每个震相进行分析是人工地震测深中尤为基础和重要的工作。不同的震相性质和差异都与测区区域构造演化密切相关、是不同地壳结构构造特征的体现(嘉世旭等、1995)。云县—宁蒗剖面提供的地震记录信息包括地壳结晶基底的回折Pg、地幔的反射PmP、折射Pn波组震相、壳内反射P1P、P2P、P3P波组震相(图3)。各震相的具体特征如下：

在3次爆破剖面上(图3)、初至波Pg和莫霍面反射波PmP震相清楚、可直接辨认。除此之外、还有3组来自壳内的反射波P1P、P2P、P3P、其能量虽较弱、但仍可辨认。另外、作为上地幔顶部的折射波Pn、在南北2炮记录中、其位于红河断裂北侧的震相记录相对清晰。

图3　实测地震记录截面及震相拾取图(折合速度 6.0km/s)Fig. 3　Seismic data from the three shots recorded along the profile(with reduced velocity of 6.0km/s). a SP1炮；b SP2炮；c SP5炮；黑色虚线表示拾取的走时

剖面上折射震相的主要特征：Pg波的可观测范围在测线南段为80km、北段可达150km；振幅、走时变化均较大、说明上地壳的速度结构较复杂、也显示了南北2段上地壳结构的明显差异。图3 SP1和SP2记录截面上Pg存在明显的走时变化、暗示了红河断裂两侧在上地壳存在较大的速度变化。Pn震相在SP1和SP5上约220km处开始出现、视速度为 8.05km/s。

壳内反射波有P1P波组：出现在Pg波之后、为来自上地壳底界面C1的反射波、能量较强、代表界面深度14～17km；P2P波组：为来自中地壳底界面C2的反射波、能量较P1P稍弱、代表界面深度21～24km；P3P波组：出现在P2P波组之后、为来自下地壳中C3界面的反射波、能量较强、代表界面深度27～31km。莫霍面反射波PmP波的特点是振幅强、容易辨认。

SP1记录截面上(图3a)、Pg波从炮点附近至100km 左右为初至震相、视速度4.0～6.0km/s。壳内反射震相P1P、P2P、P3P在此剖面50～190km 的震中距上均能较清楚地分辨；P1P波震相较弱、追踪距离较短、在100km后迅速衰减、其中红河断裂附近Pg和P1P存在明显的速度变化、暗示红河断裂错断了中上地壳。P2P和P3P震相相对较强、追踪距离较长。莫霍面反射波PmP在震中距110～230km范围内可连续追踪。上地幔顶部折射波Pn在215km 后能较好地分辨、其视速度约为 8.05km/s。

SP2记录截面上(图3b)、在南侧、Pg波从炮点附近至90km左右为初至震相；在北侧、从炮点到80km 左右为初至震相、视速度3.5～6.0km/s、其中南侧Pg在近炮点视速度低于北侧和SP1、显示了上地壳结构的横向差异。在桩号300km附近的金沙江一带、Pg走时曲线出现明显的转折变化、显示了金沙江一带上地壳结构存在较大的横向变化。壳内反射震相P1P、P2P、P3P在此剖面40～150km 的震中距上均能较清楚地分辨、P1P和P2P震相较SP1可连续追踪的距离稍长、且北支P2P的速度较高。在红河断裂附近、该记录截面上Pg和P1P同样存在明显的速度变化、进一步表明红河断裂有错断中上地壳的迹象。莫霍面反射波PmP在震中距110～120km范围出现、可连续追踪到测线两端。由于剖面长度的限制、SP2记录截面上未能观测到Pn。

SP5记录截面上(图3c)、在南侧、Pg波从炮点附近至90km左右为初至震相、在桩号300km附近金沙江一带、Pg走时曲线出现明显的转折变化、显示了金沙江一带上地壳结构的横向变化。在北侧、只有25km长的短距离接收、视速度3.5～6.0km/s。壳内反射震相P1P、P2P、P3P在此剖面50～180km 的震中距上均能较清楚地分辨、P1P波震相可在距炮点50～150km的范围内追踪、P2P和P3P震相的追踪距离也比SP1和SP2的长。莫霍面反射波PmP震相较弱、在震中距170km左右出现、可连续追踪到测线南段超过250km、这种现象可能与测线北段莫霍面埋深加深有关。上地幔顶部折射波Pn在220km 左右出现、其视速度约为 8.05km/s、能量弱于SP1。但Pn波震相主要出现在红河断裂带以北的现象在过去的研究中也有体现(阚荣举等、1986)。

3　二维地壳速度模型

通过上述主要地震震相分析、我们拾取了相应的走时、并通过试错法对所有单炮记录进行拟合、得到了一维地壳结构模型(图4)。在此基础上、我们利用射线追踪(Cerveny、1979；Cerveny、1984)计算了理论走时、并将其作为二维地壳结构模型建立的初始模型。

结合沿测线的隆起、坳陷、断裂等实际地质情况和已有的地质研究成果、我们在一维模型的基础上、建立了详细的二维初始地壳速度模型、再将沿测线的地表高程直接输入模型参与计算、然后对各炮点观测资料的反/折射震相走时、振幅通过动力学射线追踪Seis83程序(Cerveny、1979、1984)完成二维速度模型的正演拟合计算。在拟合过程中、我们遵循由浅到深依次调整拟合上地壳、中地壳、下地壳的速度结构和界面深度的方法、反复使用Seis83程序进行理论计算、最终实现理论走时与实际到时的充分拟合。本研究中、基于最终的二维速度模型得到的理论走时、振幅与观测资料相比、大多数点拟合误差在±0.1s左右。对于个别误差较大的点、有些可能是由于炮点偏离测线或测线本身的弯曲造成的；有些则是由于地质结构的复杂性造成的。综合考虑以上各种因素、最终我们在一定的误差范围内选取1个合理的速度模型。图5—7所示分别为SP1 炮、SP2炮、SP5炮相对应的射线路径及其相应的走时拟合效果。由图5—7可以看出、基于我们的二维速度模型计算得到的理论走时、各波组的振幅等特征与实测资料的记录特征达到较好的拟合。

表2　用于反演约束地壳模型的地震震相记录及参数

Table2　Number of traveltimes picked for each phase from each shot

震相SP1SP2SP5总数Pg367553164P1P165341110P2P344637117P3P405233125Pm413029100Pn30—1444

注“—”表示无此震相。

图4　基于1D速度模型正演得到的SP1、SP2、SP5炮走时拟合图Fig. 4　Computed travel time data(continuous lines)of the three shots based on 1-D velocity model.圆圈代表实测走时、曲线代表理论走时；“-”和 “+”分别代表SP2南支和北支

图5　SP1的走时拟合(a)和2D射线追踪(b)Fig. 5　Travel time fitting(a)and 2-D ray tracing(b)for SP1.圆圈代表实测走时；其他符号代表不同震相的理论走时

图6　SP2的走时拟合(a)和2D射线追踪(b)Fig. 6　Travel time fitting(a)and 2-D ray tracing(b)for SP2.圆圈代表实测走时；其他符号代表不同震相的理论走时

图7　SP5的走时拟合(a)和2D射线追踪(b)Fig. 7　Travel time fitting(a)and 2-D ray tracing(b)for SP5.圆圈代表实测走时；其他符号代表不同震相的理论走时

本研究中、共使用了660条走时记录(表2)、其中164条震中距在150km范围内的初至到时记录用于约束顶部2层的速度结构及第1个速度界面的深度；用于约束壳内不同界面C1、C2和C3深度的反射波组、P1P、P2P、P3P到时数据分别有110、117和125个；100条PmP反射到时用于约束莫霍界面的形态；44条Pn折射到时用于约束上地幔顶部的速度结构；此外、中地壳和下地壳不同层位的速度由相应的震相远端折射到时来约束。

4　地壳上地幔结构与构造特征

研究区地下构造十分复杂、破碎程度高(Huangetal.、2012)。从得到的地壳介质地震波速度模型(图8)和射线追踪图(图5-7)可以看出、该剖面地壳厚度和壳内反射射线分布等在横向上具有强烈的非均匀性、主要界限位于红河断裂、且沿剖面所跨过的断裂都在速度结构上有不同程度的体现。我们根据得到的速度结构将研究区域的地壳分为3层、其中、反射波P1P确定的C1界面以上的部分看作上地壳、由P3P波组的反射深度确定的C3界面和C1界面之间的层位称为中地壳、C3界面到PmP反射波所控制的莫霍面之间的部分称为下地壳。

图8　云县—宁蒗地震宽角折/反射剖面二维地壳P波速度结构Fig. 82-D P wave velocity structure of the crust and upper-mantle along the Yunxian-Ninglang wide-angle seismic profile.黑色虚线表示对应断裂可能的倾向、红色星形表示永胜发生过的2次7级以上大地震的位置；2条红色虚线之间为中下地壳有效分辨区段

上地壳结构：该层指自地表向下至C1界面之间的地层部分。包括基底回折波Pg控制的基底面和壳内反射波P1P控制的C1界面、共同反映了上地壳结构。由Pg波反演得到的基底面形态起伏大、深度在2.3～5.3km之间变化、宾川基底面厚度为4～5km、与白志明(2003)等宾川—遮放、洱源—江川测线所做层析成像的结果一致。基底面在剖面南段桩号130～145km和160～180km处变深4～5km、这与澜沧江断裂、无量山断裂和红河断裂形成一定程度的对应；在与中甸—思茅测线相交的景云桥地区、本文剖面下方基底面深度为3.5～4.2km、稍深于林中洋等人(1993)2.7～4.0km的结果；向剖面北段延伸、基底面在340～350km处明显变浅、埋深约2.5km左右。整个基底面到地表的速度结构形成1个强速度梯度层、由近地表面的2.9～4.0km/s到基底面之上的5.2～5.7km/s。宾川和景云桥地区该层的P波速度为 5.4km/s左右、比以往测线相同位置的Pg波速度值略小。地表速度和剖面上的覆盖层厚层相关、剖面北段基岩出露较多、地层速度较高。

由P1P反射波控制的C1界面除了在红河断裂处(桩号170～180km)深度达到16.8km外、起伏较基底面小、基本稳定在14.5～15.8km左右。该层的P波速度为6.02～6.23km/s、横向速度变化显著。红河断裂南侧、P波速度约为 6.05km/s、与中甸—思茅的速度结果一致(林中洋等、1993)；断裂北侧速度发生较大跃升、至 6.25km/s左右。这种速度的横向大幅度跃变在澜沧江断裂、无量山断裂、程海断裂处都分别有不同程度的体现。C1界面的速度跳跃为6.08～6.26km/s到6.23～6.36km/s、跳跃差在0.10～0.20km/s。

中地壳结构：该层是指C1与C3界面之间的层位、中间包括由反射P2P波控制的C2界面、深度约在22km左右。由反射波P2P和P3P的走时约束着中地壳的界面深度。限于射线覆盖、二维模型中地壳有效分辨范围在桩号130～370km之间。从界面形态看、C3界面与C1界面起伏变化基本一致。从速度结构来看、该层表现为1弱速度梯度层、上下边界速度值分别为6.20km/s和 6.55km/s、与 “滇深82”、“滇深 86”工程中的3条测线(图1)结果相同(Huetal.、1986；阚荣举等、1986; 林中洋等、1993；Zhangetal.、2009)。从横向速度变化来看、沿剖面无量山断裂以北至红河断裂以南、以及走向与剖面几乎一致的程海断裂的北端存在着低速异常、在横向上速度与前后层相差约 0.1km/s。该层埋深27.1～31.0km、较其他3条测线稍深(Huetal.、1986；林中洋等、1993；Zhangetal.、2009)；厚度为12.5～16.5km。与上地壳相比、中地壳继承了上地壳层中的低速异常、但厚度变化不如上地壳明显。

下地壳结构：从C3界面到莫霍面之间的层位为下地壳、其由壳内反射波P3P和莫霍面反射波PmP控制。限于射线覆盖、二维模型中地壳有效分辨范围在桩号150～350km之间。该层的界面形态呈明显的南浅北深的趋势、南段最浅处为45.0km、往北逐渐加深至54.0km、层厚度变化幅度较大、为16.2～23.0km；而宾川以及景云桥附近以往研究的莫霍面深度约为42～45km(阚荣举等、1986；林中洋等、1993)、本文南侧剖面结果稍深。该层的P波速度变化范围为6.49～6.90km/s、比3条相交区域的速度研究结果稍低(阚荣举等、1986；林中洋等、1993)。P波速度结构在下地壳基本呈现1个较均匀的正速度梯度层、但在红河断裂南侧横向速度存在1个弱的低速异常、而这种低速异常以更弱的方式出现在程海断裂北端。限于观测系统、下地壳射线覆盖较为稀疏、本研究所得下地壳结构复杂程度较低。然而红河断裂两侧的速度差异和程海断裂北段的低速异常仍在有效射线覆盖范围内。

上地幔顶部：本剖面SP1和SP5均有Pn记录；SP2因为剖面长度原因、未记录到Pn波。结合射线追踪、走时拟合及合成地震图显示地幔P波速度随深度增加而加大、上地幔顶部为 8.0km/s左右、比Hu等(1986)、阚荣举等(1986)等人的研究结果稍低。

5　讨论与结论

剖面地壳厚度呈现西南浅东北深的趋势、南起点云县地壳厚度为45km、北侧地壳深至54km、地壳厚度的这种变化趋势与滇西地区以往的多数地球物理探测结果一致(Huetal.、1986；林中洋等、1993；张中杰等、2005；李永华等、2009；Zhangetal.、2009)。在测线中段的宾川地区、莫霍面深度与Hu(1986)的结果一致、但稍深于白志明(2003、2004)有限差分反演的结果。从剖面结果看、上地壳和中地壳厚度变化不大、分别为14～16km和12～16km、剖面地壳增厚主要归因于下地壳的厚度增加、西南侧下地壳厚度为16km左右、而东北侧厚度高达23km。这种主要由下地壳增厚引起的地壳总厚度变化现象在已有的人工地震测深成果中屡有体现(Huetal.、1986；林中洋等、1993；张中杰等、2005)；同时张忠杰等(2005)认为下地壳增厚的原因是多重的、可能与特提斯洋壳俯冲、消减、碰撞造山和壳幔一系列复杂的相互作用有关。同时、莫霍面的连续可能在超临界震中距附近产生大振幅的SmS 震相、对于此区域的强地面震动具有放大效应(罗艳等、2010)。

速度结构结果显示、整个剖面中、速度横向变化差异最显著的地方位于红河断裂带。红河断裂是1条长期演化发展的深大断裂、至少始于前古生代、在新生代早期、受喜马拉雅运动的影响、它主要表现为左旋走滑、晚新生代时慢慢转换为右旋走滑活动(张建国、2009)。澜沧江断裂和红河断裂之间的上中地壳出现明显的低速异常、其两侧的速度差异十分明显、这与以往体波成像(Huangetal.、2012)和接收函数(Huetal.、2005；李永华等、2009)研究得到的红河断裂两侧P、S波速度甚至波速比都存在明显差异一致。红河断裂西南侧为与沉积盆地对应的低速异常区、东北侧则表现为高速异常、这可能与地表出露的玄武岩有关；白志明(2003)等也认为此处有玄武岩隆起。红河断裂带两侧的速度差约为 0.2km/s、一直延伸到下地壳、从速度差异推断红河断裂在地壳范围内倾向SW、层析成像结果表明红河断裂下方上地幔波速异常的分界有向SW移动的趋势(Huangetal.、2012)、与本文地壳内的错动倾斜趋势一致。遮放—宾川、孟连—马龙测线上红河断裂上地壳内的倾向也与本文相同、但其判断红河断裂深部为NE倾向(白志明等、2003)、这点在本次剖面速度结果中没有体现。

从速度差异判断断裂的切割深度分析、红河断裂在本条速度结构剖面上显示为1条深大断裂、至少切割至下地壳。本研究观测数据中、反映中下地壳结构的P2P和P3P在红河断裂附近走时曲线均有变化、为红河断裂切割至下地壳提供了直接证据。由于本研究观测系统较为稀疏、红河断裂两侧PmP覆盖不足、难以判定红河断裂是否错断莫霍面。在以往开展的一些人工地震测线中、徐鸣洁等(2005)根据接收函数反演结果中红河断裂下方地壳厚度存在突变推断红河断裂是陡立的超壳断裂；白志明等(2003)研究孟连—马龙DDS剖面的结果暗示红河断裂是切割地壳的深大断裂；也有层析成像结果显示上地幔存在波速差异(Huangetal.、2012)。考虑到观测系统本身的不足、本文研究地区的速度结构中红河断裂下方上地幔顶端没有明显的横向速度差异、缺乏判断其为超壳断裂的直接证据。但由于滇西地区地壳厚度自南向北逐渐加深、加之红河断裂中新世末以来存在2次右旋正断位错事件(向宏发等、2006)、且断裂中段为典型的右旋走滑变形区(向宏发等、2004)、以及滇西地区可能存在下地壳物质流动等多重地质构造因素(Roydenetal.、1997；朱守彪等、2004)、使红河断裂深部运动十分复杂。因此、即使在红河断裂下方上地幔顶部没有明显的横向速度差异、且反演结果显示断裂两侧Moho界面连续、也并不能说明红河断裂没有深切上地幔顶部这一问题。所以、红河断裂在云县—宁蒗剖面相交区域是否为超壳断裂还有待进一步验证、这也从侧面反映了红河断裂各区段地壳上地幔速度结构存在较大差异的可能。

由本文剖面研究区域最终的速度结构可以得出以下主要结论：研究区域地壳P波平均速度为6.2～6.3km/s；红河断裂两侧的地壳结构存在很大的差异、横向速度变化显著、进一步说明其为1条大地构造板块边界；地壳由西南至东北逐渐加深、且地壳增厚主要来自下地壳的增厚；测线北段金沙江与永胜地区之间地壳内有局部低速异常、这是由于程海断裂走向与测线布设方向近一致、令此段剖面所在的区域地壳破碎程度高、壳内出现的较大范围的低速异常反映了程海断裂下方极其复杂的地壳结构。

致谢本研究得到了中国地震局郑州物探中心杨卓欣研究员、潘素珍高级工程师和诸多同仁的悉心帮助和指导、以及王椿镛研究员和张先康研究员在分析解释方面给予的帮助、在此表示深深的谢意。

白志明、王椿镛. 2003. 云南地区上部地壳结构和地震构造环境的层析成像研究 [J]. 地震学报、25(2): 117—127.

BAI Zhi-ming、WANG Chun-yong. 2003. Tomographic investigation of the upper crustal structure and seismotectonic environments in Yunnan Province [J]. Acta Seismologica Sinica、25(2): 117—127(in Chinese).

白志明、王椿镛. 2004. 云南遮放—宾川和孟连—马龙宽角地震剖面的层析成像研究 [J]. 地球物理学报、47(2): 257—267.

BAI Zhi-ming、WANG Chun-yong. 2004. Tomograghy research of the Zhefang-Binchuan and Menglian-Malong wide-angle seismic profiles in Yunnan [J]. Chinese J Geophys、47(2): 257—267(in Chinese).

常利军、王椿镛、丁志峰. 2006. 云南地区 SKS 波分裂研究 [J]. 地球物理学报、49(1): 197—204.

CHANG Li-jun、WANG Chun-yong、DING Zhi-feng. 2006. A study on SKS splitting beneath the Yunnan region [J]. Chinese J Geophys、49(1): 197—204(in Chinese).

邓起东. 2007. 中国活动构造图(1︰400万)[CM]. 北京：地震出版社.

DENG Qi-dong. 2007. Active Tectonics Map of China(1︰4,000,000)[CM]. Seismological Press、Beijing(in Chinese).

嘉世旭、刘昌铨. 1995. 华北地区人工地震测深震相与地壳结构研究 [J]. 地震地质，17(2): 97—105.

JIA Shi-xu、LIU Chang-quan. 1995. Study on the seismic phase of DSS in North China [J]. Seismology and Geology、17(2): 97—105(in Chinese).

阚荣举、林中洋. 1986. 云南地壳上地幔构造的初步研究 [J]. 中国地震、2(4): 50— 61.

KAN Rong-ju、LIN Zhong-yang. 1986. A preliminary study on crustal and upper mantle structures in Yunnan [J]. Earthquake Research in China、2(4): 50— 61(in Chinese).

李永华、吴庆举、田小波、等. 2009. 用接收函数方法研究云南及其邻区地壳上地幔结构 [J]. 地球物理学报、52(1): 67—80.

LI Yong-hua、WU Qing-ju、TIAN Xiao-bo、etal. 2009. Crustal structure in the Yunnan region determined by modeling receiver functions [J]. Chinese J Geophys、52(1): 67—80(in Chinese).

林中洋、张文彬. 1993. 滇西地区地壳上地幔速度结构特征的研究 [J]. 地震学报、15(4): 427— 440.

LIN Zhong-yang、ZHANG Wen-bin. 1993. Studies of upper mantle velocity structure in Western Yunnan [J]. Acta Seismologica Sinica、15(4): 427— 440(in Chinese).

罗艳、倪四道、陈颙. 2010. 超临界距离上SmS震相对强地面振动衰减关系影响的初步探讨 [J]. 中国地震、26(2): 131—141.

LUO Yan、NI Si-dao、CHEN Yong. 2010. A case study on the effect of post-critical SmS on ground motion in Yingjiang、Yunnan Province、China [J]. Earthquake Research in China、26(2): 131—141.

谭筱红、朱志. 1999. 程海断裂新生代活动特征及性质 [J]. 云南地质科技情报、(2): 21—24.

TAN Xiao-hong、ZHU Zhi. 1999. Characteristics and nature of Cenozoic events on Chenghai Fault、Yunnan [J]. Geological Science and Technology Information、(2): 21—24(in Chinese).

王宝善、王伟涛、葛洪魁、等. 2011. 人工震源地下介质变化动态监测 [J]. 地球科学进展、(3): 249—256.

WANG Bao-shan、WANG Wei-tao、GE Hong-kui、etal. 2011. Monitoring subsurface changes with active sources [J]. Advances in Earth Science、(3): 249—256(in Chinese).

王绍晋、秦嘉政、李忠华、等. 2007. 澜沧江断裂带环境剪应力场与地震活动分析 [J]. 地震研究、30(2): 127—132.

WANG Shao-jin、QIN Jia-zheng、LI Zhong-hua、etal. 2007. Relationship between seismicity and environmental shear stress field along Lancangjiang fault zone [J]. Journal of Seismological Research、30(2): 127—132(in Chinese).

吴建平、明跃红、王椿镛. 2001. 云南数字地震台站下方的 S 波速度结构研究 [J]. 地球物理学报、44(2): 228—237.

WU Jing-ping、MING Yue-hong、WANG Chun-yong. 2001. The S-wave velocity structure beneath digital seismic stations of Yunnan Province inferred from teleseismic receiver function modeling [J]. Chinese J Geophys、44(2): 228—237(in Chinese).

向宏发、韩竹军、虢顺民、等. 2004. 红河断裂带大型右旋走滑运动定量研究的若干问题 [J]. 地球科学进展、19(增刊)：56—59.

XIANG Hong-fa、HAN Zhu-jun、GUO Shun-min、etal. 2004. Several questions about the quantitative study of the large-scale right-lateral strike-slip movement of the Red River fault zone [J]. Advances in Earth Science、19(Suppl): 56—59(in Chinese).

向宏发、万景林、韩竹军、等. 2006. 红河断裂带大型右旋走滑运动发生时代的地质分析与FT测年 [J]. 中国科学(D 辑)、36(11): 977—987.

XIANG Hong-fa、WAN jing-lin、HAN Zhu-jun、etal. 2006. Geological analysis and FT dating of the large-scale right-lateral strike-slip movement of the Red River fault zone [J]. Science in China(Ser D)、36(11)：977—987(in Chinese).

徐鸣洁、王良书、刘建华、等. 2005. 利用接收函数研究哀牢山-红河断裂带地壳上地幔特征 [J]. 中国科学(D辑)、35(8): 729—737.

XU Ming-jie、WANG Liang-shu、LIU Jian-hua、etal. 2005. Study on crust and upper mantle structure in Ailaoshan-Honghe faults zone from receiver functions [J]. Science in China(Ser D)、35(8): 729—737(in Chinese).

俞维贤、张建国、周光全、等. 2005. 2001年永胜6级地震的地表破裂与程海断裂 [J]. 地震研究、(2): 125—128.

YU Wei-xian、ZHANG Jian-guo、ZHOU Guang-quan、etal. 2005. Surface rupture of the 2001 YongshengM6 earthquake and Chenghai Fault [J]. Journal of Seismological Research、(2): 125—128(in Chinese).

张建国. 2009. 中越红河断裂活动性研究 [D]: [学位论文]. 合肥：中国科学技术大学.

ZHANG Jian-guo. 2009. Research on the activity of Red River Fault in China and Vietnam [D]. Dissertation. University of Science and Technology of China, Hefei (in Chinese).

张中杰、白志明、王椿镛、等. 2005. 冈瓦纳型和扬子型地块地壳结构: 以滇西孟连—马龙宽角反射剖面为例 [J]. 中国科学(D 辑)、34(5): 387—392.

ZHANG Zhong-jie、BAI Zhi-ming、WANG Chun-yong、etal. 2005. Crustal structure of Gangwana and Yangtze-typed blocks: An example by wide-angle seismic profile from Menglian to Malong in Western Yunnan [J]. Science in China(Ser D)、34(5): 387—392(in Chinese).

朱守彪、石耀林. 2004. 用遗传有限单元法反演川滇下地壳流动对上地壳的拖拽作用 [J]. 地球物理学报、47(2): 232—239.

ZHU Shou-biao、SHI Yao-lin. 2004. Genetic algorithm-finite element inversion of drag forces exerted by the lower crust on the upper crust in the Sichuan-Yunnan area [J]. Chinese J Geophys、47(2): 232—239(in Chinese).

Cerveny V. 1979. Ray theoretical seismograms for laterally inhomogeneous structures [J]. J Geophys、46: 335—342.

Cerveny V. 1984. SEIS83-Numerical modeling of seismic wave fields in 2-D laterally varying layered structures by the ray method [M]. In: Engdahl E R (ed). Documentation of Earthquake Algorithms. World Data Center A for Solid Earth Geophysice, Colorado, USA.

Hu Hong-xiang、Lu Han-xing、Wang Chun-yong、etal. 1986. Explosion investigation of the crustal structure in western Yunnan Province [J]. Chinese J Geophys、29(2): 133—144.

Hu Jia-fu、Su You-jin、Zhu Xiong-guan、etal. 2005. S-wave velocity and Poisson’s ratio structure of crust in Yunnan and its implication [J]. Science in China(Ser D)、48(2): 210—218.

Huang Jing、Liu Xue-jun、Su You-jing、etal. 2012. Imaging 3-D crustal P-wave velocity structure of western Yunnan with bulletin data [J]. Earthquake Science、25(2): 151—160.

Huang Zhou-chuan、Wang Liang-su、Xu Ming-jie、etal. 2007. Shear wave splitting across the Ailao Shan-Red River fault zone、SW China [J]. Geophysical Research Letters、34(20): L20301.

Kan Rong-ju、Hu Hong-xiang、Zeng Rong-sheng、etal. 1986. Crustal structure of Yunnan Province、People’s Republic of China、from seismic refraction profiles [J]. Science、234: 433— 437.

Peltzer G、Tapponnier P. 1988. Formation and evolution of strike-slip faults、rifts、and basins during the India-Asia collision: An experimental approach [J]. Journal of Geophysical Research: Solid Earth(1978-2012)、93(B12): 15085—15117.

Royden L H、Burchfiel B C、King R W、etal. 1997. Surface deformation and lower crustal flow in eastern Tibet [J]. Science、276: 788—790.

Schärer U、Tapponnier P、Lacassin R、etal. 1990. Intraplate tectonics in Asia: A precise age for large-scale Miocene movement along the Ailao Shan-Red River shear zone、China [J]. Earth and Planetary Science Letters、97(1): 65—77.

Tapponnier P、Lacassin R、Leloup P H、etal. 1990. The Ailao Shan/Red River metamorphic belt: Tertiary left-lateral shear between IndoChina and South China [J]. Nature、343: 431— 437.

Yang Y J、Ritzwoller M H、Zheng Y、etal. 2012. A synoptic view of the distribution and connectivity of the mid-crustal low velocity zone beneath Tibet [J]. Journal of Geophysical Research: Solid Earth(1978-2012)、117(B4).

Yao H J、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-I. Phase velocity maps [J]. Geophysical Journal International、166(2): 732—744.

Zhang Xi、Wang Yang-hua. 2009. Crustal and upper mantle velocity structure in Yunnan、Southwest China [J]. Tectonophysics、471(3): 171—185.

Zhao Li、Luo Yan、Liu Ting-Yu、etal. 2013. Earthquake focal mechanisms in Yunnan and their inference on the regional stress field [J]. Bulletin of the Seismological Society of America、103(4): 2498—2507.

Abstract

The Red River Fault in western Yunnan is one of the longest strike-slip faults in China and has a high seismic potential. To investigate its complicated structure、a near-NS directed 300km long wide-angle reflection/refraction seismic profile was laid out from Yunxian to Ninglang、across the Red River Fault. The 2-D velocity structure model along the profile was obtained through 1-D and 2-D analysis and fitting the observed data with combination of first-arrival traveltime tomography and forward modeling. The results indicate: In the crust、the average P-wave velocity is 6.2～6.3km/s and basically shows a positive gradient structure、but there are some low velocity anomalies at different area in upper and lower crust. Regarding the crust boundary、a relative large lateral variation exists in the depth of Moho、which goes deeper from south to north、ranging from 45km to as deep as 54km; compared to other typical continental crust、the study area demonstrates a striking thickening. It should be mentioned that the crustal thickening is mainly observed in the lower crust、while the upper and middle crust possess nearly constant thickness. We observed strong seismic velocity contrast across the Red River Fault、which emphasizes the role of the fault as an important tectonic boundary between Yangtze paraplatform and Sanjiang geosynclinal system. Along the profile、the Moho depth has no remarkable variation when crossing the Red River Fault. Combining with other study results on nearby area、it proves that there is notable heterogeneity between different parts of the Red River Fault.

CRUSTAL STRUCTURE FROM YUNXIAN-NINGLANG WIDE-ANGLE SEISMIC REFLECTION AND REFRACTION PROFILE IN NORTHWESTERN YUNNAN、CHINA

CHEN Si-wen1)WANG Bao-shan1)TIAN Xiao-feng2)WANG Fu-yun2)LIU Bao-feng2)LI Lu1)

1)KeyLaboratoryofSeismicObservationandGeophysicalImaging,InstituteofGeophysics，ChinaEarthquakeAdministration、Beijing100081、China2)GeophysicalExplorationCenter、ChinaEarthquakeAdministration、Zhengzhou450002、China

western Yunnan、two￣dimensional crustal velocity structure、wide-angle reflection/refraction、Red River Fault

10.3969/j.issn.0253-4967.2016.01.007

2014-10-09收稿、2015-10-02改回。

地震行业科研专项(201208004)、中国地震局地球物理研究所基本科研业务专项(DQJB10A01)、国家自然科学基金(41222029、41104038)、云南省陈颙院士工作站(20141C007)与新疆维吾尔自治区科研机构创新发展专项资金(2013016)共同资助。

王宝善、研究员、E-mail：wangbs@cea-igp.ac.cn。

P313.2

A

0253-4967(2016)01-0091-16

陈思文、女、1988年生、2015年毕业于中国地震局地球物理研究所、获硕士学位、主要从事地壳主动源探测方面的研究工作、E-mail：chensw1228@gmail.com。