• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.198-200
• /
• 2002

경제적이며 정확한 3차원 전자탐사 모델링을 위해 위해 Habashy et al. (1993)에 의해 제안된 국소 비선형 근사(localized nonlinear approximation)를 이용하여 전자탐사 모델링 알고리듬을 개발하였다. 전자탐사 수치모델링시 많은 계산시간 및 기억용량을 필요로 하는 Green 텐서 적분을 정확하고 빠르게 계산하기 위해, 단일 미소체를 이용한 공간파수 영역에서의 Green 텐서 적분 알고리듬을 제안하였다. 더욱이 Green 텐서의 송수신 방향 및 상반성을 고려하여 각각의 미소체에 의한 전체 미소체에의 Green 텐서 적분을 한 개의 미소체에 의한 전체 미소체에의 Green 텐서 적분 값으로 구하게 하므로 매우 적은 기억용량 만으로 Green 텐서 적분 행렬을 구성할 수 있어, 역산법에 효과적으로 적용할 수 있다. 이 수치 모델링 알고리듬을 기본으로 하여 평활화 제한을 가한 최소자승 역산 알고리듬을 개발하였다. 이 역산 알고리듬을 지표 전자탐사 및 시추공-지표 전자탐사 등에 적용하여 PC에서도 빠르게 3차원 전자탐사 역산이 수행됨을 보였다.

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.207-212
• /
• 2004

This article reviews recent developments in three-dimensional (3-D) magntotelluric (MT) imaging. The inversion of MT data is fundamentally ill-posed, and therefore the resultant solution is non-unique. A regularizing scheme must be involved to reduce the non-uniqueness while retaining certain a priori information in the solution. The standard approach to nonlinear inversion in geophysis has been the Gauss-Newton method, which solves a sequence of linearized inverse problems. When running to convergence, the algorithm minimizes an objective function over the space of models and in the sense produces an optimal solution of the inverse problem. The general usefulness of iterative, linearized inversion algorithms, however is greatly limited in 3-D MT applications by the requirement of computing the Jacobian(partial derivative, sensitivity) matrix of the forward problem. The difficulty may be relaxed using conjugate gradients(CG) methods. A linear CG technique is used to solve each step of Gauss-Newton iterations incompletely, while the method of nonlinear CG is applied directly to the minimization of the objective function. These CG techniques replace computation of jacobian matrix and solution of a large linear system with computations equivalent to only three forward problems per inversion iteration. Consequently, the algorithms are efficient in computational speed and memory requirement, making 3-D inversion feasible.

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

Non-linear elastic wavefield inversion is a powerful method for estimating elastic parameters for physical constraints that determine subsurface rock and properties. Here, I introduce six elastic-wave velocity models by reconstructing elastic-wave velocity variations from real data and a 2D elastic-wave velocity model. Reflection seismic data information is often decoupled into short and long wavelength components. The local search method has difficulty in estimating the longer wavelength velocity if the starting model is far from the true model, and source frequencies are then changed from lower to higher bands (as in the 'frequency-cascade scheme') to estimate model elastic parameters. Elastic parameters are inverted at each inversion step ('simultaneous mode') with a starting model of linear P- and S-wave velocity trends with depth. Elastic parameters are also derived by inversion in three other modes - using a P- and S-wave velocity basis $('V_P\;V_S\;mode')$; P-impedance and Poisson's ratio basis $('I_P\;Poisson\;mode')$; and P- and S-impedance $('I_P\;I_S\;mode')$. Density values are updated at each elastic inversion step under three assumptions in each mode. By evaluating the accuracy of the inversion for each parameter set for elastic models, it can be concluded that there is no specific difference between the inversion results for the $V_P\;V_S$ mode and the $I_P$ Poisson mode. The same conclusion is expected for the $I_P\;I_S$ mode, too. This gives us a sound basis for full wavelength elastic wavefield inversion.

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

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.316-322
• /
• 2002

A computationally efficient inversion scheme has been developed using the extended Born or localized nonlinear approximation to analyze electromagnetic fields obtained in a single-hole environment. The medium is assumed to be cylindrically symmetric about the borehole, and to maintain the symmetry vertical magnetic dipole source is used throughout. The efficiency and robustness of an inversion scheme is very much dependent on the proper use of Lagrange multiplier, which is often provided manually to achieve desired convergence. In this study, an automatic Lagrange multiplier selection scheme has been developed to enhance the utility of the inversion scheme in handling field data. The inversion scheme has been tested using synthetic data to show its stability and effectiveness.

• Journal of the KSME
• /
• v.34 no.7
• /
• pp.510-516
• /
• 1994

탄성파의 속도는 파동 전파 매질의 밀도와 탄선계수에 달려 있으므로, 박막이 입혀진 기판에서 전파하는 표면파에 대해서 기판과 박막의 밀도와 탄성계수 및 박막의 두께 등을 알면 전파 속 도를 계산할 수 있다. 박막의 탄성계수를 모르는 경우에는 표면파 속도를 측정하여 역으로 탄 성계수를 산출할 수 있다. 이러한 역산과정에는 일반적인 비선형 방정식의 curve-fitting에 이용될 수 있는 simplex법이 효율적으로 활용된다. 이 글에서는 표면파 속도를 측정하고 그 데이터로 부터 역산하여 박막의 탄성계수를 구하는 원리와 과정을 설명한다.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2002.09a
• /
• pp.43-50
• /
• 2002

부지고유의 강지진동모사를 위해 한반도 남부의 지각감쇠(Q)에 대한 두 가지 특성을 비선형 역산을 통해 규명하였다. 한 특성은 한반도 남부의 주요 지체구조구 및 특정 주향각에 대한 Q의 이방성이고 다른 특성은 Q의 수평방향으로의 공간적 변화 특성이다. 0의 이방성은 주향에 평행한 Q와 이에 수직한 Q의 비로 정의되었다. 사용된 지진자료는 190개 지진에 대한 3,400개 기록의 푸리에스펙트럼으로서 기록의 파선 분포는 한반도 남부 대부분을 조밀하게 덮었다. 역산 방법은 Levenberg-Marquardt 비선형역산 방법을 보다 안정화시킨 수정된 Levenberg-Marquardt 방법이 사용되었으며 기존에 도출된 광역적인 지진동파라미터를 초기해로 활용하였다. 역산 결과 국내의 중요 지체구조구 중 소백산육괴와 경상분지는 서로 다른 방향의 강한 Q 이방성을 나타내었으며 Q의 공간적 분포는 지진발생위치와 매우 밀접한 상관관계를 보여주었다

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.2
• /
• pp.291-301
• /
• 1999

A method for the synthesis of one-dimensional nonlinear distribution function is presented for the desired inverse scattering pattern. This method is based on the inverse transform of the solution of the Riccati equation derived from one-dimensional inverse scattering problem. Since the solution is analogous to the array factor or normalized space factor in collinear array antenna, the synthesis method for field pattern is applied for the construction of the involved line-source nonlinear distribution function. The suggested method is carried out under the optimization process, and is numerically verified by synthesizing the dispersive transmission line profile within the specified frequency band and control of scattered field on resistive strip.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.167-173
• /
• 1999

Since the resolution of the 2-D hole-to-hole seismic traveltime tomography is affected by the limited ray transmission angle, various methods were used to improve the resolution. Linear traveltime interpolation(LTI) ray tracing method was chosen for forward-modeling method. Inversion results using the LTI method were compared with those using the other ray tracing methods. As an inversion algorithm, SIRT method was used. In the iterative non-linear inversion method, the cost of ray tracing is quite expensive. To reduce the cost, each raypath was stored and the inversion was performed from this information. Using the proposed method, fast convergence was achieved. Inversion results are likely to be affected by the initial velocity guess, especially when the ray transmission angle was limited. To provide a good initial guess for the inversion, generalized regression neural network(GRNN) method was used. When the transmitted raypath angle is not limited or the geological model is very complex, the inversion results are not affected by initial velocity model very much. Since the raypath angles, however, are limited in most geophysical tomographic problems, the enhancement of resolution in tomography can be achieved by providing a proper initial velocity model by another inversion algorithm such as GRNN.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.