[發(fā)明專利]一種基于邊界積分反算子的鬼波壓制方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201710799722.1
• 申請(qǐng)日: 2017-09-07
• 公開（公告）號(hào): CN107450103B
• 公開（公告）日: 2020-04-28
• 發(fā)明（設(shè)計(jì)）人: 陳寶書;管西竹;汪小將;陶杰;仝中飛;楊小椿;薛冬;張連群;李松康
• 申請(qǐng)（專利權(quán)）人: 中海石油（中國(guó)）有限公司;中海石油（中國(guó)）有限公司北京研究中心
• 主分類號(hào): G01V1/36
• 分類號(hào): G01V1/36;G01V1/38
• 代理公司: 北京紀(jì)凱知識(shí)產(chǎn)權(quán)代理有限公司 11245
• 代理人: 徐寧;孫楠

本發(fā)明涉及一種基于邊界積分反算子的鬼波壓制方法，其特征在于包括以下步驟：1)根據(jù)地震波傳播的邊界積分公式以及地震深度數(shù)據(jù)，計(jì)算得到鬼波的邊界積分算子公式；2)根據(jù)步驟1)得到的鬼波邊界積分算子公式，利用鬼波算子和壓制鬼波反算子的耦合關(guān)系，得到壓制鬼波反算子的計(jì)算公式；3)根據(jù)步驟2)中得到的壓制鬼波反算子的計(jì)算公式以及含有鬼波的地震數(shù)據(jù)，計(jì)算得到壓制鬼波后的地震數(shù)據(jù)。本發(fā)明能夠根據(jù)不同采集方式獲得的地震資料，得到壓制鬼波后的地震數(shù)據(jù)，且計(jì)算方法效率高，因而本發(fā)明可以廣泛應(yīng)用于海洋地震勘探過程中。

技術(shù)領(lǐng)域

本發(fā)明涉及海洋地震勘探領(lǐng)域的地震資料處理方法，尤其涉及一種基于邊界積分反算子的鬼波壓制方法。

背景技術(shù)

在海上地震勘探過程中，由于鬼波的影響，海上拖纜地震數(shù)據(jù)采集要想獲得較低頻率的地震資料，一般采用較深震源和較大拖纜沉放深度，為了獲得較高頻率地震資料一般采用較淺震源和較小拖纜沉放深度。為了研究中深層儲(chǔ)層的地質(zhì)情況，地震勘探的目標(biāo)既要得到中深層的目標(biāo)成像，同時(shí)要求高分辨的地震資料。這就需要消除鬼波的影響，以達(dá)到擴(kuò)展地震頻帶從而獲得更好的中深目標(biāo)地層的地質(zhì)信息的目的。

隨著近幾年地震采集設(shè)備、地震處理方法的改進(jìn)，國(guó)內(nèi)外研發(fā)了很多新型的地震采集方法，通過地震采集和地震信息處理設(shè)法消除鬼波的影響。相對(duì)于常規(guī)地震資料的處理流程，寬頻地震資料處理技術(shù)存在的主要問題之一是如何壓制或削弱鬼波，從而獲得更寬頻帶的地震數(shù)據(jù)。針對(duì)不同的采集方法提出來(lái)了可以減輕或消除鬼波的影響來(lái)獲得寬頻地震數(shù)據(jù)處理技術(shù)，可以獲得對(duì)更加精細(xì)的地震資料和高分辨率的地震成果剖面，但是由于采集方式的改變，在很大程度上增加了采集成本，同時(shí)針對(duì)新的采集方式的鬼波壓制方法處理過程復(fù)雜、難度較大，不能應(yīng)用于以前大量常規(guī)采集的地震資料。

發(fā)明內(nèi)容

針對(duì)上述問題，本發(fā)明的目的是提供一種能夠壓制不同采集方式獲取的海上地震資料中的鬼波的基于邊界積分反算子的鬼波壓制方法。

為實(shí)現(xiàn)上述目的，本發(fā)明采取以下技術(shù)方案：一種基于邊界積分反算子的鬼波壓制方法，其特征在于包括以下步驟：1)根據(jù)地震波傳播的邊界積分公式以及地震深度數(shù)據(jù)，計(jì)算得到鬼波的邊界積分算子；2)根據(jù)步驟1)得到的鬼波的邊界積分算子以及鬼波算子和壓制鬼波的反算子的耦合關(guān)系，得到壓制鬼波的反算子；3)根據(jù)得到的壓制鬼波的反算子以及含有鬼波的地震數(shù)據(jù)，計(jì)算得到壓制鬼波后的地震數(shù)據(jù)。

所述步驟1)中，鬼波的邊界積分算子計(jì)算方法包括以下步驟：1.1)根據(jù)采集的地震數(shù)據(jù)中的震源深度和電纜深度信息，得到地震采集的幾何關(guān)系圖，進(jìn)而得到電纜上所接收到的地震信號(hào)與分別來(lái)自地下和來(lái)自海面反射到達(dá)電纜檢波器的地震信號(hào)的關(guān)系；1.2)根據(jù)自由界面的邊界連續(xù)條件以及步驟1.1)中地震采集的幾何關(guān)系圖，得到海水面上接收到的地震信號(hào)與分別來(lái)自地下和來(lái)自海面反射到達(dá)海水面的地震信號(hào)的關(guān)系；1.3)根據(jù)波動(dòng)方程的邊界元積分公式，得到電纜上所接收到的地震信號(hào)與以來(lái)自地下到達(dá)海水面和來(lái)自海面反射到達(dá)海水面的地震信號(hào)的邊界積分公式的關(guān)系；1.4)根據(jù)步驟1.2)和步驟1.3)得到的關(guān)系，將電纜上所接收到的地震信號(hào)表示為來(lái)自地下到達(dá)海水面的地震信號(hào)的邊界積分形式，即不含鬼波的地震信號(hào)與鬼波的邊界積分算子的表達(dá)式。

所述步驟1.1)中，電纜上所接收到的地震信號(hào)u(r)為來(lái)自地下到達(dá)電纜檢波器的信號(hào)u_t(r)與來(lái)自海面反射到達(dá)電纜檢波器的信號(hào)u_g(r)之和，即：

u(r)＝u_t(r)+u_g(r)，

式中r為電纜所在的位置矢量。

所述步驟1.2)中，海水面上所接收到的地震信號(hào)u(r₀)為來(lái)自地下到達(dá)海水面的信號(hào)u_t(r₀)與來(lái)自海面反射到達(dá)海水面的信號(hào)u_g(r₀)之和且等于0，即：