• Economic and Environmental Geology
• /
• v.18 no.3
• /
• pp.217-224
• /
• 1985

A Householder's transformation is applied to the resistivity inversion problem. The conventional resistivity inversion method is sometimes numerically unstable in interpreting a resistivity sounding data because it usually solves the normal equation derived from an observation equation. The resistivity inversion method using Householder's transformation solves the observation equation directly, so that it is numerically more stable than the conventional method. A theoretical, ill-conditioned Schlumberger sounding data was chosen to test the inversion scheme with Householder's transformation.

• 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.24 no.4
• /
• pp.149-157
• /
• 2021

Given that induced polarization (IP) and direct current (DC) resistivity surveys are similar in terms of data acquisition, most DC resistivity systems are equipped with a time-domain IP data acquisition function. In addition, the time-domain IP data include the DC resistivity values. As such, IP and DC resistivity data are intimately linked, and the inversion of IP data is a two-step process based on DC resistivity inversions. Nevertheless, IP surveys are rarely applied, in contrast to DC resistivity surveys, as proper inversion software is unavailable. In this study, through numerical modeling and inversion experiments, we analyze the problems with the conventional inverse mapping technique used to invert time-domain IP data. Furthermore, we propose a modified inverse mapping technique that can effectively suppress inversion artifacts. The performance of the technique is confirmed through inversions applied to synthetic IP data.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.11 no.4
• /
• pp.361-371
• /
• 2009

Resistivity tomography is very effective geophysical method to find out the resistivity distribution and its change in time around a tunnel. Thus, the resistivity tomogram can provide helpful information which is necessary for the effective maintenance of the tunnel. However, an air filled tunnel severely distorts tomography data, especially when the current or potential electrode is placed near the tunnel. Moreover, the distortion can often lead to misinterpretation of tomography monitoring data. To solve these problem, we developed a resistivity modeling and time-lapse inversion program which include a tunnel. In this study, using the developed program we assured that the inversion including a tunnel gives much more accurate image around a tunnel, compared with the conventional tomogram where the tunnel is not included. We also confirmed that the time-lapse inversion of resistivity monitoring data defines well resistivity changed areas around a tunnel in time.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.12a
• /
• pp.31-42
• /
• 2007

Dc resistivity monitoring has been increasingly used in order to understand the changes of subsurface conditions in terms of conductivity. The commonly adopted interpretation approach which separately inverts time-lapse data may generate inversion artifacts due to measurement error. Eventually the contaminated error amplifies the artifacts when reconstructing the difference images to quantitatively estimate the change of ground condition. In order to alleviate the problems, we defined the subsurface structure as four dimensional (4-D) space-time model and developed 4-D inversion algorithm which can calculate the reasonable subsurface structure continuously changing in time even when the material properties change during data measurements. In this paper, we discussed two case histories of resistivity monitoring to study the ground condition change when the properties of the subsurface material were artificially altered by injecting conductive materials into the ground: (1) dye tracer experiment to study the applicability of electrical resistivity tomography to monitoring of water movement in soil profile and (2) the evaluation of cement grouting performed to reinforce the ground. Through these two case histories, we demonstrated that the 4-D resistivity imaging technique is very powerful to precisely delineate the change of ground condition. Particularly owing to the 4-D inversion algorithm, we were able to reconstruct the history of the change of subsurface material property.

• The Journal of Engineering Geology
• /
• v.28 no.3
• /
• pp.367-377
• /
• 2018

This study introduces a DC resistivity inversion method that incorporates structural and inequality constraints to enhance the resolution of resistivity inversions, and presents sample inversion results with these constraints. In the constrained inversions, a base model is constructed from a layered model through interpretation of other geophysical data. Inversion tests establish that both the structural and inequality constraints produce better resistivity models than the unconstrained inversion. However, the inequality inversion not only reproduces the exact layered structure of the background, it reproduces conductive anomalies at a depth of ~ 10 m when an inexact base model of electrical resistivity is used.

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

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.180-193
• /
• 2021

The complex resistivity method is an exploration technique that can obtain various characteristic information of underground media by measuring resistivity and phase in the frequency domain, and its utilization has recently increased. In this paper, a three-dimensional inversion algorithm for the CR data was developed to increase the utilization of this method. The Poisson equation, which can be applied when the electromagnetic coupling effect is ignored, was applied to the modeling, and the inversion algorithm was developed by modifying the existing algorithm by adopting comlex variables. In order to increase the stability of the inversion, a technique was introduced to automatically adjust the Lagrangian multiplier according to the ratio of the error vector and the model update vector. Furthermore, to compensate for the loss of data due to noisy phase data, a two-step inversion method that conducts inversion iterations using only resistivity data in the beginning and both of resistivity and phase data in the second half was developed. As a result of the experiment for the synthetic data, stable inversion results were obtained, and the validity to real data was also confirmed by applying the developed 3D inversion algorithm to the analysis of field data acquired near a hydrothermal mine.