Skip to main content
Article thumbnail
Location of Repository

Inversion of 2D NMR Data

By Christopher Bose


Schlumberger Limited is a multinational company supplying oilfield and nformation services to a worldwide energy market. These services include both exploration and production tools ranging through seismic and remote sensing, well-logging and reservoir optimization. The problem described in this report is related to well-logging via Nuclear Magnetic Resonance (NMR), a relatively new and developing tool with potential to reveal a range of reservoir properties including porosity and saturation, as well as physical properties of the petroleum deposit. In order to recover this information from NMR spectra the company must have an effective, efficient and robust algorithm to perform inversion from the dataset to the unknown probability distribution on magnetic relaxation times. This ill-posed problem is encountered in diverse areas of magnetic imaging and there does not appear to be an ‘off-the-shelf’ solution which the company can apply to its problem. Company scientists have developed a sophisticated algorithm which performs well on some simple test datasets, but they are interested in knowing if there are simpler approaches which could work effectively, or if some limited but useful properties of the density are accessible with a totally different approach. Our report is organised as follows. In Section 1.2 we present a careful and complete description of the problem and the work already done by the company. In Section 1.3 we discuss Truncated Singular Value Regularisation and Tikhonov Regularisation and show how some ‘off-the-shelf’ Matlab code may be used to good effect on the test datasets provided by the company. In Section 1.4 we show that one can incorporate higher order regularisation into the company’s existing algorithm, answering one specific question raised at the beginning of the workshop. Finally, in Section 1.5 we record our unsuccessful attempt to establish an iterative algorithm for the positively constrained inversion. Finally in the last section we review our conclusions and make suggestions for future work

Topics: Energy and utilities
Year: 2003
OAI identifier:

Suggested articles


  1. (1977). Algorithms for the regularization of ill-conditioned least squares problems,
  2. (1981). Estimating solutions of the first kind integral equations with nonnegative constraints and optimal smoothing,
  3. (1998). Kronecker product and SVD approximations in image restoration,
  4. (1994). Large least squares problems involving Kronecker products,
  5. (1962). Matrix iterative analysis, Prentice-Hall,
  6. (1991). Quantitative two-dimensional time correlation relaxometry,
  7. (2001). Regularization Tools, a Matlab Package for Analysis and Solution of Discrete Ill-Posed Problems.
  8. (2002). Solving Fredholm integrals of the first kind with tensor product structure in 2 and 2.5 dimensions,
  9. (1984). The Theory of Tikhonov Regularization for Fredholm Equations of the First Kind,
  10. (1991). Topics in Matrix Analysis,

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.