[發明專利]一種處理三維地震數據的方法及裝置

在審

• 申請號: 201710341917.1
• 申請日: 2017-05-16
• 公開（公告）號: CN107272057A
• 公開（公告）日: 2017-10-20
• 發明（設計）人: 李海山;楊午陽;高建虎;張巧鳳;周春雷;楊慶;鄢高韓;杜炳毅
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產權代理有限公司11127
• 代理人: 李輝

1.一種處理三維地震數據的方法，其特征在于，該方法包括：

利用三維離散曲波變換方法對所述三維地震數據進行稀疏表示以獲取所述三維地震數據的三維稀疏表示系數；

利用硬閾值算法對所述三維稀疏表示系數進行閾值處理；

對所述閾值處理后的三維稀疏表示系數進行編碼。

2.根據權利要求1所述的方法，其特征在于，所述利用三維離散曲波變換方法對所述三維地震數據進行稀疏表示以獲取所述三維地震數據的三維稀疏表示系數的步驟包括利用下述公式來獲取所述三維稀疏表示系數：

其中，SparseCoef(j,l,k)為所述三維稀疏表示系數；Seis3D(n₁,n₂,n₃)為所獲取的三維地震數據；為三維離散曲波；為的共軛復數；n₁,n₂,n₃為三維地震數據樣點的下標，并且0≤n₁≤Nt,0≤n₂≤Nx,0≤n₃≤Ny，Nt,Nx,Ny分別為所述三維地震數據在時間方向及空間方向的采樣點數；j為尺度下標，并且1≤j≤Nj，Nj為分解的尺度數；l為分解的方向下標，并且8≤l≤Nl，Nl為分解的方向數；k＝(k₁,k₂,k₃)為所述三維稀疏表示系數的下標。

3.根據權利要求1所述的方法，其特征在于，所述利用硬閾值算法對所述三維稀疏表示系數進行閾值處理的步驟包括利用下述公式來對所述三維稀疏表示系數進行閾值處理：

其中，ThredCoef(j,l,k)為閾值處理后的三維稀疏表示系數；T為用戶設定的閾值。

4.根據權利要求1所述的方法，其特征在于，對所述閾值處理后的三維稀疏表示系數進行編碼的步驟包括利用下述公式對所述三維稀疏表示系數進行編碼：

其中，CodeCoef(m)為編碼后的所述三維稀疏表示系數；m為編碼后的所述三維稀疏表示系數的下標，其最大值為閾值處理后的非零三維稀疏表示系數的個數；ν為三維稀疏表示系數的數值。

5.根據權利要求1-4任一項所述的方法，其特征在于，該方法還包括：

對編碼后的所述三維稀疏表示系數進行存儲。

6.根據權利要求5所述的方法，其特征在于，該方法還包括：

讀取所存儲的所述三維稀疏表示系數；

對所讀取的所述三維稀疏表示系數進行解碼；

利用所述三維離散曲波變換方法的逆變換方法對解碼后的所述三維稀疏表示系數進行逆變換，以獲得壓縮恢復后的三維地震數據。

7.根據權利要求1-4任一項所述的方法，其特征在于，該方法還包括：

對編碼后的所述三維稀疏表示系數進行解碼；

利用所述三維離散曲波變換方法的逆變換方法對解碼后的所述三維稀疏表示系數進行逆變換，以獲得壓縮恢復后的三維地震數據。