• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.136-147
• /
• 2004

In this study, the characteristics of magnetotelluric (MT) responses due to a three-dimensional (3-D) body are analyzed with 3-D numerical modeling. The first model for the analysis consists of a single isolated conductive body embedded in a resistive homogeneous half-space. The second model has an additional conductive overburden while the other conditions remain the same as the first one. The analysis of apparent resistivities shows well that the 3-D effects are dominant over some frequency range for the first model. Two mechanisms, current channeling and induction, for secondary electric fields due to the conductive body are analyzed at various frequencies: at high frequencies induction is more dominant than channeling, while at low frequencies channeling is more dominant than induction. Tippers have a strong relation to the position of anomalous body and the real and imaginary parts of induction vector also indicate the position of anomalous body. off-line conductive anomaly sometimes causes severe problem in 2-D interpretation. In such case, induction vector analysis can give information on the existence and location of the anomalous body. Each parameter of the second model shows similar responses as those of the first model. The only difference is that the magnitude of all parameters is decreased and that the domain showing the 3-D effects becomes narrower. As shown in this study, the analysis of 3-D effects provides a useful and effective means to understand the 3-D subsurface structure and to interpret MT survey data.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.37-50
• /
• 2024

The magnetotelluric (MT) survey can be affected by external environmental factors. In particular, when acquiring MT data in islands, it is essential to consider the combined effect of topography and sea to understand the results and make accurate interpretations. To analyze the MT response (apparent resistivity, phase) with consideration of the effect of topography and sea, a small cone-shaped island model surrounded by deep sea was created. Two-dimensional (2-D) and three-dimensional (3-D) forward modeling were performed on the terrain model considering topography and the island model considering both topography and sea. The 2-D MT response did not reflect the topographic and sea effect of the direction orthogonal to the 2-D profile. The 3-D MT response included topographic and sea effects in all directions. The XY and YX components of the apparent resistivity were separated on undulating topography, such as a hill. A conductor at 1 km below sea level could be distinguished from topographic and sea effects in the MT response, and low resistivity anomaly was attenuated at greater depths. This study will facilitate understanding of field data measured on small islands.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2004.06a
• /
• pp.182-188
• /
• 2004

In view of the circular geometry of the Serra da Cangalha impact crater, we have carried out a 3D forward modeling computation for twenty-five MT data in order to obtain the 3D resistivity forward model for the crater region. The 3D resistivity forward model revealed a five-layer model, showing a significant reduction in the basement resistivity. We suggest that this, perhaps, could be due to the structural disturbances that have been caused by the meteorite impact on the crater about 220 million years ago resulting in brecciation, fracturing, alteration and shocked zone filled with fluids. Also, the sensitivity analysis of the 3D model chosen indicates that 3D models having a crater diameter greater than 151 are inconsistent with our data because the 3D model responses are very sensitive to changes in the diameter beyond 15 km. This analysis also reveals that, the depth limits (for the 3D body) causing the anisotropic effects seen on some of our apparent resistivity curves maximally does not extend beyond 1.2 km depth.

• Economic and Environmental Geology
• /
• v.27 no.5
• /
• pp.491-498
• /
• 1994

The integral equation method is used for magnetotelluric (MT) modeling of a finite inhomogeneity in a two-layered earth. An integral equation relates the incident plane-wave field and the scattering currents in the three-dimensional (3-D) inhomogeneity through the electric tensor Green's function appropriate to a layered earth. This paper describes an effect of overburden and basement on MT responses of 3-D body. The effect of overburden is to reduce the detectability of target, and the reduction of detectability is more apparent for conductive overburden than for resistive one. The effect of basement, on the other hand, may enhance the anomaly due to 3-D body in the upper layer. In case of the resistive basement current perturbations about the body tend to be confined to the more conductive upper layer.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.3
• /
• pp.965-974
• /
• 2013

Multiple debris flows occurred on July 27, 2012 in Mt. Umyeon, which resulted in 16 casualties and severe property demage. Accurate reproducing of the propagation and deposition of debris flow is essential for mitigating these disasters. Through applying FLO-2D model to these debris flows and comparing the results with field observations, we seek to evaluate the performance of the model and to analyse the rheological model parameters. Representative yield stress and dynamic viscosity back-calculated for the debris flows in the northern side of Mt. Umyeon are 1022 Pa and 652 $Pa{\cdot}s$, respectively. Numerical results obtained using these parameters reveal that deposition areas of debris flows in Raemian and Shindong-A regions are well reproduced in 63-85% agreement with the field observations. However, the propagation velocities of the flows are significantly underestimated, which is attributable to the inherent limitations of the model that can't take the entrainment of bed material and surface water into account. The debris flow deposition computed in Hyeongchon region where the entrainment is not significant appears to be in very good agreement with the field observation. The sensitivity study of the numerical results on model parameters shows that both sediment volume concentration and roughness coefficient significantly affect the flow thickness and velocity, which underscores the importance of careful selection of these model parameters in FLO-2D modeling.

• 한국지구과학회:학술대회논문집
• /
• 2005.09a
• /
• pp.124-131
• /
• 2005

A GDS (Geomagnetic Depth Sounding) method, one of extremely low-frequency EM methods, has been carried out to examine deep geo-electrical structures of the Korean peninsula. In this study, five additive GDS sites acquired in south-eastern area of the Korea were integrated into twelve previous GDS results. In addition, 3-D MT modeling considering the surrounding seas of the Korean peninsula was performed to evaluate sea effect at each GDS site quantitatively. As a result, Observed real induction arrows was not explained by solely sea effect, two conductive structures that are able to explain differences between observed and calculated induction arrows, was suggested. The first conductive structure is the Imjingang Belt, which is thought as a extension of Quiling-Dabie-sulu continental collision belt. The effects of the Imjingang Belt clearly appear at YIN and ICHN sites. The second one is the HCL (Highly Conductive Layer), which is considered as a conductive anomaly by mantle upwelling generated in back-basin region. The effects of the HCL are also confirmed at KZU, KMT101, 107 sites, in the south-eastern of the Korean peninsula.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• 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
• /
• v.12 no.4
• /
• pp.361-366
• /
• 2009

Layered-earth Green's functions in electormagnetic (EM) surveys play a key role in modeling the response of exploration targets. They are computed through the Hankel transforms of analytic kernels. Computational precision depends upon the choice of algebraically equivalent forms by which these kemels are expressed. Since three-dimensional (3D) modeling can require a huge number of Green's function evaluations, total computational time can be influenced by computational time for the Hankel transform evaluations. Linear digital filters have proven to be a fast and accurate method of computing these Hankel transforms. In EM modeling for 3D inversion, electric fields are generally evaluated by the secondary field formulation to avoid the singularity problem. In this study, three components of electric fields for five different sources on the surface of homogeneous half-space were derived as primary field solutions. Moreover, reflection coefficients in TE and TM modes were produced to calculate EM responses accurately for a two-layered model having a sea layer. Accurate primary fields should substantially improve accuracy and decrease computation times for Green's function-based problems like MT problems and marine EM surveys.

• Journal of the Korean Geophysical Society
• /
• v.9 no.3
• /
• pp.159-169
• /
• 2006

A GDS (Geomagnetic Depth Sounding) method, one of extremely low-frequency EM methods, has been carried out to examine conductivity anomalies in and around the Korean Peninsula. In this study, new GDS data acquired at the five sites in south-eastern area of the peninsula were incorporated into the previous GDS data. In order to quantitatively interpret observed induction arrows, the 3-D MT modeling considering the surrounding seas of the Korean Peninsula has been performed to evaluate sea effect at each GDS site. The modeling results revealed that the observed real induction arrows were not explained by solely sea effects, consequently two conductive structures that are responsible for the discrepancies between observed and calculated induction arrows were proposed. The first one is the Imjingang Belt, which is thought as an extension of Quiling-Dabie-sulu continental collision belt. The effects of the Imjingang Belt clearly appear at the site YIN and ICHN. The second one is the HCL (Highly Conductive Layer), which is considered as a conductive anomaly by mantle upwelling produced in back-basin region. The effects of the HCL are seen at the site KZU, KMT101, and KMT 107 in the south-eastern region of the Korean Peninsula.

• 한국지구과학회:학술대회논문집
• /
• 2003.09a
• /
• pp.35-36
• /
• 2003