Abstract:A new efficient random search algorithm is introduced for solving inversion problems in geophysical studies. The proposed algorithm is inherently a stochastic optimization method which is built on the concept of gradient descending and Lévy flights. Therefore, the algorithm is referred to as the Lévy gradient descent (L-GD). In which, the Lévy flights is a special class of random walk which consists of many short steps along with a few large steps. Such movements are observed in a varying range of fields, including animals’ foraging patterns, fluid dynamics, transport of light and so on. Meanwhile, the Lévy flights typically shows much higher speed in searching for sparsely located targets compared to the well-known Brownian walks, which make them preferable to drive random search algorithms. As shown in the paper, besides optimal solutions of the inverse problems, the L-GD algorithm could also produce estimations on the error distributions of the resultant model parameters. Following a detailed introduction of the methodology, parameter settings of the algorithm are discussed in length through statistical experiments. Subsequently, the proposed algorithm is evaluated using numeric tests and shows attracting properties of the global convergence and significant higher searching efficiency compared to commonly adopted stochastic optimization techniques in geophysical inversions. Moreover, the L-GD algorithm is applied to the inversions of gravity and seismic travel time data and has achieved the same accuracy as gradient-based optimization methods. Meanwhile, though error estimations generated by L-GD algorithm are essentially qualitative, they could still provide valuable information to help evaluating the resultant model parameters, which is of great importance for practical geophysical inversions.

Seismic Impedance Inversion Via L0 Gradient Minimisation

Elastic Impedance Inversion Based on Improved Particle Swarm Optimization

Reflected wave least squares reverse time migration with angle illumination compensation

Deep Learning-Based Low-Frequency Extrapolation and Impedance Inversion of Seismic Data

Seismic Acoustic Impedance Inversion Using Reweighted L1-Norm Sparse Constraint

Post‐stack impedance blocky inversion based on analytic solution of viscous acoustic wave equation

Lévy Gradient Descent: Augmented Random Search for Geophysical Inverse Problems

Three-parameter Prestack Seismic Inversion Based on L1-2 Minimization

Low-frequency constrained seismic impedance inversion combining large kernel attention and long short-term memory

High-resolution seismic acoustic impedance inversion with the sparsity-based statistical model

Pre-stack sparse envelope seismic inversion method adopting the L0-L2 norm regularization

Robust Q-compensated multidimensional impedance inversion using seislet-domain shaping regularization

Deep learning for high-resolution multichannel seismic impedance inversion

L1−2 Minimization for Exact and Stable Seismic Attenuation Compensation

High-resolution Seismic Impedance Inversion Integrating the Closed-Loop Convolutional Neural Network and Geostatistics: an Application to the Thin Interbedded Reservoir

Multi-trace Seismic Impedance Inversion with Structure-oriented Minimum Entropy Stabilizer

Joint Inversion of Gravity and Gravity Gradient Data Using Smoothed L0 Norm Regularization Algorithm with Sensitivity Matrix Compression

Global Optimizing Prestack Seismic Inversion Approach Using an Accurate Hessian Matrix Based on Exact Zoeppritz Equations

Regularized Seismic Velocity Inversion Based on Inverse Scattering Theory

Global optimization with deep-learning-based acceleration surrogate for large-scale seismic acoustic-impedance inversion

Geostatistical inversion method based on seismic waveform similarity