一種繞射點成像的方法和設(shè)備的制作方法
【專利摘要】本發(fā)明公開了一種繞射點成像的方法和設(shè)備。該方法包括：獲取當前炮集數(shù)據(jù)，并計算地震波傳播走時；選取一個未選用過的成像點作為當前繞射成像點；依據(jù)地震波傳播走時，針對當前繞射成像點的地下深度，對在各個炮檢距處的當前炮集數(shù)據(jù)進行時差校正、高通濾波及疊加，得到當前炮集數(shù)據(jù)在當前繞射成像點處的繞射點成像數(shù)據(jù)；返回執(zhí)行所述選取一個未選用過的成像點作為當前繞射成像點，直至得到當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)；對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。通過本發(fā)明實施例的技術(shù)方案，可以使繞射點處繞射波信號突出，從而利用突出繞射點信號的數(shù)據(jù)成像之后可以更容易、更準確地識別出繞射點的位置。
【專利說明】一種繞射點成像的方法和設(shè)備
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地震勘探中的信號處理技術(shù)，特別是涉及一種繞射點成像的方法和設(shè)備。
【背景技術(shù)】
[0002]在物理學(xué)中，把波繞過障礙物而傳播的現(xiàn)象稱為繞射，或衍射。在地震勘探工作中，地震波在傳播過程中，遇到一些規(guī)模較小、物性橫向變化較快、與突變相關(guān)的特殊地質(zhì)體等地層劇烈變化的地方時，會產(chǎn)生繞射現(xiàn)象而產(chǎn)生繞射波。這些地層劇烈變化的地方稱之為繞射點，例如斷層的斷點、地層尖滅點、溶洞、裂縫等都屬于能夠使地震波產(chǎn)生繞射波的繞射點。繞射點可視為一個新的震源，繞射波可視為繞射點位置的振動而產(chǎn)生的。由于繞射點所在的位置通常位于地下小尺度的地質(zhì)體，而油氣藏往往也與這些地質(zhì)體相關(guān)聯(lián)，因此，尋找地下繞射點的位置在地震勘探工作中是十分重要的。
[0003]在地震勘探工作中，通常地面某一激發(fā)點(也稱炮點)振動而在地下激發(fā)出地震波，再通過地面其他點接收到地震信號(即該炮點的炮集數(shù)據(jù))來分析地下地質(zhì)體的。其中，接收到的地震信號不僅包括繞射波，還包括有反射波等其他形式的地震波。由于地震信號中反射波的信號較強，目前最常用于分析地震信號的方法是反射波成像法。但是，正是由于繞射波的信號較弱且難以成像的很好，所以在對接收到的地震信號進行反射波成像的成像結(jié)果中，繞射點難以被發(fā)現(xiàn)。
[0004]因此，為了識別繞射點位置，現(xiàn)有技術(shù)中主要采用的是先對接收到的地震信號區(qū)分反射波與繞射波，再對繞射波進行成像處理，最終實現(xiàn)從繞射波成像數(shù)據(jù)中分析繞射點位置。但是，由于繞射波與反射波在繞射點處相切，沿繞射波軌跡的一段空間內(nèi)繞射波會被反射波干擾，這使得現(xiàn)有技術(shù)一直無法準確分離繞射波與反射波，分離出的繞射波信號一直不夠準確，從而導(dǎo)致無法準確識別繞射點位置。
【發(fā)明內(nèi)容】
[0005]本發(fā)明所要解決的技術(shù)問題是，提供一種繞射點成像的方法和設(shè)備，通過先對炮集數(shù)據(jù)進行信號處理而得到突出繞射點位置的信號數(shù)據(jù)，再對處理后的信號數(shù)據(jù)進行成像處理，從而無需再對地震信號分離反射波與繞射波，以避免按照現(xiàn)有技術(shù)中反射波與繞射波分離不準確而導(dǎo)致的無法準確識別繞射點的位置。
[0006]為解決上述技術(shù)問題，本發(fā)明提供了一種繞射點成像的方法和設(shè)備，該方法包括:
[0007]獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，并計算地震波傳播走時；
[0008]選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點；
[0009]依據(jù)所述地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距處的所述當前炮集數(shù)據(jù)進行時差校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)；[0010]返回執(zhí)行所述選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點，直至得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)；
[0011]對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
[0012]可選的，所述依據(jù)各個炮檢距的地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距的所述當前炮集數(shù)據(jù)進行動校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)，包括:
[0013]依據(jù)所述地震波傳播走時，將在各個炮檢距處的所述當前炮集數(shù)據(jù)時差校正到所述當前繞射成像點的雙程垂直旅行時，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)；
[0014]對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)；
[0015]將所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)相加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)。
[0016]可選的，還包括:
[0017]響應(yīng)于所述當前炮集數(shù)據(jù)具有多條接收線且每條接收線具有多個炮檢距的接受傳感器，選用二維濾波器執(zhí)行所述對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理。
[0018]可選的，還包括:
[0019]響應(yīng)于得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)，將所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)與所述當前繞射成像點處的總繞射點成像數(shù)據(jù)相加，重新得到所述當前繞射成像點處的總繞射點成像數(shù)據(jù)；所述當前繞射成像點處的總繞射點成像數(shù)據(jù)表示已處理過的所有炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)之和；
[0020]響應(yīng)于得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)，返回執(zhí)行所述獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，直至得到所有炮集數(shù)據(jù)處理完畢，再執(zhí)行所述對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
[0021]此外，本發(fā)明還提供了一種繞射點成像的設(shè)備，該設(shè)備包括:
[0022]炮集選取模塊，用于獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)；
[0023]走時計算模塊，用于計算地震波傳播走時；
[0024]成像點選取模塊，用于選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點；
[0025]信號處理模塊，用于依據(jù)所述地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距處的所述當前炮集數(shù)據(jù)進行時差校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)；
[0026]成像點選取觸發(fā)模塊，用于返回觸發(fā)所述成像點選取模塊再選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點，直至得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)；
[0027]成像模塊，用于對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
[0028]可選的，所述信號處理模塊包括:[0029]時差校正子模塊，用于依據(jù)所述地震波傳播走時，將在各個炮檢距處的所述當前炮集數(shù)據(jù)時差校正到所述當前繞射成像點的雙程垂直旅行時，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)；
[0030]高通濾波子模塊，用于對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)。
[0031]疊加子模塊，用于將所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)相加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)。
[0032]可選的，還包括:
[0033]濾波器確定模塊，用于響應(yīng)于所述當前炮集數(shù)據(jù)具有多條接收線且每條接收線具有多個炮檢距的接受傳感器，選用二維濾波器執(zhí)行所述對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理。
[0034]可選的，還包括:
[0035]炮集疊加模塊，用于響應(yīng)于得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)，將所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)與所述當前繞射成像點處的總繞射點成像數(shù)據(jù)相加，重新得到所述當前繞射成像點處的總繞射點成像數(shù)據(jù)；所述當前繞射成像點處的總繞射點成像數(shù)據(jù)表示已處理過的所有炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)之和；
[0036]炮集選取觸發(fā)模塊，用于響應(yīng)于得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)，返回執(zhí)行所述獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，直至得到所有炮集數(shù)據(jù)處理完畢，再執(zhí)行所述對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
[0037]與現(xiàn)有技術(shù)相比，本發(fā)明具有以下優(yōu)點:
[0038]根據(jù)本發(fā)明實施例的技術(shù)方案，先對炮集數(shù)據(jù)進行時差校正、高通濾波及疊加的信號處理，然后再對處理后的數(shù)據(jù)進行成像。一方面，在繞射點處繞射波信號的相位會發(fā)生180度反轉(zhuǎn)、振幅迅速增大，呈現(xiàn)出高頻變化的特點，而在非繞射點處繞射波信號的振幅隨空間變化則非常緩慢，呈現(xiàn)出低頻變化的特點，另一方面，反射波在空間的連續(xù)性很強、振幅隨空間的變化也呈現(xiàn)出低頻變化的特點，因此，通過對炮集數(shù)據(jù)進行高通濾波，可以壓制作為低頻信號的反射波信號以及繞射點兩邊的繞射波信號，使作為高頻信號的繞射點處繞射波信號得以突出，從而使得利用突出繞射點的數(shù)據(jù)進行成像之后，可以更容易、更準確地識別出繞射點的位置。
【專利附圖】
【附圖說明】
[0039]為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明中記載的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動的前提下，還可以根據(jù)這些附圖獲得其他的附圖。
[0040]圖1為本發(fā)明中繞射點成像的方法實施例1的流程圖；
[0041]圖2為繞射波隨空間的振幅變化示意圖；
[0042]圖3為本發(fā)明實施例中一種高通濾波后的繞射波振幅變化示例示意圖；
[0043]圖4為本發(fā)明實施例中對炮集數(shù)據(jù)的信號處理過程一實施方式的流程圖；[0044]圖5為本發(fā)明中繞射點成像的方法實施例2的流程圖；
[0045]圖6為本發(fā)明中繞射點成像的設(shè)備實施例1的結(jié)構(gòu)圖；
[0046]圖7為本發(fā)明實施例中信號處理模塊604 —實施方式的結(jié)構(gòu)圖；
[0047]圖8為本發(fā)明中繞射點成像的設(shè)備實施例2的結(jié)構(gòu)圖。
【具體實施方式】
[0048]為了使本【技術(shù)領(lǐng)域】的人員更好地理解本申請方案，下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進行清楚、完整地描述，顯然，所描述的實施例僅是本申請一部分實施例，而不是全部的實施例。基于本發(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護的范圍。
[0049]發(fā)明人經(jīng)過研究發(fā)現(xiàn)，現(xiàn)有技術(shù)中之所以先分離反射波和繞射波再從繞射波中尋找繞射點位置，原因在于地震信號中的繞射波信號較弱而使得繞射點難以被發(fā)現(xiàn)。但是，由于繞射波與反射波在繞射點處是相切的，繞射軌跡的一段空間內(nèi)繞射波會被反射波干擾，可見，反射波與繞射波實際上是難以被準確地分開的。與此同時，發(fā)明人還發(fā)現(xiàn)，在繞射點處繞射波信號的相位會發(fā)生180度反轉(zhuǎn)、振幅迅速增大，呈現(xiàn)出高頻變化的特點，而在非繞射點處繞射波信號的振幅隨空間變化則非常緩慢，呈現(xiàn)出低頻變化的特點，并且，反射波在空間的連續(xù)性很強、振幅隨空間的變化也呈現(xiàn)出低頻變化的特點。因此，地震信號中，在繞射點處的繞射波信號與繞射點兩邊的繞射波信號、整個空間的反射波信號呈現(xiàn)出不同的頻譜特性，據(jù)此可以分離出繞射點處的繞射波信號。
[0050]基于上述分析，本發(fā)明的主要思想是:先對炮集數(shù)據(jù)進行時差校正、高通濾波及疊加的信號處理，然后再對處理后的數(shù)據(jù)進行成像。這樣，通過對炮集數(shù)據(jù)進行高通濾波，使作為低頻信號的反射波信號以及繞射點兩邊的繞射波信號得以壓制，并使作為高頻信號的繞射點處繞射波信號得以突出，從而使得利用突出繞射點的數(shù)據(jù)進行成像之后，可以更容易、更準確地識別出繞射點的位置。
[0051]在介紹了本發(fā)明的主要思想之后，下面結(jié)合附圖，通過實施例來詳細說明本發(fā)明中繞射點成像的方法和設(shè)備的實現(xiàn)方式。
[0052]參見圖1，示出了本發(fā)明中繞射點成像的方法實施例1的流程圖。在本實施例中，例如具體可以包括如下步驟:
[0053]S101、獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，并計算地震波傳播走時。
[0054]其中，本實施例中所指的一個炮集數(shù)據(jù)表示的是由同一個炮點激發(fā)地震波而在各個炮檢距的接收傳感器上接收的所有數(shù)據(jù)。一般地，地震勘探時會在多個不同的炮點多次激發(fā)地震波，每個炮點激發(fā)地震波以后可以在多個不同炮檢距的接收傳感器接收傳播來的地震信號，因此，每一個炮點就會對應(yīng)一組地震信號，即該炮點的炮集數(shù)據(jù)。其中，任意一個炮點的炮集數(shù)據(jù)都可以全面地反映地下繞射點的情況，即可利用任一炮集數(shù)據(jù)進行繞射點成像?？梢岳斫獾氖?，本實施例中，可以選用未處理過的炮集數(shù)據(jù)來作為當前炮集數(shù)據(jù)對繞射點進行成像，此時，在通過疊加多次繞射點成像而放大繞射點成像結(jié)果時，就可以利用不同炮集數(shù)據(jù)對同一地質(zhì)體進行多次繞射點成像，從而減小激發(fā)位置對繞射點成像的影響。
[0055]在獲取當前炮集數(shù)據(jù)之后，由于各個炮檢距之間的距離是不同的，各個炮檢距的接收傳感器接收到的地震信號數(shù)據(jù)存在時差。為了后續(xù)對各個炮檢距處的數(shù)據(jù)進行信號處理時消除數(shù)據(jù)之間的時差，此時需要計算地震波傳播走時。例如具體地，可以計算一個地震波傳播走時表，即用表格形式表現(xiàn)的地震波傳播距離與所需時間的關(guān)系。
[0056]S102、選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點。
[0057]在繞射點成像過程中，需要在地下區(qū)域選取一個成像點假定為繞射點，然后對這個假定的繞射點進行成像，本實施例中所選取的當前繞射成像點即是這一假定的繞射點。
[0058]需要說明的是，在對一個炮集數(shù)據(jù)進行成像時，需要使成像點遍歷所需勘探的整個地下區(qū)域，為此，本實施例中需要針對整個地下區(qū)域的每個成像點而對數(shù)據(jù)進行信號處理?？梢岳斫獾氖牵捎谛枰闅v的地下區(qū)域是一個三維空間，該三維空間中的成像點位置可以是通過地面位置與地下深度位置兩者來描述。因此，在本實施例的一些實施方式中，可以把地面所在的平面劃分成大量很小的網(wǎng)格點，而在選取當前繞射成像點，可以先選擇一個未信號處理過的網(wǎng)格點，然后遍歷該網(wǎng)格點所有的地下深度位置來分別進行信號處理，再選擇一個未信號處理過網(wǎng)格點，直至遍歷所有的網(wǎng)格點為止，以此來完成遍歷地下區(qū)域所有的成像點。其中，地下深度位置可以通過成像點處的地震波傳播到地面所需的時間來表示，具體地，可以是以成像點處的地震波垂直旅行時。
[0059]S103、依據(jù)所述地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距處的所述當前炮集數(shù)據(jù)進行時差校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)。
[0060]本實施例中，對數(shù)據(jù)的信號處理至少包括時差校正、高通濾波和疊加。其中，對不同炮檢距的數(shù)據(jù)進行時差校正是為了校正各個炮檢距的炮集數(shù)據(jù)之間存在的時差，從而使得各個炮檢距的炮集數(shù)據(jù)可以進行疊加。而對不同炮檢距的數(shù)據(jù)進行疊加則是為了提高信噪比而消除繞射點兩邊的噪聲信號。
[0061]需要說明的是，通常的成像過程中，主要是通過對炮集數(shù)據(jù)進行時差校正、疊加和偏移來完成的。在沿繞射軌跡求和的疊加過程中，由于繞射點兩邊的繞射波信號振幅符號相反(如圖2所示)而互相抵消、互相壓制，使得繞射點處的繞射波信號得以保留。但是，由于通常炮集數(shù)據(jù)的采集是關(guān)于繞射點不對稱的，因此，繞射點以外的繞射波信號在成像過程中實際上并不能完全抵消，這樣會使繞射點成像信噪比降低，影響成像質(zhì)量，這也是反射成像數(shù)據(jù)中繞射信息弱的原因之一。而本實施例中，為了突出繞射點處的繞射波信號，可以先對時差校正過的炮集數(shù)據(jù)進行高通濾波，然后再對濾波后的數(shù)據(jù)進行疊加。由于繞射波信號的振幅僅在繞射點處處于高頻變化而在繞射點兩邊處于低頻變化，同時反射波的振幅隨空間變化一直處于低頻，因此，經(jīng)過高通濾波的數(shù)據(jù)中反射波信號以及繞射點兩邊的繞射波都會被壓制，而繞射點處的繞射波信號因被保留而得以突出。例如，圖3示出了一種高通濾波后的繞射波振幅變化示例示意圖，從圖3中可以看出，濾波后繞射波信號僅在繞射點處具有振幅變化。其中，本實施例中所涉及的高通濾波，表示的是具有使低頻受到壓制而使高頻相對放大或通過的信號濾波處理方式。
[0062]可以理解的是，在本實施例的一些實施方式中，可以先對炮集數(shù)據(jù)進行時差校正，然后對校正過的數(shù)據(jù)進行高通濾波，再對濾波后的數(shù)據(jù)進行疊加，從而得到成像數(shù)據(jù)。其中，時差校正具體可以是將各個炮檢距處的炮集數(shù)據(jù)都校正到當前繞射成像點的雙程垂直旅行時。具體地，如圖4所示，對炮集數(shù)據(jù)的信號處理過程，例如可以包括:[0063]S401、依據(jù)所述地震波傳播走時，將在各個炮檢距處的所述當前炮集數(shù)據(jù)時差校正到所述當前繞射成像點的雙程垂直旅行時，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)。
[0064]由于此時的信號處理是要為了針對當前繞射成像點進行成像，因此，可以將各個炮檢距處的炮集數(shù)據(jù)都校正到該當前繞射成像點的雙程垂直旅行時，從而實現(xiàn)將曲線拉平為直線?？梢岳斫獾氖?，由于各個炮檢距處所接收到的數(shù)據(jù)需要校正的時差是不同的，因此，對炮集數(shù)據(jù)的時差校正也稱為動校正。
[0065]S402、對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)。
[0066]其中，對校正數(shù)據(jù)的高通濾波處理，可以是通過預(yù)先設(shè)定為高通濾波器的數(shù)字濾波器來實現(xiàn)，該高通濾波器具體是指具有使低頻受到壓制而使高頻相對放大或通過的濾波器?？梢岳斫獾氖?，對于不同的地質(zhì)，地震信號的振幅空間變化可能有很大差別，并且，炮集數(shù)據(jù)采集的方式不同，也可能造成不同的地震信號的處理方式不同，因此，在對校正數(shù)據(jù)進行高通濾波處理時，可以依據(jù)炮集數(shù)據(jù)的情況來選擇不同設(shè)置的數(shù)字濾波器。例如，一些實施方式中，還可以響應(yīng)于所述當前炮集數(shù)據(jù)具有多條接收線且每條接收線具有多個炮檢距的接受傳感器，選用二維濾波器執(zhí)行所述對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理。
[0067]S403、將所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)相加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)。
[0068]在對炮集數(shù)據(jù)進行高通濾波之后，可以將各個炮檢距處的濾波后數(shù)據(jù)相加而得到繞射點成像數(shù)據(jù)，以提高成像數(shù)據(jù)的信噪比，使得繞射點處的信號通過疊加得以進一步突出，而使繞射點以外的噪聲信號通過疊加得以消除。
[0069]接著返回圖1。在S103執(zhí)彳丁完成以后,接著進入S104的執(zhí)打。
[0070]S104、返回執(zhí)行所述選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點，直至得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)。
[0071]在繞射點成像過程中，需要遍歷所有的成像點來獲得各個成像點的繞射點成像數(shù)據(jù)。因此，在得到所述當前炮集數(shù)據(jù)在當前繞射成像點處的繞射點成像數(shù)據(jù)之后，還需要返回S102來重新選取一個未針對該炮集數(shù)據(jù)選用過的成像點，從而得到所述當前炮集數(shù)據(jù)在重新選取的成像點處的繞射點成像數(shù)據(jù)，直至得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)為止。
[0072]需要說明的是，在前述以網(wǎng)格點及深度位置表示成像點位置的實施方式基礎(chǔ)上，可以先針對同一個網(wǎng)格點選擇不同的深度位置來獲得繞射點成像數(shù)據(jù)，在遍歷該網(wǎng)格點的所有深度位置之后，再選擇另一個網(wǎng)格點來獲得繞射點成像數(shù)據(jù)，直至遍歷所有網(wǎng)格點，從而得到各個深度位置的各個網(wǎng)格點處的繞射點成像數(shù)據(jù)。
[0073]S105、對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
[0074]可以理解的是，由于繞射點成像數(shù)據(jù)經(jīng)過了高通濾波，因此，繞射點兩邊的信號可以得到壓制，而繞射點處的信號得以突出，從而使得成像結(jié)果中的繞射點位置得以突出。
[0075]在本實施的一些實施方式中，由于單一炮點的炮集數(shù)據(jù)中繞射點處的繞射波信號往往較弱，即使經(jīng)過了前述高通濾波處理，繞射點處的繞射波信號可能依然難以與繞射點以外的噪聲信號區(qū)分。為了進一步提高繞射點處信號的信噪比，可以將多個不同炮點的炮集數(shù)據(jù)經(jīng)過高通濾波之后而在同一成像點處的成像數(shù)據(jù)進行疊加，從而使得繞射點處的信號進一步突出，繞射點以外的噪聲信號進一步消除。具體地，例如在S103執(zhí)行完成之后，還可以響應(yīng)于得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的總繞射點成像數(shù)據(jù)，將所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)與所述當前繞射成像點處的繞射點成像數(shù)據(jù)相加，重新得到所述當前繞射成像點處的總繞射點成像數(shù)據(jù)；而在S104執(zhí)行完成之后，還可以響應(yīng)于得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)，返回執(zhí)行SlOl獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，直至得到所有炮集數(shù)據(jù)處理完畢，再執(zhí)行S105對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。其中，所述當前繞射成像點處的總繞射點成像數(shù)據(jù)表示已處理過的所有炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)之和?？梢岳斫獾氖?，此時S105用于繞射點成像的數(shù)據(jù)為各個成像點的總繞射點成像數(shù)據(jù)。
[0076]本實施例的技術(shù)方案，通過對炮集數(shù)據(jù)進行高通濾波，使作為低頻信號的反射波信號以及繞射點兩邊的繞射波信號得以壓制，并使作為高頻信號的繞射點處繞射波信號得以突出，從而使得利用突出繞射點的數(shù)據(jù)進行成像之后，可以更容易、更準確地識別出繞射點的位置。
[0077]為了使本領(lǐng)域技術(shù)人員更加清楚地理解本發(fā)明的應(yīng)用方式，下面以一個應(yīng)用場景來對本發(fā)明的實施方式進行示例性地介紹。在該應(yīng)用場景中，繞射點的成像是利用多個不同炮點的炮集數(shù)據(jù)疊加而實現(xiàn)的，而成像點是由預(yù)先離散化地質(zhì)體而得到的網(wǎng)格點和深度位置來表示的。
[0078]參見圖5，示出了本發(fā)明中繞射點成像的方法實施例2的流程圖。在本實施例中，例如具體可以包括如下步驟:
[0079]S501、獲取用于繞射點成像的所有炮集數(shù)據(jù)，計算地震波傳播走時表，并選擇高通濾波器。
[0080]S502、獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)。
[0081 ] S503、選取一個未針對當前炮集數(shù)據(jù)而選用過的網(wǎng)格點作為當前網(wǎng)格點。
[0082]S504、選取一個未針對當前炮集數(shù)據(jù)及當前網(wǎng)格點而選用過的深度作為當前網(wǎng)格點的當前深度。
[0083]S505、依據(jù)地震波傳播走時，將在各個炮檢距處的當前炮集數(shù)據(jù)時差校正到在當前深度下當前網(wǎng)格點的雙程垂直旅行時，得到當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)。
[0084]S506、采用對當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理，得到當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)。
[0085]S507、將當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)相加，得到當前炮集數(shù)據(jù)在當前深度下當前網(wǎng)格點的繞射點成像數(shù)據(jù)并存儲。
[0086]S508、將當前炮集數(shù)據(jù)在當前深度下當前網(wǎng)格點的繞射點成像數(shù)據(jù)與當前深度下當前網(wǎng)格點的總繞射點成像數(shù)據(jù)相加，重新得到當前深度下當前網(wǎng)格點的總繞射點成像數(shù)據(jù)。
[0087]S509、迭代S504?S508，直至對當前網(wǎng)格點遍歷所有深度。
[0088]S510、迭代S503?S509，直至對當前炮集數(shù)據(jù)遍歷所有當前網(wǎng)格點。[0089]S511、迭代S502?S510，直至遍歷所有炮集數(shù)據(jù)。
[0090]S512、對各深度下各個網(wǎng)格點的總繞射點成像數(shù)據(jù)進行繞射點成像。
[0091]本實施例的技術(shù)方案，一方面，通過對炮集數(shù)據(jù)進行高通濾波，使作為低頻信號的反射波信號以及繞射點兩邊的繞射波信號得以壓制，并使作為高頻信號的繞射點處繞射波信號得以突出，另一方面，通過將多個不同炮點的炮集數(shù)據(jù)在同一成像點的成像數(shù)據(jù)進行了疊加，提高了繞射點處信號的信噪比，進一步突出了繞射點處的信號，從而使得利用突出繞射點的數(shù)據(jù)進行成像之后，可以更容易、更準確地識別出繞射點的位置。
[0092]對應(yīng)于方法實施例，本發(fā)明還提供了一種繞射點成像設(shè)備。
[0093]參見圖6，示出了本發(fā)明中繞射點成像的設(shè)備實施例1的結(jié)構(gòu)圖。在本實施例中，所述設(shè)備例如具體可以包括:
[0094]炮集選取模塊601，用于獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)；
[0095]走時計算模塊602,用于計算地震波傳播走時；
[0096]成像點選取模塊603，用于選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點；
[0097]信號處理模塊604，用于依據(jù)所述地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距處的所述當前炮集數(shù)據(jù)進行時差校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)；
[0098]成像點選取觸發(fā)模塊605，用于返回觸發(fā)所述成像點選取模塊再選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點，直至得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)；
[0099]成像模塊606，用于對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
[0100]參見圖7，示出了本發(fā)明實施例中信號處理模塊604—實施方式的結(jié)構(gòu)圖。在本實施方式中，所述信號處理模塊604例如可以包括:
[0101]時差校正子模塊701，用于依據(jù)所述地震波傳播走時，將在各個炮檢距處的所述當前炮集數(shù)據(jù)時差校正到所述當前繞射成像點的雙程垂直旅行時，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)；
[0102]高通濾波子模塊702，用于對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)；
[0103]疊加子模塊703，用于將所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)相力口，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)。
[0104]進一步而言，可選的，在圖7所示的實施方式基礎(chǔ)上，所述設(shè)備例如還可以包括:
[0105]濾波器確定模塊，用于響應(yīng)于所述當前炮集數(shù)據(jù)具有多條接收線且每條接收線具有多個炮檢距的接受傳感器，選用二維濾波器執(zhí)行所述對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理。
[0106]參見圖8，示出了本發(fā)明中繞射點成像的設(shè)備實施例2的結(jié)構(gòu)圖。在本實施例中，除了圖6所示的所有結(jié)構(gòu)以外，所述設(shè)備例如還可以包括:
[0107]炮集疊加模塊801，用于響應(yīng)于得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)，將所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)與所述當前繞射成像點處的總繞射點成像數(shù)據(jù)相加，重新得到所述當前繞射成像點處的總繞射點成像數(shù)據(jù)；所述當前繞射成像點處的總繞射點成像數(shù)據(jù)表示已處理過的所有炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)之和；
[0108] 炮集選取觸發(fā)模塊802，用于響應(yīng)于得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)，返回執(zhí)行所述獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，直至得到所有炮集數(shù)據(jù)處理完畢，再執(zhí)行所述對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
[0109]本發(fā)明設(shè)備實施例的技術(shù)方案，通過對炮集數(shù)據(jù)進行高通濾波，使作為低頻信號的反射波信號以及繞射點兩邊的繞射波信號得以壓制，并使作為高頻信號的繞射點處繞射波信號得以突出，從而使得利用突出繞射點的數(shù)據(jù)進行成像之后，可以更容易、更準確地識別出繞射點的位置。
[0110]需要說明的是，在本文中，諸如第一和第二等之類的關(guān)系術(shù)語僅僅用來將一個實體或者操作與另一個實體或操作區(qū)分開來，而不一定要求或者暗示這些實體或操作之間存在任何這種實際的關(guān)系或者順序。術(shù)語“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過程、方法、物品或者設(shè)備不僅包括那些要素，而且還包括沒有明確列出的其他要素，或者是還包括為這種過程、方法、物品或者設(shè)備所固有的要素。在沒有更多限制的情況下，由語句“包括一個……”限定的要素，并不排除在包括所述要素的過程、方法、物品或者設(shè)備中還存在另外的相同要素。
[0111]對于設(shè)備實施例而言，由于其基本對應(yīng)于方法實施例，所以相關(guān)之處參見方法實施例的部分說明即可。以上所描述的設(shè)備實施例僅僅是示意性的，其中所述作為分離部件說明的單元可以是或者也可以不是物理上分開的，作為單元顯示的部件可以是或者也可以不是物理單元，即可以位于一個地方，或者也可以分布到多個網(wǎng)絡(luò)單元上?？梢愿鶕?jù)實際的需要選擇其中的部分或者全部模塊來實現(xiàn)本實施例方案的目的。本領(lǐng)域普通技術(shù)人員在不付出創(chuàng)造性勞動的情況下，即可以理解并實施。
[0112]以上所述僅是本申請的【具體實施方式】，應(yīng)當指出，對于本【技術(shù)領(lǐng)域】的普通技術(shù)人員來說，在不脫離本申請原理的前提下，還可以做出若干改進和潤飾，這些改進和潤飾也應(yīng)視為本申請的保護范圍。
【權(quán)利要求】
1.一種繞射點成像的方法，其特征在于，包括: 獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，并計算地震波傳播走時； 選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點； 依據(jù)所述地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距處的所述當前炮集數(shù)據(jù)進行時差校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)； 返回執(zhí)行所述選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點，直至得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)； 對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述依據(jù)各個炮檢距的地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距的所述當前炮集數(shù)據(jù)進行動校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)，包括: 依據(jù)所述地震波傳播走時，將在各個炮檢距處的所述當前炮集數(shù)據(jù)時差校正到所述當前繞射成像點的雙程垂直旅行時，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)；對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)； 將所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)相加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)。
3.根據(jù)權(quán)利要求2所述的方法，其特征在于，還包括: 響應(yīng)于所述當前炮集數(shù)據(jù)具有多條接收線且每條接收線具有多個炮檢距的接受傳感器，選用二維濾波器執(zhí)行所述對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理。
4.根據(jù)權(quán)利要求1所述的方法，其特征在于，還包括: 響應(yīng)于得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)，將所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)與所述當前繞射成像點處的總繞射點成像數(shù)據(jù)相加，重新得到所述當前繞射成像點處的總繞射點成像數(shù)據(jù)；所述當前繞射成像點處的總繞射點成像數(shù)據(jù)表示已處理過的所有炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)之和； 響應(yīng)于得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)，返回執(zhí)行所述獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，直至得到所有炮集數(shù)據(jù)處理完畢，再執(zhí)行所述對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
5.一種繞射點成像的設(shè)備，其特征在于，包括: 炮集選取模塊，用于獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)； 走時計算模塊，用于計算地震波傳播走時； 成像點選取模塊，用于選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點； 信號處理模塊，用于依據(jù)所述地震波傳播走時，針對所述當前繞射成像點的地下深度，對在各個炮檢距處的所述當前炮集數(shù)據(jù)進行時差校正、高通濾波及疊加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)； 成像點選取觸發(fā)模塊，用于返回觸發(fā)所述成像點選取模塊再選取一個未針對所述當前炮集數(shù)據(jù)而選用過的成像點作為當前繞射成像點，直至得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)； 成像模塊，用于對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
6.根據(jù)權(quán)利要求5所述的設(shè)備，其特征在于，所述信號處理模塊包括: 時差校正子模塊，用于依據(jù)所述地震波傳播走時，將在各個炮檢距處的所述當前炮集數(shù)據(jù)時差校正到所 述當前繞射成像點的雙程垂直旅行時，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)； 高通濾波子模塊，用于對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理，得到所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)。 疊加子模塊，用于將所述當前炮集數(shù)據(jù)在各個炮檢距處的繞射點成像數(shù)據(jù)相加，得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)。
7.根據(jù)權(quán)利要求6所述的設(shè)備，其特征在于，還包括: 濾波器確定模塊，用于響應(yīng)于所述當前炮集數(shù)據(jù)具有多條接收線且每條接收線具有多個炮檢距的接受傳感器，選用二維濾波器執(zhí)行所述對所述當前炮集數(shù)據(jù)在各個炮檢距處的校正數(shù)據(jù)進行高通濾波處理。
8.根據(jù)權(quán)利要求5所述的設(shè)備，其特征在于，還包括: 炮集疊加模塊，用于響應(yīng)于得到所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)，將所述當前炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)與所述當前繞射成像點處的總繞射點成像數(shù)據(jù)相加，重新得到所述當前繞射成像點處的總繞射點成像數(shù)據(jù)；所述當前繞射成像點處的總繞射點成像數(shù)據(jù)表示已處理過的所有炮集數(shù)據(jù)在所述當前繞射成像點處的繞射點成像數(shù)據(jù)之和； 炮集選取觸發(fā)模塊，用于響應(yīng)于得到所述當前炮集數(shù)據(jù)在各個成像點處的繞射點成像數(shù)據(jù)，返回執(zhí)行所述獲取一個未處理過的炮集數(shù)據(jù)作為當前炮集數(shù)據(jù)，直至得到所有炮集數(shù)據(jù)處理完畢，再執(zhí)行所述對各個成像點處的繞射點成像數(shù)據(jù)進行繞射點成像。
【文檔編號】G01V1/28GK103941289SQ201410195993
【公開日】2014年7月23日 申請日期:2014年5月9日 優(yōu)先權(quán)日:2014年5月9日
【發(fā)明者】王真理 申請人:中國科學(xué)院地質(zhì)與地球物理研究所, 北京中科聯(lián)衡科技有限公司