CN112362553B - 一种致密砂岩微观孔隙结构表征方法 - Google Patents
一种致密砂岩微观孔隙结构表征方法 Download PDFInfo
- Publication number
- CN112362553B CN112362553B CN202011228812.3A CN202011228812A CN112362553B CN 112362553 B CN112362553 B CN 112362553B CN 202011228812 A CN202011228812 A CN 202011228812A CN 112362553 B CN112362553 B CN 112362553B
- Authority
- CN
- China
- Prior art keywords
- pore
- distribution
- radius
- throat
- mercury
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 189
- 238000012512 characterization method Methods 0.000 title claims abstract description 14
- 238000009826 distribution Methods 0.000 claims abstract description 108
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims abstract description 88
- 229910052753 mercury Inorganic materials 0.000 claims abstract description 88
- 238000003825 pressing Methods 0.000 claims abstract description 43
- 239000011435 rock Substances 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000005266 casting Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000002347 injection Methods 0.000 claims description 24
- 239000007924 injection Substances 0.000 claims description 24
- 238000005481 NMR spectroscopy Methods 0.000 claims description 22
- 238000002474 experimental method Methods 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 12
- 229920006395 saturated elastomer Polymers 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 230000001788 irregular Effects 0.000 claims description 5
- 230000035699 permeability Effects 0.000 claims description 5
- 238000005520 cutting process Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000005406 washing Methods 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 244000025254 Cannabis sativa Species 0.000 description 3
- 235000012766 Cannabis sativa ssp. sativa var. sativa Nutrition 0.000 description 3
- 235000012765 Cannabis sativa ssp. sativa var. spontanea Nutrition 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 3
- 235000009120 camo Nutrition 0.000 description 3
- 235000005607 chanvre indien Nutrition 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000002500 effect on skin Effects 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 239000011487 hemp Substances 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005587 bubbling Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000009736 wetting 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/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
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/44—Sample treatment involving radiation, e.g. heat
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
- G01N2001/2873—Cutting or cleaving
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种致密砂岩微观孔隙结构表征方法,首先,将铸体薄片划分为不同的孔隙单元,提取每一单元中孔隙的图片并统计每个孔隙的面积及周长,计算对应孔隙的形状因子,并基于不同形状的孔隙半径计算公式计算出对应的孔隙半径。然后,结合铸体薄片中的孔隙分布和恒速压汞进汞曲线的特征,同时,将高压压汞获取的喉道分布与修正后恒速压汞喉道分布相结合获取了岩心完整的喉道分布。最后,基于铸体薄片获取的形状因子分布分别提出了孔隙形状分布稳定及孔隙形状分布差异较大时岩心的孔隙半径及其分布表征方法。本发明摒弃了现有的各种微观孔隙结构表征方法的不足,获取的致密岩石完整的孔隙和喉道大小及其分布的数据更加精确。
Description
技术领域
本发明涉及非常规油气勘探技术领域,具体涉及一种致密砂岩微观孔隙结构表征方法。
背景技术
致密油是继页岩之后全球非常规油气勘探开发领域的又一大热点,已然成为我国未来重要的能源接替之一。在致密储层勘探与开发过程中,与中、高渗储层存在明显区别的是致密储层岩石颗粒细小、孔隙度和渗透率低,成岩过程中受压实作用较强,微观孔隙结构复杂,孔喉多以微纳米级为主。准确表征致密储层岩石的微观孔喉结构,弄清致密储层孔隙和喉道大小及其分布对致密储层高效开发具有十分重要的现实意义。
目前,常用的储层孔喉结构表征方法有定性和定量两个方面,针对致密储层,常用的代表性技术手段有铸体薄片、恒压压汞、恒速压汞和核磁共振。所有这些表征方法中,铸体薄片可以直接用于观察致密油储层孔喉结构,在提供直观认识的同时经过后续图像处理也可提取岩心中的相关孔喉结构信息。恒压压汞通过恒定的压力向样品中注入汞并获取毛管压力曲线,再结合Washburn方程计算获取孔喉半径分布;恒速压汞以恒定的、极低的速度注入岩心,进汞过程逼近于准静态过程,通过分析进汞过程中压力曲线的变化,可同时获得孔隙和喉道的变化特征,进而分别得到孔隙、喉道的毛细管曲线,再结合Washburn方程可以计算获取孔隙半径分布和喉道半径分布。核磁共振技术是利用氢原子核自身的磁性及其与外加磁场相互作用的原理,通过测量地层岩石孔隙流体中氢核的核磁共振弛豫信号的幅度和弛豫速率来探测地层孔隙结构的一种技术,在探测岩石孔隙中饱和含氢流体时,氢原子核弛豫信号的幅度与地层的孔隙度成正比,其弛豫速率或弛豫时间T2与孔隙大小和流体性质有关。
但是,目前这些测试方法存在以下不足:
(1)铸体薄片虽能直观的观测孔隙的大小及其连通情况,但反映的信息是二维信息,和三维空间中岩心的孔隙大小及其分布存在一定的误差。
(2)在恒速压汞实验的最初进汞阶段,进汞量的增加是由于非润湿相的汞在岩样粗糙表面的坑凹处的贴和引起的虚假进汞体积。随着压力的逐渐增大,坑凹被汞占满,此时汞还并没有真正进入孔喉***,压力也没有达到排驱压力。但在进汞量中,如把这一部分的空腔体积累计到总孔喉***的进汞量中,会造成进汞饱和度数值偏大,这一现象称为麻皮效应。受麻皮效应存在的影响,使汞饱和度偏大,所绘制的毛管曲线存在失真现象,所表征的孔隙大小及其分布存在较大的误差;此外,恒速压汞最大注汞压力为1000Psi(6.9MPa),反映的孔喉半径的下限为0.12μm,对于0.12μm以下的孔喉无法进行表征。
(3)高压压汞给出的孔喉信息虽主要反映的是喉道分布,但其中仍叠加有部分半径较小的孔隙分布,因此获取的喉道信息并不完整。
(4)此外,在常规恒速压汞和高压压汞数据处理中,将孔隙等效为球体并利用Washburn公式计算出相应的等效球体半径,用以表征岩石孔隙大小及分布,但是这对于经历了强烈压实和胶结作用、孔隙形状不规则的致密储层砂岩是不适用的,将造成计算值偏大。
(5)核磁共振技术可以获取完整的孔喉信息,但是不能直接获取孔喉半径分布。同时,将核磁共振T2弛豫时间转化为孔喉分布时,普遍将岩心孔隙形状视为定值,未考虑岩心孔隙形状差异较大时其孔隙半径的表征方法。
发明内容
针对上述问题,本发明提供一种致密砂岩微观孔隙结构表征方法,摒弃了现有的各种微观孔隙结构表征方法的不足,获取的致密砂岩储层全孔隙和喉道大小及其分布结果更加精确。
本发明采用下述的技术方案:
一种致密砂岩微观孔隙结构表征方法,包括以下步骤:
S1、选取一块致密砂岩储层的圆柱形岩心,并进行洗油和烘干,测量岩心的孔隙度和渗透率后,将岩心切割成四段;
S2、对第Ⅰ段岩心进行铸体薄片实验:将铸体薄片图像划分成不同孔隙单元,计算每个孔隙的形状因子,基于形状因子的大小判断孔隙形状,并将孔隙形状按形状分类,根据孔隙形状及不同孔隙形状的面积与其内切圆半径的关系,计算出对应孔隙形状的内切圆半径和等效半径,并绘制半径及其频率分布直方图;
S3、对第Ⅱ段岩心进行恒速压汞实验,获取修正后的恒速压汞喉道分布;
S4、对第Ⅲ段岩心进行高压压汞测试,基于高压压汞实验获取岩心完整的喉道半径及其分布;
S5、获取岩心完整的孔隙半径及其分布。
优选的,步骤S2中,所述形状因子的计算公式为:
式中,G为形状因子,无量纲;A为孔隙的截面积,μm2;P为孔隙的周长,μm。
优选的,步骤S2中,所述对应孔隙形状的内切圆半径的计算公式为:
不规则多边形:R等效=0.5L (5)
式中,A为空隙的截面积,μm2,L为孔隙长度,μm。
优选的,步骤S3中,所述获取修正后的恒速压汞喉道分布包括以下步骤:
S31、对第Ⅱ段岩心进行恒速压汞实验,获取注汞压力与注汞体积曲线;利用步骤S2中基于铸体薄片获取的孔隙半径及其分布计算出极值压力区间,并从恒速压汞注汞压力和注汞体积曲线中选出对应数据段;
S32、将曲线处理成“M”形曲线段,从曲线段中识别出孔隙和喉道对应的压力值和体积值,并绘制出修正后的孔隙毛管压力曲线、喉道毛管压力曲线和总毛管压力曲线;
S33、计算得到修正后的毛管压力曲线、孔隙、喉道半径及其分布。
优选的,所述极值压力区间的计算公式为:
所述修正后的毛管压力曲线、孔隙、喉道半径的计算公式为:
式中,Rmax为基于铸体薄片获取的孔隙半径的最大值,μm;Plb为恒速压汞曲线的注汞压力的下限值,MPa;Rmin为基于铸体薄片获取的孔隙半径的最小值,μm;Pup为恒速压汞曲线上注汞压力的上限值,MPa;Rcp为恒速压汞修正后的孔隙半径,μm;Pcp为修正后的孔隙毛管压力,MPa;Rct为恒速压汞修正后的喉道半径,μm;Pct为喉道毛管压力,MPa。
优选的,步骤S4包括以下步骤:
S41、对第Ⅲ段岩样进行高压压汞测试,获取基于高压压汞实验的孔喉半径分布;
S42、将修正后的恒速压汞喉道半径及其分布与高压压汞获取的孔喉半径及其分布绘制在同一坐标系下;
将0.12μm之前的数据选高压压汞孔喉半径及其分布,0.12μm之后的数据选择修正后的恒速压汞喉道半径及其分布;
将两个数据归一化处理后重新绘制在直角坐标系下得到岩样完整的喉道半径及其分布。
优选的,步骤S5中,所述获取岩心完整的孔隙半径及其分布包括以下步骤:
S51、对第Ⅳ段岩心饱和水,测量完全饱和水状态下的核磁共振T2谱曲线;
S52、基于步骤S2获取的孔隙形状因子及其分布判断岩心的孔隙形状分布是否稳定。
本发明的有益效果是:
1、本发明首先将铸体薄片划分为多个孔隙单元,通过获取每一个孔隙的面积、周长和长度,通过孔隙的形状因子将孔隙划分为规则的三边形、四边形、五边形和不规则多边形,根据孔隙形状的差异推导了对应孔隙的内切圆半径和等效半径的计算表达式,获取了基于铸薄片完整的孔隙大小和分别。
2、其次,将铸体薄片获取的孔隙半径分布与恒速压汞相结合,选取了恒速压汞合理的曲线区间,消除了由于麻皮效应所造成的误差,获取了修正后的孔隙和喉道大小及其分布。同时,将高压压汞获取的喉道分布与修正后恒速压汞喉道分布相结合获取了岩心完整的喉道分布。
3、最后,基于铸体薄片获取的形状因子分布提出了孔隙形状分布稳定及孔隙形状分布差异较大时岩心的孔隙半径及其分布的表征方法。本发明摒弃了现有的各种微观孔隙结构表征方法的不足,获取的致密岩石完整的孔隙和喉道大小及其分布的数据更加精确。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例的附图作简单地介绍,显而易见地,下面描述中的附图仅仅涉及本发明的一些实施例,而非对本发明的限制。
图1为本发明完全饱和水时的T2谱曲线示意图;
图2为本发明基于铸体薄片统计的孔隙形状分布示意图;
图3为本发明基于铸体薄片获取的孔隙大小及其分布示意图;
图4为本发明修正前后恒速压汞孔隙大小及其分布对比图示意图;
图5为本发明修正前后恒速压汞喉道大小及其分布示意图;
图6为本发明高压压汞喉道大小及其分布示意图;
图7为本发明完整的喉道分布示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例的附图,对本发明实施例的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于所描述的本发明的实施例,本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
除非另外定义,本公开使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本公开中使用的“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同,而不排除其他元件或者物件。“上”、“下”、“左”、“右”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。
下面结合附图和实施例对本发明进一步说明。
如图1至图7所示,一种致密砂岩微观孔隙结构表征方法,包括以下步骤:
S1、选取一块致密砂岩储层的圆柱形岩心,并进行洗油和烘干,测量岩心的孔隙度和渗透率后,将岩心切割成四段;
S2、对第Ⅰ段岩心进行铸体薄片实验:将铸体薄片图像划分成不同孔隙单元,使用高倍显微镜拍摄每一单元中所有孔隙的图片;利用ImageJ软件提取相应的孔隙轮廓、孔隙大小(孔隙截面积)A、孔隙周长P和孔隙长度L;利用以下公式计算每个孔隙的形状因子:
式中,G为形状因子,无量纲;A为孔隙的截面积,μm2;P为孔隙的周长,μm。
如表1、表2所示,基于形状因子的大小判断孔隙形状,并将孔隙形状按形状分类,根据孔隙形状及不同孔隙形状的面积与其内切圆半径的关系,计算出对应孔隙形状的内切圆半径和等效半径,并绘制半径及其频率分布直方图;
形状类型 | 所占比例 |
三边形 | 26% |
四边形 | 63% |
五边形 | 3% |
圆形 | 8% |
表1孔隙形状分类统计表
表2内切圆半径分布统计表
所述对应孔隙形状的内切圆半径的计算公式为:
不规则多边形:R等效=0.5L (5)
式中,A为空隙的截面积,μm2,L为孔隙长度,μm。
S3、对第Ⅱ段岩心进行恒速压汞实验,获取修正后的恒速压汞喉道分布;
S31、对第Ⅱ段岩心进行恒速压汞实验,记录实验过程中的注汞压力与注汞体积,利用步骤S2中基于铸体薄片获取的孔隙半径及其分布计算出极值压力区间,并从恒速压汞注汞压力和注汞体积曲线中选出对应数据段;
式中,Rmax为基于铸体薄片获取的孔隙半径的最大值,μm;Plb为恒速压汞曲线的注汞压力的下限值,MPa;Rmin为基于铸体薄片获取的孔隙半径的最小值,μm;Pup为恒速压汞曲线上注汞压力的上限值,MPa;
S32、根据恒速压汞实验原理,剔除曲线中的异常点,将曲线处理成“一升一降”的“M”形曲线段,并利用冒泡排序的算法从“M”形曲线段中识别出孔隙和喉道对应的压力值和体积值,并绘制出修正后的孔隙毛管压力曲线、喉道毛管压力曲线和总毛管压力曲线。
S33、计算以下公式计算得到修正后的毛管压力曲线、孔隙、喉道半径及其分布:
Rcp为恒速压汞修正后的孔隙半径,μm;Pcp为修正后的孔隙毛管压力,MPa;Rct为恒速压汞修正后的喉道半径,μm;Pct为喉道毛管压力,MPa。
S4、对第Ⅲ段岩心进行高压压汞测试,基于高压压汞实验获取岩心完整的喉道半径及其分布;
S41、按国家标准GB/T 29171-2012《岩石毛管压力曲线的测定》对第Ⅲ段岩样进行高压压汞测试,并获取基于高压压汞实验获取的孔喉半径分布;
S42、将修正后的恒速压汞喉道半径及其分布与高压压汞获取的孔喉半径及其分布绘制在同一坐标系下;以0.12μm为界,0.12μm之前的数据选高压压汞孔喉半径及其分布,0.12μm之后的数据选择修正后的恒速压汞喉道半径及其分布;最后,将两个数据归一化处理后重新绘制在直角坐标系下得到岩样完整的喉道半径及其分布。
S5、获取岩心完整的孔隙半径及其分布。
S51、测量第Ⅳ段岩心的干重,将岩心抽真空至133Pa,在20MPa的压力下加压饱和模拟地层水48小时后,测量饱和岩心的重量,计算岩心的有效孔隙度,当岩样的有效孔隙度与气测孔隙度的相对误差小于2%时,即满足公式8时,岩样饱和完成,否则重复上述步骤重新饱和直至岩心完全饱和;测量完全饱和水状态下的核磁共振T2谱曲线。
S5.2基于S2获取的孔隙形状因子及其分布判断岩心的孔隙形状分布是否稳定(即岩心中的孔隙形状主要以某一种形状为主);
S521、当岩心中某一种形状的孔隙分布数量占比在60%以上时,认定为稳定:将步骤S3中获取的修正后恒速压汞孔隙大小及分布和核磁共振T2谱曲线相结合,将孔隙最大值与核磁共振T2时间的最大值相对应,并以修正后的恒速压汞孔隙分布和核磁共振T2谱曲线形状为参考,利用公式9获取核磁共振T2时间与孔喉半径间的转换系数,从而将核磁共振T2谱转换为对应的孔喉分布曲线。
T2=CRNMR (9)
式中,C为转换系数,ms/μm;T2为核磁共振T2时间,ms;RNMR为核磁共振转换后孔喉半径,μm;
根据转换后的核磁共振孔喉分布曲线,绘制孔喉累计分布曲线;同时,将步骤S4中获取的完整喉道分布绘制成累计分布曲线,将两条累计分布曲线绘制在同一坐标系下,并将两条累计分布曲线相减并重新绘制后即获得岩心的完整的孔隙大小及分布。
当岩心中的孔隙的形状分布差异较大时(即,没有一种形状的孔隙数量占比超过60%):
第一、将步骤S4中获取的完整喉道分布和核磁共振T2谱曲线绘制在同一坐标系下;根据公式10中T2时间与孔喉半径间的函数关系,利用核磁共振T2谱曲线左锋标定获取弛豫率(其中的形状因子为2,喉道形状主要以长圆柱形为主),并将核磁共振T2谱曲线转化为孔喉分布曲线;
式中,ρ2为表面弛豫率,μm/ms;Fs为孔隙形状因子(长圆柱形为2、球形为3);RNMR为核磁共振转换后孔喉半径,μm;
第二、根据转换后的核磁共振孔喉分布曲线,同样绘制累计孔喉分布曲线和累计喉道分布曲线,并将两条累计分布曲线相减并重新采样后获取长圆柱形的孔隙大小及其分布;
第三、根据长圆柱形和球形形状因子的比值(2/3),将其长圆柱形孔隙大小及其分布转化为等效球体孔隙大小及其分布;
第四、根据不同截面形状下内切球与等效球体间的关系(公式11-13)分别获取不同形状下内切球半径及其分布,并基于步骤S2获取的孔隙形状因子及其分布,利用加权平均的方法获取岩心平均孔隙半径REB及其分布。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (7)
1.一种致密砂岩微观孔隙结构表征方法,其特征在于,包括以下步骤:
S1、选取一块致密砂岩储层的圆柱形岩心,并进行洗油和烘干,测量岩心的孔隙度和渗透率后,将岩心切割成四段;
S2、对第Ⅰ段岩心进行铸体薄片实验:将铸体薄片图像划分成不同孔隙单元,计算每个孔隙的形状因子,基于形状因子的大小判断孔隙形状,并将孔隙形状按形状分类,根据孔隙形状及不同孔隙形状的面积与其内切圆半径的关系,计算出对应孔隙形状的内切圆半径和等效半径,并绘制半径及其频率分布直方图;
S3、对第Ⅱ段岩心进行恒速压汞实验,获取修正后的恒速压汞喉道分布;
S4、对第Ⅲ段岩心进行高压压汞测试,基于高压压汞实验获取岩心完整的喉道半径及其分布;
S5、获取岩心完整的孔隙半径及其分布。
4.根据权利要求1所述的一种致密砂岩微观孔隙结构表征方法,其特征在于,步骤S3中,所述获取修正后的恒速压汞喉道分布包括以下步骤:
S31、对第Ⅱ段岩心进行恒速压汞实验,获取注汞压力与注汞体积曲线;利用步骤S2中基于铸体薄片获取的孔隙半径及其分布计算出极值压力区间,并从恒速压汞注汞压力和注汞体积曲线中选出对应数据段;
S32、将曲线处理成“M”形曲线段,从曲线段中识别出孔隙和喉道对应的压力值和体积值,并绘制出修正后的孔隙毛管压力曲线、喉道毛管压力曲线和总毛管压力曲线;
S33、计算得到修正后的毛管压力曲线、孔隙、喉道半径及其分布。
6.根据权利要求1所述的一种致密砂岩微观孔隙结构表征方法,其特征在于,步骤S4包括以下步骤:
S41、对第Ⅲ段岩样进行高压压汞测试,获取基于高压压汞实验的孔喉半径分布;
S42、将修正后的恒速压汞喉道半径及其分布与高压压汞获取的孔喉半径及其分布绘制在同一坐标系下;
将0.12μm之前的数据选高压压汞孔喉半径及其分布,0.12μm之后的数据选择修正后的恒速压汞喉道半径及其分布;
将两个数据归一化处理后重新绘制在直角坐标系下得到岩样完整的喉道半径及其分布。
7.根据权利要求1所述的一种致密砂岩微观孔隙结构表征方法,其特征在于,步骤S5中,所述获取岩心完整的孔隙半径及其分布包括以下步骤:
S51、对第Ⅳ段岩心饱和水,测量完全饱和水状态下的核磁共振T2谱曲线;
S52、基于步骤S2获取的孔隙形状因子及其分布判断岩心的孔隙形状分布是否稳定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011228812.3A CN112362553B (zh) | 2020-11-06 | 2020-11-06 | 一种致密砂岩微观孔隙结构表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011228812.3A CN112362553B (zh) | 2020-11-06 | 2020-11-06 | 一种致密砂岩微观孔隙结构表征方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112362553A CN112362553A (zh) | 2021-02-12 |
CN112362553B true CN112362553B (zh) | 2021-06-15 |
Family
ID=74509483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011228812.3A Active CN112362553B (zh) | 2020-11-06 | 2020-11-06 | 一种致密砂岩微观孔隙结构表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112362553B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113269712B (zh) * | 2021-04-07 | 2022-08-23 | 长江大学 | 一种基于铸体薄片构建储层岩石毛管压力曲线的方法 |
CN113239955A (zh) * | 2021-04-07 | 2021-08-10 | 长江大学 | 一种碳酸盐岩储层岩石分类方法 |
CN114965542A (zh) * | 2022-05-25 | 2022-08-30 | 西南石油大学 | 一种页岩表面弛豫率的确定方法 |
CN117345216B (zh) * | 2023-12-05 | 2024-03-15 | 西南石油大学 | 一种水侵气藏气井井周水体可动临界孔喉半径的确定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2260820A (en) * | 1991-09-20 | 1993-04-28 | Exxon Research Engineering Co | Permeability determination from NMR T2 measurements for fluids in porous media |
CN102540265A (zh) * | 2011-12-21 | 2012-07-04 | 西南石油大学 | 一种基于网络模拟的多孔介质含水饱和度计算方法 |
CN103645125A (zh) * | 2013-10-28 | 2014-03-19 | 北京大学 | 一种致密油储层渗流能力的评价方法和*** |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102565858B (zh) * | 2011-12-21 | 2014-04-16 | 西南石油大学 | 一种多孔介质含水饱和度的计算方法 |
CN103679814A (zh) * | 2012-08-31 | 2014-03-26 | 西门子公司 | 一种三角形网格模型的处理方法及装置 |
CN103278436B (zh) * | 2013-02-01 | 2015-02-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
CN104573198B (zh) * | 2014-12-23 | 2017-08-15 | 长江大学 | 基于随机分形理论的数字岩心及孔隙网络模型重构方法 |
CN105424576A (zh) * | 2015-11-12 | 2016-03-23 | 中国科学院地质与地球物理研究所 | 一种页岩孔隙分布均匀性评价方法 |
-
2020
- 2020-11-06 CN CN202011228812.3A patent/CN112362553B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2260820A (en) * | 1991-09-20 | 1993-04-28 | Exxon Research Engineering Co | Permeability determination from NMR T2 measurements for fluids in porous media |
CN102540265A (zh) * | 2011-12-21 | 2012-07-04 | 西南石油大学 | 一种基于网络模拟的多孔介质含水饱和度计算方法 |
CN103645125A (zh) * | 2013-10-28 | 2014-03-19 | 北京大学 | 一种致密油储层渗流能力的评价方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN112362553A (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112362553B (zh) | 一种致密砂岩微观孔隙结构表征方法 | |
CN105781539B (zh) | 一种致密油气储层饱和度测井计算方法 | |
CN107817199B (zh) | 一种致密砂岩多尺度孔隙模型的构建方法和应用 | |
CN103713316B (zh) | 一种基于岩石孔隙数字表征的速度预测方法及装置 | |
CN103336305B (zh) | 一种基于灰色理论划分致密砂岩储层岩石物理相的方法 | |
CN105447762B (zh) | 一种流体替换的低渗透油藏水淹信息的计算方法 | |
CN111521539B (zh) | 一种致密砂砾岩储层中不同类型孔隙结构定量表征方法 | |
CN109840497A (zh) | 一种基于深度学习的指针式水表读数检测方法 | |
WO2019205725A1 (zh) | 一种基于电成像测井分形维数的火山岩岩性识别方法 | |
CN116580021A (zh) | 一种环保混凝土降碳制品生产及质量检测方法 | |
CN108061697B (zh) | 土体三维孔隙率计算方法 | |
CN113158490A (zh) | 致密砂岩储层渗透率计算模型的建立方法 | |
CN106960433B (zh) | 一种基于图像熵以及边缘的全参考声纳图像质量评测方法 | |
CN114370269B (zh) | 深层碳酸盐岩气藏有效储层物性下限综合确定方法 | |
CN111860176A (zh) | 一种非金属夹杂全视场定量统计分布表征方法 | |
CN110208874B (zh) | 一种致密砂岩油藏有效储层识别方法 | |
CN117274263B (zh) | 一种显示器瘢痕缺陷检测方法 | |
CN110593857A (zh) | 油润湿致密砂岩核磁共振测井流体识别方法及其图版 | |
CN109856688A (zh) | 基于核磁测井双tw极化增强法的流体性质识别方法 | |
CN115935129B (zh) | 一种土壤分尺度重金属浓度值确定方法和装置 | |
CN114528729A (zh) | 一种基于多尺度耦合的潜山裂缝气藏产量预测方法 | |
CN110227968A (zh) | 基于小波变换的磁流变去除函数寄生条纹消除方法及装置 | |
CN114112856A (zh) | 一种含铀砂岩全尺寸孔径结构表征方法及*** | |
CN117974754A (zh) | 基于图像处理的储层溶洞定量识别方法、电子设备及介质 | |
CN110441204B (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 |