• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.54-63
• /
• 2015

Petrophysical parameters such as porosity and fluid saturation which provide useful information for reservoir characterization could be estimated by rock physics model (RPM) using seismic velocity and resistivity. Therefore, accurate P-wave velocity and resistivity information have to be obtained for successful estimation of the petrophysical parameters. Compared with the individual inversion of electromagnetic (EM) or seismic data, the joint inversion using both EM and seismic data together can reduce the uncertainty and gives the opportunity to use the advantages of each data. Thus, more reliable petrophysical properties could be estimated through the joint inversion. In this paper, for the successful estimation of petrophysical parameters, we proposed an effective method which applies a grid-search method to find the porosity and fluid saturation. The relations of porosity and fluid saturation with P-wave velocity and resistivity were expressed by using RPM and the improved resistivity distribution used to this study was obtained by joint inversion of seismic and EM data. When the proposed method was applied to the synthetic data which were simulated for subsea reservoir exploration, reliable petrophysical parameters were obtained. The results indicate that the proposed method can be applied for detecting a reservoir and calculating the accurate oil and gas reserves.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.185-190
• /
• 2007

In this study, we developed the 2.5D EM modeling and inversion algorithm for cross-hole source and receiver geometry. Considering the cross-hole environment, we use a VMD (vertical magnetic dipole) as a source and vertical magnetic fields as a measuring data. Developed inversion algorithm is tested for the isolated block model which has a conductive and a resistivity anomaly respectively. For the conductive anomaly, its size and resistivity are inverted well on the inversion results, while for the resistive anomaly, the location of anomalous block is shown on the inverted section, but its values are far from the exact value. Furthermore, artificial conductive anomalies are shown around the resistive anomalous zone. If we consider the inversion artifact shown in the test inversion of restive block, it is almost impossible to image the resistive zone. However, the main target of EM tomography in the engineering problem is conductive target such as fault zone, and contaminated zone etc., EM tomography algorithm can be used for detecting the anomalous zone.

• Geophysics and Geophysical Exploration
• /
• v.22 no.4
• /
• pp.186-194
• /
• 2019

Deterministic optimization, commonly used to find the geophysical inverse solutions, have its limitation that it cannot find the proper solution since it might converge into the local minimum. One of the solutions to this problem is to use global optimization based on a stochastic approach, among which a large number of particle swarm optimization (PSO) applications have been introduced. In this paper, I developed a geophysical inversion algorithm applying PSO method for the layered-earth resistivity inversion of the small-loop electromagnetic (EM) survey data and carried out numerical inversion experiments on synthetic datasets. From the results, it is confirmed that the PSO inversion algorithm could increase the inversion success rate even when attempting the inversion of small-loop EM survey data from which it might be difficult to find a best solution by applying the Gauss-Newton inversion algorithm.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.340-343
• /
• 2003

This study presents a practical procedure for the Bayesian inversion of geophysical data by Markov chain Monte Carlo (MCMC) sampling and geostatistics. We have applied geostatistical techniques for the acquisition of prior model information, and then the MCMC method was adopted to infer the characteristics of the marginal distributions of model parameters. For the Bayesian inversion of dipole-dipole array resistivity data, we have used the indicator kriging and simulation techniques to generate cumulative density functions from Schlumberger array resistivity data and well logging data, and obtained prior information by cokriging and simulations from covariogram models. The indicator approach makes it possible to incorporate non-parametric information into the probabilistic density function. We have also adopted the MCMC approach, based on Gibbs sampling, to examine the characteristics of a posteriori probability density function and the marginal distribution of each parameter. This approach provides an effective way to treat Bayesian inversion of geophysical data and reduce the non-uniqueness by incorporating various prior information.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.129-136
• /
• 2005

The damped least-squares inversion has become a most popular method in finding the solution in geophysical problems. Generally, the least-squares inversion is to minimize the object function which consists of data misfits and model constraints. Although both the data misfit and the model constraint take an important part in the least-squares inversion, most of the studies are concentrated on what kind of model constraint is imposed and how to select an optimum regularization parameter. Despite that each datum is recommended to be weighted according to its uncertainty or error in the data acquisition, the uncertainty is usually not available. Thus, the data weighting matrix is inevitably regarded as the identity matrix in the inversion. We present a new inversion scheme, in which the data weighting matrix is automatically obtained from the analysis of the data resolution matrix and its spread function. This approach, named 'extended active constraint balancing (EACB)', assigns a great weighting on the datum having a high resolution and vice versa. We demonstrate that by applying EACB to a two-dimensional resistivity tomography problem, the EACB approach helps to enhance both the resolution and the stability of the inversion process.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.217-224
• /
• 2013

We developed a time-lapse inversion using new cross-model constraints based on change ratio and resolution of model parameters. The cross-model constraint based on change ratio imposes the same penalty on the model parameters with equal change ratio. This constraint can emphasize the model parameters with significant change regardless of their increase or decrease. The resolution cross-model constraint imposes a small penalty on the model parameters with poor resolution, but a large penalty on the model parameters with good resolution. Thus, the model parameter with poor resolution can be effectively identified in the inversion result if they are significantly changed with time. Through the numerical tests for 3D resistivity monitoring data sets, the performance of these two cross-model constraints was confirmed. Finally, for the safety estimation of a sea dyke, we applied the developed time-lapse inversion to the 3D resistivity monitoring data that were acquired at a sea dike located in western coastal area of Korea. The result of time-lapse inversion suggested that there were no significant changes at the sea dike during the monitoring period.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.401-408
• /
• 2005

Recently shield TBM tunnellings are being applied to subway construction in Korean cities. Generally these kinds of tunnellings have the problems in the stability of ground such as subsidence because urban subway is constructed in the shallow depth. A sinkhole occurred on the road just above the tunnel during tunneling in Kwangju, so a survey for upper layer of the tunnel was needed. But conventional Ground Probing Radar can't be applicable due to the presence of steel-mesh screen in the shield segment, so no existent geophysical method is applicable in this site. Because the outer surface of each shield segment is electrically insulated, dipole-dipole resistivity method which is popular in engineering site investigation, was tried to this survey for the first time. Specially manufactured flexible ring-type electrodes were installed into the grouting holes at an interval of 2.4 m on the ceiling. The K-Ohm II system which has been developed by KIGAM and tested successfully in many sites, was used in this site. The system consists of 1000Volt-1Ampere constant-current transmitter, optically isolated 24 bit sigma-delta A/D conversion receiver - maximum 12 channel simultaneous measurements, and graphical automatic acquisition software for easy data quality check in real time. Borehole camera logging with circular white LED lighting was also done to investigate the state of the layer. Measured resistivity data lack of some stations due to failing opening lids of holes, shows general high-low trend well. The dipole-dipole resistivity inversion results discriminate (1) one approximately 4 meter diameter cavity (grouted but incompletely hardened, so low resistivity - less than $30{\Omega}m$), (2) weak zone (100-200${\Omega}m$), and (3) hard zone (high resistivity - more than 1000${\Omega}m$) very well for the distance of 320 meters. The 2-D inversion neglects slight absolute 3-D effect, but we can get satisfactory and useful information. Acquired resistivity section and video tapes by borehole camera logging will be reserved and reused if some problem occurs in this site in the future.

• Economic and Environmental Geology
• /
• v.14 no.3
• /
• pp.193-201
• /
• 1981

The problem of automatic inversion of apparent resistivity sounding curves resulting from horizontally layered earth models is solved using the least-squares technique. This method, which makes use of damped least-squares algorithm in conjunction with digital filtering technique, is found to be speedier and more accurate than the conventional curve-matching method. Four sounding curves were chosen to test the inversion scheme. The analysis of the theoretical sounding data associated with a three-layer model illustrates clear advantages over the conventional curve-matching method. The usefulness of the inversion method is also shown when applied to the actual field data. It was found that the best fit earth models coincide with the subsurface structures confirmed by drilling.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.236-249
• /
• 2002

We have extended the three-dimensional (3-D) resistivity imaging algorithm to cover the 3-D resistivity tomography problem, where resistivity data are acquired using electrodes installed in several boreholes as well as at the earth surface. The imaging algorithm consists of the 3-D finite element forward modeling and least-squares inversion scheme, where the ACB (Active Constraint Balancing) is adopted to enhance the resolving power of the inversion. Sensitivity analysis with numerical verifications shows that 3-D resistivity tomography is a very appealing method and can be used to get 3-D attitude of subsurface structures with very high-resolution. Moreover, we could accurately handle the topography effect, which could cause artifacts in the resistivity tomography. In the application of 3-D resistivity tomography to the real field data set acquired at the quarry mine, we could derive a very reasonable and accurate image of the subsurface.