• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.05a
• /
• pp.67-71
• /
• 2009

The high resolution seismic data processing essentially requires the latest data processing technology with a system which enables to calculate the large volume of input/output and computing of data. Nowadays, there is a high foreign dependency on the data processing and QC. In this research, The test processing had been carried out to deduce the best suited software by evaluating functionality, usability and user environment. And several attempts are made to select the most appropriated specification of hardware in order to maximize the software efficiency. The results showed that more numbers of multiple and noise were eliminated in existing prestack time migration sections, which effectively enhanced the resolution.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.163-169
• /
• 2005

Recently, the velocity stack domain is having an attention as a very useful domain for various processing in seismic data processing. In order to be used in many applications, the velocity stack should be obtained through an inversion method and the used inversion should have properties like the robustness to noise and the parsimony of velocity stack result. Iteratively Reweighted Least-Squares (IRLS) method is the one of the inversion methods that have such properties. This paper describes the theoretical background, implementation of the method, and examines the characteristics and limits of the IRLS method.

• Geophysics and Geophysical Exploration
• /
• v.20 no.4
• /
• pp.232-240
• /
• 2017

Normal moveout correction is one of the main procedures of seismic reflection data processing and a crucial pre-processing step for AVO analysis. Unfortunately, stretch phenomenon, which is the intrinsic problem of NMO correction, degrades the quality of stack section and reliability of AVO analysis. Although muting is applied to resolve this problem, it makes far-offset traces more useful to develop an advanced NMO correction technique without stretch. In this paper, easy and detailed explanations are provided on the definition and methodology of NMO correction, and then the cause of stretch is explained with its characteristics. A graphical explanation for NMO correction is given for the intuitive understanding of stretch phenomenon. Additionally, the theoretical formulation is derived to quantitatively understand the NMO correction. Through explaining the muting process to remove NMO stretch, the limitations of conventional methods are investigated and the need for a new resolution comes to discussion. We describe a stretch-free NMO correction based on inverse theory among many different stretch-free NMO corrections. Finally, the stretch-free NMO correction is verified through synthetic example and real data.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.170-176
• /
• 2005

Recently the velocity stack inversion is having many attentions as an useful way to perform various seismic data processing. In order to be used in various seismic data processing, the inversion method used should have properties such as robustness to noise and parsimony of the velocity stack result. The IRLS (Iteratively Reweighted Least-Squares) method that minimizes ${L_1}-norm$ is the one used mostly. This paper introduce another method, CGG (Conjugate Guided Gradient) method, which can be used to achieve the same goal as the IRLS method does. The CGG method is a modified CG (Conjugate Gradient) method that minimizes ${L_1}-norm$. This paper explains the CGG method and compares the result of it with the one of IRSL methods. Testing on synthetic and real data demonstrates that CGG method can be used as an inversion method f3r minimizing various residual/model norms like IRLS methods.

• Economic and Environmental Geology
• /
• v.53 no.4
• /
• pp.479-489
• /
• 2020

A fault is a geological structure that can be a migration path or a cap rock of hydrocarbon such as oil and gas, formed from source rock. The fault is one of the main targets of seismic exploration to find reservoirs in which hydrocarbon have accumulated. However, conventional fault detection methods using lateral discontinuity in seismic data such as semblance, coherence, variance, gradient magnitude and fault likelihood, have problem that professional interpreters have to invest lots of time and computational costs. Therefore, many researchers are conducting various studies to save computational costs and time for fault interpretation, and machine learning technologies attracted attention recently. Among various machine learning technologies, many researchers are conducting fault interpretation studies using the support vector machine, multi-layer perceptron, deep neural networks and convolutional neural networks algorithms. Especially, researchers use not only their own convolution networks but also proven networks in image processing to predict fault locations and fault information such as strike and dip. In this paper, by investigating and analyzing these studies, we found that the convolutional neural networks based on the U-Net from image processing is the most effective one for fault detection and interpretation. Further studies can expect better results from fault detection and interpretation using the convolutional neural networks along with transfer learning and data augmentation.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.73-85
• /
• 2006

Array observations of the vertical component of microtremors are frequently conducted to estimate a subsurface layered-earth structure on the assumption that microtremors consist predominantly of the fundamental mode Rayleigh waves. As a useful tool in the data collection, processing and analysis, the spatial autocorrelation (SPAC) method is widely used, which in practice requires a circle array consisting of M circumferential stations and one centre station (called "M-station circle array", where M is the number of stations). The present paper considers the minimum number of stations required for a circle array for efficient data collection in terms of analytical efficacy and field effort. This study first rearranges the theoretical background of the SPAC algorithm, in which the SPAC coefficient for a circle array with M infinite is solely expressed as the Bessel function, $J_0(rk)$ (r is the radius and k the wavenumber). Secondly, the SPAC coefficient including error terms independent of the microtremor energy field for an M-station circle array is analytically derived within a constraint for the wave direction across the array, and is numerically evaluated in respect of these error terms. The main results of the evaluation are: 1) that the 3-station circle array when compared with other 4-, 5-, and 9-station arrays is the most efficient and favourable for observation of microtremors if the SPAC coefficients are used up to a frequency at which the coefficient takes the first minimum value, and 2) that the Nyquist wavenumber is the most influential factor that determines the upper limit of the frequency range up to which the valid SPAC coefficient can be estimated.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.131-135
• /
• 2016

Coarray is a parallel processing technique introduced in the Fortran 2008 standard. Coarray can implement parallel processing using simple syntax. In this research, we examined applicability of Coarray to seismic parallel processing by comparing performance of seismic data processing programs using Coarray and MPI. We compared calculation time using seismic wave propagation modeling and one to one communication time using domain decomposition technique. We also compared performance of parallel reverse-time migration programs using Coarray and MPI. Test results show that the computing speed of Coarray method is similar to that of MPI. On the other hand, MPI has superior communication speed to that of Coarray.

• Geophysics and Geophysical Exploration
• /
• v.5 no.3
• /
• pp.185-192
• /
• 2002

We made a study of 3-D migration which could precisely image data of GPR (Ground Penetrating Radar) applied to NDT (Non-Destructive Test) field for the inspection of structural safety. In this study, we obtained 3-D migrated images of important targets in structuresurvey (e.g. steel pipes, cracks) by using 3-D Kirchhoff prestack depth migration scheme developed for seismic data processing. For a concrete model consisting of steel pipe and void, the targets have been well defined with opposite amplitude according to the parameters of the targets. And migrated images using Parallel-Broadside array (XX configuration) have shown higher resolution than those using Perpendicular-Broadside array (YY configuration) when steel pipes had different sizes. Therefore, it is required to analyze the migrated image of XX configuration as well as that of general YY configuration in order to get more accurate information. As the last stage, we chose a model including two steel pipes which cross each other. The upper pipe has been resolved clearly but the lower has been imaged bigger than the model size due to the high conductivity of the upper steel.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.129-136
• /
• 2017

3D seismic data processing methods such as full waveform inversion or reverse-time migration require 3D wave propagation modeling and heavy calculations. We compared efficiency and accuracy of a Xeon Phi coprocessor to those of a high-end server CPU using 3D frequency-domain wave propagation modeling. We adopted the OpenMP parallel programming to the time-domain finite difference algorithm by considering the characteristics of the Xeon Phi coprocessors. We applied the Fourier transform using a running-integration to obtain the frequency-domain wavefield. A numerical test on frequency-domain wavefield modeling was performed using the 3D SEG/EAGE salt velocity model. Consequently, we could obtain an accurate frequency-domain wavefield and attain a 1.44x speedup using the Xeon Phi coprocessor compared to the CPU.

• The Journal of the Acoustical Society of Korea
• /
• v.23 no.1
• /
• pp.8-16
• /
• 2004

An inversion method is presented for the determination of the compressional wave speed, compressional wave attenuation, thickness of the sediment layer and density as a function of depth for a horizontally stratified ocean bottom. An experiment for estimating those properties was conducted in the shallow water of South Sea in Korea. In the experiment, a light bulb implosion and the propagating sound were measured using a VLA (vertical line array). As a method for estimating the geoacoustic properties, a coherent broadband matched field processing combined with Genetic Algorithm was employed. When a time-dependent signal is very short, the Fourier transform results are not accurate, since the frequency components are not locatable in time and the windowed Fourier transform is limited by the length of the window. However, it is possible to do this using the wavelet transform a transform that yields a time-frequency representation of a signal. In this study, this transform is used to identify and extract the acoustic components from multipath time series. The inversion is formulated as an optimization problem which maximizes the cost function defined as a normalized correlation between the measured and modeled signals in the wavelet transform coefficient vector. The experiments and procedures for deploying the light bulbs and the coherent broadband inversion method are described, and the estimated geoacoustic profile in the vicinity of the VLA site is presented.