三分量地震波積分偏移方法和裝置製造方法
2023-12-04 09:11:46 1
三分量地震波積分偏移方法和裝置製造方法
【專利摘要】本發明提供了一種三分量地震波積分偏移方法和裝置，其中，該方法包括：計算三維速度模型空間中當前網格點的積分加權因子；根據當前網格點的位置，在走時表中獲取地震波從炮點到檢波點之間的全走時；根據三分量地震波的振幅，獲取該網格點到檢波器的地震波出射方向矢量；對所述方向矢量和三分量地震波的振幅矢量做點乘運算得到一個標量振幅；根據該網格點的積分加權因子對該網格點對應的標量振幅做偏移的加權運算，並將加權運算的結果與該網格點原存儲的數據進行積分求和，得到該網格點的偏移圖像。本發明解決了現有技術中三分量檢波器成像效果不好、噪聲幹擾大的技術問題，達到了有效提高三分量檢波器成像效果的技術效果。
【專利說明】H分量地震波積分偏移方法和裝置
【技術領域】
[0001] 本發明設及地震數據處理【技術領域】，特別設及一種=分量地震波積分偏移方法和 裝置。
【背景技術】
[0002] =分量檢波器是一種可W記錄地下地層物理更全面的信息的檢波器，然而，目前 的技術仍然只用單一分量地震波或將=分量轉換成單一分量進行地震偏移成像（W下簡 稱單分量偏移），前者浪費了很多有價值的信息，後者所用信息量有所增加，但丟失了=維 地震波的方向性。
[0003] 圖1所示是由震源發出的地震波照射到反射體A後，由檢波器接收，波的傳播路徑 如實線所示，在對此波進行偏移計算時按照它的走時還原反射體到它應在的位置A，但同時 也會將一虛假反射體還原到A'的位置，因A'是A相對於W震源和檢波器方向為軸的鏡像 位置，因為由震源到A'再到檢波器波的虛假走時和經A的真實走時是相同的，而單分量偏 移方法並不知道地震波震源的方向，只好將反射的能量放在真實反射體A和它的鏡像反射 體A'位置上，而A'的出現使偏移結果出現噪聲，因此，採用目前的單分量地震波偏移方法 在地質界面不存在的地方也會產生虛假的界面，難W得到好的成像效果。
【發明內容】
[0004] 本發明實施例提供了一種=分量地震波積分偏移方法，W解決現有技術中=分量 檢波器成像效果不好的技術問題，該方法包括：
[0005] 計算S維速度模型空間中當前網格點的積分加權因子；
[0006] 根據當前網格點的位置，在走時表中獲取地震波從炮點到檢波點之間的全走時；
[0007] 根據=分量地震波的振幅，獲取該網格點到檢波器的地震波出射方向矢量；
[000引對所述方向矢量和=分量地震波的振幅矢量做點乘運算得到一個標量振幅；
[0009] 根據該網格點的積分加權因子對該網格點對應的標量振幅做偏移的加權運算，並 將加權運算的結果與該網格點原存儲的數據進行積分求和，得到該網格點的偏移圖像。
[0010] 在一個實施例中，計算S維速度模型空間中當前網格點的積分加權因子，包括：
[0011] 對所述S維速度模型空間中的網格點，由孔距參數計算出相對此網格點的積分加 權因子。
[0012] 在一個實施例中，所述孔距參數按照W下方式獲取：
[0013] 獲取共炮點反射波地震數據，其中，所述地震數據為=分量地震數據；
[0014] 根據所述共炮點反射波地震數據的炮點和檢波點之間的幾何位置，計算所述共炮 點反射波地震數據的孔距參數。
[0015] 在一個實施例中，由孔距參數計算出相對此網格點的積分加權因子，包括：
[0016] 按照W下公式計算所述積分加權因子：
【權利要求】
1. 一種三分量地震波積分偏移方法，其特徵在於，包括： 計算三維速度模型空間中當前網格點的積分加權因子； 根據當前網格點的位置，在走時表中獲取地震波從炮點到檢波點之間的全走時； 根據三分量地震波的振幅，獲取該網格點到檢波器的地震波出射方向矢量； 對所述方向矢量和三分量地震波的振幅矢量做點乘運算得到一個標量振幅； 根據該網格點的積分加權因子對該網格點對應的標量振幅做偏移的加權運算，並將加 權運算的結果與該網格點原存儲的數據進行積分求和，得到該網格點的偏移圖像。
2. 如權利要求1所述的方法，其特徵在於，計算三維速度模型空間中當前網格點的積 分加權因子，包括： 對所述三維速度模型空間中的網格點，由孔距參數計算出相對此網格點的積分加權因 子。
3. 如權利要求2所述的方法，其特徵在於，所述孔距參數按照以下方式獲取： 獲取共炮點反射波地震數據，其中，所述地震數據為三分量地震數據； 根據所述共炮點反射波地震數據的炮點和檢波點之間的幾何位置，計算所述共炮點反 射波地震數據的孔距參數。
4. 如權利要求2所述的方法，其特徵在於，由孔距參數計算出相對此網格點的積分加 權因子，包括： 按照以下公式計算所述積分加權因子：
其中，11表不幾何擴散補償因子，W2表不校正反射角和偏移界面傾角的因子，W3表不有 限孔距補償因子，rg表示成像點到震源的距離，rs表示成像點到檢波器的距離，其中：
i= 1代表最近的檢波器，i= 2代表最遠的檢波器，h表示從震源到檢波器之間距離 的一半，0表示入射和反射之間的夾角，u表示反射體傾角，0g表示成像點到檢波器的入 射線和到地表的垂線之間的夾角，0 3表示成像點到震源點的入射線和到地表的垂線之間 的夾角。
5. 如權利要求1所述的方法，其特徵在於，根據該網格點的積分加權因子對該網格點 對應的標量振幅做偏移的加權運算，包括： 按照以下公式對該網格點對應的標量振幅做偏移的加權運算：
其中，M(X)表示網格點的成像點X處的三維三分量地震積分偏移， + 表示地震波經成像點X，再到檢波器的全程走時為ts(xs,x)+1^(x，xj 時的三分量地震波的振幅，及表示由X點到檢波點的地震波出射方向矢量，?表示點乘，W(Xs,XmX)表示積分加權因子，ts (Xs,X)表示震源到成像點的走時，tjx，X1^表示成像點到 檢波器的走時，\表示震源的空間坐標，X表示成像點的空間坐標，表示檢波器的空間坐
6. 如權利要求1至5中任一項所述的方法，其特徵在於，在根據當前網格點的位置，在 走時表中獲取地震波從炮點到檢波點之間的全走時之前，所述方法還包括： 確定所述三維速度模型中每一個網格點到所有炮點和所有檢波點的走時表。
7. 如權利要求6所述的方法，其特徵在於，確定所述三維速度模型中每一個網格點到 所有炮點和所有檢波點的走時表，包括： 在所述三維速度模型空間中設定一個震源點，震源點所產生的地震波向四周傳播射 線，由近到遠計算模型空間中每個網格點的射線的到時； 將確定的每個網格點的射線的到時作為該震源點的走時表。
8. -種三分量地震波積分偏移裝置，其特徵在於，包括： 積分加權因子計算單元，用於計算三維速度模型空間中當前網格點的積分加權因子； 全走時獲取單元，用於根據當前網格點的位置，在走時表中獲取地震波從炮點到檢波 點之間的全走時； 方向矢量獲取單元，用於根據三分量地震波的振幅，獲取該網格點到檢波器的地震波 出射方向矢量； 標量振幅獲取單元，用於對所述方向矢量和三分量地震波的振幅矢量做點乘運算得到 一個標量振幅； 偏移計算單元，用於根據該網格點的積分加權因子對該網格點對應的標量振幅做偏移 的加權運算，並將加權運算的結果與該網格點原存儲的數據進行積分求和，得到該網格點 的偏移圖像。
9. 如權利要求8所述的裝置，其特徵在於，所述積分加權因子計算單元具體用於對所 述三維速度模型空間中的網格點，由孔距參數計算出相對此網格點的積分加權因子，並按 照以下公式由孔距參數計算出相對此網格點的積分加權因子：
其中，11表不幾何擴散補償因子，W2表不校正反射角和偏移界面傾角的因子，W3表不有 限孔距補償因子，rg表示成像點到震源的距離，rs表示成像點到檢波器的距離，其中：
i= 1代表最近的檢波器，i= 2代表最遠的檢波器，h表示從震源到檢波器之間距離 的一半，0表示入射和反射之間的夾角，u表示反射體傾角，0g表示成像點到檢波器的入 射線和到地表的垂線之間的夾角，0 3表示成像點到震源點的入射線和到地表的垂線之間 的夾角。
10.如權利要求8所述的裝置，其特徵在於，所述偏移計算單元具體用於按照以下公式 對該網格點對應的標量振幅做偏移的加權運算：
其中，M(X)表示網格點的成像點X處的三維三分量地震積分偏移，P,.〇s,x) +ai,;))表示地震波經成像點X，再到檢波器的全程走時為ts(xs，x)+tr(x，X) 時的三分量地震波的振幅，》表示由x點到檢波點的地震波出射方向矢量，?表示點乘，W(Xs,XmX)表示積分加權因子，ts (Xs,X)表示震源到成像點的走時，tjx，X1^表示成像點到 檢波器的走時，\表示震源的空間坐標，X表示成像點的空間坐標，表示檢波器的空間坐
【文檔編號】G01V1/36GK104502972SQ201410759540
【公開日】2015年4月8日 申請日期:2014年12月11日 優先權日:2014年12月11日
【發明者】王大鵬, 餘剛, 王熙明, 李彥鵬, 張固瀾 申請人:中國石油集團東方地球物理勘探有限責任公司, 王大鵬