Advanced BEM-based methodologies to identify and simulate wave fields in complex geostructures

To enhance the applicability of BEM for geomechanical modeling numerically optimized BEM models, hybrid FEM-BEM models, and parallel computation of seismic Full Waveform Inversion (FWI) in GPU are implemented. Inverse modeling of seismic wave propagation in inhomogeneous and heterogeneous half-plane is implemented in Boundary Element Method (BEM) using Particle Swarm Optimization (PSO). The Boundary Integral Equations (BIE) based on the fundamental solutions for homogeneous elastic isotropic continuum are modified by introducing mesh-dependent variables. The variables are optimized to obtain the site-specific impedance functions. The PSO-optimized BEM models have significantly improved the efficiency of BEM for seismic wave propagation in arbitrarily inhomogeneous and heterogeneous media. Similarly, a hybrid BEM-FEM approach is developed to evaluate the seismic response of a complex poroelastic soil region containing underground structures. The far-field semi-infinite geological region is modeled via BEM, while the near-field finite geological region is modeled via FEM. The BEM region is integrated into the global FEM system using an equivalent macro-finite-element. The model describes the entire wave path from the seismic source to the local site in a single hybrid model. Additionally, the computational efficiency of time domain FWI algorithm is enhanced by parallel computation in CPU and GPU

Seismic modeling and inversion using half-precision floating-point numbers

New processors are increasingly supporting half-precision floating-point numbers, often with a significant throughput gain over single-precision operations. Seismic modeling, imaging, and inversion could benefit from such an acceleration, but it is not obvious how the accuracy of the solution can be preserved with a very narrow 16-bit representation. By scaling the finite-difference expression of the isotropic elastic wave equation, we have found that a stable solution can be obtained despite the very narrow dynamic range of the half-precision format.We develop an implementation with the CUDA platform, which, on most recent graphics processing units (GPU), is nearly twice as fast and uses half the memory of the equivalent single-precision version. The error on seismograms caused by the reduced precision is shown to correspond to a negligible fraction of the total seismic energy and is mostly incoherent with seismic phases. Finally, we find that this noise does not adversely impact full-waveform inversion nor reverse time migration, which both benefit from the higher throughput of half-precision computation

Seismic velocity estimation: A deep recurrent neural-network approach

ABSTRACT: Applying deep learning to 3D velocity model building remains a challenge due to the sheer volume of data required to train large-scale artificial neural networks. Moreover, little is known about what types of network architectures are appropriate for such a complex task. To ease the development of a deep-learning approach for seismic velocity estimation, we have evaluated a simplified surrogate problem — the estimation of the root-mean-square (rms) and interval velocity in time from common-midpoint gathers — for 1D layered velocity models. We have developed a deep neural network, whose design was inspired by the information flow found in semblance analysis. The network replaces semblance estimation by a representation built with a deep convolutional neural network, and then it performs velocity estimation automatically with recurrent neural networks. The network is trained with synthetic data to identify primary reflection events, rms velocity, and interval velocity. For a synthetic test set containing 1D layered models, we find that rms and interval velocity are accurately estimated, with an error of less than 44  m/s for the rms velocity. We apply the neural network to a real 2D marine survey and obtain accurate rms velocity predictions leading to a coherent stacked section, in addition to an estimation of the interval velocity that reproduces the main structures in the stacked section. Our results provide strong evidence that neural networks can estimate velocity from seismic data and that good performance can be achieved on real data even if the training is based on synthetics. The findings for the 1D problem suggest that deep convolutional encoders and recurrent neural networks are promising components of more complex networks that can perform 2D and 3D velocity model building