疊后信噪比估算方法及裝置制造方法
【專利摘要】一種疊后信噪比估算方法及裝置，該方法包括：采集地震數據并進行預處理，形成預處理后的地震數據；對預處理后的地震數據進行疊加，形成疊后地震數據；選擇一個時間空間窗，確定時間空間窗的中間地震道；以中間地震道為中心確定時間空間窗中第一地震道的第一樣點，選定N個樣點長的時間窗，并計算時間空間窗中剩余每一地震道的第一樣點；確定時間窗的位置，根據時間窗的位置選擇中間地震道的N個樣點長的數據，根據中間地震道的樣點數據與剩余地震道的樣點數據生成第一樣點的相關比值之和；計算時間空間窗內中間地震道的所有樣點的最大相關值之和及中間地震道的最大相關值總和；計算所述疊后的每個地震道在選擇的時間空間窗內的最大相關值總和。
【專利說明】疊后信噪比估算方法及裝置
【技術領域】
[0001]本發明本方法涉及地球物理勘探數據處理技術，特別是關于是一種疊后信噪比估算方法及裝置。
【背景技術】
[0002]在地震數據處理技術中，噪音常常影響疊加成像效果，噪聲衰減在很多地震數據處理中都變成了常規處理流程，去噪后信噪比的提升一直是一大難題。如何對野外數據采集和室內數據處理進行地震數據信噪比宏觀定性分析和數字顯示，準確地分析資料處理中的各個環節質量，一直是地震勘探研究人員想要實現的目標。
[0003]目前信噪比主要是在指定時間內，用信號與噪聲振幅的比值來衡量信號的強弱的方法。現有技術中，主要信噪比估算方法有:1)時間域疊加法；2)時間域奇異值分減法；3)模型道相關估算法。這三種信噪比估算方法，在地下構造傾角變化較大時，信噪比估算精度不高。時間域疊加和模型道相關主要是因為地下構造傾角變化時同相軸不能準確地同相疊力口。而時間域奇異值分減法本身也不適應地下構造傾角較大變化的地震數據噪聲估計。
【發明內容】
[0004]本發明提供了一種疊后信噪比估算方法及裝置，以精確地震估算地震數據的信噪比。
[0005]為了實現上述目的，本發明提供一種疊后信噪比估算方法，所述的疊后信噪比估算方法包括:
[0006]步驟1:采集地震數據；
[0007]步驟2:對所述地震數據進行預處理，形成預處理后的地震數據；
[0008]步驟3:對所述預處理后的地震數據進行疊加，形成疊后地震數據；
[0009]步驟4:在所述疊后地震數據上選擇一個時間空間窗，確定所述時間空間窗的中間地震道；
[0010]步驟5:以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的第一樣點；
[0011]步驟6:將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置；
[0012]步驟7:根據所述時間窗的位置從所述時間窗中選擇所述中間地震道的N個樣點長的數據，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和；
[0013]步驟8:重復執行步驟5至步驟7，生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和；
[0014]步驟9:重復執行步驟4至步驟8，生成所述疊后的每個地震道在選擇的時間空間窗內的最大相關值總和。
[0015]在一實施例中，對所述地震數據進行預處理，包括:對所述地震數據進行地震數據解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0016]在一實施例中，所述疊后信噪比估算方法還包括:
[0017]判斷是否對所述預處理后的地震數據進行了疊前偏移處理；
[0018]如果否，對所述預處理后的地震數據進行動校正。
[0019]在一實施例中，在所述步驟7中，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和，包括:以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，分別獲得所述剩余地震道的最大相關值，將所述剩余地震道的最大相關值按比例相加，獲得所述相關比值之和。
[0020]在一實施例中，在所述步驟8中，生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和，包括:
[0021]生成所述時間空間窗內中間地震道中除了第一個樣點外的所有樣點的最大互相關值之和；
[0022]將所述時間空間窗內中間地震道中包括第一個樣點的所有樣點的最大互相關值之和相加，形成所述時間空間窗內中間地震道最大互相關值的總和。
[0023]在一實施例中，所述疊后信噪比估算方法還包括:對每個地震道的所述最大相關值總和計算常用對數。
[0024]在一實施例中，在所述步驟4中，選擇所述時間空間窗時，時間窗中包括一個完整的同相軸，空間窗的地震道的個數為奇數，所述空間窗按所述疊后地震數據地震道順序號依次滾動。
[0025]為了實現上述目的，本發明提供一種疊后信噪比估算裝置，所述的疊后信噪比估算裝置包括:
[0026]數據采集單元，用于采集地震數據；
[0027]數據預處理單元，用于對所述地震數據進行預處理，形成預處理后的地震數據；
[0028]數據疊加單元，用于對所述預處理后的地震數據進行疊加，形成疊后地震數據；
[0029]時間窗選擇單元，用于在所述疊后地震數據上選擇一個時間空間窗，確定所述時間空間窗的中間地震道；
[0030]樣點計算單元，用于以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的樣點；
[0031]位置確定單元，用于將所述第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置；
[0032]相關比值生成單元，用于根據所述時間窗的位置從所述時間窗中選擇所述中間地震道的樣點數據，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和；
[0033]總和計算單元，用于生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和。
[0034]在一實施例中，所述數據預處理單元具體用于:對所述地震數據進行地震數據解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0035]在一實施例中，所述的疊后信噪比估算裝置還包括:
[0036]判斷單元，用于判斷是否對所述預處理后的地震數據進行了疊前偏移處理；
[0037]動校正單元，用于對所述預處理后的地震數據進行動校正。
[0038]在一實施例中，所述的相關比值生成單元具體用于:以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，分別獲得所述剩余地震道的最大相關值，將所述剩余地震道的最大相關值按比例相加，獲得所述相關比值之和。
[0039]在一實施例中，所述的總和計算單元具體用于:生成所述時間空間窗內中間地震道中除了第一個樣點外的所有樣點的最大互相關值之和；將所述時間空間窗內中間地震道中包括第一個樣點的所有樣點的最大互相關值之和相加，形成所述時間空間窗內中間地震道最大互相關值的總和。
[0040]在一實施例中，所述的疊后信噪比估算裝置還包括:對數計算單元，用于對每個地震道的所述最大相關值總和計算常用對數。
[0041 ] 在一實施例中，選擇所述時間空間窗中的時間窗中包括一個完整的同相軸，空間窗的地震道的個數為奇數。
[0042]本發明提供的疊后信噪比估算方法，為適應不同地下構造地震數據信噪比估算，更精確地震估算地震數據的信噪比，為定性分析地震數據信噪比的提高提供了有利的工具。
【專利附圖】
【附圖說明】
[0043]為了更清楚地說明本發明實施例或現有技術中的技術方案，下面將對實施例或現有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發明的一些實施例，對于本領域普通技術人員來講，在不付出創造性勞動性的前提下，還可以根據這些附圖獲得其他的附圖。
[0044]圖1為本發明實施例的疊后信噪比估算方法流程圖一；
[0045]圖2為本發明實施例的疊后信噪比估算方法流程圖二 ；
[0046]圖3為本發明實施例的疊后信噪比估算裝置結構框圖一；
[0047]圖4為本發明實施例的疊后信噪比估算裝置結構框圖二 ；
[0048]圖5為本發明實施例的疊后信噪比估算裝置結構框圖三；
[0049]圖6為模型數據和信噪比曲線顯示圖(上:模型數據，下:信噪比曲線)；
[0050]圖7為傾角模型數據和信噪比曲線顯示圖(上:模型數據，下:信噪比曲線)；
[0051]圖8為實際地震數據及信噪比曲線顯示圖(上:模型數據，下:信噪比曲線)；
[0052]圖9為低信噪比實際地震數據及信噪比曲線顯示圖(上:模型數據，下:信噪比曲線)。
【具體實施方式】
[0053]下面將結合本發明實施例中的附圖，對本發明實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本發明一部分實施例，而不是全部的實施例。基于本發明中的實施例，本領域普通技術人員在沒有做出創造性勞動前提下所獲得的所有其他實施例，都屬于本發明保護的范圍。
[0054]如圖1所示，本發明提供一種疊后信噪比估算方法，所述的疊后信噪比估算方法包括:
[0055]SlOl:采集地震數據；
[0056]S102:對所述地震數據進行預處理，形成預處理后的地震數據；
[0057]S103:對所述預處理后的地震數據進行疊加，形成疊后地震數據；
[0058]S104:在所述疊后地震數據上選擇一個時間空間窗，確定所述時間空間窗的中間地震道；
[0059]S105:以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的第一樣點；
[0060]S106:將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置；
[0061]S107:根據所述時間窗的位置從所述時間窗中選擇所述中間地震道的N個樣點長的數據，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和；
[0062]S108:重復執行S105至S107，生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和；
[0063]S109:重復執行S104至S108，生成所述疊后的每個地震道在選擇的時間空間窗內的最大相關值總和。
[0064]由圖1的流程可知，本發明的疊后信噪比估算方法，首先根據疊后地震數據，先選定需要計算靜校正的時間窗和空間窗，計算時間空間窗內的同相軸的傾角，確定相應相關時間窗長度。按計算最大相關值大小的和(按比例)求最大相關值的分貝數，分貝數值的大小直接反應信噪比的高低。本發明與現有技術中的幾種時間域計算信噪比的方法相比，更能適應傾角變化的疊后信噪比估算，信噪比估算精度更可靠。
[0065]S102中，對地震數據進行預處理可以包括:對所述地震數據進行地震數據解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0066]在一實施例中，對所述預處理后的地震數據進行疊加前，還可以判斷是否對預處理后的地震數據進行了疊前偏移處理；如果沒有對預處理后的地震數據進行疊前偏移處理，對預處理后的地震數據進行動校正。
[0067]在S104中，空間窗地震道的選擇一般為奇數，如3，5，7，9，11等。對于確定中間地震道，如果選擇11個地震道的空間窗，那么時間空間窗的中間地震道就是第6個地震道。
[0068]S105具體實施時，以S104選擇出的中間地震道為中心，確定時間空間窗中第一地震道的第一樣點，以第一樣點為中心選定N個樣點長的一段小的時間窗，N為奇數。給定一個道間時間差tl，根據采樣率Si可以計算出11 = tl/si。如果選定11個地震道空間窗，假設第6個地震道的第一個樣點值為k，那么第7個地震道的第一個樣點值就為k-ΙΙ，第8個到11個地震道依次為k-2*ll、k-3*ll、k-4*ll、k-5*ll，同樣第I至5個地震道的第一個樣點為:k+5*ll、k+4*ll、k+3*ll、k+2*ll、k+ll，這樣空間11個地震到的第一個時間窗樣點值就被確定下來。
[0069]需要注意的是，上式中tl與同相軸傾角有關系，一般都用選定時間空間窗內同相軸的時間差除以同相軸跨越的地震道數來確定，而且tl可以根據不同位置同相軸傾角的大小做相應的空間變化，上式中的加減法也與同相軸的傾角有關系，如果隨選定時間空間窗地震道道號增大，選定的同相軸時間在減小，那么就采用上段列舉的式子計算，如果隨選定時間空間窗地震道道號增大，選定的同相軸時間在增大，則需把上段的加法變成減法，以適應同相軸傾角的變化。
[0070]S106具體實施時，根據S105確定的每個地震道的第一個樣點和選定的一段小時間窗樣點長度N*si，將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點，確定所述時間窗的位置；
[0071]S107具體實施時，據S106確定的小段時間窗，選擇中間地震道對應的地震數據，以中間地震道的樣點數據與剩余地震道的樣點數據做互相關，分別獲得剩余幾個地震道的最大相關值，將剩余地震道的最大相關值按比例相加，獲得所述相關比值之和。與中心地震道越近的地震道比例系數越大，相反就越小，但是所有比例系數的和為I。
[0072]需要注意的是，如果選擇5個地震道，那么第I個地震道的比例為I/(2* (1+2))，第2個地震道為2*1/(2*(1+2)),第4個地震道為2*1/(2*(1+2)),第5個地震道為I/(2*(1+2))，I/(2*(1+2))+2*1/(2*(1+2))+2*1/(2*(1+2))+1/(2*(1+2)) = I。
[0073]S108具體實施時，包括:重復S105至S107，生成S104選定的時間空間窗內中間地震道中除了第一個樣點外的所有樣點的最大互相關值之和；
[0074]將時間空間窗內中間地震道中包括第一個樣點的所有樣點的最大互相關值之和相加，形成時間空間窗內中間地震道最大互相關值的總和，道間時間為tl。
[0075]S109具體實施時，包括:對疊加后的震數據的所有地震道都重復執行S104至S108，生成疊后的每個地震道在選擇的時間空間窗內的最大相關值總和。
[0076]在一實施例中，如圖2所示，疊后信噪比估算方法還包括:
[0077]SllO:對每個地震道的所述最大相關值總和計算常用對數。這樣就把最大相關值轉換成了分貝數，便于信噪比曲線顯示。
[0078]在一實施例中，在所述S104中，選擇所述時間空間窗時，時間窗中包括一個完整的同相軸，空間窗的地震道的個數為奇數，所述空間窗按所述疊后地震數據地震道順序號依次滾動。
[0079]如圖3所示，本發明提供一種疊后信噪比估算裝置，所述的疊后信噪比估算裝置包括:數據采集單元301,數據預處理單元302,數據疊加單元303,時間窗選擇單元304,樣點計算單元305，位置確定單元306，相關比值生成單元307以及總和計算單元308。
[0080]數據采集單元301用于采集地震數據。
[0081]數據預處理單元302用于對地震數據進行預處理，形成預處理后的地震數據。
[0082]數據疊加單元303用于對預處理后的地震數據進行疊加，形成疊后地震數據。
[0083]時間窗選擇單元304用于在疊后地震數據上選擇一個時間空間窗，確定時間空間窗的中間地震道。
[0084]樣點計算單元305用于以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的樣點。
[0085]位置確定單元306用于將所述第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點，確定所述時間窗的位置。
[0086]相關比值生成單元307用于根據所述時間窗的位置從所述時間窗中選擇所述中間地震道的樣點數據，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和。
[0087]總和計算單元308用于生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和。
[0088]在一實施例中，數據預處理單元302具體用于:對所述地震數據進行地震數據解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0089]在一實施例中，如圖4所示，疊后信噪比估算裝置還包括:判斷單元401及動校正單元402。
[0090]判斷單元401用于判斷是否對所述預處理后的地震數據進行了疊前偏移處理，動校正單元402用于對所述預處理后的地震數據進行動校正。
[0091]在一實施例中，所述的相關比值生成單元307具體用于:以中間地震道的樣點數據與剩余地震道的樣點數據做互相關，分別獲得剩余地震道的最大相關值，將剩余地震道的最大相關值按比例相加，獲得相關比值之和。
[0092]在一實施例中，所述的總和計算單元308具體用于:生成時間空間窗內中間地震道中除了第一個樣點外的所有樣點的最大互相關值之和；將時間空間窗內中間地震道中包括第一個樣點的所有樣點的最大互相關值之和相加，形成時間空間窗內中間地震道最大互相關值的總和。
[0093]在一實施例中，如圖5所示，所述的疊后信噪比估算裝置還包括:對數計算單元501，用于對每個地震道的所述最大相關值總和計算常用對數。這樣就把最大相關值轉換成了分貝數，便于信噪比曲線顯示。
[0094]在一實施例中，選擇所述時間空間窗中的時間窗中包括一個完整的同相軸，空間窗的地震道的個數為奇數。
[0095]本發明通過野外采集地震數據，對地震數據進行預處理，預處理大概包括:地震數據解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析，疊前偏移等地震數據處理，形成預處理后的地震數據，然后把預處理后的地震數據進行疊加，形成疊后地震數據，如果沒有進行疊前偏移處理，疊加前需要做動校正和地震數據切除處理。在疊加后的地震數據上選擇一個時間空間窗，時間窗的選擇標準一般必須包括一個完整的同相軸，同相軸可以不連續，空間窗地震道的選擇一般為奇數3，5，7，9，11等，空間窗可以任意選擇，不大于實際疊加剖面的道數，但是如果選擇的同相軸信噪比高，地震道選擇越少，信噪比估算精度就越高，估算出的信噪比就越可靠。時間空間窗選擇后，確定中心地震道，根據同相軸傾角時間差確定每個地震道的相關時窗，然后用中心地震道分別與相鄰地震道做互相關處理，獲得最大互相關值。把這些值乘以不同的比例值加在一起，比例值的大小與離中心地震道的距離有關，距離越近值越大，相反比例值越小。同樣的方法處理剩余的每個地震道，每個地震道都有最大相關值的總和，和的大小反應了疊后地震信噪比的高低，為了方便信噪比曲線顯示，對最大相關值的總和取對數，獲得信噪比分貝數，用分貝數去度量地震數據信噪比的高低。
[0096]本發明一種適應構造傾角變化的疊后高效信噪比估算技術首先根據疊后地震數據，先選定需要計算靜校正的時間窗和空間窗，計算時間空間窗內的同相軸的傾角，確定相應相關時間窗長度。按計算最大相關值大小的和(按比例)，然后求最大相關值的分貝數，分貝數值的大小直接反應信噪比的高低。本發明與其它幾種時間域計算信噪比的方法相t匕，更能適應傾角變化的疊后信噪比估算，信噪比估算精度更可靠。
[0097]為了更好的說明本發明實施例，下面結合兩個例子進行簡單說明。
[0098]實施例1是本發明的一種適應構造傾角變化的疊后高效信噪比估算方法，運用于模型地震數據信噪比估算。圖6是模型數據和信噪比曲線顯示示意圖，圖6包括兩個部分，圖6的上面是模型數據，下面是信噪比曲線，模型數據和信噪比曲線的橫坐標都是CMP號，每個CMP號代表一個地震道，模型數據的縱坐標是時間，信噪比曲線圖的縱坐標是分貝數，顯示的兩條信噪比曲線一條是本方法計算出的信噪比曲線，曲線不光滑，一條是平均信噪比曲線，是在本方法計算出的信噪比基礎上進行了光滑處理，便于信噪比顯示和整體觀察。圖6設計的是一個不連續的水平同相軸，左邊四個同相軸，信噪比較高，右邊一個同相軸信噪比較低，中間有部分同相軸缺失信噪比就更低，這與信噪比曲線顯示的信噪比高低相吻合，證明了本發明信噪比估算的正確性。圖7是傾角模型數據和信噪比曲線顯示，與圖6相t匕，為了驗證本方法信噪比估算能適應構造傾角變化，圖7只設計了一個傾斜同相軸，整個傾角模型的信噪比存在很小幅度的變化，這種變化是隨機的，圖7中的信噪比曲線也證實了這一點，整個數據平均信噪比曲線的分貝值沒有超過5分貝，與理論分析相吻合，這證明了本發明能適應構造傾角變化的信噪比估算。
[0099]實施例2是本發明的一種適應構造傾角變化的疊后高效信噪比估算方法，運用于實際地震數據信噪比估算。圖8是實際地震數據及信噪比曲線顯示示意圖，圖8上方是模型數據，圖8下方是信噪比曲線，估算圖8信噪比的時間窗是圖8的整個地震數據，地震數據剖面上可以看出右邊信噪比稍高，中間最低，這與圖8下方的信噪比曲線完全吻合。圖9是一個低信噪比地震數據，有一組傾斜同相軸，從整個地震數據剖面看，數據左邊信噪比高，右邊沒有同相軸，信噪比應該比較低。圖9的信噪比曲線恰好證明了我們肉眼對地震數據信噪比的分析和判斷結果。
[0100]理論分析、模型數據試算、以及實際地震數據實驗表明:本發明提高了疊后存在傾角的地震數據信噪比估算精度，為研究和定性分析地震數據處理中信噪比的提高程度提供有利的質量控制工具。
[0101]本領域內的技術人員應明白，本發明的實施例可提供為方法、系統、或計算機程序產品。因此，本發明可采用完全硬件實施例、完全軟件實施例、或結合軟件和硬件方面的實施例的形式。而且，本發明可采用在一個或多個其中包含有計算機可用程序代碼的計算機可用存儲介質(包括但不限于磁盤存儲器、CD-ROM、光學存儲器等)上實施的計算機程序產品的形式。
[0102]本發明是參照根據本發明實施例的方法、設備(系統)、和計算機程序產品的流程圖和/或方框圖來描述的。應理解可由計算機程序指令實現流程圖和/或方框圖中的每一流程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結合。可提供這些計算機程序指令到通用計算機、專用計算機、嵌入式處理機或其他可編程數據處理設備的處理器以產生一個機器，使得通過計算機或其他可編程數據處理設備的處理器執行的指令產生用于實現在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的裝置。
[0103]這些計算機程序指令也可存儲在能引導計算機或其他可編程數據處理設備以特定方式工作的計算機可讀存儲器中，使得存儲在該計算機可讀存儲器中的指令產生包括指令裝置的制造品，該指令裝置實現在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能。
[0104]這些計算機程序指令也可裝載到計算機或其他可編程數據處理設備上，使得在計算機或其他可編程設備上執行一系列操作步驟以產生計算機實現的處理，從而在計算機或其他可編程設備上執行的指令提供用于實現在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的步驟。
[0105]本發明中應用了具體實施例對本發明的原理及實施方式進行了闡述，以上實施例的說明只是用于幫助理解本發明的方法及其核心思想；同時，對于本領域的一般技術人員，依據本發明的思想，在【具體實施方式】及應用范圍上均會有改變之處，綜上所述，本說明書內容不應理解為對本發明的限制。
【權利要求】
1.一種疊后信噪比估算方法，其特征在于，所述的疊后信噪比估算方法包括: 步驟1:采集地震數據； 步驟2:對所述地震數據進行預處理，形成預處理后的地震數據； 步驟3:對所述預處理后的地震數據進行疊加，形成疊后地震數據； 步驟4:在所述疊后地震數據上選擇一個時間空間窗，確定所述時間空間窗的中間地震道; 步驟5:以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的第一樣點； 步驟6:將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置； 步驟7:根據所述時間窗的位置從所述時間窗中選擇所述中間地震道的N個樣點長的數據，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和； 步驟8:重復執行步驟5至步驟7，生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和； 步驟9:重復執行步驟4至步驟8，生成所述疊后的每個地震道在選擇的時間空間窗內的最大相關值總和。
2.根據權利要求1所述的疊后信噪比估算方法，其特征在于，對所述地震數據進行預處理，包括:對所述地震數據進行地震數據解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
3.根據權利要求1所述的疊后信噪比估算方法，其特征在于，所述疊后信噪比估算方法還包括: 判斷是否對所述預處理后的地震數據進行了疊前偏移處理； 如果否，對所述預處理后的地震數據進行動校正。
4.根據權利要求1所述的疊后信噪比估算方法，其特征在于，在所述步驟7中，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和，包括:以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，分別獲得所述剩余地震道的最大相關值，將所述剩余地震道的最大相關值按比例相加，獲得所述相關比值之和。
5.根據權利要求1所述的疊后信噪比估算方法，其特征在于，在所述步驟8中，生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和，包括: 生成所述時間空間窗內中間地震道中除了第一個樣點外的所有樣點的最大互相關值之和； 將所述時間空間窗內中間地震道中包括第一個樣點的所有樣點的最大互相關值之和相加，形成所述時間空間窗內中間地震道最大互相關值總和。
6.根據權利要求4所述的疊后信噪比估算方法，其特征在于，所述疊后信噪比估算方法還包括:對每個地震道的所述最大相關值總和計算常用對數。
7.根據權利要求1所述的疊后信噪比估算方法，其特征在于，在所述步驟4中，選擇所述時間空間窗時，時間窗中包括一個完整的同相軸，空間窗的地震道的個數為奇數，所述空間窗按所述疊后地震數據地震道順序號依次滾動。
8.一種疊后信噪比估算裝置，其特征在于，所述的疊后信噪比估算裝置包括: 數據采集單元，用于采集地震數據； 數據預處理單元，用于對所述地震數據進行預處理，形成預處理后的地震數據； 數據疊加單元，用于對所述預處理后的地震數據進行疊加，形成疊后地震數據； 時間窗選擇單元，用于在所述疊后地震數據上選擇一個時間空間窗，確定所述時間空間窗的中間地震道； 樣點計算單元，用于以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的樣點； 位置確定單兀，用于將所述第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置； 相關比值生成單元，用于根據所述時間窗的位置從所述時間窗中選擇所述中間地震道的樣點數據，以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，生成所述第一樣點的相關比值之和； 總和計算單元，用于生成所述時間空間窗內中間地震道的所有樣點的最大相關值之和及所述中間地震道的最大相關值總和。
9.根據權利要求8所述的疊后信噪比估算裝置，其特征在于，所述數據預處理單元具體用于:對所述地震數據進行地震數據解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
10.根據權利要求8所述的疊后信噪比估算裝置，其特征在于，所述的疊后信噪比估算裝置還包括: 判斷單元，用于判斷是否對所述預處理后的地震數據進行了疊前偏移處理； 動校正單元，用于對所述預處理后的地震數據進行動校正。
11.根據權利要求8所述的疊后信噪比估算裝置，其特征在于，所述的相關比值生成單元具體用于:以所述中間地震道的樣點數據與剩余地震道的樣點數據做互相關，分別獲得所述剩余地震道的最大相關值，將所述剩余地震道的最大相關值按比例相加，獲得所述相關比值之和。
12.根據權利要求8所述的疊后信噪比估算裝置，其特征在于，所述的總和計算單元具體用于:生成所述時間空間窗內中間地震道中除了第一個樣點外的所有樣點的最大互相關值之和；將所述時間空間窗內中間地震道中包括第一個樣點的所有樣點的最大互相關值之和相加，形成所述時間空間窗內中間地震道最大互相關值的總和。
13.根據權利要求12所述的疊后信噪比估算裝置，其特征在于，所述的疊后信噪比估算裝置還包括:對數計算單元，用于對每個地震道的所述最大相關值總和計算常用對數。
14.根據權利要求13所述的疊后信噪比估算裝置，其特征在于，選擇所述時間空間窗中的時間窗中包括一個完整的同相軸，空間窗的地震道的個數為奇數。
【文檔編號】G01V1/36GK104330821SQ201410404335
【公開日】2015年2月4日 申請日期:2014年8月15日 優先權日:2014年8月15日
【發明者】熊定鈺, 錢忠平, 陳海云, 代麗華, 祝寬海 申請人:中國石油天然氣集團公司, 中國石油集團東方地球物理勘探有限責任公司