CN113984589A - 一种计算岩石迂曲度和气体扩散系数的方法 - Google Patents
一种计算岩石迂曲度和气体扩散系数的方法 Download PDFInfo
- Publication number
- CN113984589A CN113984589A CN202111284159.7A CN202111284159A CN113984589A CN 113984589 A CN113984589 A CN 113984589A CN 202111284159 A CN202111284159 A CN 202111284159A CN 113984589 A CN113984589 A CN 113984589A
- Authority
- CN
- China
- Prior art keywords
- core
- pore
- gas
- rock
- diffusion coefficient
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 88
- 238000009792 diffusion process Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000011148 porous material Substances 0.000 claims abstract description 136
- 239000012530 fluid Substances 0.000 claims abstract description 46
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 30
- 238000009826 distribution Methods 0.000 claims abstract description 26
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 18
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 239000007789 gas Substances 0.000 claims description 70
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 30
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 239000003345 natural gas Substances 0.000 claims description 11
- 239000001569 carbon dioxide Substances 0.000 claims description 8
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 8
- 239000001307 helium Substances 0.000 claims description 8
- 229910052734 helium Inorganic materials 0.000 claims description 8
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims description 8
- 239000003921 oil Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012417 linear regression Methods 0.000 claims description 4
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 claims description 3
- 238000001816 cooling Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 8
- 238000002474 experimental method Methods 0.000 abstract description 6
- 238000005516 engineering process Methods 0.000 abstract description 3
- 239000003208 petroleum Substances 0.000 abstract description 2
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical class O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 8
- 239000012267 brine Substances 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 229920001971 elastomer Polymers 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 238000001069 Raman spectroscopy Methods 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000035508 accumulation Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- 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
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N2013/003—Diffusion; diffusivity between liquids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Dispersion Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出了一种计算岩石迂曲度和气体扩散系数的方法,属于地质与石油工程技术领域。该方法基于核磁共振岩心分析技术和分形理论,利用岩心T2分布对岩心复杂孔隙通道结构进行准确的定量表征,得到岩心的平均孔隙半径和孔隙数目,进而得到岩心孔隙通道迂曲度。再根据扩散迂曲度的定义式,通过岩心孔隙通道迂曲度对气体在大体积纯流体中的扩散系数进行修正,得到气体在饱和岩心中的有效扩散系数。本发明不需要通过高温高压复杂的直接测量实验即可获得岩石迂曲度和气体扩散系数,具有无损、准确、省时、省力、省钱等显著优势,完全满足岩心分析行业标准化和批量化要求,适合在工业实验室推广应用。
Description
技术领域
本发明属于地质与石油工程技术领域,涉及一种计算岩石迂曲度和气体扩散系数的方法。
背景技术
气体扩散是由于分子(原子等)的热运动而产生的传质现象,气体分子在浓度差作用下,从高浓度区通过各种介质向低浓度区自由迁移达到动态平衡的一种物理过程,只要有浓度差,就有气体分子扩散。气体扩散作用贯穿于地下气体资源(天然气、二氧化碳、氦气等)形成、运移、聚集和散失的整个地质过程。因此,扩散系数是定量评估地下气体资源的扩散损失量、盖层封闭能力和气藏保存时间的关键参数,对研究天然气藏(如常规天然气、致密天然气、煤层气、页岩气、天然气水合物等)、二氧化碳气藏、氦气藏等地下多种气体资源的成藏与开发理论具有重要指导作用。
扩散系数是表示气体扩散作用强弱的物理量,根据菲克定律,扩散系数是沿扩散方向,在单位时间每单位浓度梯度的条件下,垂直通过单位面积的气体量。由于地下岩石孔隙中通常被流体(油、气、水等)所饱和,气体(天然气、二氧化碳、氦气等)在岩石中的扩散主要是通过孔隙流体介质的扩散,与在纯流体中的扩散相比,气体在流体饱和岩石中的扩散还受到岩石内孔隙通道的连通性和弯曲程度限制,导致扩散系数严重下降,由于孔隙通道的结构复杂性和气体分散的随机性,增加了实验测量的难度。
目前,常规的岩石中气体扩散系数测量方法是用岩心套筒(橡胶或铅套筒)并施加外部围压(根据模拟地层压力设定测量实验围压)的方法将岩心(柱塞状岩石样品)圆柱侧表面包裹紧,将岩心固定在夹持器中,抽真空,根据模拟地层温度设定测量实验温度,岩心两端同时注入不同种类高压气体(比如一端是甲烷、另一端是氮气)并控制压力平衡,然后,利用气相色谱仪等仪器测量一定时间内通过岩心的气体扩散量或浓度,通过菲克定律来计算扩散系数。上述测量实验所用岩心如果太短,容易被压碎,如果太长,由于气体扩散速度很慢,实验耗时太长,因此,通常控制在1.5~2.5cm,即便如此,测量一块岩心也要5~30天(岩心越致密,测量时间越长)。显然,如此长的测量时间中,实验***可能漏气,岩心两端压力会出现不平衡,气体会从岩心与橡胶套筒之间的缝隙慢慢渗透,气相色谱仪的精度以及取样点的位置等,以上这些都会严重影响测量结果的准确性,另外,该方法将岩心作为一个整体“黑箱”模型,不能反映岩石内部孔隙通道结构对气体扩散系数的影响。而且,高温高压测量实验***造价与使用费用昂贵,耗费大量人力物力。因此,上述方法不适合批量岩心实验测量。
发明内容
为了解决上述问题,本发明提供了一种基于核磁共振T2分布计算岩心孔隙通道迂曲度,再通过岩心迂曲度计算气体在完全饱和岩心中的有效扩散系数的方法。该方法利用低场核磁共振技术测量完全被流体(水、油等含氢流体)饱和的岩心弛豫时间T2分布,进而获得岩心的孔隙分布。假设岩心中孔隙通道是由相同长度与相同孔径的平行通道组成,基于岩心T2分布和多孔介质分形理论,计算得到岩心孔隙数目与平均孔隙直径,再通过岩心孔隙体积计算孔隙平均长度,进而得到岩心的孔隙通道迂曲度。通过岩心孔隙通道迂曲度修正气体(天然气、二氧化碳、氦气等)在大体积介质(如水、油、气等)中的扩散系数,从而得到气体在完全饱和岩心中的有效扩散系数。
本发明采用的技术方案是:
一种计算岩石迂曲度和气体扩散系数的方法,包括以下步骤:
步骤1:岩心样品准备
将待测岩石制成柱塞状岩心样品,将岩心放入恒温箱内加热烘干内部水分至恒重,然后放入干燥器中冷却至室温,测量岩心直径、长度和干重。再将岩心放入流体(水、油等含氢流体)中,用真空泵抽真空至无气泡溢出,让岩心孔隙被流体完全饱和。取出已经完全饱和的岩心,并测量岩心湿重。
步骤2:利用核磁共振获取完全流体饱和岩心的弛豫时间T2分布
利用低场核磁共振CPMG脉冲序列获取步骤1得到的完全饱和岩心的弛豫时间T2分布;通过岩心的弛豫时间T2分布可以得到岩心中孔隙流体体积为Vw;另外,根据步骤1所测量的岩心直径、长度数据计算得到岩心整体体积为V,则岩心孔隙率为
步骤3:基于岩心的T2分布计算其孔隙分形维数
步骤2中已经测得完全流体饱和岩心的T2分布,用公式(1)计算出不同弛豫时间T2对应的孔隙流体体积VT2所占总孔隙流体体积Vw的百分比Sv:
当仅用一种流体饱和岩心多孔介质时,由多孔介质孔隙流体的NMR横向弛豫(T2弛豫)机制可以推导出弛豫时间T2与孔隙半径r关系如公式(2)。由岩石物理学可知,多孔介质毛细管压力和孔径间存在关系如公式(3)。因此通过公式(2)、(3)推导出多孔介质毛细管压力Pc与弛豫时间T2之间的关系,Pc与T2呈反比关系,如公式(4)所示:
其中,Fs是几何形状因子,当孔隙为球状时Fs=3,孔隙为柱状时Fs=2;ρ是孔隙表面的T2弛豫强度(μm/ms);r是孔隙半径(cm);Pc是当孔径为r时的毛细管压力(Mpa);σ为孔隙流体的表面张力(mN/m);θ为孔隙流体与孔隙表面间的接触角。
根据多孔介质分形理论,可以得出岩心毛管压力曲线的分形表达式(5),再通过公式(4)、(5)推导出岩心基于低场核磁共振T2分布的分形表达式(6):
其中,Pcmin为孔径为r时的毛细管压力,T2max为岩心中最大孔隙所对应的横向弛豫时间,Df为岩心孔隙的分形维数。Pc与T2呈反比关系,当T2为最大值T2max时Pc为最小值Pcmin。
对式(6)两边取对数得:
lgSv=(3-Df)lgT2+(Df-3)lgT2max(7)
用上式(7)直接对岩心T2分布的lgT2与lgSv数据进行线性回归分析,得到岩心孔隙的分形维数Df。
步骤4:基于T2分布计算岩心孔隙通道迂曲度
由弛豫时间T2与孔隙半径r关系公式(2),可得到弛豫时间T2与孔隙直径λ关系公式,T2与孔隙直径λ理论上呈正比关系:
基于多孔介质分形理论和岩心T2分布计算岩心的孔隙通道数目与平均孔隙直径。具体为:根据岩心的孔隙表面积与孔隙大小分布满足分形标度规律,得到孔隙直径λ大于或等于a(a为小于λmax的任意正数)的累积孔隙数目N满足如下的分形标度关系:
式中,0<Df<2(在二维空间里),或0<Df<3(在三维空间里)。
将公式(9)中的a取值为最小孔隙直径λmin,则得到岩心中的总孔隙数目为:
式中λmax、λmin分别为最大孔隙通道直径、最小孔隙通道直径。
平均孔隙直径λav采用孔隙分布概率密度函数方程公式(11)计算得到。孔隙分布概率公式表示在孔隙直径为λ的孔隙数目占总孔隙数目的百分比,推导出平均孔隙直径λav公式如公式(12)所示:
假设岩心中孔隙通道是由多条相同长度、相同孔径的平行通道组成。经过以上计算,已获得岩心的平均孔隙直径λav和总孔隙数目N,并通过岩心样品孔隙率和体积,可获取岩心孔隙体积。因此,通过岩心平均孔隙直径λav、总孔隙数目N和孔隙体积计算得到岩心孔隙通道平均长度:
因此,岩心孔隙通道迂曲度计算公式为:
步骤5:基于岩心孔隙通道迂曲度计算气体在饱和岩心中的有效扩散系数
通过查阅资料获得不同温压条件下气体(天然气、二氧化碳、氦气等)在大体积纯流体(油、气、水等)中的扩散系数,则岩心中气体有效扩散系数为:
式中,De为气体在流体饱和岩心中有效扩散系数,D为气体在大体积纯流体中的扩散系数。
本发明的有益效果:本发明所述方法基于核磁共振岩心分析技术和分形理论,利用岩心T2分布对岩心复杂孔隙通道结构进行准确的定量表征,得到岩心的平均孔隙半径和孔隙数目,进而得到岩心孔隙通道迂曲度。再根据扩散迂曲度的定义式,通过岩心孔隙通道迂曲度对气体在大体积纯流体(油、气、水等)中的扩散系数进行修正,得到气体在饱和岩心中的有效扩散系数。本发明不需要通过高温高压复杂的直接测量实验即可获得岩石迂曲度和气体扩散系数,具有无损、准确、省时、省力、省钱等显著优势,完全满足岩心分析行业标准化和批量化要求,适合在工业实验室推广应用。
附图说明
图1是方法流程框图。
图2是完全盐水饱和状态下岩心核磁共振T2分布。
图3是岩心孔隙流体体积VT2所占总孔隙体积Vw的百分比Sv分布。
图4是岩心lgT2与lgSv线性回归分析。
具体实施方式
下面结合附图和技术方案对本发明作进一步说明。
一种计算岩石迂曲度和气体扩散系数的方法,包括以下步骤:
步骤1:岩心样品准备
选取典型的露头砂岩Berea(渗透率350~450md)作为将待测岩石,将待测岩石制成柱塞状岩心,将岩心放入恒温箱内加热110℃烘干内部水分至恒重,然后放入干燥器中冷却至室温,测量岩心直径2.52cm、长度5.03cm、干重51.5g。然后将岩心放入盐水中,用真空泵抽真空至无气泡溢出,保持真空48h以上,让岩心充分饱和盐水,取出已经完全盐水饱和的岩心,并测量岩心湿重57.0g。
步骤2:获取完全水饱和状态下岩心的T2分布
利用低场核磁共振CPMG脉冲序列获取步骤1得到的完全盐水饱和岩心的弛豫时间T2分布;通过标定得到岩心总孔隙水体积为5.35cm3,另外,根据所测量的岩心直径、长度数据计算得到岩心整体体积为25.17cm3,则岩心孔隙率为0.21。
步骤3:基于T2分布计算孔隙分形维数
步骤2中已经测得完全盐水饱和状态下岩心T2分布,计算出不同T2值时对应的孔隙流体体积VT2所占总孔隙水体积Vw的百分比Sv。直接对岩心T2分布的lgT2与lgSv数据进行线性回归分析,拟合直线的斜率K=3-Df,解得岩心孔隙分形维数Df为2.37,相关系数0.99。
步骤4:基于T2数据计算岩心孔隙通道迂曲度
根据核磁共振原理,完全盐水饱和岩心的T2分布可以反映岩石内部结构,横向弛豫时间T2与孔径r理论上呈线性正比关系。因此,通过饱和岩心T2分布,可以得到岩心的孔径分布。基于多孔介质分形理论和岩心T2分布计算得到,岩心中的孔隙数目N=4.94*107,平均孔隙直径λav=4.88*10-7(m),计算的得到孔隙通道平均长度Lav=0.14(m),迂曲度为τ=2.83。此外,通过电阻率法测得该岩心的迂曲度为3.10,与本方法计算结果吻合较好。
步骤5:基于岩心孔隙通道迂曲度计算气体在饱和岩心中的有效扩散系数
根据文献资料(Guo HR,Chen Y,Lu WJ,et al.In situ Raman spectroscopicstudy of diffusion coefficients of methane in liquid water under highpressure and wide temperatures),当地下深度为3km时,压力31Mpa,温度为403K,甲烷在纯盐水中的扩散系数D=9.34*10-9(m2/s),则计算得到在该条件下甲烷在纯盐水饱和岩心中的有效扩散系数De=3.30*10-9(m2/s)。
Claims (8)
1.一种计算岩石迂曲度和气体扩散系数的方法,其特征在于,该方法包括以下步骤:
步骤一:岩心样品准备
将待测岩石制成柱塞状岩心样品,将岩心放入恒温箱内加热烘干内部水分至恒重,然后冷却至室温,测量岩心直径、长度和干重;再将岩心放入流体中,使岩心孔隙被流体完全饱和,测量岩心湿重;
步骤二:利用核磁共振获取完全流体饱和岩心的弛豫时间T2分布
利用低场核磁共振CPMG脉冲序列获取步骤一得到的完全流体饱和岩心的弛豫时间T2分布;通过岩心的弛豫时间T2分布得到岩心中孔隙流体体积,由岩心中孔隙流体体积和岩心整体体积计算得到岩心孔隙率;
步骤三:基于核磁共振理论和多孔介质分形理论,通过岩心弛豫时间T2分布计算其孔隙分形维数;
步骤四:基于T2分布计算岩心孔隙通道迂曲度
通过核磁共振T2获得岩心孔隙半径分布;基于多孔介质分形理论和岩心T2分布计算岩心的总孔隙数目与平均孔隙直径;通过岩心平均孔隙直径、总孔隙数目和孔隙体积计算得到岩心孔隙通道平均长度,进而得到岩心孔隙通道迂曲度;
步骤五:通过岩心孔隙通道迂曲度对气体在大体积纯流体中的扩散系数进行修正,得到气体在饱和岩心中的有效扩散系数。
2.根据权利要求1所述的一种计算岩石迂曲度和气体扩散系数的方法,其特征在于,所述步骤三中,计算孔隙分形维数的方法为:
基于核磁共振理论,用公式(1)计算出不同弛豫时间T2对应的孔隙流体体积VT2所占总孔隙流体体积Vw的百分比Sv:
多孔介质毛细管压力Pc与弛豫时间T2之间的关系如式(2)所示:
式中,Fs是几何形状因子,当孔隙为球状时Fs=3,孔隙为柱状时Fs=2;ρ是孔隙表面的T2弛豫强度,μm/ms;Pc是当孔径为r时的毛细管压力,Mpa;σ为孔隙流体的表面张力,mN/m;θ为孔隙流体与孔隙表面间的接触角;
根据多孔介质分形理论,得出岩心毛管压力曲线的分形表达式(3),再通过公式(2)、(3)推导出岩心基于低场核磁共振T2分布的分形表达式(4):
其中,Pcmin为孔径为r时的毛细管压力,T2max为岩心中最大孔隙所对应的横向弛豫时间,Df为岩心孔隙的分形维数;Pc与T2呈反比关系,当T2为最大值T2max时Pc为最小值Pcmin;
对式(4)两边取对数得:
lgSv=(3-Df)lgT2+(Df-3)lgT2max (5)
用式(5)直接对岩心T2分布的lgT2与lgSv数据进行线性回归分析,得到岩心孔隙的分形维数Df。
6.根据权利要求1、2或5所述的一种计算岩石迂曲度和气体扩散系数的方法,其特征在于,所述的气体包括天然气、二氧化碳或氦气;所述流体包括油、气或水。
7.根据权利要求3所述的一种计算岩石迂曲度和气体扩散系数的方法,其特征在于,所述的气体包括天然气、二氧化碳或氦气;所述流体包括油、气或水。
8.根据权利要求4所述的一种计算岩石迂曲度和气体扩散系数的方法,其特征在于,所述的气体包括天然气、二氧化碳或氦气;所述流体包括油、气或水。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111284159.7A CN113984589B (zh) | 2021-11-01 | 2021-11-01 | 一种计算岩石迂曲度和气体扩散系数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111284159.7A CN113984589B (zh) | 2021-11-01 | 2021-11-01 | 一种计算岩石迂曲度和气体扩散系数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113984589A true CN113984589A (zh) | 2022-01-28 |
CN113984589B CN113984589B (zh) | 2023-05-12 |
Family
ID=79745442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111284159.7A Active CN113984589B (zh) | 2021-11-01 | 2021-11-01 | 一种计算岩石迂曲度和气体扩散系数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113984589B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115598023A (zh) * | 2022-11-08 | 2023-01-13 | 四川中德禄江装备科技有限公司(Cn) | 一种确定松散破碎岩土体三维迂曲度方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5696448A (en) * | 1995-06-26 | 1997-12-09 | Numar Corporation | NMR system and method for formation evaluation using diffusion and relaxation log measurements |
US20080303520A1 (en) * | 2007-06-08 | 2008-12-11 | University Of New Brunswick | Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks |
CN102879306A (zh) * | 2012-10-10 | 2013-01-16 | 大连理工大学 | 一种利用磁共振成像技术检测气液扩散过程的装置及其方法 |
CN111563927A (zh) * | 2020-05-14 | 2020-08-21 | 西南石油大学 | 一种基于岩石微ct图像的孔隙迂曲度计算方法 |
US20210223192A1 (en) * | 2018-06-01 | 2021-07-22 | Board Of Regents, The University Of Texas System | Method for assessment of pore-throat size distribution and permeability in porous media |
-
2021
- 2021-11-01 CN CN202111284159.7A patent/CN113984589B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5696448A (en) * | 1995-06-26 | 1997-12-09 | Numar Corporation | NMR system and method for formation evaluation using diffusion and relaxation log measurements |
US20080303520A1 (en) * | 2007-06-08 | 2008-12-11 | University Of New Brunswick | Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks |
CN102879306A (zh) * | 2012-10-10 | 2013-01-16 | 大连理工大学 | 一种利用磁共振成像技术检测气液扩散过程的装置及其方法 |
US20210223192A1 (en) * | 2018-06-01 | 2021-07-22 | Board Of Regents, The University Of Texas System | Method for assessment of pore-throat size distribution and permeability in porous media |
CN111563927A (zh) * | 2020-05-14 | 2020-08-21 | 西南石油大学 | 一种基于岩石微ct图像的孔隙迂曲度计算方法 |
Non-Patent Citations (8)
Title |
---|
ZHAOWEN LI 等: "Experimental Study of Carbon Dioxide Diffusion in Oil-Saturated Porous Media under Reservoir Conditions", 《INDUSTRIAL & ENGINEERING CHEMISTRY RESEARCH》 * |
于水: "《墙体材料热湿物性参数研究》", 31 August 2018, 东北大学出版社 * |
刘乐乐等: "分形理论在天然气水合物研究领域的应用", 《海洋地质前沿》 * |
张超谟等: "基于核磁共振T2谱分布的储层岩石孔隙分形结构研究", 《石油天然气学报》 * |
河海大学《水利大辞典》编辑修订委员会: "《水利大辞典》", 31 October 2015, 上海辞书出版社 * |
王启立等: "《煤基多孔介质的表征、流动及应用》", 28 February 2020, 中国矿业大学出版社 * |
葛新民: "非均质碎屑岩储层孔隙结构表征及测井精细评价研究", 《中国优秀博硕士学位论文全文数据库(博士) 基础科学辑》 * |
陈振标等: "利用NMR T2 谱分布研究储层岩石孔隙分形结构", 《岩性油气藏》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115598023A (zh) * | 2022-11-08 | 2023-01-13 | 四川中德禄江装备科技有限公司(Cn) | 一种确定松散破碎岩土体三维迂曲度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113984589B (zh) | 2023-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Geffen et al. | Experimental investigation of factors affecting laboratory relative permeability measurements | |
US10288571B2 (en) | Absolute porosity and pore size determination of pore types in media with varying pore sizes | |
Shen et al. | Water imbibition and drainage of high rank coals in Qinshui Basin, China | |
CN109443867B (zh) | 一种对致密岩石的物性参数进行连续检测的方法 | |
CN107991335B (zh) | 致密砂岩水锁损害评价试验方法 | |
CN103926267A (zh) | 一种定量评价应力敏感过程中孔喉变化程度的方法 | |
CN114136861B (zh) | 一种储气库近井地带干化盐析效应实验***及评价方法 | |
Carles et al. | Low permeability measurements using steady-state and transient methods | |
US20220291157A1 (en) | Method for establishing mathematical model of relationship between spontaneous imbibition volume and time of porous medium | |
MAO et al. | An improved method of using NMR T2 distribution to evaluate pore size distribution | |
CN113431537A (zh) | 一种非稳态变流速大尺度岩心水驱气相对渗透率测试方法 | |
CN113984589B (zh) | 一种计算岩石迂曲度和气体扩散系数的方法 | |
CN209821099U (zh) | 基于核磁共振的多功能致密气储层动态参数联测装置 | |
Zhou et al. | The determination of surface relaxivity and application to coal spontaneous imbibition | |
Wang et al. | Gas-phase relative permeability characterization on tight-gas samples | |
CN113916745A (zh) | 一种无损测量水驱气藏微观孔隙结构变化规律的实验方法 | |
CN113984590B (zh) | 一种计算非均质岩石空间迂曲度和气体扩散系数分布的方法 | |
CN110672495A (zh) | 一种基于低场磁共振技术的水泥基材料水分渗透率的预测方法 | |
CN115146483A (zh) | 一种新的可动流体动用下限计算方法 | |
Sun et al. | Experimental investigation on shale gas transport characteristics in nanopores under high temperature and high pressure | |
Grattoni et al. | An improved technique for deriving drainage capillary pressure from NMR T2 distributions | |
CN109142163B (zh) | 一种饱和度与渗透率耦合检测的水锁解除能力评价方法 | |
CN113109234A (zh) | 重水核磁束缚水探测的低含油饱和度渗流规律修正方法 | |
Lombard et al. | Measurement of capillary pressure curves at reservoir conditions | |
Zhao et al. | Experimental investigation on spontaneous counter‐current imbibition in water‐wet natural reservoir sandstone core using MRI |
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 |