• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.173-182
• /
• 2009

In April 2008, KIGAM carried out an ocean-bottom seismometer (OBS) survey in the central Ulleung Basin where strong bottom simulating reflectors (BSRs) were revealed from previous surveys and some gas-hydrate samples were retrieved by direct sampling. The purpose of this survey is to estimate the velocity structure near the BSR in the gas hydrate prospect area using wide-angle seismic data recorded on the ocean-bottom seismometers. Along with the OBS survey, a 2-D seismic survey was performed whereby stratigraphic and preliminary velocity information was obtained. Two methods were applied to wide-angle data for estimating P wave velocity; one is velocity analysis in the $\tau$-p domain and the other is seismic traveltime inversion. A 1-D interval velocity profile was obtained by the first method, which was refined to layered velocity structure by the latter method. A layer stripping method was adopted for modeling and inversion. All velocity profiles at each OBS site clearly show velocity reversal at BSR depths due to the presence of gas hydrates. In addition, we could confirm high velocity in the column/chimney structure.

• Economic and Environmental Geology
• /
• v.28 no.4
• /
• pp.395-404
• /
• 1995

The geologic structure of the Cheju volcanic island has been investigated by analyzing the gravity and magnetic data. Bouguer gravity map shows apparent circular low anomalies at the central volacanic edifice, and the maximum difference of the anomaly values on the island appears to be 30 mgal. The subsurface structure of the island is modeled by three-dimensional depth inversion of gravity data by assuming the model consists of a stacked grid of rectangular prisms of volcanic rocks bounded below by basement rocks. The gravity modeling reveals that the interface between upper volvanic rocks and underlying basement warps downward under Mt. Halla with the maximum depth of 5 km. Magnetic data involve aeromagnetic and surface magnetic survey data. Both magnetic anomaly maps show characteristic features which resemble the typical pattern of total magnetic anomalies caused by a magnetic body magnetized in the direction of the geomagnetic field in the middle latitude region, though details of two maps are somewhat different. The reduced-to-pole magnetic anomaly maps reveal that main magnetic sources in the island are rift zones and the Halla volcanic edifice. The apparent magnetic boundaries inferred by the method of Cordell and Grauch (1985) are relatively well matched with known geologic boundaries such as that of Pyosunri basalt and Sihungri basalt which form the latest erupted masses. Inversion of aeromagnetic data was conducted with two variables: depth and susceptibility. The inversion results show high susceptibility bodies in rift zones along the long axis of the island, and at the central volcano. Depths to the basement are 1.5~3 km under the major axis, 1~1.5 km under the lava plateau and culminates at about 5 km under Mt. Halla. The prominent anomalies showing N-S trending appear in the eastern part of both gravity and magnetic maps. It is speculated that this trend may be associated with an undefined fault developed across the rift zones.

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

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.33-38
• /
• 2012

A resistivity method has been applied to wide range of engineering and environmental problems with the help of automatic and precise data acquisition. Thus, more accurate modeling and inversion of time-lapse monitoring data are required since resistivity monitoring has been introduced to quantitatively find out subsurface changes With respect to time. Here, we used the finite element method (FEM) for 3D resistivity modeling since the method is easy to realize complex topography and arbitrary shaped anomalous bodies. In the FEM, the linear elements, also referred to as first order elements, have certain advantages of simple formulation and narrow bandwidth of system equation. However, the linear elements show the poor accuracy and slow convergence of the solution with respect to the number of elements or nodes. To achieve the higher accuracy of finite element solution, high order elements are generally used. In this study, we developed a 3D resistivity modeling program using high order Serendipity elements. Comparing the Serendipity element solutions for a cube model with the linear element solutions, we assured that the Serendipity element solutions are more accurate than the linear element solutions in the 3D resistivity modeling.

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

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.148-154
• /
• 2004

This article reviews the development of three-dimensional (3-D) magnetotelluric (MT) modeling. The 3-D modeling of electromagnetic fields is essential in understanding the physics of MT soundings, and in implementing an inversion method to reconstruct a 3-D resistivity image. Although various numerical schemes have been developed over the last two decades, practical methods have been quite limited. However, the recent rapid improvement in computer speed and memory, as well as the advance in iterative solution algorithms for a large system of equations, makes it possible to model the MT responses of complex 3-D structures, which have been very difficult to simulate before. The use of staggered grids in finite difference method has become popular, conserving a magnetic flux and an electric current and allowing for realistic discontinuous fields. The convergence of numerical solutions has been greatly accelerated by adopting Krylov subspace methods, proper preconditioning techniques, and static divergence corrections. The vector finite-element method using edge elements is also free from the discontinuity problem, and seems a natural choice for modeling complex structures including irregular topography because its flexibility allows one to capture full geometric complexity.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.214-220
• /
• 2008

The distortion effect of electrical response for two-dimensional (2-D) DC resistivity method was verified in terms of 2-D inversion result of synthetic data obtained by three-dimensional (3-D) modeling, which is frequently applied to assess the safety of center core-type fill dam structure. The distortion effect is due to 2-D interpretation for 3-D structure. By the modeling analysis, we found that the water level is correctly described in the resistivity section around the middle part rather than each end side of the embankment due to the 3-D terrain effect, when the material of the embankment is assumed as horizontally uniform. And when we set the slope of outer rock fill part as uniform. the sharper the slope of the center core is, the more similar the resistivity section reflects. On the other hand, when the slope of the rock fill is steep, the resistivity section shows the water level at lower position than the real one, and the 3-D distortion effect at the end side of the embankment was enhanced.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.21 no.6
• /
• pp.135-146
• /
• 2013

This paper presents a feedforward control algorithm for the EGR and VGT systems of passenger car diesel engines. The air-to-fuel ratio and boost pressure are selected as control indicators and the positions of EGR valve and VGT vane are used as control inputs of the EGR and VGT controller. In order to compensate the non-linearity and coupled dynamics of the EGR and VGT systems, we have proposed a non-linear model-based feedforward control algorithm which is obtained from static model inversion approach. It is observed that the average modeling errors of the feedforward algorithm is about 2% using stationary engine experiment data of 225 operating conditions. Using a feedback controller including proportional-integral, the modeling error is compensated. Furthermore, it is validated that the proposed feedforward algorithm generates physically acceptable trajectories of the actuator and successfully tracks the desired values through engine experiments.

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

This paper presents sensitivity analysis of generalized rapid relaxation inversion (GRRI) algorithm for inverting controlled-source audio-frequency magnetotelluric (CSAMT) data. The algorithm was originally developed by modifying the RRI algorithm to recover a two-dimensional (2-D) conductivity structure of the Earth from MT data, but can be extended to include CSAMT data if it is combined with 2.5-D forward modeling. These GRRI approximate sensitivities are validated by comparison with exact 1-D and 2.5-D sensitivities. The comparison shows that the GRRI sensitivity is a good approximation to the exact sensitivity and has about half magnitude of the RRI sensitivity. Although the magnitude of the GRRI sensitivity is still slightly larger than that of the 2.5-D sensitivity, both sensitivities are broadly similar in shape when source-receiver offsets are greater than one skin depth on the Earth.