本發(fā)明涉及地震數(shù)據(jù)處理技術(shù)領(lǐng)域，尤其涉及一種采集腳印壓制方法及系統(tǒng)。

背景技術(shù)：

隨著巖性油氣藏和隱蔽性油氣藏勘探的不斷深入，若想對(duì)油氣藏和隱蔽性油氣藏的分布得出準(zhǔn)確的油藏評(píng)價(jià)結(jié)果和儲(chǔ)層預(yù)測(cè)方案，需要采集準(zhǔn)確的地震數(shù)據(jù)，因此，在采集地震數(shù)據(jù)的過(guò)程中，對(duì)地震成像分辨率和地震資料的細(xì)節(jié)刻畫(huà)提出了較高的要求。采集腳印噪聲又稱采集痕跡噪聲，是一種由于人為因素產(chǎn)生的地震噪聲，具體是指在地震資料采集與處理過(guò)程中留下的人為痕跡，表現(xiàn)方式為地震剖面或切片上不同面源位置處規(guī)則性的振幅變化假象。由于采集腳印噪聲對(duì)儲(chǔ)層預(yù)測(cè)有一定的影響，它的存在會(huì)模糊砂體邊界和構(gòu)造特征，因此，該采集腳印噪聲會(huì)影響地震成像及地震屬性分析，進(jìn)而造成一些巖性體分布的假象，對(duì)油藏評(píng)價(jià)和儲(chǔ)層預(yù)測(cè)產(chǎn)生了影響。

目前，經(jīng)典的壓制采集腳印方法是指在傳統(tǒng)道集疊加后的地震數(shù)據(jù)體上，應(yīng)用三維錐形濾波等地震處理技術(shù)來(lái)消除地震數(shù)據(jù)中以垂直和水平條紋形式出現(xiàn)的采集腳印噪聲。假設(shè)三維地震數(shù)據(jù)體的炮線和檢波線是正交的，抽取一條炮線和檢波線上的所有地震道組成一個(gè)正交子集，相當(dāng)于對(duì)地下局部三維地質(zhì)體進(jìn)行單次覆蓋地震采集，該正交子集所構(gòu)建的時(shí)空域稱為十字交叉排列域，其炮線和檢波線的交點(diǎn)稱為交叉點(diǎn)。具體的，首先將三維地震數(shù)據(jù)體抽取到十字交叉排列域，其次，利用三維錐形濾波技術(shù)來(lái)壓制地震數(shù)據(jù)采集過(guò)程中出現(xiàn)的線性噪聲，突出有效信號(hào)，進(jìn)而壓制掉一部分的采集腳印噪聲。

然而，上述三維錐形濾波技術(shù)中對(duì)于采集腳印噪聲的壓制有其局限性，無(wú)法徹底消除采集腳印噪聲的影響，而且由炮線和檢波線組成的十字交叉排列的稀疏性會(huì)影響去噪的保真度，同時(shí)會(huì)產(chǎn)生空間假頻，不符合高保真處理要求，所取得的結(jié)果也很難完全滿足振幅、頻率、相位等屬性在油氣識(shí)別中 的需要。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明提供一種采集腳印壓制方法及系統(tǒng)，以解決現(xiàn)有采集腳印壓制方法無(wú)法完整消除采集腳印的影響，且不符合高保真處理要求致使地震數(shù)據(jù)處理結(jié)果不能完全滿足振幅、頻率、相位等屬性在油氣識(shí)別中需要的問(wèn)題。

本發(fā)明提供的一種采集腳印壓制方法，包括：

將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內(nèi)，得到N個(gè)共偏移距共方位角域道集，所述共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數(shù)，N為大于或等于1的整數(shù)；

分別將所述炮線之間的距離和所述檢波線之間的距離進(jìn)行減半處理，根據(jù)每個(gè)所述共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，每個(gè)所述共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集；

對(duì)每個(gè)所述子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)，所述加密數(shù)據(jù)包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)；

將所述加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到4N個(gè)子共偏移距共方位角域道集；

對(duì)所述4N個(gè)子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

本發(fā)明還提供一種采集腳印壓制系統(tǒng)，包括：

獲取模塊，用于將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內(nèi)，得到N個(gè)共偏移距共方位角域道集，所述共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數(shù)，N為大于或等于1的整數(shù)；

處理模塊，用于分別將所述炮線之間的距離和所述檢波線之間的距離進(jìn)行減半處理，根據(jù)每個(gè)所述共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，每個(gè)所述共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集；

插值模塊，用于對(duì)每個(gè)所述子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)，所述加密數(shù)據(jù)包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)；

所述獲取模塊，還用于將所述加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到4N個(gè)子共偏移距共方位角域道集；

所述處理模塊，還用于對(duì)所述4N個(gè)子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

本發(fā)明采集腳印壓制方法及系統(tǒng)，通過(guò)將待處理的地震數(shù)據(jù)分選到十字排列域內(nèi)，并根據(jù)地震道與偏移距和方位角的關(guān)系，得到了多個(gè)共偏移距共方位角域道集，通過(guò)分別將炮線距和檢波線距進(jìn)行減半處理，且根據(jù)每個(gè)共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，可計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，進(jìn)而使每個(gè)共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集，分別對(duì)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)的加密數(shù)據(jù)，最后將該加密數(shù)據(jù)重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到四倍于共偏移距共方位角域道集的子共偏移距共方位角域道集，對(duì)所有的子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，便可得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。本發(fā)明實(shí)施例的技術(shù)方案，利用共偏移距共方位角域道集的單次覆蓋特性，通過(guò)在同一偏移距和同一方位角內(nèi)進(jìn)行插值，實(shí)現(xiàn)了加密炮線和檢波線的目的，實(shí)現(xiàn)方案簡(jiǎn)單、容易實(shí)現(xiàn)，對(duì)壓制采集腳印、淺層成像和高保真處理有重要的意義，符合高保真處理要求，所取得的結(jié)果也完全滿足振幅、頻率、相位等屬性在油氣識(shí)別中的需要。

附圖說(shuō)明

為了更清楚地說(shuō)明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作一簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1為本發(fā)明提供的采集腳印壓制方法實(shí)施例一的流程示意圖；

圖2為某一十字排列域及該十字排列域上的共偏移距共方位角域道集；

圖3為從十字排列域中抽取的某一共偏移距共方位角域道集；

圖4為炮線距和檢波線距減半前形成的十字排列域；

圖5為炮線距和檢波線距減半后形成的十字排列域；

圖6為共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值前的數(shù)據(jù)分布示意圖；

圖7為共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值后的數(shù)據(jù)分布示意圖；

圖8為原始地震數(shù)據(jù)對(duì)共偏移距共方位角域道集的劃分示意圖；

圖9為加密數(shù)據(jù)對(duì)共偏移距共方位角域道集的劃分示意圖；

圖10為原始地震數(shù)據(jù)在共偏移距共方位角域道集上形成的項(xiàng)鏈圖；

圖11為加密數(shù)據(jù)在共偏移距共方位角域道集上形成的項(xiàng)鏈圖；

圖12為原始地震數(shù)據(jù)在共偏移距共方位角域道集上的偏移距范圍；

圖13為加密數(shù)據(jù)分別在相應(yīng)的四個(gè)子共偏移距共方位角域道集上的偏移距范圍；

圖14為原始地震數(shù)據(jù)的共偏移距偏移處理成果切片得到的地震數(shù)據(jù)分析結(jié)果；

圖15為加密數(shù)據(jù)的共偏移距偏移處理成果切片得到的地震數(shù)據(jù)分析結(jié)果；

圖16為本發(fā)明提供的采集腳印壓制系統(tǒng)實(shí)施例一的結(jié)構(gòu)示意圖。

具體實(shí)施方式

為使本發(fā)明實(shí)施例的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚，下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有作出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。

三維地震勘探是一種高密度面積采集技術(shù)，是三維體積勘探。它利用炮點(diǎn)和檢波點(diǎn)網(wǎng)格的靈活組合獲得分布均勻的地下CDP點(diǎn)網(wǎng)格和確定的覆蓋次數(shù)。在三維地震數(shù)據(jù)采集中，炮檢點(diǎn)有規(guī)律地變化，會(huì)造成在疊加或偏移后的地震數(shù)據(jù)上地震波振幅出現(xiàn)有規(guī)律性的變化，在時(shí)間和深度切片上表現(xiàn)為振幅周期性變化，這種周期性振幅假象就是采集腳印。采集腳印對(duì)儲(chǔ)層預(yù)測(cè)有一定的影響，它的存在模糊了砂體邊界和構(gòu)造特征，造成一些巖性體分布 的假象，對(duì)油藏評(píng)價(jià)和儲(chǔ)層預(yù)測(cè)了產(chǎn)生影響。

觀測(cè)系統(tǒng)參數(shù)的選擇特別是炮線距和檢波線距的選擇，會(huì)影響采集腳印現(xiàn)象的產(chǎn)生，若炮線距和檢波線距過(guò)大，則會(huì)產(chǎn)生較嚴(yán)重的采集腳印，并且炮線距的影響大于檢波線距的影響。在實(shí)際應(yīng)用中，炮線距和檢波線距一般選擇為道距的4～6倍。

野外三維地震數(shù)據(jù)采集一旦完成，炮線距、檢波線距和道距就已確定。因此，采集到的地震數(shù)據(jù)中是否存在采集腳印主要取決于炮檢線距的大小。但是，一般采集地震數(shù)據(jù)時(shí)為了節(jié)省成本，通常采用較大炮檢線距，因此會(huì)出現(xiàn)采集腳印。

現(xiàn)階段，經(jīng)典的采集腳印壓制方法是在傳統(tǒng)地震道集疊加后的地震數(shù)據(jù)體上，應(yīng)用三維錐形濾波等地震處理技術(shù)，消除地震數(shù)據(jù)中以垂直和水平條紋形式出現(xiàn)的采集腳印噪聲。實(shí)際上，三維錐形濾波地震處理技術(shù)是在采集腳印已進(jìn)入成像過(guò)程中的“事后處理”思路，其對(duì)于采集腳印現(xiàn)象的壓制有一定的局限性，無(wú)法通過(guò)處理完整消除采集腳印的影響，而且不符合高保真處理要求，所取得的結(jié)果很難完全滿足振幅、頻率、相位等屬性在油氣識(shí)別中的需要。

針對(duì)現(xiàn)有采集腳印壓制方法存在的問(wèn)題和缺陷，本發(fā)明提供了一種采集腳印壓制方法及系統(tǒng)，通過(guò)在共偏移距共方位角域內(nèi)插炮檢線，將炮線距和檢波線距縮小一倍，通過(guò)在共偏移距共方位角域內(nèi)進(jìn)行高精度插值，并按照新的炮線距和檢波線距重新抽取共偏移距共方位角域道集，獲得的道集個(gè)數(shù)是原來(lái)的4倍，即覆蓋次數(shù)增大了原來(lái)的4倍，而且每個(gè)道集的偏移距范圍縮減為原來(lái)的一半。這種高密度、小線距、小偏移距范圍的共偏移距共方位角域道集在地震數(shù)據(jù)疊前偏移時(shí)能有效阻止采集腳印進(jìn)入到成像，是一種“源頭控制”的“事前壓制”思路，解決了現(xiàn)有三維錐形濾波技術(shù)存在的消除采集腳印噪聲不徹底，且不符合高保真處理要求造成的所取得的結(jié)果難以滿足振幅、頻率、相位等屬性在油氣識(shí)別中需要的問(wèn)題。

圖1為本發(fā)明提供的采集腳印壓制方法實(shí)施例一的流程示意圖。如圖1所示，本發(fā)明實(shí)施例一提供的采集腳印壓制方法，包括：

步驟101：將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內(nèi)，得到N個(gè)共偏移距共方位角域道集；

其中，共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數(shù)，N為大于或等于1的整數(shù)。

具體的，根據(jù)地震資料從觀測(cè)系統(tǒng)中抽取一條炮線和檢波線，炮線由一列炮點(diǎn)組成，檢波線由一列檢波點(diǎn)組成。炮線和檢波線交叉形成十字排列域，交叉的方式可以是正交，也可以是斜交，本發(fā)明實(shí)施例并不對(duì)此進(jìn)行限定。因此，將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到上述十字排列域內(nèi)，可以得到多個(gè)共偏移距共方位角域道集，該共偏移距共方位角域道集的數(shù)量與具有相同偏移距和相同方位角的區(qū)域被覆蓋的次數(shù)相關(guān)。

實(shí)際應(yīng)用中，覆蓋次數(shù)是指炮點(diǎn)對(duì)地下反射界面同一個(gè)點(diǎn)的追蹤次數(shù)，反映到十字排列域就是地震數(shù)據(jù)分選的每個(gè)網(wǎng)格的次數(shù)。因此，共偏移距共方位角域內(nèi)的每個(gè)地震道的數(shù)據(jù)的集合組成共偏移距共方位角域道集。

由上述可知，十字排列域包含了多個(gè)共偏移距共方位角域道集，從該十字排列域中抽取任一共偏移距共方位角域道集，該第一共偏移距共方位角域道集在空間展開(kāi)后具有單次覆蓋的特征。

本發(fā)明實(shí)施例主要是對(duì)具有單次覆蓋特征的共偏移距共方位角域道集進(jìn)行處理來(lái)減低或消除采集腳印噪聲。

步驟102：分別將炮線之間的距離和檢波線之間的距離進(jìn)行減半處理，根據(jù)每個(gè)共偏移距共方位角域道集中的地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，每個(gè)共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集；

將組成十字排列域的炮線距和檢波線距分別減半處理，并根據(jù)該共偏移距共方位角域道集中每個(gè)地震道，提取出該地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，根據(jù)該地震道的方位角及其對(duì)應(yīng)的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，可以計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)。處理完成后，可使每個(gè)共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集。

步驟103：對(duì)每個(gè)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)；

其中，該加密數(shù)據(jù)包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)。

由于同一共偏移距共方位角域道集，其具有相同的方位角，因此，在子 共偏移距共方位角域道集進(jìn)行插值得到的加密數(shù)據(jù)，其一致性顯著好于直接在共偏移距共方位角域道集內(nèi)進(jìn)行插值得到的數(shù)據(jù)，因此，按照本發(fā)明實(shí)施例提供的方法進(jìn)行數(shù)據(jù)插值的精度明顯高于常規(guī)的共偏移距域插值的精度。

步驟104：將上述加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到4N個(gè)子共偏移距共方位角域道集；

本實(shí)施例將加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，重新抽取共偏移距共方位角域道集，共獲得4N個(gè)子共偏移距共方位角域道集，該道集個(gè)數(shù)是原來(lái)的4倍，即具有相同偏移距和相同方位角區(qū)域的覆蓋次數(shù)增大了原來(lái)的4倍。這種高密度、小線距、小偏移距范圍的子共偏移距共方位角域道集，等價(jià)于炮線距和檢波線距減半的野外采集。

由步驟103可知，插值后的數(shù)據(jù)加上原始數(shù)據(jù)稱為加密數(shù)據(jù)。加密數(shù)據(jù)按新的炮線距和檢波線距重新劃分共偏移距共方位角域，再抽取共偏移距共方位角域道集。由于炮線距和檢波線距縮小為原來(lái)的一倍，原來(lái)的一個(gè)共偏移距共方位角域一化為四。

步驟105：對(duì)上述4N個(gè)子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

具體的，對(duì)上述4N個(gè)子共偏移距共方位角域道集重新進(jìn)行偏移、疊加，再提取振幅切片，進(jìn)而可得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

本發(fā)明實(shí)施例一提供的采集腳印壓制方法，通過(guò)將待處理的地震數(shù)據(jù)分選到十字排列域內(nèi)，并根據(jù)地震道與偏移距和方位角的關(guān)系，得到了多個(gè)共偏移距共方位角域道集，通過(guò)分別將炮線距和檢波線距進(jìn)行減半處理，且根據(jù)每個(gè)共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，可計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，進(jìn)而使每個(gè)共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集，分別對(duì)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)的加密數(shù)據(jù)，最后將該加密數(shù)據(jù)重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到四倍于共偏移距共方位角域道集的子共偏移距共方位角域道集，對(duì)所有的子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，便可得到壓制采集 腳印后的地震數(shù)據(jù)分析結(jié)果。本發(fā)明實(shí)施例的技術(shù)方案，利用共偏移距共方位角域道集的單次覆蓋特性，通過(guò)在同一偏移距和同一方位角內(nèi)進(jìn)行插值，實(shí)現(xiàn)了加密炮線和檢波線的目的，實(shí)現(xiàn)方案簡(jiǎn)單、容易實(shí)現(xiàn)，對(duì)壓制采集腳印、淺層成像和高保真處理有重要的意義，符合高保真處理要求，所取得的結(jié)果也完全滿足振幅、頻率、相位等屬性在油氣識(shí)別中的需要。

進(jìn)一步的，在上述實(shí)施例一提供的采集腳印壓制方法中，上述步驟102中的，根據(jù)每個(gè)共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，可具體包括以下三種情況：

情況一：通過(guò)改變共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，而不改變所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)。

具體的，可采用下述坐標(biāo)變換公式(1)得到新增加的炮點(diǎn)的坐標(biāo)：

其中，x_s和y_s分別為共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x和y分別為新增加的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，θ為共偏移距共方位角域道集中炮線的方位角，dl_s為共偏移距共方位角域道集中的炮線距。

情況二：通過(guò)改變共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，而不改變共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，計(jì)算出新增加的檢波點(diǎn)的坐標(biāo)。

詳細(xì)的，通過(guò)采用下述坐標(biāo)變換公式(2)得到新增加的檢波點(diǎn)的坐標(biāo)：

其中，x_r和y_r分別為共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x′和y′分別為新增加的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，α為共偏移距共方位角域道集中檢波線的方位角，dl_r為共偏移距共方位角域道集中的檢波線距。

情況三：通過(guò)改變共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)和 檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)和檢波點(diǎn)的坐標(biāo)。

具體的，在該情況下可具體通過(guò)情況一中的公式(1)求出新增加的炮點(diǎn)的坐標(biāo)，通過(guò)情況二中的公式(2)求出新增加的檢波點(diǎn)的坐標(biāo)。

進(jìn)一步的，在上述實(shí)施例一提供的采集腳印壓制方法中，上述步驟103，也即，對(duì)每個(gè)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)，具體包括：

利用五維插值法對(duì)每個(gè)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值處理，得到加密數(shù)據(jù)。

在地震資料數(shù)字處理中，為了保證采集到的地震數(shù)據(jù)函數(shù)在時(shí)間和空間上的連續(xù)性，往往需要對(duì)其進(jìn)行插值處理，根據(jù)插值原理，應(yīng)用較少的地震數(shù)據(jù)作為參數(shù)值，確定出任意時(shí)間和空間的函數(shù)值。由于同一共偏移距共方位角域道集內(nèi)數(shù)據(jù)具有相同的方位角，因此數(shù)據(jù)的一致性顯著好于共偏移距域道集，故數(shù)據(jù)插值的精度明顯高于常規(guī)的共偏移距域插值。

地震數(shù)據(jù)體可以實(shí)現(xiàn)五維采集，并且許多地震處理技術(shù)也可在五維空間內(nèi)實(shí)現(xiàn)。地震處理技術(shù)大都對(duì)地震數(shù)據(jù)空間分布的規(guī)則性有較嚴(yán)格的要求，但地震數(shù)據(jù)往往由于各種原因，空間采樣是不規(guī)則或缺失的，因而會(huì)造成地震數(shù)據(jù)處理結(jié)果較差。因此，本發(fā)明實(shí)施例采用不規(guī)則觀測(cè)系統(tǒng)五維插值算法對(duì)地震數(shù)據(jù)進(jìn)行插值重建，為地震數(shù)據(jù)處理提供高質(zhì)量的資料，提高了地震數(shù)據(jù)處理的效果，精度很高。

值得說(shuō)明的是，具體的插值方法可視處理系統(tǒng)的需要而定，本發(fā)明并不對(duì)此進(jìn)行限定。

下面以采集腳印壓制方法在實(shí)際中的應(yīng)用為例進(jìn)行說(shuō)明。

具體的，利用的地震數(shù)據(jù)資料具體為：炮線之間的距離為500米，檢波線之間的距離為400米，炮線和檢波線組成的十字排列域的覆蓋次數(shù)為180次，炮線的方位角為88度，檢波線的方位角為358度。

本發(fā)明該實(shí)施例的具體實(shí)現(xiàn)步驟如下：

第一：將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內(nèi)，得到180個(gè)共偏移距共方位角域道集；

在該實(shí)施例中，該十字排列域共被180個(gè)共偏移距共方位角域道集剖分，個(gè)數(shù)與由炮線和檢波線組成的十字排列域的覆蓋次數(shù)相同。

圖2為某一十字排列域及該十字排列域上的共偏移距共方位角域道集。如圖2所示，該十字排列域上共有180個(gè)不同顏色的方塊，其代表180個(gè)共偏移距共方位角域道集，表明相同偏移距和相同方位角區(qū)域被覆蓋的次數(shù)為180次。

第二：分別將炮線距和檢波線距進(jìn)行減半處理，根據(jù)每個(gè)共偏移距共方位角域道集中的地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，每個(gè)共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集；

圖3為從十字排列域中抽取的某一共偏移距共方位角域道集。如圖3所示，從上述十字排列域中抽取該共偏移距共方位角域道集在空間的展布具有單次覆蓋的特征。值得說(shuō)明的是，本發(fā)明實(shí)施例可任意抽取一個(gè)共偏移距共方位角域道集，只要其在空間的展布滿足單次覆蓋特征均可，本發(fā)明實(shí)施例并不對(duì)此進(jìn)行限定。

圖4為炮線距和檢波線距減半前形成的十字排列域；圖5為炮線距和檢波線距減半后形成的十字排列域。如圖4和圖5所示，十字排列域的橫線表示炮點(diǎn)形成的炮線，縱線表示檢波點(diǎn)形成的檢波線，從圖5可知，炮線距和檢波線距均縮小為原來(lái)的一半，由于原來(lái)的炮線之間的距離為500米，檢波線之間的距離為400米，所以，形成的新的炮線距為250米，新的檢線距為200米。

第三：利用五維插值法對(duì)每個(gè)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)；該加密數(shù)據(jù)包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)。

圖6為共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值前的數(shù)據(jù)分布示意圖；圖7為共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值后的數(shù)據(jù)分布示意圖。圖6中的空白部分為地震數(shù)據(jù)分選到共偏移距共方位角域道集后形成的空洞，也即，本步驟主要是利用五維插值法對(duì)每個(gè)子共偏移距共方位角域道集中的空洞進(jìn)行數(shù)據(jù)插值。從圖6和圖7可以看出，插值后得到的加密數(shù)據(jù)與原始地震數(shù)據(jù)的振幅、頻率、相位特性一致，包括面波的特征也很一致，這表明五維插值法的插值精度很高。

第四：將上述加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到4N個(gè)子共偏移距共方位角域道集；

由于插值后的數(shù)據(jù)加上原始數(shù)據(jù)稱為加密數(shù)據(jù)，對(duì)加密數(shù)據(jù)按新的炮線距和檢波線距重新劃分共偏移距共方位角域單元，再抽取共偏移距共方位角域道集，由于炮線距和檢波線距均縮小為了原來(lái)的一倍，所以，原來(lái)的一個(gè)共偏移距共方位角域單元一化為四。

圖8為原始地震數(shù)據(jù)對(duì)共偏移距共方位角域道集的劃分示意圖；圖9為加密數(shù)據(jù)對(duì)共偏移距共方位角域道集的劃分示意圖。圖10為原始地震數(shù)據(jù)在共偏移距共方位角域道集上形成的項(xiàng)鏈圖；圖11為加密數(shù)據(jù)在共偏移距共方位角域道集上形成的項(xiàng)鏈圖。

從圖8可知，圖8上共有180個(gè)完整的共偏移距共方位角域道集，個(gè)數(shù)與原始采集覆蓋次數(shù)相同，加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到720個(gè)子共偏移距共方位角域道集，也即，十字排列域形成的某一網(wǎng)格覆蓋次數(shù)增大為原來(lái)的4倍，具體參見(jiàn)圖9所示。類似的，圖10和圖11分別所示的原始地震數(shù)據(jù)和加密數(shù)據(jù)在共偏移距共方位角域道集上形成的項(xiàng)鏈圖也同樣說(shuō)明了上述結(jié)論。

圖12為原始地震數(shù)據(jù)在共偏移距共方位角域道集上的偏移距范圍。圖13為加密數(shù)據(jù)分別在相應(yīng)的四個(gè)子共偏移距共方位角域道集上的偏移距范圍。

參照?qǐng)D12和圖13，由于炮線距和檢波線距均縮小了一倍，因此，偏移距范圍也縮小了一倍，顯示在圖上后可知四個(gè)子共偏移距共方位角域道集的顏色更均一。所以，偏移距范圍越小，意味著偏移距分布的一致性越好，也更有利于淺層的成像和高保真處理。

第五：對(duì)上述720個(gè)子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

具體的，對(duì)上述720個(gè)子共偏移距共方位角域道集重新進(jìn)行偏移、疊加，再提取振幅切片，便可得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

圖14為原始地震數(shù)據(jù)的共偏移距偏移處理成果切片得到的地震數(shù)據(jù)分析結(jié)果，圖15為加密數(shù)據(jù)的共偏移距偏移處理成果切片得到的地震數(shù)據(jù)分析結(jié)果。從圖14可以看出，雖然對(duì)原始地震數(shù)據(jù)的處理已經(jīng)施加了常規(guī)壓制采集腳印的手段，但采集腳印還是明顯存在。從圖15可以看出，對(duì)加密數(shù)據(jù)處理后得到的地震數(shù)據(jù)分析結(jié)果顯示，采集腳印幾乎消失，效果非常好。

上述結(jié)果表明，本發(fā)明實(shí)施例利用共偏移距共方位角域道集的單次覆蓋特性，通過(guò)在同一方位角內(nèi)進(jìn)行五維插值實(shí)現(xiàn)了加密炮線和檢波線的目的，整個(gè)技術(shù)實(shí)現(xiàn)過(guò)程高效、快捷、方便，且對(duì)于壓制采集腳印、淺層成像、高保真處理都具有實(shí)用的價(jià)值。

圖16為本發(fā)明提供的采集腳印壓制系統(tǒng)實(shí)施例一的結(jié)構(gòu)示意圖。如圖16所示，本發(fā)明實(shí)施例一提供的采集腳印壓制系統(tǒng)，包括：

獲取模塊161，用于將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內(nèi)，得到N個(gè)共偏移距共方位角域道集；

其中，共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數(shù)，N為大于或等于1的整數(shù)。

處理模塊162，用于分別將炮線之間的距離和檢波線之間的距離進(jìn)行減半處理，根據(jù)每個(gè)共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，每個(gè)共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集；

插值模塊163，用于分別對(duì)四個(gè)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)；

加密數(shù)據(jù)包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)。

上述獲取模塊161，還用于將加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到4N個(gè)子共偏移距共方位角域道集；

上述處理模塊162，還用于對(duì)4N個(gè)子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

本發(fā)明實(shí)施例提供的采集腳印壓制系統(tǒng)，可以用于執(zhí)行如圖1所示采集腳印壓制方法實(shí)施例的技術(shù)方案，其實(shí)現(xiàn)原理和技術(shù)效果類似，此處不再贅述。

進(jìn)一步的，在上述實(shí)施例一提供的采集腳印壓制系統(tǒng)中，上述處理模塊162，包括：第一處理單元、第二處理單元和第三處理單元；

其中，第一處理單元，用于通過(guò)改變共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，而不改變共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn) 坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)；

具體的，該第一處理單元，具體用于采用下述坐標(biāo)變換公式(1)得到新增加的炮點(diǎn)的坐標(biāo)：

其中，x_s和y_s分別為共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x和y分別為新增加的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，θ為共偏移距共方位角域道集中炮線的方位角，dl_s為共偏移距共方位角域道集中的炮線距。

第二處理單元，用于通過(guò)改變共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，而不改變共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，計(jì)算出新增加的檢波點(diǎn)的坐標(biāo)；

具體的，該第二處理單元，具體用于采用下述坐標(biāo)變換公式(2)得到新增加的檢波點(diǎn)的坐標(biāo)：

其中，x_r和y_r分別為共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x′和y′分別為新增加的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，α為共偏移距共方位角域道集中檢波線的方位角，dl_r為共偏移距共方位角域道集中的炮線距。

第三處理單元，用于通過(guò)改變共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)和檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)和檢波點(diǎn)的坐標(biāo)。

進(jìn)一步的，在上述實(shí)施例一提供的采集腳印壓制系統(tǒng)中，上述插值模塊163，具體用于利用五維插值法對(duì)四個(gè)子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值處理，得到加密數(shù)據(jù)。

本發(fā)明實(shí)施例提供的采集腳印壓制系統(tǒng)，通過(guò)在炮線之間和檢波線之間內(nèi)插炮線和檢波線，使得組成十字排列域的炮線和檢波線加密，并利用共偏移距共方位角域道集的單次覆蓋特性，通過(guò)在同一方位角內(nèi)進(jìn)行五維高精度插值得到了加密數(shù)據(jù)，將該加密數(shù)據(jù)按照偏移距和方位角重新分選到加密炮線和加密檢波線組成的十字排列域內(nèi)，使獲得的子共偏移距共方位角域道集 為原來(lái)的四倍，整個(gè)技術(shù)實(shí)現(xiàn)過(guò)程高效、快捷、方便，最后對(duì)所有的子共偏移距共方位角域道集進(jìn)行處理，得到的地震數(shù)據(jù)分析結(jié)果，有效的壓制了采集腳印的產(chǎn)生，因此，該技術(shù)方案對(duì)壓制采集腳印、淺層成像、高保真處理均具有很高的實(shí)用價(jià)值。

最后應(yīng)說(shuō)明的是：以上各實(shí)施例僅用以說(shuō)明本發(fā)明的技術(shù)方案，而非對(duì)其限制；盡管參照前述各實(shí)施例對(duì)本發(fā)明進(jìn)行了詳細(xì)的說(shuō)明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解：其依然可以對(duì)前述各實(shí)施例所記載的技術(shù)方案進(jìn)行修改，或者對(duì)其中部分或者全部技術(shù)特征進(jìn)行等同替換；而這些修改或者替換，并不使相應(yīng)技術(shù)方案的本質(zhì)脫離本發(fā)明各實(shí)施例技術(shù)方案的范圍。