[发明专利]一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法

权利要求书

1.一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法，其特征在于，具体包括以下步骤：

步骤一：输入纵横波速度场、多分量实际观测炮记录、起伏地表高程和观测系统文件；

步骤二：根据起伏地表高程生成正交贴体网格，并将纵横波速度场变换到曲坐标系；

步骤三：计算起伏地表条件下正向延拓的棱柱波震源P波波场和S波波场；曲坐标系纵横波解构方程的棱柱波FI波场正向延拓的计算式如式(1)所示：

其中，L₀表示震源波场正向延拓算子，和分别为入射P波震源波场和入射S波震源波场，F表示震源函数；

曲坐标系纵横波解构方程的棱柱波IF波场正向延拓的计算式如式(2)所示：

其中，L₁表示震源波场到达沉积层后向上反射的波场延拓算子，和分别为反射P波震源波场和反射S波震源波场；

步骤四：计算起伏地表条件下反向延拓的棱柱波检波点P波波场和S波波场；曲坐标系纵横波解构方程的棱柱波FI波场逆时延拓的计算式如式(3)所示：

其中，表示检波点波场到达沉积层后向上反射的波场延拓算子，和分别是检波点反射P波波场和反射S波波场，d_obs为输入的观测数据；

曲坐标系纵横波解构方程的棱柱波IF波场逆时延拓的计算式如式(4)所示：

其中，表示检波点波场的逆时延拓算子，和分别是检波点入射P波波场和入射S波波场

步骤五：求取曲坐标系下棱柱波多分量成像结果；

步骤六：将曲坐标系下的棱柱波多分量成像结果变换到笛卡尔坐标系下；

步骤七：输出纵横波解构方程的弹性棱柱波逆时偏移成像。

2.如权利要求1所述的一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法，其特征在于，所述步骤五具体为：

棱柱波FI波场互相关成像条件为式(5)所示：

其中，I_p1(x)和I_s1(x)分别为棱柱波FI纵波波场和横波波场的成像结果；

棱柱波IF波场互相关成像条件为式(6)所示：

其中，I_p2(x)和I_s2(x)分别为棱柱波IF纵波波场和横波波场的成像结果，x为空间坐标，t为时间，T为总计算时间长度；

最后，多分量棱柱波的成像结果为棱柱波IF波场和棱柱波FI波场的成像结果之和为式(7)所示：

其中，I_p(x)为棱柱波纵波波场成像结果，I_s(x)为棱柱波横波波场成像结果。