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.

A Hybrid Steepest Descent Method For L-Infinity Geometry Problems

Smooth Approximation of L_infinity-Norm for Multi-view Geometry

A New Hybrid Descent Algorithm for Large-Scale Nonconvex Optimization and Application to Some Image Restoration Problems

On the Convergence of A Data-Driven Regularized Stochastic Gradient Descent for Nonlinear Ill-Posed Problems

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

Convergence of Inexact Steepest Descent Algorithm for Multiobjective Optimizations on Riemannian Manifolds Without Curvature Constraints

A Non-Local Gradient Based Approach of Infinity Laplacian with Γ-convergence

A hybrid inexact regularized Newton and negative curvature method

Linear Convergence of Inexact Descent Method and Inexact Proximal Gradient Algorithms for Lower-Order Regularization Problems

O(log T) Projections for Stochastic Optimization of Smooth and Strongly Convex Functions

Descent Methods for Vector Optimization Problems: A Majorization-Minimization Perspective

L0+L1+L2 mixed optimization: a geometric approach to seismic imaging and inversion using concepts in topology and semigroup

Accelerating Globally Optimal Consensus Maximization in Geometric Vision

Optimization Algorithms for Faster Computational Geometry.

A class of new three-term descent conjugate gradient algorithms for large-scale unconstrained optimization and applications to image restoration problems

Universal Gradient Descent Ascent Method for Nonconvex-Nonconcave Minimax Optimization

An Analysis of a Multi-Level Projected Steepest Descent Iteration for Nonlinear Inverse Problems in Banach Spaces Subject to Stability Constraints

Faster Acceleration for Steepest Descent

The Role of Level-Set Geometry on the Performance of PDHG for Conic Linear Optimization

The Stochastic Steepest Descent Method for Robust Optimization in Banach Spaces

Online Learning Guided Curvature Approximation: A Quasi-Newton Method with Global Non-Asymptotic Superlinear Convergence