CN110346258B - 一种致密岩石油相相对渗透率测定的方法 - Google Patents
一种致密岩石油相相对渗透率测定的方法 Download PDFInfo
- Publication number
- CN110346258B CN110346258B CN201910647504.5A CN201910647504A CN110346258B CN 110346258 B CN110346258 B CN 110346258B CN 201910647504 A CN201910647504 A CN 201910647504A CN 110346258 B CN110346258 B CN 110346258B
- Authority
- CN
- China
- Prior art keywords
- rock sample
- rock
- oil
- section
- water
- 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.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 165
- 230000035699 permeability Effects 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 49
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 52
- 238000002474 experimental method Methods 0.000 claims abstract description 23
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 15
- 238000010521 absorption reaction Methods 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 14
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims abstract description 8
- 229910052753 mercury Inorganic materials 0.000 claims abstract description 8
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 229910021380 Manganese Chloride Inorganic materials 0.000 claims abstract description 6
- 239000011565 manganese chloride Substances 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000005520 cutting process Methods 0.000 claims abstract description 4
- 238000003825 pressing Methods 0.000 claims abstract description 4
- 238000005303 weighing Methods 0.000 claims abstract description 4
- 230000002431 foraging effect Effects 0.000 claims abstract 2
- 239000011148 porous material Substances 0.000 claims description 47
- 239000012530 fluid Substances 0.000 claims description 23
- 230000003595 spectral effect Effects 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 9
- 238000001035 drying Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 5
- 239000002390 adhesive tape Substances 0.000 claims description 3
- 230000032683 aging Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000005311 nuclear magnetism Effects 0.000 claims description 3
- 238000009738 saturating Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000004140 cleaning Methods 0.000 claims 1
- 230000008859 change Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000005406 washing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 238000011160 research Methods 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
- 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/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- 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
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Fluid Mechanics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种致密岩石油相相对渗透率测定方法,步骤如下:S1、取致密岩石的圆柱形岩样;S2、将岩样切割为三段;S3、对第一段岩样进行油水毛管压力实验;S4、对第二段岩样进行恒速压汞实验;S5、测量第三段岩样的孔隙度φHe、渗透率k、长度L、直径D、干重m0和密度ρ;S6、将第三段岩样抽真空加压饱和MnCl2溶液;S7、将饱和完成的第三段岩样开展油驱水实验,取出岩心进行老化;S8、将老化后的第三段岩样擦去表面油后称重记为m2,同时用核磁共振仪获取岩样的T2谱曲线和核磁孔隙度;然后进行自吸水排油实验;最后进行数据处理。本发明的方法测试实验数据计量精度高,操作简单,充分利用的致密岩石的特性及相关渗流理论,具有较高的适用性。
Description
技术领域
本发明涉及油气勘探开发技术领域,特别是一种致密岩石油相相对渗透率测定方法。
背景技术
致密油作为非常规油气勘探开发领域的一大新热点,已然成为我国未来重要的能源接替之一。油相相对渗透率曲线是认识致密油气藏油水流动特征的重要参数之一,有效确定岩石油相相对渗透率对油藏的采出程度、分析油井产水规律、提高采收率以及剩余油饱和度的研究均有重要意义。因此,众学者开展了大量针对油藏相对渗透率的研究,在分析多相流体间相互作用关系的基础上,形成了多种油相相对渗透率的确定的方法,但概括起来可分为室内实验、数值模拟和经验公式。
室内实验又包含稳态法和非稳态法两种,稳态法的基本理论是一维达西渗流理论,并且忽略了毛管压力和重力作用,假设两相流体互不相溶且不可压缩。试验时在总流量不变的条件下,将油水按一定流量比例同时恒速注入岩样,当进口、出口压力及油、水流量稳定时,岩样含水饱和度不再变化,此时可利用测定岩样进口、出口压力及油、水流量,由达西定律直接计算出岩样的油、水有效渗透率及相对渗透率值。非稳态法主要是以Buckley-Leverett一维两相水驱油前缘推进理论为基础的,同样该方法也忽略毛管压力和重力作用,假设两相流体互不相溶且不可压缩,岩样任一横截面内油水饱和度是均匀的。实验时将岩心事先饱和另一种流体,用另一种流体进行驱替,在岩样出口端记录每种流体的产量和岩样两端的压力差随时间的变化,最后用JBN方法计算得到油-水相对渗透率。
数值模拟包含主要通过历史拟合或者孔隙网络模拟两种方法,历史拟合是主要是通过设定油相和水相相对渗透率与含水饱和度间存在某种函数关系,该函数关系中有一定的待定参数,给定待定参数的初始值,利用差分方法求解给定条件下岩心内流体运动的渗流模型,计算岩心出口端的累计产油量(产水量)和驱替压差,不断调节待定参数的数值使得计算出的累计产油量和驱替压差与实测值之间的离差平方和最小,从而确定最优的一组待定参数。孔隙网络模拟主要是通过构建与多孔介质相对应得孔隙网络模型老研究储层的传导性质和模拟多孔介质内多相流的宏观物理性质(如相对渗透率、毛管压力等)。
经验公式是国内外学者通过统计大量的矿场资料和室内实验的结果,建立了相对渗透率与其它物理性质间的关系,通过获取储层的其它相关物理性质求解相对渗透率。
但是,现有的技术存在以下不足:
(1)稳态法虽然测试过程和数据处理简单,但该方法中岩心内流体的运动过程并不能代表油藏中实际的流体的运动,且该方法需要较长的测试周期,对仪器和计量装置的精度要求较高。
(2)非稳态法测试周期短,驱替过程与实际油藏相符,故而被广泛应用与室内实验,但该方法忽略了毛管压力影响,认为流体在岩心内是以等饱和面推进,而致密储层发育有大量的纳米级孔喉,储层孔隙结构复杂并非均质性强,且具有显著的毛管效应,因此,对于致密储层而言,利用非稳态的驱替数据获取的相对渗透率曲线往往出现异常。
(3)历史拟合法虽原理简单,但实际计算十分繁琐,计算量较大,需要大量的生产数据和储层评价结果。孔隙网络模拟的准确与否取决于对多孔介质孔隙结构和多孔介质内多相流体间的运输规律的认知,且该方法计算复杂,对计算机性能较高。
(4)经验公式是通过统计获取,每个公式均具有对应的适用对象,且含有大量的经验参数,因此,该方法可信度较低。
综上所述,由于致密储层岩石表现出较高的毛管压力,这种现象将使得岩心具有较强的自发吸入湿相流体的能力。现有的相对渗透率测定方法不能准确的测量致密储层岩石的相对渗透率。
发明内容
本发明的目的针对现有的相对渗透率测量方法不能准确测量致密储层岩石油相相对渗透率的技术缺陷,提供一种致密岩石油相相对渗透率测定方法。
本发明提供的致密岩石油相相对渗透率测定方法,包括如下步骤:
S1、取长度7-11cm,直径2-3cm的致密岩石的圆柱形岩样,并将岩样进行洗油、洗盐、烘干等处理。
S2、将清洗烘干的岩样切割为三段,第一段和第二段长度为2-3cm,第三段长度为3-5cm,然后再将三段岩样进行烘干。
S3、按国家标准GB/T 29171-2012《岩石毛管压力曲线的测定》对第一段岩样进行油水毛管压力实验,测定岩石毛管压力曲线,毛管压力记为Pcow。
S4、对第二段岩样进行恒速压汞实验测定孔喉分布特征。
S6、测量完成后,先将第三段岩样抽真空至133Pa,在20MPa压力下加压饱和MnCl2溶液48小时后,MnCl2溶液的浓度大于50000ppm,测量岩样的质量,记为m1,计算岩样的有效孔隙体积Veff和有效孔隙度φeff,当岩样的有效孔隙度与气测孔隙度的相对误差小于2%时,即满足时,岩样饱和完成,否则按上述步骤重新饱和。计算岩样的有效孔隙体积Veff和有效孔隙度φeff的公式如下:
S7、将饱和完成的第三段岩样放入岩心加持器中,连接驱替装置,并在岩心夹持器出口端接油水分离计量管;在恒压状态下,用实验用油开展油驱水实验,待驱至含油率达到99%时停止驱替,卸载围压和内压,取出岩心放入实验用油中进行老化,老化时间不少于10d。
S8、将老化后的第三段岩样擦去表面油后称重记为m2,同时用核磁共振仪获取岩样的T2谱曲线和核磁孔隙度,其T2谱曲线的面积记为A0,核磁孔隙度记为φNMR;测量完成后,将第三段岩样周边和底边缠绕一层生胶带使其与流体隔绝,只保持上端面处于敞开状态,然后将岩样放入自吸瓶中开展自吸水排油实验,每间隔一段时间取出岩样测量其T2谱曲线和重量,T2谱曲线的面积记为Ai,重量记为m3i,同时记录每次测量的时间ti,i=1,2,3…,n-1,每次测量完成后将岩样重新封闭后放入自吸瓶中继续自吸,待自吸量连续24h不变后停止自吸,取出岩样再次测量其T2谱曲线和重量,其T2谱曲线面积记为An,重量记为m3n。
S9、数据处理,具体包括如下步骤:
S91、流量的计算:
根据自吸水排油量实验过程中一定时间间隔内的T2谱曲线面积的变化,可以计算该段时间内的自吸水排油量。首先计算某一时间间隔内的可动流体饱和度ΔSwi(小数),计算公式如下:
式中,φNMR—老化后第三段岩样核磁共振孔隙度,小数;
Vp—孔隙体积,cm3;
Ai,Ai+1—分别是自吸水排油实验中ti和ti+1时刻岩样的T2谱曲线的面积。
然后,计算某一时间间隔内自吸水排油量Voi(cm3)公式如下:
Voi=ΔSwi×Vp。
S92、压力的计算,具体如下:
(1)根据恒速压汞的孔喉分布特征和核磁共振曲线的分布特征,将每条T2谱所对应的弛豫时间转换为相应的孔喉大小,公式如下:
T2i=CirNMRi
式中:T2i—每条T2谱曲线的弛豫时间,ms;
Ci—每条T2谱曲线与孔喉半径间的转换系数,ms/μm;
rNMRi—基于核磁共振转换后的孔喉半径,μm;
(2)拟合测量的油水毛管压力曲线的函数表达式,并利用该函数表达式计算(1)中获得的核磁孔喉半径所对应的毛管压力值Pcowi(MPa),公式如下:
式中,C—基于毛管压力曲线拟合得到的参数。
式中,fi—每级孔喉所占的频率。
S93、相渗透率和相对渗透率计算:
(1)计算流体的渗流面积A(cm2):
式中,D—第三段岩样的直径,cm;
(2)某一时间间隔油相渗透率koi(mD):
式中,Voi—某一时间间隔自吸水排油量,cm3;
μow—油水粘度,mPa·s;
Δti—某一时间间隔,s;
L—第三段岩样长度,cm;
(3)相对渗透率的计算:
式中:kroi—某一时间间隔油相相对渗透率,小数;
k—岩样的渗透率,mD。
S94、含水饱和度的计算:
(1)计算束缚水饱和度Swc(小数):
Swc=100-φNMR×Vp
(2)计算某一时间间隔的含水饱和度Swi(小数):
Swi=Swc+ΔSwi。
S95、根据计算出的含水饱和度和油相相对渗透率绘制油相相对渗透率曲线。
与现有技术相比,本发明的有益之处在于:
由于致密储层岩石表现出较高的毛管压力,这种现象将使得岩心具有较强的自发吸入湿相流体的能力,因此,本发明充分利用致密岩石的该项特性,利用每个时间间隔的T2谱曲线及其变化特征,获取相应时间段内的自吸水排油量,同时,结合每个阶段的T2谱曲线和恒速压汞结果,将T2谱信息有效的转化为孔喉信息,然后,利用油水毛管压力实验结果获取了相应阶段的驱替压力,最终利用达西渗流理论获取了岩心的油相相对渗透率,同时根据自吸增水的原理获取了相应的含水饱和度。该测定方法得到的实验数据计量精度高,操作简单,充分利用的致密岩石的特性及相关渗流理论,具有较高的适用性。
本发明的其它优点、目标和特征将部分通过下面的说明体现,部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
图1、圆柱形岩样切割为三段的示意图。
图2、自吸水排油过程中各个时间段内T2谱曲线的变化。
图3、拟合测量的油水毛管压力曲线。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
一种致密岩石油相相对渗透率测定方法,包括如下步骤:
S1、取长度9cm,直径2.5cm的致密岩石的圆柱形岩样,并将岩样进行洗油、洗盐、烘干等处理。
S2、按图1所示,将清洗烘干的岩样切割为三段,第一段和第二段长度为2.5cm,第三段长度为4cm,然后再将三段岩样放入烘箱中继续烘干1天。
S3、按国家标准GB/T 29171-2012《岩石毛管压力曲线的测定》对第一段岩样进行油水毛管压力实验,测定岩石毛管压力曲线,毛管压力记为Pcow。
S4、采用美国ASPE730型自动孔隙结构测试仪,对第二段岩样进行恒速压汞实验测定孔喉分布特征。
S6、测量完成后,先将第三段岩样抽真空至133Pa,在20MPa压力下加压饱和MnCl2溶液48小时后,MnCl2溶液的浓度大于50000ppm,测量岩样的质量,记为m1,计算岩样的有效孔隙体积Veff和有效孔隙度φeff,计算岩样的有效孔隙体积Veff和有效孔隙度φeff的公式如下:
S7、将饱和完成的第三段岩样放入岩心加持器中,连接驱替装置,并在岩心夹持器出口端接油水分离计量管;在恒压状态下,用实验用油开展油驱水实验,待驱至含油率达到99%时停止驱替,卸载围压和内压,取出岩心放入实验用油中进行老化,老化时间不少于10d。
S8、将老化后的第三段岩样擦去表面油后称重记为m2,同时用核磁共振仪获取岩样的T2谱曲线和核磁孔隙度,其T2谱曲线的面积记为A0,核磁孔隙度记为φNMR;测量完成后,将第三段岩样周边和底边缠绕一层生胶带使其与流体隔绝,只保持上端面处于敞开状态,然后将岩样放入自吸瓶中开展自吸水排油实验,每间隔一段时间取出岩样测量其T2谱曲线和重量,T2谱曲线的面积记为Ai,重量记为m3i,同时记录每次测量的时间ti,i=1,2,3…,n-1,每次测量完成后将岩样重新封闭后放入自吸瓶中继续自吸,待自吸量连续24h不变后停止自吸,取出岩样再次测量其T2谱曲线和重量,其T2谱曲线面积记为An,重量记为m3n。
S9、数据处理,具体包括如下步骤:
S91、流量的计算:
如图2所示,是自吸水排油过程中各个时间段内T2谱曲线的变化。根据自吸水排油量实验过程中一定时间间隔内的T2谱曲线面积的变化,可以计算该段时间内的自吸水排油量。首先计算某一时间间隔内的可动流体饱和度ΔSwi(小数),计算公式如下:
式中,φNMR—老化后第三段岩样核磁共振孔隙度,小数;Vp—孔隙体积,cm3;Ai、Ai+1—分别是自吸水排油实验中ti和ti+1时刻岩样的T2谱曲线的面积。
然后,计算某一时间间隔内自吸水排油量Voi(cm3)公式如下:
Voi=ΔSwi×Vp。
S92、压力的计算,具体如下:
(1)根据恒速压汞的孔喉分布特征和核磁共振曲线(图2)的分布特征,将每条T2谱所对应的弛豫时间转换为相应的孔喉大小,公式如下:
T2i=CirNMRi
式中:T2i—每条T2谱曲线的弛豫时间,ms;Ci—每条T2谱曲线与孔喉半径间的转换系数,ms/μm;rNMRi—基于核磁共振转换后的孔喉半径,μm。
(2)拟合测量的油水毛管压力曲线的函数表达式(图3),并利用该函数表达式计算(1)中获得的核磁孔喉半径所对应的毛管压力值Pcowi(MPa),公式如下:
式中,C—基于毛管压力曲线拟合得到的参数。
式中,fi—每级孔喉所占的频率。
S93、相渗透率和相对渗透率计算:
(1)计算流体的渗流面积A(cm2):
式中,D—第三段岩样的直径,cm;
(2)某一时间间隔油相渗透率koi(mD):
式中,Voi—某一时间间隔自吸水排油量,cm3;μow—油水粘度,mPa·s;Δti—某一时间间隔,s;L—第三段岩样长度,cm。
(3)相对渗透率的计算:
式中:kroi—某一时间间隔油相相对渗透率,小数;k—岩样的渗透率,mD。
S94、含水饱和度的计算:
(1)计算束缚水饱和度Swc(小数):
Swc=100-φNMR×Vp。
(2)计算某一时间间隔的含水饱和度Swi(小数):
Swi=Swc+ΔSwi。
S95、根据上述公式计算出的含水饱和度和油相相对渗透率绘制油相相对渗透率曲线。
综上所述,本发明基于致密岩石的岩心具有较强的自发吸入湿相流体的能力的特点,提供了一种适用于致密岩石油相相对渗透率测定方法,该方法得到的实验数据计量精度高,操作简单,充分利用的致密岩石的特性及相关渗流理论,具有较高的适用性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (9)
1.一种致密岩石油相相对渗透率测定方法,其特征在于,包括如下步骤:
S1、取致密岩石的长度大于7cm的圆柱形岩样,并进行清洗烘干处理;
S2、将岩样切割为三段,将三段岩样再次烘干;
S3、对第一段岩样进行油水毛管压力实验,测定岩石毛管压力曲线,毛管压力记为Pcow;
S4、对第二段岩样进行恒速压汞实验测定孔喉分布特征;
S5、测量第三段岩样的孔隙度φHe、渗透率k、长度L、直径D、干重m0和密度ρ,并计算孔隙体积Vp;
S6、测量完成后,将第三段岩样抽真空加压饱和MnCl2溶液,饱和完成后测得第三段岩样质量m1;
S7、将饱和完成的第三段岩样放入岩心加持器中,连接驱替装置,并在岩心夹持器出口端接油水分离计量管;在恒压状态下,用实验用油开展油驱水实验,待驱至含油率达到99%时停止驱替,卸载围压和内压,取出岩心放入实验用油中进行老化;
S8、将老化后的第三段岩样擦去表面油后称重记为m2,同时用核磁共振仪获取岩样的T2谱曲线和核磁孔隙度,其T2谱曲线的面积记为A0,核磁孔隙度记为φNMR;测量完成后,将第三段岩样周边和底边缠绕一层生胶带使其与流体隔绝,只保持上端面处于敞开状态,然后将岩样放入自吸瓶中开展自吸水排油实验,每间隔一段时间取出岩样测量其T2谱曲线和重量,T2谱曲线的面积记为Ai,重量记为m3i,同时记录每次测量的时间ti,i=1,2,3…,n-1,每次测量完成后将岩样重新封闭后放入自吸瓶中继续自吸,待自吸量连续24h不变后停止自吸,取出岩样再次测量其T2谱曲线和重量,其T2谱曲线面积记为An,重量记为m3n;
S9、数据处理,具体包括如下步骤:
S91、流量的计算:某一时间间隔内的可动流体饱和度ΔSwi和自吸水排油量Voi计算公式分别如下:
Voi=ΔSwi×Vp
式中,φNMR—老化后第三段岩样核磁共振孔隙度;
Vp—孔隙体积,cm3;
Ai,Ai+1—分别是自吸水排油实验中ti和ti+1时刻岩样的T2谱曲线的面积;
S92、压力的计算,具体如下:
(1)根据恒速压汞的孔喉分布特征和核磁共振曲线的分布特征,将每条T2谱所对应的弛豫时间转换为相应的孔喉大小,公式如下:
T2i=CirNMRi
式中:T2i—每条T2谱曲线的弛豫时间,ms;
Ci—每条T2谱曲线与孔喉半径间的转换系数,ms/μm;
rNMRi—基于核磁共振转换后的孔喉半径,μm;
(2)拟合测量的油水毛管压力曲线的函数表达式,并利用该函数表达式计算(1)中获得的核磁孔喉半径所对应的毛管压力值Pcowi,公式如下:
式中,C—基于毛管压力曲线拟合得到的参数;
式中,fi—每级孔喉所占的频率;
S93、油相渗透率和相对渗透率计算:
(1)计算流体的渗流面积A:
式中,D第三段岩样的直径;
(2)某一时间间隔油相渗透率koi:
式中,Voi—某一时间间隔自吸水排油量,cm3;
μow—油水粘度,mPa·s;
Δti—某一时间间隔,s;
L—第三段岩样长度,cm;
(3)相对渗透率的计算:
式中:kroi—某一时间间隔相对渗透率,小数;
k—岩样的渗透率,mD;
S94、含水饱和度的计算:
(1)计算束缚水饱和度Swc:
Swc=100-φNMR×Vp
(2)计算某一时间间隔的含水饱和度Swi:
Swi=Swc+ΔSwi
S95、绘制油相相对渗透率曲线。
5.如权利要求3所述的致密岩石油相相对渗透率测定方法,其特征在于,所述MnCl2溶液的浓度大于50000ppm。
6.如权利要求1所述的致密岩石油相相对渗透率测定方法,其特征在于,所述步骤S1中,圆柱形岩样长度7-11cm,直径2-3cm。
7.如权利要求1所述的致密岩石油相相对渗透率测定方法,其特征在于,所述步骤S2中,第一段和第二段长度为2-3cm,第三段长度为3-5cm。
8.如权利要求1所述的致密岩石油相相对渗透率测定方法,其特征在于,按国家标准GB/T 29171-2012《岩石毛管压力曲线的测定》对第一段岩样进行油水毛管压力实验。
9.如权利要求1所述的致密岩石油相相对渗透率测定方法,其特征在于,步骤S7中,老化时间不少于10d。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910647504.5A CN110346258B (zh) | 2019-07-17 | 2019-07-17 | 一种致密岩石油相相对渗透率测定的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910647504.5A CN110346258B (zh) | 2019-07-17 | 2019-07-17 | 一种致密岩石油相相对渗透率测定的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110346258A CN110346258A (zh) | 2019-10-18 |
CN110346258B true CN110346258B (zh) | 2020-06-02 |
Family
ID=68176654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910647504.5A Active CN110346258B (zh) | 2019-07-17 | 2019-07-17 | 一种致密岩石油相相对渗透率测定的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110346258B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112782477B (zh) * | 2019-11-11 | 2024-05-14 | 中国石油化工股份有限公司 | 一种测定岩心不同润湿状态下电性响应特征的方法及系统 |
CN111307852B (zh) * | 2020-03-26 | 2023-11-28 | 中国石油天然气集团有限公司 | 一种混积岩的核磁及岩电实验工艺 |
CN111879678B (zh) * | 2020-07-31 | 2022-06-07 | 西南石油大学 | 一种基于自吸法的致密砂岩气水相对渗透率实验方法 |
CN114065462B (zh) * | 2020-08-05 | 2022-11-04 | 中国石油天然气股份有限公司 | 岩石孔隙结构演化分析方法及系统 |
CN112798494B (zh) * | 2021-02-01 | 2024-07-30 | 西南石油大学 | 一种长岩心干化渗流实验装置及实验方法 |
CN113075102B (zh) * | 2021-03-11 | 2022-02-01 | 中国地质大学(北京) | 一种建立多孔介质自发渗吸量与时间关系数学模型的方法 |
CN113063713B (zh) * | 2021-03-23 | 2022-04-05 | 西南石油大学 | 一种大直径长岩心渗流截面上非均匀压力分布测试方法 |
CN113075108B (zh) * | 2021-03-26 | 2022-02-08 | 西南石油大学 | 一种考虑束缚水饱和度的岩心多次应力敏感测试方法 |
CN113295592B (zh) * | 2021-05-20 | 2022-06-07 | 中国石油大学(北京) | 一种低渗透岩石的相对渗透率测定系统及方法 |
CN113295594B (zh) * | 2021-05-20 | 2022-07-12 | 中国石油大学(北京) | 一种基于分流模型与ct扫描的相对渗透率测定系统及方法 |
CN113433157B (zh) * | 2021-06-24 | 2023-12-15 | 西南石油大学 | 基于核磁共振t2谱建立随机单元等效岩心模型的方法 |
CN113640890B (zh) * | 2021-08-31 | 2022-12-02 | 成都理工大学 | 基于核磁共振测井的储层气-水相对渗透率曲线构建方法、系统 |
CN113916748B (zh) * | 2021-10-29 | 2024-02-02 | 东北石油大学 | 轻质油测页岩基质渗透率和采收率的装置和方法 |
CN115078163A (zh) * | 2022-08-05 | 2022-09-20 | 西南石油大学 | 一种页岩油储层岩石润湿性表征方法 |
CN117269000B (zh) * | 2023-11-23 | 2024-02-09 | 中国石油大学(华东) | 一种测定致密岩心两相相对渗透率的方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267721B (zh) * | 2013-05-03 | 2015-07-08 | 中国石油天然气集团公司 | 一种致密砂岩储层孔隙含水特征与赋存状态的评价方法 |
CN105353419A (zh) * | 2015-11-25 | 2016-02-24 | 中国地质大学(北京) | 一种基于储集层分类的核磁毛管压力曲线构造方法 |
US11774622B2 (en) * | 2017-11-06 | 2023-10-03 | Schlumberger Technology Corporation | System and method for reducing a feasible region of solutions in a relative permeability and capillary pressure curve |
CN108956678B (zh) * | 2018-06-11 | 2020-12-04 | 西南石油大学 | 一种基于核磁共振测井的t2谱敏感参数提取方法 |
CN109916799B (zh) * | 2019-03-22 | 2021-04-20 | 西南石油大学 | 测量非常规致密气藏自发渗吸相对渗透率的实验方法 |
-
2019
- 2019-07-17 CN CN201910647504.5A patent/CN110346258B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110346258A (zh) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110346258B (zh) | 一种致密岩石油相相对渗透率测定的方法 | |
CN109443867B (zh) | 一种对致密岩石的物性参数进行连续检测的方法 | |
US10161891B1 (en) | Method for characterizing rock physical characteristics of deeply buried carbonate rocks | |
CN109030292B (zh) | 一种致密岩石润湿性确定的新方法 | |
CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及系统 | |
CN107991335B (zh) | 致密砂岩水锁损害评价试验方法 | |
RU2525093C1 (ru) | Способ прогнозирования изменения свойств призабойной зоны пласта под воздействием бурового раствора | |
CN108827853B (zh) | 基于核磁共振的致密储层岩电测量装置及测量方法 | |
RU2455482C2 (ru) | Способ определения профиля притока флюидов и параметров околоскважинного пространства | |
CN110927035A (zh) | 一种低渗致密砂岩束缚水饱和度计算方法 | |
CN104697915A (zh) | 一种页岩微观孔隙大小及流体分布的分析方法 | |
CN103233726B (zh) | 一种水淹层饱和度、产水率测井评价模型的实验刻度方法 | |
CN109580689B (zh) | 一种核磁共振测井t2截止值的逐点计算方法 | |
CN108444890B (zh) | 测试中、高渗岩心液测渗透率的非稳态滴定装置及方法 | |
CN113431537B (zh) | 一种非稳态变流速大尺度岩心水驱气相对渗透率测试方法 | |
CN106501151A (zh) | 一种基于渗吸和离子扩散特性的页岩孔径测量装置及方法 | |
CN115078163A (zh) | 一种页岩油储层岩石润湿性表征方法 | |
CN208155804U (zh) | 基于核磁共振的致密储层岩电测量装置 | |
CN115559715A (zh) | 一种超高压低渗透气藏储层产水评价方法 | |
CN115078210A (zh) | 页岩孔隙结构测试方法 | |
CN105401937A (zh) | 一种基于孔隙结构的饱和度指数预测方法 | |
CN110672495A (zh) | 一种基于低场磁共振技术的水泥基材料水分渗透率的预测方法 | |
CN110905493A (zh) | 一种测量井下地层流体污染率的方法 | |
CN114370269A (zh) | 深层碳酸盐岩气藏有效储层物性下限综合确定方法 | |
CN111879678B (zh) | 一种基于自吸法的致密砂岩气水相对渗透率实验方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |