專利名稱：一種陡傾角地層及裂縫地震偏移方法和裝置的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及針對陡傾角地層及裂縫的地震數(shù)據(jù)疊前深度偏移成像方法，具體來說就是一種利用雙反射疊前偏移方法、實現(xiàn)針對陡傾角地層及裂縫成像的地震偏移方法和裝置。
背景技術(shù)：
長期以來，陡傾角地層及裂縫地震成像難題都是制約復(fù)雜構(gòu)造區(qū)地震勘探進程的瓶頸，尤其是裂縫成像對于天然氣勘探具有重要意義，但是傳統(tǒng)疊前深度偏移方法在解決陡傾角地層及裂縫地震成像問題時存在很多困難。傳統(tǒng)疊前深度偏移方法一般是利用一次反射波進行成像，而將其它類型的地震波視為噪音，所謂的一次反射波如圖Ia所示，其波場傳播路徑為炮點10-地下反射層12-檢波點14，如果要利用一次反射波對陡傾角地層或裂縫進行成像，就必須直接記錄到來自這些陡傾角地層或裂縫的反射波，但是這就要求地震野外采集數(shù)據(jù)的孔徑非常大，而采用常規(guī)地震勘探方法，極少有可能記錄到這些反射波。而圖Ib和圖Ic所示為雙反射波，根據(jù)兩次反射發(fā)生次序的不同，雙反射波可以分為兩大類HV型雙反射波(圖Ib所示)和VH型雙反射波(圖Ic所示)。與傳統(tǒng)的一次反射波不同，對于HV型雙反射波(圖Ib所示)而言， 其傳播路徑為炮點10-地下水平界面16-陡傾角地層或裂縫18-檢波點14 ；對于VH型雙反射波(圖Ic)而言，其傳播路徑為炮點10-陡傾角地層或裂縫18-地下水平界面16-檢波點14，基于雙反射波的上述傳播特性，能夠?qū)崿F(xiàn)針對陡傾角地層及裂縫的偏移成像，而正演模擬研究也表明了該類地震波的存在，同時，利用雙反射波對陡傾角地層及裂縫進行偏移成像時，地震采集數(shù)據(jù)的孔徑可以相當(dāng)小，一般的野外采集設(shè)計均能滿足要求。而傳統(tǒng)疊前深度偏移方法在面對陡傾角地層及裂縫地震成像問題時存在以下困難首先，傳統(tǒng)的KirchhofT積分法和單程波疊前深度偏移方法由于自身算法的限制，僅能夠?qū)σ淮畏瓷洳ㄟM行成像，無法實現(xiàn)與陡傾角地層及裂縫相關(guān)的雙反射波的偏移成像；其次，逆時偏移雖然是基于對雙程聲波方程的直接求解，能夠?qū)Πp反射波在內(nèi)的不同類型的地震波均能準(zhǔn)確成像、且無傾角限制，但是在面對實際地震資料陡傾角裂縫偏移成像問題時，仍然遇到很多困難由于逆時偏移為全波場成像，不壓制一次反射波的能量，而一次波的能量要比雙反射波強得多，從而會掩蓋與陡傾角裂縫相關(guān)的雙反射波成像能量；同時，逆時偏移需要非常大的計算量，并且對速度模型精度要求極高，但逆時偏移本身無法提供有效的速度建模方法。

發(fā)明內(nèi)容
本發(fā)明實施例是提供一種具有較高計算效率、專門針對陡傾角地層及裂縫成像的地震數(shù)據(jù)疊前偏移方法和裝置。
為達上述目的，一方面，本發(fā)明實施例提供了一種陡傾角地層及裂縫地震偏移方法，所述方法包括輸入地震數(shù)據(jù)，所述地震數(shù)據(jù)包括通過預(yù)處理得到的共炮點道集地震數(shù)據(jù)、深度域偏移速度模型以及常規(guī)疊前深度偏移結(jié)果；利用所述常規(guī)疊前深度偏移結(jié)果，在深度域定義一個基礎(chǔ)地下水平反射層；對輸入的所述共炮點道集地震數(shù)據(jù)進行重新排列，生成共檢波點道集地震數(shù)據(jù)， 所述共炮點道集地震數(shù)據(jù)用于HV型雙反射波的疊前偏移，所述共檢波點道集地震數(shù)據(jù)用于對VH型雙反射波的疊前偏移；利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共炮點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共炮點道集地震數(shù)據(jù)；利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算HV型雙反射波偏移中需要的第一旅行時信息；根據(jù)所述延拓后的共炮點道集地震數(shù)據(jù)和所述第一旅行時信息，并利用成像條件對HV型雙反射波以步長Δ ζ進行逐層成像，得到HV型雙反射波的偏移結(jié)果；利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共檢波點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共檢波點道集地震數(shù)據(jù)；利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算VH型雙反射波偏移中需要的第二旅行時信息；根據(jù)所述延拓后的共檢波點道集地震數(shù)據(jù)和所述第二旅行時信息，并利用成像條件對VH型雙反射波以步長Δ ζ進行逐層成像，得到VH型雙反射波的偏移結(jié)果；對所述HV型雙反射波的偏移結(jié)果和所述VH型雙反射波的偏移結(jié)果進行疊加成像，獲得目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果。另一方面，為達上述目的，本發(fā)明實施例提供了一種陡傾角地層及裂縫地震偏移裝置，所述裝置包括地震數(shù)據(jù)輸入單元，用于輸入地震數(shù)據(jù)，所述地震數(shù)據(jù)包括通過預(yù)處理得到的共炮點道集地震數(shù)據(jù)、深度域偏移速度模型以及常規(guī)疊前深度偏移結(jié)果；基礎(chǔ)地下水平反射層定義單元，用于利用所述常規(guī)疊前深度偏移結(jié)果，在深度域定義一個基礎(chǔ)地下水平反射層；共檢波點道集地震數(shù)據(jù)生成單元，用于對輸入的所述共炮點道集地震數(shù)據(jù)進行重新排列，生成共檢波點道集地震數(shù)據(jù)，所述共炮點道集地震數(shù)據(jù)用于HV型雙反射波的疊前偏移，所述共檢波點道集地震數(shù)據(jù)用于對VH型雙反射波的疊前偏移；共炮點道集地震數(shù)據(jù)延拓處理單元，用于利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共炮點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共炮點道集地震數(shù)據(jù)；第一旅行時信息計算單元，用于利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算HV型雙反射波偏移中需要的第一旅行時信息；HV型雙反射波的偏移結(jié)果生成單元，用于根據(jù)所述延拓后的共炮點道集地震數(shù)據(jù)和所述第一旅行時信息，并利用成像條件對HV型雙反射波以步長Δ ζ進行逐層成像，得到HV型雙反射波的偏移結(jié)果；共檢波點道集地震數(shù)據(jù)延拓處理單元，用于利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共檢波點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共檢波點道集地震數(shù)據(jù)；第二旅行時信息計算單元，用于利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算VH型雙反射波偏移中需要的第二旅行時信息；VH型雙反射波的偏移結(jié)果生成單元，用于根據(jù)所述延拓后的共檢波點道集地震數(shù)據(jù)和所述第二旅行時信息，并利用成像條件對VH型雙反射波以步長Δζ進行逐層成像，得到VH型雙反射波的偏移結(jié)果；陡傾角地層及裂縫成像結(jié)果生成單元，用于對所述HV型雙反射波的偏移結(jié)果和所述VH型雙反射波的偏移結(jié)果進行疊加成像，獲得目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果。本發(fā)明實施例提供的上述方法及裝置的有益技術(shù)效果在于本發(fā)明實施例基于雙反射疊前偏移方法，實現(xiàn)了專門針對陡傾角地層及裂縫的地震成像，該方法將與陡傾角地層及裂縫相關(guān)的雙反射波作為有效信號，對包括一次反射波在內(nèi)的其它各種類型的地震波進行壓制，從而能夠得到高質(zhì)量的陡傾角地層及裂縫成像結(jié)果。同時，本發(fā)明實施例在實現(xiàn)雙反射疊前偏移方法的過程中，采用保幅傅里葉有限差分算子進行延拓，并采用有限差分方法進行旅行時計算，因此，本發(fā)明實施例的方法及裝置具有較好的保幅性和橫向變速適應(yīng)能力。


為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn)有技術(shù)描述中所需要使用的附圖做一簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。圖Ia是現(xiàn)有技術(shù)的一次反射波的示意圖；圖Ib是現(xiàn)有技術(shù)的HV型雙反射波的示意圖；圖Ic是現(xiàn)有技術(shù)的VH型雙反射波的示意圖；圖2是本發(fā)明實施例的陡傾角地層及裂縫地震偏移方法的流程圖；圖加是本發(fā)明實施例的有限差分旅行時計算的矩形網(wǎng)格示意圖；圖3a是本發(fā)明實施例的陡傾角地層速度模型；圖北是本發(fā)明實施例的陡傾角地層模型的雙反射偏移結(jié)果；圖如是本發(fā)明實施例的裂縫速度模型；圖4b是本發(fā)明實施例的裂縫模型的雙反射偏移結(jié)果；圖5是本發(fā)明實施例的陡傾角地層及裂縫地震偏移裝置的功能框圖；圖6是本發(fā)明實施例的共炮點道集地震數(shù)據(jù)延拓處理單元508的具體功能框圖。
具體實施例方式為使本發(fā)明實施例的目的、技術(shù)方案和優(yōu)點更加清楚，下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進行清楚、完整地描述，顯然，所描述的實施例是本發(fā)明一部分實施例，而不是全部的實施例?；诒景l(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護的范圍。本發(fā)明實施例是提供一種具有較高計算效率、專門針對陡傾角地層及裂縫成像的地震數(shù)據(jù)疊前偏移新方法。圖2為本發(fā)明實施例的一種陡傾角地層及裂縫地震偏移方法的流程圖。如圖2所示，本發(fā)明實施例的雙反射疊前偏移方法，具體步驟包括202、輸入地震數(shù)據(jù)，所述地震數(shù)據(jù)包括通過對野外采集的地震數(shù)據(jù)進行預(yù)處理得到的共炮點道集地震數(shù)據(jù)、通過常規(guī)疊前深度偏移處理得到的深度域偏移速度模型以及常規(guī)疊前深度偏移結(jié)果；204、利用常規(guī)疊前深度偏移結(jié)果，在深度域定義一個基礎(chǔ)地下水平反射層；該水平反射層的深度應(yīng)該大于所要研究的陡傾角地層或裂縫的深度，并且應(yīng)是一個強反射層位，同時還應(yīng)是偏移速度模型的層位之一。另外，該水平層不一定是完全水平的反射層位，可以是局部帶有一定傾角的強反射界面。206、將步驟202中輸入的共炮點道集地震數(shù)據(jù)進行重新排列，生成共檢波點道集地震數(shù)據(jù)，這兩種道集數(shù)據(jù)分別用于不同類型雙反射波的疊前偏移，從而能夠保證程序?qū)崿F(xiàn)的一致性；在本發(fā)明實施例中，為了程序?qū)崿F(xiàn)的一致性，對于HV和VH兩種不同類型的雙反射波分別采用不同類型的地震數(shù)據(jù)進行偏移處理首先，HV型雙反射波(圖lb)的傳播路徑為炮點-地下水平界面-陡傾角地層或裂縫-檢波點，對于該類雙反射波，采用共炮點道集地震數(shù)據(jù)進行偏移處理；而VH型雙反射波(圖Ic)的傳播路徑為炮點-陡傾角地層或裂縫-地下水平界面-檢波點，對于該類雙反射波，采用共檢波點道集地震數(shù)據(jù)進行偏移處理。208、利用保幅傅里葉有限差分延拓算子對共炮點道集地震數(shù)據(jù)進行逐層延拓處理；在本發(fā)明實施例中，對于HV型雙反射波，利用共炮點道集地震數(shù)據(jù)^； 行偏移處理，首先選取延拓算子對每一個共炮點道集地震數(shù)據(jù)以步長△ ζ進行逐層延拓，延拓過程中所用到的速度模型為步驟202中輸入的深度域偏移速度模型。在本發(fā)明實施例中，延拓算子采用的是保幅傅里葉有限差分延拓算子，該算子兼有相移法和有限差分方法的優(yōu)點， 對于地下劇烈橫向變速具有很強的適應(yīng)能力，同時通過增加振幅恢復(fù)項能夠有效提高偏移結(jié)果的保幅性，并且能夠補償?shù)卣鸩▊鞑ミ^程中的能量損失，尤其是有利于改善深層成像質(zhì)量。對共炮點道集地震數(shù)據(jù)戶“Hi)利用保幅傅里葉有限差分算子進行延拓的具體實現(xiàn)過程如下所示首先，輸入共炮點道集數(shù)據(jù)的單程波場，根據(jù)公式(1)，計算得到對應(yīng)的聲壓波場ps{x,z,m) = J^lps{x,z,m)⑴其中
權(quán)利要求
1.一種陡傾角地層及裂縫地震偏移方法，其特征在于，所述方法包括輸入地震數(shù)據(jù)，所述地震數(shù)據(jù)包括通過預(yù)處理得到的共炮點道集地震數(shù)據(jù)、深度域偏移速度模型以及常規(guī)疊前深度偏移結(jié)果；利用所述常規(guī)疊前深度偏移結(jié)果，在深度域定義一個基礎(chǔ)地下水平反射層； 對輸入的所述共炮點道集地震數(shù)據(jù)進行重新排列，生成共檢波點道集地震數(shù)據(jù)，所述共炮點道集地震數(shù)據(jù)用于HV型雙反射波的疊前偏移，所述共檢波點道集地震數(shù)據(jù)用于對 VH型雙反射波的疊前偏移；利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共炮點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共炮點道集地震數(shù)據(jù)；利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算HV型雙反射波偏移中需要的第一旅行時信息；根據(jù)所述延拓后的共炮點道集地震數(shù)據(jù)和所述第一旅行時信息，并利用成像條件對HV 型雙反射波以步長Δ ζ進行逐層成像，得到HV型雙反射波的偏移結(jié)果；利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共檢波點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共檢波點道集地震數(shù)據(jù)；利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算VH型雙反射波偏移中需要的第二旅行時信息；根據(jù)所述延拓后的共檢波點道集地震數(shù)據(jù)和所述第二旅行時信息，并利用成像條件對 VH型雙反射波以步長Δ ζ進行逐層成像，得到VH型雙反射波的偏移結(jié)果；對所述HV型雙反射波的偏移結(jié)果和所述VH型雙反射波的偏移結(jié)果進行疊加成像，獲得目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果。
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述基礎(chǔ)地下水平反射層的深度大于所要研究的陡傾角地層或裂縫的深度，所述基礎(chǔ)地下水平反射層是一個強反射層位，并且是所述深度域偏移速度模型的層位之一。
3.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述基礎(chǔ)地下水平反射層是水平的反射層位，或者是局部帶有一定傾角的強反射界面。
4.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共炮點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共炮點道集地震數(shù)據(jù)包括根據(jù)輸入的共炮點道集數(shù)據(jù)的單程波場，獲得對應(yīng)的聲壓波場 根據(jù)所述聲壓波場，獲得針對背景速度的相移項； 根據(jù)所述針對背景速度的相移項，獲得針對背景速度的振幅恢復(fù)項； 根據(jù)所述針對背景速度的振幅恢復(fù)項，獲得針對變速擾動的時移項； 根據(jù)所述針對變速擾動的時移項，獲得針對變速擾動的振幅恢復(fù)項； 根據(jù)所述針對變速擾動的振幅恢復(fù)項，獲得有限差分補償項； 根據(jù)所述有限差分補償項，獲得有限差分補償?shù)恼穹謴?fù)項，其中所述有限差分補償?shù)恼穹謴?fù)項為延拓后的共炮點道集地震數(shù)據(jù)。
5.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述第一旅行時信息為=Thv= T1 (s， h)+T2(h，ο)，其中，T1 (s，h)為從炮點s到所述基礎(chǔ)地下水平反射層上反射點h的旅行時，T2(h, ο)為從所述基礎(chǔ)地下水平反射層上反射點h到地下目標(biāo)成像網(wǎng)格點ο的旅行時。
6.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述第二旅行時信息為TVH= T1 (r， h)+T2(h，o)，其中，T1 (r,h)為從檢波點r到所述基礎(chǔ)地下水平反射層上反射點h的旅行時， T2(h, ο)為從所述基礎(chǔ)地下水平反射層上反射點h到地下目標(biāo)成像網(wǎng)格點ο的旅行時。
7.根據(jù)權(quán)利要求5所述的方法，其特征在于，利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算HV型雙反射波偏移中需要的第一旅行時信息包括采用所述有限差分旅行時計算方法分別計算所述T1 (s，h)和所述T2 (h，ο)，所述有限差分進行旅行時計算方法包括如下關(guān)系式
8.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果滿足下列關(guān)系式I (nx, nz) = Ihv (nx, nz)+IVH(nx, nz)其中，I(nx，nz)表示目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果，Ihv (ηχ,ηζ)表示HV型雙反射波偏移結(jié)果，Ivh(nx, nz)表示VH型雙反射波偏移結(jié)果。
9.一種陡傾角地層及裂縫地震偏移裝置，其特征在于，所述裝置包括地震數(shù)據(jù)輸入單元，用于輸入地震數(shù)據(jù)，所述地震數(shù)據(jù)包括通過預(yù)處理得到的共炮點道集地震數(shù)據(jù)、深度域偏移速度模型以及常規(guī)疊前深度偏移結(jié)果；基礎(chǔ)地下水平反射層定義單元，用于利用所述常規(guī)疊前深度偏移結(jié)果，在深度域定義一個基礎(chǔ)地下水平反射層；共檢波點道集地震數(shù)據(jù)生成單元，用于對輸入的所述共炮點道集地震數(shù)據(jù)進行重新排列，生成共檢波點道集地震數(shù)據(jù)，所述共炮點道集地震數(shù)據(jù)用于HV型雙反射波的疊前偏移，所述共檢波點道集地震數(shù)據(jù)用于對VH型雙反射波的疊前偏移；共炮點道集地震數(shù)據(jù)延拓處理單元，用于利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共炮點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共炮點道集地震數(shù)據(jù)；第一旅行時信息計算單元，用于利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算HV型雙反射波偏移中需要的第一旅行時信息；HV型雙反射波的偏移結(jié)果生成單元，用于根據(jù)所述延拓后的共炮點道集地震數(shù)據(jù)和所述第一旅行時信息，并利用成像條件對HV型雙反射波以步長Δ ζ進行逐層成像，得到HV型雙反射波的偏移結(jié)果；共檢波點道集地震數(shù)據(jù)延拓處理單元，用于利用保幅傅里葉有限差分延拓算子并采用所述深度域偏移速度模型對所述共檢波點道集地震數(shù)據(jù)進行逐層延拓處理，得到延拓后的共檢波點道集地震數(shù)據(jù)；第二旅行時信息計算單元，用于利用有限差分旅行時計算方法并采用所述深度域偏移速度模型，同時根據(jù)所述基礎(chǔ)地下水平反射層，計算VH型雙反射波偏移中需要的第二旅行時信息；VH型雙反射波的偏移結(jié)果生成單元，用于根據(jù)所述延拓后的共檢波點道集地震數(shù)據(jù)和所述第二旅行時信息，并利用成像條件對VH型雙反射波以步長Δ ζ進行逐層成像，得到VH 型雙反射波的偏移結(jié)果；陡傾角地層及裂縫成像結(jié)果生成單元，用于對所述HV型雙反射波的偏移結(jié)果和所述 VH型雙反射波的偏移結(jié)果進行疊加成像，獲得目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果。
10.根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述基礎(chǔ)地下水平反射層的深度大于所要研究的陡傾角地層或裂縫的深度，所述基礎(chǔ)地下水平反射層是一個強反射層位，并且是所述深度域偏移速度模型的層位之一。
11.根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述基礎(chǔ)地下水平反射層是水平的反射層位，或者是局部帶有一定傾角的強反射界面。
12.根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述共炮點道集地震數(shù)據(jù)延拓處理單元具體包括第一計算模塊，用于根據(jù)輸入的共炮點道集數(shù)據(jù)的單程波場，獲得對應(yīng)的聲壓波場 第二計算模塊，用于根據(jù)所述聲壓波場，獲得針對背景速度的相移項； 第三計算模塊，用于根據(jù)所述針對背景速度的相移項，獲得針對背景速度的振幅恢復(fù)項；第四計算模塊，用于根據(jù)所述針對背景速度的振幅恢復(fù)項，獲得針對變速擾動的時移項；第五計算模塊，用于根據(jù)所述針對變速擾動的時移項，獲得針對變速擾動的振幅恢復(fù)項；第六計算模塊，用于根據(jù)所述針對變速擾動的振幅恢復(fù)項，獲得有限差分補償項； 第七計算模塊，用于根據(jù)所述有限差分補償項，獲得有限差分補償?shù)恼穹謴?fù)項，其中所述有限差分補償?shù)恼穹謴?fù)項為延拓后的共炮點道集地震數(shù)據(jù)。
13.根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述第一旅行時信息為THV= I\(s， h)+T2(h，ο)，其中，T1 (s，h)為從炮點s到所述基礎(chǔ)地下水平反射層上反射點h的旅行時， T2(h, ο)為從所述基礎(chǔ)地下水平反射層上反射點h到地下目標(biāo)成像網(wǎng)格點ο的旅行時。
14.根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述第二旅行時信息為TVH= 1\0·， h)+T2(h，o)，其中，T1 (r,h)為從檢波點r到所述基礎(chǔ)地下水平反射層上反射點h的旅行時， T2(h, ο)為從所述基礎(chǔ)地下水平反射層上反射點h到地下目標(biāo)成像網(wǎng)格點ο的旅行時。
15.根據(jù)權(quán)利要求13所述的裝置，其特征在于，所述第一旅行時信息計算單元，具體用于采用所述有限差分旅行時計算方法分別計算所述T1 (s，h)和所述T2 (h，0)，所述有限差分進行旅行時計算方法包括如下關(guān)系式
16.根據(jù)權(quán)利要求9所述的裝置，其特征在于，所述陡傾角地層及裂縫成像結(jié)果生成單元，具體用于基于如下關(guān)系式獲得所述目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果 I (nx, nz) = Ihv (nx, nz)+IVH(nx, nz)其中，I(nx，nz)表示目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果，Ihv (ηχ,ηζ)表示HV型雙反射波偏移結(jié)果，Ivh(nx, nz)表示VH型雙反射波偏移結(jié)果。
全文摘要
本發(fā)明提供一種陡傾角地層及裂縫地震偏移方法和裝置，該方法包括輸入地震數(shù)據(jù)；定義基礎(chǔ)地下水平反射層；生成共檢波點道集地震數(shù)據(jù)；利用保幅傅里葉有限差分延拓算子對共炮點道集地震數(shù)據(jù)進行逐層延拓處理；利用有限差分旅行時計算方法計算HV型雙反射波偏移中的旅行時；對HV型雙反射波進行逐層成像；利用保幅傅里葉有限差分延拓算子對共檢波點道集地震數(shù)據(jù)進行逐層延拓處理；利用有限差分旅行時計算方法計算VH型雙反射波偏移中的旅行時；對VH型雙反射波進行逐層成像；對HV和VH兩種不同雙反射波的偏移結(jié)果進行疊加成像，得到目標(biāo)區(qū)的陡傾角地層及裂縫成像結(jié)果。該方法計算效率較高，能得到高質(zhì)量的陡傾角地層及裂縫成像結(jié)果。
文檔編號G01V1/28GK102540253SQ20111045637
公開日2012年7月4日 申請日期2011年12月30日 優(yōu)先權(quán)日2011年12月30日
發(fā)明者呂彬, 王宇超, 胡自多, 蘇勤, 雍學(xué)善 申請人:中國石油天然氣股份有限公司