CN108072901B - 一种获得准确的静水压力和上覆地层压力的方法及*** - Google Patents
一种获得准确的静水压力和上覆地层压力的方法及*** Download PDFInfo
- Publication number
- CN108072901B CN108072901B CN201610987564.8A CN201610987564A CN108072901B CN 108072901 B CN108072901 B CN 108072901B CN 201610987564 A CN201610987564 A CN 201610987564A CN 108072901 B CN108072901 B CN 108072901B
- Authority
- CN
- China
- Prior art keywords
- pressure
- hydrostatic pressure
- overburden
- depth
- target layer
- 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
Links
- 230000002706 hydrostatic effect Effects 0.000 title claims abstract description 131
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000015572 biosynthetic process Effects 0.000 title claims description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 24
- 238000011160 research Methods 0.000 claims description 30
- 238000005259 measurement Methods 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 239000011435 rock Substances 0.000 claims description 15
- 239000008398 formation water Substances 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 238000009530 blood pressure measurement Methods 0.000 claims description 5
- 238000005553 drilling Methods 0.000 description 6
- 239000003208 petroleum Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Images
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. for interpretation or for event detection
- 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. for interpretation or for event detection
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
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)。
通常采用的从零时间开始积分计算静水压力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)。
为了将深度域测井数据同时间域的地震数据匹配起来,首先需要对深度域测井数据进行时深转换,通过如式(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号(其中CDP是指Common Depth Point,共深度点),纵坐标为时间,图中竖线为井上计算的静水压力和上覆地层压力变密度方式的显示,从图中可以看出,井上计算的静水压力和上覆地层压力颜色与地震剖面计算的完全融合,说明两者计算结果在井点处完全一致,说明基于地震计算的静水压力和上覆地层压力与基于井资料计算的非常吻合。
因此,本发明中的方法,消除了常规基于地震计算中,需要先进行密度预测,由于密度预测通常不准,容易照成累积误差;该方法相较常规方法能够得到更加准确的静水压力和上覆地层压力。
实施例三
如图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:基于获得的目标层的纵波阻抗信息,获得目标层准确的内静水压力和上覆地层压力;
目标层内准确的静水压力的公式如下:
目标层内准确的上覆地层压力的公式如下:
其中,PHY(x,y,t)为目标层内准确的静水压力,POV(x,y,t)为目标层内准确的上覆地层压力,Ip(x,y,t)为纵波阻抗,Vpwater代表水的速度,ρwater指的是水的密度,PHY(x,y,t0)为目标层顶界面静水压力,POV(x,y,t0)为目标层顶界面上覆地层压力。
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为深度。
4.根据权利要求1所述的获得准确的静水压力和上覆地层压力的方法,其特征在于,在所述步骤S300中,采用解释软件中的数据内插或者地质建模的模块,构建研究区目标层顶界面静水压力和目标层顶界面上覆地层压力的等效趋势面。
5.根据权利要求1所述的获得准确的静水压力和上覆地层压力的方法,其特征在于,在所述步骤S400中,采用Jason软件进行叠后波阻抗反演。
6.一种获得准确的静水压力和上覆地层压力的***,该***包括:
压力初步测定模块,所述初步测定模块基于深度域测井数据获得研究区内单井上的静水压力和上覆地层压力;
初步测定扩大模块,所述初步测定扩大模块重复初定测定模块的操作内容,得到研究区内所有井基于深度域测井数据的静水压力和上覆地层压力;
等效趋势面构建模块,所述等效趋势面构建模块基于研究区内已知的井的井震标定获得的时深关系、初步测定扩大模块获得的静水压力和上覆地层压力,构建研究区目标层顶界面静水压力和上覆地层压力的等效趋势面;
纵波阻抗信息模块,所述纵波阻抗信息模块基于叠后地震数据,针对目标层进行叠后波阻抗反演,获得目标层的纵波阻抗信息;以及
准确压力测定模块,所述准确压力测定模块基于纵波阻抗信息模块获得的目标层的纵波阻抗信息,获得目标层准确的内静水压力和上覆地层压力;
目标层内准确的静水压力的公式如下:
目标层内准确的上覆地层压力的公式如下:
其中,PHY(x,y,t)为目标层内准确的静水压力,POV(x,y,t)为目标层内准确的上覆地层压力,Ip(x,y,t)为纵波阻抗,Vpwater代表水的速度,ρwater指的是水的密度,PHY(x,y,t0)为目标层顶界面静水压力,POV(x,y,t0)为目标层顶界面上覆地层压力。
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为深度。
9.根据权利要求6所述的获得准确的静水压力和上覆地层压力的***,其特征在于,所述等效趋势面构建模块中,采用解释软件中的数据内插或者地质建模的模块,构建研究区目标层顶界面静水压力和目标层顶界面上覆地层压力的等效趋势面。
10.根据权利要求6所述的获得准确的静水压力和上覆地层压力的***,其特征在于,所述纵波阻抗信息模块,采用Jason软件进行叠后波阻抗反演。
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 CN108072901A (zh) | 2018-05-25 |
CN108072901B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927790A (zh) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | 基于叠前弹性参数的反演地层压力预测方法及*** |
CN115857011B (zh) * | 2022-12-09 | 2023-06-20 | 中国科学院地质与地球物理研究所 | 一种基于区域最优化拟合求取地震计传递函数的方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101942992B (zh) * | 2010-08-19 | 2013-11-20 | 中国石油大学(北京) | 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 |
CN104769458A (zh) * | 2014-07-15 | 2015-07-08 | 杨顺伟 | 一种基于柯西分布的叠后波阻抗反演方法 |
CN104267429B (zh) * | 2014-09-30 | 2016-11-09 | 中国石油天然气股份有限公司 | 确定地层压力的方法及装置 |
CN104483703B (zh) * | 2014-12-31 | 2017-05-10 | 中国石油天然气股份有限公司 | 一种井震联合古地层压力预测方法 |
CN105445791B (zh) * | 2015-11-25 | 2017-09-19 | 成都理工大学 | 一种基于多种地震属性的地层孔隙压力预测方法 |
CN105626056B (zh) * | 2015-12-22 | 2019-01-18 | 中国石油天然气集团公司 | 一种确定目标区三维地层孔隙压力的方法和装置 |
-
2016
- 2016-11-10 CN CN201610987564.8A patent/CN108072901B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108072901A (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3571532B1 (en) | Systematic evaluation of shale plays | |
US11209561B2 (en) | Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function | |
EP2846175B1 (en) | Seismic survey analysis | |
US7840394B2 (en) | Method for generating a 3D earth model | |
AU2008282547B2 (en) | Method for determining seismic data quality | |
AU2011356685B2 (en) | Methods and systems regarding models of underground formations | |
US10359529B2 (en) | Singularity spectrum analysis of microseismic data | |
AU2013384285B2 (en) | System, method and computer program product for predicting well production | |
CN105089615B (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
NO337423B1 (no) | Fremgangsmåte for oppdatering av en geologisk modell med seismiske data | |
WO2009062935A1 (en) | Forming a geological model | |
UA76697C2 (uk) | Спосіб прогнозування кількісних значень властивості гірської породи або текучого середовища у продуктивній товщі за сейсмічними даними | |
EP2888606A1 (en) | Systems and methods for 3d seismic data depth conversion utilizing artificial neural networks | |
GB2535752A (en) | A method of utilising impedance logs derived from VSP data for use in assessing an oilfield subterranean formation | |
Soleymani et al. | Velocity based pore pressure prediction—A case study at one of the Iranian southwest oil fields | |
CN108072901B (zh) | 一种获得准确的静水压力和上覆地层压力的方法及*** | |
US6982928B2 (en) | Seismic P-wave velocity derived from vibrator control system | |
EP3301485B1 (en) | Crosswell microseismic system | |
EP3929630B1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
CN109991663B (zh) | 工区地震速度体校正方法和装置 | |
US10209380B2 (en) | Methods and systems for juxtaposition across geological discontinuities | |
Floricich et al. | An engineering-driven approach for separating pressure and saturation using 4D seismic: Application to a Jurassic reservoir in the UK North Sea | |
US20240019599A1 (en) | Reservoir properties derived using ultra-deep resistivity inversion data | |
Ponte et al. | Well-Seismic Integration to Pore Pressure Prediction using Multivariate Geostatistics: a Case Study in a Brazilian Equatorial Margin Basin | |
CN115685345A (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 |