A matching pursuit approach to the geophysical inverse problem of seismic travel time tomography under the ray theory approximation

Naomi Schneider,Volker Michel,Karin Sigloch,Eoghan J. Totten
2023-09-01
Abstract:Seismic travel time tomography is a geophysical imaging method to infer the 3-D interior structure of the solid Earth. Most commonly formulated as a linear(ized) inverse problem, it maps differences between observed and expected wave travel times to interior regions where waves propagate faster or slower than the expected average. The Earth's interior is typically parametrized by a single kind of localized basis function. Here we present an alternative approach that uses matching pursuits on large dictionaries of basis functions. Within the past decade the (Learning) Inverse Problem Matching Pursuits ((L)IPMPs) have been developed. They combine global and local trial functions. An approximation is built in a so-called best basis, chosen iteratively from an intentionally overcomplete set or dictionary. In each iteration, the choice for the next best basis element reduces the Tikhonov-Phillips functional. This is in contrast to classical methods that use either global or local basis functions. The LIPMPs have proven its applicability in inverse problems like the downward continuation of the gravitational potential as well as the MEG-/EEG-problem from medical imaging. Here, we remodel the Learning Regularized Functional Matching Pursuit (LRFMP), which is one of the LIPMPs, for travel time tomography in a ray theoretical setting. In particular, we introduce the operator, some possible trial functions and the regularization. We show a numerical proof of concept for artificial travel time delays obtained from a contrived model for velocity differences. The corresponding code is available at <a class="link-external link-https" href="https://doi.org/10.5281/zenodo.8227888" rel="external noopener nofollow">this https URL</a> under the licence CC-BY-NC-SA 3.0 DE.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is the geophysical inverse problem in seismic wave travel - time tomography. Specifically, the researchers hope to infer the three - dimensional structure inside the Earth and identify the areas that cause wave velocity changes by analyzing the time differences of seismic waves propagating inside the Earth. ### Problem Background Seismic wave travel - time tomography is an imaging method in geophysics used to infer the three - dimensional internal structure of the solid Earth. Traditional methods usually formulate this problem as a linear (or linearized) inverse problem, that is, mapping the differences between the observed wave travel - times and the expected travel - times to certain areas inside the Earth where the wave velocities are faster or slower than the expected average values. ### Limitations of Existing Methods Currently, the parameterization of the Earth's interior usually uses a single type of local basis function, which may lead to a model that is not flexible enough to capture the complex changes inside the Earth well. In addition, due to the ill - posed nature of the inverse problem itself (i.e., the non - uniqueness and instability of the solution), traditional linearized methods face many challenges when dealing with actual data. ### New Method Proposed in the Paper To overcome the above limitations, this paper introduces a new method based on Matching Pursuit, which is especially suitable for seismic wave travel - time tomography. This method utilizes multiple basis functions in a large dictionary and constructs an approximate solution by iteratively selecting the best basis functions. This method combines the advantages of global and local trial functions and can better adapt to the changes of the Earth's internal structure at different scales. ### Core Ideas of the Method 1. **Dictionary Construction**: Construct a dictionary containing multiple types of basis functions, including orthogonal polynomials and Finite - Element - based Hat Functions (FEHFs). These basis functions can flexibly represent the internal anomalies of the Earth at different scales and shapes. 2. **Iterative Selection**: In each iteration, select the basis function from the dictionary that can minimize the Tikhonov - Phillips functional to the greatest extent and gradually construct the optimal solution. 3. **Regularization**: To deal with the ill - posed nature of the inverse problem, a regularization term in the H1(Sobolev) space is introduced to ensure the stability and physical rationality of the solution. ### Application Verification The paper verifies the effectiveness of this method through numerical experiments, using synthetic whole - mantle test structures and actual seismic wave path data. The results show that the new method can provide more accurate estimates of the Earth's internal structure while maintaining computational efficiency. ### Summary This paper aims to improve the solution method of the geophysical inverse problem in seismic wave travel - time tomography by introducing the Matching Pursuit method, in order to obtain a more accurate and flexible model of the Earth's internal structure.