Characterization of porous media by T2-T2 correlation beyond fast diffusion limit

Pore size distribution and surface relaxivity are two important properties of porous media such as rock samples and can be obtained by NMR methods. However, it is difficult to obtain these information beyond the fast diffusion limit. Here we present a new method to directly characterize the averaged pore size of a porous sample with a narrow pore size distribution. This method is based on the parallel plates pore model and the T₂-T₂ correlation sequence. The pore size (a) - surface relaxivity (ρ) correlation maps were obtained using the non-negative least squares method. Three kinds of glass bead samples were measured and the averaged pore size and surface relaxivity were extracted.

Measurement of interfacial area from NMR time dependent diffusion and relaxation measurements

The interfacial area between two immiscible phases in porous media is an important parameter for describing and predicting 2 phase flow. Although present in several models, experimental investigations are sparse due to the lack of appropriate measurement techniques. We propose two NMR techniques for the measurement of oil-water interfacial area: (i) a time dependent NMR diffusion technique applicable in static conditions, similar to those used for the measurement of the solid specific surface of a porous media, and (ii) a fast relaxation technique applicable in dynamic conditions while flowing, based on an interfacial relaxation mechanism induced by the inclusion of paramagnetic salts in the water phase. For dodecane relaxing on doped water, we found an oil interfacial relaxivity of 1.8μm/s, large enough to permit the measurement of specific interfacial surface as small as 1000cm²/cm³. We demonstrate both NMR techniques in drainage followed by imbibition, in a model porous media with a narrow pore size distribution. While flowing, we observe that the interfacial area is larger in imbibition than in drainage, implying a different organization of the oil phase. In a carbonate sample with a wide pore size distribution, we evidence the gradual invasion of the smallest pores as the oil-water pressure difference is increased.

Detection of intermolecular homonuclear dipolar coupling in organic rich shale by transverse relaxation exchange

The mechanism behind surface relaxivity within organic porosity in shales has been an unanswered question. Here, we present results that confirm the existence of intermolecular homonuclear dipolar coupling between solid and liquid phases in sedimentary organic matter. Transverse magnetization exchange measurements were performed on an organic-rich shale saturated with liquid hydrocarbon. Liquid and solid constituents were identified through both sample resaturation and through their T₁/T₂ ratios. Extensive cross peaks are observed in the T₂-T₂ exchange spectra between the solid and liquid constituents, indicating an exchange of magnetization between the two phases. This result cannot arise from physical molecular diffusion, and the dissolution energies are too high for chemical exchange, such that the magnetization exchange must arise from intermolecular homonuclear dipolar coupling. These results both confirm a possible source of surface relaxivity in organic matter and emphasize caution in the use of standard porous media interpretations of relaxation results in shales because of coupling between different magnetization environments.

Evaluation of sandstone surface relaxivity using laser-induced breakdown spectroscopy

Nuclear magnetic resonance (NMR) relaxometry is a common technique used to assess the pore size of fluid-filled porous materials in a wide variety of fields. However, the NMR signal itself only provides a relative distribution of pore size. To calculate an absolute pore size distribution from the NMR data, the material's surface relaxivity needs to be known. Here, a method is presented using laser-induced breakdown spectroscopy (LIBS) to evaluate surface relaxivity in sandstones. NMR transverse and longitudinal relaxation was measured on a set of sandstone samples and the surface relaxivity was calculated from the pore size distribution determined with MICP measurements. Using multivariate analysis, it was determined that the LIBS data can predict with good accuracy the longitudinal (R²∼0.84) and transverse (R²∼0.79) surface relaxivity. Analysis of the regression coefficients shows significant influence from several elements. Some of these are elements previously established to have an effect on surface relaxivity, such as iron and manganese, while others are not commonly associated with surface relaxivity, such as cobalt and titanium. Furthermore, LIBS provides advantages compared to current methods to calibrate surface relaxivity in terms of speed, portability, and sample size requirements. While this paper focuses on geological samples, the method could potentially be expanded to other types of porous materials.

Direct correlation of diffusion and pore size distributions with low field NMR

The time-dependent diffusion coefficient (D) is a powerful tool to probe microstructure in porous media, and can be obtained by the NMR method. In a real porous sample, molecular diffusion is very complex. Here we present a new method which directly measures the relationship between effective diffusion coefficients and pore size distributions without knowing surface relaxivity. This method is used to extract structural information and explore the relationship between D and a in porous media having broad pore size distributions. The diffusion information is encoded by the Pulsed Field Gradient (PFG) method and the pore size distributions are acquired by the Decay due to Diffusion in the Internal Field (DDIF) method. Two model samples were measured to verify this method. Restricted diffusion was analyzed, and shows that most fluid molecules experience pore wall. The D(a) curves obtained from correlation maps were fitted to the Padé approximant equation and a good agreement was found between the fitting lines and the measured data. Then a sandstone sample with unknown structure was measured. The state of confined fluids was analyzed and structural information, such as pore size distributions, were extracted. The D - T1 correlation maps were also obtained using the same method, which yielded surface relaxivities for different samples. All the experiments were conducted on 2MHz NMR equipment to obtain accurate diffusion information, where internal gradients can be neglected. This method is expected to have useful applications in the oil industry, particularly for NMR logging in the future.

Robust determination of surface relaxivity from nuclear magnetic resonance DT(2) measurements

Nuclear magnetic resonance (NMR) is a powerful tool to probe into geological materials such as hydrocarbon reservoir rocks and groundwater aquifers. It is unique in its ability to obtain in situ the fluid type and the pore size distributions (PSD). The T1 and T2 relaxation times are closely related to the pore geometry through the parameter called surface relaxivity. This parameter is critical for converting the relaxation time distribution into the PSD and so is key to accurately predicting permeability. The conventional way to determine the surface relaxivity ρ2 had required independent laboratory measurements of the pore size. Recently Zielinski et al. proposed a restricted diffusion model to extract the surface relaxivity from the NMR diffusion-T2 relaxation (DT2) measurement. Although this method significantly improved the ability to directly extract surface relaxivity from a pure NMR measurement, there are inconsistencies with their model and it relies on a number of preset parameters. Here we propose an improved signal model to incorporate a scalable LT and extend their method to extract the surface relaxivity based on analyzing multiple DT2 maps with varied diffusion observation time. With multiple diffusion observation times, the apparent diffusion coefficient correctly describes the restricted diffusion behavior in samples with wide PSDs, and the new method does not require predetermined parameters, such as the bulk diffusion coefficient and tortuosity. Laboratory experiments on glass beads packs with the beads diameter ranging from 50 μm to 500 μm are used to validate the new method. The extracted diffusion parameters are consistent with their known values and the determined surface relaxivity ρ2 agrees with the expected value within ±7%. This method is further successfully applied on a Berea sandstone core and yields surface relaxivity ρ2 consistent with the literature.