确定致密砂岩表面弛豫率的方法及装置与流程

1.本发明涉及致密砂岩油藏提高采收率技术领域，尤其涉及一种确定致密砂岩表面弛豫率的方法及装置。


背景技术：

2.致密砂岩储层孔隙类型以微-纳米级为主，流体渗流特征复杂，难以定量描述。低场核磁t2谱是一种用于测定孔隙内微观流体宏观统计规律的有效手段，弛豫时间t2与孔隙分布存在密切联系，可以一定程度上反应孔隙分布特征。
3.在均匀分布磁场中，不考虑扩散弛豫和自由弛豫的影响(相比于表面弛豫的影响可忽略)，低场核磁共振弛豫时间t2与孔隙半径可建立以下关系：
[0004][0005]
式中，t2为弛豫时间，ms；ρ为表面弛豫率，μm/s；s为岩心表面积，cm2；v为孔隙体积，cm3；r为孔隙半径，cm；c为常数，c＝1，2，3分别用于平板模型，毛细管束模型和球状模型，对于致密砂岩，一般选用毛细管束模型，即c＝2。
[0006]
因此，通过计算表面弛豫率可以将弛豫时间t2转换为孔隙半径，并利用t2谱表征孔隙分布规律。对于同一岩心样品，表面弛豫率计算方法大致可分为三类：(1)利用低场核磁测试结果直接计算，即使用cpmg或ir脉冲序列和有限的扩散脉冲序列测试(如脉冲场梯度和脉冲场梯度激发回波)这种方式；(2)结合孔隙半径(包括压汞法和氮气吸附法)与弛豫时间(低场核磁法)测试结果，计算表面弛豫率；(3)结合比表面积(包括氮气吸附、阳离子交换量和成像分析)与弛豫时间(低场核磁法)测试结果，计算表面弛豫率。
[0007]
其中，第一种方法是直接测量方法，不适用于非常规致密砂岩，因为在这类存在大量微纳米孔隙的岩心中，扩散弛豫现象影响显著，氢质子快速扩弛豫更复杂，对测试结果影响较大。第二种和第三种方法则是间接测量方法，其测定结果取决于使用何种测量手段评价致密砂岩的孔隙度，表面积和孔隙分布等。其中，适合致密砂岩孔隙分布评价的手段是高压压汞测量技术，结合高压压汞测定的孔隙半径与低场核磁测定的弛豫时间，可求解表面弛豫率。但由于高压压汞测试过程耗时耗力，操作流程复杂，且对岩心造成永久损害。


技术实现要素：

[0008]
本发明实施例提出一种确定致密砂岩表面弛豫率的方法，用以确定致密砂岩表面弛豫率，可操作性强，准确度高，该方法包括：
[0009]
获得饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；
[0010]
根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样品的拟t2截止值；
[0011]
确定不同转速对应的饱和油岩心样品的孔隙半径；
[0012]
根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率。
[0013]
本发明实施例提出一种确定致密砂岩表面弛豫率的装置，用以确定致密砂岩表面弛豫率，可操作性强，准确度高，该装置包括：
[0014]
t2谱获得模块，用于获得饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；
[0015]
拟t2截止值获得模块，用于根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样品的拟t2截止值；
[0016]
孔隙半径确定模块，用于确定不同转速对应的饱和油岩心样品的孔隙半径；
[0017]
致密砂岩表面弛豫率确定模块，用于根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率。
[0018]
本发明实施例还提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述确定致密砂岩表面弛豫率的方法。
[0019]
本发明实施例还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述确定致密砂岩表面弛豫率的方法的计算机程序。
[0020]
在本发明实施例中，获得饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样品的拟t2截止值；确定不同转速对应的饱和油岩心样品的孔隙半径；根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率。在上述过程中，直接采用低场核磁测量和高速离心处理过程获得第一核磁t2谱累积积分曲线和第二核磁t2谱累积积分曲线，从而确定饱和油岩心样品的拟t2截止值，在确定不同转速对应的饱和油岩心样品的孔隙半径后，最终获得致密砂岩表面弛豫率，整个过程避免了高压压汞测试过程，不对岩心造成永久损害，可操作性强，且准确度高。
附图说明
[0021]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本
发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
[0022]
图1为本发明实施例中确定致密砂岩表面弛豫率的方法的流程图；
[0023]
图2为本发明实施例确定致密砂岩表面弛豫率的方法的详细流程图；
[0024]
图3-图6分别为本发明实施例中致密砂岩柱塞样品a11-a14不同转速对应的第二核磁t2谱累积信号幅值的示意图；
[0025]
图7-图10分别为本发明实施例中致密砂岩柱塞样品a11-a14对应的拟合函数的示意图；
[0026]
图11-图14分别为本发明实施例中致密砂岩柱塞样品a11-a14采用不同方法获得的孔隙直径的对比示意图；
[0027]
图15为本发明实施例中确定致密砂岩表面弛豫率的装置的示意图；
[0028]
图16为本发明实施例中计算机设备的示意图。
具体实施方式
[0029]
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
[0030]
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本技术的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本技术的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
[0031]
图1为本发明实施例中确定致密砂岩表面弛豫率的方法的流程图，如图1所示，该方法包括：
[0032]
步骤101，获得饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；
[0033]
步骤102，根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样品的拟t2截止值；
[0034]
步骤103，确定不同转速对应的饱和油岩心样品的孔隙半径；
[0035]
步骤104，根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率。
[0036]
在本发明实施例提出的方法中，直接采用低场核磁测量和高速离心处理过程获得第一核磁t2谱累积积分曲线和第二核磁t2谱累积积分曲线，从而确定饱和油岩心样品的拟t2截止值，在确定不同转速对应的饱和油岩心样品的孔隙半径后，最终获得致密砂岩表面
弛豫率，整个过程避免了高压压汞测试过程，不对岩心造成永久损害，可操作性强，且准确度高。
[0037]
具体实施时，在步骤101中，进行低场核磁测量时，可以选用低场核磁装置，进行高速离心处理时，可以选用高速冷冻离心机中，在设定转速下离心2小时，获得设定转速对应的第二核磁t2谱累积积分曲线。之后，增加高速冷冻离心机的转速，重复上述步骤，直到测定的第二核磁t2谱累积积分曲线几乎不再变化(例如，t2谱累积积分面积变化小于3％)为止，这样可以得到多条不同转速对应的第二核磁t2谱累积积分曲线。
[0038]
在一实施例中，根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样品的拟t2截止值，包括：
[0039]
确定第一核磁t2谱累积积分曲线和不同转速对应的第二核磁t2谱累积积分曲线的交点对应的t2值为不同转速对应的饱和油岩心样品的拟t2截止值。
[0040]
在上述实施例中，对每条第二核磁t2谱累积积分曲线来说，将该第二核磁t2谱累积积分曲线的直线段反向延长，与第一核磁t2谱累积积分曲线相交，交点处的t2值即为拟t2截止值。上述确定拟t2截止值的方法简便，可操作性强。
[0041]
在一实施例中，确定不同转速对应的饱和油岩心样品的孔隙半径，包括：
[0042]
获得对饱和油岩心样品进行高速离心处理时的不同转速对应的离心力；
[0043]
根据不同转速对应的离心力，确定不同转速对应的饱和油岩心样品的孔隙半径。
[0044]
在上述实施例中，上述确定孔隙半径的原理时，在高速离心时，离心力等于饱和油岩心样品的毛管力，即不同转速对应的离心力为不同转速对应的毛管力，从而根据不同转速对应的毛管力，确定不同转速对应的饱和油岩心样品的孔隙半径，这种方法计算孔隙半径方便快捷，准确度也高。
[0045]
在一实施例中，采用如下公式，获得对饱和油岩心样品进行高速离心处理时的不同转速对应的离心力：
[0046][0047]
其中，p
c
为离心力，单位为mpa；l为饱和油岩心样品长度，单位为cm；r
e
为饱和油岩心样品旋转半径，单位为cm；δρ为油气两相密度差，单位为g/cm3；n为离心机转速，单位为r/min；
[0048]
采用如下公式，根据不同转速对应的离心力，确定不同转速对应的饱和油岩心样品的孔隙半径：
[0049][0050]
其中，p
ci
为毛管力，单位为mpa；σ为油气表面张力，单位为mn/m；θ为润湿角，单位为
°
；r为孔隙半径，单位为cm。
[0051]
在一实施例中，根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率，包括：
[0052]
根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定不同转速对应的表面弛豫率；
[0053]
根据不同转速对应的饱和油岩心样品的拟t2截止值和表面弛豫率，确定饱和油岩心样品的拟t2截止值和表面弛豫率的拟合函数；
[0054]
基于所述拟合函数，确定拟t2截止值为零时的表面弛豫率为饱和油岩心样品的致密砂岩表面弛豫率。
[0055]
在上述实施例中，可以通过线性拟合的方式确定多个拟t2截止值和多个表面弛豫率，形成拟合函数，另该拟合函数的拟t2截止值(即t2值坐标轴)为零时的表面弛豫率为饱和油岩心样品的致密砂岩表面弛豫率，即最终的致密砂岩表面弛豫率。
[0056]
在一实施例中，采用如下公式，根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定不同转速对应的表面弛豫率：
[0057][0058]
其中，t
2,pcutoff
为拟t2截止值，单位为ms；ρ为表面弛豫率，单位为μm/s；r为孔隙半径，单位为cm；c＝2。
[0059]
另外，还可以根据上述的致密砂岩表面弛豫率，将t2谱转化为孔隙半径分布，与平均值法和高压压汞法确定的孔隙分布结果进行对比，对比时，可以对属于同一致密砂岩柱塞大样品进行切片处理，第一份采用平均值法获得孔隙半径分布，第二份采用高压压汞法获得孔隙半径分布，第三份采用本发明提出的方法获得孔隙半径分布。通过对比可以确定本发明实施例获得的孔隙半径分布与平均值法获得的孔隙半径分布高度吻合，与高压压汞法总体趋势基本一致，因此，最后确定的表面弛豫率计算结果真实、可靠。
[0060]
采用平均值法获得孔隙半径分布后，采用平均值法计算表面弛豫率的公式如下：
[0061][0062][0063][0064]
其中，t
2lm
为弛豫时间对数平均值，单位为ms；r
p
为平均孔隙半径，单位为μm；r
i
为第i点孔隙半径，单位为μm；s
i
为第i点汞饱和度，单位为％。
[0065]
基于上述实施例，本发明提出如下一个实施例来说明确定致密砂岩表面弛豫率的方法的详细流程，图2为本发明实施例确定致密砂岩表面弛豫率的方法的详细流程图，包括：
[0066]
步骤201，获得饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；
[0067]
步骤202，确定第一核磁t2谱累积积分曲线和不同转速对应的第二核磁t2谱累积积分曲线的交点对应的t2值为不同转速对应的饱和油岩心样品的拟t2截止值；
[0068]
步骤203，获得对饱和油岩心样品进行高速离心处理时的不同转速对应的离心力；
[0069]
步骤204，根据不同转速对应的离心力，确定不同转速对应的饱和油岩心样品的孔隙半径；
[0070]
步骤205，根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定不同转速对应的表面弛豫率；
[0071]
步骤206，根据不同转速对应的饱和油岩心样品的拟t2截止值和表面弛豫率，确定饱和油岩心样品的拟t2截止值和表面弛豫率的拟合函数；
[0072]
步骤207，基于所述拟合函数，确定拟t2截止值为零时的表面弛豫率为饱和油岩心样品的致密砂岩表面弛豫率。
[0073]
当然，可以理解的是，上述详细流程还可以有其他变化例，相关变化例均应落入本发明的保护范围。
[0074]
下面给出一具体实施例，来说明本发明提出的方法的具体应用。
[0075]
选取鄂尔多斯盆地延长组主力开发层系长632四块致密砂岩柱塞样品，取芯深度2075-2225m，样品岩性主要为石英长石砂岩，粘土矿物组成包括伊利石(45.8％)，绿泥石(55.4％)和伊/蒙混层(8.1％)；
[0076]
致密砂岩柱塞样品经过洗油、烘干处理后，在致密砂岩柱塞样品端面截取一段(长度为1.8-2.0cm)用于高压压汞测试，另取一段(长度为1-1.2cm)进行接触角测试(接触角法)，剩余部分(包括a11-a14共4部分，长度为3.6-3.8cm)使用真空加压饱和装置进行抽真空饱和油处理(抽真空48小时，30mpa围压，煤油饱和5天)；
[0077]
使用micromr12-025v小尺寸核磁共振分析仪测定致密砂岩饱和油处理之后的致密砂岩柱塞样品a11-a14，获得每部分致密砂岩柱塞样品对应的第一核磁t2谱累积积分曲线，然后将每部分致密砂岩柱塞样品放置于csc-12(s)超级岩心高速冷冻离心机，在转速3000rpm/min下离心2小时，取出，测定转速3000rpm/min对应的第二核磁t2谱累积积分曲线；
[0078]
将转速增加至4000rpm，重复上述步骤，直到测定的第二核磁t2谱累积积分曲线形态几乎不再变化(t2谱累积积分面积变化小于3％)为止，图3-图6分别为本发明实施例中致密砂岩柱塞样品a11-a14不同转速对应的第二核磁t2谱累积信号幅值的示意图，其累积信号幅值即累积积分曲线，图3-图6中，按照从上到下的顺序累积积分曲线对应的转速分别为0rpm、1000rpm、3000rpm、5000rpm、6000rpm、7000rpm、8000rpm、9000rpm。
[0079]
对于每部分致密砂岩柱塞样品，对每条第二核磁t2谱累积积分曲线来说，将该第二核磁t2谱累积积分曲线的直线段反向延长，与第一核磁t2谱累积积分曲线相交，交点处的t2值即为拟t2截止值；
[0080]
根据公式(2)和公式(3)计算不同转速对应的每部分致密砂岩柱塞样品的孔隙半径。
[0081]
然后采用公式(4)计算不同转速对应的每部分致密砂岩柱塞样品的表面弛豫率，根据不同转速对应的饱和油岩心样品的拟t2截止值和表面弛豫率，确定饱和油岩心样品的拟t2截止值和表面弛豫率的拟合函数，图7-图10分别为本发明实施例中致密砂岩柱塞样品a11-a14对应的拟合函数的示意图；基于所述拟合函数，确定拟t2截止值为零时的表面弛豫率为饱和油岩心样品的致密砂岩表面弛豫率，从而可确定致密砂岩柱塞样品a11-a14的最
终的致密砂岩表面弛豫率。
[0082]
图11-图14分别为本发明实施例中致密砂岩柱塞样品a11-a14采用不同方法获得的孔隙直径的对比示意图，其中，本发明方法为图11-图14中的拟t2截止值法，表1为致密砂岩柱塞样品a11-a14采用不同方法获得的表面弛豫率的结果，两者误差在可接受的范围内。
[0083]
表1平均值法与本发明方法的表面弛豫率的结果
[0084][0085]
综上所述，在本发明实施例提出的方法中，获得饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样品的拟t2截止值；确定不同转速对应的饱和油岩心样品的孔隙半径；根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率。在上述过程中，直接采用低场核磁测量和高速离心处理过程获得第一核磁t2谱累积积分曲线和第二核磁t2谱累积积分曲线，从而确定饱和油岩心样品的拟t2截止值，在确定不同转速对应的饱和油岩心样品的孔隙半径后，最终获得致密砂岩表面弛豫率，整个过程避免了高压压汞测试过程，不对岩心造成永久损害，可操作性强，且准确度高。
[0086]
本发明实施例还提出一种确定致密砂岩表面弛豫率的装置，其原理与确定致密砂岩表面弛豫率的方法类似，这里不再赘述。
[0087]
图15为本发明实施例中确定致密砂岩表面弛豫率的装置的示意图，如图15所示，该装置包括：
[0088]
t2谱获得模块1501，用于获得饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；
[0089]
拟t2截止值获得模块1502，用于根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样
品的拟t2截止值；
[0090]
孔隙半径确定模块1503，用于确定不同转速对应的饱和油岩心样品的孔隙半径；
[0091]
致密砂岩表面弛豫率确定模块1504，用于根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率。
[0092]
在一实施例中，拟t2截止值获得模块1502具体用于：
[0093]
确定第一核磁t2谱累积积分曲线和不同转速对应的第二核磁t2谱累积积分曲线的交点对应的t2值为不同转速对应的饱和油岩心样品的拟t2截止值。
[0094]
在一实施例中，孔隙半径确定模块1503具体用于：
[0095]
获得对饱和油岩心样品进行高速离心处理时的不同转速对应的离心力；
[0096]
根据不同转速对应的离心力，确定不同转速对应的饱和油岩心样品的孔隙半径。
[0097]
在一实施例中，孔隙半径确定模块1503具体用于：
[0098]
采用如下公式，获得对饱和油岩心样品进行高速离心处理时的不同转速对应的离心力：
[0099][0100]
其中，p
c
为离心力，单位为mpa；l为饱和油岩心样品长度，单位为cm；r
e
为饱和油岩心样品旋转半径，单位为cm；δρ为油气两相密度差，单位为g/cm3；n为离心机转速，单位为r/min；
[0101]
采用如下公式，根据不同转速对应的离心力，确定不同转速对应的饱和油岩心样品的孔隙半径：
[0102][0103]
其中，p
ci
为毛管力，单位为mpa；σ为油气表面张力，单位为mn/m；θ为润湿角，单位为
°
；r为孔隙半径，单位为cm。
[0104]
在一实施例中，致密砂岩表面弛豫率确定模块1504具体用于：
[0105]
根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定不同转速对应的表面弛豫率；
[0106]
根据不同转速对应的饱和油岩心样品的拟t2截止值和表面弛豫率，确定饱和油岩心样品的拟t2截止值和表面弛豫率的拟合函数；
[0107]
基于所述拟合函数，确定拟t2截止值为零时的表面弛豫率为饱和油岩心样品的致密砂岩表面弛豫率。
[0108]
在一实施例中，致密砂岩表面弛豫率确定模块1504具体用于：
[0109]
采用如下公式，根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定不同转速对应的致密砂岩表面弛豫率：
[0110][0111]
其中，t
2,pcutoff
为拟t2截止值，单位为ms；ρ为表面弛豫率，单位为μm/s；r为孔隙半径，单位为cm；c＝2。
[0112]
综上所述，在本发明实施例提出的装置中，获得饱和油岩心样品的第一核磁t2谱
累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，所述饱和油岩心样品为进行抽真空饱和油处理后的致密砂岩柱塞样品，所述第一核磁t2谱累积积分曲线是对饱和油岩心样品进行低场核磁测量获得的；所述多条不同转速对应的第二核磁t2谱累积积分曲线是对采用不同转速进行高速离心处理后的饱和油岩心样品进行低场核磁测量获得的；根据饱和油岩心样品的第一核磁t2谱累积积分曲线和多条不同转速对应的第二核磁t2谱累积积分曲线，确定不同转速对应的饱和油岩心样品的拟t2截止值；确定不同转速对应的饱和油岩心样品的孔隙半径；根据不同转速对应的饱和油岩心样品的拟t2截止值和孔隙半径，确定致密砂岩表面弛豫率。在上述过程中，直接采用低场核磁测量和高速离心处理过程获得第一核磁t2谱累积积分曲线和第二核磁t2谱累积积分曲线，从而确定饱和油岩心样品的拟t2截止值，在确定不同转速对应的饱和油岩心样品的孔隙半径后，最终获得致密砂岩表面弛豫率，整个过程避免了高压压汞测试过程，不对岩心造成永久损害，可操作性强，且准确度高。
[0113]
本技术的实施例还提供一种计算机设备，图16为本发明实施例中计算机设备的示意图，该计算机设备能够实现上述实施例中的确定致密砂岩表面弛豫率的方法中全部步骤，所述电子设备具体包括如下内容：
[0114]
处理器(processor)1601、存储器(memory)1602、通信接口(communications interface)1603和总线1604；
[0115]
其中，所述处理器1601、存储器1602、通信接口1603通过所述总线1604完成相互间的通信；所述通信接口1603用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输；
[0116]
所述处理器1601用于调用所述存储器1602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的确定致密砂岩表面弛豫率的方法中的全部步骤。
[0117]
本技术的实施例还提供一种计算机可读存储介质，能够实现上述实施例中的确定致密砂岩表面弛豫率的方法中全部步骤，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的确定致密砂岩表面弛豫率的方法的全部步骤。
[0118]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0119]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0120]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指
令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0121]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0122]
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。