[发明专利]一种对地震波沿传播路径进行Q补偿的偏移成像方法及装置

权利要求书

1.一种对地震波沿传播路径进行Q补偿的偏移成像方法，其特征在于，所述方法包括：

获取目标地层倾角信息；

根据所述目标地层倾角信息，确定菲涅尔带的边界；

基于确定的所述菲涅尔带的边界，在S域内对所述菲涅尔带内的PP波或PS波沿传播路径进行Q补偿，得到补偿后的波场值；

根据所述菲涅尔带的大小优化偏移孔径，并基于优化后的所述偏移孔径，确定Q补偿后的波场值在时间偏移剖面上的位置；

所述获取目标地层倾角信息，包括：

获取Kirchhoff叠前时间偏移剖面；

在所述Kirchhoff叠前时间偏移剖面上拾取目标地震层位；

基于所述目标地震层位，建立多个层位格架；

根据公式1：确定与各所述层位格架对应的倾角，其中，ΔX为与所述层位格架对应的地层界面上相邻两点之间的水平距离，ΔZ为与所述层位格架对应的地层界面上相邻两点之间的垂直距离，为与所述层位格架对应的地层界面的倾角，V_ave表示所述层位格架之间的地震波平均速度，ΔT表示时间剖面上所述相邻两点之间的垂向时间间隔；

利用Kriging插值公式在相邻两个所述层位格架间进行插值得到倾角模型；

基于倾角模型，以及目标地层信息，获取目标地层倾角信息；

在基于确定的所述菲涅尔带的边界，在S域内对所述菲涅尔带内的PP波或PS波沿传播路径进行Q补偿，得到补偿后的波场值之前，所述方法还包括：

建立目标等效Q模型；

所述建立目标等效Q模型，包括：

建立常等效Q模型；

根据所述倾角模型和所述常等效Q模型，进行Kirchhoff叠前时间偏移；

将使用所述倾角模型和所述常等效Q模型获得的成像道，与使用所述倾角模型获得到的成像道相比较，判断补偿后成像道的深浅子波是否一致；

若否，则通过自上往下逐步增大的Q扫描方式对所述常等效Q模型进行更新操作，得到等效Q模型；

利用更新后的所述等效Q模型替代所述常等效Q模型，直至补偿后的成像道深浅子波一致，则输出更新的等效Q模型，并将所述等效Q模型确定为目标等效Q模型；

若是，则输出所述常等效Q模型，并将所述常等效Q模型确定为目标等效Q模型；

所述目标等效Q模型包括：Q_eff,P模型和Q_eff,PS模型，所述Q_eff,P模型为PP波的等效Q模型，所述Q_eff,PS模型为PS波的等效Q模型，所述方法还包括：

在建立目标等效Q模型时，依次建立Q_eff,P模型和Q_eff,PS模型；

所述基于确定的所述菲涅尔带的边界，在S域内对所述菲涅尔带内的PP波或PS波沿传播路径进行Q补偿，得到补偿后的波场值，包括：

根据公式6：确定P波在成像点处的等效Q值，其中，Q_eff,P表示P波在成像点处的等效Q值，Q_n,P表示P波在成像点处的层Q值，ΔZ_n表示每一层的厚度，V_n,P表示P波的层速度；

在S域内，根据公式7：对所述菲涅尔带内的PP波沿传播路径进行Q补偿，得到补偿后的波场值，其中，t_PP表示PP波在黏弹性介质中的传播时间；f_PP(x,y,ω,τ)表示菲涅尔带记录的PP波，经过t_PP时间传播衰减后，到达水平地表(x,y)处在S域内的地震响应；F_PP(x,y,t_PP)表示对f_PP(x,y,ω,τ)在S域内沿传播路径对所述PP波波场进行Q补偿后的波场值；ω表示频率；

根据公式8：确定PS波在成像点处的等效Q值，其中，Q_eff,PS表示PS波在成像点处的等效Q值，ΔZ_n表示每一层的厚度，Q_n,S表示S波的层Q值，V_n,S表示PS波的层速度；

在S域内，根据公式9：对所述菲涅尔带内的PS波沿传播路径进行Q补偿，得到补偿后的波场值，其中，t_PS表示PS波在黏弹性介质中的传播时间，f_PS(x,y,ω,τ)表示菲涅尔带内的PS波经过t_PS时间传播衰减后，到达水平地表(x,y)处在S域内的地震响应，F_PS(x,y,t_PS)表示对f_PS(x,y,ω,τ)在S域内沿传播路径对所述PS波波场进行Q补偿后的波场值。