[發(fā)明專利]一種地震數(shù)據(jù)處理方法及裝置

有效

• 申請?zhí)枺?/td>: 201810460300.6
• 申請日: 2018-05-15
• 公開（公告）號: CN108897036B
• 公開（公告）日: 2020-04-10
• 發(fā)明（設(shè)計）人: 王偉;趙玉敏;李皓;周振曉;唐博文;張文波
• 申請（專利權(quán)）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司;中國石油大學(北京)
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三友知識產(chǎn)權(quán)代理有限公司 11127
• 代理人: 李輝

本申請?zhí)峁┮环N地震數(shù)據(jù)處理方法及裝置。所述方法包括：根據(jù)采集到的地震數(shù)據(jù)，獲得反映地層結(jié)構(gòu)特征的投影誤差算子；根據(jù)所述投影誤差算子，構(gòu)建地震數(shù)據(jù)反演的目標函數(shù)，所述目標函數(shù)中包括：所述地震數(shù)據(jù)、去除噪聲的地震數(shù)據(jù)、所述投影誤差算子之間的函數(shù)關(guān)系；對所述目標函數(shù)中進行迭代反演，獲得去除噪聲的目標地震數(shù)據(jù)。利用本申請中各實施例，反演出的地震數(shù)據(jù)有效的壓制了隨機噪聲對地震信號的影響，并能夠很好的保護有效信號的振幅，提高了地震資料同相軸的連續(xù)性和信噪比，增強了地震信號對地下結(jié)構(gòu)和油氣儲層的檢測精度，為地震資料構(gòu)造解釋和儲層預測提供了高質(zhì)量的基礎(chǔ)數(shù)據(jù)。

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

本申請屬于地球物理勘探技術(shù)領(lǐng)域，尤其涉及一種地震數(shù)據(jù)處理方法及裝置。

背景技術(shù)

地震勘探是一種利用人工地震技術(shù)探測地下結(jié)構(gòu)的勘探方法。它按照一定的方式人工激發(fā)地震波，利用稱之為檢波器的裝置接收來自地下的反射信號，通過對反射信號的處理和分析探測地下結(jié)構(gòu)。檢波器在接收地震信號的同時，也接收到大量來自地下和地表的隨機噪聲，降低了地震記錄信噪比，嚴重干擾地震信號反映地下結(jié)構(gòu)的能力，如何有效的壓制隨機噪聲和恢復弱信號是地震資料處理領(lǐng)域的重要研究內(nèi)容。

預測濾波方法是比較經(jīng)典、常用的隨機噪聲衰減方法，該方法不僅能夠在時間-空間域進行，也能在頻率-空間域進行。預測濾波方法的基本思想是依據(jù)地震信號的可預測性，即地震信號可以表示為相鄰地震道與預測濾波器的褶積。然而，現(xiàn)有技術(shù)中的預測濾波隨機噪聲衰減方法存在兩個問題，第一個問題是該在去噪過程中使用了兩個相互矛盾的噪聲模型，噪聲模型假設(shè)不一致；第二個問題是預測濾波方法容易傷害有效信號，尤其是當?shù)貙咏Y(jié)構(gòu)比較復雜時，其反射特征很難滿足預測濾波方法的要求，使得預測濾波方法保幅性很差，很難達到理想的隨機噪聲衰減結(jié)果。以上問題可能會導致地震數(shù)據(jù)處理結(jié)果不準確，如：噪聲去除不完全或者將部分有效信號去除等。因此，業(yè)內(nèi)亟需一種能夠提高地震數(shù)據(jù)處理準確性，進一步提高地震數(shù)據(jù)反映地下結(jié)果能力的實施方案。

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

本申請目的在于提供一種地震數(shù)據(jù)處理方法及裝置，提高了地震數(shù)據(jù)處理結(jié)果的準確性，為后續(xù)地震資料的處理提供了準確的數(shù)據(jù)基礎(chǔ)。

一方面本申請?zhí)峁┝艘环N地震數(shù)據(jù)處理方法，包括：

根據(jù)采集到的地震數(shù)據(jù)，獲得反映地層結(jié)構(gòu)特征的投影誤差算子；

根據(jù)所述投影誤差算子，構(gòu)建地震數(shù)據(jù)反演的目標函數(shù)，所述目標函數(shù)中包括：所述地震數(shù)據(jù)、去除噪聲的地震數(shù)據(jù)、所述投影誤差算子之間的函數(shù)關(guān)系；

對所述目標函數(shù)中進行迭代反演，獲得去除噪聲的目標地震數(shù)據(jù)。

進一步地，所述方法的另一個實施例中，所述投影誤差算子包括：

上式中，

h(x,t)表示所述投影誤差算子，x表示空間方向上的參數(shù)，x＝-px,-(px-1),…,(px-1),px，t表示時間方向上的參數(shù)，t＝-pt,-(pt-1),…,(pt-1),pt，px是空間方向的算子長度，pt是時間方向的算子長度，D表示地震數(shù)據(jù)組成的列向量，D＝d(nx,nt)，nx＝1,2,…N,nt＝1,2,…T，N表示所述地震數(shù)據(jù)的總道數(shù)，T表示每道采樣點數(shù)，I表示單位矩陣，λ表示權(quán)重系數(shù)，A表示所述地震數(shù)據(jù)重新排布后的矩陣，A的第i行的元素為d(nx,nt),m≤nx≤m+2px且nx≠m+px，n≤nt≤n+2pt，m＝(i-1)|(N-2px)+1，n＝(i-1)\(N-2px)+1。

進一步地，所述方法的另一個實施例中，所述目標函數(shù)包括：J＝||d-s||²+μF(hs)

上式中，J表示所述目標函數(shù)，s表示所述去除噪聲的地震數(shù)據(jù)，d表示所述地震數(shù)據(jù)，μ表示反演權(quán)重系數(shù)，h表示所述投影誤差算子的褶積矩陣，F(xiàn)表示稀疏約束函數(shù)。