CN116165709B - 一种基于缝洞体识别的油气勘探方法、装置及设备 - Google Patents

一种基于缝洞体识别的油气勘探方法、装置及设备 Download PDF

Info

Publication number
CN116165709B
CN116165709B CN202310077353.0A CN202310077353A CN116165709B CN 116165709 B CN116165709 B CN 116165709B CN 202310077353 A CN202310077353 A CN 202310077353A CN 116165709 B CN116165709 B CN 116165709B
Authority
CN
China
Prior art keywords
target
attribute
determining
value
boundary
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
CN202310077353.0A
Other languages
English (en)
Other versions
CN116165709A (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 University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202310077353.0A priority Critical patent/CN116165709B/zh
Publication of CN116165709A publication Critical patent/CN116165709A/zh
Application granted granted Critical
Publication of CN116165709B publication Critical patent/CN116165709B/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. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/307Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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. analysis, for interpretation, for correction
    • 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. analysis, for interpretation, for correction
    • 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
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase
    • G01V2210/632Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]

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是本说明书的一个实施例提供的电子设备的结构组成示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案,下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本说明书一部分实施例,而不是全部的实施例。基于本说明书中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本说明书保护的范围。
考虑到碳酸盐岩储层是重要的石油天然气储产层之一。碳酸盐岩储产量约占全球储产量的一半以上。许多大型油气田和高产油气井的油气产量,大都是来自碳酸盐岩储集层。目前,石油天然气勘探开发热点逐渐向纵深转移,超深海相碳酸盐岩油气资源的有效勘探开发是油气上游业务可持续发展的必然选择。研究表明,走滑断裂机械破碎叠加后期岩溶改造对于碳酸盐岩储层的形成至关重要。优质储层主要发育在断裂周围的缝洞带内,为洞穴、孔洞及裂缝构成的缝洞系统,地震剖面多呈现串珠、杂乱及空白反射,非均性质性极强,预测难度大。
进一步,考虑到目前生产上主要采用叠后平面属性识别各类储层,如均方根振幅、最大波谷值、最大能量等振幅类属性识别串珠状储层。此类方法缺点有两个,一是无法识别由洞-孔-缝构成的完整缝洞体或缝洞系统,二是无法从空间立体角度描述缝洞系统,制约了钻井轨道设计,无法保障优质储层钻遇率。
针对现有方法存在的上述问题,本说明书引入一种基于缝洞体识别的油气勘探方法,从以实现对缝洞体或缝洞系统的识别,从而可以进行钻井轨道设计,提高油气勘探效率,保证优质储层钻遇率。
基于上述思路,本说明书提出一种基于缝洞体识别的油气勘探方法,首先,确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体;其次,根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体;进一步,根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值;根据所述边界属性门槛值,确定目标数据体;最后,根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探。参阅图1所示,本说明书实施例提供了一种基于缝洞体识别的油气勘探方法。具体实施时,该方法可以包括以下内容。
S101:确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体。
S102:根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体。
S103:根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体。
S104:根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值。
S105:根据所述边界属性门槛值,确定目标数据体。
S106:根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物。
S107:根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探。
在一些实施例中,上述确定目标灰岩地层的顶界值和底界值之前,在具体实施时,还可以包括:
S1:获取目标地震工区中的地震数据体,其中,所述目标地震工区中包括目标灰岩地层;
S2:根据所述地震数据体,确定目标灰岩地层的顶界值和底界值。
在一些实施例中,上述地震数据体可以在钻井之前,通过采集信号并对采集的信号进行处理获得,上述地震数据体可以为三维叠前深度偏移纯波地震数据体。上述目标灰岩地层可以为下古生界奥陶-寒武系发育而来的,可以将奥陶纪形成的地层称为奥陶系,分为下、中、上三个统。上述目标灰岩地层容易受到溶蚀形成缝洞系统或缝洞体,上述缝洞系统或缝洞体由洞、孔、缝构成。
在一些实施例中,可以利用三维叠前深度偏移纯波地震数据体,解释目标灰岩地层的顶层位和底层位,分别得到顶界值和底界值,其中,地层解释精度可以为1线×1道。
在一些实施例中,可以将顶界值向上漂移10个采样点作为计算顶界,将底界值向下漂移10个采样点作为计算底界,将计算顶界和计算底界作为上述的目标灰岩地层的顶界值和底界值,以计算顶界和计算底界为时窗,计算均方根振幅数据体。
在一些实施例中,上述根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体,在具体实施时,可以包括:
S1:根据所述均方根振幅数据体,分别确定沿线梯度数据体和沿道梯度数据体;
S2:根据所述沿线梯度数据体和所述沿道梯度数据体,确定目标缝洞体的振幅梯度属性体。
在一些实施例中,上述根据所述均方根振幅数据体,分别确定沿线梯度数据体和沿道梯度数据体,在具体实施时,可以包括:
S1:根据所述均方根振幅数据体,分别确定沿线方向的梯度值和沿道方向的梯度值;
S2:根据所述沿线方向的梯度值,确定沿线梯度数据体;
S3:根据所述沿道方向的梯度值,确定沿道梯度数据体。
在一些实施例中,上述沿线方向的梯度值可以为沿LINE方向的梯度值,上述沿道方向的梯度值可以为沿TRACE方向的梯度值,其中,LINE方向可以理解为在三维空间中(水平的x、y方向,纵向的z方向)的沿x方向,TRACE方向就是与x垂直的y方向。上述梯度值可以通过求导的方式确定,如:沿着x和y两个方向分别求导。
在一些实施例中,可以利用半径为3的高斯核函数对沿线方向的梯度值进行平滑处理得到沿线梯度数据体,可以利用半径为3的高斯核函数对沿道方向的梯度值进行平滑处理得到沿道梯度数据体。
在一些实施例中,可以利用沿线梯度数据体和沿道梯度数据体,计算综合数据体,综合数据体可以通过以下算式确定:
其中,dxy为综合数据体,dx为沿线梯度数据体,dy为沿道梯度数据体。
在一些实施例中,可以通过沿线梯度数据体和沿道梯度数据体计算出的综合数据体,确定振幅梯度属性体。如:可以利用半径为5的高斯核函数对综合数据体(dxy)做平滑处理,得到最终的振幅梯度属性体。
在一些实施例中,上述根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值,在具体实施时,可以包括:
S1:获取已钻井气测参数信息,其中,所述气测参数信息包括气测曲线;
S2:根据所述气测曲线和所述振幅梯度属性体,确定振幅梯度属性体的属性值;
S3:根据所述属性值,确定目标缝洞体的边界属性门槛值。
在一些实施例中,上述根据所述属性值,确定目标缝洞体的边界属性门槛值,在具体实施时,可以包括:
S1:对所述振幅梯度属性体的属性值进行求平均处理,得到所述属性值的平均值;
S2:将所述平均值作为目标缝洞体的边界属性门槛值。
在一些实施例中,上述已钻井可以为一口或者可以为多口已钻井,相应地,获取已钻井气测参数信息可以为获取一口已钻井的气测参数信息或者可以为获取多口已钻井中各口已钻井的气测参数信息。相应地,上述对所述振幅梯度属性体的属性值进行求平均处理可以为一口井的振幅梯度属性体的属性值或者对多口井的振幅梯度属性体的属性值进行求平均处理,可以将得到的平均值作为目标缝洞体的边界属性门槛值。
在一些实施例中,上述气测曲线可以在钻井过程中实时生成,上述气测曲线可以反应缝洞体或缝洞系统中的油气存在情况,如果缝洞体或缝洞系统中存在油气,则气测曲线会随之升高。可以通过上述气测曲线在升高或降低,确定升高或降低处对应的振幅梯度属性体的属性值。
在一些实施例中,上述根据所述属性值,确定目标缝洞体的边界属性门槛值,在具体实施时,还可以包括:
S1:确定所述振幅梯度属性体的属性值中的异常属性值;
S2:将所述异常属性值进行求平均处理,得到所述异常属性值的平均值;
S3:将所述异常属性值的平均值作为目标缝洞体的边界属性门槛值。
在一些实施例中,上述异常属性可以为在钻井过程中,钻到缝洞体边缘时,使得钻速突然增加,钻时急剧降低的异常情况。
在一些实施例中,上述根据所述边界属性门槛值,确定目标数据体,在具体实施过程中,还可以先将边界属性门槛值与上述振幅梯度属性体的属性值进行比对处理,得到比对处理结果,将比对处理结果对小于边界属性门槛值的振幅梯度属性体的属性值赋值为0,将对大于等于边界属性门槛值的振幅梯度属性体的属性值不赋值,从而得到新的数据体,作为目标数据体。通过赋值可以实现局部数据的可视化。
在一些实施例中,可以将目标数据体导入至三维可视化软件中,利用目标数据体实现对缝洞体的识别,即实现缝洞体的轮廓和填充物的识别。
在一些实施例中,参阅图2所示,上述根据所述目标轮廓和目标填充物,进行油气勘探,在具体实施时,可以包括:
S201:根据所述目标轮廓和目标填充物,确定目标灰岩地层的目标钻井轨道;
S202:根据所述目标钻井轨道,形成钻井方案;
S203:根据所述钻井方案,对目标灰岩地层进行油气勘探。
在一些实施例中,在识别出目标缝洞体的目标轮廓和目标填充物后,可以进行钻井轨道设计,形成目标钻井轨道,再根据目标钻井轨道,确定钻井方案,再基于钻井方案对储层进行开发,可以提高油层钻遇率以及油气勘探效率。
下面结合一个具体实施例对上述方法进行说明,然而值得注意的是,该具体实施例仅是为了更好地说明本申请,并不构成对本申请的不当限定。
在具体实施前,首先,获取目标地震工区中的地震数据体,其中,所述目标地震工区中包括目标灰岩地层;其次,根据所述地震数据体,确定目标灰岩地层的顶界值和底界值。在具体实施时,首先,确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体;其次,根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;然后,根据所述均方根振幅数据体,分别确定沿线方向的梯度值和沿道方向的梯度值,根据所述沿线方向的梯度值,确定沿线梯度数据体;根据所述沿道方向的梯度值,确定沿道梯度数据体;进一步,再根据所述沿线梯度数据体和所述沿道梯度数据体,确定目标缝洞体的振幅梯度属性体。进一步,获取已钻井气测参数信息,其中,所述气测参数信息包括气测曲线;根据所述气测曲线和所述振幅梯度属性体,确定振幅梯度属性体的属性值;对所述振幅梯度属性体的属性值进行求平均处理,得到所述属性值的平均值;将所述平均值作为目标缝洞体的边界属性门槛值;最后,将振幅梯度属性体的属性值与所述边界属性门槛值进行比对处理,得到比对处理结果,当对比处理结果中振幅梯度属性体的属性值小于边界属性门槛值时,将振幅梯度属性体的属性值赋值为0,当对比处理结果中振幅梯度属性体的属性值大于等于边界属性门槛值时,不对振幅梯度属性体的属性值赋值;根据赋值的情况,获得新的数据体,作为目标数据体;最后,根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;根据所述目标轮廓和所述目标填充物,确定目标灰岩地层的目标钻井轨道;根据所述目标钻井轨道,形成钻井方案;根据所述钻井方案,对目标灰岩地层进行油气勘探。通过上述方法可以能够准确识别缝洞体的轮廓和内部填充物,从而为后续设计钻井方案奠定基础,最终能够有效提高油气勘探效率。
在一个具体的场景示例中,可以应用本说明书实施例提供的方法实现缝洞系统的识别。具体实施时,参阅图3所示,可以包括以下步骤。
步骤一:获取均方根振幅数据体
1.1、利用三维叠前深度偏移纯波地震数据体,完成奥陶系灰岩顶层位和底层位解释,解释精度1线×1道。
1.2、将步骤1.1中得到的灰岩顶层位向上漂移10个采样点作为计算顶界,灰岩底层位向下漂移10个采样点作为计算底界,以此为时窗计算均方根振幅数据体。
步骤二、获取奥陶系灰岩段振幅梯度属性
2.1、利用步骤1.2中得到的均方根振幅数据体计算沿LINE方向的梯度。
2.2、利用半径为3的高斯核函数对步骤2.1得到的LINE方向的梯度数据体做平滑,得到数据体dx。
2.3、利用步骤1.2中得到的均方根振幅数据体计算沿TRACE方向的梯度。
2.4、利用半径为3的高斯核函数对步骤2.3得到的TRACE方向的梯度数据体做平滑,得到数据体dy。
2.5、将步骤2.2和步骤2.4中计算得的dx和dy代入如下公式,计算得到dxy。
2.6、利用半径为5的高斯核函数对步骤2.3得到的数据体dxy做平滑,得到最终的振幅梯度属性体。
步骤三、精细井震标定,确定振幅梯度属性门槛值
3.1、统计区块已钻井在灰岩地层钻进中的气测曲线升高或钻时曲线降低处所对应的步骤2.4得到振幅梯度属性体属性值。
3.2、将步骤3.1得到的钻时异常处的振幅梯度属性值求取平均值,作为断裂破碎带缝洞系统边界属性门槛值。
3.3、将小于步骤3.2得到的缝洞系统边界门槛值的属性值赋值为0,大于等于步骤3.2得到的缝洞系统边界门槛值的属性值不变。
3.4、将步骤3.3得到的数据体导入三维可视化软件,即实现了碳酸盐岩缝洞系统的识别。
虽然本说明书提供了如下述实施例或附图4所示的方法操作步骤或装置结构,但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中,这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时,可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。
基于上述一种基于缝洞体识别的油气勘探方法,本说明书还提出一种基于缝洞体识别的油气勘探装置的实施例。如图4所示,所述基于装置具体可以包括以下模块:
第一确定模块401,用于确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体;
第二确定模块402,用于根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体;根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值;
识别模块403,用于根据所述边界属性门槛值,确定目标数据体;根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;
钻井模块404,用于根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探。
在一些实施例中,上述第一确定模块401之前具体还可以用于获取目标地震工区中的地震数据体,其中,所述目标地震工区中包括目标灰岩地层;根据所述地震数据体,确定目标灰岩地层的顶界值和底界值。
在一些实施例中,上述第二确定模块402具体可以用于根据所述均方根振幅数据体,分别确定沿线梯度数据体和沿道梯度数据体;根据所述沿线梯度数据体和所述沿道梯度数据体,确定目标缝洞体的振幅梯度属性体。
在一些实施例中,上述第二确定模块402具体还可以用于根据所述均方根振幅数据体,分别确定沿线方向的梯度值和沿道方向的梯度值;根据所述沿线方向的梯度值,确定沿线梯度数据体;根据所述沿道方向的梯度值,确定沿道梯度数据体。
在一些实施例中,上述第二确定模块402具体还可以用于获取已钻井气测参数信息,其中,所述气测参数信息包括气测曲线;根据所述气测曲线和所述振幅梯度属性体,确定振幅梯度属性体的属性值;根据所述属性值,确定目标缝洞体的边界属性门槛值。
在一些实施例中,上述第二确定模块402具体还可以用于对所述振幅梯度属性体的属性值进行求平均处理,得到所述属性值的平均值;将所述平均值作为目标缝洞体的边界属性门槛值。
需要说明的是,上述实施例阐明的单元、装置或模块等,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本说明书时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
由上可见,基于本说明书实施例提供的一种基于缝洞体识别的油气勘探装置,可以实现对缝洞体的识别。基于识别后的结果,可以进行钻井轨道设计,提高油层钻遇率和油气勘探效率。
本说明书实施例还提供一种基于缝洞体识别的油气勘探电子设备,包括处理器以及用于存储处理器可执行指令的存储器,所述处理器具体实施时可以根据指令执行以下步骤:确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体;根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体;根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值;根据所述边界属性门槛值,确定目标数据体;根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探。
为了能够更加准确地完成上述指令,参阅图5所示,本说明书实施例还提供了另一种具体的电子设备,其中,所述电子设备包括网络通信端口501、处理器502以及存储器503,上述结构通过内部线缆相连,以便各个结构可以进行具体的数据交互。
其中,所述网络通信端口501,具体可以用于确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体。
所述处理器502,具体可以用于根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体;根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值;根据所述边界属性门槛值,确定目标数据体;根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探。
所述存储器503,具体可以用于存储相应的指令程序。
在本实施例中,所述网络通信端口501可以是与不同的通信协议进行绑定,从而可以发送或接收不同数据的虚拟端口。例如,所述网络通信端口可以是负责进行web数据通信的端口,也可以是负责进行FTP数据通信的端口,还可以是负责进行邮件数据通信的端口。此外,所述网络通信端口还可以是实体的通信接口或者通信芯片。例如,其可以为无线移动网络通信芯片,如GSM、CDMA等;其还可以为Wifi芯片;其还可以为蓝牙芯片。
在本实施例中,所述处理器502可以按任何适当的方式实现。例如,处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application SpecificIntegrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。本说明书并不作限定。
在本实施例中,所述存储器503可以包括多个层次,在数字系统中,只要能保存二进制数据的都可以是存储器;在集成电路中,一个没有实物形式的具有存储功能的电路也叫存储器,如RAM、FIFO等;在系统中,具有实物形式的存储设备也叫存储器,如内存条、TF卡等。
本说明书实施例还提供了一种基于上述缝洞体识别的油气勘探方法的计算机存储介质,所述计算机存储介质存储有计算机程序指令,在所述计算机程序指令被执行时实现:确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体;根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体;根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值;根据所述边界属性门槛值,确定目标数据体;根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探。
在本实施例中,上述存储介质包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的,用于进行网络连接通信的接口。
虽然本说明书提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。第一,第二等词语用来表示名称,而并不表示任何特定的顺序。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本说明书,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施例的描述可知,本领域的技术人员可以清楚地了解到本说明书可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本说明书的技术方案本质上可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设备等)执行本说明书各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本说明书可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本说明书,本领域普通技术人员知道,本说明书有许多变形而不脱离本说明书的精神,希望所附的权利要求包括这些变形而不脱离本说明书的精神。

Claims (9)

1.一种基于缝洞体识别的油气勘探方法,其特征在于,包括:
确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体;
根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;
根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体;
根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值;
根据所述边界属性门槛值,确定目标数据体;
根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;
根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探;
其中,所述确定目标灰岩地层的顶界值和底界值,包括:
利用地震数据体,解释目标灰岩地层的顶层位和底层位,分别得到顶界值和底界值;
将顶界值向上漂移10个采样点作为计算顶界,将底界值向下漂移10个采样点作为计算底界,将计算顶界和计算底界作为目标灰岩地层的顶界值和底界值;
其中,所述根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值,包括:
获取已钻井的气测参数信息,其中,所述气测参数信息包括气测曲线;
根据所述气测曲线和所述振幅梯度属性体,确定振幅梯度属性体的属性值;
根据所述属性值,确定目标缝洞体的边界属性门槛值;
其中,所述根据所述边界属性门槛值,确定目标数据体,包括:
将边界属性门槛值与所述振幅梯度属性体的属性值进行比对处理,得到比对处理结果,将大于等于边界属性门槛值的振幅梯度属性体的属性值不赋值,从而得到新的数据体,作为目标数据体。
2.根据权利要求1所述的方法,其特征在于,所述确定目标灰岩地层的顶界值和底界值之前,还包括:
获取目标地震工区中的地震数据体,其中,所述目标地震工区中包括目标灰岩地层。
3.根据权利要求1所述的方法,其特征在于,所述根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体,包括:
根据所述均方根振幅数据体,分别确定沿线梯度数据体和沿道梯度数据体;
根据所述沿线梯度数据体和所述沿道梯度数据体,确定目标缝洞体的振幅梯度属性体。
4.根据权利要求3所述的方法,其特征在于,所述根据所述均方根振幅数据体,分别确定沿线梯度数据体和沿道梯度数据体,包括:
根据所述均方根振幅数据体,分别确定沿线方向的梯度值和沿道方向的梯度值;
根据所述沿线方向的梯度值,确定沿线梯度数据体;
根据所述沿道方向的梯度值,确定沿道梯度数据体。
5.根据权利要求1所述的方法,其特征在于,所述根据所述属性值,确定目标缝洞体的边界属性门槛值,包括:
对所述振幅梯度属性体的属性值进行求平均处理,得到所述属性值的平均值;
将所述平均值作为目标缝洞体的边界属性门槛值。
6.根据权利要求1所述的方法,其特征在于,所述根据所述目标轮廓和目标填充物,对目标灰岩地层进行油气勘探,包括:
根据所述目标轮廓和目标填充物,确定目标灰岩地层的目标钻井轨道;
根据所述目标钻井轨道,形成钻井方案;
根据所述钻井方案,对目标灰岩地层进行油气勘探。
7.一种基于缝洞体识别的油气勘探装置,其特征在于,包括:
第一确定模块,用于确定目标灰岩地层的顶界值和底界值,其中,所述目标灰岩地层包括目标缝洞体;
第二确定模块,根据所述顶界值和底界值,确定目标缝洞体的均方根振幅数据体;根据所述均方根振幅数据体,确定目标缝洞体的振幅梯度属性体;根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值;
识别模块,用于根据所述边界属性门槛值,确定目标数据体;根据所述目标数据体对所述目标缝洞体进行识别,得到目标缝洞体的目标轮廓和目标填充物;
钻井模块,用于根据所述目标轮廓和所述目标填充物,对目标灰岩地层进行油气勘探;
其中,所述确定目标灰岩地层的顶界值和底界值,包括:
利用地震数据体,解释目标灰岩地层的顶层位和底层位,分别得到顶界值和底界值;
将顶界值向上漂移10个采样点作为计算顶界,将底界值向下漂移10个采样点作为计算底界,将计算顶界和计算底界作为目标灰岩地层的顶界值和底界值;
其中,所述根据所述振幅梯度属性体,确定目标缝洞体的边界属性门槛值,包括:
获取已钻井的气测参数信息,其中,所述气测参数信息包括气测曲线;
根据所述气测曲线和所述振幅梯度属性体,确定振幅梯度属性体的属性值;
根据所述属性值,确定目标缝洞体的边界属性门槛值;
其中,所述根据所述边界属性门槛值,确定目标数据体,包括:
将边界属性门槛值与所述振幅梯度属性体的属性值进行比对处理,得到比对处理结果,将大于等于边界属性门槛值的振幅梯度属性体的属性值不赋值,从而得到新的数据体,作为目标数据体。
8.一种基于缝洞体识别的油气勘探设备,其特征在于,包括:包括处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现权利要求1-6任一项所述方法的步骤。
9.一种计算机可读存储介质,其特征在于,其上存储有计算机指令,所述指令被处理器执行时实现权利要求1至6中任一项所述方法的步骤。
CN202310077353.0A 2023-01-19 2023-01-19 一种基于缝洞体识别的油气勘探方法、装置及设备 Active CN116165709B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310077353.0A CN116165709B (zh) 2023-01-19 2023-01-19 一种基于缝洞体识别的油气勘探方法、装置及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310077353.0A CN116165709B (zh) 2023-01-19 2023-01-19 一种基于缝洞体识别的油气勘探方法、装置及设备

Publications (2)

Publication Number Publication Date
CN116165709A CN116165709A (zh) 2023-05-26
CN116165709B true CN116165709B (zh) 2024-01-30

Family

ID=86414245

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310077353.0A Active CN116165709B (zh) 2023-01-19 2023-01-19 一种基于缝洞体识别的油气勘探方法、装置及设备

Country Status (1)

Country Link
CN (1) CN116165709B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158795A (zh) * 2015-08-27 2015-12-16 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用地层叠前纹理属性值来检测缝洞的方法
CN111290021A (zh) * 2020-03-25 2020-06-16 北京奥能恒业能源技术有限公司 基于梯度结构张量谱分解的碳酸盐岩缝洞增强识别方法
CN111413731A (zh) * 2019-01-08 2020-07-14 中国石油天然气集团有限公司 一种碳酸盐岩缝洞体的地震识别方法及装置
CN112698392A (zh) * 2020-11-20 2021-04-23 中国石油天然气股份有限公司 一种多层系储集体立体刻画与油气空间分布及定容方法
CN112698398A (zh) * 2020-11-20 2021-04-23 中国石油天然气股份有限公司 一种深大断裂体系空间刻画方法
CN114721046A (zh) * 2021-01-05 2022-07-08 中国石油天然气股份有限公司 缝洞检测方法、装置、计算机存储介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11143776B2 (en) * 2020-02-17 2021-10-12 China Petroleum & Chemical Corporation Computer-implemented method and system for small cave recognition using seismic reflection data

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158795A (zh) * 2015-08-27 2015-12-16 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用地层叠前纹理属性值来检测缝洞的方法
CN111413731A (zh) * 2019-01-08 2020-07-14 中国石油天然气集团有限公司 一种碳酸盐岩缝洞体的地震识别方法及装置
CN111290021A (zh) * 2020-03-25 2020-06-16 北京奥能恒业能源技术有限公司 基于梯度结构张量谱分解的碳酸盐岩缝洞增强识别方法
CN112698392A (zh) * 2020-11-20 2021-04-23 中国石油天然气股份有限公司 一种多层系储集体立体刻画与油气空间分布及定容方法
CN112698398A (zh) * 2020-11-20 2021-04-23 中国石油天然气股份有限公司 一种深大断裂体系空间刻画方法
CN114721046A (zh) * 2021-01-05 2022-07-08 中国石油天然气股份有限公司 缝洞检测方法、装置、计算机存储介质

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Seismic fracture characterization in tight sand reservoirs: A case study of the Xujiahe Formation, Sichuan Bbasin, China;Guangzhi Zhang等;Journal of Applied Geophysics;全文 *
四川盆地高峰场地区长兴组生物礁气藏类型研究;车琴;天然气勘探与开发;第15-16页 *
地震属性在塔河油田碳酸盐岩缝洞型油藏连通性分析及其注水开发中的应用;李小波;李新华;荣元帅;彭小龙;王可可;;油气地质与采收率(第06期);全文 *
塔里木油田奥陶系碳酸盐岩缝洞储层的定量地震描述;郑多明;李志华;赵宽志;沈凤;党青宁;高莲花;;中国石油勘探(第Z1期);全文 *

Also Published As

Publication number Publication date
CN116165709A (zh) 2023-05-26

Similar Documents

Publication Publication Date Title
CN104267431B (zh) 一种地震构造解释的小断层获取方法
US9026420B2 (en) Method and system of updating a geocellular model
CN104329079A (zh) 一种识别气测录井油气层的方法以及系统
CN111722277B (zh) 用于断溶体油气藏勘探的断溶体圈闭描述方法及装置
CN103941286A (zh) 一种基于平面导航的快速三维层位解释方法
CN104316958A (zh) 一种识别不同尺度地层断裂的相干处理方法
CN112147684A (zh) 表征同沉积断层活动强度的方法和装置
CN109598068B (zh) 古构造约束建模方法、装置和设备
CN109100803B (zh) 微断裂的确定方法和装置
CN106979002B (zh) 分支水平井预测剖面的确定方法和装置
CN116165709B (zh) 一种基于缝洞体识别的油气勘探方法、装置及设备
CN109143375A (zh) 一种确定钻井方位的方法及装置
CN108279435B (zh) 断层剖面的确定方法和装置
CN107368688B (zh) 海相单砂体的确定方法和装置
CN112698382A (zh) 小尺度断控岩溶型储层地震预测方法及装置
CN108646305B (zh) 一种界面的倾角获取方法、装置、电子设备及存储介质
CN109387873A (zh) 一种缝洞储集体反演方法及系统
Zhang et al. Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China
CN114429478A (zh) 断控岩溶体刻画方法、装置、存储介质和计算机设备
CN108343424B (zh) 钻井位置的确定方法和装置
CN108132485A (zh) 一种数据处理方法和装置
CN112394416A (zh) 非均质断裂控制储层预测方法及装置
CN113031056B (zh) 构造约束下断层封闭性分析方法及装置
CN108510586B (zh) 一种构建地质导向地层模型的方法
CN112780253B (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