[发明专利]一种叠前道集的全时间域相位一致性校正方法

说明书

技术领域

本发明属于石油物探领域，涉及一种叠前道集的全时间域相位一致性校正方法。

背景技术

提高地震剖面的信噪比和分辨率始终是地震处理解释研究的重要内容，在动静校正后的CRP道集内振幅和相位相同的假设条件下，且消除了剩余时差，地震叠加技术能明显提高资料信噪比。地震资料处理中通常用反褶积的方法来消除子波影响，其中应用范围最广的最小相位反褶积假设地震波是最小相位的，这种假设往往得不到满足时会产生预料不到的剩余相位。

人们提出了很多关于相位校正的方法，如在时间域和频率域都可以实现的常相位校正方法、变相位校正方法、分频常相位校正、时差常相位校正及加权叠加、地表一致性相位校正等。其中与本发明最相似的是变相位校正方法，变相位校正的基本思想是将整个地震记录按不同的属性（频率、时间、空间）划分为多个小窗口，然后假定每个小窗口中的相位校正量为常数，所使用的校正方法与常相位校正方法相同。

上述的相位校正方法，是在时间域同相轴消除剩余时差的基础上进行的，而已有的消除剩余时差方法对地震资料信噪比有较高的要求，在搜索邻近地震道的时候遇到断裂情况或者信号弱的地方容易失去方向。

发明内容

为克服现有的相位校正方法存在的技术缺陷，本发明公开了一种叠前道集的全时间域相位一致性校正方法。

本发明所述的一种叠前道集的全时间域相位一致性校正方法，包括如下步骤：

    S1.设定时间窗口，选择叠前道集D1中的多个道叠加形成第一模板道，所述叠加为将各道相同时间点的振幅值取平均作为第一模板道的振幅；所述时间窗口包含至少一个完整地震波形；

    S2.将第一模板道分别与叠前道集各道做互相关，找出相关度最大值位置,求取出时移量DT1；遍历叠前道集内的每一道及每一道在时间窗口内的全部时间点；

S3.利用时移量DT1将叠前道集D1拉平后得到D2；

    S4. 从D2中选择至少包含两个完整地震波性的时间窗口WT，并从D2中选择多个道叠加形成第二模板道，将时间窗口WT内全部道和第二模板道进行希尔伯特变换后，对WT中的每一道，利用下式计算相似度S；

                                                    

其中，

式中H表示希尔伯特变换，表示时间窗口WT中的道，为第二模板道进行希尔伯特变换后得到的包络，时间t1、t2表示预先定义的时间窗口WT两端节点；

S5. 对WT中每一道，选择相似度S最大的值对应的角度 作为相位修正量，移动时间窗口WT，重复步骤S4，直至遍历D2中的全部道。

优选的，所述步骤S1中选择叠前道集中的近道和/或中道叠加。

优选的，所述步骤S1中用于叠加的道占整个叠前道集数量的20%-50%。

 优选的，所述步骤S4中选择叠前道集D2中的近道和/或中道叠加。

本发明中的非物理性道集拉平消除剩余时差和分段常相位校正串行使用的方法实现CRP道集的相位一致性校正使相位校正对资料的要求降低，效率增高。优选利用近、中道叠加道做模板，可以减少远道动校拉伸对第一模板道的影响。

附图说明

图1为叠前道集的一个图形化显示形式图，

图2为本发明的一种具体实施方式示意图。

具体实施方式

下面结合附图，对本发明的具体实施方式作进一步的详细说明。

  本发明所述的一种叠前道集的全时间域相位一致性校正方法包含如下步骤，S1.设定时间窗口，选择叠前道集D1中的多个道叠加形成第一模板道，所述叠加为将各道相同时间点的振幅值取平均作为第一模板道的振幅；所述时间窗口包含至少一个完整地震波形；

    S2.将第一模板道分别与叠前道集各道做互相关，找出相关度最大值位置,求取出时移量DT1；遍历叠前道集内的每一道及每一道在时间窗口内的全部时间点；

S3.利用时移量DT1将叠前道集D1拉平后得到D2 ；

    S4.从D2中选择至少包含两个完整地震波性的时间窗口WT，并从WT中选择多个道叠加形成第二模板道，将时间窗口WT内全部道和第二模板道进行希尔伯特变换后，对WT中的每一道，利用下式计算相似度S；

     

其中

式中H表示希尔伯特变换，表示WT中的道,为第二模板道进行希尔伯特变换后得到的包络，时间t1、t2表示预先定义的时间窗口WT两端节点；