The presence of solid bitumen in a deep carbonate reservoir formation (Kuwait) occludes primary inter-crystalline porosity and thus impacts permeability. Here we describe a petrophysical model to both quantify in situ solid bitumen and estimate the effect of this bitumen on reservoir productivity. This is accomplished by adapting established nuclear magnetic resonance (NMR) workflows calibrated to laboratory NMR techniques on dolomite core samples with secondary organic matter or bitumen, primarily insoluble pyrobitumen.

The Density-NMR, "deficit porosity", method has been tested in the industry to identify the presence of solid organic (bitumen) matter. However, this method is not universally robust for quantifying bitumen saturation, as the magnitude of the "deficit porosity" is dependent on inter-echo time (TE), in situ temperature, and bitumen viscosity. In this study laboratory NMR techniques are utilized to calibrate the evaluation model. This was achieved through a combination of both one-dimensional (1D) T2 relaxation time and two-dimensional (2D) T1-T2 map measurements. The measurements were performed at 2 MHz (LF) for the purpose of comparing to NMR well logging, and at 23 MHz (HF) for higher sensitivity and better map resolutions, respectively.

The HF measurements with an inter-echo spacing time (TE) of 0.07 ms show that a solid-like component with its T2 < 0.1 ms can be detected at ambient temperature. In contrast, the LF maps or T2 distributions failed to detect such solid-like components due to the longer inter-echo spacing time (TE = 0.2 ms) used in the measurements. However, the LF signals are enhanced at the elevated temperatures, indicating contributions from the mobilized solid-like components. In order to interpret NMR logging, further data analysis was applied to the LF NMR T2 distributions at various temperatures by fitting the T2 spectra into a summation of three Gaussian functions. The results suggest that the two longer T2 components are temperature dependent following the Arrhenius relationship, while the shortest T2 component is temperature independent in the temperature range studied.

The recognition that LF NMR (comparable to the logging tool measurement) can detect the bitumen component at reservoir conditions implies that the "deficit porosity" method may underestimate the volume of bitumen in situ, as the shortest peak (at T2 ≈0.54 ms) is associated with the solid-like organic matter. Therefore, we introduce a workflow where "correct" the NMR (total) porosity for bitumen by identifying this short T2 bitumen peak in the (well-log) NMR T2 distribution and removing it by Gaussian peak fitting.

Lastly, a permeability model for bitumen prone intervals is presented based on the relative proportion of solid bitumen to original (or paleo-) inter-crystalline porosity. This model is calibrated to pyrolysis and organic petrology from the available core dataset. Pore filling solid bitumen generally reduces permeability by about one order of magnitude, relative to equivalent bitumen-free zones. The permeability reduction can be up to two orders of magnitude in the most bitumen-rich intervals.