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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.10a
• /
• pp.71-74
• /
• 2009

호주 Victoria주에서 2007년에 이어 2008년에 Victoria주 북부 지역에서 추가로 탐사 자료를 획득하였으며, 이에 대한 2차원 및 3차원 해석을 수행하였다. 새로이 얻은 자료는 이전의 측선과 나란하게 설정하여 이전 결과에서 해석하였던 전기비저항 영상의 연장성을 밝히고자 하였다. 2차원 및 3차원 역산 결과를 이 지역의 지질자료와 비료 해석한 결과 2007년 MT자료 해석으로부터 확인된 고비저항 및 저비저항대의 공통적인 특징을 새로운 측선에서도 확인할 수 있었으며, 또한 이 지역의 대규모 단층대로 영상화할 수 있었다.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.318-325
• /
• 2008

Assigning an exact boundary condition is of great importance in three-dimensional (3D) magnetotelluric (MT) modeling, in which no source is considered in a computing domain. This paper presents a 3D MT modeling algorithm utilizing a Dirichlet condition for a 2D host. To compute boundary values for a model with a 2D host, we need to conduct additional 2D MT modeling. The 2D modeling consists of transverse magnetic and electric modes, which are determined from the relationship between the polarization of plane wave and the strike direction of the 2D structure. Since the 3D MT modeling algorithm solves Maxwell's equations for electric fields using the finite difference method with a staggered grid that defines electric fields along cell edges, electric fields are calculated at the same place in the 2D modeling. The algorithm developed in this study can produce reliable MT responses for a 3D model with a 2D host.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.286-301
• /
• 2008

When magnetotelluric (MT) data are obtained in vicinity of the coast, the sea can distort observed MT responses, especially those of deep part of subsurface. We introduce an iterative method to correct the sea effect, based on the previous topographic correction method which removes the distortions due to topographic changes in seafloor MT data. The method first corrects the sea effect in observed MT impedance, and then inverts corrected responses in a model space without the sea. Due to mutual coupling between sea and subsurface structure, the correction and inversion steps are iterated until changes in each result become negligible. The method is validated for 1-D and 2-D structure using synthetic MT data produced by 3-D forward modeling including surrounding seas. In all cases, the method closely recovers the given structure after a few iterations. To test the applicability of the proposed method to field data, we generate synthetic MT data for the Jeju Island whose 1-D conductivity structure is well known, using 3-D forward modeling. The distortions due to the surrounding sea start to appear below the frequency about 1 Hz, and are relatively severe in the electrical field perpendicular to the coastline because of the location of the observation sites. The proposed method successfully eliminates the sea effect after three iterations, and both 1-D and 2-D inversion of corrected responses closely recover the given subsurface structure of the Jeju Island model.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.71-78
• /
• 1998

Inversion schemes of 2-D MT survey data generally take enormous computational time and computer memory. In addition, careful attention must be paid in handling MT data, especially in cases of TM mode, inversion results can be seriously distorted because of static effect caused by current channeling across inhomogeneous surface boundaries. There-fore inversion algorithm using the GRRI scheme for TM mode MT data was implemented. This scheme is based on a perturbation analysis with a locally 2-D analysis and local inversions were sequently performed over each divided section without additional forward modelings. The algorithm was applied to several synthetic data for the purpose of verification of its efficiency and applicability. With less computer resources than conventional schemes, it could handle static effect directly by including current channeling across inhomogeneous boundaries. Thus it is expected to be used for an useful tool such as a real-time inversion scheme in the field.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.407-415
• /
• 2010

MT soundings were carried out in 2008, in northern Victoria, Australia, as a continuing collaboration research of 2007 between Republic of Korea, Australia, and Japan. The main purpose of this research is to investigate electrical conductivity structure and thus help understanding of tectonic structure in central Victoria, which is believed to be closely linked to mineralization and magmatic processes of this region. The survey area is located in western Lachlan Fold Belts, which is the part of Tasman Fold Belts in southeastern Australia. An MT profile of 2008 is almost parallel to the one of 2007 and approximately 50 km away. The 2D inversion result of MT data also shows that the position of conductivity discontinuity near surface are well matched with the positions of major faults, such as Avoca Fault, which is the structural boundary between Stawell and Bendigo Zones, and Heathcote Fault Zone, which marks the boundary between Bendigo and Melbourne Zones. It is also confirmed from resistivity image that internal faults in Bendigo Zone are in listric form, which is implied to be formed by structural shortening during compressional orogenic activity in Silurian.

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

Two-dimensional (2-D) inversion of magnetotelluric (MT) data for two survey lines having south-north direction from Jeiu Island has been carried out. The 2-D models show a thick layer having around 10 ohm-m in the depth of a few hundred meters throughout the survey area, which can be considered as the unconsolidated sedimentary layer. And they also show a conductive anomaly at the central part of each survey lines. But unfortunately by now, we do not have any further information about the anomaly. Comparison of the 2-D inversion model using MT band only and that using both AMT and MT bands said that it is helpful for us to include AMT band as well as MT band in the inversion to interpret not only the shallow part but also the deep structures.

• Geophysics and Geophysical Exploration
• /
• v.27 no.2
• /
• pp.129-143
• /
• 2024

In magnetotelluric (MT) surveys, small inhomogeneities near the surface cause a static shift in which apparent resistivities shift regardless of frequency. As the static shift in MT data leads to errors in subsurface structure interpretation, many studies have been conducted over the past few decades to mitigate or remove the distortions it caused. The most representative method involves removing static shifts from the data before inversion. Conversely, static shifts can be corrected during inversion or included in the inversion process. In addition, other geophysical data can be used to remove static shifts. However, the correction methods are limited to one-dimensional (1D) static responses, and limitations remain in two- or three-dimensional (2D or 3D) interpretation of distorted MT data owing to static shifts. This study provides a foundation for future studies on static shift by analyzing several previously published methods.