本發(fā)明涉及一種電磁式可控震源并行震源地震勘探時壓制串?dāng)_噪聲的方法，尤其是電磁式可控震源拋物線Radon變換壓制并行震源串?dāng)_噪聲方法。
背景技術(shù)：
：：多震源激發(fā)就是在多個不同位置上激發(fā)多個震源并行工作，因此得到的地震記錄是一種由多個震源產(chǎn)生的混合地震記錄，提高并行震源地震勘探數(shù)據(jù)質(zhì)量，但多震源地震數(shù)據(jù)為并行處理問題。博士論文“電磁式可控震源地震信號檢測關(guān)鍵技術(shù)研究中”運用多方向中值濾波技術(shù)壓制并行震源帶來的干擾噪聲，但該方法往往需要多次迭代才能達(dá)到較好信噪比，數(shù)據(jù)處理速度較慢，信號細(xì)節(jié)有一定損失，CN104765069A公開的“一種壓制同步激發(fā)采集領(lǐng)炮干擾方法”提出采用alpha-trimmed矢量中值閾值濾波去除隨機(jī)噪聲，該方法算法容易引入高頻干擾，CN103675903A公開的“一種隨機(jī)去噪多震源地震波場分離方法”沒有明確指出壓制噪聲方法，CN104536034A公開的“多震源并行激發(fā)采集與混合地震記錄分離方法”適于多震源采用互異控制信號時的噪聲壓制，而當(dāng)前大部分并行震源采用相同或相似的控制信號。技術(shù)實現(xiàn)要素：：本發(fā)明的目的就在于針對上述現(xiàn)有技術(shù)的不足，提供一種電磁式可控震源拋物線Radon變換壓制并行震源串?dāng)_噪聲方法。本發(fā)明的主要思想是：多震源并行激發(fā)技術(shù)提高了生產(chǎn)效率，但采集得到的地震數(shù)據(jù)同時混合了來自多炮點的多個并行震源地震信號，這樣獲得的數(shù)據(jù)并不能直接使用，本發(fā)明是通過將多源地震記錄即共炮點道集轉(zhuǎn)換到共接收點道集后，通過Radon正變換后在Radon域里分離有效信號與串?dāng)_噪聲，去噪后再進(jìn)行Radon反變換，實現(xiàn)了多源地震記錄中包含串?dāng)_噪聲的隨機(jī)噪聲壓制。本發(fā)明是通過以下技術(shù)方案實現(xiàn)的：電磁式可控震源拋物線Radon變換壓制并行震源串?dāng)_噪聲方法，包括以下步驟：a、在勘探區(qū)域內(nèi)，采用多個震源構(gòu)成一組并行震源，采用近同時激發(fā)，令各震源并行工作，獲得共炮點地震記錄集，各地震記錄均含有來自其它炮點的串?dāng)_地震信號；b、進(jìn)行共炮點道集到共檢波點道集的轉(zhuǎn)換，此時并行震源的串?dāng)_噪聲表現(xiàn)為隨機(jī)噪聲形式；c、在共接收點道集中運用拋物線Radon正變換。如公式：U＝LT(LLT+λE)-1d(1)這里d是t-x域數(shù)據(jù)，U是τ-q域數(shù)據(jù)，其中：L=ejωq1x12ejωq1x22Lejωq1xnx2ejωq2x12ejωq2x22Lejωq2xnx2LLLLejωqnqx12ejωqnqx22Lejωqnqxnx2---(2)]]>L稱為變換算子，包括nq行，nx列，L的第i行第k列對應(yīng)的元素i＝1,2Lnq，k＝1,2Lnx，LT為L共軛轉(zhuǎn)置算子，λ為阻尼因子，一般取0.1≤λ≤1，E為單位矩陣，nx為檢波點的采集道數(shù)；nq采用如下方式定義，任取滿足的q值,其中xmax為最大偏移距，Δx為道間距，f為地震波最小頻率，令Δq為曲率采樣間隔，任取滿足的Δq值，定義,d、根據(jù)U,繪制τ-q域圖像，其有效信號呈規(guī)律的點狀或拖長點狀，其余為噪聲信號，保留U中有效信號，其余數(shù)據(jù)置0，得到U'；e、對U'進(jìn)行拋物線反Radon變換，具體如公式d'＝LU'(3)其中d'為反變換后的t-x域共檢波點地震記錄；f、對所有共檢波點道集重復(fù)步驟c～e，直至完成全部共檢波點地震記錄；g、完成所有共檢波點道集到共炮點道集轉(zhuǎn)換。有益效果：經(jīng)試驗，本發(fā)明公開的一種電磁式可控震源拋物線Radon變換壓制并行震源串?dāng)_噪聲方法能夠?qū)崿F(xiàn)在多源地震勘探中壓制隨機(jī)延遲帶來的串?dāng)_噪聲,該算法處理數(shù)據(jù)快，信號細(xì)節(jié)保護(hù)好，在不增加地震勘探數(shù)據(jù)采集成本情況下，能夠有效改善并行震源地震勘探數(shù)據(jù)質(zhì)量。附圖說明：圖1第50道檢波器對應(yīng)的共接收點地震記錄,圖2經(jīng)過拋物線Radon變換處理后的地震記錄,圖3去噪后的Radon域地震記錄,圖4拋物線Radon反變換處理后的共接收點地震記錄,具體實施方式：下面結(jié)合附圖和實施例做進(jìn)一步的詳細(xì)說明：在本實施例中使用2個震源為一組的方法進(jìn)行激發(fā)，整個測區(qū)為5000m,兩源間距為1000m，炮間距和道間距都10m，記錄時間為3s,基于拋物線Radon變換去除并行震源串行噪聲的方法，包括以下步驟：a、用兩個震源為一組的方法激發(fā)100炮得到模擬共炮點混合地震記錄，其中設(shè)置隨機(jī)激發(fā)延時；b、道集轉(zhuǎn)換。在共炮點集中由于各震源并行激發(fā)時存在隨機(jī)激發(fā)延遲，因此將地震記錄從共炮點集轉(zhuǎn)換到其它集共檢波點集后，則由于激發(fā)延遲的隨機(jī)性導(dǎo)致干擾源記錄由連續(xù)形式轉(zhuǎn)換成隨機(jī)形式，本例中以第50道檢波器為例；c、在第50道檢波器形成的共接收點道集中根據(jù)公式(1):(2)作拋物線正Radon變換,本實施例中q≤5(s/km2)，取q＝2，在xmax＝1000，Δx＝10，f＝10，λ＝0.6，nq＝400，nx＝100，Δq≤0.1(s/km2)，取Δq＝0.01；d、本實施例中t-x域中的四層雙曲線經(jīng)變換后在τ-q域中成四個點狀或拖長點狀，提取這些信號；e、將提取的信號進(jìn)行拋物線反Radon變換，根據(jù)公式：d'＝LU'(3)；f、對所有共檢波點道集重復(fù)步驟c～e，直至完成全部共檢波點地震記錄；g、完成所有共檢波點道集到共炮點道集轉(zhuǎn)換。當(dāng)前第1頁1 2 3