Multiple level-set joint inversion of traveltime and gravity data with application to ore delineation: A synthetic study

We have developed a multiple level-set method for simultaneous inversion of gravity and seismic traveltime data. The method recovers the boundaries between regions with distinct physical properties assumed constant and known, creating structurally consistent models of two subsurface properties: P-wave velocity and density. In single level-set methods, only two rock units can be considered: background and inclusion. Such methods have been applied to examples representing various geophysical scenarios, including in the context of joint inversion. In multiple level-set methods, several units can be considered, which make them far more applicable to real earth scenarios. Recently, a multiple level-set method has been proposed for inversion of magnetic data. We extend the multiple level-set formulation to joint inversion of gravity and traveltime data, improving upon previous work, and we investigate applicability of such an inversion method in ore delineation. In mineral exploration environments, traditional seismic imaging and inversion methods are challenging because of the small target size and the specific physical property contrasts involved. First-arrival seismic traveltime and gravity data complement each other, and we found that joint multiple level-set inversion is more beneficial than separate inversions, especially with limited data and slow targets. Our method is more robust than the joint inversion method based on clustering of physical properties in recovery of piecewise homogeneous models not well-constrained by the data. To justify the known property assumption, we found the degree of robustness of the multiple level-set joint inversion to uncertainties arising from incomplete knowledge of small-scale subsurface physical property variations and target composition.