[发明专利]弹性波矢量成像方法、装置、存储介质及设备

权利要求书

1.一种弹性波矢量成像方法，其特征在于，包括：

利用成像区域中各点各向异性介质的设定方向的P波速度、设定方向的S波速度以及各向异性参数，建立群角与极化角之间的转换关系；

应用弹性波方程在成像区域进行多分量地震数据的波场深度延拓，得到震源正传波场和炮记录反传波场，以及相应的坡印廷矢量；

按设定时刻间距，利用坡印廷矢量和群角与极化角之间的转换关系，对震源正传波场和炮记录反传波场进行波场分解；

利用多分量地震数据中所有炮对应的波场分解结果进行各向异性介质中弹性波矢量成像；

其中，所述按设定时刻间距，利用坡印廷矢量和群角与极化角之间的转换关系，对震源正传波场和炮记录反传波场进行波场分解，包括：

分别针对P波和SV波，利用设定时刻间距的坡印廷矢量计算得到相应的群角值；

分别针对P波和SV波，根据群角值和群角与极化角之间的转换关系，获取成像区域中各点的极化角角度；

分别针对P波和SV波，将极化角角度代入波场分解公式，计算得到伪P波和伪S波，作为波场分解结果。

2.如权利要求1所述的弹性波矢量成像方法，其特征在于，利用成像区域中各点各向异性介质的设定方向的P波速度、设定方向的S波速度以及各向异性参数，建立群角与极化角之间的转换关系，包括：

对Christoff方程求特征值，并结合设定方向的P波速度、设定方向的S波速度以及各向异性参数，分别针对P波和SV波，计算得到以相角为自变量的相速度，其中，设定方向包括对称横方向和沿对称轴方向；

分别针对P波和SV波，将相速度和相速度对相角的偏导数代入群角与相角之间的关系式，计算得到群角与相角之间的转换关系；

分别针对P波和SV波，在Christoff方程中求解极化角，并结合各向异性参数、相速度及群角与相角之间的转换关系，计算得到群角与极化角之间的转换关系。

3.如权利要求2所述的弹性波矢量成像方法，其特征在于，分别针对P波和SV波，在Christoff方程中求解极化角，并结合各向异性参数、相速度及群角与相角之间的转换关系，计算得到群角与极化角之间的转换关系，包括：

分别针对P波和SV波，在设定角度范围内按设定角度间隔取群角角度值；

分别针对P波和SV波，在Christoff方程中求解极化角，并结合各向异性参数、相速度、群角角度值及群角与相角之间的转换关系，计算得到群角和极化角之间的数值对应表，作为群角与极化角之间的转换关系。

4.如权利要求1所述的弹性波矢量成像方法，其特征在于，应用弹性波方程在成像区域进行多分量地震数据的波场深度延拓，得到震源正传波场和炮记录反传波场，以及相应的坡印廷矢量，包括：

对多分量地震数据进行预处理，得到多炮可用于逆时偏移的炮记录；

由每炮的炮记录的频宽生成宽频子波，并利用宽频子波进行数值模拟，生成成像区域的边界记录；

针对每炮的边界记录和炮记录，利用弹性波方程数值解法进行弹性波场反传播，并在不同时刻将边界记录赋值给成像区域的边界位置，得到不同时刻的震源正传波场及相应的坡印廷矢量，在不同时刻将炮记录赋值给检波点位置，得到不同时刻的炮记录反传波场及相应的坡印廷矢量。

5.如权利要求1所述的弹性波矢量成像方法，其特征在于，分别针对P波和SV波，根据群角值和群角与极化角之间的转换关系，获取极化角角度，包括：

分别针对P波和SV波，利用群角值查找群角与极化角之间的转换关系的数值对应表，得到相应的极化角值；

分别针对P波和SV波，利用得到的极化角值进行插值，得到成像区域中各点的极化角角度。