CN111435174A - 强反射地区地震资料振幅补偿方法及装置 - Google Patents
强反射地区地震资料振幅补偿方法及装置 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000015572 biosynthetic process Effects 0.000 claims description 21
- 238000005070 sampling Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 abstract description 15
- 238000010586 diagram Methods 0.000 description 14
- 230000002238 attenuated effect Effects 0.000 description 13
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 10
- 239000011435 rock Substances 0.000 description 8
- 238000009792 diffusion process Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting 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所示可以得知,本发明实施例通过获取强反射地区地震资料,将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数,根据第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数,最后根据第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿,有效补偿了由于强反射地层屏蔽作用而衰减的振幅,提高了预测强反射地层下伏奥陶系碳酸盐岩缝洞储层的准确性。
具体实施时,获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料。
发明人发现,在强反射地区,强反射地层在空间上的厚度、岩性变化剧烈,地震波在通过强反射地层时,强反射地层会在空间上对地震波产生不同程度的屏蔽作用,这使得在强反射地区地震资料中,除了存在由几何扩散、地层吸收导致的时间方向振幅衰减和由近地表结构差异造成的空间方向振幅差异,还存在由于强反射地层屏蔽作用导致的空间非一致性振幅衰减。仅采用现有的地震资料振幅补偿方法,无法补偿由于强反射地层屏蔽作用而衰减的振幅,进而在预测强反射地层下伏奥陶系碳酸盐岩缝洞储层时可能出现错误,致使油气勘探失利。本发明实施例提供的强反射地区地震资料振幅补偿方法,将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,并通过强反射地层上时窗和强反射地层下时窗的地震数据计算得到强反射地区地震资料的振幅补偿系数,从而补偿强反射地区地震资料中由于强反射地层屏蔽作用而衰减的振幅,提高了强反射地层下伏奥陶系碳酸盐岩缝洞储层预测的准确率。
实施例中,首先获取强反射地区地震资料,然后将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料。需要说明的是,划分时窗时一定要精细,强反射地层上时窗和强反射地层下时窗内不能包含强反射地层地层,但尽可能包含强反射地层地层上下所有地层。
具体实施时,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。
实施例中,首先根据强反射地层上时窗内的地震资料和时窗长度,确定强反射地层上时窗内的地震资料的第一振幅算数平均值。按如下公式确定强反射地层上时窗内的地震资料的第一振幅算数平均值:
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层上时窗内的地震资料中第i道第j样点的值,Ai为强反射地层上时窗内的地震资料的第一振幅算数平均值,N1为强反射地层上时窗内的时窗长度。
实施例中,确定强反射地层上时窗内的地震资料的第一振幅算数平均值之后,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第二振幅算数平均值。按如下公式确定强反射地层上时窗内的地震资料的第二振幅算数平均值:
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层上时窗内的地震资料中第i道第j样点的值,N1为强反射地层上时窗内的时窗长度,M为地震道个数,A为强反射地层上时窗内的地震资料的第二振幅算数平均值。
实施例中,确定强反射地层上时窗内的地震资料的第二振幅算数平均值之后,根据所述第一振幅算数平均值和第二振幅算数平均值,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。按如下公式确定强反射地层上时窗内的地震资料的第一振幅补偿系数:
其中,i为地震道序号,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,Ai为强反射地层上时窗内的地震资料的第一振幅算数平均值,A为强反射地层上时窗内的地震资料的第二振幅算数平均值。
具体实施时,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。
实施例中,首先根据强反射地层下时窗内的地震资料和时窗长度,确定强反射地层下时窗内的地震资料的第三振幅算数平均值。按如下公式确定强反射地层下时窗内的地震资料的第三振幅算数平均值:
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层下时窗内的地震资料中第i道第j样点的值,Bi为强反射地层下时窗内的地震资料的第三振幅算数平均值,N2为强反射地层下时窗内的时窗长度,N为各地震道中地震资料样点个数。
实施例中,确定强反射地层下时窗内的地震资料的第三振幅算数平均值之后,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第四振幅算数平均值。按如下公式确定强反射地层下时窗内的地震资料的第四振幅算数平均值:
其中,i为地震道序号,j为各地震道中各样点序号,Xij为强反射地层下时窗内的地震资料中第i道第j样点的值,N2为强反射地层下时窗内的时窗长度,N为各地震道中地震资料样点个数,M为地震道个数,B为强反射地层下时窗内的地震资料的第四振幅算数平均值。
实施例中,确定强反射地层下时窗内的地震资料的第四振幅算数平均值之后,根据第三振幅算数平均值和第四振幅算数平均值,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。按如下公式确定强反射地层下时窗内的地震资料的第二振幅补偿系数:
其中,i为地震道序号,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,Bi为强反射地层下时窗内的地震资料的第三振幅算数平均值,B为强反射地层下时窗内的地震资料的第四振幅算数平均值。
具体实施时,根据第一振幅补偿系数和第二振幅补偿系数,确定强反射地区地震资料的第三振幅补偿系数。
实施例中,按如下公式确定强反射地区地震资料的第三振幅补偿系数:
其中,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进一步用于,按如下公式确定强反射地区地震资料的第三振幅补偿系数:
其中,i为地震道序号,j为各地震道中各样点序号,OP3ij为强反射地区地震资料中第i道第j样点的第三振幅补偿系数,OP1i为强反射地层上时窗内的第i道地震资料的第一振幅补偿系数,OP2i为强反射地层下时窗内的第i道地震资料的第二振幅补偿系数,N为各地震道中地震资料样点个数。
一个实施例中,第三振幅补偿系数为振幅补偿因子。
综上所述,本发明实施例通过获取强反射地区地震资料,将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数,根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数,根据第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数,最后根据第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿,有效补偿了由于强反射地层屏蔽作用而衰减的振幅,提高了预测强反射地层下伏奥陶系碳酸盐岩缝洞储层的准确性。本发明实施例通过将强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,并根据强反射地层上时窗和强反射地层下时窗的地震数据计算得到强反射地区地震资料的振幅补偿系数,有效补偿了强反射地区地震资料中由于强反射地层屏蔽作用而衰减的振幅,提高了强反射地层下伏奥陶系碳酸盐岩缝洞储层预测的准确率。
本领域内的技术人员应明白,本发明的实施例可提供为方法、***、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (18)
1.一种强反射地区地震资料振幅补偿方法,其特征在于,包括:
获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料;
根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数;
根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数;
根据所述第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数;
根据所述第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿。
2.如权利要求1所述的方法,其特征在于,所述根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数,包括:
根据强反射地层上时窗内的地震资料和时窗长度,确定强反射地层上时窗内的地震资料的第一振幅算数平均值;
根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第二振幅算数平均值;
根据所述第一振幅算数平均值和第二振幅算数平均值,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。
6.如权利要求1所述的方法,其特征在于,所述根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数,包括:
根据强反射地层下时窗内的地震资料和时窗长度,确定强反射地层下时窗内的地震资料的第三振幅算数平均值;
根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第四振幅算数平均值;
根据所述第三振幅算数平均值和第四振幅算数平均值,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。
11.如权利要求1所述的方法,其特征在于,所述第三振幅补偿系数为振幅补偿因子。
12.一种强反射地区地震资料振幅补偿装置,其特征在于,包括:
时窗划分模块,用于获取强反射地区地震资料,将所述强反射地区地震资料划分为强反射地层时窗、强反射地层上时窗和强反射地层下时窗,其中所述强反射地层时窗内的地震资料是地震波经过强反射地层时产生的地震资料,所述强反射地层上时窗内的地震资料是地震波经过强反射地层上覆地层时产生的地震资料,所述强反射地层下时窗内的地震资料是地震波经过强反射地层下伏地层时产生的地震资料;
上时窗模块,用于根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第一振幅补偿系数;
下时窗模块,用于根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第二振幅补偿系数;
补偿因子模块,用于根据所述第一振幅补偿系数和第二振幅补偿系数,确定所述强反射地区地震资料的第三振幅补偿系数;
振幅补偿模块,用于根据所述第三振幅补偿系数,对强反射地区地震资料衰减的振幅进行补偿。
13.如权利要求12所述的装置,其特征在于,所述上时窗模块进一步用于:
根据强反射地层上时窗内的地震资料和时窗长度,确定强反射地层上时窗内的地震资料的第一振幅算数平均值;
根据强反射地层上时窗内的地震资料、时窗长度和地震道个数,确定强反射地层上时窗内的地震资料的第二振幅算数平均值;
根据所述第一振幅算数平均值和第二振幅算数平均值,确定强反射地层上时窗内的地震资料的第一振幅补偿系数。
14.如权利要求12所述的装置,其特征在于,所述下时窗模块进一步用于:
根据强反射地层下时窗内的地震资料和时窗长度,确定强反射地层下时窗内的地震资料的第三振幅算数平均值;
根据强反射地层下时窗内的地震资料、时窗长度和地震道个数,确定强反射地层下时窗内的地震资料的第四振幅算数平均值;
根据所述第三振幅算数平均值和第四振幅算数平均值,确定强反射地层下时窗内的地震资料的第二振幅补偿系数。
16.如权利要求12所述的装置,其特征在于,所述第三振幅补偿系数为振幅补偿因子。
17.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至11任一所述方法。
18.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至11任一所述方法的计算机程序。
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)
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)
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 |
-
2019
- 2019-01-15 CN CN201910034579.6A patent/CN111435174B/zh active Active
Patent Citations (7)
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)
Title |
---|
周星: "数据驱动下的火成岩能量屏蔽补偿方法", 《 中国石油学会2017年物探技术研讨会论文集》 * |
姚军: "三维叠后地震资料振幅补偿方法在地震沉积学研究中的应用", 《 CPS/SEG北京2018国际地球物理会议暨展览电子论文集》 * |
Cited By (3)
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 |