1.一種不規(guī)則地震數(shù)據(jù)的重建方法,其特征在于,包括以下步驟:
獲取一個炮檢距矢量面元OVT道集,所述OVT道集是預(yù)先將原始地震數(shù)據(jù)進(jìn)行重排后得到的多個OVT道集之一;
將當(dāng)前OVT道集進(jìn)行規(guī)則化處理,并將規(guī)則化處理后的OVT道集內(nèi)的地震道沿時間方向做一維傅里葉變換,獲得多個頻率分量的地震數(shù)據(jù);
對每個頻率分量的地震數(shù)據(jù)按照預(yù)設(shè)的重建規(guī)則進(jìn)行重建,獲得每個頻率分量的新地震數(shù)據(jù);
將所述每個頻率分量的新地震數(shù)據(jù)沿時間方向做一維傅里葉反變換,獲得當(dāng)前OVT道集的時域新地震數(shù)據(jù);
重復(fù)以上步驟,以獲得所有OVT道集的時域新地震數(shù)據(jù);
將所有OVT道集的時域新地震數(shù)據(jù)與所述原始地震數(shù)據(jù)進(jìn)行合并處理,獲得重建結(jié)果。
2.根據(jù)權(quán)利要求1所述的不規(guī)則地震數(shù)據(jù)的重建方法,其特征在于,所述對每個頻率分量按照預(yù)設(shè)的重建規(guī)則進(jìn)行重建,具體包括:
取出一個頻率分量的地震數(shù)據(jù)作為當(dāng)前待重建對象;
將當(dāng)前待重建對象沿空間方向做二維傅里葉變換,獲取空間域地震數(shù)據(jù);
根據(jù)預(yù)設(shè)的自適應(yīng)門檻函數(shù)確定當(dāng)前OVT道集的門檻值;
將所述空間域地震數(shù)據(jù)中,其樣點(diǎn)值小于所述門檻值的樣點(diǎn)值置零,獲得新空間域地震數(shù)據(jù);
將所述新空間域地震數(shù)據(jù)沿空間方向做二維傅里葉反變換,并用二維傅里葉反變換后獲得的變換結(jié)果替換當(dāng)前頻率分量的地震數(shù)據(jù)中相同位置處的地震數(shù)據(jù),從而獲取當(dāng)前頻率分量的新地震數(shù)據(jù);
重復(fù)上述步驟,以獲得當(dāng)前OVT道集的每個頻率分量的新地震數(shù)據(jù)。
3.根據(jù)權(quán)利要求2所述的不規(guī)則地震數(shù)據(jù)的重建方法,其特征在于,在所述獲取當(dāng)前頻率分量的新地震數(shù)據(jù)之后,還包括以下步驟:
判斷當(dāng)前獲得的當(dāng)前頻率分量的新地震數(shù)據(jù)是否達(dá)到預(yù)設(shè)要求;
如果未達(dá)到,則將當(dāng)前頻率分量的新地震數(shù)據(jù)作為當(dāng)前待重建對象,進(jìn)行迭代重建,直至迭代次數(shù)達(dá)到預(yù)設(shè)的迭代次數(shù)上限為止。
4.根據(jù)權(quán)利要求1所述的不規(guī)則地震數(shù)據(jù)的重建方法,其特征在于,所述將當(dāng)前OVT道集進(jìn)行規(guī)則化處理,包括:
對當(dāng)前OVT道集內(nèi)的地震道,根據(jù)其所屬的空間網(wǎng)格分別沿X方向順序和Y方向順序進(jìn)行重排處理;
在所述重排處理后,對于當(dāng)前OVT道集中的每個空間網(wǎng)格,判斷其內(nèi)是否有且僅有一個地震道;
如果有空間網(wǎng)格內(nèi)沒有地震道,則在該空間網(wǎng)格中生成一個其內(nèi)所有樣點(diǎn)值為零的地震道;
如果有空間網(wǎng)格內(nèi)有至少兩個地震道,則根據(jù)預(yù)設(shè)的目標(biāo)函數(shù)從中選擇一個作為該空間網(wǎng)格的地震道。
5.根據(jù)權(quán)利要求4所述的不規(guī)則地震數(shù)據(jù)的重建方法,其特征在于,所述的目標(biāo)函數(shù)包括:
其中,obj為目標(biāo)函數(shù)最小值,Off為當(dāng)前地震道的炮檢距,OffOVT為當(dāng)前OVT道集的中心炮檢距,OffInt為當(dāng)前OVT道集的炮檢距間隔,Azm為當(dāng)前地震道的方位角,AzmOVT為當(dāng)前OVT道集的中心方位角,AzmInt為當(dāng)前OVT道集的方位角間隔;X和Y分別為當(dāng)前地震道的東坐標(biāo)和北坐標(biāo),XBin和YBin分別為當(dāng)前空間網(wǎng)格的東坐標(biāo)和北坐標(biāo),DX和DY分別為當(dāng)前空間網(wǎng)格的東坐標(biāo)方向邊長和北坐標(biāo)方向邊長。
6.根據(jù)權(quán)利要求2所述的不規(guī)則地震數(shù)據(jù)的重建方法,其特征在于,所述根據(jù)預(yù)設(shè)的自適應(yīng)門檻函數(shù)確定當(dāng)前OVT道集的門檻值,包括:
確定當(dāng)前OVT道集中每個地震道的振幅譜,并將所有地震道的振幅譜逆序排列,表示為g(j),j=1,2,...,N;
根據(jù)公式獲取當(dāng)前OVT道集中所有地震道的振幅譜累積量s;
根據(jù)公式獲取當(dāng)前OVT道集的振幅譜累積量間隔sinc,其中niter是總迭代次數(shù);
根據(jù)公式thd(l)=g(k)獲取當(dāng)前OVT道集每個迭代所用的門檻值,其中k是滿足的最小值,l為迭代次數(shù)索引,l=1,2,...,niter。
7.一種不規(guī)則地震數(shù)據(jù)的重建裝置,其特征在于,包括:
OVT道集獲取模塊,用于獲取一個OVT道集,所述OVT道集是預(yù)先將原始地震數(shù)據(jù)進(jìn)行重排后得到的多個OVT道集之一;
OVT道集規(guī)則化模塊,用于將當(dāng)前OVT道集進(jìn)行規(guī)則化處理;
一維傅里葉變換模塊,用于將規(guī)則化處理后的OVT道集內(nèi)的地震道沿時間方向做一維傅里葉變換,獲得多個頻率分量的地震數(shù)據(jù);
數(shù)據(jù)重建模塊,用于對每個頻率分量的地震數(shù)據(jù)按照預(yù)設(shè)的重建規(guī)則進(jìn)行重建,獲得每個頻率分量的新地震數(shù)據(jù);
一維傅里葉反變換模塊,用于將所述每個頻率分量的新地震數(shù)據(jù)沿時間方向做一維傅里葉反變換,獲得當(dāng)前OVT道集的時域新地震數(shù)據(jù);
循環(huán)控制模塊,用于控制上述各個模塊的循環(huán),以獲得所有OVT道集的時域新地震數(shù)據(jù);
重建結(jié)果獲取模塊,用于將所有OVT道集的時域新地震數(shù)據(jù)與所述原始地震數(shù)據(jù)進(jìn)行合并處理,獲得重建結(jié)果。
8.根據(jù)權(quán)利要求7所述的不規(guī)則地震數(shù)據(jù)的重建裝置,其特征在于,所述數(shù)據(jù)重建模塊,具體包括:
重建對象確定子模塊,用于取出一個頻率分量的地震數(shù)據(jù)作為當(dāng)前待重建對象;
二維傅里葉變換子模塊,用于將當(dāng)前待重建對象沿空間方向做二維傅里葉變換,獲取空間域地震數(shù)據(jù);
門檻值獲取子模塊,用于根據(jù)預(yù)設(shè)的自適應(yīng)門檻函數(shù)確定當(dāng)前OVT道集的門檻值;
樣點(diǎn)值取舍子模塊,用于將所述空間域地震數(shù)據(jù)中,其樣點(diǎn)值小于所述門檻值的樣點(diǎn)值置零,獲得新空間域地震數(shù)據(jù);
二維傅里葉反變換子模塊,用于將所述新空間域地震數(shù)據(jù)沿空間方向做二維傅里葉反變換;
數(shù)據(jù)替換子模塊,用于將二維傅里葉反變換后獲得的變換結(jié)果替換當(dāng)前頻率分量的地震數(shù)據(jù)中相同位置處的地震數(shù)據(jù),從而獲取當(dāng)前頻率分量的新地震數(shù)據(jù);
循環(huán)控制子模塊,用于控制上述各個子模塊的循環(huán),以獲得當(dāng)前OVT道集的每個頻率分量的新地震數(shù)據(jù)。
9.根據(jù)權(quán)利要求8所述的不規(guī)則地震數(shù)據(jù)的重建裝置,其特征在于,還包括:
待重建對象迭代模塊,用于在所述獲取當(dāng)前頻率分量的新地震數(shù)據(jù)之后,判斷當(dāng)前獲得的當(dāng)前頻率分量的新地震數(shù)據(jù)是否達(dá)到預(yù)設(shè)要求;如果未達(dá)到,則將當(dāng)前頻率分量的新地震數(shù)據(jù)作為當(dāng)前待重建對象,進(jìn)行迭代重建,直至迭代次數(shù)達(dá)到預(yù)設(shè)的迭代次數(shù)上限為止。
10.根據(jù)權(quán)利要求7所述的不規(guī)則地震數(shù)據(jù)的重建裝置,其特征在于,所述將OVT道集規(guī)則化模塊,包括:
重排子模塊,用于對當(dāng)前OVT道集內(nèi)的地震道,根據(jù)其所屬的空間網(wǎng)格分別沿X方向順序和Y方向順序進(jìn)行重排處理;
判斷子模塊,用于在所述重排處理后,對于當(dāng)前OVT道集中的每個空間網(wǎng)格,判斷其內(nèi)是否有且僅有一個地震道;
生成子模塊,用于在當(dāng)前OVT道集中有空間網(wǎng)格內(nèi)沒有地震道時,在該空間網(wǎng)格中生一個其內(nèi)所有樣點(diǎn)值為零的地震道;
選擇子模塊,用于在當(dāng)前OVT道集中有空間網(wǎng)格內(nèi)有至少兩個地震道時,根據(jù)預(yù)設(shè)的目標(biāo)函數(shù)從中選擇一個作為該空間網(wǎng)格的地震道。
11.根據(jù)權(quán)利要求10所述的不規(guī)則地震數(shù)據(jù)的重建裝置,其特征在于,所述的目標(biāo)函數(shù)包括:
其中,obj為目標(biāo)函數(shù)最小值,Off為當(dāng)前地震道的炮檢距,OffOVT為當(dāng)前OVT道集的中心炮檢距,OffInt為當(dāng)前OVT道集的炮檢距間隔,Azm為當(dāng)前地震道的方位角,AzmOVT為當(dāng)前OVT道集的中心方位角,AzmInt為當(dāng)前OVT道集的方位角間隔;X和Y分別為當(dāng)前地震道的東坐標(biāo)和北坐標(biāo),XBin和YBin分別為當(dāng)前空間網(wǎng)格的東坐標(biāo)和北坐標(biāo),DX和DY分別為當(dāng)前空間網(wǎng)格的東坐標(biāo)方向邊長和北坐標(biāo)方向邊長。
12.根據(jù)權(quán)利要求8所述的不規(guī)則地震數(shù)據(jù)的重建裝置,其特征在于,所述門檻值獲取子模塊根據(jù)預(yù)設(shè)的自適應(yīng)門檻函數(shù)確定當(dāng)前OVT道集的門檻值,包括:
確定當(dāng)前OVT道集中每個地震道的振幅譜,并將所有地震道的振幅譜逆序排列,表示為g(j),j=1,2,...,N;
根據(jù)公式獲取當(dāng)前OVT道集中所有地震道的振幅譜累積量s;
根據(jù)公式獲取當(dāng)前OVT道集的振幅譜累積量間隔sinc,其中niter是總迭代次數(shù);
根據(jù)公式thd(l)=g(k)獲取當(dāng)前OVT道集每個迭代所用的門檻值,其中k是滿足的最小值,l為迭代次數(shù)索引,l=1,2,...,niter。