• Journal of the Korean earth science society
• /
• v.30 no.7
• /
• pp.834-844
• /
• 2009

Two-dimensional magnetotelluric (MT) inversion algorithm using minimum support (MS) stabilizer functional was implemented in this study to enhance the contrast of inverted images. For this implementation, this study derived a formula in discrete form for creeping model updates in the least-squares linearized inversion. A spatially varying regularization parameter determination algorithm, which is known as ACB (Active Constraint Balancing), was also adopted to stabilize the inversion process when using MS stabilizer as a model constraint. Inversion experiments for a simple isolated body model show well the feature of MS stabilizer in concentrating the anomalous body compared with the second-order derivative model constraint. This study also compared MS stabilizer and the second-order derivative model constraints for a model having multiple anomalous bodies to show the applicability of the algorithm into field data.

• The Journal of Engineering Geology
• /
• v.28 no.3
• /
• pp.367-377
• /
• 2018

This study introduces a DC resistivity inversion method that incorporates structural and inequality constraints to enhance the resolution of resistivity inversions, and presents sample inversion results with these constraints. In the constrained inversions, a base model is constructed from a layered model through interpretation of other geophysical data. Inversion tests establish that both the structural and inequality constraints produce better resistivity models than the unconstrained inversion. However, the inequality inversion not only reproduces the exact layered structure of the background, it reproduces conductive anomalies at a depth of ~ 10 m when an inexact base model of electrical resistivity is used.

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

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.8
• /
• pp.32-41
• /
• 1999

In this paper, an inversion method using chirp signals and two near field receivers is proposed. Inversion problems can be formulated into the probabilistic models composed of signals, a forward model and noise. Forward model to simulate chirp signals is chosen to be the source-wavelet-convolution planewave modeling method. The solution of the inversion problem is defined by a posteriori pdf. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the ranges for a priori uniform distribution is based. The genetic algorithm can be applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L₂norm of the difference between measured and modeled signals. The observed signals can be separated into a set of two signals reflected from the upper and lower boundaries of a sediment. The separation of signals and successive applications of the genetic algorithm optimization process reduce the search space, therefore improving the inversion results. Not only the marginal pdf but also the statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm. The examples applied here show that, for synthetic data with noise, it is possible to carry out an inversion for sedimentary layers using the proposed inversion method.

• Proceedings of the KSEEG Conference
• /
• 1999.04a
• /
• pp.248-250
• /
• 1999

• Geophysics and Geophysical Exploration
• /
• v.22 no.4
• /
• pp.202-209
• /
• 2019

In this study, an acoustic full-waveform inversion using Adam optimizer was proposed. The steepest descent method, which is commonly used for the optimization of seismic waveform inversion, is fast and easy to apply, but the inverse problem does not converge correctly. Various optimization methods suggested as alternative solutions require large calculation time though they were much more accurate than the steepest descent method. The Adam optimizer is widely used in deep learning for the optimization of learning model. It is considered as one of the most effective optimization method for diverse models. Thus, we proposed seismic full-waveform inversion algorithm using the Adam optimizer for fast and accurate convergence. To prove the performance of the suggested inversion algorithm, we compared the updated P-wave velocity model obtained using the Adam optimizer with the inversion results from the steepest descent method. As a result, we confirmed that the proposed algorithm can provide fast error convergence and precise inversion results.

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

It is difficult to obtain high-resolution images by 3D gravity inversion, because the problem is extremely underdetermined - there are too many model parameters. In order to reduce the number of model parameters we propose a 3D gravity inversion scheme utilising Euler deconvolution as a priori information. The essential point of this scheme is the reduction of the nonuniqueness of solutions by restricting the inversion space with the help of Euler deconvolution. We carry out a systematic exploration of the growing body process, but only in the restricted space within a certain radius of the Euler solutions. We have tested our method with synthetic gravity data, and also applied it to a real dataset, to delineate underground cavities in a limestone area. We found that we obtained a more reasonable subsurface density image by means of this combination between the Euler solution and the inversion process.

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.131-138
• /
• 2013

We developed a method for inverting magnetic data to image 2D susceptibility models. The major difficulty in the inversion of the potential data is the nonuniqueness. Furthermore, generally the number of inversion blocks are greater than the number of the magnetic data available, and thus the magnetic inversion leads to under-determined problem, which aggravates the nonuniqueness. When the magnetic data were inverted by the general least-squares method, the anomalous susceptibility would be concentrated near the surface in the inverted section. To overcome this nonuniqueness problem, we propose a new resolution model constraint that is calculated from the parameter resolution. The model constraint imposes large penalty on the model parameter with good resolution, on the other hand small penalty on the model parameter with poor resolution. Thus, the deep-seated model parameter, generally having poor resolution, can be effectively resolved. The developed inversion algorithm is applied to the inversion of the synthetic data for typical models of magnetic anomalies and is tested on real airborne data obtained at the Okcheon belt of Korea.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.86-92
• /
• 2006

The full-waveform inversion algorithm using normalised seismic wavefields can avoid potential inversion errors due to source estimation required in conventional full-waveform inversion methods. In this paper, we have modified the inversion scheme to install a weighted smoothness constraint for better resolution, and to implement a staged approach using normalised wavefields in order of increasing frequency instead of inverting all frequency components simultaneously. The newly developed scheme is verified by using a simple two-dimensional fault model. One of the most significant improvements is based on introducing weights in model parameters, which can be derived from integrated sensitivities. The model-parameter weighting matrix is effective in selectively relaxing the smoothness constraint and in reducing artefacts in the reconstructed image. Simultaneous multiple-frequency inversion can almost be replicated by multiple single-frequency inversions. In particular, consecutively ordered single-frequency inversion, in which lower frequencies are used first, is useful for computation efficiency.

• Geophysics and Geophysical Exploration
• /
• v.21 no.1
• /
• pp.15-25
• /
• 2018

Seismic inversion is a high-resolution tool to delineate the subsurface structures which may contain oil or gas. On the other hand, marine controlled-source electromagnetic (mCSEM) inversion can be a direct tool to indicate hydrocarbon. Thus, the joint inversion using both EM and seismic data together not only reduces the uncertainties but also takes advantage of both data simultaneously. In this paper, we have developed a simultaneous joint inversion approach for the direct estimation of reservoir petrophysical parameters, by linking electromagnetic and seismic data through rock physics model. A cross-gradient constraint is used to enhance the resolution of the inversion image and the maximum likelihood principle is applied to the relative weighting factor which controls the balance between two disparate data. By applying the developed algorithm to the synthetic model simulating the simplified gas field, we could confirm that the high-resolution images of petrophysical parameters can be obtained. However, from the other test using the synthetic model simulating an anticline reservoir, we noticed that the joint inversion produced different images depending on the model constraint used. Therefore, we modified the algorithm which has different model weighting matrix depending on the type of model parameters. Smoothness constraint and Marquardt-Levenberg constraint were applied to the water-saturation and porosity, respectively. When the improved algorithm is applied to the anticline model again, reliable porosity and water-saturation of reservoir were obtained. The inversion results indicate that the developed joint inversion algorithm can be contributed to the calculation of the accurate oil and gas reserves directly.