CN105277981A - 基于波场延拓补偿的非一致性时移地震面元匹配方法 - Google Patents

Abstract

本发明涉及地质勘探地中震资料处理方法领域，具体是一种基于波场延拓补偿的非一致性时移地震面元匹配方法。其总体方案是分析第一期资料和第二期资料的面元大小，根据两期面元位置和大小，做折中处理；对面元较小资料进行抽稀处理，使其满足中间面元大小；对面元较大资料进行道集补偿处理，使其匹配到中间面元大小。其中道集补偿处理采用三维波场延拓方法，首先对地震数据进行傅里叶变换，然后在频谱域进行炮集之间的插值，再进行傅里叶反变换得到补偿以后的炮集记录。与常规面元匹配技术相比，基于波场延拓补偿的非一致性时移地震面元匹配方法充分利用两期资料信息，能够得到较好的时移地震分析和识别效果。

Description

基于波场延拓补偿的非一致性时移地震面元匹配方法

技术领域

本发明涉及地质勘探地中震资料处理方法领域，具体是一种基于波场延拓补偿的非一致性时移地震面元匹配方法。

背景技术

随着油气藏的开发，通过对比同一地区不同时期采集的两期地震资料的时移地震技术，得到越来越多的应用。该技术对于检测油气藏的开发状态以及制定下一步滚动勘探具有重要意义。但是，由于两期地震资料采集时的非一致性，造成观测系统上存在差异，包括覆盖次数相差悬殊，面元大小和位置不同等，导致两期地震资料差异不能很好地显示单独由油气藏引起的变化。其中在采集方面，面元的大小对地震资料具有重要影响。

目前应用比较多的是把具有较小面元的第二期资料直接匹配到具有较大面元的第一期资料，这种方法能够较好地匹配数据，消除由于面元变化所引起的油气藏的差异变化，并且也能使能量达到均衡。但是由于第一期资料分辨率有限，单纯地把第二期资料匹配到第一期资料，浪费了大量的第二期资料的信息。另一种应用较多的是把第一期资料通过插值，匹配到第二期资料上，这种方法能够充分利用第二期资料的高分辨率，但是由于第二期资料分辨率比第一期资料高很多，在插值过程中容易形成假异常，不利于时移地震油气藏的分析和识别。

发明内容

本发明的目的是针对目前实际生产中的时移地震绝大多数两期观测系统都不一致，两期地震资料的面元大小、面元位置和覆盖次数等差异会引起两期地震资料差值的异常，直接影响到时移地震的精度的问题，提出一种基于波场延拓补偿的非一致性时移地震面元匹配方法，通过非一致性时移地震面元匹配使两期地震资料的面元大小、面元位置和覆盖次数达到同一个水平，既充分利用了第二期资料的地震信息，又不至于因过度插值出现假油藏差异异常。

本发明的总体技术方案是：

基于波场延拓补偿的非一致性时移地震面元匹配方法是：分析第一期资料和第二期资料的面元大小，根据两期面元位置和大小，做折中处理；对面元较小资料进行抽稀处理，使其满足中间面元大小；对面元较大资料进行道集补偿处理，使其匹配到中间面元大小。

其中道集补偿处理采用三维波场延拓方法，首先对地震数据进行傅里叶变换，然后在频谱域进行炮集之间的插值，再进行傅里叶反变换得到补偿以后的炮集记录。

本发明的具体技术方案是：

第一步：根据油气藏的构造特征和解释数据，建立典型三维模型，并通过测井数据和地震解释结果给出地层速度值；

第二步：设定两期观测系统，第一期为较大面元，覆盖次数较低，第二期为较小面元，覆盖次数较高；

第三步：采用有限差分交错网格声波方程进行正演模拟，得到正演模拟剖面；

第四步：对面元较大的一期资料的炮集数据进行傅里叶变换；

第五步：在频率域采用三维波场延拓方法，根据两炮集的数据以及观测系统中的位置，插值出两炮之间的炮集数据；

第六步：对插值出的炮集数据进行傅里叶反变化，得到时间域的炮集记录；

第七步：对面元较小的二期资料的炮集数据进行抽稀处理，达到与插值后的一期资料相同的面元和覆盖次数；

第八步：比较插值后的一期资料和抽稀后的二期资料，采用剖面相似系数来比较插值出的数据与理论数据之间的误差；

第九步：根据比较结果，优化基于三维波场延拓原理的插值方法。

本发明的技术方案解决了目前实际生产中的时移地震绝大多数两期观测系统都不一致，因两期地震资料的面元大小、面元位置和覆盖次数等差异会引起两期地震资料差值的异常，而影响到时移地震分析和识别精度的问题，通过非一致性时移地震面元匹配技术可以使两期地震资料的面元大小、面元位置和覆盖次数达到同一个水平，既充分利用了第二期资料的地震信息，又不至于因过度插值出现假油藏差异异常。与常规面元匹配技术相比，基于波场延拓补偿的非一致性时移地震面元匹配方法充分利用两期资料信息，能够得到较好的时移地震分析和识别效果。

附图说明

图1一种基于波场延拓补偿的非一致性时移地震面元匹配方法实施例的流程图

图2本发明一种实施例的第一期4炮地震记录

图3本发明一种实施例的面元匹配后炮点位置

图4本发明一种实施例的第二期25炮地震记录

图5本发明一种实施例的第2炮相似系数曲线

图6本发明一种实施例的第4炮相似系数曲线

图7本发明一种实施例的第5炮相似系数曲线

图8本发明一种实施例的第6炮相似系数曲线

图9本发明一种实施例的第8炮相似系数曲线。

具体实施方式

分析第一期资料和第二期资料的面元大小，根据两者面元位置和大小，做折中处理。对第二期资料进行抽稀处理，使其满足中间面元大小；对第一期资料进行道集补偿处理，使其匹配到中间面元大小。这里补偿技术采用三维波场延拓方法，首先对地震数据进行傅里叶变换，然后在频谱域进行炮集之间的插值，再进行傅里叶反变换得到补偿以后的炮集记录。

结合附图1，具体实现方式如下：

第一步：根据油气藏的构造特征和解释数据，建立典型三维模型，并通过测井数据和地震解释结果给出地层速度值；

第二步：设计两期观测系统，第一期为较大面元，覆盖次数较低，第二次为较小面元，覆盖次数较高；

第三步：采用有限差分交错网格声波方程进行正演模拟，得到正演模拟剖面；

第四步：对面元较大的一期资料的炮集数据进行傅里叶变换；

第五步：在频率域采用三维波场延拓原理，根据两炮集的数据以及观测系统中的位置，插值出两炮之间的炮集数据；

第六步：对插值出的炮集数据进行傅里叶反变化，得到时间域的炮集记录；

第七步：对面元较小的二期资料的炮集数据进行抽稀处理，达到与插值后的一期资料相同的面元和覆盖次数；

第八步：比较插值后的一期资料和抽稀后的二期资料，采用剖面相似系数来比较插值出的数据与理论数据之间的误差；

第九步：根据比较结果，优化基于三维波场延拓原理的插值方法。

结合附图2-9，通过一个应用实施例进一步说明。

以一个正方形研究工区为例，第一期资料共有4炮（图2所示），面元大小为2.5m*100m，第二期资料共有25炮（图4所示），面元大小为2.5m*25m。根据前面的方案，这里对第一期资料进行补偿，每两个炮之间插出一炮的数据，这样第一期资料就变成了9炮（图3所示），面元2.5m*50m，其中第2、4、5、6、8炮为补偿出来的炮集。

补偿后的新资料可以与第二期资料中相对应的9炮数据进行相同面元下的数据对比和分析。

如图5-9所示，为了展示面元匹配的效果，这里给出补偿炮与第二期原始炮的相似系数曲线。

通过以上实施例和附图所示进一步表明，与常规面元匹配技术相比，本发明的基于波场延拓补偿的非一致性时移地震面元匹配方法充分利用两期资料信息，能够得到较好的时移地震效果。

Claims (3)

1.基于波场延拓补偿的非一致性时移地震面元匹配方法，其特征是：分析第一期资料和第二期资料的面元大小，根据两期面元位置和大小，做折中处理；对面元较小资料进行抽稀处理，使其满足中间面元大小；对面元较大资料进行道集补偿处理，使其匹配到中间面元大小。

2.根据权利要求1所述的基于波场延拓补偿的非一致性时移地震面元匹配方法，其特征是：道集补偿处理采用三维波场延拓方法，首先对地震数据进行傅里叶变换，然后在频谱域进行炮集之间的插值，再进行傅里叶反变换得到补偿以后的炮集记录。

3.根据权利要求1或2所述的基于波场延拓补偿的非一致性时移地震面元匹配方法，其特征是：

第一步：根据油气藏的构造特征和解释数据，建立典型三维模型，并通过测井数据和地震解释结果给出地层速度值；

第二步：设定两期观测系统，第一期为较大面元，覆盖次数较低，第二期为较小面元，覆盖次数较高；

第三步：采用有限差分交错网格声波方程进行正演模拟，得到正演模拟剖面；

第四步：对面元较大的一期资料的炮集数据进行傅里叶变换；

第五步：在频率域采用三维波场延拓方法，根据两炮集的数据以及观测系统中的位置，插值出两炮之间的炮集数据；

第六步：对插值出的炮集数据进行傅里叶反变化，得到时间域的炮集记录；

第七步：对面元较小的二期资料的炮集数据进行抽稀处理，达到与插值后的一期资料相同的面元和覆盖次数；

第八步：比较插值后的一期资料和抽稀后的二期资料，采用剖面相似系数来比较插值出的数据与理论数据之间的误差；

第九步：根据比较结果，优化基于三维波场延拓原理的插值方法。