一种岩心核磁共振T2谱转化孔径分布方法及系统技术方案

本申请涉及一种岩心核磁共振T2谱转化孔径分布方法及系统。该方法包括：对岩心样品孔径分布特征进行CT扫描，获得孔径累积分布曲线；对CT扫描后的岩心样品进行原油的抽真空饱和，待岩心充分饱和完毕后将岩心取出；将取出后的饱和充分的岩心样品进行核磁共振扫描，获得T2谱分布曲线，进而获得不同初始值Ts下的T2谱累积分布曲线；根据T2谱累积分布曲线，获得T2谱累积分布图版；校对并匹配孔径累积分布与T2谱图版。通过改变T2谱累积分布曲线的初始值Ts以匹配CT扫描结果的孔径累积分布曲线，操作更加方便快捷，并且拟合匹配效果较好，转化精度较高。度较高。度较高。

【技术实现步骤摘要】
一种岩心核磁共振T2谱转化孔径分布方法及系统


[0001]本申请涉及岩石的孔隙结构领域，具体地涉及一种基于CT扫描技术的岩心核磁共振T2谱转化孔径分布方法及系统。

技术介绍

[0002]随着能源需求的日益增长，致密油等非常规能源引起了更加广泛的关注和重视，逐渐在我国能源供给方面发挥着越来越重要的作用。不同于常规储层，致密等非常规储层的孔隙结构较为复杂，主要体现在孔隙连通性差和孔径分布差异大等方面。孔隙结构影响着致密储层油气的赋存、运移及储层质量的优劣，同时对产能具有重要影响。因此为了更好的开发致密等非常规储层，准确的描述岩石的孔隙结构显得尤为关键。
[0003]核磁共振技术由于其快速无损的特点，在研究岩土材料工程特性和微观结构特征的关系领域发挥着越来越重要的作用。现有研究核磁共振T2谱弛豫时间转化孔径分布的方式是结合压汞测试获得的孔径分布特征和T2谱分布特征，利用最小二乘法等相关算法使得两者误差最小，从而获得孔径大小和T2谱弛豫时间两者间的转化系数，进而将核磁共振的T2谱弛豫时间转化为孔隙半径。但是针对致密等孔隙连通性较差的岩心，高压压汞测试不能够识别因为连通性较差使得汞不能进入的孔隙，同时测试样品较小缺乏一定的代表性，因此会对孔径分布结果造成较大的影响，进而会对孔径转化造成较大的误差。

技术实现思路

[0004]针对上述问题，本申请的目的是提供一种基于CT扫描技术的岩心核磁共振T2谱转化孔径分布方法，该方法可以较好的实现T2谱弛豫时间向孔隙半径的转化，同时对致密等连通性较差的岩心孔喉结构识别能力更强，测试样品具有代表性，拟合效果更好，大大提高了孔径分布的转化效果，为后续岩石孔径分布的研究提供了有效途径。
[0005]为实现上述目的，本申请采取以下技术方案：
[0006]第一方面，本专利技术提供的一种岩心核磁共振T2谱转化孔径分布方法，包括：
[0007]对岩心样品孔径特征进行CT扫描获得孔径累积分布曲线；
[0008]对CT扫描后的岩心样品进行原油的抽真空饱和，待岩心充分饱和完毕后将岩心取出；
[0009]将取出后的饱和充分的岩心样品进行核磁共振扫描获得T2谱累积分布曲线；
[0010]根据T2谱累积分布曲线，获得不同初始值下的T2谱累积分布曲线；
[0011]根据T2谱累积分布曲线，获得T2谱累积分布图版；
[0012]校对并匹配孔径累积分布与T2谱累积分布图版；
[0013]根据匹配好的孔径累积分布曲线和T2谱累积分布曲线建立T2谱和孔径大小的关系。
[0014]对岩心样品孔径特征进行CT扫描获得孔径累积分布曲线包括：对干燥洁净的柱塞状岩心样品进行微纳米CT扫描，对柱塞状岩心样品沿程的切片扫描结果进行统计分析以获
得不同孔径的占比，最后获得孔径累积分布曲线。
[0015]将取出后的饱和充分的岩心样品进行核磁共振扫描获得T2谱累积分布曲线包括：获得岩心样品的弛豫时间T2和信号幅度的关系，弛豫时间T2的数值在一定程度上反应孔径的大小，信号幅度大小反应了对应弛豫时间T2值的孔径中所含油样的信号强度，信号幅度与弛豫时间T2围成的面积大小可反映所含油样的多少，以便获得不同弛豫时间T2数值范围内的峰面积占比，进而获得T2谱累积分布曲线。
[0016]获得T2谱累积分布图版包括：将获得了T2谱从0ms开始的累积分布曲线，分别将T2谱的起始值Ts设置为5ms、10ms、15ms、20ms，然后分别作出不同弛豫时间T2起始值下的累积分布曲线，绘制在同一坐标系下，得到不同弛豫时间T2起始值下的T2谱累积分布图版。
[0017]校对孔径累积分布与T2谱图版包括：
[0018]根据关系式(1)获得横向弛豫时间T2和孔隙半径r存在一定的转化关系：
[0019][0020]其中
[0021]V为孔隙体积；
[0022]S为孔隙表面积；
[0023]ρ
T
为弛豫率；
[0024]F
s
为孔隙形状因子；
[0025]r为孔隙半径。
[0026]将CT扫描结果获得的孔径累积分布曲线绘制在T2谱累积分布图版中，通过左右平移孔径累积分布曲线来匹配不同起始值的T2谱累积分布曲线，最终获得匹配好的结果。
[0027]孔径累积分布曲线满足关系式(2)：
[0028]Q1＝f1(r)
ꢀꢀꢀ
(2)
[0029]T2谱累积分布曲线满足关系式(3)：
[0030]Q2＝f2(T2)
ꢀꢀꢀ
(3)
[0031]其中
[0032]Q1为孔径累积分布变化量；
[0033]Q2为T2谱累积分布变化量；
[0034]f1为拟合的孔径累积变化量与孔隙半径的函数关系；
[0035]f2为拟合的T2谱信号幅度累计变化量与横向弛豫时间的函数关系；
[0036]根据Q1＝Q2可得f1(r)＝f2(T2)，
[0037]从而建立T2谱中不同横向弛豫时间与孔隙半径的关系。
[0038]第二方面，本专利技术提供的一种岩心核磁共振T2谱转化孔径分布系统，包括：
[0039]第三方面，本专利技术提供的一种计算机存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现本专利技术第一方面所述的岩心核磁共振T2谱转化孔径分布方法的步骤。
[0040]第四方面，本专利技术提供的一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现本专利技术第一方面所述的岩心核磁共振T2谱转化孔径分布方法的步骤。
[0041]本申请由于采取以上技术方案，其具有以下优点：
[0042]本申请提供了核磁共振中弛豫时间T2向对应岩心孔隙半径转化的方法，不同于传统通过压汞测试获取孔径转化系数的方法，CT扫描技术在识别复杂孔隙结构更直观、更准确且样品更具有代表性，特别针对存在的微小孔隙及盲端孔隙，从而能够获得更加准确的孔径累积分布特征，保证了孔径分布转化的准确性，为致密页岩等非常规储层原油在岩心中分布特征的定量研究奠定了基础；
[0043]通过改变T2谱累积分布曲线的初始值Ts以匹配CT扫描结果的孔径累积分布曲线，操作更加方便快捷，并且拟合匹配效果较好，转化精度较高；
[0044]随着纳米CT扫描技术的发展，能够将T2谱转化成孔隙半径的下限持续增大。
附图说明
[0045]通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本申请的限制。在整个附图中，用相同的附图标记表示相同的部件。在附图中：
[0046]图1为岩心样品中存在的不同孔隙结构类型示意图；
[0047]图2为充分饱和后的岩心样品的核磁共振T2谱结果；
[0048]图3为根据核磁共振T2谱获得的不同起始本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种岩心核磁共振T2谱转化孔径分布方法，其特征在于，包括：对岩心样品孔径特征进行CT扫描获得孔径累积分布曲线；对CT扫描后的岩心样品进行原油的抽真空饱和，待岩心充分饱和完毕后将岩心取出；将取出后的饱和充分的岩心样品进行核磁共振扫描获得T2谱累积分布曲线；根据T2谱累积分布曲线，获得不同初始值下的T2谱累积分布曲线；根据T2谱累积分布曲线，获得T2谱累积分布图版；校对并匹配孔径累积分布与T2谱累积分布图版；根据匹配好的孔径累积分布曲线和T2谱累积分布曲线建立T2谱和孔径大小的关系。2.根据权利要求1所述的岩心核磁共振T2谱转化孔径分布方法，其特征在于，所述对岩心样品孔径特征进行CT扫描获得孔径累积分布曲线包括：对干燥洁净的柱塞状岩心样品进行微纳米CT扫描，对柱塞状岩心样品沿程的切片扫描结果进行统计分析以获得不同孔径的占比，最后获得孔径累积分布曲线。3.根据权利要求1所述的岩心核磁共振T2谱转化孔径分布方法，其特征在于，所述将取出后的饱和充分的岩心样品进行核磁共振扫描获得T2谱累积分布曲线包括：获得岩心样品的弛豫时间T2和信号幅度的关系，弛豫时间T2的数值在一定程度上反应孔径的大小，信号幅度大小反应了对应弛豫时间T2值的孔径中所含油样的信号强度，信号幅度与弛豫时间T2围成的面积大小可反映所含油样的多少，以便获得不同弛豫时间T2数值范围内的峰面积占比，进而获得T2谱累积分布曲线。4.根据权利要求1所述的岩心核磁共振T2谱转化孔径分布方法，其特征在于，所述获得T2谱累积分布图版包括：将获得了T2谱从0ms开始的累积分布曲线，分别将T2谱的起始值Ts设置为5ms、10ms、15ms、20ms，然后分别作出不同弛豫时间T2起始值下的累积分布曲线，绘制在同一坐标系下，得到不同弛豫时间T2起始值下的T2谱累积分布图版。5.根据权利要求1所述的岩心核磁共振T2谱转化孔径分布方法，其特征在于，所述校对孔径累积分布与T2谱图版包括：根据关系式(1)获得横向弛豫时间T2和孔隙半径r存在一定的转化关系：其中V为孔隙体积；S为孔隙表面积；ρ
T
为弛豫率；F

【专利技术属性】
技术研发人员：黄世军，杨明洋，赵凤兰，陈新阳，孙浩月，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：