2. 中国地震局地壳应力研究所, 北京 100085;
3. 四川省地震局测绘工程院, 四川雅安 625000
2. Institute of Crustal Dynamics, China Earthquake Administration, Beijing 100085, China;
3. Survey Engineering Institute, Earthquake Administration of Sichuan Province, Ya'an Sichuan 625000, China
鲜水河断裂带地处青藏断块区东北边缘, 南北地震带的中段, 其北为巴颜喀拉块体, 其南为川滇块体, 是一条左行走滑的全新世活动断裂带[1], 全长300余公里, 总体走向N40 ° -50 °W, 具有向北东方向的微弯曲, 倾向主要为NE, 小部分为SW, 倾角较陡, 一般在70 °左右, 其北西段是由老乾宁向北西延伸的部分, 形迹相对简单和连续, 属于“单一型”几何结构[2], 内部可细分为左阶斜列的炉霍段、道孚段、乾宁段(图 1), 其中炉霍段1973年2月6日发生过炉霍7.6级地震, 道孚段1981年1月24日发生过道孚6.9级地震; 南东段(老乾宁以南)则相对复杂, 由雅拉河、康定(色拉哈)、折多塘、磨西四条分支断层组成.20世纪70年代开始, 我国地震工作者沿该断裂带的主干断裂进行跨断层形变观测, 早期为水准和短基线测量, 至80年代末期, 又相继布设了跨断层定点连续蠕变测点, 积累了大量宝贵的断层形变资料.近十多年来, 鲜水河周边地区发生多次强震, 其中包括1997年玛尼MS7.5级地震, 2001年昆仑山MS8.1级大地震, 2008年汶川MS8.0级大地震, 这几次地震均发生在巴颜喀拉块体的边界断裂带上, 鲜水河断裂带作为巴颜喀喇块体的西南边界, 在这一时期发生了什么样的变化, 对于区域动力学研究意义重大, 需要通过断层形变观测数据深入分析.
鉴于鲜水河断裂带在青藏高原内部变形中所起的关键作用, 前人的研究涉及地震地质条件、地壳形变观测、地震危险性等方面.在空间展布上, 该断裂带具有雁列组合和断层弯曲的几何特征[3].滑动速率方面, 由地貌断错分析得到的近代地质滑动速率介于5~15mm / a [4-7], GP S观测到的现今滑动速率为6~11mm / a [8-10], 近场的短基线测量显示其北段为1.04~3.82mm / a, 中南段最低0.08mm / a [11].诸多研究显示以老乾宁为分界点(图 1), 断裂带的北西段与南东段在地貌、深部结构、滑动速率及地震活动性方面具有分段性, 北西段强震活动频度高, 走滑速率相对较大[12]; 南东段断裂出现分支, 地貌上表现为隆起区[13], 走向发生顺时针偏转, 地震频度相对北西段低, 走滑速率小于北西段.岩石的破碎程度[14]、断层泥的粒径分布[15]、地壳厚度以及重磁异常[14]的分布也显示其分段差异明显.此外, 也有学者从短周期蠕滑事件频次差异[16]、断层蠕变的分维值[17]的角度探讨了该断裂带的活动特征.
如上所述, 前人已对鲜水河断裂带的分段性进行过研究, 但大多是基于断层几何方面的讨论, 对断层形变观测反映的分段活动关注较少.近十几年来, 该断裂带周边强震频发, 对其断层活动可能存在影响.2008年, 巴颜喀拉块体东南边界的龙门山断裂曾发生逆冲兼右旋走滑的汶川8.0级大地震[18], 它是巴颜喀拉地块整体S EE向运动在龙门山推覆构造带上弹性加载以至释放的结果[19], 然而, 鲜水河断裂带长期以来均为左行走滑, 在运动学上与巴颜喀拉块体S EE向的运动相互矛盾, 汶川地震前鲜水河断裂带理应出现右行走滑或左行走滑的减弱, 这样才能为汶川地震积累弹性应变能, 这个疑问有待实际的断层形变资料去证实.以往的研究显示, 不同观测方法得到的鲜水河断裂带走滑速率差异显著, 那么, 断裂带两侧的变形场是如何分布的?
针对上述问题, 本文利用鲜水河断裂带1990年1月-2009年12月的蠕变与短基线数据, 采用小波变换与断层运动学分析方法, 获取构造活动产生的断层形变速率.结合近场断层形变测量与GP S资料分析该断裂带的分段活动特征及时空演化.
2 鲜水河断裂带断层形变观测概况 2.1 蠕变观测断层蠕变仪是一种监测断层活动的有力工具[20-22], 常见杆式、线式和梁式三种[23], 本质上均为断层近场的基线测量.鲜水河断裂带的6个蠕变测点由北向南依次为虾拉沱、恰叫、沟普、龙灯坝、老乾宁和紫马垮(图 1), 采用DS J水平蠕变仪, 产出数据为日值, 目前已经积累了20余年的历史数据. DS J仪器应用比较法原理, 采用超低膨胀系数、高抗拉张强度的含Nb铟瓦合金丝, 以确定张力下形成的弦长为长度基准, 观测两个被测参考点之间水平距离的微量变化, 仪器跨距一般为3~50m, 机械式仪器灵敏度为0.002mm, 光电式为0.001mm, 量程为30mm(可扩展), 具有温度自动补偿、抗干扰能力强、长期稳定性高等优点.测线布设方式是一条与断层走向斜交(所夹锐角介于30 °~35 °), 另一条与断层走向近正交.
本研究使用1990年1月-2009年12月的日值数据, 部分台站由于仪器故障等原因存在缺值, 缺值统计情况及其他基础资料见表 1.
鲜水河断裂带8个短基线观测场地由北向南分别为侏倭、格篓、虚墟、虾拉沱、沟普、龙灯坝、老乾宁和折多塘(图 1), 测线边长介于72~288m之间, 观测周期基本为每年6次, 极少的时段为每年3、4次, 测量精度优于10 -6, 由于篇幅所限, 场地基础资料可参阅文献[24].便于分析, 统一使用1990年后的数据(龙灯坝短基线场地上, 两条斜交测线的观测结果相互矛盾而无法解释, 故未予采用).短基线观测场地中, 虾拉沱、沟普、老乾宁三个场地同时具有蠕变观测, 可以进行对比分析.
3 数据处理方法近地表的跨断层形变观测受季节、气象等环境干扰明显, 数据一般会呈现年周期变化, 并叠加有高频噪声和随机干扰, 研究显示季节变化、降雨都能够触发或影响断层的蠕滑行为[25-26], 在计算断层形变速率时应尽量消除或减小这些干扰的影响.以往研究中, 主要采用拟合的方法提取数据的趋势变化, 消除年变后计算断层形变速率[24, 27].本文将采用适应性更强的小波变换方法提取长周期趋势形变, 利用趋势变化计算断层运动速率.
3.1 小波变换提取趋势形变小波变换是20世纪发展起来的一种能够进行多分辨率分析的理论, 可以提供信号随时间变化的频谱特征[28].首先对原始数据进行二进小波变换, 根据每层的频率范围和变化特征, 确定年周期变化和趋势变化所在的频段, 将年变和其他高频信息所在频段的小波系数置零, 最后进行重构便得到包含构造活动信息的长周期趋势形变.
针对蠕变仪数据, 选用d b 6二进制小波, 为避免过度分解, 根据数据长度, 计算其最大分解层数为9.蠕变仪日值数据的频率范围介于0~1. 154× 10 -5Hz (对应周期1~∞天), 小波分析结果显示8、9层上的细节部分显示出明显的周期性变化, 其频带为2. 261×10 -8~9. 042×10 -8Hz(128~512天), 因此这部分信号一定包括了周期为365天的年变成分, 9层近似部分则显示为明显的趋势性变化, 频率范围介于0~2.261×10 -8Hz(512~∞天), 说明这部分为信号的长周期趋势项.将所有细节部分的小波系数置零, 剔除相对高频成分, 利用9层近似部分重构便得到趋势形变.
对于短基线数据, 由于观测周期不均匀, 对原始数据进行三次样条插值处理, 插值后每年得到间隔均匀(两月)的6个值.小波基仍选用d b 6小波, 最大分解层数为3, 其三层近似部分的频率范围是0~2.411×10 -8Hz (480~∞天), 因此利用该层的小波系数进行重构可滤除年周期变化, 进而得到趋势变化.上述方法剔除了资料中的年周期变化成分, 以及其他相对高频的细节部分, 这样会损失数据的高频信息, 但对于提取长周期的断层趋势形变来说是合理、有效的.
3.2 计算断层运动参数采用上节中的方法, 提取各条测线的趋势变化, 对其做一阶差分, 得到测线长度的日变化量或两月变化量, 在此基础上根据断层与测线相对位置的几何关系计算断层走滑、张压运动速率(对于缺值的时段进行了线性插值处理).
在小变形情况下, 近似认为断层两盘均为刚体, 并假设地表测值在一定程度上反映了断层整体运动的相对变化.由于断层垂直方向的位移是两条测线伸缩量的高阶分量, 所以这里仅考虑水平方向上的位移.设断层一盘相对另一盘在水平面上的错动矢量为V(Δ S, ΔT), 其中Δ S是平行断层的水平扭错分量, ΔT是垂直断层的张压分量.欲得到断层错动矢量V(Δ S, ΔT), 至少需要两条与断层交角不同的跨断层基线测量, 设两条测线与断层夹角分别为α1、α2 (测线逆时针旋转至断层走向的角度), 测线长度的变化量分别为Δ L 1、Δ L 2, 根据几何关系有如下关系式:
(1) |
由(1)式可得:
(2) |
规定测线长度变化量Δ L 1、Δ L 2以伸长为正; 水平扭错分量Δ S以左行为正; 水平张压分量ΔT以压缩为正.
折多塘场地只有一条斜交测线, 因此没有计算其张压运动分量, 只近似按下式计算水平扭错分量:
(3) |
其中Δ S含义同上, Δ L为该条测线的伸缩量, α为测线与断层走向所夹锐角.
最终, 使用蠕变仪数据得到的Δ S、ΔT为日速率, 短基线得到的Δ S、ΔT为两月(60天)速率.
4 处理结果 4.1 测线长度该变量及其趋势变化特征各条测线长度改变量(即原始观测数据)及采用小波变换提取的趋势变化如图 2所示,曲线上升表示测线伸长,下降表示缩短.可以看出,采用小波变换提取断层形变的趋势变化,对原始数据拟合程度高,无“过拟合”现象,能够较好地反映其时变特性.
从虾拉沱、沟普场地的蠕变和短基线观测来看,主要反映断层走滑运动的斜交测线长度改变量比直交测线大一个数量级,表明鲜水河断裂带以水平走滑为主.
断层蠕变数据中(图 2右列曲线),虾拉沱、沟普斜交测线累计伸缩量较大,趋势显著,相比之下近直交测线的变化不明显;其他测点的伸缩量均较小,趋势也比较波动.
短基线方面(图 2左列曲线),侏倭、虾拉沱、沟普场地的趋势较为单调,累计伸缩量相对蠕变观测大很多.例如,虾拉沱蠕变仪斜交测线约为15mm,短基线则达到44mm之多.相比之下,其他场地的长度改变量不显著,趋势转折较频繁.
在虾拉沱、沟普、老乾宁三个场地上,同时具有蠕变和短基线观测,图 3是两者观测结果的对比(注意图中两种观测的纵轴数值范围不同).由于沟普、老乾宁场地上蠕变和短基线观测的斜交测线布设方式不同,为便于比较,将短基线的斜交测线结果进行了翻转.虾拉沱和沟普场地上,两种观测结果的趋势基本一致,斜交测线均表现为比较单调的上升趋势,近直交测线长度改变量较小,趋势变化不显著.老乾宁场地上,两种观测反映出的断层活动均比较微弱,两者的趋势变化存在较大差异,其中,斜交测线大致在2000年以后呈现相反的趋势;直交测线的变化趋势也有所不同,蠕变仪的直交测线在2001、2002年前后由升转降,而短基线直交测线的趋势变化不明显.
20年来鲜水河断裂带各测点累计走滑量与张压量的空间分布见图 4.
走滑运动方面, 蠕变仪和短基线都显示鲜水河断裂带基本为左行走滑, 且短基线的活动量大于蠕变仪.蠕变仪中累计走滑量最大的是虾拉沱, 达到15.58mm, 其次为沟普的9.46mm, 其他测点的量值都比较小, 其中恰叫和紫马垮两点比较特殊, 显示微弱的右行.短基线中累计走滑量最大的也是虾拉沱, 达到50. 08 mm, 其次为侏倭和沟普, 分别为20.11mm、18.22mm, 其他场地的走滑都很微弱, 其中折多塘观测至今的终态走滑量为右行.两种观测手段均显示, 炉霍、道孚断层的走滑量较大, 而乾宁、折多塘断层活动微弱.
张压运动方面, 两种手段的活动量相差不多.蠕变方面, 虾拉沱、恰叫、沟普显示张性活动, 龙灯坝、老乾宁均为压性, 且龙灯坝较强, 达到3.75mm, 最南端的紫马垮为张性.短基线结果显示, 侏倭、虾拉沱场地为张性活动, 累计拉张量分别为-4.21mm、-4.66mm, 沟普的短基线则为较强的压性, 累计张压量达到5.91mm.以上结果表明, 炉霍、道孚断层表现为张性活动(GP S观测也显示甘孜-炉霍段具有较大的拉张分量[29]), 而乾宁断层为压性活动, 最南端的紫马垮为张性.
对比同一场地的蠕变和短基线结果, 虾拉沱场地上, 两者反映的断层活动性质相同, 均为张性左行走滑, 但活动量上短基线大于蠕变观测; 沟普场地上, 蠕变与短基线均显示左行走滑, 且短基线活动量较大, 但两者显示的张压活动相反, 蠕变为微小的张性, 短基线则为较强的压性.老乾宁场地上两种观测的活动量都很小, 难以说明问题.虽然存在部分矛盾, 但总体上两种手段所反映的断层活动特征较为一致, 均显示出北西强南东弱的活动特点, 主要差别在于活动量的大小.
4.3 断层走滑随时间的变化根据各个测点的累计走滑量, 得到20年来蠕变观测的平均走滑速率介于0.01~0.78mm / a, 短基线的平均走滑速率在0.02~2.46mm / a之间.
下面以观测时间内周边强震的发生时间为节点, 计算出震间时段的平均走滑速率(图 5), 选取的三次强震分别为1997年11月8日玛尼7. 5级(E 1)、2001年11月14日昆仑山8.1级(E 2)、2008年5月12日汶川8.0级(E 3线观测中(图 5 b), 除格篓外, 其他场地也为衰减态势, 说明鲜水河断裂带在炉霍、道孚地震后, 断层面逐步愈合, 向闭锁的方向发展.
图 6给出了各测点走滑速率的年均值变化, 除上述的整体衰减特征外, 沟普以北的场地走滑速率年均值的变化范围较小(短基线表现更明显), 说明炉霍、道孚断层蠕滑比较稳定; 以南波动范围较大, 变化频繁, 表示乾宁、折多塘断层蠕滑状态复杂.
综上所述, 二十余年的断层蠕变与短基线观测显示, 鲜水河断裂带的走滑呈衰减趋势, 分段性显著, 以道孚县为界, 以北的炉霍、道孚断层与以南的乾宁、折多塘断层在活动量、活动性质上差异明显, 炉霍、道孚断层走滑量大且活动方式稳定, 显示张性和左行走滑; 乾宁、折多塘断层活动微弱, 走滑量小且活动状态复杂, 其中乾宁断层为压性和左行走滑, 折多塘断层为微弱的右行走滑, 最南端紫马垮场地为张性和右行, 且以张性为主.
4.4 断层形变的时空演化为更直观地展示断层形变的时空演化过程, 以时间为X轴, 场地的纬度为Y轴, 采用克里金(K r i g i n g)插值, 做出蠕变和短基线观测得到的走滑、张压运动速率等值线图(图 7).
概略的来看, 大部分时间断裂带的走滑都是左行的(图 7(a、b), 色标为红色), 但在某些时段存在逆向走滑(右行走滑), 比较明显的两个时间段为2001年昆仑山地震和2008年汶川地震前后(色标为蓝色), 右行走滑较明显的场地包括恰叫、龙灯坝、老乾宁、折多塘.2001年昆仑山8.1级地震前后, 老乾宁、折多塘场地出现了明显的右行错动, 这在以往的研究中也有所论述[30].汶川地震前, 从蠕变结果来看(图 7 b), 沟普场地从2007年开始左行蠕滑速率开始减小, 之后减小区域向南扩展.短基线结果显示(图 7 a), 沟普及以北的场地左行走滑从2006年中开始持续减弱; 折多塘场地在2006年底至2007年底发生过逆向走滑; 老乾宁场地在2007年中也开始出现右行走滑, 且右行区域有扩大趋势.由此可见, 汶川地震前2年左右, 炉霍、道孚断层表现为左行走滑减弱, 乾宁、折多塘断层在2007年确实出现过右行走滑.
其他时间里, 个别场地也发生过逆向走滑, 但不是很突出, 可靠性及原因有待深入分析.2009年9月, 沟普、龙灯坝的蠕变测点由之前的左行变为右行, 并有扩大趋势; 同期的短基线观测也显示沟普、折多塘场地相继发生了逆向走滑.因此, 逆向走滑区域至少已扩展到乾宁断层和折多塘断层,值得关注.
对比蠕变仪和短基线所反映的断层走滑随时间的变化,虾拉沱场地表现较为一致,观测时段内均为较强的左行走滑;沟普的总体变化趋势有一定的对应性,2009年前基本为左行,之后变为右行,并且都存在多次走滑速率强弱的交替变化;老乾宁测点的一致性较差,两者所反映的旋性有时相反.
张压运动速率的时空变化如图 7(c、d)所示.空间上来看, 蠕变观测显示道孚断层为张性活动, 乾宁断层则以压性为主(图 7 d); 短基线显示的张压活动分为两个阶段(图 7 c), 2003年前为第一个阶段, 炉霍、道孚断层张性活动占主导, 乾宁断层以压性为主, 2003年后为第二个阶段, 其北西段基本都为张性活动.汶川地震后, 恰叫、沟普蠕变仪的张性活动比较醒目.
对比两种观测所反映的断层张压活动, 两者在2003年之前均显示沟普以北为张性, 以南为压性, 但2003年之后, 短基线显示乾宁断层为张性, 蠕变却显示压性.相比之下, 两种手段所反映的走滑运动一致性要优于张压运动, 这可能是因为鲜水河断裂带以左行走滑为主, 走滑运动分量更显著所致.
总结起来, 在走滑运动方面, 虽然鲜水河断裂带以左行走滑为主, 但某些部位存在逆向走滑现象, 比较突出的为乾宁、折多塘断层.汶川地震前2年, 炉霍、道孚断层左行走滑减弱, 乾宁、折多塘断层在2007年出现过右行走滑.张压运动方面, 炉霍、道孚断层以张性为主, 乾宁断层以压性为主.
5 结论(1)不同段落断层活动方式存在差异性.鲜水河断裂带断层活动具有明显的分段性, 以道孚县为界, 以北的炉霍、道孚断层走滑量较大且活动方式稳定, 显示张性和左行走滑; 以南的乾宁、折多塘断层活动微弱, 层为压性和左行走滑, 折多塘断层为微弱的右行走滑.
(2)不同时期断层走滑方式存在交替性.鲜水河断裂带虽以左行走滑错动为主, 但汶川地震前在一些断层段出现过逆向走滑现象.汶川地震前2年, 炉霍、道孚断层左行走滑减弱, 乾宁、折多塘断层在2007年出现过右行走滑, 至2009年底, 逆向走滑区域保持扩展态势.
(3)不同测点间距得到的断层错动速率不一致.随远离断层距离的增加, 平行断层的位移速率逐渐增大, 测点间距18. 7~65. 1 m的蠕变速率为0.01~0.78mm / a; 测点间距72~288m的短基线测量为0.02~2.46mm / a, 点距十几至几十公里的GP S观测为6~11mm / a [8-10]; 地质滑动速率为5~15mm / a [4-7].
6 讨论 6.1 鲜水河断裂带分段活动的影响因素对于鲜水河断裂带分段活动的成因, 有学者认为该断裂带在南东段走向发生偏转, 来自左旋运动的水平侧压增加, 导致断裂南西盘南东向水平滑移受阻, 并转换为垂直运动分量[6]; 此外, 南东段断裂发生分支, 可能对两侧地块的相对运动起到分配作用, 使其走滑速率小于北西段[4].但是, 上述研究的分段位置划在老乾宁附近, 而在本文中, 断层形变观测显示的分段点在道孚县附近, 以北的炉霍、道孚断层与以南的乾宁、折多塘断层在活动量、活动性质上差异明显.
道孚断层与乾宁断层在道孚县附近以左阶雁列衔接, 同时走向发生约15 °的偏转(图 1), 1981年道孚6.9级地震就发生在这个部位, 说明该拐折点是一个较强的障碍体, 上述断层活动差异可能与这种几何结构的改变有关.根据断层滑动资料[31]和震源机制解[32-33], 鲜水河地区主压应力轴方位角介于255 °~295 °之间, 在不考虑其他因素的情况下, 炉霍、道孚断层与应力主轴的夹角要小于乾宁断层与应力主轴的夹角, 故而前两者走滑运动较强, 后者走滑较弱并显示压扭性, 实测的断层形变资料基本证实了这一点.构造物理实验中, 小角度拐折的存在可使其两侧的断层交替成为主导活动构造, 强度较高的非连续部位的变形对断层带的变形有一定的控制作用[ 34], 因此在道孚县附近的拐折构造可能是影响两侧断层活动差异的因素之一.
在巴颜喀拉块体内部, 与龙门山断裂带近平行的龙日坝断裂带, 将巴颜喀拉块体划分为西部的阿坝和东部的龙门山两个次级块体(图 1左下缩略图), 龙日坝断裂带吸收、转换了巴颜喀拉块体的部分S E向运动, 两个次级块体在滑移方向和速率上均存在一定的差异[ 10, 35].龙日坝断裂与鲜水河断裂带大致交汇于炉霍北东[35], 接近于道孚县附近的拐折部位, 阿坝块体与龙门山块体的差异运动, 可能导致它们与川滇块体之间的相对运动有所不同, 反映在断层形变上, 可能以道孚县为界, 两侧断层活动呈现差异.另外, 龙门山地区的GP S观测显示, 在道孚县南西可能存在一个具有右旋运动的地块边界[ 36], 其可靠性及潜在的影响还需进一步的证实.
总之, 本文显示的分段活动可能与断层几何及巴颜喀拉块体内部次级块体的差异运动有关.
6.2 汶川地震前后鲜水河断裂带的逆向走滑汶川地震发生于巴颜喀拉块体的东南边界, 即挤压性质的龙门山逆冲断裂带.由4.4节的描述可知, 汶川地震前, 鲜水河断裂带的炉霍、道孚断层左行走滑减弱, 乾宁、折多塘断层在2007年出现右行走滑, 即以道孚县为分段点, 两侧断层活动的变化模式有所不同.由于鲜水河断裂带的左行走滑是巴颜喀拉与川滇块体南东向的差异运动所致, 因此上述变化说明, 在汶川地震前巴颜喀拉块体东部的南东向运动加剧, 这在跨龙门山断裂的短水准观测中有直接的体现, 五个测点中有三个出现压性增强或由张性转压性的活动[37].
汶川地震后,震区的震后形变GP S连续观测表明,大部分观测点的震后位移方向与同震位移方向相近,表明震后形变在某种程度上是同震形变的延续,其形变源有某种相似性[38].2009年9月,沟普、龙灯坝的蠕变测点为明显的右行走滑,沟普、折多塘的短基线场地也显示断层旋性发生了改变,两种手段同时观测到如此大范围的逆向走滑现象,说明这并非局部行为.更重要的是,逆向走滑区在道孚县附近便不再向北发展,而龙日坝断裂与鲜水河断裂带恰好在此交汇,而且GP S观测发现跨越龙日坝断裂带100km的距离上,在1997-2007年存在约1.8mm / a的水平缩短作用[39],这可能说明由于龙日坝断裂的存在,鲜水河断裂带以道孚县为界,两侧的断层段分属不同的次级块体,因此其断层活动在震前和震后的变化都具有分段性.
综合上述信息,鲜水河断裂带在汶川地震前后出现逆向走滑的可能机理为:汶川地震前巴颜喀拉块体南东向运动增强,龙门山断裂带挤压加剧,由于龙日坝断裂的影响,道孚县以北的炉霍、道孚断层表现为左行走滑减弱,以南的乾宁、折多塘断层表现为逆向走滑;震后,龙门山断裂带积累的弹性应变能释放,表现出与同震形变方向相同的震后松弛形变,震后松弛的范围可能只牵涉龙门山地块,或者以龙门山地块为主,而阿坝地块由于龙日坝断裂的影响,震后形变不如龙门山地块显著,因此,这种变形模式导致地震后道孚县以南的乾宁、折多塘断层发生右行走滑,而以北的炉霍、道孚断层仍为左行.
6.3 断层形变观测与其他观测手段的对比与分析本文的计算结果显示, 蠕变观测得到的鲜水河断裂带走滑速率在0.01~0.78mm / a区间, 短基线测量介于0.02~2.46mm / a, 均小于文献[11]的结果(使用1980-2001年的短基线数据, 本文则为1990年1月-2009年12月), 这种差异应是分析时段的不同造成的, 断层形变随构造活动不断发展变化, 如前所述, 鲜水河断裂带的走滑运动整体上是衰减的, 因此最近20年的平均走滑速率会有所减小.
相比地质滑动速率(5~15 mm / a) [4-7]和GP S观测(6~11mm / a) [8-10], 蠕变和短基线测量显得都很微小.我们选取相同时段的蠕变仪、短基线、GP S观测结果进行比较.布设在道孚的七美(Q I ME)和瓦日(WAR I)跨断层GP S连续观测点, 点距为16.7km (图 1), 与沟普跨断层测量场地的位置接近.2006年2月-2007年11月期间, GP S连续观测显示道孚处平均左旋走滑速率为(6.3±0.3)mm / a [8], 而同期的蠕变与短基线观测均不到2 mm / a.我们将不同手段得到的走滑速率进行对数曲线拟合, 为了更贴合实际, 将断层两侧≥100km的远场走滑速率限定为GP S观测的较大值11mm / a(假设断层两侧变形场反对称分布, 断层一侧平行断层的位移速率则为5.5mm / a), 所得的拟合方程具有统计学意义(图 8 a), 之后由拟合方程得到离开断层的距离与位移速率之间的关系(表 2), 可以看出距离断层越远, 平行断层的位移速率越大, 但变化率逐渐减小, 当离开断层的距离超过30km左右时, 位移速率已达到4.45mm / a, 占总位移速率5.5mm / a的80%以上, 此后位移速率随距离增加缓慢.
我们将不同手段得到的走滑速率除以各自的点距(基线长度), 得到单位距离上(1km)的位移速率, 由于断层的走滑量相对于归分析发现(图 8 b), 断裂带中心的变形强度最大, 远离断层时视剪应变率急剧衰减, 其变形模式更接近剪切不连续的脆性断裂[40].
由上述分析可知, 不同测点间距的断层错动速率不统一, 在相同时段内, 随测点间距的增加, 平行断层的位移速率按对数函数增长, 视剪应变率按幂函数衰减, 超过一定的断层影响区将趋向稳定, 因此测量的空间“尺子”大小不同, 就会导致观测到的形变信息成分不同.近场观测的测线长度从几十米到几百米, 很明显这样的跨距不可能穿越整个变形场, 而只能反映出断裂带的一部分变形.对于远场的GP S观测, 其点距可达几十公里, 可能包含了断层的错动、断层两侧块体的黏弹性或塑性变形, 更能反映块体之间的相对运动.在一个大地震旋回中, 断层形变的发展阶段常常可划分为震前期-同震期-震后期-间震期[40], 断层近、远场的形变速率差异在不同阶段有所不同, 在震后滑动阶段, 断层面近乎开放, 近场滑动速率较高, 与远场速率差异小[41]; 而在震前闭锁阶段, 断层近场的相对运动微弱, 与远场变形速率的差异就会很大.因此, 除了要考虑测量空间尺度的不同, 还要看断层形变处于大地震旋回的哪个阶段, 由目前沟普场地的形变观测情况来看, 近、远场位移速率相差悬殊, 断层的应变积累程度应较高.
地质资料显示的断层滑动速率, 其时间跨度要比现今的形变测量大很多, 是历史上多次震时快速错动与震间缓慢滑动累加后的平均, 而且震时的突然错动应该占主要成分, 代表了断层活动的长期背景, 它与现今断层形变是不同时间尺度的观测, 几十年的形变观测相对于漫长的地质时期过于短暂, 获得的信息量还很有限[42], 它们之间的差异说明断层错动速率在时间上不是常数.
本文的结果只是基于观测资料的初步分析, 还存在一些问题和无法解释的现象, 例如, 同一场地的断层蠕变与短基线观测在某些时段变化不同步, GP S观测并未发现鲜水河断裂带的逆向走滑等, 需要做深入的思考和研究.
致谢感谢两位匿名审稿专家的建设性修改意见, 使本文在原有水平上有很大提高.论文撰写过程中, 与中国地震局地质研究所刘力强研究员、刘培洵副研究员,四川省地震局测绘工程院苏琴高级工程师进行过多次有益交流,在此一并表示感谢!
[1] | Allen C R, Lou Z L, Qian H, et al. Field study of a highly active fault zone:The Xianshuihe fault of Southwestern China. Geological Society of America Bulletin , 1991, 103(9): 1178-1199. DOI:10.1130/0016-7606(1991)103<1178:FSOAHA>2.3.CO;2 |
[2] | 潘懋, 梁海华, 蔡永恩, 等. 中国川西地区鲜水河断裂和则木河断裂几何学、运动学特征及地震活动性对比研究. 中国地震 , 1994, 10(1): 28–37. Pan M, Liang H H, Cai Y E, et al. Comparision between the Xianshuihe fault zone and Zemuhe fault zone, western Sichuan province, China. Earthquake Research in China (in Chinese) (in Chinese) , 1994, 10(1): 28-37. |
[3] | 闻学泽, AllenC R, 罗灼礼, 等. 鲜水河全新世断裂带的分段性、几何特征及其地震构造意义. 地震学报 , 1989, 11(4): 362–372. Wen X Z, Allen C R, Luo Z L, et al. Segmentation, geometric features, and their seismotectonic implications for the Holocene Xianshuihe fault zone. Acta Seismologica Sinica (in Chinese) (in Chinese) , 1989, 11(4): 362-372. |
[4] | 钱洪, 艾伦CR, 罗灼礼, 等. 全新世以来鲜水河断裂的活动特征. 中国地震 , 1988, 4(2): 9–18. Qian H, Allen C R, Luo Z L, et al. The active characteristics of Xianshuihe fault in Holocene. Earthquake Research in China (in Chinese) (in Chinese) , 1988, 4(2): 9-18. |
[5] | 闻学泽, 徐锡伟, 郑荣章, 等. 甘孜-玉树断裂的平均滑动速率与近代大地震破裂. 中国科学(D辑:地球科学) , 2003, 33(Z1): 199–208. Wen X Z, Xu X W, Zheng R Z, et al. The average slip rate and modern strong earthquake rupture of Ganzi-Xianshuihe fault. Science in China (Series D) (in Chinese) (in Chinese) , 2003, 33(Z1): 199-208. |
[6] | 徐锡伟, 闻学泽, 郑荣章, 等. 川滇地区活动块体最新构造变动样式及其动力来源. 中国科学(D辑:地球科学) , 2003, 33(Z1): 151–162. Xu X W, Wen X Z, Zheng R Z, et al. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region. Science in China (Series D) (in Chinese) , 2003, 33(Z1): 151-162. |
[7] | Holt W E, Ni J F, Wallace T C, et al. The active tectonics of the Eastern Himilayan syntaxis and surrounding regions. Journal of Geophysical Research , 1991, 96(B9): 14595-146321. DOI:10.1029/91JB01021 |
[8] | 王敏, 沈正康, 甘卫军, 等. GPS连续监测鲜水河断裂形变场动态演化. 中国科学(D辑:地球科学) , 2008, 38(5): 575–581. Wang M, Shen Z K, Gan W J, et al. Deformation dynamic evolution of Xianshuihe fault by GPS continuous observation. Science in China (Series D) (in Chinese) (in Chinese) , 2008, 38(5): 575-581. |
[9] | Chen Z, Burchfiel B C, Liu Y, et al. Global positioning system measurements from Eastern Tibet and their implications for India/Eurasia intercontinental deformation. J. Geophys. Res. , 2000, 105(B7): 16215-16227. DOI:10.1029/2000JB900092 |
[10] | Shen Z K, Lu J N, Wang M, et al. Contemporary crustal deformation around the Southeast borderland of the Tibetan plateau. J. Geophys. Res. , 2005, 110(B11): B11409. DOI:10.1029/2004JB003421 |
[11] | 程万正, 刁桂苓, 吕弋培, 等. 川滇地块的震源力学机制、运动速率和活动方式. 地震地质 , 2003, 25(1): 71–87. Cheng W Z, Diao G L, Lü Y P, et al. Focal mechanisms, displacement rate and mode of motion of the Sichuan-Yunnan block. Seismology and Geology (in Chinese) (in Chinese) , 2003, 25(1): 71-87. |
[12] | 龙德雄, 邓天岗. 鲜水河断裂带断层运动分段特征的初步研究. 地震研究 , 1986, 9(5): 591–600. Long D X, Deng T G. A preliminary study on segment feature of fault motion along Xianshuihe Fault belt. Journal of Seismological Research (in Chinese) (in Chinese) , 1986, 9(5): 591-600. |
[13] | 闻学泽, 白兰香. 鲜水河活动断裂带形变组合与运动特征的研究. 中国地震 , 1985, 1(4): 53–59. Wen X Z, Bai L X. Study on the deformation composition and the motion feature of the Xianshuihe active fault zone. Earthquake Research in China (in Chinese) (in Chinese) , 1985, 1(4): 53-59. |
[14] | 张家涛, 姚光亮. 鲜水河断裂带的分段活动特征. 地壳形变与地震 , 1990, 10(3): 52–60. Zhang J T, Yao G L. The characteristic of segmented activity of Xianshuihe fault zone. Crustal Deformation and Earthquake (in Chinese) (in Chinese) , 1990, 10(3): 52-60. |
[15] | 杨主恩, 郭芳, 李铁明, 等. 鲜水河断裂西北段的断层泥特征及其地震地质意义. 地震地质 , 1999, 21(1): 21–28. Yang Z E, Guo F, Li T M, et al. The characteristics of fault gouges from the Xianshuihe fault zone and its geological implications. Seismology and Geology (in Chinese) (in Chinese) , 1999, 21(1): 21-28. |
[16] | 周硕愚, 施顺英, 宋永厚, 等. 鲜水河断裂带现今断层运动短周期事件的发现与初步研究. 地壳形变与地震 , 1996, 16(4): 1–8. Zhou S Y, Shi S Y, Song Y H, et al. Discovery and preliminary research for short-period phenomena of present-day fault movement on Xianshuihe fault zone. Crustal Deformation and Earthquake (in Chinese) (in Chinese) , 1996, 16(4): 1-8. |
[17] | 吕弋培, 王兰, 苏琴. 断层蠕变的分维变化及其与地震之关系. 地壳形变与地震 , 1999, 19(4): 62–67. Lü Y P, Wang L, Su Q. Changes of fractal dimension of fault creep and their relations to earthquake. Crustal Deformation and Earthquake (in Chinese) (in Chinese) , 1999, 19(4): 62-67. |
[18] | 张国民, 邵志刚, 刘桂萍, 等. 汶川8.0级地震构造动力成因分析. 地震 , 2009, 29(1): 164–173. Zhang G M, Shao Z G, Liu G P, et al. Analysis on the geodynamic cause of the Wenchuan 8.0 earthquake.. Earthquake (in Chinese) , 2009, 29(1): 164-173. |
[19] | 徐锡伟, 闻学泽, 叶建青, 等. 汶川MS8.0地震地表破裂带及其发震构造. 地震地质 , 2008, 30(3): 597–629. Xu X W, Wen X Z, Ye J Q, et al. The MS8.0 Wenchuan earthquake surface ruptures and its seismogenic structure. Seismology and Geology (in Chinese) , 2008, 30(3): 597-629. |
[20] | Nason R D, Phillippsborn F R, Yamashita P A. Catalog of creep-meter measurements in central California from 1968 to 1972. U S Geol. Surv. Open File Rep , 1974. |
[21] | Goulty N R, Davis P M, Gilman R, et al. Meteorological noise in wire strainmeter data from parfield, California. Bull. Seismol. Soc. Am. , 1979, 69(5): 1983-1988. |
[22] | Lee J C, Jeng F S, Chu H T, et al. A rod-type creepmeter for measurement of displacement in active fault zone. Earth Planets Space , 2000, 52: 321-328. DOI:10.1186/BF03351643 |
[23] | King C Y, Nason R D, Burford R O. Coseismic steps recorded on creep meters along the San Andreas Fault. J. Geophys. Res. , 1977, 82(11): 1655-1662. DOI:10.1029/JB082i011p01655 |
[24] | 杨永林, 苏琴. 鲜水河断裂带现今活动特征研究. 大地测量与地球动力学 , 2007, 27(6): 22–27. Yang Y L, Su Q. On present activity characteristics of Xianshuihe fault Zone. Journal of Geodesy and Geodynamics (in Chinese) (in Chinese) , 2007, 27(6): 22-27. |
[25] | Roeloffs E. Creep rate changes at parkfield, California 1996-1999:Seasonal, precipitation induced, and tectonic. J. Geophys. Res. , 2001, 106(B8): 16525-16547. DOI:10.1029/2001JB000352 |
[26] | Lee J C, Angelier J, Chu H T, et al. Active fault creep variations at Chihshang, Taiwan, revealed by creep meter, monitoring, 1998-2001. J. Geophys. Res. , 2003, 108(B11): 2528. DOI:10.1029/2003JB002394 |
[27] | 吕弋培, 李铁明, 廖华. 鲜水河断裂带跨断层形变测量及其地震学意义. 地震地质 , 1997, 19(4): 333–340. Lü Y P, Li T M, Liao H. Across-fault deformation measurement along the Xianshuihe fault zone and its seismological implications. Seismology and Geology (in Chinese) (in Chinese) , 1997, 19(4): 333-340. |
[28] | Mallat S. A theory for multiresolution signal decomposition the wavelet representation. IEEE Transactions on PAMI , 1989, 11(7): 674-693. DOI:10.1109/34.192463 |
[29] | 王阎昭, 王恩宁, 沈正康, 等. 基于GPS资料约束反演川滇地区主要断裂现今活动速率. 中国科学(D辑:地球科学) , 2008, 38(5): 582–597. Wang Y Z, Wang E N, Shen Z K, et al. Constrained inversion on major fault slip rate of Sichuan-Yunnan region based on GPS data. Science in China (Series D) (in Chinese) (in Chinese) , 2008, 38(5): 582-597. |
[30] | 苏琴, 杨永林, 张永久, 等. 昆仑山口西8.1级地震对鲜水河断裂带断层活动的影响. 地震研究 , 2009, 32(2): 169–174. Su Q, Yang Y L, Zhang Y J, et al. Influence of the Kunlun Mountain MS8.1 earthquake to the Xianshuihe fault activity. Journal of Seismological Research (in Chinese) , 2009, 32(2): 169-174. |
[31] | 谢富仁, 李宏. 利用断层滑动资料确定鲜水河断裂带现代构造应力的方向和大小. 地震学报 , 1995, 17(2): 164–171. Xie F R, Li H. Magnitude and direction of modern tectonic stress field of Xianshuihe fault zone using fault slip data. Acta Seismologica Sinica (in Chinese) (in Chinese) , 1995, 17(2): 164-171. |
[32] | 孙建中, 施顺英, 周硕愚, 等. 利用地震矩张量反演鲜水河断裂带现今运动学特征. 地壳形变与地震 , 1994, 14(4): 9–15. Sun J Z, Shi S Y, Zhou S Y, et al. Inversion of the present kinematic characteristics of Xianshuihe fault zone from seismic moment tensor. Crustal Deformation and Earthquake (in Chinese) (in Chinese) , 1994, 14(4): 9-15. |
[33] | 刘春, 邵辉成, 石军, 等. 利用地震矩张量初步分析断裂带的运动学特征. 地震地质 , 2005, 27(3): 429–436. Liu C, Shao H C, Shi J, et al. Preliminary analysis of kinematic features of a fault zone by seismic moment tensors. Seismology and Geology (in Chinese) (in Chinese) , 2005, 27(3): 429-436. |
[34] | 马瑾, 马胜利, 刘力强, 等. 断层几何结构与物理场的演化及失稳特征. 地震学报 , 1996, 18(2): 200–207. Ma J, Ma S L, Liu L Q, et al. Geometrical textures of faults, evolution of physical field and instability characteristics. Acta Seismologica Sinica (in Chinese) (in Chinese) , 1996, 18(2): 200-207. |
[35] | 徐锡伟, 闻学泽, 陈桂华, 等. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义. 中国科学(D辑:地球科学) , 2008, 38(5): 529–542. Xu X W, Wen X Z, Chen G H, et al. Discovery of the Longriba faults, eastern part of the Bayankela tectonic block and its geodynamic implications. Science in China (Series D) (in Chinese) , 2008, 38(5): 529-542. |
[36] | Zhang K L, Wei D P. Implications from the kinematic pattern of the Longmenshan Region. Tectonophysics , 2011, 504(1-4): 57-64. DOI:10.1016/j.tecto.2011.01.019 |
[37] | 焦青, 杨选辉, 许丽卿, 等. 汶川8.0级地震前后龙门山断裂活动特征浅析. 大地测量与地球动力学 , 2008, 28(4): 7–11. Jiao Q, Yang X H, Xu L Q, et al. Preliminary study on motion characteristics of Longmenshan fault before and after MS8.0 Wenchuan earthquake. Journal of Geodesy and Geodynamics (in Chinese) , 2008, 28(4): 7-11. |
[38] | 王敏, 沈正康, 万永革等. 2008年汶川MW7.9地震GPS震后形变场监测与龙门山断裂带震后滑移及其邻域流变学响应研究.//中国地球物理2010--中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集. 2010. |
[39] | 杜方, 闻学泽, 张培震, 等. 2008年汶川8.0级地震前横跨龙门山断裂带的震间形变. 地球物理学报 , 2009, 52(11): 2729–2738. Du F, Wen X Z, Zhang P Z, et al. Interseismic deformation across the Longmenshan fault zone before the 2008M8.0 Wenchuan earthquake. Chinese J. Geophys. (in Chinese) , 2009, 52(11): 2729-2738. |
[40] | Scholz C H. The Mechanics of Earthquakes and Faulting. (2nd ed). New York: Canbridge University Press, 2002 . |
[41] | 杜方, 闻学泽, 张培震. 鲜水河断裂带炉霍段的震后滑动与形变. 地球物理学报 , 2010, 53(11): 2355–2366. Du F, Wen X Z, Zhang P Z. Post-seismic slip and deformation on the Luhuo segment of the Xianshuihe fault zone. Chinese J. Geophys. (in Chinese) , 2010, 53(11): 2355-2366. |
[42] | 车兆宏, 张艳梅. 南北地震带中南段断层现今活动. 地震 , 2001, 21(3): 31–38. Che Z H, Zhang Y M. The recent fault activity along the middle and southern segments of South-North Seismic Belt. Earthquake (in Chinese) (in Chinese) , 2001, 21(3): 31-38. |