專利名稱：基于分?jǐn)?shù)階傅里葉變換的相干噪聲抑制方法
技術(shù)領(lǐng)域：
本發(fā)明屬于石油地震勘探數(shù)據(jù)處理領(lǐng)域，涉及用于消除噪聲干擾的技術(shù)，更具體地說，涉及一種用于利用分?jǐn)?shù)階傅里葉變換來抑制相干噪聲的方法。
背景技術(shù)：
利用地下介質(zhì)彈性和密度的差異，通過觀測和分析大地對人工激發(fā)地震波的響應(yīng)，推斷地下巖層的性質(zhì)和形態(tài)的地球物理勘探方法稱為地震勘探。地震勘探是鉆探前勘測石油與天然氣資源的重要手段，在煤田和工程地質(zhì)勘查、區(qū)域地質(zhì)研究和地殼研究等方面，也得到廣泛應(yīng)用。在地震勘探中，最為基礎(chǔ)和重要的工作是如何采集到野外地震數(shù)據(jù)并對其進(jìn)行有效的處理和解釋。為了能夠獲得有效的地震數(shù)據(jù)，需要在野外設(shè)置大量的監(jiān)測儀器，其中，每一個監(jiān)測儀器采集到的數(shù)據(jù)稱作一道地震數(shù)據(jù)。在采集到的地震資料中，不可避免地存在著一定數(shù)量的噪聲干擾。對于山地地震資料而言，其噪聲來源十分復(fù)雜，噪聲對有效波有較強(qiáng)的干擾，尤其是地震資料中普遍存在著相當(dāng)嚴(yán)重的線性干擾，如淺層強(qiáng)能量的地表多次折射波和面波，中層、深層的強(qiáng)聲波干擾，以及頻帶很寬的地表直達(dá)波等。一般情況下，這些干擾破壞了有效反射信號，嚴(yán)重時覆蓋了整個地震記錄，完全淹沒了有效信號，大大降低了地震資料的信噪比。對這些干擾，常規(guī)的處理方法是采用f_k濾波法、τ -ρ變換法、中值濾波法、速度陷波法等。這些方法雖然對相干干擾的消除有一定的作用，但它們本身都存在一定的局限性。例如，f_k濾波技術(shù)要求均勻的空間采樣，而陸地?cái)?shù)據(jù)尤其3D數(shù)據(jù)常常難以滿足這一苛刻條件；經(jīng)典的τ -ρ變換也存在同樣的問題；速度陷波濾波器則要求滿足噪聲為線狀噪聲的這一條件。這些方法都對不同道間的地震數(shù)據(jù)的振幅變化相當(dāng)敏感，即使細(xì)微的變化也會使濾波的效果變得極差。所以如果這些相干噪聲不能得到很好的抑制，則最終會影響地震記錄的疊加成像質(zhì)量及地震剖面的橫向分辨率。

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種能夠基于分?jǐn)?shù)階傅里葉變換來抑制地震資料中的相干噪聲的方法，從而提高在較復(fù)雜的地質(zhì)環(huán)境中獲得的地震數(shù)據(jù)(尤其是低信噪比的地震數(shù)據(jù))的處理質(zhì)量。根據(jù)本發(fā)明的一方面，提供一種用于在監(jiān)測的地震數(shù)據(jù)中抑制相干噪聲的方法， 所述方法包括在監(jiān)測的地震數(shù)據(jù)中估計(jì)相干噪聲的頻帶和視速度范圍，并提取所述頻帶內(nèi)的地震數(shù)據(jù)；在視速度范圍內(nèi)，針對提取的所述頻帶內(nèi)的地震數(shù)據(jù)來確定與視速度對應(yīng)的視速度方向；基于所述視速度和對應(yīng)的視速度方向來確定視傾角；根據(jù)視傾角對提取的所述頻帶內(nèi)的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換；對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測，以得到預(yù)測的相干噪聲；以及從監(jiān)測的地震數(shù)據(jù)中消除預(yù)測的相干噪聲。估計(jì)相干噪聲的頻帶和視速度范圍的步驟可包括人工估計(jì)監(jiān)測的地震數(shù)據(jù)中相干噪聲的頻帶，并估計(jì)出相干噪聲的視速度范圍。
提取所述頻帶內(nèi)的地震數(shù)據(jù)的步驟可包括在估計(jì)出相干噪聲的頻帶之后，構(gòu)建相應(yīng)的帶通濾波器，應(yīng)用所述帶通濾波器對監(jiān)測的地震數(shù)據(jù)進(jìn)行濾波，以得到濾波后的所述頻帶內(nèi)的地震數(shù)據(jù)。確定視速度方向的步驟可包括在估計(jì)出的視速度范圍內(nèi)，選取一個或多個視速度來對提取的所述頻帶內(nèi)的地震數(shù)據(jù)依次進(jìn)行掃描，并由此確定與所述一個或多個視速度對應(yīng)的視速度方向。當(dāng)對提取的所述頻帶內(nèi)的地震數(shù)據(jù)依次進(jìn)行掃描時，沿著各個方向?qū)λ鲱l帶內(nèi)的地震數(shù)據(jù)進(jìn)行采樣，從而在每個方向上得到一組采樣數(shù)據(jù)，然后對每個方向上的一組采樣數(shù)據(jù)相加求和，將絕對值最大的求和值所對應(yīng)的方向確定為相應(yīng)視速度的視速度方向。將視速度的方向確定為視傾角的方向。確定視傾角的步驟可包括對于二維疊前地震數(shù)據(jù)，基于視速度和監(jiān)測的地震數(shù)據(jù)中的地震道間距來確定視傾角的大小。視傾角的大小=(1000X地震道間距)/視速度，其中，地震道間距的單位為米，視速度的單位為米/秒。確定視傾角的步驟包括對于三維疊前地震數(shù)據(jù)，基于視速度和監(jiān)測的地震數(shù)據(jù)中的偏移距來確定視傾角的大小。根據(jù)視傾角對提取的所述頻帶內(nèi)的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換的步驟可包括將提取的所述頻帶內(nèi)的地震數(shù)據(jù)與視傾角下的分?jǐn)?shù)傅里葉變換核矩陣相乘，以得到分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)。對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測的步驟可包括采用F-X域?yàn)V波算法對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測。對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測的步驟可包括利用在所述頻帶外的頻段求取的預(yù)測因子對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測。所述頻帶外的頻段是指相干噪聲較少的頻段。所述相干噪聲較少的頻段為中頻段。利用在所述頻帶外的頻段求取的預(yù)測因子對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測的步驟可包括用求取的預(yù)測因子與所述頻帶內(nèi)的地震數(shù)據(jù)進(jìn)行褶積，從而得到相應(yīng)的相干預(yù)測濾波結(jié)果。用求取的預(yù)測因子與所述頻帶內(nèi)的地震數(shù)據(jù)進(jìn)行褶積的步驟可包括根據(jù)中頻段與相干噪聲的所述頻帶之間的頻率關(guān)系，用求取的預(yù)測因子與所述頻帶內(nèi)的地震數(shù)據(jù)進(jìn)行褶積。從監(jiān)測的地震數(shù)據(jù)中消除預(yù)測的相干噪聲的步驟可包括將預(yù)測的相干噪聲進(jìn)行分?jǐn)?shù)階傅里葉逆變換，以得到時空域中的相干噪聲，然后從監(jiān)測到的地震數(shù)據(jù)中減去時空域中的相干噪聲，從而得到抑制相干噪聲后的有效地震數(shù)據(jù)。


通過下面結(jié)合附圖進(jìn)行的對實(shí)施例的描述，本發(fā)明的上述和/或其它目的和優(yōu)點(diǎn)將會變得更加清楚，其中圖1是示出根據(jù)本發(fā)明示例性實(shí)施例的相干噪聲抑制方法的流程圖2A示出根據(jù)本發(fā)明示例性實(shí)施例監(jiān)測到的地震數(shù)據(jù)的示意圖；圖2B示出根據(jù)本發(fā)明示例性實(shí)施例過濾出的地震數(shù)據(jù)的示意圖；圖2C示出根據(jù)本發(fā)明示例性實(shí)施例對過濾出的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換的示意圖；圖2D示出根據(jù)本發(fā)明示例性實(shí)施例對過濾出的地震數(shù)據(jù)執(zhí)行相干預(yù)測的示意圖；圖2E示出根據(jù)本發(fā)明示例性實(shí)施例抑制相干噪聲后的地震數(shù)據(jù)的示圖；圖3A到圖3C示出根據(jù)本發(fā)明示例性實(shí)施例針對疊前二維數(shù)據(jù)來抑制相干噪聲的示圖；以及圖4A到圖4C示出根據(jù)本發(fā)明示例性實(shí)施例針對疊前三維數(shù)據(jù)來抑制相干噪聲的示圖。
具體實(shí)施例方式現(xiàn)將詳細(xì)描述本發(fā)明的示例性實(shí)施例，所述實(shí)施例的示例在附圖中示出，其中，相同的標(biāo)號指示相同的部分。以下將通過參照附圖來說明所述實(shí)施例，以便解釋本發(fā)明。圖1是示出根據(jù)本發(fā)明示例性實(shí)施例的相干噪聲抑制方法的流程圖。如圖1所示， 根據(jù)本發(fā)明示例性實(shí)施例的相干噪聲抑制方法通過基于相干噪聲的視速度及其方向來確定其視傾角，并按照視傾角將估計(jì)的相干噪聲數(shù)據(jù)轉(zhuǎn)換到分?jǐn)?shù)域中，從而在分?jǐn)?shù)域執(zhí)行相干噪聲預(yù)測。上述方法能夠在較低信噪比的地震資料中有效地預(yù)測出分?jǐn)?shù)域的相干噪聲， 在將分?jǐn)?shù)域的相干噪聲變換到時空域之后，從原始的地震資料中減去變換到時空域的預(yù)測相干噪聲，以起到抑制相干噪聲的作用。具體說來，在步驟S100，在監(jiān)測的地震數(shù)據(jù)中估計(jì)相干噪聲的頻帶和視速度范圍， 并提取所述頻帶內(nèi)的地震數(shù)據(jù)。圖2A示出根據(jù)本發(fā)明示例性實(shí)施例監(jiān)測到的地震數(shù)據(jù)的示意圖。如圖2A所示， 在由多個監(jiān)測儀器監(jiān)測到的具有多個地震道的地震數(shù)據(jù)中，必然存在著一定數(shù)量的相干干擾。相干干擾如地表折射多次波、面波和聲波及直達(dá)波等，由于傳播機(jī)理各不相同，使得它們相互之間以及與有效反射波之間存在諸多方面的差異，其顯著特征之一就是頻譜的衰減規(guī)律有較大不同，也就是說，不同的干擾波，其優(yōu)勢頻帶各不相同。因此，針對不同的線性干擾，在各自的優(yōu)勢頻帶范圍內(nèi)進(jìn)行識別，不僅可以更準(zhǔn)確地識別出線性干擾波，而且可以保證其它頻段內(nèi)的有效信號不受影響，使處理結(jié)果具有更高的保真度，而相干噪聲一般集中于地震數(shù)據(jù)所有頻率中的低頻部分。因此，作為示例，可人工識別地震數(shù)據(jù)中相干噪聲的頻帶，并同時估計(jì)出相干噪聲的視速度范圍。本領(lǐng)域技術(shù)人員應(yīng)理解本發(fā)明并不受限于任何具體的用于估計(jì)相干噪聲的頻帶和視速度范圍的方式，任何能夠有助于從監(jiān)測到的地震數(shù)據(jù)中分離出相干噪聲并估計(jì)出相干噪聲的大致視速度范圍的手段均可等同地應(yīng)用于本發(fā)明。在估計(jì)出相干噪聲的頻帶之后，構(gòu)建相應(yīng)的帶通濾波器，應(yīng)用所述帶通濾波器對地震數(shù)據(jù)進(jìn)行濾波，以得到濾波后的地震數(shù)據(jù)。圖2B示出根據(jù)本發(fā)明示例性實(shí)施例過濾出的地震數(shù)據(jù)的示意圖。如圖2B所示，濾波后的地震數(shù)據(jù)中主要包括相干噪聲，但同時也含有一定數(shù)量的有效波。為了進(jìn)一步得到其中的有效波，有必要對濾波后的地震數(shù)據(jù)進(jìn)行進(jìn)一步的處理。
然后，在步驟S200，在視速度范圍內(nèi)，針對提取的所述頻帶內(nèi)的地震數(shù)據(jù)來確定與視速度對應(yīng)的視速度方向。作為示例，在估計(jì)出的視速度范圍內(nèi)，可根據(jù)實(shí)際需要選取一個或多個視速度 (例如，按照預(yù)定的速度間隔來選取一個或多個典型的視速度)來對過濾出的地震數(shù)據(jù)依次進(jìn)行掃描，并由此確定與所述一個或多個視速度對應(yīng)的視速度方向。在利用每個視速度進(jìn)行掃描時，可沿著各個方向?qū)λ鲞^濾出的地震數(shù)據(jù)進(jìn)行采樣，這樣，每個方向上可得到一組采樣數(shù)據(jù)，然后對每個方向上的一組采樣數(shù)據(jù)相加求和，根據(jù)相干噪聲的特性可知絕對值最大的求和值所對應(yīng)的方向則為相應(yīng)視速度的視速度方向，也是相應(yīng)視速度下的相干噪聲的視速度方向。在得到了相干噪聲的視速度和視速度方向之后，可在步驟S300基于所述視速度和對應(yīng)的視速度方向來確定視傾角，其中，視傾角的方向即為視速度的方向。關(guān)于如何確定視傾角的大小，作為示例，對于二維疊前數(shù)據(jù)而言，可認(rèn)為其中的相干噪聲滿足以下兩個條件(1)干擾波沿同相軸呈現(xiàn)線性關(guān)系；(2)干擾噪聲的視速度不同于有效波的視速度。因此，在這種情況下，可基于道間距和視速度來確定視傾角的大小。作為優(yōu)選而非限定性的計(jì)算方式，可根據(jù)下面的等式來確定視傾角的大小a a = (1000 X地震道間距)/視速度。其中，地震道間距的單位為m，視速度的單位為m/s。作為另一示例，對于三維疊前數(shù)據(jù)而言，由于在進(jìn)行三維數(shù)據(jù)的野外采集時，炮點(diǎn)不在檢波線上，因此，地震道間距是變化的，以致在較遠(yuǎn)的偏移距上，面波和直達(dá)波不呈線性變化。在這種情況下，可基于偏移距和視速度來確定視傾角的大小。在步驟S300確定視傾角之后，本發(fā)明的示例性實(shí)施例可在步驟S400根據(jù)視傾角對提取的所述頻帶內(nèi)的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換，從而得到分?jǐn)?shù)域的相應(yīng)數(shù)據(jù)。圖 2C示出根據(jù)本發(fā)明示例性實(shí)施例對濾波后的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換的示意圖。如圖2C所示，當(dāng)根據(jù)在步驟S300確定的視傾角對在步驟SlOO濾波后的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換(將所述地震數(shù)據(jù)與視傾角下的分?jǐn)?shù)傅里葉變換核矩陣相乘以得到分?jǐn)?shù)域下的對應(yīng)數(shù)據(jù))時，可看出所述變換處理相當(dāng)于將地震數(shù)據(jù)在時空域中的坐標(biāo)軸沿視傾角的方向旋轉(zhuǎn)相應(yīng)的角度a。通過分?jǐn)?shù)階傅里葉變換，可以看出分?jǐn)?shù)域中的相干噪聲在視傾角的方向上呈現(xiàn)良好的相關(guān)性，而與有效波的相關(guān)性則比較差。因此，本發(fā)明的示例性實(shí)施例在步驟S500對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測，以得到預(yù)測的相干噪聲。根據(jù)分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)的特性(即，相干噪聲呈現(xiàn)良好的相關(guān)性，而又與有效波不相關(guān)，其中，隨機(jī)噪聲不可預(yù)測)，可通過執(zhí)行相干預(yù)測來得到相干噪聲。圖2D示出根據(jù)本發(fā)明示例性實(shí)施例對濾波后的地震數(shù)據(jù)執(zhí)行相干預(yù)測的示意圖。以下，首先結(jié)合最簡單的情況來說明相干預(yù)測的基本概念假設(shè)對于一組信道時差為At1的相鄰地震道，其中，第一地震道上的地震波數(shù)據(jù)可表示為W1 (t)，則在第一地震道上經(jīng)過時間、的反射波為W1 (t-、)，對其進(jìn)行傅里葉變換， 可得到，對于某一個頻率f，則相應(yīng)地構(gòu)成了一個復(fù)數(shù)序列Wi,W, (f)e-i2^ , 體‘、,W 人 f]e_n7tfh...
其中是對反射波W1(It1)進(jìn)行傅里葉變換的輸出結(jié)果，f為頻率，、為時間，對應(yīng)的離散Z變換形式為
權(quán)利要求
1.一種用于在監(jiān)測的地震數(shù)據(jù)中抑制相干噪聲的方法，所述方法包括在監(jiān)測的地震數(shù)據(jù)中估計(jì)相干噪聲的頻帶和視速度范圍，并提取所述頻帶內(nèi)的地震數(shù)據(jù)；在視速度范圍內(nèi)，針對提取的所述頻帶內(nèi)的地震數(shù)據(jù)來確定與視速度對應(yīng)的視速度方向；基于所述視速度和對應(yīng)的視速度方向來確定視傾角； 根據(jù)視傾角對提取的所述頻帶內(nèi)的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換； 對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測，以得到預(yù)測的相干噪聲；以及從監(jiān)測的地震數(shù)據(jù)中消除預(yù)測的相干噪聲。
2.如權(quán)利要求1所述的方法，其中，估計(jì)相干噪聲的頻帶和視速度范圍的步驟包括人工估計(jì)監(jiān)測的地震數(shù)據(jù)中相干噪聲的頻帶，并估計(jì)出相干噪聲的視速度范圍。
3.如權(quán)利要求1所述的方法，其中，提取所述頻帶內(nèi)的地震數(shù)據(jù)的步驟包括在估計(jì)出相干噪聲的頻帶之后，構(gòu)建相應(yīng)的帶通濾波器，應(yīng)用所述帶通濾波器對監(jiān)測的地震數(shù)據(jù)進(jìn)行濾波，以得到濾波后的所述頻帶內(nèi)的地震數(shù)據(jù)。
4.如權(quán)利要求1所述的方法，其中，確定視速度方向的步驟包括在估計(jì)出的視速度范圍內(nèi)，選取一個或多個視速度來對提取的所述頻帶內(nèi)的地震數(shù)據(jù)依次進(jìn)行掃描，并由此確定與所述一個或多個視速度對應(yīng)的視速度方向。
5.如權(quán)利要求4所述的方法，其中，當(dāng)對提取的所述頻帶內(nèi)的地震數(shù)據(jù)依次進(jìn)行掃描時，沿著各個方向?qū)λ鲱l帶內(nèi)的地震數(shù)據(jù)進(jìn)行采樣，從而在每個方向上得到一組采樣數(shù)據(jù)，然后對每個方向上的一組采樣數(shù)據(jù)相加求和，將絕對值最大的求和值所對應(yīng)的方向確定為相應(yīng)視速度的視速度方向。
6.如權(quán)利要求1所述的方法，其中，將視速度的方向確定為視傾角的方向。
7.如權(quán)利要求6所述的方法，其中，確定視傾角的步驟包括對于二維疊前地震數(shù)據(jù)， 基于視速度和監(jiān)測的地震數(shù)據(jù)中的地震道間距來確定視傾角的大小。
8.如權(quán)利要求7所述的方法，其中，視傾角的大小=(1000X地震道間距)/視速度，其中，地震道間距的單位為米，視速度的單位為米/秒。
9.如權(quán)利要求6所述的方法，其中，確定視傾角的步驟包括對于三維疊前地震數(shù)據(jù)， 基于視速度和監(jiān)測的地震數(shù)據(jù)中的偏移距來確定視傾角的大小。
10.如權(quán)利要求1所述的方法，其中，根據(jù)視傾角對提取的所述頻帶內(nèi)的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換的步驟包括將提取的所述頻帶內(nèi)的地震數(shù)據(jù)與視傾角下的分?jǐn)?shù)傅里葉變換核矩陣相乘，以得到分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)。
11.如權(quán)利要求1所述的方法，其中，對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測的步驟包括采用F-X域?yàn)V波算法對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測。
12.如權(quán)利要求1所述的方法，其中，對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測的步驟包括利用在所述頻帶外的頻段求取的預(yù)測因子對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測。
13.如權(quán)利要求12所述的方法，其中，所述頻帶外的頻段是指相干噪聲較少的頻段。
14.如權(quán)利要求13所述的方法，其中，所述相干噪聲較少的頻段為中頻段。
15.如權(quán)利要求14所述的方法，其中，利用在所述頻帶外的頻段求取的預(yù)測因子對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測的步驟包括用求取的預(yù)測因子與所述頻帶內(nèi)的地震數(shù)據(jù)進(jìn)行褶積，從而得到相應(yīng)的相干預(yù)測濾波結(jié)果。
16.如權(quán)利要求15所述的方法，其中，用求取的預(yù)測因子與所述頻帶內(nèi)的地震數(shù)據(jù)進(jìn)行褶積的步驟包括根據(jù)中頻段與相干噪聲的所述頻帶之間的頻率關(guān)系，用求取的預(yù)測因子與所述頻帶內(nèi)的地震數(shù)據(jù)進(jìn)行褶積。
17.如權(quán)利要求1所述的方法，其中，從監(jiān)測的地震數(shù)據(jù)中消除預(yù)測的相干噪聲的步驟包括將預(yù)測的相干噪聲進(jìn)行分?jǐn)?shù)階傅里葉逆變換，以得到時空域中的相干噪聲，然后從監(jiān)測到的地震數(shù)據(jù)中減去時空域中的相干噪聲，從而得到抑制相干噪聲后的有效地震數(shù)據(jù)。
全文摘要
提供一種基于分?jǐn)?shù)階傅里葉變換的相干噪聲抑制方法，所述方法包括在監(jiān)測的地震數(shù)據(jù)中估計(jì)相干噪聲的頻帶和視速度范圍，并提取所述頻帶內(nèi)的地震數(shù)據(jù)；在視速度范圍內(nèi)，針對提取的所述頻帶內(nèi)的地震數(shù)據(jù)來確定與視速度對應(yīng)的視速度方向；基于所述視速度和對應(yīng)的視速度方向來確定視傾角；根據(jù)視傾角對提取的所述頻帶內(nèi)的地震數(shù)據(jù)執(zhí)行分?jǐn)?shù)階傅里葉變換；對分?jǐn)?shù)階傅里葉變換后的地震數(shù)據(jù)執(zhí)行相干預(yù)測，以得到預(yù)測的相干噪聲；以及從監(jiān)測的地震數(shù)據(jù)中消除預(yù)測的相干噪聲。
文檔編號G01V1/36GK102323618SQ201110130279
公開日2012年1月18日 申請日期2011年5月19日 優(yōu)先權(quán)日2011年5月19日
發(fā)明者何光明, 呂文彪, 張華 , 曹中林, 王玉華, 羅紅明, 陳愛萍 申請人:中國石油集團(tuán)川慶鉆探工程有限公司