S U M M A R Y Non-uniqueness in the geophysical inverse problem is well recognized and so too is the ability to obtain solutions with different character by altering the form of the regularization function. Of particular note is the use of p norms with p ∈ [0, 2] which gives rise to sparse or smooth models. Most algorithms are designed to implement a single p norm for the entire model domain. This is not adequate when the fundamental character of the model changes throughout the volume of interest. In such cases we require a generalized regularization function where each sub-volume of the model domain has penalties on smallness and roughness and its own suite of p parameters. Solving the inverse problem using mixed p norms in the regularization (especially for p < 1) is computationally challenging. We use the Lawson formulation for the p norm and solve the optimization problem with Iterative Reweighted Least Squares. The algorithm has two stages; we first solve the l2-norm problem and then we switch to the desired suite of p norms; there is one value of p for each term in the objective function. To handle the large changes in numerical values of the regularization function when p values are changed, and to ensure that each component of the regularization is contributing to the final solution, we successively rescale the gradients in our Gauss–Newton solution. An indicator function allows us to evaluate our success in finding a solution in which components of the objective function have been equally influential. We use our algorithm to generate an ensemble of solutions with mixed p norms. This illuminates some of the non-uniqueness in the inverse problem and helps prevent overinterpretation that can occur by having only one solution. In addition, we use this ensemble to estimate the suite of p values that can be used in a final inversion. First, the most common features of our ensemble are extracted using principal component analysis and edge detection procedures; this provides a reference model. A correlation of each member of the ensemble with the reference model, carried out in a windowed domain, then yields a set of p values for each model cell. The efficacy of our technique is illustrated on a synthetic 2-D cross-well example. We then apply our technique to the field example that motivated this research, the 3-D inversion of magnetic data at a kimberlite site in Canada. Since the final regularization terms have different sets of p values in different regions of model space we are able to recover compact regions associated with the kimberlite intrusions, continuous linear features with sharp edges that are associated with dykes and a background that is relatively smooth. The result has a geologic character that would not have been achievable without the use of spatially variable mixed norms.

Dominique Fournier

Douglas W Oldenburg

Geophysical Journal International

《国际地球物理杂志》发表关于理论、计算、应用和观测地球物理学各个方面的高质量研究论文、特快信函、特邀评论论文和书评。

The Geophysical Journal International publishes high-quality research papers, rapid communications, invited reviews, and book reviews covering all aspects of theoretical, computational, applied, and observational geophysics.

《Geophysical Journal International》发表高质量的研究论文、特快信函、特邀评论和书评，涵盖理论、计算、应用和观测地球物理学的各个方面。

Category	Quartile
地球科学	3
地球科学, 地球化学与地球物理	3

Inversion using spatially variable mixed ℓp norms

ivySCI AI Smartly Parses PDF, Answers Researchers' Questions, and Helps You Understand Papers in Seconds

Journal Info

Category	Quartile
GEOCHEMISTRY & GEOPHYSICS	2