Low-frequency seismic data are crucial for convergence of full-waveform inversion (FWI) to reliable subsurface properties. However, it is challenging to acquire field data with an appropriate signal-to-noise ratio in the low-frequency part of the spectrum. We have extrapolated low-frequency data from the respective higher frequency components of the seismic wavefield by using deep learning. Through wavenumber analysis, we find that extrapolation per shot gather has broader applicability than per-trace extrapolation. We numerically simulate marine seismic surveys for random subsurface models and train a deep convolutional neural network to derive a mapping between high and low frequencies. The trained network is then tested on sections from the BP and SEAM Phase I benchmark models. Our results indicate that we are able to recover 0.25 Hz data from the 2 to 4.5 Hz frequencies. We also determine that the extrapolated data are accurate enough for FWI application.
Building realistic and reliable models of the subsurface is the primary goal of seismic imaging. We have constructed an ensemble of convolutional neural networks (CNNs) to build velocity models directly from the data. Most other approaches attempt to map full data into 2D labels. We exploit the regularity of seismic acquisition and train CNNs to map gathers of neighboring common midpoints (CMPs) to vertical 1D velocity logs. This allows us to integrate well-log data into the inversion, simplify the mapping by using the 1D labels, and accommodate larger dips relative to using single CMP inputs. We dynamically generate the training data in parallel with training the CNNs, which reduces overfitting. Data generation and training of CNNs is more computationally expensive than conventional full-waveform inversion (FWI). However, once the network is trained, data sets with similar acquisition parameters can be inverted much faster than with FWI. The multiCMP CNN ensemble is tested on multiple realistic synthetic models, performs well, and was combined with FWI for even better performance.
23 June 2022 in IEEE Transactions on Geoscience and Remote Sensing
Low-frequency (LF) signal content in seismic data as well as a realistic initial model are key ingredients for robust and efficient full-waveform inversions (FWIs). However, acquiring LF data is challenging in practice for active seismic surveys. Data-driven solutions show promise to extrapolate LF data given a high-frequency counterpart. While being established for synthetic acoustic examples, the application of bandwidth extrapolation to field datasets remains nontrivial. Rather than aiming to reach superior accuracy in bandwidth extrapolation, we propose to jointly reconstruct LF data and a smooth background subsurface model within a multitask deep learning framework. We automatically balance data, model, and trace-wise correlation loss terms in the objective functional and show that this approach improves the extrapolation capability of the network. We also design a pipeline for generating synthetic data suitable for field data applications. Finally, we apply the same trained network to synthetic and real marine streamer datasets and run an elastic FWI from the extrapolated dataset.
Full-waveform inversion (FWI) attempts to resolve an ill-posed nonlinear optimization problem to retrieve the unknown subsurface model parameters from seismic data. In general, FWI fails to obtain an adequate representation of models with large high-velocity structures over a wide region, such as salt bodies and the sediments beneath them, in the absence of low frequencies in the recorded seismic signal, due to nonlinearity and nonuniqueness. We alleviate the ill posedness of FWI associated with data sets affected by salt bodies using model regularization. We have split the optimization problem into two parts: First, we minimize the data misfit and the total variation in the model, seeking to achieve an inverted model with sharp interfaces; and second, we minimize sharp velocity drops with depth in the model. Unlike conventional industrial salt flooding, our technique requires minimal human intervention and no information about the top of the salt. Those features are demonstrated on data sets of the BP 2004 and Sigsbee2A models, synthesized from a Ricker wavelet of dominant frequency 5.5 Hz and minimum frequency 3 Hz. We initiate the inversion process with a simple model in which the velocity increases linearly with depth. The model is well-retrieved when the same constant density acoustic code is used to simulate the observed data, which is still one of the most common FWI tests. Moreover, our technique allows us to reconstruct a reasonable depiction of the salt structure from the data synthesized independently with the BP 2004 model with variable density. In the Sigsbee2A model, we manage to even capture some of the fine layering beneath the salt. In addition, we evaluate the versatility of our method on a field data set from the Gulf of Mexico.
14 March 2019 in Journal of Geophysical Research: Solid Earth
Full‐waveform inversion (FWI) optimizes the subsurface properties of geophysical earth models in such a way that the modeled data, based on these subsurface properties, match the observed data. The anisotropic properties, whether monoclinic, orthorhombic, triclinic, or vertical transversally isotropic (VTI), of the subsurface, be it a fractured reservoir or the core‐mantle boundary, are necessary to describe the observed wave phenomena. There are no principal limitations on the complexity of the anisotropy that can be inverted using FWI. However, the question remains – what kind of anisotropic descriptions of the elastic properties of the earth can or cannot be inverted reliably from seismic waveforms? We reveal the resolution that can be achieved through reconstructions of each elastic parameter by building vertical resolution patterns from the scattering radiation patterns of body waves. A visual analysis of these patterns indicates "tradeoffs", i.e., perturbations of parameters that have the same reflection‐based scattering patterns as other perturbations. Each tradeoff leads to an apparent ambiguity in the inversion, which must be addressed by additional assumptions, constraints, or regularizations. For orthorhombic media, we find the exact tradeoffs that exist between parameters. Our parameterization isolates the VTI parameters and, therefore, we also obtain tradeoffs for every scattering mode for VTI media. We discover that only monoclinic parameters are recoverable from the first order scattering of monotypic waves. We summarize the tradeoffs found in tables for easy reference. This paper is intended to be useful for researchers setting up anisotropic FWI problems, interpreting or controlling the quality of such inversions.