46] 圖7a是單點(diǎn)模型根據(jù)現(xiàn)有技術(shù)中等效偏移距偏移方法的成像結(jié)果；
[0047] 圖7b是單點(diǎn)模型根據(jù)本發(fā)明實(shí)施例的方法的成像結(jié)果；
[0048] 圖8a是水平界面模型根據(jù)現(xiàn)有技術(shù)中等效偏移距偏移方法的成像結(jié)果；
[0049] 圖8b是水平界面模型根據(jù)本發(fā)明實(shí)施例的方法的成像結(jié)果；
[0050] 圖9a是斷層模型根據(jù)現(xiàn)有技術(shù)中等效偏移距偏移方法的成像結(jié)果；
[0051] 圖9b是斷層模型根據(jù)本發(fā)明實(shí)施例的方法的成像結(jié)果；。
【具體實(shí)施方式】
[0052]本發(fā)明的實(shí)施例提供一種依據(jù)等效偏移距的疊前時(shí)間偏移成像方法，既利用近散 射點(diǎn)的能量，又利用遠(yuǎn)散射點(diǎn)處的能量，可以避免現(xiàn)有技術(shù)中只利用部分散射波能量成像 的不足，使成像過(guò)程中噪音變小，從而使得成像效果更加真實(shí)。
[0053] 實(shí)施例
[0054] 圖1是本實(shí)施例的基于等效偏移距的疊前時(shí)間偏移成像方法的流程圖。具體步驟 如圖所示：
[0055] 首先，在步驟S101中，抽取每個(gè)CSP道集并獲取均方根速度分布。
[0056] 隨后，執(zhí)行步驟S102,根據(jù)克?；舴蚱圃韺?duì)每個(gè)CSP道集進(jìn)行等效偏移距偏 移，即基于目標(biāo)散射點(diǎn)位置處的均方根速度對(duì)每個(gè)CSP道集進(jìn)行等效偏移距偏移，對(duì)目標(biāo) 散射點(diǎn)及其附近相對(duì)的較遠(yuǎn)散射點(diǎn)的散射波能量偏移收斂形成當(dāng)前散射點(diǎn)及其附近相對(duì) 較遠(yuǎn)散射點(diǎn)的偏移成像結(jié)果。如此以來(lái)，對(duì)每個(gè)CSP道集可以得到包含目標(biāo)散射點(diǎn)和附近 散射點(diǎn)的偏移成像結(jié)果的等效偏移道集。
[0057] 接下來(lái)，執(zhí)行步驟S103,從等效偏移道集中抽取與目標(biāo)散射點(diǎn)的水平面投影坐標(biāo) 相同的道形成共等效點(diǎn)道集。其中，將目標(biāo)散射點(diǎn)的水平面投影坐標(biāo)作為目標(biāo)散射點(diǎn)的等 效點(diǎn)的物理坐標(biāo)。對(duì)于偏移之后的CSP道集，按照目標(biāo)散射點(diǎn)的等效點(diǎn)的物理坐標(biāo)重行選 排，抽取與等效點(diǎn)的物理坐標(biāo)相同的道形成共等效點(diǎn)道集。
[0058] 最后，執(zhí)行步驟S104,對(duì)每個(gè)目標(biāo)散射點(diǎn)的共等效點(diǎn)道集進(jìn)行水平疊加，獲得全部 散射點(diǎn)的由近散射波能量和遠(yuǎn)散射波能量形成的成像結(jié)果。具體來(lái)說(shuō)，使得目標(biāo)散射點(diǎn)在 其CSP道集中的成像點(diǎn)的能量與目標(biāo)散射點(diǎn)在其他CSP道集中的來(lái)自同一成像點(diǎn)的能量疊 加，得到由近散射波能量和遠(yuǎn)散射波能量形成的全部散射點(diǎn)的成像結(jié)果。
[0059] 以下參照?qǐng)D2對(duì)步驟S102中對(duì)其中一個(gè)CSP道集進(jìn)行等效偏移距偏移的具體方 法進(jìn)行詳細(xì)說(shuō)明。
[0060] 首先，在步驟S201中，選擇CSP道集中的一道做為零偏移距道并選定零偏移距旅 行時(shí)間t。。
[0061] 然后，執(zhí)行步驟S202,首先，獲取指定的零偏移距道中時(shí)間t。對(duì)應(yīng)的目標(biāo)散射點(diǎn)處 的均方根速度v"s ;然后計(jì)算零偏移距道與其所在CSP道集中其他道，S卩非零偏移距道的水 平距離δ，以及每個(gè)δ匕對(duì)應(yīng)的走時(shí)h，其中i= 1，2……，N，N為覆蓋次數(shù)。其中，按照 下述表達(dá)式計(jì)算走時(shí)h
[0062] Cl)
[0063] ?女丨、本，彳土ZFW5山3中，分別計(jì)算非零偏移距道的每一走時(shí)h處的采樣點(diǎn)值的等 效偏移數(shù)值，即將CSP道集中非零偏移距道的走時(shí)h的采樣點(diǎn)值乘以?xún)A斜因子cosΘ，再乘 以球面擴(kuò)散因子R(v"s，r)，得到采樣點(diǎn)值的等效偏移數(shù)值cosΘR(v"s，r)Pin(tl)。
[0064] 其中，非零偏移距道在水平面的投影為走時(shí)h的采樣點(diǎn)的等效點(diǎn)。Θ是零偏移距 旅行時(shí)間t。對(duì)應(yīng)的目標(biāo)散射點(diǎn)與δ匕對(duì)應(yīng)的非零偏移距道在水平面投影的連線與垂直方 向的夾角，也就是目標(biāo)散射點(diǎn)與采樣點(diǎn)的等效點(diǎn)的連線與垂直方向的夾角。r是目標(biāo)散射點(diǎn) 與采樣點(diǎn)的等效點(diǎn)的距離，Vnis是目標(biāo)散射點(diǎn)處的均方根速度。PJh)為每個(gè)走時(shí)h處的 采樣點(diǎn)值。
[0065] 球面擴(kuò)散因子在二維偏移中表為-
在三維偏移中表為：
[0066] 在步驟S204中，將CSP道集中非零偏移距道的采樣點(diǎn)值的等效偏移數(shù)值疊加并輸 出到選定的零偏移距旅行時(shí)間t。對(duì)應(yīng)的目標(biāo)散射點(diǎn)。用下式表示：
[0067]
C2)
[0068] 其中，(t。）為零偏移距旅行時(shí)間對(duì)應(yīng)的目標(biāo)散射點(diǎn)的波場(chǎng)數(shù)值，Pin (tj為每個(gè) 走時(shí)t處的采樣點(diǎn)值，ΛX為CSP道集的道間距。
[0069] 接下來(lái)，在步驟S205中判斷是否已經(jīng)遍歷全部零偏移距旅行時(shí)間t。。如果是，則 執(zhí)行步驟S206,輸出CSP道集中其他非零偏移距道對(duì)零偏移距道的疊加道，如果否，則執(zhí)行 步驟S201，重新設(shè)定下一個(gè)零偏移距旅行時(shí)間t。。
[0070] 這樣可以遍歷該道中全部零偏移距旅行時(shí)間t。，由零偏移距道中全部目標(biāo)散射點(diǎn) 的等效偏移數(shù)值疊加結(jié)果得到零偏移距道的疊加道。
[0071] 隨后，在步驟S207中判斷是否已經(jīng)遍歷當(dāng)前CSP道集中所有的道。如果是，則執(zhí) 行步驟S208,輸出一個(gè)CSP道集中的全部疊加道，將全部疊加道組合得到CSP道集的疊加道 集；如果否，則執(zhí)行步驟S201，重新選擇CSP道集中的下一個(gè)道做為零偏移距道。
[0072] 然后，在步驟S209中，對(duì)疊加道組合得到CSP道集的疊加道集以濾波因子進(jìn)行褶 積形成等效偏移距偏移道集，表示為：
[0073]
C3)
[0074] 其中，Ρ_表示疊加道集，P(t)為濾波因子，m表示CSP道集編號(hào)。
[0075] Σ表示遍歷所指定的零偏移距道中的全部零偏移距旅行時(shí)間t。，計(jì)算并組合目 標(biāo)散射點(diǎn)的波場(chǎng)數(shù)值ρ_α。）。Σ表示遍歷當(dāng)前csp道集中所有的道，計(jì)算并組合疊加道， 得到CSP道集的疊加道集。運(yùn)算符號(hào)*表示褶積。
[0076] 以下結(jié)合圖3和圖4對(duì)本實(shí)施例中等效偏移距偏移的效果進(jìn)行詳細(xì)說(shuō)明。圖3是 CSP道集的克?；舴蚱圃硎疽鈭D。
[0077] 如圖3所示，選擇CSP道集中的第j道（j= 1，2，…N)作為零偏移距道，設(shè)定一 個(gè)雙程旅行時(shí)間t。。A點(diǎn)為所在CSP道集第j道上一個(gè)目標(biāo)散射點(diǎn)，圖中0A距離為該目標(biāo) 散射點(diǎn)的零偏移距旅行時(shí)間t。。第j道到其他每一道的水平距離為δhyN表示CSP道集 的覆蓋次數(shù)。v"s為目標(biāo)散射點(diǎn)A處的均方根速度，第i道在水平面的投影B為走時(shí)t的 采樣點(diǎn)的等效點(diǎn)。根據(jù)表達(dá)式（1)計(jì)算水平距離Sh(圖3中的0B)對(duì)應(yīng)的走時(shí)h。
[0078] 在圖3中，Θ是〇A與AB的夾角，也就是目標(biāo)散射點(diǎn)與采樣點(diǎn)的等效點(diǎn)的連線與 垂直方向的夾角。r表示AB的距離?？捎杀磉_(dá)式（2)計(jì)算目標(biāo)散射點(diǎn)A的波場(chǎng)數(shù)值，進(jìn)而 根據(jù)表達(dá)式（3)計(jì)算CSP道集的等效偏移距偏移道集。
[0079] 圖3中顯示為虛線的雙曲線表示走時(shí)h的軌跡，代表CSP道集進(jìn)行克?；舴蚱?時(shí)的零偏移距動(dòng)校正軌跡。
[0080] 圖4a是一個(gè)傾斜界面的CSP道集進(jìn)行克?；舴蚱魄暗某上裥Ч疽鈭D，其中傾 斜界面的圖像并不收斂，在傾斜界面的下方還呈現(xiàn)一條由近散射點(diǎn)的能量形成的雙曲線。 圖4b為偏移后的成像示意圖，可以看出經(jīng)過(guò)克?；舴蚱坪?，CSP道集上的每個(gè)散射點(diǎn)能 量都會(huì)被收斂到相應(yīng)的一點(diǎn)，最終形成一個(gè)清晰的傾斜界面。
[0081] 以下根據(jù)圖5和圖6說(shuō)明在步驟S103中如何從等效偏移道集中抽取與目標(biāo)散射 點(diǎn)的水平面投影坐標(biāo)相同的道形成共等效點(diǎn)道集。
[0082] 圖5是一個(gè)傾斜界面上一個(gè)點(diǎn)在相鄰多個(gè)共散射點(diǎn)道集的成像效果示意圖。如圖 5所示，在傾斜界面上的A點(diǎn)經(jīng)過(guò)偏移后在CSP1道集、CSP2道集、CSP3道集上都有成像點(diǎn)， 并且A點(diǎn)在CSP1道集、CSP2道集、CSP3道集上的成像點(diǎn)的