[發明專利]一種利用GPU提高波場延拓計算效率的方法

在審

• 申請號: 201210315284.4
• 申請日: 2012-08-30
• 公開（公告）號: CN103675895A
• 公開（公告）日: 2014-03-26
• 發明（設計）人: 孔祥寧;張慧宇;段心標;徐兆濤;張兵;孫武亮
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 郭韞

1.一種利用GPU提高波場延拓計算效率的方法，其特征在于：所述方法首先擴大震源波場和記錄波場的數據網格，然后利用GPU進行多線程并行計算，并應用CUDA提供的二維傅里葉變換函數實現波場延拓。

2.根據權利要求1所述的利用GPU提高波場延拓計算效率的方法，其特征在于：所述方法包括以下步驟：

(1)擴大震源波場和記錄波場的數據網格得到處理后的震源波場和記錄波場；

(2)應用裂步傅里葉法對所述處理后的震源波場和記錄波場進行波場延拓；

(3)設延拓深度iz＝1；

(4)利用CUDA提供的二維傅里葉變換函數對震源波場進行二維正傅里葉變換；

(5)利用CUDA提供的二維傅里葉變換函數對記錄波場進行二維正傅里葉變換；

(6)參考速度場的相移計算：利用步驟用(4)和步驟(5)的結果進行相移計算，參考速度場是平均速度，即輸入速度場的平均值；

(7)利用CUDA提供的二維傅里葉變換函數對震源波場進行二維反傅里葉變換；

(8)利用CUDA提供的二維傅里葉變換函數對記錄波場進行二維反傅里葉變換；

(9)擾動速度場的波場計算與相關成像：利用步驟(7)和步驟(8)的結果進行相移計算，擾動速度場是輸入速度場與參考速度場的差值；

(10)判斷iz+1是否大于NZ，所述NZ＝延拓深度/延拓步長，如果是，則轉入步驟(11)，如果否，則返回步驟(4)；

(11)結束。

3.根據權利要求2所述的利用GPU提高波場延拓計算效率的方法，其特征在于：所述步驟(1)是這樣實現的：對震源波場和記錄波場的長度進行鑲邊處理，即當震源波場和記錄波場的長度小于1000時，將其長度鑲邊處理為128的倍數，對于不足的數據處填充零；當震源波場和記錄波場的長度大于1000而小于2048時，將其長度鑲邊處理成256的倍數，對于不足的數據處填充零。

4.根據權利要求2所述的利用GPU提高波場延拓計算效率的方法，其特征在于：所述步驟(2)至(11)都是在GPU上進行的。