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

• Journal of the Korean earth science society
• /
• v.30 no.6
• /
• pp.699-708
• /
• 2009

Recently as the interest in the development of domestic ore deposits has increased, we can easily find some studies on exploration geophysics-based ore-deposit survey in literature. Based on the fact that mineralized zone are generally more conductive than surrounding media, electrical resistivity survey among several geophysical surveys has been applied to investigate metallic ore deposits. Most of them are grounded on 2-D survey. However, 2-D inversion may lead to some misinterpretation for 3-D geological structures. In this study, we investigate the feasibility of the 3-D electrical resistivity survey to 3-D vein-type ore deposits. We first simulate 2-D dipole-dipole survey data for survey lines normal to the strike and 3-D pole-pole survey data, and then perform 3-D inversion. For 3-D ore-body structures, we assume a width-varying dyke, a wedge-shaped, and a fault model. The 3-D inversion results are compared to 2-D inversion results. By comparing 3-D inversion results for 2-D dipole-dipole survey data to 3-D inversion results for 3-D pole-pole survey data, we could note that the 2-D dipole-dipole survey data yield better inversion results than the 3-D pole-pole data, which is due to the main characteristic of the pole-pole array. From these results, we are convinced that if we have certain information on the direction of the strike, it would be desirable to apply 2-D dipole-diple survey for the survey lines normal to the strike. However, in most cases, we do not have any information on the direction of the strike, because we already developed the ore deposit with the outcrops and the remaining ore deposits are buried under the surface. In that case, performing 3-D pole-pole electrical resistivity survey would be a reasonable choice to obtain more accurate interpretation on ore body structure in spite of low resolution of pole-pole array.

• 한국지구물리탐사학회:학술대회논문집
• /
• 1999.08a
• /
• pp.21-36
• /
• 1999

Recently resistivity survey is widely used for site investigations in the field of civil engineering. Since such application area requires accurate interpretation tools especially in the area of complicated geology and rough terrain topography, we developed a three-dimensional (3-D) resistivity inversion code, which can reconstruct real earth structures. Furthermore, the inversion code gives resolution-enhanced images by applying the ACB(Active Constraint Balancing) method. With the help of this inversion code, 3-D resistivity survey is now used as new techniques for site investigations in civil engineering problem. By imaging the 3-D resistivity distribution, we could get useful informations such as depth distribution of basement rock, distribution of weak zone, fractures and cavities which is crucial to civil engineers.

• The Journal of Engineering Geology
• /
• v.17 no.2 s.52
• /
• pp.271-277
• /
• 2007

To find out the anomalous zone in cut slope composed of phyllite and shist, we performed resistivity tomography using a pole-dipole way. The electrical distribution that propagates from a current source in lower part of slope is measured by a potential electrode in upper part of slope. Apparent resistivity data are inverted with an iterative regularized inversion method to reconstruct 3D resistivity image. By comparing with the resistivity images in relation to each section, the images of anomalous zone correspond to their positions represented in cut slope. Therefore, the application of resistivity tomography in cut slope is useful to recognize the extension of anomalous zone.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.

• Geomechanics and Engineering
• /
• v.35 no.2
• /
• pp.109-119
• /
• 2023

Ground subsidence in urban areas due to excessive development and degraded underground facilities is a serious problem. Geophysical surveys have been conducted to estimate the distribution and scale of cavities and subsidence. In this study, electrical resistivity tomography (ERT) was performed near an area of road subsidence in an urban area. The subsidence arose due to groundwater leakage that carried soil into a neighboring excavation site. The ERT survey line was located between the main subsidence area and an excavation site. Because ERT data are affected by rapid topographic changes and surrounding structures, the influence of the excavation site on the data was analyzed through field-scale numerical modeling. The effect of an excavation should be considered when interpreting ERT data because it can lead to wrong anomalous results. A method for performing 2D inversion after correcting resistivity data for the effect of the excavation site was proposed. This method was initially tested using a field-scale numerical model that included the excavation site and subsurface anomaly, which was a loosened zone, and was then applied to field data. In addition, ERT data were interpreted using an existing in-house 3D algorithm, which considered the effect of excavation sites. The inversion results demonstrated that conductive anomalies in the loosened zone were greater compared to the inversion that did not consider the effects of excavation.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.236-249
• /
• 2002

We have extended the three-dimensional (3-D) resistivity imaging algorithm to cover the 3-D resistivity tomography problem, where resistivity data are acquired using electrodes installed in several boreholes as well as at the earth surface. The imaging algorithm consists of the 3-D finite element forward modeling and least-squares inversion scheme, where the ACB (Active Constraint Balancing) is adopted to enhance the resolving power of the inversion. Sensitivity analysis with numerical verifications shows that 3-D resistivity tomography is a very appealing method and can be used to get 3-D attitude of subsurface structures with very high-resolution. Moreover, we could accurately handle the topography effect, which could cause artifacts in the resistivity tomography. In the application of 3-D resistivity tomography to the real field data set acquired at the quarry mine, we could derive a very reasonable and accurate image of the subsurface.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.217-224
• /
• 2013

We developed a time-lapse inversion using new cross-model constraints based on change ratio and resolution of model parameters. The cross-model constraint based on change ratio imposes the same penalty on the model parameters with equal change ratio. This constraint can emphasize the model parameters with significant change regardless of their increase or decrease. The resolution cross-model constraint imposes a small penalty on the model parameters with poor resolution, but a large penalty on the model parameters with good resolution. Thus, the model parameter with poor resolution can be effectively identified in the inversion result if they are significantly changed with time. Through the numerical tests for 3D resistivity monitoring data sets, the performance of these two cross-model constraints was confirmed. Finally, for the safety estimation of a sea dyke, we applied the developed time-lapse inversion to the 3D resistivity monitoring data that were acquired at a sea dike located in western coastal area of Korea. The result of time-lapse inversion suggested that there were no significant changes at the sea dike during the monitoring period.

• Economic and Environmental Geology
• /
• v.39 no.4 s.179
• /
• pp.473-483
• /
• 2006

Although application of electrical methods in Korea began with observation of self potentials before World War II, the methods were developed slowly by the beginning of 1980's when a major burst of development activity took place. DC resistivity methods are applied in Korea more to geotechnical problems rather than to environmental ones unlike other developed countries. As with every other branch of technology, the evolving speed of the silicon chip and of streaming data to hard disk has revolutionized data collection and noise reduction processing. The last two decades saw major advances in data collection, processing, and interpretation of electrical data. Development of smooth-model two-dimensional (2D) resistivity inversion is one of the most visible changes to geophysical interpretation of the last 40 years and is now routinely applied to apparent resistivity data. The ability to represent resistivities in section rather than pseudosection view has revolutionized interpretation. Although calculation of sensitivities for general electromagnetic problems require numerous forward modelings, DC resistivity methods can enjoy computational efficiencies if sources and receivers occupy the same position, and previously intractable 3D inversion is now becoming available.