CN105700017A - 一种确定油气分布数据的方法和装置 - Google Patents

一种确定油气分布数据的方法和装置 Download PDF

Info

Publication number
CN105700017A
CN105700017A CN201610159384.0A CN201610159384A CN105700017A CN 105700017 A CN105700017 A CN 105700017A CN 201610159384 A CN201610159384 A CN 201610159384A CN 105700017 A CN105700017 A CN 105700017A
Authority
CN
China
Prior art keywords
pore pressure
formation pore
formation
interval velocity
seismic
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
CN201610159384.0A
Other languages
English (en)
Other versions
CN105700017B (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 Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201610159384.0A priority Critical patent/CN105700017B/zh
Publication of CN105700017A publication Critical patent/CN105700017A/zh
Application granted granted Critical
Publication of CN105700017B publication Critical patent/CN105700017B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters
    • G01V2210/6248Pore pressure

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (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

一种确定油气分布数据的方法和装置
技术领域
本发明涉及地球物理勘探技术领域,尤其涉及一种确定油气分布数据的方法和装置。
背景技术
在地球物理勘探技术领域,地层孔隙压力是油气勘探开发工作中的基础数据之一。地层孔隙压力不仅能为钻井工程提供参考数据,关系钻井安全,同时由于欠压实、构造运动特别是生烃作用和流体充注的影响可以形成地层异常高压,多个目标区已表明油气的产出与异常高压相关,因此,准确的确定出地层孔隙压力对油气勘探开发有着十分重要的意义。
现有技术中确定的目标开发区域的油气分布数据,需要先确定目标区的空间地层孔隙压力,一般的,空间地层孔隙压力可以利用测井数据和地震数据结合高压低速理论确定。但现有技术中依据高压低速理论即地层孔隙压力和地震层速度呈负相关关系,而砂砾岩储集层的速度高于泥岩地层的速度,由此计算得到的储集层的地层孔隙压力纵向上始终为低压,与实际情况不符,无法真实反映储集层的地层孔隙压力,导致确定出的油气分布数据精度低,无法有效为下一步进行油气勘探预测工作提供帮助。
发明内容
本申请的目的是提供一种确定油气分布数据的方法和装置,可以更准确的确定出目标储集层的油气分布数据,从而可以为油气勘探预测工作提供数据支持。
本申请提供的确定油气分布数据的方法和装置是这样实现的:
一种确定油气分布数据的方法,所述方法包括:
基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;
利用所述目标储集层的地震波阻抗反演体结合加德纳公式计算得到所述目标储集层的第一地震层速度体;
利用所述目的层的地震层速度计算得到目的层的埋藏深度速度数据体;
利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;
基于确定出的所述目标区的地震层速度与地层孔隙压力关系,由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体;
确定单井地层孔隙压力曲线数据,利用所述单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力;
根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据。
一种确定油气分布数据的装置,所述装置包括:
地震波阻抗反演体确定模块,用于基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;
第一计算模块,用于利用所述目标储集层的地震波阻抗反演体结合加德纳公式计算得到所述目标储集层的第一地震层速度体;
第二计算模块,用于利用所述目的层的地震层速度计算得到目的层的埋藏深度速度数据体;
校正处理模块,用于利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;
地层孔隙压力体确定模块,用于基于确定出的所述目标区的地震层速度与地层孔隙压力关系,由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体;
压力曲线数据确定模块,用于确定单井地层孔隙压力曲线数据;
空间地层孔隙压力确定模块,用于利用所述单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力;
油气分布数据确定模块,用于根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据。
本申请基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;可以实现分岩性确定地震波阻抗反演体。然后,结合加德纳公式计算得到所述目标储集层的第一地震层速度体;利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;并由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体,保证利用单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力的准确性;最后根据所述空间地层孔隙压力中地层孔隙压力分布情况数据可以准确的确定出所述目标区的目标储集层的油气分布数据,实现通过圈定异常高压区进行油气预测。与现有技术相比,利用本申请实施例提供的技术方案可以更准确的确定出目标储集层的油气分布数据,从而可以为油气勘探预测工作提供数据支持。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请提供的确定油气分布数据的方法的一种实施例的流程图;
图2是本申请提供的确定出所述目标区的单井地层孔隙压力计算公式的一种实施例的流程示意图;
图3是分别利用伊顿法、比值法和等效深度法求取的地层孔隙压力与实测地层孔隙压力一种实施例的比较示意图;
图4是采用不同伊顿指数的伊顿法确定出的地层孔隙压力与实测地层孔隙压力的一种实施例的比较示意图;
图5是本申请提供确定所述目标区的目标储集层的油气分布数据一种实施例的流程示意图;
图6是本申请确定油气分布数据的装置一种实施例的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面以几个具体的例子详细说明本申请实施例的具体实现。
以下首先介绍本申请一种确定油气分布数据的方法的一种实施例。图1是本申请提供的确定油气分布数据的方法的一种实施例的流程图,本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不表示唯一的执行顺序。在实际中的系统或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。具体的,如图1所示,所述方法可以包括:
S110:基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体。
本申请实施例中可以基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体。具体的,所述目标区的地震波阻抗反演体可以包括但不限于递归反演、稀疏脉冲反演和基于模型的反演等方法反演获得的。所述目标储集层可以包括所述目的层中的储集层。
在实际应用中,不同区域对应的储集层也不一样。以上述准噶尔盆地西北缘玛南斜坡区为目标区为例,所述目标储集层可以是砂砾岩储集层,一般埋藏深度为3000m至3500m。具体的,所述岩性敏感参数可以包括所述目标储集层与其他地层差异较大的岩性参数。具体的,比如砂砾岩储集层的岩性敏感参数可以包括波阻抗和电阻率数据。
S120:利用所述目标储集层的地震波阻抗反演体结合加德纳公式计算得到所述目标储集层的第一地震层速度体。
本申请实施例中,在步骤S110获得所述目标储集层的地震波阻抗反演体之后,可以利用所述目标储集层的地震波阻抗反演体结合加德纳公式计算得到所述目标储集层的第一地震层速度体。具体的,所述加德纳(Gardner)公式可以如下所示:
ρ=a(VP)b
上式中,ρ表示地层密度,单位为g/m3;VP表示第一地震层速度体,所述第一地震层速度体可以表示所述目的层的地震层速度进行分岩性处理后的目标储集层的地震层速度体,单位为m/s;a、b为常数,可以根据目标区具体情况确定。一般的a=0.31,b=1/4。
S130:利用所述目的层的地震层速度计算得到目的层的埋藏深度速度数据体。
本申请实施例中,可以利用所述目的层的地震层速度计算得到目的层的埋藏深度速度数据体。具体的,可以根据地震层速度体与埋深关系计算得到的与埋藏深度相对应的地震层速度体,所述相对应的地震层速度体作为目的层的埋藏深度速度数据体。
进一步的,所述目的层的地震层速度可以利用迪克斯公式结合叠加速度计算得到。
S140:利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体。
本申请实施例中,在步骤S130获得埋藏深度速度数据体之后,可以利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体。具体的,所述埋藏深度校正处理包括采用如下计算公式:
vih=viL-(vif-vimin)
上式中,vih表示第二地震层速度体,所述第二地震层速度体表示所述第一地震层速度体进行埋藏深度校正处理后的地震层速度体;viL表示第一地震层速度体;vif表示埋藏深度速度数据体,所述埋藏深度速度数据体根据地震层速度体与埋深关系计算得到的与埋藏深度相对应的地震层速度体;vimin表示所述目的层顶界面的地震层速度。
S150:基于确定出的所述目标区的地震层速度与地层孔隙压力关系,由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体。
本申请实施例中,可以基于确定出的所述目标区的地震层速度与地层孔隙压力关系,由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体。具体的,所述确定出的所述目标区的地震层速度与地层孔隙压力关系可以包括:
利用所述第二地震层速度体中的参考井的第二地震层速度和所述参考井的单井地层孔隙压力确定出的所述目标区地震层速度与地层孔隙压力关系。
具体的,所述参考井可以包括所述目标区中获得第二地震层速度和单井地层孔隙压力的井。
S160:确定单井地层孔隙压力曲线数据,利用所述单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力。
本申请实施例中,可以确定单井地层孔隙压力曲线数据,利用所述单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力。这里通过对目标储集层的地层孔隙压力体进行井约束的校正处理之后,有效提高了目标储集层的空间地层孔隙压力的准确性。
具体的,所述确定单井地层孔隙压力曲线数据可以包括:
基于目标区参考井的地层孔隙压力进行多方法多参数对比拟合分析处理,确定出所述目标区的单井地层孔隙压力计算公式。
利用所述单井地层孔隙压力计算公式计算得到所述目标区的单井地层孔隙压力。
根据所述目标区的单井地层孔隙压力确定出单井地层孔隙压力曲线数据。
具体的,图2是本申请提供的确定出所述目标区的单井地层孔隙压力计算公式的一种实施例的流程示意图,具体的,所述基于目标区参考井的地层孔隙压力进行多方法多参数对比拟合分析处理,确定出所述目标区的单井地层孔隙压力计算公式可以包括:
S161:利用多种地层孔隙压力计算方法对所述目标区参考井的声波时差测井数据和密度测井数据进行地层孔隙压力计算。
具体的,在实际应用中,所述多种地层孔隙压力计算方法可以包括但不限于伊顿法、比值法和等效深度法等。
S162:将所述地层孔隙压力计算的计算结果与所述参考井的地层孔隙压力进行比对拟合处理,确定出所述比对拟合处理的最大拟合度值所对应地层孔隙压力计算方法。
具体的,所述拟合度值可以包括利用一定算法计算得到的反映地层孔隙压力之间的拟合程度的数值。一般的拟合程度越高,所述拟合度值越大。
在一个具体的实施例中,以准噶尔盆地西北缘玛南斜坡区为目标区,分别利用伊顿法、比值法和等效深度法对所述目标区参考井的声波时差测井数据和密度测井数据求取上覆地层压力和泥岩正常压实趋势线后得到的地层孔隙压力。如图3所示,图3是分别利用伊顿法、比值法和等效深度法求取的地层孔隙压力与实测地层孔隙压力一种实施例的比较示意图。其中,310、330、350实测地层孔隙压力对应的曲线,320、340、360依次为伊顿法、比值法和等效深度法对应求得的地层孔隙压力所对应的曲线。从图中可见,针对所述目标区,伊顿法求取地层孔隙压力与实测的地层孔隙压力之间的拟合相对较好。相应的,可以确定出所述比对拟合处理的最大拟合度值所对应地层孔隙压力计算方法为伊顿法。
S163:对所述最大拟合度值所对应地层孔隙压力计算方法进行参数调节处理,确定出所述参数调节处理对应地层孔隙压力与所述参考井的地层孔隙压力进行比对拟合处理的最大拟合度值所对应参数。
在一个具体的实施例中,以上述准噶尔盆地西北缘玛南斜坡区为目标区为例,可以对伊顿法中的参数进行参数调节处理,相应的,这里的参数可以是伊顿指数。具体的,假设将所述伊顿指数依次设置为0.3、0.65、1.2后求取相应的地层孔隙压力。如图4所示,图4是采用不同伊顿指数的伊顿法确定出的地层孔隙压力与实测地层孔隙压力的一种实施例的比较示意图,其中,410、430、450实测地层孔隙压力对应的曲线,420、440、460依次为伊顿法的伊顿指数为0.3、0.65、1.2对应求得的地层孔隙压力所对应的曲线。从图中可见,伊顿指数为0.65时,伊顿法求取地层孔隙压力与实测的地层孔隙压力之间的拟合相对较好。
S164:当所述最大拟合度值所对应地层孔隙压力计算方法的参数为所述所对应参数时,将所述地层孔隙压力计算方法所对应计算公式作为所述单井地层孔隙压力计算公式。
具体的,以上述准噶尔盆地西北缘玛南斜坡区为目标区为例,所述单井地层孔隙压力计算公式可以如下所示:
P = P o - ( P o - P w ) ( V i V n o r m a l ) 0.65
上式中,P表示目标区的地层孔隙压力,单位为MPa;Vi表示地震层速度,单位为m/s;Po表示上覆地层压力,单位为MPa;Pw表示静水压力,单位为MPa;Vnormal表示正常压实趋势速度,单位为m/s;0.65表示伊顿参数。
因此,本申请实施例进行多方法多参数对比拟合分析处理之后,可以准确的确定出所述目标区的单井地层孔隙压力计算公式,可以有效保证了所述目标区的单井地层孔隙压力的准确性,同时为后续空间地层孔隙压力的确定提供数据支持。
S170:根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据。
本申请实施例中,在步骤S160得到所述目标储集层的空间地层孔隙压力之后,可以根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据。如图5所示,图5是本申请提供确定所述目标区的目标储集层的油气分布数据一种实施例的流程示意图,具体的,可以包括:
S171:判断所述空间地层孔隙压力中地层孔隙压力数值大于预设地层孔隙压力数值的百分比,将所述百分比作为所述目标区的目标储集层的含油气百分比。
具体的,所述预设地层孔隙压力数值可以预先设置,优选的,所述方法还可以包括:
根据所述目标区的单井地层孔隙压力确定含油气储集层的地层孔隙压力数值,将所述地层孔隙压力数值作为预设地层孔隙压力数值。
S172:根据所述含油气百分比和所述目的层的厚度确定所述油气分布范围和油气厚度,将所述油气分布范围和油气厚度作为油气分布数据。
由此可见,本申请一种确定油气分布数据的方法的实施例基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;可以实现分岩性确定地震波阻抗反演体。然后,结合加德纳公式计算得到所述目标储集层的第一地震层速度体;利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;并由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体,保证利用单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力的准确性;最后根据所述空间地层孔隙压力中地层孔隙压力分布情况数据可以准确的确定出所述目标区的目标储集层的油气分布数据,实现通过圈定异常高压区进行油气预测。与现有技术相比,利用本申请实施例提供的技术方案可以更准确的确定出目标储集层的油气分布数据,从而可以为油气勘探预测工作提供数据支持。
本申请另一方面还提供一种确定油气分布数据的装置,如图6所示,图6是本申请确定油气分布数据的装置一种实施例的结构示意图。所述装置600可以包括:
地震波阻抗反演体确定模块610,可以用于基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;
第一计算模块620,可以用于利用所述目标储集层的地震波阻抗反演体结合加德纳公式计算得到所述目标储集层的第一地震层速度体;
第二计算模块630,可以用于利用所述目的层的地震层速度计算得到目的层的埋藏深度速度数据体;
校正处理模块640,可以用于利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;
地层孔隙压力体确定模块650,可以用于基于确定出的所述目标区的地震层速度与地层孔隙压力关系,由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体;
压力曲线数据确定模块660,可以用于确定单井地层孔隙压力曲线数据;
空间地层孔隙压力确定模块670,可以用于利用所述单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力;
油气分布数据确定模块680,可以用于根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据。
在一个优选的实施例中,所述油气分布数据确定模块680可以包括:
第一数据处理单元,可以用于判断所述空间地层孔隙压力中地层孔隙压力数值大于预设地层孔隙压力数值的百分比,将所述百分比作为所述目标区的目标储集层的含油气百分比;
第二数据处理单元,可以用于根据所述含油气百分比和所述目的层的厚度确定所述油气分布范围和油气厚度,将所述油气分布范围和油气厚度作为油气分布数据。
在一个优选的实施例中,所述装置600还可以包括:
预设地层孔隙压力数值确定模块,可以用于根据所述目标区的单井地层孔隙压力确定含油气储集层的地层孔隙压力数值,将所述地层孔隙压力数值作为预设地层孔隙压力数值。
在一个优选的实施例中,所述压力曲线数据确定模块660可以包括:
压力计算公式确定单元,可以用于基于目标区参考井的地层孔隙压力进行多方法多参数对比拟合分析处理,确定出所述目标区的单井地层孔隙压力计算公式;
第一计算单元,可以用于利用所述单井地层孔隙压力计算公式计算得到所述目标区的单井地层孔隙压力;
压力曲线数据确定单元,可以用于根据所述目标区的单井地层孔隙压力确定出单井地层孔隙压力曲线数据。
在一个优选的实施例中,所述压力计算公式确定单元可以包括:
第二计算单元,可以用于利用多种地层孔隙压力计算方法对所述目标区参考井的声波时差测井数据和密度测井数据进行地层孔隙压力计算;
拟合处理单元,可以用于将所述地层孔隙压力计算的计算结果与所述参考井的地层孔隙压力进行比对拟合处理,确定出所述比对拟合处理的最大拟合度值所对应地层孔隙压力计算方法;
参数确定单元,可以用于对所述最大拟合度值所对应地层孔隙压力计算方法进行参数调节处理,确定出所述参数调节处理对应地层孔隙压力与所述参考井的地层孔隙压力进行比对拟合处理的最大拟合度值所对应参数;
压力计算公式设定单元,可以用于当所述最大拟合度值所对应地层孔隙压力计算方法的参数为所述所对应参数时,将所述地层孔隙压力计算方法所对应计算公式作为所述单井地层孔隙压力计算公式。
在一个优选的实施例中,所述埋藏深度校正处理可以包括采用如下计算公式:
vih=viL-(vif-vimin)
上式中,vih表示第二地震层速度体,所述第二地震层速度体表示所述第一地震层速度体进行埋藏深度校正处理后的地震层速度体;viL表示第一地震层速度体;vif表示埋藏深度速度数据体,所述埋藏深度速度数据体根据地震层速度体与埋深关系计算得到的与埋藏深度相对应的地震层速度体,vimin表示所述目的层顶界面的地震层速度。
在一个优选的实施例中,所述确定出的所述目标区的地震层速度与地层孔隙压力关系可以包括:
利用所述第二地震层速度体中的参考井的第二地震层速度和所述参考井的单井地层孔隙压力确定出的所述目标区地震层速度与地层孔隙压力关系。
由此可见,本申请一种确定油气分布数据的方法和装置的实施例基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;可以实现分岩性确定地震波阻抗反演体。然后,结合加德纳公式计算得到所述目标储集层的第一地震层速度体;利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;并由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体,保证利用单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力的准确性;最后根据所述空间地层孔隙压力中地层孔隙压力分布情况数据可以准确的确定出所述目标区的目标储集层的油气分布数据,实现通过圈定异常高压区进行油气预测。与现有技术相比,利用本申请实施例提供的技术方案可以更准确的确定出目标储集层的油气分布数据,从而可以为油气勘探预测工作提供数据支持。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (14)

1.一种确定油气分布数据的方法,其特征在于,所述方法包括:
基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;
利用所述目标储集层的地震波阻抗反演体结合加德纳公式计算得到所述目标储集层的第一地震层速度体;
利用所述目的层的地震层速度计算得到目的层的埋藏深度速度数据体;
利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;
基于确定出的所述目标区的地震层速度与地层孔隙压力关系,由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体;
确定单井地层孔隙压力曲线数据,利用所述单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力;
根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据。
2.根据权利要求1所述的方法,其特征在于,所述根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据包括:
判断所述空间地层孔隙压力中地层孔隙压力数值大于预设地层孔隙压力数值的百分比,将所述百分比作为所述目标区的目标储集层的含油气百分比;
根据所述含油气百分比和所述目的层的厚度确定所述油气分布范围和油气厚度,将所述油气分布范围和油气厚度作为油气分布数据。
3.根据权利要求2所述的方法,其特征在于,所述方法还包括:
根据所述目标区的单井地层孔隙压力确定含油气储集层的地层孔隙压力数值,将所述地层孔隙压力数值作为预设地层孔隙压力数值。
4.根据权利要求1至3任意一项所述的方法,其特征在于,所述确定单井地层孔隙压力曲线数据包括:
基于目标区参考井的地层孔隙压力进行多方法多参数对比拟合分析处理,确定出所述目标区的单井地层孔隙压力计算公式;
利用所述单井地层孔隙压力计算公式计算得到所述目标区的单井地层孔隙压力;
根据所述目标区的单井地层孔隙压力确定出单井地层孔隙压力曲线数据。
5.根据权利要求4所述的方法,其特征在于,所述基于目标区参考井的地层孔隙压力进行多方法多参数对比拟合分析处理,确定出所述目标区的单井地层孔隙压力计算公式包括:
利用多种地层孔隙压力计算方法对所述目标区参考井的声波时差测井数据和密度测井数据进行地层孔隙压力计算;
将所述地层孔隙压力计算的计算结果与所述参考井的地层孔隙压力进行比对拟合处理,确定出所述比对拟合处理的最大拟合度值所对应地层孔隙压力计算方法;
对所述最大拟合度值所对应地层孔隙压力计算方法进行参数调节处理,确定出所述参数调节处理对应地层孔隙压力与所述参考井的地层孔隙压力进行比对拟合处理的最大拟合度值所对应参数;
当所述最大拟合度值所对应地层孔隙压力计算方法的参数为所述所对应参数时,将所述地层孔隙压力计算方法所对应计算公式作为所述单井地层孔隙压力计算公式。
6.根据权利要求1至3任意一项所述的方法,其特征在于,所述埋藏深度校正处理包括采用如下计算公式:
vih=viL-(vif-vimin)
上式中,vih表示第二地震层速度体,所述第二地震层速度体表示所述第一地震层速度体进行埋藏深度校正处理后的地震层速度体;viL表示第一地震层速度体;vif表示埋藏深度速度数据体,所述埋藏深度速度数据体根据地震层速度体与埋深关系计算得到的与埋藏深度相对应的地震层速度体,vimin表示所述目的层顶界面的地震层速度。
7.根据权利要求1至3任意一项所述的方法,其特征在于,所述确定出的所述目标区的地震层速度与地层孔隙压力关系包括:
利用所述第二地震层速度体中的参考井的第二地震层速度和所述参考井的单井地层孔隙压力确定出的所述目标区地震层速度与地层孔隙压力关系。
8.一种确定油气分布数据的装置,其特征在于,所述装置包括:
地震波阻抗反演体确定模块,用于基于目标区目的层的岩性敏感参数从所述目标区的地震波阻抗反演体中确定出目标储集层的地震波阻抗反演体;
第一计算模块,用于利用所述目标储集层的地震波阻抗反演体结合加德纳公式计算得到所述目标储集层的第一地震层速度体;
第二计算模块,用于利用所述目的层的地震层速度计算得到目的层的埋藏深度速度数据体;
校正处理模块,用于利用所述埋藏深度速度数据体对所述第一地震层速度体进行埋藏深度校正处理,得到第二地震层速度体;
地层孔隙压力体确定模块,用于基于确定出的所述目标区的地震层速度与地层孔隙压力关系,由所述第二地震层速度体确定出所述目标区的目标储集层的地层孔隙压力体;
压力曲线数据确定模块,用于确定单井地层孔隙压力曲线数据;
空间地层孔隙压力确定模块,用于利用所述单井地层孔隙压力曲线数据对所述目标储集层的地层孔隙压力体进行校正处理,得到所述目标储集层的空间地层孔隙压力;
油气分布数据确定模块,用于根据所述空间地层孔隙压力中地层孔隙压力分布情况数据确定所述目标区的目标储集层的油气分布数据。
9.根据权利要求8所述的装置,其特征在于,所述油气分布数据确定模块包括:
第一数据处理单元,用于判断所述空间地层孔隙压力中地层孔隙压力数值大于预设地层孔隙压力数值的百分比,将所述百分比作为所述目标区的目标储集层的含油气百分比;
第二数据处理单元,用于根据所述含油气百分比和所述目的层的厚度确定所述油气分布范围和油气厚度,将所述油气分布范围和油气厚度作为油气分布数据。
10.根据权利要求9所述的装置,其特征在于,所述装置还包括:
预设地层孔隙压力数值确定模块,用于根据所述目标区的单井地层孔隙压力确定含油气储集层的地层孔隙压力数值,将所述地层孔隙压力数值作为预设地层孔隙压力数值。
11.根据权利要求8至10任意一项所述的装置,其特征在于,所述压力曲线数据确定模块包括:
压力计算公式确定单元,用于基于目标区参考井的地层孔隙压力进行多方法多参数对比拟合分析处理,确定出所述目标区的单井地层孔隙压力计算公式;
第一计算单元,用于利用所述单井地层孔隙压力计算公式计算得到所述目标区的单井地层孔隙压力;
压力曲线数据确定单元,用于根据所述目标区的单井地层孔隙压力确定出单井地层孔隙压力曲线数据。
12.根据权利要求11所述的装置,其特征在于,所述压力计算公式确定单元包括:
第二计算单元,用于利用多种地层孔隙压力计算方法对所述目标区参考井的声波时差测井数据和密度测井数据进行地层孔隙压力计算;
拟合处理单元,用于将所述地层孔隙压力计算的计算结果与所述参考井的地层孔隙压力进行比对拟合处理,确定出所述比对拟合处理的最大拟合度值所对应地层孔隙压力计算方法;
参数确定单元,用于对所述最大拟合度值所对应地层孔隙压力计算方法进行参数调节处理,确定出所述参数调节处理对应地层孔隙压力与所述参考井的地层孔隙压力进行比对拟合处理的最大拟合度值所对应参数;
压力计算公式设定单元,用于当所述最大拟合度值所对应地层孔隙压力计算方法的参数为所述所对应参数时,将所述地层孔隙压力计算方法所对应计算公式作为所述单井地层孔隙压力计算公式。
13.根据权利要求8至10任意一项所述的装置,其特征在于,所述埋藏深度校正处理包括采用如下计算公式:
vih=viL-(vif-vimin)
上式中,vih表示第二地震层速度体,所述第二地震层速度体表示所述第一地震层速度体进行埋藏深度校正处理后的地震层速度体;viL表示第一地震层速度体;vif表示埋藏深度速度数据体,所述埋藏深度速度数据体根据地震层速度体与埋深关系计算得到的与埋藏深度相对应的地震层速度体,vimin表示所述目的层顶界面的地震层速度。
14.根据权利要求8至10任意一项所述的装置,其特征在于,所述确定出的所述目标区的地震层速度与地层孔隙压力关系包括:
利用所述第二地震层速度体中的参考井的第二地震层速度和所述参考井的单井地层孔隙压力确定出的所述目标区地震层速度与地层孔隙压力关系。
CN201610159384.0A 2016-03-21 2016-03-21 一种确定油气分布数据的方法和装置 Active CN105700017B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610159384.0A CN105700017B (zh) 2016-03-21 2016-03-21 一种确定油气分布数据的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610159384.0A CN105700017B (zh) 2016-03-21 2016-03-21 一种确定油气分布数据的方法和装置

Publications (2)

Publication Number Publication Date
CN105700017A true CN105700017A (zh) 2016-06-22
CN105700017B CN105700017B (zh) 2018-06-01

Family

ID=56232150

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610159384.0A Active CN105700017B (zh) 2016-03-21 2016-03-21 一种确定油气分布数据的方法和装置

Country Status (1)

Country Link
CN (1) CN105700017B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109509111A (zh) * 2017-09-15 2019-03-22 中国石油化工股份有限公司 探井地层压力的预测方法及系统
CN109736784A (zh) * 2018-04-27 2019-05-10 长江大学 沉积岩地层孔隙压力预测计算方法
WO2020072367A1 (en) * 2018-10-04 2020-04-09 Schlumberger Technology Corporation Seismic pore-pressure prediction using prestack seismic inversion
CN111060986A (zh) * 2019-10-18 2020-04-24 中国石油化工股份有限公司 一种地层压力预测方法及岩性油藏评价方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120143510A1 (en) * 2007-05-25 2012-06-07 Aftab Alam High resolution attributes for seismic data processing and interpretation
CN103777245A (zh) * 2012-10-17 2014-05-07 中国石油化工股份有限公司 基于地震资料的油气成藏条件定量评价方法
CN104483703A (zh) * 2014-12-31 2015-04-01 中国石油天然气股份有限公司 一种井震联合古地层压力预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120143510A1 (en) * 2007-05-25 2012-06-07 Aftab Alam High resolution attributes for seismic data processing and interpretation
CN103777245A (zh) * 2012-10-17 2014-05-07 中国石油化工股份有限公司 基于地震资料的油气成藏条件定量评价方法
CN104483703A (zh) * 2014-12-31 2015-04-01 中国石油天然气股份有限公司 一种井震联合古地层压力预测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘强国 等: "《录井方法与原理》", 30 November 2011, 石油工业出版社 *
徐宝荣 等: "异常地层压力预测技术在准噶尔盆地的应用", 《新疆石油地质》 *
牛丽娟: "压力敏感性对低渗透油藏弹性产能影响", 《科学技术与工程》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109509111A (zh) * 2017-09-15 2019-03-22 中国石油化工股份有限公司 探井地层压力的预测方法及系统
CN109509111B (zh) * 2017-09-15 2021-09-28 中国石油化工股份有限公司 探井地层压力的预测方法及系统
CN109736784A (zh) * 2018-04-27 2019-05-10 长江大学 沉积岩地层孔隙压力预测计算方法
CN109736784B (zh) * 2018-04-27 2019-11-05 长江大学 沉积岩地层孔隙压力预测计算方法
WO2020072367A1 (en) * 2018-10-04 2020-04-09 Schlumberger Technology Corporation Seismic pore-pressure prediction using prestack seismic inversion
CN111060986A (zh) * 2019-10-18 2020-04-24 中国石油化工股份有限公司 一种地层压力预测方法及岩性油藏评价方法

Also Published As

Publication number Publication date
CN105700017B (zh) 2018-06-01

Similar Documents

Publication Publication Date Title
CN106951660B (zh) 一种海相碎屑岩水平井储层测井解释方法及装置
CN104047597B (zh) 含油气泥页岩地层测井曲线标准化方法
Vernik et al. Velocity anisotropy in shales: A petrophysical study
CN104267429A (zh) 确定地层压力的方法及装置
CN105700017A (zh) 一种确定油气分布数据的方法和装置
EA020635B1 (ru) Способ сейсмической разведки для поиска углеводородов с применением поля средних скоростей сейсмической волны
CN105549079A (zh) 一种地球物理参数的全波形反演模型的建立方法和装置
CN105068117A (zh) 用于裂缝性介质的avo反演方法、装置和设备
CN105093311A (zh) 叠合盆地多期差异隆升区地层剥蚀厚度测量方法
CN105626056A (zh) 一种确定目标区三维地层孔隙压力的方法和装置
CN105242318A (zh) 一种确定砂体连通关系的方法及装置
CN108873065B (zh) 砂岩优质储层预测方法及装置
Zhang et al. Simultaneous inversion for microseismic event location and velocity model in Vaca Muerta Formation
Gillespie et al. Hierarchical kink band development in the Appalachian Plateau decollement sheet
CN108629459B (zh) 储层含烃孔隙的检测方法及装置
Martorana et al. Seismo-stratigraphic model of “La Bandita” area in the Palermo Plain (Sicily, Italy) through HVSR inversion constrained by stratigraphic data
CN107728205A (zh) 一种地层压力预测方法
Rodríguez-Pradilla Microseismic monitoring of a hydraulic-fracturing operation in a CBM reservoir: Case study in the Cerrejón Formation, Cesar-Ranchería Basin, Colombia
Marín-Moreno et al. Inverse modelling and seismic data constraints on overpressure generation by disequilibrium compaction and aquathermal pressuring: application to the Eastern Black Sea Basin
Pochat et al. A simple method of determining sand/shale ratios from seismic analysis of growth faults: an example from upper Oligocene to lower Miocene Niger Delta deposits
CN106869915B (zh) 一种水平井井间隔夹层预测方法及装置
CN115857047B (zh) 一种地震储层综合预测方法
Huang et al. Fast-forward modeling of compressional arrival slowness logs in high-angle and horizontal wells
CN109991663B (zh) 工区地震速度体校正方法和装置
CN111435174A (zh) 强反射地区地震资料振幅补偿方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant