CN113790996A - 基于离心法测量岩心束缚流体饱和度的方法 - Google Patents
基于离心法测量岩心束缚流体饱和度的方法 Download PDFInfo
- Publication number
- CN113790996A CN113790996A CN202111027008.3A CN202111027008A CN113790996A CN 113790996 A CN113790996 A CN 113790996A CN 202111027008 A CN202111027008 A CN 202111027008A CN 113790996 A CN113790996 A CN 113790996A
- Authority
- CN
- China
- Prior art keywords
- saturation
- rock sample
- bound
- fluid
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 164
- 239000011435 rock Substances 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000010438 heat treatment Methods 0.000 claims abstract description 45
- 239000011148 porous material Substances 0.000 claims abstract description 36
- 238000009795 derivation Methods 0.000 claims abstract description 14
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 10
- 238000009738 saturating Methods 0.000 claims abstract description 4
- 239000004927 clay Substances 0.000 claims description 16
- 238000005303 weighing Methods 0.000 claims description 13
- 238000005119 centrifugation Methods 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 238000001035 drying Methods 0.000 claims description 5
- 238000003672 processing method Methods 0.000 abstract description 2
- 239000000243 solution Substances 0.000 description 5
- 238000005481 NMR spectroscopy Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000007864 aqueous solution Substances 0.000 description 2
- 239000002734 clay mineral Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000010792 warming Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000018044 dehydration Effects 0.000 description 1
- 238000006297 dehydration reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/02—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by absorbing or adsorbing components of a material and determining change of weight of the adsorbent, e.g. determining moisture content
- G01N5/025—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by absorbing or adsorbing components of a material and determining change of weight of the adsorbent, e.g. determining moisture content for determining moisture content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/04—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder
- G01N5/045—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder for determining moisture content
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明公开了一种基于离心法测量岩心束缚流体饱和度的方法,包括如下步骤,S1、得到岩样干重;S2、将岩样加压饱和,得到饱和岩样的质量;S3、采用离心法对岩样做脱水处理;S4对岩样做热处理,记录每次增温后的岩样质量;S5、计算岩样的孔隙束缚流体饱和度;S6、对Swi‑Ti关系曲线中的束缚流体饱和度Swi求导;S7、对Swi’‑Ti关系曲线中的Swi的一阶导数Swi’求导,绘制Swi”‑Ti关系曲线;S8、根据Swi”‑Ti关系曲线,得到两种类型束缚流体的截止温度;S9、计算得到两种束缚流体的饱和度;本申请处理方法简单,易于操作,并且两种束缚流体的区间表现明显,划分准确性高。
Description
技术领域
本发明涉及勘探开发技术领域,具体涉及一种基于离心法测量岩心束缚流体饱和度的方法。
背景技术
对于典型的致密岩石,孔隙空间充满自由流体FF(free fluid)、毛管束缚流体CAF(capillary bound fluid)和粘土束缚流体CBF(clay bound fluid)。自由流体是以游离态的形式存在,可以在孔道中自由流动的流体,相应的存储空间又称为自由孔隙。在油气藏成储的过程中,受孔隙壁表面润湿性差异和细小孔道毛管力的影响,油气在运移过程中无法驱替走所有的水,这些不能被驱替出孔隙的水就成了孔道中的束缚流体;其中分布和残存在岩石中细小非粘土矿物接触位置、微细孔隙中、颗粒表面和角隅处的称为毛管束缚流体,对应的孔隙空间称为毛管束缚孔隙;残留在粘土矿物颗粒表面的束缚流体称为粘土束缚流体,对应的孔隙空间称为粘土束缚孔隙。不同油藏由于岩石及流体性质不同,油气运移条件有差异,束缚水饱和度的大小差别较大。
低场核磁共振(NMR,nuclear magnetic resonance)测量得到的横向弛豫时间(T2)可以表征岩石孔隙中流体的分布状态,然而,如何根据NMR测量得到的T2分布确定T2截止值(T2C)以准确界定自由流体和束缚流体含量是目前研究的难点问题。
目前,T2截止值的确定主要包括两种模型:第一种是常规单T2c模型,该方法通常根据核磁共振T2谱曲线的凹面最低点处对应的值划分,将小于该低点处对应的值部分划分为束缚流体,高于该低点对应的部分划分为自由流体。但是单T2C不是划分自由流体与束缚流体的准确值,目前在物理学上也很难找到一个弛豫时间阈值,高于这个阈值的弛豫时间对应的孔隙流体能够自由产出而无残余流体,低于这个阈值的弛豫时间所对应的孔隙流体就难以产出。所以用单T2截止值来划分孔隙类型会引入许多不确定性。第二种是双T2c模型,由于单T2截止值的局限性,后人提出了双T2截止值模型来划分孔隙类型。早期研究通常采用经验T2C值来快速划分孔隙类型,如致密砂岩中常取T2C1为10ms、T2C2为3ms,划分孔隙为自由孔隙、毛管束缚孔隙、粘土束缚孔隙三类;每类孔隙被相应的流体(自由流体FF、毛管束缚流体CAF、粘土束缚流体CBF)所占据。但是使用经验T2C值无法正确表征所有岩石类型的孔隙大小,对于某些岩类的孔隙划分可能产生很大的误差。早期研究中有学者采用离心方法来确定两个T2C的值,但实验证明,离心所驱替出的流体通常只是自由流体,孔隙中仍然存在大量的毛管束缚流体和粘土束缚流体;采用多级离心确定双T2C来划分孔隙结构仍无法准确划分自由流体、毛管束缚流体和粘土束缚流体。因此,目前亟需一种能够准确划分及计算三种流体饱和度的方法。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于离心处理测量岩心束缚流体饱和度的方法。
本发明所采用的技术方案是:
基于离心法测量岩心束缚流体饱和度的方法,包括如下步骤,
S1、烘干岩样,称量并得到岩样干重md;
S2、将岩样置于地层水溶液中加压饱和,使岩样孔隙中完全充满地层水,称量并得到饱和岩样的质量ms;
S3、采用离心法对岩样做脱水处理,称量并得到离心后样品的质量mc,并根据该离心后样品的质量计算自由流体的饱和度;
S4、采用等时差梯度增温方法对离心后的岩样做热处理,增温次数为i,记录每次增温后的温度,并称量每次增温后的岩样质量,分别得到第i次增温后的热处理的温度Ti以及第i次增温后的岩样质量mi;
S5、根据岩样干重以及每次增温后的岩样质量,计算每次增温后岩样的孔隙束缚流体饱和度,得到第i次增温后的束缚流体饱和度Swi,并绘制束缚流体饱和度Swi随温度Ti变化的关系曲线,即Swi-Ti关系曲线;
S6、对Swi-Ti关系曲线中的每个束缚流体饱和度Swi求导,得到每个束缚流体饱和度Swi的一阶导数Swi’,并绘制Swi的一阶导数Swi’随时间Ti变化的关系曲线,即Swi’-Ti关系曲线;
S7、对Swi’-Ti关系曲线中的Swi的一阶导数Swi’求导,得到每个束缚流体饱和度Swi的二阶导数Swi”,并绘制Swi的二阶导数Swi”随时间Ti变化的关系曲线,即Swi”-Ti关系曲线;
S8、根据Swi”-Ti关系曲线特征的差异性,得到两种束缚流体的截止温度;
S9、根据步骤S4得到的Swi-Ti关系曲线以及步骤S7得到的两种束缚流体的截止温度,分别计算得到两种束缚流体的饱和度。
进一步,所述步骤S3中,所述根据该离心后样品的质量计算自由流体的饱和度,具体采用如下公式:
SFF=(mC-md)/(ms-md)*100%;
其中,SFF为自由流体饱和度,mC为离心后样品的质量,md为岩样干重,ms为饱和岩样的质量。
进一步,所述步骤S4中,所述采用等时差梯度增温方法对岩样做热处理,具体采用如下控温方程:
Ti=T0+i*ΔT;
其中,T0为起始温度,Δt为加热时长,ΔT为升温梯度,i为增温次数,Ti为第i次增温后的热处理的温度。
进一步,所述步骤S5中,根据所述岩样质量和岩样干重,计算每次热处理后岩样的孔隙束缚流体饱和度Swi,具体采用如下公式:
Swi=(mi-md)/(ms-md)*100%;
其中,Swi为第i次增温后的束缚流体饱和度,mi为第i次增温后的岩样质量,md为岩样干重,ms为饱和岩样的质量。
进一步,所述步骤S6中,对Swi-Ti关系曲线中的每个束缚流体饱和度Swi求导,得到每个束缚流体饱和度Swi的一阶导数Swi’,具体采用如下公式:
Swi’=(Swi-Sw(i-1))/ΔT*100%;
其中,i>=1。
进一步,所述步骤S7中,对Swi’-Ti关系曲线中的Swi的一阶导数Swi’求导,得到每个束缚流体饱和度Swi的二阶导数Swi”,具体采用如下公式:
Swi”=(Swi’-Sw(i-1)’)/ΔT*100%;
其中,i>=2。
进一步,所述两种束缚流体分别毛管束缚流体和粘土束缚流体。
进一步,所述步骤S8中,所述根据Swi”-Ti关系曲线特征的差异性,得到两种束缚流体的截止温度,具体为,
根据Swi”-Ti关系曲线呈现出的明显不同的两个区间,得到对应两个区间的两种束缚流体,两个区间的交点分别为两种束缚流体的分界点,该两种束缚流体的分界点对应的温度Ti为两种束缚流体的截止温度。
本申请的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本申请的实践了解到。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中,类似的元件或部分一般由类似的附图标记标识。附图中,各元件或部分并不一定按照实际的比例绘制。
图1为本申请实施例所提供的束缚流体饱和度Swi随温度Ti变化的关系曲线图(Swi-Ti关系曲线);
图2为本申请实施例所提供的Swi的一阶导数Swi’随时间Ti变化的关系曲线图(Swi’-Ti关系曲线);
图3为本申请实施例所提供的Swi的二阶导数Swi”随时间Ti变化的关系曲线图(Swi”-Ti关系曲线)。
具体实施方式
下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案,因此只作为示例,而不能以此来限制本发明的保护范围。
需要注意的是,除非另有说明,本申请使用的技术术语或者科学术语应当为本发明所属领域技术人员所理解的通常意义。
本申请基于离心法测量岩心束缚流体饱和度的方法,包括如下步骤,
S1、烘干岩样,称量并得到岩样干重md;具体为,将岩样放置于烘箱中烘干,烘干温度为200℃,烘干时间为24小时。
S2、配置模拟地层水溶液,将岩样置于该地层水溶液中加压饱和,使岩样孔隙中完全充满地层水,称量并得到饱和岩样的质量ms。
S3、采用离心法对岩样做脱水处理,称量并得到离心后样品的质量mc,并根据该离心后样品的质量计算自由流体的饱和度。
采用离心处理驱替出的水均为自由流体,而毛管束缚流体和粘土束缚流体均留存在岩样中,无法采用离心法驱替出。
通过称量离心后岩样的质量mc,根据如下公式可计算出自由流体饱和度:
SFF=(mC-md)/(ms-md)*100%;
其中,SFF为自由流体饱和度,mC为离心后样品的质量,md为岩样干重,ms为饱和岩样的质量。
则束缚流体饱和度SCAF+CBF=100-SFF;其中,SCAF-CBF为束缚流体饱和度,SFF为自由流体饱和度。
S4、采用等时差梯度增温方法对岩样做热处理,增温次数为i,记录每次增温后的热处理的温度,并称量每次增温后的岩样质量,分别得到第i次增温后的热处理的温度Ti以及第i次增温后的岩样质量mi。
由于离心法无法驱替出束缚流体,因此,在本实施例中,采用等时差梯度增温方法对离心后的岩样进行热处理,以驱替出岩样中的束缚流体。
在采用等时差梯度增温方法对岩样做热处理时,具体采用如下控温方程:
Ti=T0+i*ΔT;
其中,T0为起始温度,Δt为加热时长,ΔT为升温梯度,i为增温次数,Ti为第i次增温后的热处理的温度。
本实施例中,T0取20℃,Δt取20分钟,ΔT取20℃,Ti的最大值限定为200℃。
在一种示例性实施例中,第1次增温,将岩样置于烘箱中,设置温度为40℃(T1为40℃),加热20分钟后,测量岩样质量,得到第1次增温后的岩样质量m1。
第2次增温,再次将岩样置于烘箱中,设置温度为60℃(T2为60℃),加热20分钟后,测量岩样质量,得到第2次增温后的岩样质量m2。
第3次增温,再次将岩样置于烘箱中,设置温度为80℃(T3为80℃),加热20分钟后,测量岩样质量,得到第3次增温后的岩样质量m3;
按照上述增温方法,根据控温方程,将升温梯度ΔT设置为20℃,加热时长Δt设置为20分钟,对岩样进行热处理,并记录每次热处理的温度以及增温后的岩样质量,直至热处理的温度达到200℃时停止。
S4、根据岩样干重以及每次增温后的岩样质量,计算每次增温后岩样的孔隙束缚流体饱和度,得到第i次增温后的束缚流体饱和度Swi,并绘制束缚流体饱和度Swi随温度Ti变化的关系曲线,即Swi-Ti关系曲线(参见图1)。
在根据岩样质量和岩样干重,计算每次热处理后岩样的孔隙束缚流体饱和度Swi,具体采用如下公式:
Swi=(mi-md)/(ms-md)*100%;
其中,Swi为第i次增温后的束缚流体饱和度,mi为第i次增温后的岩样质量,md为岩样干重,ms为饱和岩样的质量;当i=0时,mi=mc,此时Sw0=SFF。
S5、对Swi-Ti关系曲线中的每个束缚流体饱和度Swi求导,得到每个束缚流体饱和度Swi的一阶导数Swi’,并绘制Swi的一阶导数Swi’随时间Ti变化的关系曲线,即Swi’-Ti关系曲线(参见图2)。
对Swi-Ti关系曲线中的每个束缚流体饱和度Swi求导,得到每个束缚流体饱和度Swi的一阶导数Swi’,具体采用如下公式:
Swi’=(Swi-Sw(i-1))/ΔT*100%;
其中,Swi’为Swi的一阶导数,Swi为第i次增温后的束缚流体饱和度,Sw(i-1)为第i-1次增温后的束缚流体饱和度,i>=1。
S6、对Swi’-Ti关系曲线中的Swi的一阶导数Swi’求导,得到每个束缚流体饱和度Swi的二阶导数Swi”,并绘制Swi的二阶导数Swi”随时间Ti变化的关系曲线,即Swi”-Ti关系曲线(参见图3)。
对Swi’-Ti关系曲线中的Swi的一阶导数Swi’求导,得到每个饱和度Swi的二阶导数Swi”,具体采用如下公式:
Swi”=(Swi’-Sw(i-1)’)/ΔT*100%;
其中,Swi”为Swi的二阶导数,Swi’为Swi的一阶导数,即第i次增温后的束缚流体饱和度的导数,Sw(i-1)’为Swi-1的一阶导数,即第i-1次增温后的束缚流体饱和度的导数,i>=2。
S7、根据Swi”-Ti关系曲线特征的差异性,得到两种束缚流体的截止温度。两种束缚流体分别为毛管束缚流体和粘土束缚流体。
根据Swi”-Ti关系曲线特征的差异性,得到两种束缚流体的截止温度,具体为,根据Swi”-Ti关系曲线呈现出的明显不同的两个区间,得到对应两个区间的两种束缚流体,两个区间的交点分别为两种束缚流体的分界点,该两种束缚流体的分界点对应的温度Ti为两种束缚流体的截止温度。
在本实施例中,Swi”-Ti关系曲线如图3所示,由图3中,可以看出Swi”-Ti关系曲线具有明显不同的两个区间,根据同一种流体斜率接近的特点,可以很容易的将该两个区间对应于两种束缚流体,从左至右分别为毛管束缚流体和粘土束缚流体;对应区间的交点即分别为毛管束缚流体和粘土束缚流体的分界点,即TCAF-CBF为毛管束缚流体和粘土束缚流体的截止温度。
S8、根据步骤S4得到的Swi-Ti关系曲线以及步骤S7得到的两种束缚流体的截止温度,分别计算得到两种束缚流体的饱和度。
具体为,将TCAF-CBF代入Swi-Ti关系曲线中(参见图1),计算得到粘土束缚流体饱和度SCBF,根据离心法的质量直接计算得到的自由流体饱和度SFF,毛管束缚流体饱和度SCAF=100-SCBF-SFFF。
本申请利用不同温度阶段脱出水的类型不同,在曲线上的不同温度阶段呈现不同的曲线形态,首先采用离心法以脱水的方式驱替出岩样中的自由流体,并直接计算出自由流体的饱和度,再采用等时差梯度增温方法对岩样做热处理,能够得到每次热处理后的岩样质量,从而计算出每次热处理后的束缚流体饱和度,并根据每次热处理后的束缚流体饱和度绘制束缚流体饱和度Swi随温度Ti变化的关系曲线(Swi-Ti关系曲线),通过对该关系曲线二次求导,可以得到斜率变化趋势明显不同的两个区间,根据同一种流体斜率接近的特点可以很容易的得到两种不同类型束缚流体的区间以及截止值,从而可以更加简单的得到岩样中束缚流体饱和度,处理方法简单,易于操作,并且两种束缚流体的区间表现明显,划分准确性高。
在本申请中,除非另有明确的规定和限定,术语“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
本发明的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、系统和技术,以便不模糊对本说明书的理解。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、系统、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、系统、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (8)
1.基于离心法测量岩心束缚流体饱和度的方法,其特征在于,包括如下步骤,
S1、烘干岩样,称量并得到岩样干重md;
S2、将岩样置于地层水溶液中加压饱和,使岩样孔隙中完全充满地层水,称量并得到饱和岩样的质量ms;
S3、采用离心法对岩样做脱水处理,称量并得到离心后样品的质量mc,并根据该离心后样品的质量计算自由流体的饱和度;
S4、采用等时差梯度增温方法对离心后的岩样做热处理,增温次数为i,记录每次增温后的温度,并称量每次增温后的岩样质量,分别得到第i次增温后的热处理的温度Ti以及第i次增温后的岩样质量mi;
S5、根据岩样干重以及每次增温后的岩样质量,计算每次增温后岩样的孔隙束缚流体饱和度,得到第i次增温后的束缚流体饱和度Swi,并绘制束缚流体饱和度Swi随温度Ti变化的关系曲线,即Swi-Ti关系曲线;
S6、对Swi-Ti关系曲线中的每个束缚流体饱和度Swi求导,得到每个束缚流体饱和度Swi的一阶导数Swi’,并绘制Swi的一阶导数Swi’随时间Ti变化的关系曲线,即Swi’-Ti关系曲线;
S7、对Swi’-Ti关系曲线中的Swi的一阶导数Swi’求导,得到每个束缚流体饱和度Swi的二阶导数Swi”,并绘制Swi的二阶导数Swi”随时间Ti变化的关系曲线,即Swi”-Ti关系曲线;
S8、根据Swi”-Ti关系曲线特征的差异性,得到两种束缚流体的截止温度;
S9、根据步骤S4得到的Swi-Ti关系曲线以及步骤S7得到的两种束缚流体的截止温度,分别计算得到两种束缚流体的饱和度。
2.根据权利要求1所述的基于离心法测量岩心束缚流体饱和度的方法,其特征在于,所述步骤S3中,所述根据该离心后样品的质量计算自由流体的饱和度,具体采用如下公式:
SFF=(mC-md)/(ms-md)*100%;
其中,SFF为自由流体饱和度,mC为离心后样品的质量,md为岩样干重,ms为饱和岩样的质量。
3.根据权利要求1所述的基于离心法测量岩心束缚流体饱和度的方法,其特征在于,所述步骤S4中,所述采用等时差梯度增温方法对岩样做热处理,具体采用如下控温方程:
Ti=T0+i*ΔT;
其中,T0为起始温度,Δt为加热时长,ΔT为升温梯度,i为增温次数,Ti为第i次增温后的热处理的温度。
4.根据权利要求2所述的基于离心法测量岩心束缚流体饱和度的方法,其特征在于,所述步骤S5中,根据所述岩样质量和岩样干重,计算每次热处理后岩样的孔隙束缚流体饱和度Swi,具体采用如下公式:
Swi=(mi-md)/(ms-md)*100%;
其中,Swi为第i次增温后的束缚流体饱和度,mi为第i次增温后的岩样质量,md为岩样干重,ms为饱和岩样的质量。
5.根据权利要求4所述的基于离心法测量岩心束缚流体饱和度的方法,其特征在于,所述步骤S6中,对Swi-Ti关系曲线中的每个束缚流体饱和度Swi求导,得到每个束缚流体饱和度Swi的一阶导数Swi’,具体采用如下公式:
Swi’=(Swi-Sw(i-1))/ΔT*100%;
其中,i>=1。
6.根据权利要求5所述的基于离心法测量岩心束缚流体饱和度的方法,其特征在于,所述步骤S7中,对Swi’-Ti关系曲线中的Swi的一阶导数Swi’求导,得到每个束缚流体饱和度Swi的二阶导数Swi”,具体采用如下公式:
Swi”=(Swi’-Sw(i-1)’)/ΔT*100%;
其中,i>=2。
7.根据权利要求1所述的基于离心法测量岩心束缚流体饱和度的方法,其特征在于,所述两种束缚流体分别毛管束缚流体和粘土束缚流体。
8.根据权利要求7所述的基于离心法测量岩心束缚流体饱和度的方法,其特征在于,所述步骤S8中,所述根据Swi”-Ti关系曲线特征的差异性,得到两种束缚流体的截止温度,具体为,
根据Swi”-Ti关系曲线呈现出的明显不同的两个区间,得到对应两个区间的两种束缚流体,两个区间的交点分别为两种束缚流体的分界点,该两种束缚流体的分界点对应的温度Ti为两种束缚流体的截止温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111027008.3A CN113790996A (zh) | 2021-09-02 | 2021-09-02 | 基于离心法测量岩心束缚流体饱和度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111027008.3A CN113790996A (zh) | 2021-09-02 | 2021-09-02 | 基于离心法测量岩心束缚流体饱和度的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113790996A true CN113790996A (zh) | 2021-12-14 |
Family
ID=78879475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111027008.3A Pending CN113790996A (zh) | 2021-09-02 | 2021-09-02 | 基于离心法测量岩心束缚流体饱和度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113790996A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5299453A (en) * | 1993-01-28 | 1994-04-05 | Mobil Oil Corporation | Method for determining oil and water saturation of core samples at overburden pressure |
CN104697915A (zh) * | 2015-03-20 | 2015-06-10 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩微观孔隙大小及流体分布的分析方法 |
CN106525688A (zh) * | 2016-11-21 | 2017-03-22 | 中国石油大学(华东) | 一种饱和页岩孔隙流体分离及饱和度计算的实验方法 |
CN109725016A (zh) * | 2018-11-29 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种用于含重油、沥青质岩心的核磁共振实验室测量方法 |
CN110618170A (zh) * | 2019-09-26 | 2019-12-27 | 中国石油天然气股份有限公司 | 岩土的饱和度指数的确定方法及装置 |
CN111537543A (zh) * | 2020-06-03 | 2020-08-14 | 中国矿业大学 | 低场核磁共振测定页岩黏土与脆性矿物相对含量的方法 |
CN112540096A (zh) * | 2020-11-27 | 2021-03-23 | 武汉大学 | 饱和冻结岩石未冻束缚水和未冻自由水含量的获取方法 |
-
2021
- 2021-09-02 CN CN202111027008.3A patent/CN113790996A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5299453A (en) * | 1993-01-28 | 1994-04-05 | Mobil Oil Corporation | Method for determining oil and water saturation of core samples at overburden pressure |
CN104697915A (zh) * | 2015-03-20 | 2015-06-10 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩微观孔隙大小及流体分布的分析方法 |
CN106525688A (zh) * | 2016-11-21 | 2017-03-22 | 中国石油大学(华东) | 一种饱和页岩孔隙流体分离及饱和度计算的实验方法 |
CN109725016A (zh) * | 2018-11-29 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种用于含重油、沥青质岩心的核磁共振实验室测量方法 |
CN110618170A (zh) * | 2019-09-26 | 2019-12-27 | 中国石油天然气股份有限公司 | 岩土的饱和度指数的确定方法及装置 |
CN111537543A (zh) * | 2020-06-03 | 2020-08-14 | 中国矿业大学 | 低场核磁共振测定页岩黏土与脆性矿物相对含量的方法 |
CN112540096A (zh) * | 2020-11-27 | 2021-03-23 | 武汉大学 | 饱和冻结岩石未冻束缚水和未冻自由水含量的获取方法 |
Non-Patent Citations (3)
Title |
---|
孙玉立: "基于NMR的致密砂岩N2-CO2吞吐增采机理研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅰ辑(月刊)》 * |
杨明 等: "不同阶煤孔隙结构与流体特性的核磁共振试验研究", 《中国安全科学学报》 * |
陈琼: "黏土吸附结合水动力学模型及机理研究", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅱ辑 (月刊)》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107991335B (zh) | 致密砂岩水锁损害评价试验方法 | |
CN111795993B (zh) | 一种高温高压下岩石热物性瞬态测试系统及方法 | |
CN113790996A (zh) | 基于离心法测量岩心束缚流体饱和度的方法 | |
Kumar et al. | Modelling of breather for transformer health assessment | |
Grunewald et al. | Towards an engineering model of material characteristics for input to ham transport simulations-Part 1: an approach | |
CN112986091B (zh) | 基于电阻率测试的层理缝发育页岩的含油饱和度评价方法 | |
CN113790997A (zh) | 基于梯度热处理技术测量岩心束缚流体饱和度的方法 | |
CN113791102B (zh) | 基于核磁共振双t2截止值确定岩心束缚流体分布的方法 | |
Salager et al. | Effect of temperature on water retention phenomena in deformable soils: theoretical and experimental aspects | |
CN115470722B (zh) | 一种获取多孔介质型碎片床干涸热流密度的方法及系统 | |
Tao et al. | Effects of Contact Angle on the Hysteresis Effect of Soil‐Water Characteristic Curves during Dry‐Wet Cycles | |
CN112049623A (zh) | 一种中基性火山岩束缚水出水的判定方法和装置 | |
CN112613174B (zh) | 考虑多重吸附机制的页岩吸附甲烷能力评估方法 | |
Guimarães et al. | A tool to predict water absorption in porous building materials | |
RU2535527C1 (ru) | Способ определения количественного состава многокомпонентной среды (варианты) | |
Berger et al. | On The Integrationof Hygrothermal Bridges Into Whole Building Ham Modeling | |
Harun et al. | Drying comparison of nonhygroscopic and hygroscopic materials | |
Harun et al. | Drying simulation of ceramic shell build up process | |
CN113657054B (zh) | 乳化堵塞损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 | |
US20170254764A1 (en) | Pressurizable dilatometer system and apparatus | |
Knight et al. | A laboratory procedure for estimating irreducible water saturation from cuttings | |
Plaisted | Characterization of soil unsaturated flow properties using steady state centrifuge methods | |
CN116449102A (zh) | 一种岩石电阻率的函数描述方法 | |
Salinas et al. | Three-dimensional simulation of isothermal wood drying of radiata pine using effective diffusion coefficients | |
Lv et al. | Time‐Dependent Creep Damage Model of Frozen Soil with Unsteady Fractional Order under Temperature and Stress Coupling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |