CN118091083A - 基于三参数的干燥岩石体积模量计算方法、装置、设备及介质 - Google Patents

基于三参数的干燥岩石体积模量计算方法、装置、设备及介质 Download PDF

Info

Publication number
CN118091083A
CN118091083A CN202410220370.XA CN202410220370A CN118091083A CN 118091083 A CN118091083 A CN 118091083A CN 202410220370 A CN202410220370 A CN 202410220370A CN 118091083 A CN118091083 A CN 118091083A
Authority
CN
China
Prior art keywords
dry rock
bulk modulus
rock sample
physical characteristic
dry
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
CN202410220370.XA
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.)
CNOOC International Energy Services Beijing Ltd
Original Assignee
CNOOC International Energy Services Beijing 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 CNOOC International Energy Services Beijing Ltd filed Critical CNOOC International Energy Services Beijing Ltd
Priority to CN202410220370.XA priority Critical patent/CN118091083A/zh
Publication of CN118091083A publication Critical patent/CN118091083A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

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

Abstract

本发明公开了一种基于三参数的干燥岩石体积模量计算方法、装置、设备及介质。该方法包括：确定目标干燥岩石的第一物理特性参数；根据所述目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量。该方法能够解决现有干燥岩石体积模量计算方法适应性差，从而使得在此基础上形成的地震波速度计算不够准确的问题，提高了干燥岩石体积模量的计算精度。

Description

基于三参数的干燥岩石体积模量计算方法、装置、设备及介质
技术领域
本发明涉及勘探开发技术领域，尤其涉及一种基于三参数的干燥岩石体积模量计算方法、装置、设备及介质。
背景技术
目前随着油气藏勘探开发的不断深化，一方面油气藏勘探逐渐从构造油气藏向隐蔽油气藏等转变，导致油气藏的发现越来越困难，另一方面，油气藏开发中剩余油分布等方面的认识越来越困难，要较好地解决这两个方面面临的困难，需要加大岩石物理理论适应性的研究，特别是适用于一个地区干燥岩石体积模量的研究，干燥岩石体积模量的合理计算往往是岩石物理是否适用一个地区的关键。
但目前因干燥岩石体积模量计算方法适应性差，使得基于现有干燥岩石体积模量计算方法基础上形成的地震波速度计算不够准确，制约了对油气藏勘探开发的正确认识，影响了合适的勘探开发方案部署。
发明内容
本发明提供了一种基于三参数的干燥岩石体积模量计算方法、装置、设备及介质，以解决干燥岩石体积模量计算方法适应性差的问题。
根据本发明的一方面，提供了一种基于三参数的干燥岩石体积模量计算方法，包括：
确定目标干燥岩石的第一物理特性参数，所述第一物理特性参数为对目标干燥岩石进行数据采集得到的数据，所述第一物理特性参数至少包括：目标干燥岩石的密度、目标干燥岩石的孔隙度及有效压力，所述有效压力为单位面积内施加到目标干燥岩石上的压力；
根据所述目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量，所述关系模型为根据干燥岩石样本的第二物理特性参数与干燥岩石样本的干燥岩石体积模量构建的关系式，所述干燥岩石样本为与目标干燥岩石处于同一地区且具有相同物理特性的岩石样本，所述第二物理特性参数为对干燥岩石样本进行数据采集得到的数据。
根据本发明的另一方面，提供了一种基于三参数的干燥岩石体积模量计算装置，包括：
物理参数确定模块，用于确定目标干燥岩石的第一物理特性参数，所述第一物理特性参数为对目标干燥岩石进行数据采集得到的数据，所述第一物理特性参数至少包括：目标干燥岩石的密度、目标干燥岩石的孔隙度及有效压力，所述有效压力为单位面积内施加到目标干燥岩石上的压力；
干燥岩石体积模量确定模块，用于根据所述目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量，所述关系模型为根据干燥岩石样本的第二物理特性参数与干燥岩石样本的干燥岩石体积模量构建的关系式，所述干燥岩石样本为与目标干燥岩石处于同一地区且具有相同物理特性的岩石样本，所述第二物理特性参数为对干燥岩石样本进行数据采集得到的数据。
根据本发明的另一方面，提供了一种电子设备，所述电子设备包括：
至少一个处理器；以及与所述至少一个处理器通信连接的存储器；
其中，所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行本发明任一实施例所述的基于三参数的干燥岩石体积模量计算方法。
根据本发明的另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现本发明任一实施例所述的基于三参数的干燥岩石体积模量计算方法。
本发明实施例的技术方案，根据目标干燥岩石的第一物理特性参数，根据目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量。该方法能够解决现有干燥岩石体积模量计算方法适应性差，从而使得在此基础上形成的地震波速度计算不够准确的问题，提高了干燥岩石体积模量的计算精度。
应当理解，本部分所描述的内容并非旨在标识本发明的实施例的关键或重要特征，也不用于限制本发明的范围。本发明的其它特征将通过以下的说明书而变得容易理解。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种基于三参数的干燥岩石体积模量计算方法的流程图；
图2为本发明实施例提供的一种有效压力为42.76MPa时不同干燥岩石样品干燥岩石体积模量与孔隙度的交会图；
图3为本发明实施例提供的一种有效压力为30.34MPa时不同干燥岩石样品干燥岩石体积模量与孔隙度的交会图；
图4为本发明实施例提供的一种有效压力为42.76MPa时不同干燥岩石样品干燥岩石体积模量与密度的交会图；
图5为本发明实施例提供的一种有效压力为30.34MPa时不同干燥岩石样品干燥岩石体积模量与密度的交会图；
图6为本发明实施例提供的一种干燥岩石样品的干燥岩石体积模量与有效压力的交会图；
图7为本发明实施例提供的一种第一干燥岩石体积模量与实际干燥岩石体积模量的相关性示意图；
图8为本发明实施例提供的一种基于三参数的干燥岩石体积模量计算装置的结构示意图；
图9为实现本发明实施例的基于三参数的干燥岩石体积模量计算的电子设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
图1为本发明实施例提供的一种基于三参数的干燥岩石体积模量计算的流程图，本实施例可适用于获取的干燥岩石体积模量的情况，该方法可以由基于三参数的干燥岩石体积模量计算装置来执行，该基于三参数的干燥岩石体积模量计算装置可以采用硬件和/或软件的形式实现，该基于三参数的干燥岩石体积模量计算装置可配置于任何具有网络通信功能的电子设备中。如图1所示，该方法包括：
S110、确定目标干燥岩石的第一物理特性参数。
其中，第一物理特性参数为对目标干燥岩石进行数据采集得到的数据，第一物理特性参数至少包括：目标干燥岩石的密度、目标干燥岩石的孔隙度及有效压力，有效压力为单位面积内施加到目标干燥岩石上的压力。
示例性的，从S油田中采集目标干燥岩石，对目标干燥岩石经过去油及烘干等处理，测量目标干燥岩石的质量和体积，根据质量和体积得到目标干燥岩石的密度；利用体积法获得干燥岩石的孔隙度；测量获得单位面积内施加到目标干燥岩石上的有效压力。
其中，体积法测定干燥岩石孔隙度的过程为：对样品进行烘干，测量其重量m1；然后让样品饱和水，再次测量其重量m2；接着测量干燥岩石的孔隙度Por＝(m2-m1)/(岩样体积*水的密度)。
S120、根据目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量。
其中，关系模型为根据干燥岩石样本的第二物理特性参数与干燥岩石样本的干燥岩石体积模量构建的关系式，干燥岩石样本为与目标干燥岩石处于同一地区且具有相同物理特性的岩石样本，第二物理特性参数为对干燥岩石样本进行数据采集得到的数据。
将目标干燥岩石的第一物理特性参数输入至关系模型中，经过计算得到第一干燥岩石体积模量。
可选的，关系模型的构建过程，包括步骤A1-A2：
步骤A1、获取干燥岩石样本的第二物理特性参数及第二干燥岩石体积模量数据。
其中，第二干燥岩石体积模量数据为干燥岩石样本在不同有效压力下的干燥岩石体积模量数据。
示例性的，从S油田中采集干燥岩石样本，对干燥岩石样本经过烘干及去油处理，测量干燥岩石样本的质量和体积，根据质量和体积得到目标干燥岩石的密度；利用体积法获得干燥岩石的孔隙度；对干燥岩石样本进行加压测试，得到第二干燥岩石体积模量数据。
可选的，获取干燥岩石样本的第二物理特性参数及第二干燥岩石体积模量数据，包括步骤B1-B4：
步骤B1、对目标岩样进行烘干及去油处理后，进行编号得到干燥岩石样本。
示例性的，将S油田储层获得的N块目标岩样进行烘干及去油处理后，进行编号得到干燥岩石样本i，其中，i＝1,2,…,N。
步骤B2、根据干燥岩石样本对应的质量及体积确定干燥岩石样本的密度。
示例性的，利用n块干燥岩石样本对应的质量和体积获得干燥岩石密度，记为ρi，其中，i＝1,2,…,n。
示例性的，假设共有9块目标岩样，目标岩样进行烘干及去油处理后，利用9块干燥岩石样本对应的质量和体积获得干燥岩石密度，记为ρi，其中，i＝1,2,…,9，如表1所示。
表1 9块干燥岩石样本对应的密度数据
步骤B3、对干燥岩石样本进行孔隙度测量，得到干燥岩石样本的孔隙度。
示例性的，对n块干燥岩石样本利用体积法获得干燥岩石的孔隙度，记为其中，i＝1,2,…,n。
示例性的，对9块干燥岩石样本利用体积法获得干燥岩石的孔隙度，记为其中，i＝1,2,…,9，如表2所示。
表2 9块干燥岩石样本的孔隙度数据
样品名称 孔隙度
S1 0.27
S2 0.22
S3 0.30
S4 0.22
S5 0.20
S6 0.10
S7 0.28
S8 0.30
S9 0.28
步骤B4、对干燥岩石样本进行加压测试，得到第二干燥岩石体积模量数据。
示例性的，对n块干燥岩石样品分别进行m次加压测试，获得n块干燥岩石样品在m次不同有效压力下的干燥岩石体积模量数据，获得n块岩石样品分别在m次不同有效压力下的第二干燥岩石体积模量数据,记为其中，i＝1,2,…,n；j＝1,2,…,m。
示例性的，对9块干燥岩石样本分别进行6次有效压力下的获得干燥岩石体积模量，获得岩石样品分别在6次不同有效压力下的第二干燥岩石体积模量数据,记为其中，i＝1,2,…,n；j＝1,2,…,m。
表3 9块岩样6次不同有效压力下测量的第二干燥岩石体积模量数据
步骤A2、根据干燥岩石样本的第二物理特性参数构建干燥岩石样本的关系模型。
根据干燥岩石样本的第二物理特性参数之间的关系构建干燥岩石样本的关系模型。
可选的，根据干燥岩石样本的第二物理特性参数构建干燥岩石样本的关系模型，包括步骤C1-C4：
步骤C1、根据干燥岩石样本的孔隙度及第二干燥岩石体积模量构建交会图，得到第一近似关系。
其中，第一近似关系为干燥岩石样本的孔隙度及第二干燥岩石体积模量的近似关系。
示例性的，将干燥岩石样本在一定加压测试对应的有效压力下干燥岩石样本的第二干燥岩石体积模量与干燥岩石孔隙度数据生成交会图，对交会图进行模拟分析，分析出干燥岩石体积模量与干燥岩石孔隙度之间存在的近似线性关系，如图2，图3所示。图2提供了干燥岩石第二干燥岩石体积模量在有效压力为42.76MPa时，不同干燥岩石样品第二干燥岩石体积模量与孔隙度的交会图，图中圆圈表示不同样品在不同孔隙度下对应的第二干燥岩石体积模量，绿线是拟合的第二干燥岩石体积模量与孔隙度的关系，其关系具有典型的线性特征，对于其它压力，干燥岩石体积模量与孔隙度也具有这一线性特征。如图2展示出在有效压力为30.34MPa时不同干燥岩石样品第二干燥岩石体积模量与孔隙度的交会图。
步骤C2、根据所述干燥岩石样本的密度及所述第二干燥岩石体积模量构建交会图，得到第二近似关系。
其中，所述第二近似关系为所述干燥岩石样本的密度及所述第二干燥岩石体积模量的近似关系。
示例性的，将一定有效压力下干燥岩石第二干燥岩石体积模量与干燥岩石密度数据生成交会图，对交会图进行分析，分析出干燥岩石第二干燥岩石体积模量与干燥岩石密度之间存在的近似线性关系，如图4，图5所示。图4提供了干燥岩石第二干燥岩石体积模量在有效压力为42.76MPa时，不同干燥岩石样品第二干燥岩石体积模量与密度的交会图，图中圆圈表示不同样品在不同密度对应的第二干燥岩石体积模量，绿线是拟合的第二干燥岩石体积模量与密度的关系，其关系具有典型的线性特征，对于其它压力，第二干燥岩石体积模量与密度也具有这一线性特征。如图5展示出在有效压力为30.34MPa时不同干燥岩石样品第二干燥岩石体积模量与密度的交会图。
步骤C3、根据第二干燥岩石体积模量及有效压力构建交会图，得到第三近似关系。
其中，第三近似关系为第二干燥岩石体积模量及有效压力的近似关系。
示例性的，将一定孔隙度、密度的干燥岩石体积模量与有效压力数据生成交会图，对交会图进行分析，分析出干燥岩石体积模量与有效压力之间近似幂指数关系，如图6所示。
步骤C4、根据第一近似关系、第二近似关系及第三近似关系构建干燥岩石样本的关系模型。
对第一近似关系、第二近似关系及第三近似关系进行分析，在此基础上构建干燥岩石样本的关系模型。
可选的，干燥岩石样本的关系模型，包括：
通过如下公式确定干燥岩石样本的关系模型：
其中，为干燥岩石样本的孔隙度；ρ干燥岩石样本密度，单位为：g/cm3；P为有效压力，单位为：MPa；a0、a1、a2、a3、a4、a5、a6为拟合系数。
可选的，干燥岩石样本的关系模型确定过程，包括：
采用非线性最小二乘法确定干燥岩石样本的关系模型中的拟合参数。
利用Matlab软件中提供的非线性最小二乘算法－列文伯格-马夸尔特法(Levenberg-Marquardt)来确定干燥岩石样本的关系模型中的拟合参数。干燥岩石样本的关系模型中的拟合参数，如表4所示。
表4干燥岩石样本的关系模型中的拟合参数
a0 a1 a2 a3 a4 a5 a6
-3.93 7.11 2.24 -3.55 16.67 0.12 -2.89
其中，Levenberg-Marquardt法为非线性回归中回归参数最小二乘估计的一种估计方法。最优化是寻找使得函数值最小的参数向量。
可选的，根据所述目标干燥岩石的物理特性参数采用关系模型确定第一干燥岩石体积模量后，包括步骤D1-D2：
步骤D1、根据第一干燥岩石体积模量与实际干燥岩石体积模量进行对比，得到第一干燥岩石体积模量与实际干燥岩石体积模量的相关性。
将第一干燥岩石体积模量与实际干燥岩石体积模量进行对比，根据对比的得到的关系曲线确定第一干燥岩石体积模量与实际干燥岩石体积模量的相关性。
步骤D2、根据相关性对关系模型进行优化。
根据第一干燥岩石体积模量与实际干燥岩石体积模量的相关性的大小，对关系模型的拟合参数进行优化，直至相关性满足预设需求。
示例性的，如图7所示为第一干燥岩石体积模量与实际干燥岩石体积模量对比图，从图中可以看出第一干燥岩石体积模量与实际干燥岩石体积模量具有很好的相关性，第一干燥岩石体积模量与实际干燥岩石体积模量的相关性达到了0.9502。
本实施例的技术方案，根据目标干燥岩石的第一物理特性参数，根据目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量。该方法能够解决现有干燥岩石体积模量计算方法适应性差，从而使得在此基础上形成的地震波速度计算不够准确的问题，提高了干燥岩石体积模量的计算精度。
图8为本发明实施例提供的一种基于三参数的干燥岩石体积模量计算装置的结构示意图。本实施例可适用于获取干燥岩石的干燥岩石体积模量的情况，该基于三参数的干燥岩石体积模量计算装置可以采用硬件和/或软件的形式实现，该基于三参数的干燥岩石体积模量计算装置可配置于任何具有网络通信功能的电子设备中。如图8所示，该装置包括：物理参数确定模块210及干燥岩石体积模量确定模块220，其中：
物理参数确定模块210：用于确定目标干燥岩石的第一物理特性参数，所述第一物理特性参数为对目标干燥岩石进行数据采集得到的数据，所述第一物理特性参数至少包括：目标干燥岩石的密度、目标干燥岩石的孔隙度及有效压力，所述有效压力为单位面积内施加到目标干燥岩石上的压力；
干燥岩石体积模量确定模块220：用于根据所述目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量，所述关系模型为根据干燥岩石样本的第二物理特性参数与干燥岩石样本的干燥岩石体积模量构建的关系式，所述干燥岩石样本为与目标干燥岩石处于同一地区且具有相同物理特性的岩石样本，所述第二物理特性参数为对干燥岩石样本进行数据采集得到的数据。
可选的，干燥岩石体积模量确定模块220，包括：
数据获取单元：用于获取干燥岩石样本的第二物理特性参数及第二干燥岩石体积模量数据，所述第二干燥岩石体积模量数据为干燥岩石样本在不同有效下的干燥岩石体积模量数据；
关系模型确定单元：用于根据所述干燥岩石样本的第二物理特性参数构建干燥岩石样本的关系模型。
可选的，数据获取单元，包括：
干燥岩石样本确定子单元：用于对目标岩样进行烘干及去油处理后，进行编号得到干燥岩石样本；
密度确定子单元：用于根据所述干燥岩石样本对应的质量及体积确定所述干燥岩石样本的密度；
孔隙度确定子单元：用于对所述干燥岩石样本进行孔隙度测量，得到干燥岩石样本的孔隙度；
第二干燥岩石体积模量确定子单元：用于对所述干燥岩石样本进行加压测试，得到第二干燥岩石体积模量数据。
可选的，关系模型确定单元，包括：
第一近似关系确定子单元：用于根据所述干燥岩石样本的孔隙度及所述第二干燥岩石体积模量构建交会图，得到第一近似关系，所述第一近似关系为所述干燥岩石样本的孔隙度及所述第二干燥岩石体积模量的近似关系；
第二近似关系确定子单元：用于根根据所述干燥岩石样本的密度及所述第二干燥岩石体积模量构建交会图，得到第二近似关系，所述第二近似关系为所述干燥岩石样本的密度及所述第二干燥岩石体积模量的近似关系；
第三近似关系确定子单元：用于根据所述第二干燥岩石体积模量及所述有效压力构建交会图，得到第三近似关系，所述第三近似关系为所述第二干燥岩石体积模量及有效压力的近似关系；
关系模型确定子单元：根据所述第一近似关系、所述第二近似关系及所述第三近似关系构建干燥岩石样本的关系模型。
可选的，关系模型确定子单元，具体用于：
通过如下公式确定干燥岩石样本的关系模型：
其中，为干燥岩石样本的孔隙度；ρ干燥岩石样本密度，单位为：g/cm3；P为有效压力，单位为：MPa；a0、a1、a2、a3、a4、a5、a6为拟合系数。
可选的，关系模型确定子单元，具体用于：
采用非线性最小二乘法确定所述干燥岩石样本的关系模型中的拟合参数。
可选的，干燥岩石体积模量确定模块220，包括：
相关性确定单元：用于根据所述第一干燥岩石体积模量与实际干燥岩石体积模量进行对比，得到所述第一干燥岩石体积模量与实际干燥岩石体积模量的相关性；
优化单元：用于根据所述相关性对所述关系模型进行优化。
本发明实施例中所提供的基于三参数的干燥岩石体积模量计算装置可执行上述本发明任意实施例中所提供的基于三参数的干燥岩石体积模量计算方法，具备执行该基于三参数的干燥岩石体积模量计算方法相应的功能和有益效果，详细过程参见前述实施例中基于三参数的干燥岩石体积模量计算方法的相关操作。
图9为实现本发明实施例的基于三参数的干燥岩石体积模量计算的电子设备的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备(如头盔、眼镜、手表等)和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。
如图9所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器(ROM)12、随机访问存储器(RAM)13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器(ROM)12中的计算机程序或者从存储单元18加载到随机访问存储器(RAM)13中的计算机程序，来执行各种适当的动作和处理。在RAM 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、ROM 12以及RAM 13通过总线14彼此相连。输入/输出(I/O)接口15也连接至总线14。
电子设备10中的多个部件连接至I/O接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如基于三参数的干燥岩石体积模量计算方法。
在一些实施例中，基于三参数的干燥岩石体积模量计算方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由ROM 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到RAM 13并由处理器11执行时，可以执行上文描述的基于三参数的干燥岩石体积模量计算方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行基于三参数的干燥岩石体积模量计算方法。
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、现场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、负载可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
用于实施本发明的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本发明的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
为了提供与用户的交互，可以在电子设备上实施此处描述的系统和技术，该电子设备具有：用于向用户显示信息的显示装置(例如，CRT(阴极射线管)或者LCD(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(LAN)、广域网(WAN)、区块链网络和互联网。
计算系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与VPS服务中，存在的管理难度大，业务扩展性弱的缺陷。
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。
上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。

Claims (10)

1.一种基于三参数的干燥岩石体积模量计算方法，其特征在于，包括：
确定目标干燥岩石的第一物理特性参数及有效压力，所述第一物理特性参数为对目标干燥岩石进行数据采集得到的数据，所述第一物理特性参数至少包括：目标干燥岩石的密度及目标干燥岩石的孔隙度，所述有效压力为单位面积内施加到目标干燥岩石上的压力；
根据所述目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量，所述关系模型为根据干燥岩石样本的第二物理特性参数与干燥岩石样本的干燥岩石体积模量构建的关系式，所述干燥岩石样本为与目标干燥岩石处于同一地区且具有相同物理特性的岩石样本，所述第二物理特性参数为对干燥岩石样本进行数据采集得到的数据。
2.根据权利要求1所述的方法，其特征在于，所述关系模型的构建过程，包括：
获取干燥岩石样本的第二物理特性参数及第二干燥岩石体积模量数据，所述第二干燥岩石体积模量数据为干燥岩石样本在不同有效下的干燥岩石体积模量数据；
根据所述干燥岩石样本的第二物理特性参数构建干燥岩石样本的关系模型。
3.根据权利要求2所述的方法，其特征在于，获取干燥岩石样本的第二物理特性参数及第二干燥岩石体积模量数据，包括：
对目标岩样进行烘干及去油处理后，进行编号得到干燥岩石样本；
根据所述干燥岩石样本对应的质量及体积确定所述干燥岩石样本的密度；
对所述干燥岩石样本进行孔隙度测量，得到干燥岩石样本的孔隙度；
对所述干燥岩石样本进行加压测试，得到第二干燥岩石体积模量数据。
4.根据权利要求2所述的方法，其特征在于，根据所述干燥岩石样本的第二物理特性参数构建干燥岩石样本的关系模型，包括：
根据所述干燥岩石样本的孔隙度及所述第二干燥岩石体积模量构建交会图，得到第一近似关系，所述第一近似关系为所述干燥岩石样本的孔隙度及所述第二干燥岩石体积模量的近似关系；
根据所述干燥岩石样本的密度及所述第二干燥岩石体积模量构建交会图，得到第二近似关系，所述第二近似关系为所述干燥岩石样本的密度及所述第二干燥岩石体积模量的近似关系；
根据所述第二干燥岩石体积模量及所述有效压力构建交会图，得到第三近似关系，所述第三近似关系为所述第二干燥岩石体积模量及有效压力的近似关系；
根据所述第一近似关系、所述第二近似关系及所述第三近似关系构建干燥岩石样本的关系模型。
5.根据权利要求4所述的方法，其特征在于，所述干燥岩石样本的关系模型，包括：
通过如下公式确定干燥岩石样本的关系模型：
其中，为干燥岩石样本的孔隙度；ρ干燥岩石样本密度，单位为：g/cm3；P为有效压力，单位为：MPa；a0、a1、a2、a3、a4、a5、a6为拟合系数。
6.根据权利要求5所述的方法，其特征在于，所述干燥岩石样本的关系模型确定过程，包括：
采用非线性最小二乘法确定所述干燥岩石样本的关系模型中的拟合参数。
7.根据权利要求1所述的方法，其特征在于，根据所述目标干燥岩石的物理特性参数采用关系模型确定第一干燥岩石体积模量后，包括：
根据所述第一干燥岩石体积模量与实际干燥岩石体积模量进行对比，得到所述第一干燥岩石体积模量与实际干燥岩石体积模量的相关性；
根据所述相关性对所述关系模型进行优化。
8.一种基于三参数的干燥岩石体积模量计算装置，其特征在于，包括：
物理参数确定模块，用于确定目标干燥岩石的第一物理特性参数，所述第一物理特性参数为对目标干燥岩石进行数据采集得到的数据，所述第一物理特性参数至少包括：目标干燥岩石的密度、目标干燥岩石的孔隙度及有效压力，所述有效压力为单位面积内施加到目标干燥岩石上的压力；
干燥岩石体积模量确定模块，用于根据所述目标干燥岩石的第一物理特性参数采用关系模型确定第一干燥岩石体积模量，所述关系模型为根据干燥岩石样本的第二物理特性参数与干燥岩石样本的干燥岩石体积模量构建的关系式，所述干燥岩石样本为与目标干燥岩石处于同一地区且具有相同物理特性的岩石样本，所述第二物理特性参数为对干燥岩石样本进行数据采集得到的数据。
9.一种电子设备，其特征在于，所述电子设备包括：
至少一个处理器；以及与所述至少一个处理器通信连接的存储器；
其中，所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-7中任一项所述的基于三参数的干燥岩石体积模量计算方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现权利要求1-7中任一项所述的基于三参数的干燥岩石体积模量计算方法。
CN202410220370.XA 2024-02-28 2024-02-28 基于三参数的干燥岩石体积模量计算方法、装置、设备及介质 Pending CN118091083A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410220370.XA CN118091083A (zh) 2024-02-28 2024-02-28 基于三参数的干燥岩石体积模量计算方法、装置、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410220370.XA CN118091083A (zh) 2024-02-28 2024-02-28 基于三参数的干燥岩石体积模量计算方法、装置、设备及介质

Publications (1)

Publication Number Publication Date
CN118091083A true CN118091083A (zh) 2024-05-28

Family

ID=91152929

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410220370.XA Pending CN118091083A (zh) 2024-02-28 2024-02-28 基于三参数的干燥岩石体积模量计算方法、装置、设备及介质

Country Status (1)

Country Link
CN (1) CN118091083A (zh)

Similar Documents

Publication Publication Date Title
CN111159904A (zh) 一种用于地热资源量估算评价的方法
WO2019223585A1 (zh) 一种地震震源反演不确定性分析方法、存储介质及服务器
CN110857625B (zh) 碳酸盐岩有封闭水体缝洞单元储量获取方法、装置及设备
CN118091083A (zh) 基于三参数的干燥岩石体积模量计算方法、装置、设备及介质
CN111027892A (zh) 气井的储层渗透率的确定方法、装置和服务器
CN115984501A (zh) 一种地质三维模型建立方法、装置、设备及存储介质
CN115721276A (zh) 一种门静脉压力梯度无创测量方法及装置
CN115907655A (zh) 一种人防建筑模型的审查方法、装置、设备及介质
CN118655169B (zh) 基于高斯函数仿射变换的流体识别方法及装置
CN118351989A (zh) 气藏有效砂体结构模型构建方法、装置、设备及介质
CN116933968A (zh) 一种含油气饱和度的测量方法、装置、电子设备以及介质
CN116738517B (zh) 一种储层预测方法、装置、电子设备及介质
CN117669119A (zh) 各向异性岩石水平主应力确定方法、装置、设备及介质
CN118294337A (zh) 一种孔喉分布确定方法、装置、电子设备及存储介质
CN118443514A (zh) 非密闭取心样品含油饱和度确定方法、装置、设备及介质
CN117972469A (zh) 一种火成岩分类识别方法、装置、设备及存储介质
CN118443905A (zh) 非密闭取心样品含气饱和度确定方法、装置、设备及介质
CN111898785B (zh) 基于页岩的断裂韧性空间分布特征预测方法及系统
CN118313102A (zh) 页岩气藏最终可采储量的预测方法、装置、设备及介质
CN118057390A (zh) 油藏数值模拟拟合精度的确定方法、装置、设备及介质
CN118607684A (zh) 一种抽油机井的产液量计量方法、装置、设备及介质
CN118730843A (zh) 预测岩石孔隙结构方法、装置、电子设备和存储介质
CN116070406A (zh) 页岩气井剩余潜力的评估方法、装置、电子设备及介质
CN116201524A (zh) 断块型圈闭体积校正系数确定方法、装置、设备及介质
CN118191935A (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