一種多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地球物理領(lǐng)域，尤其涉及一種多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法。
【背景技術(shù)】
[0002]目前國(guó)內(nèi)大規(guī)模工業(yè)化應(yīng)用的采集多震源混合數(shù)據(jù)分離方法主要是在偽分離道集的共偏移距道集進(jìn)行濾波去噪，得到分離數(shù)據(jù)。
[0003]然而，由于上述過(guò)程是利用有效信號(hào)的統(tǒng)計(jì)特性來(lái)壓制混疊噪音，因此，對(duì)于混疊信噪比較低的多源混疊數(shù)據(jù)，數(shù)據(jù)分離的質(zhì)量較差進(jìn)而導(dǎo)致利用常規(guī)方法進(jìn)行偏移成像的結(jié)果較差。

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

[0004]本發(fā)明提供一種多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法，去除混疊噪音的影響，能夠提高數(shù)據(jù)分離的質(zhì)量。
[0005]本發(fā)明解決技術(shù)問(wèn)題采用如下技術(shù)方案:一種多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法，包括:S10、獲取時(shí)域混疊的二維炮記錄m;S20、根據(jù)所述二維炮記錄m和預(yù)設(shè)的從共炮點(diǎn)道集到共中心點(diǎn)道集的轉(zhuǎn)換函數(shù)Ψ，獲取混疊數(shù)據(jù)在共中心點(diǎn)道集b;S30、根據(jù)預(yù)設(shè)一致性濾波器Π和所述混疊數(shù)據(jù)在共中心點(diǎn)道集b獲取分離的有效信號(hào)S和分離的鄰炮干擾N;S40、根據(jù)所述共中心點(diǎn)道集的轉(zhuǎn)換函數(shù)Ψ將所述分離的有效信號(hào)S和分離的鄰炮干擾N分別轉(zhuǎn)換到共炮點(diǎn)道集中，得到道集轉(zhuǎn)換后分離的有效信號(hào)S'和道集轉(zhuǎn)換后分離的鄰炮干擾N';S50、根據(jù)所述道集轉(zhuǎn)換后分離的有效信號(hào)S'和道集轉(zhuǎn)換后分離的鄰炮干擾N’，分別獲取分離的鄰炮干擾中泄露的有效信號(hào)NS、分離的鄰炮干擾中提純的鄰炮干擾NN和分離的有效信號(hào)中提純的有效信號(hào)SS以及分離的有效信號(hào)中殘留的鄰炮干擾SN;S60、根據(jù)所述分離的鄰炮干擾中泄露的有效信號(hào)NS、分離的鄰炮干擾中提純的鄰炮干擾NN、分離的有效信號(hào)中提純的有效信號(hào)SS以及分離的有效信號(hào)中殘留的鄰炮干擾SN，對(duì)應(yīng)重新組合，構(gòu)建最終分離的有效信號(hào)SI和最終分離的鄰炮干擾NI。
[0006]可選的，本發(fā)明實(shí)施例提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中，所述S20，包括:將所述二維炮記錄m乘以所述轉(zhuǎn)換函數(shù)Ψ，得到所述混疊數(shù)據(jù)在共中心點(diǎn)道集b。
[0007]可選的，本發(fā)明實(shí)施例提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中，所述S30，包括S301、將所述預(yù)設(shè)一致性濾波器Π乘以所述混疊數(shù)據(jù)在共中心點(diǎn)道集b，得到分離的有效信號(hào)S;S302、將所述混疊數(shù)據(jù)在共中心點(diǎn)道集b減去所述分離的有效信號(hào)S，得到所述分離的鄰炮干擾N。
[0008]可選的，本發(fā)明實(shí)施例提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中，所述S40，包括:S401、根據(jù)所述轉(zhuǎn)換函數(shù)Ψ獲取從共中心點(diǎn)道集到共炮點(diǎn)道集的轉(zhuǎn)換函數(shù)WT;S402、將所述轉(zhuǎn)換函數(shù)屯了乘以所述分離的有效信號(hào)S，得到道集轉(zhuǎn)換后分離的有效信號(hào)V ;S403、將所述轉(zhuǎn)換函數(shù)$7乘以所述分離的鄰炮干擾N，得到道集轉(zhuǎn)換后分離的鄰炮干擾N'。
[0009]可選的，本發(fā)明實(shí)施例提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中，所述S50，包括:S501、將所述預(yù)設(shè)一致性濾波器Π乘以所述道集轉(zhuǎn)換后分離的鄰炮干擾N ’，得到所述分離的鄰炮干擾中提純的鄰炮干擾NN; S502、將所述道集轉(zhuǎn)換后分離的鄰炮干擾N'減去所述分離的鄰炮干擾中提純的鄰炮干擾NN，得到所述分離的鄰炮干擾中泄露的有效信號(hào)NS;S503、將所述預(yù)設(shè)一致性濾波器Π乘以所述道集轉(zhuǎn)換后分離的有效信號(hào)SS得到所述分離的有效信號(hào)中提純的有效信號(hào)SS;S504、將所述道集轉(zhuǎn)換后分離的有效信號(hào)S'減去所述分離的有效信號(hào)中提純的有效信號(hào)SS，得到所述分離的有效信號(hào)中殘留的鄰炮干擾SN0
[0010]可選的，本發(fā)明實(shí)施例提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中，所述S60，包括:S601、將所述分離的有效信號(hào)中提純的有效信號(hào)SS加上所述分離的鄰炮干擾中泄露的有效信號(hào)NS，得到最終分離的有效信號(hào)SI ;S602、將所述分離的有效信號(hào)中殘留的鄰炮干擾SN加上所述分離的鄰炮干擾中提純的鄰炮干擾NN，得到最終分離的鄰炮干擾
Nl0
[0011 ]可選的，本發(fā)明實(shí)施例提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中，所述一致性濾波器為多維矢量中值一致性濾波器。
[0012]本發(fā)明具有如下有益效果:采用多域多次分離的辦法去除混疊噪音的影響，能夠提高數(shù)據(jù)分離的質(zhì)量;上述過(guò)程能夠?yàn)槭褂贸Ｒ?guī)的地震方法進(jìn)行偏移成像提供數(shù)據(jù)，提高偏移的質(zhì)量。本發(fā)明實(shí)施例提供的技術(shù)方案解決了現(xiàn)有技術(shù)由于利用有效信號(hào)的統(tǒng)計(jì)特性來(lái)壓制混疊噪音，因此，解決對(duì)于混疊信噪比比較低的多源混疊炮數(shù)據(jù)，數(shù)據(jù)分離的質(zhì)量較差，進(jìn)而導(dǎo)致利用常規(guī)方法進(jìn)行偏移成像的結(jié)果較差的問(wèn)題。
【附圖說(shuō)明】
[0013]圖1為本發(fā)明實(shí)施例1提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法的流程圖；
[0014]圖2為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中二維炮記錄的不意圖；
[0015]圖3為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中混疊數(shù)據(jù)在共中心點(diǎn)道集b的示意圖；
[0016]圖4為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中分離的有效信號(hào)S的示意圖；
[0017]圖5為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中分離的鄰炮干擾N的不意圖；
[0018]圖6為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中道集轉(zhuǎn)換后在共炮點(diǎn)道集中的有效信號(hào)S'的示意圖；
[0019]圖7為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中道集轉(zhuǎn)換后在共炮點(diǎn)道集中的鄰炮干擾N'的示意圖；
[0020]圖8為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中最終分離的有效信號(hào)Si的示意圖；
[0021]圖9為圖1所示的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法中重新構(gòu)建的鄰炮干擾NI的示意圖。
【具體實(shí)施方式】
[0022]下面結(jié)合實(shí)施例及附圖對(duì)本發(fā)明的技術(shù)方案作進(jìn)一步闡述。
[0023]實(shí)施例1
[0024]如圖1所示，本發(fā)明實(shí)施例提供的多域多次分離的高效采集多震源混合數(shù)據(jù)分離方法，包括:
[0025]步驟101，獲取時(shí)域混疊的二維炮記錄m。
[0026]步驟102，根據(jù)所述二維炮記錄m和預(yù)設(shè)的從共炮點(diǎn)道集到共中心點(diǎn)道集的轉(zhuǎn)換函數(shù)Ψ，獲取混疊數(shù)據(jù)在共中心點(diǎn)道集b。
[0027]其中，通過(guò)步驟102獲取共中心點(diǎn)道集b的過(guò)程包括:將所述二維炮記錄m乘以所述轉(zhuǎn)換函數(shù)Ψ，得到所述混疊數(shù)據(jù)在共中心點(diǎn)道集b。
[0028]即混疊數(shù)據(jù)在共中心點(diǎn)的道集b= Wm 公式一
[0029]步驟103，根據(jù)預(yù)設(shè)一致性濾波器Π和所述混疊數(shù)據(jù)在共中心點(diǎn)道集b獲取分離的有效信號(hào)S和分離的鄰炮干擾N。
[0030]其中，通過(guò)步驟103獲取分離的有效信號(hào)S和分離的鄰炮干擾N的過(guò)程包括:將所述預(yù)設(shè)一致性濾波器Π乘以所述混疊數(shù)據(jù)在共中心點(diǎn)道集b，得到分離的有效信號(hào)S;將所述混疊數(shù)據(jù)在共中心點(diǎn)道集b減去所述分離的有效信號(hào)S，得到所述分離的鄰炮干擾N。該一致性濾波器可以為一致性濾波器，具體可以為多維矢量中值一致性濾波器。
[0031 ]即分離的有效信號(hào)S= IIb 公式二
[0032]分離的鄰炮干擾N=b_S
[0033]步驟104，根據(jù)所述共中心點(diǎn)道集的轉(zhuǎn)換函數(shù)Ψ將所述分離的有效信號(hào)S和分離的鄰炮干擾N分別轉(zhuǎn)換到共炮點(diǎn)道集中，得到道集轉(zhuǎn)換后分離的有效信號(hào)S'和道集轉(zhuǎn)換后分離的鄰炮干擾N'。
[0034]在本實(shí)施例中，通過(guò)步驟104獲取道集轉(zhuǎn)換后分離的有效信號(hào)S'和道集轉(zhuǎn)換后分離的鄰炮干擾K的過(guò)程包括:根據(jù)所述轉(zhuǎn)換函數(shù)Ψ獲取從共中心點(diǎn)道集到共炮點(diǎn)道集的轉(zhuǎn)換函數(shù)Ψτ;將所述轉(zhuǎn)換函數(shù)Ψτ乘以所述分離的有效信號(hào)S，得到道集轉(zhuǎn)換后分離的有效信號(hào)S';將所述轉(zhuǎn)換函數(shù)$7乘以所述分離的鄰炮干擾N，得到道集轉(zhuǎn)換后分離的鄰炮干擾K。
[0035]即道集轉(zhuǎn)換后分離的有效信號(hào)S7=WtS公式三
[0036]道集轉(zhuǎn)換后分離的鄰炮干擾V =ΨΤΝ
[0037]步驟105，根據(jù)所述道集轉(zhuǎn)換后分離的有效信號(hào)S'和道集轉(zhuǎn)換后分離的鄰炮干擾N丨，分別獲取分離的鄰炮干擾中泄露的有效信號(hào)NS、分離的鄰炮干擾中提純的鄰炮干擾NN和分離的有效信號(hào)中提純的有效信號(hào)SS以及分離的有效信號(hào)中