• Journal of Soil and Groundwater Environment
• /
• v.17 no.5
• /
• pp.56-67
• /
• 2012

In the present study, an enhanced subsurface prediction algorithm based on a non-parametric geostatistical model and a history matching technique through Gibbs sampler is developed and the iterative prediction improvement procedure is proposed. The developed model is applied to a simple two-dimensional synthetic case where domain is composed of three different hydrogeologic media with $500m{\times}40m$ scale. In the application, it is assumed that there are 4 independent pumping tests performed at different vertical interval and the history curves are acquired through numerical modeling. With two hypothetical borehole information and pumping test data, the proposed prediction model is applied iteratively and continuous improvements of the predictions with reduced uncertainties of the media distribution are observed. From the results and the qualitative/quantitative analysis, it is concluded that the proposed model is good for the subsurface prediction improvements where the history data is available as a supportive information. Once the proposed model be a matured technique, it is believed that the model can be applied to many groundwater, geothermal, gas and oil problems with conventional fluid flow simulators. However, the overall development is still in its preliminary step and further considerations needs to be incorporated to be a viable and practical prediction technique including multi-dimensional verifications, global optimization, etc. which have not been resolved in the present study.

• Tunnel and Underground Space
• /
• v.32 no.5
• /
• pp.285-297
• /
• 2022

We report results of Extended Leak-Off Test (XLOT) conducted in a large diameter borehole, which is drilled for installation of deep borehole geophysical monitoring system to monitor micro-earthquakes and fault behavior of major fault zones in the southeastern Korean Peninsula. The borehole was planned to secure a final diameter of 200 mm (or more) at a depth of ~1 km, with 12" diameter wellbore to intermediate depths, and 7-7/8" (~200 mm) to the bottom hole depth. We drilled first the 12" borehole to approximately 504 m deep and installed American Petroleum Institute standard 8-5/8" casing, then annulus between the casing and bedrock was fully cemented. XLOT was carried out for several purposes such as confirming casing and cementing integrity, measuring rock stress states. To that end, we drilled additional 4 m long open hole interval to directly inject water and pressurize into the rock mass using the upper API casings. During the XLOT, flow rates and interval pressures were recorded in real time. Based on the logs we tried to analyze hydraulic conductivity of the test interval.

• Economic and Environmental Geology
• /
• v.38 no.2 s.171
• /
• pp.143-153
• /
• 2005

Euiseong subbasin, included in the Kyungsang Basin, was created by the result of volcanic activity in the late Cretaceous, and contacts with Milyang and Youngyang subbasins by Palgongsan and Andong faults, respectively. In this study, geophysical survey is implemented fur investigating surface and subsurface geologic structure in Euiseong subbasin which composed with the complex of volcanic and plutonic rocks. To understand surface geologic feature, IRS satellite image and DEM(Digital Terrain Map) are used for analyzing lineament and its density. The numbers of lineaments show major trend in $N55^{\circ}\~65^{\circ}W$, and aspects of lineament lengths show major trend in $N55^{\circ}\~65^{\circ}W$ and N-S directions. 13 delineate subsurface density discontinuity; Power spectrum analysis was implemented for gravity anomaly data, resulting $4-5{\cal}km$ depth of basin basement and $0.5-0.6{\cal}km$ depth of shallow discontinuity. From the result of power spectrum analysis, 2.5-D modelings were implemented along two profiles of A-A' and B-B', and they show subsurface geology in detail. Analytic signal method for detecting boundaries of magnetic basements show 0.001-130 nT/m values, and high energy area show good correspondency with the boundaries of Palgongsan granite and caldera areas in Euiseong subbasin.

• Economic and Environmental Geology
• /
• v.27 no.6
• /
• pp.563-569
• /
• 1994

Resistivity Image Profiling has been applied to a geological survey in the Keoje-do area. Survey lines are located near the KD-02 and KD-06 wells in the area, where we have already sampled all rock cores and carried out several kinds of geophysical logs. In each site a resistivity cross-section is obtained by inverting pole-pole apparent resistivities. Comparing the reconstructed resistivity section with the rock cores and logging data obtained in the well provides a detailed picture of subsurface geology. The geology of KD-02 site is composed of conglomerate, sandstone and shale with fractures. The reconstructed resistivity image is useful for assessing the grade of weathering of these rocks. The KD-06 site is mainly underlain by granitic rocks, and its fresh basement can be delineated by resistivities over $1000{\Omega}{\cdot}m$.

• Economic and Environmental Geology
• /
• v.30 no.1
• /
• pp.25-33
• /
• 1997

The purpose of this study is to examine the feasibility of using stable isotopes as a hydrologic tracer, and to elucidate the groundwater circulation system and the source of S component dissolved in thermal water of the Chonju Jukrim thermal spring district based on the O, H and S isotopic variabilities of environmental materials including bedrock, rainwater, surface water, shallow subsurface water and thermal spring water. The ${\delta}^{18}O$ and ${\delta}D$ of subsurface waters and surface water show highly restricted range and plotted on the same meteoric water line as a ${\delta}D=8{\delta}^{18}O+19$ line, and derivate from the mean annual isotopic composition of the rain water but are analogous to those of rain waters precipitated during winter season, indicating that ground waters are originated from the meteoric water and are strongly affected by the seasonal variation of air mass. Thermal spring waters are more depleted in ${\delta}^{18}O$ and ${\delta}D$ than those of shallow ground water and surface water. It can be explained by the difference of recharge area. The hydrochemical properties of subsurface waters and surface water devide into two groups: $Ca(HCO_3)_2$ type including shallow subsurface water and surface water, and $Na(HCO_3)$ type of thermal spring waters. The ${\delta}^{34}S$ values of thermal spring water show very high positive and quitely distinct from those of shallow subsurface water and surface water that are similar to those of bed rocks, indicating that sulfate dissolved in thermal spring water has not only a terrigenic origin, but also originates partially from the foreign source containing very heavy ${\delta}^{34}S$ component such as an ancient sea water. However, the presence of $H_2S$ can not be ignore the affact of the isotopic fractionation to explaine the heavy ${\delta}^{34}S$ of thermal spring water. Overall, the Oxygen and Hydrogen stable isotopes can identify the source and the circulation system of the natural waters and the S-isotopes can provide a crucial clue on tracing the dissolved material transports in the circulation system of the natural water.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• Economic and Environmental Geology
• /
• v.21 no.1
• /
• pp.107-113
• /
• 1988

We present a precedure for interpreting dipole-dipole apparent resistivity data. The procedure is constructed by combining a forward two-dimensional finite element modeling and an inverse technique with Householder's transformation. In the interpretation, subsurface structure is divided into some blocks with constant resistivities. Our inversion technique is tested on synthetic and field data. We found that geologic constraint is required for successful interpretation.

• Economic and Environmental Geology
• /
• v.25 no.2
• /
• pp.167-177
• /
• 1992

The gravity measurement has been conducted at 53 and 34 stations with an interval of 1~1.5 km along the national roads of about 47 km and 34 km running from Duksungri to Yangpori and from Angangri to Byungpori, Kyungsangbookdo, respectively. The subsurface geology and geologic structure of Tertiary Pohang and Janggi basins along two survey lines are interpreted quantitatively by applying Fourier series and Talwani methods for Bouguer gravity anomaly. The depths of Conrad discontinuity vary from 11.8 to 12.5 km and 11.5 to 13.2 km along the survey lines between Duksungri and Yangpori, and Angangri and Byungpori, respectively. The depths of pre-Cambrian Gneiss complex underneath Kyungsang Supergroup vary from 3.8 to 4.2 km and 3.8 to 4.6 km along the survey lines between Duksungri and Yangpori, and Angangri and Byungpori, respectively. Massive granite bodies which are not exposed along the survey line between Duksungri and Yangpori are distributed on a large scale at the subsurface between Duksungri and Ochun, and Daegokri and Yangpori. Along the survey line between Angangri and Byungpori, it is exposed at Angangri, and extends underneath Chungrimdong, Pohang city. Andesite is distributed on a small scale underneath Pohang city and Ochun. The thicknesses of Tertiary Yonil and Janggi Groups are 0.2~0.9 km and 0.1~0.5 km, respectively. The Tuffaceous rocks which are the lowest formation of Tertiary sedimentary rocks are distributed with the thickness of 0.2 km at the surface and between Kyungsang Supergroup and Yonil or Janggi Groups. The Yonil and Janggi Groups are in fault contact by a fault running through Ochun and Chungrimdong, Pohang city. Two other faults are newly found near Heunghae-eup and Hyungsan river.

• The Journal of Engineering Geology
• /
• v.12 no.1
• /
• pp.23-34
• /
• 2002

Borehole image analysis has been carried out to obtain the detailed geological data by approach of direct observation. Direct application of borehole image analysis inevitably gives rise to a few of restriction of data acquisition due to the limited information within narrow borehole space. Considering the apparent dip of discontinuity surface depending upon the direction, the visualized program of two-dimensional subsurface discontinuities is coded. Borehole image analysis can compensate the distribution of subsurface discontinuity extending into the expected area of investigation. In order to draw subsurface profile in the proposed area of subsurface construction, visualized program is coded as a window GUI (Graphic User Interface) using Fortran and Visual Basic Programming languages. It is to open publicly for the usage of whoever is in want. Discontinuity distribution map is visualized along the Proposed line of tunnel in the Janggye-ri area, Jangsu-gun. Using the visualized program, the limited information from borehole spatially applies into analysis of overall subsurface structures, and the distributional characteristics of discontinuity anticipate at the proposed area. In addition, spacing and extension of joint and depth of discontinuity effecting tunnel safety can be visualized along the direction of the proposed tunnel. These lines of visualization apply design and construction of fundanmental structures.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.457-465
• /
• 2013

We integrated and correlated datasets from surface and subsurface geophysics, drilling cores, and engineering geology to identify geological interfaces and characterize the joints and fracture zones within the rock mass. The regional geometry of a geologically weak zone was investigated via a fence projection of electrical resistivity data and a borehole image-processing system. Subsurface discontinuities and intensive fracture zones within the rock mass are delineated by cross-hole seismic tomography and analyses of dip directions in rose diagrams. The dynamic elastic modulus is studied in terms of the P-wave velocity and Poisson's ratio. Subsurface discontinuities, which are conventionally identified using the N value and from core samples, can now be identified from anomalous reflection coefficients (i.e., acoustic impedance contrast) calculated using a pair of well logs, comprising seismic velocity from suspension-PS logging and density from logging. Intensive fracture zones identified in the synthetic seismogram are matched to core loss zones in the drilling core data and to a high concentration of joints in the borehole imaging system. The upper boundaries of fracture zones are correlated to strongly negative amplitude in the synthetic trace, which is constructed by convolution of the optimal Ricker wavelet with a reflection coefficient. The standard deviations of dynamic elastic moduli are higher for fracture zones than for acompact rock mass, due to the wide range of velocities resulting from the large numbers of joints and fractures within the zone.