[发明专利]三分量地震波积分偏移方法和装置

权利要求书

1.一种三分量地震波积分偏移方法，其特征在于，包括：

计算三维速度模型空间中当前网格点的积分加权因子；

根据当前网格点的位置，在走时表中获取地震波从炮点到检波点之间的全走时；

根据三分量地震波的振幅，获取该网格点到检波器的地震波出射方向矢量；

对所述方向矢量和三分量地震波的振幅矢量做点乘运算得到一个标量振幅；

根据该网格点的积分加权因子对该网格点对应的标量振幅做偏移的加权运算，并将加权运算的结果与该网格点原存储的数据进行积分求和，得到该网格点的偏移图像。

2.如权利要求1所述的方法，其特征在于，计算三维速度模型空间中当前网格点的积分加权因子，包括：

对所述三维速度模型空间中的网格点，由孔距参数计算出相对此网格点的积分加权因子。

3.如权利要求2所述的方法，其特征在于，所述孔距参数按照以下方式获取：

获取共炮点反射波地震数据，其中，所述地震数据为三分量地震数据；

根据所述共炮点反射波地震数据的炮点和检波点之间的几何位置，计算所述共炮点反射波地震数据的孔距参数。

4.如权利要求2所述的方法，其特征在于，由孔距参数计算出相对此网格点的积分加权因子，包括：

按照以下公式计算所述积分加权因子：

W1=zrg+rs(rgrs)3/2]]>

W2=(1-α^2β^2)(1+α^2)1/2(1+β^2)1/2]]>

W3=1β2(α^)-β1(α^)]]>

其中，W₁表示几何扩散补偿因子，W₂表示校正反射角和偏移界面倾角的因子，W₃表示有限孔距补偿因子，r_g表示成像点到震源的距离，r_s表示成像点到检波器的距离，其中：

α^=tan(υ)=tan((θg-θs)/2)]]>

β^=tan(θ)=tan((θg+θs)/2)]]>

βi=(1+α^)2+4α^2hi2/z2-(1+α^)22α^2hi/z,i=1,2]]>

i＝1代表最近的检波器，i＝2代表最远的检波器，h表示从震源到检波器之间距离的一半，θ表示入射和反射之间的夹角，υ表示反射体倾角，θ_g表示成像点到检波器的入射线和到地表的垂线之间的夹角，θ_s表示成像点到震源点的入射线和到地表的垂线之间的夹角。