• Geophysics and Geophysical Exploration
• /
• v.22 no.3
• /
• pp.160-167
• /
• 2019

Resistivity inversion result may be distorted if the seepage line fluctuation within central core with the change of reservoir water level as well as the conductivity of the reservoir water is not taken into consideration because the reservoir dike is composed of three-dimensional (3D) resistivity structure. Consequently, to accurately analyze the resistivity changes inside the reservoir dike according to the change of reservoir water level, 3D electrical resistivity modeling for the 2D survey line considering topography and physical properties of dam components was carried out. In addition, 2D inversion was performed with the simulated 2D resistivity data for a given 3D model in order to compare it with the inversion result of real field data. For 283 reservoirs in Korea, 2D inversion results for the simulated 2D data and field 2D resistivity data were compared. Finally, the reservoirs with an inversion ratio of 50% or less were selected as reservoirs that require further precise investigation.

• Geophysics and Geophysical Exploration
• /
• v.18 no.1
• /
• pp.14-20
• /
• 2015

Three-dimensional (3D) resistivity method is an effective tool in the engineering site survey because it can provide a 3D resistivity distribution of the site. In this study, we tried to find out faults, fractures and coal seams that can cause the collapse of the tunnel. We carried out 2D resistivity survey along 5 parallel lines and 11 cross lines and merged all the apparent resistivity data for 3D inversion. Finally, from the 3D resistivity image and drilling data we presented the 3D distribution of faults, fractures and coal seams that are considered the main cause of the tunnel collapse.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.55-58
• /
• 2005

Low resistivity zone is observed at the lower part of a CFRD (Concrete Face Rockfill Dam). Generally, CFRD tends not to have any saturated zone within the body, but the result of resistivity survey shows that it is possible for the dam to be saturated under 20m depth with water. The level of reservoir was under 10m from the crest. We suspect that this result may come from the wrong 2D inversion process ignoring the 3D geometry of dams. For the analysis of possibility of distortion by different geometry, we perform the 3D forward modeling for the dam and apply the 2D inversion process. And then we check the point of traditional interpretation of resistivity data. By the analysis, it is found that the result of 2D inversion process of 3D geometry of dams, seems to have deep relation with the reservoir level, and the complex 3D structure hide some internal electrical anomaly of dams from resistivity information.

• Journal of the Korean Geophysical Society
• /
• v.8 no.4
• /
• pp.211-214
• /
• 2005

Low resistivity zone is observed at the lower part of a CFRD (Concrete Face Rockfill Dam). Generally, CFRD tends not to have any saturated zone within the body, but the result of resistivity survey shows that it is possible for the dam to be saturated under 20m depth with water. The level of reservoir was under 10 m from the crest. We suspect that this result may come from the wrong 2D inversion process ignoring the 3D geometry of dams. For the analysis of possibility of distortion by different geometry, we perform the 3D forward modeling for the dam and apply the 2D inversion process. And then we check the point of traditional interpretation of resistivity data. By the analysis, it is found that the result of 2D inversion process of 3D geometry of dams, seems to have deep relation with the reservoir level, and the complex 3D structure hide some internal electrical anomaly of dams from resistivity information.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.199-206
• /
• 2003

Characteristics of the static shift are discussed by comparing the three-dimensional MT inversion with/without static shift parameterization. The galvanic distortion by small-scale shallow feature often leads severe distortion in inverted resistivity structures. The new inversion algorithm is applied to four numerical data sets contaminated by different amount of static shift. In real field data interpretations, we generally do not have any a-priori information about how much the data contains the static shift. In this study, we developed an algorithm for finding both Lagrangian multiplier for smoothness and the trade-off parameter for static shift, simultaneously in 3-D MT inversion. Applications of this inversion routine for the numerical data sets showed quite reasonable estimation of static shift parameters without any a-priori information. The inversion scheme is successfully applied to all the four data sets, even when the static shift does not obey the Gaussian distribution. Allowing the static shift parameters have non-zero degree of freedom to the inversion, we could get more accurate block resistivities as well as static shifts in the data. When inversion does not consider the static shift as inversion parameters (conventional MT inversion), the block resistivities on the surface are modified considerably to match possible static shift. The inhomogeneous blocks on the surface can generate the static shift at low frequencies. By those mechanisms, the conventional 3-D MT inversion can reconstruct the resistivity structures to some extent in the deeper parts even when moderate static shifts are in the data. As frequency increased, however, the galvanic distortion is not frequency independent any more, and thus the conventional inversion failed to fit the apparent resistivity and phase, especially when strong static shift is added. Even in such case, however, reasonable estimation of block resistivity as well as static shift parameters were obtained by 3-D MT inversion with static shift parameterization.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.104-113
• /
• 2014

Electric resistivity monitoring was applied to evaluate the soft ground in reclaimed land in order to figure out the applicability of physical prospecting. For this, electrical resistivity monitoring data were acquired for total three months and analyzed those data with independent inversion, time-lapse inversion, and 4D inversion methods. The result was compared for various inversion methods so as to figure out what showed the soft soil most properly. Moreover, drilling and CPT(Cone Penetration Test) data were also used in order to find out if each of those inversion methods could distinguish either bed rock or the soft soil clearly. And according to the result, time-lapse inversion showed less inversion artifacts than independent inversion, so it could indicate the soft soil better. If data gained for a longer period than three months are used, 4D inversion has been found to be a more efficient analysis method than the time-lapse inversion method. Electrical resistivity monitoring on the soft soil has been found to be a useful method that can analyze the spatio-temporal electric state of the ground serially.

• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.139-149
• /
• 2018

The resistivity monitoring is a practical method to resolve changes in resistivity of underground structures over time. With the advance of sophisticated automatic data acquisition system and rapid data communication technology, resistivity monitoring has been widely applied to understand spatio-temporal changes of subsurface. In this study, a new 4D inversion algorithm is developed, which can effectively emphasize significant changes of underground resistivity with time. To overcome the overly smoothing problem in 4D inversion, the Lagrangian multipliers in the space-domain and time-domain are determined automatically so that the proportion of the model constraints to the misfit roughness remains constant throughout entire inversion process. Furthermore, a focusing model constraint is added to emphasize significant spatio-temporal changes. The performance of the developed algorithm is demonstrated by the numerical experiments using the synthetic data set for a time-lapse model.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.167-172
• /
• 2007

The frequency-domain small-loop electromagnetic (EM) instruments are increasingly used for shallow environmental and geotechnical surveys because of their portability and speed. However, it is well known that the data quality is generally so poor that quantitative interpretation of the data is not justified in many cases. We present an inversion method that allows the correction for the calibration errors and also constructs multidimensional resistivity models. The key point in this method is that the data are collected at least at two different heights. The forward modeling used in the inversion is based on an efficient 3-D finite-difference method, and its solution was checked against 2-D finite-element solution. The synthetic and real data examples demonstrate that the joint inversion recovers reliable resistivity models from multi-frequency data severely contaminated by the calibration errors.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.