지구물리와물리탐사 (Geophysics and Geophysical Exploration)

  • 제5권4호
  • /
  • Pages.291-298
  • /
  • 2002
  • /
  • 1229-1064(pISSN)
  • /
  • 2384-051X(eISSN)
한국지구물리·물리탐사학회 (Korean Society of Earth and Exploration Geophysicists)
권호 표지

지형을 포함한 MT 탐사 자료의 2차원 모델링과 역산

Two-dimensional Modeling and Inversion of MT Data Including Topography

  • 이성곤 (한국지질자원연구원 탐사개발연구부) ;
  • 송윤호 (한국지질자원연구원 탐사개발연구부) ;
  • 김정호 (한국지질자원연구원 탐사개발연구부) ;
  • 정승환 (한국지질자원연구원 탐사개발연구부)
  • Lee Seong Kon (Geophysical Exploration and Mining Division, Korea Institute of Geoscience and Mineral Resource) ;
  • Song Yoonho (Geophysical Exploration and Mining Division, Korea Institute of Geoscience and Mineral Resource) ;
  • Kim Jung-Ho (Geophysical Exploration and Mining Division, Korea Institute of Geoscience and Mineral Resource) ;
  • Chung Seung-Hwan (Geophysical Exploration and Mining Division, Korea Institute of Geoscience and Mineral Resource)
  • 발행 : 2002.11.01
PDF 다운로드
⟨ 이전 논문 다음 논문 ⟩

초록

본 연구에서는 지형을 포함한 2차원 MT 역산 알고리듬을 개발하였다. 역산 과정시 필요한 모델 반응 계산을 위하여 유한요소법을 이용하였다. 공기와 지표면의 경계를 기준으로 고도에 따라 각 요소들의 절점을 수직으로 이동시킴으로써 추가적인 계산시간의 증가없이 간편히 지형을 구현하였다. 역산에서는 공간적인 함수로서 라그랑지 곱수를 결정하는 알고리듬을 채택하여 역산의 분해능과 안정성을 높이고자 하였다. 수치모델 실험을 통하여 TM과 TE모드 자료의 지형효과를 고찰하였고, 수치 자료의 역산을 통하여 지형을 포함한 역산의 타당성을 살펴보았다. 또한 현장 자료에 대하여 적용하여 본 연구에서 개발된 지형을 포함한 MT자료 역산 알고리듬의 적용성을 확인하였다.

We have developed a two-dimensional (2-D) magnetotelluric (MT) inversion algorithm, which can include topographic effects in inversion. We use the finite element method (FEM) to incorporate topography into forward calculation. Topography is implemented simply by moving nodes of rectangular elements in z-direction according to the elevation of air-earth interface. In the inversion process, we adopt a spatially variable Lagrangian multiplier algorithm in the smoothness-constrained least-squares inversion. The inversion algorithm developed in this study reconstructs subsurface resistivity structure quite well when topography variation exists. Also, it turns out to be effective in both resolution and stability from a model study and field data application.

키워드

참고문헌

  1. 김희준, 1990, 지자기 지전류 자료의 2차원 역산: 지질학회지, 26, 442-449
  2. 손정술, 1996, 유한요소법을 이용한 MT 탐사자료의 지형보정: 공학석사학위논문, 서울대학교
  3. 송윤호, 조인기 , 정승환, 김정호, 전정수, 1992, 제주도에서의 MT 탐사: 전기, 전자탐사 연구 KR-92-1G-5, 한국자원연구소, 33-69
  4. 송윤호, 조인기, 정승환, 김정호, 전정수, 1993, 2차원 MT 탐사-EMAP 탐사 : 전기, 전자탐사 연구 KR-93-1G-7, 한국자원연구소, 37-86
  5. 신동훈, 2000, 지형효과를 고려한 2차원 지자기지전류 역산: 공학석사학위논문, 서울대학교
  6. 이성곤, 현병구, 정승환, 송윤호, 1994, 2차원 MT 모델링과 공간 필터링을 이용한 정적효과 제거에 관한 연구: 한국자원공학회지, 31, 444-456
  7. 정승환,오재호 외 16인, 1994, MT 및 DGS 탐사에 의한 경상분지내 지열 및 탄화수소 탐사: 한국자원연구소 KR-94(T)-18, 과학기술처
  8. 조인기, 승윤호, 정승환, 김정호, 전정수, 1993, 원거리 기준점 MT 탐사: 전기, 전자탐사 연구 KR-93-1G-7, 한국자원연구소, 5-34
  9. Morrison, H.F., and Nichols, E. A., 1996, Continuous impedance profiling for mineral exploration: 66th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1286-1289
  10. Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P., 1992, Numerical recipes in FORTRAN: the art of scientific computing: Cambridge Univ. Press
  11. Rodi, W. L., 1976, A technique for improving the accuracy of finite element solutions for magnetotelluric data: Geophys. J. Roy. Astr. Soc., 44, 483-506
  12. Sasaki, Y., 1989, Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data: Geophysics, 54, 254-262
  13. Torres-Verdin, C., and Bostick, F. X. Jr., 1992, Principles of spatial surface electric field filtering in magnetotellurics: Electromagnetic array profiling (EMAP): Geophysics, 57, 603-622
  14. Yi, M.J., and Kim, J.H., 1998, Enhancing the resolving power of the least-squares inversion with active constraint balancing: 68th Ann. Internat. Mtg., Soc. ExpI. Geophys., Expanded Abstracts, 485-488
  15. Wannamaker, P. E., Stodt, J. A., and Rijo, L., 1986, Two-dimen-sional topographic responses in magnetotellurics modeled using finite elements: Geophysics, 51, 2131-2144