CN108072901A - 一种获得准确的静水压力和上覆地层压力的方法及系统 - Google Patents

一种获得准确的静水压力和上覆地层压力的方法及系统 Download PDF

Info

Publication number
CN108072901A
CN108072901A CN201610987564.8A CN201610987564A CN108072901A CN 108072901 A CN108072901 A CN 108072901A CN 201610987564 A CN201610987564 A CN 201610987564A CN 108072901 A CN108072901 A CN 108072901A
Authority
CN
China
Prior art keywords
pressure
overlying formation
hydrostatic pressure
formation pressure
depth
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
CN201610987564.8A
Other languages
English (en)
Other versions
CN108072901B (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 Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201610987564.8A priority Critical patent/CN108072901B/zh
Publication of CN108072901A publication Critical patent/CN108072901A/zh
Application granted granted Critical
Publication of CN108072901B publication Critical patent/CN108072901B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface

Abstract

本发明中的获得准确的静水压力和上覆地层压力的方法及系统，基于井震标定获得的时深关系，结合单井计算的静水压力和上覆地层压力，利用商业软件地质建模的思路构建出目标层顶界面静水压力和上覆地层压力的等效趋势面。然后重新推导了时间域静水压力和上覆地层压力的计算公式，通过叠后波阻抗反演结果能够获得目标层内静水压力和上覆地层压力。本发明所提供的获得准确的静水压力和上覆地层压力的方法及系统可以获得准确的静水压力和上覆地层压力，从而解决了现有技术由于外界各种因素的影响，而难以获得准确的静水压力和上覆地层压力计算结果的技术问题。

Description

一种获得准确的静水压力和上覆地层压力的方法及系统
技术领域
本发明涉及油气地球物理勘探领域，具体的是涉及地层压力预测中的一种获得准确的静水压力和上覆地层压力的方法。
背景技术
地层压力是石油勘探开发工作中的基础数据之一，对异常地层压力(特别是异常高压)的研究，愈来愈引起人们的注意。在石油勘探过程中，地层压力为油气的分布、运移、聚集提供重要信息；在石油钻井工程中，地层压力不仅是确定钻井液密度、井身结构的依据，而且关系到能否安全、快速、经济钻井，甚至影响到钻井的成败。
由于地震资料具有对地层分辨率高、覆盖面积广、资料丰富等优点，地震资料常被用来进行地层压力的钻前预测。常用的地层压力钻前预测方法为Eaton法和Fillippone法。
其中Eaton法的地层压力预测公式表示为：
P＝POV-(POV-PHY)(v/vnct)n
式中，POV为上覆地层压力，PHY为静水压力。
Fillippone法的地层压力预测公式表示为：
式中，POV也是上覆地层压力。
由上述内容可以看出，无论采用何种方法上覆地层压力POV和静水压力PHY的准确估算对于获得可靠的地层压力预测结果都非常关键。
静水压力是指与岩石表面及地表连通的开放体系下的水柱压力。静水压力相当于目的层到水源水柱的垂直高度。静水压力的计算公式为：
PHY＝0.098ρw·h
上式中PHY为静水压力，单位为kg/cm2；ρw为地层水的相对密度，单位是kg/cm3；h是水柱高度，单位为m。
上覆地层压力是指覆盖在某一深度地层以上的地层基岩和岩石孔隙中流体的总重量所造成的对这个地层的压力。上覆地层压力的计算公式为：
POV＝0.098ρb·h
上式中POV为上覆地层压力，单位为kg/cm2；ρb为岩石的体密度，单位是kg/cm3；h是地层的深度，单位为m。对于整个沉积层段来说，沉积岩的体密度不是固定不变的。
从上面静水压力和上覆地层压力的计算公式可以看出，获取岩石的体密度和地层水密度以及获得地层精确深度是计算上覆地层压力和静水压力的关键。但事实上，由于地表起伏的影响，我们很难获得准确的地层深度。另外，通常测井工程师并不是从起始深度开始测量，造成密度测井曲线在起始深度段缺失。这些因素都会给我们获得准确的上覆地层压力和静水压力造成困难。
发明内容
本发明提供一种获得准确的静水压力和上覆地层压力的方法，以解决现有技术中常规地层压力预测中难以获得准确的上覆地层压力和静水压力的技术问题。
针对上述技术问题，本发明提供一种获得准确的静水压力和上覆地层压力的方法，其特征在于，包括以下步骤：
S100：基于深度域测井数据获得研究区内单井上的静水压力和上覆地层压力；
S200：重复步骤S100，得到研究区内所有井的静水压力和上覆地层压力
S300：基于研究区内已知的井的井震标定获得的时深关系、步骤S200中获得的静水压力和上覆地层压力，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面；
S400：基于叠后地震数据，针对目标层进行叠后地震反演叠后波阻抗反演，获得目标层的纵波阻抗信息；
S500：基于获得的目标层的纵波阻抗信息，获得目标层准确的内静水压力和上覆地层压力。
在一个实施例中，在所述步骤S100中，静水压力和上覆地层压力根据下述方法获得：
静水压力计算公式如下：
PHY(z)＝0.098ρw·z
其中PHY(z)为静水压力，ρw为地层水密度，z为深度；
上覆地层压力计算公式如下：
POV(z)＝0.098ρb·z
其中POV(z)为上覆地层压力，ρb为岩石体密度，z为深度；
其特征在于，所述岩石体密度有以下算法获得：
ρb＝ρmudline+A*ZB
其中ρmudline为近地表地层密度，A为第一拟合系数，B为第二拟合系数，z为深度。
在一个实施例中，在所述步骤S300中，所述的时深关系的关系式如下：
其中t为时间；vp为纵波速度；z为深度。
在一个实施例中，在所述步骤S300中，采用解释软件中的数据内插或者地质建模的模块，构建研究区目标层顶界面静水压力和目标层顶界面上覆地层压力的等效趋势面。
在一个实施例中，在所述步骤S400中，采用Jason软件进行叠后波阻抗反演。
本发明的另一方面，还提供一种获得准确的静水压力和上覆地层压力的系统，该系统包括：
压力初步测定模块，所述初步测定模块基于深度域测井数据获得研究区内单井上的静水压力和上覆地层压力；
初步测定扩大模块，所述初步测定扩大模块重复初定测定模块的操作内容，得到研究区内所有井基于深度域测井数据的静水压力和上覆地层压力；
等效趋势面构建模块，所述等势面构建模块基于研究区内已知的井的井震标定获得的时深关系、初步测定扩大模块获得的静水压力和上覆地层压力，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面；
纵波阻抗信息模块，所述纵波阻抗信息模块基于叠后地震数据，针对目标层进行叠后波阻抗反演，获得目标层的纵波阻抗信息；以及
准确压力测定模块，所述准确压力测定模块基于纵波阻抗信息模块获得的目标层的纵波阻抗信息，获得目标层准确的内静水压力和上覆地层压力。
在一个实施例中，所述压力初步测定模块中的静水压力和上覆地层压力根据下述方法获得：
静水压力计算公式如下：
PHY(z)＝0.098ρw·z
其中PHY(z)为静水压力，ρw为地层水密度，z为深度；
上覆地层压力计算公式如下：
POV(z)＝0.098ρb·z
其中POV(z)为上覆地层压力，ρb为岩石体密度，z为深度；
其特征在于，所述岩石体密度有以下算法获得：
ρb＝ρmudline+A*ZB
其中ρmudline为近地表地层密度，A为第一拟合系数，B为第二拟合系数，Z为深度。
在一个实施例中，所述等效趋势面构建模块中的时深关系的关系式如下：
其中t为时间，vp为纵波速度，z为深度。
在一个实施例中，所述等效趋势面构建模块中，采用解释软件中的数据内插或者地质建模的模块，构建研究区目标层顶界面静水压力和目标层顶界面上覆地层压力的等效趋势面。
在一个实施例中，所述纵波阻抗信息模块，采用Jason软件进行叠后波阻抗反演。
本发明中的获得准确的静水压力和上覆地层压力的方法及系统，基于井震标定获得的时深关系，结合单井计算的静水压力和上覆地层压力，利用商业软件地质建模的思路构建出目标层顶界面静水压力和上覆地层压力的等效趋势面。然后重新推导了时间域静水压力和上覆地层压力的计算公式，通过叠后波阻抗反演结果能够获得目标层内静水压力和上覆地层压力。本发明所提供的获得准确的静水压力和上覆地层压力的方法，解决了在现有技术中，静水压力和上覆地层压力的获得过程中，由于地表起伏的影响，我们很难获得准确的地层深度；通常测井工程师并不是从起始深度开始测量，造成密度测井曲线在起始深度段缺失；常规基于密度资料从零时间开始积分的方法容易造成累积误差，难以获得准确的静水压力和上覆地层压力计算结果等技术问题。根据本发明所提供的方法和系统可以获得更加精确的静水压力和上覆地层压力，为进一步的获得可靠地层压力预测结果打好基础。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
图1是本发明实施例一中的获得准确的静水压力和上覆地层压力的方法的流程图；
图2是本发明实施例二中工区目的层段静水压力预测结果图；
图3是本发明实施例二中工区目的层段上覆地层压力预测结果图；
图4是本发明实施例三中的获得准确的静水压力和上覆地层压力的系统的结构框图。
在附图中，相同的部件使用相同的附图标记。附图并未按照实际的比例。
具体实施方式
下面将结合附图对本发明作进一步说明。
本发明基于深度域测井数据Well(z)(这里z代表深度域，单位为m或者ft)和时间域地震数据Seis(x，y，t)(这里x，y代表每个地震道对应的空间位置，可以是x，y坐标，也可以是纵横测线号，t代表时间域，单位为s或者ms)，通过以下步骤实现对于上覆地层压力POV(x，y，t)和静水压力PHY(x，y，t)准确预测：
(i)基于深度域测井数据计算研究区内，单井上的静水压力PHY(z)和上覆地层压力POV(z)。
PHY(z)＝0.098ρw·z (1)
其中ρw为地层水密度，z为深度；
POV(z)＝0.098ρb·z (2)
其中ρb为岩石体密度，z为深度；
通常测井工程师并不是从起始深度开始测量，造成密度测井曲线在起始深度段缺失，本发明中，采用Amoco公式，如公式(3)计算密度，
ρb＝ρmudline+A*ZB (3)
其中参数ρmudline，A，B由用户自己给出。
(ii)重复步骤(i)得到研究区内所有井的静水压力PHY(z)和上覆地层压力POV(z)。
(iii)基于研究区已知井井震标定获得的时深关系(基于常规地震解释软件可以获得，比如Jason、HRS等软件均可以获得)和步骤(ii)得到的静水压力PHY(z)和上覆地层压力POV(z)，利用解释软件中用于数据内插或者地质建模的模块，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面PHY(x，y，t0)和POV(x，y，t0)。其中t0为目标层顶界面层位时间属性。
在本实施例中采用的是Jason软件Model generator模块获得。输入数据为步骤步骤(ii)得到的静水压力和上覆地层压力、工区地层格架(包含目标层顶界面层位数据horizon(x，y，t0))，采用克里金插值算法获取工区究区目标层顶界面静水压力和上覆地层压力的等效趋势面PHY(x，y，t0)和POV(x，y，t0)。具体获得方法，参考Jason软件Modelgenerator模块使用说明书。
通常从零时间开始积分计算静水压力PHY(x，y，t)和上覆地层压力POV(x，y，t)的方法，容易造成误差累积。通过本发明步骤(iii)获得目标层顶界面静水压力PHY(x，y，t0)和上覆地层压力POV(x，y，t0)的等效趋势面之后，可以保证在目标层顶界面获得完全准确的静水压力和上覆地层压力。
因此，我们的问题变成针对目标层内的静水压力PHY(x，y，t)和上覆地层压力POV(x，y，t)的精确计算。
(iv)基于叠后地震数据体Seis(x，y，t)，通过叠后波阻抗反演，获取纵波阻抗信息Ip(x，y，t)。
(v)基于获得的波阻抗以及基于步骤(i)至步骤(iii)中的其他数据，进行目标层内静水压力PHY(x，y，t)和上覆地层压力POV(x，y，t)精确计算。
目前，基于叠后地震数据通过地震反演的方法获取纵波阻抗的方法有很多种，商业软件业比较成熟，比如Jason，Strata等软件。通过选择合适的子波、低频模型以及合适的反演参数，针对目标层开展叠后地震反演，都能够获得目标层内纵向分辨率较高的纵波阻抗Ip(x，y，t)。在本实施例中采用的是Jason软件InverTrace-Plus模块获得。采用稀疏脉冲反演算法获取工区纵波阻抗信息Ip(x，y，t)。具体获得方法，参考Jason软件InverTrace-Plus模块使用说明书。
为了将深度域测井数据同时间域的地震数据匹配起来，首先需要对深度域测井数据进行时深转换，通过如式(4)所示，将深度域测井数据转换到与地震匹配的时间域数据
深度域内计算静水压力的公式为式(5)，深度域内计算上覆地层压力的公式为式(6)：
将时深转换关系式(4)，带入式(5)和式(6)，
得到时间域计算静水压力PHY(t)的公式为式(7)，和时间域上覆地层压力POV(t)计算公式为式(8)：
将目标层内叠后地震反演获得的纵波阻抗Ip(x，y，t)带入式(7)和式(8)，即可以得到目标层内精确的静水压力PHY(x，y，t)的公式为式(9)和上覆地层压力POV(x，y，t)的公式为式(10)：
其中：是Vpwater代表水的速度，ρwater指的水的数据
实施例一
如图1所示，为本实施例中获得准确的静水压力和上覆地层压力的方法的流程图，本实施例中的获得准确的静水压力和上覆地层压力的方法，包括以下步骤：
S100：基于深度域测井数据获得研究区内单井上的静水压力和上覆地层压力；
S200：重复步骤S100，得到研究区内所有井的静水压力和上覆地层压力
S300：基于研究区内已知的井的井震标定获得的时深关系、步骤S200中获得的静水压力和上覆地层压力，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面；
S400：基于叠后地震数据，针对目标层进行叠后波阻抗反演，获得目标层的纵波阻抗信息；
S500：基于获得的目标层的纵波阻抗信息，获得目标层准确的内静水压力和上覆地层压力。
实施例二
如图2所示，根据本发明的实施例获取目标层内静水压力和上覆地层压力的方法包括以下步骤：
步骤S100，基于深度域测井数据计算单井上的静水压力PHY(z)和上覆地层压力POV(z)：
PHY(z)＝0.098ρw·z，
其中ρw为地层水密度，本实施例中给定位1.08g/cm3；
POV(z)＝0.098ρb·z，
其中ρb为岩石体密度，
通常测井工程师并不是从起始深度开始测量，造成密度测井曲线在起始深度段缺失，采用Amoco公式计算密度，ρb＝ρmudline+A*ZB，其中参数ρmudline，A，B由分别为1.95，0.156，0.159，在有实测密度曲线的部分，用实测密度曲线计算上覆地层压力。
步骤S200，重复S100计算不同井的静水压力和上覆地层压力。
步骤S300，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面，需要用到步骤S100中获得单井计算的静水压力和上覆地层压力、S200中获得多井计算的静水压力和上覆地层压力以及井震标定的时深关系，利用jason软件EarthModel模块地质建模的思路构建出目标层顶界面静水压力和上覆地层压力的等效趋势面PHY(x，y，t0)和POV(x，y，t0)。
步骤S400，基于叠后地震数据体Seis(x，y，t)，通过Jason软件开展叠后波阻抗反演，获取纵波阻抗信息Ip(x，y，t)；
步骤S500，基于获得的波阻抗进行目标层内静水压力PHY(x，y，t)和上覆地层压力POV(x，y，t)精确计算。
根据上述的计算得到如图2和图3所示的结果：
其中图2为静水压力计算结果，图3为上覆地层压力计算结果，从图中可以看出，基于地震计算的静水压力和上覆地层压力与基于井资料计算的非常吻合。
图2和图3中，横坐标为CDP号，纵坐标为时间，图中竖线为井上计算的静水压力和上覆地层压力变密度方式的显示，从图中可以看出，井上计算的静水压力和上覆地层压力颜色与地震剖面计算的完全融合，说明两者计算结果在井点处完全一致，说明基于地震计算的静水压力和上覆地层压力与基于井资料计算的非常吻合；
因此，本发明中的方法，消除了常规基于地震计算中，需要先进行密度预测，由于密度预测通常不准，容易照成累积误差的缺陷；该方法相较常规方法能够得到更加准确的静水压力和上覆地层压力。
实施例三
如图4所示为本发明中的一种获得准确的静水压力和上覆地层压力的系统的结构框图，该系统包括：
压力初步测定模块10，所述初步测定模块10基于深度域测井数据获得研究区内单井上的静水压力和上覆地层压力；
初步测定扩大模块20，所述初步测定扩大模块20重复初定测定模块的操作内容，得到研究区内所有井基于深度域测井数据的静水压力和上覆地层压力；
等效趋势面构建模块30，所述等势面构建模块30基于研究区内已知的井的井震标定获得的时深关系、初步测定扩大模块获得的静水压力和上覆地层压力，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面；
纵波阻抗信息模块40，所述纵波阻抗信息模块40基于叠后地震数据，针对目标层进行叠后波阻抗反演获得目标层的纵波阻抗信息；以及
准确压力测定模块50，所述准确压力测定模块50基于纵波阻抗信息模块获得的目标层的纵波阻抗信息，获得目标层准确的内静水压力和上覆地层压力。
其中压力初步测定模块10、初步测定扩大模块20、等效趋势面构建模块30、纵波阻抗信息模块40和准确压力测定模块50，可以分别执行实施例一和实施二中的步骤S100、S200、S300、S400和S500的流程，在此不再赘述。
本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (10)

1.一种获得准确的静水压力和上覆地层压力的方法，其特征在于，包括以下步骤：
S100：基于深度域测井数据获得研究区内单井上的静水压力和上覆地层压力；
S200：重复步骤S100，得到研究区内所有井基于深度域测井数据的静水压力和上覆地层压力
S300：基于研究区内已知的井的井震标定获得的时深关系、步骤S200中获得的静水压力和上覆地层压力，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面；
S400：基于叠后地震数据，针对目标层进行叠后波阻抗反演，获得目标层的纵波阻抗信息；
S500：基于获得的目标层的纵波阻抗信息，获得目标层准确的内静水压力和上覆地层压力。
2.根据权利要求1所述的获得准确的静水压力和上覆地层压力的方法，在所述步骤S100中，静水压力和上覆地层压力根据下述方法获得：
静水压力计算公式如下：
PHY(z)＝0.098ρw·z
其中PHY(z)为静水压力，ρw为地层水密度，z为深度；
上覆地层压力计算公式如下：
POV(z)＝0.098ρb·z
其中POV(z)为上覆地层压力，ρb为岩石体密度，z为深度；
其特征在于，所述岩石体密度有以下算法获得：
ρb＝ρmudline+A*ZB
其中ρmudline为近地表地层密度，A为第一拟合系数，B为第二拟合系数，z为深度。
3.根据权利要求1所述的获得准确的静水压力和上覆地层压力的方法，其特征在于，在所述步骤S300中，所述的时深关系的关系式如下：
<mrow> <mi>t</mi> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>z</mi> </msubsup> <mfrac> <mrow> <mn>2</mn> <mi>d</mi> <mi>z</mi> </mrow> <mrow> <mi>v</mi> <mi>p</mi> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
其中t为时间，vp为纵波速度，z为深度。
4.根据权利要求1所述的获得准确的静水压力和上覆地层压力的方法，其特征在于，在所述步骤S300中，采用解释软件中的数据内插或者地质建模的模块，构建研究区目标层顶界面静水压力和目标层顶界面上覆地层压力的等效趋势面。
5.根据权利要求1所述的获得准确的静水压力和上覆地层压力的方法，其特征在于，在所述步骤S400中，采用Jason软件进行叠后波阻抗反演。
6.一种获得准确的静水压力和上覆地层压力的系统，该系统包括：
压力初步测定模块，所述初步测定模块基于深度域测井数据获得研究区内单井上的静水压力和上覆地层压力；
初步测定扩大模块，所述初步测定扩大模块重复初定测定模块的操作内容，得到研究区内所有井基于深度域测井数据的静水压力和上覆地层压力；
等效趋势面构建模块，所述等势面构建模块基于研究区内已知的井的井震标定获得的时深关系、初步测定扩大模块获得的静水压力和上覆地层压力，构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面；
纵波阻抗信息模块，所述纵波阻抗信息模块基于叠后地震数据，针对目标层进行叠后波阻抗反演，获得目标层的纵波阻抗信息；以及
准确压力测定模块，所述准确压力测定模块基于纵波阻抗信息模块获得的目标层的纵波阻抗信息，获得目标层准确的内静水压力和上覆地层压力。
7.根据权利要求6所述的获得准确的静水压力和上覆地层压力的系统，其特征在于，所述压力初步测定模块中的静水压力和上覆地层压力根据下述方法获得：
静水压力计算公式如下：
PHY(z)＝0.098ρw·z
其中PHY(z)为静水压力，ρw为地层水密度，z为深度；
上覆地层压力计算公式如下：
POV(z)＝0.098ρb·z
其中POV(z)为上覆地层压力，ρb为岩石体密度，z为深度；
其特征在于，所述岩石体密度有以下算法获得：
ρb＝ρmudline+A*ZB
其中ρmudline为近地表地层密度，A为第一拟合系数，B为第二拟合系数，z为深度。
8.根据权利要求6所述的获得准确的静水压力和上覆地层压力的系统，其特征在于，所述等效趋势面构建模块中的时深关系的关系式如下：
<mrow> <mi>t</mi> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>z</mi> </msubsup> <mfrac> <mrow> <mn>2</mn> <mi>d</mi> <mi>z</mi> </mrow> <mrow> <mi>v</mi> <mi>p</mi> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
其中t为时间，vp为速度，z为深度。
9.根据权利要求6所述的获得准确的静水压力和上覆地层压力的系统，其特征在于，所述等效趋势面构建模块中，采用解释软件中的数据内插或者地质建模的模块，构建研究区目标层顶界面静水压力和目标层顶界面上覆地层压力的等效趋势面。
10.根据权利要求6所述的获得准确的静水压力和上覆地层压力的系统，其特征在于，所述纵波阻抗信息模块，采用Jason软件进行叠后波阻抗反演。
CN201610987564.8A 2016-11-10 2016-11-10 一种获得准确的静水压力和上覆地层压力的方法及系统 Active CN108072901B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610987564.8A CN108072901B (zh) 2016-11-10 2016-11-10 一种获得准确的静水压力和上覆地层压力的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610987564.8A CN108072901B (zh) 2016-11-10 2016-11-10 一种获得准确的静水压力和上覆地层压力的方法及系统

Publications (2)

Publication Number Publication Date
CN108072901A true CN108072901A (zh) 2018-05-25
CN108072901B CN108072901B (zh) 2020-12-01

Family

ID=62153345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610987564.8A Active CN108072901B (zh) 2016-11-10 2016-11-10 一种获得准确的静水压力和上覆地层压力的方法及系统

Country Status (1)

Country Link
CN (1) CN108072901B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110927790A (zh) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 基于叠前弹性参数的反演地层压力预测方法及系统
CN115857011A (zh) * 2022-12-09 2023-03-28 中国科学院地质与地球物理研究所 一种基于区域最优化拟合求取地震计传递函数的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101942992A (zh) * 2010-08-19 2011-01-12 中国石油大学(北京) 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法
CN104267429A (zh) * 2014-09-30 2015-01-07 中国石油天然气股份有限公司 确定地层压力的方法及装置
CN104483703A (zh) * 2014-12-31 2015-04-01 中国石油天然气股份有限公司 一种井震联合古地层压力预测方法
WO2016008105A1 (zh) * 2014-07-15 2016-01-21 杨顺伟 一种基于柯西分布的叠后波阻抗反演方法
CN105445791A (zh) * 2015-11-25 2016-03-30 成都理工大学 一种基于多种地震属性的地层孔隙压力预测方法
CN105626056A (zh) * 2015-12-22 2016-06-01 中国石油天然气集团公司 一种确定目标区三维地层孔隙压力的方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101942992A (zh) * 2010-08-19 2011-01-12 中国石油大学(北京) 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法
WO2016008105A1 (zh) * 2014-07-15 2016-01-21 杨顺伟 一种基于柯西分布的叠后波阻抗反演方法
CN104267429A (zh) * 2014-09-30 2015-01-07 中国石油天然气股份有限公司 确定地层压力的方法及装置
CN104483703A (zh) * 2014-12-31 2015-04-01 中国石油天然气股份有限公司 一种井震联合古地层压力预测方法
CN105445791A (zh) * 2015-11-25 2016-03-30 成都理工大学 一种基于多种地震属性的地层孔隙压力预测方法
CN105626056A (zh) * 2015-12-22 2016-06-01 中国石油天然气集团公司 一种确定目标区三维地层孔隙压力的方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王玉柱 等: "华南古生界页岩储层压力预测方法及其应用研究", 《煤炭学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110927790A (zh) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 基于叠前弹性参数的反演地层压力预测方法及系统
CN115857011A (zh) * 2022-12-09 2023-03-28 中国科学院地质与地球物理研究所 一种基于区域最优化拟合求取地震计传递函数的方法

Also Published As

Publication number Publication date
CN108072901B (zh) 2020-12-01

Similar Documents

Publication Publication Date Title
US11519262B2 (en) Systematic evaluation of shale plays
US10386513B2 (en) Seismic survey analysis
US9217802B2 (en) Seismic image enhancement
CN106837297A (zh) 一种识别井间连通性及油水动态预测的方法
US10359529B2 (en) Singularity spectrum analysis of microseismic data
NO20160203A1 (en) Stratigraphic function
EP3400463A1 (en) Amplitude inversion on partitioned depth image gathers using point spread function
US20130345985A1 (en) Seismic orthogonal decomposition attribute
CN105386756B (zh) 一种应用应变量计算脆性地层孔隙度的方法
CN104656136A (zh) 基于实际模型地震模拟指导下的油气储层低频阴影识别技术
CN109667573A (zh) 三维页岩储层孔隙压力预测方法、装置和电子设备
CN105938203A (zh) 一种储层特性的检测方法及装置
Li et al. Geostatistical models for shales in distributary channel point bars (Ferron Sandstone, Utah): From ground-penetrating radar data to three-dimensional flow modeling
CN106285642A (zh) 一种基于地震资料的地层孔隙压力预测方法
Yusof et al. Uncertainty analysis of hydrocarbon in place calculation using 3D seismic and well data during appraisal stage–Case study of Goldie Field, offshore Sarawak
CN108072901A (zh) 一种获得准确的静水压力和上覆地层压力的方法及系统
US11892579B2 (en) Crosswell microseismic system
Du et al. Integrated shale gas reservoir modeling
GB2523460A (en) Singularity spectrum analysis of microseismic data
Abbaszadeh et al. Development of an integrated reservoir model for a naturally fractured volcanic reservoir in China
Babasafari et al. Practical workflows for monitoring saturation and pressure changes from 4D seismic data: A case study of Malay Basin
CN107976711A (zh) 一种地震信息定量预测烃源岩有机碳质量分数方法及装置
Alam et al. Interpretation of geologic facies in seismic volume using key rock elastic properties and high-definition facies templates
CN104316959B (zh) 一种基于等效流体声波阻抗的流体识别方法
Doyle et al. The" Gypsyy" Field Research Program in Integrated Reservoir Characterization

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