ALBERT

Beschreibung: In most full-waveform inversion (FWI) problems, sufficient prior information is available to constrain the velocity of certain parts of the model, e.g., the water column or, in some cases, near-surface velocities. We take advantage of this situation and develop a fast Schur-complement-based forward modeling and inversion approach by partitioning the velocity model into two parts. The first part consists of the constrained zone that does not change during the inversion, whereas the second part is the anomalous zone to be updated during the inversion. For this decomposition, we partially factorize the governing system of linear equations by computing a Schur complement for the anomalous zone. The Schur complement system is then solved to compute the fields in the anomalous zone, which are then back substituted to compute the fields in the constrained region. For each successive modeling steps with new anomalous zone velocities, the corresponding Schur complement is easily computed using simple algebra. Because the anomalous part of the model is comparatively smaller than the whole model, considerable computational savings can be achieved using our Schur approach. Additionally, we showed that the Schur complement method maintains the accuracy of standard frequency-domain finite difference formulations, but this comes at a slightly higher peak memory requirement. Our FWI workflow shows reduced runtime by 15%–57% depending upon the depth of the water column without losing any accuracy compared to the standard method.

Beschreibung: A good starting model is imperative in full-waveform inversion (FWI) because it solves a least-squares inversion problem using a local gradient-based optimization method. A suboptimal starting model can result in cycle skipping leading to poor convergence and incorrect estimation of subsurface properties. This problem is especially crucial for salt models because the strong velocity contrasts create substantial time shifts in the modeled seismogram. Incorrect estimation of salt bodies leads to velocity inaccuracies in the sediments because the least-squares gradient aims to reduce traveltime differences without considering the sharp velocity jump between sediments and salt. We have developed a technique to estimate velocity models containing salt bodies using a combination of global and local optimization techniques. To stabilize the global optimization algorithm and keep it computationally tractable, we reduce the number of model parameters by using sparse parameterization formulations. The sparse formulation represents sediments using a set of interfaces and velocities across them, whereas a set of ellipses represents the salt body. We use very fast simulated annealing (VFSA) to minimize the misfit between the observed and synthetic data and estimate an optimal model in the sparsely parameterized space. The VFSA inverted model is then used as a starting model in FWI in which the sediments and salt body are updated in the least-squares sense. We partition model updates into sediment and salt updates in which the sediments are updated like conventional FWI, whereas the shape of the salt is updated by taking the zero crossing of an evolving level set surface. Our algorithm is tested on two 2D synthetic salt models, namely, the Sigsbee 2A model and a modified SEG Advanced Modeling Program (SEAM) Phase I model while fixing the top of the salt. We determine the efficiency of the VFSA inversion and imaging improvements from the level set FWI approach and evaluate a few sources of uncertainty in the estimation of salt shapes.

Beschreibung: Migration techniques are an integral part of seismic imaging workflows. Least-squares reverse time migration (LSRTM) overcomes some of the shortcomings of conventional migration algorithms by compensating for illumination and removing sampling artifacts to increase spatial resolution. However, the computational cost associated with iterative LSRTM is high and convergence can be slow in complex media. We implement pre-stack LSRTM in a deep learning framework and adopt strategies from the data science domain to accelerate convergence. The proposed hybrid framework leverages the existing physics-based models and machine learning optimizers to achieve better and cheaper solutions. Using a time-domain formulation, we show that mini-batch gradients can reduce the computation cost by using a subset of total shots for each iteration. Mini-batch approach does not only reduce source cross-talk but also is less memory intensive. Combining mini-batch gradients with deep learning optimizers and loss functions can improve the efficiency of LSRTM. Deep learning optimizers such as the adaptive moment estimation are generally well suited for noisy and sparse data. We compare different optimizers and demonstrate their efficacy in mitigating migration artifacts. To accelerate the inversion, we adopt the regularised Huber loss function in conjunction. We apply these techniques to 2D Marmousi and 3D SEG/EAGE salt models and show improvements over conventional LSRTM baselines. The proposed approach achieves higher spatial resolution in less computation time measured by various qualitative and quantitative evaluation metrics.