CN113688515B - 非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法 - Google Patents

非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法 Download PDF

Info

Publication number
CN113688515B
CN113688515B CN202110947752.9A CN202110947752A CN113688515B CN 113688515 B CN113688515 B CN 113688515B CN 202110947752 A CN202110947752 A CN 202110947752A CN 113688515 B CN113688515 B CN 113688515B
Authority
CN
China
Prior art keywords
equivalent
tight sandstone
solid matrix
uniform
fluid
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
CN202110947752.9A
Other languages
English (en)
Other versions
CN113688515A (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.)
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC Research Institute 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 China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN202110947752.9A priority Critical patent/CN113688515B/zh
Publication of CN113688515A publication Critical patent/CN113688515A/zh
Application granted granted Critical
Publication of CN113688515B publication Critical patent/CN113688515B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Public Health (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明公开了一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法，包括：获取非均匀致密砂岩地层固体基质的等效弹性模量数据；计算非均匀致密砂岩地层流体混合物的等效体积模量数据；建立非均匀致密砂岩地层的流体饱和的岩石物理模型；获得非均匀致密砂岩地层孔隙结构参数；计算非均匀致密砂岩地层干岩石骨架弹性模量。本发明公开的一种非均匀致密砂岩地层干岩石骨架弹性模量的计算方法，相比常规经验公式方法，不仅更适用于具有非均匀孔隙结构的致密砂岩地层，而且还提高对非均匀致密砂岩地层岩石骨架模量描述的准确性。

Description

非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法
技术领域
本发明涉及石油天然气地球物理技术领域，具体涉及一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法。
背景技术
致密砂岩受成岩过程中早期岩石作用和矿物胶结作用的影响，具有低孔、低渗及发育有微裂缝的地质特征，孔隙空间以次生孔隙为主，包含少量原生残余粒间孔，还发育有少量微裂缝，复杂的孔隙结构使储层表现出很强的非均质性。
干岩石骨架等效弹性模量对饱和岩石等效模量的预测有重要的影响，目前，干岩石骨架等效弹性模量可以通过岩石物理实验室直接测量获取，但是对于原位岩石性质的研究，尤其是干岩石骨架等效弹性模量的计算，较为困难，常用的方法为经验模型法。经验模型法一般是在大量数据分析的基础上建立干岩石骨架等效弹性模量与矿物弹性模量和孔隙度的公式表达式关系，通过公式表达式对干岩石骨架等效弹性模量进行计算。然而，由于这些经验岩石骨架模型的建立基于一定假设或某一类特定岩石，因此，不同经验模型对地质条件和岩石孔隙度的应用条件是有差异的，具有局限性。
常规油气致密砂岩地层干岩石骨架等效弹性模量与岩石的基质矿物、孔隙度、非均匀孔隙结构等多种因素有关，常用经验模型不能充分描述地层的孔隙特征等因素对干岩石骨架等效弹性模量的影响，导致非均匀致密砂岩地层岩石骨架等效弹性模量计算的准确性低。
发明内容
本发明的目的在于提供一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法，用以解决现有的非均匀致密砂岩地层岩石骨架等效弹性模量计算准确性低的问题。
本发明提供一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法，包括：
获取非均匀致密砂岩地层固体基质的等效弹性模量数据；
计算非均匀致密砂岩地层流体混合物的等效体积模量数据；
根据所述非均匀致密砂岩地层固体基质的等效弹性模量数据和所述非均匀致密砂岩地层流体混合物的等效体积模量数据，结合非均匀致密砂岩地层固体基质、流体混合物和非均匀孔隙结构，建立非均匀致密砂岩地层的流体饱和的岩石物理模型；
根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型，进行非均匀致密砂岩地层孔隙结构参数反演，获得非均匀致密砂岩地层孔隙结构参数；
根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型与所述非均匀致密砂岩地层孔隙结构参数，计算非均匀致密砂岩地层干岩石骨架弹性模量。
进一步地，所述非均匀致密砂岩地层固体基质的等效弹性模量数据包括非均匀致密砂岩地层岩石固体基质的等效体积模量、等效剪切模量、等效密度和等效体积，所述获取非均匀致密砂岩地层固体基质的等效弹性模量数据包括以下步骤：
步骤A1：计算非均匀致密砂岩地层岩石固体基质的等效体积模量；
步骤A2：计算非均匀致密砂岩地层岩石固体基质的等效剪切模量；
步骤A3：分别计算非均匀致密砂岩地层固体基质的等效密度和等效体积。
进一步地，所述步骤A3包括以下步骤：
步骤A31：计算非均匀致密砂岩地层岩石固体基质等效密度，其计算公式为：
Figure BDA0003217326080000021
式中，ρo表示岩石固体基质等效密度；k表示从第1到m种固体基质成分，m表示固体基质成分个数；ρk表示构成固体基质各成分的密度；vk表示构成固体基质各成分的体积含量；
步骤A32：计算非均匀致密砂岩地层岩石固体基质等效体积，其计算公式为：
Figure BDA0003217326080000022
式中，vo表示岩石固体基质等效体积含量；vk表示构成固体基质各成分的体积含量；k表示从第1到m种固体基质成分，m表示固体基质成分个数。
进一步地，所述计算非均匀致密砂岩地层流体混合物的等效体积模量数据的步骤中，假定非均匀致密砂岩地层混合流体的等效剪切模量μf＝0，计算非均匀致密砂岩地层混合流体的等效模量数据，包括非均匀致密砂岩地层混合流体的等效体积模量、等效密度和等效体积，
其中，非均匀致密砂岩地层混合流体的等效体积模量的计算公式为：
Figure BDA0003217326080000023
式中，Kf表示混合流体的等效体积模量；i表示从第1到n种流体成分，n表示流体成分数目；Ki表示构成流体各成分的体积模量；fi表示各构成成分的体积含量；
非均匀致密砂岩地层混合流体的等效密度的计算公式为：
Figure BDA0003217326080000024
式中，ρf表示非均匀致密砂岩地层混合流体的等效密度；i表示从第1到n种流体成分，n表示流体成分数目；fi表示各构成成分的体积含量；ρi表示构成流体各成分的密度；
非均匀致密砂岩地层混合流体的等效体积的计算公式为：
Figure BDA0003217326080000031
式中，f表示流体混合物的等效体积含量；i表示从第1到n种流体成分，n表示流体成分数目；fi表示各构成成分的体积含量。
进一步地，所述建立非均匀致密砂岩地层的流体饱和的岩石物理模型包括以下步骤：
步骤C1：引入软孔比例参数fsoft描述致密砂岩的非均匀孔隙结构；
步骤C2：计算非均匀致密砂岩的孔隙几何因子；
步骤C3：分别计算非均匀致密砂岩的等效体积模量和等效剪切模量；
步骤C4：计算非均匀致密砂岩的密度、纵波速度和横波速度。
本发明还涉及一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算装置，包括：
第一处理单元，用于获取非均匀致密砂岩地层固体基质的等效弹性模量数据；
第二处理单元，用于计算非均匀致密砂岩地层流体混合物的等效体积模量数据；
第三处理单元，用于根据所述非均匀致密砂岩地层固体基质的等效弹性模量数据和所述非均匀致密砂岩地层流体混合物的等效体积模量数据，结合非均匀致密砂岩地层固体基质、流体混合物和非均匀孔隙结构，建立非均匀致密砂岩地层的流体饱和的岩石物理模型；
第四处理单元，用于根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型，进行非均匀致密砂岩地层孔隙结构参数反演，获得非均匀致密砂岩地层孔隙结构参数；
第五处理单元，用于根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型与所述非均匀致密砂岩地层孔隙结构参数，计算非均匀致密砂岩地层干岩石骨架弹性模量。
本发明还涉及一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述方法的步骤。
本发明还涉及一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法的步骤。
与现有技术相比，本发明的有益效果是：
本发明公开了一种非均匀致密砂岩地层干岩石骨架弹性模量的计算方法，通过岩石物理建模、基于模型的孔隙结构参数反演、非均匀致密砂岩地层岩石骨架模量计算等步骤，提供了一种适用于非均匀致密砂岩储层、且能够提高岩石骨架模量计算准确度的方法。根据建立的非均匀致密砂岩岩石物理模型计算岩石骨架模量的方法考虑了矿物基质、孔隙度、孔隙结构等因素，相比常规经验公式方法，不仅更适用于具有非均匀孔隙结构的致密砂岩地层，而且还提高对非均匀致密砂岩地层岩石骨架模量描述的准确性。
具体实施方式
实施例1
实施例1提供一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法，该方法包括以下步骤：
步骤A：获取非均匀致密砂岩地层固体基质的等效弹性模量数据，具体包括以下步骤：
首先，获取非均匀致密砂岩地层固体基质的常规测井数据，包括岩石固体基质中的矿物组分及其含量；获取各成分对应的矿物弹性模量参数，包括体积模量与剪切模量；根据岩石固体基质中的矿物组分及其含量、矿物弹性模量参数，计算非均匀致密砂岩地层固体基质的等效弹性模量等数据，该非均匀致密砂岩地层固体基质的等效弹性模量数据包括非均匀致密砂岩地层岩石固体基质的等效体积模量、等效剪切模量、等效密度和等效体积，具体包括以下步骤：
步骤A1：计算非均匀致密砂岩地层岩石固体基质的等效体积模量，其计算公式如式1至式4所示：
KHS＝(KHS++KHS-)/2 (式1)
KHS+＝Λ(μmax) (式2)
KHS-＝Λ(μmin) (式3)
Figure BDA0003217326080000041
式中，KHS表示岩石固体基质的等效体积模量；
KHS-表示岩石固体基质等效体积模量的下限值，μmin表示岩石固体基质各成分剪切模量的最小值；
KHS+表示岩石固体基质等效体积模量的上限值，μmax表示岩石固体基质各成分剪切模量的最大值；
Λ(l)为表示岩石固体基质等效体积模量上限值或下限值的函数，l表示岩石固体基质各成分剪切模量的最大值μmax或最小值μmin，K(r)表示岩石固体基质第r个组成成分的体积模量，尖括号＜＞表示对岩石固体基质各组成成分按其体积含量求加权平均。
步骤A2：计算非均匀致密砂岩地层岩石固体基质等效剪切模量，其计算公式如式5至式9所示：
μHS＝(μHS++μHS-)/2 (式5)
μHS+＝Γ(ζ(Kmax,μmax)) (式6)
μHS-＝Γ(ζ(Kmin,μmin)) (式7)
Figure BDA0003217326080000051
Figure BDA0003217326080000052
式中，μHS表示岩石固体基质的等效剪切模量；
μHS+表示岩石固体基质等效剪切模量的上限值，Kmax表示岩石固体基质各成分体积模量的最大值，μmax是岩石固体基质各成分剪切模量的最大值；
μHS-表示岩石固体基质等效剪切模量的下限值，Kmin是岩石固体基质各成分体积模量的最小值，μmin是岩石固体基质各成分剪切模量的最小值；
Γ(z)为表示岩石固体基质等效剪切模量上限值μHS+或下限值μHS-的函数，μ(r)表示岩石固体基质第r个组成成分的剪切模量；z为岩石固体基质等效剪切模量上限值μHS+或下限值μHS-的中间变量；
K表示岩石固体基质各成分体积模量的最大值Kmax或最小值Kmin；
μ表示岩石固体基质各成分剪切模量的最大值μmax或最小值μmin。
步骤A3：分别计算非均匀致密砂岩地层固体基质的等效密度和等效体积，其计算公式如式10至式11所示：
步骤A31：计算非均匀致密砂岩地层岩石固体基质等效密度，其计算公式为：
Figure BDA0003217326080000053
式中，ρo表示岩石固体基质等效密度；k表示从第1到m种固体基质成分，m表示固体基质成分个数；ρk表示构成固体基质各成分的密度；vk表示构成固体基质各成分的体积含量。
步骤A32：计算非均匀致密砂岩地层岩石固体基质等效体积，其计算公式为：
Figure BDA0003217326080000054
式中，vo表示岩石固体基质等效体积含量；vk表示构成固体基质各成分的体积含量；k表示从第1到m种固体基质成分，m表示固体基质成分个数。
步骤B：假定非均匀致密砂岩地层混合流体的等效剪切模量μf＝0，计算非均匀致密砂岩地层混合流体的等效模量数据，该等效模量数据包括非均匀致密砂岩地层混合流体的等效体积模量、等效密度和等效体积，具体包括以下步骤：
首先，获取非均匀致密砂岩地层流体混合物的常规测井数据，包括地层流体混合物的流体组分及其含量；获取各成分对应的流体弹性模量参数，包括体积模量与剪切模量；根据非均匀致密砂岩地层流体混合物流体组分及其含量、流体弹性模量参数，计算非均匀致密砂岩地层流体混合物的等效弹性模量等数据，具体包括以下步骤：
其中，非均匀致密砂岩地层混合流体的等效体积模量的计算公式为：
Figure BDA0003217326080000061
式中，Kf表示混合流体的等效体积模量；i表示从第1到n种流体成分，n表示流体成分数目；Ki表示构成流体各成分的体积模量；fi表示各构成成分的体积含量。
非均匀致密砂岩地层混合流体的等效密度的计算公式为：
Figure BDA0003217326080000062
式中，ρf表示非均匀致密砂岩地层混合流体的等效密度；i表示从第1到n种流体成分，n表示流体成分数目；fi表示各构成成分的体积含量；ρi表示构成流体各成分的密度。
非均匀致密砂岩地层混合流体的等效体积含量的计算公式为：
Figure BDA0003217326080000063
式中，f表示流体混合物的等效体积含量；i表示从第1到n种流体成分，n表示流体成分数目；fi表示各构成成分的体积含量。
步骤C：利用步骤A获得的非均匀致密砂岩地层固体基质的等效弹性模量数据和步骤B获得的非均匀致密砂岩地层流体混合物的等效体积模量数据，结合非均匀致密砂岩地层固体基质、流体混合物和非均匀孔隙结构，建立非均匀致密砂岩地层的流体饱和的岩石物理模型，包括以下步骤：
步骤C1：引入软孔比例参数fsoft描述致密砂岩的非均匀孔隙结构；
设岩石的总孔隙空间为总孔隙度φtotal，高纵横比的微观孔隙为硬孔隙，其孔隙度为φstiff，低纵横比的微裂隙为软孔隙，其孔隙度为φsoft，岩石的总孔隙空间等效为硬孔隙与软孔隙这两部分的组合，其表达式为：
φtotal＝φstiff+φsoft (式15)
假设将硬孔隙等效为球形孔隙，将软孔隙等效为微裂缝，即纵横比αstiff＝1的球形孔隙为硬孔隙，纵横比为αsoft＝0.01的微裂缝为软孔隙，岩石物理模型中以总孔隙空间中的软孔比例fsoft＝φsoft/φtotal表征孔隙结构特征；将混合流体填充入硬孔隙、软孔隙中，硬孔隙的体积模量Kstiff、软孔隙的体积模量Ksoft与岩石总体积模量Kf相等，即Kstiff＝Ksoft＝Kf，硬孔隙的剪切模量μstiff、软孔隙的剪切模量μsoft与岩石总剪切模量μf相等，即μstiff＝μsoft＝μf。
步骤C2：计算非均匀致密砂岩的孔隙几何因子；
假设将硬孔隙等效为球形孔隙，将软孔隙等效为微裂缝，即纵横比αstiff＝1的球形孔隙为硬孔隙，纵横比为αsoft＝0.01的微裂缝为软孔隙，分别计算硬孔隙和软孔隙的几何因子，而计算致密砂岩等效体积模量所需的岩石固体基质的几何因子PHS＝1，计算致密砂岩等效剪切模量所需的岩石固体基质的几何因子QHS＝1；
非均匀致密砂岩的孔隙几何因子的计算公式如下：
Figure BDA0003217326080000071
Figure BDA0003217326080000072
Figure BDA0003217326080000073
Figure BDA0003217326080000074
Figure BDA0003217326080000075
Figure BDA0003217326080000076
式中，ξ表示孔隙几何因子的中间变量；
μHS表示岩石固体基质的等效剪切模量；
KHS表示岩石固体基质的等效体积模量；
β表示软孔隙几何因子所需的中间变量；
Pstiff表示致密砂岩体积模量所需的硬孔隙几何因子；
Kstiff表示硬孔隙的体积模量；
Qstiff表示致密砂岩剪切模量所需的硬孔隙几何因子；
Psoft表示致密砂岩体积模量所需的软孔隙几何因子；
Qsoft表示致密砂岩剪切模量所需的软孔隙几何因子；
μsoft表示软孔隙的剪切模量；
αsoft表示软孔隙的纵横比；
Ksoft表示软孔隙的体积模量。
非均匀致密砂岩地层的流体饱和的岩石物理模型包括非均匀致密砂岩的等效体积模量和等效剪切模量、非均匀致密砂岩的密度和纵波速度和横波速度。
步骤C3：分别计算非均匀致密砂岩的等效体积模量和等效剪切模量，包括以下步骤：
步骤C31：计算致密砂岩的等效体积模量，其计算公式为：
Figure BDA0003217326080000081
式中，角标j表示第j相等效致密砂岩组分，N表示等效致密砂岩组分数目；xj为岩石固体基质、硬孔隙和软孔隙各自的体积含量；Kj表示岩石固体基质等效体积模量KHS以及硬孔隙的等效体积模量Kstiff、软孔隙的等效体积模量Ksoft，
Figure BDA0003217326080000082
表示致密砂岩的等效体积模量，Pj表示计算致密砂岩等效体积模量所需的岩石固体基质几何因子PHS、硬孔隙几何因子Pstiff和软孔隙几何因子Psoft。
步骤C32：计算致密砂岩的等效剪切模量，其计算公式为：
Figure BDA0003217326080000083
式中，μj表示岩石固体基质等效剪切模量μHS以及硬孔隙的等效剪切模量μstiff、软孔隙的等效剪切模量μsoft，
Figure BDA0003217326080000084
表示致密砂岩的等效剪切模量，Qj表示计算致密砂岩等效剪切模量所需的岩石固体基质几何因子QHS、硬孔隙几何因子Qstiff和软孔隙几何因子Qsoft。
步骤C4：计算非均匀致密砂岩的密度、纵波速度和横波速度，其计算公式为：
Figure BDA0003217326080000085
Figure BDA0003217326080000086
ρ*＝ρov0+ρff(式26)
式中，Vp表示致密砂岩的纵波速度，
Figure BDA0003217326080000087
表示致密砂岩的等效体积模量，
Figure BDA0003217326080000088
表示致密砂岩的等效剪切模量，ρ*表示致密砂岩的密度；Vs表示致密砂岩的横波速度；ρo表示岩石固体基质等效密度，v0表示岩石固体基质等效体积含量；ρf表示流体混合物的等效密度，f表示流体混合物的等效体积含量。
步骤D：根据所述步骤C建立的非均匀致密砂岩地层的流体饱和的岩石物理模型，进行非均匀致密砂岩地层孔隙结构参数反演，获得非均匀致密砂岩地层孔隙结构参数。
软孔比例参数fsoft为未知参数，为求取该参数并将其应用于非均匀致密砂岩地层干岩石骨架弹性模量的计算，需进行软孔比例参数fsoft的反演。非均匀致密砂岩储层孔隙结构参数的反演依托于模型的正演，基于步骤C中建立的非均匀致密砂岩岩石物理模型，将不同软孔比例fsoft作为拟合参数输入岩石物理模型中，计算得到不同软孔比例fsoft所对应不同的纵波速度和横波速度的正演结果。
通过寻找使得计算的纵波速度和井中实测的纵波速度达到最佳拟合时对应的软孔比例fsoft参数，输出该最佳拟合时对应的软孔比例fsoft参数、纵波速度和横波速度，在缺少井中横波实测数据的情况下可作为对横波速度的估计。最佳拟合时对应的软孔比例fsoft参数，即为反演出的软孔比例的最佳结果。
步骤E：根据所述步骤C建立的非均匀致密砂岩地层的流体饱和的岩石物理模型与所述步骤D获得的非均匀致密砂岩地层孔隙结构参数，计算非均匀致密砂岩地层干岩石骨架弹性模量。
干岩石骨架模量可以通过将流体饱和岩石中孔隙空间填充物即流体的体积模量和剪切模量设为零来模拟。根据步骤D中非均匀致密砂岩地层孔隙结构参数，如软孔比例，在软孔比例所表示的孔隙结构的基础上，将孔隙空间填充物的弹性模量设为零，根据步骤C中建立的非均匀致密砂岩地层的流体饱和的岩石物理模型，计算干岩石骨架弹性模量。
虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims (8)

1.一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法，其特征在于，包括：
获取非均匀致密砂岩地层固体基质的等效弹性模量数据；
计算非均匀致密砂岩地层流体混合物的等效体积模量数据；
根据所述非均匀致密砂岩地层固体基质的等效弹性模量数据和所述非均匀致密砂岩地层流体混合物的等效体积模量数据，结合非均匀致密砂岩地层固体基质、流体混合物和非均匀孔隙结构，建立非均匀致密砂岩地层的流体饱和的岩石物理模型；
根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型，进行非均匀致密砂岩地层孔隙结构参数反演，获得非均匀致密砂岩地层孔隙结构参数；
根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型与所述非均匀致密砂岩地层孔隙结构参数，计算非均匀致密砂岩地层干岩石骨架弹性模量。
2.根据权利要求1所述的方法，其特征在于，所述非均匀致密砂岩地层固体基质的等效弹性模量数据包括非均匀致密砂岩地层岩石固体基质的等效体积模量、等效剪切模量、等效密度和等效体积，所述获取非均匀致密砂岩地层固体基质的等效弹性模量数据包括以下步骤：
步骤A1：计算非均匀致密砂岩地层岩石固体基质的等效体积模量；
步骤A2：计算非均匀致密砂岩地层岩石固体基质的等效剪切模量；
步骤A3：分别计算非均匀致密砂岩地层固体基质的等效密度和等效体积。
3.根据权利要求2所述的方法，其特征在于，所述步骤A3包括以下步骤：
步骤A31：计算非均匀致密砂岩地层岩石固体基质等效密度，其计算公式为：
Figure FDA0003217326070000011
式中，ρo表示岩石固体基质等效密度；k表示从第1到m种固体基质成分，m表示固体基质成分个数；ρk表示构成固体基质各成分的密度；vk表示构成固体基质各成分的体积含量；
步骤A32：计算非均匀致密砂岩地层岩石固体基质等效体积，其计算公式为：
Figure FDA0003217326070000012
式中，vo表示岩石固体基质等效体积含量；vk表示构成固体基质各成分的体积含量；k表示从第1到m种固体基质成分，m表示固体基质成分个数。
4.根据权利要求1所述的方法，其特征在于，在所述计算非均匀致密砂岩地层流体混合物的等效体积模量数据的步骤中，假定非均匀致密砂岩地层混合流体的等效剪切模量μf＝0，计算非均匀致密砂岩地层混合流体的等效模量数据，该等效模量数据包括非均匀致密砂岩地层混合流体的等效体积模量、等效密度和等效体积，
其中，非均匀致密砂岩地层混合流体的等效体积模量的计算公式为：
Figure FDA0003217326070000021
式中，Kf表示混合流体的等效体积模量；i表示从第1到n种流体成分，n表示流体成分数目；Ki表示构成流体各成分的体积模量；fi表示各构成成分的体积含量；
非均匀致密砂岩地层混合流体的等效密度的计算公式为：
Figure FDA0003217326070000022
式中，ρf表示非均匀致密砂岩地层混合流体的等效密度；i表示从第1到n种流体成分，n表示流体成分数目；fi表示各构成成分的体积含量；ρi表示构成流体各成分的密度；
非均匀致密砂岩地层混合流体的等效体积的计算公式为：
Figure FDA0003217326070000023
式中，f表示流体混合物的等效体积含量；i表示从第1到n种流体成分，n表示流体成分数目；fi表示各构成成分的体积含量。
5.根据权利要求1所述的方法，其特征在于，所述建立非均匀致密砂岩地层的流体饱和的岩石物理模型包括以下步骤：
步骤C1：引入软孔比例参数fsoft描述致密砂岩的非均匀孔隙结构；
步骤C2：计算非均匀致密砂岩的孔隙几何因子；
步骤C3：分别计算非均匀致密砂岩的等效体积模量和等效剪切模量；
步骤C4：计算非均匀致密砂岩的密度、纵波速度和横波速度。
6.一种非均匀致密砂岩地层干岩石骨架等效弹性模量的计算装置，其特征在于，包括：
第一处理单元，用于获取非均匀致密砂岩地层固体基质的等效弹性模量数据；
第二处理单元，用于计算非均匀致密砂岩地层流体混合物的等效体积模量数据；
第三处理单元，用于根据所述非均匀致密砂岩地层固体基质的等效弹性模量数据和所述非均匀致密砂岩地层流体混合物的等效体积模量数据，结合非均匀致密砂岩地层固体基质、流体混合物和非均匀孔隙结构，建立非均匀致密砂岩地层的流体饱和的岩石物理模型；
第四处理单元，用于根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型，进行非均匀致密砂岩地层孔隙结构参数反演，获得非均匀致密砂岩地层孔隙结构参数；
第五处理单元，用于根据所述非均匀致密砂岩地层的流体饱和的岩石物理模型与所述非均匀致密砂岩地层孔隙结构参数，计算非均匀致密砂岩地层干岩石骨架弹性模量。
7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-5中任一项所述的方法的步骤。
8.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-5任一所述的方法的步骤。
CN202110947752.9A 2021-08-18 2021-08-18 非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法 Active CN113688515B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110947752.9A CN113688515B (zh) 2021-08-18 2021-08-18 非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110947752.9A CN113688515B (zh) 2021-08-18 2021-08-18 非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法

Publications (2)

Publication Number Publication Date
CN113688515A CN113688515A (zh) 2021-11-23
CN113688515B true CN113688515B (zh) 2022-08-23

Family

ID=78580674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110947752.9A Active CN113688515B (zh) 2021-08-18 2021-08-18 非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法

Country Status (1)

Country Link
CN (1) CN113688515B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110824556A (zh) * 2019-10-22 2020-02-21 中国石油天然气股份有限公司 一种非常规致密砂岩储层的岩石物理模型建立方法及应用
CN111090125A (zh) * 2019-12-30 2020-05-01 中国石油大学(北京) 确定致密砂岩储层的弹性模量的方法和装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110824556A (zh) * 2019-10-22 2020-02-21 中国石油天然气股份有限公司 一种非常规致密砂岩储层的岩石物理模型建立方法及应用
CN111090125A (zh) * 2019-12-30 2020-05-01 中国石油大学(北京) 确定致密砂岩储层的弹性模量的方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
含不连通孔隙的致密砂岩储层岩石弹性模量预测方法;刘倩等;《石油物探》;20151125(第06期);全文 *

Also Published As

Publication number Publication date
CN113688515A (zh) 2021-11-23

Similar Documents

Publication Publication Date Title
US7676349B2 (en) Integrated anisotropic rock physics model
CN112133377B (zh) 一种海域天然气水合物赋存状态的判别方法及系统
CN109655903B (zh) 页岩层横波速度预测方法及系统
CN105487122B (zh) 用于确定储层流体识别因子的敏感性的方法及系统
US10386512B2 (en) System and method for petro-elastic modeling
CN110348590A (zh) 一种页岩吸附气含量计算方法
CN111090125A (zh) 确定致密砂岩储层的弹性模量的方法和装置
CN105425280A (zh) 一种矿物模量和孔隙结构预测方法
CN110344822B (zh) 含钙质致密储层含水饱和度的确定方法
CN113688515B (zh) 非均匀致密砂岩地层干岩石骨架等效弹性模量的计算方法
CN112346130A (zh) 一种富有机质岩石横波速度预测方法、存储介质及系统
CN112147689A (zh) 一种碳酸盐岩储层解释方法及解释量版建立方法
CN112946739B (zh) 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法
CN113031109B (zh) 岩石物理建模方法及装置
CN116930023A (zh) 致密砂岩相控分类孔隙度测井精细解释方法和装置
Jensen et al. Estimation of elastic moduli of mixed porous clay composites
CN113917531B (zh) 基于扰动体积模量的压力系数预测方法及电子设备
CN115079261B (zh) 一种基于多参数定量解释模板的致密砂岩气储层评价方法
CN115248462B (zh) 储层可动性数据的预测方法、系统、存储介质及电子设备
CN112558152B (zh) 一种提高测井横波速度预测精度的方法及系统
CN117238393A (zh) 一种赋存水合物的沉积物的声学特征模拟方法和装置
CN118688854A (zh) 地层破裂压力预测方法及装置
CN114428311A (zh) 一种储层参数直接反演方法及系统
CN118626763A (zh) 基于三维岩石物理模版的物性参数概率化反演方法及装置
Bała The influence of pore media distribution on the elastic parameters of rocks in Miocene sediments (Carpathian Foredeep, Poland, B deposit)

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