基于核磁共振双T2截止值确定岩心束缚流体分布的方法与流程

基于核磁共振双t2截止值确定岩心束缚流体分布的方法
技术领域
1.本发明涉及勘探开发技术领域，具体涉及一种基于核磁共振双t2截止值确定岩心束缚流体分布的方法。


背景技术：

2.对于典型的致密岩石，孔隙空间充满自由流体ff(free fluid)、毛管束缚流体caf(capillarybound fluid)和粘土束缚流体cbf(claybound fluid)。自由流体是以游离态的形式存在，可以在孔道中自由流动的流体，相应的存储空间又称为自由孔隙。在油气藏成储的过程中，受孔隙壁表面润湿性差异和细小孔道毛管力的影响，油气在运移过程中无法驱替走所有的水，这些不能被驱替出孔隙的水就成了孔道中的束缚流体；其中分布和残存在岩石中细小非粘土矿物接触位置、微细孔隙中、颗粒表面和角隅处的称为毛管束缚流体，对应的孔隙空间称为毛管束缚孔隙；残留在粘土矿物颗粒表面的束缚流体称为粘土束缚流体，对应的孔隙空间称为粘土束缚孔隙。不同油藏由于岩石及流体性质不同，油气运移条件有差异，束缚水饱和度的大小差别较大。
3.低场核磁共振(nmr，nuclearmagnetic resonance)测量得到的横向弛豫时间(t2)可以表征岩石孔隙中流体的分布状态，然而，如何根据nmr测量得到的t2分布确定t2截止值(t
2c
)以准确界定自由流体和束缚流体含量是目前研究的难点问题。
4.目前，t2截止值的确定主要包括两种模型：第一种是常规单t
2c
模型，该方法通常根据核磁共振t2谱曲线的凹面最低点处对应的值划分，将小于该低点处对应的值部分划分为束缚流体，高于该低点对应的值部分划分为自由流体。但是单t
2c
不是划分自由流体与束缚流体的准确值，目前在物理学上也很难找到一个弛豫时间阈值，高于这个阈值的弛豫时间对应的孔隙流体能够自由产出而无残余流体，低于这个阈值的弛豫时间所对应的孔隙流体就难以产出。所以用单t2截止值来划分孔隙类型都会引入许多不确定性。第二种是双t
2c
模型，由于单t2截止值的局限性，后人提出了双t2截止值模型来划分孔隙类型。早期研究通常采用经验t
2c
值来快速划分孔隙类型，如致密砂岩中常取t
2c1
为10ms、t
2c2
为3ms，划分孔隙为自由孔隙、毛管束缚孔隙、粘土束缚孔隙三类；每类孔隙被相应的流体(自由流体ff、毛管束缚流体caf、粘土束缚流体cbf)所占据。
5.但是使用经验双t
2c
值无法正确表征所有岩石类型的孔隙大小，对于某些岩类的孔隙划分可能产生很大的误差。早期研究中有学者采用离心方法来确定两个t
2c
的值，但实验证明，离心所驱替出的流体通常只是自由流体可动部分，孔隙中仍然存在大量的毛管束缚流体和粘土束缚流体；采用多级离心确定双t
2c
来划分孔隙结构仍无法准确划分自由流体、毛管束缚流体和粘土束缚流体。因此，目前亟需一种能够准确计算三种流体饱和度的方法。
6.上述方法既无法准确得出三种流体的饱和度，又无法确定三种流体在孔隙中的分布状，因此，目前亟需一种能够准确计算三种流体饱和度，同时能够确定三种流体在孔隙中的分布的方法。


技术实现要素：

7.针对现有技术中的缺陷，本发明的目的是提供一种基于核磁共振双t2截止值确定岩心束缚流体分布的方法。
8.本发明所采用的技术方案是：
9.基于核磁共振双t2截止值确定岩心束缚流体分布的方法，包括如下步骤，
10.s1、获取岩样干重m
d
和饱和岩样的质量m
s
；
11.s2、采用等时差梯度增温方法对岩样做热处理，增温次数为i，记录每次增温后的温度，并称量每次增温后的岩样质量，分别得到第i次增温后的热处理的温度t
i
以及第i次增温后的岩样质量m
i
；
12.s3、根据岩样干重以及每次增温后的岩样质量，计算每次增温后岩样的孔隙流体饱和度，得到第i次增温后的流体饱和度sw
i
，并绘制流体饱和度sw
i
随温度t
i
变化的关系曲线，即sw
i
‑
t
i
关系曲线；
13.s4、对sw
i
‑
t
i
关系曲线中的每个流体饱和度sw
i
求导，得到每个流体饱和度sw
i
的一阶导数sw
i’，并绘制sw
i
的一阶导数sw
i’随时间t
i
变化的关系曲线，即sw
i
’‑
t
i
关系曲线；
14.s5、对sw
i
’‑
t
i
关系曲线中的sw
i
的一阶导数sw
i’求导，得到每个饱和度sw
i
的二阶导数sw
i”，并绘制sw
i
的二阶导数sw
i”随时间t
i
变化的关系曲线，即sw
i
”‑
t
i
关系曲线；
15.s6、根据sw
i
”‑
t
i
关系曲线特征的差异性，得到三种流体的截止温度；
16.s7、重新烘干并饱和岩样，测量饱和岩样的核磁t2谱，记作t2
ff
；
17.s8、根据步骤s6得到的三种流体的截止温度，将饱和岩样依次置于对应的三种流体的截止温度下进行热处理，并测量每次热处理后的核磁t2谱，分别记作t2
caf
和t2
cbf
；
18.s9、根据所述t2
ff
、t2
caf
和t2
cbf
求取双t2截止值t2c1和t2c2；
19.s10、根据所述t2c1和t2c2，将t2谱划分为自由孔隙、毛管束缚孔隙和粘土束缚孔隙。
20.进一步，所述步骤s2中，所述采用等时差梯度增温方法对岩样做热处理，具体采用如下控温方程：
21.t
i
＝t0 i*δt；
22.其中，t0为起始温度，δt为加热时长，δt为升温梯度，i为增温次数，t
i
为第i次增温后的热处理的温度。
23.进一步，所述步骤s3中，根据所述岩样质量和岩样干重，计算每次热处理后岩样的孔隙流体饱和度sw
i
，具体采用如下公式：
24.sw
i
＝(m
i
‑
m
d
)/(m
s
‑
m
d
)*100％；
25.其中，sw
i
为第i次增温后的流体饱和度，m
i
为第i次增温后的岩样质量，m
d
为岩样干重，m
s
为饱和岩样的质量。
26.进一步，所述步骤s4中，对sw
i
‑
t
i
关系曲线中的每个流体饱和度sw
i
求导，得到每个流体饱和度sw
i
的一阶导数sw
i’，具体采用如下公式：
27.sw
i’＝(sw
i
‑
sw
(i
‑
1)
)/δt*100％；
28.其中，i>＝1。
29.进一步，所述步骤s5中，对sw
i
’‑
t
i
关系曲线中的sw
i
的一阶导数sw
i’求导，得到每个饱和度sw
i
的二阶导数sw
i”，具体采用如下公式：
30.sw
i”＝(sw
i
’‑
sw
(i
‑
1)’)/δt*100％；
31.其中，i>＝2。
32.进一步，所述三种流体分别为自由流体、毛管束缚流体和粘土束缚流体。
33.进一步，所述步骤s6中，所述根据sw
i
”‑
t
i
关系曲线特征的差异性，得到三种流体的截止温度，具体为，
34.根据sw
i
”‑
t
i
关系曲线呈现出的明显不同的三个区间，得到对应三个区间的三种流体，三个区间的交点分别为三种流体的分界点，该三种流体的分界点对应的温度t
i
为三种流体的截止温度；其中自由流体和毛管束缚流体的截止温度记作t
ff
‑
caf
，毛管束缚流体和粘土束缚流体的截止温度记作t
caf
‑
cbf
。
35.进一步，所述步骤s8具体包括如下步骤：
36.s801、将饱和岩样置于烘箱内，阶梯增温至截止温度t
ff
‑
caf
，并在该截止温度t
ff
‑
caf
加热时长为δt，取出并测量核磁t2谱，此时岩样中剩余流体的t2谱反映总束缚流体的分布，记作t2
caf
；
37.s801、再次将岩样置于烘箱内，阶梯增温至截止温度t
caf
‑
cbf
，并在该截止温度t
caf
‑
cbf
加热时长为δt，取出并测量核磁t2谱，此时样品中剩余流体的t2谱反映粘土束缚流体的分布，记作t2
cbf
。
38.进一步，所述步骤s9中，根据所述t2
ff
、t2
caf
和t2
cbf
求取双t2截止值t2c1和t2c2，具体包括如下步骤：
39.s901、分别绘制t2
ff
、t2
caf
和t2
cbf
三条曲线的累积曲线c
ff
、c
caf
和c
cbf
，t2信号累积的最大值分别记为c
ff(max)
、c
caf(max)
和c
cbf(max)
；
40.s902、根据步骤s901绘制的累积曲线做两条投影线，第一条投影线起自t2
caf
累积谱的最大值c
caf(max)
，终点为该第一条投影线与t2
ff
累积曲线的交点；第二条投影线起自t2
cbf
累积谱的最大值c
cbf(max)
，终点为该第二条投影线与t2
ff
累积曲线的交点；
41.s903、分别过第一条投影线与t2
ff
累积曲线的交点和第二条投影线与t2
ff
累积曲线的交点作t2谱横坐标的垂向投影线，与横坐标的交点从右往左依次为t2c1和t2c2，其中t2c1为自由流体t2截止值，t2c2为毛管束缚流体截止值。
42.进一步，所述步骤s10中，根据所述t2c1和t2c2，将t2谱划分为自由孔隙、毛管束缚孔隙和粘土束缚孔隙，具体为，
43.按t2谱，将t2c1右侧部分的孔隙划分为自由孔隙；将位于t2c1与t2c2之间的部分划分为毛管束缚孔隙；将t2c2左侧部分划分为粘土束缚孔隙。
44.本技术的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本技术的实践了解到。
附图说明
45.为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中，类似的元件或部分一般由类似的附图标记标识。附图中，各元件或部分并不一定按照实际的比例绘制。
46.图1为本技术实施例所提供的流体饱和度sw
i
随温度t
i
变化的关系曲线图(sw
i
‑
t
i
关系曲线)；
47.图2为本技术实施例所提供的sw
i
的一阶导数sw
i’随时间t
i
变化的关系曲线图(sw
i
’‑
t
i
关系曲线)；
48.图3为本技术实施例所提供的sw
i
的二阶导数sw
i”随时间t
i
变化的关系曲线图(sw
i
”‑
t
i
关系曲线)；
49.图4为本技术实施例所提供的t2谱示意图；
50.图5为本技术实施例所提供的累积曲线示意图；
51.图6为本技术实施例所提供的累积曲线的两条投影线示意图；
52.图7为本技术实施例所提供的投影线与累积曲线的交点的垂向投影线示意图；
53.图8为本技术实施例所提供的束缚流体分布示意图。
具体实施方式
54.下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案，因此只作为示例，而不能以此来限制本发明的保护范围。
55.需要注意的是，除非另有说明，本技术使用的技术术语或者科学术语应当为本发明所属领域技术人员所理解的通常意义。
56.本技术基于核磁共振双t2截止值确定岩心束缚流体分布的方法，包括如下步骤，
57.s1、烘干岩样，称量并得到岩样干重m
d
；具体为，将岩样放置于烘箱中烘干，烘干温度为200℃，烘干时间为24小时。
58.配置模拟地层水溶液，将岩样置于该地层水溶液中加压饱和，使岩样孔隙中完全充满地层水，称量并得到饱和岩样的质量m
s
。
59.s2、采用等时差梯度增温方法对岩样做热处理，增温次数为i，记录每次增温后的热处理的温度，并称量每次增温后的岩样质量，分别得到第i次增温后的热处理的温度t
i
以及第i次增温后的岩样质量m
i
。
60.在采用等时差梯度增温方法对岩样做热处理时，具体采用如下控温方程：
61.t
i
＝t0 i*δt；
62.其中，t0为起始温度，δt为加热时长，δt为升温梯度，i为增温次数，t
i
为第i次增温后的热处理的温度。
63.本实施例中，t0取20℃，δt取20分钟，δt取20℃，t
i
的最大值限定为200℃。
64.在一种示例性实施例中，第1次增温，将岩样置于烘箱中，设置温度为40℃(t1为40℃)，加热20分钟后，测量岩样质量，得到第1次增温后的岩样质量m1。
65.第2次增温，再次将岩样置于烘箱中，设置温度为60℃(t2为60℃)，加热20分钟后，测量岩样质量，得到第2次增温后的岩样质量m2。
66.第3次增温，再次将岩样置于烘箱中，设置温度为80℃(t3为80℃)，加热20分钟后，测量岩样质量，得到第3次增温后的岩样质量m3；
67.按照上述增温方法，根据控温方程，将升温梯度δt设置为20℃，加热时长δt设置为20分钟，对岩样进行热处理，并记录每次热处理的温度以及增温后的岩样质量，直至热处理的温度达到200℃时停止。
68.s3、根据岩样干重以及每次增温后的岩样质量，计算每次增温后岩样的孔隙流体
饱和度，得到第i次增温后的流体饱和度sw
i
，并绘制流体饱和度sw
i
随温度t
i
变化的关系曲线，即sw
i
‑
t
i
关系曲线(参见图1)。
69.在根据岩样质量和岩样干重，计算每次热处理后岩样的孔隙流体饱和度sw
i
，具体采用如下公式：
70.sw
i
＝(m
i
‑
m
d
)/(m
s
‑
m
d
)*100％；
71.其中，sw
i
为第i次增温后的流体饱和度，m
i
为第i次增温后的岩样质量，m
d
为岩样干重，m
s
为饱和岩样的质量；当i＝0时，m
i
＝m
s
，此时sw0＝100％。
72.s4、对sw
i
‑
t
i
关系曲线中的每个流体饱和度sw
i
求导，得到每个流体饱和度sw
i
的一阶导数sw
i’，并绘制sw
i
的一阶导数sw
i’随时间t
i
变化的关系曲线，即sw
i
’‑
t
i
关系曲线(参见图2)。
73.对sw
i
‑
t
i
关系曲线中的每个流体饱和度sw
i
求导，得到每个流体饱和度sw
i
的一阶导数sw
i’，具体采用如下公式：
74.sw
i’＝(sw
i
‑
sw
(i
‑
1)
)/δt*100％；
75.其中，sw
i’为sw
i
的一阶导数，sw
i
为第i次增温后的流体饱和度，sw
(i
‑
1)
为第i
‑
1次增温后的流体饱和度，i>＝1。
76.s5、对sw
i
’‑
t
i
关系曲线中的sw
i
的一阶导数sw
i’求导，得到每个饱和度sw
i
的二阶导数sw
i”，并绘制sw
i
的二阶导数sw
i”随时间t
i
变化的关系曲线，即sw
i
”‑
t
i
关系曲线(参见图3)。
77.对sw
i
’‑
t
i
关系曲线中的sw
i
的一阶导数sw
i’求导，得到每个饱和度sw
i
的二阶导数sw
i”，具体采用如下公式：
78.sw
i”＝(sw
i
’‑
sw
(i
‑
1)’)/δt*100％；
79.其中，sw
i”为sw
i
的二阶导数，sw
i’为sw
i
的一阶导数，即第i次增温后的流体饱和度的导数，sw
(i
‑
1)’为sw
i
‑1的一阶导数，即第i
‑
1次增温后的流体饱和度的导数，i>＝2。
80.s6、根据swi
”‑
t
i
关系曲线特征的差异性，得到三种流体的截止温度。三种流体分别为自由流体、毛管束缚流体和粘土束缚流体。
81.根据swi
”‑
t
i
关系曲线特征的差异性，得到三种流体的截止温度，具体为，根据swi
”‑
t
i
关系曲线呈现出的明显不同的三个区间，得到对应三个区间的三种流体，三个区间的交点分别为三种流体的分界点，该三种流体的分界点对应的温度t
i
为三种流体的截止温度。
82.在本实施例中，swi
”‑
t
i
关系曲线如图3所示，由图3中，可以看出swi
”‑
t
i
关系曲线具有明显不同的三个区间，根据同一种流体斜率接近的特点，可以很容易的将该三个区间对应于三种流体，从左至右分别为自由流体、毛管束缚流体和粘土束缚流体；对应区间的交点即分别为自由流体、毛管束缚流体和粘土束缚流体的分界点。其中，t
ff
‑
caf
为自由流体和毛管束缚流体的截止温度，t
caf
‑
cbf
为毛管束缚流体和粘土束缚流体的截止温度。
83.s7、重新烘干并饱和岩样，测量饱和岩样的核磁t2谱，记作t2
ff
。
84.具体的，将岩样放置于烘箱中烘干，烘干温度为200℃，烘干时间为24小时。
85.配置模拟地层水溶液，将岩样置于该地层水溶液中加压饱和，使岩样孔隙中完全充满地层水，测量饱和岩样的核磁t2谱，记作t2
ff
。
86.s8、根据步骤s6得到的三种流体的截止温度，将饱和岩样依次置于对应的三种流
体的截止温度下进行热处理，并测量每次热处理后的核磁t2谱，分别记作t2
caf
和t2
cbf
。
87.具体的，s801、将饱和岩样置于烘箱内，阶梯增温至截止温度t
ff
‑
caf
，并在该截止温度t
ff
‑
caf
下加热时长为δt，取出并测量核磁t2谱，此时岩样中剩余流体的t2谱反映总束缚流体的分布，记作t2
caf
；
88.s802、再次将岩样置于烘箱内，阶梯增温至截止温度t
caf
‑
cbf
，并在该截止温度t
caf
‑
cbf
下加热时长为δt，取出并测量核磁t2谱，此时样品中剩余流体的t2谱反映粘土束缚流体的分布，记作t2
cbf
。
89.具体在阶梯增温时，可按步骤s2的等时差梯度增温方法重新进行增温热处理。即t0取20℃，每次增温20℃,加热时长δt取20分钟，直至温度达到t
ff
‑
caf
或t
caf
‑
cbf
。
90.在一种示例性实施例中，假设t
ff
‑
caf
为75℃，t
caf
‑
cbf
为110℃；将饱和岩样置于烘箱内，加热温度依次设为40℃、60℃和75℃，每次加热20分钟，加热完成后，取出测量核磁t2谱；再次将岩样置于烘箱内，加热温度依次设为80℃、100℃和110℃，每次加热20分钟，加热完成后，取出并测量核磁t2谱。
91.s9、根据所述t2
ff
、t2
caf
和t2
cbf
求取双t2截止值t2c1和t2c2；具体包括如下步骤：
92.s901、如图4和图5所示，分别绘制t2
ff
、t2
caf
和t2
cbf
三条曲线的累积曲线c
ff
、c
caf
和c
cbf
，t2信号累积的最大值分别记为c
ff(max)
、c
caf(max)
和c
cbf(max)
。
93.其中，c
ff(max)
＝sum(t2
ff(1)
,t2
ff(2)
,
……
，t2
ff(x)
)；
94.c
caf(max)
＝sum(t2
caf(1)
,t2
caf(2)
,
……
，t2
caf(x)
)；
95.c
cbf(max)
＝sum(t2
cbf(1)
,t2
cbf(2)
,
……
，t2
cbf(x)
)；
96.其中，x为核磁t2测试的布点数，t2
ff(x)
表示t2
ff
横坐标上第x个点对应的纵坐标幅度值，t2
caf(x)
表示t2
caf
横坐标上第x个点对应的纵坐标幅度值，t2
cbf(x))
表示t2
cbf
横坐标上第x个点对应的纵坐标幅度值。
97.s902、如图6所示，根据步骤s901绘制的累积曲线做两条投影线，第一条投影线起自t2
caf
累积谱的最大值c
caf(max)
，终点为该第一条投影线与t2
ff
累积曲线的交点(标记为
①
)；第二条投影线起自t2
cbf
累积谱的最大值c
cbf(max)
，终点为该第二条投影线与t2
ff
累积曲线的交点(标记为
②
)。
98.s903、如图7所示，分别过第一条投影线与t2
ff
累积曲线的交点和第二条投影线与t2
ff
累积曲线的交点作t2谱横坐标的垂向投影线(标记为和
③
和
④
)，与横坐标的交点从右往左依次为t2c1和t2c2，其中t2c1为自由流体t2截止值，t2c2为毛管束缚流体截止值。
99.s10、如图8所示，根据所述t2c1和t2c2，将t2谱划分为自由孔隙、毛管束缚孔隙和粘土束缚孔隙。
100.具体的，按t2谱，将t2c1右侧部分(即t2>t2c1)的孔隙划分为自由孔隙；将位于t2c1与t2c2之间的部分划分为毛管束缚孔隙；将t2c2左侧部分(t2<t2c2)划分为粘土束缚孔隙。
101.本技术利用不同温度阶段脱出水的类型不同，在曲线上的不同温度阶段呈现不同的曲线形态，采用等时差梯度增温方法对岩样做热处理，能够得到每次热处理后的岩样质量，从而计算出每次热处理后的流体饱和度，并根据每次热处理后的流体饱和度绘制流体饱和度sw
i
随温度t
i
变化的关系曲线(sw
i
‑
t
i
关系曲线)，通过对该关系曲线二次求导，可以得到斜率变化趋势明显不同的三个区间，根据同一种流体斜率接近的特点可以很容易的得
到三种不同流体的区间以及截止温度；根据该截止温度，再次进行热处理，可以分别得到岩样在该截止温度下的核磁t2谱，根据该核磁t2谱即可获得该岩样的双t2截止值以及三种流体的分布；本技术处理方法简单，易于操作，并且三种流体的区间表现明显，划分准确性高，利用三种流体的截止温度，可以更准确的得到双t2截止值，从而准确的划分流体分布。
102.在本技术中，除非另有明确的规定和限定，术语“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
103.本发明的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、系统和技术，以便不模糊对本说明书的理解。
104.在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、系统、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、系统、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
105.最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。