Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging and inversion

摘要

Linearized inversion of surface seismic data for a model of the earth's subsurface requires estimating the sensitivity of the seismic response to perturbations in the earth's subsurface. This sensitivity, or Jacobian, matrix is usually quite expensive to estimate for all but the sim- plest model parameterizations. We exploit the numerical structure of the finite-element method, modern sparse matrix technology, and source-receiver reciprocity to de- velop an algorithm that explicitly calculates the Jacobian matrix at only the cost of a forward model solution.