오일러 디컨벌루션을 사전정보로 이용한 3 차원 중력 역산

3D gravity inversion with Euler deconvolution as a priori information

  • 임형래 (한국지질자원연구원 지질지구정보부) ;
  • 박영수 (한국지질자원연구원 지질지구정보부) ;
  • 임무택 (한국지질자원연구원 지질지구정보부) ;
  • 구성본 (한국지질자원연구원 지질지구정보부) ;
  • 권병두 (서울대학교 사범대학 지구과학교육과)
  • Rim, Hyoung-Rae (Korea Institute of Geoscience and Mineral Resources) ;
  • Park, Yeong-Sue (Korea Institute of Geoscience and Mineral Resources) ;
  • Lim, Mu-Taek (Korea Institute of Geoscience and Mineral Resources) ;
  • Koo, Sung-Bon (Korea Institute of Geoscience and Mineral Resources) ;
  • Kwon, Byung-Doo (Department of Earth Sciences Education, Seoul National University)
  • 발행 : 2007.02.28

초록

고해상도를 가지는 지하 밀도 영상을 얻기 위한 3 차원 중력 역산은 모델 변수들이 급격하게 많아지는 문제가 발생한다. 이 논문에서는 모델 변수들의 수를 줄이기 위해서 오일러 디컨벌루션의 해를 사전정보로 활용하는 3 차원 중력역산을 제안하였다. 이 논문에서 고안한 역산 알고리즘의 핵심은 오일러 디컨벌루션의 해가 얻어진 주위로 역산 공간을 제한하여 역산 해의 비유일성을 줄인 점이다. 먼저 중력 자료에 대한 3 차원 오일러 디컨별루션의 해를 구하고, 오일러 디컨벌루션의 해가 나타나는 주위에서만 3 차원 확장 탐색 역산을 수행하여 지하 밀도 영상을 구하였다. 이 3 차원 중력 역산 방법은 합성 모델에 적용하여 그 성능을 검증하였고, 석회암 지대에 존재하는 공동의 분포를 밝히기 위한 고정밀 중력탐사 자료 역산에도 적용하였다. 결과적으로, 오일러 디컨벌루션의 해를 사전정보로 이용한 역산을 이용하여 분해능이 향상된 고해상도의 지하 멸도 영상을 구할 수 있었다.

It is difficult to obtain high-resolution images by 3D gravity inversion, because the problem is extremely underdetermined - there are too many model parameters. In order to reduce the number of model parameters we propose a 3D gravity inversion scheme utilising Euler deconvolution as a priori information. The essential point of this scheme is the reduction of the nonuniqueness of solutions by restricting the inversion space with the help of Euler deconvolution. We carry out a systematic exploration of the growing body process, but only in the restricted space within a certain radius of the Euler solutions. We have tested our method with synthetic gravity data, and also applied it to a real dataset, to delineate underground cavities in a limestone area. We found that we obtained a more reasonable subsurface density image by means of this combination between the Euler solution and the inversion process.

키워드

참고문헌

  1. Bishop, I., Styles, P., Emsley, S. J., and Ferguson, N. S., 1997, The detection of cavities using microgravity technique; case histories from mining and karstic environments: in Modern Geophysics in Engineering Geology: McCann, D. M., Eddleston, M., Fenning, P. J., and Reeves, G. M. (eds). The Geological Society, 153-166
  2. Butler, D. K., 1984, Microgravimetric and gravity gradient techniques for detection of subsurface cavities: Geophysics 49, 1084-1096. doi: 10.1190/1.1441723
  3. Camacho, A. G., Montesinos, F. G., and Vieira, R., 2000, Gravity inversion by means of growing bodies: Geophysics 65, 95-101. doi: 10.1190/1.1444729
  4. Colley, G. C., 1963, The detection of caves by gravity measurements: Geophysical Prospecting 11, 1-9. doi: 10.1111/j.1365-2478.1963. tb02019.x
  5. Last, B. J., and Kubik, K., 1983, Compact gravity inversion: Geophysics 48, 713-721. doi: 10.1190/1.1441501
  6. Oldenburg, D. W., 1974, The inversion and interpretation of gravity anomalies: Geophysics 39, 526-536. doi: 10.1190/1.1440444
  7. Pedersen, L. B., 1979, Constrained inversion of potential field data: Geophysical Prospecting 27, 726-748. doi: 10.1111/j.1365-2478. 1979.tb00993.x
  8. Rim, H., Park, Y.-S., Lim, M., Koo, S. B., and Kwon, B., 2005, Application of microgravity survey for imaging the density distribution of the rock fill dam itself: 67th European Association of Geoscientists & Engineers Conference and Exhibition, Expanded Abstracts, P354
  9. Styles, P., and Thomas, E., 2001, The use of microgravity for characterization of karstic cavities on Grand Bahama, Bahamas: in Geotechnical and Environmental Applications of Karstic Geology and Hydrology: Beck, B. F. and Herring, J. G. (eds). Balkema, 389-394
  10. Thompson, D. T., 1982, EULDPH: A new technique for making computerassisted depth estimates from magnetic data: Geophysics 47, 31-37. doi: 10.1190/1.1441278