一種地震隨機(jī)噪聲衰減方法
【專(zhuān)利摘要】本發(fā)明提供了一種地震隨機(jī)噪聲衰減方法，屬于地震數(shù)據(jù)處理領(lǐng)域。本方法包括：(1)對(duì)含噪數(shù)據(jù)進(jìn)行行和列方向的平移得到新的二維數(shù)組；(2)對(duì)步驟(1)得到的二維數(shù)組進(jìn)行Curvelet變換得到曲波系數(shù)；(3)對(duì)步驟(2)得到的曲波系數(shù)進(jìn)行閾值去噪得到新的曲波系數(shù)；(4)對(duì)步驟(3)得到的新的曲波系數(shù)進(jìn)行Curvelet逆變換得到新的二維數(shù)組；(5)對(duì)步驟(4)得到的新的二維數(shù)組進(jìn)行行和列方向的逆向平移，恢復(fù)數(shù)據(jù)原先的狀態(tài)，即得到一次平移的去噪結(jié)果；(6)判斷是否達(dá)到設(shè)定的循環(huán)處理次數(shù)，如果達(dá)到，則處理完成，轉(zhuǎn)入步驟(7)；如果否，則返回步驟(1)；(7)對(duì)所有去噪結(jié)果進(jìn)行線性平均得到最終去噪結(jié)果。
【專(zhuān)利說(shuō)明】-種地震隨機(jī)噪聲衰減方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于地震數(shù)據(jù)處理領(lǐng)域，具體涉及一種地震隨機(jī)噪聲衰減方法，利用函數(shù) 變換方法進(jìn)行地震數(shù)據(jù)處理，將多分辨、帶通、具有方向性的函數(shù)分析方法和一種循環(huán)平移 的技術(shù)結(jié)合應(yīng)用于地震隨機(jī)噪聲衰減。

【背景技術(shù)】
[0002] 從化rvelet變換提出W來(lái)，國(guó)內(nèi)外的學(xué)者做了大量的工作，主要是按照提出者的 思路進(jìn)行了一些方法的試驗(yàn)；在圖像去噪、圖像增強(qiáng)、圖像融合、圖像恢復(fù)等方面有了一些 初步的應(yīng)用；在地震資料處理方面主要是地震數(shù)據(jù)的去噪方面做了一些探索性的研究。
[0003] 去噪在地震資料處理中占有非常重要的地位，地震記錄的頻帶很寬，在有效波的 頻帶范圍內(nèi)可能包含有多種類(lèi)型的噪聲。地震記錄信噪比的高低，將直接影響地震資料的 可靠性、參數(shù)提取的精度及提高分辨率的效果等?；痳velet變換的多尺度和高度各向異性 的特點(diǎn)，在地震數(shù)據(jù)處理領(lǐng)域有著巨大的發(fā)展?jié)摿Α?br> [0004] 隨著勘探程度的不斷提高與勘探領(lǐng)域的不斷延伸，地震勘探難度越來(lái)越大，復(fù)雜 的地表和地質(zhì)條件對(duì)深層地震信號(hào)產(chǎn)生各種干擾和崎變，在地震記錄上表現(xiàn)為信噪比降 低、噪聲類(lèi)型復(fù)雜、有效信號(hào)弱的特點(diǎn)，該嚴(yán)重制約著地震資料后續(xù)的高分辨率處理、目標(biāo) 準(zhǔn)確成像W及屬性提取和反演。因此，地震噪聲衰減仍舊是地震在資料處理的關(guān)鍵問(wèn)題。 目前，化rvelet變換應(yīng)用于地震去噪的方法主要是闊值去噪，通過(guò)對(duì)含噪數(shù)據(jù)的多尺度 化rvelet變換，利用噪聲和有效信號(hào)在不同方向和尺度的分布，通過(guò)設(shè)定一定的闊值條件， 將噪聲和有效信號(hào)進(jìn)行分離，進(jìn)而達(dá)到去噪效果。
[0005] 在地震勘探中，隨機(jī)噪聲是一種頻帶較寬的干擾波，常規(guī)的去噪方法效果不理想。


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

[0006] 本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題，提供一種地震隨機(jī)噪聲衰減 方法，其是基于化rvelet變換的切cle spinning(循環(huán)平移方法）實(shí)現(xiàn)對(duì)圖像去噪的。
[0007] 本發(fā)明是通過(guò)W下技術(shù)方案實(shí)現(xiàn)的：
[0008] 一種地震隨機(jī)噪聲衰減方法，包括：
[0009] (1)對(duì)含噪數(shù)據(jù)進(jìn)行行和列方向的平移得到新的二維數(shù)組；
[0010] (2)對(duì)步驟（1)得到的二維數(shù)組進(jìn)行化rvelet變換得到曲波系數(shù)；州rvelet變 換是現(xiàn)有的技術(shù)，調(diào)用變換函數(shù)即可（請(qǐng)參考文獻(xiàn)（Candes E，Demanet以Donoho D，et al. Fast discrete curvelet transforms[J]. Multiscale Modeling&Simulation, 2006, 5(3) ：861-899)；
[0011] (3)對(duì)步驟（2)得到的曲波系數(shù)進(jìn)行闊值去噪得到新的曲波系數(shù)；闊值去噪是一 種現(xiàn)有的技術(shù)，對(duì)curvelet變換后得到的曲波系數(shù)采用統(tǒng)一的闊值進(jìn)行取舍，大于闊值的 系數(shù)保留，小于闊值的去除置零，得到新的曲波系數(shù)；
[001引 （4)對(duì)步驟（3)得到的新的曲波系數(shù)進(jìn)行化rvelet逆變換得到新的二維 數(shù)組；州rvelet逆變換也是現(xiàn)有的技術(shù)，調(diào)用函數(shù)即可（請(qǐng)參考文獻(xiàn)（Candes E， Demanet L, Donoho D, et al.Fast discrete curvelet transforms[J]. Multiscale Modelin泌Simulation, 2006,5(3) :861_899);
[0013] (5)對(duì)步驟（4)得到的新的二維數(shù)組進(jìn)行行和列方向的逆向平移，恢復(fù)數(shù)據(jù)原先 的狀態(tài)，即得到一次平移的去噪結(jié)果；該也是cycle spinning技術(shù)的其中一部分；
[0014] (6)判斷是否達(dá)到設(shè)定的循環(huán)處理次數(shù)，如果達(dá)到，則處理完成，轉(zhuǎn)入步驟（7);如 果否，則返回步驟（1);
[0015] 該里是對(duì)原來(lái)的含噪數(shù)據(jù)重復(fù)步驟（1)到步驟化）的處理，該樣做的目的就是讓 噪聲出現(xiàn)在不同的地方進(jìn)行反復(fù)的去噪，得到更好的去噪效果。
[0016] (7)對(duì)所有去噪結(jié)果進(jìn)行線性平均得到最終去噪結(jié)果。
[0017] 所述步驟（1)是該樣實(shí)現(xiàn)的：
[0018] 采用循環(huán)平移方法（即cycle spinning)，即通過(guò)一個(gè)平移函數(shù)對(duì)由含噪數(shù)據(jù)組 成的二維矩陣沿行和列方向進(jìn)行平移，形成新的二維數(shù)組。每次移動(dòng)1行1列。
[0019] 所述步驟巧）中所述一次平移的去噪結(jié)果如下：
[0020]

【權(quán)利要求】
1. 一種地震隨機(jī)噪聲衰減方法，其特征在于；所述方法包括： (1) 對(duì)含噪數(shù)據(jù)進(jìn)行行和列方向的平移得到新的二維數(shù)組； (2) 對(duì)步驟（1)得到的二維數(shù)組進(jìn)行化rvelet變換得到曲波系數(shù)； (3) 對(duì)步驟（2)得到的曲波系數(shù)進(jìn)行闊值去噪得到新的曲波系數(shù)； (4) 對(duì)步驟（3)得到的新的曲波系數(shù)進(jìn)行化rvelet逆變換得到新的二維數(shù)組； (5) 對(duì)步驟（4)得到的新的二維數(shù)組進(jìn)行行和列方向的逆向平移，恢復(fù)數(shù)據(jù)原先的狀 態(tài)，即得到一次平移的去噪結(jié)果； (6) 判斷是否達(dá)到設(shè)定的循環(huán)處理次數(shù)，如果達(dá)到，則處理完成，轉(zhuǎn)入步驟（7);如果 否，則返回步驟（1); (7) 對(duì)所有去噪結(jié)果進(jìn)行線性平均得到最終去噪結(jié)果。
2. 根據(jù)權(quán)利要求1所述的地震隨機(jī)噪聲衰減方法，其特征在于：所述步驟（1)是該樣 實(shí)現(xiàn)的： 采用循環(huán)平移方法，即通過(guò)一個(gè)平移函數(shù)對(duì)由含噪數(shù)據(jù)組成的二維矩陣沿行和列方向 進(jìn)行平移，形成新的二維數(shù)組。
3. 根據(jù)權(quán)利要求2所述的地震隨機(jī)噪聲衰減方法，其特征在于；所述步驟巧）中所述 一次平移的去噪結(jié)果如下：
其中，S為循環(huán)平移算子，下標(biāo)i，j和-i，-j分別為行和列方向上的平移量，C為 化rvelet變換算子，C-i為逆化rvelet變換算子，A為闊值算子。
4. 根據(jù)權(quán)利要求3所述的地震隨機(jī)噪聲衰減方法，其特征在于；所述步驟化）中的所 述循環(huán)處理次數(shù)為1到含噪數(shù)據(jù)二維數(shù)組的行數(shù)X列數(shù)之間，即最小平移次數(shù)為1，最大平 移次數(shù)為行數(shù)X列數(shù)。
5. 根據(jù)權(quán)利要求4所述的地震隨機(jī)噪聲衰減方法，其特征在于：所述步驟（7)是該樣 實(shí)現(xiàn)的： 利用下面的線性平均公式求得最終去噪結(jié)果：
其中，町，馬分別表示行和列方向上的最大平移量。

【文檔編號(hào)】G01V1/36GK104422961SQ201310409633
【公開(kāi)日】2015年3月18日 申請(qǐng)日期:2013年9月10日 優(yōu)先權(quán)日:2013年9月10日
【發(fā)明者】薛詩(shī)桂 申請(qǐng)人:中國(guó)石油化工股份有限公司, 中國(guó)石油化工股份有限公司石油物探技術(shù)研究院