本發(fā)明屬于地震數(shù)據(jù)處理,更具體地,涉及一種地表一致性反褶積方法、裝置、電子設(shè)備及存儲介質(zhì)。
背景技術(shù):
1、我國西部地區(qū)存在許多巨厚沙漠探區(qū)(如塔里木盆地、騰格里沙漠),由于這些工區(qū)地表條件屬典型的風(fēng)積沙漠地貌,沙丘厚度大。由于可控震源技術(shù)的高效安全環(huán)保等特性,可控震源技術(shù)被國內(nèi)各大油田公司應(yīng)用到這種巨厚沙漠區(qū)。但是沙丘厚度大并且地表松散,因此地表吸收衰減嚴重,能量下傳受限,在采集資料上根本看不到勘探目的層的信號,因此無法利用目的層完成地表一致性處理,常規(guī)地表一致性處理流程無法解決地表一致性的問題。
2、與目的層信號極弱相比,可控震源的回轉(zhuǎn)波能量強,信噪比高,沒有虛假信息干擾(如沒有相干噪音及折射多次波),并且對近地表變化敏感,因此回轉(zhuǎn)波是進行近地表效應(yīng)校正的最佳選擇。
3、公開于本發(fā)明背景技術(shù)部分的信息僅僅旨在加深對本發(fā)明的一般背景技術(shù)的理解,而不應(yīng)當(dāng)被視為承認或以任何形式暗示該信息構(gòu)成已為本領(lǐng)域技術(shù)人員所公知的現(xiàn)有技術(shù)。
技術(shù)實現(xiàn)思路
1、本發(fā)明的目的是提出一種地表一致性反褶積方法、裝置、電子設(shè)備及存儲介質(zhì),實現(xiàn)充分利用回轉(zhuǎn)波信噪比高、對近地表變化敏感的特點,有效解決巨厚沙漠地區(qū)近地表效應(yīng)引起的子波不一致性。
2、為實現(xiàn)上述目的,本發(fā)明提出了一種地表一致性反褶積方法、裝置、電子設(shè)備及存儲介質(zhì)。
3、根據(jù)本發(fā)明的第一方面,提出了一種地表一致性反褶積方法,包括:
4、獲取地震回轉(zhuǎn)波的初至數(shù)據(jù);
5、基于所述初至數(shù)據(jù),抽取xyo道集;
6、其中,x代表縱向,y代表橫向,o代表偏移距;
7、基于所述xyo道集建立回轉(zhuǎn)波地表一致性褶積模型;
8、基于所述回轉(zhuǎn)波地表一致性褶積模型計算所述xyo道集的炮點地表一致性分量和檢波點地表一致性分量;
9、基于所述炮點地表一致性分量和所述檢波點地表一致性分量計算所述xyo道集的地表一致性反褶積數(shù)據(jù)。
10、可選地,所述回轉(zhuǎn)波地表一致性褶積模型的表達式為:
11、時間域:pb(t)=w(t)*si(t)*rj(t)*q(t);
12、頻率域:
13、其中,si(t)為炮點的地表一致性效應(yīng),rj(t)為檢波點的地表一致性效應(yīng),和分別是si(t)和rj(t)的傅里葉變換,w(t)為xyo道集的平均響應(yīng),為w(t)的傅里葉變換,w(t)可以用xyo道集的疊加道w′(t)代替,q(t)為與地下構(gòu)造相關(guān)的殘余項,是q(t)的傅里葉變換,ω為角頻率,b為與炮點和檢波點索引對唯一對應(yīng)的索引號,(i,j)→b,i=1,…,ns,j=1,…,nr,ns為炮點唯一索引號的總數(shù),nr為檢波點唯一索引號的總數(shù)。
14、可選地,計算所述xyo道集的炮點地表一致性分量和檢波點地表一致性分量的表達式為:
15、
16、其中,ω=2πf,ω為角頻率,為頻率為f時索引b對應(yīng)地震道的值,為頻率為f時xyo道集的平均響應(yīng)分量,為頻率為f時炮點地表一致性分量,為頻率為f時檢波點地表一致性分量,q(f)為頻率為f時與地下構(gòu)造相關(guān)的地表一致性分量。
17、可選地,計算所述xyo道集的地表一致性反褶積數(shù)據(jù)的表達式為:
18、
19、其中,為所述xyo道集的地表一致性反褶積數(shù)據(jù),∈為防止被零除的微小正數(shù)。
20、可選地,,根據(jù)所述索引號是炮點和檢波點索引對唯一對應(yīng)的索引號,選擇對應(yīng)的所述炮點一致性分量和所述檢波點地表一致性分量進行計算,獲取所述xyo道集的地表一致性反褶積數(shù)據(jù)。
21、根據(jù)本發(fā)明的第二方面,提出了一種地表一致性反褶積裝置,包括:
22、獲取模塊,用于獲取獲取地震回轉(zhuǎn)波的初至數(shù)據(jù);
23、抽取模塊,用于基于所述初至數(shù)據(jù),抽取xyo道集;
24、其中,x代表縱向,y代表橫向,o代表偏移距;
25、建立模塊,用于基于所述xyo道集建立回轉(zhuǎn)波地表一致性褶積模型;
26、第一計算模塊,用于基于所述回轉(zhuǎn)波地表一致性褶積模型計算所述xyo道集的炮點地表一致性分量和檢波點地表一致性分量;
27、第二計算模塊,用于基于所述炮點地表一致性分量和所述檢波點地表一致性分量計算所述xyo道集的地表一致性反褶積數(shù)據(jù)。
28、可選地,所述回轉(zhuǎn)波地表一致性褶積模型的表達式為:
29、時間域:pb(t)=w(t)*si(t)*rj(t)*q(t);
30、頻率域:
31、其中,si(t)為炮點的地表一致性效應(yīng),rj(t)為檢波點的地表一致性效應(yīng),和分別是si(t)和rj(t)的傅里葉變換,w(t)為xyo道集的平均響應(yīng),為w(t)的傅里葉變換,w(t)可以用xyo道集的疊加道w′(t)代替,q(t)為與地下構(gòu)造相關(guān)的殘余項,是q(t)的傅里葉變換,ω為角頻率,b為與炮點和檢波點索引對唯一對應(yīng)的索引號,(i,j)→b,i=1,…,ns,j=1,…,nr,ns為炮點唯一索引號的總數(shù),nr為檢波點唯一索引號的總數(shù)。
32、可選地,計算所述xyo道集的炮點地表一致性分量和檢波點地表一致性分量的表達式為:
33、
34、其中,ω=2πf,ω為角頻率,為頻率為f時索引b對應(yīng)地震道的值,為頻率為f時xyo道集的平均響應(yīng)分量,為頻率為f時炮點地表一致性分量,為頻率為f時檢波點地表一致性分量,q(f)為頻率為f時與地下構(gòu)造相關(guān)的地表一致性分量。
35、根據(jù)本發(fā)明的第三方面,提出了一種電子設(shè)備,所述電子設(shè)備包括:
36、至少一個處理器;以及,
37、與所述至少一個處理器通信連接的存儲器;其中,
38、所述存儲器存儲有可被所述至少一個處理器執(zhí)行的指令,所述指令被所述至少一個處理器執(zhí)行,以使所述至少一個處理器能夠執(zhí)行權(quán)利要求1-5任一所述的地表一致性反褶積方法。
39、根據(jù)本發(fā)明的第四方面,提出了一種非暫態(tài)計算機可讀存儲介質(zhì),該非暫態(tài)計算機可讀存儲介質(zhì)存儲計算機指令,該計算機指令用于使計算機執(zhí)行第一方面任一所述的地表一致性反褶積方法。
40、本發(fā)明的有益效果在于:本發(fā)明利用回轉(zhuǎn)波建立地表一致性褶積模型,將炮點和檢波點處的近地表效應(yīng)分離出來,進而計算炮點地表一致性分量和檢波點地表一致性分量,根據(jù)炮點地表一致性分量和檢波點地表一致性分量實現(xiàn)地表一致性反褶積;本發(fā)明充分利用了回轉(zhuǎn)波信噪比高、對近地表變化敏感的特點,有效解決巨厚沙漠地區(qū)近地表效應(yīng)引起的子波不一致性。
41、本發(fā)明的系統(tǒng)具有其它的特性和優(yōu)點,這些特性和優(yōu)點從并入本文中的附圖和隨后的具體實施方式中將是顯而易見的,或者將在并入本文中的附圖和隨后的具體實施方式中進行詳細陳述,這些附圖和具體實施方式共同用于解釋本發(fā)明的特定原理。
1.一種地表一致性反褶積方法,其特征在于,包括:
2.根據(jù)權(quán)利要求1所述的地表一致性反褶積方法,其特征在于,所述回轉(zhuǎn)波地表一致性褶積模型的表達式為:
3.根據(jù)權(quán)利要求2所述的地表一致性反褶積方法,其特征在于,計算所述xyo道集的炮點地表一致性分量和檢波點地表一致性分量的表達式為:
4.根據(jù)權(quán)利要求3所述的地表一致性反褶積方法,其特征在于,計算所述xyo道集的地表一致性反褶積數(shù)據(jù)的表達式為:
5.根據(jù)權(quán)利要求4所述的地表一致性反褶積方法,其特征在于,根據(jù)所述索引號是炮點和檢波點索引對唯一對應(yīng)的索引號,選擇對應(yīng)的所述炮點一致性分量和所述檢波點地表一致性分量進行計算,獲取所述xyo道集的地表一致性反褶積數(shù)據(jù)。
6.一種電子設(shè)備,其特征在于,所述電子設(shè)備包括:
7.一種非暫態(tài)計算機可讀存儲介質(zhì),其特征在于,該非暫態(tài)計算機可讀存儲介質(zhì)存儲計算機指令,該計算機指令用于使計算機執(zhí)行權(quán)利要求1-5任一所述的地表一致性反褶積方法。
8.一種地表一致性反褶積裝置,其特征在于,包括:
9.根據(jù)權(quán)利要求8所述的地表一致性反褶積裝置,其特征在于,所述回轉(zhuǎn)波地表一致性褶積模型的表達式為:
10.根據(jù)權(quán)利要求8所述的地表一致性反褶積裝置,其特征在于,計算所述xyo道集的炮點地表一致性分量和檢波點地表一致性分量的表達式為: