• Economic and Environmental Geology
• /
• v.43 no.3
• /
• pp.259-267
• /
• 2010

Resistivity logging instruments were designed to measure electrical resistivity of formation, which can be directly interpreted to provide water-saturation profile. Short and long normal logging measurements are made under groundwater level. In some investigation sites, groundwater level reaches to a depth of a few meters. It has come to attention that the proximity of groundwater level might distort short and long normal logging readings, when the measurements are made near groundwater level, owing to the proximity of an insulating air. This study investigates the effects of the proximity of groundwater level (and also the proximity of earth surface) on the normal by simulating normal logging measurements near groundwater level. In the simulation, we consider all the details of real logging situation, i.e., the presence of wellbore, the tool mandrel with current and potential electrodes, and currentreturn and reference-potential electrodes. We also model the air to include the earth’'s surface in the simulation rather than the customary choice of imposing a boundary condition. To obtain apparent resistivity, we compute the voltage, i.e., potential difference between monitoring and reference electrodes. For the simulation, we use a twodimensional, goal-oriented and high-order self-adaptive hp finite element refinement strategy (h denotes the element size and p the polynomial order of approximation within each element) to obtain accurate simulation results. Numerical results indicate that distortion on the normal logging is greater when the reference potential electrode is closer to the borehole and distortions on long normal logging are larger than those on short normal logging.

• Journal of the Korean Geophysical Society
• /
• v.6 no.3
• /
• pp.131-137
• /
• 2003

In this paper the liquefaction of sand bed under oscillating water pressure are treated as a basic study of the prevention works against the scouring around the hydraulic structures. The results of the former resurch show that the occurrence of the liquefaction depends on both properties of the oscillating water pressure and of the sand layer. Considering the latter properties, that is , the resistivity against the liquefaction increases with the increase of the permeability of the sand bed, we propose the displcement method as one of the prevention works, which is a method to displace the upper layer of the sand bed by the sand with large permeability. The effects of this method are investigated theoretically and experimentally. By the experimental study, it is shown that the proposed displacement method has the apparent effect to prevent the liquefaction. The experimental results are explained fairly well by the theoretical analysis based on the theory of the flow through the elastic porous media.

• Journal of the Korean earth science society
• /
• v.29 no.6
• /
• pp.447-453
• /
• 2008

Recently, the electrical DC survey has frequently been performed to assess the safety of embankment. This study showed that the damaged section of embankment could be appropriately detected by the survey only when the three dimensional effect was correctly considered. The shape of the three dimensional embankment was numerically implemented, and a proper modeling was performed to confirm the effect by analyzing the apparent and inverted true resistivity. Then, the field work was carried out. The three dimensional interpretation distinguished the erroneous weak zones from the geometrical artifact, and the embankment was ensured as safe both by the additional survey performed in rainy season and the partial excavation for direct observation.

• Journal of Soil and Groundwater Environment
• /
• v.7 no.3
• /
• pp.3-17
• /
• 2002

Data analysis of groundwater monitoring wells and geostatistical methods are used to identify the local characteristics of sea water intrusion and the range of sea water intrusion at the southeastern coastal area of Busan, Korea. Rainfall and groundwater level of two monitoring wells show a linear correlation because of the direct groundwater recharge by the precipitation. However, rainfall and electric conductivity have the inverse relationship because of the increase of groundwater. Electric conductivity rapidly increased at 24m depth and exceeded 20,000$\mu\textrm{s}$/cm near 26m depth in the monitoring wells. The variations of groundwater level and electric conductivity show that the interface between sea water and fresh water tends to move upward when groundwater level goes down. In the cross correlation analysis, groundwater level versus rainfall represents the largest cross correlation coefficient in 0 time lag but the cross correlation coefficient of electric conductivity versus rainfall is the largest when the time lag is 9 days. This suggests that the fluctuations of groundwater level respond to rainfall in a short time, but the interface between sea water and fresh water respond very slow to rainfall. Horizontal extents of sea water intrusion are estimated to 14 m from the east of Line 1, and 25 m from the southeast end of Line 2 in the inversion of dipole-dipole profiling data of two survey lines. The data of VES by the Schulumberger array in May and July show lognormal distributions. In the kriged apparent resistivity and earth resistivity distributions, the resistivities of July are increased comparing to those of May. This reflects that the concentration of sea water in aquifer is reduced due to the increased groundwater recharge from the rainfall in June and July. In analyzing the vertical and horizontal apparent resistivities and earth resistivity distributions, the geostatistical methods are very useful to identify the variations of earth resistivity distributions at the coastal area.

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.239-245
• /
• 2009

To verify the reinforcing effect of grouting materials composed of colloid cement and ordinary portland cement on the water leakage region in a small scale dike, we performed a tubecasing method and applied surface electrical resistivity survey including electrical resistivity tomography (ERT) to find resistivity variation before and after grouting. Hydraulic conductivities after grouting show 10 times lower than those of before grouting. These variation indicates that the cement grout blocks the leakage pathway effectively. As the results of dipole-dipole resistivity survey along the dike, resistivity distribution after grouting did not represent noticeable spatial variation in time. Resistivity monitoring results at the dike with vertical electrical sounding (VES) showed that the region of decreasing apparent resistivity was occupied by the grout after grouting. Predicted resistivities from the inversion of ERT data well matched with results of VES at the same regions. From the ERT using check holes to inspect the effect of grouting, we could find that the ERT is quite effective to identify spatially the grout region in a dike.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.322-331
• /
• 2007

In this study, electrical resistivity tomography (ERT) were conducted to find the mineralized zone at the Geumpung mine in Dojeon-ri, Susan-myeon, Jecheon-si, Chungcheongbuk-do. The deviation of the inclined borehole was measured to obtain the exact positions of the electrodes for correcting apparent resistivity values from ERT. Geophysical loggings such as resistivity and natural gamma were conducted to obtain the properties of the material near the borehole. Measurements of the physical properties of the cores, such as porosity, water content, density, susceptibility, resistivity were performed to analyze the correlation between physical properties and resistivity. Grade analysis for core sample was also conducted to identify relationship between grade and resistivity. Rock property analysis shows that the resistivity is more dominated by susceptibility and grade than by porosity and water content in the mineralized zone. The results of ERT are well consistent with geophysical logging data and geologic column. So ERT is powerful method to identify conductive mineralized zone.

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.1-6
• /
• 2000

The 1 $km^2$ area studied is located in Sukchun-ri, Hwasung-koon, the southern part of Kyeonggi-do. Even though this site has been known as a contaminated area caused by seawater intrusions, geophysical and geochemical surveys have never been carried out at the site to determine the extent of the seawater contamination and to investigate whether the seawater intrusion is in progress. The purpose of this study is to determine the extent of seawater contamination and a preferred channel of the seawater intrusion using geophysical methods such as DC resistivity surveys with Schlumberger array and a dipole-dipole array. In order to determine whether the seawater intrusion is in progress in the area, DC resistivity monitoring with Schlumberger array was performed. According to the resistivity map obtained from the inversion of the resistivity data measured with Schlumberger array, the study area is divided into two districts as relatively lowly resistive (less than 30 ohm-m) and highly resistive (more than 30 ohm-m) areas. The distribution of the lowly resistive area is consistent with the distribution of the layer composed of clay minerals, and the resistivity of this layer decreases slowly as approaching to the old seashore. Hydrogeological analysis shows that the clay layer within a distance of about 200 m from the seashore has been already contaminated by sea-water and its electric conductivity is 8 times higher than that of the sand layer covered by the clay layer. According to the results of the 2-dimensional DC resistivity surveys with a dipole-dipole array, there are two preferred channels of the seawater intrusion in the site, and both the channels are in the NW-SE direction from the old seashore. The lowly resistive zone in the southern channel extends to a depth of 80 m. The DC resistivity monitoring with Schlumberger array was carried out along the preferred channel which has the low resistivity Bone (fracture zone) that extended to a depth of 80 m. The time series of apparent resistivity, measured at a distance of 260 m from the old coast line, fluctuates with a period of 12 hours. From these observations, it can be concluded that the seawater intrusion caused by tidal action is still in progress along the fractured zone interpreted by the DC resistivity surveys with a dipole-dipole array.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.37-50
• /
• 2024

The magnetotelluric (MT) survey can be affected by external environmental factors. In particular, when acquiring MT data in islands, it is essential to consider the combined effect of topography and sea to understand the results and make accurate interpretations. To analyze the MT response (apparent resistivity, phase) with consideration of the effect of topography and sea, a small cone-shaped island model surrounded by deep sea was created. Two-dimensional (2-D) and three-dimensional (3-D) forward modeling were performed on the terrain model considering topography and the island model considering both topography and sea. The 2-D MT response did not reflect the topographic and sea effect of the direction orthogonal to the 2-D profile. The 3-D MT response included topographic and sea effects in all directions. The XY and YX components of the apparent resistivity were separated on undulating topography, such as a hill. A conductor at 1 km below sea level could be distinguished from topographic and sea effects in the MT response, and low resistivity anomaly was attenuated at greater depths. This study will facilitate understanding of field data measured on small islands.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.390-397
• /
• 2003

A three-dimensional (3D) inversion technique has been developed for interpretation of magnetotelluric (MT) data. The inversion method is based on the linearized least-squares (Gauss-Newton) method with smoothness regularization. In addition to the underground 3D resistivity distribution, static shifts are also treated as unknown parameters in the inversion. The forward modeling is by the staggered-grid finite difference method. A Bayesian criterion ABle is applied to search the optimum trade-off among the minimization of the data misfit, model roughness and static shifts. The method has been applied to several MT datasets obtained at geothermal fields in Japan and other Asian countries. In this paper, two examples will be discussed: one is the data at the Ogiri geothermal area, southwestern Japan, and the other is at the Pohang low-enthalpy geothermal field, southeastern Korea. The inversion of the Ogiri data has been performed stably, resulting in a good fitting between the observed and computed apparent resistivities and phases. The recovered 3D resistivity structure is generally similar to the two-dimensional (2D) inversion models, although the deeper portion of the 3D model seems to be more realistic than that of the 2D model. The 3D model is also in a good agreement with the geological model of the geothermal reservoirs. 3D interpretation of the Pohang MT data is still preliminary. Although the fitting to the observed data is very good, the preliminary 3D model is not reliable enough because the station coverage is not sufficient for a 3D inversion.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.275-284
• /
• 2007

For precise interpretation of magnetotelluric (MT) data distorted by irregular surface terrain, topography effects are investigated by computing apparent resistivities, phases, tippers and induction vectors for a three-dimensional (3D) hill-and-valley model. To compute MT responses for the 3D surface topography model, we use a 3D MT modeling algorithm based on an edge finite-element method which is free from vector parasites. Distortions on the apparent resistivity and phase are mainly caused by distorted currents that flow along surface topography. The distribution of tipper amplitudes over both hill and valley are the same, while the tipper points toward the center of hill and the base of the valley. The real part of induction vector also points in the same direction as that of tipper, while the imaginary part in the opposite direction.