云南地处印度板块与欧亚板块中国大陆碰撞带东缘,地震活动频度高、震级大、分布广(皇甫岗等,2010),是研究地壳运动与应变积累及其与强震关系的热点区域(魏文薪等,2012; 周海涛等,2017)。1982年初以来,先后在川滇菱形块体西南、东南边界带上建立了多个跨断层形变监测场地,用于监测大地形变,分析强震孕育过程中断层形变的发展演变过程。随着观测数据的不断积累,对跨断层数据的处理与分析研究也逐步深入。张超等(1981,1987)详细推导了利用跨断层基线、水准数据计算断层运动走滑、张压、垂直升降三分量的计算公式,在此基础上,多位学者分析了断层运动时间演化特征及其与地震活动的关系。薄万举等(1997)利用信息合成的方法提取了断层的趋势及短期变化信息; 李腊月等(2013,2014)利用主成分分析的方法综合分析了首都圈跨断层形变资料; 焦青和周俊萍(1998)利用跨断层资料分析了断层位移速率的动态演化特征与强震的关系; 李瑞莎等(2009)论证了川滇地区断层活动与地震的关系。另外,对预测指标的构建和应用研究也越来越多元和系统化,张晶等(2011)采用断层三维立体模型计算了断层三维运动参数,进而提出了断层活动协调比的概念; 李瑞莎等(2013)根据断层协调比的概念,讨论了起点变化对协调比计算结果的影响,提出了移动窗计算断层活动协调比的方法以及离散度评定指标; 张希等(2012)基于灰色关联度指标,研究了鲜水河断裂活动特性与大震的关系。

以往的研究大都基于单个指标进行单一尺度的区域地震危险性分析,本文对前人提出的多种跨断层异常分析方法进行了整合,并重点对反应断层应变状态变化的断层活动协调比指标进行了改进研究,在此基础上梳理了云南地区最新的跨断层数据,首先利用断层活动异常频次(刘强等,2010)分析整个云南地区的强震背景; 并重点计算了滇东南区域的断层运动年速率变化、断层运动主成分信息合成、断层运动协调比指标,并结合震例对该区域的强震危险性进行了综合分析。

近年来,云南地区多个跨断层场地先后出现了大幅度的拉张、压缩变化,断层活动加剧,异常变化显著。针对多个场地异常发展变化的群体性、同步性特征与强震孕育的关系,刘强等(2010)提出了断层形变异常参数,其计算公式如下:

R_i,j=((m_i,j-m_i-1,j)/(t_i,j-t_i-1,j))/(1/n∑ⁿ_i=1〖JB<2|〗(m_i,j-m_i-1,j)/(t_i,j-t_i-1,j)〖JB>2|〗+(1/n∑ⁿ_i=1(〖JB<2|〗(m_i,j-m_i-1,j)/(t_i,j-t_i-1,j)〖JB>2|〗-1/n∑ⁿ_i=1〖JB<2|〗(m_i,j-m_i-1,j)/(t_i,j-t_i-1,j)〖JB>2|〗)²)^1/2)(1)

式中:n为所用资料的总年数减去1; i表示年(i=1,2,…,n); j表示月(j=1,2,…,12); m_{i-1, j}和m_{i, j}分别为第i-1年和第i年的第j个月的观测值; t_{i-1, j}和t_{i, j}为第i-1年和第i年的第j月观测时间值; R_{i, j}为第i年第j个月的异常参数值。

笔者选取了持续时间较长、观测较为完整的10个场地的水准观测数据(表1),计算了各测项的异常参数值,并以 〖JB<1|〗R_{i, j}〖JB>1|〗>1 作为异常阈值统计了断层活动异常频次变化(图1)。从图1可以看出:自1982年观测以来,云南及周边地区有4组6.5级以上强震发生,均发生在异常频次较高的时段后(5组高频异常后2年左右有4组对应了该区域的成组强震),目前异常频次较高且持续时间较长,整个云南地区的强震危险性背景较强。

表1 云南地区跨断层场地概况(据刘强等,2010)
Tab.1 The overview of cross fault sites in Yunnan(based on Liu et al,2010)

图1 1982年2月—2018年8月云南地区断层活动异常频次与MS≥6.5强震对应关系
Fig.1 The relationship between the anomaly frequency of fault activity and MS≥6.5 strong earthquakes from Feb.,1982 to Aug.,2018

如图2所示,滇东南地区存在多项显著的跨断层观测异常,例如石屏基线2-1、2-3自2010年起呈现持续拉张变化,目前转平且伴有短期压缩变化; 石屏水准2-1、2-3自2013年起出现持续正断活动,目前速率变缓且伴有短期逆断活动变化; 建水基线2-1在大幅拉张后出现回复趋势; 宜良基线1-2出现趋势压性异常并在低位震荡变化; 峨山水准1-2在正断背景下出现显著的逆断变化。

为了消除年变的影响,进一步分析各场地速率变化异常特征,计算上述异常项的年滑动速率(每个月一个值)及年平均滑动速率(每年一个值),如图3所示。从图中可见,建水基线2-1拉张年平均速率,以及石屏基线2-3压缩年平均速率持续增大,表明楚雄—建水断裂应变积累不断增强,其他场地变化不明显。

图2 滇东南地区显著跨断层形变异常
Fig.2 The significant fault deformation anomalies in Southeast Yunnan

图3 滇东南地区显著跨断层异常场地运动速率
Fig.3 The fault movement rates of significant across fault sites in Southeast Yunnan

跨断层观测数据既包含断层活动的构造信息,又包含季节、气象和地下水资料等非构造的干扰信息,而这些干扰信息往往具有一定的相关性。主成分分析方法(PCA)能将多个具有相关性的指标(对应于多个观测台站数据),分解成几个相互独立的综合指标(信息不重叠,每个主成分可以单独分析),进而通过贡献率加权重组,得到能够反映区域构造运动特征的时间序列,在此基础上更客观地分析其变化和地震的关系。

在使用PCA对断层运动数据进行信息合成之前,首先利用同台的水准、基线数据及其与断层的夹角信息计算断层运动三分量,即走滑ΔS、张压ΔP和垂直升降量ΔH(张超,1981)。选取图4矩形框中的楚雄、峨山、建水、石屏、通海、寻甸、宜良7个台站,分别计算其三维运动分量。图5结果显示,整个研究区域的运动具有较好的一致性,以右旋张性运动为主,正断运动显著。

分别将垂直升降、走滑、张压分量代入PCA模型,计算了断层运动沿3个正交方向上的PCA运动合成,其中垂直升降量的PCA信息合成结果与云南地区的强震对应结果较好。如图6所示,

图4 滇东南地区断层运动PCA合成台站及震例选取范围
Fig.4 The station and range of earthquake selection used for PCA in Southeast Yunnan

图5 滇东南地区7个场地断层运动分解结果
Fig.5 The result of fault movement decomposition of seven sites in Southeast Yunnan

2007年6月3日云南宁洱6.4级地震、2014年10月7日云南景谷6.6级地震前,综合指标出现明显的台阶,具体表现为“正断活动加剧—趋于平稳—发震”的形态演化特征,目前的计算结果显示该类变化不突出,该区域断层正断活动加剧异常正在演化过程中。

图6 滇东南地区垂向运动PCA合成
Fig.6 The PCA synthesis of vertical movement in Southeastern Yunnan

云南地区的断层活动异常频次计算结果表明,频次高值与本地及周边强震(M≥6.5)具有较好的对应关系,近年来该区域断层运动异常显著,具有较强的强震发震背景。2018年8月13,14日通海5.0级地震发生后,通过对滇东南各个场地的跨断层观测进行三维运动分解、PCA信息合成,笔者发现其计算结果与原始曲线变化具有较好一致性,均表明该区域自2012年以来,多个场地均出现显著的拉张变化,区域正断运动背景显著; 在此背景下,各异常台站的运动速率的不同,显示出同一构造区内不同断层段的运动差异,这种运动差异有利于应变积累,可能对地震具有一定的触发作用。此外,改进断层运动协调比的计算结果显示,石屏、宜良台出现了显著的应变状态变化。

2018年8月13,14日通海5.0级地震前,滇东南地区存在多项显著的跨断层观测趋势异常,短期变化不显著,尽管震前距离震中较近的石屏台、宜良台出现了显著的断层运动协调比异常变化,但结合以往震例分析结果表明,各项异常仍在演化过程中,且尚未达到以往6级以上地震前的异常幅度和持续时间。综上,滇东南地区的断层运动异常幅度较大、持续时间较长,2次通海5.0级地震震级较小,不能完全对应给该区域的多项异常,滇东南地区的应变积累还在持续,区域地震危险性仍在增强。