南北地震带是纵穿中国大陆中部的一条构造带，从北端鄂尔多斯地块西缘起，经过秦岭、龙门山，沿鲜水河-安宁河-小江断裂带进入缅甸境内，是青藏高原与鄂尔多斯高原、四川盆地及华南地块的分隔带^[1]，历史上曾孕育和发生过多次大地震。2013-07-22甘肃岷县-漳县交界发生的M6.6地震，位于南北地震带北段，介于西秦岭北缘断裂带和东昆仑断裂带之间，是继2008年汶川特大地震及2010年玉树地震之后，发生在南北地震带上的又一次大地震^[2]。本文通过分布在南北地震带上观测质量较好的10个洞体应变台2013年观测资料，分析地震前后的变化特征，提取震前异常信息，并结合其他学者研究结果^[3-4]，在南北地震带这一构造运动背景上，初步探讨这次地震的演化过程。

1 数据处理

本研究采用分布在南北地震带上的10个洞体应变台观测资料(图 1), 各台距震中的距离在表 1中列出。这些台站的观测仪器为应变伸缩仪，观测值为地表的相对伸缩量，伸长为正值，缩短为负值。时间段为2013年，数据采样率为每分钟一个样点。首先进行数据预处理，包括个别缺值内插，以及去趋势性变化。分别对每个台站的北南向和东西向进行多尺度小波分解，得到细节分量的第14层分量。小波分解的目的是为了尽可能去除高频干扰和长周期潮汐变化对结果的影响。以每个分量的细节分量的第14层结果作为输入数据计算，每天得到一个形变分量的范数。

范数是在赋范线性空间中定义，且满足一定条件(非负性、齐次性、三角不等式和相容性)的函数。常见范数有1-范数、2-范数和∞-范数。本研究采用2-范数：

2-范数也叫Euclid范数，得到的是向量矩阵的最大奇异值。洞体应变观测值为地表的相对伸缩量。从每天的观测值得到一个相对伸缩量的范数，在时序上表述了地表伸缩量的最大变化率，是对当地构造运动状态的一种描述。

2 结果与讨论

每天相对伸缩量的范数时序表达如图 2。为了方便讨论，给出一条参考线，范数为10 000，对范数大于该值的变化进行重点讨论。

从图 2可以看到，应变范数在背景变化的基础上在地震前各台都有所增大，其中变幅较大的台站依次为昭通台、昆明台、永胜台和宕昌台。有学者研究了川滇地区应变与地震的关系，认为应变持续加速变化超过相应阈值后可能有M5.8以上地震发生，这与本研究结果是相吻合的^[5]。本研究中，各台呈现出一个共同的特征, 即除昭通台以外，其他各台东西分量的范数变化远大于南北分量。从异常变化出现的时间可见，最早出现异常的是昭通台，异常起始于2013-06初，之后依次为永胜台、昆明台和宕昌台，起始时间分别为6月中旬和7月初。地震之后，最大变化幅度出现的时间显示出从震中开始由北及南逐渐滞后的现象，而处在震中北方一侧的兰州台、合作台地震引起的幅度变化则极其微弱。综合各台变化幅度及异常起始时间以及各台的地理分布，似乎显示出以昭通地区为中心向外辐射的基本态势。如果沿南北地震带考察应变范数的变化特征可以看到，在地震之前大致呈现由南及北的异常变化过程，而地震之后则显示出由北及南的推进次序。一些文献从不同时间尺度研究了南北地震带地震发生的时空关系，认为存在南北两端的跳跃现象^[6-8]，一定程度上揭示了南北两端地震的相互关系，这些观点与本文结果相一致。

文献[9]利用292个GPS站点多年的观测结果，结合本地区大地震的震源机制解、主干活动断裂分布及其他研究结果，将青藏高原东南缘划分为11个子块体，并给出除去整体运动后的运动特征图，见图 3。

图 3中，虚线表示子板块边界，粗箭头方向示意了除去整体运动后块体运动的方向^[2]。另外，给出了每个台站形变范数时序图，如小框图所示，左图为北南向，右图为东西向。为了便于比较各台站的变化幅值，每个台站两个分量的幅值分别以昭通台相应分量的最大值为参考进行归一化，岷县-漳县地震的时序位置用小箭头标出。

由图 3可以看出，这一区域的运动特征大体可以概括为3种情形，第1种是藏东子块体、香格里拉子块体、雅江子块体、滇中子块体、景谷子块体、保山子块体和西盟子块体的整体顺时针旋转运动；第2种是阿坝子块体的北东向推进；第3种是南北地震带东侧自南向北的运动。小庙台、昭通台、永胜台和昆明台在距震中760~1 100 km的范围，肃南台距离震中约650 km。从图 2可以看到，在地震前存在两种现象：一是小庙台、昭通台、昆明台、永胜台形变范数远大于其他台站，甚至大于距震中最近的宕昌台，特别是北南向(可见本区域内最大的形变异常不是在震中区)；二是异常出现的时序依次为小庙台、永胜台、昭通台、昆明台和宕昌台，是一个由南及北的过程。结合文献[9]的结果推测，2013-07-22岷县-漳县M6.6地震的发生与南北地震带东侧自南向北的运动过程密切相关，这一过程的动力源从本文掌握的资料看可能来自昭通-小庙这一区域，这一异常过程向北传递至甘东南地区触发了本次地震。另外可以看出，对于2013-04芦山M7.0地震，在小庙台和同仁台有较明显的震前异常反映，其余多台未见明显变化。由上述现象得到两点基本认识：1)这两次地震具有不同的动力来源；2)岷县-漳县地震之前存在一个大体自南向北的动力学过程，最大幅度的形变并非地震发生的充要条件。来自文献[9]的GPS观测结果(图 3)较好地解释了各个台形变范数变化幅度大小的顺次关系及异常出现时间的时序关系。换句话说，也支持了本研究的上述观点。

3 结语

1) 本研究通过小波分解和范数计算方法较好地获得了2013-07-22岷县-漳县M6.6地震前的异常信息。最早异常出现在震前50 d，最近台站的异常出现在震前20 d。

2) 2013-07-22岷县-漳县M6.6地震的发生与南北地震带东侧自南向北的运动过程密切相关，这一过程的动力源从本文掌握的资料看应该来自小庙-昭通这一区域；而震后的形变传递过程是一个反向过程。

3) 洞体应变的范数计算对于形变异常的提取是一种简单而有效的处理方法

4) 由于目前形变观测还未达到成场化观测，限于观测台站和资料及相应震例的数量，本研究未涉及对异常的判定指标及与地震位置和大小的定量关系的研究，这是今后进一步工作的一个方向。