Numerical Solution of the 3-D Travel Time Tomography Problem

Michael V. Klibanova Jingzhi Li,Wenlong Zhang

DOI: https://doi.org/10.48550/arXiv.2209.09420

2022-09-20

Abstract:The first numerical solution of the 3-D travel time tomography problem is presented. The globally convergent convexification numerical method is applied.

Numerical Analysis

What problem does this paper attempt to address?

Convexification for the 3D Problem of Travel Time Tomography

Michael V. Klibanov,Jingzhi Li,Vladimir G. Romanov,Zhipeng Yang

2024-09-21

Abstract:The travel time tomography problem is a coefficient inverse problem for the eikonal equation. This problem has well known applications in seismic. The eikonal equation is considered here in the circular cylinder, where point sources run along its axis and measurements of travel times are conductes on the whole surface of this cylinder. A new version of the globally convergent convexification numerical method for this problem is developed. Results of numerical studies are presented.

Numerical Analysis
Numerical inversion of 3D geodesic X-ray transform arising from traveltime tomography

Tak Shing Au Yeung,Eric T. Chung,Gunther Uhlmann

DOI: https://doi.org/10.48550/arXiv.1804.10006

2018-04-26

Abstract:In this paper, we consider the inverse problem of determining an unknown function defined in three space dimensions from its geodesic X-ray transform. The standard X-ray transform is defined on the Euclidean metric and is given by the integration of a function along straight lines. The geodesic X-ray transform is the generalization of the standard X-ray transform in Riemannian manifolds and is defined by integration of a function along geodesics. This paper is motivated by Uhlmann and Vasy's theoretical reconstruction algorithm for geodesic X-ray transform and mathematical formulation for traveltime tomography to develop a novel numerical algorithm for the stated goal. Our numerical scheme is based on a Neumann series approximation and a layer stripping approach. In particular, we will first reconstruct the unknown function by using a convergent Neumann series for each small neighborhood near the boundary. Once the solution is constructed on a layer near the boundary, we repeat the same procedure for the next layer, and continue this process until the unknown function is recovered on the whole domain. One main advantage of our approach is that the reconstruction is localized, and is therefore very efficient, compared with other global approaches for which the reconstructions are performed on the whole domain. We illustrate the performance of our method by showing some test cases including the Marmousi model. Finally, we apply this method to a travel time tomography in 3D, in which the inversion of the geodesic X-ray transform is one important step, and present several numerical results to validate the scheme.

Numerical Analysis
Study of Travel Time Tomography by Iterative Tikhonov Regularization Inversion

WANG Zhen-yu,LIU Guo-hua

DOI: https://doi.org/10.3785/j.issn.1008-973x.2005.02.019

2005-01-01

Abstract:Ill-posed problem is the bottleneck to solution of inversion. The singular value decomposition (SVD) technique was applied to analyze the characteristic of inversion equations. The iterative Tikhonov regularization method was studied to solve 2D acoustic travel time tomography. It was a quasi-linear inversion algorithm, and the L-curve method was proposed for finding the optimal regularization parameter to avoid aimless random trial calculations. Four convergence criteria were also suggested for guiding the inversion iterations. The tomography results of the numerical simulation model with an abnormal core showed that the solutions had optimal balance between the solution norm and the corresponding residual norm. Compared with the conventional simultaneous iterative reconstruction techniques (SIRTs), the iterative Tikhonov regularization method has better precision and convergence efficiency.
Implementation strategy of 3D reverse time migration on GPU/CPU clusters

Shouwei Liu,Huazhong Wang,Shengchang Chen,Xiangning Kong

DOI: https://doi.org/10.6038/cjg20131023

2013-01-01

Abstract:Prestack reverse time migration(RTM) is currently the most precise seismic imaging method. It needs appropriate implementation strategies and efficient computation platforms, because of the needs for large amount of calculation and large storage. Based on high order finite difference RTM, this paper focuses on solving the questions of the lack of video memory and artificial boundary, which must be solved for Graphic Processing Unit (GPU) RTM. Using the domain decomposition technique, we can efficiently realize three dimensional (3D) RTM of any production scale. It eliminates the limitation of GPU memory size. The control equation of the traditional perfectly matched layer (PML) boundary condition is different to that of the inner area. It is not suitable for GPU calculation. This paper realizes the Nearly PML (NPML) boundary condition on GPU. It makes high order finite difference calculation eliminate the needs of conditional branches. And the amount of auxiliary wavefields in boundary area is small. Therefore, wavefields propagation using NPML is efficient on GPU. 3D theoretical data and real field data RTM results prove the correctness of this paper's methods.
Study on Inversion Methods for Travel Time Tomography Based on Generalized Inverse Theory

王振宇,刘国华,梁国钱

DOI: https://doi.org/10.3785/j.issn.1008-973x.2005.01.001

2005-01-01

Abstract:Conventional inversion methods of computer tomography (IMCT) have low convergence rates, and cannot deal with prior information and incomplete projection problems. To improve their robustness, firstly, a new IMCT based on generalized inverse theory was proposed for solving the inversion equation system. Then, a multi-objective optimization problem was constructed and resolved using the ideal point method to avoid another bottleneck in the conventional IMCT, i.e., finding an optimum damping coefficient. Two numerical examples were modeled to validate the above-devised two measures. It was shown that the optimal methodology is insensitive to the initial values of damping coefficient, and that the proposed IMCT is more effective and accurate than the conventional ones.
3D Wave-Equation First-Arrival Tomography Using Monochromatic Traveltime Sensitivity Kernel

Feng Bo,Xu RongWei,Wang HuaZhong,Luo Fei

DOI: https://doi.org/10.1190/segam2021-3594587.1

2023-01-01

Abstract:Seismic traveltime tomography is a key technique for near-surface velocity model building. Specifically, the wave-equation tomography can correctly handle the first-order diffraction phenomena when the scale of velocity heterogeneity is comparable with the incident wavelength. Although wave-equation-based traveltime inversion can theoretically produce a velocity model with higher resolution compared with ray tomography, the computational cost is much higher which prohibits its implementation for industry-scaled 3-D problems. To address this problem, by introducing the random boundary condition and the discrete Fourier integration, we propose a practical 3-D waveequation first -arrival tomography using the monochromatic gradient inversion scheme. The advantages of the proposed method are: in terms of computational cost and storage space, the proposed method only needs very few memory spaces (i. e., compared with wave-equation forward modeling, only two additional monochromatic wavefields and one image volume are required), and meanwhile, the computational cost is decreased at least by 1/3 (i. e., compared with the wavefield reconstruction scheme, only two wavefield propagations are needed, and no special treatment for numerical boundaries). Besides, the introduction of random boundary conditions further simplifies the modeling algorithm which in turn makes it well suited for GPU (Graphics Processing Unit) acceleration. Theoretically, considering that the major area that may influence the first -order diffraction of seismic waves is usually referred to as the first Fresnel zone, which has a very similar shape to the monochromatic sensitivity kernel calculated with the dominant frequency. On one hand, the monochromatic sensitivity kernel within the first Fresnel zone will be coherently increased because the scattered wavefield caused by random boundaries becomes less spatial correlative. On the other hand, the multi-fold seismic acquisition system makes the sensitivity kernels outside the first Fresnel zone further canceled. Numerical example demonstrates that the smoothed monochromatic gradient contains sufficient low-wavenumber components which are adequate for background velocity inversion. Moreover, a 3-D synthetic data test demonstrates that a high-resolution near-surface velocity model can be inverted using wave-equation traveltime tomography with affordable computer resources. Compared with the traditional wave-equation tomography methods, the proposed method not only reduces the memory requirement remarkably but also significantly increases the computational efficiency. Therefore, it is very promising for dealing with industrial-sized 3-D inverse problems.
Traveltime calculations for qP, qSV, and qSH waves in 3D tilted transversely isotropic media using fast sweeping method with a Newton method local solver

Yongming Lu,Wei Zhang

DOI: https://doi.org/10.22541/essoar.167336515.56790498/v1

2023-01-01

Abstract:Efficient and accurate traveltime calculations of seismic waves have important applications in tomography, prestack migration, earthquake location, etc. Anisotropy significantly affects the traveltimes of seismic wave. For high-resolution imaging and inversion, it is necessary to consider anisotropy in the traveltime calculations. The fast sweeping method (FSM) does not need to track and store the minimum traveltime point of wavefront, which has important applications in computing the anisotropic first-arrival traveltime. The conventional method that solves a transformed traveltime quartic equation combined with FSM is suitable for general anisotropic media. However, this method involves finding the intervals of roots and solving the quartic equations using bisection algorithm at each iteration, leading to high computational cost and instability for the 3D problems. In our previous work, for the vertical transversely isotropic (VTI) case, we developed an FSM to compute the qP-wave first-arrival traveltimes by analytically solving the simplified quadratic slowness equation in a specific triangular-pyramid stencil. This method greatly improves the computational efficiency. However, for the qP and qSV waves, analytically solving the slowness equation cannot be extended to tilted transversely isotropic (TTI) media. To address this problem, we introduced the Newton method in the triangular-pyramid local solver to quickly solve the TTI slowness equation. For the qSH wave, its slowness equation is quadratic and simple to solve. The proposed method provides an efficient procedure for the traveltime calculations of qP, qSV, and qSH waves in 3D general TTI media. Numerical examples have verified the efficiency and accuracy of the proposed method.
Wave-equation First-Arrival Tomography for 3-D Near-Surface Velocity Inversion

B. Feng,R. Wu,H. Wang,C. Zhang,R. Xu

DOI: https://doi.org/10.3997/2214-4609.202113178

2021-01-01

Abstract:Summary In wave-equation traveltime inversion, the major computation comes from the calculation of the functional gradient or Fréchet derivative. Considering that the goal of traveltime inversion is to recover the large-scale velocity model (i.e., the low-wavenumber component), we propose a computationally economical method to calculate the functional gradient using the monochromatic traveltime sensitivity kernel, which considerably reduces the computational cost as well as memory requirement. 3-D field test demonstrates that a high-resolution near-surface velocity model can be inverted using wave-equation traveltime tomography with affordable computer resources.
An Iterative Kalman Filter-based Method for Traveltime Tomography

Yun Zhang,Xihai Li,Tian Liu,Jihao Liu,Chaoying Bai

DOI: https://doi.org/10.1093/jge/gxae048

2024-04-18

Journal of Geophysics and Engineering

Abstract:In addition to directly determining the inverted velocity, it is crucial to assess the uncertainty (or confidence) of the relevant velocity imaging results. This is especially true when dealing with time-lapse seismic data, which play a crucial role in detecting changes in fluid movement in petroleum reservoirs. In this study, we developed an iterative Kalman filter-based method to monitor time-lapse changes in the subsurface using traveltime tomography. Our approach considers the uncertainty of the results. We successfully verified the validity of our proposed method through synthetic crosshole and time-lapse tests. Our iterative Kalman filter-based method provides reasonable reconstruction results suitable for accessing uncertainties, and the algorithm can be used to monitor changes in the subsurface medium.

geochemistry & geophysics
Adjoint-state Traveltime Tomography for Azimuthally Anisotropic Media in Spherical Coordinates

Jing Chen,Guoxu Chen,Masaru Nagaso,Ping Tong

DOI: https://doi.org/10.1093/gji/ggad093

IF: 3.352

2023-01-01

Geophysical Journal International

Abstract:SUMMARYTong has proposed an adjoint-state traveltime tomography method to determine velocity heterogeneity and azimuthal anisotropy. This method, however, ignores the Earth’s curvature when deriving the eikonal equation for azimuthally anisotropic media. Thus, further coordinate transformation or approximation is required to ensure the accuracy of traveltime prediction in large-scale tomography. To address this problem, we derive the eikonal equation for azimuthally anisotropic media in spherical coordinates, which naturally considers the Earth’s curvature. Another key ingredient is the forward modelling algorithm, whose accuracy and efficiency dominate the numerical error and computational cost of the inversion. In this study, we apply a modified fast sweeping method to solve the eikonal equation in spherical coordinates. Two approaches, including the third-order weighted essentially non-oscillatory approximation and multiplicative factorization technique, are applied to improve the accuracy. According to the numerical experiments, this new eikonal solver achieves a second-order accuracy and is about two orders of magnitude more accurate than the commonly used first-order fast sweeping method with similar runtime. Taking advantage of the two improvements, we develop a novel eikonal equation-based adjoint-state traveltime tomography method for azimuthally anisotropic media in spherical coordinates. This method is applicable for large-scale tomography, and its performance is verified by a synthetic checkerboard test and a practical seismic tomographic inversion in central California near Parkfield.
Travel time tomography by ray tracing using the fast sweeping method

DOI: https://doi.org/10.1007/s11770-024-1095-6

2024-05-17

Applied Geophysics

Abstract:This paper presents a gradient-descent travel time tomography method for solving the acoustic-type velocity model inversion problem. Similarly to the adjoint-state method, the proposed method is based on the Eikonal equation, enabling simultaneous calculation of contributions from all common-source receivers to the gradient. This overcomes the inefficiency inherent in conventional travel time tomography methods, which rely on a two-point ray tracing process. By directly calculating Fréchet derivatives, our method avoids the complex derivation processes associated with the adjoint-state method. The key to calculating the Fréchet derivatives is to calculate a so-called ray-path term. Consequently, compared to the adjoint-state method, the proposed method can explicitly obtain the ray paths, resulting in a more concise and intuitive derivation process. Furthermore, our method retains the benefits of the adjoint-state method, such as speed, low memory usage, and robustness. This paper focuses on elucidating the principles and algorithms for calculating the ray-path term based on the fast sweeping method. The algorithms could be further speeded up by using parallel computational techniques. Synthetic tests demonstrate that our proposed travel time tomographic method accurately calculates ray paths, regardless of the complexity of the model and recording geometry.

geochemistry & geophysics
Template-based CT Reconstruction with Optimal Transport and Total Generalized Variation

Yiming Gao,Zhengmeng Jin,Xu Li

DOI: https://doi.org/10.1088/1361-6420/aceb17

IF: 2.408

2023-01-01

Inverse Problems

Abstract:x-ray computed tomography (CT) has been widely used in clinical diagnosis as a modality of medical imaging. To decrease the radiation dose patients suffering from, sparse-view CT has gained much attention in medical imaging field. In this paper, we propose to design a variational model based on dynamic optimal transportation and total generalized variation for CT reconstruction. This is a joint task involving inverse problem and template registration. The final state image of the optimal transport problem is unknown and needs to be reconstructed through CT inversion, while the given initial state can be regarded as a template to provide some structural information for the final one. Moreover, the existence and stability of minimizers to our proposed model are shown in continuous space. In discretization with the continuity equation, we utilize the well-known staggered grid in fluid mechanics and develop a first-order algorithm based on primal-dual method for numerically solving the proposed model. Finally, numerical experiments for sparse-view CT reconstruction are exhibited to show the performance of our proposed model in recovering images with high quality and structure preservation.
Wave-equation Based Travel-time Tomography for First Arrival

B. Feng,H.Z. Wang,X. Hu

DOI: https://doi.org/10.3997/2214-4609.20141575

2014-01-01

Proceedings

Abstract:Summary The conventional travel-time sensitivity kernel (TSK) is derived using the first-order Born approximation. However, the first-order Rytov approximation is more suitable for the first arrival tomography due to its better behaving for the forward scattering. In this paper, we present a wave-equation based travel-time tomography method, which is derived from the first-order Rytov-approximation to the wave equation. The travel-time sensitivity kernel can be calculated using the well-known RTM-like algorithm. The L2-norm of travel-time residuals is minimized by the Gauss-Newton algorithm. The presented method is applied to the first arrival travel-time inversion. Synthetic examples show that the long-wavelength part of the velocity model can be recovered well.
Microlocal Methods for The Elastic Travel Time Tomography Problem for Transversely Isotropic Media

Yuzhou Zou

DOI: https://doi.org/10.48550/arXiv.2112.14455

2021-12-29

Abstract:We study the travel time tomography problem for transversely isotropic media using the artificial boundary method used in works in X-ray tomography and boundary rigidity problems. The analysis involved requires a modification of the scattering pseudodifferential calculus in order to provide parametrices for operators which are degenerately elliptic; such a calculus is constructed in this paper and used to solve the tomography problem of interest.

Analysis of PDEs,Differential Geometry
Calculating Traveltimes in 2D General Tilted Transversely Isotropic Media Using Fast Sweeping Method

Yongming Lu,Wei Zhang,Jianfeng Zhang

DOI: https://doi.org/10.1016/j.cageo.2024.105724

2024-01-01

Abstract:Traveltime calculations play an important role in the field of exploration seismology, such as traveltime tomography and seismic imaging and so on. Seismic anisotropy poses a challenge for traveltime calculation, because anisotropic eikonal solvers are more complex than the isotropic counter part. To solve the eikonal equations in 2D tilted transversely isotropic (TTI) media, we have developed a fast algorithm combine with fast sweeping method to compute the first arrival traveltimes of quasi-P (qP)-, quasi-SV (qSV)-, and quasi-SH(qSH)-waves. For the qP- and qSV-waves, we analyzed the quartic coupled slowness surface equation derived from the Christoffel equation. Then, we constructed a local solver to relate traveltime and slowness. We found that in the local solver, one component of the slowness vector is known and the corresponding slowness equation is monotonic. This provides a strong basis for the fast iterative algorithm we proposed, where we use the Newton method to solve the qP- and qSV-wave slowness equation to determine the related traveltimes. For the qSH wave, the slowness equation is quadratic and simple to solve. Numerical experiments demonstrate that the proposed method can obtain accurate traveltimes for simple and complicated 2D TTI models.
3D Weak-Dispersion Reverse Time Migration with a Stereomodeling Method

Li Jingshuang,Dinghui Yang,Xueyuan Huang,Mike Fehler

DOI: https://doi.org/10.1190/segam2013-0996.1

2013-01-01

Abstract:The finite difference method has been widely used in seismic modeling and reverse time migration. However, it generally has two issues: large computational cost and numerical dispersion. Recently, a nearly-analytic discrete operator was developed to approximate the partial differential operators. Based on this spatial discretization, many weak-dispersion and efficient StereoModeling methods have been developed, which are found to be superior to conventional algorithms in suppressing numerical dispersion. In this paper, we generalize one StereoModeling method, the nearly-analytic central difference method (NACD), from 2D to 3D and apply it to 3D reverse-time migration. Numerical results show that the NACD can be used effectively as a new tool for seismic modeling and migration. The reverse time migration (RTM) results for the 3D SEG/EAGE Phase A classic dataset 1 show that the NACD can get a much better image than the Lax-Wendroff correction (LWC) method particularly when using a coarse grid size.
Stability and statistical inversion of travel time tomography

Ashwin Tarikere Ashok Kumar Nag,Hanming Zhou

DOI: https://doi.org/10.1088/1361-6420/ad4911

IF: 2.408

2024-05-10

Inverse Problems

Abstract:In this paper, we consider the travel time tomography problem for conformal metrics on a bounded domain, which seeks to determine the conformal factor of the metric from the lengths of geodesics joining boundary points. We establish forward and inverse stability estimates for simple conformal metrics under some a priori conditions. We then apply the stability estimates to show the consistency of a Bayesian statistical inversion technique for travel time tomography with discrete, noisy measurements.

mathematics, applied,physics, mathematical
A fast local level set adjoint state method for first arrival transmission traveltime tomography with discontinuous slowness

wenbin li,shingyu leung

DOI: https://doi.org/10.1093/gji/ggt244

IF: 3.352

2013-01-01

Geophysical Journal International

Abstract:We propose a numerical algorithm for solving first arrival transmission traveltime tomography problems where the underlying slowness is piecewise continuous. The idea is based upon our previously efficient approach for smooth slowness inversion (Leung & Qian) using the fast sweeping method and the adjoint state method. In this work, we further incorporate the level set method to implicitly represe...
Mathematical Modeling and Solution of the Tomography Problem in Domains with Reflecting Obstacles

Kamen Lozev

DOI: https://doi.org/10.48550/arXiv.1209.0130

2012-09-02

Abstract:This work develops new numerical methods for the solution of the tomography problem in domains with reflecting obstacles. We compare the solution's performance for Lambertian reflection, for classical tomography with ubroken rays and for specular reflection. Our numerical method using Lambertian reflection improves the solution's accuracy by an order of magnitude compared to classical tomography with ubroken rays and for tomography in the presence of a specularly reflecting obstacle the numerical method improves the solution's accuracy approximately by a factor of three times. We present efficient new algorithms for the solution's software implementation and analyze the solution's performance and effectiveness.

Numerical Analysis
Practical wave-equation first-arrival tomography: Generalized Rytov approximation, double-difference measurement, and an efficient Gauss-Newton algorithm

Bo Feng,Ru-Shan Wu,Huazhong Wang

DOI: https://doi.org/10.1190/SEGAM2020-3427790.1

2020-01-01

SEG Technical Program Expanded Abstracts

Abstract:PreviousNext No AccessSEG Technical Program Expanded Abstracts 2020Practical wave-equation first-arrival tomography: Generalized Rytov approximation, double-difference measurement, and an efficient Gauss-Newton algorithmAuthors: Bo FengRu-Shan WuHuazhong WangBo FengTongji UniversitySearch for more papers by this author, Ru-Shan WuUniversity of CaliforniaSearch for more papers by this author, and Huazhong WangTongji UniversitySearch for more papers by this authorhttps://doi.org/10.1190/segam2020-3427790.1 SectionsAboutPDF/ePub ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InRedditEmail AbstractIn the conventional finite-frequency traveltime tomography, the traveltime Fréchet derivative is usually derived using the Born- or the Rytov approximation, which requires that the magnitude of velocity perturbations should be small enough to satisfy the weak-scattering assumption. For large-scale and strong velocity perturbations, the generalized Rytov approximation (GRA) has been proposed, which can predict the accumulated phase-change of the forward scattered wave, making it more suited for traveltime inverse problems. To mitigate the source uncertainties, we introduce the doubledifference (DD) approach to calculate the traveltime shift. By combining GRA with the DD approach, we propose a new finite-frequency traveltime sensitivity kernel, i.e. the GRA-based double-difference traveltime sensitivity kernel, which inherits the advantage of both the DD measurement and the GRA method. To solve the associated traveltime inverse problem, we use a matrix-free Gauss-Newton algorithm to accelerate the convergence rate, where the Hessian-vector product is computed approximately using the ray approximation, making it more computationally efficient. Therefore, the proposed method is quite promising for 3-D problems.Presentation Date: Monday, October 12, 2020Session Start Time: 1:50 PMPresentation Time: 2:15 PMLocation: Poster Station 8Presentation Type: PosterKeywords: tomography, inversion, traveltimePermalink: https://doi.org/10.1190/segam2020-3427790.1FiguresReferencesRelatedDetails SEG Technical Program Expanded Abstracts 2020ISSN (print):1052-3812 ISSN (online):1949-4645Copyright: 2020 Pages: 3887 publication data© 2020 Published in electronic format with permission by the Society of Exploration GeophysicistsPublisher:Society of Exploration Geophysicists HistoryPublished Online: 30 Sep 2020 CITATION INFORMATION Bo Feng, Ru-Shan Wu, and Huazhong Wang, (2020), "Practical wave-equation first-arrival tomography: Generalized Rytov approximation, double-difference measurement, and an efficient Gauss-Newton algorithm," SEG Technical Program Expanded Abstracts : 3664-3668. https://doi.org/10.1190/segam2020-3427790.1 Plain-Language Summary KeywordstomographyinversiontraveltimePDF DownloadLoading ...

Numerical Solution of the 3-D Travel Time Tomography Problem

Convexification for the 3D Problem of Travel Time Tomography

Numerical inversion of 3D geodesic X-ray transform arising from traveltime tomography

Study of Travel Time Tomography by Iterative Tikhonov Regularization Inversion

Implementation strategy of 3D reverse time migration on GPU/CPU clusters

Study on Inversion Methods for Travel Time Tomography Based on Generalized Inverse Theory

3D Wave-Equation First-Arrival Tomography Using Monochromatic Traveltime Sensitivity Kernel

Traveltime calculations for qP, qSV, and qSH waves in 3D tilted transversely isotropic media using fast sweeping method with a Newton method local solver

Wave-equation First-Arrival Tomography for 3-D Near-Surface Velocity Inversion

An Iterative Kalman Filter-based Method for Traveltime Tomography

Adjoint-state Traveltime Tomography for Azimuthally Anisotropic Media in Spherical Coordinates

Travel time tomography by ray tracing using the fast sweeping method

Template-based CT Reconstruction with Optimal Transport and Total Generalized Variation

Wave-equation Based Travel-time Tomography for First Arrival

Microlocal Methods for The Elastic Travel Time Tomography Problem for Transversely Isotropic Media

Calculating Traveltimes in 2D General Tilted Transversely Isotropic Media Using Fast Sweeping Method

3D Weak-Dispersion Reverse Time Migration with a Stereomodeling Method

Stability and statistical inversion of travel time tomography

A fast local level set adjoint state method for first arrival transmission traveltime tomography with discontinuous slowness

Mathematical Modeling and Solution of the Tomography Problem in Domains with Reflecting Obstacles

Practical wave-equation first-arrival tomography: Generalized Rytov approximation, double-difference measurement, and an efficient Gauss-Newton algorithm