疊后信噪比估算方法及裝置制造方法
【專利摘要】一種疊后信噪比估算方法及裝置，該方法包括：采集地震數(shù)據(jù)并進行預(yù)處理，形成預(yù)處理后的地震數(shù)據(jù)；對預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)；選擇一個時間空間窗，確定時間空間窗的中間地震道；以中間地震道為中心確定時間空間窗中第一地震道的第一樣點，選定N個樣點長的時間窗，并計算時間空間窗中剩余每一地震道的第一樣點；確定時間窗的位置，根據(jù)時間窗的位置選擇中間地震道的N個樣點長的數(shù)據(jù)，根據(jù)中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)生成第一樣點的相關(guān)比值之和；計算時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及中間地震道的最大相關(guān)值總和；計算所述疊后的每個地震道在選擇的時間空間窗內(nèi)的最大相關(guān)值總和。
【專利說明】疊后信噪比估算方法及裝置

【技術(shù)領(lǐng)域】
[0001]本發(fā)明本方法涉及地球物理勘探數(shù)據(jù)處理技術(shù)，特別是關(guān)于是一種疊后信噪比估算方法及裝置。

【背景技術(shù)】
[0002]在地震數(shù)據(jù)處理技術(shù)中，噪音常常影響疊加成像效果，噪聲衰減在很多地震數(shù)據(jù)處理中都變成了常規(guī)處理流程，去噪后信噪比的提升一直是一大難題。如何對野外數(shù)據(jù)采集和室內(nèi)數(shù)據(jù)處理進行地震數(shù)據(jù)信噪比宏觀定性分析和數(shù)字顯示，準(zhǔn)確地分析資料處理中的各個環(huán)節(jié)質(zhì)量，一直是地震勘探研究人員想要實現(xiàn)的目標(biāo)。
[0003]目前信噪比主要是在指定時間內(nèi)，用信號與噪聲振幅的比值來衡量信號的強弱的方法。現(xiàn)有技術(shù)中，主要信噪比估算方法有:1)時間域疊加法；2)時間域奇異值分減法；3)模型道相關(guān)估算法。這三種信噪比估算方法，在地下構(gòu)造傾角變化較大時，信噪比估算精度不高。時間域疊加和模型道相關(guān)主要是因為地下構(gòu)造傾角變化時同相軸不能準(zhǔn)確地同相疊力口。而時間域奇異值分減法本身也不適應(yīng)地下構(gòu)造傾角較大變化的地震數(shù)據(jù)噪聲估計。


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

[0004]本發(fā)明提供了一種疊后信噪比估算方法及裝置，以精確地震估算地震數(shù)據(jù)的信噪比。
[0005]為了實現(xiàn)上述目的，本發(fā)明提供一種疊后信噪比估算方法，所述的疊后信噪比估算方法包括:
[0006]步驟1:采集地震數(shù)據(jù)；
[0007]步驟2:對所述地震數(shù)據(jù)進行預(yù)處理，形成預(yù)處理后的地震數(shù)據(jù)；
[0008]步驟3:對所述預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)；
[0009]步驟4:在所述疊后地震數(shù)據(jù)上選擇一個時間空間窗，確定所述時間空間窗的中間地震道；
[0010]步驟5:以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的第一樣點；
[0011]步驟6:將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置；
[0012]步驟7:根據(jù)所述時間窗的位置從所述時間窗中選擇所述中間地震道的N個樣點長的數(shù)據(jù)，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和；
[0013]步驟8:重復(fù)執(zhí)行步驟5至步驟7，生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和；
[0014]步驟9:重復(fù)執(zhí)行步驟4至步驟8，生成所述疊后的每個地震道在選擇的時間空間窗內(nèi)的最大相關(guān)值總和。
[0015]在一實施例中，對所述地震數(shù)據(jù)進行預(yù)處理，包括:對所述地震數(shù)據(jù)進行地震數(shù)據(jù)解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0016]在一實施例中，所述疊后信噪比估算方法還包括:
[0017]判斷是否對所述預(yù)處理后的地震數(shù)據(jù)進行了疊前偏移處理；
[0018]如果否，對所述預(yù)處理后的地震數(shù)據(jù)進行動校正。
[0019]在一實施例中，在所述步驟7中，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和，包括:以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，分別獲得所述剩余地震道的最大相關(guān)值，將所述剩余地震道的最大相關(guān)值按比例相加，獲得所述相關(guān)比值之和。
[0020]在一實施例中，在所述步驟8中，生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和，包括:
[0021]生成所述時間空間窗內(nèi)中間地震道中除了第一個樣點外的所有樣點的最大互相關(guān)值之和；
[0022]將所述時間空間窗內(nèi)中間地震道中包括第一個樣點的所有樣點的最大互相關(guān)值之和相加，形成所述時間空間窗內(nèi)中間地震道最大互相關(guān)值的總和。
[0023]在一實施例中，所述疊后信噪比估算方法還包括:對每個地震道的所述最大相關(guān)值總和計算常用對數(shù)。
[0024]在一實施例中，在所述步驟4中，選擇所述時間空間窗時，時間窗中包括一個完整的同相軸，空間窗的地震道的個數(shù)為奇數(shù)，所述空間窗按所述疊后地震數(shù)據(jù)地震道順序號依次滾動。
[0025]為了實現(xiàn)上述目的，本發(fā)明提供一種疊后信噪比估算裝置，所述的疊后信噪比估算裝置包括:
[0026]數(shù)據(jù)采集單元，用于采集地震數(shù)據(jù)；
[0027]數(shù)據(jù)預(yù)處理單元，用于對所述地震數(shù)據(jù)進行預(yù)處理，形成預(yù)處理后的地震數(shù)據(jù)；
[0028]數(shù)據(jù)疊加單元，用于對所述預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)；
[0029]時間窗選擇單元，用于在所述疊后地震數(shù)據(jù)上選擇一個時間空間窗，確定所述時間空間窗的中間地震道；
[0030]樣點計算單元，用于以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的樣點；
[0031]位置確定單元，用于將所述第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置；
[0032]相關(guān)比值生成單元，用于根據(jù)所述時間窗的位置從所述時間窗中選擇所述中間地震道的樣點數(shù)據(jù)，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和；
[0033]總和計算單元，用于生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和。
[0034]在一實施例中，所述數(shù)據(jù)預(yù)處理單元具體用于:對所述地震數(shù)據(jù)進行地震數(shù)據(jù)解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0035]在一實施例中，所述的疊后信噪比估算裝置還包括:
[0036]判斷單元，用于判斷是否對所述預(yù)處理后的地震數(shù)據(jù)進行了疊前偏移處理；
[0037]動校正單元，用于對所述預(yù)處理后的地震數(shù)據(jù)進行動校正。
[0038]在一實施例中，所述的相關(guān)比值生成單元具體用于:以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，分別獲得所述剩余地震道的最大相關(guān)值，將所述剩余地震道的最大相關(guān)值按比例相加，獲得所述相關(guān)比值之和。
[0039]在一實施例中，所述的總和計算單元具體用于:生成所述時間空間窗內(nèi)中間地震道中除了第一個樣點外的所有樣點的最大互相關(guān)值之和；將所述時間空間窗內(nèi)中間地震道中包括第一個樣點的所有樣點的最大互相關(guān)值之和相加，形成所述時間空間窗內(nèi)中間地震道最大互相關(guān)值的總和。
[0040]在一實施例中，所述的疊后信噪比估算裝置還包括:對數(shù)計算單元，用于對每個地震道的所述最大相關(guān)值總和計算常用對數(shù)。
[0041 ] 在一實施例中，選擇所述時間空間窗中的時間窗中包括一個完整的同相軸，空間窗的地震道的個數(shù)為奇數(shù)。
[0042]本發(fā)明提供的疊后信噪比估算方法，為適應(yīng)不同地下構(gòu)造地震數(shù)據(jù)信噪比估算，更精確地震估算地震數(shù)據(jù)的信噪比，為定性分析地震數(shù)據(jù)信噪比的提高提供了有利的工具。

【專利附圖】

【附圖說明】
[0043]為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。
[0044]圖1為本發(fā)明實施例的疊后信噪比估算方法流程圖一；
[0045]圖2為本發(fā)明實施例的疊后信噪比估算方法流程圖二 ；
[0046]圖3為本發(fā)明實施例的疊后信噪比估算裝置結(jié)構(gòu)框圖一；
[0047]圖4為本發(fā)明實施例的疊后信噪比估算裝置結(jié)構(gòu)框圖二 ；
[0048]圖5為本發(fā)明實施例的疊后信噪比估算裝置結(jié)構(gòu)框圖三；
[0049]圖6為模型數(shù)據(jù)和信噪比曲線顯示圖(上:模型數(shù)據(jù)，下:信噪比曲線)；
[0050]圖7為傾角模型數(shù)據(jù)和信噪比曲線顯示圖(上:模型數(shù)據(jù)，下:信噪比曲線)；
[0051]圖8為實際地震數(shù)據(jù)及信噪比曲線顯示圖(上:模型數(shù)據(jù)，下:信噪比曲線)；
[0052]圖9為低信噪比實際地震數(shù)據(jù)及信噪比曲線顯示圖(上:模型數(shù)據(jù)，下:信噪比曲線)。

【具體實施方式】
[0053]下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例?；诒景l(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護的范圍。
[0054]如圖1所示，本發(fā)明提供一種疊后信噪比估算方法，所述的疊后信噪比估算方法包括:
[0055]SlOl:采集地震數(shù)據(jù)；
[0056]S102:對所述地震數(shù)據(jù)進行預(yù)處理，形成預(yù)處理后的地震數(shù)據(jù)；
[0057]S103:對所述預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)；
[0058]S104:在所述疊后地震數(shù)據(jù)上選擇一個時間空間窗，確定所述時間空間窗的中間地震道；
[0059]S105:以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的第一樣點；
[0060]S106:將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置；
[0061]S107:根據(jù)所述時間窗的位置從所述時間窗中選擇所述中間地震道的N個樣點長的數(shù)據(jù)，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和；
[0062]S108:重復(fù)執(zhí)行S105至S107，生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和；
[0063]S109:重復(fù)執(zhí)行S104至S108，生成所述疊后的每個地震道在選擇的時間空間窗內(nèi)的最大相關(guān)值總和。
[0064]由圖1的流程可知，本發(fā)明的疊后信噪比估算方法，首先根據(jù)疊后地震數(shù)據(jù)，先選定需要計算靜校正的時間窗和空間窗，計算時間空間窗內(nèi)的同相軸的傾角，確定相應(yīng)相關(guān)時間窗長度。按計算最大相關(guān)值大小的和(按比例)求最大相關(guān)值的分貝數(shù)，分貝數(shù)值的大小直接反應(yīng)信噪比的高低。本發(fā)明與現(xiàn)有技術(shù)中的幾種時間域計算信噪比的方法相比，更能適應(yīng)傾角變化的疊后信噪比估算，信噪比估算精度更可靠。
[0065]S102中，對地震數(shù)據(jù)進行預(yù)處理可以包括:對所述地震數(shù)據(jù)進行地震數(shù)據(jù)解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0066]在一實施例中，對所述預(yù)處理后的地震數(shù)據(jù)進行疊加前，還可以判斷是否對預(yù)處理后的地震數(shù)據(jù)進行了疊前偏移處理；如果沒有對預(yù)處理后的地震數(shù)據(jù)進行疊前偏移處理，對預(yù)處理后的地震數(shù)據(jù)進行動校正。
[0067]在S104中，空間窗地震道的選擇一般為奇數(shù)，如3，5，7，9，11等。對于確定中間地震道，如果選擇11個地震道的空間窗，那么時間空間窗的中間地震道就是第6個地震道。
[0068]S105具體實施時，以S104選擇出的中間地震道為中心，確定時間空間窗中第一地震道的第一樣點，以第一樣點為中心選定N個樣點長的一段小的時間窗，N為奇數(shù)。給定一個道間時間差tl，根據(jù)采樣率Si可以計算出11 = tl/si。如果選定11個地震道空間窗，假設(shè)第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與同相軸傾角有關(guān)系，一般都用選定時間空間窗內(nèi)同相軸的時間差除以同相軸跨越的地震道數(shù)來確定，而且tl可以根據(jù)不同位置同相軸傾角的大小做相應(yīng)的空間變化，上式中的加減法也與同相軸的傾角有關(guān)系，如果隨選定時間空間窗地震道道號增大，選定的同相軸時間在減小，那么就采用上段列舉的式子計算，如果隨選定時間空間窗地震道道號增大，選定的同相軸時間在增大，則需把上段的加法變成減法，以適應(yīng)同相軸傾角的變化。
[0070]S106具體實施時，根據(jù)S105確定的每個地震道的第一個樣點和選定的一段小時間窗樣點長度N*si，將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點，確定所述時間窗的位置；
[0071]S107具體實施時，據(jù)S106確定的小段時間窗，選擇中間地震道對應(yīng)的地震數(shù)據(jù)，以中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，分別獲得剩余幾個地震道的最大相關(guān)值，將剩余地震道的最大相關(guān)值按比例相加，獲得所述相關(guān)比值之和。與中心地震道越近的地震道比例系數(shù)越大，相反就越小，但是所有比例系數(shù)的和為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具體實施時，包括:重復(fù)S105至S107，生成S104選定的時間空間窗內(nèi)中間地震道中除了第一個樣點外的所有樣點的最大互相關(guān)值之和；
[0074]將時間空間窗內(nèi)中間地震道中包括第一個樣點的所有樣點的最大互相關(guān)值之和相加，形成時間空間窗內(nèi)中間地震道最大互相關(guān)值的總和，道間時間為tl。
[0075]S109具體實施時，包括:對疊加后的震數(shù)據(jù)的所有地震道都重復(fù)執(zhí)行S104至S108，生成疊后的每個地震道在選擇的時間空間窗內(nèi)的最大相關(guān)值總和。
[0076]在一實施例中，如圖2所示，疊后信噪比估算方法還包括:
[0077]SllO:對每個地震道的所述最大相關(guān)值總和計算常用對數(shù)。這樣就把最大相關(guān)值轉(zhuǎn)換成了分貝數(shù)，便于信噪比曲線顯示。
[0078]在一實施例中，在所述S104中，選擇所述時間空間窗時，時間窗中包括一個完整的同相軸，空間窗的地震道的個數(shù)為奇數(shù)，所述空間窗按所述疊后地震數(shù)據(jù)地震道順序號依次滾動。
[0079]如圖3所示，本發(fā)明提供一種疊后信噪比估算裝置，所述的疊后信噪比估算裝置包括:數(shù)據(jù)采集單元301,數(shù)據(jù)預(yù)處理單元302,數(shù)據(jù)疊加單元303,時間窗選擇單元304,樣點計算單元305，位置確定單元306，相關(guān)比值生成單元307以及總和計算單元308。
[0080]數(shù)據(jù)采集單元301用于采集地震數(shù)據(jù)。
[0081]數(shù)據(jù)預(yù)處理單元302用于對地震數(shù)據(jù)進行預(yù)處理，形成預(yù)處理后的地震數(shù)據(jù)。
[0082]數(shù)據(jù)疊加單元303用于對預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)。
[0083]時間窗選擇單元304用于在疊后地震數(shù)據(jù)上選擇一個時間空間窗，確定時間空間窗的中間地震道。
[0084]樣點計算單元305用于以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的樣點。
[0085]位置確定單元306用于將所述第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點，確定所述時間窗的位置。
[0086]相關(guān)比值生成單元307用于根據(jù)所述時間窗的位置從所述時間窗中選擇所述中間地震道的樣點數(shù)據(jù)，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和。
[0087]總和計算單元308用于生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和。
[0088]在一實施例中，數(shù)據(jù)預(yù)處理單元302具體用于:對所述地震數(shù)據(jù)進行地震數(shù)據(jù)解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
[0089]在一實施例中，如圖4所示，疊后信噪比估算裝置還包括:判斷單元401及動校正單元402。
[0090]判斷單元401用于判斷是否對所述預(yù)處理后的地震數(shù)據(jù)進行了疊前偏移處理，動校正單元402用于對所述預(yù)處理后的地震數(shù)據(jù)進行動校正。
[0091]在一實施例中，所述的相關(guān)比值生成單元307具體用于:以中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，分別獲得剩余地震道的最大相關(guān)值，將剩余地震道的最大相關(guān)值按比例相加，獲得相關(guān)比值之和。
[0092]在一實施例中，所述的總和計算單元308具體用于:生成時間空間窗內(nèi)中間地震道中除了第一個樣點外的所有樣點的最大互相關(guān)值之和；將時間空間窗內(nèi)中間地震道中包括第一個樣點的所有樣點的最大互相關(guān)值之和相加，形成時間空間窗內(nèi)中間地震道最大互相關(guān)值的總和。
[0093]在一實施例中，如圖5所示，所述的疊后信噪比估算裝置還包括:對數(shù)計算單元501，用于對每個地震道的所述最大相關(guān)值總和計算常用對數(shù)。這樣就把最大相關(guān)值轉(zhuǎn)換成了分貝數(shù)，便于信噪比曲線顯示。
[0094]在一實施例中，選擇所述時間空間窗中的時間窗中包括一個完整的同相軸，空間窗的地震道的個數(shù)為奇數(shù)。
[0095]本發(fā)明通過野外采集地震數(shù)據(jù)，對地震數(shù)據(jù)進行預(yù)處理，預(yù)處理大概包括:地震數(shù)據(jù)解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析，疊前偏移等地震數(shù)據(jù)處理，形成預(yù)處理后的地震數(shù)據(jù)，然后把預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)，如果沒有進行疊前偏移處理，疊加前需要做動校正和地震數(shù)據(jù)切除處理。在疊加后的地震數(shù)據(jù)上選擇一個時間空間窗，時間窗的選擇標(biāo)準(zhǔn)一般必須包括一個完整的同相軸，同相軸可以不連續(xù)，空間窗地震道的選擇一般為奇數(shù)3，5，7，9，11等，空間窗可以任意選擇，不大于實際疊加剖面的道數(shù)，但是如果選擇的同相軸信噪比高，地震道選擇越少，信噪比估算精度就越高，估算出的信噪比就越可靠。時間空間窗選擇后，確定中心地震道，根據(jù)同相軸傾角時間差確定每個地震道的相關(guān)時窗，然后用中心地震道分別與相鄰地震道做互相關(guān)處理，獲得最大互相關(guān)值。把這些值乘以不同的比例值加在一起，比例值的大小與離中心地震道的距離有關(guān)，距離越近值越大，相反比例值越小。同樣的方法處理剩余的每個地震道，每個地震道都有最大相關(guān)值的總和，和的大小反應(yīng)了疊后地震信噪比的高低，為了方便信噪比曲線顯示，對最大相關(guān)值的總和取對數(shù)，獲得信噪比分貝數(shù)，用分貝數(shù)去度量地震數(shù)據(jù)信噪比的高低。
[0096]本發(fā)明一種適應(yīng)構(gòu)造傾角變化的疊后高效信噪比估算技術(shù)首先根據(jù)疊后地震數(shù)據(jù)，先選定需要計算靜校正的時間窗和空間窗，計算時間空間窗內(nèi)的同相軸的傾角，確定相應(yīng)相關(guān)時間窗長度。按計算最大相關(guān)值大小的和(按比例)，然后求最大相關(guān)值的分貝數(shù)，分貝數(shù)值的大小直接反應(yīng)信噪比的高低。本發(fā)明與其它幾種時間域計算信噪比的方法相t匕，更能適應(yīng)傾角變化的疊后信噪比估算，信噪比估算精度更可靠。
[0097]為了更好的說明本發(fā)明實施例，下面結(jié)合兩個例子進行簡單說明。
[0098]實施例1是本發(fā)明的一種適應(yīng)構(gòu)造傾角變化的疊后高效信噪比估算方法，運用于模型地震數(shù)據(jù)信噪比估算。圖6是模型數(shù)據(jù)和信噪比曲線顯示示意圖，圖6包括兩個部分，圖6的上面是模型數(shù)據(jù)，下面是信噪比曲線，模型數(shù)據(jù)和信噪比曲線的橫坐標(biāo)都是CMP號，每個CMP號代表一個地震道，模型數(shù)據(jù)的縱坐標(biāo)是時間，信噪比曲線圖的縱坐標(biāo)是分貝數(shù)，顯示的兩條信噪比曲線一條是本方法計算出的信噪比曲線，曲線不光滑，一條是平均信噪比曲線，是在本方法計算出的信噪比基礎(chǔ)上進行了光滑處理，便于信噪比顯示和整體觀察。圖6設(shè)計的是一個不連續(xù)的水平同相軸，左邊四個同相軸，信噪比較高，右邊一個同相軸信噪比較低，中間有部分同相軸缺失信噪比就更低，這與信噪比曲線顯示的信噪比高低相吻合，證明了本發(fā)明信噪比估算的正確性。圖7是傾角模型數(shù)據(jù)和信噪比曲線顯示，與圖6相t匕，為了驗證本方法信噪比估算能適應(yīng)構(gòu)造傾角變化，圖7只設(shè)計了一個傾斜同相軸，整個傾角模型的信噪比存在很小幅度的變化，這種變化是隨機的，圖7中的信噪比曲線也證實了這一點，整個數(shù)據(jù)平均信噪比曲線的分貝值沒有超過5分貝，與理論分析相吻合，這證明了本發(fā)明能適應(yīng)構(gòu)造傾角變化的信噪比估算。
[0099]實施例2是本發(fā)明的一種適應(yīng)構(gòu)造傾角變化的疊后高效信噪比估算方法，運用于實際地震數(shù)據(jù)信噪比估算。圖8是實際地震數(shù)據(jù)及信噪比曲線顯示示意圖，圖8上方是模型數(shù)據(jù)，圖8下方是信噪比曲線，估算圖8信噪比的時間窗是圖8的整個地震數(shù)據(jù)，地震數(shù)據(jù)剖面上可以看出右邊信噪比稍高，中間最低，這與圖8下方的信噪比曲線完全吻合。圖9是一個低信噪比地震數(shù)據(jù)，有一組傾斜同相軸，從整個地震數(shù)據(jù)剖面看，數(shù)據(jù)左邊信噪比高，右邊沒有同相軸，信噪比應(yīng)該比較低。圖9的信噪比曲線恰好證明了我們?nèi)庋蹖Φ卣饠?shù)據(jù)信噪比的分析和判斷結(jié)果。
[0100]理論分析、模型數(shù)據(jù)試算、以及實際地震數(shù)據(jù)實驗表明:本發(fā)明提高了疊后存在傾角的地震數(shù)據(jù)信噪比估算精度，為研究和定性分析地震數(shù)據(jù)處理中信噪比的提高程度提供有利的質(zhì)量控制工具。
[0101]本領(lǐng)域內(nèi)的技術(shù)人員應(yīng)明白，本發(fā)明的實施例可提供為方法、系統(tǒng)、或計算機程序產(chǎn)品。因此，本發(fā)明可采用完全硬件實施例、完全軟件實施例、或結(jié)合軟件和硬件方面的實施例的形式。而且，本發(fā)明可采用在一個或多個其中包含有計算機可用程序代碼的計算機可用存儲介質(zhì)(包括但不限于磁盤存儲器、CD-ROM、光學(xué)存儲器等)上實施的計算機程序產(chǎn)品的形式。
[0102]本發(fā)明是參照根據(jù)本發(fā)明實施例的方法、設(shè)備(系統(tǒng))、和計算機程序產(chǎn)品的流程圖和/或方框圖來描述的。應(yīng)理解可由計算機程序指令實現(xiàn)流程圖和/或方框圖中的每一流程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結(jié)合?？商峁┻@些計算機程序指令到通用計算機、專用計算機、嵌入式處理機或其他可編程數(shù)據(jù)處理設(shè)備的處理器以產(chǎn)生一個機器，使得通過計算機或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的裝置。
[0103]這些計算機程序指令也可存儲在能引導(dǎo)計算機或其他可編程數(shù)據(jù)處理設(shè)備以特定方式工作的計算機可讀存儲器中，使得存儲在該計算機可讀存儲器中的指令產(chǎn)生包括指令裝置的制造品，該指令裝置實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能。
[0104]這些計算機程序指令也可裝載到計算機或其他可編程數(shù)據(jù)處理設(shè)備上，使得在計算機或其他可編程設(shè)備上執(zhí)行一系列操作步驟以產(chǎn)生計算機實現(xiàn)的處理，從而在計算機或其他可編程設(shè)備上執(zhí)行的指令提供用于實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的步驟。
[0105]本發(fā)明中應(yīng)用了具體實施例對本發(fā)明的原理及實施方式進行了闡述，以上實施例的說明只是用于幫助理解本發(fā)明的方法及其核心思想；同時，對于本領(lǐng)域的一般技術(shù)人員，依據(jù)本發(fā)明的思想，在【具體實施方式】及應(yīng)用范圍上均會有改變之處，綜上所述，本說明書內(nèi)容不應(yīng)理解為對本發(fā)明的限制。
【權(quán)利要求】
1.一種疊后信噪比估算方法，其特征在于，所述的疊后信噪比估算方法包括: 步驟1:采集地震數(shù)據(jù)； 步驟2:對所述地震數(shù)據(jù)進行預(yù)處理，形成預(yù)處理后的地震數(shù)據(jù)； 步驟3:對所述預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)； 步驟4:在所述疊后地震數(shù)據(jù)上選擇一個時間空間窗，確定所述時間空間窗的中間地震道; 步驟5:以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的第一樣點； 步驟6:將所述每一地震道第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置； 步驟7:根據(jù)所述時間窗的位置從所述時間窗中選擇所述中間地震道的N個樣點長的數(shù)據(jù)，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和； 步驟8:重復(fù)執(zhí)行步驟5至步驟7，生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和； 步驟9:重復(fù)執(zhí)行步驟4至步驟8，生成所述疊后的每個地震道在選擇的時間空間窗內(nèi)的最大相關(guān)值總和。
2.根據(jù)權(quán)利要求1所述的疊后信噪比估算方法，其特征在于，對所述地震數(shù)據(jù)進行預(yù)處理，包括:對所述地震數(shù)據(jù)進行地震數(shù)據(jù)解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
3.根據(jù)權(quán)利要求1所述的疊后信噪比估算方法，其特征在于，所述疊后信噪比估算方法還包括: 判斷是否對所述預(yù)處理后的地震數(shù)據(jù)進行了疊前偏移處理； 如果否，對所述預(yù)處理后的地震數(shù)據(jù)進行動校正。
4.根據(jù)權(quán)利要求1所述的疊后信噪比估算方法，其特征在于，在所述步驟7中，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和，包括:以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，分別獲得所述剩余地震道的最大相關(guān)值，將所述剩余地震道的最大相關(guān)值按比例相加，獲得所述相關(guān)比值之和。
5.根據(jù)權(quán)利要求1所述的疊后信噪比估算方法，其特征在于，在所述步驟8中，生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和，包括: 生成所述時間空間窗內(nèi)中間地震道中除了第一個樣點外的所有樣點的最大互相關(guān)值之和； 將所述時間空間窗內(nèi)中間地震道中包括第一個樣點的所有樣點的最大互相關(guān)值之和相加，形成所述時間空間窗內(nèi)中間地震道最大互相關(guān)值總和。
6.根據(jù)權(quán)利要求4所述的疊后信噪比估算方法，其特征在于，所述疊后信噪比估算方法還包括:對每個地震道的所述最大相關(guān)值總和計算常用對數(shù)。
7.根據(jù)權(quán)利要求1所述的疊后信噪比估算方法，其特征在于，在所述步驟4中，選擇所述時間空間窗時，時間窗中包括一個完整的同相軸，空間窗的地震道的個數(shù)為奇數(shù)，所述空間窗按所述疊后地震數(shù)據(jù)地震道順序號依次滾動。
8.一種疊后信噪比估算裝置，其特征在于，所述的疊后信噪比估算裝置包括: 數(shù)據(jù)采集單元，用于采集地震數(shù)據(jù)； 數(shù)據(jù)預(yù)處理單元，用于對所述地震數(shù)據(jù)進行預(yù)處理，形成預(yù)處理后的地震數(shù)據(jù)； 數(shù)據(jù)疊加單元，用于對所述預(yù)處理后的地震數(shù)據(jù)進行疊加，形成疊后地震數(shù)據(jù)； 時間窗選擇單元，用于在所述疊后地震數(shù)據(jù)上選擇一個時間空間窗，確定所述時間空間窗的中間地震道； 樣點計算單元，用于以所述中間地震道為中心，確定所述時間空間窗中第一地震道的第一樣點，以所述第一樣點為中心選定N個樣點長的時間窗，并計算所述時間空間窗中剩余每一地震道的樣點； 位置確定單兀，用于將所述第一個樣點作為所述時間窗的中心樣點向上向下分別移動N/2個樣點,確定所述時間窗的位置； 相關(guān)比值生成單元，用于根據(jù)所述時間窗的位置從所述時間窗中選擇所述中間地震道的樣點數(shù)據(jù)，以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，生成所述第一樣點的相關(guān)比值之和； 總和計算單元，用于生成所述時間空間窗內(nèi)中間地震道的所有樣點的最大相關(guān)值之和及所述中間地震道的最大相關(guān)值總和。
9.根據(jù)權(quán)利要求8所述的疊后信噪比估算裝置，其特征在于，所述數(shù)據(jù)預(yù)處理單元具體用于:對所述地震數(shù)據(jù)進行地震數(shù)據(jù)解編、靜校正、噪聲衰減、振幅補償、反褶積、速度分析及疊前偏移處理。
10.根據(jù)權(quán)利要求8所述的疊后信噪比估算裝置，其特征在于，所述的疊后信噪比估算裝置還包括: 判斷單元，用于判斷是否對所述預(yù)處理后的地震數(shù)據(jù)進行了疊前偏移處理； 動校正單元，用于對所述預(yù)處理后的地震數(shù)據(jù)進行動校正。
11.根據(jù)權(quán)利要求8所述的疊后信噪比估算裝置，其特征在于，所述的相關(guān)比值生成單元具體用于:以所述中間地震道的樣點數(shù)據(jù)與剩余地震道的樣點數(shù)據(jù)做互相關(guān)，分別獲得所述剩余地震道的最大相關(guān)值，將所述剩余地震道的最大相關(guān)值按比例相加，獲得所述相關(guān)比值之和。
12.根據(jù)權(quán)利要求8所述的疊后信噪比估算裝置，其特征在于，所述的總和計算單元具體用于:生成所述時間空間窗內(nèi)中間地震道中除了第一個樣點外的所有樣點的最大互相關(guān)值之和；將所述時間空間窗內(nèi)中間地震道中包括第一個樣點的所有樣點的最大互相關(guān)值之和相加，形成所述時間空間窗內(nèi)中間地震道最大互相關(guān)值的總和。
13.根據(jù)權(quán)利要求12所述的疊后信噪比估算裝置，其特征在于，所述的疊后信噪比估算裝置還包括:對數(shù)計算單元，用于對每個地震道的所述最大相關(guān)值總和計算常用對數(shù)。
14.根據(jù)權(quán)利要求13所述的疊后信噪比估算裝置，其特征在于，選擇所述時間空間窗中的時間窗中包括一個完整的同相軸，空間窗的地震道的個數(shù)為奇數(shù)。
【文檔編號】G01V1/36GK104330821SQ201410404335
【公開日】2015年2月4日 申請日期:2014年8月15日 優(yōu)先權(quán)日:2014年8月15日
【發(fā)明者】熊定鈺, 錢忠平, 陳海云, 代麗華, 祝寬海 申請人:中國石油天然氣集團公司, 中國石油集團東方地球物理勘探有限責(zé)任公司