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

说明书

本发明涉公开了一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法，具体涉及石油地球物理勘探领域。该方法具体包括输入纵横波速度场、多分量实际观测炮记录、起伏地表高程和观测系统文件；根据起伏地表高程生成正交贴体网格，并将纵横波速度场变换到曲坐标系；计算起伏地表条件下正向延拓的棱柱波震源P波波场和S波波场；计算起伏地表条件下反向延拓的棱柱波检波点P波波场和S波波场；求取曲坐标系下棱柱波多分量成像结果；将曲坐标系下的棱柱波多分量成像结果变换到笛卡尔坐标系下；输出纵横波解构方程的弹性棱柱波逆时偏移成像。

技术领域

本发明涉及石油地球物理勘探领域，具体涉及一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法。

背景技术

剧烈的山前带起伏地表给地震成像带来巨大困难。常规地震勘探应用反射波进行成像，但是无法对高陡构造进行成像，因此，利用棱柱波信息可对高陡构造进行成像，常规棱柱波成像技术都是基于声波近似假设对单分量地震记录进行成像的，无法充分利用横波信息。

发明内容

本发明的目的是针对上述不足，提出了一种在校正起伏地表影响的前提下，充分利用分离的纵横波信息，提高高陡构造的成像精度，压制纵横波串扰噪音的基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法。

本发明具体采用如下技术方案：

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

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

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

步骤三：计算起伏地表条件下正向延拓的棱柱波震源P波波场和S波波场；

步骤四：计算起伏地表条件下反向延拓的棱柱波检波点P波波场和S波波场；

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

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

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

优选地，所述步骤三具体为：

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

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

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

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

优选地，所述步骤四具体为：

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

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

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

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

优选地，所述步骤五具体为：