• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.58-68
• /
• 2011

As part of basic studies of monitoring carbon dioxide ($CO_2$) storage using electrical and seismic surveys, laboratory experiments have been conducted to measure resistivity and P-wave velocity changes due to the injection of $CO_2$ into water-saturated sandstone. The rock sample used is a cylinder of Berea sandstone. $CO_2$ was injected under supercritical conditions (10 MPa, $40^{\circ}C$). The experimental results show that resistivity increases monotonously throughout the injection period, while P-wave velocity and amplitude decrease drastically due to the supercritical $CO_2$ injection. A reconstructed P-wave velocity tomogram clearly images $CO_2$ migration in the sandstone sample. Both resistivity and seismic velocity are useful for monitoring $CO_2$ behaviour. P-wave velocity, however, is less sensitive than resistivity when the $CO_2$ saturation is greater than ~20%. The result indicates that the saturation estimation from resistivity can effectively complement the difficulty of $CO_2$ saturation estimations from seismic velocity variations. By combining resistivity and seismic velocity we were able to estimate $CO_2$ saturation distribution and the injected $CO_2$ behaviour in our sample.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.33-42
• /
• 1999

AVO and complex trace analyses mainly used to characterize natural gas reservoir were tested in this paper for a possible application to detection of major geological discontinuities such as fracture zones. The test data used in this study were calculated by utilizing a viscoelastic numerical program which was based on the generalized Maxwell body for a horizontal fracture model. In AVO analysis of a horizontal fracture zone, p-wave reflection appears to be variant depending upon the acoustic-impedence contrast and the offset distance. The fracture zone is also effectively clarified both in gradient stack and range-limited stack in which fracture zone reflection is attenuated with the increasing offset distance. In complex attribute plots (instantaneous amplitude, frequency, and phase), the top and bottom of the fracture Tone are characterized by a zone of strong amplitudes and an event of the same phase. Low frequency characteristics appear at the fracture zone and the underneath. Amplitude attenuation and waveform dispersion are dependent on Q-contrast between the fracture zone and the surrounding media. They were properly compensated by optimum inverse Q-filtering.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.17-25
• /
• 1999

A traveltime tomography has been carried out by transforming electromagnetic data in frequency domain to wave-like domain. The transform uniquely relates a field satisfying a diffusion equation to an integral of the corresponding wavefield. But direct transform of frequency domain magnetic fields to wave-field domain is ill-posed problem because the kernel of the integral transform is highly damped. In this study, instead of solving such an unstable problem, it is assumed that wave-fields in transformed domain can be approximated by sum of ray series. And for further simplicity, reflection and refraction energy compared to that of direct wave is weak enough to be neglected. Then first arrival can be approximated by calculating the traveltime of direct wave only. But these assumptions are valid when the conductivity contrast between background medium and the target anomalous body is low enough. So this approach can only be applied to the models with low conductivity contrast. To verify the algorithm, traveltime calculated by this approach was compared to that of direct transform method and exact traveltime, calculated analytically, for homogeneous whole space. The error in first arrival picked by this study was less than that of direct transformation method, especially when the number of frequency samples is less than 10, or when the data are noisy. Layered earth model with varying conductivity contrasts and inclined dyke model have been successfully imaged by applying nonlinear traveltime tomography in 30 iterations within three CPU minutes on a IBM Pentium Pro 200 MHz.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.8-16
• /
• 1999

In terms of high resolution, seismic reflection survey is by far the most significant geophysical method applied to define subsurface structure. In shallow seismic reflection survey, it is, however, difficult to obtain high resolution image due to both the wave attenuation in the unconsolidated layer and the existence of source-generated surface waves Therefore, when collecting data, it is imperative to select proper equipments and choose optimum field data acquisition parameters for acquiring high S/N data. In this survey, a small size hammer was used as a low energy source and 40-Hz vertical geophones were used as receivers. Trigger signal was obtained from the hammer starter attached in the aluminum plate and thus it was possible to control the source onset time for the vertical stack. During the field work, a modified standard CMP technique was introduced to achieve the many-fold CMP data effectively. Data processing was conducted by the 'Seismic Unix' which is mounted on PC with a Linux operating system. The main distinctions were the emphasis and detail placed on near-surface velocity analysis and the extra care exercised in muting.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.4
• /
• pp.228-237
• /
• 2015

Field observations were conducted to collect hydrodynamic and morphological data, which are needed to account for mechanisms of bathymetry changes caused by physical forcings, in Haeundae beach. In order to quantitatively describe characteristics of wave transformations and current patterns in space in winter and summer, in-situ sensors for measuring waves and current profiles were installed at three locations in the cross-shore direction and also three locations in the along-shore direction. As for the results of wave measurements, waves with main direction from the east dominate in winter while waves are incident from the S and the ESE in summer. Analysis of current data reveals that currents over the study domain are considerably influenced by a pattern of tidal motions, thereby, mainly oscillating in the direction of tidal currents, i.e., east-west directions, in both winter and summer. Currents tend to be influenced by local bathymetry in the shallow water region, with the direction changed along the depth contours and the magnitude reduced as they approach the shoreline. The results analysed from the hydrodynamic data through this study can be further combined with the morphological and bathymetry data, leading to the quantification of seasonal sediment transport rates and sand budget changes.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.10 no.2
• /
• pp.118-124
• /
• 2007

The oceans could absorb almost all the anthropogenic carbon dioxide the mankind has been producing eventually, but in the nature the air-sea $CO_2$ exchange occurs very slowly and to lower the atmospheric $CO_2$ concentration substantially $CO_2$ must be injected to the interior of the ocean directly. If we inject $CO_2$ collected at the major $CO_2$ sources into the international waters in the Philippine Sea or east of Japan, we could store the $CO_2$ in the oceans effectively for a few hundred years. When $CO_2$ is dissolved into the water, PH drops. The creatures adapted to the deep oceans where environment is very stable could be affected by even a small change in pH significantly. If, therefore, we are to inject $CO_2$ into the oceans, we must assess the effect of $CO_2$ injection in the marine ecosystem beforehand. Only when the damage to the marine ecosystem is smaller than the benefit from the $CO_2$ injection, $CO_2$ injection is effective.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.29 no.5
• /
• pp.543-550
• /
• 2011

The gravity anomaly is a basic geophysical data applied in various fields such as geophysics, geodesy and national defense. In general, the gravity anomaly is used through a interpolation process based on the constructed database. The gravity variation, however, is appeared in various shapes depending on the topography and the density of the underground structures. Therefore, the interpolation could lead to a large differences if the gravity fields do not satisfy the assumptions on the signal behavior like linear or a certain degree polynomials. Furthermore, the interpolation does not reflect the physical characteristics of the gravity such as the harmonic condition. In this study, the gravity modeling using the plane Fourier series and radial basis functions are performed to overcome the problems in the usual interpolation. The results of the modeling is analyzed for the case of the gravity referenced navigation focused on the signal characteristics. Based on the study, it was found that the results from modeling are not much different to that from the interpolation in a smoothly varied area. In case of the highly varied area, however, a large differences are appeared among the three methods. Especially, the Fourier series shows the most smooth variations in the modeled gravity values while the highest variations appeared in the interpolation. Applying to the gravity referenced navigation, it was found that the modeling is more effective in calculation cost. It is considered that the results from this study provides a basis on effective modeling of the gravity fields in terms of the signal characteristics and resolution for various application fields.

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

In order to investigate the velocity structure of the southern part of the Korean peninsula, exploded seismic signals were recorded for 120 s along a 294-km WNW-ESE line and 150 s along a 335-km NNW-SSE line in 2002 and 2004, respectively. Velocity tomograms were derived from inverting first arrival times. One-dimensional velocity models derived by joint analyses of teleseismic receiver functions and surface wave dispersion at several stations near the profiles were uesd to build initial models. The raypaths indicate several midcrust interfaces including ones at approximate depths of 2.0 and 14.9 km with refraction velocities of approximately 6.0 and 7.1 km/s, respectively. The deepest significant interface varies in depth from 30.8 km to 36.1 km. The critically refracting velocity varies from 7.8 to 8.1 km/s along this interface which may correspond to the Moho discontinuity. The velocity tomograms show (1) existence of a low-velocity zone centered at 6-7 km depth under the Okchon fold belt, (2) extension of the Yeongdon fault down to greater than 10 km, and (3) existence of high-velocity materials under the Gyeongsan basin whose thickness is less than 4.2 km.

• Geophysics and Geophysical Exploration
• /
• v.17 no.3
• /
• pp.137-146
• /
• 2014

This study presents a procedure for discrimination of artificial events from earthquakes occurred in and around the Korean Peninsula using data set in the Wonju KSRS seismograph network, Korea. Two training sets representing natural and artificial earthquakes were constructed with 150 and 56 events, respectively, with high signal to noise ratio. A frequency band, Pg(4-6 Hz)/Lg(5-7 Hz), which is optimal for the discrimination of seismic sources was derived from the two-dimensional grid of Pg/Lg spectral amplitude ratio. The corrections for the effects of earthquake magnitude and hypocentral distance were carried out for improvement of discrimination capability. For correcting the effect of magnitude dependence due to the inverse proportionality of corner frequency to seismic moment, the Brune's source spectrum was subtracted from the observation spectrum. The spectrum was corrected using the optimal damping coefficient to remove damping effect with the hypocentral distance. The effect of locally varying spectrum ratio was cancelled correcting variation of wave propagation along the ray path. The performance in discrimination between training sets of natural and artificial events was compared using the Mahalanobis distance in each step of correction. The procedure of magnitude, distance, and path corrections show clear improvements of the discrimination results with increasing Mahalanobis distance, from 1.98 to 3.01, between two training sets.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.257-267
• /
• 2013

3-D ore modeling was performed to understand the configuration of ore bodies by integrated analysis of borehole and geophysical data in iron-mine area. Five representative indices of rocks were designated, which were obtained from geological survey and borehole. The five indices of rocks were geostatistically simulated by Sequential Indicator Simulation method to delineate boundary of the ore bodies. And Ordinary Kriging and Sequential Gaussian Simulation was applied to make secondary information using resistivity data from magnetotellurics and DC resistivity survey, and this information was used for simple kriging with local varying means, one of integrated kriging techniques. From the correlation analysis between each properties, it was found that high grade of ore is characterized by increased density, whereas the electrical resistivity decreases. With the integrated results of geophysical and borehole data, it was also found that the real configuration of ore body was similar to the modeled result and information about ore grade in 3-D space was obtained.