Compound Regularization of Full-waveform Inversion for Imaging Piecewise Media

Abstract

The nonlinear and ill-posed nature of full waveform inversion (FWI) requires us to use sophisticated regularization techniques to solve it. In most applications, the model parameters may be described by physical properties (e.g., wave speeds, density, attenuation, anisotropic parameters) which are piecewise functions of space. Compound regularizations are thus necessary to reconstruct properly such parameters by FWI. We consider different implementations of compound regularizations in the wavefield reconstruction inversion (WRI) method, a formulation of FWI that extends its search space and prevent the so-called cycle skipping pathology. Our hybrid regularizations rely on Tikhonov and total variation (TV) functionals, from which we build two classes of hybrid regularizers: the first class is simply obtained by a convex combination (CC) of the two functionals, while the second relies on their infimal convolution (IC). In the former class, the model of parameters is required to simultaneously satisfy different priors, while in the latter the model is broken into its basic components, each satisfying a distinct prior (e.g. smooth, piecewise constant, piecewise linear). We implement these types of compound regularizations in the WRI optimization problem using the alternating direction method of multipliers (ADMM). Then, we assess our regularized WRI in the framework of seismic imaging applications. Using a wide range of subsurface models, we conclude that compound regularizer based on IC leads to the lowest error in the parameter reconstruction compared to that obtained with the CC counterpart and the Tikhonov and TV regularizers when used independently.

References in corpus (4)

• A penalty method for PDE-constrained optimization in inverse problems
• Improving full waveform inversion by wavefield reconstruction with the alternating direction method of multipliers
• Priors with Coupled First and Second Order Differences for Manifold-Valued Image Processing
• ADMM-based multi-parameter wavefield reconstruction inversion in VTI acoustic media with TV regularization

Cited by in corpus (13)

• Full Waveform Inversion with Adaptive Regularization
• Extended full waveform inversion in the time domain by the augmented Lagrangian method
• ADMM-based multi-parameter wavefield reconstruction inversion in VTI acoustic media with TV regularization
• Complex-valued Imaging with Total Variation Regularization: An Application to Full-Waveform Inversion in Visco-acoustic Media
• Automatic Differentiation-based Full Waveform Inversion with Flexible Workflows
• Robust Wavefield Inversion via Phase Retrieval
• Efficient extended-search space full-waveform inversion with unknown source signatures
• Optimal Space-Variant Anisotropic Tikhonov Regularization for Full Waveform Inversion of Sparse Data
• Inexact Augmented Lagrangian Method-Based Full-waveform Inversion with Randomized Singular Value Decomposition
• Improved full-waveform inversion for seismic data in the presence of noise based on the K-support norm
• Simultaneous Automatic Picking and Manual Picking Refinement for First-Break
• ADMM-based full-waveform inversion for microseismic imaging
• Randomized source sketching for full waveform inversion