利用Lg波计算地震应力降的方法及在青藏高原东缘典型震例中的应用

地震应力降标志震源破裂过程中作用在断层系统上的应力释放水平,是刻画震源物理过程和预测震源辐射特性的重要参数.地震应力降主要受控于构造环境、震源机制和地震类型等.观测得到的应力降还受到观测频带的影响,所以它的绝对值难以准确测定.采用体波、面波以及零频观测(大地测量、GPS、InSAR)等不同类型的数据获得的应力降也存在差异.对于地震波数据,通常采用间接方法去除传播路径中的衰减效应并获得震源激发谱,进而估计应力降.随着宽频带地壳Lg波衰减模型的建立,可以构建直接校正路径衰减的方法,从而获得对Lg波震源激发函数的准确估计.进而通过对观测和理论震源谱的拟合获得整个地震序列中各个地震的标量地震矩、拐角频率和高频下降率等震源参数,再根据断层模型计算应力降.本文分别以青藏高原东缘典型的构造地震2017年Ms7.0九寨沟地震和潜在的工业注水诱发地震2019年Ms6.0长宁地震为例,计算两个地震序列应力降的时-空变化过程,探索构造地震与诱发地震之间的潜在的物理差异. 2017年九寨沟地震的主震应力降为27 MPa,余震震级和应力降均呈快速下降趋势. 2019年长宁地震的主震应力降为32 MPa,余震序列的应力降以起伏形式缓慢下降,并曾出现过与主震应力降数值相当的余震.这一现象主要来源于长宁地区长期的工业注水干扰了区域应力场,提升了断层内的流体静压力并降低了断层滑动的驱动应力.这一过程需经过震后较长的时间才能恢复.上述两个地震序列具有相似的主震应力降但属于完全不同的类型,说明在该区域无法通过单独观测主震应力降来区分构造与诱发地震,而通过研究整个序列中应力降的时-空发展过程则有可能揭示出与此有关的进一步信息.