• Journal of the Korean Society for Nondestructive Testing
• /
• v.25 no.1
• /
• pp.9-14
• /
• 2005

The angular dependence(profile) of backward radiated ultrasound was measured for glass specimens with random surface roughness using ultrasonic goniometer that ran changes the incident angle continuously. It was concluded that the roughened region had greater acoustic impedance than the unperturbed region. The comparison of backward radiations showed that the amplitude of peak and the area of radiation profile were increased with surface roughness. It was suggested from the sensitive dependence of the profile area that the profile of backward radiation could be applied to in the nondestructive evaluation of sulfate region. Inclined C-scan technique with the transducer inclined at Rayleigh angle showed the reverse of luminosity and the high signal to noise ratio so that it provided high resolution.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.26-32
• /
• 1999

Among the various seismic data processing sequences, the velocity analysis is the most time consuming and man-hour intensive processing steps. For the production seismic data processing, a good velocity analysis tool as well as the high performance computer is required. The tool must give fast and accurate velocity analysis. There are two different approches in the velocity analysis, batch and interactive. In the batch processing, a velocity plot is made at every analysis point. Generally, the plot consisted of a semblance contour, super gather, and a stack pannel. The interpreter chooses the velocity function by analyzing the velocity plot. The technique is highly dependent on the interpreters skill and requires human efforts. As the high speed graphic workstations are becoming more popular, various interactive velocity analysis programs are developed. Although, the programs enabled faster picking of the velocity nodes using mouse, the main improvement of these programs is simply the replacement of the paper plot by the graphic screen. The velocity spectrum is highly sensitive to the presence of the noise, especially the coherent noise often found in the shallow region of the marine seismic data. For the accurate velocity analysis, these noise must be removed before the spectrum is computed. Also, the velocity analysis must be carried out by carefully choosing the location of the analysis point and accuarate computation of the spectrum. The analyzed velocity function must be verified by the mute and stack, and the sequence must be repeated most time. Therefore an iterative, interactive, and unified velocity analysis tool is highly required. An interactive velocity analysis program, xva(X-Window based Velocity Analysis) was invented. The program handles all processes required in the velocity analysis such as composing the super gather, computing the velocity spectrum, NMO correction, mute, and stack. Most of the parameter changes give the final stack via a few mouse clicks thereby enabling the iterative and interactive processing. A simple trace indexing scheme is introduced and a program to nike the index of the Geobit seismic disk file was invented. The index is used to reference the original input, i.e., CDP sort, directly A transformation techinique of the mute function between the T-X domain and NMOC domain is introduced and adopted to the program. The result of the transform is simliar to the remove-NMO technique in suppressing the shallow noise such as direct wave and refracted wave. However, it has two improvements, i.e., no interpolation error and very high speed computing time. By the introduction of the technique, the mute times can be easily designed from the NMOC domain and applied to the super gather in the T-X domain, thereby producing more accurate velocity spectrum interactively. The xva program consists of 28 files, 12,029 lines, 34,990 words and 304,073 characters. The program references Geobit utility libraries and can be installed under Geobit preinstalled environment. The program runs on X-Window/Motif environment. The program menu is designed according to the Motif style guide. A brief usage of the program has been discussed. The program allows fast and accurate seismic velocity analysis, which is necessary computing the AVO (Amplitude Versus Offset) based DHI (Direct Hydrocarn Indicator), and making the high quality seismic sections.

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

EZTOMO is a crosswell seismic tomography software system. The system has capability of event picking, raytracing, inversion, error analysis, and visualization of the processing results. Waveform of the first arrival signal has been utilized to select the event of the first motion, and uncertainty measured in estimation of the first breaks has been utilized to improve the inversion process.

• Geophysics and Geophysical Exploration
• /
• v.26 no.2
• /
• pp.37-51
• /
• 2023

The ground roll is the most common coherent noise in land seismic data and has an amplitude much larger than the reflection event we usually want to obtain. Therefore, ground roll suppression is a crucial step in seismic data processing. Several techniques, such as f-k filtering and curvelet transform, have been developed to suppress the ground roll. However, the existing methods still require improvements in suppression performance and efficiency. Various studies on the suppression of ground roll in seismic data have recently been conducted using deep learning methods developed for image processing. In this paper, we introduce three models (DnCNN (De-noiseCNN), pix2pix, and CycleGAN), based on convolutional neural network (CNN) or conditional generative adversarial network (cGAN), for ground roll suppression and explain them in detail through numerical examples. Common shot gathers from the same field were divided into training and test datasets to compare the algorithms. We trained the models using the training data and evaluated their performances using the test data. When training these models with field data, ground roll removed data are required; therefore, the ground roll is suppressed by f-k filtering and used as the ground-truth data. To evaluate the performance of the deep learning models and compare the training results, we utilized quantitative indicators such as the correlation coefficient and structural similarity index measure (SSIM) based on the similarity to the ground-truth data. The DnCNN model exhibited the best performance, and we confirmed that other models could also be applied to suppress the ground roll.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.176-184
• /
• 2017

To interpret the seismic image precisely, random noises should be suppressed and the continuity of the image should be enhanced by using the appropriate smoothing techniques. Structure-Oriented Filter-Edge Preserving (SOF-EP) technique is one of the methods, that have been actively researched and used until now, to efficiently smooth seismic data while preserving the continuity of signal. This technique is based on the principle that diffusion occurs from large amplitude to small one. In a continuous structure such as a horizontal layer, diffusion or smoothing is operated along the layer, thereby increasing the continuity of layers and eliminating random noise. In addition, diffusion or smoothing across boundaries at discontinuous structures such as faults can be avoided by employing the continuity decision factor. Accordingly, the precision of the smoothing technique can be improved. However, in the case of the structure-oriented semblance technique, which has been used to calculate the continuity factor, it takes lots of time depending on the size of the filter and data. In this study, we first implemented the SOF-EP method and confirmed its effectiveness by applying it step by step to the field data. Next, we proposed and applied the corner response method which can efficiently calculate the continuity decision factor instead of structure-oriented semblance. As a result, we could confirm that the computation time can be reduced by about 6,000 times or more by applying the corner response method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.227-232
• /
• 2005

Least squares (${\ell}^2-norm$) solutions of seismic inversion tend to be very sensitive to data points with large errors. The ${\ell}^p-norm$ minimization for $1{\le}p<2$ gives more robust solutions, but usually with higher computational cost. Iteratively reweighted least squares (IRLS) gives efficient approximate solutions of these ${\ell}^p-norm$ problems. I propose a simple way to implement the IRLS method for a hybrid ${\ell}^1/{\ell}^2$ minimization problem that behaves as ${\ell}^2$ fit for small residual and ${\ell}^1$ fit for large residuals. Synthetic and a field-data examples demonstrates the improvement of the hybrid method over least squares when there are outliers in the data.

• Economic and Environmental Geology
• /
• v.48 no.1
• /
• pp.25-39
• /
• 2015

$CO_2$ geological storage plays an important role in reduction of greenhouse gas emissions, but there is a lack of research for CCS demonstration. To achieve the goal of CCS, storing $CO_2$ safely and permanently in underground geological formations, it is essential to understand the characteristics of them, such as total storage capacity, stability, etc. and establish an injection strategy. We perform the impedance inversion for the seismic data acquired from the Ulleung Basin in 2012. To review the possibility of $CO_2$ storage, we also construct porosity models and extract attributes of the prospects from the seismic data. To improve the quality of seismic data, amplitude preserved processing methods, SWD(Shallow Water Demultiple), SRME(Surface Related Multiple Elimination) and Radon Demultiple, are applied. Three well log data are also analysed, and the log correlations of each well are 0.648, 0.574 and 0.342, respectively. All wells are used in building the low-frequency model to generate more robust initial model. Simultaneous pre-stack inversion is performed on all of the 2D profiles and inverted P-impedance, S-impedance and Vp/Vs ratio are generated from the inversion process. With the porosity profiles generated from the seismic inversion process, the porous and non-porous zones can be identified for the purpose of the $CO_2$ sequestration initiative. More detailed characterization of the geological storage and the simulation of $CO_2$ migration might be an essential for the CCS demonstration.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.25-30
• /
• 1998

Time-lapse crosswell seismic data, recorded before and after the cavity filling, showed that the filling increased the velocity at a known cavity zone in an old mine site in Inchon area. The seismic response depicted on the tomogram and in conjunction with the geologic data from drillings imply that the size of the cavity may be either small or filled by debris. In this study, I attempted to evaluate the filling effect by analyzing velocity measured from the time-lapse tomograms. The data acquired by a downhole airgun and 24-channel hydrophone system revealed that there exists measurable amounts of source statics. I presented a methodology to estimate the source statics. The procedure for this method is: 1) examine the source firing-time for each source, and remove the effect of irregular firing time, and 2) estimate the residual statics caused by inaccurate source positioning. This proposed multi-step inversion may reduce high frequency numerical noise and enhance the resolution at the zone of interest. The multi-step inversion with different starting models successfully shows the subtle velocity changes at the small cavity zone. The inversion procedure is: 1) conduct an inversion using regular sized cells, and generate an image of gross velocity structure by applying a 2-D median filter on the resulting tomogram, and 2) construct the starting velocity model by modifying the final velocity model from the first phase. The model was modified so that the zone of interest consists of small-sized grids. The final velocity model developed from the baseline survey was as a starting velocity model on the monitor inversion. Since we expected a velocity change only in the cavity zone, in the monitor inversion, we can significantly reduce the number of model parameters by fixing the model out-side the cavity zone equal to the baseline model.

• The Journal of Engineering Geology
• /
• v.29 no.2
• /
• pp.153-162
• /
• 2019

Rapid variations in the geometry (i.e., thickness) of the refractor and low velocities affect greatly the imaging of the reflectors of land seismic data. Conventional solutions to obtain the weathering models utilizes first break picking method, which requires time consuming steps and causes the human error in picking the first arrivals. A new interferometric approach (interferometric refraction statics, IRS) which utilizes the first arrival signal (S/N enhanced by refraction convolution stack) instead of first break picking, is tested in this study to the synthetic data from the velocity structure provided by surface geophysics (refraction, MASW) and borehole geophysics (tomography, SPS logging) for the Cheongju granitic bodies. The results of IRS approach are found to be better than the ones from conventional first break picking in terms of continuities and horizontal resolution of the reflectors. The unresolved long-wavelength statics in brute stack are much removed by IRS weathering correction and the overlying refractors (the base of shallow weathering zone) are incidentally delineated in the refraction convolution stack.

• Geophysics and Geophysical Exploration
• /
• v.25 no.1
• /
• pp.14-25
• /
• 2022

The passive surface wave method measures seismic signals from ambient noises or vibrations of natural phenomena without using an artificial source. Since passive sources are usually in lower frequencies than artificial ones being able to ensure the information on deeper geological structures, the passive surface wave method can investigate deeper geological structures. In the passive method, frequency dispersion curves are obtained after data acquisition, and the dispersion curves are analyzed by assuming 1D-layered earth, which is like the method of active surface wave survey. However, when computing dispersion curves, the passive method first obtains and analyzes coherence curves of received signals from a set of receivers based on spatial autocorrelation. In this review, we explain how passive surface wave methods measure signals, and make data processing and interpretation, before analyzing field application cases.