CN111721683A - 岩样微观孔隙结构定量检测方法及装置 - Google Patents

岩样微观孔隙结构定量检测方法及装置 Download PDF

Info

Publication number
CN111721683A
CN111721683A CN201910206682.4A CN201910206682A CN111721683A CN 111721683 A CN111721683 A CN 111721683A CN 201910206682 A CN201910206682 A CN 201910206682A CN 111721683 A CN111721683 A CN 111721683A
Authority
CN
China
Prior art keywords
curve
rock sample
magnetic resonance
nuclear magnetic
mercury intrusion
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.)
Pending
Application number
CN201910206682.4A
Other languages
English (en)
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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201910206682.4A priority Critical patent/CN111721683A/zh
Publication of CN111721683A publication Critical patent/CN111721683A/zh
Pending legal-status Critical Current

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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Dispersion Chemistry (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种岩样微观孔隙结构定量检测方法及装置,该方法包括:获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线;根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系;获取检测到的所述储层的待测岩样的核磁共振T2谱曲线;根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。本发明解决了现有的微观孔隙结构测试方法无法实现对复杂岩样的微观孔隙结构进行定量的检测的技术问题。

Description

岩样微观孔隙结构定量检测方法及装置
技术领域
本发明涉及岩样微观孔隙结构分析技术,具体而言,涉及一种岩样微观孔隙结构定量检测方法及装置。
背景技术
储层微观孔隙结构决定着油气藏渗流能力和产油气能力,尤其对于低渗致密油气藏来说,微观孔隙结构对油气藏开发影响非常大,决定着储层开发潜力,是开展储层评价重要参数。目前,针对储层微观孔隙结构研究主要有压汞法(恒速压汞、高压压汞)、核磁共振法、CT扫描法和铸体薄片法等。其中压汞法用于检测出岩样的压汞曲线,进而根据压汞曲线可定量描述岩心微观孔隙结构,但压汞法只适合尺寸规则的标准岩样,例如直径2.5cm,厚度0.5cm的标准岩样,不适合形状复杂、不规则的复杂岩样(对于裂缝孔隙型碳酸盐岩、砂砾岩和疏松砂岩等,很难取得规则岩样,实验岩样难以制作)。核磁共振法对岩样形状要求较少,但目前只能定性描述岩心微观孔隙结构,将孔隙分为细孔隙、小孔隙、中孔隙、大孔隙和洞。CT扫描和铸体薄片直观展示岩心微观孔隙结构,但缺乏识别精度和数字重建困难。现有的微观孔隙结构测试方法无法实现对复杂岩样的微观孔隙结构进行定量的检测。
发明内容
本发明的主要目的在于提供一种岩样微观孔隙结构定量检测方法及装置,以解决现有的微观孔隙结构测试方法无法实现对复杂岩样的微观孔隙结构进行定量的检测的技术问题。
为了实现上述目的,根据本发明的一个方面,提供了一种岩样微观孔隙结构定量检测方法,该方法包括:获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线;
根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系;
获取检测到的所述储层的待测岩样的核磁共振T2谱曲线;
根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
可选的,所述获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线,包括:
获取检测到的所述储层的第一标准岩样的核磁共振T2谱曲线;以及
获取检测到的所述储层的第二标准岩样的压汞曲线。
可选的,所述获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线,包括:
分别获取所述储层的多块标准岩样的核磁共振T2谱曲线和压汞曲线。
可选的,所述根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系,包括:
根据每块标准岩样的核磁共振T2谱曲线和压汞曲线确定每块标准岩样对应的核磁共振T2谱曲线与压汞曲线的对应关系;
将确定的多个对应关系取平均值,得到所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
可选的,所述储层的核磁共振T2谱曲线与压汞曲线的对应关系包括:T2弛豫时间与喉道半径的对应关系。
为了实现上述目的,根据本发明的另一方面,提供了一种岩样微观孔隙结构定量检测装置,该装置包括:标准岩样曲线获取单元,用于获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线;
曲线对应关系确定单元,用于根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系;
待测岩样曲线获取单元,用于获取检测到的所述储层的待测岩样的核磁共振T2谱曲线;
微观孔隙结构定量确定单元,用于根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
可选的,所述标准岩样曲线获取单元包括:
第一曲线获取模块,用于获取检测到的所述储层的第一标准岩样的核磁共振T2谱曲线;
第二曲线获取模块,用于获取检测到的所述储层的第二标准岩样的压汞曲线。
可选的,所述标准岩样曲线获取单元包括:
多岩样曲线获取模块,用于分别获取所述储层的多块标准岩样的核磁共振T2谱曲线和压汞曲线。
可选的,所述曲线对应关系确定单元包括:
对应关系确定模块,用于根据每块标准岩样的核磁共振T2谱曲线和压汞曲线确定每块标准岩样对应的核磁共振T2谱曲线与压汞曲线的对应关系;
计算模块,用于将确定的多个对应关系取平均值,得到所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
可选的,所述储层的核磁共振T2谱曲线与压汞曲线的对应关系包括:T2弛豫时间与喉道半径的对应关系。
为了实现上述目的,根据本发明的另一方面,提供了另一种岩样微观孔隙结构定量检测装置,该装置包括:ISCO泵、中间容器、核磁共振仪、高压压汞仪以及处理器;
所述中间容器,用于盛装地层水、标准岩样以及与所述标准岩样处于同一储层的待测岩样;
所述ISCO泵,用于给所述中间容器加压,使所述标准岩样和所述待测岩样含水饱和;
所述核磁共振仪,用于分别检测含水饱和的所述标准岩样及检测含水饱和的所述待测岩样的核磁共振T2谱曲线;
高压压汞仪,用于检测含水饱和的所述标准岩样的压汞曲线;
所述处理器,用于根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定储层的核磁共振T2谱曲线与压汞曲线的对应关系;根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
为了实现上述目的,根据本发明的另一方面,还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述岩样微观孔隙结构定量检测方法中的步骤。
为了实现上述目的,根据本发明的另一方面,还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序在计算机处理器中执行时实现上述岩样微观孔隙结构定量检测方法中的步骤。
本发明的有益效果为:本发明提供的岩样微观孔隙结构定量检测方法在对复杂岩样进行微观孔隙结构定量测试时只需进行核磁共振T2谱测试,具有实验精度高,实验岩样不需特殊处理,操作简单等优点。此外,本发明的岩样微观孔隙结构定量检测方法实验成本低且快捷,测试一块复杂岩样核磁共振T2谱曲线所需时间不到5分钟,极大地提高了复杂岩样微观孔隙结构定量检查效率,适合现场、大批量微观孔隙结构测试,为开展储层评价提供基础参数与技术支撑。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1是本发明实施例岩样微观孔隙结构定量检测方法的流程图;
图2是本发明实施例获取检测的标准岩样的曲线的流程图;
图3是本发明实施例确定曲线对应关系的方法的流程图;
图4是本发明实施例岩样微观孔隙结构定量检测装置的结构框图;
图5是本发明实施例标准岩样曲线获取单元的结构框图;
图6是本发明实施例曲线对应关系确定单元的结构框图;
图7是本发明实施例标准岩样核磁共振T2谱曲线与压汞曲线对比图;
图8是本发明实施例核磁共振法与压汞法测试结果对比图;
图9是本发明实施例计算机设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
本发明的目的是提供一种在实验室中数据测量精确、易于操作的、快速的复杂岩样微观孔隙结构定量测试方法与装置,以测试分析储层微观孔隙结构,为油气藏储层评价提供技术支撑。
图1是本发明实施例岩样微观孔隙结构定量检测方法的流程图,如图1所示,本发明实施例的岩样微观孔隙结构定量检测方法包括步骤S101至步骤S104。
步骤S101,获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线。
在本发明的实施例中,岩样的核磁共振T2谱曲线和压汞曲线可以在实验室中通过检测设备得出,该检测设备包括:ISCO泵、中间容器、核磁共振仪和高压压汞仪。其中,中间容器主要用于盛装地层水、复杂岩样和标准岩样;ISCO泵主要用于给盛装地层水、复杂岩样和标准岩样的中间容器加压(压力值ISCO泵上显示),让标准岩样和复杂岩样充分饱和;核磁共振仪主要用于测量饱和地层水后的标准岩样和复杂岩样核磁共振T2谱,获得标准岩样和复杂岩样核磁共振T2谱曲线;高压压汞仪主要用于对标准岩样进行压汞测试,获取标准岩样的压汞曲线。
图7是本发明实施例标准岩样核磁共振T2谱曲线与压汞曲线对比图,在图7中同时展示了检测到的砂岩的标准岩样的核磁共振T2谱曲线和压汞曲线。在本发明的实施例中,该核磁共振T2谱曲线可以为核磁共振T2弛豫时间累计曲线,该核磁共振T2弛豫时间累计曲线的横坐标为T2弛豫时间(单位ms),纵坐标为累计幅度。在本发明的实施例中,该压汞曲线可以为压汞喉道半径累计曲线,该压汞喉道半径累计曲线的横坐标为喉道半径(单位μm),纵坐标为累计幅度。从图7中可以看出砂岩核磁共振T2弛豫时间累计曲线与压汞喉道半径累计曲线存在很好的一致性。
步骤S102,根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
在本发明的实施例中,可以将标准岩样的核磁共振T2谱曲线与压汞曲线进行对比,并建立该两种曲线的换算关系,进而根据换算关系就可以实现两种曲线之间相互转换。在本发明实施例中,该核磁共振T2谱曲线可以为核磁共振T2弛豫时间累计曲线,该核磁共振T2弛豫时间累计曲线的横坐标为T2弛豫时间(单位ms),纵坐标为累计幅度;该压汞曲线可以为压汞喉道半径累计曲线,该压汞喉道半径累计曲线的横坐标为喉道半径(单位μm),纵坐标为累计幅度。由于核磁共振T2弛豫时间累计曲线和压汞喉道半径累计曲线的纵坐标均为累计幅度,因此核磁共振T2谱曲线与压汞曲线的对应关系也可以为T2弛豫时间与喉道半径的对应关系。
在本发明的实施例中,储层的核磁共振T2谱曲线与压汞曲线的对应关系可以由计算机对两曲线进行分析对比得出。在本发明的实施例中,储层的核磁共振T2谱曲线与压汞曲线的对应关系也可以通过分析软件得出,例如Matlab等。
在如图7所示的实施例中,通过分析,砂岩T2弛豫时间(单位ms)与喉道半径R(单位μm)满足如下关系:
Figure BDA0001999200010000061
在得到T2弛豫时间与喉道半径的对应关系后,就可以根据该对应关系实现核磁共振T2谱曲线与压汞曲线之间相互转换。
步骤S103,获取检测到的所述储层的待测岩样的核磁共振T2谱曲线。
在本发明的实施例中,该待测岩样为与上述标准岩样位于同一储层的复杂岩样。该复杂岩样的核磁共振T2谱曲线可以由以下步骤检测得出:通过ISCO泵给盛装地层水和复杂岩样的中间容器加压(压力值ISCO泵上显示),让复杂岩样充分饱和;通过核磁共振仪测量饱和地层水后的复杂岩样的核磁共振T2谱,获得复杂岩样核磁共振T2谱曲线。
步骤S104,根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
在本发明实施例中,在得到复杂岩样的核磁共振T2谱曲线后,可以根据在上述步骤S102中得出的储层的核磁共振T2谱曲线与压汞曲线的对应关系,进而通过数值反演得到复杂岩样对应的压汞曲线。在本发明的实施例中,在得到复杂岩样对应的压汞曲线后,就可以根据现有技术的方法,确定该复杂岩样定量的微观孔隙结构。
图8是本发明实施例核磁共振法与压汞法测试结果对比图,图8中的压汞法曲线为采用本发明方法得出的复杂岩样的压汞曲线。如图8所示,通过对比本发明方法的测试结果与核磁共振法的测试结果,可以看出两种方法测试结构高度一致,相关性好,因此也可以利用核磁共振快速、定量测试岩样微观孔隙结构,获取喉道分布曲线,避免压汞,既环保,又提高测试效率,尤其对于复杂不规则岩样微观孔隙结构,形状要求无法满足测试要求,本文所述复杂岩样微观孔隙结构快速、定量测试方法可满足测试要求,使之微观孔隙结构定量测试成为可能。
由以上描述可以看出,本发明提供的岩样微观孔隙结构定量检测方法在对复杂岩样进行微观孔隙结构定量测试时只需进行核磁共振T2谱测试,具有实验精度高,实验岩样不需特殊处理,操作简单等优点。此外,本发明的岩样微观孔隙结构定量检测方法实验成本低且快捷,测试一块复杂岩样核磁共振T2谱曲线所需时间不到5分钟,极大地提高了复杂岩样微观孔隙结构定量检查效率,适合现场、大批量微观孔隙结构测试,为开展储层评价提供基础参数与技术支撑。
图2是本发明实施例获取检测的标准岩样的曲线的流程图,如图2所示,本发明实施例的获取检测的标准岩样的曲线的方法包括步骤S201和步骤S202。
步骤S201,获取检测到的所述储层的第一标准岩样的核磁共振T2谱曲线。
步骤S201,获取检测到的所述储层的第二标准岩样的压汞曲线。
在本发明的实施例中,获取的标准岩样的核磁共振T2谱曲线和压汞曲线,可以为检测得到的同一块标准岩样的核磁共振T2谱曲线和压汞曲线。在本发明的其他实施例中,为了便于实验也可以采用两块标准岩样分别检测的方法,一块标准岩样用于测量核磁共振T2谱曲线,另一块用于测量压汞曲线。当然该两块标准岩样为同一储层的标准岩样,此外该两块标准岩样的尺寸可以为同样大小或不同大小。
图3是本发明实施例确定曲线对应关系的方法的流程图,如图3所示,本发明实施例的确定曲线对应关系的方法包括步骤S301和步骤S303。
步骤S301,分别获取所述储层的多块标准岩样的核磁共振T2谱曲线和压汞曲线。
步骤S302,根据每块标准岩样的核磁共振T2谱曲线和压汞曲线确定每块标准岩样对应的核磁共振T2谱曲线与压汞曲线的对应关系。
步骤S303,将确定的多个对应关系取平均值,得到所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
在本发明的实施例中,根据标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系时,为了提高该对应关系的准确性,也可以设置多块标准岩样,分别检测出每块标准岩样的核磁共振T2谱曲线和压汞曲线,进而根据每块标准岩样的两种曲线来确定储层的核磁共振T2谱曲线与压汞曲线的对应关系。
在本申请实施例中,上述根据多块标准岩样的两种曲线来确定储层的核磁共振T2谱曲线与压汞曲线的对应关系具体可以为,分别根据每块标准岩样的两种曲线求出该块标准岩样对应的所述对应关系,进而根据每块标准岩样求出的对应关系确定出最终较为准确的储层的核磁共振T2谱曲线与压汞曲线的对应关系。在本发明的实施例中所述对应关系可以为T2弛豫时间与喉道半径之间的对应系数,入上述步骤S102中的公式。根据每块标准岩样求出的对应关系确定出最终较为准确的对应关系可以为,根据每块标准岩样求出的T2弛豫时间与喉道半径之间的对应系数,进而通过求平均值得出最终较为准确的储层的核磁共振T2谱曲线与压汞曲线的对应关系。
本发明还公开了一种在实验室环境中应用本发明的岩样微观孔隙结构定量检测方法定量确定复杂岩样微观孔隙结构的实验方法。该实验的步骤具体包括:
1、根据高压压汞仪要求,将岩心样本(直径2.5cm,长度5cm左右)切割成两块标准岩样,一块直径2.5cm、长度0.5cm,另一块直径2.5cm、长度4cm左右的岩样(切割过程中会有部分长度损失),其中直径2.5cm,长度0.5cm的标准岩样用于进行高压压汞测试,获取压汞曲线;
2、将切割后的直径2.5cm,长度4cm左右的标准岩样进行真空、利用ISCO泵加压饱和地层水,并进行核磁共振测试,获取核磁共振T2谱曲线;
3、在电脑上对压汞测试获得的压汞曲线(喉道半径累计分布曲线)与核磁共振获得的核磁共振T2谱曲线(T2谱累计分布曲线)进行比对;建立压汞曲线与核磁共振T2谱曲线的换算关系式;
4、将复杂岩样进行抽真空、利用ISCO泵加压饱和地层水,并进行核磁共振测试,获取复杂岩样核磁共振T2谱曲线;
5、利用步骤3获得的换算关系式和步骤4获得的复杂岩样核磁共振T2谱曲线,数值计算确定复杂岩样压汞曲线,定量确定复杂岩样微观孔隙结构,结束实验。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
基于同一发明构思,本发明实施例还提供了一种岩样微观孔隙结构定量检测装置,可以用于实现上述实施例所描述的岩样微观孔隙结构定量检测方法,如下面的实施例所述。由于岩样微观孔隙结构定量检测装置解决问题的原理与岩样微观孔隙结构定量检测方法相似,因此岩样微观孔隙结构定量检测装置的实施例可以参见岩样微观孔隙结构定量检测方法的实施例,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图4是本发明实施例岩样微观孔隙结构定量检测装置的结构框图,如图4所示,本发明实施例岩样微观孔隙结构定量检测装置包括:标准岩样曲线获取单元1、曲线对应关系确定单元2、待测岩样曲线获取单元3以及微观孔隙结构定量确定单元4。
标准岩样曲线获取单元1,用于获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线。
在本发明的实施例中,岩样的核磁共振T2谱曲线和压汞曲线可以在实验室中通过检测设备得出,该检测设备包括:ISCO泵、中间容器、核磁共振仪和高压压汞仪。其中,中间容器主要用于盛装地层水、复杂岩样和标准岩样;ISCO泵主要用于给盛装地层水、复杂岩样和标准岩样的中间容器加压(压力值ISCO泵上显示),让标准岩样和复杂岩样充分饱和;核磁共振仪主要用于测量饱和地层水后的标准岩样和复杂岩样核磁共振T2谱,获得标准岩样和复杂岩样核磁共振T2谱曲线;高压压汞仪主要用于对标准岩样进行压汞测试,获取标准岩样的压汞曲线。
曲线对应关系确定单元2,根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
在本发明的实施例中,可以将标准岩样的核磁共振T2谱曲线与压汞曲线进行对比,并建立该两种曲线的换算关系,进而根据换算关系就可以实现两种曲线之间相互转换。在本发明实施例中,该核磁共振T2谱曲线可以为核磁共振T2弛豫时间累计曲线,该核磁共振T2弛豫时间累计曲线的横坐标为T2弛豫时间(单位ms),纵坐标为累计幅度;该压汞曲线可以为压汞喉道半径累计曲线,该压汞喉道半径累计曲线的横坐标为喉道半径(单位μm),纵坐标为累计幅度。由于核磁共振T2弛豫时间累计曲线和压汞喉道半径累计曲线的纵坐标均为累计幅度,因此核磁共振T2谱曲线与压汞曲线的对应关系也可以为T2弛豫时间与喉道半径的对应关系。
在本发明的实施例中,储层的核磁共振T2谱曲线与压汞曲线的对应关系可以由计算机对两曲线进行分析对比得出。在本发明的实施例中,储层的核磁共振T2谱曲线与压汞曲线的对应关系也可以通过分析软件得出,例如Matlab等。
在如图7所示的实施例中,通过分析,砂岩T2弛豫时间(单位ms)与喉道半径R(单位μm)满足如下关系:
Figure BDA0001999200010000101
在得到T2弛豫时间与喉道半径的对应关系后,就可以根据该对应关系实现核磁共振T2谱曲线与压汞曲线之间相互转换。
待测岩样曲线获取单元3,用于获取检测到的所述储层的复杂岩样的核磁共振T2谱曲线。
在本发明的实施例中,该待测岩样为与上述标准岩样位于同一储层的复杂岩样。复杂岩样的核磁共振T2谱曲线可以由以下步骤检测得出:通过ISCO泵给盛装地层水和复杂岩样的中间容器加压(压力值ISCO泵上显示),让复杂岩样充分饱和;通过核磁共振仪测量饱和地层水后的复杂岩样的核磁共振T2谱,获得复杂岩样核磁共振T2谱曲线。
微观孔隙结构定量确定单元4,用于根据所述复杂岩样的核磁共振T2谱曲线以及所述对应关系确定出所述复杂岩样的压汞曲线,以根据该压汞曲线定量确定所述复杂岩样的微观孔隙结构。
在本发明实施例中,在得到复杂岩样的核磁共振T2谱曲线后,可以根据得出的储层的核磁共振T2谱曲线与压汞曲线的对应关系,进而通过数值反演得到复杂岩样对应的压汞曲线。在本发明的实施例中,在得到复杂岩样对应的压汞曲线后,就可以根据现有技术的方法,确定该复杂岩样定量的微观孔隙结构。
图5是本发明实施例标准岩样曲线获取单元的结构框图,如图5所示,本发明实施例的标准岩样曲线获取单元包括:第一曲线获取模块101和第二曲线获取模块102。
第一曲线获取模块101,用于获取检测到的所述储层的第一标准岩样的核磁共振T2谱曲线。
第二曲线获取模块102,用于获取检测到的所述储层的第二标准岩样的压汞曲线。
在本发明的实施例中,获取的标准岩样的核磁共振T2谱曲线和压汞曲线,可以为检测得到的同一块标准岩样的核磁共振T2谱曲线和压汞曲线。在本发明的其他实施例中,为了便于实验也可以采用两块标准岩样分别检测的方法,一块标准岩样用于测量核磁共振T2谱曲线,另一块用于测量压汞曲线。当然该两块标准岩样为同一储层的标准岩样,此外该两块标准岩样的尺寸可以为同样大小或不同大小。
如图5所示,本发明实施例的标准岩样曲线获取单元还包括:多岩样曲线获取模块103,用于分别获取所述储层的多块标准岩样的核磁共振T2谱曲线和压汞曲线。
图6是本发明实施例曲线对应关系确定单元的结构框图,如图6所示,本发明实施例曲线对应关系确定单元包括:对应关系确定模块201和计算模块202。
对应关系确定模块201,用于根据每块标准岩样的核磁共振T2谱曲线和压汞曲线确定每块标准岩样对应的核磁共振T2谱曲线与压汞曲线的对应关系。
计算模块202,用于将确定的多个对应关系取平均值,得到所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
在本发明的实施例中,根据标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系时,为了提高该对应关系的准确性,也可以设置多块标准岩样,分别检测出每块标准岩样的核磁共振T2谱曲线和压汞曲线,进而根据每块标准岩样的两种曲线来确定储层的核磁共振T2谱曲线与压汞曲线的对应关系。
在本申请实施例中,上述根据多块标准岩样的两种曲线来确定储层的核磁共振T2谱曲线与压汞曲线的对应关系具体可以为,分别根据每块标准岩样的两种曲线求出该块标准岩样对应的所述对应关系,进而根据每块标准岩样求出的对应关系确定出最终较为准确的储层的核磁共振T2谱曲线与压汞曲线的对应关系。在本发明的实施例中所述对应关系可以为T2弛豫时间与喉道半径之间的对应系数,如上述步骤S102中的公式。根据每块标准岩样求出的对应关系确定出最终较为准确的对应关系可以为,根据每块标准岩样求出的T2弛豫时间与喉道半径之间的对应系数,进而通过求平均值得出最终较为准确的储层的核磁共振T2谱曲线与压汞曲线的对应关系。
由以上描述可以看出,本发明提供的岩样微观孔隙结构定量检测装置在对复杂岩样进行微观孔隙结构定量测试时只需核磁共振T2谱曲线,具有实验精度高,实验岩样不需特殊处理,操作简单等优点。此外,本发明的岩样微观孔隙结构定量检测装置成本低且快捷,测试一块复杂岩样核磁共振T2谱曲线所需时间不到5分钟,极大地提高了复杂岩样微观孔隙结构定量检查效率,适合现场、大批量微观孔隙结构测试,为开展储层评价提供基础参数与技术支撑。
为了实现上述目的,根据本申请的另一方面,还提供了另一种岩样微观孔隙结构定量检测装置,包括:ISCO泵、中间容器、核磁共振仪、高压压汞仪以及处理器;
所述中间容器,用于盛装地层水、标准岩样以及与所述标准岩样处于同一储层的待测岩样;
所述ISCO泵,用于给所述中间容器加压,使所述标准岩样和所述待测岩样含水饱和;
所述核磁共振仪,用于分别检测含水饱和的所述标准岩样及检测含水饱和的所述待测岩样的核磁共振T2谱曲线;
高压压汞仪,用于检测含水饱和的所述标准岩样的压汞曲线;
所述处理器,用于根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定储层的核磁共振T2谱曲线与压汞曲线的对应关系;根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
为了实现上述目的,根据本申请的另一方面,还提供了一种计算机设备。如图9所示,该计算机设备包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述岩样微观孔隙结构定量检测方法中的步骤。
处理器可以为中央处理器(Central Processing Unit,CPU)。处理器还可以为其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片,或者上述各类芯片的组合。
存储器作为一种非暂态计算机可读存储介质,可用于存储非暂态软件程序、非暂态计算机可执行程序以及单元,如本发明上述方法实施例中对应的程序单元。处理器通过运行存储在存储器中的非暂态软件程序、指令以及模块,从而执行处理器的各种功能应用以及作品数据处理,即实现上述方法实施例中的方法。
存储器可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储处理器所创建的数据等。此外,存储器可以包括高速随机存取存储器,还可以包括非暂态存储器,例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中,存储器可选包括相对于处理器远程设置的存储器,这些远程存储器可以通过网络连接至处理器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
所述一个或者多个单元存储在所述存储器中,当被所述处理器执行时,执行上述实施例中的方法。
上述计算机设备具体细节可以对应参阅上述实施例中对应的相关描述和效果进行理解,此处不再赘述。
为了实现上述目的,根据本申请的另一方面,还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序在计算机处理器中执行时实现上述岩样微观孔隙结构定量检测方法中的步骤。本领域技术人员可以理解,实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)、随机存储记忆体(RandomAccessMemory,RAM)、快闪存储器(Flash Memory)、硬盘(HardDisk Drive,缩写:HDD)或固态硬盘(Solid-State Drive,SSD)等;所述存储介质还可以包括上述种类的存储器的组合。显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (13)

1.一种岩样微观孔隙结构定量检测方法,其特征在于,包括:
获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线;
根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系;
获取检测到的所述储层的待测岩样的核磁共振T2谱曲线;
根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
2.根据权利要求1所述的岩样微观孔隙结构定量检测方法,其特征在于,所述获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线,包括:
获取检测到的所述储层的第一标准岩样的核磁共振T2谱曲线;以及
获取检测到的所述储层的第二标准岩样的压汞曲线。
3.根据权利要求1所述的岩样微观孔隙结构定量检测方法,其特征在于,所述获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线,包括:
分别获取所述储层的多块标准岩样的核磁共振T2谱曲线和压汞曲线。
4.根据权利要求3所述的岩样微观孔隙结构定量检测方法,其特征在于,所述根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系,包括:
根据每块标准岩样的核磁共振T2谱曲线和压汞曲线确定每块标准岩样对应的核磁共振T2谱曲线与压汞曲线的对应关系;
将确定的多个对应关系取平均值,得到所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
5.根据权利要求1所述的岩样微观孔隙结构定量检测方法,其特征在于,所述储层的核磁共振T2谱曲线与压汞曲线的对应关系包括:T2弛豫时间与喉道半径的对应关系。
6.一种岩样微观孔隙结构定量检测装置,其特征在于,包括:
标准岩样曲线获取单元,用于获取检测到的储层的标准岩样的核磁共振T2谱曲线和压汞曲线;
曲线对应关系确定单元,用于根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定所述储层的核磁共振T2谱曲线与压汞曲线的对应关系;
待测岩样曲线获取单元,用于获取检测到的所述储层的待测岩样的核磁共振T2谱曲线;
微观孔隙结构定量确定单元,用于根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
7.根据权利要求6所述的岩样微观孔隙结构定量检测装置,其特征在于,所述标准岩样曲线获取单元包括:
第一曲线获取模块,用于获取检测到的所述储层的第一标准岩样的核磁共振T2谱曲线;
第二曲线获取模块,用于获取检测到的所述储层的第二标准岩样的压汞曲线。
8.根据权利要求6所述的岩样微观孔隙结构定量检测装置,其特征在于,所述标准岩样曲线获取单元包括:
多岩样曲线获取模块,用于分别获取所述储层的多块标准岩样的核磁共振T2谱曲线和压汞曲线。
9.根据权利要求8所述的岩样微观孔隙结构定量检测装置,其特征在于,所述曲线对应关系确定单元包括:
对应关系确定模块,用于根据每块标准岩样的核磁共振T2谱曲线和压汞曲线确定每块标准岩样对应的核磁共振T2谱曲线与压汞曲线的对应关系;
计算模块,用于将确定的多个对应关系取平均值,得到所述储层的核磁共振T2谱曲线与压汞曲线的对应关系。
10.根据权利要求6所述的岩样微观孔隙结构定量检测装置,其特征在于,所述储层的核磁共振T2谱曲线与压汞曲线的对应关系包括:T2弛豫时间与喉道半径的对应关系。
11.一种岩样微观孔隙结构定量检测装置,其特征在于,包括:ISCO泵、中间容器、核磁共振仪、高压压汞仪以及处理器;
所述中间容器,用于盛装地层水、标准岩样以及与所述标准岩样处于同一储层的待测岩样;
所述ISCO泵,用于给所述中间容器加压,使所述标准岩样和所述待测岩样含水饱和;
所述核磁共振仪,用于分别检测含水饱和的所述标准岩样及检测含水饱和的所述待测岩样的核磁共振T2谱曲线;
高压压汞仪,用于检测含水饱和的所述标准岩样的压汞曲线;
所述处理器,用于根据所述标准岩样的核磁共振T2谱曲线和压汞曲线确定储层的核磁共振T2谱曲线与压汞曲线的对应关系;根据所述待测岩样的核磁共振T2谱曲线以及所述对应关系确定出所述待测岩样的压汞曲线,以根据该压汞曲线定量确定所述待测岩样的微观孔隙结构。
12.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至5任一项方法中的步骤。
13.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序在计算机处理器中执行时实现如权利要求1至5任意一项方法中的步骤。
CN201910206682.4A 2019-03-19 2019-03-19 岩样微观孔隙结构定量检测方法及装置 Pending CN111721683A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910206682.4A CN111721683A (zh) 2019-03-19 2019-03-19 岩样微观孔隙结构定量检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910206682.4A CN111721683A (zh) 2019-03-19 2019-03-19 岩样微观孔隙结构定量检测方法及装置

Publications (1)

Publication Number Publication Date
CN111721683A true CN111721683A (zh) 2020-09-29

Family

ID=72562335

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910206682.4A Pending CN111721683A (zh) 2019-03-19 2019-03-19 岩样微观孔隙结构定量检测方法及装置

Country Status (1)

Country Link
CN (1) CN111721683A (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6008645A (en) * 1997-03-11 1999-12-28 Conoco Inc. Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions
US20120275658A1 (en) * 2011-02-28 2012-11-01 Hurley Neil F Petrographic image analysis for determining capillary pressure in porous media
CN103279647A (zh) * 2013-05-03 2013-09-04 中国石油天然气集团公司 一种基于孔隙特征参数计算致密砂岩储层渗透率的方法
CN104634718A (zh) * 2015-03-05 2015-05-20 中国石油大学(华东) 应用核磁共振表征致密砂岩孔径分布的标定方法
CN108268712A (zh) * 2018-01-05 2018-07-10 中国石油大学(北京) 核磁共振确定孔隙介质毛管压力的方法和装置
CN108956416A (zh) * 2018-06-06 2018-12-07 中国地质大学(北京) 一种基于Matlab分析致密砂岩储层孔隙表征的方法
CN109030311A (zh) * 2018-07-16 2018-12-18 西南石油大学 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法
CN109060623A (zh) * 2018-06-30 2018-12-21 西南石油大学 致密砂岩孔隙大小分布联合表征方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6008645A (en) * 1997-03-11 1999-12-28 Conoco Inc. Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions
US20120275658A1 (en) * 2011-02-28 2012-11-01 Hurley Neil F Petrographic image analysis for determining capillary pressure in porous media
CN103279647A (zh) * 2013-05-03 2013-09-04 中国石油天然气集团公司 一种基于孔隙特征参数计算致密砂岩储层渗透率的方法
CN104634718A (zh) * 2015-03-05 2015-05-20 中国石油大学(华东) 应用核磁共振表征致密砂岩孔径分布的标定方法
CN108268712A (zh) * 2018-01-05 2018-07-10 中国石油大学(北京) 核磁共振确定孔隙介质毛管压力的方法和装置
CN108956416A (zh) * 2018-06-06 2018-12-07 中国地质大学(北京) 一种基于Matlab分析致密砂岩储层孔隙表征的方法
CN109060623A (zh) * 2018-06-30 2018-12-21 西南石油大学 致密砂岩孔隙大小分布联合表征方法
CN109030311A (zh) * 2018-07-16 2018-12-18 西南石油大学 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
唐巨鹏 等: "煤层气赋存运移的核磁共振成像理论及应用", 东北大学出版社 *
李天降等: "核磁共振与压汞法的孔隙结构一致性研究", 《天然气工业》 *
陈欢庆等: "储层孔隙结构研究进展", 《天然气地球科学》 *

Similar Documents

Publication Publication Date Title
EP3593159B1 (en) Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr
CN104697915B (zh) 一种页岩微观孔隙大小及流体分布的分析方法
CN109932297B (zh) 一种致密砂岩储层渗透率的计算方法
EP3088871B1 (en) Raman spectrum detection method
CN109100812B (zh) 基于核磁共振的岩石孔隙分形维数评价方法及装置
CN110927035A (zh) 一种低渗致密砂岩束缚水饱和度计算方法
CN108444881B (zh) 一种适用于陆相泥页岩微纳米尺度储集空间的表征方法
CN109340586B (zh) 一种供水管道泄露的检测方法及系统
US10534055B2 (en) NMR method for determining non-oil volume of a rock sample
RU2013141026A (ru) Способ и устройство изготовления эталонов для неразрушающего контроля пористости
CN105676309A (zh) 裂缝性油藏逆向渗吸采收率预测方法
CN111827968A (zh) 一种基于核磁共振测井的储层非均质性评价方法及装置
US20140203807A1 (en) Methods for performing nmr measurements on porous media
CN109612900A (zh) 一种储层岩石渗透率预测方法、装置及存储介质
CN111721683A (zh) 岩样微观孔隙结构定量检测方法及装置
CN116593373A (zh) 一种卷烟孔隙率及孔隙分布的检测方法及系统
CN111220639B (zh) 基于核磁共振的气驱水时岩心含气饱和度测定方法及装置
CN107843531A (zh) 低渗透低粘油油层孔隙结构的核磁共振参数表征方法
CN106837318B (zh) 岩层稠油含量获取方法和装置
US11566992B2 (en) Method and device for characterizing microbial carbonate pores, and server
WO2013169137A1 (en) Method and apparatus for analyzing a drill core sample
WO2022204297A1 (en) The combined discrete gaussian analysis of micp and nmr t2 distributions of multi-modal carbonate rocks
CN113376069A (zh) 致密储层物性下限的确定方法及装置
CN113253354A (zh) 一种核磁t2谱截止值的确定方法、装置、介质及电子设备
CN114065462B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200929