本申請涉及石油地震勘探中時移地震數(shù)據(jù)處理技術(shù)領(lǐng)域，特別涉及一種地震數(shù)據(jù)時差校正的方法及系統(tǒng)。

背景技術(shù)：

四維地震技勘探術(shù)，也稱時移地震勘探技術(shù)，是利用針對同一勘探區(qū)域的兩期或多期地震響應(yīng)的差異信息進行油藏動態(tài)特征研究的技術(shù)。該技術(shù)通過分析由于開采因素所造成的儲層內(nèi)孔隙壓力、孔隙度、流體成分、流體飽和度及溫度等油藏特性的變化，從而了解儲層中流體流動狀況以及流體成分的變化情況，為制定或更新油氣開發(fā)方案提供有效依據(jù)。如果兩期時移地震數(shù)據(jù)之間的地震響應(yīng)差異僅包括油藏變化產(chǎn)生的期望振幅差異，時移地震勘探技術(shù)即可通過兩期時移地震數(shù)據(jù)相減以消除由非儲層巖石的構(gòu)造、巖性等形成的共同地震響應(yīng)，進而從相減得到的期望振幅差異來研究儲層內(nèi)部流體性質(zhì)的變化。在時移地震勘探技術(shù)實施過程中，由于受到施工環(huán)境變化、采集儀器和采集參數(shù)差異等因素的影響，所述地震響應(yīng)差異往往既包括油藏變化產(chǎn)生的期望振幅差異，又包括由兩期地震數(shù)據(jù)時差等非油藏變化產(chǎn)生的非期望振幅差異。因此，需要消除兩期時移地震數(shù)據(jù)時差，以獲得無時差影響的期望振幅差異。

針對上述兩期時移地震數(shù)據(jù)時差，目前常用的一種時差校正方法是小時窗互相關(guān)方法，其方法是利用互相關(guān)方法求取基礎(chǔ)數(shù)據(jù)和監(jiān)測數(shù)據(jù)在相同小時窗(比如時窗的窗長為8毫秒至64毫秒的時窗)位置處的時間樣點的互相關(guān)極值，將求取的互相關(guān)極值作為所述小時窗中心點對應(yīng)的時間樣點的時差量，將小時窗移動一個時間樣點，計算移動后的小時窗中心點對應(yīng)的時間樣點的時差量，直至完成監(jiān)測地震數(shù)據(jù)中所有時間樣點的時差量的計算，利用計算的時差量對監(jiān)測地震數(shù)據(jù)進行時差校正。

時移地震數(shù)據(jù)時差產(chǎn)生的原因包括采集環(huán)境和采集參數(shù)變化引起的采集相關(guān)時差，也包括上覆儲層由于流體變化引入的局部時差。

發(fā)明人發(fā)現(xiàn)現(xiàn)有技術(shù)中至少存在如下問題：我國油田往往含多個薄油層，開發(fā)方式多采用多層聯(lián)采，上覆儲層流體變化將引起下層地震反射信號的到達(dá)時間發(fā)生變化，形成與儲層相關(guān)的兩期時移地震數(shù)據(jù)局部時差。目前這種上覆儲層流體變化引起的局部時差問題并沒有引起工業(yè)界重視，但這部分局部時差對薄油層的兩期時移地震數(shù)據(jù)的振幅差異影響較大。針對上述上覆儲層流體變化引起的局部時差問題，需要選取時窗的窗長較小的時窗計算時差量，但小時窗互相關(guān)方法中所采用的互相關(guān)方法隨時窗的窗長減小而計算精度降低。因此，小時窗互相關(guān)方法對上覆儲層流體變化引起的兩期時移地震數(shù)據(jù)局部時差的計算精度較低。

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

本申請實施例的目的是提供一種地震時差校正的方法及系統(tǒng)，以提高地震數(shù)據(jù)時差校正的精度。

為解決上述技術(shù)問題，本申請實施例提供一種地震數(shù)據(jù)時差校正的方法及系統(tǒng)是這樣實現(xiàn)的：

一種地震數(shù)據(jù)時差校正的方法，包括：

獲取目標(biāo)區(qū)域第一時間對應(yīng)的基礎(chǔ)地震數(shù)據(jù)，和第二時間對應(yīng)的監(jiān)測地震數(shù)據(jù)；

分別獲取基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)中相同位置對應(yīng)的第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)，分別選取第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點；

預(yù)設(shè)多個第一時差量，計算與每個所述第一時差量分別對應(yīng)的第一時窗位置處的時間樣點的第一方差值；

根據(jù)所述第一方差值確定最小方差值，將所述最小方差值對應(yīng)的第二時差量作為第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點對應(yīng)的目標(biāo)時差量；

根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到目標(biāo)地震數(shù)據(jù)。

優(yōu)選方案中，所述計算與每個所述第一時差量分別對應(yīng)的第一時窗位置處的時間樣點的第一方差值，包括：

獲取第一地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間對應(yīng)的第一地震道數(shù)據(jù)的第一振幅值；

根據(jù)第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間減去第一時差量，得到第二地震道數(shù)據(jù)中第一時窗位置處的第二采樣時間，獲取第二采樣時間對應(yīng)的第二地震道數(shù)據(jù)的第二振幅值；

根據(jù)第一振幅值和第二振幅值，計算第一時差量對應(yīng)的第一時窗位置處的時間樣點的第一方差值。

優(yōu)選方案中，所述根據(jù)第一振幅值和第二振幅值，計算第一時差量對應(yīng)的第一時窗位置處的時間樣點的第一方差值，采用下述公式實現(xiàn)：

公式中，第一振幅值_i表示第一地震道數(shù)據(jù)中第一時窗位置處的第i個時間樣點的第一振幅值，第二振幅值_i為第二地震道數(shù)據(jù)中第一時窗位置處的第i個時間樣點的第二振幅值。

優(yōu)選方案中，所述預(yù)設(shè)多個第一時差量，包括：

選取第一時差量的第一時差量范圍和第一時差增量；

根據(jù)第一時差量范圍和第一時差增量，得到多個第一時差量。

優(yōu)選方案中，第一時差增量為所述監(jiān)測地震數(shù)據(jù)的采樣間隔的四分之一。

優(yōu)選方案中，所述方法還包括：檢驗第一時窗的窗長、第一時差量范圍中最大第一時差量和最小第一時差量，以及第一時差增量的選取是否合適。

優(yōu)選方案中，所述檢驗第一時窗的窗長、第一時差量范圍中最大第一時差量和最小第一時差量，以及第一時差增量的選取是否合適，包括：

獲取所述監(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時差量，設(shè)定所述監(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時差量構(gòu)成的曲線的平滑閾值；

當(dāng)所述監(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時量差構(gòu)成的曲線的平滑度大于平滑閾值，第一時窗的窗長的選取不合適；

根據(jù)所述監(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正；

將所述基礎(chǔ)地震數(shù)據(jù)與所述監(jiān)測地震數(shù)據(jù)相減的地震數(shù)據(jù)在第一主測線或第一聯(lián)絡(luò)測線上的非儲層變化引起的校正前地震響應(yīng)，同所述基礎(chǔ)地震數(shù)據(jù)與所述經(jīng)過時差校正后的所述監(jiān)測地震數(shù)據(jù)相減的地震數(shù)據(jù)在第一主測線或第一聯(lián)絡(luò)測線上的非儲層變化引起的校正后地震響應(yīng)對比，設(shè)定所述校正后地震響應(yīng)對比所述校正前地震響應(yīng)的地震響應(yīng)減弱閾值；

當(dāng)所述校正后地震響應(yīng)對比所述校正前地震響應(yīng)的地震響應(yīng)減弱度小于所述地震響應(yīng)減弱閾值，第一時差量范圍中最大第一時差量、第一時差量范圍中最小第一時差量或第一時差增量的選取不合適。

優(yōu)選方案中，所述根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到目標(biāo)地震數(shù)據(jù)，包括：

根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到校正地震數(shù)據(jù)；

對所述校正地震數(shù)據(jù)進行辛格函數(shù)插值，得到目標(biāo)地震數(shù)據(jù)。

一種地震數(shù)據(jù)時差校正的系統(tǒng)，包括：地震數(shù)據(jù)獲取單元、時間樣點獲取單元、方差值計算單元、目標(biāo)時差量獲取單元和目標(biāo)地震數(shù)據(jù)獲取單元；其中，

所述地震數(shù)據(jù)獲取單元，用于獲取目標(biāo)區(qū)域第一時間對應(yīng)的基礎(chǔ)地震數(shù)據(jù)，和第二時間對應(yīng)的監(jiān)測地震數(shù)據(jù)；

所述時間樣點獲取單元，用于分別獲取基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)中相同位置對應(yīng)的第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)，分別選取第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點；

所述方差值計算單元，用于預(yù)設(shè)多個第一時差量，計算與每個所述第一時差量分別對應(yīng)的第一時窗位置處的時間樣點的第一方差值；

所述目標(biāo)時差量獲取單元，用于根據(jù)所述第一方差值確定最小方差值，將所述最小方差值對應(yīng)的第二時差量作為第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點對應(yīng)的目標(biāo)時差量；

所述目標(biāo)地震數(shù)據(jù)獲取單元，用于根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到目標(biāo)地震數(shù)據(jù)。

優(yōu)選方案中，所述方差值計算單元包括：第一振幅值獲取單元、第二振幅值獲取單元和第一方差值計算單元；其中，

所述第一振幅值獲取單元，用于獲取第一地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間對應(yīng)的第一地震道數(shù)據(jù)的第一振幅值；

所述第二振幅值獲取單元，用于根據(jù)第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間減去第一時差量，得到第二地震道數(shù)據(jù)中第一時窗位置處的第二采樣時間，獲取第二采樣時間對應(yīng)的第二地震道數(shù)據(jù)的第二振幅值；

所述第一方差值計算單元，用于根據(jù)第一振幅值和第二振幅值，計算第一時差量對應(yīng)的第一時窗位置處的時間樣點的第一方差值。

優(yōu)選方案中，所述系統(tǒng)還包括：檢驗單元；其中，

所述檢驗單元，用于檢驗第一時窗的窗長、第一時差量范圍中最大第一時差量和最小第一時差量，以及第一時差增量的選取是否合適。

優(yōu)選方案中，所述目標(biāo)地震數(shù)據(jù)獲取單元包括：校正地震數(shù)據(jù)獲取單元和函數(shù)插值單元；其中，

所述校正地震數(shù)據(jù)獲取單元，用于根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到校正地震數(shù)據(jù)；

所述函數(shù)插值單元，用于對所述校正地震數(shù)據(jù)進行辛格函數(shù)插值，得到目標(biāo)地震數(shù)據(jù)。

本申請?zhí)峁┝艘环N地震數(shù)據(jù)時差校正的方法及系統(tǒng)，所述方法基于小時窗方差的方式獲取并校正兩期地震時差，不僅能夠校正各類與采集相關(guān)的較大的時差量，也可以校正上覆儲層流體引起的局部較小的時差量，消除了非儲層變化引起的非期望振幅差異，顯著提高了兩期地震期望振幅差異精度。

附圖說明

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

圖1是本申請一種地震數(shù)據(jù)時差校正的方法實施例的流程圖；

圖2是本申請實施例中Marmousi2模型基礎(chǔ)地震數(shù)據(jù)的振幅剖面、監(jiān)測地震數(shù)據(jù)的振幅剖面以及基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)的振幅差剖面的示意圖；

圖3是本申請實施例中Marmousi2模型基礎(chǔ)地震數(shù)據(jù)與監(jiān)測地震數(shù)據(jù)之間的目標(biāo)時差量剖面以及基礎(chǔ)地震數(shù)據(jù)與目標(biāo)地震數(shù)據(jù)的振幅差剖面的示意圖；

圖4是本申請實施例中基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)的振幅差在目的儲層頂?shù)难貙忧衅@示的示意圖；

圖5是本申請實施例中基礎(chǔ)地震數(shù)據(jù)和目標(biāo)地震數(shù)據(jù)的振幅差在目的儲層頂?shù)难貙忧衅@示的示意圖；

圖6是本申請一種地震數(shù)據(jù)時差校正的系統(tǒng)實施例的組成結(jié)構(gòu)圖；

圖7是本申請地震數(shù)據(jù)時差校正的系統(tǒng)實施例中方差值計算單元的組成結(jié)構(gòu)圖；

圖8是本申請地震數(shù)據(jù)時差校正的系統(tǒng)實施例中目標(biāo)地震數(shù)據(jù)獲取單元的組成結(jié)構(gòu)圖。

具體實施方式

本申請實施例提供一種地震數(shù)據(jù)時差校正的方法及系統(tǒng)。

為了使本技術(shù)領(lǐng)域的人員更好地理解本申請中的技術(shù)方案，下面將結(jié)合本申請實施例中的附圖，對本申請實施例中的技術(shù)方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本申請一部分實施例，而不是全部的實施例?；诒旧暾堉械膶嵤├绢I(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都應(yīng)當(dāng)屬于本申請保護的范圍。

圖1是本申請一種地震數(shù)據(jù)時差校正的方法實施例的流程圖。如圖1所示，所述地震數(shù)據(jù)時差校正的方法，包括以下步驟。

步驟S101：獲取目標(biāo)區(qū)域第一時間對應(yīng)的基礎(chǔ)地震數(shù)據(jù)，和第二時間對應(yīng)的監(jiān)測地震數(shù)據(jù)。

具體地，可以采用Marmousi2模型模擬目標(biāo)區(qū)域第一時間對應(yīng)的基礎(chǔ)地震數(shù)據(jù)的振幅剖面，和第二時間對應(yīng)的監(jiān)測地震數(shù)據(jù)的振幅剖面。其中Marmousi2模型是基于Marmousi模型增加了一些儲層并改進后的模型；Marmousi模型最早是由法國石油研究院(IFP)所屬協(xié)會根據(jù)西非海上油田制作的，這種模型及其聲波有限差分合成數(shù)據(jù)被應(yīng)用于許多地球物理科研項目，例如，石油地震勘探數(shù)據(jù)處理等科研項目；Marmousi模型處理地震數(shù)據(jù)的過程是地質(zhì)人員根據(jù)實際地震數(shù)據(jù)畫出地質(zhì)模型，如泥巖地層，砂巖地層，斷層等，再根據(jù)測井速度對地質(zhì)模型填上速度，得到速度模型，對速度模型做正演模擬，獲得了合成的地震數(shù)據(jù)。第一時間對應(yīng)的基礎(chǔ)地震數(shù)據(jù)和第二時間對應(yīng)的監(jiān)測地震數(shù)據(jù)分別為同一目標(biāo)區(qū)域在不同時間的地震數(shù)據(jù)。這樣可以為后續(xù)時差校正的步驟做好地震數(shù)據(jù)的準(zhǔn)備。

例如，圖2是本申請實施例中Marmousi2模型基礎(chǔ)地震數(shù)據(jù)的振幅剖面、監(jiān)測地震數(shù)據(jù)的振幅剖面以及基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)的振幅差剖面的示意圖。圖2中(a)圖為基礎(chǔ)地震數(shù)據(jù)的振幅剖面圖，圖2中(b)圖為監(jiān)測地震數(shù)據(jù)的振幅剖面圖，圖2中(c)圖為基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)的振幅差剖面圖。圖2中(a)圖和(b)圖的箭頭所示為實際流體發(fā)生變化的儲層地震數(shù)據(jù)的振幅剖面。圖2中(c)圖中基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)的振幅差可以包括：采集參數(shù)和采集環(huán)境引起的采集相關(guān)時差導(dǎo)致的振幅差，以及上覆儲層流體變化引起的局部時差導(dǎo)致的振幅差?？梢栽谒霰O(jiān)測地震數(shù)據(jù)中設(shè)置采集參數(shù)和采集環(huán)境引起的采集相關(guān)時差為2毫秒(ms)。

步驟S102：分別獲取基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)中相同位置對應(yīng)的第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)，分別選取第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點。

具體地，根據(jù)所述基礎(chǔ)地震數(shù)據(jù)和所述監(jiān)測地震數(shù)據(jù)，可以分別獲取所述基礎(chǔ)地震數(shù)據(jù)和所述監(jiān)測地震數(shù)據(jù)中相同位置對應(yīng)的第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)。可以分別選取第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點。其中，第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)可以分別為所述基礎(chǔ)地震數(shù)據(jù)和所述監(jiān)測地震數(shù)據(jù)中任意相同位置對應(yīng)的地震道數(shù)據(jù)。第一時窗位置可以為第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)中任意時窗位置，其中任意時窗位置對應(yīng)的時窗的窗長可以為32ms。這里選取第一時窗位置處的時間樣點為后續(xù)計算第一時窗位置處的時間樣點的時差做準(zhǔn)備，同時也體現(xiàn)了后續(xù)計算的第一時窗位置處的時間樣點的時差是隨時間、空間變化的。

步驟S103：預(yù)設(shè)多個第一時差量，計算與每個所述第一時差量分別對應(yīng)的第一時窗位置處的時間樣點的第一方差值。

具體地，可以選取第一時差量的第一時差量范圍為-4ms～4ms?？梢赃x取第一時差量的第一時差增量為所述監(jiān)測地震數(shù)據(jù)的采樣間隔的四分之一，即為0.25ms；其中，所述監(jiān)測地震數(shù)據(jù)的采樣間隔時間可以為1ms。根據(jù)第一時差量范圍和第一時差增量，可以得到多個第一時差量，即為-4ms、-3.75ms、-3.5ms、…0ms、…3.5ms、3.75ms、4ms。其中省略號部分為按照所述時差增量依次增加后的第一時差量。

根據(jù)所述選取的第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點，可以獲取第一地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間對應(yīng)的第一地震道數(shù)據(jù)的第一振幅值。其中第一采樣時間為第一地震道數(shù)據(jù)中第一時窗位置處的任意時間樣點對應(yīng)的采樣時間。根據(jù)第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間減去第一時差量，可以得到第二地震道數(shù)據(jù)中的第二采樣時間?？梢垣@取第二采樣時間對應(yīng)的第二地震道數(shù)據(jù)中的時間樣點的第二振幅值。第一振幅值和第二振幅值分別為第一地震道數(shù)據(jù)中第一采樣時間和第二地震道數(shù)據(jù)中第二采樣時間對應(yīng)的振幅值。

根據(jù)第一振幅值和第二振幅值，可以采用下述公式計算第一時差量對應(yīng)的第一時窗位置處的時間樣點的第一方差值：

公式中，第一振幅值_i表示第一地震道數(shù)據(jù)中第一時窗位置處的第i個時間樣點的第一振幅值，第二振幅值_i為第二地震道數(shù)據(jù)中第一時窗位置處的第i個時間樣點的第二振幅值。第一方差值可以為第一時差量范圍-4ms～4ms內(nèi)不同時差量對應(yīng)的所述計算的多個方差值。這個步驟可以通過所述計算的方差值體現(xiàn)不同時差量的真實性，從而為后續(xù)步驟中確定準(zhǔn)確的時差量做準(zhǔn)備。

步驟S104：根據(jù)所述第一方差值確定最小方差值，將所述最小方差值對應(yīng)的第二時差量作為第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點對應(yīng)的目標(biāo)時差量。

具體地，根據(jù)所述計算的每一個第一時差量對應(yīng)的第一時窗位置處的時間樣點的第一方差值，可以比較所述第一方差值的大小，可以確定最小方差值。最小方差值可以表明其對應(yīng)的時差量最能真實體現(xiàn)第二地震道數(shù)據(jù)中第一時窗的中間位置處的時間樣點的時差?？梢詫⑺鲎钚》讲钪祵?yīng)的第二時差量作為第二地震道數(shù)據(jù)中第一時窗的中間位置處的時間樣點對應(yīng)的目標(biāo)時差量。其中，選取第一時窗位置處的時間樣點個數(shù)為奇數(shù)。第二時差量為最小方差值對應(yīng)的時差量。

可以將第一時窗的位置移動一個時間樣點。采用上述步驟可以計算移動后的第一時窗的中間位置處的時間樣點對應(yīng)的目標(biāo)時差量，直至完成所述監(jiān)測地震數(shù)據(jù)中第二地震道中每一個時間樣點對應(yīng)的目標(biāo)時差量。從而可以完成所述監(jiān)測地震數(shù)據(jù)中每一個地震道中每一個時間樣點對應(yīng)的目標(biāo)時差量。不僅能夠得到各類與采集相關(guān)的較大的時差量，也可以得到上覆儲層流體引起的局部較小的時差量。

例如，圖3中(a)圖為圖2中(a)圖中的基礎(chǔ)地震數(shù)據(jù)和圖2中(b)圖中的監(jiān)測地震數(shù)據(jù)之間的目標(biāo)時差量剖面圖。如圖3中(a)圖所示，可以看出本發(fā)明不僅可以準(zhǔn)確計算人為設(shè)置的2ms系統(tǒng)時差，也可以準(zhǔn)確得到上覆儲層流體變化引入的局部時差。

在另一個實施方式中，還可以檢驗第一時窗的窗長、第一時差量范圍中最大第一時差量和最小第一時差量，以及第一時差增量的選取是否合適。

具體地，可以獲取所述監(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時差量?？梢詾樗霰O(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時差量構(gòu)成的曲線設(shè)定一平滑閾值。當(dāng)所述監(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時量差構(gòu)成的曲線的平滑度大于平滑閾值時，第一時窗的窗長的選取可能不合適，可以重新選擇所述監(jiān)測地震數(shù)據(jù)中第二時窗位置的第二時窗的窗長。所述平滑閾值的取值范圍可以為大于或等于2毫秒。第一主測線和第一聯(lián)絡(luò)測線可以分別為所述監(jiān)測地震數(shù)據(jù)中任意主測線和任意聯(lián)絡(luò)測線。所述監(jiān)測地震數(shù)據(jù)可以包括三個維度，分別為：主測線維度、聯(lián)絡(luò)測線維度和時間維度；其中，所述時間維度、主測線維度和所述聯(lián)絡(luò)測線維度兩兩相互垂直。

根據(jù)所述監(jiān)測地震數(shù)據(jù)中第一主測線或第一聯(lián)絡(luò)測線上的時間樣點的目標(biāo)時差量可以對所述監(jiān)測地震數(shù)據(jù)進行時差校正。將所述基礎(chǔ)地震數(shù)據(jù)與所述監(jiān)測地震數(shù)據(jù)相減的地震數(shù)據(jù)在第一主測線或第一聯(lián)絡(luò)測線上的非儲層變化引起的校正前地震響應(yīng)，可以同所述基礎(chǔ)地震數(shù)據(jù)與所述經(jīng)過時差校正后的監(jiān)測地震數(shù)據(jù)相減的地震數(shù)據(jù)在第一主測線或第一聯(lián)絡(luò)測線上的非儲層變化引起的校正后地震響應(yīng)對比?？梢詾樗鲂Ｕ蟮卣痦憫?yīng)對比所述校正前地震響應(yīng)的減弱程度設(shè)定一地震響應(yīng)減弱閾值。當(dāng)所述校正后地震響應(yīng)對比所述校正前地震響應(yīng)的地震響應(yīng)減弱度小于所述地震響應(yīng)減弱閾值時，第一時差量范圍中最大第一時差量、第一時差量范圍中最小第一時差量或第一時差增量的選取可能不合適，可以重新選擇所述監(jiān)測地震數(shù)據(jù)中第二時差量范圍中最大第二時差量、第二時差量范圍中最小第二時差量或第二時差增量。所述地震響應(yīng)減弱閾值的取值范圍可以為大于或等于90％。

步驟S105：根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到目標(biāo)地震數(shù)據(jù)。

具體地，根據(jù)所述監(jiān)測地震數(shù)據(jù)中每一個地震道中每一個時間樣點對應(yīng)的目標(biāo)時差量，可以對所述監(jiān)測地震數(shù)據(jù)進行時差校正?？梢缘玫叫Ｕ卣饠?shù)據(jù)。可以對所述校正地震數(shù)據(jù)進行辛格函數(shù)(Sinc)插值，可以得到目標(biāo)地震數(shù)據(jù)。Sinc函數(shù)是一個簡單的最小相位的子波序列；由于所述時差校正后的校正地震數(shù)據(jù)采樣間隔小于所述監(jiān)測地震數(shù)據(jù)采樣間隔且采樣時間序列可能不連續(xù)，可以采用Sinc函數(shù)的八點截斷實現(xiàn)所述校正地震數(shù)據(jù)的垂向插值。所述Sinc函數(shù)可以采用下述公式表示：

$Sinc\left(t\right)=\frac{sin\left(t\right)}{t}$

公式中，t表示為所述校正地震數(shù)據(jù)中的采樣時間，Sinc(t)表示為在所述校正地震數(shù)據(jù)中的采樣時間t處插入的振幅值。

例如，圖3中(b)圖為圖2中(a)圖中的基礎(chǔ)地震數(shù)據(jù)和2中(b)圖中的目標(biāo)地震數(shù)據(jù)的振幅差剖面圖，圖3中(b)圖的箭頭所示為實際儲層變化的期望振幅差異剖面?？梢耘c圖2中(c)圖中的基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)的振幅差比較。如圖3中(b)圖所示，可以看出由所述采集腳印相關(guān)時差和上覆儲層流體引起的時差導(dǎo)致的非儲層振幅差異可以得到消除，實際儲層變化的期望振幅差異可以得到凸顯。

圖4是本申請實施例中基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)的振幅差在目的儲層頂?shù)难貙忧衅@示的示意圖。圖5是本申請實施例中基礎(chǔ)地震數(shù)據(jù)和目標(biāo)地震數(shù)據(jù)的振幅差在目的儲層頂?shù)难貙忧衅@示的示意圖。圖4和圖5的橢圓形圈內(nèi)區(qū)域部分包含了井眼處實際儲層變化引起的振幅差。如圖4和圖5所示，可以看出由所述采集腳印相關(guān)時差和上覆儲層流體引起的時差導(dǎo)致的非儲層振幅差異可以得到有效消除。橢圓形圈內(nèi)區(qū)域部分所示的井眼處實際儲層變化引起的振幅差異更為突出。

所述地震數(shù)據(jù)時差校正的方法實施例，基于小時窗方差的方式獲取并校正基礎(chǔ)地震數(shù)據(jù)和監(jiān)督地震數(shù)據(jù)之間時差量，不僅能夠校正各類與采集相關(guān)的較大的時差量，也可以校正上覆儲層流體引起的局部較小的時差量，消除了非儲層變化引起的非期望振幅差異，顯著提高了兩期地震期望振幅差異精度。

圖6是本申請一種地震數(shù)據(jù)時差校正的系統(tǒng)實施例的組成結(jié)構(gòu)圖。如圖6所示，所述地震數(shù)據(jù)時差校正的系統(tǒng)可以包括：地震數(shù)據(jù)獲取單元100、時間樣點獲取單元200、方差值計算單元300、目標(biāo)時差量獲取單元400和目標(biāo)地震數(shù)據(jù)獲取單元500。

所述地震數(shù)據(jù)獲取單元100，可以用于獲取目標(biāo)區(qū)域第一時間對應(yīng)的基礎(chǔ)地震數(shù)據(jù)，和第二時間對應(yīng)的監(jiān)測地震數(shù)據(jù)。

所述時間樣點獲取單元200，可以用于分別獲取基礎(chǔ)地震數(shù)據(jù)和監(jiān)測地震數(shù)據(jù)中相同位置對應(yīng)的第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)?？梢苑謩e選取第一地震道數(shù)據(jù)和第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點。

所述方差值計算單元300，可以用于預(yù)設(shè)多個第一時差量，可以計算與每個所述第一時差量分別對應(yīng)的第一時窗位置處的時間樣點的第一方差值。

所述目標(biāo)時差量獲取單元400，可以用于根據(jù)所述第一方差值確定最小方差值，將所述最小方差值對應(yīng)的第二時差量作為第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點對應(yīng)的目標(biāo)時差量。

所述目標(biāo)地震數(shù)據(jù)獲取單元500，可以用于根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到目標(biāo)地震數(shù)據(jù)。

在另一種實施方式中，所述地震數(shù)據(jù)時差校正的系統(tǒng)還可以包括：檢驗單元。

所述檢驗單元，可以用于檢驗第一時窗的窗長、第一時差量范圍中最大第一時差量和最小第一時差量，以及第一時差增量的選取是否合適。

圖7是本申請地震數(shù)據(jù)時差校正的系統(tǒng)實施例中方差值計算單元的組成結(jié)構(gòu)圖。如圖7所示，圖6中方差值計算單元300包括：第一振幅值獲取單元310、第二振幅值獲取單元320和第一方差值計算單元330。

所述第一振幅值獲取單元310，可以用于獲取第一地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間對應(yīng)的第一地震道數(shù)據(jù)的第一振幅值。

所述第二振幅值獲取單元320，可以用于根據(jù)第二地震道數(shù)據(jù)中第一時窗位置處的時間樣點的第一采樣時間減去第一時差量，得到第二地震道數(shù)據(jù)中第一時窗位置處的第二采樣時間。可以獲取第二采樣時間對應(yīng)的第二地震道數(shù)據(jù)的第二振幅值。

所述第一方差值計算單元330，可以用于根據(jù)第一振幅值和第二振幅值，可以計算第一時差量對應(yīng)的第一時窗位置處的時間樣點的第一方差值。

圖8是本申請地震數(shù)據(jù)時差校正的系統(tǒng)實施例中目標(biāo)地震數(shù)據(jù)獲取單元的組成結(jié)構(gòu)圖。如圖8所示，圖6中目標(biāo)地震數(shù)據(jù)獲取單元500包括：校正地震數(shù)據(jù)獲取單元510和函數(shù)插值單元520。

所述校正地震數(shù)據(jù)獲取單元510，可以用于根據(jù)所述目標(biāo)時差量對所述監(jiān)測地震數(shù)據(jù)進行時差校正，得到校正地震數(shù)據(jù)。

所述函數(shù)插值單元520，可以用于對所述校正地震數(shù)據(jù)進行辛格函數(shù)插值，得到目標(biāo)地震數(shù)據(jù)。

所述地震數(shù)據(jù)時差校正的系統(tǒng)實施例與所述地震數(shù)據(jù)時差校正的方法實施例相對應(yīng)，可以實現(xiàn)校正各類與采集相關(guān)的較大的時差量，也可以實現(xiàn)校正上覆儲層流體引起的局部較小的時差量，消除了非儲層變化引起的非期望振幅差異，顯著提高了兩期地震期望振幅差異精度。

在20世紀(jì)90年代，對于一個技術(shù)的改進可以很明顯地區(qū)分是硬件上的改進(例如，對二極管、晶體管、開關(guān)等電路結(jié)構(gòu)的改進)還是軟件上的改進(對于方法流程的改進)。然而，隨著技術(shù)的發(fā)展，當(dāng)今的很多方法流程的改進已經(jīng)可以視為硬件電路結(jié)構(gòu)的直接改進。設(shè)計人員幾乎都通過將改進的方法流程編程到硬件電路中來得到相應(yīng)的硬件電路結(jié)構(gòu)。因此，不能說一個方法流程的改進就不能用硬件實體模塊來實現(xiàn)。例如，可編程邏輯器件(Programmable Logic Device,PLD)(例如現(xiàn)場可編程門陣列(Field Programmable Gate Array，F(xiàn)PGA))就是這樣一種集成電路，其邏輯功能由用戶對器件編程來確定。由設(shè)計人員自行編程來把一個數(shù)字系統(tǒng)“集成”在一片PLD上，而不需要請芯片制造廠商來設(shè)計和制作專用的集成電路芯片2。而且，如今，取代手工地制作集成電路芯片，這種編程也多半改用“邏輯編譯器(logic compiler)”軟件來實現(xiàn)，它與程序開發(fā)撰寫時所用的軟件編譯器相類似，而要編譯之前的原始代碼也得用特定的編程語言來撰寫，此稱之為硬件描述語言(Hardware Description Language，HDL)，而HDL也并非僅有一種，而是有許多種，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等，目前最普遍使用的是VHDL(Very-High-Speed Integrated Circuit Hardware Description Language)與Verilog2。本領(lǐng)域技術(shù)人員也應(yīng)該清楚，只需要將方法流程用上述幾種硬件描述語言稍作邏輯編程并編程到集成電路中，就可以很容易得到實現(xiàn)該邏輯方法流程的硬件電路。

控制器可以按任何適當(dāng)?shù)姆绞綄崿F(xiàn)，例如，控制器可以采取例如微處理器或處理器以及存儲可由該(微)處理器執(zhí)行的計算機可讀程序代碼(例如軟件或固件)的計算機可讀介質(zhì)、邏輯門、開關(guān)、專用集成電路(Application Specific Integrated Circuit，ASIC)、可編程邏輯控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存儲器控制器還可以被實現(xiàn)為存儲器的控制邏輯的一部分。

本領(lǐng)域技術(shù)人員也知道，除了以純計算機可讀程序代碼方式實現(xiàn)控制器以外，完全可以通過將方法步驟進行邏輯編程來使得控制器以邏輯門、開關(guān)、專用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來實現(xiàn)相同功能。因此這種控制器可以被認(rèn)為是一種硬件部件，而對其內(nèi)包括的用于實現(xiàn)各種功能的裝置也可以視為硬件部件內(nèi)的結(jié)構(gòu)?；蛘呱踔?，可以將用于實現(xiàn)各種功能的裝置視為既可以是實現(xiàn)方法的軟件模塊又可以是硬件部件內(nèi)的結(jié)構(gòu)。

上述實施例闡明的系統(tǒng)、裝置、模塊或單元，具體可以由計算機芯片或?qū)嶓w實現(xiàn)，或者由具有某種功能的產(chǎn)品來實現(xiàn)。

為了描述的方便，描述以上裝置時以功能分為各種單元分別描述。當(dāng)然，在實施本申請時可以把各單元的功能在同一個或多個軟件和/或硬件中實現(xiàn)。

通過以上的實施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請可借助軟件加必需的通用硬件平臺的方式來實現(xiàn)?；谶@樣的理解，本申請的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，在一個典型的配置中，計算設(shè)備包括一個或多個處理器(CPU)、輸入/輸出接口、網(wǎng)絡(luò)接口和內(nèi)存。該計算機軟件產(chǎn)品可以包括若干指令用以使得一臺計算機設(shè)備(可以是個人計算機，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請各個實施例或者實施例的某些部分所述的方法。該計算機軟件產(chǎn)品可以存儲在內(nèi)存中，內(nèi)存可能包括計算機可讀介質(zhì)中的非永久性存儲器，隨機存取存儲器(RAM)和/或非易失性內(nèi)存等形式，如只讀存儲器(ROM)或閃存(flash RAM)。內(nèi)存是計算機可讀介質(zhì)的示例。計算機可讀介質(zhì)包括永久性和非永久性、可移動和非可移動媒體可以由任何方法或技術(shù)來實現(xiàn)信息存儲。信息可以是計算機可讀指令、數(shù)據(jù)結(jié)構(gòu)、程序的模塊或其他數(shù)據(jù)。計算機的存儲介質(zhì)的例子包括，但不限于相變內(nèi)存(PRAM)、靜態(tài)隨機存取存儲器(SRAM)、動態(tài)隨機存取存儲器(DRAM)、其他類型的隨機存取存儲器(RAM)、只讀存儲器(ROM)、電可擦除可編程只讀存儲器(EEPROM)、快閃記憶體或其他內(nèi)存技術(shù)、只讀光盤只讀存儲器(CD-ROM)、數(shù)字多功能光盤(DVD)或其他光學(xué)存儲、磁盒式磁帶，磁帶磁磁盤存儲或其他磁性存儲設(shè)備或任何其他非傳輸介質(zhì)，可用于存儲可以被計算設(shè)備訪問的信息。按照本文中的界定，計算機可讀介質(zhì)不包括短暫電腦可讀媒體(transitory media)，如調(diào)制的數(shù)據(jù)信號和載波。

本說明書中的各個實施例均采用遞進的方式描述，各個實施例之間相同相似的部分互相參見即可，每個實施例重點說明的都是與其他實施例的不同之處。尤其，對于系統(tǒng)實施例而言，由于其基本相似于方法實施例，所以描述的比較簡單，相關(guān)之處參見方法實施例的部分說明即可。

本申請可用于眾多通用或?qū)Ｓ玫挠嬎銠C系統(tǒng)環(huán)境或配置中。例如：個人計算機、服務(wù)器計算機、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的消費電子設(shè)備、網(wǎng)絡(luò)PC、小型計算機、大型計算機、包括以上任何系統(tǒng)或設(shè)備的分布式計算環(huán)境等等。

本申請可以在由計算機執(zhí)行的計算機可執(zhí)行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類型的例程、程序、對象、組件、數(shù)據(jù)結(jié)構(gòu)等等。也可以在分布式計算環(huán)境中實踐本申請，在這些分布式計算環(huán)境中，由通過通信網(wǎng)絡(luò)而被連接的遠(yuǎn)程處理設(shè)備來執(zhí)行任務(wù)。在分布式計算環(huán)境中，程序模塊可以位于包括存儲設(shè)備在內(nèi)的本地和遠(yuǎn)程計算機存儲介質(zhì)中。

雖然通過實施例描繪了本申請，本領(lǐng)域普通技術(shù)人員知道，本申請有許多變形和變化而不脫離本申請的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請的精神。