Mini-batch optimized full waveform inversion with geological constrained gradient filtering

Abstract High computation cost and generating solutions without geological sense have hindered the wide application of Full Waveform Inversion (FWI). Source encoding technique is a way to dramatically reduce the cost of FWI but subject to fix-spread acquisition setup requirement and slow convergence for the suppression of cross-talk. Traditionally, gradient regularization or preconditioning is applied to mitigate the ill-posedness. An isotropic smoothing filter applied on gradients generally gives non-geological inversion results, and could also introduce artifacts. In this work, we propose to address both the efficiency and ill-posedness of FWI by a geological constrained mini-batch gradient optimization method. The mini-batch gradient descent optimization is adopted to reduce the computation time by choosing a subset of entire shots for each iteration. By jointly applying the structure-oriented smoothing to the mini-batch gradient, the inversion converges faster and gives results with more geological meaning. Stylized Marmousi model is used to show the performance of the proposed method on realistic synthetic model.

[1]  A. Tarantola Inversion of seismic reflection data in the acoustic approximation , 1984 .

[2]  Tariq Alkhalifah,et al.  Scattering-angle based filtering of the waveform inversion gradients , 2014 .

[3]  Max Deffenbaugh,et al.  Efficient seismic forward modeling using simultaneous random sources and sparsity , 2008 .

[4]  Yi Luo,et al.  Full-traveltime inversion , 2016 .

[5]  Gao Jing Random seismic noise suppression via structure-adaptive median filter , 2012 .

[6]  Jinghuai Gao,et al.  Multimutation Differential Evolution Algorithm and Its Application to Seismic Inversion , 2016, IEEE Transactions on Geoscience and Remote Sensing.

[7]  Felix J. Herrmann,et al.  Fast waveform inversion without source‐encoding , 2013 .

[8]  Curtis C. Ober,et al.  Faster shot-record depth migrations using phase encoding , 1998 .

[9]  Thomas A. Dickens,et al.  Encoding Multiple Shot Gathers In Prestack Migration , 2000 .

[10]  Changsoo Shin,et al.  An algorithm adapting encoded simultaneous-source full-waveform inversion to marine-streamer acquisition data , 2014 .

[11]  Tariq Alkhalifah,et al.  The optimized gradient method for full waveform inversion and its spectral implementation , 2016 .

[12]  Dave Hale Structure-oriented Bilateral Filtering of Seismic Images , 2011 .

[13]  Gerard T. Schuster,et al.  Least-squares reverse time migration of marine data with frequency-selection encoding , 2013 .

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

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

[16]  Ru-Shan Wu,et al.  Seismic envelope inversion and modulation signal model , 2014 .

[17]  Mrinal K. Sen,et al.  Estimating a starting model for full-waveform inversion using a global optimization method , 2016 .

[18]  M. Sacchi,et al.  Comparison of multifrequency selection strategies for simultaneous-source full-waveform inversion , 2014 .

[19]  Taejong Kim,et al.  High-resolution model building with multistage full-waveform inversion for narrow-azimuth acquisition data , 2016 .

[20]  T. Leeuwen,et al.  Fast randomized full-waveform inversion with compressive sensing , 2012 .

[21]  Gijs C. Fehmers,et al.  Fast structural interpretation with structure-oriented filtering , 2002 .

[22]  Jean Virieux,et al.  An efficient frequency-domain full waveform inversion method using simultaneous encoded sources , 2011 .

[23]  J. Krebs,et al.  Fast full-wavefield seismic inversion using encoded sources , 2009 .

[24]  Mark W. Schmidt,et al.  Hybrid Deterministic-Stochastic Methods for Data Fitting , 2011, SIAM J. Sci. Comput..

[25]  Antoine Guitton,et al.  Constrained full-waveform inversion by model reparameterization1 , 2012 .

[26]  Louis A. Romero,et al.  Phase encoding of shot records in prestack migration , 2000 .

[27]  F. Herrmann,et al.  A new optimization approach for source-encoding full-waveform inversion , 2013 .

[28]  R. Pratt,et al.  Waveform-based simulated annealing of crosshole transmission data: A semi-global method for estimating seismic anisotropy , 2014 .