• Journal of the Korean Geophysical Society
• /
• v.4 no.1
• /
• pp.35-46
• /
• 2001

한반도 중부에 위치한 대전 지진관측소(TJN) 하부의 세부 지각구조를 밝혀내기 위하여 수신함수를 이용한 선형화된 역산(linearized inversion) 방법을 적용하였다. 본 방법의 비단일해(nonuniqueness)와 초기 모델 의존성의 문제를 해결하기 위해 근사 초기 속도 모델로부터 72개의 서로 다른 초기 모델을 구하여 역산을 수행한 후 결과모델들의 평균 속도 모델을 제시하는 방법을 사용하였다. 역산 결과 총 72개의 모델 중 뚜렷한 지각-맨틀 경계를 보이는 43개의 모델만이 조건에 만족하는 결과를 나타내었다. 모든 모델에서 속도 구조는 전체적으로 깊이에 따라 속도의 불연속면이나 급격한 증가없이 연속적인 변화를 하며, 모호면의 깊이는 30～32.5 km의 범위로 나타났다. 평균적인 하부 지각의 속도는 6.5 km/s, 상부 맨틀의 속도는 7.8 km/s로 뚜렷한 속도 변화를 보였다. 결과 모델 군은 중부지각(mid-crust)에서의 속도를 기준으로 약한 저속도층을 나타내는 군과 상대적으로 일정한 속도를 가지는 군으로 구분되었다. 단지 지진파형의 비교만으로 두 모델군 중 합당한 모델군의 선택은 불가능하였다. 따라서 수신 함수를 이용하여 연구 지역의 신뢰할 만한 지각 구조를 구하기 위해서는 그 지역에 대한 지질학적, 지구물리학적 추가정보와의 동반 해석이 요구된다.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.220-226
• /
• 2011

In the elastic wave equations, both horizontal and vertical displacements are defined. Since we can measure both the horizontal and vertical displacements in field acquisition, these displacements compose a displacement vector. In this study, we propose a frequency-domain elastic waveform inversion technique taking advantage of the magnitudes of displacement vectors to define objective function. When we apply this displacement-vector objective function to the frequency-domain waveform inversion, the inversion process naturally incorporates the back-propagation algorithm. Through the inversion examples with the Marmousi model and the SEG/EAGE salt model, we could note that the RMS error of the solution obtained by our algorithm decreased more stably than that of the conventional method. Particularly, the density of the Marmousi model and the low-velocity sub-salt zone of the SEG/EAGE salt model were successfully recovered. Since the gradient direction obtained from the proposed objective function is numerically unstable, we need additional study to stabilize the gradient direction. In order to perform the waveform inversion using the displacementvector objective function, it is necessary to acquire multi-component data. Hence, more rigorous study should be continued for the multi-component land acquisition or OBC (Ocean Bottom Cable) multi-component survey.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.214-219
• /
• 2011

We perform the frequency-domain waveform inversion based on the residual-selection strategy. In the residual-selection strategy, we classify time-domain residual wavefields into several groups according to the order of absolute amplitudes. Because the residual wavefields are normalized after regularization of the gradient directions within each group, the residual-selection strategy plays a role in enhancing the small-amplitude wavefields, which contributes to improving the deep parts of inverted subsurface images. After classifying residuals in the time domain, they are transformed to the frequency domain. Waveform inversion is performed in the frequency domain using the back-propagation technique which has been popularly used in reverse-time migration. The residual-selection strategy is applied to the SEG/EAGE salt and IFP Marmousi models. Numerical results show that the residual-selection strategy yields better results than the conventional frequency-domain waveform inversion.

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

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

In the geophysical monitoring to understand the change of subsurface material properties with time, the time-invariant static subsurface model is commonly adopted to reconstruct a time-lapse image. This assumption of static model, however, can be invalid particularly when fluid migrates very quickly in highly permeable medium in the brine injection experiment. In such case, the resultant subsurface images may be severely distorted. In order to alleviate this problem, we develop a new least-squares inversion algorithm under the assumption that the subsurface model will change continuously in time. Instead of sampling a time-space model into numerous space models with a regular time interval, a few reference models in space domain at different times pre-selected are used to describe the subsurface structure continuously changing in time; the material property at a certain space coordinate are assumed to change linearly in time. Consequently, finding a space-time model can be simplified into obtaining several reference space models. In order to stabilize iterative inversion and to calculate meaningful subsurface images varying with time, the regularization along time axis is introduced assuming that the subsurface model will not change significantly during the data acquisition. The performance of the proposed algorithm is demonstrated by the numerical experiments using the synthetic data of crosshole dc resistivity tomography.

• Proceedings of the Acoustical Society of Korea Conference
• /
• autumn
• /
• pp.385-388
• /
• 2004

본 논문에서는 천해 해양환경에서 이동 음원과 수직선배열을 이용한 실 해상 실험 자료를 이용하여 지음향 역산 및 정합장처리를 수행하였다. 협대역 신호를 이용하여 음원의 이동에 따른 전달손실을 계산하였고, 이를 음파전달모델을 이용한 결과와 비교하는 방식으로 역산의 목적함수를 구성함으로써 해저면 지음향 인자들의 평균적인 값을 역산할 수 있었다. 역산 결과로 얻은 지음향 인자들을 사용하여 정합장처리를 통해 음원의 이동경로를 추적함으로써 역산 결과의 타당성을 검증하였고, 이를 기존 지질 자료를 이용한 정합장처리 결과와 비교하였다.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.180-186
• /
• 2003

We have developed an one-dimensional (ID) inversion program that can invert multiple frequency small-loop EM data from horizontal coplanar (HCP) and vertical coplanar (VCP) configurations. The inverse problem is solved using least-squares method with active constraint balancing (ACB) method and Jacobian matrix is calculated analytically. Tests using synthetic data from simple ID models indicate that conductivity and depth of each layer can be estimated properly when both real and imaginary data are used together.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.38-49
• /
• 2020

Full-waveform inversion (FWI) is an optimization process of fitting observed and modeled data to reconstruct high-resolution subsurface physical models. In acoustic FWI (AFWI), pressure data acquired using a marine streamer has mainly been used to reconstruct the subsurface P-wave velocity models. With recent advances in marine seismic-acquisition techniques, acquiring multi-component data in marine environments have become increasingly common. Thus, AFWI strategies must be developed to effectively use marine multi-component data. Herein, we proposed an AFWI strategy using horizontal and vertical particle-acceleration data. By analyzing the modeled acoustic data and conducting sensitivity kernel analysis, we first investigated the characteristics of each data component using AFWI. Common-shot gathers show that direct, diving, and reflection waves appearing in the pressure data are separated in each component of the particle-acceleration data. Sensitivity kernel analyses show that the horizontal particle-acceleration wavefields typically contribute to the recovery of the long-wavelength structures in the shallow part of the model, and the vertical particle-acceleration wavefields are generally required to reconstruct long- and short-wavelength structures in the deep parts and over the whole area of a given model. Finally, we present a sequential-inversion strategy for using the particle-acceleration wavefields. We believe that this approach can be used to reconstruct a reasonable P-wave velocity model, even when the pressure data is not available.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.326-334
• /
• 2008

The resistivity method has been used to image the electrical properties of the subsurface. Especially, this method has become suitable for monitoring since data could be rapidly and automatically acquired. In this study, we developed a time-lapse inversion algorithm for the interpretation of resistivity monitoring data. The developed inversion algorithm imposes a big penalty on the model parameter with small change, while a minimal penalty on the model parameter with large change compared to the reference model. Through the numerical experiments, we can ensure that the time-lapse inversion result shows more accurate and focused image where model parameters have changed. Also, applying the timelapse inversion method to the leakage detection of an embankment dam, we can confirm that there are three major leakage zones, but they have not changed over time.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.164-170
• /
• 2021

Frequency-domain and time-domain induced polarization methods can provide spectral information about subsurface media. Analysis of spectral characteristics has been studied mainly in the frequency-domain, however, time-domain induced polarization research has recently become popular. In this study, assuming a homogeneous half-space model, an inversion method was developed to extract Cole-Cole parameters from the measured secondary potential or electrical resistivity. Since the Cole-Cole parameters of chargeability, time constant, and frequency index are not independent of each other, various problems, such as slow convergence rate, initial model problem, local minimum problem, and divergence, frequently occur when conventional nonlinear inversion is applied. In this study, we developed an effective inversion method using the initial model close to the true model by introducing a grid search method. Finally, the validity of the developed inversion method was verified using inversion experiments.