CN106960264B - 低渗透砂砾岩地层压力预测方法及装置 - Google Patents

低渗透砂砾岩地层压力预测方法及装置 Download PDF

Info

Publication number
CN106960264B
CN106960264B CN201710217188.9A CN201710217188A CN106960264B CN 106960264 B CN106960264 B CN 106960264B CN 201710217188 A CN201710217188 A CN 201710217188A CN 106960264 B CN106960264 B CN 106960264B
Authority
CN
China
Prior art keywords
longitudinal wave
pressure
wave velocity
stress
core
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
CN201710217188.9A
Other languages
English (en)
Other versions
CN106960264A (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.)
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 CN201710217188.9A priority Critical patent/CN106960264B/zh
Publication of CN106960264A publication Critical patent/CN106960264A/zh
Application granted granted Critical
Publication of CN106960264B publication Critical patent/CN106960264B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Agronomy & Crop Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Husbandry (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本申请实施例提供了一种低渗透砂砾岩地层压力预测方法及装置，该方法包括：在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度；获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率；并获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率；根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数；根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力。本申请实施例可提高低渗透砂砾岩地层压力预测的准确度。

Description

低渗透砂砾岩地层压力预测方法及装置
技术领域
本申请涉及油气勘探开发技术领域，尤其是涉及一种低渗透砂砾岩地层压力预测方法及装置。
背景技术
有效应力是岩石力学的一个基本概念，最早是由Terzaghi提出的，他把有效应力定义为上层总压力与孔隙流体压力的差值：Pd＝Pc-Pp。其中，Pd表示Terzaghi有效应力，即压差；Pc是围压；Pp是孔隙流体压力。有效应力把围压和孔隙流体压力这两个独立的变量简化为一个变量，从而方便了对多孔材料性质的分析，围压和孔隙流体压力变化引起的任意物理性质的变化都可以用有效应力来描述。因此，确定有效应力对于油气藏岩石的渗透率和孔隙度的影响意义重大。
研究表明，有效应力与材料性质具有一一对应关系，只要有效应力保持不变，无论Pc和Pp如何变化，性质Q都不会改变。假定函数Q(Pc，Pp)足够光滑，根据二元函数的全微分形式：
Figure BDA0001262584860000011
整理得到：
Figure BDA0001262584860000012
上式可写为：
Figure BDA0001262584860000013
其中dPeff＝dPc-nQdPp，式中的系数nQ就称为性质Q的有效应力系数。有效应力系数反映了性质Q对Pc和Pp变化的相对敏感程度。Robin指出，不同性质的材料具有不同的有效应力系数。
在油气勘探开发工程中，对于低渗透砂砾岩储层，由于其结构的特殊性，目前尚无可行的低渗透砂砾岩储层的地层压力的预测方案。
发明内容
本申请实施例的目的在于提供一种低渗透砂砾岩地层压力预测方法及装置，以实现低渗透砂砾岩地层压力预测。
为达到上述目的，一方面，本申请实施例提供了一种低渗透砂砾岩地层压力预测方法，包括以下步骤：
在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度；
获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率；并获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率；
根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数；
根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力。
本申请实施例的低渗透砂砾岩地层压力预测方法，所述根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数，包括：
根据公式
Figure BDA0001262584860000021
计算所述岩心的有效应力系数；
其中，n为所述岩心的有效应力系数，
Figure BDA0001262584860000022
为第一变化率，
Figure BDA0001262584860000023
为第二变化率，Vp1为第一纵波速度，Pf为孔隙压力，Vp2为第二纵波速度，Pd为差应力。
本申请实施例的低渗透砂砾岩地层压力预测方法，所述根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力，包括：
根据公式Pf＝Pov-nσ计算所述岩心所对应储层的地层压力；
其中，pf为地层压力，Pov为上覆岩层压力，n为所述岩心的有效应力系数，σ为骨架有效应力。
本申请实施例的低渗透砂砾岩地层压力预测方法，所述获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率，包括：
确定在所述设定孔隙压力下，各个指定差应力时的第一纵波速度与有效应力的第一关系曲线；
根据所述第一关系曲线确定所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率。
本申请实施例的低渗透砂砾岩地层压力预测方法，所述获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率，包括：
确定在所述设定差应力下，各个指定孔隙压力时的第二纵波速度与有效应力的第二关系曲线；
根据所述第二关系曲线确定所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率。
另一方面，本申请实施例还提供了一种低渗透砂砾岩地层压力预测装置，包括：
纵波速度获取模块，用于在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度；
变化率获取模块，用于获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率；并获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率；
有效应力系数获取模块，用于根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数；
地层压力获取模块，用于根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力。
本申请实施例的低渗透砂砾岩地层压力预测装置，所述根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数，包括：
根据公式
Figure BDA0001262584860000031
计算所述岩心的有效应力系数；
其中，n为所述岩心的有效应力系数，
Figure BDA0001262584860000032
为第一变化率，
Figure BDA0001262584860000033
为第二变化率，Vp1为第一纵波速度，Pf为孔隙压力，Vp2为第二纵波速度，Pd为差应力。
本申请实施例的低渗透砂砾岩地层压力预测装置，所述根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力，包括：
根据公式Pf＝Pov-nσ计算所述岩心所对应储层的地层压力；
其中，pf为地层压力，Pov为上覆岩层压力，n为所述岩心的有效应力系数，σ为骨架有效应力。
本申请实施例的低渗透砂砾岩地层压力预测装置，所述获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率，包括：
确定在所述设定孔隙压力下，各个指定差应力时的第一纵波速度与有效应力的第一关系曲线；
根据所述第一关系曲线确定所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率。
本申请实施例的低渗透砂砾岩地层压力预测装置，所述获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率，包括：
确定在所述设定差应力下，各个指定孔隙压力时的第二纵波速度与有效应力的第二关系曲线；
根据所述第二关系曲线确定所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率。
本申请实施例在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取岩心在各个指定孔隙压力时的第二纵波速度；然后根据第一纵波速度在设定孔隙压力下随差应力的第一变化率，以及第二纵波速度在设定差应力下随孔隙压力的第二变化率确定岩心的有效应力系数；最后根据岩心的有效应力系数确定岩心所对应储层的地层压力，从而提供了一种低渗透砂砾岩储层的地层压力的预测方案。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本申请一实施例的低渗透砂砾岩地层压力预测方法流程图；
图2为本申请一实施例中低渗透砂砾岩实验测得的纵波速度与有效压力的关系；
图3为本申请一实施例中两种有效应力系数与孔隙度之间的关系；
图4a为利用现有技术的Biot系数得到的压力系数与埋深的关系；
图4b为利用本申请一实施例的有效应力系数得到的压力系数与埋深的关系；
图5a为利用现有技术的Biot系数得到的压力系数与实测压力系数的关系；
图5b为利用本申请一实施例的有效应力系数得到压力系数与实测压力系数的关系；
图6为本申请一实施例的低渗透砂砾岩地层压力预测装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
参考图1所示，本申请实施例的低渗透砂砾岩地层压力预测方法包括以下步骤：
S101、在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度。
本申请实施例中，所述岩心为目标地层的岩心。所述在设定孔隙压力下即保持孔隙压力在某一固定值下。同样，在设定差应力下即保持差应力在某一固定值下。其中，所述差应力为岩心围压与孔隙压力的差值。
本申请实施例中，所述在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度可以为通过岩石物理实验测量差应力分别为5MPa、10MPa、20MPa、30MPa、40MPa和50MPa等时的纵波速度。同样，所述在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度可以为通过岩石物理实验测量孔隙压力分别为5MPa、10MPa、20MPa、30MPa、40MPa和50MPa等时的纵波速度。
S102、获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率；并获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率。
本申请实施例中，所述获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率可以包括以下步骤：
确定在所述设定孔隙压力下，各个指定差应力时的第一纵波速度与有效应力的第一关系曲线；在上述获取岩心在各个指定差应力时的第一纵波速度后，可以记录测量值，并将其绘制在有效应力与纵波速度的交会图中，从而可以得到第一关系曲线，如图2中左侧虚线所示；
根据所述第一关系曲线确定所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率。当第一关系曲线确定后，根据所述第一关系曲线可以确定所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率。
本申请实施例中，所述获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率，可以包括：
确定在所述设定差应力下，各个指定孔隙压力时的第二纵波速度与有效应力的第二关系曲线；在上述获取所述岩心在各个指定孔隙压力时的第二纵波速度后，可以记录测量值，并将其绘制在有效应力与纵波速度的交会图中，从而可以得到第二关系曲线，如图2中右侧虚线所示；
根据所述第二关系曲线确定所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率。当第二关系曲线确定后，根据所述第二关系曲线可以确定所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率。
S103、根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数。
所述根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数，包括：
本申请实施例中，根据公式
Figure BDA0001262584860000061
计算所述岩心的有效应力系数；
其中，n为所述岩心的有效应力系数，
Figure BDA0001262584860000062
为第一变化率，
Figure BDA0001262584860000063
为第二变化率，Vp1为第一纵波速度，Pf为孔隙压力，Vp2为第二纵波速度，Pd为差应力。
从公式
Figure BDA0001262584860000064
可以看出，本申请实施例中，有效应力系数在计算时分子为差压力一定时速度随孔隙压力的变化率，分母为孔隙压力一定时速度随差应力的变化率。这两种变化率通过岩石物理实验可求得。本申请的发明人进一步研究发现：低渗透砂砾岩有效应力系数与常规的Biot有效应力系数相比有以下几点差异：
①同等条件下本申请实施例的有效应力系数高于Biot有效应力系数；
②两种有效应力系数与孔隙度之间有明显的相关性，主要表现为随孔隙度增大有效应力系数均增大且接近于1；但在低孔低渗地区两者的差异较大，特别是孔隙度在8％～10％之间时，Biot有效应力系数为0.3～0.4之间而本申请实施例的低渗透砂砾岩有效应力系数在0.7～0.8之间。
③如图3所示，两类有效应力系数的差异随孔隙度的增大而变小，即在常规储层中，两类有效应力系数是可以通用的。
S104、根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力。
本申请实施例中，所述根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力，包括：
根据公式Pf＝Pov-nσ计算所述岩心所对应储层的地层压力；
其中，pf为地层压力，Pov为上覆岩层压力，n为所述岩心的有效应力系数，σ为骨架有效应力。
本申请实施例中，根据有效应力定理，利用低渗透砂砾岩储层的有效应力系数可以确定地层压力，然后通过地层压力与静水压力的比值确定异常压力地层段。
为了验证本申请实施例的效果，本申请发明人在准噶尔盆地西北缘玛北地区，分别利用现有技术的Biot系数得到的压力系数与埋深的关系(如图4a所示)、压力系数与实测压力系数(如图5a所示)；利用本申请请实施例的有效应力系数得到的压力系数与埋深的关系(如图4b所示)、压力系数与实测压力系数(如图4b所示)，对比图4a和图4b，以及图5a和图5b可知，利用本申请一实施例的有效应力系数得到压力系数与地质规律吻合的更好，主要表现为：
①用单井计算的地层压力系数随埋深的增大而增大的规律。
②新的地层压力预测模型预测的压力系数与实测的线性关系更好，相关度更高。
③利用新模型计算的压力系数与甜点储层预测的结果吻合较好，且能够体现出在甜点储层发育区随着埋深的增大地层压力系数增大的特征。
虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
参考图6所示，为本申请实施例的低渗透砂砾岩地层压力预测装置，其可以包括：
纵波速度获取模块61，可以用于在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度；
变化率获取模块62，可以用于获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率；并获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率；
有效应力系数获取模块63，可以用于根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数；
地层压力获取模块64，可以用于根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力。
本申请实施例的低渗透砂砾岩地层压力预测装置与上述图1所示的低渗透砂砾岩地层压力预测方法对应，因此，有关于本申请实施例的低渗透砂砾岩地层压力预测装置的细节，请参见上述图1所示的低渗透砂砾岩地层压力预测方法，在此不再赘述。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (8)

1.一种低渗透砂砾岩地层压力预测方法，其特征在于，包括以下步骤：
在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度；
获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率；并获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率；
根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数；
根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力；
其中，所述根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数，包括：
根据公式
Figure FDA0002380848340000011
计算所述岩心的有效应力系数；
其中，n为所述岩心的有效应力系数，
Figure FDA0002380848340000012
为第一变化率，
Figure FDA0002380848340000013
为第二变化率，Vp1为第一纵波速度，Pf为孔隙压力，Vp2为第二纵波速度，Pd为差应力。
2.根据权利要求1所述的低渗透砂砾岩地层压力预测方法，其特征在于，所述根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力，包括：
根据公式Pf＝Pov-nσ计算所述岩心所对应储层的地层压力；
其中，pf为地层压力，Pov为上覆岩层压力，n为所述岩心的有效应力系数，σ为骨架有效应力。
3.根据权利要求1所述的低渗透砂砾岩地层压力预测方法，其特征在于，所述获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率，包括：
确定在所述设定孔隙压力下，各个指定差应力时的第一纵波速度与有效应力的第一关系曲线；
根据所述第一关系曲线确定所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率。
4.根据权利要求1所述的低渗透砂砾岩地层压力预测方法，其特征在于，所述获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率，包括：
确定在所述设定差应力下，各个指定孔隙压力时的第二纵波速度与有效应力的第二关系曲线；
根据所述第二关系曲线确定所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率。
5.一种低渗透砂砾岩地层压力预测装置，其特征在于，包括：
纵波速度获取模块，用于在设定孔隙压力下，获取岩心在各个指定差应力时的第一纵波速度；并在设定差应力下，获取所述岩心在各个指定孔隙压力时的第二纵波速度；
变化率获取模块，用于获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率；并获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率；
有效应力系数获取模块，用于根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数；
地层压力获取模块，用于根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力；
其中，所述根据所述第一变化率和所述第二变化率确定所述岩心的有效应力系数，包括：
根据公式
Figure FDA0002380848340000021
计算所述岩心的有效应力系数；
其中，n为所述岩心的有效应力系数，
Figure FDA0002380848340000022
为第一变化率，
Figure FDA0002380848340000023
为第二变化率，Vp1为第一纵波速度，Pf为孔隙压力，Vp2为第二纵波速度，Pd为差应力。
6.根据权利要求5所述的低渗透砂砾岩地层压力预测装置，其特征在于，所述根据所述岩心的有效应力系数确定所述岩心所对应储层的地层压力，包括：
根据公式Pf＝Pov-nσ计算所述岩心所对应储层的地层压力；
其中，pf为地层压力，Pov为上覆岩层压力，n为所述岩心的有效应力系数，σ为骨架有效应力。
7.根据权利要求5所述的低渗透砂砾岩地层压力预测装置，其特征在于，所述获取所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率，包括：
确定在所述设定孔隙压力下，各个指定差应力时的第一纵波速度与有效应力的第一关系曲线；
根据所述第一关系曲线确定所述第一纵波速度在所述设定孔隙压力下随差应力的第一变化率。
8.根据权利要求5所述的低渗透砂砾岩地层压力预测装置，其特征在于，所述获取所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率，包括：
确定在所述设定差应力下，各个指定孔隙压力时的第二纵波速度与有效应力的第二关系曲线；
根据所述第二关系曲线确定所述第二纵波速度在所述设定差应力下随孔隙压力的第二变化率。
CN201710217188.9A 2017-04-05 2017-04-05 低渗透砂砾岩地层压力预测方法及装置 Active CN106960264B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710217188.9A CN106960264B (zh) 2017-04-05 2017-04-05 低渗透砂砾岩地层压力预测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710217188.9A CN106960264B (zh) 2017-04-05 2017-04-05 低渗透砂砾岩地层压力预测方法及装置

Publications (2)

Publication Number Publication Date
CN106960264A CN106960264A (zh) 2017-07-18
CN106960264B true CN106960264B (zh) 2020-11-06

Family

ID=59483274

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710217188.9A Active CN106960264B (zh) 2017-04-05 2017-04-05 低渗透砂砾岩地层压力预测方法及装置

Country Status (1)

Country Link
CN (1) CN106960264B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020131039A1 (en) * 2018-12-18 2020-06-25 Halliburton Energy Services, Inc. Determining when applied stress to a core rock sample has equilibrated in the core rock sample
CN111101930B (zh) * 2019-12-19 2022-05-20 西南石油大学 一种气藏开发模式下单井开采提产潜力评价方法
CN114076976B (zh) * 2020-08-14 2024-03-26 中国石油天然气股份有限公司 一种利用纵横波速度联合预测有效应力系数的方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104698492A (zh) * 2013-12-09 2015-06-10 中国石油天然气股份有限公司 一种计算异常地层压力的方法
CN105468886A (zh) * 2014-09-04 2016-04-06 中国石油化工股份有限公司 基于岩石物理学参数计算地层压力的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104698492A (zh) * 2013-12-09 2015-06-10 中国石油天然气股份有限公司 一种计算异常地层压力的方法
CN105468886A (zh) * 2014-09-04 2016-04-06 中国石油化工股份有限公司 基于岩石物理学参数计算地层压力的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
THE EFFECTIVE PRESSURE LAW FOR PERMEABILITY DURING PORE PRESSURE AND CONFINING PRESSURE CYCLING OF SEVERAL CRYSTALLINE ROCKS;Yves Bernabe;《JOURNAL OF GEOPHYSICAL RESEARCH》;19870110;第92卷;第649-657页 *
基于有效应力法的碳酸盐岩地层孔隙压力测井计算;夏宏泉等;《钻采工艺》;20050525;第28卷(第3期);第28-30页 *
川西坳陷深层岩石力学特征及其影响因素;高艳霞等;《油气地质与采收率》;20071115;第14卷(第6期);第23-25页 *
常规方法预测碳酸盐岩地层压力的偏差分析;王子振等;《中国石油大学学报( 自然科学版)》;20141020;第38卷(第5期);第96-101页 *
特低渗透砂岩有效应力系数测定;高涛等;《西安石油大学学报( 自然科学版)》;20140725;第29卷(第4期);第46-49、54页 *

Also Published As

Publication number Publication date
CN106960264A (zh) 2017-07-18

Similar Documents

Publication Publication Date Title
CN109034447B (zh) 一种水驱油藏含水上升率变化的预测方法及装置
Ennis-King et al. Onset of convection in anisotropic porous media subject to a rapid change in boundary conditions
CN106896406B (zh) 基于阻抗域属性预测优质储层的方法和装置
CN106960264B (zh) 低渗透砂砾岩地层压力预测方法及装置
CN106991223B (zh) 确定多层油藏水驱波及系数的方法及装置
CN107451311B (zh) 一种油水相对渗透率曲线计算方法及装置
CN108005644B (zh) 一种倾斜煤层动态渗透率预测方法及装置
CN108681793B (zh) 深层油藏采油指数预测方法及装置
CN107656312B (zh) 基于分方位角叠加的河道砂体预测方法及装置
CN108873065A (zh) 砂岩优质储层预测方法及装置
CN108319743B (zh) 古老海相碳酸盐岩油气资源丰度预测方法及装置
CN106772599B (zh) 一种计算地层横波速度的方法及装置
CN107977732B (zh) 边水气藏水平井的见水时间预测方法及装置
Barkve et al. Analysis of reinfiltration in fractured porous media
Park et al. Use of ensemble Kalman filter with 3-dimensional reservoir characterization during waterflooding
CN106908842B (zh) 一种确定单井沉积微相曲线的方法及装置
CN109063228B (zh) 一种水驱油藏含水上升率变化的确定方法及装置
CN109033508B (zh) 一种油藏含水率及含水上升率的确定方法和装置
CN108920781B (zh) 一种油藏含水率及含水上升率的确定方法和装置
CN108956943B (zh) 煤岩割理压缩系数测定方法及装置
CN114139386B (zh) 一种基于多层合采的连通性评估方法及装置
CN111856571B (zh) 测井频段速度向地震频段速度转换方法、装置及设备
CN111381279A (zh) 储层孔隙度定量预测方法及装置
Ahmadinia et al. A history matching approach to estimate caprock morphology parameters for CO2 storage in saline aquifers
CN111856567B (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