Three-dimensional gravity inversion using automatic regional-residual components separation with application to the Iwaki area in Fukushima Prefecture, Japan

Abstract

We have developed a simple method for analysing gravity data to define 3D subsurface structures consistent with geological maps and reflection seismic survey results. This new method makes it possible to automatically separate gravity anomalies caused by structures that are the target of analysis (residual gravity) and not subject to gravity anomaly analysis (trend surfaces or regional gravity variations). This separation is achieved by assuming that the gravity anomaly trend surface is extremely smooth prior to analysis. In addition, this analysis can be carried out automatically, so it is possible to use a boot-strap simulation to evaluate the effect that gravity data accuracy and survey point distribution have on the analysis results. Using this newly developed method, we conducted a 3D subsurface structural analysis based on newly obtained gravity survey data in an area around Iwaki City, Fukushima Prefecture. The results of the 3D structural analysis revealed sharp depressions in the basement rock and basin-shaped underground structures along the Yunodake and Yamada faults. Finally, we were able to estimate the influence of survey point distribution and abrupt changes in the structure on the analysis results using a boot-strap simulation of the 3D analysis.