• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.191-201
• /
• 1999

The interpretation of resistivity data has, so far, mainly been made under the assumption that the earth is of relatively simple structure and then using one or two-dimensional inversion scheme. Since real earth structure and topography are fully three-dimensional and very complicated In nature, however, such assumptions often lead to misinterpretation of the earth structures. In such situations, three-dimensional inversion is probably the only way to get correct image of the earth. In this study, we have developed a three-dimensional inversion code using the finite element solution for the forward problem. The forward modeling algorithm simulates the real field situation with irregular topography. The inverse problem is solved iteratively using the least-squares method with smoothness constraint. Our inversion scheme employs ACB (Active Constraint Balancing) to enhance the resolving power of the inversion. Including Irregular surface topography in the inversion, we can accurately define the earth structures without artifact in the numerical tests. We could get reasonable image of earth structure by Inverting the real field data sets taken over highway bridge construction site.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.139-148
• /
• 2021

Negative apparent chargeability data in time-domain induced polarization (IP) survey are very often considered as noise and rejected before the inversion stage. Negative IP data can, however, occur naturally as a consequence of the distribution of chargeable zones in the ground. In some cases, the negative apparent chargeability values may account for most of the data measured. Negative IP data are caused by the geometry of chargeable zones and electrode positions. Negative apparent chargeability data appear most frequently when a dipole-dipole array is used. In this study, the effect of negative apparent chargeability data on inversion results is analyzed through the numerical 2D time-domain IP modeling and nonlinear inversion. The results demonstrate clearly that negative apparent chargeability data have to be included in the inversion as they contain important information on the distribution of subsurface chargeability.

• The Korean Journal of Petroleum Geology
• /
• v.13 no.1
• /
• pp.17-23
• /
• 2007

Gravity characteristics are investigated in the vicinity of the DASAN scientific station, located at the Svalbard Archipelago, the Arctic using ArcGP data. Boundary effects of free-air gravity anomalies, which appeared generally at the continental margin, are erased after Bouguer correction was applied. Complete Bouguer anomalies produced after terrain correction by GrOPO30 show that gravity anomalies increase from continent to marine. This phenomena seem to be related to the rise of Moho discontinuity. The cut-off frequency of 0.16 was decided after power spectrum analysis and the gravity anomalies were divided into two parts. Residual anomalies in high frequency part show that characteristics of high values along the faults and of low values related to thick sediments in the continent. Characteristic is low values from basement subsidence of continental slope or thick sediments in the marine. The undulation of Moho discontinuity from 3-D inversion modeling show typical characteristics of continental margin that become higher from Svalbard archipelago to Knipovich ridge bordering Eurasian plate.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.83-86
• /
• 2006

From numerical modeling test $-{\rho}a$ by one steel casing borehole near resistivity survey line can be acquired. Negative apparent resistivities even in the flat area are surely subsurface information. Inversion technique for those need to be developed in the near future.

• Geophysics and Geophysical Exploration
• /
• v.1 no.3
• /
• pp.155-160
• /
• 1998

EM tomography technique has been developed. The algorithm used the extended Born approximations for forward modeling and reconstructed a conductivity image by a smoothness constraint least squares inversion method. Observed data, the vertical components of secondary magnetic fields, were simulated with the 3-D integral equation code. The results showed that the location of anomalous body could be imaged very well, but conductivity of the body was lower than real one and the vertical resolution was much higher than the horizontal resolution.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.308-313
• /
• 2007

The small-loop electromagnetic technique has been used successfully for many geophysical qualitative investigations, particularly for shallow engineering and environmental surveys. Recently, various geophysical imaging methods based on numerical modeling and inversion have been tried in order to get more quantitative subsurface structure. However, conventional 2.5D small loop EM modeling takes a lot of time because responses should be calculated for several wave numbers and transformed into space domain. In this study, we developed a 2.5D HCP small loop EM modeling algorithm using extended Born approximation, which does not require transformation. Also, we checked its validity by comparison with other numerical results.

• Journal of the Korean Geophysical Society
• /
• v.10 no.1
• /
• pp.13-25
• /
• 2007

We investigated the undulation of Moho depth and the crustal structure of the continental margin in the East Sea along the Korea Peninsula from inversion and modelling using potential data and previous seismic results. Free-air gravity anomalies generally reflect topography effect. Bouguer gravity anomalies increase toward the Ulleung Basin, indicating that Moho depth is shallower under the Ulleung Basin. Positive magnetic anomalies exist along the continental margin and decrease toward the Ulleung Basin. In analytic signal, the small anomaly in the Hupo Bank infers that the Hupo Bank is uplifted by igneous intrusion and the strong anomaly on the continental slope denotes existence of SDR(seaward dipping reflectors), which are in accordance with the location of SDR detected in previous seismic studies. The inversion result of Bouguer gravity anomaly and the 2-dimensional gravity modelling indicate that the undulation of Moho depth shallows from the continental shelf toward the Ulleung Basin. This is in good agreement with the Moho depth calculated by the previous seismic velocity model using ocean bottom seismometer(OBS). The 2-dimensional gravity modelling infers magmatic underplating zone under the lower continental crust on the continental margin of the East Sea, indicating the possible rifiting of the continental margin.

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.119-130
• /
• 2013

We developed a method for inverting magnetic data to recover the 3D susceptibility models. The major difficulty in the inversion of the potential data is the non-uniqueness and the vast computing time. The insufficient number of data compared with that of inversion blocks intensifies the non-uniqueness problem. Furthermore, there is poor depth resolution inherent in magnetic data. To overcome this non-uniqueness problem, we propose a resolution model constraint that imposes large penalty on the model parameter with good resolution; on the other hand, small penalty on the model parameter with poor resolution. Using this model constraint, the model parameter with a poor resolution can be effectively resolved. Moreover, the wavelet transform and parallel solving were introduced to save the computing time. Through the wavelet transform, a large system matrix was transformed to a sparse matrix and solved by a parallel linear equation solver. This procedure is able to enormously save the computing time for the 3D inversion of magnetic data. The developed inversion algorithm is applied to the inversion of the synthetic data for typical models of magnetic anomalies and real airborne data obtained at the Geumsan area of Korea.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.19-28
• /
• 2007

In this study, we propose a joint inversion method, using genetic algorithms, to estimate an S-wave velocity structure for deep sedimentary layers from receiver functions and surface-wave phase velocity observed at several sites. The method takes layer continuity over a target area into consideration by assuming that each layer has uniform physical properties, especially an S-wave velocity, at all the sites in a target area in order to invert datasets acquired at different sites simultaneously. Numerical experiments with synthetic data indicate that the proposed method is effective in reducing uncertainty in deep structure parameters when modelling only surface-wave dispersion data over a limited period range. We then apply the method to receiver functions derived from earthquake records at one site and two datasets of Rayleigh-wave phase velocity obtained from microtremor array surveys performed in central Tokyo, Japan. The estimated subsurface structure is in good agreement with the results of previous seismic refraction surveys and deep borehole data. We also conclude that the proposed method can provide a more accurate and reliable model than individual inversions of either receiver function data only or surface-wave dispersion data only.