Determining bed boundaries from inversion of EM logging data using generalmeasures of model structure and data misfit

In this paper we propose a new method to locate bed boundaries by carrying out a 1-D nonlinear inversion of electromagnetic (EM) logging data. We firs t solve for 1-D resistivity structure in which the earth is modeled using l ayers of constant thickness. This thickness is determined based upon the to ol resolution and the desired resolution from the user. We use general meas ures of data misfit and model structure in the inversion to construct piece wise-constant models through the iteratively reweighted least-squares (IRLS ) procedure, and we minimize a generic global model objective function subj ect to data constraints. The general measure includes the traditional l(p) norm as a special case. The recovered piecewise-constant models simulate tr aditional models comprised of a few uniform layers and hence permit easy de termination of bed boundaries. The generic model objective function allows the incorporation of prior geological information and provides measures of the closeness to the reference model and the amount of structure in the rec overed model. We develop an efficient scheme to incorporate the IRLS proced ure into the search for the appropriate regularization parameter at each it eration. The iterations required by the IRLS algorithm are combined with th ose for the Linearized, iterative inversion procedure. The applicability of the algorithm is demonstrated using synthetic and field data examples.