[发明专利]一种保持低频信息的谱延拓方法

说明书

本发明提供了一种保持低频信息的谱延拓方法。所述方法采用数据分层、频率分段、子波时变延拓的思路；求取分层数据振幅谱，再求取分段频率振幅包络拟合值，再求取输入子波补偿因子，能量补偿后的输入子波变换为时间域子波；同理，利用解释合成地震记录数据，求取希望中、高频谱延拓的输出子波并变换到时间域，最后利用输入、希望子波进行反褶积得到分辨率提高的剖面，并对每层子波振幅谱进行过度衰减处理，得到保持低频信息的高分辨率资料，为后续岩性解释提供丰富低频信息的高品质资料。

技术领域

本发明涉及石油勘探地震资料处理领域，更具体地讲，涉及一种保持低频信息的谱延拓方法。

背景技术

高分辨率处理技术一直是岩性油气藏勘探的关键技术，随着众多地球物理学家对低频地震异常与油气藏密切联系的逐渐重视，因为在山地地震勘探资料中，反射波低频部分的变化和异常记载了反映地下岩石孔隙、流体等丰富信息，保持低频信息的潜在意义和技术应用价值越来越受到国内外地球物理人员的重视，所以如何在高分辨率处理中保留更丰富的低频信息是目前进行高分辨率处理技术面临的难题。

常规高分辨率处理技术大都以求取地震资料反褶积算子，通过算子来对剖面能量的再分配，往往是以损失部分低频信息来换回高频信息的拓展，无法做到二者兼顾，从而给后续的岩性反演解释带来一定的困难。

发明内容

针对现有技术中存在的不足，本发明的目的之一在于解决上述现有技术中存在的一个或多个问题。例如，本发明能够获得保持低频信息的高分辨率资料的方法。

为了实现上述目的，本发明的提供了一种保持低频信息的谱延拓方法，所述方法包括以下步骤：A、对叠后地震数据进行分层；B、确定第j层地震数据振幅谱的最低频率和最高频率，并根据第j层地震数据振幅谱的最低频率和最高频率求取俞氏子波频率域下的第二振幅谱；C、合成单道合成记录数据并转换到频率域下，得到第三振幅谱；D、将第j层地震数据的第一振幅谱、第二振幅谱和第三振幅谱在其各自最低频率和最高频率之间分成n个频率段的振幅数据和E、利用曲线拟合计算每个频率段的第一振幅的振幅包络值第二振幅的振幅包络值和第三振幅的振幅包络值F、利用第一振幅的振幅包络值和第二振幅的振幅包络值求取补偿因子a，利用补偿因子a对第二振幅谱进行能量补偿，将能量补偿后的子波变换为时间域子波G、利用第二振幅的振幅包络值和第三振幅的振幅包络值求取中频段和高频段的补偿因子b，利用补偿因子b对能量补偿后的第二振幅谱进行中频段和高频段能量补偿，然后将进行中频段和高频段能量补偿后的子波变换为时间域子波H、利用子波与第j层地震数据反褶积得反射系数，利用反射系数与子波褶积得保持低频信息的高分辨率剖面；I、重复所述步骤B至H，直到每层地震数据均得到保持低频信息的高分辨率剖面；其中，i＝1,…,n；j＝1,…,N，N为叠后地震数据分层后的层数。

根据本发明保持低频信息的谱延拓方法的一个实施例，在所述步骤A中，对叠后地震数据按浅层、中层和深层进行分层。

根据本发明保持低频信息的谱延拓方法的一个实施例，所述步骤B包括根据俞氏子波公式(1)求取俞氏子波频率域下的第二振幅谱，

其中，f_l为第j层地震数据振幅谱的最低频率，f_h为第j层地震数据振幅谱的最高频率。

根据本发明保持低频信息的谱延拓方法的一个实施例，所述步骤C利用测井声波曲线合成并标定的单道测井反射系数曲线，利用勘探目的层所需主频确定子波主频，与反射系数褶积合成单道合成记录数据，再将合成记录数据变换到频率域，得到第三振幅谱，其中，目的层用户可自己定义，可把当前处理层作为目的层，比如第j层。

根据本发明保持低频信息的谱延拓方法的一个实施例，所述步骤D通过设置n个窄带滤波器将第j层数据分成n个频率段的振幅数据。

根据本发明保持低频信息的谱延拓方法的一个实施例，所述步骤E包括：