本發(fā)明一實施例中的縱波終疊剖面和轉(zhuǎn)換波終疊剖面；
[0041] 圖5顯示了圖4經(jīng)過層位自動匹配處理的縱波終疊剖面和轉(zhuǎn)換波終疊剖面W及相 應(yīng)的縱橫波速度比剖面；
[0042] 圖6顯示了本發(fā)明另一實施例中某工區(qū)L385線通過速度交互分析獲得的縱波偏 移速度模型；
[0043] 圖7顯示了L385線縱橫波速度比剖面；
[0044]圖8顯示了根據(jù)圖7所示的縱橫波速度比剖面換算出的橫波偏移速度模型；
[0045] 圖9顯示了L385線轉(zhuǎn)換波層位一致性疊前時間偏移剖面和縱波疊前時間偏移剖 面的比較結(jié)果；
[0046] 圖10顯示了通過不同的方法獲得的兩種轉(zhuǎn)換波疊前時間偏移剖面的對比圖。
【具體實施方式】
[0047]為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚，W下結(jié)合附圖對本發(fā)明作進(jìn)一步 地詳細(xì)說明。圖3顯示了本發(fā)明的橫波偏移速度建模方法的一種實施方案的流程圖。需要 說明的是，圖3顯示的流程圖僅用于解釋本發(fā)明，但不構(gòu)成對本發(fā)明的限制。
[0048]S100,分別對縱波地震資料和轉(zhuǎn)換波地震資料進(jìn)行疊加處理，獲取縱波和轉(zhuǎn)換波 的等效速度終疊剖面，并建立縱波偏移速度模型W及獲取縱波疊前時間偏移剖面。
[0049]S200,對縱波終疊剖面和轉(zhuǎn)換波終疊剖面進(jìn)行層位匹配處理，求取兩個終疊剖面 之間的時移序列，根據(jù)兩個終疊剖面之間的時移序列反演縱橫波速度比。
[0050] 縱波終疊剖面和轉(zhuǎn)換波終疊剖面可W視為兩組二維圖像。在本實施例中，優(yōu)選地 根據(jù)動態(tài)圖像變形DIW(DynamicImageWa巧ing)最優(yōu)化技術(shù)，W送兩組二維圖像最小平方 誤差為目標(biāo)，求取兩個終疊剖面之間的時移序列。目P，按照下式求解兩個時間序列的圖像 f[i]與g[i]之間的時移序列u[i]的最優(yōu)化問題：
[0051]
[0052] 式（1)中定義的函數(shù)D表示一種距離函數(shù)，是與圖像誤差函數(shù)e有關(guān)的累計誤差 函數(shù)。
[0053]
(2)
[0054] 式（2)中l(wèi)[i]為延遲量，表示每個樣點上的時移大小。
[0055] 時移序列U山滿足W下約束條件：
[005引 |u山I《1 做
[0057] 而圖像的誤差函數(shù)e可表述為兩組圖像差的平方：
[0058] e[i，1]三（f[i]-g[i+l])2 (4)
[0059] 根據(jù)求取的時移序列，就可W完成縱波終疊剖面和轉(zhuǎn)換波終疊剖面的反射波組的 同相軸的層位匹配。也即將來自地下同一反射層位的縱波和轉(zhuǎn)換波的同相軸匹配在縱波時 間上，注冊為同層反射波組。
[0060] 然后根據(jù)兩個終疊剖面之間的時移序列，通過下式求取縱橫波速度比：
[006。
（5)
[0062] 式巧）中Y表示縱橫波速度比，Vp，Vg分別表示縱波速度和橫波速度，tpp表示縱 波雙程旅行時，化/dtpp表示時移函數(shù)U對縱波雙程旅行時的偏導(dǎo)數(shù)。式（5)描述了時移序 列對縱波時間偏導(dǎo)數(shù)與縱橫波速度比之間的關(guān)系。其中，偏導(dǎo)數(shù)對應(yīng)時間應(yīng)變極限，C是比 例因子，取值范圍大約在2. 0左右。例如當(dāng)C取2. 0,時間應(yīng)變極限為10%。
[0063] 優(yōu)選地，在本實施例中，還W縱波速度場的時間-速度對（已知）為控制點，對上 述縱橫波速度比做層內(nèi)的趨勢平滑處理，W獲得每層穩(wěn)定變化的縱橫波速度比。
[0064] S300,根據(jù)縱波偏移速度模型和縱橫波速度比，建立橫波偏移速度模型。
[0065] 根據(jù)已知的縱波偏移速度模型和平滑處理后的縱橫波速度比，通過下式換算出橫 波偏移速度，建立橫波偏移速度模型：
[0066]
( 6)
[0067] S400,根據(jù)縱波偏移速度模型和橫波偏移速度模型，對轉(zhuǎn)換波地震資料進(jìn)行疊前 時間偏移處理，獲取其時間與縱波旅行時一致的轉(zhuǎn)換波疊前時間偏移剖面。
[0068] S500,對縱波疊前時間偏移剖面和轉(zhuǎn)換波疊前時間偏移剖面再次進(jìn)行層位匹配處 理，求取兩個偏移剖面之間的時移序列。
[0069] 與步驟S200類似，在層位匹配處理過程中，也采用動態(tài)圖像變形DIW值ynamiC ImageWarping)最優(yōu)化技術(shù)，計算縱波疊前時間偏移剖面與轉(zhuǎn)換波疊前時間偏移剖面之間 的時移序列，完成兩個偏移剖面的層位自動匹配。
[0070]S600,判斷兩個偏移剖面之間的時移序列是否大于給定的精度闊值：
[0071] 如果是，進(jìn)入步驟S700;
[007引如果否，進(jìn)入步驟S900;
[0073]S700,根據(jù)兩個偏移剖面之間的時移序列反演縱橫波速度比微量變化，根據(jù)縱橫 波速度比微量變化修正縱橫波速度比。
[0074] 首先根據(jù)縱波疊前時間偏移剖面與轉(zhuǎn)換波疊前時間偏移剖面之間的時移序列，通 過下式求取縱橫波速度比微量變化：
[0075]
( 7)
[0076]式中，ΔY(tPP)表示縱橫波速度比微量變化，化'Mtpp表示兩個偏移剖面之間 時移函數(shù)U'對縱波雙程旅行時tpp的偏導(dǎo)數(shù)，C'是比例因子，取值范圍大約在1.0左右。
[0077] 優(yōu)選地，在本實施例中，還W縱波速度場的時間-速度對（已知）為控制點，對縱 橫波速度比微量變化做層內(nèi)的趨勢平滑處理，獲得每層穩(wěn)定變化的縱橫波速度比微量變 化。
[0078] 最后根據(jù)平滑處理后的縱橫波速度比微量變化，通過下式修正縱橫波速度比：
[0079] Y(tpp) =Y(tpp) ±ΔY(tpp) 做
[0080]S800,根據(jù)縱波偏移速度模型和修正后的縱橫波速度比，更新橫波偏移速度模型， 返回步驟S400。
[0081]S900,輸出橫波偏移速度模型。
[008引理論上，上述實施例中，經(jīng)過步驟S400獲得的轉(zhuǎn)換波疊前時間偏移剖面的時間應(yīng) 當(dāng)與縱波疊前時間偏移剖面的時間完全一致。但在實際應(yīng)用時，送兩個偏移剖面之間可能 仍會存在一定的時移序列，也即剩余時移序列。送種剩余時移序列是由于建立的橫波偏移 速度模型不合適而產(chǎn)生的。因此本發(fā)明提出在步驟S300建立橫波偏移速度模型之后，進(jìn)一 步執(zhí)行步驟S400~步驟S800,W通過偏移剖面的層位匹配來修正不合適的橫波偏移速度 模型，使其達(dá)到成像速度模型所需要的精度。具體來說，應(yīng)用剩余時移序列反演縱橫波速度 比的微量變化，根據(jù)縱橫波速度比的微量變化修正縱橫波速度比，然后根據(jù)修正的縱橫波 速度比更新橫波偏移速度模型。
[0083] -般而言，通過一次迭代修正過程就可W使橫波偏移速度模型達(dá)到成像速度模型 所需要的精度，當(dāng)然也不排除其他的情況。例如不進(jìn)行迭代修正，直接輸出步驟S400獲得 的橫波偏移速度模型。又或者反復(fù)執(zhí)行步驟S400~步驟S800,進(jìn)行多次迭代，直至橫波偏 移速度模型滿足給定的精度要求。
[0084] 實施例一
[0085] 圖4顯示了一實施例中的縱波終疊剖面和轉(zhuǎn)換波終疊剖面。其中，圖4(a)為縱波 終疊剖面，記錄時長為5砂，圖4(b)為轉(zhuǎn)換波終疊剖面，記錄時長為8砂。
[0086] 圖5顯示了經(jīng)過層位自動匹配處理的縱波終疊剖面和轉(zhuǎn)換波終疊剖面，W及根據(jù) 時移量獲得的縱橫波速度比剖面。其中，圖5 (a)為縱波終疊剖面，記錄時長為5砂，圖5化） 為壓縮后的轉(zhuǎn)換波終疊剖面，記錄時長為5砂，圖5(c)為根據(jù)時移量獲得的縱橫波速度比 剖面，記錄時長為5砂。
[0087] 實施例二
[0088] 圖6顯示了某工區(qū)L385線通過速度交