• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.227-241
• /
• 2022

Full waveform inversion (FWI) in the field of seismic data processing is an inversion technique that is used to estimate the velocity model of the subsurface for oil and gas exploration. Recently, deep learning (DL) technology has been increasingly used for seismic data processing, and its combination with FWI has attracted remarkable research efforts. For example, DL-based data processing techniques have been utilized for preprocessing input data for FWI, enabling the direct implementation of FWI through DL technology. DL-based FWI can be divided into the following methods: pure data-based, physics-based neural network, encoder-decoder, reparameterized FWI, and physics-informed neural network. In this review, we describe the theory and characteristics of the methods by systematizing them in the order of advancements. In the early days of DL-based FWI, the DL model predicted the velocity model by preparing a large training data set to adopt faithfully the basic principles of data science and apply a pure data-based prediction model. The current research trend is to supplement the shortcomings of the pure data-based approach using the loss function consisting of seismic data or physical information from the wave equation itself in deep neural networks. Based on these developments, DL-based FWI has evolved to not require a large amount of learning data, alleviating the cycle-skipping problem, which is an intrinsic limitation of FWI, and reducing computation times dramatically. The value of DL-based FWI is expected to increase continually in the processing of seismic data.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.86-92
• /
• 2006

The full-waveform inversion algorithm using normalised seismic wavefields can avoid potential inversion errors due to source estimation required in conventional full-waveform inversion methods. In this paper, we have modified the inversion scheme to install a weighted smoothness constraint for better resolution, and to implement a staged approach using normalised wavefields in order of increasing frequency instead of inverting all frequency components simultaneously. The newly developed scheme is verified by using a simple two-dimensional fault model. One of the most significant improvements is based on introducing weights in model parameters, which can be derived from integrated sensitivities. The model-parameter weighting matrix is effective in selectively relaxing the smoothness constraint and in reducing artefacts in the reconstructed image. Simultaneous multiple-frequency inversion can almost be replicated by multiple single-frequency inversions. In particular, consecutively ordered single-frequency inversion, in which lower frequencies are used first, is useful for computation efficiency.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.49-56
• /
• 1998

Seismic tomography has been widely used as high resolution subsurface imaging techniques in engineering applications. Although most of the techniques have been using travel time inversion, waveform method is being driven forward owing to the progress of computational environments. Although full-waveform inversion method has been known as the best method in terms of model resolving power without high-frequency restriction and weak scattering approximation, it has practical disadvantage that it is apt to get stuck in local minimum if the initial guess is far from the actual model and it consumes so much time to calculate. In this study, 2-D full-waveform inversion algorithm in acoustic medium is developed, which uses result of traveltime tomography as initial model. From the application on synthetic data, it is proved that this approach can efficiently reduce the problem of conventional approaches: our algorithm shows much faster convergence rate and improvement of model resolution. Result of application on physical modeling data also shows much improvement. It is expected that this algorithm can be applicable to real data.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.129-135
• /
• 2012

The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.214-216
• /
• 2022

Distributed acoustic sensing (DAS), an increasingly growing acquisition technique in the oil and gas exploration and seismology fields, has been used to record seismic signals using optical cables as receivers. With the development of imaging methods for DAS data, full waveform inversion (FWI) is been applied to DAS data to obtain high-resolution property models such as P- and S-velocity. However, because the DAS systems measure strain from the phase distortion between two points along optical cables, DAS data must be transformed from strain to particle velocity for FWI algorithms. In this study, a plane-wave FWI algorithm based on the relationship between strain and horizontal particle velocity in the plane-wave assumption is proposed to apply FWI to DAS data. Under the plane-wave assumption, strain equals the horizontal particle velocity, which is scaled by the velocity at the receiver position. This relationship was confirmed using a numerical experiment. Furthermore, 4-layer and modified Marmousi-2 velocity models were used to verify the applicability of the proposed FWI algorithm in various survey environments. The proposed FWI was implemented in land and marine survey environments and provided high-resolution P- and S-velocity models.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.203-213
• /
• 2022

To understand an underlying geological structure via seismic imaging, the velocity information of the subsurface medium is crucial. Although the full-waveform inversion (FWI) method is considered useful for estimating subsurface velocity models, 3D FWI needs a lot-of computing power and time. Herein, we compare the calculation efficiency and accuracy of frequency-domain 2D and 3D acoustic FWIs. Thereafter, we demonstrate that the artifacts from 2D approximation can be partially suppressed via frequency-domain 2D FWI by employing diffraction angle filtering (DAF). By applying DAF, which employs only big reflection angle components, the impact of noise and out-of-plane reflections can be reduced. Additionally, it is anticipated that the DAF can create long-wavelength velocity structures for 3D FWI and migration.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.242-251
• /
• 2022

Automatic differentiation automatically calculates the derivatives of a function using the chain rule once the forward operation of a function is defined. Given the recent development of computing libraries that support automatic differentiation, many researchers have adopted automatic differentiation techniques to solve geophysical inverse problems. We analyzed the advantages, disadvantages, and performances of automatic differentiation techniques using the gradient calculations of seismic full waveform inversion objective functions. The gradients of objective functions can be expressed as multiplications of the derivatives of the model parameters, wavefields, and objective functions using the chain rule. Using numerical examples, we demonstrated the speed of analytic differentiation and the convenience of complex gradient calculations for automatic differentiation. We calculated derivatives of model parameters and objective functions using automatic differentiation and derivatives of wavefields using analytic differentiation.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.33-45
• /
• 2002

It was necessary to devise new techniques to overcome and enhance the resolution limits of traveltime tomography. Waveform inversion has been one of the methods for giving very high resolution result. High resolution image could be acquired because waveform inversion used not only phase but amplitude. But waveform inversion was much time consuming Job because forward and backward modeling was needed at each iteration step. Velocity-stress method was used for effective modeling. Resolution limits of imaging methods such as travel time inversion, acoustic and elastic waveform inversion were investigated with numerical models. it was investigated that Resolution limit of waveform inversion was similar tn resolution limit of migration derived by Schuster. Horizontal resolution limit could be improved with increased coverage by adding VSP data in cross hole that had insufficient coverage. Also, waveform inversion was applied to realistic models to evaluate applicability and using initial guess of travel time tomograms to reduce non-linearity of waveform inversion showed that the better reconstructed image could be acquired.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2011.10a
• /
• pp.49-53
• /
• 2011

• Geophysics and Geophysical Exploration
• /
• v.24 no.2
• /
• pp.53-66
• /
• 2021

Velocity model building is an essential procedure in seismic data processing. Conventional techniques, such as traveltime tomography or velocity analysis take longer computational time to predict a single velocity model and the quality of the inversion results is highly dependent on human expertise. Full-waveform inversions also depend on an accurate initial model. Recently, deep neural network techniques are gaining widespread acceptance due to an increase in their integration to solving complex and nonlinear problems. This study investigated cases of seismic velocity model building using deep neural network techniques by classifying items according to the neural networks used in each study. We also included cases of generating training synthetic velocity models. Deep neural networks automatically optimize model parameters by training neural networks from large amounts of data. Thus, less human interaction is involved in the quality of the inversion results compared to that of conventional techniques and the computational cost of predicting a single velocity model after training is negligible. Additionally, unlike full-waveform inversions, the initial velocity model is not required. Several studies have demonstrated that deep neural network techniques achieve outstanding performance not only in computational cost but also in inversion results. Based on the research results, we analyzed and discussed the characteristics of deep neural network techniques for building velocity models.