6

An Efficient Implementation of Least-Squares Reverse Time Migration Based on Stable Pseudo-Acoustic Tti Wave Equation

Conventional tools for seismic imaging normally ignore the anisotropy of the media to produce images of the subsurface and such omission reduces the quality of the images. Therefore, it is necessary to develop methods which account for subsurface …

Solving one-step wave extrapolation matrix method using Krylov methods for matrix functions

We investigate the use of a restarted Krylov subspace method for computing the action of a function of large sparse matrix on a vector to numerically solve the one-step wave extrapolation matrix (OSEM) method. The OSEM method belongs to a category …

Seismic modeling from scratch using Devito: A demonstration with a typical Brazilian pre-salt model

Generating seismic modeling software is a prime activity in exploration seismology as it plays an important role in almost all its aspects such as data acquisition, processing, interpretation, and reservoir characterization. The need to write high …

Using of Fluidfft for the Optimal Choice of HPC FFT Library: A 3D Seismic Modeling Example

In this work, we illustrate a practical application of the use of Fluidfft, a newly developed common API (C++ and Python) for performing Fast Fourier Transforms (FFT) with different FFT libraries. Fluidfft is an open-source software aimed at filling …

A Comparison Between Two Elitism-Based Global Optimization Methods for Estimating a Starting Velocity Model for FWI

The performance of two global optimization methods (GOMs) with embedded elitism strategies for solving the 2D acoustic seismic inverse problem is compared here using synthetic acoustic data of the Marmousi model. A real-coded elitist genetic …

Towards automatically building starting models for full-waveform inversion using global optimization methods: A PSO approach via DEAP+ Devito

In this work, we illustrate an example in estimating the macromodel of velocities in the subsurface through the use of global optimization methods (GOMs). The optimization problem is solved using DEAP (Distributed Evolutionary Algorithms in Python) …

Lanczos Bidiagonalization Method for Parallel 3-D Gravity Inversion-Application to Basement Relief Definition

It is present an efficient parallel algorithm for the inversion of 3-D gravity data, which goal is to estimate the depth of a sedimentary basin in which the density contrast varies parabolically with depth. The efficiency of the gravity inversion …