TY - JOUR

T1 - Joint two-dimensional dc resistivity and seismic traveltime inversion with cross-gradients constraints.

AU - Meju, Maxwell

AU - Gallardo, Luis A.

N1 - We produce the world's first objective way of mathematically coupling mechanical and electrical properties of the Earth in multi-physics based imaging, enabling the fusion of information taken from different measurement platforms. Our algorithm is now used in other research fields. First author was a PhD student, supervised by second author. RAE_import_type : Journal article RAE_uoa_type : Earth Systems and Environmental Sciences

PY - 2004/3/1

Y1 - 2004/3/1

N2 - It is now common practice to perform collocated DC resistivity and seismic refraction surveys that complement each other in the search for more accurate characterization of the subsurface. Although conventional separate DC resistivity and seismic models can be diagnostic, we posit that better results can be derived from jointly estimated models. We make the assumption that both methods must be sensing the same underlying geology and have developed an innovative resistivity-velocity cross-gradients relationship to evaluate the structural features common to both methods. The cross-gradients function is incorporated as a constraint in a nonlinear least squares problem formulation, which is solved using the Lagrange multiplier method. The resultant iterative two-dimensional (2-D) joint inversion scheme is successfully applied to synthetic data (serving as validation tests here) and to field data from collocated DC resistivity and seismic refraction profiling experiments and also compared to conventional separate inversion results. The joint inversion results are shown to be superior to those from separate 2-D inversions of the respective data sets, since our algorithm leads to resistivity and velocity models with remarkable structural agreement.

AB - It is now common practice to perform collocated DC resistivity and seismic refraction surveys that complement each other in the search for more accurate characterization of the subsurface. Although conventional separate DC resistivity and seismic models can be diagnostic, we posit that better results can be derived from jointly estimated models. We make the assumption that both methods must be sensing the same underlying geology and have developed an innovative resistivity-velocity cross-gradients relationship to evaluate the structural features common to both methods. The cross-gradients function is incorporated as a constraint in a nonlinear least squares problem formulation, which is solved using the Lagrange multiplier method. The resultant iterative two-dimensional (2-D) joint inversion scheme is successfully applied to synthetic data (serving as validation tests here) and to field data from collocated DC resistivity and seismic refraction profiling experiments and also compared to conventional separate inversion results. The joint inversion results are shown to be superior to those from separate 2-D inversions of the respective data sets, since our algorithm leads to resistivity and velocity models with remarkable structural agreement.

U2 - 10.1029/2003JB002716

DO - 10.1029/2003JB002716

M3 - Journal article

VL - 109

SP - B03311

JO - Journal of Geophysical Research: Solid Earth

JF - Journal of Geophysical Research: Solid Earth

IS - (B3)

ER -