• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.191-196
• /
• 2004

Seismic waveform inversion can be solved by using the classical Gauss-Newton method, which needs to construct the huge Hessian by the directly computed Jacobian. The property of Hessian mainly depends upon a source and receiver aperture, a velocity model, an illumination Bone and a frequency content of source wavelet. In this paper, we try to invert the Marmousi seismic data by controlling the huge Hessian appearing in the Gauss-Newton method. Wemake the two kinds of he approximate Hessian. One is the banded Hessian and the other is the approximate Hessian with automatic gain function. One is that the 1st updated velocity model from the banded Hessian is nearly the same of the result from the full approximate Hessian. The other is that the stability using the automatic gain function is more improved than that without automatic gain control.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.12a
• /
• pp.85-85
• /
• 2007

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

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

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2006.03a
• /
• pp.25-32
• /
• 2006

We present a imaging method of seismic sources by time reversal propagation of seismic waves. Time-reversal wave propagation is actively used in medical imaging, non destructive testing and waveform tomography. Time-reversal wave propagation is based on the time-reversal invariance and the spatial reciprocity of the wave equation. A signal is recorded by an array of receivers, time-reversed and then back-propagated into the medium. The time-reversed signal propagates back into the same medium and the energy refocuses back at the source location. The increasing power of computers and numerical methods makes it possible to simulate more accurately the propagation of seismic waves in heterogenous media. In this work, a staggered-grid finite-difference solution of the elastic wave equation is employed for the wave propagation simulation. With numerical experiments, we show that the time-reversal imaging will enable us to explore the spatio-temporal history of complex earthquake.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.57 no.3
• /
• pp.55-63
• /
• 2015

Embankment of agricultural reservoir started by four major rivers project. Most agricultural reservoirs must insure the agricultural water, they need must be ensured stability of embankment. Recently, there is a growing interest in seismic stability of structure by earthquake. Results of evaluation of the structural stability through seismic vibration test and numerical analysis, maximum displacement and the maximum acceleration is a similar trends. Appeared by increasing occurrence of the value of the displacement and acceleration of the structure with the result long period wave type in accordance with the seismic wave in the case of seismic waves, which shows the results of similar tendency as long period wave type consists of waveform seismic acceleration. Model test and numerical analysis results with in order to increase embankment agricultural reservoir, the displacement was found to ensure it is displayed within one percentage structural stability of the embankment.

• The Journal of the Acoustical Society of Korea
• /
• v.39 no.6
• /
• pp.592-599
• /
• 2020

This paper proposes a Convolutional Recurrent Neural Net (CRNN) structure that can simultaneously reflect both static and dynamic characteristics of seismic waveforms for various earthquake events classification. Addressing various earthquake events, including not only micro-earthquakes and artificial-earthquakes but also macro-earthquakes, requires both effective feature extraction and a classifier that can discriminate seismic waveform under noisy environment. First, we extract the static characteristics of seismic waveform through an attention-based convolution layer. Then, the extracted feature-map is sequentially injected as input to a multi-input single-output Long Short-Term Memory (LSTM) network structure to extract the dynamic characteristic for various seismic event classifications. Subsequently, we perform earthquake events classification through two fully connected layers and softmax function. Representative experimental results using domestic and foreign earthquake database show that the proposed model provides an effective structure for various earthquake events classification.

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