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

• Kyungpook Mathematical Journal
• /
• v.52 no.4
• /
• pp.383-395
• /
• 2012

In this paper, we use a computational algorithm for the inversion of the one and two-dimensional $\mathcal{L}_2$-transform based on the Bromwich's integral and the Fourier series. The new inversion formula can evaluate the inverse of the $\mathcal{L}_2$-transform with considerable accuracy over a wide range of values of the independent variable and can be devised for the functions which are not Laplace transformable and have damping motion in small interval near origin.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.291-298
• /
• 2002

We have developed a two-dimensional (2-D) magnetotelluric (MT) inversion algorithm, which can include topographic effects in inversion. We use the finite element method (FEM) to incorporate topography into forward calculation. Topography is implemented simply by moving nodes of rectangular elements in z-direction according to the elevation of air-earth interface. In the inversion process, we adopt a spatially variable Lagrangian multiplier algorithm in the smoothness-constrained least-squares inversion. The inversion algorithm developed in this study reconstructs subsurface resistivity structure quite well when topography variation exists. Also, it turns out to be effective in both resolution and stability from a model study and field data application.

• Proceedings of the KSEG Conference
• /
• 2001.03a
• /
• pp.67-71
• /
• 2001

• Atmosphere
• /
• v.18 no.3
• /
• pp.207-224
• /
• 2008

Severe downslope windstorms observed at Gangneung, Korea in the springtime during the last 30 years are studied to understand their generation mechanisms. 92 severe wind cases are selected for which the maximum instantaneous wind speeds exceed two standard deviation of total mean plus ($18.7ms^{-1}$). They are categorized into the three mechanisms (hydraulic jump, partial reflection, and critical-level reflection) proposed in previous studies based on the flow condition, which is calculated using the wind and temperature profile observed at one upstream rawinsonde station, Osan. Among the three, partial reflection is found to be the most frequent mechanism for the last 30 years (1976 - 2005). To understand the role of inversion in generating severe downslope windstorms, horizontal velocity perturbation was calculated analytically for the atmosphere with an inversion layer. It turned out that the intensity of downslope wind was increased by inversion layer of specific heights, which are well matched with the observations. For better understanding the generation mechanisms, two-dimensional numerical simulations are conducted for the 92 severe wind cases using the ARPS model. In most simulations, surface wind speed exceeds the value of the severe-wind criterion, and each simulated case can be explained by its own generation mechanism. However, in most simulations, the simulated surface wind speed is larger than the observed, due to ignoring the flow-splitting effect in the two-dimensional framework.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.21-26
• /
• 2005

Two-dimensional (2-D) inversion of magnetotelluric (MT) data for two survey lines having south-north direction from Jeiu Island has been carried out. The 2-D models show a thick layer having around 10 ohm-m in the depth of a few hundred meters throughout the survey area, which can be considered as the unconsolidated sedimentary layer. And they also show a conductive anomaly at the central part of each survey lines. But unfortunately by now, we do not have any further information about the anomaly. Comparison of the 2-D inversion model using MT band only and that using both AMT and MT bands said that it is helpful for us to include AMT band as well as MT band in the inversion to interpret not only the shallow part but also the deep structures.

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

• Economic and Environmental Geology
• /
• v.21 no.1
• /
• pp.107-113
• /
• 1988

We present a precedure for interpreting dipole-dipole apparent resistivity data. The procedure is constructed by combining a forward two-dimensional finite element modeling and an inverse technique with Householder's transformation. In the interpretation, subsurface structure is divided into some blocks with constant resistivities. Our inversion technique is tested on synthetic and field data. We found that geologic constraint is required for successful interpretation.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.60-68
• /
• 2007

Conventional analysis of spectral induced polarization (SIP) data consists of measuring impedances over a range of frequencies, followed by spectral analysis to estimate spectral parameters. For the quantitative and accurate estimation of subsurface SIP parameter distribution, however, a sophisticated and stable inversion technique is required. In this study, we have developed a two-step inversion approach to obtain the two-dimensional distribution of SIP parameters. In the first inversion step, all the SIP data measured over a range of frequencies are simultaneously inverted, adopting cross regularisation of model complex resistivities at each frequency. The cross regularisation makes it possible to enhance the noise characteristics of the inversion by imposing a strong assumption, that complex resistivities should show similar characteristics over a range of frequencies. In numerical experiments, we could verify that our inversion approach successfully reduced inversion artefacts. As a second step, we have also developed an inversion algorithm to obtain SIP parameters based on the Cole-Cole model, in which frequency-dependent complex resistivities from the first step are inverted to obtain a two-dimensional distribution of SIP parameters. In numerical tests, the SIP parameter images showed a fairly good match with the exact model, which suggests that SIP imaging can provide a very useful subsurface image to complement resistivity.

• Economic and Environmental Geology
• /
• v.22 no.3
• /
• pp.267-276
• /
• 1989

This paper presents results of interpretaton of gravity data by iterative nonlinear inversion methods. The gravity data are obtained by a theoretical formula for two-dimensional 2-layer structure. Depths to the basement of the structure are determined from the gravity data by four interative inversion methods. The four inversion methods used here are the Gradient, Gauss-Newton, Newton-Raphson, and Full Newton methods. Inversions are performed by using different initial guesses of depth for the over-determined, even-determined, and under-determined cases. This study shows that the depth can be determined well by all of the methods and most efficiently by the Newton-Raphson method.