ADMM-based full-waveform inversion for microseismic imaging

Full waveform inversion (FWI) is beginning to be used to characterize weak seismic events at different scales, an example of which is microseismic event (MSE) characterization. However, FWI with unknown sources is a severely underdetermined optimization problem, and hence requires strong prior information about the sources and/or the velocity model. The frequency-domain wavefield reconstruction inversion method (WRI) has shown promising results to mitigate the nonlinearity of the FWI objective function that is generated by cycle-skipping. WRI relies on the reconstruction of data-assimilated wavefields, which approach the true wavefields near the receivers, a helpful feature when the source is added as an additional optimization variable. We present an adaptation of a recently proposed version of WRI based on the alternating direction method of multipliers (ADMM) that first finds the location of the MSEs and then reconstructs the wavefields and the source signatures jointly. Finally, the subsurface model is updated to focus the MSEs at their true location. The method does not require prior knowledge of the number of MSEs. The inversion is stabilized by sparsifying regularizations separately tailored to the source location and velocity model subproblems. The method is tested on the Marmousi model using one MSE and two clusters of MSEs with two different initial velocity models, an accurate one and a rough one, as well as with added noise. In all cases, the method accurately locates the MSEs and recovers their source signatures.

[1]  L. Eisner,et al.  Recent Advances and Challenges of Waveform‐Based Seismic Location Methods at Multiple Scales , 2020, Reviews of Geophysics.

[2]  Gregory C. Beroza,et al.  Reverse-Time Migration for Microseismic Sources Using the Geometric Mean as an Imaging Condition , 2015 .

[3]  Tingting Wu,et al.  AN OPTIMAL 9-POINT FINITE DIFFERENCE SCHEME FOR THE HELMHOLTZ EQUATION WITH PML , 2012 .

[4]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[5]  Hicks,et al.  Gauss–Newton and full Newton methods in frequency–space seismic waveform inversion , 1998 .

[6]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[7]  T. Zhu,et al.  Full-waveform inversion of passive seismic data for sources and velocities , 2016 .

[8]  Felix J. Herrmann,et al.  Mitigating local minima in full-waveform inversion by expanding the search space , 2013 .

[9]  Felix J. Herrmann,et al.  Source Estimation for Wavefield Reconstruction Inversion , 2015 .

[10]  A. Owen A robust hybrid of lasso and ridge regression , 2006 .

[11]  E. Tessmer,et al.  Reverse modelling for seismic event characterization , 2005 .

[12]  Felix J. Herrmann,et al.  Source Estimation for Frequency-domain FWI with Robust Penalties , 2012 .

[13]  Ali Gholami,et al.  Improving full-waveform inversion by wavefield reconstruction with the alternating direction method of multipliers , 2018, GEOPHYSICS.

[14]  Maksymilian Dryja,et al.  An Additive Variant of the Schwarz Alternating Method for the Case of Many Subregions , 2018 .

[15]  M. Maharramov,et al.  Total-variation minimization with bound constraints , 2015, 1505.05694.

[16]  Ali Gholami,et al.  Compound Regularization of Full-Waveform Inversion for Imaging Piecewise Media , 2019, IEEE Transactions on Geoscience and Remote Sensing.

[17]  T. Alkhalifah,et al.  Microseismic event estimation and velocity analysis based on a source-focusing function , 2019, GEOPHYSICS.

[18]  R. Pratt Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

[19]  George A. McMechan,et al.  Determination of source parameters by wavefield extrapolation , 1982 .

[20]  S. Minkoff,et al.  Microseismic event estimation in noisy data via full waveform inversion , 2015 .

[21]  Full-waveform inversion for microseismic events using sparsity constraints , 2017 .

[22]  James Rickett The variable projection method for waveform inversion with an unknown source function , 2013 .

[23]  Gary Martin,et al.  Marmousi2 An elastic upgrade for Marmousi , 2006 .

[24]  E. Haber,et al.  On optimization techniques for solving nonlinear inverse problems , 2000 .

[25]  René-Édouard Plessix,et al.  A parametrization study for surface seismic full waveform inversion in an acoustic vertical transversely isotropic medium , 2011 .

[26]  Ki Ha Lee,et al.  Source-independent full waveform inversion of seismic data , 2003 .

[27]  Michael B. Smith,et al.  Hydraulic Fracturing: History of an Enduring Technology , 2010 .

[28]  S. Operto,et al.  Full waveform inversion by proximal Newton method using adaptive regularization , 2020, Geophysical Journal International.

[29]  T. Leeuwen,et al.  A penalty method for PDE-constrained optimization in inverse problems , 2015, 1504.02249.

[30]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[31]  C. Thurber,et al.  Advances in seismic event location , 2000 .

[32]  Stephen P. Boyd,et al.  Proximal Algorithms , 2013, Found. Trends Optim..

[33]  I. Tsvankin,et al.  Gradient calculation for waveform inversion of microseismic data in VTI media 185 x 1 x 3 θ , 2014 .

[34]  S. Minkoff,et al.  A self-adjoint velocity-stress full-waveform inversion approach to microseismic source estimation , 2018, GEOPHYSICS.

[35]  Rongrong Wang,et al.  Fast sparsity-promoting microseismic source estimation , 2018, Geophysical Journal International.

[36]  N. Anders Petersson,et al.  Source Estimation by Full Wave Form Inversion , 2014, J. Sci. Comput..

[37]  J. Maddison,et al.  Full wavefield inversion of ambient seismic noise , 2018 .

[38]  Aria Abubakar,et al.  Application of the finite-difference contrast-source inversion algorithm to seismic full-waveform data , 2009 .

[39]  H S Aghamiry,et al.  Efficient extended-search space full-waveform inversion with unknown source signatures , 2021 .

[40]  S. Operto,et al.  Implementing bound constraints and total-variation regularization in extended full waveform inversion with the alternating direction method of multiplier: application to large contrast media , 2019, Geophysical Journal International.