• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.4
• /
• pp.221-230
• /
• 2021

This paper proposes full waveform inversion (FWI) for estimating the physical properties of a layered half-space. An FWI solution is obtained using a genetic algorithm (GA), which is a well-known global optimization approach. The dynamic responses of a layered half-space subjected to a harmonic vertical disk load are measured and compared with those calculated using the estimated physical properties. The responses are calculated using the thin-layer method, which is accurate and efficient for layered media. Subsequently, a numerical model is constructed for a layered half-space using mid-point integrated finite elements and perfectly matched discrete layers. An objective function of the global optimization problem is defined as the L2-norm of the difference between the observed and estimated responses. A GA is used to minimize the objective function and obtain a solution for the FWI. The accuracy of the proposed approach is applied to various problems involving layered half-spaces. The results verify that the proposed FWI based on a GA is suitable for estimating the material properties of a layered half-space, even when the measured responses include measurement noise.

• 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.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.2
• /
• pp.5-14
• /
• 2024

Reinforced concrete bridge decks are the first to be damaged by vehicle loads and rain infiltration. Concrete deterioration primarily occurs owing to the corrosion of rebars and other metal components by chlorides used for snow and ice melting. The structural condition and concrete deterioration of the bridge decks within the pavement were evaluated using ground-penetrating radar (GPR) survey data. To evaluate concrete deterioration in bridges, it is necessary to develop GPR data analysis techniques to accurately identify deteriorated locations and rebar positions. GPR exploration involves the acquisition of reflection and diffraction wave signals due to differences in radar wave propagation velocity in geotechnical media. Therefore, a full-waveform inversion (FWI) method was developed to evaluate the deterioration of reinforced concrete bridge decks by estimating the radar wave propagation velocity in geotechnical media using GPR data. Numerical experiments using a GPR velocity model confirmed the deterioration phenomena of bridge decks, such as concrete delamination and rebar corrosion, verifying the applicability of the developed technology. Moreover, using the synthetic GPR data, FWI facilitates the determination of rebar positions and concrete deterioration locations using inverted velocity images.

• 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.23 no.1
• /
• pp.38-49
• /
• 2020

Full-waveform inversion (FWI) is an optimization process of fitting observed and modeled data to reconstruct high-resolution subsurface physical models. In acoustic FWI (AFWI), pressure data acquired using a marine streamer has mainly been used to reconstruct the subsurface P-wave velocity models. With recent advances in marine seismic-acquisition techniques, acquiring multi-component data in marine environments have become increasingly common. Thus, AFWI strategies must be developed to effectively use marine multi-component data. Herein, we proposed an AFWI strategy using horizontal and vertical particle-acceleration data. By analyzing the modeled acoustic data and conducting sensitivity kernel analysis, we first investigated the characteristics of each data component using AFWI. Common-shot gathers show that direct, diving, and reflection waves appearing in the pressure data are separated in each component of the particle-acceleration data. Sensitivity kernel analyses show that the horizontal particle-acceleration wavefields typically contribute to the recovery of the long-wavelength structures in the shallow part of the model, and the vertical particle-acceleration wavefields are generally required to reconstruct long- and short-wavelength structures in the deep parts and over the whole area of a given model. Finally, we present a sequential-inversion strategy for using the particle-acceleration wavefields. We believe that this approach can be used to reconstruct a reasonable P-wave velocity model, even when the pressure data is not available.

• Journal of the Korean earth science society
• /
• v.37 no.5
• /
• pp.277-285
• /
• 2016

We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.3
• /
• pp.203-211
• /
• 2023

A Markov chain Monte Carlo (MCMC) simulation is proposed for probabilistic full waveform inversion (FWI) in a layered half-space. Dynamic responses on the half-space surface are estimated using the thin-layer method when a harmonic vertical force is applied. Subsequently, a posterior probability distribution function and the corresponding objective function are formulated to minimize the difference between estimations and observed data as well as that of model parameters from prior information. Based on the gradient of the objective function, a proposal distribution and an acceptance probability for MCMC samples are proposed. The proposed MCMC simulation is applied to several layered half-space examples. It is demonstrated that the proposed MCMC simulation for probabilistic FWI can estimate probabilistic material properties such as the shear-wave velocities of a layered half-space.