• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.383-392
• /
• 2007

Karstic features and mining-related cavities not only lead to severe restrictions in land utilizations, but also constitute serious concern about geohazard and groundwater contamination. A microgravity survey was applied for detecting, mapping and monitoring karstic cavities in the test site at Muan prepared by KIGAM. The gravity data were collected using an AutoGrav CG-3 gravimeter at about 800 stations by 5 m interval along paddy paths. The density distribution beneath the profiles was drawn by two dimensional inversion based on the minimum support stabilizing functional, which generated better focused images of density discontinuities. We also imaged three dimensional density distribution by growing body inversion with solution from Euler deconvolution as a priori information. The density image showed that the cavities were dissolved, enlarged and connected into a cavity network system, which was supported by drill hole logs. A time-lapse microgravity was executed on the road in the test site for monitoring the change of the subsurface density distribution before and after grouting. The data were adjusted for reducing the effects due to the different condition of each survey, and inverted to density distributions. They show the change of density structure during the lapsed time, which implies the effects of grouting. This case history at the Muan test site showed that the microgravity with accuracy and precision of ${\mu}Gal$ is an effective and practical tool for detecting, mapping and monitoring the subsurface cavities.

• Journal of the Korean Geophysical Society
• /
• v.2 no.1
• /
• pp.39-44
• /
• 1999

High-resolution images are drawn from existing seismic data which were originally obtained by Korea Ocean Research & Development Institute (KORDI) during 1994-1997 for deep seismic studies on the East Sea of Korea. These images are analyzed for mapping Quaternary faults and near-bottom gas pockets. First 12 channels are selected from shot gathers for reprocessing. The processing sequence adopted for high-resolution seismic images comprises data copy, trace editing, true amplitude recovery, common-midpoint sorting, initial muting, prestack deconvolution, bandpass filtering, stacking, highpass filtering, poststack deconvolution, f-x migration, and automatic gain control (AGC). Among these processing steps, predictive deconvolution, highpass filtering, and short window AGC are the most significant in enhancement of resolution. More than 200 Quaternanry faults are interpreted on the migrated sections in the shallow depths beneath the seafloor. Although numerous faults are found mostly at the western continental slope and boundaries of the Ulleung Basin, significant amount of the faults are also indicated within the basin. Many of these faults are believed to be formed with reactivation of basement, from geotectonic activities including volcanism, and often originated in Tertiary, indicating that the tectonic regime of the East Sea might be unstable. Existence of shallow gas pockets casts real hazardous warnings to deep-sea drillings and/or to underwater constructions such as inter-island cables and gas pipelines. On the other hand, discovery of these gas pockets heightens the interests in developing natural resources in the East Sea. Reprocessed seismic sections, however, show no typical seismic characteristics for gas hydrates such as bottom-simulating reflectors in the western continental slope and ocean floor.

• The Korean Journal of Petroleum Geology
• /
• v.2 no.2 s.3
• /
• pp.43-50
• /
• 1994

The geologic structure of Gongju Basin, which is a Cretaceous sedimentary basin located on the boundary of Gyeonggi Massif and Ogcheon Belt, is modeled by using gravity data and interpreted in relation with basin forming tectonism. The electrical survey with dipole-dipole array was also conducted to uncover the development of fractures in the two fault zones which form the boundaries of the basin. In the process of gravity data reduction, the terrain correction was performed by using the conic prism model, which showed better results specially for topography having a steep slope. The gravity model of the geologic structure of Gongju basin is obtained by forward modeling based on the surface geology and density inversion. It reveals that the width of the basin at its central part is about $4{\cal}km$ and about $2.5{\cal}km$ at the southern part. The depth of crystalline basement beneath sedimentary rocks of the basin is about $700{\~}400{\cal}m$ below the sea level and it is thinner in the center than in margin. The fault of the southeastern boundary appears more clearly than that of the northwestern boundary, and its fracture zone may extended to the depth of more than $1{\cal}km$. Therefore, it is thought that the tectonic movement along the fault in the southeastern boundary was much stronger. These results coincide with the appearance of broad low resistivity anomaly at the southeastern boundary of the basin in the resistivity section. The fracture zones having low density are also recognized inside the basin from the gravity model. The swelling feature of basement and the fractures in sedimentary rocks of the basin suggest that the compressional tectonic stress had also involved after the deposition of the Cretaceous sediments.

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

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2010

Seismic reflection surveying is one of the most widely used and effective techniques for coal seam structure delineation and risk mitigation for underground longwall mining. However, the ability of the method can be compromised by the presence of volcanic cover. This problem arises within parts of the Bowen and Sydney Basins of Australia and seismic surveying can be unsuccessful. As a consequence, such areas are less attractive for coal mining. Techniques to improve the success of seismic surveying over basalt flows are needed. In this paper, we use elastic wave-equation-based forward modelling techniques to investigate the effects and characteristics of seismic wave propagation under different settings involving changes in basalt properties, its thickness, lateral extent, relative position to the shot position and various forms of inhomogeneity. The modelling results suggests that: 1) basalts with high impedance contrasts and multiple flows generate strong multiples and weak reflectors; 2) thin basalts have less effect than thick basalts; 3) partial basalt cover has less effect than full basalt cover; 4) low frequency seismic waves (especially at large offsets) have better penetration through the basalt than high frequency waves; and 5) the deeper the coal seams are below basalts of limited extent, the less influence the basalts will have on the wave propagation. In addition to providing insights into the issues that arise when seismic surveying under basalts, these observations suggest that careful management of seismic noise and the acquisition of long-offset seismic data with low-frequency geophones have the potential to improve the seismic results.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.80-87
• /
• 2010

This study was conducted to develop a reservoir modelling workflow to reproduce the heterogeneous distribution of effective permeability that impacts on the performance of SAGD (Steam Assisted Gravity Drainage), the in-situ bitumen recovery technique in the Athabasca Oil Sands. Lithologic facies distribution is the main cause of the heterogeneity in bitumen reservoirs in the study area. The target formation consists of sand with mudstone facies in a fluvial-to-estuary channel system, where the mudstone interrupts fluid flow and reduces effective permeability. In this study, the lithologic facies is classified into three classes having different characteristics of effective permeability, depending on the shapes of mudstones. The reservoir modelling workflow of this study consists of two main modules; facies modelling and permeability modelling. The facies modelling provides an identification of the three lithologic facies, using a stochastic approach, which mainly control the effective permeability. The permeability modelling populates mudstone volume fraction first, then transforms it into effective permeability. A series of flow simulations applied to mini-models of the lithologic facies obtains the transformation functions of the mudstone volume fraction into the effective permeability. Seismic data contribute to the facies modelling via providing prior probability of facies, which is incorporated in the facies models by geostatistical techniques. In particular, this study employs a probabilistic neural network utilising multiple seismic attributes in facies prediction that improves the prior probability of facies. The result of using the improved prior probability in facies modelling is compared to the conventional method using a single seismic attribute to demonstrate the improvement in the facies discrimination. Using P-wave velocity in combination with density in the multiple seismic attributes is the essence of the improved facies discrimination. This paper also discusses sand matrix porosity that makes P-wave velocity differ between the different facies in the study area, where the sand matrix porosity is uniquely evaluated using log-derived porosity, P-wave velocity and photographically-predicted mudstone volume.

• Korean Journal of Remote Sensing
• /
• v.28 no.5
• /
• pp.477-487
• /
• 2012

Sea surface winds in the sea off the east coast of Korea were derived from L-band ALOS (Advanced Land Observing Satellite) PALSAR (Phased Array type L-band Synthetic Aperture Radar) data and their characteristics of errors were analyzed. We could retrieve high-resolution wind vectors off the east coast of Korea including the coastal region, which has been substantially unavailable from satellite scatterometers. Retrieved SAR-wind speeds showed a good agreement with in-situ buoy measurement by showing relatively small an root-mean-square (RMS) error of 0.67 m/s. Comparisons of the wind vectors from SAR and scatterometer presented RMS errors of 2.16 m/s and $19.24^{\circ}$, 3.62 m/s and $28.02^{\circ}$ for L-band GMF (Geophysical Model Function) algorithm 2009 and 2007, respectively, which tended to be somewhat higher than the expected limit of satellite scatterometer winds errors. L-band SAR-derived wind field exhibited the characteristic dependence on wind direction and incidence angle. The previous version (L-band GMF 2007) revealed large errors at small incidence angles of less than $21^{\circ}$. By contrast, the L-band GMF 2009, which improved the effect of incidence angle on the model function by considering a quadratic function instead of a linear relationship, greatly enhanced the quality of wind speed from 6.80 m/s to 1.14 m/s at small incident angles. This study addressed that the causes of wind retrieval errors should be intensively studied for diverse applications of L-band SAR-derived winds, especially in terms of the effects of wind direction and incidence angle, and other potential error sources.

• Economic and Environmental Geology
• /
• v.49 no.6
• /
• pp.469-477
• /
• 2016

To evaluate the $CO_2$ storage capacity for the reservoir rock, the laboratory scale technique to measure the amount of $scCO_2$, replacing pore water of the reservior rock after the $CO_2$ injection was developed in this study. Laboratory experiments were performed to measure the $scCO_2$ displacement efficiency of the conglomerate and the sandstone in Janggi basin, which are classified as available $CO_2$ storage rocks in Korea. The high pressurized stainless steel cell containing two different walls was designed and undisturbed rock cores acquired from the deep drilling site around Janggi basin were used for the experiments. From the lab experiments, the average $scCO_2$ displacement efficiency of the conglomerate and the sandstone in Janggi basin was measured at 31.2% and 14.4%, respectively, which can be used to evaluate the feasibility of the Janggi basin as a $scCO_2$ storage site in Korea. Assuming that the effective radius of the $CO_2$ storage formations is 250 m and the average thickness of the conglomerate and the sandstone formation under 800 m in depth is 50 m each (from data of the drilling profile and the geophysical survey), the $scCO_2$ storage capacity of the reservoir rocks around the probable $scCO_2$ injection site in Janggi basin was calculated at 264,592 metric ton, demonstrating that the conglomerate and the sandstone formations in Janggi basin have a great potential for use as a pilot scale test site for the $CO_2$ storage in Korea.

• Economic and Environmental Geology
• /
• v.38 no.6 s.175
• /
• pp.643-656
• /
• 2005

The East Sea, a marginal sea or back-arc basin, consists of Japan Basin, Yamato Basin, and Ulleung Basin and is surrounded by the Pacific Plate and Philippine Sea Plate. Ulleung Basin locates in the southwestern part of the East Sea and shows the depth of 1,500 m in average and 2,500 m in maximum, connecting to the Japan Basin along 2,000 m contour. The slope of the seafloor is greater in the western side of the basin than in the southern and the eastern side. The crustal thickness of the Ulleung Basin from the OBS tends to get thicker toward the north and the west side and the sediment thickness of the Ulleung Basin is getting thicker toward the southeast side and reaches up to 12 km. The crustal type of the Ulleung Basin was variously suggested as like as a rifted continental crust, an extended continental crust, and an incipient oceanic trust. The origin of the crustal formation and the Ulleung Basin, however, is still controversial. Based on the bathymetry and gravtiy anomaly data for this study, the axis of the Ulleung Basin shows that the basin develops along the axis trending NW-SE direction and reveals a general symmetry of the bathymetry. And also the free-air gravity anomalies show a very similar pattern to the bathymetry of the basin. The sediment thickness is relatively thicker in the southeastern side of the basin than in the northwestern side. Although the crustal age of the Ulleung Basin is supposed to be younger than them of the Japan Basin and the Yamato Basin, the free-air gravity anomalies of the Ulleung Basin ranging -40 to 50 mGals are lower than the other basins, which suggests that the densities of crust and sediment of the Ulleng Basin are lower than the Japan Basin and the Yamato Basin.

• Geophysics and Geophysical Exploration
• /
• v.13 no.3
• /
• pp.243-248
• /
• 2010

Two suspicious events, which were claimed as underground nuclear tests by North Korea, were detected in the northern Korean Peninsula on October 9, 2006 and May 25, 2009. The KIGAM and Korea-China Joint seismic stations are distributed uniformly along the boundaries between North Korea and adjacent countries. In this study, the data from broadband stations with the distance of 200 to 550 km from the test site are used to analyze and compare two nuclear tests of North Korea. By comparing the time differences of the Pn-wave arrival times of 1st and 2nd tests at multiple stations, the relative locations of two test sites could be calculated precisely. From the geometrical calculation with the velocity of Pn wave $V_{Pn}$ = 8 km/s, the 2nd test site is estimated to move in the WNW direction from 1st one with the distance of 2 km. Body wave magnitude, mb of the 2nd test, which was announced officially as the network average of 4.5, varies widely with the directional location of stations from 4.1 to 5.2. The magnitude obtained from Lg wave, $m_b$(Lg), shows less variation between 4.3 to 4.7 with the average of 4.6. The moving-window spectra of time traces of 1st and 2nd tests show very similar pattern with different scale level. In addition, the corner frequencies of P wave of 1st and 2nd tests at each station show no or negligible difference. This indicates the burial depths of two tests might be very similar. The relative yield amount of the 2nd test is estimated 8 times larger than that of the 1st from the weighted average of ground-velocity amplitude ratios.