[發明專利]一種起伏海面的虛反射壓制方法和裝置

有效

• 申請號: 201611087106.5
• 申請日: 2016-12-01
• 公開（公告）號: CN106646596B
• 公開（公告）日: 2019-08-06
• 發明（設計）人: 馬光凱;周錚錚;錢忠平;趙波;任曉喬;尹天奎
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 李輝

1.一種起伏海面的虛反射壓制方法，其特征在于，所述方法包括：

基于起伏海面反射系數計算公式，對獲取的頻率波數域的地震數據對應的角頻率和波數、與所述頻率波數域的地震數據相對應的海浪高度、以及相對應的海水速度進行計算得到反射系數，所述起伏海面反射系數計算公式以預設參數擬合獲得；

利用所述頻率波數域的地震數據對應的角頻率計算得到白噪聲系數；

利用所述反射系數、所述白噪聲系數分別計算得到預設數量的拖纜掃描深度所對應的頻率波數域的虛反射壓制算子；

將所述頻率波數域的虛反射壓制算子分別與頻率波數域的地震數據相乘得到預設數量的頻率波數域壓制數據；

將所述預設數量的頻率波數域壓制數據進行時間空間參量的二維傅里葉反變換，得到預設數量的時間空間域壓制數據；

以預設的時窗間隔劃分所述時間空間域壓制數據，分別計算同一時窗內所述預設數量的拖纜掃描深度所對應的頻域數據的絕對值的和，將絕對值的和最小所對應的拖纜掃描深度的時間空間域壓制數據作為相應時窗的壓制虛反射后地震數據。

2.根據權利要求1所述的方法，其特征在于，所述起伏海面反射系數計算公式包括如下計算公式：

上式中，r表示反射系數；k_x表示頻率波數域地震數據在x方向對應的波數；ω表示頻率波數域的地震數據對應的角頻率；h表示海浪高度；c表示海水速度；e表示自然常數；π表示圓周率。

3.根據權利要求1或2任意一項所述的方法，其特征在于，所述預設數量的拖纜掃描深度包括：

根據預設掃描范圍、掃描間隔和獲取的檢波器深度確定出的拖纜掃描深度。

4.根據權利要求1或2任意一項所述的方法，其特征在于，所述利用所述頻率波數域的地震數據對應的角頻率計算得到白噪聲系數包括采用下述公式計算：

上式中，white表示白噪聲系數；ω表示頻率波數域的地震數據對應的角頻率。

5.根據權利要求1或2任意一項所述的方法，其特征在于，所述利用所述反射系數、所述白噪聲系數分別計算得到預設數量的拖纜掃描深度所對應的頻率波數域的虛反射壓制算子包括采用下述公式計算：

上式中，A(ω,k_x)_z表示拖纜掃描深度為z時頻率波數域的虛反射壓制算子，z＝(i-1)*d_z，i表示掃描序號，i取值范圍為[0，n]；n表示掃描次數，z₁表示獲取的拖纜掃描的起始深度；z₂表示獲取的拖纜掃描的終止深度；d_z表示預設掃描間隔；white表示拖纜掃描深度為z時的白噪聲系數；G(ω,k_x,r)_z表示拖纜掃描深度為z時包括虛反射的地震波傳播算子，表示G(ω,k_x,r)_z的共軛；r表示反射系數；k_x表示頻率波數域地震數據在x方向對應的波數；ω表示頻率波數域的地震數據對應的角頻率。

6.根據權利要求1或2任意一項所述的方法，其特征在于，所述獲取的頻率波數域的地震數據對應的角頻率和波數包括：

對采集的地震數據中單炮地震數據分別進行時間空間參量的二維傅里葉變換后得到的頻率波數域的地震數據對應的角頻率和波數。

7.根據權利要求1或2任意一項所述的方法，其特征在于，所述預設數量的拖纜掃描深度所對應的頻域數據包括：

對所述預設數量的時間空間域壓制數據進行一維傅里葉變換后得到的數據。