Regularized extremal bounds analysis (REBA): An approach to quantifying uncertainty in nonlinear geophysical inverse problems

[1] Geophysical measurements are band-limited in nature, contain noise, and bear a nonlinear relationship to the subterranean features being sampled. These result in uncertainty in data interpretation and necessitate a regularization approach. Although model uncertainty and non-uniqueness can be reduced by combining measurements of fundamentally different physical attributes of a subsurface target under investigation or by using available a priori information about the target, quantifying non-uniqueness remains a difficult task in nonlinear geophysical inversion. This paper develops a new theoretical framework for regularized extremal bounds analysis (REBA) of model uncertainty using the most-squares formalism. The nonlinear most-squares method allows for combining data constraints and their associated uncertainties in an objective manner to determine the model bounds. By making the assumption that the relevant bounding models must sample the same underlying geology, it is proposed that structural dissimilarity between the extreme models can be quantified using the cross-products of the gradients of their property fields and should serve as an objective measure of interpretational uncertainty in multidimensional inversion.

[1]  N. Linde,et al.  Local earthquake (LE) tomography with joint inversion for P‐ and S‐wave velocities using structural constraints , 2006 .

[2]  M. A. Meju,et al.  Iterative most-squares inversion: application to magnetotelluric data , 1992 .

[3]  Thomas Kalscheuer,et al.  A non-linear truncated SVD variance and resolution analysis of two-dimensional magnetotelluric models , 2007 .

[4]  D. Jackson The use of a priori data to resolve non‐uniqueness in linear inversion , 1979 .

[5]  Max A. Meju,et al.  Biased estimation: a simple framework for inversion and uncertainty analysis with prior information , 1994 .

[6]  A. Binley,et al.  Improved hydrogeophysical characterization using joint inversion of cross‐hole electrical resistance and ground‐penetrating radar traveltime data , 2006 .

[7]  S. Hubbard,et al.  Joint inversion of crosshole radar and seismic traveltimes , 2008 .

[8]  M. Meju,et al.  Characterization of heterogeneous near‐surface materials by joint 2D inversion of dc resistivity and seismic data , 2003 .

[9]  Roel Snieder Uncertainty Estimation in Inverse Problems; an Uncertain Proposition , 2004 .

[10]  D. W. Vasco,et al.  Invariance, groups, and non-uniqueness; the discrete case , 2005 .

[11]  V. Sakkas,et al.  Heterogeneous crust and upper mantle across southern Kenya and the relationship to surface deformation as inferred from magnetotelluric imaging , 2007 .

[12]  M. Meju,et al.  Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints , 2004 .

[13]  Max A. Meju,et al.  Joint two‐dimensional cross‐gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification , 2007 .