一種新的石油地震勘探數(shù)據(jù)去噪方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于地震數(shù)據(jù)勘探技術(shù)領(lǐng)域，涉及一種新的石油地震勘探數(shù)據(jù)去噪方法。
【背景技術(shù)】
[0002] 石油勘探地震數(shù)據(jù)去噪，是石油地震勘探數(shù)據(jù)處理中非常重要的一個環(huán)節(jié)。去噪 方法一般是根據(jù)不同的變換域中信號和噪聲的差異進行。常用的變換域有f_k域，f-x域， τ-P域和頻域等。如果噪聲和信號在所述變換域中差異不明顯，或者無差異，則去噪無法進 行。例如，經(jīng)過常用的低通濾波器（可以選用Butterworth濾波器實現(xiàn)）后，噪聲和信號在 頻域?qū)o法再進一步分離，但數(shù)據(jù)中仍然存在噪聲，只不過此時噪聲和信號（數(shù)據(jù)）具有相 同的振幅譜。

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

[0003] 本發(fā)明的目的在于提供一種新的石油地震勘探數(shù)據(jù)去噪方法，解決了現(xiàn)有技術(shù)常 用的低通濾波器噪聲和信號在頻域?qū)o法再進一步分離，數(shù)據(jù)中仍然存在噪聲的問題。
[0004] 本發(fā)明所采用的技術(shù)方案是按照以下步驟進行：
[0005] 步驟1 :首先采集得到的石油地震勘探數(shù)據(jù)，經(jīng)過本領(lǐng)域常規(guī)的預處理；
[0006] 步驟2 :然后利用快速傅立葉變換分析地震數(shù)據(jù)的頻譜，確定有用信號的頻譜范 圍，進而明確低通濾波器的截止頻率，利用此截止頻率對地震數(shù)據(jù)進行低通濾波，濾除地震 數(shù)據(jù)的高頻噪聲；
[0007] 步驟3 :對于濾除高頻噪聲之后的地震數(shù)據(jù)，送入分數(shù)傅立葉變換（FrFT)域進行 變換，在進行FrFT變換時，首先按照一定的時間間隔掃描將其分為不同的階數(shù)，在分別計 算不同階數(shù)所對應的FrFT變換；
[0008] 步驟4 :根據(jù)掃描所得到的FrFT域變換，選擇出信號最為集中的變換數(shù)據(jù)和其所 對應的階數(shù)Btl;
[0009] 步驟5 :將步驟3中進行FrFT域變換后的數(shù)據(jù)，再次進行反變換，反變換的階數(shù) 為-?，經(jīng)此變換之后的數(shù)據(jù)即為所求數(shù)據(jù)。
[0010] 進一步，所述步驟2中的低通濾波方式采用巴特沃斯濾波器實現(xiàn)。
[0011] 進一步，所述步驟3中FrFT變換的階數(shù)范圍為[0,4]，時間間隔小于0.5秒。
[0012] 本發(fā)明的有益效果是石油勘探數(shù)據(jù)中的噪聲的進一步濾除。
【附圖說明】
[0013] 圖1是本發(fā)明方法流程示意圖；
[0014] 圖2是Ricker子波的FrFT變換圖。
【具體實施方式】
[0015] 下面結(jié)合附圖和【具體實施方式】對本發(fā)明進行詳細說明。
[0016] 如圖1所示，本發(fā)明方法按照以下步驟進行：
[0017] 步驟1 :首先在野外采集得到的石油地震勘探數(shù)據(jù)，經(jīng)過本領(lǐng)域常規(guī)的預處理（預 處理可采用數(shù)據(jù)解編、道編輯等處理方法）；
[0018] 步驟2 :利用FFT (快速傅立葉變換）分析地震數(shù)據(jù)的頻譜，確定有用信號的頻譜 范圍，進而明確低通濾波器的截止頻率，利用此截止頻率對地震數(shù)據(jù)進行低通濾波（可以 利用巴特沃斯濾波器實現(xiàn)，也可以選擇其他的實現(xiàn)形式），濾除地震數(shù)據(jù)的高頻噪聲；
[0019] 步驟3 :對于濾除高頻噪聲之后的地震數(shù)據(jù)，送入分數(shù)傅立葉變換（FrFT)域進行 變換，F(xiàn)rFT變換的階數(shù)范圍為[0,4]。在進行FrFT變換時，首先進行變換掃描，即在[0,4] 范圍內(nèi)按照一定的時間間隔（時間間隔優(yōu)選小于0.5秒）分別計算不同階數(shù)所對應的FrFT 變換，其計算公式為：
【主權(quán)項】
1. 一種新的石油地震勘探數(shù)據(jù)去噪方法，其特征在于，該新的石油地震勘探數(shù)據(jù)去噪 方法包括以下步驟： 步驟一：首先采集得到的石油地震勘探數(shù)據(jù)，預處理； 步驟二：然后利用快速傅立葉變換分析地震數(shù)據(jù)的頻譜，確定有用信號的頻譜范圍，進 而明確低通濾波器的截止頻率，利用此截止頻率對地震數(shù)據(jù)進行低通濾波，濾除地震數(shù)據(jù) 的尚頻噪聲； 步驟三：對于濾除高頻噪聲之后的地震數(shù)據(jù)，送入分數(shù)傅立葉變換FrFT域進行變換， 在進行FrFT變換時，首先按照時間間隔掃描將其分為不同的階數(shù)，在分別計算不同階數(shù)所 對應的FrFT變換； 步驟四：根據(jù)掃描所得到的FrFT域變換，選擇出信號最為集中的變換數(shù)據(jù)和其所對應 的階數(shù)aQ; 步驟五：將步驟三中進行FrFT域變換后的數(shù)據(jù)，再次進行反變換，反變換的階數(shù) 為-%，經(jīng)此變換之后的數(shù)據(jù)即為所求數(shù)據(jù)。
2. 如權(quán)利要求1所述新的石油地震勘探數(shù)據(jù)去噪方法，其特征在于：所述步驟二中的 低通濾波方式采用巴特沃斯濾波器實現(xiàn)。
3. 如權(quán)利要求1所述新的石油地震勘探數(shù)據(jù)去噪方法，其特征在于：步驟三中FrFT變 換的階數(shù)范圍為[〇,4]，時間間隔小于0.5秒。
【專利摘要】本發(fā)明公開了一種新的石油地震勘探數(shù)據(jù)去噪方法，采集得到的石油地震勘探數(shù)據(jù)，利用快速傅立葉變換分析地震數(shù)據(jù)的頻譜，確定有用信號的頻譜范圍，進而明確低通濾波器的截止頻率，利用此截止頻率對地震數(shù)據(jù)進行低通濾波，濾除地震數(shù)據(jù)的高頻噪聲；對于濾除高頻噪聲之后的地震數(shù)據(jù)，送入FrFT域進行變換，在進行FrFT變換時，根據(jù)掃描所得到的FrFT域變換，選擇出信號最為集中的變換數(shù)據(jù)和其所對應的階數(shù)a0；將FrFT域變換后的數(shù)據(jù)，再次進行反變換，反變換的階數(shù)為-a0，經(jīng)此變換之后的數(shù)據(jù)即為所求數(shù)據(jù)。發(fā)明的有益效果是石油勘探數(shù)據(jù)中的噪聲的進一步濾除。
【IPC分類】G01V1-36
【公開號】CN104597502
【申請?zhí)枴緾N201410735187
【發(fā)明人】翟明岳
【申請人】翟明岳
【公開日】2015年5月6日
【申請日】2014年12月8日