[發明專利]三分量地震波積分偏移方法和裝置

有效

• 申請號: 201410759540.8
• 申請日: 2015-08-03
• 公開（公告）號: CN104502972A
• 公開（公告）日: 2015-07-29
• 發明（設計）人: 王大鵬;余剛;王熙明;李彥鵬;張固瀾
• 申請（專利權）人: 中國石油集團東方地球物理勘探有限責任公司;王大鵬
• 主分類號: G01V1/36
• 分類號: G01V1/36;G01V1/30;G06F19/00
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 王天堯

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表示成像點到震源點的入射線和到地表的垂線之間的夾角。