地震定位与地震破裂面的ssa研究.docxVIP

地震定位与地震破裂面的ssa研究 互相关进行的地震定位法 精确确定震源位置和发震时间是地震研究中最重要的问题。 通常的方法是计算各个台站、 各个震相的理论到时和观测到时之差(走时残差)的最小值而确定震源位置的, 尽管理论地震学和计算机硬件都有了很大的发展, 但是大多数台网的定位依然离不开震相拾取的方法。 读取到时的方法至少有几个方面的缺点: 第一, 信噪比较低时到时的震相容易误认或者根本找不到; 第二, 当相同震源不同位置和时间但是相距很近的多事件发生时, 各个事件独立震相的提取是非常困难的; 第三, 台站数量较多时, 人们往往要付出很多时间对台站进行选取, 不便于大批量的计算。 总之, 震相拾取的方法在单独事件容易分离的情况和具有清晰的到时并且观测台站具有较低的背景噪声的情况下才能理想地定位。 通过使用互相关方法确定地震对的到时差, 利用双差法或主事件定位法, 对于地震丛集事件的定位得到了很大的改进,, 然而, 互相关方法会把原来并不相关的波形段进行关联。 通过波形的反演可以较好地改进定位的时间和空间位置, 其假设条件是完整的震源结构可以从不同方向和不同距离遭受干扰的各个台站得到恢复。 然而, 波形反演方法仅仅应用于非常有限的一些大地震(MW7), 这很大原因是因为近场地壳介质复杂性不能不加以考虑, 于是很难精确计算在近场范围内的短周期理论地震图。 另外, 波形反演方法要求预先假设断层的规模和产状, 但是对于很小的地震通常不可能获得这些相关的数据。 震源扫描算法(Source-Scanning Algorithm, SSA)是一种新的可以对震源分布进行成像的方法, 这种方法不需要预先知道断层产状和规模, 它可以方便地利用包括相对振幅和到时等波形数据, 检测到给定的时间和位置是否有地震震源。 通过系统地扫描一个试探性的位置范围和起震时刻, 可以发现整个震源序列分布而不需要拾取到时又无需计算理论地震图。 震源扫描算法2004年由Kao和Shan提出, 他们在对2003年3月北Cascadia俯冲带的慢地震研究中取得了比较满意的效果。 当时, 在一个小时的时间里发生了几个“tremor”事件, 震动延续的事件很长, 振幅与噪声水平接近, 因而不可能去辨别同一震源的震相, 传统的定位方法根本无法应用, 由于对震源的未知, 波形反演方法也不能实施。 最后是通过SSA方法根据相应的“亮度”函数扫描, 找出了“tremor”的震源位置和发生时间。 1 人工震源区的处理方法 震源扫描算法是在整个时空寻找可能震源的方法, 其基本原理可以表示为图1所示: 假设一个事件由N个台站记录到(图1中的A、B、C点), 首先分别归一化每个数字地震台站的记录un, 然后计算某个点(η)某个时刻(τ)的“亮度”函数 br(η,τ)=1N∑n=1N|un(τ+tηn)|(1)br(η,τ)=1Ν∑n=1Ν|un(τ+tηn)|(1) 式中,un为归一化地震记录,tη n为从点η处到台站n计算的某个最大震相的走时, 如果所有的最大振幅都是由点η和时间τ产生, 那么br(η,τ)=1。 同样,br(η,τ)=0.1意味着此处只有10%的可能性为震源。 在通常情况下,br(η,τ)=1的情况很少见, 可以认为整个时空中“亮度”函数取最大时是震源发生的位置和事件发生的时刻。 由于震源区速度模型不太清楚会造成走时计算有一些偏差, 所以必须把数据中相临近的点考虑进去, 不过, 可以通过光滑的方法来预先处理记录, 得到每个台站光滑的“亮度”函数序列 brn=smooth(|un|)(2)brn=smooth(|un|)(2) smooth表示一种光滑方法, 这样, 计算对于每点的亮度函数时只需要将各个台站相对应亮度函数相加平均即可。 由于每点的亮度函数是依次排列的, 这样下一段时间点的亮度函数可以一次求出, 这样做也免除了重复计算亮度函数, 减小了计算量, 从而有效地提高了扫描速度, 通过程序的改进, 可以使速度提高312倍。 改善算法以后的对于某一空间点(η)的某时刻(τ)的亮度函数可以表示为 br(η,τ)=1N∑n=1Nbrn(τ+tηn)(3)br(η,τ)=1Ν∑n=1Νbrn(τ+tηn)(3) 2 波场特征记录 为了检查SSA的稳定性和对相近震源的分辩能力, 笔者做了以下的震源验证试验。 假设区域有9个台站STA1～STA9, 台站之间经度和纬度的间距均为0.5°, 在某一时刻同时发生了“五角星”所示的5个事件(图2), 5相邻事件的经度和纬度之差为0.01°, 用一个半周期为0.2 s的半正弦波来模拟事件发出的信号。 在9个台站将记到由5个台网内的事件叠加在一起的9条波形记录(图3), 假设波速为4 km, 扫描时间间隔为0.02 s, 扫描的间距经度和纬