17] 對單炮地震記錄G' = (g，Jmxn進行K-L變換Y = UTG，，其中，Y = (Yij)ixnS 的K-L變換，yij為Y的第i個地震道、第j個采樣點，i = 1，j = 1，2,. . .，n，U % U的 轉置矩陣；
[0018] 步驟6)由步驟5)的Y = UV，利用K-L反變換G" = UY來重構地滾波模型 G" = (g" ij)mXn，G' = (Sij)mxn與重建地滾波G" = (g" UOmxn的最小平方誤差為：ε = |G' -G" |2-0,其中，g" u為重建的地滾波模型G"的第i個地震道、第j個采樣點，i =I, 2,. . . , m ； j = 1, 2,. . . , η ；
[0019] 步驟7)將重建的地滾波模型G" = (g" υ_)πΧη對應步驟4)進行反線性動校及靜 校正量反時移-Δ ^(1 = 1，2，...，m)，得到地滾波最終模型G" ' =(g" ' υ)πΧη4" ' υ 為地滾波最終模型G"'的第i個地震道、第j個采樣點，i = 1，2,. . .，m ;j = 1，2,. . .，η ;
[0020] 步驟8)得到壓制地滾波后的單炮地震記錄X" = (χ〃 υ)πΧη= (Xij)mxn-(g〃 ' ij)mXn，x" ijS單炮地震記錄X"的第i個地震道、第j個采樣點，i = l，2，...，m;j = 1，2，· · ·，η ;
[0021] 步驟9)重復步驟3)至步驟8)，去除單炮地震記錄X' = (X' Jmxn上不同頻率 范圍、不同視速度的地滾波，得到最終壓制地滾波后的單炮地震記錄X" '= (X" ' υ)πΧη， 其中X" ^ ^為單炮地震記錄X" ^的第i個地震道、第j個采樣點，i = 1，2,...，m ;j = I，2，· · ·，η ;
[0022] 步驟10)對單炮地震記錄X" ' = (x〃 ' υ)πΧη?行偏移距幾何反變換，得到壓制 地滾波后的最終單炮地震記錄X" " = (X" " Jmxn，其中X" " ijS單炮地震記錄X"" 的第i個地震道、第j個采樣點，i = 1，2,. . .，m ;j = 1，2,. . .，n，完成了地滾波的壓制。
[0023] 步驟1)所述的對單炮地震記錄X = (Xij)mxn進行偏移距幾何變換，按如下步驟進 行：
[0024] (1)拾取單炮地震記錄X = (Xij)mxn的初至時間FBT i，i = 1，…，m ;
[0025] (2)由單炮地震記錄X = (Xij)mxn的初至時間FBT i，生成炮檢距呈線性變化的單炮 地震記錄初至時間FBT' i，i = 1，· · ·，m ;
[0026] (3)對單炮地震記錄X= (Xij)mxJt地震道時移校正，變換為t = (X' Jmxn，時 移量為 ATi=FBTi-FBIw i。
[0027] 步驟4)中所述的由G= (gij)mXn得到G' = (g' Jmxn具體步驟如下：
[0028] (1)對單炮地震記錄G = (gij)mXn采用步驟2)得到的視速度做線性動校正，使單 炮地震記錄G = (gij)mXn上的地滾波同相軸排齊；
[0029] (2)單炮地震記錄G= (Sij)nixn上的地滾波同相軸排齊后，再結合互相關確定地震 道間的靜校正時移量Δ τ i，i = 1，2, ...，m，用靜校正時移量Λ τ 齊拉平地滾波，得到 單炮地震記錄C = (g'
[0030] 步驟5)中求單炮地震記錄G' = (g' υ)πΧη的正交投影向量矩陣U= (Uij)mxi的 具體步驟如下：
[0031] (1)對G'進行主成分分析，求G'的協方差矩陣Ce, =E[G' (G' )Τ]，其中E表 示求數學期望，Ce, e Rnxn為實對稱矩陣，（G" ) τ為G"的轉置矩陣；
[0032] (2)對心進行奇異值分解，得到奇異值（^=〇' = 1，...，111)，其中1'為矩陣(^ 的奇異值個數，則Ce,的特征值矩陣即對角陣
，e Rmxm，必然存在 一個正交特征向量矩陣W =(u，JmxmW eRmXm，滿足（U，）TCe, W =A、（W )1? U'的轉置矩陣；
[0033] (3)求正交投影向量矩陣U = (Uij)mXi，用地滾波第一個主分量代表地滾波，根據 上步的推導，則r = 1，對應對角陣為A = σ pA e R1X1，則對應得到正交投影向量矩陣U = (Uij) mX 1。
[0034] 步驟7)所述的由G" = (g" Jmxn得到G" = (g" Jmxn具體步驟如下：
[0035] (1)將 G" = (g" ij^xn進行靜校正量反時移-Δ τ p i = 1，2, · · ·，m ;
[0036] (2)在上步的基礎上，完成線性反動校，得到地滾波最終模型G" '= (g" ' y mXn0
[0037] 步驟10)所述對單炮地震記錄X" '= (X" ' Jmxn進行偏移距幾何反變換，時移 量為-ATi=-(FBTi-FBT' 得到壓制地滾波后的最終單炮地震記錄X" " = (X" " y mXn0
[0038] 本發(fā)明與現有技術相比具有以下優(yōu)點：
[0039] (1)現有技術針對疊前疊后傾斜相干噪聲，采用了道間時差排齊同相軸的傾斜 K-L變換，壓制傾斜噪音能取得較好的效果，但不適合復雜地區(qū)非縱測線上雙曲線化地滾波 的壓制；本發(fā)明采用偏移距幾何變換法將雙曲線化地滾波線性化，滿足了傾斜K-L變換的 適用條件，避免了用直線段代替曲線弧的分區(qū)分塊的傾斜K-L變換產生的誤差，提高了地 滾波模型預測的精度。使得噪音去除更徹底，從而突出有效信號能量。
[0040] (2)與現有技術相比，本方法針對地滾波在頻率、相位、速度、振幅以及空間假頻上 具有局部變化快的特點，本發(fā)明采用了分頻處理，自動調節(jié)追蹤地滾波在不同地質條件下 的變化特征，從而實現波場分離與去噪，能夠提高中生界致密油目的層的連續(xù)性，因此該方 法具有振幅保真特性。
[0041] 下面將結合附圖做進一步詳細說明。
【附圖說明】
[0042] 圖1是本發(fā)明的流程示意圖；
[0043] 圖2 (a)為去噪前的非縱測線單炮記錄；
[0044] 圖2 (b)為非線性地滾波采用偏移距幾何變換法線性化后的單炮記錄；
[0045] 圖2 (c)利用K-L變換壓制線性地滾波后的單炮記錄；
[0046] 圖2 (d)為非線性地滾波壓制后再經偏移距幾何反變換后的單炮記錄；
[0047] 圖3(a)為鄂爾多斯盆地隴東地區(qū)黃土塬非縱測線單炮記錄去噪前的非縱測線單 炮記錄；
[0048] 圖3 (b)為現有方法采用K-L變換壓制地滾波后的單炮記錄；
[0049] 圖3 (c)為現有方法采用K-L變換壓制地滾波后的噪音單炮記錄；
[0050] 圖4為采用本發(fā)明對鄂爾多斯盆地隴東地區(qū)黃土塬非縱測線單炮記錄處理后的 效果；
[0051] 圖4(a)利用偏移距幾何變換法結合K-L變換壓制地滾波后的單炮記錄；
[0052] 圖4 (b)為利用偏移距幾何變換法結合K-L變換壓制地滾波后的噪音單炮記錄；
[0053] 圖5 (a)為去噪前的非縱測線單炮記錄頻譜；
[0054] 圖5 (b)為現有方法采用K-L變換壓制地滾波后的單炮記錄頻譜；
[0055] 圖5 (c)利用本發(fā)明壓制地滾波后的單炮記錄頻譜；
[0056] 圖6 (a)為現有方法采用K-L變換壓制地滾波后的疊加剖面；
[0057] 圖6(b)為采用本發(fā)明壓制地滾波后的疊加剖面。
【具體實施方式】
[0058] 實施例1 :
[0059] 鄂爾多斯盆地巨厚黃土塬區(qū)針對中生界致密油目的層實施的非縱測線采集方式， 進一步提升了地震資料的品質，但單炮記錄上仍存在較強的地滾波，且600米的非縱距使 得地滾波呈現雙曲線化的噪音形態(tài)，使得K-L變換壓制地滾波的適用條件受到限制。本發(fā) 明采用偏移距幾何變換法，對單炮記錄上的雙曲線化地滾波變換為隨偏移距線性變化的地 滾波，然后采用K-L變頻奇異值分解去除近炮點的頻散地滾波，再進行幾何反變換將其返 回至原始狀態(tài)，從而達到在不傷害有效信號的前提下，有效去除地滾波的目的。
[0060] 本實施例提供了一種如圖1所示的地滾波壓制方法，包括如下步驟：
[0061] 步驟1)對單炮地震記