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

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

Tunnel detection is known to be one of the challenging problems in geophysical society. Among various geophysical methods, we tried to examine the applicability of electrical resistivity tomography (ERT) method to detect empty tunnel. In this study, we analyzed the ERT data acquired at the test site for tunnel detection. The inversion results have shown reasonable image of the tunnel although the resolution is quite poor. Moreover, we could obtain the three-dimensional attitude of tunnel through 3-D ERT imaging. Therefore, we expect that ERT can make contribution to the tunnel detection problem and further research effort such as fusion of geophysical methods will provide more reliable tunnel detection capability.

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

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.207-215
• /
• 2015

Electrical resistivity surveys were conducted at two subsidence areas near and at limestone mine sites, respectively, in order to estimate their causes of subsidence and the regions of potential occurrence. In addition, the linkages of mine development with these subsidences were investigated by the rock engineering analysis. Two study areas have different geological setting. One study site is the land subsidence area, which contains clay and sandy soil near the limestone mine, The other study site is the land subsidence area located just above the mine, which is expected to be relevant to the limestone mine. As results of two-dimensional (2D) electrical resistivity surveys at the sites 1 and 2, low resistivity zones, which are 70 ~ 120 ohm-m and 20 ~ 50 ohm-m, respectively, were found under the subsidence zones. For the study site 1, the possibility of subsidence was confirmed by using three-dimensional (3D) inversion performed with 2D resistivity profiles. For the study site 2, the cause of the subsidence and the possibility of subsidence occurrence were confirmed by rock engineering computation with regard to measurement line 7 in which low resistivity accompanied by subsidence area was observed.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.44-49
• /
• 2006

In this study, the resolution capabilities of electrical resistivity tomography (ERT) in the monitoring of $CO_2$ injection are investigated. The pole-pole and bipole-bipole electrode configuration types are used between two uncased boreholes straddling the $CO_2$ plume. Forward responses for an initial pre-injection model and three models for subsequent stages of $CO_2$ injection are calculated for the two different electrode configuration types, noise is added and the theoretical data are inverted with both L1- and L2-norm optimisation. The results show that $CO_2$ volumes over a certain threshold can be detected with confidence. The L1-norm proved superior to the L2-norm in most instances. Normalisation of the inverted models with the pre-injection inverse model gives good images of the regions of changing resistivity, and an integrated measure of the total change in resistivity proves to be a valid measure of the total injected volume.

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

• 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.17 no.2
• /
• pp.95-103
• /
• 2014

This paper presents a practical inversion method for recovering a three-dimensional (3D) resistivity model and static shifts simultaneously. Although this method is based on a Gauss-Newton approach that requires a sensitivity matrix, the computer time can be greatly reduced by implementing a simple and effective procedure for updating the sensitivity matrix using the Broyden's algorithm. In this research, we examine the approximate inversion procedure and the weighting factor ${\beta}$ for static shifts through inversion experiments using synthetic MT data. In methods using the full sensitivity matrix constructed only once in the iteration process, a procedure using the full sensitivity in the earlier stage is useful to produce the smallest rms data misfit. The choice of ${\beta}$ is not critical below some threshold value. Synthetic examples demonstrate that the method proposed in this paper is effective in reconstructing a 3D resistivity structure from static-shifted MT data.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.41-46
• /
• 2008

After choosing twenty sluice culverts located at national rivers by detailed appearance evaluation, hydraulic response test and acoustic televiewer in addition to surface geophysical explorations were carried out to investigate the cavities which exist at the boundary of sluice culvert and ground. Tapping, GPR, 3-D resistivity surveys and acoustic televiewer imagery obtained successful results in limited conditions according to characteristics of the sites. On the other hand, hydraulic response test successfully detected cavity existing at the lower part of sluice culvert and its continuity. However, hydraulic response test can not supply the quantitative informations on the scale and the shape of cavity. Thus it is very important for accurate investigation to make comprehensive interpretation after supplying various surveys proper to site characteristics.