一種基于區域約束的分頻段能量地震面波處理方法
【專利摘要】本發明提供了一種基于區域約束的分頻段能量地震面波處理方法，屬于油氣地震勘探領域。本方法包括：(1)輸入地震數據，所述地震數據是需要處理的地震采集單炮數據；(2)設計時窗；(3)在步驟(2)設計的時窗內對所述地震數據進行分頻處理；(4)在不同時窗、不同的頻帶內計算時窗內的平均振幅值；(5)設計振幅門檻值；(6)強幅分析；(7)面波區域約束判別與處理；(8)效果分析，如果合適，則轉入步驟(9)，如果不合適，則返回步驟(5)；(9)將經過處理后的地震數據再通過小波反變換后形成的數據輸出。本發明有效的保護了全局低頻信息，同時在面波區域內更加有效的消除面波，提高了信噪比。
【專利說明】一種基于區域約束的分頻段能量地震面波處理方法
【技術領域】
[0001] 本發明屬于油氣地震勘探領域，具體涉及一種基于區域約束的分頻段能量地震面 波處理方法。
【背景技術】
[0002] 面波是陸上地震勘探中最常見的干擾波之一，與有效信號相比它的能量往往強幾 倍甚至幾十倍，視速度低能量強，頻率主要分布在IOHz以內。隨著傳播距離增大面波振幅 延續時間變長，產生頻散形成掃帚狀分布。面波壓制處理方法的選擇和應用對地震資料的 處理有很大的影響。
[0003] 地震資料面波處理方法主要有：
[0004] 1、帶通濾波法。依據面波的低頻在IOHz以內，設計一個高通濾波器，濾波器的低 頻為10Hz，對數據進行濾波處理，此時高于這個頻率的成份通過，低于這個頻率的成份被濾 掉了，實現面波處理。這一方法簡單，但方法本身在濾去面波的同時，也將低于IOHz的有效 成分給濾掉了，損失了一些有用的信息。
[0005] 2、FXCNS壓制面波法。FXCNS預測相干噪音是在給定頻率的數據上進行的，根據最 小平方誤差準則，分別估計每個頻率的相干噪音。利用頻率和速度兩個參數限制噪音范圍， 對指定范圍以外的信號影響較小，能壓制空間采樣不規則噪音，適用范圍廣。缺點是在垂直 測線方向接收點非常稀疏，不利于相干噪音預測和壓制。
[0006] 3、FK壓制面波方法。該方法是對地震數據進行二維傅里葉變換得到F-K譜，在這 個譜內面波低頻與有效波是分開的，根據這一特性在F-K譜內將面波譜能量進行消除處理 后，再進行二維傅里葉反變換得到消除面波的數據。當數據采樣過大時FK變換易產生假 頻，反變換后信號產生失真。
[0007] 4、自適應面波衰減法。利用時頻分析的方法，根據面波和反射波在頻率分布特征、 空間分布范圍、能量等方面的差異，首先檢測出面波在時間和空間上的分布范圍，然后根據 面波的頻率特征，對其進行分頻加權壓制。缺陷是單炮之間的主頻存在較大差異時，應采用 不同的參數進行壓制，否則對有效信號會產生較大損害。
[0008] 5、基于頻率衰減特性面波壓制法。在地震記錄上確定面波帶，針對不同的炮檢距 選用不同的低切頻，在面波帶里進行高通濾波。面波的壓制在面波帶里進行，能夠較有效地 衰減面波帶里的面波能量。該方法實現起來難度很大形成不了生產能力。
[0009] 以上這些方法具有各自的優缺點，基于區域約束的分頻段能量地震面波處理方法 更具優勢，能夠擬補其它方法的不足，并未見相應的研究文獻。
[0010] 目前消除面波的方法大多都沒有考慮在去掉面波的同時也要保留面波頻率內的 有效地震成分。多數方法對有效信號產生一些改變或被消除，不利于后期的精確處理成像。 尤其是地震資料處理的重點由疊后轉到疊前，疊前AVO處理及屬性分析對地震數據的保幅 性要求高，在面波壓制過程中確保面波去除的同時，盡可能不損害有效信息，為后續的保幅 處理、屬性提取等巖性研究工作打下良好資料基礎。
【發明內容】
[0011] 本發明的目的在于解決上述現有技術中存在的難題，提供一種基于區域約束的分 頻段能量地震面波處理方法，能夠在有效消除地震面波的同時保留面波頻率內的有效地震 成分，采取區域約束機制不破壞面波區域外的接近面波頻率的有效信息，利于后期的成像 精度。
[0012] 本發明是通過以下技術方案實現的：
[0013] 一種基于區域約束的分頻段能量地震面波處理方法，包括：
[0014] (1)輸入地震數據，所述地震數據是需要處理的地震采集單炮數據；
[0015] (2)設計時窗；
[0016] (3)在步驟（2)設計的時窗內對所述地震數據進行分頻處理；
[0017] (4)在不同時窗、不同的頻帶內計算時窗內的平均振幅值；
[0018] (5)設計振幅門檻值；
[0019] (6)強幅分析；
[0020] (7)面波區域約束判別與處理；
[0021] (8)效果分析，如果合適，則轉入步驟（9)，如果不合適，則返回步驟（5);
[0022] (9)將經過處理后的地震數據再通過小波反變換后形成的數據輸出。
[0023] 所述步驟（2)是這樣實現的：
[0024] 根據面波起始時間作為分頻開始時間，設計接近一個反射子波長度的時窗，所述 時窗的長度為400-500ms。
[0025] 所述步驟（3)是這樣實現的：
【權利要求】
1. 一種基于區域約束的分頻段能量地震面波處理方法，其特征在于：所述方法包括： (1) 輸入地震數據，所述地震數據是需要處理的地震采集單炮數據； (2) 設計時窗； (3) 在步驟（2)設計的時窗內對所述地震數據進行分頻處理； (4) 在不同時窗、不同的頻帶內計算時窗內的平均振幅值； (5) 設計振幅門檻值； (6) 強幅分析； (7) 面波區域約束判別與處理； (8) 效果分析，如果合適，則轉入步驟（9)，如果不合適，則返回步驟（5); (9) 將經過處理后的地震數據再通過小波反變換后形成的數據輸出。
2. 根據權利要求1所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述步驟（2)是這樣實現的： 根據面波起始時間作為分頻開始時間，設計接近一個反射子波長度的時窗，所述時窗 的長度為400-500ms。
3. 根據權利要求2所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述步驟（3)是這樣實現的： 在時窗內利用下面的小波變換公式對地震數據進行分頻處理，將低頻面波變換出來： Wf(s,u) = ￡f(x)4sy/(s(x-u))dx ,
4. 根據權利要求3所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述步驟（4)是這樣實現的： 利用下式進行振幅強度分析，求取時窗內的平均振幅：
上式中mftk為時窗內所含有的道數，Eftk是頻帶f時窗t內k道的平均振幅，Aiftk是頻 帶f?時窗t內k道第i個采樣點的振幅。
5. 根據權利要求4所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述步驟（5)是這樣實現的： 利用平均振幅對所有采樣值進行統計分析，選取能夠包含80%采樣點數的振幅值作為 振幅門檻值；所述統計分析是采用累計疊加統計分析法，對采樣點逐點進行統計來實現的。
6. 根據權利要求5所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述步驟（6)是這樣實現的： 首先在時窗內進行平均振幅計算，利用平均振幅對所有采樣值進行統計分析，選取能 夠包含80%采樣點數的振幅值作為振幅門檻值，將大于這個門檻值的振幅作為強振幅。
7. 根據權利要求6所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述步驟（7)是這樣實現的： 分析步驟（1)輸入的地震數據中的單炮，確定能夠代表絕大多數的單炮，對該單炮進 行分析得到其面波分布，然后設計一個包含其面波的區域，將該區域作為約束區域； 在所述時窗內逐點提取振幅強度值，根據振幅強度分別采用壓制因子公式計算，在所 述約束區域的約束下對地震數據進行處理更新，具體如下： 對超過振幅門檻值的樣點，通過壓制因子公式計算得到新的振幅值，并用該新的振幅 值替換原來的振幅值。
8. 根據權利要求7所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述壓制因子公式如下： 設頻帶f?時窗t內的壓制因子Hftk為：
式中：ct為壓制系數，Rftk為頻帶f時窗t內的振幅門濫值。
9. 根據權利要求8所述的基于區域約束的分頻段能量地震面波處理方法，其特征在 于：所述步驟（8)中的所述效果分析具體如下： 對經過步驟（7)處理后的面波，在面波區域內去掉面波后，如果能夠保留有效波，并且 區域內的有效波與周邊的有效波特征基本一致，則判斷為合適，否則為不合適；所述有效波 特征包括頻率、振幅及雙曲特征。
10. 根據權利要求9所述的基于區域約束的分頻段能量地震面波處理方法，其特征 在于：所述步驟（9)中的小波反變換具體如下：采用下面的小波反變換公式對數據進行處 理：
上式中：f(x)為輸入的地震數據，V (X)為基本小波，a為尺度因子，b為平移因子，s 為尺度因子的倒數，u等于b, Wf (s, u)為變換后的小波，f (t)為變換后的地震數據，C為矩 陣。
【文檔編號】G01V1/28GK104345341SQ201310342646
【公開日】2015年2月11日 申請日期:2013年8月8日 優先權日:2013年8月8日
【發明者】郭愷 申請人:中國石油化工股份有限公司, 中國石油化工股份有限公司石油物探技術研究院