Inversion for Non-Smooth Models with Physical Bounds

Geological processes produce structures at multiple scales. A discontinuity in the subsurface can occur due to layering, tectonic activities such as faulting, folding and fractures. Traditional approaches to invert geophysical data employ smoothness constraints. Such methods produce smooth models and thefore sharp contrasts in the medium such as lithological boundaries are not easily discernible. The methods that are able to produce non-smooth models, can help interpret the geological discontinuity. In this paper we examine various approaches to obtain non-smooth models from a finite set of noisy data. Broadly they can be categorized into approaches: (1) imposing non-smooth regularization in the inverse problem and (2) solve the inverse problem in a domain that provides multi-scale resolution, such as wavelet domain. In addition to applying non-smooth constraints, we further constrain the inverse problem to obtain models within prescribed physical bounds. The optimization with non-smooth regularization and physical bounds is solved using an interior point method. We demonstrate the applicability and usefulness of these methods with realistic synthetic examples and provide a field example from crosswell radar data