CN106525684A - 一种基于孔喉结构的致密砂岩克氏渗透率的校正方法 - Google Patents
一种基于孔喉结构的致密砂岩克氏渗透率的校正方法 Download PDFInfo
- Publication number
- CN106525684A CN106525684A CN201610961286.9A CN201610961286A CN106525684A CN 106525684 A CN106525684 A CN 106525684A CN 201610961286 A CN201610961286 A CN 201610961286A CN 106525684 A CN106525684 A CN 106525684A
- Authority
- CN
- China
- Prior art keywords
- permeability
- pressure
- pore
- throat
- gas
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 179
- 238000000034 method Methods 0.000 title claims abstract description 54
- 239000011148 porous material Substances 0.000 claims abstract description 134
- 239000011435 rock Substances 0.000 claims abstract description 54
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims abstract description 25
- 229910052753 mercury Inorganic materials 0.000 claims abstract description 25
- 238000012937 correction Methods 0.000 claims abstract description 12
- 238000002474 experimental method Methods 0.000 claims abstract description 10
- 238000006073 displacement reaction Methods 0.000 claims abstract description 7
- 238000002347 injection Methods 0.000 claims abstract description 7
- 239000007924 injection Substances 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims description 9
- 238000000926 separation method Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000009736 wetting Methods 0.000 claims description 2
- 239000002655 kraft paper Substances 0.000 claims 2
- 239000007789 gas Substances 0.000 description 75
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 13
- 238000005259 measurement Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 6
- 229910052757 nitrogen Inorganic materials 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000005406 washing Methods 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000012267 brine Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000001307 helium Substances 0.000 description 1
- 229910052734 helium Inorganic materials 0.000 description 1
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 1
- 238000000329 molecular dynamics simulation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
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
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种基于孔喉结构的致密砂岩克氏渗透率的校正方法。该方法包括:制备岩心;选取不同压力点对岩心气测渗透率;根据岩心的气测渗透率得到岩心的克氏渗透率;对岩心进行高压压汞实验,得到岩心的排驱压力、饱和度中值压力、饱和度中值半径、孔喉半径平均值、最大进汞饱和度;根据高压压汞实验得到的实验参数对孔喉结构进行分类;对每一类的孔喉结构的岩心,进行气测渗透率与克氏渗透率的线性拟合,得到拟合参数,并将拟合参数的结果与平均孔隙压力进行指数拟合,得到克氏渗透率与平均孔隙压力和气测渗透率的函数关系式,完成基于孔喉结构的致密砂岩克氏渗透率的校正。
Description
技术领域
本发明涉及一种基于孔喉结构的致密砂岩克氏渗透率的校正方法,属于储层岩石物理性质评价领域。
背景技术
渗透率是表征储层允许通过的能力,其大小是对储层分类及评价的一个重要参数。目前,研究者经常采用气测渗透率来表征储层的渗流能力,但是,由于在致密储层中气体滑脱效应的存在,导致采用不同气体以及在不同压力下测得的渗透率都不一样,对于研究致密储层带来很大的麻烦。
针对稳态法气体渗透性测试中存在的问题,Jones(Jones,S.,1972.A rapidunsteady-state Klinkenberg permeameter.Soc.Pet.Eng.J.383-397.)提出了一种简单、非稳态装置和合适的理论来测量计算岩心的克氏渗透率、滑脱因子和湍流系数。但是,这种计算方法是基于线性拟合来找到一个非线性误差函数的收敛点,缺乏理论依据和收敛条件。Fernando(FA Pazos,A Bhaya,ALM Compan.Calculation of Klinkenbergpermeability,slip factor and turbulence factor of core plugs via nonlinearregression.Journal of Petroleum Science&Engineering.2009,67(3):159-167.)等针对Jones方法的不足,提出了一种简单的算法,该方法基于非线性拟合来计算未知参数,并且相比Jones方法有着理论依据的优势和较少的收敛条件。Firouzi等利用非平衡态分子动力学模拟滑脱效应研究氦气在页岩中的输运和渗透率参数。AI-Bulushi等(IR AI-Bulushi,RS AI-Maamari,OB Wilson.Brine versus Klinkenberg corrected gaspermeability correlation for Shuaiba carbonate formation[J].Journal ofPetroleum Science&Engineering.2012,s 92-93(4):24-29.)针对Shuaiba地层,建立了液测渗透率与经过克林肯贝格校正的气测渗透率之间的相关性,利用该相关性可以在已知经过校正的气测渗透率后计算得到液测渗透率。黄建章(黄建章,冯建明,陈心胜.获得克氏渗透率常规方法的简化[J].石油勘探与开发,1994,21(4):54-58.)等通过利用模型简化、公式推导、数据拟合得出一个数学公式,在已知某一平均压力p和此压力下的气测渗透率以及孔隙度下,便可计算得出克氏渗透率,大大简化了操作步骤,但是由于在推导公式过程中对岩心孔喉进行了简单假设,因此与实际存在一定的误差。
申请号为201510117076.7的中国专利申请,公开了一种致密砂岩气体渗透率测试装置及测试方法。该方法主要基于考虑气体滑脱效应的非稳态渗流理论的脉冲衰减法渗透率测试技术,通过测试岩心出口端的压力降落曲线便可同时获得克氏渗透率和滑脱因子,相比常规的五点法回归求取克氏渗透率,耗时短,但是该方法在拟合滑脱因子b与克氏渗透率k∞之间的关系采用经验公式,同时拟合数据量较少,容易产生较大的误差。
申请号为201510152794.8的中国专利申请提供了一种基于岩心气测渗透率测量装置的校正方法。该方法是利用岩心出口端流量与进出口稳定压力平方差的数据拟合出函数关系,利用该函数关系,同时结合达西定律,即可计算出极限压力下对应的进口压力和气测渗透率,最终建立起气测渗透率与平均压力的倒数的关系函数,该分段函数哦水平段截距即为克氏渗透率。该方法可准确获取岩心的克氏渗透率。但是该方法在测量过程中要求测量十组以上的数据,耗时多,同时要求净出口压差小并在小范围内波动,在实际测量中不容易实现。
综上所述,提供一种可以快速而准确的求得岩心的克氏渗透率的方法成为了本领域亟待解决的问题。
发明内容
为了解决上述技术问题,本发明的目的在于提供一种可以将岩心孔喉结构进行分类,研究不同孔喉结构类型下克氏渗透率与气测渗透率和平均孔隙压力的关系的克氏渗透率的校正方法,该方法快速、准确。
为了实现上述目的,本发明提供了一种基于孔喉结构的致密砂岩克氏渗透率的校正方法,该校正方法包括以下步骤:
制备岩心;
选取不同压力点对岩心气测渗透率;
根据岩心的气测渗透率得到岩心的克氏渗透率;
对岩心进行高压压汞实验,得到岩心的排驱压力、饱和度中值压力、饱和度中值半径、孔喉半径平均值、最大进汞饱和度;
在每一平均孔隙压力下,对每一类孔喉结构的岩心进行气测渗透率与克氏渗透率的线性拟合,得到拟合参数,并将不同平均孔隙压力下的拟合参数结果与平均孔隙压力进行指数函数拟合,得到克氏渗透率与平均孔隙压力和气测渗透率的函数关系式,完成基于孔喉结构的致密砂岩克氏渗透率的校正。
对于某一已知孔喉结构的岩心,将其某个平均孔隙压力下的气测渗透率代入上述的函数关系式中,即获得某一平均压力下的气测渗透率下的岩心的克氏渗透率。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,根据岩心的气测渗透率得到岩心的克氏渗透率时,根据如下公式进行计算:
其中,为平均压力时的气测渗透率,单位为μm2(D);K∞为克氏渗透率,单位为μm2(D);b为滑脱因子,单位为MPa;为平均孔隙压力,单位为MPa。
根据本发明的具体实施方式,将不同平均孔隙压力下的气测渗透率与平均压力的倒数进行线性拟合,得到的公式中截距即为所求岩心的克氏渗透率,同时利用斜率可求得岩心的滑脱因子b。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,进行高压压汞实验时的最大进汞压力为100MPa-200MPa。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,根据如下公式计算岩心的气测渗透率;
其中,为平均压力时的气测渗透率,单位为μm2(D),Q0为大气压下气体体积流量,单位为cm3/s;p0为大气压,单位为10-1MPa;p1为入口绝对压力,单位为10-1MPa;p2为出口绝对压力,单位为10-1MPa;μ为气体粘度,单位为mPa·s;L为岩心长度,单位为cm;A为岩心端口截面积,单位为cm2。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,排驱压力是指非润湿相开始连续进入岩心最大喉道时所对应的毛管压力。
根据本发明的具体实施方式,排驱压力是通过如下步骤得到的:
在半对数坐标中沿着毛管压力曲线平坦部分的第一个拐点做切线,切线延长与纵坐标轴相交的压力点即为排驱压力。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,饱和度中值压力是指进汞饱和度为50%时所对应的毛管压力。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,饱和度中值半径是指与饱和度中值压力相对应的喉道半径。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,喉道半径平均值是表示岩石平均孔喉半径大小的参数。根据如下公式进行计算:
其中,为平均孔喉半径,单位为μm;ri为第i个压力点所对应的喉道半径,单位为μm;si为第i个压力点所对应的累积进汞饱和度,i=1,2,3,4…,n。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,根据高压压汞实验得到的孔喉分选系数对孔喉结构进行分类,具体按照如下标准进行:
孔喉分选系数为1.1≦S<1.6,属于I类孔喉结构;
孔喉分选系数为1.6≦S<2.1,属于II类孔喉结构;
孔喉分选系数为2.1≦S<2.6,属于Ⅲ类孔喉结构。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,克氏渗透率与平均孔隙压力和气测渗透率的函数关系式按照如下步骤得到:
在每一平均孔隙压力下,对每一类孔喉结构的岩心进行气测渗透率与克氏渗透率的线性拟合k∞=akg+b,得到拟合参数a和b,并将不同平均孔隙压力下的拟合参数结果与平均孔隙压力进行指数函数拟合a(or b)=mpav -n,得到克氏渗透率与平均孔隙压力和气测渗透率的函数关系式
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,优选地,最大进汞饱和度是指最高实验压力时的汞饱和度值。
根据本发明的具体实施方式,对岩心进行气测渗透率时采用的装置包括:氮气瓶,流量控制器,第一压力传感器,第二压力传感器,第三压力传感器,手摇泵和气体流量计;
氮气瓶与流量控制器连通,氮气瓶与流量控制器之间设置有减压阀;流量控制器与岩心夹持器的入口连通;流量控制器与岩心夹持器的入口之间设置有单向阀和第一压力传感器;
岩心夹持器的围压口与手摇泵连通;岩心夹持器的围压口与手摇泵之间设置有第三压力传感器;
岩心夹持器的出口与气体流量计连通,岩心夹持器的出口与气体流量计之间设置有第二压力传感器,第二压力传感器和气体流量计之间设置有回压阀。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法中,气测岩心的渗透率参考《SYT 5336-2006岩心分析方法》执行。对岩心进行高压压汞实验,具体操作分析参考《SY/T 5346-2005岩石毛管压力曲线的测定》执行。
本发明提供的基于孔喉结构的致密砂岩克氏渗透率的校正方法,根据致密砂岩的孔喉结构特征,将岩心孔喉结构进行分类,研究不同孔喉结构类型下克氏渗透率与气测渗透率和平均孔隙压力的关系。
致密砂岩存在气体滑脱效应,孔喉结构不同,其滑脱效应程度也不相同,本发明的基于孔喉结构的致密砂岩克氏渗透率校正方法可用于测定致密砂岩的克氏渗透率,该方法主要考虑了岩心孔喉结构对滑脱效应产生的影响,针对不同的孔喉结构,研究了其克氏渗透率与气测渗透率和平均孔隙压力之间的关系。利用该方法,在已知岩心孔喉结构的基础上,代入某一平均孔隙压力及该压力下的气测渗透率即可求得岩心克氏渗透率,可以对储层进行物性评价提供较为快速准确的评价依据。
附图说明
图1为气测岩心渗透率的装置示意图;
图2为对岩心孔喉结构分类图;
图3a为Ⅰ类孔喉结构下的不同孔隙压力的气测渗透率与克氏渗透率的线性拟合曲线;
图3b为Ⅱ类孔喉结构下的不同孔隙压力的气测渗透率与克氏渗透率的线性拟合曲线;
图3c为Ⅲ类孔喉结构下的不同孔隙压力的气测渗透率与克氏渗透率的线性拟合曲线;
图4为线性拟合参数结果与平均孔隙压力的指数拟合;
图5a为Ⅰ类孔喉结构下的平均孔隙压力为0.2MPa-1MPa的气测渗透率与克氏渗透率的关系曲线;
图5b为Ⅱ类孔喉结构下的平均孔隙压力为0.2MPa-1MPa的气测渗透率与克氏渗透率的关系曲线;
图5c为Ⅲ类孔喉结构下的平均孔隙压力为0.2MPa-1MPa的气测渗透率与克氏渗透率的关系曲线。
主要附图符号说明
1氮气瓶 2减压阀 3流量控制器 4单向阀 501第一压力传感器 502第二压力传感器 503第三压力传感器 6岩心夹持器 7回压阀 8第一手摇泵 9气体流量计
具体实施方式
为了对本发明技术特征、目的和有益效果有更加清楚的理解,现对本发明的技术方案进行以下详细说明,但不能理解为对本发明的可实施范围的限定。
实施例
本实施例提供了一种基于孔喉结构的致密砂岩克氏渗透率的校正方法,该校正方法具体包括以下步骤:
一、制备岩心;
钻取天然岩心并进行编号;切岩心并磨平端面,保证岩心直径25mm,长度在60-70mm之间;使用索氏抽提器进行为期一至两个月的洗油,期间使用紫外分析仪检查洗油效果,如岩心未洗净,则其浸泡在365nm紫外光下会发出白色荧光,在这种情况下继续进行洗油,直至洗净为止;
烘干岩样至恒重,温度控制在104℃,烘干时间不小于8h,8h后每1h称量一次,两次称量的差值小于10mg;
用游标卡尺测量岩心长度;
二、选取不同压力点对岩心气测渗透率,具体选取如表1所示的数据,结果如表1所示。
表1
采用如图1所示的装置气测岩心的渗透率,该装置包括:氮气瓶1,流量控制器3,第一压力传感器501,第二压力传感器502,第三压力传感器503,第一手摇泵8和气体流量计9;
氮气瓶1与流量控制器3连通,氮气瓶1与流量控制器3之间设置有减压阀2;流量控制器3与岩心夹持器6的入口连通;流量控制器3与岩心夹持器6的入口之间设置有单向阀4和第一压力传感器501;
岩心夹持器6的围压口与第一手摇泵8连通;岩心夹持器6的围压口与第一手摇泵8之间设置有第三压力传感器503;
岩心夹持器6的出口与气体流量计9连通,岩心夹持器6的出口与气体流量计9之间设置有第二压力传感器502,第二压力传感器502和气体流量计9之间设置有回压阀7。
利用上述装置气测岩心渗透率的具体步骤包括:
装好设备;
加围压,加进口压力,等压力稳定(30min)后,采用气体流量计测量出口流量;测定岩心在入口压力、出口压力,围压分别为表1中的压力下测得渗透率;
使用下列公式计算气测渗透率;
其中,为平均压力为时的气测渗透率,单位为μm2(D),Q0为大气压下气体体积流量,单位为cm3/s;p0为大气压,单位为10-1MPa;p1为入口绝对压力,单位为10-1MPa;p2为出口绝对压力,单位为10-1MPa;μ为气体粘度,单位为mPa·s;L为岩心长度,单位为cm;A为岩心端口截面积,单位为cm2。
三、根据岩心的气测渗透率得到岩心的克氏渗透率;
以岩心的平均孔隙压力的倒数为横坐标,对应的测得的气测渗透率为纵坐标画散点图,采用公式对其进行线性拟合,得到的截距即为该岩心的克氏渗透率。
其中,为平均压力为时的气测渗透率,单位为μm2(D);K∞为克氏渗透率,单位为μm2(D);b为滑脱因子,单位为MPa;为平均孔隙压力,单位为MPa。表2为选取的30块岩心的孔隙度和各个压力下的气测渗透率以及计算得到的克氏渗透率。
表2
四、对岩心进行高压压汞实验,得到岩心的排驱压力、饱和度中值压力、饱和度中值半径、孔喉半径平均值、最大进汞饱和度,结果如表3所示。
表3
五、根据高压压汞实验得到的实验参数对孔喉结构进行分类,分类结果如图2所示。
六、对每一类的孔喉结构的岩心,采用公式k∞=akg+b进行气测渗透率与克氏渗透率的线性拟合,得到拟合参数(斜率a和截距b),并将所述拟合参数的结果与平均孔隙压力采用公式指数进行拟合,得到克氏渗透率与平均孔隙压力和气测渗透率的函数关系式(表4),完成基于孔喉结构的致密砂岩克氏渗透率的校正;
对于某一已知孔喉结构的岩心,将其某个平均孔隙压力下的气测渗透率代入上述的函数关系式中,即获得某一平均压力下的气测渗透率下的岩心的克氏渗透率。表5为针对三类孔喉结构选取了6块岩心,分别测出其在不同平孔隙压力下的气测渗透率和克氏渗透率,并利用公式计算出不同岩心的克氏渗透率,通过对比实验测得的克氏渗透率和利用本实施例的方法得到的克氏渗透率。
表4克氏渗透率与气测渗透率和平均孔隙压力的公式
表5不同方法求得的克氏渗透率对比表
通过表5可知,发现利用本实施例的方法得到的克氏渗透率与实验测得的克氏渗透率很接近,误差较小,而且平均孔隙压力越大,误差越小。
图3a为Ⅰ类孔喉结构下的不同孔隙压力的气测渗透率与克氏渗透率的线性拟合曲线;图3b为Ⅱ类孔喉结构下的不同孔隙压力的气测渗透率与克氏渗透率的线性拟合曲线;图3c为Ⅲ类孔喉结构下的不同孔隙压力的气测渗透率与克氏渗透率的线性拟合曲线;图4为线性拟合参数结果与平均孔隙压力的指数拟合;图5a为Ⅰ类孔喉结构下的平均孔隙压力为0.2MPa-1MPa的气测渗透率与克氏渗透率的关系曲线;图5b表为Ⅱ类孔喉结构下的平均孔隙压力为0.2MPa-1MPa的气测渗透率与克氏渗透率的关系曲线;图5c表为Ⅲ类孔喉结构下的平均孔隙压力为0.2MPa-1MPa的气测渗透率与克氏渗透率的关系曲线。
以上实施例说明,本发明的基于孔喉结构的致密砂岩克氏渗透率校正方法可用于测定致密砂岩的克氏渗透率,利用该方法,在已知岩心孔喉结构的基础上,代入某一平均孔隙压力及该压力下的气测渗透率即可求得岩心克氏渗透率,可以对储层进行物性评价提供较为快速准确的评价依据。
Claims (10)
1.一种基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,该校正方法包括以下步骤:
制备岩心;
选取不同压力点对岩心气测渗透率;
根据岩心的气测渗透率得到岩心的克氏渗透率;
对岩心进行高压压汞实验,得到岩心的排驱压力、饱和度中值压力、饱和度中值半径、孔喉半径平均值、最大进汞饱和度;
根据高压压汞实验得到的实验参数对孔喉结构进行分类;
在每一平均孔隙压力下,对每一类孔喉结构的岩心进行气测渗透率与克氏渗透率的线性拟合,得到拟合参数,并将不同平均孔隙压力下的拟合参数结果与平均孔隙压力进行指数函数拟合,得到克氏渗透率与平均孔隙压力和气测渗透率的函数关系式,完成基于孔喉结构的致密砂岩克氏渗透率的校正。
2.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,根据岩心的气测渗透率得到岩心的克氏渗透率时,根据如下公式进行计算:
其中,为平均压力为时的气测渗透率,单位为μm2(D);K∞为克氏渗透率,单位为μm2(D);b为滑脱因子,单位为MPa;为平均孔隙压力,单位为MPa。
3.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,进行高压压汞实验时的最大进汞压力为100MPa-200MPa。
4.根据权利要求1或2所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,根据如下公式计算岩心的气测渗透率;
其中,为平均压力时的气测渗透率,单位为μm2(D),Q0为大气压下气体体积流量,单位为cm3/s;p0为大气压,单位为10-1MPa;p1为入口绝对压力,单位为10- 1MPa;p2为出口绝对压力,单位为10-1MPa;μ为气体粘度,单位为mPa·s;L为岩心长度,单位为cm;A为岩心端口截面积,单位为cm2。
5.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,所述排驱压力是指非润湿相开始连续进入岩心最大喉道时所对应的毛管压力;所述饱和度中值压力是指进汞饱和度为50%时所对应的毛管压力。
6.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,所述饱和度中值半径是指与饱和度中值压力相对应的喉道半径。
7.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,所述喉道半径平均值根据如下公式进行计算:
其中,为平均孔喉半径,单位为μm;ri为第i个压力点所对应的喉道半径,单位为μm;si为第i个压力点所对应的累积进汞饱和度,i=1,2,3,4…,n。
8.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,根据高压压汞实验得到的孔喉分选系数对孔喉结构进行分类,具体按照如下标准进行:
孔喉分选系数为1.1≦S<1.6,属于I类孔喉结构;
孔喉分选系数为1.6≦S<2.1,属于II类孔喉结构;
孔喉分选系数为2.1≦S<2.6,属于Ⅲ类孔喉结构。
9.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,克氏渗透率与平均孔隙压力和气测渗透率的函数关系式按照如下步骤得到:
在每一平均孔隙压力下,对每一类孔喉结构的岩心进行气测渗透率与克氏渗透率的线性拟合k∞=akg+b,得到拟合参数a和b,并将不同平均孔隙压力下的拟合参数结果与平均孔隙压力进行指数函数拟合a(or b)=mpav -n,得到克氏渗透率与平均孔隙压力和气测渗透率的函数关系式
10.根据权利要求1所述的基于孔喉结构的致密砂岩克氏渗透率的校正方法,其特征在于,所述最大进汞饱和度是指最高实验压力时的汞饱和度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610961286.9A CN106525684B (zh) | 2016-10-27 | 2016-10-27 | 一种基于孔喉结构的致密砂岩克氏渗透率的校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610961286.9A CN106525684B (zh) | 2016-10-27 | 2016-10-27 | 一种基于孔喉结构的致密砂岩克氏渗透率的校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106525684A true CN106525684A (zh) | 2017-03-22 |
CN106525684B CN106525684B (zh) | 2019-04-16 |
Family
ID=58326238
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610961286.9A Active CN106525684B (zh) | 2016-10-27 | 2016-10-27 | 一种基于孔喉结构的致密砂岩克氏渗透率的校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106525684B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107144566A (zh) * | 2017-05-31 | 2017-09-08 | 中国石油大学(华东) | 利用化学‑沉积相表征细粒沉积岩非均质性的方法 |
CN109085104A (zh) * | 2018-07-23 | 2018-12-25 | 中国石油大学(北京) | 一种致密气藏的气体相对渗透率的确定方法 |
CN109490168A (zh) * | 2018-11-30 | 2019-03-19 | 重庆科技学院 | 页岩气藏不同传质扩散机理对储层渗流能力贡献率的定量评价方法 |
CN109827889A (zh) * | 2019-04-09 | 2019-05-31 | 北京艾迪佳业技术开发有限公司 | 一种克氏渗透率测量系统 |
CN109900610A (zh) * | 2017-12-07 | 2019-06-18 | 核工业北京地质研究院 | 一种多孔材料的联合测试方法 |
CN111398116A (zh) * | 2020-03-04 | 2020-07-10 | 中国地质大学(武汉) | 一种利用特定方向压汞法表征页岩各向异性的方法 |
CN111485850A (zh) * | 2019-01-09 | 2020-08-04 | 中国石油天然气股份有限公司 | 基于大孔道数据的油井堵水方法及装置 |
CN112924357A (zh) * | 2021-01-29 | 2021-06-08 | 西南石油大学 | 一种地层压力下致密岩石孔渗联测装置及方法 |
CN116087475A (zh) * | 2022-12-01 | 2023-05-09 | 中国石油大学(北京) | 用于分析二氧化碳驱原油组分变化的填砂管装置及其填制方法和分析方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0567239A1 (en) * | 1992-04-24 | 1993-10-27 | Core Holdings B.V. | Method of and apparatus for obtaining permeability related attributes |
CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
CN103900942A (zh) * | 2014-01-10 | 2014-07-02 | 中国石油大学(华东) | 基于高压压汞分析的储层微观孔喉参数的连续表征方法 |
CN104568694A (zh) * | 2014-12-19 | 2015-04-29 | 西南石油大学 | 一种致密岩心气-水相对渗透率的测试方法 |
CN104713812A (zh) * | 2015-04-01 | 2015-06-17 | 西南石油大学 | 一种基于岩心气测渗透率测量装置的校正方法 |
CN104833789A (zh) * | 2015-04-09 | 2015-08-12 | 中国石油大学(北京) | 利用致密砂岩微观孔隙结构确定气水关系的方法 |
CN105334149A (zh) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | 一种致密储层微观孔隙结构评价与储层分类方法 |
CN105784965A (zh) * | 2014-12-24 | 2016-07-20 | 中国石油天然气股份有限公司 | 基于模拟实验的地史演化测试方法和系统 |
CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
-
2016
- 2016-10-27 CN CN201610961286.9A patent/CN106525684B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0567239A1 (en) * | 1992-04-24 | 1993-10-27 | Core Holdings B.V. | Method of and apparatus for obtaining permeability related attributes |
CN103278436A (zh) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
CN103900942A (zh) * | 2014-01-10 | 2014-07-02 | 中国石油大学(华东) | 基于高压压汞分析的储层微观孔喉参数的连续表征方法 |
CN104568694A (zh) * | 2014-12-19 | 2015-04-29 | 西南石油大学 | 一种致密岩心气-水相对渗透率的测试方法 |
CN105784965A (zh) * | 2014-12-24 | 2016-07-20 | 中国石油天然气股份有限公司 | 基于模拟实验的地史演化测试方法和系统 |
CN104713812A (zh) * | 2015-04-01 | 2015-06-17 | 西南石油大学 | 一种基于岩心气测渗透率测量装置的校正方法 |
CN104833789A (zh) * | 2015-04-09 | 2015-08-12 | 中国石油大学(北京) | 利用致密砂岩微观孔隙结构确定气水关系的方法 |
CN105334149A (zh) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | 一种致密储层微观孔隙结构评价与储层分类方法 |
CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107144566A (zh) * | 2017-05-31 | 2017-09-08 | 中国石油大学(华东) | 利用化学‑沉积相表征细粒沉积岩非均质性的方法 |
CN109900610A (zh) * | 2017-12-07 | 2019-06-18 | 核工业北京地质研究院 | 一种多孔材料的联合测试方法 |
CN109900610B (zh) * | 2017-12-07 | 2021-07-20 | 核工业北京地质研究院 | 一种多孔材料的联合测试方法 |
CN109085104A (zh) * | 2018-07-23 | 2018-12-25 | 中国石油大学(北京) | 一种致密气藏的气体相对渗透率的确定方法 |
CN109085104B (zh) * | 2018-07-23 | 2020-06-12 | 中国石油大学(北京) | 一种致密气藏的气体相对渗透率的确定方法 |
CN109490168A (zh) * | 2018-11-30 | 2019-03-19 | 重庆科技学院 | 页岩气藏不同传质扩散机理对储层渗流能力贡献率的定量评价方法 |
CN109490168B (zh) * | 2018-11-30 | 2019-06-18 | 重庆科技学院 | 页岩气藏不同传质扩散机理对储层渗流能力贡献率的定量评价方法 |
CN111485850A (zh) * | 2019-01-09 | 2020-08-04 | 中国石油天然气股份有限公司 | 基于大孔道数据的油井堵水方法及装置 |
CN109827889A (zh) * | 2019-04-09 | 2019-05-31 | 北京艾迪佳业技术开发有限公司 | 一种克氏渗透率测量系统 |
CN111398116A (zh) * | 2020-03-04 | 2020-07-10 | 中国地质大学(武汉) | 一种利用特定方向压汞法表征页岩各向异性的方法 |
CN112924357A (zh) * | 2021-01-29 | 2021-06-08 | 西南石油大学 | 一种地层压力下致密岩石孔渗联测装置及方法 |
CN112924357B (zh) * | 2021-01-29 | 2022-02-01 | 西南石油大学 | 一种地层压力下致密岩石孔渗联测装置及方法 |
CN116087475A (zh) * | 2022-12-01 | 2023-05-09 | 中国石油大学(北京) | 用于分析二氧化碳驱原油组分变化的填砂管装置及其填制方法和分析方法 |
CN116087475B (zh) * | 2022-12-01 | 2024-07-26 | 中国石油大学(北京) | 用于分析二氧化碳驱原油组分变化的填砂管装置及其填制方法和分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106525684B (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106525684B (zh) | 一种基于孔喉结构的致密砂岩克氏渗透率的校正方法 | |
CN104568694B (zh) | 一种致密岩心气‑水相对渗透率的测试方法 | |
CN110296921B (zh) | 储层条件下稳态法页岩气体渗透率的测试装置及测试方法 | |
CN110598167B (zh) | 低渗透油藏油水相对渗透率实验数据的处理方法 | |
CN108169099A (zh) | 一种基于核磁共振的页岩气储层孔隙结构定量计算方法 | |
CN107622139B (zh) | 裂缝渗透率的计算方法 | |
CN106932323B (zh) | 一种页岩气藏气体有效孔隙度反演方法 | |
CN107525720A (zh) | 一种测试致密储层敏感性的装置及方法 | |
CN104374683A (zh) | 一种岩心孔隙压缩系数测试装置及其测试方法 | |
CN104297130B (zh) | 低渗透率储层的固有渗透率解释方法及系统 | |
CN210051673U (zh) | Grt-1型全自动储层岩石渗透率测量装置 | |
CN113075108B (zh) | 一种考虑束缚水饱和度的岩心多次应力敏感测试方法 | |
CN108133086A (zh) | 一种应力敏感储层中产水压裂气井裂缝半长反演方法 | |
Mathur et al. | Comparison of steady-state and transient methods for measuring shale permeability | |
CN109580454B (zh) | 一种用压力振荡法测试致密储层流体敏感性的方法 | |
Cui et al. | Permeability of the Montney Formation in the Western Canada Sedimentary Basin: insights from different laboratory measurements | |
CN115046897A (zh) | 基质-裂缝间非稳态传质效率及形状因子确定方法及装置 | |
CN105334293B (zh) | 多重介质中致密气的流态识别方法及其装置 | |
CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
Yanjie et al. | Lower limit of tight oil flowing porosity: Application of high-pressure mercury intrusion in the fourth Member of Cretaceous Quantou Formation in southern Songliao Basin, NE China | |
RU2542998C1 (ru) | Способ лабораторного определения анизотропии абсолютной газопроницаемости на полноразмерном керне | |
CN112129682A (zh) | 一种基于边界层动态变化的低渗透油藏驱替介质优选方法 | |
CN115559715A (zh) | 一种超高压低渗透气藏储层产水评价方法 | |
CN114577677B (zh) | 基于储层非均质性的致密油藏非线性渗流表征方法 | |
CN204575498U (zh) | 一种测量低渗透岩心的气测渗透率装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |