We describe a general global inversion methodology of multi-dimensional NMR logging for pore fluid typing and quantification in petroleum exploration. Although higher dimensions are theoretically possible, for practical reasons, we limit our discussion of proton density distributions as a function of two (2D) or three (3D) independent variables. The 2D can be diffusion coefficient and T-2 relaxation time (D-T-2), and the 3D can be diffusion coeflicient, T-2, and T-1 relaxation times (D-T-2-T-1) of the saturating fluids in rocks. Using the contrast between the diffusion coefficients of fluids (oil and water). the oil and water phases within the rocks call be clearly identified. This 2D or 3D proton density distribution function can be obtained from either two-window or regular type multiple CPMG echo trains encoded with diffusion, T-1, and T-2 relaxation by varying echo spacing and wait time. From this 2D/3D proton density distribution function, not only the saturations of water and oil can be determined, the viscosity of the oil and the gas-oil ratio can also be estimated based on a previously experimentally determined D-T-2 relationship. (C) 2004 Elsevier Inc. All rights reserved.