FDM을 이용한 MT 탐사의 3차원 모형 반응 연구

An Interpretive Analysis of Magnetotelluric Response for a Three-dimensional Body Using FDM

  • 한누리 (서울대학교 공과대학 지구환경시스템공학부 물리탐사연구실) ;
  • 이성곤 (한국지질자원연구운) ;
  • 송윤호 (한국지질자원연구운) ;
  • 서정희 (서울대학교 공과대학 지구환경시스템공학부 물리탐사연구실)
  • 발행 : 2004.05.01

초록

본 논문에서는 자기지전류 탐사법에 대한 3차원 이상체의 모형 반응 결과를 분석하였다. 분석에 이용된 3차원 모형은 비전도성 배경 매질 내에 전도성 고립이상체가 존재하는 모형과, 같은 모형에 전도성 표토층을 추가한 모형이다. 지하의 전도성 이상체에 의한 겉보기 전기비저항 이상에는 뚜렷한 주파수 의존성이 존재하며, 이는 이상체 주변에 발생되는 전류 집중 및 유도 전류 거동의 주파수 의존성 때문이다. 지표에서 tipper와 induction vector의 반응에도 뚜렷한 주파수 의존성이 존재하였으며, 그 외에도 tipper와 induction vector의 반응은 이상체의 위치와 직접적인 연관성이 존재하였다. 또한, 2차원 탐사를 가정하여 한 측선에서 겉보기 전기비저항과 위상, induction vector를 분석하면, 3차원 이상체의 존재 여부 및 위치 파악에 도움을 줄 것으로 판단되었다. 전도성 표토층을 추가한 모형의 반응은 대체적으로 고립이상체 모형의 반응과 비슷한 양상을 나타내었지만, 전도성 표토층의 영향으로 이상체에 의한 반응의 크기가 감소하였으며, 이상체의 영향이 나타나는 영역도 크게 감소하였다. 위와 같이 본 연구에서 논의된 3차원 모형 반응에 대한 분석 결과 및 그 방법은 지하 구조의 3차원 반응의 이해와 MT 탐사 자료의 해석에 효과적인 지침이 될 것으로 생각된다.

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.

키워드

참고문헌

  1. 송윤호, 박영수, 임무택, 조성준, 이성곤, 손정술, 김형찬, 기원서, 고희재 , 임성근, 정지민, 2003, 지열에너지 개발을 위한 심부물리탐사, KR-2002-B-14, 한국지질자원연구원
  2. 양준모, 오석훈, 이덕기, 권병두, 윤용훈, 2002, 3차원 전도체의 공간적 위치 및 크기에 따른 차이지시자의 특성 연구, 지구물리, 5, 305-319
  3. 이성곤, 1993, 2차원 MT 모델링과 공간 필터링을 이용한 정적효과 제거에 관한 연구, 공학석사 학위논문, 서울대학교
  4. Berdichevsky, M. N., and Dmitriev, V. I., 2002, Magnetotellurics in the context of the theory of ill-posed problems, Soc. Explor. Geovhys
  5. Dosso, H. W., and Chen, J., 2000, Analogue model study of EM induction in elongated conductors-2D and 3D induction arrow responses, Earth Planets Space, 52, 355-360
  6. Mackie, R. L., Madden, T. R., and Wannamaker, P. E., 1993, Three-dimensional magnetotelluric modeling using difference equations-Theory and comparisons to integral equation solutions, Geophysics, 58, 215-226
  7. Mackie, R. L., Smith, J. T., and Madden, T. R., 1994, Three-dimensional electromagnetic modeling using finite difference equations - The magnetotelluric example, Radio Science, 29, 923-935
  8. Park, S. K., 1985, Distortion of magnetotelluric sounding curves by three-dimensional structures, Geophysics, 50, 785-797
  9. Park, S. K., Orange, A. S., and Madden, T. R., 1983. Effects of three-dimensional structure on magnetotelluric sounding curves, Geophysics, 48, 1402-1405
  10. Sasaki, Y., 1989, Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data, Geophysics, 54, 254-262
  11. Vozoff, K., 1991, The magnetotelluric method, in Misac N. Nabighian, Ed., Electromagnetic methods in applied geophysics Vol. II, Soc. Explor. Geophys
  12. Yee, K. S., 1966, Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media, IEEE Trans. Ant. Prop., AP-14, 302-307