[发明专利]成像方法、成像装置以及计算机存储介质

权利要求书

1.一种成像方法，其特征在于，包括：

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

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

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

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

从所述多个方向波波场中选择至少两个方向波波场，并对选择出的方向波波场进行归一化互相关矢量成像；其中，所述归一化互相关矢量成像的成像条件包括:

其中，m代表空间位置，θ代表ADCIGs中的离散角度，θ_p代表实际计算的角度，σ代表高斯函数的方差，代表炮点处左下行纵波波场，代表检波点处右上行纵波波场，代表检波点处右上行横波波场，表示炮点处左下行纵波和检波点处右上行纵波的成像结果，表示炮点处左下行纵波和检波点处右上行横波的成像结果。

2.如权利要求1所述的成像方法，其特征在于，所述对分离出的纵波波场和或横波波场在上下左右四个方向进行分解，包括：

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

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

3.如权利要求2所述的成像方法，其特征在于，所述对分离出的纵波波场进行上下方向的方向波分解，包括：

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

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

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

4.如权利要求3所述的成像方法，其特征在于，对上行方向波波场进行左右方向的方向波分解，包括：

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

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

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

5.如权利要求3所述的成像方法，其特征在于，对下行方向波波场进行左右方向的方向波分解，包括：

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

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

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

6.如权利要求2所述的成像方法，其特征在于，所述对分离出的横波波场进行上下方向的方向波分解，包括：

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

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

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