[發明專利]一種利用GPU進行逆時偏移提取角度道集的方法

有效

• 申請號: 201210184592.8
• 申請日: 2012-06-06
• 公開（公告）號: CN103472481A
• 公開（公告）日: 2013-12-25
• 發明（設計）人: 李敏;李振華;任鵬
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/32
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 郭韞

1.一種利用GPU進行逆時偏移角度道集提取的方法，其特征在于：所述方法包括以下步驟：

(1)炮點波場正向延拓；

(2)并行執行逆時偏移和角度信息統計；

(3)角度道集提取。

2.根據權利要求1所述的利用GPU進行逆時偏移角度道集提取的方法，其特征在于：所述步驟(1)包括以下步驟：

(A1)在炮點位置放入人工子波震源；

(A2)利用有限差分方法對步驟(A1)中的人工子波震源激發出的地震波的波場沿時間方向進行正向延拓，得到不同時刻的炮點波場；

所述步驟(A2)是在GPU上完成的。

3.根據權利要求2所述的利用GPU進行逆時偏移角度道集提取的方法，其特征在于：所述步驟(2)中，所述逆時偏移包括以下步驟：

(B1)將步驟(A2)中得到的不同時刻的炮點波場存儲在波場數據緩沖池中；

(B2)將從步驟(A2)得到的炮點波場沿時間方向正向延拓至最大時刻，得到炮點最大時刻的波場；正向延拓是指沿時間遞增的方向延拓波場；

(B3)從野外采集的地震數據中讀取檢波點最大時刻的波場，然后從最大時刻開始，將炮點波場和檢波點波場沿時間方向反向延拓，得到炮點和檢波點不同時刻的波場；反向延拓是指沿時間遞減的方向延拓波場；

(B4)將每個時刻的炮點波場與檢波點波場相乘得到成像結果；

(B5)每個炮點對應一個單炮，將每個單炮對應的步驟(B4)中的所述成像結果作為單炮偏移成像結果單獨存儲在磁盤上；

(B6)每個從野外采集的地震數據包含很多單炮，利用每個單炮對應的步驟(B4)中的所述成像結果做多炮疊加得到多炮疊加的偏移成像結果；所述多炮疊加是指把多個單炮成像結果針對地下相同位置處的成像值進行累加。

4.根據權利要求3所述的利用GPU進行逆時偏移角度道集提取的方法，其特征在于：所述步驟(2)中，所述角度信息統計包括以下步驟：

(C1)角度信息統計：將所述波場數據緩沖池中的炮點波場按照時間順序存儲，并且提取每個時刻每個空間位置的振幅和傳播角度；

(C2)完成步驟(A2)后，將步驟(C1)中得到的炮點波場在每個空間位置最大振幅時刻的傳播角度作為角度信息數據體存儲存儲在磁盤上；

(C3)在所述步驟(B6)得到的所述多炮疊加的偏移成像結果上，提取地層傾角作為地層傾角信息數據體存儲在磁盤上。

5.根據權利要求4所述的利用GPU進行逆時偏移角度道集提取的方法，其特征在于：所述步驟(B2)至(B5)與步驟(C1)至(C2)并行執行，步驟(B6)執行完畢后再執行步驟(C3)。

6.根據權利要求5所述的利用GPU進行逆時偏移角度道集提取的方法，其特征在于：所述步驟(3)包括：

(D1)從磁盤上讀取步驟(C2)和(C3)中得到的角度信息數據體和地層傾角信息數據體；

(D2)根據所述角度信息數據體和地層傾角信息數據體計算入射角：

入射角θ等于炮點波場在每一空間位置處達到最大振幅時刻的傳播角度β與地層傾角α一半的差，公式表達為：θ＝β-0.5α；

(D3)從磁盤上讀取步驟(B5)輸出的單炮偏移成像結果，然后將單炮偏移成像結果輸出到對應于步驟(D2)得到的入射角的道集中，得到最終的角度道集數據體。

7.根據權利要求3所述的利用GPU進行逆時偏移角度道集提取的方法，其特征在于：所述數據緩沖池中的基本數據結構包括一個浮點型數據、一個整型數據和一個浮點型數組，其中浮點型數據存儲當前時刻值、整型數據存儲所述浮點數組的長度、浮點型數組存儲當前時刻的波場。