Convergent scattering series solution of the inhomogeneous Helmholtz equation via renormalization group and homotopy continuation approaches

Abstract The inhomogeneous Helmholtz equation can be represented by an equivalent integral equation of the Lippmann-Schwinger type. Since it is very costly to solve the Lippmann-Schwinger equation exactly via matrix inversion, researchers often try to use the popular Born series, which represents a physics based iterative solution. However, the popular Born series is only guaranteed to converge if the contrast volume is sufficiently small. We here use renormalization group (RG) and homotopy continuation (HC) approaches to derive novel scattering series solutions of the Lippmann-Schwinger equation which are guaranteed to converge for arbitrary large contrast volumes (perturbations). Our RG approach is based on the use of an auxillary set of scale-dependent scattering potentials and Green's operators which gradually evolve toward the real physical scattering potential and the physical Green's operator, respectively. We show that the auxillary Green's operators satisfies the renormalization group law and we derive the corresponding renormalization group flow equation by using familiar “scaling” arguments. The HC approach comes from pure mathematics and is based on the deformation of a relatively simple reference equation (initial condition) into a more complicated target equation (the Lippmann-Schwinger equation for the physical scattering potential), via a parameter that varies continuously between 0 and 1. Our HC approach is more general than the RG approach since it involves a convergence control operator, but the two approaches are nevertheless closely connected. We show that the RG flow and HC equations can be implemented very efficiently using the Fast Fourier Transform (FFT) algorithm. Also, we present a very efficient iterative scheme based on a combination of the zeroth-order deformation equation with a floating initial. We illustrate our results with a numerical example related with seismic wavefield modeling in strongly scattering media.

[1]  D. Kouri,et al.  Inverse scattering theory: renormalization of the Lippmann-Schwinger equation for acoustic scattering in one dimension. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  global sci.  Target-Oriented Inversion of Time-Lapse Seismic Waveform Data , 2020 .

[3]  J. Carcione,et al.  Seismic modelingSeismic modeling , 2002 .

[4]  William H. Butler,et al.  Multiple scattering in solids , 1999, Graduate texts in comtemporary physics.

[5]  Renormalization in few-body nuclear physics , 2002 .

[6]  Accelerating the T‐matrix approach to seismic full‐waveform inversion by domain decomposition , 2018, Geophysical Prospecting.

[7]  A. Kuvshinov,et al.  Electromagnetic field scattering in a heterogeneous Earth: A solution to the forward problem , 1995 .

[8]  Mark Spencer Haynes,et al.  Full-wave Nonlinear Inverse Scattering for Acoustic and Electromagnetic Breast Imaging. , 2012 .

[9]  Ana B. Ramirez,et al.  Sparsity constrained contrast source inversion. , 2016, The Journal of the Acoustical Society of America.

[10]  Shijun Liao,et al.  Advances in the homotopy analysis method , 2014 .

[11]  Weng Cho Chew,et al.  An iterative solution of the two‐dimensional electromagnetic inverse scattering problem , 1989, Int. J. Imaging Syst. Technol..

[12]  M. Gell-Mann,et al.  QUANTUM ELECTRODYNAMICS AT SMALL DISTANCES , 1954 .

[13]  Bjørn Ursin,et al.  Full waveform inversion in the frequency domain using direct iterative T-matrix methods , 2015 .

[14]  Mrinal K. Sen,et al.  Frequency-domain full waveform inversion with a scattering-integral approach and its sensitivity analysis , 2013 .

[15]  Ru-Shan Wu,et al.  One-way and one-return approximations (de wolf approximation) for fast elastic wave modeling in complex media , 2007 .

[16]  Arthur B. Weglein,et al.  Inverse scattering series and seismic exploration , 2003 .

[17]  D. P. Datta,et al.  Comparative Study of Homotopy Analysis and Renormalization Group Methods on Rayleigh and Van der Pol Equations , 2014, 1405.6615.

[18]  L. T. Watson,et al.  Globally convergent homotopy methods: a tutorial , 1989, Conference on Numerical Ordinary Differential Equations.

[19]  Ru-Shan Wu,et al.  Non-linear partial derivative and its De Wolf approximation for non-linear seismic inversion , 2014 .

[20]  T. Jordan,et al.  FAST TRACK PAPER: Full three-dimensional tomography: a comparison between the scattering-integral and adjoint-wavefield methods , 2007 .

[21]  Weng Cho Chew,et al.  Integral Equation Methods for Electromagnetic and Elastic Waves , 2007, Synthesis Lectures on Computational Electromagnetics.

[22]  Renormalized scattering series for frequency-domain waveform modelling of strong velocity contrasts , 2016 .

[23]  Ru-Shan Wu,et al.  Seismic envelope inversion and renormalization group theory: Nonlinear scale separation and slow dynamics , 2016 .

[24]  Edward Roy Pike,et al.  Scattering: Scattering and Inverse Scattering in Pure and Applied Science , 2002 .

[25]  Eleftherios Kirkinis,et al.  The Renormalization Group: A Perturbation Method for the Graduate Curriculum , 2012, SIAM Rev..

[26]  L. Ryder,et al.  Quantum Field Theory , 2001, Foundations of Modern Physics.

[27]  Kristopher A. Innanen,et al.  Born series forward modelling of seismic primary and multiple reflections: an inverse scattering shortcut , 2009 .

[28]  D. A. de Wolf.  Correction to "Renormalization of EM fields in application to large-angle scattering from randomly continuous media and sparse particle distributions" , 1985 .

[29]  Kjersti Solberg Eikrem,et al.  Iterated extended Kalman filter method for time‐lapse seismic full‐waveform inversion , 2018, Geophysical Prospecting.

[30]  P. M. Berg,et al.  A contrast source inversion method , 1997 .

[31]  R. Kress,et al.  Integral equation methods in scattering theory , 1983 .

[32]  H. Munthe-Kaas.  Runge-Kutta methods on Lie groups , 1998 .

[33]  J. Yorke,et al.  The homotopy continuation method: numerically implementable topological procedures , 1978 .

[34]  D. A. de Wolf,et al.  Electromagnetic reflection from an extended turbulent medium: Cumulative forward-scatter single-backscatter approximation , 1971 .

[35]  B. Han,et al.  A homotopy method for the inversion of a two-dimensional acoustic wave equation , 2005 .

[36]  E. Kirkinis.  Renormalization group interpretation of the Born and Rytov approximations. , 2008, Journal of the Optical Society of America. A, Optics, image science, and vision.

[37]  Aria Abubakar,et al.  Three-dimensional visco-acoustic modeling using a renormalized integral equation iterative solver , 2013, J. Comput. Phys..

[38]  W. Chew,et al.  Reconstruction of two-dimensional permittivity distribution using the distorted Born iterative method. , 1990, IEEE transactions on medical imaging.

[39]  Oono,et al.  Renormalization group and singular perturbations: Multiple scales, boundary layers, and reductive perturbation theory. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[40]  Morten Jakobsen,et al.  T-matrix approach to seismic forward modelling in the acoustic approximation , 2012, Studia Geophysica et Geodaetica.

[41]  Gerwin Osnabrugge,et al.  A convergent Born series for solving the inhomogeneous Helmholtz equation in arbitrarily large media , 2016, J. Comput. Phys..

[42]  P. Morse,et al.  Methods of theoretical physics , 1955 .

[43]  P. M. Berg,et al.  Towards high resolution quantitative subsurface models by full waveform inversion , 2013 .

[44]  Bertrand Delamotte,et al.  A Hint of renormalization , 2002, hep-th/0212049.

[46]  K. Wilson.  Renormalization Group and Critical Phenomena. I. Renormalization Group and the Kadanoff Scaling Picture , 1971 .