基于疊后反射波能量統(tǒng)計(jì)的地表一致性振幅補(bǔ)償方法
【專利摘要】本發(fā)明提供了一種基于多域疊加反射波能量統(tǒng)計(jì)的地震資料的地表一致性振幅補(bǔ)償方法，利用共炮、共檢、共偏移距等多域疊加后的反射波能量統(tǒng)計(jì)，通過迭代得到地表一致性振幅補(bǔ)償因子，應(yīng)用于疊前地震記錄，能有效的消除低信噪比、復(fù)雜地表的條件下，振幅受地表非一致因素的影響，方法實(shí)現(xiàn)簡單、可靠且適應(yīng)性強(qiáng)。
【專利說明】基于疊后反射波能量統(tǒng)計(jì)的地表一致性振幅補(bǔ)償方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地震資料處理技術(shù)，是一種基于多域疊加反射波能量統(tǒng)計(jì)的地震資料的地表一致性振幅補(bǔ)償方法。
【背景技術(shù)】
[0002]地表的不一致性是造成復(fù)雜地區(qū)地震資料品質(zhì)降低的重要原因之一。通常激發(fā)條件(包括藥量)和接收條件的差異會導(dǎo)致各道反射能量不同，需要經(jīng)地表一致性處理的目的就是消除地表(包括近地表)的不一致性對地震記錄的影響。地表一致性振幅補(bǔ)償要消除的是地表的不一致性引起的能量變化，是一個(gè)相對振幅保持振幅補(bǔ)償方法，對某些資料例如海洋資料不是必要的，但對于復(fù)雜地區(qū)的地震資料則是非做不可。地表一致性振幅補(bǔ)償是泰勒等人(1981)提出的一種振幅校正技術(shù)，將地表?xiàng)l件變化而引起的反射記錄崎變歸結(jié)為炮點(diǎn)響應(yīng)、接收點(diǎn)響應(yīng)、炮檢距響應(yīng)及共中心點(diǎn)響應(yīng)的綜合反映，并提出地表一致性校正模型。其基本思想是把地震波的振幅分解為與地面及地表有關(guān)的分量和與界面有關(guān)的分量，從而消除地表及近地表因素對振幅的影響。地表一致性振幅補(bǔ)償?shù)哪Ｐ秃芏?，但都是在疊前道集中求出全測線各炮集、各接收點(diǎn)集、各炮檢距道集上的補(bǔ)償系數(shù)，而對于信噪比低的地震資料，在疊前準(zhǔn)確的統(tǒng)計(jì)出各個(gè)補(bǔ)償系數(shù)有些困難，相對疊前反射信號，疊后反射同相軸有更高的信噪比，因此，我們提出了基于疊后反射波能量統(tǒng)計(jì)的地表一致性振幅補(bǔ)償方法。

【發(fā)明內(nèi)容】

[0003]本發(fā)明是提供了一種基于多域疊后反射波能量統(tǒng)計(jì)的，能有效減小復(fù)雜地表?xiàng)l件下激發(fā)和接收等因素對振幅影響的低信噪比地震資料的地表一致性振幅補(bǔ)償方法。
[0004]本發(fā)明通過以下技術(shù)方案實(shí)現(xiàn):
[0005]一種基于疊后反射波能量統(tǒng)計(jì)的地表一致性振幅補(bǔ)償方法，具體步驟包括:
[0006]I)在野外采集地震數(shù)據(jù)，經(jīng)預(yù)處理和速度分析，定義振幅分析時(shí)窗，地震數(shù)據(jù)(DATAl)；
[0007]步驟I)所述的預(yù)處理包括加載觀測系統(tǒng)，幾何擴(kuò)散振幅補(bǔ)償，靜校正。
[0008]步驟I)所述的速度分析是通過疊加速度分析分別由縱波資料或轉(zhuǎn)換波資料得到縱波均方根速度或轉(zhuǎn)換波等效均方根速度。
[0009]步驟I)所述的定義振幅分析時(shí)窗是指按一個(gè)層位的時(shí)間(T0)，長度為A T秒的時(shí)窗，這里的時(shí)間(TO)是零炮間距時(shí)間。
[0010]2)確定對振幅影響由大到小的數(shù)據(jù)類型，振幅門檻值(ATH)，迭代次數(shù)(M)；
[0011]步驟2)所述的數(shù)據(jù)類型是炮點(diǎn)，檢波點(diǎn)，中心點(diǎn)，偏移距。
[0012]步驟2)所述的對振幅影響由大到小的數(shù)據(jù)類型，是炮號、檢波點(diǎn)號、CMP號和炮檢距中的N個(gè)選項(xiàng)，可以是I個(gè)選項(xiàng)，2個(gè)選項(xiàng)，3個(gè)選項(xiàng)，最多4個(gè)選項(xiàng)，兩項(xiàng)以上，數(shù)據(jù)類型
不重復(fù)。[0013]步驟2)所述的對振幅影響由大到小的數(shù)據(jù)類型N個(gè)選項(xiàng)是這樣確定的:對振幅影響最大的數(shù)據(jù)類型選項(xiàng)，即第I選項(xiàng)，按道頭字選擇，缺省為炮集數(shù)據(jù)，即首先計(jì)算炮集中各道的平均絕對振幅或均方根振幅。選擇數(shù)據(jù)類型時(shí)，要考慮其地表不一致性對振幅的影響，哪種類型數(shù)據(jù)對振幅的影響最大，就首先選哪種類型數(shù)據(jù)，可選擇炮號(shot number,缺省)、檢波點(diǎn)號(receiver number)、CMP號(CMP number)和炮檢距(offset);對振幅影響次大的數(shù)據(jù)類型選項(xiàng)，即第2選項(xiàng),可選擇炮號(shot number)、檢波點(diǎn)號(receivernumber,缺省)、CMP號(CMP number)和炮檢距(offset);對振幅影響占第三位的數(shù)據(jù)類型選項(xiàng)，即第3選項(xiàng)，可選擇炮號(shot number)、檢波點(diǎn)號(receiver number)、CMP號(CMPnumber,缺省)和炮檢距(offset);對振幅影響占第四位的數(shù)據(jù)類型選項(xiàng)，即第4選項(xiàng),可選擇炮號(shot number)、檢波點(diǎn)號(receiver number)、CMP 號(CMP number)和炮檢距(offset，缺省)；選項(xiàng)順序用符號I代替，1=1，2，-N；
[0014]3)按第I選項(xiàng)道集輸入地震數(shù)據(jù)，根據(jù)選定振幅分析時(shí)窗時(shí)間，確定常速動校正速度，并根據(jù)常速動校正速度確定偏移距(offset)的時(shí)窗的時(shí)間，之后在選定時(shí)窗范圍進(jìn)行第I選項(xiàng)相干疊加，計(jì)算均方根振幅或平均絕對振幅(A1),并將時(shí)窗范圍內(nèi)的數(shù)據(jù)
(DATA2)存儲起來，所有的第I選項(xiàng)完成后可以得到全工區(qū)平均振幅(3 第I選項(xiàng)每個(gè)點(diǎn)
的振幅(A1)與全工區(qū)平均振幅(3 )之比，就得到全工區(qū)各第I選項(xiàng)點(diǎn)的補(bǔ)償系數(shù)(C1)；
[0015]4)讀取(I)選定時(shí)窗范圍內(nèi)存儲的數(shù)據(jù)(DATA2)，依次按第I選項(xiàng)道集抽取數(shù)據(jù)，對每道應(yīng)用其補(bǔ)償系數(shù)C1-Ch進(jìn)行共第I選項(xiàng)點(diǎn)的相干疊加，計(jì)算均方根振幅或平均絕對
振幅(A1),每個(gè)第I選項(xiàng)點(diǎn)的振幅(A1)與全工區(qū)平均振幅(i)之比，就得到全工區(qū)各第I
選項(xiàng)點(diǎn)的補(bǔ)償系數(shù)(C1);
[0016]步驟4)所述的依次按第I選項(xiàng)是指按1=2，…，N選項(xiàng)，先從1=2，直到I=N ；
[0017]步驟4)所述的對每道應(yīng)用其補(bǔ)償系數(shù)C1-Ch，是指應(yīng)用第I選項(xiàng)順序之前的其它選項(xiàng)的補(bǔ)償系數(shù)，比如1=3，每道應(yīng)用其補(bǔ)償系數(shù)CpC2，比如1=4，每道應(yīng)用其補(bǔ)償系數(shù)Cpc2、c3
[0018]5)讀取(I)選定時(shí)窗范圍內(nèi)存儲的數(shù)據(jù)(DATA2)，依次按第I選項(xiàng)的道集抽取數(shù)據(jù)，對每道應(yīng)用其余N-1項(xiàng)的補(bǔ)償系數(shù)進(jìn)行第I選項(xiàng)點(diǎn)相干疊加，計(jì)算均方根振幅或平均絕
對振幅(A1),每個(gè)第I選項(xiàng)點(diǎn)的振幅(A1)與全工區(qū)平均振幅(i )之比，就得到全測線新的
第I選項(xiàng)點(diǎn)的補(bǔ)償系數(shù)(C1);
[0019]步驟5)所述的依次按第I選項(xiàng)的道集抽取數(shù)據(jù)，是指按1=1，…，N選項(xiàng)，先從1=1，直到I=N ；
[0020]步驟5)所述的對每道應(yīng)用其余N-1項(xiàng)的補(bǔ)償系數(shù)，是指對每道應(yīng)用其對應(yīng)的補(bǔ)償系數(shù)C1,…，C1^1,…，CI+1，Cn ；比如1=3，每道應(yīng)用其補(bǔ)償系數(shù)Cp C2，C4 ;
[0021]步驟5)所述的進(jìn)行第I選項(xiàng)相干疊加，是指第I選項(xiàng)對應(yīng)的每個(gè)道集內(nèi)的道進(jìn)行相干疊加，每個(gè)道集內(nèi)的道計(jì)算它的振幅AA，計(jì)算模型振幅Am，若AA>Am*ATH，則該道不參加該道集的相干疊加；
[0022]步驟5)所述的模型振幅
【權(quán)利要求】
1.一種基于疊后反射波能量統(tǒng)計(jì)的地表一致性振幅補(bǔ)償方法，具體步驟包括: 1)在野外采集地震數(shù)據(jù)，經(jīng)預(yù)處理和速度分析，定義振幅分析時(shí)窗、地震數(shù)據(jù)DATA1， 步驟I)所述的預(yù)處理包括加載觀測系統(tǒng)，幾何擴(kuò)散振幅補(bǔ)償，靜校正， 步驟I)所述的速度分析是通過疊加速度分析分別由縱波資料或轉(zhuǎn)換波資料得到縱波均方根速度或轉(zhuǎn)換波等效均方根速度， 步驟I)所述的定義振幅分析時(shí)窗是指按一個(gè)層位的時(shí)間TO、長度為A T秒的時(shí)窗，這里的時(shí)間TO是零炮間距時(shí)間； 2)確定對振幅影響由大到小的數(shù)據(jù)類型、振幅門檻值A(chǔ)TH、迭代次數(shù)M， 步驟2)所述的對振幅影響由大到小的數(shù)據(jù)類型是炮號、檢波點(diǎn)號、CMP號和炮檢距中的至少兩個(gè)選項(xiàng)，所述選項(xiàng)用符號I代替，其按對振幅影響由大到小順序依次編號為第1、.2、...N 選項(xiàng)，SP 1=1，2，…N，N ≤ 4 ; 3)按第I選項(xiàng)的道集輸入地震數(shù)據(jù)，根據(jù)選定振幅分析時(shí)窗時(shí)間，確定常速動校正速度，并根據(jù)常速動校正速度確定炮檢距的時(shí)窗的時(shí)間，之后在選定時(shí)窗范圍進(jìn)行第I選項(xiàng)相干疊加，計(jì)算均方根振幅或平均絕對振幅A1,并將時(shí)窗范圍內(nèi)的數(shù)據(jù)DATA2存儲起來，所有的第I選項(xiàng)完成后可以得到全工區(qū)平均振幅2 ,第I選項(xiàng)每個(gè)點(diǎn)的振幅A1與全工區(qū)平均振幅:之比，就得到全工區(qū)各第I選項(xiàng)點(diǎn)的補(bǔ)償系數(shù)C1 ； 4)讀取選定時(shí)窗范圍內(nèi)存儲的數(shù)據(jù)DATA2，依次按第I選項(xiàng)的道集抽取數(shù)據(jù)，對每道應(yīng)用其補(bǔ)償系數(shù)C1至CV1進(jìn)行共第I選項(xiàng)點(diǎn)的相干疊加，計(jì)算均方根振幅或平均絕對振幅A1,每個(gè)第I選項(xiàng)點(diǎn)的振幅A1與全工區(qū)平均振幅:比，就得到全工區(qū)各第I選項(xiàng)點(diǎn)的補(bǔ)償系數(shù)C1, 步驟4)所述的依次按第I選項(xiàng)是指按1=2，…,N選項(xiàng)順序,先從1=2開始，直到I=N ； 5)讀取選定時(shí)窗范圍內(nèi)存儲的數(shù)據(jù)DATA2，依次按第I選項(xiàng)的道集抽取數(shù)據(jù)，對每道應(yīng)用其余N-1項(xiàng)的補(bǔ)償系數(shù)進(jìn)行第I選項(xiàng)點(diǎn)相干疊加，計(jì)算均方根振幅或平均絕對振幅A1，每個(gè)第I選項(xiàng)點(diǎn)的振幅A1與全工區(qū)平均振幅I之比，就得到全測線新的第I選項(xiàng)點(diǎn)的補(bǔ)償系數(shù)C1, 步驟5)所述的依次按第I選項(xiàng)的道集抽取數(shù)據(jù)，是指按1=1，…，N選項(xiàng)順序，先從I=I開始，直到I=N ； 步驟5)所述的對每道應(yīng)用其余N-1項(xiàng)的補(bǔ)償系數(shù)，是指對每道應(yīng)用其對應(yīng)的補(bǔ)償系數(shù)。，…，Ch,…，CI+1，CN， 步驟5)所述的進(jìn)行第I選項(xiàng)點(diǎn)相干疊加，是指第I選項(xiàng)點(diǎn)對應(yīng)的每個(gè)道集內(nèi)的道進(jìn)行相干疊加，每個(gè)道集內(nèi)的道計(jì)算其振幅AA，計(jì)算模型振幅Am，若AA>Am*ATH，則該道不參加該道集的相干疊加，所述的模型振幅
【文檔編號】G01V1/28GK103645499SQ201310553765
【公開日】2014年3月19日 申請日期:2013年11月8日 優(yōu)先權(quán)日:2013年11月8日
【發(fā)明者】陳海峰, 李向陽, 狄?guī)妥? 魏建新, 蔡志光, 陳雙全 申請人:中國石油大學(xué)(北京)