用震源扫描算法(SSA)进行微震的定位.pdfVIP

第26 卷 第 3 期 地 震 Vol.26 No.3 2006 年 7 月 EARTHCUAKE Jul. 2006 用震源扫描算法 ! SSA 进行微震的定位 李文军 ! 陈棋福 中国地震局地震预测研究所 北京 100036 摘要 # 引入了一种新的方法 震源扫描算法 Source Scannin g Al g orithm SSA 来对地震定 位和地震破裂面进行研究 这种方法充分利用了数字地震波形资料 在不用精确抬取到时和计 算理论地震图的情况下达到比较理想的定位效果 首先 文中阐述了 2004 年由 Kao 等人发展 的 SSA 方法的原理 其次对扫描算法的分辨率进行了理论上的模拟 用 5 点试验证明了 SSA 的很强的分辩能力 最后将 SSA 应用于南北地震带微震的定位试验中去 发现它能够在没有 人工干预的情况下较好地定出地震的位置 理论模拟和实际定位结果都证明 SSA 是普通微 小地震定位的一种好工具 关键词 # 震源扫描算法 地震定位 南北地震带 中图分类号 # P315.6 文献标识码 # A 文章编号 # 1000-3274 2006 03-0107-09 引言 精确地确定震源位置和发震时刻是地震研究中的一个最基本问题 通常的方法是计算 各个台站 各个震相的理论到时和观测到时之差 走时残差 的最小值而确定震源位置的 尽管理论地震学和计算机硬件都有了很大的发展 但是大多数台网的定位依然离不开震相 抬取的方法 读取到时的方法至少有几个方面的缺点 第一 信噪比较低时到时的震相容易误认或 者根本找不到 第二 当相同震源不同位置和时间但是相距很近的多事件发生时 各个事 件独立震相的提取是非常困难的 第三 台站数量较多时 人们往往要付出很多时间对台 站进行选取 不便于大批量的计算 总之 震相抬取的方法在单独事件容易分离的情况和 具有清晰的到时并且观测台站具有较低的背景噪声的情况下才能理想地定位 通过使用互相关方法确定地震对的到时差 利用双差法或主事件定位法 对于地震丛 集事件的定位得到了很大的改进 1 2 然而 互相关方法会把原来并不相关的波形段进行 关联 通过波形的反演可以较好地改进定位的时间和空间位置 其假设条件是完整的震源 结构可以从不同方向和不同距离遭受干扰的各个台站得到恢复 3 然而 波形反演方法仅 仅应用于非常有限的一些大地震 M W # 7 这很大原因是因为近场地壳介质复杂性不能不 收稿日期 # 修改回日期 #基金项目 # 地震科学联合基金课题 106010 作者简介 # 李文军 1969 男 湖南怀化人 2005 年获固体地球物理学博士学位 主要从事地震学等研究 加以考虑 9 于是很难精确计算在近场范围内的短周期理论地震图 O 另外 9 波形反演方法要 求预先假设断层的规模和产状 9 但是对于很小的地震通常不可能获得这些相关的数据 O 震源扫描算法 ( SourceScannin g Al g orithm 9 SSA 是一种新的可以对震源分布进行成 像的方法 9 这种方法不需要预先知道断层产状和规模 9 它可以方便地利用包括相对振幅和 到时等波形数据 [ 4 ] 9 检测到给定的时间和位置是否有地震震源 O 通过系统地扫描一个试探 性的位置范围和起震时刻 9 可以发现整个震源序列分布而不需要抬取到时又无需计算理论 地震图 O 震源扫描算法 2004 年由 Kao 和 Shan 提出 9 他们在对 2003 年 3 月北 Cascadia 俯冲带 的慢地震研究中取得了比较满意的效果 O 当时 9 在一个小时的时间里发生了几个 6 tremor 事件 9 震动延续的事件很长 9 振幅与噪声水平接近 9 因而不可能去辨别同一震源的震相 9 传统的定位方法根本无法应用 9 由于对震源的未知 9 波形反演方法也不能实施 O 最后是通 过 SSA 方法根据相应的 6 亮度 函数扫描 9 找出了 6 tremor 的震源位置和发生时间 O 1 震源扫描算法 ( SSA 的原理 震源扫描算法是在整个时空寻找可能震源的方法 9 其基本原理可以表示为图 1 所示 : 图 1 震源扫描算法原理示意图 ( 据 Kao 和 San 9 2004 空间上的点 ! 在 时刻的 6 亮度 可以由相应台站的理论到时 ( 即 加上相应的走时 t a ! ~ t b ! ~ t c ! 计算而 来 9