跨断层监测作为地震监测的重要手段,在中国有近半个世纪的观测历史,跨断层资料为断层活动与地震关系的研究奠定了良好的基础。目前,对于跨断层资料的分析方法有很多,有的注重断层活动状态、不同场地活动关联性的分析,如张希等(2012,2018,2020)在跨断层资料的分析研究方面做了大量工作,提出了灰色关联度、特征强度等一系列方法用于异常的识别,且部分方法对断层活动状态有较好的描述,从一定程度上能够解释断层活动与地震之间的关系(田晓等,2021)。但受场地观测流程复杂、复测困难等因素影响,云南地区的跨断层观测多以单断层的监测为主,且多数场地分布在主干断裂的某一分支断裂,难以完全跨过主干断裂,导致分析断层活动机理与地震之间的关系比较困难。例如楚雄场地观测到的是长期左旋走滑活动,与该区域断层右旋走滑的机理不一致,不能完全代表主干断裂的活动特征。另外一些学者的方法主要是基于数学方法对观测时序进行异常分离与识别(薄万举等,2001; 贾鹏等,2016; 张晶等,2020; 李瑞莎等,2020,2016),但异常的识别依赖于正常背景状态的明确,这些方法能够对异常进行识别,但无法从机理上给予更多的解释。

协调比的计算通过描述断层自身各分量的协调性来表征断层周边应变积累程度,避免了断层活动与板块运动之间机理解释的矛盾性,同时也能够反映出断层异常状态与区域构造活动之间的相关性,特别适合云南地区跨断层观测场地的现状。在实际使用断层协调比描述断层协调性时,存在计算起点不同、结果出现显著差异的问题,导致其在会商应用时一直存在困难。李瑞莎等(2013)提出了改进办法,在计算断层滑动量时不再采用统一起点,而采用相邻两期的观测数据差分来计算断层的活动速率,该方法较为有效地避免了因起点不同而引起的计算结果的差异问题。从计算结果的物理意义看,李瑞莎等(2013)认为该协调比结果是断层的“瞬态协调比”,可称其为短期协调比,经过检验,断层协调比结果虽然不再受计算起点影响,但是不能较好地反映断层协调比异常的累积状态。李瑞莎等(2013)使用滑动窗计算协调比的方差来反映断层一段时间内协调比的离散程度,从一定程度上兼顾了协调比异常状态的持续时间的因素。该方法对原有协调比进行了较好的改进,利用方差虽然在反映偏离离散程度方面能够在一定程度上表达协调比偏离背景的异常频度,但在反映异常偏离累积程度及幅度、异常出现同向偏离等方面仍存在欠缺。

因此,本文提出了扰动协调比的计算方法,对“瞬态协调比”的合成算法进行了改进,在继承该方法不受起点影响的特性的同时,进一步体现了协调比持续时间和幅度等信息的合成,在表达断层异常状态的持续时间方面更具优势,以期尽可能充分显示协调比异常程度,将其并应用于震情会商中,以云南地区11处跨断层资料为例,针对川滇地区M≥6.5强震合成得到了多个预测指标。

张晶等(2011)提出了协调比方法,该方法在区分正常动态和异常变化上有显著的优势,能够较为清晰地判别断层的异常活动情况。该方法的理论基础是,无论断层处于纯走滑还是走滑兼倾滑的无障碍蠕动状态,断层3个运动分量的比值处于相对稳定的恒定值; 而当断层存在显著应变积累时,断层3个运动分量出现不协调变化,导致其比值偏离背景值,从而有效地区分正常状态和异常。许多学者针对断层协调比与地震的关系开展了大量研究工作(徐东卓等,2021; 李智蓉等,2020; 李瑞莎等,2012; 李文静,2014; 张立成,薄万举,2021; 周海涛等,2018; 苑争一等,2018)。

计算断层协调比时,首先根据断层观测值计算断层活动的三维分量,断层活动性各分量和断层及测线之间的关系如下(薄万举等,1998):设测线与断层夹角为α(α定义为测线顺时针转动到与断层线重合时所转过的角度),测线伸长量为ΔL(缩短符号为负),上盘逆时针相对下盘的扭错量为a、水平拉张量为b、垂直上升量为c、水准变化量为Δh(a>0表示左旋、b>0表示张性、c>0表示上盘上升,由上盘到下盘为往测方向),则在二条不同的测线上,同时施测水准和基线,可唯一确定断层运动参数:

定义f₁、f₂、f₃为断层活动协调比,用于描述断层3个分量两两之间的协调性,并认为当断层处于无障碍的蠕滑状态时,这3个分量为恒定值,可表示为:

从机理上来看,断层“瞬态协调比”异常存在两种可能的原因:一是按照断层协调比的定义,断层上出现显著的应变积累,导致断层活动状态发生改变,从而出现协调比异常,直到地震发生后应力得到释放为止,这类异常对应变积累区附近场地影响可能更为突出,持续时间贯穿应变积累到发震的整个过程,在此过程中计算得到的瞬态协调比将可能持续反复出现异常; 二是协调比基本原理是断层分量之间的比值,受周边其他强震影响,当周边其他区域强震发生前发震断层由稳态进入亚失稳状态时,可能导致区域应力出现调整,引进应力状态发生改变,而跨断层场地位于敏感带,容易受到这种应力状态调整的影响,多个场地出现准同步的异常变化(张希等,2020),从而导致断层协调比出现异常,可认为是强震亚失稳状态前出现的区域应变场扰动导致的异常。无论是哪种原因导致的协调比异常,瞬态协调比都可能反复出现,但可能在出现时间和幅度上存在差异。因此,为了更好地表示瞬态协调比反复出现的频次和幅度,本文提出“扰动协调比”的概念,即对于任意协调比参量时间序列,将瞬态协调比归算为一定时间窗的协调比累积值。假设现有某一跨断层场地通过相邻期观测数据计算得到的协调比分量时序为f_i(i=1,…,n),对于时序中任意时间点的值f_k,假设时间窗w内,f_k前有m个值,计算其扰动协调比为:

该归化过程总体来说使用f分量一定窗长(w)内的协调比的绝对值进行累加, 为遗忘因子,距离当前值时间越远的量对当前值贡献越小,越近则越大,从而形成了考虑影响衰减因素的协调比累积结果计算方法,最终得到的参量反映了断层协调比异常状态在特定时间窗内的异常程度,即为扰动协调比。

利用本文方法计算得到峨山场地年窗长的应力场扰动协调比结果(图1)。从图1可以看出,峨山场地对川滇菱形块体边界带,甚至龙门山断裂带附近的强震均有显著响应,且在强震前1～3年就已经出现,持续时间较长。因此,扰动协调比在断层状态描述方面有着显著优势,在窗长恒定时,其结果具备唯一性,与协调比计算的起点无关(图1c、d中为不同起点)。从物理意义上来说,协调比的扰动现象反映的是强震前区域应力场出现的扰动现象,这种扰动在前兆观测方面可能表现为大范围的前兆异常项数增多现象,而在断层上的观测则表现为断层活动状态偏离常态。

图1 峨山场地年窗长区域应力场扰动协调比时序图(绿色为准确预测地震,蓝色为漏报地震,下同)
Fig.1 Time-series diagram of coordination ratio disturbance in annual window length area of Eshan site (Green stands for accurately predicted earthquakes; blue stands for missed earthquakes,they have the same meaning in Fig.3 and Fig.4)

本文在瞬态协调比合成算法的基础上提出扰动协调比计算方法,以云南地区11处跨断层资料为例,针对1981—2021年川滇地区6.5级以上强震样本开展了震例总结和综合预测指标的提取,对预测指标进行效能评估,得到以下结论:

(1)扰动协调比计算方法在跨断层场地异常活动状态的提取中有一定优势,与周边强震存在较好的相关性,说明震前由于区域应力场的变化,导致断层的活动出现背离背景的现象,利用该现象开展强震的预测具有一定的意义,能够较好地识别震前断层出现的协调比异常状态,对于强震可能发生时间可提供中期尺度的判定依据。

(2)扰动协调比计算方法在异常信息的累积方面做出了进一步的优化,从而有效体现了一定时间尺度内断层的异常活动状态的持续效应,提升了异常信息的信噪比,有利于后续的震例总结,在震前异常的提取方面有较好优势。

但从目前计算的结果看,扰动协调比仍然存在震后效应和震前异常无法完全区分的问题。如2005年以后各个合成指标均出现不同程度的异常,特别在全部扰动协调比分量合成结果中更加明显,而该异常从大区域的地震活动看,可能与2004年印尼9.3级地震震后影响有关。因此,在扰动协调比信息合成过程中,地震样本的选取尤为重要,从扰动协调比的基本原理来说,应当和一些震级较大的地震相关性更好,因为此类地震造成的区域应变场的扰动范围可能更大,群体异常出现的几率更高。

本文研究中地震的选取不一定合理,但在一定程度上说明了在强震前由于区域应力场的扰动,存在导致跨断层场地断层协调比背离正常背景的现象,这种现象持续时间1～3 a不等,通过扰动协调比的计算方法能够较好地提取到这些异常现象。通过进一步震例的总结和筛选,可提取出符合本区域的扰动协调比综合预测指标,辅助会商决策。