CN113984589B - 一种计算岩石迂曲度和气体扩散系数的方法 - Google Patents

一种计算岩石迂曲度和气体扩散系数的方法 Download PDF

Info

Publication number
CN113984589B
CN113984589B CN202111284159.7A CN202111284159A CN113984589B CN 113984589 B CN113984589 B CN 113984589B CN 202111284159 A CN202111284159 A CN 202111284159A CN 113984589 B CN113984589 B CN 113984589B
Authority
CN
China
Prior art keywords
core
pore
gas
tortuosity
fluid
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
Application number
CN202111284159.7A
Other languages
English (en)
Other versions
CN113984589A (zh
Inventor
赵越超
宋永臣
黄明星
毛晓坤
张毅
杨明军
郑嘉男
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202111284159.7A priority Critical patent/CN113984589B/zh
Publication of CN113984589A publication Critical patent/CN113984589A/zh
Application granted granted Critical
Publication of CN113984589B publication Critical patent/CN113984589B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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，则岩心孔隙率为
Figure BDA0003332365390000031
步骤3：基于岩心的T2分布计算其孔隙分形维数
步骤2中已经测得完全流体饱和岩心的T2分布，用公式(1)计算出不同弛豫时间T2对应的孔隙流体体积VT2所占总孔隙流体体积Vw的百分比Sv：
Figure BDA0003332365390000032
当仅用一种流体饱和岩心多孔介质时，由多孔介质孔隙流体的NMR横向弛豫(T2弛豫)机制可以推导出弛豫时间T2与孔隙半径r关系如公式(2)。由岩石物理学可知，多孔介质毛细管压力和孔径间存在关系如公式(3)。因此通过公式(2)、(3)推导出多孔介质毛细管压力Pc与弛豫时间T2之间的关系，Pc与T2呈反比关系，如公式(4)所示：
Figure BDA0003332365390000033
Figure BDA0003332365390000034
Figure BDA0003332365390000035
其中，Fs是几何形状因子，当孔隙为球状时Fs＝3，孔隙为柱状时Fs＝2；ρ是孔隙表面的T2弛豫强度(μm/ms)；r是孔隙半径(cm)；Pc是当孔径为r时的毛细管压力(Mpa)；σ为孔隙流体的表面张力(mN/m)；θ为孔隙流体与孔隙表面间的接触角。
根据多孔介质分形理论，可以得出岩心毛管压力曲线的分形表达式(5)，再通过公式(4)、(5)推导出岩心基于低场核磁共振T2分布的分形表达式(6)：
Figure BDA0003332365390000041
Figure BDA0003332365390000042
其中，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与孔隙直径λ理论上呈正比关系：
Figure BDA0003332365390000043
式中，λ为孔隙直径，
Figure BDA0003332365390000044
为转换系数，通过T2分布可以获取岩心孔隙直径λ分布。
基于多孔介质分形理论和岩心T2分布计算岩心的孔隙通道数目与平均孔隙直径。具体为：根据岩心的孔隙表面积与孔隙大小分布满足分形标度规律，得到孔隙直径λ大于或等于a(a为小于λmax的任意正数)的累积孔隙数目N满足如下的分形标度关系：
Figure BDA0003332365390000051
式中，0<Df<2(在二维空间里)，或0<Df<3(在三维空间里)。
将公式(9)中的a取值为最小孔隙直径λmin，则得到岩心中的总孔隙数目为：
Figure BDA0003332365390000052
式中λmax、λmin分别为最大孔隙通道直径、最小孔隙通道直径。
平均孔隙直径λav采用孔隙分布概率密度函数方程公式(11)计算得到。孔隙分布概率公式表示在孔隙直径为λ的孔隙数目占总孔隙数目的百分比，推导出平均孔隙直径λav公式如公式(12)所示：
Figure BDA0003332365390000053
Figure BDA0003332365390000054
由于
Figure BDA0003332365390000055
远大于
Figure BDA0003332365390000056
方程(12)化简为：
Figure BDA0003332365390000057
假设岩心中孔隙通道是由多条相同长度、相同孔径的平行通道组成。经过以上计算，已获得岩心的平均孔隙直径λav和总孔隙数目N，并通过岩心样品孔隙率和体积，可获取岩心孔隙体积。因此，通过岩心平均孔隙直径λav、总孔隙数目N和孔隙体积计算得到岩心孔隙通道平均长度：
Figure BDA0003332365390000058
因此，岩心孔隙通道迂曲度计算公式为：
Figure BDA0003332365390000059
式中，Lav为岩心孔隙通道平均长度，L为岩心几何长度，
Figure BDA00033323653900000510
为岩心孔隙率，A为岩心横截面积，τ为岩心孔隙通道迂曲度。
步骤5：基于岩心孔隙通道迂曲度计算气体在饱和岩心中的有效扩散系数
通过查阅资料获得不同温压条件下气体(天然气、二氧化碳、氦气等)在大体积纯流体(油、气、水等)中的扩散系数，则岩心中气体有效扩散系数为：
Figure BDA0003332365390000061
式中，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 (3)

1.一种计算岩石迂曲度和气体扩散系数的方法，其特征在于，该方法包括以下步骤：
步骤一：岩心样品准备
将待测岩石制成柱塞状岩心样品，将岩心放入恒温箱内加热烘干内部水分至恒重，然后冷却至室温，测量岩心直径、长度和干重；再将岩心放入流体中，使岩心孔隙被流体完全饱和，测量岩心湿重；
步骤二：利用核磁共振获取完全流体饱和岩心的弛豫时间T2分布
利用低场核磁共振CPMG脉冲序列获取步骤一得到的完全流体饱和岩心的弛豫时间T2分布；通过岩心的弛豫时间T2分布得到岩心中孔隙流体体积，由岩心中孔隙流体体积和岩心整体体积计算得到岩心孔隙率；
步骤三：基于核磁共振理论和多孔介质分形理论，通过岩心弛豫时间T2分布计算其孔隙分形维数；计算孔隙分形维数的方法为：
基于核磁共振理论，用公式(1)计算出不同弛豫时间T2对应的孔隙流体体积VT2所占总孔隙流体体积Vw的百分比Sv：
Figure FDA0004062691630000011
多孔介质毛细管压力Pc与弛豫时间T2之间的关系如式(2)所示：
Figure FDA0004062691630000012
式中，Fs是几何形状因子，当孔隙为球状时Fs＝3，孔隙为柱状时Fs＝2；ρ是孔隙表面的T2弛豫强度，μm/ms；Pc是当孔径为r时的毛细管压力，Mpa；σ为孔隙流体的表面张力，mN/m；θ为孔隙流体与孔隙表面间的接触角；
根据多孔介质分形理论，得出岩心毛管压力曲线的分形表达式(3)，再通过公式(2)、(3)推导出岩心基于低场核磁共振T2分布的分形表达式(4)：
Figure FDA0004062691630000013
Figure FDA0004062691630000021
其中，Pcmin为孔径为r时的毛细管压力，T2max为岩心中最大孔隙所对应的横向弛豫时间，Df为岩心孔隙的分形维数；Pc与T2呈反比关系，当T2为最大值T2max时Pc为最小值Pcmin；
对式(4)两边取对数得：
lgSv＝(3-Df)lgT2+(Df-3)lgT2max (5)
用式(5)直接对岩心T2分布的lgT2与lgSv数据进行线性回归分析，得到岩心孔隙的分形维数Df；
步骤四：基于T2分布计算岩心孔隙通道迂曲度
通过核磁共振T2获得岩心孔隙半径分布；基于多孔介质分形理论和岩心T2分布计算岩心的总孔隙数目与平均孔隙直径；通过岩心平均孔隙直径、总孔隙数目和孔隙体积计算得到岩心孔隙通道平均长度，进而得到岩心孔隙通道迂曲度；计算岩心孔隙通道迂曲度的方法为：
设岩心中孔隙通道是由多条相同长度、相同孔径的平行通道组成；岩心多孔介质的总孔隙通道数目N与平均孔隙直径λav分别为：
Figure FDA0004062691630000022
Figure FDA0004062691630000023
式中，λmax、λmin分别为最大孔隙通道直径、最小孔隙通道直径；
再通过岩心平均孔隙直径λav、总孔隙数目N和孔隙体积计算得到岩心孔隙通道平均长度Lav：
Figure FDA0004062691630000024
最后，得到岩心孔隙通道迂曲度τ：
Figure FDA0004062691630000025
式中，Lav为岩心孔隙通道平均长度，L为岩心几何长度，Df为岩心孔隙分形维数，
Figure FDA0004062691630000032
为岩心孔隙率，A为岩心横截面积，τ为岩心孔隙通道迂曲度；
步骤五：通过岩心孔隙通道迂曲度对气体在大体积纯流体中的扩散系数进行修正，得到气体在饱和岩心中的有效扩散系数。
2.根据权利要求1所述的一种计算岩石迂曲度和气体扩散系数的方法，其特征在于，所述步骤五中，饱和岩心中的有效扩散系数的计算方法为：
通过不同温压条件下气体在大体积纯流体中的扩散系数，通过岩心孔隙通道迂曲度对气体在大体积纯流体中的扩散系数进行修正，得到气体在饱和岩心中的有效扩散系数：
Figure FDA0004062691630000031
式中，De为气体在流体饱和岩心中有效扩散系数，D为气体在大体积纯流体中的扩散系数。
3.根据权利要求1或2所述的一种计算岩石迂曲度和气体扩散系数的方法，其特征在于，所述的气体包括天然气、二氧化碳或氦气；所述流体包括油、气或水。
CN202111284159.7A 2021-11-01 2021-11-01 一种计算岩石迂曲度和气体扩散系数的方法 Active CN113984589B (zh)

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 CN113984589A (zh) 2022-01-28
CN113984589B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115598023B (zh) * 2022-11-08 2023-04-07 昆明理工大学 一种确定松散破碎岩土体三维迂曲度方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
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
US7567079B2 (en) * 2007-06-08 2009-07-28 University Of New Brunswick Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks
CN102879306B (zh) * 2012-10-10 2014-12-10 大连理工大学 一种利用磁共振成像技术检测气液扩散过程的方法
US12025576B2 (en) * 2018-06-01 2024-07-02 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图像的孔隙迂曲度计算方法

Also Published As

Publication number Publication date
CN113984589A (zh) 2022-01-28

Similar Documents

Publication Publication Date Title
CN107991335B (zh) 致密砂岩水锁损害评价试验方法
Geffen et al. Experimental investigation of factors affecting laboratory relative permeability measurements
CN105910971B (zh) 富有机质致密岩心气体渗透率和扩散系数的联测方法
CN110595953B (zh) 一种页岩混合润湿性的实验测试装置及方法
Shen et al. Water imbibition and drainage of high rank coals in Qinshui Basin, China
CN109001243B (zh) 一种采用低场核磁共振评价煤的动态水锁效应的方法与装置
CN109187615A (zh) 一种地层压力条件下岩石纳米孔径分布测量装置与方法
CN103926267A (zh) 一种定量评价应力敏感过程中孔喉变化程度的方法
WO2018164953A1 (en) Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr
CN108827853B (zh) 基于核磁共振的致密储层岩电测量装置及测量方法
CN109959672B (zh) 一种二氧化碳置换页岩甲烷的定量评价方法及其应用
CN105910941B (zh) 基于压力板仪的冻土未冻水含量的测试方法
CN113984589B (zh) 一种计算岩石迂曲度和气体扩散系数的方法
CN113484216B (zh) 一种评估致密砂岩气藏水相返排率及合理返排压差的方法
Carles et al. Low permeability measurements using steady-state and transient methods
MAO et al. An improved method of using NMR T2 distribution to evaluate pore size distribution
CN113866069A (zh) 一种页岩岩心渗透率实验装置和方法
CN110672495A (zh) 一种基于低场磁共振技术的水泥基材料水分渗透率的预测方法
CN209821099U (zh) 基于核磁共振的多功能致密气储层动态参数联测装置
Schneebeli et al. Measurements of water potential and water content in unsaturated cystalline rock
CN117825430A (zh) 一种评价天然气藏不同生产制度下流体动用规律的方法
CN113984590B (zh) 一种计算非均质岩石空间迂曲度和气体扩散系数分布的方法
CN109115657B (zh) 一种饱和度与渗透率耦合检测的水锁解除能力评价装置
CN116754600A (zh) 一种模拟地层流体环境下导热系数校正装置
Sun et al. Experimental investigation on shale gas transport characteristics in nanopores under high temperature and high pressure

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