최근에 수행되고 있는 중력 탐사는 고분해능의 중력계와 GPS(Global positioning system)를 통한 정밀한 측정과 측지가 이루어지고 있다. 중력탐사에서 모델링과 역산의 기술은 많은 발전이 있어왔지만, 중력자료처리는 거의 변화가 없었다. 통상적인 정밀한 중력 자료 보정을 통한 부우게이상은 측정점의 고도에서 기준면까지의 물질의 영향을 일정한 밀도를 이용해 제거해 버리기 때문에 측정점 바로 하부의 이상체에 의한 영향을 상당히 왜곡시키게 된다. 본 연구에서는 탐사 지역의 지형을 DEM(Digital Elevation Map) 자료와 Multiquadric equation을 이용하여 실제 지형과 유사한 Multiquadric surface를 자동적으로 구성하고, 이를 블록화 함으로써 보정의 대상이었던 기준면 상부에 대한 밀도를 탐사 지역의 지질 정보와 지형을 포함하는 역산을 통해 수치적으로 계산하였다. 이러한 지형을 포함한 역산 방법을 3차 원중력지형역산(3DGTI; 3-D Gravity Terrain Inversion)이라 한다. 이 연구의 효율성을 검증하기 위하여 주변암과 밀도차가 존재하는 관입지역에 대한 모델을 구성하고 적용한 결과 기존의 부게 보정 방법을 적용한 부게 이상도에 비해 자료의 왜곡이 감소하는 효과를 얻을 수 있었다. 이를 통하여 지형 역산을 통한 객관적인 부게 밀도의 결정과 부게 보정시 실제의 수평적인 밀도 변화를 반영함으로써 기존의 문제점을 보완하였다. 게다가, 3DGTI로부터 얻어진 밀도분포는 지형의 윤곽을 그대로 표현하고 있어서 보다 실질적인 지질을 보여준다고 하겠다. 이 방법을 화강암체가 관입하고 있는 마산·창원 일대에서의 중력 탐사 자료에 적용해본 결과 기존 방법보다 관입 화강암체의 위치와 그 규모를 알아내는데 더 효과적이었다. 따라서, 수평적인 밀도 변화가 뚜렷하게 존재하는 지역의 경우, 새로운 중력 자료 처리 방법이 기존의 부게 보정에서 발생하였던 문제점을 해결함으로써 천부의 분해능을 높이고, 심부의 밀도 분포도 좀더 정확하게 계산할 수 있으리라 생각된다.

Bayesian inversion for gravity and resistivity data was performed to investigate the cavity structure appearing as a lava tunnel in Cheju Island, Korea. Dipole-dipole DC resistivity data were proposed for a prior information of gravity data and we applied the geostatistical techniques such as kriging and simulation algorithms to provide a prior model information and covariance matrix in data domain. The inverted resistivity section gave the indicator variogram modeling for each threshold and it provided spatial uncertainty to give a prior PDF by sequential indicator simulations. We also presented a more objective way to make data covariance matrix that reflects the state of the achieved field data by geostatistical technique, cross-validation. Then Gaussian approximation was adopted for the inference of characteristics of the marginal distributions of model parameters and Broyden update for simple calculation of sensitivity matrix and SVD was applied. Generally cavity investigation by geophysical exploration is difficult and success is hard to be achieved. However, this exotic multiple interpretations showed remarkable improvement and stability for interpretation when compared to data-fit alone results, and suggested the possibility of diverse application for Bayesian inversion in geophysical inverse problem.