基于疊后反射波能量統計的地表一致性振幅補償方法
【專利摘要】本發明提供了一種基于多域疊加反射波能量統計的地震資料的地表一致性振幅補償方法，利用共炮、共檢、共偏移距等多域疊加后的反射波能量統計，通過迭代得到地表一致性振幅補償因子，應用于疊前地震記錄，能有效的消除低信噪比、復雜地表的條件下，振幅受地表非一致因素的影響，方法實現簡單、可靠且適應性強。
【專利說明】基于疊后反射波能量統計的地表一致性振幅補償方法
【技術領域】
[0001]本發明涉及地震資料處理技術，是一種基于多域疊加反射波能量統計的地震資料的地表一致性振幅補償方法。
【背景技術】
[0002]地表的不一致性是造成復雜地區地震資料品質降低的重要原因之一。通常激發條件(包括藥量)和接收條件的差異會導致各道反射能量不同，需要經地表一致性處理的目的就是消除地表(包括近地表)的不一致性對地震記錄的影響。地表一致性振幅補償要消除的是地表的不一致性引起的能量變化，是一個相對振幅保持振幅補償方法，對某些資料例如海洋資料不是必要的，但對于復雜地區的地震資料則是非做不可。地表一致性振幅補償是泰勒等人(1981)提出的一種振幅校正技術，將地表條件變化而引起的反射記錄崎變歸結為炮點響應、接收點響應、炮檢距響應及共中心點響應的綜合反映，并提出地表一致性校正模型。其基本思想是把地震波的振幅分解為與地面及地表有關的分量和與界面有關的分量，從而消除地表及近地表因素對振幅的影響。地表一致性振幅補償的模型很多，但都是在疊前道集中求出全測線各炮集、各接收點集、各炮檢距道集上的補償系數，而對于信噪比低的地震資料，在疊前準確的統計出各個補償系數有些困難，相對疊前反射信號，疊后反射同相軸有更高的信噪比，因此，我們提出了基于疊后反射波能量統計的地表一致性振幅補償方法。
【發明內容】
[0003]本發明是提供了一種基于多域疊后反射波能量統計的，能有效減小復雜地表條件下激發和接收等因素對振幅影響的低信噪比地震資料的地表一致性振幅補償方法。
[0004]本發明通過以下技術方案實現:
[0005]一種基于疊后反射波能量統計的地表一致性振幅補償方法，具體步驟包括:
[0006]I)在野外采集地震數據，經預處理和速度分析，定義振幅分析時窗，地震數據(DATAl)；
[0007]步驟I)所述的預處理包括加載觀測系統，幾何擴散振幅補償，靜校正。
[0008]步驟I)所述的速度分析是通過疊加速度分析分別由縱波資料或轉換波資料得到縱波均方根速度或轉換波等效均方根速度。
[0009]步驟I)所述的定義振幅分析時窗是指按一個層位的時間(T0)，長度為A T秒的時窗，這里的時間(TO)是零炮間距時間。
[0010]2)確定對振幅影響由大到小的數據類型，振幅門檻值(ATH)，迭代次數(M)；
[0011]步驟2)所述的數據類型是炮點，檢波點，中心點，偏移距。
[0012]步驟2)所述的對振幅影響由大到小的數據類型，是炮號、檢波點號、CMP號和炮檢距中的N個選項，可以是I個選項，2個選項，3個選項，最多4個選項，兩項以上，數據類型
不重復。[0013]步驟2)所述的對振幅影響由大到小的數據類型N個選項是這樣確定的:對振幅影響最大的數據類型選項，即第I選項，按道頭字選擇，缺省為炮集數據，即首先計算炮集中各道的平均絕對振幅或均方根振幅。選擇數據類型時，要考慮其地表不一致性對振幅的影響，哪種類型數據對振幅的影響最大，就首先選哪種類型數據，可選擇炮號(shot number,缺省)、檢波點號(receiver number)、CMP號(CMP number)和炮檢距(offset);對振幅影響次大的數據類型選項，即第2選項,可選擇炮號(shot number)、檢波點號(receivernumber,缺省)、CMP號(CMP number)和炮檢距(offset);對振幅影響占第三位的數據類型選項，即第3選項，可選擇炮號(shot number)、檢波點號(receiver number)、CMP號(CMPnumber,缺省)和炮檢距(offset);對振幅影響占第四位的數據類型選項，即第4選項,可選擇炮號(shot number)、檢波點號(receiver number)、CMP 號(CMP number)和炮檢距(offset，缺省)；選項順序用符號I代替，1=1，2，-N；
[0014]3)按第I選項道集輸入地震數據，根據選定振幅分析時窗時間，確定常速動校正速度，并根據常速動校正速度確定偏移距(offset)的時窗的時間，之后在選定時窗范圍進行第I選項相干疊加，計算均方根振幅或平均絕對振幅(A1),并將時窗范圍內的數據
(DATA2)存儲起來，所有的第I選項完成后可以得到全工區平均振幅(3 第I選項每個點
的振幅(A1)與全工區平均振幅(3 )之比，就得到全工區各第I選項點的補償系數(C1)；
[0015]4)讀取(I)選定時窗范圍內存儲的數據(DATA2)，依次按第I選項道集抽取數據，對每道應用其補償系數C1-Ch進行共第I選項點的相干疊加，計算均方根振幅或平均絕對
振幅(A1),每個第I選項點的振幅(A1)與全工區平均振幅(i)之比，就得到全工區各第I
選項點的補償系數(C1);
[0016]步驟4)所述的依次按第I選項是指按1=2，…，N選項，先從1=2，直到I=N ；
[0017]步驟4)所述的對每道應用其補償系數C1-Ch，是指應用第I選項順序之前的其它選項的補償系數，比如1=3，每道應用其補償系數CpC2，比如1=4，每道應用其補償系數Cpc2、c3
[0018]5)讀取(I)選定時窗范圍內存儲的數據(DATA2)，依次按第I選項的道集抽取數據，對每道應用其余N-1項的補償系數進行第I選項點相干疊加，計算均方根振幅或平均絕
對振幅(A1),每個第I選項點的振幅(A1)與全工區平均振幅(i )之比，就得到全測線新的
第I選項點的補償系數(C1);
[0019]步驟5)所述的依次按第I選項的道集抽取數據，是指按1=1，…，N選項，先從1=1，直到I=N ；
[0020]步驟5)所述的對每道應用其余N-1項的補償系數，是指對每道應用其對應的補償系數C1,…，C1^1,…，CI+1，Cn ；比如1=3，每道應用其補償系數Cp C2，C4 ;
[0021]步驟5)所述的進行第I選項相干疊加，是指第I選項對應的每個道集內的道進行相干疊加，每個道集內的道計算它的振幅AA，計算模型振幅Am，若AA>Am*ATH，則該道不參加該道集的相干疊加；
[0022]步驟5)所述的模型振幅
【權利要求】
1.一種基于疊后反射波能量統計的地表一致性振幅補償方法，具體步驟包括: 1)在野外采集地震數據，經預處理和速度分析，定義振幅分析時窗、地震數據DATA1， 步驟I)所述的預處理包括加載觀測系統，幾何擴散振幅補償，靜校正， 步驟I)所述的速度分析是通過疊加速度分析分別由縱波資料或轉換波資料得到縱波均方根速度或轉換波等效均方根速度， 步驟I)所述的定義振幅分析時窗是指按一個層位的時間TO、長度為A T秒的時窗，這里的時間TO是零炮間距時間； 2)確定對振幅影響由大到小的數據類型、振幅門檻值ATH、迭代次數M， 步驟2)所述的對振幅影響由大到小的數據類型是炮號、檢波點號、CMP號和炮檢距中的至少兩個選項，所述選項用符號I代替，其按對振幅影響由大到小順序依次編號為第1、.2、...N 選項，SP 1=1，2，…N，N ≤ 4 ; 3)按第I選項的道集輸入地震數據，根據選定振幅分析時窗時間，確定常速動校正速度，并根據常速動校正速度確定炮檢距的時窗的時間，之后在選定時窗范圍進行第I選項相干疊加，計算均方根振幅或平均絕對振幅A1,并將時窗范圍內的數據DATA2存儲起來，所有的第I選項完成后可以得到全工區平均振幅2 ,第I選項每個點的振幅A1與全工區平均振幅:之比，就得到全工區各第I選項點的補償系數C1 ； 4)讀取選定時窗范圍內存儲的數據DATA2，依次按第I選項的道集抽取數據，對每道應用其補償系數C1至CV1進行共第I選項點的相干疊加，計算均方根振幅或平均絕對振幅A1,每個第I選項點的振幅A1與全工區平均振幅:比，就得到全工區各第I選項點的補償系數C1, 步驟4)所述的依次按第I選項是指按1=2，…,N選項順序,先從1=2開始，直到I=N ； 5)讀取選定時窗范圍內存儲的數據DATA2，依次按第I選項的道集抽取數據，對每道應用其余N-1項的補償系數進行第I選項點相干疊加，計算均方根振幅或平均絕對振幅A1，每個第I選項點的振幅A1與全工區平均振幅I之比，就得到全測線新的第I選項點的補償系數C1, 步驟5)所述的依次按第I選項的道集抽取數據，是指按1=1，…，N選項順序，先從I=I開始，直到I=N ； 步驟5)所述的對每道應用其余N-1項的補償系數，是指對每道應用其對應的補償系數。，…，Ch,…，CI+1，CN， 步驟5)所述的進行第I選項點相干疊加，是指第I選項點對應的每個道集內的道進行相干疊加，每個道集內的道計算其振幅AA，計算模型振幅Am，若AA>Am*ATH，則該道不參加該道集的相干疊加，所述的模型振幅
【文檔編號】G01V1/28GK103645499SQ201310553765
【公開日】2014年3月19日 申請日期:2013年11月8日 優先權日:2013年11月8日
【發明者】陳海峰, 李向陽, 狄幫讓, 魏建新, 蔡志光, 陳雙全 申請人:中國石油大學(北京)
專業數控銑床產品供應商
20年生產經驗，實力工廠，精工品質，質量有保障
