致密岩石核磁共振T2截止值的拟合计算方法与流程

致密岩石核磁共振t2截止值的拟合计算方法
技术领域
1.本发明涉及油气田开发地质领域，特别是涉及到一种致密岩石核磁共振t2截止值的拟合计算方法。


背景技术：

2.核磁共振技术是表征储层孔隙结构并获取岩石可动流体饱和度的一项重要技术，具有安全无损快捷的优点，且可通过核磁共振测井应用于非取心井中。该技术通过t2弛豫时间获取岩样的孔隙结构及可动流体信息，理论上t2分布与孔隙大小存在线性关系，也就是说，较大的t2弛豫时间对应较大的孔隙半径。
3.确定样品的t2截止值是核磁共振实验中的一项关键技术。核磁共振t2截止值是区分可动流体与束缚流体的界限，饱和样品t2谱中大于t2截止值的部分为可动流体，反之则为束缚流体。这意味着t2截止值的确定直接关系到对储层可动流体的评价，进而影响到对储层评级及经济价值的评价结果。
4.目前确定t2截止值的方法主要包括经验值法、毛管压力法、形态分析法以及离心法等。早期研究中多采用固定的t2截止值，例如碳酸盐岩样品t2截止值一般取100ms，碎屑岩样品t2截止值一般取33ms，这种方法被称为经验值法。经验值法过程简单，但现有研究认为不同样品的t2截止值变化范围很大，对不同地区不同层位的储层使用固定t2截止值会导致较大的误差，这一现象在非常规储层中尤为明显。
5.为获取不同样品的t2截止值，一些研究借助压汞实验获取的毛管压力曲线确定t2截止值，即毛管压力法。这种方法中，先利用同一样品的压汞毛管压力曲线确定样品束缚水饱和度，再将束缚水饱和度与孔隙度相乘得到束缚流体含量，最后在饱和样品t2谱中寻找到该束缚流体含量对应的t2值，即为t2截止值。毛管压力法较为准确，但存在压力点确定问题。过大的压力点会导致将部分束缚流体计算为可动流体，过小的压力点则会导致确定的束缚水饱和度较高。另外，压汞实验会对岩心造成污染，该方法也无法应用在核磁共振测井中。
6.形态分析法则根据饱和t2谱的形态特征人工确定或根据经验公式计算每个样品的t2截止值，这种方法不需要进行岩心实验，可直接应用于核磁共振测井中，但该方法结果受人工对t2谱分类的影响，准确性较低，同时存在工作量较大的问题。
7.离心法中，首先对饱含流体样品进行核磁共振测试，获取饱和样品的t2分布谱，然后对样品在一定条件下进行离心并对离心后样品进行核磁共振测试，获取离心样品的t2分布谱，最后通过对比饱含流体及离心后样品的t2谱曲线获取t2截止值。该方法结果最为准确，也是目前应用最为广泛的方法。但离心法需要进行大量的离心操作，耗费的时间和人力成本较大，同时也存在结果无法直接应用于核磁共振测井这一问题。
8.在申请号：cn201410555235.7的中国专利申请中，涉及到一种核磁共振测井t2谱中t2截止值的确定方法，尤其是一种依据核磁共振测井t2谱各时间点t2谱强度所反映的孔隙结构特征构建t2截止值的方法。本发明能节省原本大量岩心的实验分析标定t2截止值的
时间和成本；同时，也不用建立在t2谱分类基础上经验公式求取t2截止值，实现的t2截止值动态自动求取，该方法实现了科学简便快速求取t2截止值。
9.在申请号：cn201811203489.7的中国专利申请中，涉及到一种核磁共振测井t2截止值的逐点计算方法，对目标工区的岩样核磁共振实验数据，分别计算核磁共振t2谱的截止时间几何平均值gc和伪几何平均值g'，构建gc、g'数据矩阵后，由最小二乘法算法确定目标工区的伪t2截止值t2'c。由t2'c计算核磁共振t2谱的伪截止时间几何平均值g'c，确定t2截止值t2cutoff与g'c之间的拟合关系式。计算目标工区的实测核磁共振测井资料的g'c，再利用前述确定拟合关系式，实现核磁共振测井t2截止值的逐点计算方法。
10.在申请号：cn200510072232.9的中国专利申请中，涉及到一种核磁共振测井t2谱t2截止值的确定方法，包括以下步骤：步骤1、对核磁共振测井t2谱按形态特征进行分类；步骤2、根据不同类型的测井t2谱拟合得出其离心谱；步骤3、将离心谱进行面积积分，得出离心谱的面积；步骤4、计算测井t2谱的面积；步骤5、将测井t2谱的面积进行累加，当累加面积等于或与离心谱面积最相近时，此时对应的t2值即为t2截止值。
11.以上现有技术与本发明并不相同，不能解决核磁共振t2截止值确定及在核磁共振测井上应用的问题，我们发明了一种核磁共振t2截止值的拟合计算方法。


技术实现要素：

12.本发明的目的是提供一种致密岩石核磁共振t2截止值的拟合计算方法，以减少离心过程耗费的时间和人力成本，同时该方法可直接应用于核磁共振测井，为油气田的勘探与开发提供技术支撑。
13.本发明的目的可通过如下技术措施来实现：致密岩石核磁共振t2截止值的拟合计算方法，该致密岩石核磁共振t2截止值的拟合计算方法包括：步骤1，将样品进行预处理；步骤2，开展常规物性测试实验，获取样品的孔隙度和渗透率；步骤3，开展饱和样品核磁共振实验，获取饱和样品t2分布数据；步骤4，开展离心后样品核磁共振试验，获取离心样品t2分布数据；步骤5，建立核磁共振t2截止值拟合公式。
14.本发明的目的还可通过如下技术措施来实现：
15.在步骤1中，对样品进行预处理，包括钻取柱状岩样、洗油和烘干。
16.在步骤1中，切取长度和直径均为2.5cm的柱状样品并进行洗油处理，之后将样品置于烘箱中至样品恒重备用。
17.在步骤1中，洗油处理采用蒸馏抽提法，溶剂选用65％三氯甲烷和35％甲醇的混合溶液，为防止粘土矿物转化，烘干温度不超过63℃。
18.在步骤2中，对样品进行常规物性测试，确定样品的孔隙度和渗透率，常规储层样品渗透率测试采用稳态法，非常规储层样品渗透率测试采用脉冲衰减法。
19.在步骤3中，开展饱和样品核磁共振测试，获取饱和样品t2分布数据，其中饱和流体根据地层水的矿化度和类型配备。
20.在步骤4中，开展离心样品核磁共振测试，获取离心样品t2分布数据，其中常规储层离心力使用0.689mpa、100psi，非常规储层离心力根据样品情况确定，要求可动水含量随离心力增大不再明显减少。
21.在步骤5中，利用离心样品t2分布数据计算离心后样品的累积孔隙度，并在饱和样
截止值，计算过程依据《岩样核磁共振参数实验室测量规范》sy/t 6490-2014标准规定执行，具体过程为：利用离心样品t2分布数据计算离心后样品的累积孔隙度，并在饱和样品t2谱上寻找一点，使其左侧累积孔隙度与离心样品总累积孔隙度相等，该点对应t2值即为t2截止值。之后利用步骤2获取的岩心孔隙度和渗透率，利用数学回归分析方法，得出t2截止值与储层常规物性参数的转换关系，即核磁共振t2截止值的拟合计算公式。核磁共振t2截止值拟合计算公式见公式(i)。
[0038][0039]
其中t
2cutoff
为t2截止值，ms；为孔隙度，％；k为渗透率，md；a、b和c为常数，通过回归求得。
[0040]
为更加明确本发明效果，在应用本发明的实施例中，以胜利油田某区块超低渗-特低渗储层为具体实例详细说明具体过程。
[0041]
本发明中实例样品取自渤海湾盆地惠民凹陷中央隆起带某区块沙河街组基山砂岩储层中，储集层孔隙度为3.6％-18.1％，空气渗透率为0.009-5.464md。表1为该实例储层样品孔隙度、渗透率、实测t2截止值及计算t2截止值数据。
[0042]
表1胜利油田某区块基山砂岩储层核磁共振t2截止值拟合数据表
[0043]
序号孔隙度(％)渗透率(md)实测t2截止值计算t2截止值111.680.27110.9323.06217.063.03642.3734.44315.332.58953.6736.2949.890.02625.665.46517.121.85626.5929.9368.840.07633.4717.11716.221.26928.2928.2489.660.13013.1620.31915.711.34736.1229.73109.430.1037.8018.661113.710.42010.9223.111218.052.69128.2831.481317.871.76620.6428.07148.650.0796.0917.791510.870.62535.2632.04166.860.07022.6920.13173.650.00911.288.02188.840.0605.5714.93
[0044]
步骤1，首先在岩性均匀部分钻取直径和长度均为2.5cm的规则柱状岩样，并利用三氯甲烷(65)和甲醇(35)的混合溶液对样品进行洗油处理，样品清洗过程采用蒸馏抽提法，最大限度的去除样品中的残留油和地层水。之后将柱状样品置于烘箱中进行烘干至样品恒重备用(每2小时使用高精度天平称重一次，连续3次测量样品质量变化量小于0.01g)，
以去除样品中的残留水和易挥发杂质，为防止粘土矿物发生转化，烘干温度为63℃。
[0045]
步骤2，开展常规物性测试工作，获取样品的孔隙度和渗透率。由于样品渗透率较低，属于超低渗-特低渗储层，渗透率测试选用脉冲衰减法，实验在coretest公司的ap-608自动渗透率-孔隙度测试仪上进行，围压选用3.45mpa(500psi)。
[0046]
步骤3，开展饱和样品核磁共振实验，实验在纽迈公司的mesomr23-060h-i型核磁共振仪上进行，获取饱和样品t2分布数据(图2)。其中饱和流体按照该实例储层地层水的矿化度和类型配备，为保证较高的信噪比，实验中使用的参数主要为：回波间隔(te)为0.1ms，完全恢复时间(tw)为3s，扫描次数(ns)为128，回波个数(nech)为8192。
[0047]
步骤4，对完成步骤3中核磁共振测试的饱和样品进行离心，离心过程在卢湘仪公司的gl-21m高速冷冻离心机上进行。由于该实例储层属于非常规储层，经多次测量确定离心机转速为10000r/min，对应的离心力为3.785mpa(549psi)，离心时长为150分钟。对离心后样品开展核磁共振实验，实验在纽迈公司的mesomr23-060h-i型核磁共振仪上进行，获取离心样品t2分布数据(图2)，实验中参数同步骤3中饱和样品核磁共振测试。
[0048]
步骤5中，利用步骤3和4中获取的饱和样品和离心样品t2分布数据计算实测t2截止值，具体过程为：利用离心样品t2分布数据计算离心后样品的累积孔隙度，并在饱和样品t2谱上寻找一点，使其左侧累积孔隙度与离心样品总累积孔隙度相等，该点对应t2值即为t2截止值(图3)。之后利用步骤2获取的岩心孔隙度和渗透率，利用数学回归分析方法，得出该实例储层核磁共振t2截止值的拟合计算公式：
[0049][0050]
其中t
2cutoff
为t2截止值，ms；为孔隙度，％；k为渗透率，md。
[0051]
图4为该实例储层的核磁共振t2截止值拟合效果图，图5为该实例储层的计算t2截止值与实测t2截止值对比图，结果显示该方法应用效果较好。
[0052]
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域技术人员来说，其依然可以对前述实施例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。除说明书所述的技术特征外，均为本专业技术人员的已知技术。