Fast Calculation of the Sensitivity Matrix for Responses to the Earth's Conductivity: General Strategy and Examples

作者：Pankratov O V^*; Kuvshinov A V

来源：Izvestiya-Physics of the Solid Earth, 2010, 46(9): 788-804.

DOI：10.1134/S1069351310090089

Modern electromagnetic studies of the Earth deal with three-dimensional data sets and their inversion in the class of three-dimensional models. Due to the large scale of the 3D inverse problem, it is usually solved using gradient-type iteration techniques. These techniques require repeated calculations of the gradient of the penalty i.e., the sum of the misfit and the regularization term) with respect to the model parameters that describe the spatial distribution of conductivity. Since there are generally no analytical tools for calculating the misfit gradient, its computation requires the application of some numerical techniques. Unfortunately, it is impossible to perform mass calculations of the misfit gradient even using modern computing if one applies conventional numerical differentiation. In reality, in this case the number of calls of the forward problem is proportional to the number of the sought parameters, which is very large in case of a 3D inversion. However, there is a much more efficient method for calculation of the misfit gradient, namely, the adjoint-field technique. At present, this technique is widely used in the numerical schemes of the 3D inversion of electromagnetic and other data. With this technique we can calculate the misfit gradient with the same computational burden as that required for the solution of a single additional forward problem. Although this technique is widely used, we failed to find any comprehensive description of this technique in the literature that would allow one to apply this approach for any particular scheme of sounding, profiling, or transillumination as easily as a general formula is applied for a given case. In the present paper, we obtain general formulas for the quick calculation of the derivatives of the frequency-domain responses and the derivatives of the misfit function with respect to the variations in conductivity. We also show how the derived general formulas are transformed into particular expressions when we consider the response functions of the conventional and generalized magnetotelluric methods. In addition, we demonstrate how to apply our general methodology to the inversion of magnetovariational induction vectors (tippers), the phase tensor, and multipoint transfer functions, namely, the horizontal magnetic and telluric tensors.

• 出版日期2010-9