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

AVO analysis was conducted on hydrocarbon-bearing structures by applying the crossplot and offset-coordinate amplitude polynomial techniques. To evaluate the applicability of the AVO analysis, it was conducted on synthetic data that were generated with an anticline model, and field data from the hydrocarbon-bearing Colony Sand bed in Canada. Analysis of synthetic data from the anticline model demonstrates that the crossplot method yields zero-offset reflection amplitude and amplitude variation with negative values for the upper interface of the hydrocarbon-bearing layer. The crossplot values are clustered in the third quadrant. The results of AVO analysis based on the coefficients of the amplitude polynomial are similar to those from the crossplots. These well correlated results of AVO analysis on field and synthetic data suggest that both methods successfully investigate the characteristics of the reflections from the upper interface of a hydrocarbon-bearing layer. Analysis based on the incident-angle equation facilitates the application of various interpretation methods. However, it requires the conversion of seismic data to an incident angle gather. By contrast, analysis using coefficients of the amplitude polynomial is cost-effective because it allows examining amplitude variation with offset without involving the conversion process. However, it warrants further investigation into versatile application. The two different techniques can be complement each other effectively as AVO-analysis tools for the detection of hydrocarbon reservoirs.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• KEPCO Journal on Electric Power and Energy
• /
• v.3 no.1
• /
• pp.29-34
• /
• 2017

Leaching kinetics for extracting yttrium from the coal fly ash was investigated in the presence of sulfuric acid during extraction. The leaching kinetics of yttrium were conducted at reactant densities of 5~1,000 g coal fly ash per L of $1.0{\sim}10.0N\;H_2SO_4$, agitation speed of 250 rpm and temperature ranging from 30 to $90^{\circ}C$. As a result, the leaching kinetic model was determined in a two-step model based on the shrinking core model with spherical particles. The first step was proceeded by chemical reaction at ash surface, and the second step was proceeded by ash layer diffusion because the leaching conversion of yttrium by the first chemical reaction increases with increased the time irrelevant to the temperature whereas it increases with increased the leaching temperature. The activation energy of the first chemical leaching step was determined to be $1.163kJmol^{-1}$. After the first chemical reaction, the activation energy of ash layer diffusion leaching was derived to be $41.540kJmol^{-1}$. The optimum conditions for leaching the yttrium metal of 60 % were found to be the slurry density of 250 g fly ash per L of $H_2SO_4$, solvent concentration of $2.0N\;H_2SO_4$, second step leaching of temperatures of $30^{\circ}C$ for 3 hours and then $90^{\circ}C$ for 3 hours at agitation rate of 250 rpm.

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

Two-dimensional (2-D) inversion of magnetotelluric (MT) data for two survey lines having south-north direction from Jeiu Island has been carried out. The 2-D models show a thick layer having around 10 ohm-m in the depth of a few hundred meters throughout the survey area, which can be considered as the unconsolidated sedimentary layer. And they also show a conductive anomaly at the central part of each survey lines. But unfortunately by now, we do not have any further information about the anomaly. Comparison of the 2-D inversion model using MT band only and that using both AMT and MT bands said that it is helpful for us to include AMT band as well as MT band in the inversion to interpret not only the shallow part but also the deep structures.

• Journal of Environmental Science International
• /
• v.20 no.2
• /
• pp.275-289
• /
• 2011

In order to clarify the relation between sea breeze penetration and Planetary Boundary Layer development in southeastern part of the Korean Peninsula, several numerical assessments were carried out using atmospheric numerical model WRF(Weather Research and Forecasting). Compared with onset time of sea breeze at eastern coast area(Uljin), the time at southern coast region(Masan) with complex costal line tend to delay for several hours. The penetration patterns of sea breeze between two coastal regions are some different due to the shape of their coastal line and back ground topography. Intensified valley wind due to high topography of lee side of Uljin can help penetration of sea breeze at early time. So penetration of sea breeze at early time often prevent PBL to develop at Uljin and lower PBL height last for a day time. But because of late penetration of sea breeze at Masna, PBL Height dramatically decrease after 1500LST. The distribution of front genesis function based on the heat and momentum variation are explained obviously the sea breeze penetration patterns and agreed well with the PBL height distribution.

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

The crustal structure of Korean Peninsula is investigated by analyzing phase velocity dispersion data of Rayleigh wave. Earthquakes recorded by three component seismographs during 1999 - 2004 in South Korea are used in this study. The fundamental mode signals of Rayleigh waves are obtained from vertical components of seismograms by multiple filter technique method and phase match filter method. Velocity dispersion curves of surface waves for 14 propagation paths on the great circle are computed from the fundamental mode signals on the great circle path by two-station method. Treating the shear velocity of each layer as an independent parameter, phase velocities of Rayleigh wave are inverted. The result models are regarded as average structure for surface wave propagation paths respectively. All the results can be explained by an earth model of the Korean Peninsula comprising crust of shear-wave velocity increasing from 2.8 to 3.25 km/sec from top to 33 km depth and uppermost mantle of shear-wave velocity between 4.55 and 4.67 km/sec.

• The Journal of Engineering Geology
• /
• v.27 no.3
• /
• pp.305-312
• /
• 2017

The most effective way to distinguish subsurface interfaces that produce various geophysical responses is through the integration of multiple geophysical methods, with each method detecting both a complementary and unique set of distinct physical properties relating to the subsurface. In this study, shallow seismic reflection (SSR) and ground penetrating radar (GPR) surveys were conducted at the Cheongju-Gadeok site of the Korea National Groundwater Monitoring Network to map the water table, which was measured at 12 m depth during the geophysical surveys. The water table proved to be a good target reflector in both datasets, as the abrupt transition from the overlying unsaturated weathered rock to the underlying saturated weathered rock yielded large acoustic impedance and dielectric constant contrasts. The two datasets were depth converted and integrated into a single section, with the SSR and GPR surveys conducted to ensure subsurface imaging at approximately the same wavelength. The GPR data provided detailed information on the upper ~15 m of the section, whereas the SSR data imaged structures at depths of 10-45 m. The integrated section thus captured the full depth coverage of the sandy clay, water table, weathered rock, soft rock, and hard rock structures, which correlated well with local drillcore and water table observations. Incorporation of these two geophysical datasets yielded a synthetic section that resembled a simplified aquifer model, with the best-fitting seismic velocity, dielectric constant, and porosity of the saturated weathered layer being $v_{seismic}=1000m/s$, ${\varepsilon}_r=16$, and ${\phi}=0.32$, respectively.

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

Two-dimensional velocity tomograms of P- and S-waves were obtained by inverting traveltimes of first arrivals. The two sections of shear-wave velocity show similar features as a whole, with smaller values on the section from surface wave dispersions. Difficulties in picking SH-wave phases due to noise and later arrivals than P waves and PS converted waves are experienced. In addition, a flat layer model based on the surface wave inversion prohibits applications of the method where sgear wave velocities vary strongly in the lateral direction.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.