基于梯度热处理技术测量岩心束缚流体饱和度的方法与流程

1.本发明涉及勘探开发技术领域，具体涉及一种基于梯度热处理技术测量岩心束缚流体饱和度的方法。


背景技术：

2.对于典型的致密岩石，孔隙空间充满自由流体ff(free fluid)、毛管束缚流体caf(capillary bound fluid)和粘土束缚流体cbf(clay bound fluid)。自由流体是以游离态的形式存在，可以在孔道中自由流动的流体，相应的存储空间又称为自由孔隙。在油气藏成储的过程中，受孔隙壁表面润湿性差异和细小孔道毛管力的影响，油气在运移过程中无法驱替走所有的水，这些不能被驱替出孔隙的水就成了孔道中的束缚流体；其中分布和残存在岩石中细小非粘土矿物接触位置、微细孔隙中、颗粒表面和角隅处的称为毛管束缚流体，对应的孔隙空间称为毛管束缚孔隙；残留在粘土矿物颗粒表面的束缚流体称为粘土束缚流体，对应的孔隙空间称为粘土束缚孔隙。不同油藏由于岩石及流体性质不同，油气运移条件有差异，束缚水饱和度的大小差别较大。
3.低场核磁共振(nmr，nuclear magnetic 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)所占据。但是使用经验双t
2c
值无法正确表征所有岩石类型的孔隙大小，对于某些岩类的孔隙划分可能产生很大的误差。早期研究中有学者采用离心方法来确定两个t
2c
的值，但实验证明，离心所驱替出的流体通常只是自由流体，孔隙中仍然存在大量的毛管束缚流体和粘土束缚流体；采用多级离心确定双t
2c
来划分孔隙结构仍无法准确划分自由流体、毛管束缚流体和粘土束缚流体。因此，目前亟需一种能够准确划分及计算三种流体饱和度的方法。


技术实现要素：

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