CN110987761B - 一种定量计算岩心孔隙连通性的方法 - Google Patents

一种定量计算岩心孔隙连通性的方法 Download PDF

Info

Publication number
CN110987761B
CN110987761B CN201911289477.5A CN201911289477A CN110987761B CN 110987761 B CN110987761 B CN 110987761B CN 201911289477 A CN201911289477 A CN 201911289477A CN 110987761 B CN110987761 B CN 110987761B
Authority
CN
China
Prior art keywords
nuclear magnetic
pore
centrifugation
core
magnetic resonance
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
CN201911289477.5A
Other languages
English (en)
Other versions
CN110987761A (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.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and 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 Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201911289477.5A priority Critical patent/CN110987761B/zh
Publication of CN110987761A publication Critical patent/CN110987761A/zh
Application granted granted Critical
Publication of CN110987761B publication Critical patent/CN110987761B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Abstract

本发明提供一种可以定量对不同孔径的孔隙之间的连通关系进行测量的计算岩心孔隙连通性的方法,包括以下步骤:S1:饱和水岩心的核磁共振实验;S2:对岩心进行多次离心,对每一次离心后的岩心均开展核磁共振实验;S3:计算离心力对应的孔隙半径;S4:定量计算孔隙连通性。本发明采用核磁共振与多次离心相结合的方法解决了孔隙连通性难以定量计算的难题。

Description

一种定量计算岩心孔隙连通性的方法
技术领域
本发明属于油气藏开发领域,具体涉及一种定量计算岩心孔隙连通性的方法。
背景技术
孔隙连通性是油气藏开发中的重要参数。它反映了油气藏的储集能力和渗透能力,影响了油气的流动性、油气井产能、油气藏采收率等,对生产制度的制定和开发方案的编制具有重要的指导意义。
分析孔隙连通性的常规实验方法有:CT扫描、扫描电镜、压汞法、核磁共振等,但CT扫描和扫描电镜只能对单个或者几个孔隙的通性的好坏进行微观、定性地分析,难以对进行宏观、定量的计算与评价。并且,如果视域相对较大,则难以获得微小孔隙的图像,如果想要获得微小孔隙的图像,则难以将微小视域的图像模型外推到整个储集层。压汞法、核磁共振法能够获取岩心的孔径分布,但难以对不同孔径的孔隙之间的连通关系进行测量。
发明内容
为了解决上述问题,本发明提供一种可以定量对不同孔径的孔隙之间的连通关系进行测量的计算岩心孔隙连通性的方法,包括以下步骤:
S1:饱和水岩心的核磁共振实验;
S2:对岩心进行多次离心,对每一次离心后的岩心均开展核磁共振实验;
S3:计算离心力对应的孔隙半径;
S4:定量计算孔隙连通性。
进一步的,所述饱和水岩心的核磁共振实验包括:
将岩心洗净烘干并抽真空后在饱和水状态下开展核磁共振实验,测试并绘制核磁T2谱曲线,计算核磁T2谱曲线与T2时间轴围成的面积为S。
所述岩心多次离心后的核磁共振实验包括:
以离心力F1对饱和水的岩心进行第一次离心,再对第一次离心后岩心开展核磁共振实验,测试核磁T2谱,并计算出T2曲线与T2时间轴围成的面积为S1
以大于离心力F1的离心力F2进行第二次离心,离心结束后,对第二次离心后的岩心开展核磁共振实验,测试核磁T2谱,并计算出T2曲线与T2时间轴围成的面积为S2
岩心多次离心后的核磁共振实验包括:
对岩心开展n次离心,第n次离心的离心力为Fn,并且,Fn>Fn-1,对离心n次后的岩心开展核磁共振实验,测试核磁T2谱,并计算出核磁T2谱与T2时间轴围成的面积为Sn
进一步的,所述计算离心力对应的孔隙半径包括:
采用以下公式获取孔隙半径:
Figure GDA0003413722380000021
式中,F为离心力,σ为气水界面张力,θ为润湿角,r为孔隙半径。
所述定量计算孔隙连通性包括:
在核磁共振和n次离心实验和m次离心实验后,第m次离心对应的离心力为Fm,第m次离心对应的孔隙半径为rm,第m次离心对应核磁T2谱与T2时间轴围成的面积为Sm,第n次离心对应的离心力为Fn,第n次离心对应的孔隙半径为rn,第n次离心对应核磁T2谱与T2时间轴围成的面积为Sn,n大于m,采用以下公式获取rn的连通关系:
比rn大的孔隙与rn相连通的占比为:
Figure GDA0003413722380000031
孔径在rn~rm之间的孔隙与rn连通的占比为:
Figure GDA0003413722380000032
比rn小的孔隙与rn相连通的占比为:
Figure GDA0003413722380000033
本发明的有益效果是:本发明通过对饱和水的岩心进行离心,对流体在孔隙与孔隙间的流动来反映孔隙与孔隙之间的连通情况。每次离心后开展核磁共振实验,则能够对对应尺寸孔隙间的连通情况进行定量计算。采用核磁共振与多次离心相结合的方法解决了孔隙连通性难以定量计算的难题。
附图说明
图1本发明流程图。
图2岩心多次离心后测试的核磁T2谱示意图。
具体实施方式
如图1所示,本发明提供一种可以定量对不同孔径的孔隙之间的连通关系进行测量的计算岩心孔隙连通性的方法,包括以下步骤:
S1:饱和水岩心的核磁共振实验;
S2:对岩心进行多次离心,对每一次离心后的岩心均开展核磁共振实验;
S3:计算离心力对应的孔隙半径;
S4:定量计算孔隙连通性。
下面对各步骤进行具体说明。
饱和水岩心的核磁共振实验步骤包括,将岩心洗净烘干、抽真空后饱和水。在饱和水状态下,开展核磁共振实验,测试并绘制核磁T2谱曲线,并计算核磁T2谱曲线与X轴(T2时间轴)围成的面积为S。
岩心多次离心后的核磁共振实验包括,以一个较小的离心力F1对饱和水的岩心进行离心(记为离心1),再对离心1后的岩心开展核磁共振实验,测试核磁T2谱,并计算出T2曲线与T2时间轴围成的面积为S1;随后对离心1后的岩心以离心力F2进行离心(记为离心2),并且F2>F1,离心结束后,对离心2后的岩心开展核磁共振实验,测试核磁T2谱,并计算出T2曲线与T2时间轴围成的面积为S2。根据需要,可对岩心开展n次离心,第n次离心(记为离心n)的离心力为Fn,并且,Fn>Fn-1,对离心n次后的岩心开展核磁共振实验,测试核磁T2谱,并计算出核磁T2谱与T2时间轴围成的面积为Sn
计算离心力对应的孔隙半径步骤包括,在离心时,当离心力大于或等于某孔隙对应的毛管压力时,该孔隙内的流体才可能被离心出岩心。因此,离心力等于某孔径对应的毛管压力是该孔径内流体被离心出岩心的临界条件。根据毛管压力与孔径的关系式,可以推导离心力与孔隙半径的关系式:
Figure GDA0003413722380000041
式中,F为离心力,pc为毛管压力,σ为气水界面张力,θ为润湿角,r为孔隙半径。
因此,离心力F对应的孔隙半径r可表示为:
Figure GDA0003413722380000042
采用式(2)计算F1、F2、……、Fn对应的孔径分别为:r1、r2、……、rn。其中:F1<F2<……<Fn,r1>r2>……>rn
定量计算孔隙连通性步骤包括,计算每次离心后,核磁T2谱曲线与X轴(T2时间轴)围成的面积。以4次离心为例(如图1所示),分别计算S、S1、S2、S3、S4。采用式(2)计算S、S1、S2、S3、S4对应的孔隙半径分别为r1、r2、r3、r4
离心时,孔隙p(孔径为rp)内的流体被离心出岩心需要同时满足以下两个条件:①离心力大于或等于孔隙p内的毛管压力;②孔隙p与孔径大于rp的孔隙连通。根据这一原理,通过离心后流体滞留在岩心内的比例(通过S反映),即可计算各孔隙间的连通性。
因此,r1~r4的孔隙连通性可通过下列方式计算:
(1)r1的连通关系:
比r1大的孔隙与r1连通的占比为:
Figure GDA0003413722380000051
比r1小的孔隙与r1连通的占比为:
Figure GDA0003413722380000052
(2)r2的连通关系:
比r2大的孔隙与r2连通的总占比为:
Figure GDA0003413722380000053
其中,孔径在r2~r1之间的孔隙与r2连通的占比为:
Figure GDA0003413722380000054
比r2小的孔隙与r2连通的占比为:
Figure GDA0003413722380000055
(3)r3的连通关系:
比r3大的孔隙与r3连通的占比为:
Figure GDA0003413722380000056
其中,孔径在r3~r2之间的孔隙与r3连通的占比为:
Figure GDA0003413722380000057
比r3小的孔隙与r3连通的占比为:
Figure GDA0003413722380000058
(4)r4的连通关系:
比r4大的孔隙与r4连通的占比为:
Figure GDA0003413722380000059
其中,孔径在r4~r3之间的孔隙与r4连通的占比为:
Figure GDA00034137223800000510
比r4小的孔隙与r4连通的占比为:
Figure GDA00034137223800000511
因此,对一般规律进行总结。在核磁共振和多次离心实验后,离心m对应的离心力为Fm,对应的孔隙半径为rm,其核磁T2谱与T2时间轴围成的面积为Sm。离心n(n>m),对应的离心力为Fn(Fn>Fm),对应的孔隙半径为rn(rn<rm),其核磁T2谱与T2时间轴围成的面积为Sn(Sn<Sm)。
则rn的连通关系为:
比rn大的孔隙与rn相连通的占比为:
Figure GDA0003413722380000061
其中,孔径在rn~rm之间的孔隙与rn连通的占比为:
Figure GDA0003413722380000062
比rn小的孔隙与rn相连通的占比为:
Figure GDA0003413722380000063
因此,孔隙连通性可由上述通式计算得出。
虽然已经参考优选实施例对本发明进行了描述,但在不脱离本发明的范围的情况下,可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是,只要不存在结构冲突,各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例,而是包括落入权利要求的范围内的所有技术方案。

Claims (3)

1.一种定量计算岩心孔隙连通性的方法,包括以下步骤:
S1:饱和水岩心的核磁共振实验;
S2:对岩心进行多次离心,对每一次离心后的岩心均开展核磁共振实验;包括:
以离心力F1对饱和水的岩心进行第一次离心,再对第一次离心后岩心开展核磁共振实验,测试核磁T2谱,并计算出T2曲线与T2时间轴围成的面积为S1
以大于离心力F1的离心力F2进行第二次离心,离心结束后,对第二次离心后的岩心开展核磁共振实验,测试核磁T2谱,并计算出T2曲线与T2时间轴围成的面积为S2
对岩心开展n次离心,第n次离心的离心力为Fn,并且,Fn>Fn-1,对离心n次后的岩心开展核磁共振实验,测试核磁T2谱,并计算出核磁T2谱与T2时间轴围成的面积为Sn
S3:计算离心力对应的孔隙半径;
S4:定量计算孔隙连通性,包括:
在核磁共振和n次离心实验和m次离心实验后,第m次离心对应的离心力为Fm,第m次离心对应的孔隙半径为rm,第m次离心对应核磁T2谱与T2时间轴围成的面积为Sm,第n次离心对应的离心力为Fn,第n次离心对应的孔隙半径为rn,第n次离心对应核磁T2谱与T2时间轴围成的面积为Sn,n大于m,采用以下公式获取rn的连通关系:
比rn大的孔隙与rn相连通的占比为:
Figure FDA0003413722370000011
孔径在rn~rm之间的孔隙与rn连通的占比为:
Figure FDA0003413722370000012
比rn小的孔隙与rn相连通的占比为:
Figure FDA0003413722370000021
2.如权利要求1所述的一种定量计算岩心孔隙连通性的方法,其特征在于,所述饱和水岩心的核磁共振实验包括:
将岩心洗净烘干并抽真空后在饱和水状态下开展核磁共振实验,测试并绘制核磁T2谱曲线,计算核磁T2谱曲线与T2时间轴围成的面积为S。
3.如权利要求1所述的一种定量计算岩心孔隙连通性的方法,其特征在于,所述计算离心力对应的孔隙半径包括:
采用以下公式获取孔隙半径:
Figure FDA0003413722370000022
式中,F为离心力,σ为气水界面张力,θ为润湿角,r为孔隙半径。
CN201911289477.5A 2019-12-13 2019-12-13 一种定量计算岩心孔隙连通性的方法 Active CN110987761B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911289477.5A CN110987761B (zh) 2019-12-13 2019-12-13 一种定量计算岩心孔隙连通性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911289477.5A CN110987761B (zh) 2019-12-13 2019-12-13 一种定量计算岩心孔隙连通性的方法

Publications (2)

Publication Number Publication Date
CN110987761A CN110987761A (zh) 2020-04-10
CN110987761B true CN110987761B (zh) 2022-03-18

Family

ID=70094040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911289477.5A Active CN110987761B (zh) 2019-12-13 2019-12-13 一种定量计算岩心孔隙连通性的方法

Country Status (1)

Country Link
CN (1) CN110987761B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112964616B (zh) * 2021-02-05 2022-05-10 福州大学 一种确定水泥基材料不同尺度孔隙连通性的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104697915A (zh) * 2015-03-20 2015-06-10 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 一种页岩微观孔隙大小及流体分布的分析方法
CN104819923A (zh) * 2015-05-17 2015-08-05 西南石油大学 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法
CN105866009A (zh) * 2016-05-30 2016-08-17 中国石油大学(北京) 一种计算致密油储层有效孔隙度的方法及装置
CN106249306A (zh) * 2016-10-12 2016-12-21 贵州大学 基于核磁共振的页岩孔隙结构检测方法
CN106442600A (zh) * 2016-11-23 2017-02-22 中国石油大学(华东) 页岩束缚水含量确定方法
CN109030311A (zh) * 2018-07-16 2018-12-18 西南石油大学 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2489005B (en) * 2011-03-15 2015-11-25 Schlumberger Holdings Examination of porosity by NMR and intrusion porosimetry
US10156531B2 (en) * 2015-10-27 2018-12-18 Saudi Arabian Oil Company Measuring connectivity between different pore types in porous media

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104697915A (zh) * 2015-03-20 2015-06-10 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 一种页岩微观孔隙大小及流体分布的分析方法
CN104819923A (zh) * 2015-05-17 2015-08-05 西南石油大学 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法
CN105866009A (zh) * 2016-05-30 2016-08-17 中国石油大学(北京) 一种计算致密油储层有效孔隙度的方法及装置
CN106249306A (zh) * 2016-10-12 2016-12-21 贵州大学 基于核磁共振的页岩孔隙结构检测方法
CN106442600A (zh) * 2016-11-23 2017-02-22 中国石油大学(华东) 页岩束缚水含量确定方法
CN109030311A (zh) * 2018-07-16 2018-12-18 西南石油大学 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Comparative analysis on watdifferent reservoir rocks by nuclear magnetic resonanceer movability in pores of;Song Li 等;《ENERGY EXPLORATION & EXPLOITATION》;20151231;第33卷(第5期);第689-705页 *

Also Published As

Publication number Publication date
CN110987761A (zh) 2020-04-10

Similar Documents

Publication Publication Date Title
CN106979917B (zh) 一种页岩气储层孔隙结构的表征方法以及页岩气储层的评价方法
CN108507920B (zh) 一种油水两相相对渗透率的模拟方法
CN110987761B (zh) 一种定量计算岩心孔隙连通性的方法
CN111128209B (zh) 一种基于混合掩蔽学习目标的语音增强方法
CN1946920A (zh) 用于储层流体表征的井下流体分析用方法和装置
CN108982287B (zh) 一种确定页岩气储层游离吸附比的方法和装置
CN109187614B (zh) 基于核磁共振和质谱的代谢组学数据融合方法及其应用
CA2652336A1 (en) Methods and apparatus for analysis of downhole asphaltene gradients and applications thereof
CN106525688A (zh) 一种饱和页岩孔隙流体分离及饱和度计算的实验方法
CN110849785B (zh) 一种利用多次压汞实验表征页岩孔隙连通性的方法
CN107843611A (zh) 低渗透砂岩储层可动流体饱和度核磁共振参数表征新方法
CN114577677B (zh) 基于储层非均质性的致密油藏非线性渗流表征方法
CN111860863B (zh) 一种基于协方差矩阵的主动学习方法
CN107133879A (zh) 一种筛选相似油田的方法
CN116383573B (zh) 一种基于多区相变传质渗流耦合的凝析气产能评价方法
CN115755195A (zh) 一种储层岩石模拟生产条件下的润湿性表征方法
CN106285658B (zh) 一种毛管力曲线校正的方法及装置
CN113253354A (zh) 一种核磁t2谱截止值的确定方法、装置、介质及电子设备
CN109459402B (zh) 一种光谱法水质cod参数预测方法
CN116183458A (zh) 页岩油有效孔隙度测定方法
US20220058303A1 (en) Method and system for quantitatively evaluating kerogen swelling oil in shale
CN108572129B (zh) 致密油有效储层孔隙阈值界定方法及系统
CN115201244A (zh) 致密低渗气藏气水两相微观渗流评价方法
CN113449261B (zh) 页岩气储层开采性分析方法及装置
CN113176186A (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