[发明专利]一种碳酸盐岩储层的逆时偏移方法

摘要

本发明公开了一种碳酸盐岩储层的逆时偏移方法，该方法采用针对碳酸盐岩的缝洞孔隙构造常伴随有由于构造运动所形成的断裂和大角度的形态扭曲，运用逆时偏移方法对其进行精确刻画，本发明有限差分方法求解波动方程进行地震波场模拟，波场传播至边界处的处理方法采用PML吸收边界条件，震源波场恢复部分采用逆PML波场恢复技术，计算部分采用GPU并行加速技术。本发明采用PML吸收边界条件，比随机边界条件更能够有效的消除边界反射干扰，采用逆PML波场恢复技术，相比于传统的波场存储策略，在不增加计算量的前提下，大幅减少了存储量的需求，而且波场恢复效果良好，无振幅衰减，采用GPU并行加速技术，提高了运算效率。

主权项

一种碳酸盐岩储层的逆时偏移方法，其特征在于，所述的碳酸盐岩储层的逆时偏移方法采用有限差分方法求解波动方程进行波场模拟，波场传播至边界处的处理方法采用PML吸收边界条件，震源波场恢复部分采用逆PML波场恢复技术，计算部分采用GPU并行加速技术；该方法采用有限差分方法求解波动方程进行波场模拟，具体方法为：用U表示某一时刻t二维空间任一点(x，z)处的位移，二维声波方程可表示如下：∂2U∂x2+∂2U∂z2=1v2∂2U∂t2]]>式中，U为声波波场，v为速度场，时间导数采用二阶中心差分、空间导数为2N阶差分精度的二维声波方程的高阶有限差分格式为：Un+1=(x,z)=2Un(x,z)-Un-1(x,z)+dt2v2{Cx2[Un(x,z)]+Cz2[Un(x,z)]}]]>式中，Cx2[Un(x,z)]=Σl=-NNalΔx2U(x+l,z)]]>Cz2[Un(x,z)]=Σl=-NNalΔz2U(x,z+l)]]>其中，al为有限差分系数，即al=(-1)l+1l2(-1)l+1Πi=1,i≠lNi2Πi=1l-1(l2-i2)Πi=l+1N(i2-l2).]]>