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

Korean Society of Earth and Exploration Geophysicists (한국지구물리·물리탐사학회)
권호 표지
DOI QR코드

DOI QR Code

On the Efficient Three-Dimensional Inversion of Static Shifted MT Data

정적효과를 포함한 자기지전류 자료의 효율적인 3차원 역산에 관하여

  • Jang, Hannuree (Energy Resources Institute, Pukyong National University) ;
  • Jang, Hangilro (Department of Energy Resources Engineering, Pukyong National University) ;
  • Kim, Hee Joon (Department of Energy Resources Engineering, Pukyong National University)
  • 장한누리 (부경대학교 에너지자원연구소) ;
  • 장한길로 (부경대학교 에너지자원공학과) ;
  • 김희준 (부경대학교 에너지자원공학과)
  • Received : 2013.05.10
  • Accepted : 2014.03.18
  • Published : 2014.05.31
Download PDF
⟨ Previous Next ⟩

Abstract

This paper presents a practical inversion method for recovering a three-dimensional (3D) resistivity model and static shifts simultaneously. Although this method is based on a Gauss-Newton approach that requires a sensitivity matrix, the computer time can be greatly reduced by implementing a simple and effective procedure for updating the sensitivity matrix using the Broyden's algorithm. In this research, we examine the approximate inversion procedure and the weighting factor ${\beta}$ for static shifts through inversion experiments using synthetic MT data. In methods using the full sensitivity matrix constructed only once in the iteration process, a procedure using the full sensitivity in the earlier stage is useful to produce the smallest rms data misfit. The choice of ${\beta}$ is not critical below some threshold value. Synthetic examples demonstrate that the method proposed in this paper is effective in reconstructing a 3D resistivity structure from static-shifted MT data.

MT 자료에서 3차원 전기비저항 구조와 정적효과를 동시에 구하기 위한 실용적인 역산법을 소개한다. 이 방법은 감도행렬이 필요한 Gauss-Newton법을 기반으로 하고 반복과정에서 Broyden의 방식으로 감도를 수정하는 것을 기본으로 하고 있다. 이 논문에서는 합성 MT 자료에 대한 역산실험을 통해 근사역산법의 성능과 정적효과에 대한 가중치에 대해 검토하였다. 해석적으로 구해지는 초기감도를 Broyden의 방식으로 수정하는 역산법은 초기감도를 끝까지 쓰는 역산법보다 자료오차를 줄이는데 효과적이었다. 그리고 완전한 감도행렬을 반복 중간에서 단 한번만 사용하는 근사역산법으로서는 반복 전반부에서 완전한 감도를 사용할 때 자료오차를 가장 많이 줄이는 것으로 나타났다. 정적효과에 대한 가중치는 어느 특정 한계값 이하로 선택하면 최종 자료오차에는 결정적인 영향을 주지 않는다. 합성 MT 자료에 대한 실험 결과 이 역산법은 정적효과가 포함된 MT 자료로부터 3차원 전기비저항 구조를 재현하는데 효과적임을 확인하였다.

Keywords

References

  1. deGroot-Hedlin, C., 1991, Removal of static shift in two dimensions by regularized inversion, Geophysics, 56, 2102-2106. https://doi.org/10.1190/1.1443022
  2. Groom, R. W., and Bailey, R. C., 1989, Decomposition of magnetotelluric impedance tensors in the presence of local three-dimensional galvanic distortion, J. Geophys. Res., 94, 1913-1925. https://doi.org/10.1029/JB094iB02p01913
  3. Han, N., Nam, M. J., Kim, H. J., Lee, T. J., Song, Y., and Suh, J. H., 2007, An Efficient 3D Inversion of MT Data Using Approximate Sensitivities, Mulli-Tamsa, 10(4), 259-267.
  4. Han, N., Nam, M. J., Kim, H. J., Lee, T. J., Song, Y., and Suh, J. H., 2008, Efficient three-dimensional inversion of magnetotelluric data using approximate sensitivities, Geophys. J. Int., 175, 477-485. https://doi.org/10.1111/j.1365-246X.2008.03894.x
  5. Jones, A. G., 1988, Static shift of magnetotelluric data and its removal in a sedimentary basin environment, Geophysics, 53, 967-978. https://doi.org/10.1190/1.1442533
  6. Lee, S. K., Kim, H. J., Song, Y., and Lee, C.-K., 2009, MT2DInvMatlab-A program in MATLAB and FORTRAN for two-dimensional magnetotelluric inversion, Computers & Geosciences, 35, 1722-1734. https://doi.org/10.1016/j.cageo.2008.10.010
  7. Lee, T. J., Uchida, T., Sasaki, Y., and Song, Y., 2003, Characteristics of static shift in 3-D MT inversion, Mulli-Tamsa, 6, 199-206.
  8. Lee, T. J., Song, Y., and Uchida, T., 2007, Three-dimensional magnetotelluric surveys for geothermal development in Pohang, Korea, Exploration Geophysics, 38, 44-49. https://doi.org/10.1071/EG07010
  9. Loke, M. H., and Barker, R. D., 1996, Practical technique for 3D resistivity surveys and data inversion, Geophys. Prosp., 44, 499-523. https://doi.org/10.1111/j.1365-2478.1996.tb00162.x
  10. McGillivray, P. R., Oldenburg, D. W., Ellis, R. G., and Habashy, T. M., 1994, Calculation of sensitivities for the frequency-domain electromagnetic problem, Geophys. J. Int., 116, 1-4. https://doi.org/10.1111/j.1365-246X.1994.tb02121.x
  11. Ogawa, Y., and Uchida, T., 1996, A two-dimensional magnetotelluric inversion assuming Gaussian static shift, Geophys. J. Int., 126, 69-76. https://doi.org/10.1111/j.1365-246X.1996.tb05267.x
  12. Sasaki, Y., 1999, Three-dimensional frequency-domain electromagnetic modeling using the finite-difference method, Butsuri-Tansa, 52, 421-431. (Japanese with English abstract)
  13. Sasaki, Y., 2004, Three-dimensional inversion of static-shifted magnetotelluric data, Earth Planets Space, 56, 239-248. https://doi.org/10.1186/BF03353406
  14. Sasaki, Y., and Meju, M. A., 2006, Three-dimensional joint inversion for magnetotelluric resistivity and static shift distributions in complex media, J. Geophys. Res., 111, B05151.
  15. Siripunvaraporn, W., Egbert, G., Lenbury, Y., and Uyeshima, M., 2005, Three-dimensional magnetotelluric inversion: data-space method, Physics of the Earth and Planetary Interiors, 150, 3-14. https://doi.org/10.1016/j.pepi.2004.08.023
  16. Smith, J. T., 1996, Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator, Geophysics, 61, 1319-1324. https://doi.org/10.1190/1.1444055
  17. Sternberg, B. K., Washburne, J. C., and Pellerin, L., 1988, Correction for the static shift in magnetotellurics using transient electromagnetic soundings, Geophysics, 53, 1459-1468. https://doi.org/10.1190/1.1442426
  18. Takasugi, S., Tanaka, N., Kawakami, N., and Muramatsu, S., 1992, High spatial resolution of the resistivity structure revealed by a dense network MT measurement - A case study in the Minamikayabe area, Hokkaido, Japan, J. Geomag. Geoelectr., 44, 289-308. https://doi.org/10.5636/jgg.44.289
  19. Torres-Verdin, C., and Bostick, F. X., 1992, Principles of spatial surface electric fields filtering in magnetotellurics: Electromagnetic array profiling (EMAP), Geophysics, 57, 603-622. https://doi.org/10.1190/1.1443273
  20. Uchida, T., Lee, T. J., Honda, M., and Andan, A., 2002, 2-D and 3-D interpretation of magnetotelluric data in the Bajawa geothermal field, central Flores, Indonesia, Bull. Geol. Surv. Japan, 53, 265-283. https://doi.org/10.9795/bullgsj.53.265
  21. Unsworth, M. J., Travis, B. J., and Chave, A. D., 1993, Electromagnetic induction by a finite electric dipole source over a 2-D earth, Geophysics, 58, 198-214. https://doi.org/10.1190/1.1443406
  22. Yee, K. S., 1966, Numerical solution of initial boundary value problems involving Maxwell's equation in isotropic media, IEEE Trans. Anten. Prop., AP-14, 302-307.