本發明專利技術涉及一種地震記錄面波干擾壓制方法及裝置，獲取地震記錄xi(t)；根據所述地震記錄xi(t)獲取面波干擾高截頻fg；根據所述面波干擾高截頻fg確定變化曲線q(t)；對地震記錄xi(t)做高截頻為fg的低通濾波，得到面波干擾頻帶的地震記錄yi(t)；并對地震記錄xi(t)做低截頻為fg高截頻為2fg的帶通濾波，得到非面波干擾頻段的地震記錄zi(t)；根據面波干擾頻帶的地震記錄yi(t)獲取振幅變化曲線py(t)；并根據非面波干擾頻段的地震記錄zi(t)獲取振幅變化曲線pz(t)；根據變化曲線q(t)、振幅變化曲線py(t)和振幅變化曲線pz(t)對面波干擾頻段yi(t)進行調制獲取干擾頻帶地震記錄根據干擾頻帶地震記錄面波干擾頻段yi(t)以及地震記錄xi(t)獲取面波干擾壓制之后的地震記錄
【技術實現步驟摘要】
本專利技術涉及地震記錄處理領域，特別涉及一種地震記錄面波干擾壓制方法及裝 置。
技術介紹
地震勘探是一種利用人工地震技術探測地下結構的勘探方法。它按照一定的方式人工激發地震波，利用稱之為檢波器的裝置接收來自地下的反射信號，通過對反射信號的處理和分析探測地下結構。檢波器裝置在接收地震信號的同時，也接收了多種類型的噪聲干擾，如何有效地壓制噪聲，恢復被噪聲污染的地震信號，是勘探地震資料處理工作的重要研究內容。在各類噪聲干擾中，面波是最為普遍、同時也是最難壓制的噪聲干擾。與有效信號相比，面波干擾具有低頻、低速、頻散、強能量等特點，在地震記錄上呈掃帚狀分布，其能量往往是有效信號的上百倍，在面波污染的地震記錄上，幾乎見不到有效信號的影子。低通濾波是最為經典的面波壓制方法，它利用面波干擾與有效信號在頻率上的差異利用低通濾波消除面波干擾。由于該方法在消除面波干擾的同時，也消除了有效信號的低頻分量，不利于地震資料高精度保幅處理，所以，目前該方法在工業界很少使用。f-k濾波、Radon變換濾波、頻率空間域預測濾波等是另外一類面波壓制的方法。盡管采用的數學方法不同，但這類方法本質上都是利用了面波干擾的空間相干性，其壓制效果在很大程度上依賴于面波干擾在線性方向的相關程度。如果面波干擾具有數學意義上完整的相關性，即使利用簡單的f-k濾波，其壓制噪聲、恢復弱信號的能力也是非常滿意的。然而，由于接收條件和傳播路徑的差異，以及面波固有的頻散特性，面波干擾更多地表現為視覺上的相干性，并不具備完全的數學意義上的相干性。因此，直接利用空間相干類濾波方法往往達不到理想的面波壓制效果。針對傅立葉變換存在的問題，近十年來，先后發展了基于小波變換、廣義S變換、curvelet變換的面波壓制方法，這類方法的共同特點是，采用了最新的數學變換方法，但實現和操作比較復雜，不利于工業化推廣應用，目前尚未見到大規模工業應用的報到。錐形濾波(也稱FKK濾波，又稱十字交叉濾波)是工業界地震資料處理系統一種新的三維地震資料面波壓制方法，該方法從本質上將仍然屬于f-k濾波，只是由二維擴展到三維，且對地震記錄進行了分選和重排，使之更加適應三維f-k濾波的要求。該方法較二維相干噪聲壓制方法具有更強的面波壓制能力，但依然無法擺脫面波頻散造成的視線性問題，處理之后的地震記錄依然殘存有較強的面波能量。綜上所述，在保持地震信號低頻分量的同時，有效地壓制面波干擾對地震信號的影響，增強地震記錄的信噪比，是目前高精度地震勘探迫切需要解決的技術問題。
技術實現思路
本專利技術的目的是針對上述問題，提出一種地震記錄面波干擾壓制方法及裝置，可以壓制面波干擾對地震記錄有效信號的影響，有效地恢復被面波污染的地震信號，提高地震記錄信噪比，增強地震信號反映地下結構的能力，提高地震勘探精度。為實現上述目的，本專利技術提供了一種地震記錄面波干擾壓制方法，包括獲取地震記錄XiU),其中，i = 1,2,…η, η為地震道數，t是時間，單位是毫秒；根據所述地震記錄Xi (t)獲取面波干擾高截頻fg，單位為赫茲；根據所述面波干擾高截頻fg確定變化曲線q(t)；對地震記錄Xi (t)做高截頻為fg的低通濾波，得到面波干擾頻帶的地震記錄Ji (t)i = 1,2, *··η ;并對地震記錄Xi (t)做低截頻為fg高截頻為2fg的帶通濾波,得到非面波干擾頻段的地震記錄Zi (t), i = 1,2,…η;根據面波干擾頻帶的地震記錄yi(t)獲取振幅變化曲線py(t);并根據非面波干擾頻段的地震記錄Zi (t)獲取振幅變化曲線Pz (t);根據變化曲線q(t)、振幅變化曲線Py(t)和振幅變化曲線Pz(t)對面波干擾頻段Yi (t)進行調制獲取干擾頻帶地震記錄3^0;根據干擾頻帶地震記錄又.(O、面波干擾頻段Ii (t)以及地震記錄Xi (t)獲取面波干擾壓制之后的地震記錄SG)O可選的，在本專利技術一實施例中，所述根據所述地震記錄Xi (t)獲取面波干擾高截頻fg包括對地震記錄Xi⑴作傅里葉變換獲取Xi (f);根據傅里葉變換后的地震記錄Xi (f)計算振幅譜Ai (f)；計算η道地震記錄的平均振幅譜A (f)；對平均振幅譜A (f)進行低通濾波，得到平均振幅譜A (f)的變化趨勢B (f)；獲取平均振幅譜的變化趨勢B (f)的第一極值點所對應的頻率fd ；根據頻率fd確定面波干擾高截頻fg。可選的，在本專利技術一實施例中，所述平均振幅譜A(f)的表達式為權利要求1.一種地震記錄面波干擾壓制方法，其特征在于，包括獲取地震記錄Xi (t),其中，i = 1,2,…η, η為地震道數；根據所述地震記錄Xi (t)獲取面波干擾高截頻fg ；根據所述面波干擾高截頻fg確定變化曲線q(t)；對地震記錄Xi (t)做高截頻為fg的低通濾波，得到面波干擾頻帶的地震記錄Ji (t)，i=1,2, *··η ;并對地震記錄Xi (t)做低截頻為fg高截頻為2fg的帶通濾波,得到非面波干擾頻段的地震記錄Zi (t) i = 1,2,…η;根據面波干擾頻帶的地震記錄yi(t)獲取振幅變化曲線py(t);并根據非面波干擾頻段的地震記錄zdt)獲取振幅變化曲線pz(t)；根據變化曲線q(t)、振幅變化曲線py(t)和振幅變化曲線？2(0對面波干擾頻段yi(t)進行調制，獲取干擾頻帶地震記錄滅(O .根據干擾頻帶地震記錄又O)、面波干擾頻段yjt)以及地震記錄Xi (t)獲取面波干擾壓制之后的地震記錄S (O。2.根據權利要求I所述的方法，其特征在于，所述根據地震記錄Xi(t)獲取面波干擾高截頻fg包括對地震記錄Xi (t)作傅里葉變換獲取Xi (f)；根據傅里葉變換后的地震記錄Xi (f)計算振幅譜Ai (f)；計算η道地震記錄的平均振幅譜A (f)；對平均振幅譜A (f)進行低通濾波，得到平均振幅譜A (f)的變化趨勢B (f)；獲取平均振幅譜的變化趨勢B (f)的第一極值點所對應的頻率fd ；根據頻率fd確定面波干擾高截頻fg。3.根據權利要求2所述的方法，其特征在于，所述平均振幅譜A(f)的表達式為4.根據權利要求2所述的方法，其特征在于，所述平均振幅譜A(f)的變化趨勢B(f)的表達式為5.根據權利要求2所述的方法，其特征在于，所述面波干擾高截頻fg的表達式為fg = I· 5fdo6.根據權利要求I所述的方法，其特征在于，所述根據面波干擾高截頻fg確定不同頻段能量之比變化曲線q(t)包括對η個地震道中炮檢距最大的m個地震道所對應的地震記錄Xi (t)進行高截頻為fg的低通濾波，得到O. O到fg頻段的地震記錄(t), i = n-m+1, n-m+2,----, η ;并對η個地震道中炮檢距最大的m個地震道所對應的地震記錄Xi⑴進行低截頻為fg、高截頻為2fg的帶通濾波，得到fg到2fg頻段的地震記錄7.根據權利要求6所述的方法，其特征在于，所述變化曲線q(t)的表達式為8.根據權利要求I所述的方法，其特征在于，所述振幅變化曲線py(t)的表達式為9.根據權利要求I或8所述的方法，其特征在于，所述振幅變化曲線pz(t)的表達式為10.根據權利要求I所述的方法，其特征在于，所述干擾頻帶地震記錄又(O的表達式為11.根據權利要求I所述的方本文檔來自技高網...
【技術保護點】
一種地震記錄面波干擾壓制方法，其特征在于，包括：獲取地震記錄xi(t)，其中，i＝1，2，…n，n為地震道數；根據所述地震記錄xi(t)獲取面波干擾高截頻fg；根據所述面波干擾高截頻fg確定變化曲線q(t)；對地震記錄xi(t)做高截頻為fg的低通濾波，得到面波干擾頻帶的地震記錄yi(t)，i＝1，2，…n；并對地震記錄xi(t)做低截頻為fg高截頻為2fg的帶通濾波，得到非面波干擾頻段的地震記錄zi(t)i＝1，2，…n；根據面波干擾頻帶的地震記錄yi(t)獲取振幅變化曲線py(t)；并根據非面波干擾頻段的地震記錄zi(t)獲取振幅變化曲線pz(t)；根據變化曲線q(t)、振幅變化曲線py(t)和振幅變化曲線pz(t)對面波干擾頻段yi(t)進行調制，獲取干擾頻帶地震記錄根據干擾頻帶地震記錄面波干擾頻段yi(t)以及地震記錄xi(t)獲取面波干擾壓制之后的地震記錄FDA00002315757000011.tif,FDA00002315757000012.tif,FDA00002315757000013.tif
【技術特征摘要】
【專利技術屬性】
技術研發人員：李國發，秦德海，付立新，祝文亮，岳英，翟桐立，
申請(專利權)人：中國石油大學北京，
類型：發明
國別省市：
還沒有人留言評論。發表了對其他瀏覽者有用的留言會獲得科技券。