水驱油过程的核磁共振二维谱研究

岩心水驱油过程中油水分布状况是岩心多孔介质的重要性质. 水驱油过程的研究是进一步进行提高采收率研究的基础. 核磁共振扩散-弛豫二维谱提供了岩心中流体性质的多方面信息,与核磁共振一维弛豫谱相比极大地提高了区分油水的能力. 该文通过2组岩心水驱油实验,从不同含油饱和度的扩散-弛豫二维谱中提取出水的一维弛豫谱,在原油粘度比较高的情况下获得了驱替过程中油水在不同孔隙中的分布状况以及润湿性等信息, 解决了单独用一维弛豫谱方法难以区分油水的问题. 该文的研究方法对油田提高采收率的研究有比较大的参考价值.     

岩心孔隙介质中流体的核磁共振弛豫

弛豫时间是核磁共振研究中的一个重要参数,岩心孔隙介质流体的弛豫过程是自由流体弛豫机制、表面弛豫机制和流体的扩散弛豫机制共同作用的结果,它包含了丰富的孔隙和流体本身的信息. 弛豫时间和自扩散系数的测量及对弛豫时间的分析是核磁共振技术应用于岩心分析和石油勘测的重要内容.     

低场核磁共振扩散-横向弛豫二维反演算法

针对低场核磁共振一维反演中无法分辨一维谱中重叠组分和目前报道的扩散-横向弛豫二维反演算法计算量大、计算耗时长的问题,提出了一种计算量小、计算效率高、耗时短的扩散-横向弛豫二维反演算法. 首先对扩散系数D-和横向弛豫时间T₂进行布点;其次根据信号采集条件计算出两个核心矩阵,并分别进行奇异值分解;然后,由所采集信号计算出两个核心矩阵的奇异值截断值,分别对两个核心矩阵的奇异值矩阵进行截断并求其逆矩阵;最后计算出初始反演结果,并添加非负约束经过多次迭代得到最终反演结果. 实验结果证明,提出的扩散-横向弛豫二维反演算法在不影响反演结果准确性的基础上,能极大提高计算效率.     

核磁共振二维谱反演

核磁共振二维谱包含扩散系数 D 和横向弛豫时间 T₂ 的信息,利用核磁共振仪来测量多孔介质中物质的信息,根据其弛豫时间和扩散系数的差别来区分不同物质;利用全局反演方法,提出了核磁共振二维谱反演的物理模型和数学模型;介绍了传统的奇异值分解（SVD）和改进的奇异值分解反演算法;采用改进的奇异值分解法对核磁共振二维谱进行反演,其反演算法具有计算速度快和二维谱分布连续等优点. 它适合于信噪比较高的数据反演,当原始数据信噪比SNR≥100时,可以对二维谱图进行定量分析;当60≤SNR<100时,可以对二维谱图进行定性分析. 核磁共振二维谱可以一次性直接区分油和水,为核磁共振测井提供了新的科学技术.     

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

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

一种基于组分补偿的二维核磁共振测井数据高精度处理方法

二维核磁共振技术能够对储层中各类含氢流体进行无损、快速、定量的测量和表征，但受限于采集方式和参数，核磁共振设备在对页岩油等致密储层中的有机质、沥青等超快弛豫组分进行检测时，经常出现由于信号采集不完整所导致的二维谱中流体组分缺失或不准的问题.本文提出了基于超快弛豫组分补偿技术的T₂-T₁二维谱高精度反演方法，该方法将一维核磁共振前端信号补偿技术进行推广，通过在二维核磁数据反演前对回波数据进行组分补偿，能够有效解决二维核磁共振测井前端信号漏失的问题.实验及测井数据的应用表明，该方法在页岩油等富含快弛豫组分信号的储层中，可以得到更加精准和完整的储层信息.     

水分子在微孔隙介质中的受限扩散模拟

在饱和多相流体孔隙介质中,利用核磁共振测得的扩散系数可以区分油气水,并用于饱和度等岩石物理参数的求取,但在微小孔隙(20μm)中,流体的表观扩散系数会偏离其自由扩散系数.通过有限差分和随机游走模拟方法对水分子在单孔隙及多孔隙介质中的扩散进行数值模拟,考察水分子扩散受限的程度.结合核磁共振在测井及实验室岩石物理应用的技术,进一步分析了受限扩散对求取油水饱和度等岩石物理参数造成的影响.同时,分析总结了两种模拟方法的使用条件及优缺点.     

水分子在微孔隙介质中的受限扩散模拟

在饱和多相流体孔隙介质中,利用核磁共振测得的扩散系数可以区分油气水,并用于饱和度等岩石物理参数的求取,但在微小孔隙（< 20 μm）中,流体的表观扩散系数会偏离其自由扩散系数．通过有限差分和随机游走模拟方法对水分子在单孔隙及多孔隙介质中的扩散进行数值模拟,考察水分子扩散受限的程度．结合核磁共振在测井及实验室岩石物理应用的技术,进一步分析了受限扩散对求取油水饱和度等岩石物理参数造成的影响．同时,分析总结了两种模拟方法的使用条件及优缺点．     

基于二维核磁共振弛豫谱的雷四段孔隙度计算方法

基于常规测井资料的体积模型和一维核磁共振（NMR）测井资料的固定截止值法,可以计算地层孔隙度,但是扩径段孔隙度计算结果偏大.本文通过川西海相雷四段岩心样品NMR实验和实测二维NMR测井资料,开展扩径段的T₂、T₁孔隙度计算方法研究.首先总结不同孔径流体在T₂谱和T₁谱上的响应特征,分析钻井液流体峰截止值的分布范围和影响因素,建立钻井液流体弛豫时间截止值的计算模型;然后根据未扩径段和扩径段的粘土束缚水弛豫时间截止值确定原则,确立变粘土束缚水弛豫时间截止值的有效孔隙度计算方法.多口井的应用效果表明,基于T₁谱的孔隙度计算结果精度更高、定量分析误差小,能有效解决扩径段孔隙度计算结果偏大的问题,满足储层评价的要求.     

川西气田白云岩储层二维核磁共振测井气水识别方法

川西气田雷口坡组四段白云岩储层具有良好的天然气勘探开发前景,但是储层分布具有较强的非均质性,利用常规测井和一维核磁共振测井方法判断储层流体性质存在多解性．通过分析川西雷四段白云岩储层地质特征与测井响应特征,开展岩心二维核磁共振实验,明确了实验室分析与测井过程中环境、参数等影响因素,建立了岩心分析结果的校正方法,明确了钻井液、束缚流体、可动水与天然气信号在二维核磁共振测井中谱分布区间,建立了基于T₂-T₁、T₁/T₂（R）的白云岩储层二维核磁共振测井气水识别图版．利用图版对川西雷四段储层流体性质开展评价,二维核磁共振测井解释结论得到了实钻测试结果的验证,该气水识别方法填补了常规测井、一维核磁共振测井在评价储层气水关系中的缺陷,可有效解决白云岩储层流体性质判别难题．     

A novel improved method for analysis of 2D diffusion-relaxation data--2D PARAFAC-Laplace decomposition

This paper demonstrates how the multi-linear PARAFAC model can with advantage be used to decompose 2D diffusion-relaxation correlation NMR spectra prior to 2D-Laplace inversion to the T(2)-D domain. The decomposition is advantageous for better interpretation of the complex correlation maps as well as for the quantification of extracted T(2)-D components. To demonstrate the new method seventeen mixtures of wheat flour, starch, gluten, oil and water were prepared and measured with a 300 MHz nuclear magnetic resonance (NMR) spectrometer using a pulsed gradient stimulated echo (PGSTE) pulse sequence followed by a Carr-Purcell-Meiboom-Gill (CPMG) pulse echo train. By varying the gradient strength, 2D diffusion-relaxation data were recorded for each sample. From these double exponentially decaying relaxation data the PARAFAC algorithm extracted two unique diffusion-relaxation components, explaining 99.8% of the variation in the data set. These two components were subsequently transformed to the T(2)-D domain using 2D-inverse Laplace transformation and quantitatively assigned to the oil and water components of the samples. The oil component was one distinct distribution with peak intensity at D=3 x 10(-12) m(2) s(-1) and T(2)=180 ms. The water component consisted of two broad populations of water molecules with diffusion coefficients and relaxation times centered around correlation pairs: D=10(-9) m(2) s(-1), T(2)=10 ms and D=3 x 10(-13) m(2) s(-1), T(2)=13 ms. Small spurious peaks observed in the inverse Laplace transformation of original complex data were effectively filtered by the PARAFAC decomposition and thus considered artefacts from the complex Laplace transformation. The oil-to-water ratio determined by PARAFAC followed by 2D-Laplace inversion was perfectly correlated with known oil-to-water ratio of the samples. The new method of using PARAFAC prior to the 2D-Laplace inversion proved to have superior potential in analysis of diffusion-relaxation spectra, as it improves not only the interpretation, but also the quantification.     

Quantitative evaluation of porous media wettability using NMR relaxometry

We propose a new method to determine wettability indices from NMR relaxometry. The new method uses the sensitivity of low field NMR relaxometry to the fluid distribution in oil-water saturated porous media. The model is based on the existence of a surface relaxivity for both oil and water, allowing the determination of the amount of surface wetted either by oil or by water. The proposed NMR wettability index requires the measurement of relaxation time distribution at four different saturation states. At the irreducible water saturation, we determine the dominant relaxation time of oil in the presence of a small amount of water, and at the oil residual saturation, we determine the dominant relaxation time of water in the presence of a small amount of oil. At 100% water and 100% oil saturation, we determine the surface relaxivity ratio. The interaction of oil with the surface is also evidenced by the comparison of the spin-lattice (T1) and spin-locking (T1rho) relaxation times. The new NMR index agrees with standard wettability measurements based on drainage-imbibition capillary pressure curves (USBM test) in the range [-0.3-1].     

Diffusion-relaxation distribution functions of miscible fluids measured in grossly inhomogeneous fields

In this paper, the recently developed technique of measuring diffusion-relaxation distribution functions in grossly inhomogeneous fields was applied to the characterization of bulk fluids. The measurements were conducted in a static stray field of a superconducting magnet at a Larmor frequency of 5 MHz, relevant for ex situ nuclear magnetic resonance applications in oil-well logging or materials testing. Results are presented for samples of water and hydrocarbon oils, including a number of crude oils. The measurements indicate strong correlations between diffusion coefficients and relaxation times that are sample specific. The diffusion-relaxation correlation function provides information on the correlation between the rotational and the translational diffusion coefficients of each component of the fluid.     

应用核磁共振进行聚驱后泡沫驱油特性研究

采用并联岩心进行聚驱后泡沫驱油实验,利用核磁共振技术,对驱替后岩心的不同直径孔隙内的流体分布进行了研究,得到了水驱、聚合物驱、泡沫驱替阶段的驱出油的孔径范围以及剩余油分布.实验结果表明,与水驱和聚合物驱相比,泡沫驱增大了波及的孔径的范围.泡沫可以在水驱和聚驱易发生窜流的大孔径通道形成封堵,从而波及到了水驱和聚驱未波及到的孔径,不但大幅度的提高了低渗岩心的采收率,也驱出了部分高渗岩心的小孔径的油.     

Exploring the separate NMR responses from crude oil and water in rock cores

In analysis of transverse relaxation time (T ₂) curves in a Carr-Purcell-Meiboom-Gill (CPMG) experiment in a multicomponent system originating from measurements of oil and water in rock cores, where internal magnetic field gradients broaden the line widths significantly, there is very little direct information to be extracted of the different components contributing to the totalT ₂ relaxation time curve. From the study of rock cores saturated with different amounts of crude oil and water, we show that with an optimised experimental setup it is possible to extract information from the nuclear magnetic resonance response that is not resolved by any other methods. This setup combines pulsed field gradient methods with the CPMG experiment utilizing data from both rock cores and bulk oil and water. Then it becomes feasible to separate the signals from oil and water where the two-dimensional inverse Laplace transform ordinarily seems to fail.     

Diffusion-relaxation distribution functions of sedimentary rocks in different saturation states

We present diffusion-relaxation distribution functions measured on four rock cores that were prepared in a succession of different saturation states of brine and crude oil. The measurements were performed in a static gradient field at a Larmor frequency of 1.76 MHz. The diffusion-relaxation distribution functions clearly separate the contributions from the two fluid phases. The results can be used to identify the wetting and non-wetting phase, to infer fluid properties of the phases, and to obtain additional information on the geometrical arrangement of the phases. We also observe effects due to restricted diffusion and susceptibility induced internal gradients.     

Rapid encoding of T(1) with spectral resolution in n-dimensional relaxation correlations

Nuclear magnetic resonance (NMR) T(1) relaxation times have been encoded in the second dimension of two-dimensional relaxation correlation and exchange experiments using a rapid "double-shot"T(1) pulse sequence. This technique also retains chemical shift information (delta) for short T(2)( *) materials. In this way, a spectral dimension can be incorporated into a T(2)-T(1)-delta correlation without an increase in experimental time compared to the conventional, chemically insensitive T(1)-T(2) correlation. Here, the T(2)-T(1)-delta pulse sequence is used to unambiguously identify oil and water fractions in a permeable rock. A novel T(1)-T(1)-(delta) relaxation exchange measurement is also introduced and used to observe diffusive exchange of water in cellulose fibres.     

分子动力学模拟研究方解石表面润湿性反转机理

利用分子动力学模拟技术从分子尺度探究方解石表面润湿性反转机理.首先,研究方解石表面润湿性反转过程;而后,从原油分子-方解石表面与原油分子-原油分子/水分子相互作用两个方面系统揭示方解石表面润湿性反转机理.结果:(1)水分子能够驱离方解石表面弱吸附的非极性分子造成润湿性的改变,但不能驱离强吸附的极性分子使润湿性反转难以实现;(2)原油分子极性越强与方解石表面相互作用越强,极性分子与方解石表面之间主要为静电力,非极性分子与方解石表面之间主要为范德华力;(3)原油分子极性越相近分子之间的相互作用越强,分子极性相差越大分子之间的相互作用越弱.非极性分子之间主要是范德华力,极性分子之间主要是静电力;(4)原油分子在方解石表面和水分子的共同作用下形成乙酸-吡啶-水-甲苯-己烷的稳定吸附序列.本研究为靶向提高采收率技术的设计与应用提供理论基础.