当前，地震震源参数快速反演仍是一个具有挑战性的科学研究问题，主要有两方面原因：（1）具有复杂破裂过程的大震激发的高频波形较难模拟；（2）地震台网覆盖不均匀、台站仪器参数不准确以及波形记录信噪比低等因素导致地震矩张量自动反演实现起来困难。为了解决上述问题，中科院地质地球所地球与行星物理院重点实验室地震学研究团队的赵旭博士后、姚振兴研究员，联合法国斯特拉斯堡大学的合作者，进一步发展了W震相（W-phase, 图1）波形快速反演地震震源参数方法。研究了基于近震W震相的新时窗和合适滤波频段，运用并行技术，提升计算效率。
图1 2016年印尼Mw7.8地震GNI台站垂直向记录的W震相波形。波形附近两个红色圆圈表示选用的时窗。黑线为观测记录，红线为理论波形
研究结果表明（图2和图3）：相对于原方法（Kanamori and Rivera，2008），新方法能在震后4～7 分钟内准确测定M5~9级地震质心矩张量，将地震震级下限从Mw6.5降至Mw5.0，测定时间仅为原来1/5。
图2 基于近震W震相方法自动反演地震（M ≥ 5.0）质心矩张量。研究的地震事件共有147个，含2011年11月至2015年11月中国及其邻区地震（M ≥ 5.0）、2008汶川Mw7.9地震和2011年日本Mw9.1地震。震源球颜色根据相对于美国哥伦比亚大学Global CMT计算的角距离Φ大小表示（见色标）。选用台站共约850个（见图右上角）
图3 自动反演的矩震级（纵轴）同Global CMT给出的矩震级（横轴）比较。两者成良好的对角线分布，吻合较好
上述研究成果发表在国际地球物理学期刊Journal of Geophysical Research: Solid Earth 上（Zhao et al. Regional W-phase source inversion formoderate to large earthquakes in China and neighboring areas. Journal of Geophysical Research: Solid Earth. 2017,122: 10052-10068)。
原文链接
