SeismoClip:服务于地震波形反演的交互式截窗程序

SeismoClip:An interactive time window selection tool for waveform inversion

  • 摘要: 地球内部介质的多尺度变化特征蕴含了板块构造动力学和地质演化的重要信息,因此获取可靠的地下介质参数模型是地震学研究的核心目标之一。相较于基于高频近似射线理论且使用震相走时信息少的成像方法,波形反演方法是利用地震波形记录中不同频率的相位和振幅等动力学信息来约束介质参数的,是目前构建高分辨度地下介质模型最具潜力的成像方法。从本质上分析,波形反演是非线性,在迭代过程中不加筛选的使用所有波形记录,容易导致迭代难以收敛或陷入局部极小值。所以,波形反演成功的关键之一是合理选取用于约束目标结构的波形数据时间窗。精准界定该时间窗口对提升数据有效信息的利用率、降低噪声干扰以及增强反演的准确性和稳定性具有至关重要的作用。为了高效便捷地提取有效波形信息,本研究开发了一个开源波形数据自动筛选程序SeismoClip。SeismoClip是一个基于MATLAB的交互式截窗程序,为波形反演而设计,它不仅可以作为独立的截窗工具,也可以作为其它自动截窗程序的后续处理工具。该程序通过提供阵列分析、多频段滤波、长短时窗比曲线等多种辅助分析工具及图件,大幅提高了截窗的准确度。此外,SeismoClip程序的设计非常注重用户体验,具备高度可定制性,支持用户根据具体需求添加或调整程序功能。

     

    Abstract:
    The multi-scale structure of the Earth’s interior contains important information on plate tectonics and geological processes, making the development of reliable models a critical objective in seismological research. Compared to conventional travel time tomography based on high-frequency ray theory and limited seismic phase arrival time information, waveform inversion utilizes the dynamical properties in the waveforms such as the frequency-dependent phases and amplitudes to constrain model parameters, thus representing the most promising approach for constructing high-resolution structural models. Waveform inversion is inherently nonlinear, thus blind use of all records without proper screening during iterative inversion leads to difficulties in convergence or getting trapped into local minima. Therefore, one of the key factors for successful waveform inversions is the selection of time windows of signals to constrain the intended structural targets. Proper selection of time windows is crucial for enhancing the utilization of effective information, reducing noise interference, and improving the accuracy and stability of inversions.
    To efficiently and conveniently extract effective waveform information, we develop the SeismoClip, an open-source waveform data automatic selection program. SeismoClip is a MATLAB-based interactive tool, designed for waveform inversion. It can be used not only as an independent window selector but also as a post-processing tool for other automatic window selection programs. Various auxiliary functions, such as array analysis, multi-frequency-band filtering, and short-time-average/long-time-average (STA/LTA) analysis, are also employed to significantly improve the accuracy and robustness of window selection performance. Furthermore, SeismoClip emphasizes user experience with extensive customizability, allowing researchers to make modifications to adapt program functions according to specific research needs. The core function of SeismoClip is to interactively process and select time windows from the input synthetic and observed seismic waveforms, and generate the corresponding window files. The processing flow is mainly divided into four steps: data processing, array analysis, interactive window selection, and automated adjustment.
    In recent years, the widespread deployment of seismic stations has led to the availability of huge volumes of waveform records from densely distributed seismic arrays. The waveform similarity of the same seismic phase between stations and the variation of phase arrival times with epicentral distance can help us determine a reasonable interval for the seismic phase window. Therefore, we have added an array-based module to the SeismoClip to assist in the analysis and extraction of phase windows. Since recorded waveforms are affected by noise and thus have a lower signal-to-noise ratios compared to theoretical seismograms, phase information in records can easily be obscured by noise. Therefore, SeismoClip performs array analysis based on the STA/LTA curve calculated from noise-free synthetic seismograms generated from the same source parameters.
    During the interactive selection process, SeismoClip displays both observed and synthetic waveforms in array form for reference. The interactive selection process is applied to single-channel waveforms. SeismoClip offers a variety of interactive selection interface options to accommodate different user needs, a including multi-frequency-band filtering single-channel selection interface, an interface highlighting weak-amplitude body-wave phases, and an interface based on existing selection results. This allows users to precisely extract seismic phase windows based on waveform data. For the interactively selected windows, SeismoClip performs a series of automated adjustments under the constraints of existing window selections to achieve a better fit. Users can choose either directly use the results of the interactive selection or allow the program to make automatic adjustments.
    In this article, we compare the selection results of SeismoClip with those from the automated time-window selection algorithm FLEXWIN, and discuss the distribution of sensitivity kernels corresponding to different time-windows to highlight the importance of precise time-window selection. Comparative analysis with real data shows that the time windows selected by SeismoClip are generally more accurate and reasonable than those selected by automated selection, particularly in selecting windows from waveforms with low signal-to-noise ratios.
    Overall, SeismoClip represents as a convenient and effective time-window selection tool for waveform inversion, based on array analysis of specific seismic phase tracking and multi-frequency-band analysis. Moreover, by real-time evaluation of the similarity and cross-correlation travel time differences between the observed and synthetic waveforms within the selected time window, it is possible to effectively avoid cycle skipping in waveform inversions. In cases where waveform features are not prominent, SeismoClip can be used to carefully select the time window to obtain appropriate sensitivity kernels and enhance the stability and accuracy of the inversions. For users utilizing an automated time-selection algorithm, we suggest starting with less restrictive parameters to conduct an initial window selection, and then utilizing the“interactive selection interface based on the eaisting selection results”of SeismoClip to review and refine the initial selection results.

     

/

返回文章
返回