[发明专利]弹性波场的方向波分解方法、装置以及计算机存储介质

说明书

本申请公开了一种弹性波场的方向波分解方法、装置以及计算机存储介质。该方法包括：对输入的波场进行基于希尔伯特变换的时间复数域拓展，得到时间复数域拓展波场；并对时间复数域拓展波场进行弹性波数值模拟；基于预设的变异函数的IDW插值算法，对弹性波数值模拟后的时间复数域拓展波场进行纵横波解耦分离，得到纵波时间复数域波场及横波时间复数域波场；基于希尔伯特变换的空间复数域拓展，对纵波时间复数域波场及横波时间复数域波场进行前后左右上下6个方向的方向波分解，得到多个方向波波场。本申请能对三维弹性波场进行方向波分解，提高计算效率，能利用分解后的各方向波计算坡印廷矢量，进而获取更准确的角度域成像道集，提高成像质量。

技术领域

本申请涉及地震勘探领域，尤其是涉及一种弹性波场的方向波分解方法、装置以及计算机存储介质。

背景技术

弹性波逆时偏移中成像噪音包括耦合噪音，低频噪音和成像假象噪音。耦合噪音主要来自于纵波和横波耦合在一起成像造成的噪音。低频噪音和成像假象噪音主要来源于特定方向的方向波互相关成像，比如，I_ldld,I_lulu,I_rdrd,I_ruru...这些相同下标的组合是产生成像低频噪声的主要来源；又比如，I_luld,I_lurd,I_ruld,I_rurd...这些炮点波场上行波与检波点波场下行波，当存在强速度变化时，由于后向散射的存在，会导致一次波假象，但该类组合在成盐丘底部构造的像时，又是准确的。

另外，如果地下反射层倾角较大，采用炮点左行和检波点右行的组合，可以取得较好的结果。因此，为了压制低频噪音和成像假象噪音，需要对成像前的波场进行方向的分解。Fei等(2010)提出可以将上行和下行波场分离，然后对不同传播方向的波场分别应用成像条件，组合成像。

Fei等(2015)提出成像假象的概念，他们认为Liu(2011)的做法虽然能够较好地压制低频噪音，但是还是忽略了当速度变化剧烈时，后向散射造成的成像假象，即由炮点前向波场和检波点后向波场互相关造成的假象I_ud，称之为一次波假象(primary false image)，并由此提出了de-primary RTM算法。王一博等(2016)在Fei等(2015)工作的基础上，推导出了标量波方程上下左右4个象限方向波的分离。

现有技术中，上下左右方向波分解可以较好地压制由同方向传播的方向波相关产生的低频噪音，并且一定程度上压制由炮点前向波场和检波点后向波场互相关造成的假象I_ud噪音，但对于三维情况，成像噪音的产生的原因更加复杂，只进行上下左右方向波的分解，不能满足达到预设的成像质量的要求。

发明内容

本发明的目的在于提供了一种弹性波场的方向波分解方法、装置以及计算机存储介质，以解决上述的问题。

为了实现上述目的，本发明提供了一种弹性波场的方向波分解方法、装置以及计算机存储介质，该弹性波场的方向波分解方法包括：

对输入的波场进行基于希尔伯特变换的时间复数域拓展，得到时间复数域拓展波场；

对所述时间复数域拓展波场进行弹性波数值模拟；

基于预设的变异函数的IDW插值算法，对所述弹性波数值模拟之后的时间复数域拓展波场进行纵横波解耦分离，得到纵波时间复数域波场以及横波时间复数域波场；

基于所述希尔伯特变换的空间复数域拓展，对所述纵波时间复数域波场和所述横波时间复数域波场进行前后左右上下6个方向的方向波分解，得到多个方向波波场。

在一可选的实施方式中，所述基于希尔伯特变换的空间复数域拓展，对所述纵波时间复数域波场和所述横波时间复数域波场进行前后左右上下6个方向的方向波分解，得到多个方向波波场，包括：