自適應(yīng)低頻異常振幅壓制的方法
【專利摘要】一種自適應(yīng)低頻異常振幅壓制的方法，該方法包括：得到頻率域的地震數(shù)據(jù)；從頻率域的地震數(shù)據(jù)分離出低頻段的地震數(shù)據(jù)；從分離出的低頻段的地震數(shù)據(jù)得到時(shí)間域的低頻地震數(shù)據(jù)；從時(shí)間域的低頻地震數(shù)據(jù)提取異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)；根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值；根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)；從單炮地震數(shù)據(jù)中減去低頻異常振幅分布區(qū)數(shù)據(jù)與低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)之差，得到低頻異常振幅壓制后的地震數(shù)據(jù)。根據(jù)本發(fā)明，可以自動去除采集的地震數(shù)據(jù)中的低頻異常振幅干擾，同時(shí)提高地震數(shù)據(jù)處理效率。
【專利說明】自適應(yīng)低頻異常振幅壓制的方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于石油地震資料勘探領(lǐng)域，更具體地講，涉及一種自動、有效地去除采集的地震數(shù)據(jù)中包含的低頻異常振幅干擾，同時(shí)提高地震數(shù)據(jù)處理效率的自適應(yīng)低頻異常振幅壓制方法。
【背景技術(shù)】
[0002]低頻異常振幅是石油地震資料勘探中廣泛存在的一種干擾波，它主要集中在地震記錄的近偏移距附近，具有低頻、低速、衰減較慢的特點(diǎn)，強(qiáng)低頻異常振幅的存在，將嚴(yán)重地影響反褶積和疊加效果，對保幅處理尤為不利。在地震資料處理中低頻異常振幅壓制是提高地震資料信噪比和分辨率的重要步驟，因此在地震資料處理中如何去除或提取低頻異常振幅就成了目前地震勘探數(shù)據(jù)處理的主要問題之一。
[0003]目前在地震資料處理過程中大多采用高通濾波的方法來壓制低頻異常振幅，主要是利用低頻異常振幅和反射波的頻率差異及速度差異進(jìn)行濾波，如一維濾波和頻率波數(shù)濾波，但因低頻異常振幅的頻帶往往和反射波中的低頻帶攪在一起，即便是頻率波數(shù)域?yàn)V波也不能只對與有效信號相同頻帶的噪聲進(jìn)行壓制，其實(shí)在濾波過程的同時(shí)也對有效波的低頻成分造成了損害，而且斜坡長的濾波響應(yīng)在反褶積或譜白化處理后，低頻異常振幅往往又被突出出來。因此，需要一種能夠有效的去除采集的地震數(shù)據(jù)中包含的低頻異常振幅干擾，同時(shí)提高地震數(shù)據(jù)處理效率的地震數(shù)據(jù)處理方法。

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

[0004]本發(fā)明的目的在于提供一種能夠自動去除采集的地震數(shù)據(jù)中包含的低頻異常振幅干擾，同時(shí)提高地震數(shù)據(jù)處理效率的自適應(yīng)低頻異常振幅壓制方法，從而滿足日漸興起的高分辨率勘探的精度要求。
[0005]本發(fā)明的一方面提供一種自適應(yīng)低頻異常振幅壓制的方法，包括:對單炮地震數(shù)據(jù)的每道進(jìn)行傅立葉正變換，得到頻率域的地震數(shù)據(jù)；從頻率域的地震數(shù)據(jù)分離出低頻段的地震數(shù)據(jù)；將分離出的低頻段的地震數(shù)據(jù)進(jìn)行傅立葉反變換，得到時(shí)間域的低頻地震數(shù)據(jù)；從時(shí)間域的低頻地震數(shù)據(jù)提取低頻異常振幅分布區(qū)數(shù)據(jù)、淺層折射區(qū)數(shù)據(jù)以及無干擾分布區(qū)數(shù)據(jù)；根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值；根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)；從單炮地震數(shù)據(jù)中減去低頻異常振幅分布區(qū)數(shù)據(jù)與低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)之差，得到低頻異常振幅壓制后的地震數(shù)據(jù)。
[0006]可選地，所述低頻段的頻率范圍為0-20HZ。
[0007]可選地，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值的步驟包括:對提取的無干擾分布區(qū)數(shù)據(jù)進(jìn)行傅立葉正變換，得到頻率域的無干擾分布區(qū)數(shù)據(jù)；計(jì)算頻率域的無干擾分布區(qū)數(shù)據(jù)的低頻異常振幅閥值。
[0008]可選地，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值的步驟包括:將無干擾分布區(qū)數(shù)據(jù)沿著時(shí)間的方向劃分為多個(gè)時(shí)間窗口 ；對每個(gè)時(shí)間窗口的數(shù)據(jù)進(jìn)行傅立葉正變換，得到相應(yīng)的頻率域的地震數(shù)據(jù)；將得到的每個(gè)時(shí)間窗口的頻率域的地震數(shù)據(jù)劃分為多個(gè)頻段；計(jì)算每個(gè)時(shí)間窗口的每個(gè)頻段的低頻異常振幅閥值。
[0009]可選地，根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)的步驟包括:對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，將提取的低頻異常振幅分布區(qū)數(shù)據(jù)旋轉(zhuǎn)變換到分?jǐn)?shù)域的一個(gè)時(shí)頻平面；從旋轉(zhuǎn)變換得到的數(shù)據(jù)查找低頻段的振幅大于低頻異常振幅閥值的區(qū)域；在查找到的區(qū)域填充預(yù)定值，從而獲取時(shí)頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)；對獲取的時(shí)頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉反變換，得到時(shí)間域的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)。
[0010]可選地，對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，將提取的低頻異常振幅分布區(qū)數(shù)據(jù)旋轉(zhuǎn)變換到分?jǐn)?shù)域的一個(gè)時(shí)頻平面的步驟包括:將低頻異常振幅分布區(qū)數(shù)據(jù)沿著時(shí)間的方向劃分為多個(gè)時(shí)間窗口 ；對每個(gè)時(shí)間窗口的數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，得到相應(yīng)的分?jǐn)?shù)域的地震數(shù)據(jù)。
[0011]可選地，從旋轉(zhuǎn)變換得到的數(shù)據(jù)查找低頻段的振幅大于低頻異常振幅閥值的區(qū)域的步驟包括:將每個(gè)時(shí)間窗口的分?jǐn)?shù)域的地震數(shù)據(jù)劃分為多個(gè)頻段；將每個(gè)時(shí)間窗口的每個(gè)頻段的振幅值與每個(gè)時(shí)間窗口的每個(gè)頻段的低頻異常振幅閥值進(jìn)行比較，在每個(gè)時(shí)間窗口的每個(gè)頻段中查找振幅值大于對應(yīng)的低頻異常振幅閥值的區(qū)域。
[0012]可選地，所述預(yù)定值為零。
[0013]可選地，在二維分?jǐn)?shù)階傅立葉變換中的旋轉(zhuǎn)角度為P π/2，其中，P為二維分?jǐn)?shù)階傅立葉變換的階數(shù)。
[0014]可選地，低頻異常振幅分布區(qū)數(shù)據(jù)是時(shí)間域的低頻地震數(shù)據(jù)中低頻異常振幅集中分布的區(qū)域中的地震數(shù)據(jù)，無干擾分布區(qū)數(shù)據(jù)是時(shí)間域的低頻地震數(shù)據(jù)中不存在低頻異常振幅干擾的區(qū)域中的地震數(shù)據(jù)。
[0015]根據(jù)本發(fā)明的自適應(yīng)低頻異常振幅壓制的方法，可以自動去除采集的地震數(shù)據(jù)中的低頻異常振幅干擾，同時(shí)提高地震數(shù)據(jù)處理效率。將在接下來的描述中部分闡述本發(fā)明總體構(gòu)思另外的方面和/或優(yōu)點(diǎn)，還有一部分通過描述將是清楚的，或者可以經(jīng)過本發(fā)明總體構(gòu)思的實(shí)施而得知。
【專利附圖】

【附圖說明】
[0016]通過下面結(jié)合附圖進(jìn)行的詳細(xì)描述，本發(fā)明的上述和其它目的、特點(diǎn)和優(yōu)點(diǎn)將會變得更加清楚，其中:
[0017]圖1是示出根據(jù)本發(fā)明實(shí)施例的自適應(yīng)低頻異常振幅壓制的方法的流程圖；
[0018]圖2是示出從地震數(shù)據(jù)提取低頻異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)的示意圖。
【具體實(shí)施方式】
[0019]提供參照附圖的以下描述以幫助對由權(quán)利要求及其等同物限定的本發(fā)明的實(shí)施例的全面理解。包括各種特定細(xì)節(jié)以幫助理解，但這些細(xì)節(jié)僅被視為是示例性的。因此，本領(lǐng)域的普通技術(shù)人員將認(rèn)識到在不脫離本發(fā)明的范圍和精神的情況下，可對描述于此的實(shí)施例進(jìn)行各種改變和修改。此外，為了清楚和簡潔，省略對公知的功能和結(jié)構(gòu)的描述。
[0020]本發(fā)明的自適應(yīng)低頻異常振幅壓制的方法包括:對單炮地震數(shù)據(jù)的每道進(jìn)行傅立葉正變換，得到頻率域的地震數(shù)據(jù)；從頻率域的地震數(shù)據(jù)分離出低頻段的地震數(shù)據(jù)；將分離出的低頻段的地震數(shù)據(jù)進(jìn)行傅立葉反變換，得到時(shí)間域的低頻地震數(shù)據(jù)；從時(shí)間域的低頻地震數(shù)據(jù)提取低頻異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)；根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值；根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)；從單炮地震數(shù)據(jù)中減去低頻異常振幅分布區(qū)數(shù)據(jù)與低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)之差，得到低頻異常振幅壓制后的地震數(shù)據(jù)。
[0021]圖1是示出根據(jù)本發(fā)明的實(shí)施例的自適應(yīng)低頻異常振幅壓制的方法的流程圖。
[0022]在步驟101中，對單炮地震數(shù)據(jù)的每道進(jìn)行傅立葉正變換，得到頻率域的地震數(shù)據(jù)。將地震數(shù)據(jù)由時(shí)間域變換到頻率域的主要目的是能夠?qū)⒌卣饠?shù)據(jù)按照頻率進(jìn)行劃分。
[0023]每道地震數(shù)據(jù)記為A (Xi，t)，Xi為第i道地震數(shù)據(jù)的地震信號，i表示地震數(shù)據(jù)的道數(shù)，i e [1，m]，且m為大于O的整數(shù)，t表示時(shí)間，對每道地震數(shù)據(jù)進(jìn)行傅立葉變換可表示如下，
[0024]Y(xi； o)=fft(A(Xi, t)) (I)
[0025]式(I)中，fft(A(Xi,t))表示對 A (Xi, t)進(jìn)行傅立葉變換，Y (Xi, ω)是 A(x”t)對應(yīng)的頻率域的地震數(shù)據(jù)，ω為角頻率，co=2 3if，f為頻率。
[0026]在步驟102中，從頻率域的地震數(shù)據(jù)分離出低頻段的地震數(shù)據(jù)。優(yōu)選地，低頻段為低頻異常振幅主要集中的頻段。例如，低頻段的頻率范圍為0-20HZ。需要壓制地震數(shù)據(jù)中的低頻異常振幅干擾只需針對地震數(shù)據(jù)的低頻成分即可。
[0027]在步驟103中，將分離出的低頻段的地震數(shù)據(jù)進(jìn)行傅立葉反變換，得到時(shí)間域的低頻地震數(shù)據(jù)。將得到的時(shí)間域的每道低頻地震數(shù)據(jù)記為B(Xi，t)。換言之，B(Xi，t)只是地震數(shù)據(jù)A(Xi，t)中的一部分，即低頻異常振幅存在的部分。
[0028]在步驟104中，從時(shí)間域的低頻地震數(shù)據(jù)提取低頻異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)。低頻異常振幅分布區(qū)數(shù)據(jù)是時(shí)間域的低頻地震數(shù)據(jù)中低頻異常振幅集中分布的區(qū)域中的地震數(shù)據(jù)，無干擾分布區(qū)數(shù)據(jù)是時(shí)間域的低頻地震數(shù)據(jù)中不存在低頻異常振幅干擾的區(qū)域中的地震數(shù)據(jù)。
[0029]時(shí)間域的低頻地震數(shù)據(jù)可被分為低頻異常振幅分布區(qū)數(shù)據(jù)、無干擾分布區(qū)數(shù)據(jù)和折射波區(qū)數(shù)據(jù)，在自適應(yīng)低頻異常振幅壓制方法中，只有低頻異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)參與壓制低頻異常振幅的運(yùn)算，折射波分布區(qū)數(shù)據(jù)并不參與運(yùn)算。
[0030]圖2是示出從時(shí)間域的低頻地震數(shù)據(jù)提取低頻異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)的示意圖。
[0031]如圖2所示，橫坐標(biāo)為道號，縱坐標(biāo)為時(shí)間。時(shí)間域的低頻地震數(shù)據(jù)可被劃分為折射波區(qū)域、無干擾區(qū)域和低頻異常振幅分布區(qū)域，圖中從上到下的區(qū)域依次為折射波區(qū)域
I1、無干擾區(qū)域I2和低頻異常振幅分布區(qū)域13。
[0032]低頻異常振幅分布區(qū)域I3的地震數(shù)據(jù)可稱為低頻異常振幅分布區(qū)數(shù)據(jù)，記為C(Xi, t)，無干擾區(qū)域I2的地震數(shù)據(jù)可稱為無干擾分布區(qū)數(shù)據(jù)，記為C1 (Xi, t)。換言之，從B (Xi, t)中提取出低頻異常振幅分布區(qū)數(shù)據(jù)C(Xi，t)和無干擾分布區(qū)數(shù)據(jù)C1Ui, t)。
[0033]折射波區(qū)域I1的上邊界線L1為初至波線(地震發(fā)生后，地震觀測點(diǎn)最先接收到的波稱為初至波)。為此，首先得到檢測單炮地震數(shù)據(jù)的地區(qū)的替換速度，再將單炮地震數(shù)據(jù)的每道的偏移距除以替換速度，從而針對每道計(jì)算得到一個(gè)對應(yīng)的時(shí)間值Tli (i表示道號)，這樣可以得到多個(gè)點(diǎn)(i，Tli),連接這些點(diǎn)獲得初至線。
[0034]可通過如下方式確定折射波區(qū)域I1與無干擾區(qū)域I2的分界線L2 (即，折射波區(qū)域I1的下邊界線，無干擾區(qū)域I2的上邊界線):根據(jù)單炮地震數(shù)據(jù)估算最小的折射波速度，再將單炮地震數(shù)據(jù)中每道的偏移距除以折射波速度，從而針對每道計(jì)算得到一個(gè)對應(yīng)的時(shí)間值T2i (i表示道號)，這樣可以得到多個(gè)點(diǎn)(i，T2i)，連接這些點(diǎn)獲得邊界線L2。
[0035]可通過如下方式確定無干擾區(qū)域I2與低頻異常振幅分布區(qū)域I3的分界線L3 (SP，無干擾區(qū)域I2的下邊界線，低頻異常振幅分布區(qū)域I3的上邊界線):首先根據(jù)單炮地震數(shù)據(jù)估算面波的速度(低頻異常振幅主要集中在近偏移距，在地震記錄上表現(xiàn)形式也就是面波)，再將單炮地震數(shù)據(jù)中每道的偏移距除以面波的速度，從而針對每道計(jì)算得到一個(gè)對應(yīng)的時(shí)間值T3i (i表示道號)，這樣可以得到多個(gè)點(diǎn)(i，T3i)，連接這些點(diǎn)獲得邊界線L3。
[0036]在步驟105中，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值。
[0037]具體地講，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值的步驟為:
[0038](a)對提取的無干擾分布區(qū)數(shù)據(jù)進(jìn)行傅立葉正變換，得到頻率域的無干擾分布區(qū)數(shù)據(jù)。具體地講，首先將無干擾分布區(qū)數(shù)據(jù)C1(Xpt)沿著時(shí)間的方向劃分為多個(gè)窗口(例如，0-500ms，500-1500ms, 1500-2500ms, 2500-4000ms,以此類推)，其次，對每一個(gè)時(shí)間窗口的地震數(shù)據(jù)進(jìn)行傅立葉正變換，得到相應(yīng)的頻率域的地震數(shù)據(jù)。
[0039](b)計(jì)算頻率域的無干擾分布區(qū)數(shù)據(jù)的低頻異常振幅閥值。具體地講，首先將得到的每個(gè)時(shí)間窗口的頻率域的地震數(shù)據(jù)劃分為多個(gè)頻段(例如，l-5Hz，5-10Hz，10-15HZ，15-20HZ,以此類推)，其次，統(tǒng)計(jì)每個(gè)時(shí)間窗口的每個(gè)頻段的數(shù)據(jù)的最大振幅值A(chǔ)m，計(jì)算每個(gè)時(shí)間窗口每個(gè)頻段的低頻異常振幅閥值Ατ。
[0040]A1=FXAm (2)
[0041]其中，F(xiàn)是權(quán)系數(shù)因子，0〈F〈100。
[0042]計(jì)算的低頻異常振幅閥值即是低頻異常振幅振幅的最小門檻值，從而實(shí)現(xiàn)低頻異常振幅的自適應(yīng)壓制，權(quán)系數(shù)因子的目的在于調(diào)節(jié)門檻值以便最佳的壓制低頻異常振幅。
[0043]在步驟106中，根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)。
[0044]具體地講，根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)的步驟為:
[0045](a)對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，將提取的低頻異常振幅分布區(qū)數(shù)據(jù)旋轉(zhuǎn)變換到分?jǐn)?shù)域的一個(gè)時(shí)頻平面。具體地說，首先將低頻異常振幅分布區(qū)數(shù)據(jù)C(Xi，t)沿著時(shí)間的方向劃分為多個(gè)時(shí)窗，對低頻異常振幅分布區(qū)數(shù)據(jù)劃分的時(shí)窗必須與步驟105中無干擾分布區(qū)數(shù)據(jù)劃分的時(shí)窗相同(例如，0-500ms，500-1500ms，1500-2500ms, 2500-4000ms,以此類推)，其次，對每個(gè)時(shí)窗的數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，通過適當(dāng)?shù)男D(zhuǎn)角度，把時(shí)空域上的地震數(shù)據(jù)變換到視速度方向上的分?jǐn)?shù)域中。對低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換可表示如下，
[0046]
【權(quán)利要求】
1.一種自適應(yīng)低頻異常振幅壓制的方法，所述方法包括以下步驟: 對單炮地震數(shù)據(jù)的每道進(jìn)行傅立葉正變換，得到頻率域的地震數(shù)據(jù)； 從頻率域的地震數(shù)據(jù)分離出低頻段的地震數(shù)據(jù)； 將分離出的低頻段的地震數(shù)據(jù)進(jìn)行傅立葉反變換，得到時(shí)間域的低頻地震數(shù)據(jù)； 從時(shí)間域的低頻地震數(shù)據(jù)提取異常振幅分布區(qū)數(shù)據(jù)和無干擾分布區(qū)數(shù)據(jù)； 根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值； 根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)； 從單炮地震數(shù)據(jù)中減去低頻異常振幅分布區(qū)數(shù)據(jù)與低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)之差，得到低頻異常振幅壓制后的地震數(shù)據(jù)。
2.如權(quán)利要求1所述的方法，其中，所述低頻段的頻率范圍為0-20HZ。
3.如權(quán)利要求1所述的方法，其中，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值的步驟包括: (a)對提取的無干擾分布區(qū)數(shù)據(jù)進(jìn)行傅立葉正變換，得到頻率域的無干擾分布區(qū)數(shù)據(jù)； (b)計(jì)算頻率域的無干擾分布區(qū)數(shù)據(jù)的低頻異常振幅閥值。
4.如權(quán)利要求3所述`的方法，其中，根據(jù)提取的無干擾分布區(qū)數(shù)據(jù)獲取低頻異常振幅閥值的步驟包括: 將無干擾分布區(qū)數(shù)據(jù)沿著時(shí)間的方向劃分為多個(gè)時(shí)間窗口； 對每個(gè)時(shí)間窗口的數(shù)據(jù)進(jìn)行傅立葉正變換，得到相應(yīng)的頻率域的地震數(shù)據(jù)； 將得到的每個(gè)時(shí)間窗口的頻率域的地震數(shù)據(jù)劃分為多個(gè)頻段； 計(jì)算每個(gè)時(shí)間窗口的每個(gè)頻段的低頻異常振幅閥值。
5.如權(quán)利要求1所述的方法，其中，根據(jù)低頻異常振幅閥值對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行處理，得到低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)的步驟包括: (a)對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，將提取的低頻異常振幅分布區(qū)數(shù)據(jù)旋轉(zhuǎn)變換到分?jǐn)?shù)域的一個(gè)時(shí)頻平面； (b)從旋轉(zhuǎn)變換得到的數(shù)據(jù)查找低頻段的振幅大于低頻異常振幅閥值的區(qū)域； (C)在查找到的區(qū)域填充預(yù)定值，從而獲取時(shí)頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)； (d)對獲取的時(shí)頻平面上的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉反變換，得到時(shí)間域的低頻異常振幅壓制后的低頻異常振幅分布區(qū)數(shù)據(jù)。
6.如權(quán)利要求5所述的方法，其中，對提取的低頻異常振幅分布區(qū)數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，將提取的低頻異常振幅分布區(qū)數(shù)據(jù)旋轉(zhuǎn)變換到分?jǐn)?shù)域的一個(gè)時(shí)頻平面的步驟包括: 將低頻異常振幅分布區(qū)數(shù)據(jù)沿著時(shí)間的方向劃分為多個(gè)時(shí)間窗口； 對每個(gè)時(shí)間窗口的數(shù)據(jù)進(jìn)行二維分?jǐn)?shù)階傅立葉變換，得到相應(yīng)的分?jǐn)?shù)域的地震數(shù)據(jù)。
7.如權(quán)利要求6所述的方法，其中，從旋轉(zhuǎn)變換得到的數(shù)據(jù)查找低頻段的振幅大于低頻異常振幅閥值的區(qū)域的步驟包括: 將每個(gè)時(shí)間窗口的分?jǐn)?shù)域的地震數(shù)據(jù)劃分為多個(gè)頻段；將每個(gè)時(shí)間窗口的每個(gè)頻段的振幅值與每個(gè)時(shí)間窗口的每個(gè)頻段的低頻異常振幅閥值進(jìn)行比較，在每個(gè)時(shí)間窗口的每個(gè)頻段中查找振幅值大于對應(yīng)的低頻異常振幅閥值的區(qū)域。
8.如權(quán)利要求5所述的方法，其中，所述預(yù)定值為零。
9.如權(quán)利要求5所述的方法，其中，在二維分?jǐn)?shù)階傅立葉變換中的旋轉(zhuǎn)角度為P/2，其中，P為二維分?jǐn)?shù)階傅立葉變換的階數(shù)。
10.如權(quán)利要求1所述的方法，其中，低頻異常振幅分布區(qū)數(shù)據(jù)是時(shí)間域的低頻地震數(shù)據(jù)中低頻異常振幅集中分布的區(qū)域中的地震數(shù)據(jù)，無干擾分布區(qū)數(shù)據(jù)是時(shí)間域的低頻地震數(shù)據(jù)中不存在低頻異常振幅干擾的區(qū)`域中的地震數(shù)據(jù)。
【文檔編號】G01V1/28GK103869361SQ201410083083
【公開日】2014年6月18日 申請日期:2014年3月7日 優(yōu)先權(quán)日:2014年3月7日
【發(fā)明者】曹中林, 何光明, 陳愛萍, 呂文彪, 周強(qiáng), 張華 , 李振 申請人:中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司