• The Journal of Engineering Geology
• /
• v.25 no.3
• /
• pp.349-356
• /
• 2015

Although geophysical methods are useful and generally provide valuable information about the subsurface, it is important to recognize their limitations. A common limitation is the lack of sufficient contrast in physical properties between different layers. Thus, multiple methods are commonly used to best constrain the physical properties of different layers and interpret each section individually. Ground penetrating radar (GPR) and shallow seismic reflection (SSR) methods, used for shallow and very shallow subsurface imaging, respond to dielectric and velocity contrasts between layers, respectively. In this study, we merged GPR and SSR data from a test site within the Cheongui granitic mass, where the water table is ~3 m deep all year. We interpreted the data in combination with field observations and existing data from drill cores and well logs. GPR and SSR reflections from the tops of the sand layer, water table, and weathered and soft rocks are successfully mapped in a single section, and they correlate well with electrical resistivity data and SPS (suspension PS) well-logging profiles. In addition, subsurface interfaces in the integrated section correlate well with S-wave velocity structures from multi-channel analysis shear wave (MASW) data, a method that was recently developed to enhance lateral resolution on the basis of CMP (common midpoint) cross-correlation (CMPCC) analysis.

• Journal of Soil and Groundwater Environment
• /
• v.27 no.spc
• /
• pp.19-33
• /
• 2022

When subsurface is polluted, contaminants tend to migrate through groundwater flow path. The groundwater flow path is highly dependent upon underground geological structures in the contaminated area. Geophysical survey is an useful tool to identify subsurface geological structure. In addition, geophysical logging in a borehole precisely provides detailed information about geological characteristics in vicinity of the borehole, including fractures, lithology, and groundwater level. In this work, surface seismic refraction and electrical resistivity surveys were conducted in a test site located in Namyangju city, South Korea, along with well logging tests in five boreholes installed in the site. Geophysical data and well logging data were collected and processed to construct an 3D geological map in the site.

• 한국지구물리탐사학회:학술대회논문집
• /
• 1999.08a
• /
• pp.92-116
• /
• 1999

Tomographic approaches to image underground structure using electrical properties, can be divided into DC resistivity, electromagnetic, and radar tomography, based on the operating frequency. DC resistivity and radar tomography methods have been recently applied to site investigation for engineering purpose in Korea. This paper review these two tomography methods, through the case histories acquired in Korea. As another method of borehole radar survey, borehole radar reflection method is included, and its inherent problem and solution are discussed, how to find the azimuth angle of reflector using direction-finding-antenna. Since the velocity anisotropy of radar wave has been commonly encountered in field data, anisotropic radar tomography is discussed in this paper. In DC resistivity tomography, two subjects are focussed, electrode arrays, and borehole effect owing to the conductive fluid in borehole. Using the numerical modeling data, various kinds of electrode ways are compared, and borehole effect is illustrated. Most of the case histories presented in this paper are compared with known geology, core logging data, and/or Televiewer images.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.30-36
• /
• 2006

Japan's first pilot-scale $CO_2$ sequestration experiment has been conducted in Nagaoka, where 10400 t of $CO_2$ have been injected in an onshore aquifer at a depth of about 1100 m. Among various measurements conducted at the site for monitoring the injected $CO_2$, we conducted time-lapse crosswell seismic tomography between two observation wells to determine the distribution of $CO_2$ in the aquifer by the change of P-wave velocities. This paper reports the results of the crosswell seismic tomography conducted at the site. The crosswell seismic tomography measurements were carried out three times; once before the injection as a baseline survey, and twice during the injection as monitoring surveys. The velocity tomograms resulting from the monitoring surveys were compared to the baseline survey tomogram, and velocity difference tomograms were generated. The velocity difference tomograms showed that velocity had decreased in a part of the aquifer around the injection well, where the injected $CO_2$ was supposed to be distributed. We also found that the area in which velocity had decreased was expanding in the formation up-dip direction, as increasing amounts of $CO_2$ were injected. The maximum velocity reductions observed were 3.0% after 3200 t of $CO_2$ had been injected, and 3.5% after injection of 6200 t of $CO_2$. Although seismic tomography could map the area of velocity decrease due to $CO_2$ injection, we observed some contradictions with the results of time-lapse sonic logging, and with the geological condition of the cap rock. To investigate these contradictions, we conducted numerical experiments simulating the test site. As a result, we found that part of the velocity distribution displayed in the tomograms was affected by artefacts or ghosts caused by the source-receiver geometry for the crosswell tomography in this particular site. The maximum velocity decrease obtained by tomography (3.5%) was much smaller than that observed by sonic logging (more than 20%). The numerical experiment results showed that only 5.5% velocity reduction might be observed, although the model was given a 20% velocity reduction zone. Judging from this result, the actual velocity reduction can be more than 3.5%, the value we obtained from the field data reconstruction. Further studies are needed to obtain more accurate velocity values that are comparable to those obtained by sonic logging.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2000.09a
• /
• pp.54-69
• /
• 2000

There are groundwater trouble by high-salinity yield inducing sea-water intrusion in Cheju Island. It is used groundwater-GIS(Well-lnfo) in the maintenance and management of groundwater in Cheju Island to grasp groundwater trouble area and cause of high-salinity yield. For 16 wells certain to yield high-salinity, we logged specific electrical conductivity(EC) and tried to get hold of freshwater and saltwater relationship. As result of distribution of $Cl^-$ by depth, it is showed up groundwater trouble by high-salinity yield in the east coastal area and the partly north coastal area. The reason of high-salinity groundwater yield are low-groundwater level by the structure of geology and low-hydraulic gradient etc. There is necessity for management to development and use of groundwater in the high-salinity area, special management area.

• 한국지구물리탐사학회:학술대회논문집
• /
• 1999.08a
• /
• pp.137-155
• /
• 1999

Recently, Televiewer(Borehole Acoustic Scanner(Televiewer)) has come to be widely used specially for the general engineering construction design. The Televiewer tool using a focussed acoustic beam is to detect the amplitude and traveltime of each reflected acoustic signal at the wall, resulting in the amplitude- and traveltime image respectively. Fractures can be well detected, because they easily scatter the acoustic energy due to the highly narrow beam. In addition, the drilling work will rough the borehole wall so that the acoustic energy can be scattered simply due to the roughness of the wall. Thus, the amplitude level can be directed associated with the elastic properties(impedance) and the hardness of the rock as well. Meanwhile, the traveltime image provides an information about the borehole shape and can be converted to a high precision 3D caliper log(max. 288 arms). In this paper, based on the high resolution of Televiewer images, general evaluation methods are illustrated to derive very reliable rock parameters.

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

• Economic and Environmental Geology
• /
• v.43 no.6
• /
• pp.615-623
• /
• 2010

Seismic refraction survey is a geophysical method that delineates subsurface velocity structure using direct wave and critically refracted wave. The generalized reciprocal method(GRM) is an inversion technique which uses travel-time data from several forward and reverse shots and which can provide the geometry of irregular inclined refractors and structures underlain by hidden layer such as low velocity zone and thin layer. In this study, a simple Excel-GRM routine was tested for fast mapping of the interface between weathering layer and bedrock during the survey, with employing a pair of forward and reverse shots. This routine was proved to control the maximum dip of approximately $30^{\circ}C$ and maximum velocity contrast of 0.6, based on the panel tests in terms of dipping angle and velocity contrast for the two-layer inclined models. In contrast with conventional operation of five to seven shots with sufficient offset distance and indoor data analysis thereafter, this routine was performed in the field shortly after data acquisition. Depth to the bedrock provided by Excel-GRM, during the field survey for Cheongju granite area, correlates well with the elevation of the surface of soft rock from the drill core and SPS logging data. This cost-effective routine developed for quickly delineating the bedrock surface in the field survey will be readily applicable to mapping of weathering zone in narrow zone with small variation of elevation of bedrock.

• Tunnel and Underground Space
• /
• v.26 no.3
• /
• pp.166-180
• /
• 2016

The present study numerically simulated the CO2 injection into the saline aquifer of CO2CRC Otway pilot project and the resulting hydrological-mechanical coupled process in the storage site by TOUGH-FLAC simulator. A three-dimensional numerical model was generated using the stochastic geological model which was established based on well log and core data. It was estimated that the CO2 injection of 30,000t over a period of 200 days increased the pressure near the injection point by 0.5 MPa at the most. The pressure increased rapidly and tended to approach a certain value at an early stage of the injection. The hydrological and mechanical behavior observed from the CO2 flow, effective stress change and stress-strength ratio revealed that the CO2 injection into the saline aquifer under the given condition would not have significant effects on the mechanical safety of the storage site and the hydrological state around the adjacent fault.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.242-249
• /
• 2008

Recently, AVO analysis has been widely used in oil exploration with seismic subsurface section as a direct indicator of the existence of the gas. In the case of the deep reservoirs like the gas reservoirs in the East-sea, it is often difficult to observe AVO responses in CMP gathers even though the bright spots are shown in the stacked section. Because the reservoir becomes more consolidated as its depth deepens, P-wave velocity does not decrease significantly when the pore fluid is replaced by the gas. Thus the difference in Poisson's ratio, which is a key factor for AVO response, between the reservoir and the layer above it does not increase significantly. In this study, we analyzed the effects of Poisson's ratio difference on AVO response with a variety of Poisson's ratios for the upper and lower layers. The results show that, as the difference in Poisson's ratio between the upper and lower layers decreases, the change in the reflection amplitude with incidence angle decreases and AVO responses become insignificant. To consider the limitation of AVO responses shown in the gas reservoir in East-sea, the velocity model was made by simulation Gorae V structure with seismic data and well logs. The results of comparing AVO responses observed from the synthetic data with theoretical AVO responses calculated by using material properties show that the amount of the change in reflection amplitude with increasing incident angle is very small when the difference in Poisson's ratio between the upper and lower layers is small. In addition, the characteristics of AVO responses were concealed by noise or amplitude distortion arisen during preprocessing. To overcome such limitations of AVO analysis of the data from deep reservoirs, we need to acquire precisely reflection amplltudes In data acquisition stage and use processing tools which preserve reflection amplitude in data processing stage.