• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.268-274
• /
• 2007

Traditional two-dimensional (2D) interpretation of magnetotelluric (MT) data utilizes only transverse magnetic (TM)-mode data, because 2D inversion of transverse electric (TE)-mode data results in spurious features when 3D structures exist in the subsurface. The application of a 3D inversion algorithm to a single MT profile can reduce contamination due to off-profile anomalies and help us to incorporate TE-mode data in the interpretation. In this study, we conduct 2D and 3D inversions of MT data observed along two lines in Jeju Island. First, we invert apparent resistivities and phases in the TM and TE modes separately. Then, we perform 2D joint inversion of both TM- and TE-mode data and 3D inversion of both Zxy- and Zyx-mode data corresponding to TE- and TM-mode data in 2D. The resistivity images derived from all four data show that the geoelectrical structure in Jeju Island is a three-layered earth with the resistive-conductive-resistive stratigraphy within a depth of 5 km. The 3D inversion does not produce clear anomalies in the reconstructed profile image, while all of 2D do. This attributed to the possibility that 2D inversion results are distorted by exiting off-profile 3D anomalies in Jeju. With 3D inversion of 2D profile MT data, we can deduce more reliable results that are not seriously distorted by off-profile 3D anomalies.

• Geophysics and Geophysical Exploration
• /
• v.16 no.1
• /
• pp.27-35
• /
• 2013

In this study, we developed the trial product on 3D magnetic inversion tentatively named 'KMag3D'. Also, we briefly introduced its own function and graphic user interface on which especially focused through the development in the form of user manual. KMag3D is consisted of two fundamental frame for the 3D magnetic inversion. First, algebraic reconstruction technique was selected as a 3D inversion algorithm instead of least square method conventionally used in various magnetic inversion. By comparison, it was turned out that algebraic reconstruction algorithm was more effective and economic than that of least squares in aspect of both computation time and memory. Second, for the effective determination of the 3D initial and a-priori information model required in the execution of our algorithm, we proposed the practical technique based on the assemblage of 2D forward modeling and inversion results for individual user-selected 2D profiles. And in succession, initial and a-priori information model were constructed by appropriate interpolation along the strke direction. From this, we concluded that our technique is both suitable and very practical for the application of 3D magentic inversion problem.

• 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.20 no.1
• /
• pp.33-42
• /
• 2017

Recently, the grounded electrical-source airborne transient electromagnetic (GREATEM) system with high power source was introduced to achieve deeper investigation depth and to overcome high noise level. Although the GREATEM is a transient electromagnetic system using a long grounded wire as the transmitter, GREATEM data have been interpreted with 1D earth models because 2D or 3D modeling and inversion of vast airborne data are complicated and expensive to calculate. Generally, 1D inversion is subsequently applied to every survey point and combining 1D images together forms the stitched conductivity-depth image. However, the stitched models often result in abrupt variations in neighboring models. To overcome this problem, laterally constrained inversion (LCI) has been developed in inversion of ATEM data, which can yield layered sections with lateral smooth transitions. In this study, we analysed the GREATEM data through 1D numerical modeling for a curved grounded wire source. Furthermore, we developed a laterally constrained inversion scheme for continuous GREATEM data based on a layered earth model. All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. Applying the developed LCI technique to the GREATEM data, it was confirmed that the laterally constrained inversion can provide laterally smooth model sections that reflect the layering of the survey area effectively.

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

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

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

• Economic and Environmental Geology
• /
• v.17 no.3
• /
• pp.163-169
• /
• 1984

A linear inversion of heat flow values using heat production data with reliable value is studied in this work. To evaluate 2-D problem, a thin vertical sheet model is considered. Making use of a relation based on potential theory, a new relation between $q_{rad}$ and $A_0$ is derived. The forward calculations with noise and without noise are shown. The inversion of random search is comparable to that of ridge regression method. The agreements between the computed and best fit after inversion suggest the importance of random search method in the inversion technique.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.171-179
• /
• 2021

We outline a process for estimating Cole-Cole parameters from time-domain induced polarization (IP) data. The IP transients are all inverted to 2D Cole-Cole earth models that include resistivity, chargeability, relaxation time, and the frequency exponent. Our inversion algorithm consists of two stages. We first convert the measured voltage decay curves into time series of current-on time apparent resistivity to circumvent the negative chargeability problem. As a first step, a 4D inversion recovers the resistivity model at each time channel that increases monotonically with time. The desired intrinsic Cole-Cole parameters are then recovered by inverting the resistivity time series of each inversion block. In the second step, the Cole-Cole parameters can be estimated readily by setting the initial model close to the true value through a grid search method. Finally, through inversion procedures applied to synthetic data sets, we demonstrate that our algorithm can image the Cole-Cole earth models effectively.