[发明专利]确定地震波场的坡印廷矢量的方法、装置以及计算机存储介质

权利要求书

1.一种确定地震波场的坡印廷矢量的方法，其特征在于，包括：

构建地下炮点正传波场以及检波点反传波场；

判断传播所述正传波场以及所述反传波场的介质类型；

当所述介质为各向同性介质时，根据预设的方程解耦算法，对所述正传波场或所述反传波场进行纵横波场分离，以及当所述介质为非均匀各向异性时，根据预设的变异函数IDW插值算法，对所述正传波场或反传波场进行纵横波场分离；

对分离出的纵波波场或横波波场在上下左右四个方向进行分解，得到所述纵波波场或所述横波波场的多个方向波波场；

从所述多个方向波波场中选择出多个方向波波场，并计算选择出的所述多个方向波波场的坡印廷矢量；其中，

所述当所述介质为非均匀各向异性时，根据预设的变异函数IDW插值算法，对所述正传波场或反传波场进行纵横波场分离，包括：

当在非均匀各向异性介质下时，利用傅里叶变换分别将各时刻的所述正传波场或反传波场由空间域变换至波数域；

选取N个参考模型，根据所述参考模型计算分离算子，并在波数域利用自褶积组合窗函数对所述分离算子进行截断优化；

在波数域利用所述分离算子对所述参考模型下的所述正传波场或反传波场进行分离，得到所述参考模型下的分离后的所述纵横波波场；

将所述参考模型下的分离后的纵横波波场由波数域变换至空间域，N个所述参考模型对应得到N个所述参考模型下的分离后的纵横波波场，N为正整数；

在空间域，利用变异函数改进IDW插值算法并根据改进后的所述IDW插值算法计算N个参考模型的权重系数，对N个所述参考模型下的分离后的纵横波波场进行加权插值处理，得到最终分离后的纵波波场和横波波场。

2.如权利要求1所述的确定地震波场的坡印廷矢量的方法，其特征在于，所述对分离出的纵波波场或横波波场在上下左右四个方向进行分解包括：

对分离出的纵波波场和横波波场进行上下方向的方向波分离，得到对应的上行方向波波场和下行方向波波场；

对上行方向波波场和下行方向波波场进行左右方向的方向波分离，得到对应的左上行、左下行方向波波场，和右上行、右下行方向波波场。

3.如权利要求2所述的确定地震波场的坡印廷矢量的方法，其特征在于，所述对分离出的纵波波场进行上下方向的方向波分解，包括：

将分离出的纵波波场进行时间复数域拓展；

在Z方向对时间复数域拓展之后的纵波波场进行傅里叶变换，得到波数域纵波波场；

从所述波数域纵波波场中选择波数大于零的部分，波数小于零的部分置为零，进行Z方向上的反傅里叶变换，得到纵波波场的下行方向波波场；并从所述波数域纵波波场中选择波数小于零的部分，波数大于零的部分置为零，进行Z方向上的反傅里叶变换，得到所述纵波波场的上行方向波波场。

4.如权利要求3所述的确定地震波场的坡印廷矢量的方法，其特征在于，所述对上行方向波波场进行左右方向的方向波分解，包括：

将所述纵波波场的上行方向波波场进行时间复数域拓展；

在X方向对时间复数域拓展之后的所述纵波波场的上行方向波波场进行傅里叶变换，得到波数域纵波波场；

从所述波数域纵波波场的上行方向波波场中选择波数大于零的部分，波数小于零的部分置为零，进行X方向上的反傅里叶变换，得到所述纵波波场的右上行方向波波场；并从所述波数域纵波波场中选择波数小于零的部分，波数大于零的部分置为零，进行X方向上的反傅里叶变换，得到所述纵波波场的左上行方向波波场。

5.如权利要求3所述的确定地震波场的坡印廷矢量的方法，其特征在于，对下行方向波波场进行左右方向的方向波分解，包括：

将所述纵波波场的下行方向波波场进行时间复数域拓展；

在X方向对时间复数域拓展之后的所述纵波波场的下行方向波波场进行傅里叶变换，得到波数域纵波波场；

从所述波数域纵波波场的下行方向波波场中选择波数大于零的部分，波数小于零的部分置为零，进行X方向下的反傅里叶变换，得到所述纵波波场的右下行方向波波场；并从所述波数域纵波波场中选择波数小于零的部分，波数大于零的部分置为零，进行X方向下的反傅里叶变换，得到所述纵波波场的左下行方向波波场。