• Geophysics and Geophysical Exploration
• /
• 제20권4호
• /
• pp.226-231
• /
• 2017

Since microgravity survey aims to delineate subsurface density structures in small scale, it requires inversion method, which is able to resolve small scale structures. It can be achieved by adopting a stabilizing functional which separates density boundary distinctly, which is different concept from general inversion routines. We composed Matlab-based interactive two-dimensional microgravity data inversion package containing several kinds of inversion routines with different stabilizing functional, for handling various geologic conditions and survey purposes. Different kinds of inversion routines in the package were verified and examined with representative synthetic data sets generated by numerical modeling. In addition, we applied the developed package to a real microgravity survey data.

• The Journal of Engineering Geology
• /
• 제14권1호
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.459-463
• /
• 2003

Waveform inversion requires extracting a reliable low frequency content of seismic data for estimating of the low wave number velocity model. The low frequency content of the seismic data is usually discarded or neglected because of the band-limited response of the source and the receivers. In this study, however small the spectral of the low frequency seismic data is, we assume that it is possible to extract a reliable phase information of the low frequency from the seismic data and use it in waveform inversion. To this end, we exploit the frequency domain finite element modeling and source-receiver reciprocity to calculate the $Frech\`{e}t$ derivative of the phase of the seismic data with respect to the earth model parameter such as velocity, and then apply a damped least squares method to invert the phase of the seismic data. Through numerical example, we will attempt to demonstrate the feasibility of our method in estimating the correct velocity model for prestack depth migration.

• Proceedings of the KSEEG Conference
• /
• 대한자원환경지질학회 2002년도 춘계 공동학술발표회
• /
• pp.159-162
• /
• 2002

We tried to obtain an initial velocity model for prestack depth migration via waveform inversion. For application of any field data we chose a smooth background layered velocity model (v=v0 + k x z) as an initial velocity model. Newton type waveform inversion needs to invert huge Hessian matrix. In order to compute full Hessian matrix arising from full aperture data and full illumination zone, we meet insurmountable difficulties of paying astronomical computing cost. For the layered media, approximate Hessian emerging from single shot aperture data can be used repeatedly for split spread source configuration. In our work of using this Hessian characteristic of layered media we attempted to obtain the approximate velocity model as close as possible to the true velocity model in first iteration.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• 제24권4호
• /
• pp.533-542
• /
• 2000

Analytical experiments to determine the line-of-sight temperature distribution is conducted by using spectral radiation intensities. For this study, fourteen narrow bands of $25cm^{-1}$ interval in $CO_2\;4.3{\mu}m$ band ($2,050cm^{-1}$ to $2375cm^{-1}$) are selected. The applied system is a one-dimensional gas slab filled with 100% $CO_2$ gas at 1 atm. Two types of temperature profile are tested; parabolic and boundary layer types. Three kinds of radiation calculation are used in the iteration procedure for the temperature inversion; LBL(Line by Line), SNB(Statistical Narrow Band) and WNB(WSGGM. based Narrow Band) models. The LBL solution shows perfect agreement while some error of temperature prediction is caused by radiation modeling error when using SNB and WNB models. The inversion result shows that the WNB model may be used more accurately in spectral remote sensing techniques than the traditional SNB model.

• Journal of the Korean earth science society
• /
• 제23권1호
• /
• pp.15-29
• /
• 2002

Bayesian inversion for gravity and resistivity data was performed to investigate the cavity structure appearing as a lava tunnel in Cheju Island, Korea. Dipole-dipole DC resistivity data were proposed for a prior information of gravity data and we applied the geostatistical techniques such as kriging and simulation algorithms to provide a prior model information and covariance matrix in data domain. The inverted resistivity section gave the indicator variogram modeling for each threshold and it provided spatial uncertainty to give a prior PDF by sequential indicator simulations. We also presented a more objective way to make data covariance matrix that reflects the state of the achieved field data by geostatistical technique, cross-validation. Then Gaussian approximation was adopted for the inference of characteristics of the marginal distributions of model parameters and Broyden update for simple calculation of sensitivity matrix and SVD was applied. Generally cavity investigation by geophysical exploration is difficult and success is hard to be achieved. However, this exotic multiple interpretations showed remarkable improvement and stability for interpretation when compared to data-fit alone results, and suggested the possibility of diverse application for Bayesian inversion in geophysical inverse problem.

• Geomechanics and Engineering
• /
• 제35권2호
• /
• pp.109-119
• /
• 2023

Ground subsidence in urban areas due to excessive development and degraded underground facilities is a serious problem. Geophysical surveys have been conducted to estimate the distribution and scale of cavities and subsidence. In this study, electrical resistivity tomography (ERT) was performed near an area of road subsidence in an urban area. The subsidence arose due to groundwater leakage that carried soil into a neighboring excavation site. The ERT survey line was located between the main subsidence area and an excavation site. Because ERT data are affected by rapid topographic changes and surrounding structures, the influence of the excavation site on the data was analyzed through field-scale numerical modeling. The effect of an excavation should be considered when interpreting ERT data because it can lead to wrong anomalous results. A method for performing 2D inversion after correcting resistivity data for the effect of the excavation site was proposed. This method was initially tested using a field-scale numerical model that included the excavation site and subsurface anomaly, which was a loosened zone, and was then applied to field data. In addition, ERT data were interpreted using an existing in-house 3D algorithm, which considered the effect of excavation sites. The inversion results demonstrated that conductive anomalies in the loosened zone were greater compared to the inversion that did not consider the effects of excavation.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• 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.

• Geophysics and Geophysical Exploration
• /
• 제22권2호
• /
• pp.49-55
• /
• 2019

In this study, we presented a reference-shot subset method for stable convergence of full waveform inversion using a cyclic-shot subsampling technique. Full waveform inversion needs repetitive modeling of wave propagation and thus its calculation time increases as the number of sources increases. In order to reduce the computation time, we can use a cyclic-shot subsampling method; however, it makes the cost function oscillate in the early stage of the inversion and causes a problem in applying the convergence criteria. We introduced a method in which the cost function is calculated using a fixed reference-shot subset while updating the model parameters using the cyclic-shot subsampling method. Through the examples of full waveform inversion using the Marmousi velocity model, we confirmed that the convergence of cost function becomes stable even under the cyclic-shot subsampling method if using a reference-shot subset.

• Geophysics and Geophysical Exploration
• /
• 제24권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.