A new iterative solver for the time-harmonic wave equation

The time-harmonic wave equation, also known as the Helmholtz equation, is obtained if the constant-density acoustic wave equation is transformed from the time domain to the frequency domain. Its discretization results in a large, sparse, linear system of equations. In two dimensions, this system can be solved efficiently by a direct method. In...

Resistivity imaging with controlled-source electromagnetic data: Depth and data weighting

We discuss some computational aspects of resistivity imaging by inversion of offshore controlled-source electromagnetic data. We adopt the classic approach to imaging by formulating it as an inverse problem. A weighted least-squares functional measures the misfit between synthetic and observed data. Its minimization by a quasi-Newton algorithm...

A fast 3-D free-surface topography method for acoustic full-waveform inversion

We propose a finite-difference scheme for the simulation of seismic waves interacting with 3-D freesurface topography. The intended application is velocity model building by acoustic full-waveform inversion (FWI). The scheme follows an immersed boundary approach for wave equations in the firstorder stress-velocity formulation, discretized on a...

Acoustic VTI Full-waveform Inversion with 3-D Free-surface Topography

In land applications, topography variation may impact the imaging if not taken into account. With low-frequency and wide-aperture data, the long-to-intermediate wavelength components of the velocity model can be recovered by full-waveform inversion. Standard static corrections to handle the topography do not work satisfactorily on long-offset...