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.

The discrete analogue of the method of quickest descent for an inverse acoustic problem in case of a smooth source

The Discrete Adjoint Method: Efficient Derivatives for Functions of Discrete Sequences

Gradient-type Approaches to Inverse and Ill-Posed Problems of Mathematical Physics

A geometric integration approach to smooth optimisation: Foundations of the discrete gradient method

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

Numerical Analysis for a Hyperbolic PDE-Constrained Optimization Problem in Acoustic Full Waveform Inversion

Accelerated optimization algorithms and ordinary differential equations: the convex non Euclidean case

On the equivalence between SOR-type methods for linear systems and discrete gradient methods for gradient systems

Convergence Rate Analysis of Continuous- and Discrete-Time Smoothing Gradient Algorithms

An Optimization Approach to Multi-Dimensional Time Domain Acoustic Inverse Problems

A discrete adjoint gradient approach for equality and inequality constraints in dynamics

First Order Methods beyond Convexity and Lipschitz Gradient Continuity with Applications to Quadratic Inverse Problems

Computational enhancement of discrete gradient method

On Some Versions of Subspace Optimization Methods with Inexact Gradient Information

An Optimization Approach to a Three-Dimensional Acoustic Inverse Problem in the Time Domain

Automatically Differentiable Higher-Order Parabolic Equation for Real-Time Underwater Sound Speed Profile Sensing

Direct Runge-Kutta Discretization Achieves Acceleration.

Higher order discrete gradient method by the discontinuous Galerkin time-stepping method

Acoustic shape optimization using energy stable curvilinear finite differences

A Homotopy Method for the Inversion of a Two-Dimensional Acoustic Wave Equation

Faster Gradient Descent and the Efficient Recovery of Images