CN103675904B - 一种井震匹配解释性目标处理方法及装置 - Google Patents
一种井震匹配解释性目标处理方法及装置 Download PDFInfo
- Publication number
- CN103675904B CN103675904B CN201210332849.XA CN201210332849A CN103675904B CN 103675904 B CN103675904 B CN 103675904B CN 201210332849 A CN201210332849 A CN 201210332849A CN 103675904 B CN103675904 B CN 103675904B
- Authority
- CN
- China
- Prior art keywords
- seismic
- well
- frequency range
- volume
- frequency domain
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种井震匹配解释性目标处理方法及装置,该方法包括:进行基础数据准备:首先解释目的层的顶底层位,以目的层为中心沿层选取时窗,分析该时窗内数据体的有效频带,在有效频带基础上选取常相位子波;然后利用常相位子波及优选的井制作合成地震记录;设计多个频段的三角滤波门;进行比例系数的求取及应用:在目标处理的时窗内,在井点处求取每个频段的合成地震记录与实际地震数据的比例系数;进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;并应用求取的比例系数得到空间统一比例系数井震匹配处理后的地震数据体;进行匹配效果监控,最终获得预期的井震匹配目标处理数据体。本发明可提高地震资料的分辨能力。
Description
技术领域
本发明涉及油藏地球物理技术领域,尤其涉及一种井震匹配解释性目标处理方法及装置。
背景技术
地震成果数据体在有效频带内,经常有优势频率成分非常突出的情况,这样其它频率成分的贡献相对减小,导致地震资料延续相位增多,分辨能力下降。这种现象源于地震勘探过程中的多个环节,首先地震成果数据体是地下共反射点反射波多次叠加的结果,由于一些原因,比如静校正误差、速度参数误差、各向异性参数误差等,都可以导致这一多次叠加不能全部同相叠加,这种不完全同相叠加的主要效果是对高频成分的强烈衰减,频率越高衰减越大;另一方面,地震反射波在整个地震勘探过程中,经历过多个带通滤波,比如检波器的频率特性、检波器组合接收、炮点组合激发、用优势频带进行速度分析和剩余静校正量求取,它们的实际效果都有一定程度的突出优势频率成分、压制其它成分的作用。
综上可见,如何恢复被压制频率成分在地震资料中的作用,以提高提高地震资料的分辨能力,从而为油藏地球物理及油藏地质研究提供高精度的地震资料,是目前亟待解决的一个技术问题。
发明内容
本发明实施例提供一种井震匹配解释性目标处理方法及装置,以提高提高地震资料的分辨能力,从而为油藏地球物理及油藏地质研究提供高精度的地震资料。
一方面,本发明实施例提供了一种井震匹配解释性目标处理方法,所述井震匹配解释性目标处理方法,包括:
a、进行基础数据准备:首先解释目的层的顶底层位,以目的层为中心,沿层选取时间窗口,简称时窗,作为目标处理求取比例系数的数据体时窗,分析该时窗内数据体的有效频带,在有效频带基础上选取常相位子波;然后利用常相位子波及优选的井制作合成地震记录;
b、进行滤波门设计:依据地震数据体目的层有效频带,设计多个频段的三角滤波门,其中,前一个滤波门的后斜坡与后一个滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择滤波门的个数;
c、进行比例系数的求取及应用:在目标处理的时窗内,在井点处求取每个频段的合成地震记录与实际地震数据的比例系数;将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;并应用求取的比例系数得到空间统一比例系数井震匹配处理后的地震数据体;
d、进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,继续上述的步骤,不断优化三角滤波门设计及分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
可选的,在本发明一实施例中,所述步骤a中,在多井情况下,进一步优选合成地震记录与井旁地震道相关性较好的井,制作多井合成地震记录,生成相应的合成地震记录数据体。
可选的,在本发明一实施例中,所述步骤c中,在时间域进行比例系数求取及应用,包括:
(3.1)将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;
(3.2)在目标处理的时窗内,在井点处求取每个频段合成地震记录子体和相应地震子体的地震属性,并计算其地震属性的比值,然后在多个地震属性比值的基础上,优选多井不同频段比例系数组合,同时用多井平均求取不同频段统一比例系数组合;其中,所述地震属性包括:每个频段合成地震记录子体和相应地震子体的均方根振幅、总振幅、总能量;
(3.3)利用上述优选的多井不同频段比例系数组合,进行空间合理插值得到多井不同频段比例系数数据体,用该比例系数数据体分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到多井空间变比例系数井震匹配处理后的地震数据体;
(3.4)用上述多井平均求取的不同频段统一比例系数组合,分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间统一比例系数井震匹配处理后的地震数据体。
可选的,在本发明一实施例中,所述步骤c中,在频率域进行比例系数求取及应用,包括:
(4.1)在目标时窗内,通过傅里叶变换将合成地震记录数据体和地震数据体,转换为频率域的地震数据体,即频域合成地震记录数据体与频域地震数据体;
(4.2)将频域合成地震记录数据体与频域地震数据体分别按照设计的三角滤波门进行滤波,分别得到多个频段的频域合成地震记录子体和频域地震子体;
(4.3)在井点处求取每个频段的频域合成地震记录子体与频域地震子体的累加振幅值,并计算其累加振幅的比值,将其作为频率域多井不同频段的比例系数组合,同时用多井平均求取频率域不同频段统一比例系数组合;
(4.4)利用上述频率域多井不同频段的比例系数组合,进行空间合理插值得到频率域多井不同频段的比例系数数据体,用该比例系数数据体分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域多井空间变比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体;
(4.5)用上述多井平均求取频率域不同频段统一比例系数组合,分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间统一比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体。
可选的,在本发明一实施例中,所述步骤d中,对井震匹配处理后的地震数据体进行分析检查,包括:井旁地震道与合成地震记录相关性分析、信噪比分析、目标时窗内井震振幅谱相对高低关系分析。
另一方面,本发明实施例提供了一种井震匹配解释性目标处理装置,所述井震匹配解释性目标处理装置,包括:
基础数据准备单元,用于进行基础数据准备:首先解释目的层的顶底层位,以目的层为中心,沿层选取时间窗口,简称时窗,作为目标处理求取比例系数的数据体时窗,分析该时窗内数据体的有效频带,在有效频带基础上选取常相位子波;然后利用常相位子波及优选的井制作合成地震记录;
滤波门设计单元,用于进行滤波门设计:依据地震数据体目的层有效频带,设计多个频段的三角滤波门,其中,前一个滤波门的后斜坡与后一个滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择滤波门的个数;
比例系数处理单元,用于进行比例系数的求取及应用:在目标处理的时窗内,在井点处求取每个频段的合成地震记录与实际地震数据的比例系数;将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;并应用求取的比例系数得到空间统一比例系数井震匹配处理后的地震数据体;
匹配效果监控单元,用于进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,继续上述的步骤,不断优化三角滤波门设计及分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
可选的,在本发明一实施例中,所述基础数据准备单元,进一步用于在多井情况下,进一步优选合成地震记录与井旁地震道相关性较好的井,制作多井合成地震记录,生成相应的合成地震记录数据体。
可选的,在本发明一实施例中,所述比例系数处理单元,进一步用于在时间域进行比例系数求取及应用,包括:
(8.1)将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;
(8.2)在目标处理的时窗内,在井点处求取每个频段合成地震记录子体和相应地震子体的地震属性,并计算其地震属性的比值,然后在多个地震属性比值的基础上,优选多井不同频段比例系数组合,同时用多井平均求取不同频段统一比例系数组合;其中,所述地震属性包括:每个频段合成地震记录子体和相应地震子体的均方根振幅、总振幅、总能量;
(8.3)利用上述优选的多井不同频段比例系数组合,进行空间合理插值得到多井不同频段比例系数数据体,用该比例系数数据体分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到多井空间变比例系数井震匹配处理后的地震数据体;
(8.4)用上述多井平均求取的不同频段统一比例系数组合,分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间统一比例系数井震匹配处理后的地震数据体。
可选的,在本发明一实施例中,所述比例系数处理单元,进一步用于在频率域进行比例系数求取及应用,包括:
(9.1)在目标时窗内,通过傅里叶变换将合成地震记录数据体和地震数据体,转换为频率域的地震数据体,即频域合成地震记录数据体与频域地震数据体;
(9.2)将频域合成地震记录数据体与频域地震数据体分别按照设计的三角滤波门进行滤波,分别得到多个频段的频域合成地震记录子体和频域地震子体;
(9.3)在井点处求取每个频段的频域合成地震记录子体与频域地震子体的累加振幅值,并计算其累加振幅的比值,将其作为频率域多井不同频段的比例系数组合,同时用多井平均求取频率域不同频段统一比例系数组合;
(9.4)利用上述频率域多井不同频段的比例系数组合,进行空间合理插值得到频率域多井不同频段的比例系数数据体,用该比例系数数据体分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域多井空间变比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体;
(9.5)用上述多井平均求取频率域不同频段统一比例系数组合,分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间统一比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体。
可选的,在本发明一实施例中,所述匹配效果监控单元,进一步用于对井震匹配处理后的地震数据体进行分析检查,包括:井旁地震道与合成地震记录相关性分析、信噪比分析、目标时窗内井震振幅谱相对高低关系分析。
上述技术方案具有如下有益效果:在地震数据体上,以目的层为中心,沿层选取合适的时间窗口,在该时窗的有效频带内,选取一个合适的常相位子波,制作井的合成地震记录,将井旁的振幅谱向合成地震记录振幅谱靠拢,得到井震匹配参数,将该参数应用到整个地震数据体上,用来弥补优势频率之外的频率成分受到的压制,减少地震资料延续相位,提高地震资料的分辨能力,同时提高了地震资料井旁道与合成地震记录的吻合程度,以利于提高波阻抗反演的精度,为油藏地球物理及油藏地质研究提供高精度的地震资料。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一种井震匹配解释性目标处理方法流程图;
图2为本发明实施例一种井震匹配解释性目标处理装置结构示意图;
图3为本发明实施例井震匹配处理原理示意图;
图4a为本发明实施例等间距三角滤波门示意图;
图4b为本发明实施例倍频程三角滤波门示意图;
图5为本发明应用实例井震匹配处理前后效果综合对比示意图;
图6为本发明应用实例多井合成地震记录数据体示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种井震匹配解释性目标处理方法及装置,以提高提高地震资料的分辨能力,从而为油藏地球物理及油藏地质研究提供高精度的地震资料。
如图1所示,为本发明实施例一种井震匹配解释性目标处理方法流程图,所述井震匹配解释性目标处理方法,包括:
a、进行基础数据准备:首先解释目的层的顶底层位,以目的层为中心,沿层选取时间窗口,简称时窗,作为目标处理求取比例系数的数据体时窗,分析该时窗内数据体的有效频带,在有效频带基础上选取常相位子波;然后利用常相位子波及优选的井制作合成地震记录;
b、进行滤波门设计:依据地震数据体目的层有效频带,设计多个频段的三角滤波门,其中,前一个滤波门的后斜坡与后一个滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择滤波门的个数;
c、进行比例系数的求取及应用:在目标处理的时窗内,在井点处求取每个频段的合成地震记录与实际地震数据的比例系数;将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;并应用求取的比例系数得到空间统一比例系数井震匹配处理后的地震数据体;
d、进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,继续上述的步骤,不断优化三角滤波门设计及分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
可选的,所述步骤a中,在多井情况下,进一步优选合成地震记录与井旁地震道相关性较好的井,制作多井合成地震记录,生成相应的合成地震记录数据体。
如图3所示,为本发明实施例井震匹配处理原理示意图,如图4a所示,为本发明实施例等间距三角滤波门示意图。可选的,所述步骤c中,在时间域进行比例系数求取及应用,包括:
(3.1)将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;
(3.2)在目标处理的时窗内,在井点处求取每个频段合成地震记录子体和相应地震子体的地震属性,并计算其地震属性的比值,然后在多个地震属性比值的基础上,优选多井不同频段比例系数组合,同时用多井平均求取不同频段统一比例系数组合;其中,所述地震属性包括:每个频段合成地震记录子体和相应地震子体的均方根振幅、总振幅、总能量;
(3.3)利用上述优选的多井不同频段比例系数组合,进行空间合理插值得到多井不同频段比例系数数据体,用该比例系数数据体分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到多井空间变比例系数井震匹配处理后的地震数据体;
(3.4)用上述多井平均求取的不同频段统一比例系数组合,分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间统一比例系数井震匹配处理后的地震数据体。
如图4b所示,为本发明实施例倍频程三角滤波门示意图。可选的,所述步骤c中,在频率域进行比例系数求取及应用,包括:
(4.1)在目标时窗内,通过傅里叶变换将合成地震记录数据体和地震数据体,转换为频率域的地震数据体,即频域合成地震记录数据体与频域地震数据体;
(4.2)将频域合成地震记录数据体与频域地震数据体分别按照设计的三角滤波门进行滤波,分别得到多个频段的频域合成地震记录子体和频域地震子体;
(4.3)在井点处求取每个频段的频域合成地震记录子体与频域地震子体的累加振幅值,并计算其累加振幅的比值,将其作为频率域多井不同频段的比例系数组合,同时用多井平均求取频率域不同频段统一比例系数组合;
(4.4)利用上述频率域多井不同频段的比例系数组合,进行空间合理插值得到频率域多井不同频段的比例系数数据体,用该比例系数数据体分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域多井空间变比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体;
(4.5)用上述多井平均求取频率域不同频段统一比例系数组合,分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间统一比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体。
可选的,在本发明一实施例中,所述步骤d中,对井震匹配处理后的地震数据体进行分析检查,包括:井旁地震道与合成地震记录相关性分析、信噪比分析、目标时窗内井震振幅谱相对高低关系分析。
对应于上述方法实施例,如图2所示,为本发明实施例一种井震匹配解释性目标处理装置结构示意图,所述井震匹配解释性目标处理装置,包括:
基础数据准备单元21,用于进行基础数据准备:首先解释目的层的顶底层位,以目的层为中心,沿层选取时间窗口,简称时窗,作为目标处理求取比例系数的数据体时窗,分析该时窗内数据体的有效频带,在有效频带基础上选取常相位子波;然后利用常相位子波及优选的井制作合成地震记录;
滤波门设计单元22,用于进行滤波门设计:依据地震数据体目的层有效频带,设计多个频段的三角滤波门,其中,前一个滤波门的后斜坡与后一个滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择滤波门的个数;
比例系数处理单元23,用于进行比例系数的求取及应用:在目标处理的时窗内,在井点处求取每个频段的合成地震记录与实际地震数据的比例系数;将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;并应用求取的比例系数得到空间统一比例系数井震匹配处理后的地震数据体;
匹配效果监控单元24,用于进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,继续上述的步骤,不断优化三角滤波门设计及分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
可选的,所述基础数据准备单元21,进一步用于在多井情况下,进一步优选合成地震记录与井旁地震道相关性较好的井,制作多井合成地震记录,生成相应的合成地震记录数据体。
可选的,所述比例系数处理单元23,进一步用于在时间域进行比例系数求取及应用,包括:
(8.1)将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;
(8.2)在目标处理的时窗内,在井点处求取每个频段合成地震记录子体和相应地震子体的地震属性,并计算其地震属性的比值,然后在多个地震属性比值的基础上,优选多井不同频段比例系数组合,同时用多井平均求取不同频段统一比例系数组合;其中,所述地震属性包括:每个频段合成地震记录子体和相应地震子体的均方根振幅、总振幅、总能量;
(8.3)利用上述优选的多井不同频段比例系数组合,进行空间合理插值得到多井不同频段比例系数数据体,用该比例系数数据体分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到多井空间变比例系数井震匹配处理后的地震数据体;
(8.4)用上述多井平均求取的不同频段统一比例系数组合,分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间统一比例系数井震匹配处理后的地震数据体。
可选的,所述比例系数处理单元23,进一步用于在频率域进行比例系数求取及应用,包括:
(9.1)在目标时窗内,通过傅里叶变换将合成地震记录数据体和地震数据体,转换为频率域的地震数据体,即频域合成地震记录数据体与频域地震数据体;
(9.2)将频域合成地震记录数据体与频域地震数据体分别按照设计的三角滤波门进行滤波,分别得到多个频段的频域合成地震记录子体和频域地震子体;
(9.3)在井点处求取每个频段的频域合成地震记录子体与频域地震子体的累加振幅值,并计算其累加振幅的比值,将其作为频率域多井不同频段的比例系数组合,同时用多井平均求取频率域不同频段统一比例系数组合;
(9.4)利用上述频率域多井不同频段的比例系数组合,进行空间合理插值得到频率域多井不同频段的比例系数数据体,用该比例系数数据体分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域多井空间变比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体;
(9.5)用上述多井平均求取频率域不同频段统一比例系数组合,分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间统一比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域井震匹配处理后的地震数据体。
可选的,所述匹配效果监控单元24,进一步用于对井震匹配处理后的地震数据体进行分析检查,包括:井旁地震道与合成地震记录相关性分析、信噪比分析、目标时窗内井震振幅谱相对高低关系分析。
本发明实施例上述方法或装置技术方案具有如下有益效果:在地震数据体上,以目的层为中心,沿层选取合适的时间窗口,在该时窗的有效频带内,选取一个合适的常相位子波,制作井的合成地震记录,将井旁的振幅谱向合成地震记录振幅谱靠拢,得到井震匹配参数,将该参数应用到整个地震数据体上,用来弥补优势频率之外的频率成分受到的压制,减少地震资料延续相位,提高地震资料的分辨能力,同时提高了地震资料井旁道与合成地震记录的吻合程度,以利于提高波阻抗反演的精度,为油藏地球物理及油藏地质研究提供高精度的地震资料。
以下本发明应用实例以时间域求取比例系数为例作如下叙述:
(1)解释研究区目的层的顶底层位,即萨零油组顶界和高四油组底界,在萨零油组至高四油组目的层范围内,分析地震数据体的有效频带,本次实例工区中地震数据体的有效频带为15-75Hz,在该有效频带基础上选取合适的常相位子波,本实例选择75度相位的带通子波:即5-15-70-90Hz频段的带通子波;优选声波、密度测井曲线段较长且测井资料质量较高的井,识别并剔除测井曲线异常值,在研究区中最终优选了合成记录与井旁地震道相关性较好的5口井;然后利用上述常相位带通子波及优选的5口井制作高质量的合成地震记录数据体;如图6所示,为本发明应用实例多井合成地震记录数据体示意图,黑白背景为合成记录数据体,井旁黑白柱为井波阻抗曲线的变密度显示。
(2)按照地震数据体有效频带设计三角滤波门(参见图4a),本发明应用实例三角滤波门按照等间距设计:即0-20-40Hz、20-40-60Hz、40-60-80Hz、60-80-100Hz,总共4个三角滤波门;
(3)将合成地震记录数据体和地震数据体分别按照上述的4个三角滤波门进行分频滤波,分别得到4个频段的合成地震记录子体和地震子体;
(4)在萨零油组至高四油组目的层之间,分别在5口井点处提取每个合成地震记录子体和相应地震子体上的均方根振幅、总振幅、总能量等地震属性,并计算其地震属性的比值,然后在多个地震属性比值的基础上,优选了均方根振幅属性计算的比例系数作为5口井不同频段最佳比例系数组合,同时求取每个频段5口井比例系数的均值,作为不同频段的统一比例系数组合;
(5)利用上述优选的5口井不同频段最佳比例系数组合,进行空间插值得到5口井不同频段比例系数数据体,用该比例系数数据体分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间变比例系数井震匹配处理后的地震数据体;
(6)用上述多井平均求取的不同频段的统一比例系数组合,分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间统一比例系数井震匹配处理后的地震数据体;
(7)在目的层时窗内,利用傅里叶变换,分别计算井震匹配处理后的地震数据体与合成地震记录数据体的振幅谱,检查这两个数据体频谱的相对高低关系,同时分析井旁地震道与合成地震记录相关性、信噪比等,从空间变比例系数和统一比例系数井震匹配处理数据体中,优选井震最佳匹配处理地震数据体,即地震井旁道振幅谱向合成地震记录振幅谱形态逼近、波形逼近,达到井震最佳匹配的效果;如图5所示,为本发明应用实例井震匹配处理前后效果综合对比示意图,上半部分的为三个地震剖面对比,下半部分为目的层的井旁地震道频谱对比。
(8)如果没有达到井震最佳匹配的满意效果,继续上述的步骤,不断优化三角滤波门设计及分频段地震子体比例系数估算等环节,最终获得预期的满意的井震匹配目标处理数据体。
本发明应用实例通过井震匹配处理,地震资料的振幅谱向合成地震记录的振幅谱形态逼近、波形逼近,达到了井震最佳匹配的效果,在地震资料的有效频带内,提高了地震资料的分辨能力,同时提高了地震资料井旁道与合成地震记录的吻合程度,以利于提高波阻抗反演的精度,为油藏地球物理及油藏地质研究提供高精度的地震资料。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrativelogical block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrativecomponents),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列(FPGA)或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种井震匹配解释性目标处理方法,其特征在于,所述井震匹配解释性目标处理方法,包括:
a、进行基础数据准备:针对目的层段,分析地震数据体的有效频带,在该有效频带基础上选取常相位子波;然后利用常相位子波及测井资料制作合成地震记录,生成相应的合成地震记录数据体;
b、进行三角滤波门设计:依据地震数据体目的层段的有效频带,设计多个频段的三角滤波门,其中,前一个三角滤波门的后斜坡与后一个三角滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择三角滤波门的个数;
c、在时间域进行比例系数求取及应用:
(2.1)将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;
(2.2)针对目的层段,在井点处求取每个频段合成地震记录子体和地震子体的地震属性,并计算其地震属性的比值,然后在多个地震属性比值的基础上,选择多井不同频段空间变比例系数,同时用多井平均求取不同频段空间统一比例系数;其中,所述地震属性包括:每个频段合成地震记录子体和地震子体的均方根振幅、总振幅、总能量;
(2.3)利用上述选择的多井不同频段空间变比例系数,进行空间合理插值得到多井不同频段空间变比例系数数据体,用该空间变比例系数数据体分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间变比例系数井震匹配处理后的地震数据体;
(2.4)用上述多井平均求取的不同频段空间统一比例系数,分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间统一比例系数井震匹配处理后的地震数据体;
d、进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,继续步骤b和步骤c,不断优化三角滤波门设计及分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
2.如权利要求1所述井震匹配解释性目标处理方法,其特征在于,所述步骤d中,
对井震匹配处理后的地震数据体进行分析检查,包括:对目的层段的井震匹配处理后的地震数据体与合成地震记录数据体的振幅谱相对高低关系进行分析,对井旁地震道与合成地震记录的相关性和信噪比进行分析。
3.一种井震匹配解释性目标处理方法,其特征在于,所述井震匹配解释性目标处理方法,包括:
a、进行基础数据准备:针对目的层段,分析地震数据体的有效频带,在该有效频带基础上选取常相位子波;然后利用常相位子波及测井资料制作合成地震记录,生成相应的合成地震记录数据体;
b、进行三角滤波门设计:依据地震数据体目的层段的有效频带,设计多个频段的三角滤波门,其中,前一个三角滤波门的后斜坡与后一个三角滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择三角滤波门的个数;
c、在频率域进行比例系数求取及应用,包括:
(3.1)针对目的层段,通过傅里叶变换将合成地震记录数据体和地震数据体,转换为频率域的地震数据体,即频域合成地震记录数据体与频域地震数据体;
(3.2)将频域合成地震记录数据体与频域地震数据体分别按照设计的三角滤波门进行滤波,分别得到多个频段的频域合成地震记录子体和频域地震子体;
(3.3)在井点处求取每个频段的频域合成地震记录子体与频域地震子体的累加振幅值,并计算其累加振幅的比值,将其作为频率域多井不同频段的空间变比例系数,同时用多井平均求取频率域不同频段空间统一比例系数;
(3.4)利用上述频率域多井不同频段的空间变比例系数,进行空间合理插值得到频率域多井不同频段的空间变比例系数数据体,用该空间变比例系数数据体分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间变比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域空间变比例系数井震匹配处理后的地震数据体;
(3.5)用上述多井平均求取频率域不同频段空间统一比例系数,分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间统一比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域空间统一比例系数井震匹配处理后的地震数据体
d、进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,继续步骤b和步骤c,不断优化三角滤波门设计及分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
4.如权利要求3中所述井震匹配解释性目标处理方法,其特征在于,所述步骤d中,
对井震匹配处理后的地震数据体进行分析检查,包括:对目的层段的井震匹配处理后的地震数据体与合成地震记录数据体的振幅谱相对高低关系进行分析,对井旁地震道与合成地震记录的相关性和信噪比进行分析。
5.一种井震匹配解释性目标处理装置,其特征在于,所述井震匹配解释性目标处理装置,包括:
基础数据准备单元,用于进行基础数据准备:针对目的层段,分析地震数据体的有效频带,在该有效频带基础上选取常相位子波;然后利用常相位子波及测井资料制作合成地震记录,生成相应的合成地震记录数据体;
滤波门设计单元,用于进行三角滤波门设计:依据地震数据体目的层有效频带,设计多个频段的三角滤波门,其中,前一个三角滤波门的后斜坡与后一个三角滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择三角滤波门的个数;
比例系数处理单元,用于在时间域进行比例系数求取及应用,包括:
(6.1)将合成地震记录数据体和地震数据体分别按照设计的三角滤波门进行分频滤波,分别得到多个频段的合成地震记录子体和地震子体;
(6.2)针对目的层段,在井点处求取每个频段合成地震记录子体和相应地震子体的地震属性,并计算其地震属性的比值,然后在多个地震属性比值的基础上,选择多井不同频段空间变比例系数,同时用多井平均求取不同频段空间统一比例系数;其中,所述地震属性包括:每个频段合成地震记录子体和相应地震子体的均方根振幅、总振幅、总能量;
(6.3)利用上述选择的多井不同频段空间变比例系数,进行空间合理插值得到多井不同频段空间变比例系数数据体,用该空间变比例系数数据体分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间变比例系数井震匹配处理后的地震数据体;
(6.4)用上述多井平均求取的不同频段空间统一比例系数,分别乘以相应频段的地震子体,然后将处理后的不同频段地震子体相加,得到空间统一比例系数井震匹配处理后的地震数据体;
匹配效果监控单元,用于进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,滤波门设计单元不断优化三角滤波门设计,比例系数处理单元继续分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
6.如权利要求5所述井震匹配解释性目标处理装置,其特征在于,所述匹配效果监控单元,进一步用于对井震匹配处理后的地震数据体进行分析检查,包括:对目的层段的井震匹配处理后的地震数据体与合成地震记录数据体的振幅谱相对高低关系进行分析,对井旁地震道与合成地震记录的相关性和信噪比进行分析。
7.一种井震匹配解释性目标处理装置,其特征在于,所述井震匹配解释性目标处理装置,包括:
基础数据准备单元,用于进行基础数据准备:针对目的层段,分析地震数据体的有效频带,在该有效频带基础上选取常相位子波;然后利用常相位子波及测井资料制作合成地震记录,生成相应的合成地震记录数据体;
滤波门设计单元,用于进行三角滤波门设计:依据地震数据体目的层有效频带,设计多个频段的三角滤波门,其中,前一个三角滤波门的后斜坡与后一个三角滤波门的前斜坡宽度保证相等,各三角滤波门在斜坡中点处相互重叠,其间距是等间距或者是倍频程间距;设计的三角滤波门频段要将地震有效频带包含在内,依据频谱形态的复杂程度选择三角滤波门的个数;
比例系数处理单元,用于在频率域进行比例系数求取及应用,包括:
(7.1)针对目的层段,通过傅里叶变换将合成地震记录数据体和地震数据体,转换为频率域的地震数据体,即频域合成地震记录数据体与频域地震数据体;
(7.2)将频域合成地震记录数据体与频域地震数据体分别按照设计的三角滤波门进行滤波,分别得到多个频段的频域合成地震记录子体和频域地震子体;
(7.3)在井点处求取每个频段的频域合成地震记录子体与频域地震子体的累加振幅值,并计算其累加振幅的比值,将其作为频率域多井不同频段的空间变比例系数,同时用多井平均求取频率域不同频段空间统一比例系数;
(7.4)利用上述频率域多井不同频段的空间变比例系数,进行空间合理插值得到频率域多井不同频段的比例系数空间变比例系数数据体,用该空间变比例系数数据体分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间变比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域空间变比例系数井震匹配处理后的地震数据体;
(7.5)用上述多井平均求取频率域不同频段空间统一比例系数,分别乘以相应频段的频域地震子体,然后将处理后的不同频段的频域地震子体相加,得到频率域空间统一比例系数井震匹配处理后的地震数据体;在此基础上,应用反傅里叶变换,得到相应的时间域空间统一比例系数井震匹配处理后的地震数据体;
匹配效果监控单元,用于进行匹配效果监控:对井震匹配处理后的地震数据体进行分析检查,如果没有达到预期井震最佳匹配的满意效果,滤波门设计单元不断优化三角滤波门设计,比例系数处理单元继续分频段地震子体比例系数估算,最终获得预期的井震匹配目标处理数据体。
8.如权利要求7所述井震匹配解释性目标处理装置,其特征在于,所述匹配效果监控单元,进一步用于对井震匹配处理后的地震数据体进行分析检查,包括:对目的层段的井震匹配处理后的地震数据体与合成地震记录数据体的振幅谱相对高低关系进行分析,对井旁地震道与合成地震记录的相关性和信噪比进行分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210332849.XA CN103675904B (zh) | 2012-09-10 | 2012-09-10 | 一种井震匹配解释性目标处理方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210332849.XA CN103675904B (zh) | 2012-09-10 | 2012-09-10 | 一种井震匹配解释性目标处理方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675904A CN103675904A (zh) | 2014-03-26 |
CN103675904B true CN103675904B (zh) | 2016-08-17 |
Family
ID=50314023
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210332849.XA Active CN103675904B (zh) | 2012-09-10 | 2012-09-10 | 一种井震匹配解释性目标处理方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675904B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536042B (zh) * | 2014-12-24 | 2017-01-25 | 中国石油天然气股份有限公司 | 一种二维叠后地震资料振幅补偿方法及装置 |
CN104698502A (zh) * | 2015-03-22 | 2015-06-10 | 山西煤炭进出口集团科技发展有限公司 | 由井出发的瞬变电磁勘探地形校正处理方法 |
CN104820243B (zh) * | 2015-04-20 | 2017-04-12 | 中国石油天然气股份有限公司 | 一种三维叠后地震资料振幅补偿方法及装置 |
CN106610505B (zh) * | 2016-12-29 | 2019-03-22 | 中国石油大学(华东) | 一种基于dtw和aba联合的井震资料匹配方法 |
CN107356970B (zh) * | 2017-09-11 | 2020-05-01 | 中国石油大学(华东) | 一种高精度的井震数据匹配方法 |
CN109212603A (zh) * | 2018-11-08 | 2019-01-15 | 北京邮电大学 | 一种基于gs变换滤波和emd去噪的地震信号检测算法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9206681D0 (en) * | 1992-03-27 | 1992-05-13 | Geco As | A method of making a synthetic seismogram |
CN1412575A (zh) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | 一种基于多相介质理论的油气检测方法 |
US7447114B2 (en) * | 2005-08-31 | 2008-11-04 | Geotrace Technologies, Inc. | Non-linear seismic trace matching to well logs |
CN102169190B (zh) * | 2011-01-06 | 2012-10-24 | 中国科学院地质与地球物理研究所 | 一种调制补充子空间的井约束叠前弹性参数反演方法 |
-
2012
- 2012-09-10 CN CN201210332849.XA patent/CN103675904B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103675904A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675904B (zh) | 一种井震匹配解释性目标处理方法及装置 | |
CN104950334B (zh) | 一种预测储层分布的方法及装置 | |
CN104395779B (zh) | 用于估计和减弱地震数据中的噪声的系统和方法 | |
CN111123354B (zh) | 基于频变反射振幅衰减预测致密气层的方法及设备 | |
CN105425289B (zh) | 确定低频波阻抗的方法和装置 | |
CN106249299B (zh) | 强反射屏蔽下薄层弱反射地震能量恢复方法及装置 | |
CN105863628B (zh) | 一种油田开发阶段储层精细预测方法 | |
CN109001813B (zh) | 一种压制多次波的方法、装置及系统 | |
CN105301636B (zh) | 速度模型的建立方法和装置 | |
CN108445538A (zh) | 基于反射地震资料建立深度域层q模型的方法和系统 | |
CN107315194A (zh) | 地震叠前时间偏移速度的确定方法和装置 | |
EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
CN106707342B (zh) | 共炮点道集多级面波压制方法和装置 | |
CN104820243B (zh) | 一种三维叠后地震资料振幅补偿方法及装置 | |
CN106646603B (zh) | 一种实际时移地震资料处理差异的可靠性判断方法 | |
CN107643539A (zh) | 一种基于煤层地震响应特征分析剥离强屏蔽层的方法 | |
CN106033125A (zh) | 压制叠前大角度道集干涉的提频方法 | |
Wei et al. | Deep learning-based P-and S-wave separation for multicomponent vertical seismic profiling | |
CN105954803A (zh) | 叠后地震反演方法及装置 | |
CN106125139B (zh) | 一种三维地震数据处理方法及系统 | |
CN109471162A (zh) | 层间多次波处理方法、系统、电子设备及可读介质 | |
CN108427140B (zh) | 一种用于小尺度缝洞储集体地震识别的方法 | |
CN108254785A (zh) | 储层确定方法和装置 | |
CN108693558B (zh) | 地震数据处理方法和装置 | |
CN109856672A (zh) | 基于深度波数谱的瞬变波包提取方法、存储介质与终端 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |