• Journal of the Korean Geophysical Society
• /
• v.9 no.2
• /
• pp.129-134
• /
• 2006

Recently, spillways are need to control stable water level for supporting main dams because of floods by unusual change of weather such as Typhoon Rusa. This study has been focused on the amount of leakage through the rock mass distributed fractures and joints under the opened emergency spillway. It is very important to evaluate the amount of leakage as these affect stability of spillway by interaction between effective stress and pore pressure. The commercial program MAFIC has been used for analyzing groundwater flow in fractured rock mass. The results showed that the values of range, average and deviation of leakage were 2.85∼3.79×10-1, 3.32×10-1 and 1.70×10-2 m3/day/m2 respectively. Secondary, we have estimated the effect of grouting after the transmissivity(Tf) of joint 1 as main pathway of leakage known from above results was changed from 1.78×10-7 to 1.59×10-9 m2/s. The results showed that the values of range, average and deviation of leakage were 7.80×10-4∼1.53×10-3, 1.18×10-3 and 1.32×10-4 m3/day/m2 respectively. As the result, the amount of leakage after grouting has been decreased by a ratio of 1 to 277.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.350-360
• /
• 2008

Electrical resistivity is one of physical property of the earth and measured by electrical resistivity survey, electrical resistivity logging and laboratory test. Recently, electrical resistivity is widely used in determination of rock quality in support pattern design of road and railway tunnel construction sites. To get more reliable rock quality data from electrical resistivity, it needs a lot of test and study on correlation of resistivity and rock quality. Firstly, we did rock property test in laboratory, such as P wave velocity, Young's modulus, uniaxial compressive strength (UCS) and electrical resistivity. We correlate each test results and we found out that electrical resistivity has highly related to P wave velocity, Young's modulus and UCS. Next, we accomplished electrical resistivity survey in field site and carried out electrical resistivity logging at in-situ area. We also performed rock classification, such as RQD, RMR and Q-system and we correlate electrical resistivity to RMR data. We found out that electrical resistivity logging data are highly correlate to RMR. Also we found out that electrical resistivity survey data are lower than electrical resistivity logging data when there are faults or fractures. And it cause electrical resistivity survey data to lowly correlate to RMR.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.286-301
• /
• 2008

When magnetotelluric (MT) data are obtained in vicinity of the coast, the sea can distort observed MT responses, especially those of deep part of subsurface. We introduce an iterative method to correct the sea effect, based on the previous topographic correction method which removes the distortions due to topographic changes in seafloor MT data. The method first corrects the sea effect in observed MT impedance, and then inverts corrected responses in a model space without the sea. Due to mutual coupling between sea and subsurface structure, the correction and inversion steps are iterated until changes in each result become negligible. The method is validated for 1-D and 2-D structure using synthetic MT data produced by 3-D forward modeling including surrounding seas. In all cases, the method closely recovers the given structure after a few iterations. To test the applicability of the proposed method to field data, we generate synthetic MT data for the Jeju Island whose 1-D conductivity structure is well known, using 3-D forward modeling. The distortions due to the surrounding sea start to appear below the frequency about 1 Hz, and are relatively severe in the electrical field perpendicular to the coastline because of the location of the observation sites. The proposed method successfully eliminates the sea effect after three iterations, and both 1-D and 2-D inversion of corrected responses closely recover the given subsurface structure of the Jeju Island model.

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

In order to reduce leakage from a reservoir, a large amount of cement milk (grout) was injected from boreholes drilled around the shores of the reservoir, and monitored to establish the infiltration of cement milk into the bedrock under the reservoir. From laboratory tests using rock core samples, it was revealed that the resistivity of cement milk is much lower than that of the groundwater at this location. Therefore, it was expected that the resistivity of the zones filled with cement milk would be significantly reduced. Geophysical surveys are expected to be suitable methods to check the effectiveness of grouting in improving the water-retaining performance of a reservoir. DC electrical surveys (seven in total) and two Controlled Source Audio-frequency Magneto-Telluric (CSAMT) surveys were conducted along survey lines in the reservoir to monitor the infiltration of cement milk during the grouting. Extremely low resistivity zones ($10\;{\Omega}m$ or less) were observed in resistivity sections obtained by 2D inversion. The zones are inferred to be fractured zones filled with cement milk. In sections showing the rate of change of resistivity, three zones that showed significant change showed gradual expansion to deeper parts as the grouting progressed. These zones correspond to highly permeable zones detected by Lugeon tests at grout boreholes. We have confirmed that it is possible to measure the resistivity change by DC electrical and CSAMT surveys from the surface of the reservoir. It seems that such monitoring results could be reflected in future grouting plans.

• Geophysics and Geophysical Exploration
• /
• v.2 no.2
• /
• pp.96-103
• /
• 1999

Most of saline water inousions have been diagnosed by geophysical or geochemical approach independently. The objective of this study is to provide the effective method to detect the saline water intrusion on the ground water in the vicinity of seashore using these two methods. Schulumberger sounding, frequency domain electromagnetic sounding and geochemical analysis of ground water were carried out to explore saline water intrusion. Schulumberger sounding was implemented in dry surface condition before irrigation water was introduced into the field, while electromagnetic sounding was carried out in wet ground condition after the irrigation. The purpose of duplicated measurements on the equivalent spot at different times was to investigate the variation of anomaly zone depending on the amount of ground water. It was possible to discriminate the anomalous zone due to high water saturation from the low electric resistivity zone by high salt concentrations through this way. For the verification of the geophysical result, the ground water samples in the study area were collected and analysed at the 23 points near the measuring spots. The groundwater at the spot nearest to the sea water intrusion identified by geophysical method indicates higher salinity than the standard limit concentration for agricultural irrigation water (250 mg/1). Isotope analysis of $D({^2}H)$ vs. is ${^18}O$ and PCA analysis were used to discriminate the anthropogenic pollution from those of high salinity from sea water intrusion.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.164-175
• /
• 2011

Using natural electromagnetic (EM) fields at low frequencies, magnetotelluric (MT) surveys can investigate conductivity structures of the deep subsurface and thus are used to explore geothermal energy resources and investigate proper sites for not only geological $CO_2$ sequestration but also enhanced geothermal system (EGS). Moreover, marine MT data can be used for better interpretation of marine controlled-source EM data. In the interpretation of MT data, MT modeling schemes are important. This study improves a three dimensional (3D) MT modeling algorithm which uses edge finite elements. The algorithm computes magnetic fields by solving an integral form of Faraday's law of induction based on a finite difference (FD) strategy. However, the FD strategy limits the algorithm in computing vertical magnetic fields for a topographic model. The improved algorithm solves the differential form of Faraday's law of induction by making derivatives of electric fields, which are represented as a sum of basis functions multiplied by corresponding weightings. In numerical tests, vertical magnetic fields for topographic models using the improved algorithm overcome the limitation of the old algorithm. This study recomputes induction vectors and tippers for a 3D hill and valley model which were used for computation of the responses using the old algorithm.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.156-162
• /
• 2005

Particularly in research related to seawater intrusion the change of fluid electrical conductivity is one of major concerns, and effective monitoring can help to optimize a water pumping performance in coastal areas. Special considerations should be given to the mounting of sensors at proper depth during the monitoring design since the vertical distribution of fluid electrical conductivity is sensitive to the characteristics of seawater intrusion zone. This tells us the multi-channel electrical conductivity monitoring is of paramount consequence. It, however, is a rare event when this approach becomes routinely available in that commonly used commercial stand-alone type sensors are very expensive and inadequate for a long term monitoring of electrical conductivity or water level due to their restricted storage and difficulty of real-time control. For this reason, we have developed a real-time monitoring system that could meet these requirements. This system is user friendly, cost-effective, and easy to control measurement parameters - sampling interval, acquisition range, and others. And this devised system has been utilized for the electrical conductivity monitoring in boreholes, Yeonggwang-gun, Korea. Monitoring has been consecutively executed for 24 hours, and the responses of electrical conductivity at some channels have been regularly increased or decreased while pumping up water. It, with well logging data implemented before/after pumping water, verifies that electrical conductivity changes in the specified depths originate from fluid movements through sand layer or permeable fractured rock. Eventually, the multi-channel electrical conductivity monitoring system makes an effective key to secure groundwater resources in coastal areas.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.137-144
• /
• 2005

In this study, we offered the results of geophysical and geochemical survey on the municipal waste disposal area to delineate the size and extent of leachate contamination. Preliminary to intensive geochemical investigation, we performed two geophysical methods to characterize the survey area. Electromagnetic (EM) and magnetic method were used far site investigation. From the EM method, we can get the information of soil conductivity directly related to the leachate of the contaminations and from magnetic anomalies we can find the boundary of landfill which is not identified on the surface due to soil capping. The results of geophysical survey were well matched to those of geochemical method carried out inside and near the landfill. Electric conductivity (EC) of the groundwater sampled from low resistivity anomaly region of EM result was higher than background value and the border estimated from the magnetic survey showed good agreement with that estimated from the soil gas detection survey.

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