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

在審

• 申請?zhí)枺?/td>: 201910543575.0
• 申請日: 2019-06-21
• 公開（公告）號: CN110244361A
• 公開（公告）日: 2019-09-17
• 發(fā)明（設計）人: 王本鋒;耿建華
• 申請（專利權）人: 同濟大學
• 主分類號: G01V1/36
• 分類號: G01V1/36;G01V1/28
• 代理公司: 上海科盛知識產權代理有限公司 31225
• 代理人: 丁云

本發(fā)明涉及一種基于插值與去混疊同時處理的地震數據分離方法及系統(tǒng)，該方法包括如下步驟：(1)基于觀測地震數據提取偽分離地震數據；(2)基于Curvelet變換對偽分離地震數據進行插值與去混疊同時處理迭代更新直至迭代結束，輸出主震源、副震源數據，完成地震數據分離。與現有技術相比，本發(fā)明避免了數據不規(guī)則性對去混疊精度的影響及混疊噪聲對插值重建效果的影響，得到去混疊后的完整地震數據，有利于提高后續(xù)常規(guī)地震數據處理及偏移成像的精度。

技術領域

本發(fā)明涉及一種地震數據分離方法及系統(tǒng)，尤其是涉及一種基于插值與去混疊同時處理的地震數據分離方法及系統(tǒng)。

背景技術

隨著高效地震數據采集的深入，稀疏采集、混疊采集扮演著重要的角色。而采集得到的混疊、不規(guī)則地震數據，對傳統(tǒng)的數據處理、反演及偏移等方法帶來巨大挑戰(zhàn)，因此插值技術及去混疊技術成為高效采集地震數據處理的重要環(huán)節(jié)。插值方法一般可分為四類：基于波動方程的插值方法、基于預測濾波的插值方法、基于稀疏變換的插值方法及基于降秩的插值方法。隨著壓縮感知理論的完善，基于稀疏變換的插值重建方法扮演了越來越重要的角色。去混疊方法包括基于濾波的方法和基于反演的方法，濾波類方法利用信號相干性及混疊噪聲隨機性的特性，壓制噪聲、提取信號；反演類的方法利用信號的相干性、稀疏性等約束，進行混疊數據的有效分離。而地震數據的不規(guī)則性會影響去混疊的精度，混疊噪聲的存在會影響插值重建的效果，因此需研究插值去混疊同時處理方法。針對不規(guī)則的混疊數據，基于濾波的方法可得到完整的去混疊數據，但算法參數設置相對困難。因此需要一種地震數據插值去混疊同時處理的方法得到完整的去混疊地震數據。

發(fā)明內容

本發(fā)明的目的就是為了克服上述現有技術存在的缺陷而提供一種基于插值與去混疊同時處理的地震數據分離方法及系統(tǒng)。

本發(fā)明的目的可以通過以下技術方案來實現：

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

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

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

步驟(1)具體為：

(11)確定混疊算子Γ；

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

所述的混疊算子Γ具體確定為：

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

步驟(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為最大迭代次數。