CN110263496A - 岩芯特征尺寸确定方法及装置 - Google Patents

岩芯特征尺寸确定方法及装置 Download PDF

Info

Publication number
CN110263496A
CN110263496A CN201910651912.8A CN201910651912A CN110263496A CN 110263496 A CN110263496 A CN 110263496A CN 201910651912 A CN201910651912 A CN 201910651912A CN 110263496 A CN110263496 A CN 110263496A
Authority
CN
China
Prior art keywords
core
size
compression strength
conglomerate
diameter
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
Application number
CN201910651912.8A
Other languages
English (en)
Other versions
CN110263496B (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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201910651912.8A priority Critical patent/CN110263496B/zh
Publication of CN110263496A publication Critical patent/CN110263496A/zh
Application granted granted Critical
Publication of CN110263496B publication Critical patent/CN110263496B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Pathology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本申请提供一种岩芯特征尺寸确定方法及装置,应用于数据处理领域。获取多个不同直径的岩芯数值模型,并通过模拟获得抗压强度,其中,该岩芯数值模型中所有砾石的尺寸为预设圆形固定粒径;根据不同直径的岩芯数值模型的抗压强度,建立岩芯直径与抗压强度关系的拟合函数;通过该拟合函数计算出该圆形固定粒径下岩芯的特征尺寸。获取该特征尺寸下的岩芯样本实体对应的特征抗压强度。如此,在以后的砾岩地层取芯和野外砾岩钻芯取样中,按照该特征尺寸进行钻芯取样,并进行岩石力学实验,得到的岩石抗压强度与砾岩岩样尺寸无关,能够代表实际地层中特定粒径下的抗压强度,有效地减少地质人员的实验工作量,并得到更为可靠的岩石力学参数。

Description

岩芯特征尺寸确定方法及装置
技术领域
本申请涉及数据处理领域,具体而言,涉及一种岩芯特征尺寸确定方法及装置。
背景技术
地质层砾岩的力学参数是工程施工中的重要基础资料,常规获取地质层力学参数的方法通过对储层进行取芯,然后在全直径岩芯上钻取岩样,对岩样进行岩石力学试验,获取储层岩石的强度参数和弹性参数。这种获取地质层砾岩的力学参数的方法对于均质性较好的砂岩和泥岩是可行的,但是对于具有复杂结构的砾岩具有很大的局限性,主要原因是砾岩中的砾石导致砾岩的强度具有较强的尺寸效应,通过小尺寸岩样得到的力学参数不仅与实际地质层的力学参数不符,而且获得的力学参数还具有非常强的离散性。
发明内容
为了克服现有技术中的至少一个不足,本申请的目的之一在于提供一种岩芯特征尺寸确定方法,应用于数据处理设备,所述方法包括:
获取多个模拟抗压强度,所述多个模拟抗压强度由不同直径的岩芯数值模型模拟获得,所述岩芯数值模型中所有砾石的尺寸均为预设特定粒径;
对所述岩芯数值模型的直径以及模拟抗压强度进行拟合,获得反映岩芯直径与抗压强度关系的拟合函数;
将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸;
获取所述特征尺寸的第一岩芯样本实体对应的特征抗压强度,所述特征抗压强度用于代表所述第一岩芯样本实体对应砾岩的抗压强度。
可选地,所述方法还包括:
获取所述砾岩中各砾石的粒径累计分布频率;
将所述分布频率大于预设阈值的砾石对应的粒径作为所述特定粒径。
可选地,所述获取所述砾岩中各砾石的粒径累计分布频率的步骤包括:
获取所述砾岩的第二岩芯样本实体的截面图,所述截面图为垂直于所述第二岩芯样本实体径向的抛光图像;
对所述截面图进行图像处理,获得所述截面图中各个砾石的粒径;
统计所述砾岩中各尺寸的砾石的粒径占比,获得所述各砾石的粒径累计分布频率。
可选地,所述拟合函数的形式如下:
σ=ae-bD+c;
其中,σ为砾岩的单轴抗压强度;D为岩心样本实体的直径;a,b,c为拟合参数。
可选地,所述将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸的步骤包括:
将所述拟合函数的斜率的绝对值与所述预设斜率进行比较,获得所述砾岩的特征尺寸,所述特征尺寸的形式如下:
其中,γ为所述预设斜率。
可选地,所述预设斜率为0.01。
本申请实施例的另一目的在于提供一种岩芯特征尺寸确定装置,应用于数据处理设备,所述岩芯特征尺寸确定装置包括第一获取模块、函数拟合模块、尺寸确定模块和第二获取模块;
所述第一获取模块用于获取多个模拟抗压强度,所述多个模拟抗压强度由不同直径的岩芯数值模型模拟获得,所述岩芯数值模型中所有砾石的尺寸为预设特定粒径;
所述函数拟合模块用于对所述岩芯数值模型的直径以及模拟抗压强度进行拟合,获得反映岩芯直径与抗压强度关系的拟合函数;
所述尺寸确定模块用于将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸;
所述第二获取模块用于获取所述特征尺寸的第一岩芯样本实体对应的特征抗压强度,所述特征抗压强度用于代表所述第一岩芯样本实体对应砾岩的抗压强度。
可选地,所述岩芯特征尺寸确定装置还包括第三获取模块和长度确定模块;
所述第三获取模块用于获取所述砾岩中各砾石的粒径累计分布频率;
所述长度确定模块用于将将所述分布频率大于预设阈值的砾石对应的粒径作为所述特定粒径。
可选地,所述岩芯特征尺寸确定装置还包括第四获取模块、图像处理模块和粒径统计模块;
所述第四获取模块用于获取所述砾岩的第二岩芯样本实体的截面图,所述截面图为所述第二岩芯样本实体的抛光图像;
所述图像处理模块用于对所述截面图进行图像处理,获得所述截面图中各个砾石粒径;
所述粒径统计模块用于统计所述砾岩中各尺寸的砾石的粒径占比,获得所述各砾石的粒径累计分布频率。
可选地,所述拟合函数的形式如下:
σ=ae-bD+c;
其中,σ为砾岩的单轴抗压强度;D为岩心样本实体的直径;a,b,c为拟合参数。
相对于现有技术而言,本申请具有以下有益效果:
本申请实施例提供一种岩芯特征尺寸确定方法及装置,应用于数据处理领域。获取多个不同直径的岩芯数值模型,并通过模拟获得抗压强度,其中,该岩芯数值模型中所有砾石的尺寸为预设特定粒径;根据不同直径的岩芯数值模型的抗压强度,建立岩芯直径与抗压强度关系的拟合函数;通过该拟合函数计算出该特定粒径下岩芯的特征尺寸;获取该特征尺寸下的岩芯样本实体对应的特征抗压强度。如此,在以后的砾岩地层取芯和野外砾岩钻芯取样中,按照该特征尺寸进行钻芯取样,并进行岩石力学实验,得到的岩石抗压强度与砾岩岩样尺寸无关,能够代表实际地层中特定粒径下的抗压强度,有效地减少地质人员的实验工作量,并得到更为可靠的岩石力学参数。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本申请的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本申请实施例提供的数据处理设备的硬件结构图;
图2为本申请实施例提供的岩芯特征尺寸的确定方法的步骤流程图;
图3为本申请实施例提供的粒径分布统计示意图;
图4为本申请实施例提供的砾岩的特征尺寸图;
图5为本申请实施例提供的岩芯特征尺寸确定装置的示意图之一;
图6为本申请实施例提供的岩芯特征尺寸确定装置的示意图之二。
图标:100-数据处理设备;110-岩芯特征尺寸确定装置;120-存储器;130-处理器;1101-第一获取模块;1102-函数拟合模块;1103-尺寸确定模块;1104-第二获取模块;1105-第三获取模块;1106-长度确定模块;1107-第四获取模块;1108-图像处理模块;1109-粒径统计模块。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
在本申请的描述中,需要说明的是,术语“第一”、“第二”、“第三”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
请参照图1,图1是本申请实施例提供的数据处理设备100的硬件结构图,所述数据处理设备100包括岩芯特征尺寸确定装置110、存储器120和处理器130。
所述存储器120和处理器130各元件相互之间直接或间接地电性连接,以实现数据的传输或交互。例如,这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述岩芯特征尺寸确定装置110包括至少一个可以软件或固件(firmware)的形式存储于所述存储器120中或固化在所述数据处理设备的操作系统(operating system,OS)中的软件功能模块。所述处理器130用于执行所述存储器120中存储的可执行模块,例如所述岩芯特征尺寸确定装置110所包括的软件功能模块及计算机程序等。
其中,所述存储器120可以是,但不限于,随机存取存储器(Random AccessMemory,RAM),只读存储器(Read Only Memory,ROM),可编程只读存储器(ProgrammableRead-Only Memory,PROM),可擦除只读存储器(Erasable Programmable Read-OnlyMemory,EPROM),电可擦除只读存储器(Electric Erasable Programmable Read-OnlyMemory,EEPROM)等。其中,存储器120用于存储程序,所述处理器130在接收到执行指令后,执行所述程序。
所述处理器130可能是一种集成电路芯片,具有信号的处理能力。上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
请参照图2,图2为应用于图1所示的数据处理设备100的岩芯特征尺寸确定方法的流程图,通过该方法确定出岩芯的特征尺寸,获取该特征尺寸的第一岩芯样本对应的特征抗压强度,该特征抗压强度与砾岩样品尺寸无关,能够代表实际砾岩地层中特定粒径下的抗压强度。以下将对所述方法包括各个步骤进行详细阐述。
步骤S100,获取多个模拟抗压强度,所述多个模拟抗压强度由不同直径的岩芯数值模型模拟获得,所述岩芯数值模型中所有砾石的尺寸为预设特定粒径。
该数据处理设备100获取砾岩的第二岩芯样本实体对应的物理参数和力学参数,其中,物理参数包括砾石含量以及砾石的粒径分布,力学参数包括砾石与基质的压入硬度比、砾石与基质的压入模量比、单轴抗压强度以及弹性模量。
该数据处理设备100根据第二岩芯样本实体对应的物理参数和力学参数,生成该第二岩芯样本实体对应的数值模型,并通过该数值模型模拟出多个不同直径的虚拟岩芯,以及其对应的模拟抗压强度。
值得说明的是,该数值模型中所有砾石的尺寸为预设特定粒径。由于砾岩中的砾石大小以及砾石含量对砾岩的抗压强度有密切的关系,该数值模型中砾石的粒径以及砾石含量很大程度上决定了该数值模型的准确程度。目前钻取岩芯的设备尺寸有限,不可能获得很大直径的岩芯;这就导致钻取的岩芯具有很大局限性,不同位置钻取的岩芯其砾石含量以及砾径分布具有很大的差别。
本实施例中,为了生成该数值模型,将该第二岩芯样本实体进行相应的处理,获得所述特定粒径。
获取所述砾岩的第二岩芯样本实体的截面图,所述截面图为所述第二岩芯样本实体的抛光图像。
在一种可能的实施方式中,该岩芯样本实体为圆柱结构,该截面图为垂直于第二岩芯样本实体径向的横截面的抛光图像。
在另一种可能的实施方式中,该岩芯样本实体为圆柱结构,该截面图为沿该第二岩芯样本实体径向方向的纵切面的抛光图像。
该数据处理设备100对该界面图进行图像处理,获得该图像中各个砾石的粒径,统计所述砾岩中各尺寸的砾石的粒径占比,获得所述各砾石的粒径累计分布频率。
在一种可能的实施方式中,该数据处理设备100根据砾岩中各尺寸的砾石含量占比,将占比大于第一预设阈值的砾石对应的直径作为该预设直径长度。
请参照图3,在另一种可能的实施方式中,该数据处理设备100根据砾岩中各尺寸的砾石含量占比,按照砾石直径由小到大的方式,累计各直径砾石的含量占比之和,将所述占比和大于第二阈值的砾石直径作为该特定粒径。其中,1~2mm的粒径分布频率为9%,2~4mm的粒径分布频率为32%,4~8mm的粒径分布频率为42%,8~16的粒径分布频率为8%,16~32mm的粒径分布频率为7%,32~64mm的粒径分布频率为2%。该第二阈值为80%,由于1~2mm、2~4mm以及4~8mm的累计粒径分布平率之和为83%,大于80%,则将8mm作为该预设特定粒径。
如此,将数值模型中的砾石的粒径设置为该预设特定粒径,能最大限度的保证最后得到的特征尺寸具有更强的适用性。即如果砾岩中实际砾石的粒径小于预设该预设特定粒径,通过该预设特定粒径下得到的特征尺寸大小的第一岩芯样本进行测试获得的特征抗压强度,依然能够保证获得的该特征抗压强度具有较强的可靠性。
步骤S200,对所述岩芯数值模型的直径以及模拟抗压强度进行拟合,获得反映岩芯直径与抗压强度关系的拟合函数。
该数据处理设备100通过非线性最小二乘法,对虚拟岩芯的直径以及该直径对应的模拟抗压强度进行拟合,获得反映岩芯直径与抗压强度关系的拟合函数。
在一种可能的示例中,该拟合函数的形式如下:
σ=ae-bD+c;
其中,σ为砾岩的单轴抗压强度;D为岩心样本实体的直径;a,b,c为拟合参数。
步骤S300,将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸。
设该预设斜率为γ,获得该拟合函数的导数,将该拟合函数的导数与预设斜率进行比较,获得该砾岩的特征尺寸,其中,该特征尺寸的形式如下:
在一种可能的示例中,该预设斜率为0.01,该预设斜率根据实际的工程需求进行选取。
步骤S400,获取所述特征尺寸的第一岩芯样本实体对应的特征抗压强度,所述特征抗压强度用于代表所述第一岩芯样本实体对应砾岩的抗压强度。
请参照图4,图4为本实施例提供的不同特定粒径下,砾岩的特征尺寸图。
请参照图5,本实施例还提供一种岩芯特征尺寸确定装置110,岩芯特征尺寸确定装置110包括至少一个可以软件形式存储于存储器120中的功能模块。从功能上划分,岩芯特征尺寸确定装置110包括第一获取模块1101、函数拟合模块1102、尺寸确定模块1103和第二获取模块1104。
该第一获取模块1101用于获取多个模拟抗压强度,所述多个模拟抗压强度由不同直径的岩芯数值模型模拟获得,所述岩芯数值模型中所有砾石的尺寸均为预设特定粒径。
在本实施例中,该第一获取模块1101用于执行图2中的步骤S100,关于该第一获取模块1101的详细描述可以参考步骤S100的详细描述。
该函数拟合模块1102用于对所述岩芯数值模型的直径以及模拟抗压强度进行拟合,获得反映岩芯直径与抗压强度关系的拟合函数。
在本实施例中,该函数拟合模块1102用于执行图2中的步骤S200,关于该函数拟合模块1102的详细描述,可以参考步骤S200的详细描述。
该尺寸确定模块1103用于将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸。
在本实施例中,该尺寸确定模块1103用于执行图2中的步骤S300,关于该尺寸确定模块1103的详细描述可以参考步骤S300的详细描述。
该第二获取模块1104用于获取所述特征尺寸的第一岩芯样本实体对应的特征抗压强度,所述特征抗压强度用于代表所述第一岩芯样本实体对应砾岩的抗压强度。
在本实施例中,该第二获取模块1104用于执行图2中的步骤S400,关于该第二获取模块1104的详细描述,可以参考步骤S400的详细描述。
可选地,请参照图6,该所述岩芯特征尺寸确定装置110还包括第三获取模块1105和长度确定模块1106。
该第三获取模块1105用于获取所述砾岩中各砾石的粒径累计分布频率;
该长度确定模块1106用于将将所述分布频率大于预设阈值的砾石对应的粒径作为所述特定粒径。
可选地,请再次参照图6,所述岩芯特征尺寸确定装置110还包括第四获取模块1107、图像处理模块1108和粒径统计模块1109。
该第四获取模块1107用于获取所述砾岩的第二岩芯样本实体的截面图,所述截面图为所述第二岩芯样本实体的抛光图像;
该图像处理模块1108用于对所述截面图进行图像处理,获得所述截面图中各个砾石的粒径;
该粒径统计模块1109用于统计所述砾岩中各尺寸的砾石的粒径占比,获得所述各砾石的粒径累计分布频率。
可选地,所述拟合函数的形式如下:
σ=ae-bD+c;
其中,σ为砾岩的单轴抗压强度;D为岩心样本实体的直径;a,b,c为拟合参数。
综上所述,本申请实施例提供一种岩芯特征尺寸确定方法及装置,应用于数据处理领域。获取多个不同直径的岩芯数值模型,并通过模拟获得抗压强度,其中,该岩芯数值模型中所有砾石的尺寸为预设特定粒径;根据不同直径的岩芯数值模型的抗压强度,建立岩芯直径与抗压强度关系的拟合函数;通过该拟合函数计算出该特定粒径下岩芯的特征尺寸,获取该特征尺寸下的岩芯样本实体对应的特征抗压强度。如此,在以后的砾岩地层取芯和野外砾岩钻芯取样中,按照该特征尺寸进行钻芯取样,并进行岩石力学实验,得到的岩石抗压强度与砾岩岩样尺寸无关,能够代表实际地层中特定粒径下的抗压强度,有效地减少地质人员的实验工作量,并得到更为可靠的岩石力学参数。
在本申请所提供的实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述,仅为本申请的各种实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种岩芯特征尺寸确定方法,其特征在于,应用于数据处理设备,所述方法包括:
获取多个模拟抗压强度,所述多个模拟抗压强度由不同直径的岩芯数值模型模拟获得,所述岩芯数值模型中所有砾石的尺寸均为预设特定粒径;
对所述岩芯数值模型的直径以及模拟抗压强度进行拟合,获得反映岩芯直径与抗压强度关系的拟合函数;
将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸;
获取所述特征尺寸的第一岩芯样本实体对应的特征抗压强度,所述特征抗压强度用于代表所述第一岩芯样本实体对应砾岩的抗压强度。
2.根据权利要求1所述的岩芯特征尺寸确定方法,其特征在于,所述方法还包括:
获取所述砾岩中各砾石的粒径累计分布频率;
将所述分布频率大于预设阈值的砾石对应的粒径作为所述特定粒径。
3.根据权利要求2所述的岩芯特征尺寸确定方法,其特征在于,所述获取所述砾岩中各砾石的粒径累计分布频率的步骤包括:
获取所述砾岩的第二岩芯样本实体的截面图,所述截面图为垂直于所述第二岩芯样本实体径向的抛光图像;
对所述截面图进行图像处理,获得所述截面图中各个砾石的粒径;
统计所述砾岩中各尺寸的砾石的粒径占比,获得所述各砾石的粒径累计分布频率。
4.根据权利要求1所述的岩芯特征尺寸确定方法,其特征在于,所述拟合函数的形式如下:
σ=ae-bD+c;
其中,σ为砾岩的单轴抗压强度;D为岩心样本实体的直径;a,b,c为拟合参数。
5.根据权利要求4所述的岩芯特征尺寸确定方法,其特征在于,所述将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸的步骤包括:
将所述拟合函数的斜率的绝对值与所述预设斜率进行比较,获得所述砾岩的特征尺寸,所述特征尺寸的形式如下:
其中,γ为所述预设斜率。
6.根据权利要求5所述的岩芯特征尺寸确定方法,其特征在于,所述预设斜率为0.01。
7.一种岩芯特征尺寸确定装置,其特征在于,应用于数据处理设备,所述岩芯特征尺寸确定装置包括第一获取模块、函数拟合模块、尺寸确定模块和第二获取模块;
所述第一获取模块用于获取多个模拟抗压强度,所述多个模拟抗压强度由不同直径的岩芯数值模型模拟获得,所述岩芯数值模型中所有砾石的尺寸均为预设特定粒径;
所述函数拟合模块用于对所述岩芯数值模型的直径以及模拟抗压强度进行拟合,获得反映岩芯直径与抗压强度关系的拟合函数;
所述尺寸确定模块用于将所述拟合函数的斜率与预设斜率相等时对应的岩芯直径作为特征尺寸;
所述第二获取模块用于获取所述特征尺寸的第一岩芯样本实体对应的特征抗压强度,所述特征抗压强度用于代表所述第一岩芯样本实体对应砾岩的抗压强度。
8.根据权利要求7所述岩芯特征尺寸确定装置,其特征在于,所述岩芯特征尺寸确定装置还包括第三获取模块和长度确定模块;
所述第三获取模块用于获取所述砾岩中各砾石的粒径累计分布频率;
所述长度确定模块用于将将所述分布频率大于预设阈值的砾石对应的粒径作为所述特定粒径。
9.根据权利要求8所述岩芯特征尺寸确定装置,其特征在于,所述岩芯特征尺寸确定装置还包括第四获取模块、图像处理模块和粒径统计模块;
所述第四获取模块用于获取所述砾岩的第二岩芯样本实体的截面图,所述截面图为所述第二岩芯样本实体的抛光图像;
所述图像处理模块用于对所述截面图进行图像处理,获得所述截面图中各个砾石粒径;
所述粒径统计模块用于统计所述砾岩中各尺寸的砾石的粒径占比,获得所述各砾石的粒径累计分布频率。
10.根据权利要求7所述的岩芯特征尺寸确定装置,其特征在于,所述拟合函数的形式如下:
σ=ae-bD+c;
其中,σ为砾岩的单轴抗压强度;D为岩心样本实体的直径;a,b,c为拟合参数。
CN201910651912.8A 2019-07-18 2019-07-18 岩芯特征尺寸确定方法及装置 Active CN110263496B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910651912.8A CN110263496B (zh) 2019-07-18 2019-07-18 岩芯特征尺寸确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910651912.8A CN110263496B (zh) 2019-07-18 2019-07-18 岩芯特征尺寸确定方法及装置

Publications (2)

Publication Number Publication Date
CN110263496A true CN110263496A (zh) 2019-09-20
CN110263496B CN110263496B (zh) 2021-11-19

Family

ID=67926964

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910651912.8A Active CN110263496B (zh) 2019-07-18 2019-07-18 岩芯特征尺寸确定方法及装置

Country Status (1)

Country Link
CN (1) CN110263496B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113313131A (zh) * 2021-07-29 2021-08-27 四川省冶勘设计集团有限公司 基于图像处理的数字化岩芯辨识方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0274931B1 (fr) * 1986-12-17 1990-10-31 Societe Nationale D'etude Et De Construction De Moteurs D'aviation "Snecma" Dispositif de commande d'aubes à calage variable de redresseur de turbomachine
CN104977210A (zh) * 2015-07-01 2015-10-14 山东科技大学 一种受不同采动影响坚硬岩体Hoek-Brown参数m、s的求取方法
CN105021458A (zh) * 2015-07-14 2015-11-04 中国石油大学(华东) 一种含油泥页岩杨氏模量的定量评价方法
CN107704718A (zh) * 2017-11-27 2018-02-16 中南大学 一种计算岩石材料在压缩试验峰值强度点处弹性应变能密度的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0274931B1 (fr) * 1986-12-17 1990-10-31 Societe Nationale D'etude Et De Construction De Moteurs D'aviation "Snecma" Dispositif de commande d'aubes à calage variable de redresseur de turbomachine
CN104977210A (zh) * 2015-07-01 2015-10-14 山东科技大学 一种受不同采动影响坚硬岩体Hoek-Brown参数m、s的求取方法
CN105021458A (zh) * 2015-07-14 2015-11-04 中国石油大学(华东) 一种含油泥页岩杨氏模量的定量评价方法
CN107704718A (zh) * 2017-11-27 2018-02-16 中南大学 一种计算岩石材料在压缩试验峰值强度点处弹性应变能密度的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIXI LIANG等: "Experimental study on the wettability and adsorption characteristics of Longmaxi Formation shale in the Sichuan Basin, China", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 *
李士斌: "深井岩石破碎规律及破碎的分形机理研究", 《中国博士学位论文全文数据库工程科技I辑》 *
高艳华等: "基于持续屈服节理模型的节理直剪数值试验", 《中南大学学报(自然科学版)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113313131A (zh) * 2021-07-29 2021-08-27 四川省冶勘设计集团有限公司 基于图像处理的数字化岩芯辨识方法及系统
CN113313131B (zh) * 2021-07-29 2021-09-21 四川省冶勘设计集团有限公司 基于图像处理的数字化岩芯辨识方法及系统

Also Published As

Publication number Publication date
CN110263496B (zh) 2021-11-19

Similar Documents

Publication Publication Date Title
CN103548032B (zh) 使用微震事件数据的水力裂缝表征的系统和方法
Zhang et al. Size effect of RQD and generalized representative volume elements: a case study on an underground excavation in Baihetan dam, Southwest China
WO2003067025A2 (en) Interpretation and design of hydraulic fracturing treatments
EP3394647B1 (en) A methodology for building realistic numerical forward stratigraphic models in data sparse environment
CN110532707B (zh) 焖井时间确定方法及设备
Gong et al. AE energy evolution during CJB fracture affected by rock heterogeneity and column irregularity under lateral pressure
CN109635244A (zh) 岩石可钻性预测方法、系统、存储介质及电子终端
CN111985081B (zh) 一种测井曲线构建方法、系统、设备及可读存储介质
CN107861162B (zh) 基于微电极测井资料的天然裂缝识别方法及系统
CN114547953B (zh) 一种基于优化设计图版的压裂施工参数优化方法及系统
CN105283867A (zh) 用于基于平均有效裂缝长度的分布优化现有井和设计新井的系统和方法
CN110263496A (zh) 岩芯特征尺寸确定方法及装置
CN110232208A (zh) 用于识别优势渗流段的储层参数模型构建方法及设备
CN113221347B (zh) 一种井壁稳定性钻井优化方法、装置及设备
Zhu et al. Multicluster fractures propagation during temporary plugging fracturing in naturally fractured reservoirs integrated with dynamic perforation erosion
Wu et al. Integrating complex fracture modeling and EDFM to optimize well spacing in shale oil reservoirs
Wijaya et al. Effect of choke management in optimizing shale-oil production with models of different recovery driving mechanisms
Saeidi et al. Assessment of rock mass erosion in unlined spillways using developed vulnerability and fragility functions
Yang et al. Intelligent rating method of tunnel surrounding rock based on one-dimensional convolutional neural network
Pankaj et al. Back to basics: revisiting proppant and fluid selection for unconventional reservoirs using a new 2D slurry transport model
CN113622901B (zh) 页岩油储层品质分析方法、装置、设备及介质
Moradi et al. Learning from Behavioral Frac Maps: A Montney Case Study in Integration of Modern Microseismic and Production Data Analyses
AU2012396846B2 (en) System, method and computer program product for evaluating and ranking geobodies using a Euler Characteristic
CN109100802A (zh) 储层确定方法和装置
CN113622904A (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