• 한국신재생에너지학회:학술대회논문집
• /
• 2010.11a
• /
• pp.121.2-121.2
• /
• 2010

제주도는 지질학적으로 제4기에 형성된 화산섬으로 지금까지 고온의 지열징후는 보고된 바 없으나, 남한에서 가장 최근까지 화산활동이 있었던 것으로 기록되어 있어 화산활동과 관련된 심부 지열자원 부존 가능성은 아직 열려있다고 할 수 있다. 본 연구에서는 제주도에서 지열부존 가능성을 타진하고 제주도 심부 지질구조 파악을 목적으로 2차원 및 3차원 자기지전류 (MT) 탐사를 수행하였다. 탐사는 중산간지역에서 한라산을 중심으로 동, 서, 남, 북의 4방향 4측선과 제주 서부지역에 남북방향의 1측선을의 총 5개 측선에 대해 수행하였으며, 이에 대한 MT 탐사자료의 2차원 및 3차원 역산을 통하여 한라산 하부 및 주변의 심부 지질구조를 파악하고자 하였다. 역산 해석 결과는 천부 구조는 기존 시추조사 결과 밝혀진 층서구조의 형태를 잘 나타내어 획득된 자료의 신뢰도가 높음을 지시하였다. 즉, 제주도 최 상부를 피복하고 있는 현무암 등의 화산암류는 고비저항(수백 ohm-m)으로, 그 하부의 해성 미고결퇴적층(U층 및 서귀포층)은 저비저항으로, 그리고 최하부의 응회암이나 화강암으로 구성된 기반암은 1,000 ohm-m 이상의 고비저항 층으로 잘 구분되어 나타났다. 특히, 제주도에서 특징적으로 해수면 하부 수십 ~ 수백 m에 존재하는 것으로 알려진 미고결퇴적층이 10 ohm-m 내외로 측선 전반에 걸쳐 나타났다. 이는 기존의 시추결과에서 미고결 퇴적층이 제주도 전역에 걸쳐 해수면 하부 100 m 내외의 심도에서 관찰되는 것과 일치하는 결과이다. 기반암 하부에서는 특징적으로 모든 측선의 중앙부에서 저비저항 이상대가 영상화되었으며 이는 2차원 역산과 3차원 역산해석에서 공통적으로 나타났다. 특히, 3차원 해석에서는 이러한 저비저항 이상대가 한라산 정상에서 서북쪽 부근에 나타나는데 이는 과거의 화산활동과 관련된 지질학적인 구조에 의한 영향일 가능성과 측선의 양단과 중앙에서 주변 바다의 영향이 다르게 나타나기 때문일 가능성으로 볼 수 있다. 즉, 전자는 심부에 발단된 각각의 파쇄대가 모든 측선의 중앙부에서 교차하거나 이를 통한 한라산의 생성과정과 연관된 지질학적인 구조일 가능성을 의미한다. 만약 한라산을 형성한 화성활동의 영향이 아직 지하 심부에 남아있다면 지열수의 부존 혹은 마그마의 통로가 되었을 단층의 영향으로 한라산 하부에 저비저항 이상대로 나타날 가능성이 높다. 그러나 후자에 의한 가능성도 배제할 수는 없으므로 향후 주변바다에 대한 영향을 고려한 3차원 역산해석이나 심부시추 등을 통한 상세한 지질조사 등 추가적인 연구가 이루어져야 할 것으로 판단된다.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.184-196
• /
• 2008

We perform MT soundings in Bendigo, the northern part of Victoria, Australia, to investigate the deep subsurface geologic structure. The primary purpose of this survey is to figure out whether the discontinuity such as faults extends northward. The time series of MT signal were measured over 11 days at 71 measurement stations together with at remote reference, which help enhance the quality of impedance estimation and its interpretation. The impedances are estimated by robust processing using remote reference technique and then inverted with 2D MT 2D inversion. We can see that known faults are clearly imaged in MT 2D inversion. Comparing resistivity images from MT 2D inversion with interpreted boundary from reflection seismic exploration, two interpretations match well each other.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.89-97
• /
• 2007

A three-dimensional (3D) magnetotelluric (MT) survey has been carried out to delineate subsurface structures and possible fractures, for development of low-temperature geothermal resources in Pohang, Korea. Quite good quality MT data could be obtained throughout the survey region by locating the remote reference in Kyushu, Japan, which is ${\sim}480\;km$ from the centre of the field site. 3D modelling and inversion are performed taking into account the sea effect in MT measurements near the seashore. The nearby sea in the Pohang area affects MT data at frequencies below $1\;Hz{\sim}0.2\;Hz$, depending on the distance from the seashore. The most severe sea effects were observed in the south-east parts of the survey area, closer to Youngil Bay. 3D inversion with and without the seawater constraint showed very similar results at shallow depths, roughly down to 2 km. At greater depths, however, a strong sea effect seems to form a fictitious conductive structure in ordinary 3D inversion, especially in the south-eastern part of the survey region. Comparison between drilling results and the resistivity profiles from inversions showed that five layered structures can be distinguished the subsurface beneath the target area. They are: (a) semi-consolidated mudstones with resistivity less than $10\;{\Omega}m$, which are ${\sim}300\;m$ thick in the northern part and ${\sim}600\;m$ thick in the southern part of the survey area; (b) occasional occurrence of trachybasalt and lapilli tuff within the mudstone layer has resistivity of a few tens of${\Omega}m$, (c) intrusive rhyolite ${\sim}400\;m$ thick has resistivity of several hundreds of ${\Omega}m$, (d) alternating sandstone and mudstone down to 1.5 km depth shows resistivity of ${\sim}100\;{\Omega}m$, (e) a conductive structure was found at a depth of ${\sim}3\;km$, but more geological and geophysical study should be carried out to identify this structure.

• 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.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.10 no.4
• /
• pp.345-352
• /
• 2007

Two-dimensional (2D) MT surveys near the Century mine in Australia have been performed with very far remote reference in Esashi, Japan (RR_ESS) as well as Gregory Downs (RR_GREG), which are roughly 6,400 km and 80 km apart from the survey area, respectively. Good quality of MT data could be obtained by remote reference processing with RR_GREG, while the coherency of magnetic fields between field sites and RR_ESS was not sufficient to be used as remote data. Both 2D and 3D inversion of 2D profile data represented the general geological structure beneath the survey area. The main target of the survey, Termite Range Fault, appeared as a boundary between a conductive block to the north and a resistive block to the south in the reconstructed resistivity section, and is inclined slightly to the north-east direction.

• Journal of the Korean earth science society
• /
• v.32 no.6
• /
• pp.548-559
• /
• 2011

3-D Multi-geophysical surveys were carried out around the Hwasan caldera at the Euisung Sub-basin. To overcome the limitations of resolutions in previous studies, dense gravity data and magnetotelluric (MT) data were obtained and analyzed. In this study, the independent inversion models from gravity and MT data were integrated using correlation and classification approaches for 3-D imaging of the geologic structures. A Structure Index (SI) method was proposed and applied to the integration and classification analyses. This method consists of Type Angle (TA) and Type Intensity (TI) values, which are estimated by the spatial correlation and abnormality of the physical properties. The SI method allowed the classification analysis to be effectively performed. Major findings are as follows: 1) pyroclastic rocks around the central area of the Hwasan caldera with lower density and resistivity than those of neighboring regions extended to a depth of around 1 km, 2) intrusive igneous rocks with high resistivity and density were imaged around the ring fault boundary, and 3) a basement structure with low resistivity and high density, at a depth of 3-5 km, was inferred by the SI analysis.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.299-308
• /
• 2009

A magentotelluric (MT) survey at the Uiseong area has been performed for the site investigation of pilot scale $CO_2$ sequestration. The purpose of the MT survey is to delineate deeply extended fracture systems that can act as a leakage path of injected $CO_2$ Plume. Since the target area is extremely noisy in electromagentic sense, low frequency data below 1 Hz cannot be used for inversion. Two- and three-dimensional interpretation of the MT data showed a very clear conductive anomaly, which has the direction of $N55\sim65^{\circ}W$ and is extended roughly down to 1.6 km. It have the same direction with the strike-slip faults, the Gaeum and Geumcheon Faults. On the contrary, the eastern part of the survey area shows relatively homogeneous to the depth of 2 km though some small fractures at shallow depths can be found. Test drilling and high-definition borehole surveys should be followed at the eastern part of the survey area and hydraulic fracturing is required for injection of $CO_2$, because mean porosity of the sedimetary rock in the area is only 1.47%.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.145-155
• /
• 2005

A two-dimensional (2-D) interpretation of MT data has been performed for the purpose of fracture detection for geothermal development. Remote stations have been operated in Kyushu, Japan (480 km apart) as well as in Korea (60 km and 165 km apart in 2002 and 2003 data set, respectively). Apparent resistivity and phase curves calculated by remote processing with the Japan remote data showed enough quality for 2-D inversion for the whole frequency range. Remote reference processing with Korea remote reference data also showed quite good continuity in apparent resistivity and phase curves except some noisy frequency bands; around the power frequency, 60 Hz, and around the dead band $10^{-1}Hz\;Hz\;\~1\;Hz$, where the natural EM signal is known to be very weak. Even though the subsurface showed severe three-dimensional (3-D) characteristics in the survey area so that 2-D inversion by itself could not give enough information for deep geological structures, the 2-D inversion for the 5 survey lines showed several common features. The conductive semi-consolidate mudstone layer is dipping from north to south (about 500 m depth on the south and 200 m on the north most part of the survey area). The boundary between the low (L-2) and high (H-2) resistivity anomalies can be thought as a major fault with strike $N15^{\circ}E$, passing through the sites 206, 112 and 414. The shallow (< 1 km) conductive anomalies (L-4) seem to be fracture zones having strike E-W (at site 105) and $N60^{\circ}W$ (at site 434). And there exists a conductive layer in the western and west-southern part of the survey area in the depth below $2\~3\;km$, for which further investigation is to be needed.

• Geophysics and Geophysical Exploration
• /
• v.17 no.1
• /
• pp.28-33
• /
• 2014

We investigate the spatial distribution of highly conductive layer using the one-dimensional inversions of the new magnetotelluric (MT) measurements obtained at the mid-mountain (400 ~ 900 m in elevation) western area of Jeju Island and the previous MT data over Jeju Island, Korea. The conductive layer indicates the sedimentary layer comprised of Seoguipo Fomation and U Formation. There is a definite positive correlation between the top of conductive layer and the earth surface in elevation. On the contrary, the bottom of conductive layer has a negative correlation with the surface elevation. In other words, the conductive layer has a shape of convex lens, which is thickest in the central part. The basement beneath the conductive layer could be concave in the central part of Jeju Island. A kriging considering the correlation between the layer boundary and the surface elevation provides a reliable geoelectric structure model of Jeju Island. However, further studies, i.e. three-dimensional modeling and interpretation integrated with other geophysical or logging data, are required to reveal the possible presence of three-dimensional conductive body near the subsurface vent of Mt. Halla and the causes of the bias in the depths of layer estimated from MT and core log data.