基于应变释放加速的未来强震危险性分析与预测.docxVIP

基于应变释放加速的未来强震危险性分析与预测 主震前不同时段、不同空间与应变释放的震例分析 地震专家长期以来一直注意到，在一些强震发生之前，可以观察到地震活动和地震活动所代表的变形速度的“加速”现象。其中，随着地震活动的变形释放和加快，还有许多方法可以用菲比安顿率进行定量描述。它还开发了一种基于“加速”进程的计算方法，该过程被称为“破坏时间分析”。该方法通过强震前的应变释放加速（Accelerating Moment Release,AMR）现象进行强震发生时间的预测，相应的分析方法也已在一些国际震例的回溯性检验中显示了较好的预报有效性。国内也有一些学者对“破裂时间分析”方法进行改进与应用研究。其中，针对中国大陆二维空间中AMR现象的普遍性研究认为：在空间上具有AMR加速异常现象的震例似乎分布在几个特定的地块或边界上；对于不同时间和地点的强震，甚至同一地点的不同时间的地震，出现AMR异常的时空范围不相同，具有复杂性，即表征应变释放程度的m值在时空的二维空间呈现复杂分布。考虑到在AMR方法的实际应用中，特别是在预测效果上，存在如何选取合理的空间扫描范围等问题，通过震例的回顾性研究是探索相关问题、评价AMR方法用于预测的有效性、进而开展应用研究的一种重要途径。 1976年唐山7.8级地震、1998年张北6.2级地震、2006年文安5.1级地震是华北北部地区有较好地震监测能力以来发生的三次重要的破坏性地震。对这三次事件开展回顾性的震例研究，有助于华北及首都圈地区的地震预测探索。本文通过分析这三次地震之前的应变释放现象在不同时间尺度、不同空间范围所呈现的特征，归纳出主震前特定时间段能表征应变释放（AMR）加速特征的优势空间扫描半径；进而围绕晋冀蒙交界地区（40°～41°N,113°～115°E）和首都圈地区（38.5°～41°N,114°～120°E）进行计算与分析，以期将地震活动定量分析的AMR方法应用于这两个地区的未来强震危险性分析与预测。 1 方法和数据 1.1 声发射应变加速释放特征 利用断裂力学和裂隙扩展理论得到、经Bufe和Varnes整理、Brehm和Braile改进的“破裂时间”方程为 式中，tf是未来主震的发震时刻，Ω为地震活动性参量（能量E、标量地震矩或频次），t是地震前兆事件的时间，k是表征区域地震活动水平的常数。m1为加速释放，m=1表征释放为稳定的线性；m1为减速释放。Bufe和Varnes曾提出，对一给定地区，m值可由以前发生的主震的前震加速序列来估算。式（1）中能量E与震级MS的转换关系为 结合断裂力学以及地震矩与断裂长度的关系: 得到地震矩M0与震级MS的转换关系为 其中，M0的单位是N·m。式（4）中A,B为未知常量。 针对式（1）中m值的意义，蒋海昆等在实验室中考察了不同温度、压力条件下岩石变形过程中的声发射数据。他们观察到当围压固定时，m值随温度升高逐渐变大，声发射应变从加速释放逐渐过渡到匀速释放，这意味着不同温度条件下岩石变形过程中内部微破裂形式的差异，可能导致应变释放类型的较大差异；在地壳浅部的温压条件下，岩石破坏前会显示一定的微破裂加速释放特征，此时的m值小于1；在渐进式破坏区段，应变释放呈逐渐减弱的减速释放态势，m明显大于1；在深部温压条件下，应变释放加速特征明显，m值明显较低；完整岩样破裂前声发射应变加速释放特征显著，而宏观剪裂面的黏滑之前，声发射应变基本上匀速释放。假定临界区域是圆形，Bowman等给出一种确定临界区域的准则：对于临界区域最后的估计是该区域内的数据能被幂律很好地拟合。 选用累积贝尼奥夫应变（CBS）作为地震释放的量度，采用公式（5）先计算出地震矩，然后取地震矩平方根的累积作为变量，其正比于累积贝尼奥夫应变。为了量化贝尼奥夫应变的加速程度，Bowman定义了曲率参数c,c是幂律拟合均方差与线性拟合均方差之比，其中最小c值所对应的半径即为临界半径。经验关系为 式中r为临界半径，M为相应的地震震级，c,d是常数。 根据以上的相关理论与方法，在震例的回顾性研究中，将已发生的主震震中位置作为圆心，使用不同的半径连续计算出式（6）的幂律拟合残差与线性拟合残差的最小比值c和临界半径。 1.2 主要震例的预处理 Brehm和Braile在研究新马德里地震带时，针对应变释放加速地震序列的判别，提出了选用地震目录资料的两个准则：主震前的事件在适当的时间段内是完整的（即震级-频度关系是线性的），其完整记录必须包括震级比主震震级小2个震级单位的所有事件。本研究遵循这两个准则。 本研究从中国地震台网地震目录中挑选出1970年1月至2012年12月河北及邻区（112°～121°E,34°～43°N）范围的地震目录，并对挑选出的地震目录进行完整性分析。结果显示2002年以来ML≥1.6地震是完整的，