• Geosystem Engineering
• /
• v.21 no.5
• /
• pp.251-261
• /
• 2018

With the urbanization in recent years, the power line interference noise in electromagnetic signal is increasing day by day, and has gradually become an unavoidable component of noises in magnetotelluric signal detection. Therefore, a kind of power line interference noise elimination method based on independent component analysis in wavelet domain for magnetotelluric signal is put forward in this paper. The method first uses wavelet decomposition to change single-channel signal into multi-channel signal, and then takes advantage of blind source separation principle of independent component analysis to eliminate power line interference noise. There is no need to choose the layer number of wavelet decomposition and the wavelet base of wavelet decomposition according to the observed signal. On the treatment effect, it is better than the previous power line interference removal method based on independent component analysis. Through the de-noising processing to actual magnetotelluric measuring data, it is shown that this method makes both the apparent resistivity curve near 50 Hz and the phase curve near 50 Hz become smoother and steadier than before processing, i.e., it effectively eliminates the power line interference noise.

• Journal of the Korean earth science society
• /
• v.30 no.7
• /
• pp.834-844
• /
• 2009

Two-dimensional magnetotelluric (MT) inversion algorithm using minimum support (MS) stabilizer functional was implemented in this study to enhance the contrast of inverted images. For this implementation, this study derived a formula in discrete form for creeping model updates in the least-squares linearized inversion. A spatially varying regularization parameter determination algorithm, which is known as ACB (Active Constraint Balancing), was also adopted to stabilize the inversion process when using MS stabilizer as a model constraint. Inversion experiments for a simple isolated body model show well the feature of MS stabilizer in concentrating the anomalous body compared with the second-order derivative model constraint. This study also compared MS stabilizer and the second-order derivative model constraints for a model having multiple anomalous bodies to show the applicability of the algorithm into field data.

• Journal of the Korean Society for Geothermal and Hydrothermal Energy
• /
• v.4 no.1
• /
• pp.1-10
• /
• 2008

A magnetotelluric survey and geological survey have been performed to delineate deeply extended fracture systems at the geothermal field in Seokmo Island, Inchon, Korea. One borehole(BH-1) succeed to meet a large fracture system at the depth of 750 m where approximately $72^{\circ}C$ geothermal water is overflowing, while the other borehole(BH-2), which is about 200 m eastward from BH-1, failed to develop enough geothermal water even at the depth of 1,200m. Though there have been so many electric noise sources around the survey region, good quality of MT data above 1 Hz could be obtained with careful installation, remote reference processing. Inversion of MT data for two lines roughly perpendicular to the two major lineaments respectively show that the two lineaments are related to the fracture systems that are extended at least down to 1.5 km depth and inclined eastwards. From the interpretation, additional drilling for BH-2 is recommended and finally meet the fracture systems at the depth of 1,280 m and resulted in overflow of large amount of geothermal water of temperature $69.4^{\circ}C$ from BH-2.

• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.261-270
• /
• 2006

Two-dimensional (2-D) inversion of magnetotelluric (MT) data for two survey lines having south-north direction from Jeju Island has been carried out. Broad band MT sounding curves with good quality could be gathered by performing audio-frequency magnetotelluric (AMT) survey during the MT survey and by operating the remote reference in Kyushu Island, Japan. Comparison of the 2-D inversion model using MT band only and that using both AMT and MT bands for the field data as well as for the data from numerical 2-D modeling said that high frequency information from AMT survey can be useful for interpreting not only the shallow part but also the deep structures, especially when the formation is resistive. The 2-D inversion models of field data 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. It can be either the effect of the surrounding sea water, or the structures due to ancient volcanic events. But unfortunately by now, we do not have any further information about the anomaly.

• Journal of the Korean earth science society
• /
• v.30 no.5
• /
• pp.550-564
• /
• 2009

When magnetotelluric (MT) data are obtained in the vicinity of the coast, the surrounding seas make it difficult to interpret subsurface structure, especially the deep part of the subsurface. We introduce an iterative method to correct the sea effect, based on the previous topographic correction method that removes the distortion due to topographic changes in seafloor MT data. The method first corrects the sea effect in observed MT impedance, and then inverts corrected response in a model space without the sea. Due to mutual coupling between the sea and the subsurface structure, the correction and inversion steps are iterated until the changes in each result become negligible. The method is tested for 1- and 2-D structures using synthetic MT data produced by 3-D forward modeling including surrounding seas. In all cases, the method closely recovers the true structure assumed to generate synthetic responses after a few iterations.

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

An efficient three-dimensional (3D) inversion of magnetotelluric (MT) data can be carried out by using approximate sensitivities or avoiding the calculation of a full sensitivity matrix. In this paper, we propose approximate sensitivities for efficient 3D MT inversion based on the Gauss-Newton method and test and compare four kinds of sensitivities. Applying the four sensitivities to both synthetic and field data shows that the effects of sensitivities are highly dependent on data and thus applying various combinations of sensitivities is recommended for efficient inversion and good images.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.10 no.2
• /
• pp.147-153
• /
• 2007

The conjugate gradient (CG) method is one of the most efficient algorithms for solving a linear system of equations. In addition to being used as a linear equation solver, it can be applied to a least-squares problem. When the CG method is applied to large-scale three-dimensional inversion of magnetotelluric data, two approaches have been pursued; one is the linear CG inversion in which each step of the Gauss-Newton iteration is incompletely solved using a truncated CG technique, and the other is referred to as the nonlinear CG inversion in which CG is directly applied to the minimization of objective functional for a nonlinear inverse problem. In each procedure we only need to compute the effect of the sensitivity matrix or its transpose multiplying an arbitrary vector, significantly reducing the computational requirements needed to do large-scale inversion.

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.232-240
• /
• 2015

Magnetotelluric (MT) survey investigates electrical structure of subsurface by measuring natural electromagnetic fields on the earth surface. For the accurate interpretation of MT data, the precise three-dimensional (3-D) modeling algorithm is prerequisite. Since MT responses are affected by electrical anisotropy of medium, the modeling algorithm has to incorporate the electrical anisotropy especially when analyzing time-lapse MT data sets, for monitoring engineered geothermal system (EGS) reservoir, because changes in different-vintage MT-data sets are small. This study developed a MT modeling algorithm for the simulation MT responses in the presence of electrical anisotropy by improving a pre-existing staggered-grid finite-difference MT modeling algorithm. After verifying the developed algorithm, we analyzed the effect of vertical transversely isotropic (VTI) anisotropy on MT responses. In addition, we are planning to extend the applicability of the developed algorithm which can simulate not only the horizontal transversely isotropic (HTI) anisotropy, but also the tiled transversely isotropic (TTI) anisotropy.

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

A multi-geophysical surveys were carried out at Hwasan caldera which is located in Euisung Sub-basin. In order to overcome the limitation of the previous studies, dense gravity data and magnetotelluric (MT) data were obtained and integrated. In this study, the independent inversion models from gravity and MT method were integrated using a correlation and classification approaches to map geologic structure. The results of integration analysis indicated followings; 1) pyroclastic rocks around the central area of Hwasan caldera have lower density and resistivity when compared with those of neighborhood regions and are extended to around 1 km in depth, 2) the high resistivity and density intrusive igneous rocks are imaged around the ring fault boundary, and 3) the basement structure, which has low resistivity and high density, 5 km deep inferred by integration analysis. Also, for integration analysis, we suggested Structure Index method. This method is analyzed using Type Angle and Type Intensity, which are calculated by the spatial correlation of the physical properties. In this study, we can perform the integration analysis effectively using Structure Index method.