• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.225-233
• /
• 2013

The accuracy of the automatic NMO velocity analysis, which is used for an effective and objective NMO velocity analysis, is highly affected by the velocity resolution of the velocity spectrum. In this study, we have developed an automatic NMO velocity algorithm, where the velocity spectra are created using high-resolution bootstrapped differential semblance (BDS), and the velocity analysis on CMP gathers is performed in parallel with MPI. We also compared the velocity models from the developed automatic NMO velocity algorithm with high-resolution BDS to those from BDS. To verify the developed automatic velocity analysis module we created synthetic seismic data from a velocity model including horizon layers. We confirmed that the developed automatic velocity analysis module estimated velocity more accurately. In addition, NMO velocity which yielded a CMP stacked section, where the coherency of the events were improved, was estimated when the developed module was applied to a marine field data set.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.93-101
• /
• 2006

In this study, we propose a joint inversion method, using genetic algorithms, to determine the shear-wave velocity structure of deep sedimentary layers from receiver functions and surface-wave phase velocity. Numerical experiments with synthetic data indicate that the proposed method can avoid the trade-off between shear-wave velocity and thickness that arises when inverting the receiver function only, and the uncertainty in deep structure from surface-wave phase velocity inversion alone. We apply the method to receiver functions obtained from earthquake records with epicentral distances of about 100 km, and Rayleigh-wave phase velocities obtained from a microtremor array survey in the Kanto Plain, Japan. The estimated subsurface structure is in good agreement with the previous results of seismic refraction surveys and deep borehole data.

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

• Journal of the Korean Geophysical Society
• /
• v.5 no.1
• /
• pp.29-39
• /
• 2002

Geophysical data sets from the Chojeong area in the Chungbok-Do are compositely studied in terms of multi-attribute interpretations for the subsurface mappings of shallow fracture zones, associated with groundwater reservoir. Utilizing a GIS software, the attribute data were implemented to a database; a lineament from the satellite image, electrical resistivities and its standard deviation, radioactivity, seismic velocity, and bedrock depth. In an attempt to interpret 1-D electrical sounding data in 3-D views, 1-D data are firstly performed horizontal and vertical inter- and extrapolation. Reconstruction of a resistivity volume is found to be an effective scheme for subsurface mapping of shallow fracture zones. Shallow fracture zones are located in the southeastern part of the study area, which are commonly correlated with the various exploration data.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.227-233
• /
• 2011

Attention has been focused on geothermal energy as an alternative energy because it is continuously operable without external supply. Most of geothermal anomalies in Korea are related to deep circulation of groundwater through a fracture system in granite area. Therefore it is very important to understand the distribution of the fracture system which is the main channel of ground water. In this research, we constructed the velocity models with a fracture system and the layered sediments, respectively, and generated synthetic data sets with them to verify the presented vertical seismic profiling (VSP) preprocessing scheme. We compared the results from conventional VSP preprocessing flow to those from VSP preprocessing flow considering fracture system. We noticed that the preprocessing flow considering fracture system retains more sufficient signal including down-going wave than conventional preprocessing. In addition, we applied 3D VSP prestack phase screen migration to the preprocessed reversed VSP (RVSP) data from Seokmo Island so that we were able to image fracture structure of the geothermal site in Seokmo Island.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.60-66
• /
• 2006

Three differing sandstones, two synthetic and one field sample, have been tested ultrasonically under a range of confining pressures and pore pressures representative of in-situ reservoir pressures. These sandstones include: a synthetic sandstone with calcite intergranular cement produced using the CSIRO Calcite In-situ Precipitation Process (CIPS); a synthetic sandstone with silica intergranular cement; and a core sample from the Otway Basin Waarre Formation, Boggy Creek 1 well, from the target lithology for a trial $CO_2$ pilot project. Initial testing was carried on the cores at "room-dried" conditions, with confining pressures up to 65 MPa in steps of 5 MPa. All cores were then flooded with $CO_2$, initially in the gas phase at 6 MPa, $22^{\circ}C$, then with liquid-phase $CO_2$ at a temperature of $22^{\circ}C$ and pressures from 7 MPa to 17 MPa in steps of 5 MPa. Confining pressures varied from 10 MPa to 65 MPa. Ultrasonic waveforms for both P- and S-waves were recorded at each effective pressure increment. Velocity versus effective pressure responses were calculated from the experimental data for both P- and S-waves. Attenuations $(1/Q_p)$ were calculated from the waveform data using spectral ratio methods. Theoretical calculations of velocity as a function of effective pressure for each sandstone were made using the $CO_2$ pressure-density and $CO_2$ bulk modulus-pressure phase diagrams and Gassmann effective medium theory. Flooding the cores with gaseous phase $CO_2$ produced negligible change in velocity-effective stress relationships compared to the dry state (air saturated). Flooding with liquid-phase $CO_2$ at various pore pressures lowered velocities by approximately 8% on average compared to the air-saturated state. Attenuations increased with liquid-phase $CO_2$ flooding compared to the air-saturated case. Experimental data agreed with the Gassmann calculations at high effective pressures. The "critical" effective pressure, at which agreement with theory occurred, varied with sandstone type. Discrepancies are thought to be due to differing micro-crack populations in the microstructure of each sandstone type. The agreement with theory at high effective pressures is significant and gives some confidence in predicting seismic behaviour under field conditions when $CO_2$ is injected.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.211-216
• /
• 2013

In general, smoothing filters regularize functions by reducing differences between adjacent values. The smoothing filters, therefore, can regularize inverse solutions and produce more accurate subsurface structure when we apply it to full waveform inversion. If we apply a smoothing filter with a constant coefficient to subsurface image or velocity model, it will make layer interfaces and fault structures vague because it does not consider any information of geologic structures and variations of velocity. In this study, we develop a selective smoothing regularization technique, which adapts smoothing coefficients according to inversion iteration, to solve the weakness of smoothing regularization with a constant coefficient. First, we determine appropriate frequencies and analyze the corresponding wavenumber coverage. Then, we define effective maximum wavenumber as 99 percentile of wavenumber spectrum in order to choose smoothing coefficients which can effectively limit the wavenumber coverage. By adapting the chosen smoothing coefficients according to the iteration, we can implement multi-scale full waveform inversion while inverting multi-frequency components simultaneously. Through the successful inversion example on a salt model with high-contrast velocity structures, we can note that our method effectively regularizes the inverse solution. We also verify that our scheme is applicable to field data through the numerical example to the synthetic data containing random noise.

• Geophysics and Geophysical Exploration
• /
• v.23 no.3
• /
• pp.117-130
• /
• 2020

Noises can distort acquired geophysical data, leading to their misinterpretation. Potential noises sources include anthropogenic activity, natural phenomena, and instrument noises. Conventional denoising methods such as wavelet transform and filtering techniques, are based on subjective human investigation, which is computationally inefficient and time-consuming. Recently, many researchers attempted to implement neural networks to efficiently remove noise from geophysical data. This study aims to review and analyze different types of neural networks, such as artificial neural networks, convolutional neural networks, autoencoders, residual networks, and wavelet neural networks, which are implemented to remove different types of noises including seismic, transient electromagnetic, ground-penetrating radar, and magnetotelluric surveys. The review analyzes and summarizes the key challenges in the removal of noise from geophysical data using neural network, while proposes and explains solutions to the challenges. The analysis support that the advancement in neural networks can be powerful denoising tools for geophysical data.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.19-28
• /
• 2007

In this study, we propose a joint inversion method, using genetic algorithms, to estimate an S-wave velocity structure for deep sedimentary layers from receiver functions and surface-wave phase velocity observed at several sites. The method takes layer continuity over a target area into consideration by assuming that each layer has uniform physical properties, especially an S-wave velocity, at all the sites in a target area in order to invert datasets acquired at different sites simultaneously. Numerical experiments with synthetic data indicate that the proposed method is effective in reducing uncertainty in deep structure parameters when modelling only surface-wave dispersion data over a limited period range. We then apply the method to receiver functions derived from earthquake records at one site and two datasets of Rayleigh-wave phase velocity obtained from microtremor array surveys performed in central Tokyo, Japan. The estimated subsurface structure is in good agreement with the results of previous seismic refraction surveys and deep borehole data. We also conclude that the proposed method can provide a more accurate and reliable model than individual inversions of either receiver function data only or surface-wave dispersion data only.