• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.37-44
• /
• 2016

$CO_2$-EOR (Carbon Dioxide-Enhanced Oil Recovery), one of the enhanced oil recovery methods, helps to not only enhance the production of oil, but also store carbon dioxide in underground. However, if micro fractures occur when during the injection of $CO_2$, it is difficult to make permanent storage of $CO_2$ in reservoir and can cause contamination of groundwater and soil. Therefore, in this study, we performed microseismic monitoring to investigate the occurrence of fractures during the $CO_2$ injection at the Meruap oil reservoir, Indonesia. To pick the first arrivals of microseismic events, Improved MER (Modified Energy Ratio) method was used. After picking the first arrivals, hodogram analysis was carried out by using the data recorded at three component geophones to calculate the back azimuth of events. Finally, locations of microseismic events were decided by using the results of first arrival picking and hodogram analysis. Estimated locations showed that all microseismic events were occurred at surface and any fracture did not occur around the reservoir. Moreover, by analyzing noise characteristic, we confirmed that almost of picked first arrivals were due to the repetitive mechanical noise.

• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.11-19
• /
• 2016

Two representative residual static methods of traveltime decomposition and stack-power maximization are discussed in terms of application to land seismic data. For the model data with synthetic shot/receiver statics (time shift) applied and random noises added, continuities of reflection event are much improved by stack-power maximization method, resulting the derived time-shifts approximately equal to the synthetic statics. Optimal parameters (maximum allowable shift, correlation window, iteration number) for residual statics are effectively chosen with diagnostic displays of CSP (common shot point) stack and CRP (common receiver point) stack as well as CMP gather. In addition to removal of long-wavelength time shift by refraction statics, prior to residual statics, processing steps of f-k filter, predictive deconvolution and time variant spectral whitening are employed to attenuate noises and thereby to minimize the error during the correlation process. The reflectors including horizontal layer of reservoir are more clearly shown in the variable-density section through repicking the velocities after residual statics and inverse NMO correction.

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.192-198
• /
• 2009

Side scan sonar and SBP (sub-bottom profiler) play a very important role in the survey for seafloor imaging and sub-bottom profiling. In this study, we have acquired side scan sonar and SBP data from the artificial reef area. We applied digital image processing techniques to side scan sonar data in order to improve an image quality. For the enhancement of data quality and image resolution, we applied the typical seismic data processing sequence including gain recovery, muting, spectrum analysis, predictive deconvolution, migration to SBP data. We could easily estimate if artificial reef structures were settled properly and their distribution on the seafloor from the integrated interpretation of side scan sonar and SBP data. From the sampling analysis of seabed sediments, texture filtering of side scan sonar data and SBP data interpretation, we could evaluate the sediment type, distribution and thickness of seafloor sediments in detail.

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.173-182
• /
• 2009

In April 2008, KIGAM carried out an ocean-bottom seismometer (OBS) survey in the central Ulleung Basin where strong bottom simulating reflectors (BSRs) were revealed from previous surveys and some gas-hydrate samples were retrieved by direct sampling. The purpose of this survey is to estimate the velocity structure near the BSR in the gas hydrate prospect area using wide-angle seismic data recorded on the ocean-bottom seismometers. Along with the OBS survey, a 2-D seismic survey was performed whereby stratigraphic and preliminary velocity information was obtained. Two methods were applied to wide-angle data for estimating P wave velocity; one is velocity analysis in the $\tau$-p domain and the other is seismic traveltime inversion. A 1-D interval velocity profile was obtained by the first method, which was refined to layered velocity structure by the latter method. A layer stripping method was adopted for modeling and inversion. All velocity profiles at each OBS site clearly show velocity reversal at BSR depths due to the presence of gas hydrates. In addition, we could confirm high velocity in the column/chimney structure.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.41-51
• /
• 2008

To construct in situ velocity-porosity relationships for oceanic basalt, considering crack features, P- and S-wave velocity measurements on basaltic samples obtained from the eastern flank of the Juan de Fuca Ridge were carried out under confining pressures up to 40 MPa. Assuming that the changes in velocities with confining pressures are originated by micro-crack closure, we estimated micro-crack aspect ratio spectra using the Kuster-$Toks{\ddot{o}}z$ theory. The result demonstrates that the normalised aspect ratio spectra of the different samples have similar characteristics. From the normalised aspect ratio spectrum, we then constructed theoretical velocity-porosity relationships by calculating an aspect ratio spectrum for each porosity. In addition, by considering micro-crack closure due to confining pressure, a velocity-porosity relationship as a function of confining pressure could be obtained. The theoretical relationships that take into account the aspect ratio spectra are consistent with the observed relationships for over 100 discrete samples measured at atmospheric pressure, and the commonly observed pressure dependent relationships for a wide porosity range. The agreement between the laboratory-derived data and theoretically estimated values demonstrates that the velocity-porosity relationships of the basaltic samples obtained from the eastern flank of the Juan de Fuca Ridge, and their pressure dependence, can be described by the crack features (i.e. normalised aspect ratio spectra) and crack closure.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.15-25
• /
• 2008

Ocean-bottom seismometer (OBS) positions are one of the key parameters in an OBS-airgun seismic survey for crustal structure study. To improve the quality of these parameters, we have developed a new method of determining OBS positions, using airgun shot data and bathymetric data in addition to available distance measurements by acoustic transponders. The traveltimes of direct water waves emitted by airgun shots and recorded by OBSs are used as important information for determining OBS locations, in cases where there are few acoustic transponder data (<3 sites). The new method consists of two steps. A global search is performed as the first step, to find nodes of the bathymetric grid that are the closest to explaining the observed direct water-wave traveltimes from airgun shots, and acoustic ranging using a transponder system. The use of precise 2D bathymetric data is most important if the bottom topography near the OBS is extremely rough. The locations of the nodes obtained by the first step are used as initial values for the second step, to avoid falling into local convergence minima. In the second step, a non-linear inverse method is executed. If the OBS internal clock shows large drift, a secondary correction for the OBS internal clock is obtained, as well as the OBS location, as final results by this method. We discuss the error and the influence of each measurement used in the determination of OBS location.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.137-145
• /
• 2017

High-grade limestone applied to various chemical industries is abundant within upper Pungchon formation in Taebaeksan basin, South Korea. Geophysical exploration is one of the most efficient methods to investigate subsurface geological structure in an extensive area. Since the geophysical exploration for the high-grade limestone has rarely been conducted in Korea, its appropriate strategy has not been set up yet. In this study, we focused on to suggest the reasonable strategy and accumulate geophysical databases which are essential for interpreting geophysical images by characterizing laboratory physical properties of in-situ rocks. Hence, rocks were obtained from drilled cores consisting of lower Hwajeol formation, Pungchon formation, and dykes in Jeongseon area, Gangwon province. Geophysical laboratory experiments and petrography of the rocks were conducted. Since susceptibility values of the rocks in Pungchon Formation were obviously lower than those of upper Hwajeol and dykes, it is considered that the lithological boundaries could be distinguished by magnetic survey. In addition, electrical properties of the rocks in middle Pungchon formation were relatively different compared with those of upper/lower Pungchon formations. Thus, induced polarization is shown to be able to detect the high-grade limestone in upper Pungchon formation.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.185-192
• /
• 2017

A sequence of earthquakes with the main shock $M_L$ 5.8 occurred on September 12 2016 in the Gyeongju area. The main shock was the largest earthquakes in the southern part of the Korean peninsula since the instrumental seismic observation began in the peninsula in 1905 and clearly demonstrated that the Yangsan fault is seismically active. The mean focal depth of the foreshock, main shock, and aftershock of the Gyeongju earthquakes estimated by the crustal model of single layer of the Korean peninsula without the Conrad discontinuity turns out to be 12.9 km, which is 2.8 km lower than that estimated based on the IASP91 reference model with the Conrad discontinuity. The distribution of the historical and instrumental earthquakes in the Gyeongju area indicates that the Yangsan fault system comprising the main Yangsan fault and its subsidiary faults is a large fracture zone. The epicenters of the Gyeongju earthquakes show that a few faults of the Yangsan fault system are involved in the release of the strain energy accumulated in the area. That the major earthquakes of Gyeongju earthquakes occurred not on the surface but below 10 km depth suggests the necessity of the study of the distribution of deep active faults of the Yangsan fault system. The magnitude of maximum earthquake of the Gyeongju area estimated based on the earthquake data of the area turns out to be 7.3. The recurrence intervals of the earthquakes over magnitudes 5.0, 6.0 and 7.0 based on the earthquake data since 1978, which is the most complete data in the peninsula, are estimated as 80, 670, and 5,900 years, respectively. The September 2016 Gyeongju earthquakes are basically intraplate earthquakes not related to the Great East Japan earthquake of March 11 2011 which is interplate earthquake.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.163-175
• /
• 2017

In land seismic exploration, irregular surface topography and weathering layer in near surface distorts the reflected signals of data. Therefore, typical land seismic data should be compensated for this distortion by static correction. To perform the static correction, near-surface velocity is required, which can be obtained by seismic refraction survey. However, land seismic data is often acquired in a limited form of geometry depending on the equipment availability, accessibility condition, and permission for the survey site. In this situation, refraction analysis should be performed using reflection data because it is impossible to acquire refraction-oriented data due to limited source and receiver geometry. In this study, we aimed to analyze the reliability of the results obtained by refraction traveltime tomography when using reflection data with a limited number of sources and receivers from irregular surface topography. By comparing the inversion result from irregular topography with that from flat surface, we found that the surface topography affects the reliability of the inversion results to some degree. We also found that the number of sources has little effect on the inversion results unless the number of sources are very small. On the other hand, we observed that velocity distortion occurred in the overlapped part of receiver arrays when using a limited number of receivers, and therefore suggested the size of the least overlapping ratio to avoid the velocity distortion. Finally, we performed numerical tests for the model which simulates the surface topography and acquisition geometry of the survey region and verified the reliability analysis of inversion results. We identified reliable areas and suspicious area of the inverted velocity model by applying the analysis results to field data.