Tortuosity measurement and the effects of finite pulse widths on xenon gas diffusion NMR studies of porous media.

We have extended the utility of NMR as a technique to probe porous media structure over length scales of approximately 100-2000 microm by using the spin 1/2 noble gas 129Xe imbibed into the system's pore space. Such length scales are much greater than can be probed with NMR diffusion studies of water-saturated porous media. We utilized Pulsed Gradient Spin Echo NMR measurements of the time-dependent diffusion coefficient, D(t), of the xenon gas filling the pore space to study further the measurements of both the pore surface-area-to-volume ratio, S/V(p), and the tortuosity (pore connectivity) of the medium. In uniform-size glass bead packs, we observed D(t) decreasing with increasing t, reaching an observed asymptote of approximately 0.62-0.65D(0), that could be measured over diffusion distances extending over multiple bead diameters. Measurements of D(t)/D(0) at differing gas pressures showed this tortuosity limit was not affected by changing the characteristic diffusion length of the spins during the diffusion encoding gradient pulse. This was not the case at the short time limit, where D(t)/D(0) was noticeably affected by the gas pressure in the sample. Increasing the gas pressure, and hence reducing D(0) and the diffusion during the gradient pulse served to reduce the previously observed deviation of D(t)/D(0) from the S/V(p) relation. The Pade approximation is used to interpolate between the long and short time limits in D(t). While the short time D(t) points lay above the interpolation line in the case of small beads, due to diffusion during the gradient pulse on the order of the pore size, it was also noted that the experimental D(t) data fell below the Pade line in the case of large beads, most likely due to finite size effects.     

2D NMR技术在石油测井中的应用

近几年,2D NMR技术得到迅速发展,特别是在核磁共振测井领域. 该文将主要介绍2D NMR技术的脉冲序列、弛豫原理以及2D NMR技术在石油测井中应用. 2D NMR技术是在梯度场的作用下,利用一系列回波时间间隔不同的CPMG脉冲进行测量,利用二维的数学反演得到2D NMR. 2D NMR技术可以直接测量自扩散系数、弛豫时间、原油粘度、含油饱和度、可动水饱和度、孔隙度、渗透率等地层流体性质和岩石物性参数. 从2D NMR谱上,可以直观的区分油、气、水,判断储层润湿性,确定内部磁场梯度等. 2D NMR技术为识别流体类型提供了新方法.     

基于ADI-FDTD模拟的岩石核磁共振横向弛豫特征分析

将复杂的骨架-孔隙系统抽象成等效双孔介质,根据Bloch方程构建数学模型,用交替隐式时域有限差分(ADI-FDTD)和联合反演迭代法(SIRT)进行横向宏观磁化矢量的数值模拟与核磁共振T₂谱的反演,定量研究扩散系数、弛豫速率、孔隙组分比和孔隙宽度对核磁响应的影响.结果表明:横向宏观磁化矢量衰减速率与扩散系数和微孔隙分量成正比,与孔隙宽度成反比,与表面弛豫速率基本无关.当扩散系数较大、孔隙宽度较小时,核磁共振T₂谱难以直观反映孔隙组分及孔隙结构.应用核磁共振评价孔隙结构时需特别注意扩散系数和孔隙尺寸的影响.     

NMR in well logging and hydrocarbon exploration

Nuclear magnetic resonance (NMR) has become a versatile tool for the evaluation of underground hydrocarbon reservoirs. Formation attributes such as rock porosity and rock pore size distributions, as well as the relative concentrations of water, oil and gas, can be inferred from subsurface NMR. The hydrogen NMR signal encodes porosity as amplitude, pore sizes as relaxation times and fluid properties as a mixture of relaxation and diffusion rates. The paper describes the basic operating principles for NMR on cable (wireline), NMR on a drill string (logging-while-drilling) and NMR for downhole fluid sampling. The geometry of the borehole requires a magnet that projects its field into the surrounding rock, implying a grossly inhomogeneous field distribution. Experience shows that even under these circumstances, saturation-recovery and Carr-Purcell-Meiboom-Gill sequences can work well and yield meaningfulT ₁ andT ₂ information.     

基于低场核磁共振的抚顺油页岩孔隙连通性演化研究

油页岩原位注热开采过程中，储层内部孔隙结构的连通性直接影响载热介质的流动行为和传热效率，同时对油气产物的扩散和流动行为起控制作用.本文利用低场核磁共振（LF NMR）技术，考察了不同热解终温（23~650℃）处理时，饱和水及束缚水状态下抚顺油页岩的T₂谱，分析了可动流体T₂截止值、束缚流体孔隙度、饱和流体孔隙度、渗透率等NMR孔隙参数，定量研究了随热解终温升高，抚顺油页岩孔隙结构的连通性演化规律.研究结果表明热解终温对抚顺油页岩孔隙连通性及渗透率的变化起控制作用，且可动流体孔隙度对总孔隙度的增加起主要促进作用，这说明热解终温升高加大了渗透率及油气产物的输运能力.本文为深入认识油页岩原位热解过程中孔隙结构的演化提供了依据.     

The narrow pulse approximation and long length scale determination in xenon gas diffusion NMR studies of model porous media

We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Padé approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Padé interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Padé length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter.     

A global inversion method for multi-dimensional NMR logging

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 coefficient, 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 can 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.     

In situ fluid typing and quantification with 1D and 2D NMR logging

In situ nuclear magnetic resonance (NMR) fluid typing has recently gained momentum due to data acquisition and inversion algorithm enhancement of NMR logging tools. T(2) distributions derived from NMR logging contain information on bulk fluids and pore size distributions. However, the accuracy of fluid typing is greatly overshadowed by the overlap between T(2) peaks arising from different fluids with similar apparent T(2) relaxation times. Nevertheless, the shapes of T(2) distributions from different fluid components are often different and can be predetermined. Inversion with predetermined T(2) distributions allows us to perform fluid component decomposition to yield individual fluid volume ratios. Another effective method for in situ fluid typing is two-dimensional (2D) NMR logging, which results in proton population distribution as a function of T(2) relaxation time and fluid diffusion coefficient (or T(1) relaxation time). Since diffusion coefficients (or T(1) relaxation time) for different fluid components can be very different, it is relatively easy to separate oil (especially heavy oil) from water signal in a 2D NMR map and to perform accurate fluid typing. Combining NMR logging with resistivity and/or neutron/density logs provides a third method for in situ fluid typing. We shall describe these techniques with field examples.     

Tracking pore to pore exchange using relaxation exchange spectroscopy

We observe the movement of water over time between pores of differing sizes in Castlegate sandstone. To achieve this, we perform an NMR transverse relaxation exchange experiment for several mixing times. The resulting data are converted to 2D T2 distributions using a 2D inverse Laplace transform (ILT). We show for the first time that quantitative analysis of ILT distributions enables one to extract characteristic times for different pores sizes. This information is potentially useful for permeability determination as well as better understanding of exchange between specific pore subpopulations.     

Xenon NMR measurements of permeability and tortuosity in reservoir rocks

In this work we present measurements of permeability, effective porosity and tortuosity on a variety of rock samples using NMR/MRI of thermal and laser-polarized gas. Permeability and effective porosity are measured simultaneously using MRI to monitor the inflow of laser-polarized xenon into the rock core. Tortuosity is determined from measurements of the time-dependent diffusion coefficient using thermal xenon in sealed samples. The initial results from a limited number of rocks indicate inverse correlations between tortuosity and both effective porosity and permeability. Further studies to widen the number of types of rocks studied may eventually aid in explaining the poorly understood connection between permeability and tortuosity of rock cores.     

Estimation of Permeability by Integrating Nuclear Magnetic Resonance (NMR) Logs with Mercury Injection Capillary Pressure (MICP) Data in Tight Gas Sands

It has been a great challenge to determine permeability in tight gas sands due to the generally poor correlation between porosity and permeability. The Schlumberger Doll Research (SDR) and Timur–Coates permeability models, which have been derived for use with nuclear magnetic resonance (NMR) data, also lose their roles. In this study, based on the analysis of the mercury injection experiment data for 20 core plugs, which were drilled from tight gas sands in the Xujiahe Formation of central Sichuan basin, Southwest China, two empirical correlations between the pore structure index ( $ \sqrt {{K \mathord{\left/ {\vphantom {K \varphi }} \right. \kern-0em} \varphi }} $ , defined by the square root of the ratio of rock permeability and porosity) and the R ₃₅ (the pore throat radius corresponding to 35.0 % of mercury injection saturation), the pore structure index and the Swanson parameter have been developed. To consecutively estimate permeability in field applications, based on the study of experimental NMR measurements for 36 core samples, two effective statistical models, which can be used to derive the Swanson parameter and R ₃₅ from the NMR T ₂ logarithmic mean value, have been established. These procedures carried out on the experimental data set can be extended to reservoir conditions to estimate consecutive formation permeability along the intervals with which NMR logs were acquired. The processing results of several field examples using the proposed technique show that the classification scale models are effective only in tight gas reservoirs, whereas the SDR and Timur–Coates models are inapplicable. The R ₃₅-based model is of significance in thin sands with high porosity and high permeability, but the predicted permeability curves in tight gas sands are slightly lower. In tight gas and thin sands, the Swanson parameter model is all credible.     

活性炭吸附甲醛的格子Boltzmann模拟

基于孔隙尺度,结合活性炭与甲醛的真实物性参数,利用格子Boltzmann方法,选取热质耦合的LBGK模型对填充有球形活性炭的方腔内部双扩散混合对流、流固共轭传热及吸附特性进行数值模拟.分别采用二维D2Q9模型描述速度温度场,D2Q5模型描述浓度场,研究活性炭颗粒直径、孔隙率以及颗粒的排列方式对整个动态吸附性能的影响.结...     

Feasibility of diffusion-NMR surface-to-volume measurements tested by calculations and computer simulations

It has been demonstrated previously that the surface-to-volume ratio S/V can be determined from the derivative of the time-dependent diffusion coefficient D(t), in the limit t --> 0. Several questions arise concerning the practicality of determining S/V by NMR. In particular, how large are the errors generated by (1) working outside the t --> 0 limit and (2) measuring D outside the b --> 0 limit, both for narrow and full-width gradient pulses? Here b is gamma2G2delta2Delta for narrow pulses and gamma2G2t3/12 for broad pulses. These questions are addressed by random-walk computer simulations and numerical calculations in geometries relevant to small-airways of lung. The results demonstrate that one can work well outside the t --> 0 and b --> 0 limits, provided 10-20% accuracy in the measured S/V is sufficient. Emphasis is placed on the useful range of times t for which NMR determinations of lung S/V are feasible.     

页岩气中CH4-C2H6吸附-扩散分子动力学研究

本文采用巨正则蒙特卡洛(GCMC)和分子动力学(MD)模拟方法,对比分析了不同温度、压力和孔径对二元气体(CH4-C2H6)在K-伊利石中的吸附-扩散的影响.结果表明,在低压条件下,K-伊利石对C2H6的吸附能力大于CH4, C2H6优先吸附在K-伊利石孔隙表面.热力学因子随着孔径的增加而减小,C2H6的热力学因子大于...     

Concentration-dependent self-diffusion of adsorbates in mesoporous materials

The pulsed-field gradient NMR method has been applied to study self-diffusion of liquids in mesoporous materials with different pore sizes and morphologies as a function of pore loading. It is found that the effective diffusivities of adsorbate molecules in mesopores at partial loadings are related to two mechanisms, the Knudsen diffusion through the gaseous phase in the pore space and the diffusion within the layer of molecules adsorbed on the pore walls. The relative contributions of these modes, which are determined by the details of the interphase equilibrium, change with variation of the pore loading, leading to a complex behavior of the effective self-diffusion coefficient. The impact of the pore size and the adsorbate-surface interaction on self-diffusion is elucidated. Possible reasons for an experimentally obtained hysteresis in the diffusivities measured on adsorption and desorption in mesopores are discussed.     

129Xe NMR study of Xe adsorption on multiwall carbon nanotubes

¹²⁹Xe NMR spectroscopy has been used to study the adsorption of Xe on multi-wall carbon nanotubes (MWCNT). The results obtained have shown the ¹²⁹Xe NMR ability to probe the intercrystalline (aggregate) and the inner porosity of CNT. In particular, the effects on porosity of tubes openings by hydrogen exposure and of ball milling were examined. Dramatic changes observed in the ¹²⁹Xe NMR spectra after moderate ball milling of MWCNTs were attributed to the destruction of the initial intercrystalline pore structure and to the Xe access inside the nanotubes. To examine the exchange dynamics the mixture of as-made and milled MWCNTs was studied with one- and two-dimensional (1D and 2D) ¹²⁹Xe NMR. The exchange between the interior of milled nanotubes and the aggregate pores of as-made MWCNTs was fast on the NMR acquisition time scale. The Xenon exchange between the interior of the as-made MWCNTs and the large aggregate pores occurred on a longer time scale of 10 ms, as was established by 2D ¹²⁹Xe NMR exchange spectroscopy. Variable temperature ¹²⁹Xe NMR data were also discussed and analyzed in terms of the fast exchange approximation.     

注气驱替煤层气数值模拟

建立一种注气驱替煤层气的双重介质数学模型,考虑注气驱替煤层气中多组分气体渗流、吸附/解吸、扩散及孔隙度和渗透率敏感性等多物理场的耦合,用有限差分方法对数学模型数值求解,并验证模型的有效性.通过实例模拟,对比分析注入CO₂、N₂和70%N₂+30%CO₂提高采收率、产出组分动态及储层物性.计算表明:注入CO₂、N₂及其混合气体均能显著提高煤层气藏采收率,提高程度可达28%以上;N₂在注入后迅速突破产出,而CO₂突破时间较晚;注气后基质的收缩/膨胀效应在孔隙度和渗透率敏感性中起主导作用,CO₂吸附造成孔隙度和渗透率减小,而N₂吸附将增大孔隙度和渗透率.     

NMR study of the vapor phase contribution to diffusion in partially filled silica glasses with nanometer and micrometer pores

The contribution of the vapor phase to molecular diffusion in porous silica glasses with nanometer (Vycor) and micrometer (VitraPor#5) pores partially filled with water (polar) or cyclohexane (nonpolar) was investigated with the aid of field-gradient NMR diffusometry. Due to the vapor phase, the effective diffusion coefficient of cyclohexane filling micrometer pores (VitraPor#5) increased up to 10 times relative to the value in bulk liquid upon reduction of the pore space filling factor. On the other hand, the effective diffusion coefficient of water first decreases and then increases when the liquid content is reduced. The dependence of the effective diffusion coefficient on the pore filling factor is strongly related to the pore dimension. A general two-phase exchange model is presented that is well accounting for all experimental diffusion features.     

Improved Halbach sensor for NMR scanning of drill cores

A lightweight Halbach magnet system for use in nuclear magnetic resonance (NMR) studies on drill cores was designed and built. It features an improved homogeneous magnetic field with a strength of 0.22 T and a maximum accessible sensitive volume. Additionally, it is furnished with a sliding table for automatic scans of cylindrical samples. This device is optimized for nondestructive online measurements of porosity and pore size distributions of water-saturated full cylindrical and split semicylindrical drill cores of different diameters. The porosity of core plugs with diameters from 20 to 80 mm can be measured routinely using exchangeable radiofrequency coils. Advanced NMR techniques that provide 2D T(1)-T(2) correlations with an average measurement time of 30 min and permeability estimates can be performed with a special insert suitable for small core plugs with diameter and length of 20 mm.