CN111435174A - 强反射地区地震资料振幅补偿方法及装置 - Google Patents

强反射地区地震资料振幅补偿方法及装置 Download PDF

Info

Publication number
CN111435174A
CN111435174A CN201910034579.6A CN201910034579A CN111435174A CN 111435174 A CN111435174 A CN 111435174A CN 201910034579 A CN201910034579 A CN 201910034579A CN 111435174 A CN111435174 A CN 111435174A
Authority
CN
China
Prior art keywords
time window
seismic data
strong reflection
seismic
reflection stratum
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
CN201910034579.6A
Other languages
English (en)
Other versions
CN111435174B (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 CN201910034579.6A priority Critical patent/CN111435174B/zh
Publication of CN111435174A publication Critical patent/CN111435174A/zh
Application granted granted Critical
Publication of CN111435174B publication Critical patent/CN111435174B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking

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所示,该方法可以包括:
步骤101、获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料;
步骤102、根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数;
步骤103、根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数;
步骤104、根据所述第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数;
步骤105、根据所述第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿。
由图1所示可以得知,本发明实施例通过获取强反射地区地震资料,将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数,根据第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数,最后根据第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿,有效补偿了由于强反射地层屏蔽作用而衰减的振幅,提高了预测强反射地层下伏奥陶系碳酸盐岩缝洞储层的准确性。
具体实施时,获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料。
发明人发现,在强反射地区,强反射地层在空间上的厚度、岩性变化剧烈,地震波在通过强反射地层时,强反射地层会在空间上对地震波产生不同程度的屏蔽作用,这使得在强反射地区地震资料中,除了存在由几何扩散、地层吸收导致的时间方向振幅衰减和由近地表结构差异造成的空间方向振幅差异,还存在由于强反射地层屏蔽作用导致的空间非一致性振幅衰减。仅采用现有的地震资料振幅补偿方法,无法补偿由于强反射地层屏蔽作用而衰减的振幅,进而在预测强反射地层下伏奥陶系碳酸盐岩缝洞储层时可能出现错误,致使油气勘探失利。本发明实施例提供的强反射地区地震资料振幅补偿方法,将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,并通过强反射地层上时窗和强反射地层下时窗的地震数据计算得到强反射地区地震资料的振幅补偿系数,从而补偿强反射地区地震资料中由于强反射地层屏蔽作用而衰减的振幅,提高了强反射地层下伏奥陶系碳酸盐岩缝洞储层预测的准确率。
实施例中,首先获取强反射地区地震资料,然后将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料。需要说明的是,划分时窗时一定要精细,强反射地层上时窗和强反射地层下时窗内不能包含强反射地层地层,但尽可能包含强反射地层地层上下所有地层。
具体实施时,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。
实施例中,首先根据强反射地层上时窗内的地震资料和时窗长度,确定强反射地层上时窗内的地震资料的第一振幅算数平均值。按如下公式确定强反射地层上时窗内的地震资料的第一振幅算数平均值:
Figure BDA0001945429280000051
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层上时窗内的地震资料中第i道第j样点的值,Ai为强反射地层上时窗内的地震资料的第一振幅算数平均值,N1为强反射地层上时窗内的时窗长度。
实施例中,确定强反射地层上时窗内的地震资料的第一振幅算数平均值之后,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第二振幅算数平均值。按如下公式确定强反射地层上时窗内的地震资料的第二振幅算数平均值:
Figure BDA0001945429280000052
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层上时窗内的地震资料中第i道第j样点的值,N1为强反射地层上时窗内的时窗长度,M为地震道个数,A为强反射地层上时窗内的地震资料的第二振幅算数平均值。
实施例中,确定强反射地层上时窗内的地震资料的第二振幅算数平均值之后,根据所述第一振幅算数平均值和第二振幅算数平均值,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。按如下公式确定强反射地层上时窗内的地震资料的第一振幅补偿系数:
Figure BDA0001945429280000061
其中,i为地震道序号,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,Ai为强反射地层上时窗内的地震资料的第一振幅算数平均值,A为强反射地层上时窗内的地震资料的第二振幅算数平均值。
具体实施时,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。
实施例中,首先根据强反射地层下时窗内的地震资料和时窗长度,确定强反射地层下时窗内的地震资料的第三振幅算数平均值。按如下公式确定强反射地层下时窗内的地震资料的第三振幅算数平均值:
Figure BDA0001945429280000062
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层下时窗内的地震资料中第i道第j样点的值,Bi为强反射地层下时窗内的地震资料的第三振幅算数平均值,N2为强反射地层下时窗内的时窗长度,N为各地震道中地震资料样点个数。
实施例中,确定强反射地层下时窗内的地震资料的第三振幅算数平均值之后,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第四振幅算数平均值。按如下公式确定强反射地层下时窗内的地震资料的第四振幅算数平均值:
Figure BDA0001945429280000063
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层下时窗内的地震资料中第i道第j样点的值,N2为强反射地层下时窗内的时窗长度,N为各地震道中地震资料样点个数,M为地震道个数,B为强反射地层下时窗内的地震资料的第四振幅算数平均值。
实施例中,确定强反射地层下时窗内的地震资料的第四振幅算数平均值之后,根据第三振幅算数平均值和第四振幅算数平均值,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。按如下公式确定强反射地层下时窗内的地震资料的第二振幅补偿系数:
Figure BDA0001945429280000071
其中,i为地震道序号,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,Bi为强反射地层下时窗内的地震资料的第三振幅算数平均值,B为强反射地层下时窗内的地震资料的第四振幅算数平均值。
具体实施时,根据第一振幅补偿系数和第二振幅补偿系数,确定强反射地区地震资料的第三振幅补偿系数。
实施例中,按如下公式确定强反射地区地震资料的第三振幅补偿系数:
Figure BDA0001945429280000072
其中,i为地震道序号,j为各地震道中各样点序号,OP3ij为强反射地区地震资料中第i道第j样点的第三振幅补偿系数,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,N为各地震道中地震资料样点个数。
具体实施时,根据所述第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿。
实施例中,第三振幅补偿系数为振幅补偿因子,首先对强反射地区地震资料进行动校正,得到动校正后的共中心点(Common Middle Point,CMP)道集,然后将第三振幅补偿因子与动校正后的CMP道集的振幅进行相乘,从而补偿强反射地区地震资料衰减的振幅。
下面给出一个具体实施例,说明本发明实施例中强反射地区地震资料振幅补偿方法的具体应用。以原有振幅补偿技术处理为基础,在叠后数据上将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,划分方法如图2所示,然后计算强反射地区地震资料中的每一道的振幅算数平均值,计算出的间隔25道的振幅的算数平均值曲线如图3中的A所示,单位为分贝,可以看出对强反射地层下时窗的地震数据振幅有明显屏蔽作用,表现为强反射地层下时窗地震数据部分道的振幅算数平均值曲线有不同程度变弱,整体曲线有发散现象。计算强反射地区地震资料中的所有道的振幅算数平均值,表示整体数据的振幅衰减规律,用于期望输出,计算结果如图3中的B所示。根据计算出的强反射地区地震资料中的每一道的振幅算数平均值和所有道的振幅算数平均值,计算强反射地区地震资料的振幅补偿因子,如图3中的C所示。最后用计算出的振幅补偿因子与原有补偿技术处理后的动校后的CMP道集做积,即完成基于强反射地层的振幅补偿,结果如图3中的D所示,通过对比图3中的A和图3中的D可以看出,强反射地层下时窗部分道的振幅算数平均值曲线得到较好补偿,整体曲线变得更集中。图4是利用强反射地区地震资料振幅补偿方法进行振幅补偿的结果图,图4中的A是强反射地区地震资料振幅补偿前的偏移剖面,图4中的B是强反射地区地震资料振幅补偿后的偏移剖面,通过对比可以看出,应用了强反射地区地震资料振幅补偿方法的偏移剖面,强反射地区地震资料中由于强反射地层屏蔽作用而衰减的振幅得到较好的补偿,其振幅关系较好的反映了地下地层的真实反射特征。
基于同一发明构思,本发明实施例还提供了一种强反射地区地震资料振幅补偿装置,如下面的实施例所述。由于这些解决问题的原理与强反射地区地震资料振幅补偿方法相似,因此装置的实施可以参见方法的实施,重复之处不再赘述。
图5为本发明实施例中强反射地区地震资料振幅补偿装置的结构图,如图5所示,该装置包括:
时窗划分模块501,用于获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料;
上时窗模块502,用于根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数;
下时窗模块503,用于根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数;
补偿因子模块504,用于根据所述第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数;
振幅补偿模块505,用于根据所述第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿。
一个实施例中,上时窗模块502进一步用于:
根据强反射地层上时窗内的地震资料和时窗长度,确定强反射地层上时窗内的地震资料的第一振幅算数平均值;
根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第二振幅算数平均值;
根据所述第一振幅算数平均值和第二振幅算数平均值,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。
一个实施例中,下时窗模块503进一步用于:
根据强反射地层下时窗内的地震资料和时窗长度,确定强反射地层下时窗内的地震资料的第三振幅算数平均值;
根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第四振幅算数平均值;
根据所述第三振幅算数平均值和第四振幅算数平均值,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。
一个实施例中,补偿因子模块504进一步用于,按如下公式确定强反射地区地震资料的第三振幅补偿系数:
Figure BDA0001945429280000091
其中,i为地震道序号,j为各地震道中各样点序号,OP3ij为强反射地区地震资料中第i道第j样点的第三振幅补偿系数,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,N为各地震道中地震资料样点个数。
一个实施例中,第三振幅补偿系数为振幅补偿因子。
综上所述,本发明实施例通过获取强反射地区地震资料,将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数,根据第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数,最后根据第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿,有效补偿了由于强反射地层屏蔽作用而衰减的振幅,提高了预测强反射地层下伏奥陶系碳酸盐岩缝洞储层的准确性。本发明实施例通过将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,并根据强反射地层上时窗和强反射地层下时窗的地震数据计算得到强反射地区地震资料的振幅补偿系数,有效补偿了强反射地区地震资料中由于强反射地层屏蔽作用而衰减的振幅,提高了强反射地层下伏奥陶系碳酸盐岩缝洞储层预测的准确率。
本领域内的技术人员应明白,本发明的实施例可提供为方法、***、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (18)

1.一种强反射地区地震资料振幅补偿方法,其特征在于,包括:
获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料;
根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数;
根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数;
根据所述第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数;
根据所述第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿。
2.如权利要求1所述的方法,其特征在于,所述根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数,包括:
根据强反射地层上时窗内的地震资料和时窗长度,确定强反射地层上时窗内的地震资料的第一振幅算数平均值;
根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第二振幅算数平均值;
根据所述第一振幅算数平均值和第二振幅算数平均值,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。
3.如权利要求2所述的方法,其特征在于,按如下公式确定强反射地层上时窗内的地震资料的第一振幅算数平均值:
Figure FDA0001945429270000011
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层上时窗内的地震资料中第i道第j样点的值,Ai为强反射地层上时窗内的地震资料的第一振幅算数平均值,N1为强反射地层上时窗内的时窗长度。
4.如权利要求2所述的方法,其特征在于,按如下公式确定强反射地层上时窗内的地震资料的第二振幅算数平均值:
Figure FDA0001945429270000021
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层上时窗内的地震资料中第i道第j样点的值,N1为强反射地层上时窗内的时窗长度,M为地震道个数,A为强反射地层上时窗内的地震资料的第二振幅算数平均值。
5.如权利要求2所述的方法,其特征在于,按如下公式确定强反射地层上时窗内的地震资料的第一振幅补偿系数:
Figure FDA0001945429270000022
其中,i为地震道序号,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,Ai为强反射地层上时窗内的地震资料的第一振幅算数平均值,A为强反射地层上时窗内的地震资料的第二振幅算数平均值。
6.如权利要求1所述的方法,其特征在于,所述根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数,包括:
根据强反射地层下时窗内的地震资料和时窗长度,确定强反射地层下时窗内的地震资料的第三振幅算数平均值;
根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第四振幅算数平均值;
根据所述第三振幅算数平均值和第四振幅算数平均值,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。
7.如权利要求6所述的方法,其特征在于,按如下公式确定强反射地层下时窗内的地震资料的第三振幅算数平均值:
Figure FDA0001945429270000023
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层下时窗内的地震资料中第i道第j样点的值,Bi为强反射地层下时窗内的地震资料的第三振幅算数平均值,N2为强反射地层下时窗内的时窗长度,N为各地震道中地震资料样点个数。
8.如权利要求6所述的方法,其特征在于,按如下公式确定强反射地层下时窗内的地震资料的第四振幅算数平均值:
Figure FDA0001945429270000031
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层下时窗内的地震资料中第i道第j样点的值,N2为强反射地层下时窗内的时窗长度,N为各地震道中地震资料样点个数,M为地震道个数,B为强反射地层下时窗内的地震资料的第四振幅算数平均值。
9.如权利要求6所述的方法,其特征在于,按如下公式确定强反射地层下时窗内的地震资料的第二振幅补偿系数:
Figure FDA0001945429270000032
其中,i为地震道序号,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,Bi为强反射地层下时窗内的地震资料的第三振幅算数平均值,B为强反射地层下时窗内的地震资料的第四振幅算数平均值。
10.如权利要求1所述的方法,其特征在于,按如下公式确定所述强反射地区地震资料的第三振幅补偿系数:
Figure FDA0001945429270000033
其中,i为地震道序号,j为各地震道中各样点序号,OP3ij为强反射地区地震资料中第i道第j样点的第三振幅补偿系数,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,N为各地震道中地震资料样点个数。
11.如权利要求1所述的方法,其特征在于,所述第三振幅补偿系数为振幅补偿因子。
12.一种强反射地区地震资料振幅补偿装置,其特征在于,包括:
时窗划分模块,用于获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料;
上时窗模块,用于根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数;
下时窗模块,用于根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数;
补偿因子模块,用于根据所述第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数;
振幅补偿模块,用于根据所述第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿。
13.如权利要求12所述的装置,其特征在于,所述上时窗模块进一步用于:
根据强反射地层上时窗内的地震资料和时窗长度,确定强反射地层上时窗内的地震资料的第一振幅算数平均值;
根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第二振幅算数平均值;
根据所述第一振幅算数平均值和第二振幅算数平均值,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。
14.如权利要求12所述的装置,其特征在于,所述下时窗模块进一步用于:
根据强反射地层下时窗内的地震资料和时窗长度,确定强反射地层下时窗内的地震资料的第三振幅算数平均值;
根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第四振幅算数平均值;
根据所述第三振幅算数平均值和第四振幅算数平均值,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。
15.如权利要求12所述的装置,其特征在于,所述补偿因子模块进一步用于,按如下公式确定所述强反射地区地震资料的第三振幅补偿系数:
Figure FDA0001945429270000041
其中,i为地震道序号,j为各地震道中各样点序号,OP3ij为强反射地区地震资料中第i道第j样点的第三振幅补偿系数,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,N为各地震道中地震资料样点个数。
16.如权利要求12所述的装置,其特征在于,所述第三振幅补偿系数为振幅补偿因子。
17.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至11任一所述方法。
18.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至11任一所述方法的计算机程序。
CN201910034579.6A 2019-01-15 2019-01-15 强反射地区地震资料振幅补偿方法及装置 Active CN111435174B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910034579.6A CN111435174B (zh) 2019-01-15 2019-01-15 强反射地区地震资料振幅补偿方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910034579.6A CN111435174B (zh) 2019-01-15 2019-01-15 强反射地区地震资料振幅补偿方法及装置

Publications (2)

Publication Number Publication Date
CN111435174A true CN111435174A (zh) 2020-07-21
CN111435174B CN111435174B (zh) 2023-04-25

Family

ID=71580682

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910034579.6A Active CN111435174B (zh) 2019-01-15 2019-01-15 强反射地区地震资料振幅补偿方法及装置

Country Status (1)

Country Link
CN (1) CN111435174B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112379429A (zh) * 2020-11-02 2021-02-19 中国石油天然气集团有限公司 地震数据的振幅补偿方法及装置
CN113572714A (zh) * 2021-06-29 2021-10-29 紫光展锐(重庆)科技有限公司 信号补偿方法、装置、设备及存储介质、芯片、模组设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6757216B1 (en) * 2003-05-15 2004-06-29 Exxonmobil Upstream Research Company Method for post processing compensation of amplitude for misaligned and misstacked offset seismic data
CN104375188A (zh) * 2014-11-19 2015-02-25 中国石油天然气集团公司 一种地震波透射衰减补偿方法及装置
CN104536042A (zh) * 2014-12-24 2015-04-22 中国石油天然气股份有限公司 一种二维叠后地震资料振幅补偿方法及装置
CN104820243A (zh) * 2015-04-20 2015-08-05 中国石油天然气股份有限公司 一种三维叠后地震资料振幅补偿方法及装置
CN106569277A (zh) * 2015-10-12 2017-04-19 中国石油化工股份有限公司 地震属性叠后振幅补偿方法
CN108181649A (zh) * 2017-11-23 2018-06-19 中国石油天然气集团公司 一种介质分界面振幅补偿方法及装置
US20180321405A1 (en) * 2017-05-03 2018-11-08 Saudi Arabian Oil Company Refraction-based surface-consistent amplitude compensation and deconvolution

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6757216B1 (en) * 2003-05-15 2004-06-29 Exxonmobil Upstream Research Company Method for post processing compensation of amplitude for misaligned and misstacked offset seismic data
CN104375188A (zh) * 2014-11-19 2015-02-25 中国石油天然气集团公司 一种地震波透射衰减补偿方法及装置
CN104536042A (zh) * 2014-12-24 2015-04-22 中国石油天然气股份有限公司 一种二维叠后地震资料振幅补偿方法及装置
CN104820243A (zh) * 2015-04-20 2015-08-05 中国石油天然气股份有限公司 一种三维叠后地震资料振幅补偿方法及装置
CN106569277A (zh) * 2015-10-12 2017-04-19 中国石油化工股份有限公司 地震属性叠后振幅补偿方法
US20180321405A1 (en) * 2017-05-03 2018-11-08 Saudi Arabian Oil Company Refraction-based surface-consistent amplitude compensation and deconvolution
CN108181649A (zh) * 2017-11-23 2018-06-19 中国石油天然气集团公司 一种介质分界面振幅补偿方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周星: "数据驱动下的火成岩能量屏蔽补偿方法", 《 中国石油学会2017年物探技术研讨会论文集》 *
姚军: "三维叠后地震资料振幅补偿方法在地震沉积学研究中的应用", 《 CPS/SEG北京2018国际地球物理会议暨展览电子论文集》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112379429A (zh) * 2020-11-02 2021-02-19 中国石油天然气集团有限公司 地震数据的振幅补偿方法及装置
CN113572714A (zh) * 2021-06-29 2021-10-29 紫光展锐(重庆)科技有限公司 信号补偿方法、装置、设备及存储介质、芯片、模组设备
CN113572714B (zh) * 2021-06-29 2022-09-30 紫光展锐(重庆)科技有限公司 信号补偿方法、装置、设备及存储介质、芯片、模组设备

Also Published As

Publication number Publication date
CN111435174B (zh) 2023-04-25

Similar Documents

Publication Publication Date Title
Eide et al. Seismic interpretation of sill complexes in sedimentary basins: implications for the sub-sill imaging problem
EP3259619B1 (en) Method and system of processing seismic data by providing surface aperture common image gathers
US20170184748A1 (en) A method and a computing system for seismic imaging a geological formation
EA020635B1 (ru) Способ сейсмической разведки для поиска углеводородов с применением поля средних скоростей сейсмической волны
CN109188520B (zh) 薄储层厚度预测方法及装置
CN105425289A (zh) 确定低频波阻抗的方法和装置
CN108957532B (zh) 储层预测方法及装置
CN111435174A (zh) 强反射地区地震资料振幅补偿方法及装置
Farrugia et al. Validation of 1D earthquake site characterization methods with observed earthquake site amplification in Alberta, Canada
CN106574980A (zh) 用于地下地质体的岩石性质估计的***和方法
US7382683B1 (en) Computing an absorption parameter for a mode-converted seismic wave
CN111381280A (zh) 预测储层含烃饱和度的方法及装置
CN109991663B (zh) 工区地震速度体校正方法和装置
CN112198547A (zh) 深层或超深层地震资料处理方法及装置
CN103454679B (zh) 一种地震反射同相轴等时性的确定方法及***
US20230228894A1 (en) Target-oriented seismic acquisition method and apparatus, medium and device
US20210231820A1 (en) Computer-implemented method for determining a velocity image of a domain of the subsurface structural geology in an oil and gas reservoir
CN113138412A (zh) 深层页岩气孔隙度地震预测方法及装置
CN112147700A (zh) 速度异常区的低频模型构建方法及***
CN113296153B (zh) 轴对称介质各向异性参数确定方法及装置
NO20190217A1 (en) Correcting biases in microseismic-event data
CN112711071B (zh) 地层倾角校正方法及装置
CN113009579B (zh) 地震数据反演方法及装置
CN107589457A (zh) 测井资料尺度粗化的方法和装置
CN114396257B (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