CN113848158B - 二维大岩石模型孔隙度分布测试方法及装置 - Google Patents
二维大岩石模型孔隙度分布测试方法及装置 Download PDFInfo
- Publication number
- CN113848158B CN113848158B CN202010594983.1A CN202010594983A CN113848158B CN 113848158 B CN113848158 B CN 113848158B CN 202010594983 A CN202010594983 A CN 202010594983A CN 113848158 B CN113848158 B CN 113848158B
- Authority
- CN
- China
- Prior art keywords
- dimensional large
- large rock
- rock model
- different positions
- propagation speed
- 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
Links
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/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
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0232—Glass, ceramics, concrete or stone
Abstract
本发明提供了一种二维大岩石模型孔隙度分布测试方法及装置,该方法包括:测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。本发明可以测试二维大岩石模型孔隙度分布,准确度高。
Description
技术领域
本发明涉及石油开发实验技术领域,尤其涉及一种二维大岩石模型孔隙度分布测试方法及装置。
背景技术
石油开发实验技术领域中,二维大岩石模型物理模拟驱油实验可以在室内最大限度地模拟实际油藏条件,利用相似准则转换动态参数,评价不同油藏类型、注入介质、注采方式、井网、井距等条件下的驱油效率和驱替特征,为油田现场注采方案的设计提供重要的参考依据。
孔隙度作为岩石模型重要的基础参数,主要有2个用途:1)确定岩石内部容纳流体的能力,计算原油储量;2)与渗透率共同描述油藏的非均质程度。以往实验中,通常的做法是:在同一块岩石上钻取3-5块柱状岩心,利用其孔隙度的平均值代表整个二维模型的综合孔隙度,虽然柱状岩心和二维模型取自同一块岩石,但物性特征存在较大差异,孔隙度代表性差;另外,二维模型的孔隙度分布近似均质,导致储量计算出现偏差,并无法判断模型的非均质性。
综上所述,目前缺乏一种准确高的二维大岩石模型孔隙度分布测试方法。
发明内容
本发明实施例提出一种二维大岩石模型孔隙度分布测试方法,用以测试二维大岩石模型孔隙度分布,准确度高,该方法包括:
测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;
根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;
根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。
本发明实施例提出一种二维大岩石模型孔隙度分布测试装置,用以测试二维大岩石模型孔隙度分布,准确度高,该装置包括:
首波时间确定模块,用于测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;
岩心孔隙度值确定模块,用于根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;
孔隙度连续分布数据获得模块,用于根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。
本发明实施例还提出了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述二维大岩石模型孔隙度分布测试方法。
本发明实施例还提出了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述二维大岩石模型孔隙度分布测试方法的计算机程序。
在本发明实施例中,测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。在上述过程中,采用的是超声波测量原理,即通过常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间,来确定二维大岩石模型的不同位置处的岩心孔隙度值,相比现有技术中采用孔隙度的平均值代表整个二维模型的综合孔隙度的方法,本发明实施例的方法的准确度更高,可获得全面的孔隙度连续分布数据,且采用超声波测量原理,实现了对二维大岩石模型的无损测试。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中二维大岩石模型孔隙度分布测试方法的流程图;
图2为本发明实施例中超声波首波时间测试原理示意图;
图3为本发明实施例二维大岩石模型孔隙度分布测试方法的详细流程图;
图4为本发明实施例中二维大岩石模型测试点示意图;
图5为本发明实施例绘制的孔隙度连续分布图;
图6为本发明实施例中二维大岩石模型孔隙度分布测试装置的示意图;
图7为本发明实施例中计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
在本说明书的描述中,所使用的“包含”、“包括”、“具有”、“含有”等,均为开放性的用语,即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施,其中的步骤顺序不作限定,可根据需要作适当调整。
图1为本发明实施例中二维大岩石模型孔隙度分布测试方法的流程图,如图1所示,该方法包括:
步骤101,测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;
步骤102,根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;
步骤103,根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。
在本发明实施例提出的方法中,采用的是超声波测量原理,即通过常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间,来确定二维大岩石模型的不同位置处的岩心孔隙度值,相比现有技术中采用孔隙度的平均值代表整个二维模型的综合孔隙度的方法,本发明实施例的方法的准确度更高,可获得全面的孔隙度连续分布数据,且采用超声波测量原理,实现了对二维大岩石模型的无损测试。
具体实施时,图2为本发明实施例中超声波首波时间测试原理示意图,发射探头和接收探头在二维大岩石模型垂直上下表面的位置相同,该位置处的模型厚度为D,常态条件下,二维大岩石模型中饱和空气,测试超声波穿过模型的某个位置的首波时间为t。
在一实施例中,根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值,包括:
根据常态条件下的环境温度,确定超声波在空气中的传播速度;
根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度;
根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值。
在一实施例中,采用如下公式,根据常态条件下的环境温度,确定超声波在空气中的传播速度:
υa=340+0.6×(T-15) (1)
其中,υa为超声波在空气中的传播速度,m/s;T为常态条件下的环境温度,℃。
在上述实施例中,常态条件下的环境温度即测试时的环境温度,由于空气中超声波传播速度对温度较为敏感,所以温度的影响不能忽略,这样计算的超声波在空气中的传播速度更准确。
在一实施例中,采用如下公式,根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度:
其中,υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;tr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的时长,s;Di为二维大岩石模型的位置i处的厚度,cm。
在上述实施例中,抽真空后的二维大岩石模型即岩石骨架模型,超声波在不同类型岩石骨架中的传播速度差异较大,因此需要对二维大岩石模型抽真空后进行标定,上述即标定的过程。
在一实施例中,采用如下公式,根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值:
其中,φi为二维大岩石模型的位置i处的岩心孔隙度值,%;Di为二维大岩石模型的位置i处的厚度,cm;t为二维大岩石模型的位置i处的首波时间,s;υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;υa为超声波在空气中的传播速度,m/s。
在一实施例中,根据二维大岩石模型的不同位置处的孔隙度值,获得二维大岩石模型的孔隙度连续分布数据,包括:
对不同位置处的孔隙度值进行差值处理,获得二维大岩石模型的孔隙度连续分布数据。
在上述实施例中,二维大岩石模型的孔隙度连续分布数据可以采用分布图的形式展示,计算的二维大岩石模型的位置越多,分布图越精确。进行差值处理时,可以采用现有的软件,例如surfur12软件等,提高效率。
基于上述实施例,本发明提出如下一个实施例来说明二维大岩石模型孔隙度分布测试方法的详细流程,图3为本发明实施例二维大岩石模型孔隙度分布测试方法的详细流程图,包括:
步骤301,测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;
步骤302,根据常态条件下的环境温度,确定超声波在空气中的传播速度;
步骤303,根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度;
步骤304,根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值;
步骤305,对不同位置处的孔隙度值进行差值处理,获得二维大岩石模型的孔隙度连续分布数据。
当然,可以理解的是,上述详细流程还可以有其他变化例,相关变化例均应落入本发明的保护范围。
下面给出一个具体实施例,来说明本发明的方法的具体应用。
首先,获得一个二维大岩石模型,模型尺寸(长×宽×厚):30cm×30cm×4.5cm,等间隔6cm测试孔隙度,全模型共计36个测试点,图4为本发明实施例中二维大岩石模型测试点示意图,一个测试点对应一个位置,可适当增加测试点数,以提高孔隙度分布测试精度。
对于每个测试点(即位置),测试常态条件下超声波垂直穿过二维大岩石模型的该位置的首波时间,例如一个某一个位置处为1.5×10-5s。
然后,对二维大岩石模型抽真空不低于2小时,测试得到超声波穿过抽真空后的二维大岩石模型的时长,则利用公式(2)计算得到超声波穿过抽真空后的二维大岩石模型的不同位置处的二维大岩石模型的传播速度,例如某一个位置处为3000m/s。
常态条件下的环境温度为30℃,利用公式(1)计算超声波在空气中的传播速度为349m/s。
按顺序完成图4所示36个测试点的孔隙度值计算。每个测试点首波时间可以测试3次之后取平均值,耗时约1分钟,36个测试点的测试时间基本会控制在1小时内,这也是该方法的另外一个优点,效率高。
本次测试利用surfur12软件绘制孔隙度连续分布图,在绘制前,首先对不同位置处的孔隙度值进行差值处理,图5为本发明实施例绘制的孔隙度连续分布图,可以直观看到模型孔隙度分布的非均质性,并可进一步利用含油饱和度值计算模型原油储量。
综上所述,在本发明实施例提出的方法中,测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。在上述过程中,采用的是超声波测量原理,即通过常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间,来确定二维大岩石模型的不同位置处的岩心孔隙度值,相比现有技术中采用孔隙度的平均值代表整个二维模型的综合孔隙度的方法,本发明实施例的方法的准确度更高,可获得全面的孔隙度连续分布数据,且采用超声波测量原理,实现了对二维大岩石模型的无损测试。
本发明实施例还提出一种二维大岩石模型孔隙度分布测试装置,其原理与二维大岩石模型孔隙度分布测试方法类似,这里不再赘述,图6为本发明实施例中二维大岩石模型孔隙度分布测试装置的示意图,该装置包括:
首波时间确定模块601,用于测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;
岩心孔隙度值确定模块602,用于根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;
孔隙度连续分布数据获得模块603,用于根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。
在一实施例中,岩心孔隙度值确定模块602具体用于:
根据常态条件下的环境温度,确定超声波在空气中的传播速度;
根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度;
根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值。
在一实施例中,岩心孔隙度值确定模块602具体用于:
采用如下公式,根据常态条件下的环境温度,确定超声波在空气中的传播速度:
υa=340+0.6×(T-15)
其中,υa为超声波在空气中的传播速度,m/s;T为常态条件下的环境温度,℃。
在一实施例中,岩心孔隙度值确定模块602具体用于:
采用如下公式,根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度:
其中,υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;tr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的时长,s;Di为二维大岩石模型的位置i处的厚度,cm。
在一实施例中,岩心孔隙度值确定模块602具体用于:
采用如下公式,根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值:
其中,φi为二维大岩石模型的位置i处的岩心孔隙度值,%;Di为二维大岩石模型的位置i处的厚度,cm;t为二维大岩石模型的位置i处的首波时间,s;υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;υa为超声波在空气中的传播速度,m/s。
在一实施例中,孔隙度连续分布数据获得模块603具体用于:
对不同位置处的孔隙度值进行差值处理,获得二维大岩石模型的孔隙度连续分布数据。
综上所述,在本发明实施例提出的装置中,测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据。在上述过程中,采用的是超声波测量原理,即通过常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间,来确定二维大岩石模型的不同位置处的岩心孔隙度值,相比现有技术中采用孔隙度的平均值代表整个二维模型的综合孔隙度的方法,本发明的准确度更高,可获得全面的孔隙度连续分布数据,且采用超声波测量原理,实现了对二维大岩石模型的无损测试。
本申请的实施例还提供一种计算机设备,图7为本发明实施例中计算机设备的示意图,该计算机设备能够实现上述实施例中的二维大岩石模型孔隙度分布测试方法中全部步骤,所述电子设备具体包括如下内容:
处理器(processor)701、存储器(memory)702、通信接口(CommunicationsInterface)703和总线704;
其中,所述处理器701、存储器702、通信接口703通过所述总线704完成相互间的通信;所述通信接口703用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输;
所述处理器701用于调用所述存储器702中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的二维大岩石模型孔隙度分布测试方法中的全部步骤。
本申请的实施例还提供一种计算机可读存储介质,能够实现上述实施例中的二维大岩石模型孔隙度分布测试方法中全部步骤,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的二维大岩石模型孔隙度分布测试方法的全部步骤。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种二维大岩石模型孔隙度分布测试方法,其特征在于,包括:
测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;
根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;
根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据;
根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值,包括:
根据常态条件下的环境温度,确定超声波在空气中的传播速度;
根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度;
根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值;
采用如下公式,根据常态条件下的环境温度,确定超声波在空气中的传播速度:
υa=340+0.6×(T-15)
其中,υa为超声波在空气中的传播速度,m/s;T为常态条件下的环境温度,℃;
采用如下公式,根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度:
其中,υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;tr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的时长,s;Di为二维大岩石模型的位置i处的厚度,m;
采用如下公式,根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值:
其中,φi为二维大岩石模型的位置i处的岩心孔隙度值,%;Di为二维大岩石模型的位置i处的厚度,m;t为二维大岩石模型的位置i处的首波时间,s;υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;υa为超声波在空气中的传播速度,m/s。
2.如权利要求1所述的二维大岩石模型孔隙度分布测试方法,其特征在于,根据二维大岩石模型的不同位置处的孔隙度值,获得二维大岩石模型的孔隙度连续分布数据,包括:
对不同位置处的孔隙度值进行差值处理,获得二维大岩石模型的孔隙度连续分布数据。
3.一种二维大岩石模型孔隙度分布测试装置,其特征在于,包括:
首波时间确定模块,用于测试常态条件下超声波垂直穿过二维大岩石模型的不同位置处的首波时间;
岩心孔隙度值确定模块,用于根据不同位置处的首波时间,确定二维大岩石模型的不同位置处的岩心孔隙度值;
孔隙度连续分布数据获得模块,用于根据二维大岩石模型的不同位置处的岩心孔隙度值,获得二维大岩石模型的孔隙度连续分布数据;
岩心孔隙度值确定模块具体用于:
根据常态条件下的环境温度,确定超声波在空气中的传播速度;
根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度;
根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值;
岩心孔隙度值确定模块具体用于:
采用如下公式,根据常态条件下的环境温度,确定超声波在空气中的传播速度:
υa=340+0.6×(T-15)
其中,υa为超声波在空气中的传播速度,m/s;T为常态条件下的环境温度,℃;
岩心孔隙度值确定模块具体用于:
采用如下公式,根据超声波穿过抽真空后的不同位置处的二维大岩石模型的时长,获得超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度:
其中,υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;tr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的时长,s;Di为二维大岩石模型的位置i处的厚度,m;
岩心孔隙度值确定模块具体用于:
采用如下公式,根据不同位置处的首波时间、超声波在空气中的传播速度和超声波穿过抽真空后的不同位置处的二维大岩石模型的传播速度,确定二维大岩石模型的不同位置处的岩心孔隙度值:
其中,φi为二维大岩石模型的位置i处的岩心孔隙度值,%;Di为二维大岩石模型的位置i处的厚度,m;t为二维大岩石模型的位置i处的首波时间,s;υr,i为超声波穿过抽真空后的二维大岩石模型的位置i处的二维大岩石模型的传播速度,m/s;υa为超声波在空气中的传播速度,m/s。
4.如权利要求3所述的二维大岩石模型孔隙度分布测试装置,其特征在于,孔隙度连续分布数据获得模块具体用于:
对不同位置处的孔隙度值进行差值处理,获得二维大岩石模型的孔隙度连续分布数据。
5.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至2任一项所述方法。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至2任一项所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010594983.1A CN113848158B (zh) | 2020-06-28 | 2020-06-28 | 二维大岩石模型孔隙度分布测试方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010594983.1A CN113848158B (zh) | 2020-06-28 | 2020-06-28 | 二维大岩石模型孔隙度分布测试方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113848158A CN113848158A (zh) | 2021-12-28 |
CN113848158B true CN113848158B (zh) | 2023-09-26 |
Family
ID=78972270
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010594983.1A Active CN113848158B (zh) | 2020-06-28 | 2020-06-28 | 二维大岩石模型孔隙度分布测试方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113848158B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1610433A1 (ru) * | 1988-03-09 | 1990-11-30 | Всесоюзный научно-исследовательский и проектно-конструкторский институт геофизических методов исследований, испытания и контроля нефтегазоразведочных скважин | Способ определени скорости распространени поперечных ультразвуковых волн через керн |
GB9520105D0 (en) * | 1994-09-30 | 1995-12-06 | Core Holdings Bv | Contact acoustic porosity sensor and method for determining continuous rock porosity using same |
CN102262041A (zh) * | 2011-04-20 | 2011-11-30 | 中国石油天然气股份有限公司 | 一种基于多谱孔隙结构分析的饱和度确定方法 |
CN106483056A (zh) * | 2016-09-28 | 2017-03-08 | 西南石油大学 | 一种基于纵波波速的页岩孔隙度测量方法及测量装置 |
CN107449826A (zh) * | 2017-08-17 | 2017-12-08 | 中国石油大学(华东) | 流体分布检测方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150023564A1 (en) * | 2013-07-18 | 2015-01-22 | Chevron U.S.A. Inc. | System and method for estimating porosity distribution in subterranean reservoirs |
EP3121625A1 (en) * | 2015-07-20 | 2017-01-25 | CGG Services SA | Predicting mechanical and elastic rock properties of the subsurface |
-
2020
- 2020-06-28 CN CN202010594983.1A patent/CN113848158B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1610433A1 (ru) * | 1988-03-09 | 1990-11-30 | Всесоюзный научно-исследовательский и проектно-конструкторский институт геофизических методов исследований, испытания и контроля нефтегазоразведочных скважин | Способ определени скорости распространени поперечных ультразвуковых волн через керн |
GB9520105D0 (en) * | 1994-09-30 | 1995-12-06 | Core Holdings Bv | Contact acoustic porosity sensor and method for determining continuous rock porosity using same |
CN102262041A (zh) * | 2011-04-20 | 2011-11-30 | 中国石油天然气股份有限公司 | 一种基于多谱孔隙结构分析的饱和度确定方法 |
CN106483056A (zh) * | 2016-09-28 | 2017-03-08 | 西南石油大学 | 一种基于纵波波速的页岩孔隙度测量方法及测量装置 |
CN107449826A (zh) * | 2017-08-17 | 2017-12-08 | 中国石油大学(华东) | 流体分布检测方法及装置 |
Non-Patent Citations (2)
Title |
---|
一种考虑成岩作用的砂岩孔隙度解释方法;耿斌 等;石油学报(第05期);P44-47 * |
油藏条件下CO2驱岩心内流体动态分布检测技术研究;车新跃;中国优秀硕士学位论文全文数据库 工程科技I辑(第7期);第二章和第三章和第四章 * |
Also Published As
Publication number | Publication date |
---|---|
CN113848158A (zh) | 2021-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108387933B (zh) | 一种确定地层品质因子的方法、装置及系统 | |
CN105891089B (zh) | 储层渗透率确定方法及装置 | |
RU2012152063A (ru) | Способ q томографии | |
CN112946783B (zh) | 一种水合物饱和度确定方法、装置及设备 | |
CN103884633A (zh) | 一种确定岩石渗透率的方法及装置 | |
CN106383365B (zh) | 一种利用图版校正火成岩核磁共振孔隙度的方法 | |
WO2017023789A1 (en) | Method for estimating permeability of fraclured rock formations from induced slow fluid pressure waves | |
CN105628486B (zh) | 测量页岩岩石力学性质的方法及装置 | |
CN106767671B (zh) | 基于三维电子罗盘的地质结构面产状计算方法 | |
CN107356523B (zh) | 基于离散切线刚度估计的实时混合模拟试验反馈力修正方法 | |
CN106556877B (zh) | 一种地磁通化方法及装置 | |
CN103414451A (zh) | 一种应用于飞行器姿态估计的扩展卡尔曼滤波方法 | |
CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
CN105223616A (zh) | 一种页岩储层的孔隙纵横比反演方法 | |
CN113848158B (zh) | 二维大岩石模型孔隙度分布测试方法及装置 | |
CN102087107B (zh) | 绳系多传感器协同优化近海测波浮标及其滤波融合方法 | |
CN111221038B (zh) | 薄储层厚度定量预测的方法和装置 | |
CN109800508A (zh) | 嵌岩桩桩端的空洞顶板厚度的计算方法及终端设备 | |
CN112419493B (zh) | 页岩储层三维属性模型建立方法及装置 | |
CN109031437A (zh) | 随钻核磁共振自旋回波校正方法与装置 | |
CN112925018B (zh) | 一种异常压力成因判别方法及装置 | |
CN104635246A (zh) | 卫星导航信号的动态范围检测系统及检测方法 | |
CN115730175A (zh) | 一种海洋和航空重磁数据的修正方法和装置 | |
CN110287545A (zh) | 一种计算油藏内部任意部位原始含水饱和度的方法及系统 | |
CN116973972B (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 |