《T1与T2的比值在孔隙性评估中的作用要点》由会员分享,可在线阅读,更多相关《T1与T2的比值在孔隙性评估中的作用要点(18页珍藏版)》请在金锄头文库上搜索。
1、T1/T2比在孔隙性评估中的作用选自2008年第49届美国测井分析家年会论文集摘要核磁共振测井数据往往通过一个伴随一系列短脉冲的长CPMG序列获得。这 些短脉冲有一个顺序的20ms的等待时间,它们主要用于提高短T2组分的孔隙度 精度。然后再使用联合反演的方法以得到准确的核磁共振总孔隙度。这个联合反 演的技术还使我们能够优化T1/T2的比值。由于常规CPMG回波序列的等待时间 (1-10s)和它的短脉冲回波序列的等待时间差异非常大,并且数值模拟也表明 T1/T2的比值强烈的依赖于噪声概况,因此它也可以对核磁共振总孔隙度产生巨 大的影响,尤其是对像碳酸盐岩这样的长T2驰豫时间的地层。我们还发现对于
2、 那些有大量顺磁杂质的地层,优化T1/T2比的范围应当从常用的1至3提高到1 至20。为了进一步研究这个问题,我们使用了 T1-T2二维核磁共振技术对岩心 进行了核磁共振测量并且发现T1/T2的比值强烈的依赖于T2驰豫时间。我们已 经将这个依赖于T1/T2比的功能加入到反演代码当中,并且大大提高了在碳酸盐 岩地层和含有大量顺磁杂质地层中确定核磁共振总孔隙度的精度。引言核磁共振测井通常用于测量孔隙度、束缚水饱和度、渗透率。核磁共振测井 得到的 CPMG (Carr and Purcell 1954; Meiboom and Gill 1958)回波序列数据 被拉普拉斯变换倒置后,就可以得到岩石物
3、理解释所需要的核磁共振T2分布。 T2分布的质量和准确性依赖于两个主要的数据采集参数:回波间隔(TE)和两 个相邻CPMG回波序列之间的等待时间(WT)。最低回波间隔主要由所使用的核磁 共振测井仪器的频率所决定。所有商用核磁共振测井仪器的最低回波间隔都是 0.2ms或者更长。当顺磁杂质存在时,孔隙流体和颗粒表面磁化率的差异将会伴 随外部磁场的应用在孔隙空间引发静磁场梯度。这种梯度强度的变化范围取决于 顺磁杂质的数量及其分布。随着外部磁场梯度,这种诱发磁场梯度增强扩散从而 影响到T2的分布。极化等待时间的长短受测井速度、磁铁长度和NMR仪器的操作拉莫尔频率值 所限制。对于一个典型的测井速度,等待
4、时间也就几秒钟。短的等待时间可能致 使T1强烈依赖于极化因素1-exp(-WT/T1)。因此,T1分布在获得准确的T2分 布和核磁总孔隙度中起很重要的作用。因为测量T1分布是很耗费时间的,所以 我们常引入T1/T2比值来将对T1的依赖变换为对T2的依赖。这种T1/T2比可以 从核磁共振分析中确定。克莱因伯格等人用2MHz的核磁共振光谱仪研究了 105 块岩石样品,发现在回波间隔TE=0.2ms时,T1/T2比在1到3之间变化并且其 平均值为1.65。假设一种流体的T2值为1秒且等待时间WT也为1秒,利用公 式1-exp(-WT/T1)计算极化因素,当T1/T2=1时,极化因素为0.63,当T1
5、/T2=3 时,极化因素为0.28。因此,如果从NMR T2测井获得的没有经过极化校正的表 观孔隙度为20pu,那么,当T1/T2=1时,极化校正后的孔隙度可能为31.6pu, 当T1/T2=3时,极化校正后的孔隙度可能是70.5pu。很明显,T1/T2比值在孔隙 度评价上的作用是非常重要的。为了利于T2反演而优化T1/T2比,现代NMR数据常利用一个长CPMG序列和 跟随它的一系列短脉冲获得(例如Enhanced Precision Mode in CMR(Schlumberger) or PRIME of MRIL (Halliburton) with short bursts (Chen
6、 et al. 2000)。这些短脉冲有20ms的顺序等待时间,他们主要用于提高短T2 组分的孔隙度精度。这些数据集通常用一种联合反演方法处理,该方法可以在1 到3的范围内优化T1/T2比从而获得一个准确的NMR总孔隙度。这个反演方法可 以在大多数纯砂岩地层中获得很好的效果。然而,当粘土的T2组分变的过短而 不能被检测到,或者地层中含有大量的顺磁杂质原子,这种做法不能准确的估计 孔隙度。此外,对于含有长T2驰豫时间的轻烃液体的碳酸盐岩,这种方法也可 能确定出不准确的孔隙度。更确切的说,数值模拟表明,当T2过长时,联合反 演T1/T2比强烈的取决于噪音概况。为了在联合反演中改善T1/T2比的测定
7、,我们进行了 NMR岩心测量来显示二 维T1-T2分布,以便每个T2分量的T1特性都可以得到明确的界定。我们的目标 是确定怎样处理和解释受到强烈的内部场梯度或在例如碳酸盐岩这样有长T2分 量存在的影响下获得的NMR测井数据。T1/T2比对孔隙度的影响我们进行数值模拟来说明T1/T2比对NMR总孔隙度的影响。合成数据串包含 一个主要的CPMG回波序列和一个短的回波脉冲。主要的回波序列有600个回波, 每个回波有一个TE=0.2ms的回波间隔和一个WT=3000ms的等待时间。短脉冲有 30个回波,每个回波的TE=0.2ms,WT=20ms。原始回波数据是由一个T2分布模 型产生的。这个T2分布模
8、型带有一个可以从短驰豫时间变到长驰豫时间的单独 的T2值、一个基于上面的固定T1/T2=r=1的采集参数并且可以叠加高斯噪声。 该模型的孔隙度假设为100pu。因此第i个回波幅度可由以下得出:&二(1-仑项以)疽& +(1)然后,我们在联合SVD反演算法的基础上,通过T2反演出原始回波数据。4=力必1-疗心*” +与.村,=+与)=1公式中匕是反演矩阵,r是考虑到主回波序列和短脉冲的极化因素而用于反演 矩阵的T1/T2比。在此研究中用的r= T1/T2是1,3,5或7。图一显示了一个 T1/T2比的函数的拟合误差结果,该函数用于带有变化范围从1ms到2000ms的 独立T2值的T2分布模型,其
9、中分别假定T1/T2比是1(A)或者5(B)。我们注 意到当T2为10ms,100ms,甚至1000ms时,选择二次方程式最小值的优化方案 将恰好是预测的T1/T2=1或者5,这是正确的。然而,当T2是1ms或2s时,拟 合误差具有像高斯噪声一样的规律并且还稍微有些随机性。因此,选择二次方程 式最小值的优化方案不一定会产生正确的T1/T2比。对于1ms这样的短T2值, 一个不正确的T1/T2比是无关紧要的,因为它对极化因子产生很小的影响。然而, 对于2s这样的长T2值,一个不正确的T1/T2比将导致在孔隙度计算上的巨大误 差。接着,我们的讨论将集中在对于长T2组分的核磁共振反演处理上。图二显示
10、了当单一 T2驰豫时间是2s并且噪声水平是1pu(A)和4pu(B) 时,通过公式(1)计算的合成原始回波数据。图3A和3B显示了当噪声水平分 别是1pu(A)和4pu(B)时,来自于联合反演的T2分布。T1/T2比为1,3,5 和7的拟合曲线在图二中分别用深绿色、蓝绿色、紫红色和红色的线表示,而这 些回波序列(无噪声)的拟合曲线是通过对T2分布的计算得到的。这四条曲线 彼此间基本上难以区分(这就是为什么我们在图二中只看到一条红线)。不同 T1/T2比的拟合误差之间的差异比噪声水平本身要小十倍多。然而反演所给出的 视孔隙度的误差可能超过300% (例如在1pu的噪声下,T1/T2=1时,孔隙度
11、为 1.0;但是T1/T2=7时,孔隙度却为4.2)。这个数字例子清楚的表明了 T1/T2比 对于核磁共振总孔隙度有巨大的影响。为了阐明问题的统计特性,我们以相同的采集参数和模型(T2=2s,T1/T2=1) 为基础,生成10000列带有随即噪声的回波序列。图四显示了在噪声是1pu,T1/T2 比分别为1,3, 5和7时的拟合误差的频率图。由于用于产生原始回波数据的白 噪声的性质,所以各种情况下的拟合误差都呈一种正态分布。每种分布的平均值 都接近于原始数据的噪声水平,但也会随着T1/T2比而增长。正态分布的宽度大 概是噪声水平的10%并且不随着T1/T2比的不同而改变。这个特性表明当T1/T2
12、=1 时比高T1/T2值时更容易找到更多的解决方法。为了得到T1/T2比和有效解决方 法的数量之间的关系,在假定解决方法有效的情况下,我们可以使用拟合误差截 止法。假定误差截止值就像图四中深蓝色线标示的那样比噪声水平高10%或者是 1.1pu,那么当T1/T2比分别为1,3, 5和7的时候,有效解决方法的总数量将 分别占解决方法总数量的93.8%,57.7%,25.6%和10.1%。为了进一步分析有效方法的数量与T1/T2比及噪声水平之间的关系,我们在 不同的噪声水平下(1pu和4pu)对很多组T1/T2比进行了相同的模拟。结果如 图五所示。在低噪声水平的情况下,有效方法的数量随T1/T2比的
13、增加而快速减 少,但是在高噪声水平的情况下,这种减少却相当缓慢。同时,平均视孔隙度几 乎随T1/T2比呈线性增长,并且这种关系不会随噪声水平的不同而发生很大的改 变。但是孔隙度分布的宽度(在这里不会被显示)的确会随着噪声水平或T1/T2 比的增加而增加。我们同样注意到了拟合误差分布的中心会随着T1/T2比的增加 而移离噪声水平处,尤其是在低噪声水平的情况下(见图五C)。同样的,视孔 隙度的平均值也会随着孔隙度分布中的T1/T2比而发生移动(见图五D)。我们已经表明的一个最糟糕的例子就是视NMR总孔隙度和T1/T2比有强烈的 依赖关系。随着模型的T2驰豫时间的减少,情况会变的更好一点。然而像2s
14、 这样的长T2驰豫时间并不是不可能出现的情况。在许多油田,T2驰豫时间的高 峰值都会达到2s左右。T1/T2比的一丝轻微的不确定都可能造成孔隙度上的重 大的不确定。此外,随着T2驰豫时间的增加和表面驰豫机制的效果的降低,T1/T2 比将会减少到1。以后,我们应当在基于岩心测量的基础上,在T2反演中使用 一个和T1/T2比相关的函数。这种新方法已经被证实了可以正确的预测孔隙度, 并且可以避免由于在长T2驰豫时间中不正确的T1/T2比所引起的巨大的孔隙度 波动。高T1/T2比的现场实例数值模拟已经表明对于碳酸盐岩这样的长T2驰豫时间地层,不准确的T1/T2 比会对孔隙度评价产生巨大的影响,现在我们
15、看一下另外一个还有大量顺磁物质 地层的例子,对于这样的地层,当前的在1到3的范围内优化T1/T2比的方法是 更本不够的。图六显示的NMR孔隙度来自于一个取得EPM (Enhanced Precision Mode) 模式的CMR测井仪。T1/T2比通过联合反演算法在1到3的范围内进行了优化。 地层中包含了大量的可以在孔隙内引起巨大内部磁场梯度的顺磁物质。这致使 T2驰豫率得到很大的加强,并且直接影响到了用于T2反演中的T1/T2比。图六 中的第三道显示了核磁共振得到的孔隙度与中子-密度交会图确定的总孔隙度之 间的对比。在第三道中NMR所测孔隙度与交会图确定的孔隙度相比显著减少的部 分在图中用橙
16、色阴影标示了出来。第四道显示了通过CMR仪的T2反演所得到的 最佳T1/T2比(蓝色实线)和通过实验室岩心测量得到的T1/T2比(红色小点)。 后者要求使用T1、T2各自的对数平均数。第四道的两种测井方法间的差异表明 了用于常规反演中的T1/T2比的范围是肯定不恰当的。实验室测量法为了完全的研究问题,我们对这种地层的17个岩心样品进行了实验室核磁 共振T1、T2值的测量。我们也对完全饱含盐水的岩心样品使用2MHz的Maran Ultra共振仪进行二维T1-T2图的测量。对于每个岩心塞,设定三个不同的回波 间隔:0.21ms、0.6ms、1.2ms,对于每个回波间隔进行10个不同的等待时间(0.3, 1,3,10,30,100,300,1000,3000 和 10000ms)的测量,这样就会得到 10 个回波序列。每个回波间隔的数据产生一个相应的二维T1-T2图。回波间隔为 0.21ms的二维T1-T2图被认为是最准确的,因为它丢失的信号最少。它很接近 于