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

• 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.17 no.2
• /
• pp.95-103
• /
• 2014

This paper presents a practical inversion method for recovering a three-dimensional (3D) resistivity model and static shifts simultaneously. Although this method is based on a Gauss-Newton approach that requires a sensitivity matrix, the computer time can be greatly reduced by implementing a simple and effective procedure for updating the sensitivity matrix using the Broyden's algorithm. In this research, we examine the approximate inversion procedure and the weighting factor ${\beta}$ for static shifts through inversion experiments using synthetic MT data. In methods using the full sensitivity matrix constructed only once in the iteration process, a procedure using the full sensitivity in the earlier stage is useful to produce the smallest rms data misfit. The choice of ${\beta}$ is not critical below some threshold value. Synthetic examples demonstrate that the method proposed in this paper is effective in reconstructing a 3D resistivity structure from static-shifted MT data.

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

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

The advent of towed geo-electrical array surveying on water and land has resulted in datasets of magnitude approaching that of airborne electromagnetic surveying and most suited to 1D inversion. Robustness and complete automation is essential if processing and reliable interpretation of such data is to be viable. Sharp boundaries such as river beds and the top of saline aquifers must be resolved so use of smoothness constraints must be minimised. Suitable inversion algorithms must intelligently handle low signal-to-noise ratio data if conductive basement, that attenuates signal, is not to be misrepresented. A noise-level aware inversion algorithm that operates with one elastic thickness layer per electrode configuration has been coded. The noise-level aware inversion identifies if conductive basement has attenuated signal levels so that they are below noise level, and models conductive basement where appropriate. Layers in the initial models are distributed to span the effective depths of each of the geo-electric array quadrupoles. The algorithm works optimally on data collected using geo-electric arrays with an approximately exponential distribution of quadrupole effective depths. Inversion of data from arrays with linear electrodes, used to reduce contact resistance, and capacitive-line antennae is plausible. This paper demonstrates the effectiveness of the algorithm using theoretical examples and an example from a salt interception scheme on the Murray River, Australia.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.268-274
• /
• 2007

Traditional two-dimensional (2D) interpretation of magnetotelluric (MT) data utilizes only transverse magnetic (TM)-mode data, because 2D inversion of transverse electric (TE)-mode data results in spurious features when 3D structures exist in the subsurface. The application of a 3D inversion algorithm to a single MT profile can reduce contamination due to off-profile anomalies and help us to incorporate TE-mode data in the interpretation. In this study, we conduct 2D and 3D inversions of MT data observed along two lines in Jeju Island. First, we invert apparent resistivities and phases in the TM and TE modes separately. Then, we perform 2D joint inversion of both TM- and TE-mode data and 3D inversion of both Zxy- and Zyx-mode data corresponding to TE- and TM-mode data in 2D. The resistivity images derived from all four data show that the geoelectrical structure in Jeju Island is a three-layered earth with the resistive-conductive-resistive stratigraphy within a depth of 5 km. The 3D inversion does not produce clear anomalies in the reconstructed profile image, while all of 2D do. This attributed to the possibility that 2D inversion results are distorted by exiting off-profile 3D anomalies in Jeju. With 3D inversion of 2D profile MT data, we can deduce more reliable results that are not seriously distorted by off-profile 3D anomalies.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.55-58
• /
• 2005

Low resistivity zone is observed at the lower part of a CFRD (Concrete Face Rockfill Dam). Generally, CFRD tends not to have any saturated zone within the body, but the result of resistivity survey shows that it is possible for the dam to be saturated under 20m depth with water. The level of reservoir was under 10m from the crest. We suspect that this result may come from the wrong 2D inversion process ignoring the 3D geometry of dams. For the analysis of possibility of distortion by different geometry, we perform the 3D forward modeling for the dam and apply the 2D inversion process. And then we check the point of traditional interpretation of resistivity data. By the analysis, it is found that the result of 2D inversion process of 3D geometry of dams, seems to have deep relation with the reservoir level, and the complex 3D structure hide some internal electrical anomaly of dams from resistivity information.

• Journal of the Korean Geophysical Society
• /
• v.8 no.4
• /
• pp.211-214
• /
• 2005

Low resistivity zone is observed at the lower part of a CFRD (Concrete Face Rockfill Dam). Generally, CFRD tends not to have any saturated zone within the body, but the result of resistivity survey shows that it is possible for the dam to be saturated under 20m depth with water. The level of reservoir was under 10 m from the crest. We suspect that this result may come from the wrong 2D inversion process ignoring the 3D geometry of dams. For the analysis of possibility of distortion by different geometry, we perform the 3D forward modeling for the dam and apply the 2D inversion process. And then we check the point of traditional interpretation of resistivity data. By the analysis, it is found that the result of 2D inversion process of 3D geometry of dams, seems to have deep relation with the reservoir level, and the complex 3D structure hide some internal electrical anomaly of dams from resistivity information.

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

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.250-262
• /
• 2008

In order to assess the feasibility of the dipole-dipole electric method to the investigation of metallic ore deposit, both field data simulation and inversion are carried out for several simplified ore deposit models. Our interest is in a vein-type model, because most of the ore deposits (more than 70%) exist in a vein type in Korea. Based on the fact that the width of the vein-type ore deposits ranges from tens of centimeters to 2m, we change the width and the material property of the vein, and we use 40m-electrode spacing for our test. For the vein-type model with too small width, the low resistivity zone is not detected, even though the resistivity of the vein amounts to 1/300 of that of the surrounding rock. Considering a wide electrode interval and cell size used in the inversion, it is natural that the size of the low resistivity zone is overestimated. We also perform field data simulation and inversion for a vein-type model with surrounding hydrothermal alteration zones, which is a typical structure in an epithermal ore deposits. In the model, the material properties are assumed on the basis of resistivity values directly observed in a mine originated from an epithermal ore deposits. From this simulation, we can also note that the high resistivity value of the vein does not affect the results when the width of the vein is narrow. This indicates that our main target should be surrounding hydrothermal alteration zones rather than veins in field survey. From these results, we can summarize that when the vein is placed at the deep part and the difference of resistivity values between the vein and the surrounding rock is not large enough, we cannot detect low resistivity zone and interpret the subsurface structures incorrectly using the electric method performed at the surface. Although this work is a little simple, it can be used as references for field survey design and field data Interpretation. If we perform field data simulation and inversion for a number of models and provide some references, they will be helpful in real field survey and interpretation.