專利名稱：振幅保存疊前偏移方法
技術(shù)領(lǐng)域：
一般來說，本發(fā)明涉及地下地質(zhì)構(gòu)造的地震測量領(lǐng)域。更具體地說，本發(fā)明涉及基本上保存地震數(shù)據(jù)中的事件振幅的用于地震數(shù)據(jù)的疊前偏移的方法。
背景技術(shù)：
地震測量用于評估地下地質(zhì)構(gòu)造的結(jié)構(gòu)、成分和流體含量。地震測量的特定應(yīng)用是推斷有用物質(zhì)、如石油在地下地質(zhì)構(gòu)造中的存在。一般來說，地震測量包括在地表或地表附近部署地震傳感器陣列，并且還在地表或地表附近部署地震能量源或此類源的陣列。地震能量源受促動并且地震能量從源發(fā)出，一般向下通過地下穿行，直至到達(dá)地下的一個(gè)或多個(gè)聲阻抗邊界。地震能量從一個(gè)或多個(gè)阻抗邊界反射，然后向上穿行直至被一個(gè)或多個(gè)傳感器檢測。地下地質(zhì)結(jié)構(gòu)和成分從包括檢測能量的其它屬性在內(nèi)，地震能量的穿行時(shí)間以及就從地震源發(fā)出的能量而言地震能量的各種頻率分量的振幅和相位來推斷。
為了從在地表或地表附近測量的從地表或地表附近的源位置的地震穿行時(shí)間推斷地下地質(zhì)構(gòu)造的結(jié)構(gòu)，必須確定地震能量穿過各種構(gòu)造的速度。地質(zhì)構(gòu)造的速度可視地下深度(縱向)和地理位置(橫向)而有所不同。但是，地震數(shù)據(jù)僅按時(shí)間記錄。本領(lǐng)域中用于估算地質(zhì)構(gòu)造的縱向和橫向速度的已知方法依靠關(guān)于地震能量從地震源穿行到部署在地表的各種地震接收器時(shí)地震能量的穿行路徑幾何形狀的推論。
為了讓從地震數(shù)據(jù)產(chǎn)生的圖像準(zhǔn)確對應(yīng)于地下結(jié)構(gòu)的空間分布和地下地質(zhì)的成分變化，對地震數(shù)據(jù)執(zhí)行通常稱作“時(shí)間偏移”和“深度偏移”的技術(shù)。偏移是進(jìn)行地震數(shù)據(jù)中的反射事件以對應(yīng)于其真實(shí)空間位置的過程。時(shí)間偏移按時(shí)間定位反射事件并且不校正橫向速度梯度的所有效果。深度偏移嘗試按反射事件在地下地質(zhì)中的真實(shí)深度位置定位反射事件。
一般來說，有兩類時(shí)間偏移技術(shù)。一種時(shí)間偏移稱作疊前時(shí)間偏移，并且公開于例如Sun，C.，Martinez，R.，Amplitude preserving3D pre-stack Kirchhoff time migration for V(z)and VTI media，72ndAnnual International Meeting，Society of Exploration Geophysicists，Expanded Abstracts，pp.1224-1227(2002)中。疊前偏移在各個(gè)地震數(shù)據(jù)記錄(“軌跡”)上執(zhí)行并且可能是計(jì)算上昂貴的。
另一種時(shí)間偏移稱為疊后偏移，代表將各個(gè)數(shù)據(jù)記錄的數(shù)字求和(“堆疊”)以提高信噪比的對地震數(shù)據(jù)執(zhí)行的偏移技術(shù)。疊前偏移通常產(chǎn)生更好的圖像，雖然它在計(jì)算上比疊后偏移昂貴。
在Sun等人的上述文件中公開的疊前時(shí)間偏移方法是稱為“振幅保存”的類型之一，意味著在偏移過程中地震數(shù)據(jù)中事件的振幅基本上被保存。除其它好處之外，振幅保存偏移還在稱為“振幅與偏移量”(AVO)的處理技術(shù)中提供更大的準(zhǔn)確性。振幅保存偏移技術(shù)的主要目的在于使反射振幅與各種地下邊界處的反射率成比例。因此，振幅保存偏移的特征在于，恢復(fù)反射地震能量中由通過土地傳播時(shí)地震能量的幾何擴(kuò)展導(dǎo)致的振幅損失，使得振幅至少在某種程度上對應(yīng)于各種地下邊界的反射率。
疊前偏移圖像通過映射輸入地震數(shù)據(jù)軌跡的起作用的數(shù)據(jù)樣本至圖像軌跡的目標(biāo)樣本來計(jì)算，受到起作用的數(shù)據(jù)樣本時(shí)間(索引到地震源的促動)與地震能量通過圖像點(diǎn)從地震源至接收器的穿行時(shí)間相同的約束。然后，對來自所有輸入地震數(shù)據(jù)軌跡的起作用的數(shù)據(jù)樣本求和。所謂的“真振幅”偏移，即恢復(fù)由幾何擴(kuò)展導(dǎo)致的反射振幅損失以使得振幅對應(yīng)于地下反射體的反射率，是通過輸入數(shù)據(jù)樣本的加權(quán)求和而不是簡單求和來實(shí)現(xiàn)的。也就是說，起作用的數(shù)據(jù)樣本的振幅代表乘以加權(quán)因子的輸入數(shù)據(jù)樣本的和。加權(quán)因子是多個(gè)不同變量的函數(shù)，變量包括穿行時(shí)間以及對應(yīng)于每個(gè)數(shù)據(jù)樣本軌跡(稱為“采集幾何”)的源和接收器的位置。本領(lǐng)域中已知幾種不同技術(shù)可用于確定“真振幅”權(quán)函數(shù)。例如，請參見Schleicher，J.，Tygel，M.，and Hubral，P.，1993，3-D true amplitude finite-offsetmigration；Geophysics，58，1112-1126或Bleistein，N.，Cohen，J.K.，Stockwell，J.W.，Mathematics of multidimensional seismic imaging，migration，and inversion，Springer，2001。用于最常見的地震數(shù)據(jù)采集幾何的真振幅權(quán)函數(shù)，如共射點(diǎn)或共偏移幾何，包括包含起飛和出射角的余弦的平方根的項(xiàng)，以及用于源、接收器和圖像點(diǎn)的各種可能組合，甚至用于相對簡單的v(z)介質(zhì)情況的面內(nèi)和面外擴(kuò)展項(xiàng)。例如，請參見Zhang，Y.，Gray，S.，Young，J.，Exact and approximate weightsfor Kirchhoff migration，70thAnnual Internat.Mtg.，Soc.Expl.Geophys.Expanded Abstracts，1036-1039(2000)。因?yàn)闄?quán)重的生成需要大量的計(jì)算并且權(quán)重必須在偏移期間存儲和內(nèi)插，所以權(quán)函數(shù)的計(jì)算可以左右偏移過程的總計(jì)算成本。
本領(lǐng)域中已知以各種方法近似計(jì)算權(quán)函數(shù)以避免增大的輸入/輸出和計(jì)算開支。例如，請參見Dellinger，J.A.，Gray，S.H.，Murphy，G且，and Etgen，J.T.，Efficient 2.5-D true-amplitude migration；Geophysics，65，943-950(2000)。通常此類估算使用地質(zhì)構(gòu)造的恒定速度近似值進(jìn)行。但是，已經(jīng)確定使用恒定速度近似值估算的權(quán)函數(shù)經(jīng)常無法為在v(z)(速度隨深度變化)介質(zhì)中處理所謂的“振幅與偏移量”(“AVO”)提供滿意的結(jié)果。需要提供與本領(lǐng)域已知的“確切”真振幅技術(shù)所提供的相似的振幅質(zhì)量，而不會大量增加計(jì)算開支的疊前偏移技術(shù)。

發(fā)明內(nèi)容
本發(fā)明的一個(gè)方面是用于振幅保存Kirchhoff時(shí)間偏移的方法。根據(jù)本發(fā)明這一方面的方法包括使用恒定速度近似值計(jì)算幾何擴(kuò)展項(xiàng)，同時(shí)根據(jù)隨穿行時(shí)間變化的速度的模型準(zhǔn)確計(jì)算起飛和出射角。
本發(fā)明的另一方面是存儲在計(jì)算機(jī)可讀介質(zhì)中的計(jì)算機(jī)程序。該程序包括邏輯，可用于使可編程計(jì)算機(jī)執(zhí)行包括使用恒定速度近似值計(jì)算幾何擴(kuò)展項(xiàng)，同時(shí)根據(jù)隨穿行時(shí)間變化的速度的模型準(zhǔn)確計(jì)算起飛和出射角的步驟。
通過以下描述和所附權(quán)利要求，本發(fā)明的其它方面和優(yōu)點(diǎn)將會十分明顯。


圖1顯示根據(jù)本發(fā)明的方法的示范實(shí)施例的流程圖。
圖2顯示用于模擬偏移過程的示范地質(zhì)模型。
圖3和圖4顯示根據(jù)本發(fā)明從使用分析上精確的權(quán)重和估算的權(quán)重的偏移進(jìn)行的振幅與偏移量(“AVO”)分析。
圖5顯示存儲了程序的可編程計(jì)算機(jī)和計(jì)算機(jī)可讀介質(zhì)。
具體實(shí)施例方式根據(jù)本發(fā)明的方法在Kirchhoff疊前時(shí)間偏移的上下文中可以更好地理解。Kirchhoff疊前時(shí)間偏移基于標(biāo)量波動方程的Kirchhoff定積分解。Kirchhoff定積分解是惠更斯原理的數(shù)學(xué)表達(dá)式?；莞乖碚f明波場，在體積中的任何內(nèi)點(diǎn)，可通過波場邊界值的表面積分重構(gòu)，按對波算子的格林函數(shù)脈沖響應(yīng)加權(quán)。
Kirchhoff疊前時(shí)間偏移可以按以下表達(dá)式使用由(xs，ys)表示的射點(diǎn)(地震源)位置、由(xr，yr)表示的地震接收器位置以及由(Ti，xi，yi)表示的時(shí)間和空間中的圖像點(diǎn)位置在數(shù)學(xué)上表達(dá)
I(Ti，xi，yi)＝ ∫W(xi，yi，Ti，xs，ys，xr，yr)D(xs，ys，xr，yr，Ts+Tr)dxsdxrdysdyr(1) 其中，I(Ti，xi，yi)代表處于由(xi，yi)表示的等效表面(地理)位置的圖像(振幅)，Ti代表圖像時(shí)間，W(xi，yi，Ti，xs，ys，xr，yr)是加權(quán)函數(shù)，D(xs，ys，xr，yr，Ts+Tr)代表輸入(地震軌跡)數(shù)據(jù)的時(shí)間導(dǎo)數(shù)，Ttot＝TS+Tr是地震能量通過圖像點(diǎn)(xi，yi)從每個(gè)源位置(xs，ys)穿行至每個(gè)地震接收器位置(xr，yr)的總穿行時(shí)間。Tr和Ts分別代表地震能量從圖像點(diǎn)到接收器位置的穿行時(shí)間以及從源位置到圖像點(diǎn)的穿行時(shí)間。為每個(gè)源和接收器位置計(jì)算圖像點(diǎn)的穿行時(shí)間Ttot的優(yōu)選公式如下所示，可用于各向同性[v(z)]介質(zhì)和/或弱VTI(橫向上具有對稱垂直軸的各向同性)介質(zhì)Ttot=Tr+Ts+ΔT(VTI)]]>=Tr0(1+12CCc4xr6Tr02)+Ts0(1+12CCc4xs6Ts02)+ΔT(VTI)---(2)]]> 其中，CC為常數(shù)，xs和xr代表源和接收器偏移量，即源位置和圖像位置以及接收器位置和圖像位置之間的橫向距離，Tro和Tso分別代表包含圖像點(diǎn)和接收器以及圖像點(diǎn)和源之間偏移量的最多四次項(xiàng)的單程穿行時(shí)間。從等式(1)和(2)的上述參數(shù)的計(jì)算在Sun，C.，Martinez，R.，Amplitude preserving 3D pre-stack Kirchhoff timemigration for V(z)and VTI media，72nd Annual International Meeting，Society of Exploration Geophysicists，Expanded Abstracts，pp.1224-1227(2002)中說明。等式(2)中的ΔT(VTI)代表“直射線”(地震能量從源到圖像點(diǎn)和從圖像點(diǎn)到接收器點(diǎn)沿直射線路徑穿行)和“曲射線”(地震能量由于地下速度不同所導(dǎo)致的折射引起的沿曲路徑穿行)之間的穿行時(shí)間差異，其中介質(zhì)為VTI。如Sun和Martinez出版物中所述，量ΔT(VTI)可以使用以下表達(dá)式計(jì)算ΔT(VTI)=C(VTI)x4s2T20+x2sV2rms+C(VTI)x4r2T20+x2rV2rms---(3)]]> 其中量C(VTI)可由以下表達(dá)式確定C(VTI)=2ηV2rms[T02V2rms+(1+2η)x2]---(4)]]> 以上等式中的η代表各向異性參數(shù)。
如本文背景部分所述，計(jì)算振幅保存偏移的加權(quán)因子在以分析上精確的方式執(zhí)行時(shí)是計(jì)算上昂貴的運(yùn)算。如Sun和Martinez出版物中所公開，可以近似計(jì)算振幅保存權(quán)重。但是，使用本領(lǐng)域中已知方法的此類近似權(quán)重已證實(shí)具有準(zhǔn)確度限制。
一般來說，根據(jù)本發(fā)明的方法包括通過使用恒定速度近似值以及使用采用隨穿行時(shí)間不斷變化的速度模型即時(shí)準(zhǔn)確計(jì)算的起飛和出射角估算幾何擴(kuò)展項(xiàng)，用于以上v(z)(速度與深度相關(guān))介質(zhì)中的Kirchhoff時(shí)間偏移等式的所謂“真振幅”權(quán)重的近似。通過準(zhǔn)確計(jì)算起飛和出射角，估算的權(quán)重可用于產(chǎn)生比僅為權(quán)重計(jì)算和起飛及出射角計(jì)算使用恒定速度近似值的方式更準(zhǔn)確的圖像。即時(shí)計(jì)算權(quán)重減少CPU時(shí)間的高需求和大存儲空間需求，因此使得該方法十分有效。
根據(jù)本發(fā)明的方法的理論可說明如下。如Zhang等人出版物中所述，本文中背景部分所引用，在具有共偏移量幾何的v(z)(速度與深度相關(guān))介質(zhì)中的“真振幅”權(quán)函數(shù)(如Kirchhoff偏移等式(1)中使用)由以下表達(dá)式給定 w=cosαs0cosαr0v0[(ψsψr+ψrψs)(σsσr+σrσs)+sin2γ2cos2θL]---(5)]]> 其中γ是源投影和到地表的接收器射線之間的角度，L由以下表達(dá)式給定 L=(cosαs+cosαr)(ψsψrσsσr+σsσrψsψr)]]>-(1+cosαscosαr)(ψsσrψrσs+ψrσsψsσr)]]>(6) 在上述等式(5)和(6)中，αs代表沿地震能量射線路徑從地震源相對于垂直的角度，ψs和σs是從地震源位置到圖像點(diǎn)的面內(nèi)和面外幾何擴(kuò)展項(xiàng)。每個(gè)地震接收器位置的等效定義適用于αr、ψr和σr。即使忽略等式(5)中的最后一項(xiàng)Lsin2γ/(2cos2θ)(可能因?yàn)樗姆冉?jīng)常不大)，在以上形式的等式(6)中包括面內(nèi)和面外擴(kuò)展項(xiàng)仍然是計(jì)算上昂貴的。計(jì)算開支的產(chǎn)生是因?yàn)槊鎯?nèi)和面外擴(kuò)展項(xiàng)必須提前計(jì)算并且存儲在表格(涉及昂貴的臨時(shí)或永久計(jì)算機(jī)存儲空間)中并且在偏移期間內(nèi)插。
對于恒定速度介質(zhì)，具有由v0代表的速度，但是，振幅權(quán)函數(shù)可簡化為以下表達(dá)式 w=tm2v0(tstr+trts)(1tr+1ts)---(7)]]> 在上述表達(dá)式中，tm為雙程圖像時(shí)間，ts和tr分別是從源位置到圖像點(diǎn)以及從圖像點(diǎn)到接收器位置的穿行時(shí)間。
使用如上述等式(7)中計(jì)算的恒定速度振幅權(quán)函數(shù)，但是，加權(quán)因子通常被低估，主要是因?yàn)樵刺幍钠痫w角αs0和接收器處的出射角αr0由于假定的直射線路徑被高估。直射線路徑是在假定介質(zhì)具有恒定速度時(shí)假定的。
為了提高在使用此類近似權(quán)重時(shí)偏移結(jié)果的準(zhǔn)確性，在根據(jù)本發(fā)明的方法中，使用等式(9)和(10)為隨時(shí)間或深度變化的速度模型準(zhǔn)確計(jì)算起飛和出射角，如下所述。然后，可以使用恒定速度近似值計(jì)算等式(5)的擴(kuò)展項(xiàng)的比率。最后，等式(5)中的最后項(xiàng)可以忽略，因?yàn)樗ǔ７葮O小。使用這些近似值，振幅權(quán)函數(shù)可以表達(dá)為 w=cosαs0cosαr0v0tstr(tstr+trts)(1tr+1ts)---(8)]]> 因?yàn)樗俣燃俣榕c徑向(橫向或表面位置)坐標(biāo)無關(guān)，所以慢矢量的徑向分量為v(z)介質(zhì)保存(斯涅爾定律)。所以，源位置的起飛角的余弦可以通過以下等式準(zhǔn)確計(jì)算 cosαs0=1-(∂ts∂ρs)2v02---(9)]]> 同樣，出射角的余弦可以通過以下表達(dá)式計(jì)算 cosαr0=1-(∂tr∂ρr)2v02]]>(10) 在上述等式(9)中，ts是從源位置到圖像點(diǎn)的穿行時(shí)間，ρs是從v(z)介質(zhì)中的源的徑向的坐標(biāo)，v0是在表面的速度。從圖像點(diǎn)到接收器(r)位置的對應(yīng)項(xiàng)在等式(10)中。因?yàn)榇┬袝r(shí)間在分析上計(jì)算，所以等式(9)和(10)中的導(dǎo)數(shù)也會在分析上計(jì)算。很明顯，等式(8)的權(quán)函數(shù)在計(jì)算上比等式(7)更加昂貴地執(zhí)行，但獲得的結(jié)果可遠(yuǎn)遠(yuǎn)好于使用等式(7)中的恒定速度近似值獲得的結(jié)果。
計(jì)算完權(quán)函數(shù)，可以計(jì)算總穿行時(shí)間，例如，使用表達(dá)式，如用于VTI介質(zhì)中的曲射線的等式(2)。然后，穿行時(shí)間和權(quán)函數(shù)可用于使用Kirchhoff時(shí)間偏移等式(1)生成圖像。
根據(jù)本發(fā)明的方法的一個(gè)實(shí)施例如圖1中的流程圖所示。在10，對地震能量的雙程穿行時(shí)間，生成地震速度的模型。使用在10生成的模型，在12，可以使用恒定速度近似值估算權(quán)重。在14，使用公式計(jì)算起飛和出射角。在16，權(quán)重和起飛及出射角用于生成偏移圖像(即輸入地震數(shù)據(jù)樣本根據(jù)等式(1)映射至圖像樣本)。
另一方面，本發(fā)明涉及存儲在計(jì)算機(jī)可讀介質(zhì)中的計(jì)算機(jī)程序。參照圖5，如參照圖1所述的以上過程可以體現(xiàn)為計(jì)算機(jī)可讀代碼，存儲在計(jì)算機(jī)可讀介質(zhì)中，比如組成通用可編程計(jì)算機(jī)的一部分的軟盤68、CD-ROM 70或磁硬盤66。如本領(lǐng)域中已知的，計(jì)算機(jī)包括中央處理單元60、用戶輸入裝置、如鍵盤62和用戶顯示器64、如平板LCD顯示器或陰極射線管顯示器。根據(jù)本發(fā)明的這一方面，計(jì)算機(jī)可讀介質(zhì)包括可用于使計(jì)算機(jī)按以上所述和參照圖1說明的來執(zhí)行動作的邏輯。
三維合成模型，代表頁巖中嵌入的兩個(gè)充滿海水的砂層，用于測試上述實(shí)施例。圖2中顯示模型沿一面的圖形表示。水層在由參考編號1和2表示的層邊界之間表示。具有不同速度的頁巖層表示為布置在邊界2和3、4和5、5和6、6和7之間。每個(gè)頁巖層建模為具有不同速度，以便在每個(gè)此類邊界處提供地震反射體。另一個(gè)頁巖層顯示在邊界8下方，但是該模型在邊界8以下未提供任何反射體。兩個(gè)砂層中最上面的一個(gè)顯示為布置在由編號3和4表示的邊界之間。最上面的砂層是“平坦的”(即結(jié)構(gòu)傾斜為零)，下層傾斜砂層，在標(biāo)示為7、8的邊界之間表示，被建模為具有10度的傾斜。模型具有25×12.5米的箱柜空格。建模中使用的雙源容器具有相隔100米的8個(gè)長條。每個(gè)長條有313個(gè)接收器，組間隔為25m，偏移范圍在200m和8000m之間。
常見圖像采集，包含100m至8000m偏移范圍，偏移增量為100m，使用振幅保存Kirchhoff時(shí)間偏移來偏移，如上所述。偏移使用先有技術(shù)權(quán)函數(shù)近似方法、分析上精確的權(quán)函數(shù)方法以及根據(jù)本發(fā)明計(jì)算權(quán)重的方法的實(shí)施例執(zhí)行，如上所述。然后在傾斜砂層的頂部(邊界7)和底部(邊界8)執(zhí)行振幅與偏移量(“AVO”)分析。結(jié)果，分別顯示在圖3和圖4中，表明等式(8)中的振幅權(quán)函數(shù)與先有技術(shù)近似值等式(7)相比大大提高了振幅準(zhǔn)確性，特別是在遠(yuǎn)偏移處。在圖3和圖4中，從使用分析上精確的權(quán)重的偏移進(jìn)行的AVO分析如曲線20所示。使用先有技術(shù)權(quán)重近似值的偏移地震數(shù)據(jù)的AVO分析如曲線22所示。最后，使用本發(fā)明方法計(jì)算權(quán)重的偏移地震數(shù)據(jù)的AVO分析如曲線24所示。
與計(jì)算權(quán)重的精確、確定性方法相比，根據(jù)本發(fā)明的方法和計(jì)算機(jī)程序可以以更少的計(jì)算成本為振幅保存Kirchhoff疊前時(shí)間偏移提供準(zhǔn)確的加權(quán)因子。
雖然本發(fā)明針對數(shù)量有限的實(shí)施例進(jìn)行了說明，但是得益于此公開的本領(lǐng)域的技術(shù)人員將知道，只要不背離本文公開的發(fā)明范圍，可以設(shè)計(jì)其它實(shí)施例。因此，本發(fā)明的范圍應(yīng)當(dāng)僅受所附權(quán)利要求的限制。
權(quán)利要求
1.一種用于振幅保存Kirchhoff時(shí)間偏移的方法，包括使用恒定速度近似值計(jì)算幾何擴(kuò)展項(xiàng)的估算值；以及使用隨偏移時(shí)間和深度中至少一項(xiàng)而變化的速度的模型計(jì)算起飛和出射角。
2.如權(quán)利要求1所述的方法，其特征在于，還包括使用估算的權(quán)重和計(jì)算的起飛及出射角在至少一個(gè)樣本時(shí)間生成圖像。
3.如權(quán)利要求2所述的方法，其特征在于，生成圖像包括為多個(gè)地震源和地震接收器位置計(jì)算到至少一個(gè)圖像點(diǎn)的穿行時(shí)間，并且根據(jù)穿行時(shí)間和權(quán)重計(jì)算Kirchhoff定積分。
4.一種存儲在計(jì)算機(jī)可讀介質(zhì)中的計(jì)算機(jī)程序，所述程序具有可用于使可編程計(jì)算機(jī)執(zhí)行以下步驟的邏輯使用恒定速度近似值為振幅保存Kirchhoff時(shí)間偏移計(jì)算幾何擴(kuò)展項(xiàng)的估算值；以及使用隨偏移時(shí)間和深度中至少一項(xiàng)而變化的速度的模型計(jì)算起飛和出射角。
5.如權(quán)利要求4所述的程序，其特征在于，還包括可用于使計(jì)算機(jī)執(zhí)行使用估算的權(quán)重和計(jì)算的起飛及出射角在至少一個(gè)樣本時(shí)間生成圖像的邏輯。
6.如權(quán)利要求5所述的程序，其特征在于，生成圖像包括為多個(gè)地震源和地震接收器位置計(jì)算到至少一個(gè)圖像點(diǎn)的穿行時(shí)間，并且根據(jù)穿行時(shí)間和權(quán)重計(jì)算Kirchhoff定積分。
全文摘要
公開了一種用于振幅保存Kirchhoff時(shí)間偏移的方法。所述方法包括使用恒定速度近似值計(jì)算幾何擴(kuò)展項(xiàng)的估算值(12)。起飛和出射角根據(jù)隨穿行時(shí)間變化的速度模型準(zhǔn)確計(jì)算(14)。使用估算的權(quán)重和計(jì)算的起飛及出射角可以在至少一個(gè)穿行時(shí)間計(jì)算(16)圖像。
文檔編號G01V1/00GK101036069SQ200580033585
公開日2007年9月12日 申請日期2005年10月3日 優(yōu)先權(quán)日2004年10月4日
發(fā)明者李性馥, D·G·金 申請人:Pgs美洲公司