[發明專利]基于插值與去混疊同時處理的地震數據分離方法及系統

在審

• 申請號: 201910543575.0
• 申請日: 2019-06-21
• 公開（公告）號: CN110244361A
• 公開（公告）日: 2019-09-17
• 發明（設計）人: 王本鋒;耿建華
• 申請（專利權）人: 同濟大學
• 主分類號: G01V1/36
• 分類號: G01V1/36;G01V1/28
• 代理公司: 上?？剖⒅R產權代理有限公司 31225
• 代理人: 丁云

1.一種基于插值與去混疊同時處理的地震數據分離方法，其特征在于，該方法包括如下步驟：

(1)基于觀測地震數據提取偽分離地震數據；

(2)基于Curvelet變換對偽分離地震數據進行插值與去混疊同時處理迭代更新直至迭代結束，輸出主震源、副震源數據，完成地震數據分離；

步驟(1)具體為：

(11)確定混疊算子Γ；

(12)獲取偽分離地震數據：F_t^HΓ^HF_td，其中，F_t、F_t^H分別為時間方向的一維正、反傅里葉變換，Γ^H為混疊算子Γ的共軛算子，d為觀測地震數據；

步驟(2)具體為：

(21)賦值n＝1，其中，為主震源第n次迭代值，為副震源第n次迭代值；

(22)確定去混疊項和插值項并更新主震源和副震源數據：

其中，和為去混疊項，和為插值項，為主震源第n+1次迭代值，為副震源第n+1次迭代值，C、C^H為Curvelet正、反變換，T_λ為硬閾值函數，R為采樣算子，且R滿足可交換性，即RF_t＝F_tR，RΓ＝ΓR；

(23)判斷n是否小于N，若是則賦值n＝n+1并返回步驟(22)，否則輸出主震源和副震源數據，其中，N為最大迭代次數。

2.根據權利要求1所述的一種基于插值與去混疊同時處理的地震數據分離方法，其特征在于，所述的混疊算子Γ具體確定為：

其中，e為自然常數，i為虛數單位，ω為圓頻率，t_k為副震源中第k炮的時間延遲，k＝1,...,M，M為副震源中炮的數量。

3.一種基于插值與去混疊同時處理的地震數據分離系統，該系統包括：

地震數據提取模塊：該模塊基于觀測地震數據提取偽分離地震數據；

插值與去混疊分離模塊：該模塊對偽分離地震數據進行插值與去混疊并迭代更新，得到分離的主震源、副震源數據；

輸出模塊：該模塊輸出分離的主震源、副震源數據；

所述的地震數據提取模塊獲取的偽分離地震數據為F_t^HΓ^HF_td，其中，F_t、F_t^H分別為時間方向的一維正、反傅里葉變換，Γ^H為混疊算子Γ的共軛算子，d為觀測地震數據；

所述的插值與去混疊分離模塊執行以下步驟完成主震源、副震源數據分離：

(a)賦值n＝1，其中，為主震源第n次迭代值，為副震源第n次迭代值；

(b)確定去混疊項和插值項并更新主震源和副震源數據：

其中，和為去混疊項，和為插值項，為主震源第n+1次迭代值，為副震源第n+1次迭代值，C、C^H為Curvelet正、反變換，T_λ為硬閾值函數，R為采樣算子，且R滿足可交換性，即RF_t＝F_tR，RΓ＝ΓR；

(c)判斷n是否小于N，若是則賦值n＝n+1并返回步驟(b)，否則輸出主震源和副震源數據，其中，N為最大迭代次數。

4.根據權利要求3所述的一種基于插值與去混疊同時處理的地震數據分離系統，其特征在于，所述的混疊算子Γ具體確定為：

其中，e為自然常數，i為虛數單位，ω為圓頻率，t_k為副震源中第k炮的時間延遲，k＝1,...,M，M為副震源中炮的數量。