CN114428356A - 相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备 - Google Patents
相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备 Download PDFInfo
- Publication number
- CN114428356A CN114428356A CN202011012403.XA CN202011012403A CN114428356A CN 114428356 A CN114428356 A CN 114428356A CN 202011012403 A CN202011012403 A CN 202011012403A CN 114428356 A CN114428356 A CN 114428356A
- Authority
- CN
- China
- Prior art keywords
- phase
- seismic
- longitudinal wave
- solution
- wave impedance
- 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.)
- Pending
Links
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
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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
- G01V2210/624—Reservoir parameters
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备,该方法包括:基于原始地震数据,分别获取断溶体地震相和裂缝地震相;建立碳酸盐岩直板低频模型;基于直板低频模型,获取串珠地震相;基于直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型;基于相控低频模型,获得纵波阻抗数据体;基于纵波阻抗数据体,获取有效储集体的孔隙度体,从而刻画碳酸盐岩断溶体内部结构。本发明对断溶体内部储层展布、发育规模和连通性直观刻画,并对断溶体内部各个储集体的规模进行准确预测,对井位部署提供可靠依据,提高勘探开发钻井的成功率,提升勘探开发效益,且该发明不受地域限制,适用范围广泛。
Description
技术领域
本发明属于油气勘探技术领域,具体涉及一种相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备。
背景技术
2015年鲁新便等提出断溶体圈闭的概念,主要指岩溶水沿断裂方向向下渗入或上涌对破碎带溶蚀改造,形成不同三维空间的岩溶缝洞系统,在上覆泥灰岩或侧向致密灰岩的遮挡下形成的特殊圈闭类型。塔里木盆地塔河外围斜坡区受上奥陶统地层覆盖的影响,岩溶缝洞体的发育与断裂带、构造变形之间有较好的匹配关系,具有明显的断控岩溶特征,以深大断裂带为核心发生溶蚀扩大,使之成为覆盖区中—下奥陶统岩溶缝洞体发育的有利区带。在断溶体油藏模式的指导下,随着勘探开发的深入,发现了顺北大型断溶体类型油气田。目前,顺北油气田已发现落实了与断溶体发育有关的18条走滑断裂带,资源规模达17×108t油当量,展示了顺北地区良好的勘探开发前景。
前人研究认为,断溶体内部主要储集空间为与走滑断裂相关的洞穴、构造高角度缝和沿缝溶蚀的空洞。构造破裂作用是断溶体储层发育的主控因素,溶蚀改造作用进一步促进储集体发育(焦方正2018)。“串珠”强能量地震相与洞穴储层相关性强,表现为低纵波阻抗;溶蚀孔洞表现为较强能量或空白反射地震相,对应较低纵波阻抗;裂缝型储层发育段与AFE、蚂蚁体等地震属性具有较好的相关性。
断溶体是缝洞结合的地质体,刻画内容包括断溶体边界空间展布的外部形态,以及断溶体内洞穴、溶蚀孔洞和连接孔洞的裂缝的内部结构。断溶体内部储层发育受多种因素控制,地震波场复杂,准确定量刻画断溶体内部储层难度大。邓光校等(2017)提出的发明专利《碳酸盐岩断溶体的刻画方法》,主要是利用叠后常规反演得到的地下纵波阻抗数据为背景,把代表裂缝发育的蚂蚁体属性进行镂空雕刻处理,然后把镂空处理后的蚂蚁体属性叠置在地下纵波阻抗之上,实现储集体与断裂的刻画。该方法用纵波阻抗属性刻画断溶体的储层,蚂蚁体属性刻画裂缝,只是利用地震属性刻画了断溶体内部发育孔洞和裂缝储层的内部结构,不能刻画断溶体边界的外部形态。
文欢等(2017)提出的发明专利《碳酸盐岩断溶体内部结构的表征方法》,主要是利用叠后常规反演得到地下纵波阻抗数据,确定代表缝洞体储层的纵波阻抗数据阈值,对地下纵波阻抗数据进行镂空雕刻处理;计算地震资料的张量属性,并对张量属性进行空间平滑处理,确定代表断溶体破碎区域的张量属性阈值,得到断溶体的轮廓;以所述断溶体的轮廓作为边界,将所述张量属性及所述地下纵波阻抗数据显示在轮廓内,实现断溶体的外部轮廓刻画及内部结构表征。该方法用纵波阻抗属性刻画断溶体内部孔洞型储层,用张量属性刻画断溶体边界的外部形态,不能刻画断溶体内部裂缝的空间展布。
通过上述技术背景调研表明,前人刻画断溶体的方法主要把反映断溶体外部形态和内部结构的张量、纵波阻抗和蚂蚁体三种地震属性中的其中两种属性叠置在一起,不能有效刻画断溶体的外部形态和内部结构。另外,假如把这三种地震属性叠置在一起刻画断溶体,也只能定性反映断溶体的外部形态和内部结构,不能进行定量描述断溶体的空间展布和有效储集空间的规模。在断溶体内部结构中,裂缝不仅能沟通孔洞,而且也是重要的储集空间。张量、纵波阻抗和蚂蚁体三种地震属性中,只有纵波阻抗地震属性与孔隙度有直接的关系,常规的储层反演得到的纵波阻抗,不能有效反映裂缝的空间展布。因此,特别需要一种方法能定量刻画断溶体内部缝洞结构规模及其展布。
发明内容
本发明的目的是提出一种相控反演刻画碳酸盐岩断溶体内部结构的方法,解决目前不能定量刻画断溶体内部缝洞结构规模及其展布的问题。
有鉴于此,本发明提供了一种相控反演刻画碳酸盐岩断溶体内部结构的方法、装置、电子设备及介质,至少解决目前没有形成统一的可操作性强的精度高的饱和度评价方法的问题。
第一方面,本发明提供一种相控反演刻画碳酸盐岩断溶体内部结构的方法,包括:基于原始地震数据,分别获取断溶体地震相和裂缝地震相;建立碳酸盐岩直板低频模型;基于所述直板低频模型,获取串珠地震相;基于所述直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型;基于所述相控低频模型,获得纵波阻抗数据体;基于所述纵波阻抗数据体,获取有效储集体的孔隙度体,从而刻画碳酸盐岩断溶体内部结构。
可选的,根据以下步骤获取断溶体地震相:对所述原始地震数据进行构造导向滤波处理;对构造导向滤波处理后的地震数据进行边缘检测处理;对边缘检测处理后的地震数据进行体元密度增强处理,获得体元密度数据及体元密度增强处理后的地震数据;将体元密度数据大于断溶体边界阈值的体元密度增强处理后的地震数据作为断溶体地震相。
可选的,根据以下步骤获取裂缝地震相:对所述原始地震数据进行构造导向滤波处理;对构造导向滤波处理后的原始地震数据进行AFE计算,获得AFE属性数据体;将AFE属性值小于裂缝发育的AFE属性阈值的AFE属性数据体作为裂缝地震相。
可选的,设定碳酸盐岩地层的纵波阻抗背景值,基于所述碳酸盐岩地层的纵波阻抗背景值建立碳酸盐岩直板低频模型。
可选的,所述基于所述直板低频模型,获取串珠地震相包括:基于所述直板低频模型,对原始地震数据进行约束稀疏脉冲反演,获得地下纵波阻抗数据;将纵波阻抗值小于洞穴储层的纵波阻抗阈值的地下纵波阻抗数据作为串珠地震相。
可选的,所述基于所述直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型包括:将所述断溶体地震相嵌入所述直板低频模型中;向所述直板低频模型中的断溶体内输入小于所述碳酸盐岩地层的纵波阻抗背景值,将所述串珠地震相嵌入所述断溶体地震相中;将所述裂缝地震相嵌入所述断溶体地震相中,获得相控低频模型。
可选的,所述将所述裂缝地震相嵌入所述断溶体地震相中包括:利用钻探结果,在嵌入串珠地震相后的断溶体地震相中输入裂缝地震相对应的纵波阻抗值,将所述裂缝地震相嵌入断溶体地震相中。
可选的,所述基于所述相控低频模型,获得纵波阻抗数据体包括:对所述相控低频模型进行反演,获得相控低频模型的地下纵波阻抗数据体;在所述相控低频模型的地下纵波阻抗数据体中,将高于碳酸盐岩地层的纵波阻抗背景值的纵波阻抗背景值用所述碳酸盐岩地层的纵波阻抗背景值替换;在断溶体内,将高于断溶体纵波阻抗背景值的背景值用所述断溶体纵波阻抗背景值替换,获得新的低频模型;对新的低频模型进行约束稀疏脉冲反演,获得纵波阻抗数据体。
可选的,基于所述纵波阻抗数据体,获取有效储集体的孔隙度体包括:将已知钻井的纵波阻抗数据与孔隙度进行拟合,获得纵波阻抗数据与孔隙度的关系;基于纵波阻抗数据与孔隙度的关系,将所述纵波阻抗数据体转换为孔隙度体:基于所述孔隙度体和有效储集体的孔隙度阈值,获取有效储集体的孔隙度体。
第二方面,本发明还提供一种电子设备,所述电子设备包括:存储器,存储有可执行指令;处理器,所述处理器运行所述存储器中的所述可执行指令,以实现上述相控反演刻画碳酸盐岩断溶体内部结构的方法。
本发明的有益效果在于:本发明的相控反演刻画碳酸盐岩断溶体内部结构的方法利用反映断溶体轮廓、断溶体内“串珠”强能量和裂缝的地震相,建立相控低频模型,对相控低频模型进行叠后反演,得到纵波阻抗数据,进一步获得有效储集体的孔隙度体,对断溶体内部储层展布、发育规模和连通性直观刻画,并对断溶体内部各个储集体的规模进行准确预测,对井位部署提供可靠依据,提高勘探开发钻井的成功率,提升勘探开发效益,且该发明不受地域限制,适用范围广泛。
本发明具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的,或者将在并入本文中的附图和随后的具体实施例中进行详细陈述,这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法流程图。
图2示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井断溶体剖面图。
图3示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井AFE属性剖面图。
图4示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井地震与断溶体边界剖面图。
图5示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井叠后约束稀疏脉冲反演阻抗剖面图。
图6示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过过某井串珠地震相剖面图。
图7示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井断溶体相、裂缝相、串珠地震相剖面剖面。
图8示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井相控反演纵波阻抗剖面图。
图9示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法的某断溶体内部雕刻圈闭厚度分布图。
图10示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法的某断溶体内部雕刻SHB5CX、SHB5-3井所钻圈闭厚度分布图。
图11示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过SHB5CX、SHB5-3连井相控反演纵波阻抗剖面图。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
本发明提供一种相控反演刻画碳酸盐岩断溶体内部结构的方法,包括:基于原始地震数据,分别获取断溶体地震相和裂缝地震相;建立碳酸盐岩直板低频模型;基于直板低频模型,获取串珠地震相;基于直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型;基于相控低频模型,获得纵波阻抗数据体;基于纵波阻抗数据体,获取有效储集体的孔隙度体,从而刻画碳酸盐岩断溶体内部结构。
具体的,分别获取断溶体地震相、裂缝和串珠地震相,将断溶体地震相、裂缝地震相和串珠地震相嵌入到直板低频模型,获得相控低频模型,基于相控低频模型,获得纵波阻抗数据体,进一步获得有效储集体的孔隙度体,有效储集体的孔隙度体刻画了碳酸盐岩断溶体内部结构。
根据示例性的实施方式,相控反演刻画碳酸盐岩断溶体内部结构的方法利用反映断溶体轮廓、断溶体内“串珠”强能量和裂缝的地震相,建立相控低频模型,对相控低频模型进行叠后反演,得到纵波阻抗数据,进一步获得有效储集体的孔隙度体,对断溶体内部储层展布、发育规模和连通性直观刻画,并对断溶体内部各个储集体的规模进行准确预测,对井位部署提供可靠依据,提高勘探开发钻井的成功率,提升勘探开发效益,且该发明不受地域限制,适用范围广泛。
作为可选方案,根据以下步骤获取断溶体地震相:对原始地震数据进行构造导向滤波处理;对构造导向滤波处理后的地震数据进行边缘检测处理;对边缘检测处理后的地震数据进行体元密度增强处理,获得体元密度数据及体元密度增强处理后的地震数据;将体元密度数据大于断溶体边界阈值的体元密度增强处理后的地震数据作为断溶体地震相。
具体的,先对原始地震资料进行构造导向滤波处理,然后对构造导向滤波处理后的地震资料进行边缘检测处理,然后对边缘检测处理后的地震资料进行体元密度增强处理,体元密度数据及体元密度增强处理后的地震数据。
碳酸盐岩和碎屑岩的分界面,是个标志地震反射层。对该地震反射层进行大平滑处理,该层就接近水平了,平滑后的地震反射层与原始的地震反射层交叉部位对应的体元密度值数据体的值(体元密度数据体的值是从小到大连续分布的),作为断溶体边界的门槛值,大于这个门槛值的数据就是断溶体在空间上所包含的数据,因此把这个门槛值作为断溶体边界值,大于这个门槛值就是属于断溶体内部的数据,小于这个门槛值就不是断溶体的包含的数据。
作为可选方案,根据以下步骤获取裂缝地震相:对原始地震数据进行构造导向滤波处理;对构造导向滤波处理后的原始地震数据进行AFE计算,获得AFE属性数据体;将AFE属性值小于裂缝发育的AFE属性阈值的AFE属性数据体作为裂缝地震相。
对原始地震资料进行构造导向滤波处理,对构造导向滤波处理后的地震资料进行进行AFE计算,得到AFE属性数据体,AFE属性反映裂缝;将AFE属性值小于裂缝发育的AFE属性阈值的AFE属性数据体作为裂缝地震相。利用实际钻井结果,确定裂缝发育的AFE属性阈值,即AFE的门槛值。
钻探结果表明,钻遇裂缝时表现为钻井液漏失,在成像测井中也可以看到裂缝表现的地球物理特征,通过钻探证实裂缝在某一深度发育,通过精细井震标定,建立准确的时-深关系,确定实际钻到的裂缝在时间域地震剖面上的纵向位置,计算得到的反映裂缝的AFE属性数据是从原始的地震资料得到的,在纵向时间上没有改变,和原始的地震资料在纵向时间上保持一致。AFE属性数据值也是从小到大连续分布的,前面通过井震标定确定了实际钻到裂缝在时间域空间上的发育位置,在AFE数据体中对应的值就可以确定裂缝对应的AFE值,该值做为裂缝发育的AFE属性的门槛值(AFE属性值是从小到大分布的),大于这个门槛值就不发育裂缝,小于这个门槛值发育裂缝,因此小于门槛值的AFE属性体就是裂缝发育的裂缝相。
作为可选方案,设定碳酸盐岩地层的纵波阻抗背景值,基于碳酸盐岩地层的纵波阻抗背景值建立碳酸盐岩直板低频模型。
作为可选方案,基于直板低频模型,获取串珠地震相包括:基于直板低频模型,对原始地震数据进行约束稀疏脉冲反演,获得地下纵波阻抗数据;将纵波阻抗值小于洞穴储层的纵波阻抗阈值的地下纵波阻抗数据作为串珠地震相。
具体的,利用直板低频模型,对原始地震资料进行约束稀疏脉冲反演,得到地下纵波阻抗数据,给定洞穴储层的纵波阻抗门槛值,得到反映洞穴储层“串珠”地震相。
在碳酸盐岩地层中,孔隙度大于5%就是好储层,对应洞穴型储层,通过实际的钻井资料和测井资料,把孔隙度大于5%的层段找出来,通过井震标定,确定孔隙度大于5%的发育层段对应地震剖面的位置,该位置就对应洞穴型储层,表现为串珠状地震反射,孔隙度为5%对应的纵波阻抗值就是串珠地震相的门槛值,小于这个门槛值的都是洞穴层串珠地震相。
作为可选方案,基于直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型包括:将断溶体地震相嵌入直板低频模型中;向直板低频模型中的断溶体内输入小于碳酸盐岩地层的纵波阻抗背景值,将串珠地震相嵌入断溶体地震相中;将裂缝地震相嵌入断溶体地震相中,获得相控低频模型。
具体的,把断溶体地震相嵌入直板低频模型中,给断溶体内部输入低于碳酸盐岩地层的的纵波阻抗背景值,把刻画的“串珠”地震相嵌入断溶体中,利用钻探结果,输入裂缝相的纵波阻抗值,嵌入断溶体中建立相控低频模型。
作为可选方案,将裂缝地震相嵌入断溶体地震相中包括:利用钻探结果,在嵌入串珠地震相后的断溶体地震相中输入裂缝地震相对应的纵波阻抗值,将裂缝地震相嵌入断溶体地震相中。
裂缝是有孔隙度的,通过工区已钻井的的纵波阻抗和孔隙度的交汇统计分析,确定反映裂缝孔隙度的阻抗的最大门槛值。在这里裂缝的纵波阻抗门槛值是小于断溶体的门槛值的。裂缝地震相是空间分布范围的,把得到的裂缝地震相的空间分布都替换为裂缝的纵波阻抗门槛值,就得到反映裂缝相的纵波阻抗低频模型。
作为可选方案,基于相控低频模型,获得纵波阻抗数据体包括:对相控低频模型进行反演,获得相控低频模型的地下纵波阻抗数据体;在相控低频模型的地下纵波阻抗数据体中,将高于碳酸盐岩地层的纵波阻抗背景值的纵波阻抗背景值用碳酸盐岩地层的纵波阻抗背景值替换;在断溶体内,将高于断溶体纵波阻抗背景值的背景值用断溶体纵波阻抗背景值替换,获得新的低频模型;对新的低频模型进行约束稀疏脉冲反演,获得纵波阻抗数据体。
具体的,利用相控低频模型,开展约束稀疏脉冲反演,得到地下纵波阻抗数据;在地下纵波阻抗数据体中,把高于碳酸盐岩地层的纵波阻抗背景值的用于碳酸盐岩地层的纵波阻抗背景值替换,把断溶体内高于断溶体纵波阻抗背景值的用断溶体纵波阻抗背景值替换,就得到新的低频模型,用新的低频模型进行第二次约束稀疏脉冲反演,得到纵波阻抗体。
作为可选方案,基于纵波阻抗数据体,获取有效储集体的孔隙度体包括:将已知钻井的纵波阻抗数据与孔隙度进行拟合,获得纵波阻抗数据与孔隙度的关系;基于纵波阻抗数据与孔隙度的关系,将纵波阻抗数据体转换为孔隙度体:基于孔隙度体和有效储集体的孔隙度阈值,获取有效储集体的孔隙度体。
具体的,工区内碳酸盐岩有效储集体的孔隙度门槛值是已知的,高于这个孔隙度就是有效储集体。通过已钻井碳酸盐岩层段的有效孔隙度和纵波阻抗交汇得到有效孔隙度对应的纵波阻抗门槛值,小于这个门槛值的纵波阻抗就是有效储集体对应的纵波阻抗。
通过已钻井纵波阻抗与孔隙度拟合,得到孔隙度对应的纵波阻抗值,依据该孔隙度与纵波阻抗拟合曲线,把孔隙度和纵波阻抗联系起来,就可以把纵波阻抗体转换为孔隙度体。
基于有效储集体的孔隙度阈值,去除孔隙度体中非有效储集体的部分,获取有效储集体的孔隙度体。
本发明还提供一种电子设备,电子设备包括:存储器,存储有可执行指令;处理器,处理器运行存储器中的可执行指令,以实现上述相控反演刻画碳酸盐岩断溶体内部结构的方法。
实施例一
图1示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法的流程图。图2示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井断溶体剖面图。图3示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井AFE属性剖面图。图4示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井地震与断溶体边界剖面图。图5示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井叠后约束稀疏脉冲反演阻抗剖面图。图6示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过过某井串珠地震相剖面图。图7示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井断溶体相、裂缝相、串珠地震相剖面剖面。图8示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过某井相控反演纵波阻抗剖面图。图9示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法的某断溶体内部雕刻圈闭厚度分布图。图10示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法的某断溶体内部雕刻SHB5CX、SHB5-3井所钻圈闭厚度分布图。图11示出了根据本发明的一个实施例的一种相控反演刻画碳酸盐岩断溶体内部结构的方法过SHB5CX、SHB5-3连井相控反演纵波阻抗剖面图。
结合图1、图2、图3、图4、图5、图6、图7、图8、图9、图10和图11所示,该相控反演刻画碳酸盐岩断溶体内部结构的方法,包括:
步骤1:基于原始地震数据,分别获取断溶体地震相和裂缝地震相;
其中,根据以下步骤获取断溶体地震相:对原始地震数据进行构造导向滤波处理;对构造导向滤波处理后的地震数据进行边缘检测处理;对边缘检测处理后的地震数据进行体元密度增强处理,获得体元密度数据及体元密度增强处理后的地震数据;将体元密度数据大于断溶体边界阈值的体元密度增强处理后的地震数据作为断溶体地震相。
其中,根据以下步骤获取裂缝地震相:对原始地震数据进行构造导向滤波处理;对构造导向滤波处理后的原始地震数据进行AFE计算,获得AFE属性数据体;将AFE属性值小于裂缝发育的AFE属性阈值的AFE属性数据体作为裂缝地震相。
步骤2:建立碳酸盐岩直板低频模型;
其中,设定碳酸盐岩地层的纵波阻抗背景值,基于碳酸盐岩地层的纵波阻抗背景值建立碳酸盐岩直板低频模型。
步骤3:基于直板低频模型,获取串珠地震相;
其中,基于直板低频模型,获取串珠地震相包括:基于直板低频模型,对原始地震数据进行约束稀疏脉冲反演,获得地下纵波阻抗数据;将纵波阻抗值小于洞穴储层的纵波阻抗阈值的地下纵波阻抗数据作为串珠地震相。
步骤4:基于直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型;
其中,基于直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型包括:将断溶体地震相嵌入直板低频模型中;向直板低频模型中的断溶体内输入小于碳酸盐岩地层的纵波阻抗背景值,将串珠地震相嵌入断溶体地震相中;将裂缝地震相嵌入断溶体地震相中,获得相控低频模型。
其中,将裂缝地震相嵌入断溶体地震相中包括:利用钻探结果,在嵌入串珠地震相后的断溶体地震相中输入裂缝地震相对应的纵波阻抗值,将裂缝地震相嵌入断溶体地震相中。
步骤5:基于相控低频模型,获得纵波阻抗数据体;
其中,基于相控低频模型,获得纵波阻抗数据体包括:对相控低频模型进行反演,获得相控低频模型的地下纵波阻抗数据体;在相控低频模型的地下纵波阻抗数据体中,将高于碳酸盐岩地层的纵波阻抗背景值的纵波阻抗背景值用碳酸盐岩地层的纵波阻抗背景值替换;在断溶体内,将高于断溶体纵波阻抗背景值的背景值用断溶体纵波阻抗背景值替换,获得新的低频模型;对新的低频模型进行约束稀疏脉冲反演,获得纵波阻抗数据体。
步骤6:基于纵波阻抗数据体,获取有效储集体的孔隙度体,从而刻画碳酸盐岩断溶体内部结构。
其中,基于纵波阻抗数据体,获取有效储集体的孔隙度体包括:将已知钻井的纵波阻抗数据与孔隙度进行拟合,获得纵波阻抗数据与孔隙度的关系;基于纵波阻抗数据与孔隙度的关系,将纵波阻抗数据体转换为孔隙度体:基于孔隙度体和有效储集体的孔隙度阈值,获取有效储集体的孔隙度体。
在图2中,图中地震反射层与大平滑后的层交叉部位的体元密度属性值附近即为刻画断溶体边界的阈值。在图3中,为通过AFE计算反映裂缝的地震相。在图8中,直观反映了断溶体的形状和断溶体内部缝洞纵向、横向的连通性。颜色越亮代表阻抗值越低,反映储层物性越好。
将本发明应用于塔里木盆地目前较为关注的以断溶体为勘探对象的区块,某井侧钻至7945.09m,失返性漏失,放空进尺2.92m,漏失泥浆482.01m3。放空段对应图中的该方法预测的低纵波阻抗区,如见图8所示,实钻与预测结果吻合。对该方法刻画断溶体内部结构结果看出如图9的蓝色圆圈内,俯瞰SHB5CX和SHB5-3井所钻圈闭是多个圈闭连在一起,实际上SHB5CX和SHB5-3井所钻圈闭在空间上是分开的,如图10的蓝色圆圈内,用该方法预测这两口钻井的连井剖面的结果表明也是不连通的,如图11所示,这两口井后期的实际产能表明这两口井是不连通,预测结果与实际生产情况吻合。
实施例二
本公开提供一种电子设备包括,该电子设备包括:存储器,存储有可执行指令;处理器,处理器运行存储器中的可执行指令,以实现上述相控反演刻画碳酸盐岩断溶体内部结构的方法。
根据本公开实施例的电子设备包括存储器和处理器。
该存储器用于存储非暂时性计算机可读指令。具体地,存储器可以包括一个或多个计算机程序产品,该计算机程序产品可以包括各种形式的计算机可读存储介质,例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。
该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元,并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中,该处理器用于运行该存储器中存储的该计算机可读指令。
本领域技术人员应能理解,为了解决如何获取良好用户体验效果的技术问题,本实施例中也可以包括诸如通信总线、接口等公知的结构,这些公知的结构也应包含在本公开的保护范围之内。
有关本实施例的详细说明可以参考前述各实施例中的相应说明,在此不再赘述。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,包括:
基于原始地震数据,分别获取断溶体地震相和裂缝地震相;
建立碳酸盐岩直板低频模型;
基于所述直板低频模型,获取串珠地震相;
基于所述直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型;
基于所述相控低频模型,获得纵波阻抗数据体;
基于所述纵波阻抗数据体,获取有效储集体的孔隙度体,从而刻画碳酸盐岩断溶体内部结构。
2.根据权利要求1所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,根据以下步骤获取断溶体地震相:
对所述原始地震数据进行构造导向滤波处理;
对构造导向滤波处理后的地震数据进行边缘检测处理;
对边缘检测处理后的地震数据进行体元密度增强处理,获得体元密度数据及体元密度增强处理后的地震数据;
将体元密度数据大于断溶体边界阈值的体元密度增强处理后的地震数据作为断溶体地震相。
3.根据权利要求2所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,根据以下步骤获取裂缝地震相:
对所述原始地震数据进行构造导向滤波处理;
对构造导向滤波处理后的原始地震数据进行AFE计算,获得AFE属性数据体;
将AFE属性值小于裂缝发育的AFE属性阈值的AFE属性数据体作为裂缝地震相。
4.根据权利要求3所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,设定碳酸盐岩地层的纵波阻抗背景值,基于所述碳酸盐岩地层的纵波阻抗背景值建立碳酸盐岩直板低频模型。
5.根据权利要求4所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,所述基于所述直板低频模型,获取串珠地震相包括:
基于所述直板低频模型,对原始地震数据进行约束稀疏脉冲反演,获得地下纵波阻抗数据;
将纵波阻抗值小于洞穴储层的纵波阻抗阈值的地下纵波阻抗数据作为串珠地震相。
6.根据权利要求5所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,所述基于所述直板低频模型、断溶体地震相、裂缝地震相和串珠地震相,获得相控低频模型包括:
将所述断溶体地震相嵌入所述直板低频模型中;
向所述直板低频模型中的断溶体内输入小于所述碳酸盐岩地层的纵波阻抗背景值,将所述串珠地震相嵌入所述断溶体地震相中;
将所述裂缝地震相嵌入所述断溶体地震相中,获得相控低频模型。
7.根据权利要求6所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,所述将所述裂缝地震相嵌入所述断溶体地震相中包括:
利用钻探结果,在嵌入串珠地震相后的断溶体地震相中输入裂缝地震相对应的纵波阻抗值,将所述裂缝地震相嵌入断溶体地震相中。
8.根据权利要求7所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,所述基于所述相控低频模型,获得纵波阻抗数据体包括:
对所述相控低频模型进行反演,获得相控低频模型的地下纵波阻抗数据体;
在所述相控低频模型的地下纵波阻抗数据体中,将高于碳酸盐岩地层的纵波阻抗背景值的背景值用所述碳酸盐岩地层的纵波阻抗背景值替换;
在断溶体内,将高于断溶体纵波阻抗背景值的背景值用所述断溶体纵波阻抗背景值替换,获得新的低频模型;
对新的低频模型进行约束稀疏脉冲反演,获得纵波阻抗数据体。
9.根据权利要求8所述的相控反演刻画碳酸盐岩断溶体内部结构的方法,其特征在于,基于所述纵波阻抗数据体,获取有效储集体的孔隙度体包括:
将已知钻井的纵波阻抗数据与孔隙度进行拟合,获得纵波阻抗数据与孔隙度的关系;
基于纵波阻抗数据与孔隙度的关系,将所述纵波阻抗数据体转换为孔隙度体:
基于所述孔隙度体和有效储集体的孔隙度阈值,获取有效储集体的孔隙度体。
10.一种电子设备,其特征在于,所述电子设备包括:
存储器,存储有可执行指令;
处理器,所述处理器运行所述存储器中的所述可执行指令,以实现根据权利要求1-9中任一项所述的相控反演刻画碳酸盐岩断溶体内部结构的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011012403.XA CN114428356A (zh) | 2020-09-23 | 2020-09-23 | 相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011012403.XA CN114428356A (zh) | 2020-09-23 | 2020-09-23 | 相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114428356A true CN114428356A (zh) | 2022-05-03 |
Family
ID=81309643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011012403.XA Pending CN114428356A (zh) | 2020-09-23 | 2020-09-23 | 相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114428356A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115421214A (zh) * | 2022-10-10 | 2022-12-02 | 成都理工大学 | 基于断裂识别和易溶岩相组合的断溶体刻画方法 |
CN116299701A (zh) * | 2023-02-23 | 2023-06-23 | 北京阳光吉澳能源技术有限公司 | 断溶体储层识别方法、装置及电子设备 |
WO2024103475A1 (zh) * | 2022-11-14 | 2024-05-23 | 中国石油天然气股份有限公司 | 一种碳酸盐岩油气藏定量表征方法与系统 |
-
2020
- 2020-09-23 CN CN202011012403.XA patent/CN114428356A/zh active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115421214A (zh) * | 2022-10-10 | 2022-12-02 | 成都理工大学 | 基于断裂识别和易溶岩相组合的断溶体刻画方法 |
CN115421214B (zh) * | 2022-10-10 | 2023-10-20 | 成都理工大学 | 基于断裂识别和易溶岩相组合的断溶体刻画方法 |
WO2024103475A1 (zh) * | 2022-11-14 | 2024-05-23 | 中国石油天然气股份有限公司 | 一种碳酸盐岩油气藏定量表征方法与系统 |
CN116299701A (zh) * | 2023-02-23 | 2023-06-23 | 北京阳光吉澳能源技术有限公司 | 断溶体储层识别方法、装置及电子设备 |
CN116299701B (zh) * | 2023-02-23 | 2024-02-13 | 北京阳光吉澳能源技术有限公司 | 断溶体储层识别方法、装置及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114428356A (zh) | 相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备 | |
RU2669948C2 (ru) | Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности | |
CN109116428B (zh) | 一种缝洞型碳酸盐岩油藏不确定性建模方法及其装置 | |
CN110308487B (zh) | 一种断溶体型油藏定量表征方法 | |
CA2823989C (en) | Method and system of updating a geocellular model | |
CN111722277B (zh) | 用于断溶体油气藏勘探的断溶体圈闭描述方法及装置 | |
US11280179B2 (en) | Evaluation of pressure-transient behavior of wells | |
US10380793B2 (en) | Geobody surface reconstruction | |
US20220163692A1 (en) | Modeling and simulating faults in subterranean formations | |
US11294095B2 (en) | Reservoir simulations with fracture networks | |
CN115166853A (zh) | 页岩气藏天然裂缝模型建立方法、装置、电子设备及介质 | |
US9988881B2 (en) | Surface representation for modeling geological surfaces | |
EP3338115A1 (en) | Reservoir simulations with fracture networks | |
CN106869915B (zh) | 一种水平井井间隔夹层预测方法及装置 | |
CN113409463B (zh) | 一种包括尖灭处理的三维地质模型构建方法及装置 | |
CN116559938A (zh) | 一种油气藏裂缝模型建立方法及电子设备 | |
CN115407407A (zh) | 针对碳酸盐岩古溶洞及其充填的三维地质模型构建方法 | |
CN114545498A (zh) | 一种古岩溶缝洞结构的综合刻画方法 | |
CN112394416A (zh) | 非均质断裂控制储层预测方法及装置 | |
CN112782780A (zh) | 一种基于岩石物理相的储层评价方法、装置及设备 | |
CN112180466B (zh) | 一种河道砂体识别方法及装置 | |
Zhang et al. | Analysis of the Lost Circulation Problem. | |
CN107870366B (zh) | 开发动态硬约束的砂砾岩连通体地震预测方法 | |
CN116931075A (zh) | 一种基于剩余阻抗刻画断溶体储层的反演方法 | |
CN116125555A (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 |