兩界面之間的多次波的疊后預(yù)測方法和裝置的制造方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地球物理勘探技術(shù)領(lǐng)域，特別涉及一種兩界面之間的多次波的疊后預(yù) 測方法和裝置。
【背景技術(shù)】
[0002] 在一些特殊的巖體與圍巖之間的界面往往會產(chǎn)生較強的多次波或?qū)娱g多次波，例 如：碳酸鹽巖、火成巖和煤層等發(fā)育的探區(qū)就很容易產(chǎn)生較強的多次波或?qū)娱g多次波，這會 對地震資料的速度分析、成像和反演產(chǎn)生比較大的影響。
[0003] 為了減少這種影響，一般需要對多次波進行壓制，目前對于信噪比較高的地震資 料，壓制多次波一般采用疊前方法，而對于信噪比低的地震資料，壓制多次波一般采用疊后 方法。進一步的，考慮到海上地震資料信噪比較高，陸上低信噪比地震資料信噪比低，因此 壓制海上地震資料的多次波一般采用疊前方法，壓制陸上地震資料多次波一般采用疊后方 法，下面對疊后方法進行具體描述：
[0004] 疊后壓制多次波的方法主要采用波場延拓的方法，即給定地下一個或多個產(chǎn)生多 次波的界面及對應(yīng)的疊加速度，然后通過一次下行波波場延拓和一次上行波波場延拓來預(yù) 測與該界面有關(guān)的全程多次波和微屈多次波。對于兩個界面之間的層間多次波預(yù)測，理論 上可以通過六次波場延拓來實現(xiàn)，然而，這種實現(xiàn)方式的缺點是需要兩個界面之間的層速 度，而由疊加速度轉(zhuǎn)化為層速度的過程會使得層速度的結(jié)果存在較大的誤差，且，較多次的 波場延拓會使得計算效率較低。
[0005] 針對如何高效準確地預(yù)測兩個界面之間的層間多次波的問題，目前尚未提出有效 的解決方案。

【發(fā)明內(nèi)容】

[0006] 本發(fā)明實施例提供了一種兩界面之間的多次波的疊后預(yù)測方法，以解決現(xiàn)有技術(shù) 中兩界面的層間多次波的預(yù)測的準確性不高，預(yù)測效率低下的技術(shù)問題，該方法包括：
[0007] 獲取疊加剖面和與所述疊加剖面對應(yīng)的疊加速度場；
[0008] 獲取給定的第一時深和第二時深，其中，第二時深大于第一時深，第一時深對應(yīng)的 地下反射界面為第一界面，第二時深對應(yīng)的地下反射界面為第二界面；
[0009] 以所述第一界面將所述疊加剖面分為第一部分數(shù)據(jù)和第二部分數(shù)據(jù)，其中，所述 第一部分數(shù)據(jù)為所述第一界面以上部分的數(shù)據(jù)，所述第二部分數(shù)據(jù)為所述第一界面以下部 分的數(shù)據(jù)；
[0010] 從所述疊加速度場中查找出所述第一時深對應(yīng)的第一疊加速度和所述第二時深 對應(yīng)的第二疊加速度；
[0011] 根據(jù)所述第一時深和所述第一疊加速度計算得到所述第一界面的第一深度，根據(jù) 所述第二時深和所述第二疊加速度計算得到所述第二界面的第二深度；
[0012] 根據(jù)所述第一深度和所述第一疊加速度，采用下行波方程將所述疊加剖面從基準 面延拓到所述第一界面，得到第一波場，再采用上行波方程，將所述第一波場從所述第一界 面延拓到所述基準面，得到所述第一界面的全程和微屈多次波；
[0013] 根據(jù)所述第二深度和所述第二疊加速度，采用下行波方程將所述疊加剖面從基準 面延拓到所述第二界面，得到第二波場，再采用上行波方程，將所述第二波場從所述第二界 面延拓到所述基準面，得到所述第二界面的全程和微屈多次波；
[0014] 根據(jù)所述第一界面的全程和微屈多次波計算得到所述第一界面多次波的地震道 的能量，根據(jù)所述第二界面的全程和微屈多次波計算得到所述第二界面多次波的地震道的 能量；
[0015] 根據(jù)所述第一界面的全程和微屈多次波和所述第二界面的全程和微屈多次波，計 算得到所述第一界面和所述第二界面的多次波之間的互相關(guān)剖面；
[0016] 根據(jù)所述第一界面多次波的地震道的能量、所述第二界面多次波的地震道的能量 和所述互相關(guān)剖面，計算得到歸一化的互相關(guān)剖面；
[0017] 根據(jù)所述歸一化的相關(guān)剖面和所述第二部分數(shù)據(jù)，計算得到所述第一界面和所述 第二界面之間的層間多次波。
[0018] 在一個實施方式中，按照以下公式計算所述第一深度和第二深度：
[0019] Z1 (x, y) = T1 (x, y) V1 (x, y)
[0020] Z2 (x, y) = T2 (x, y) V2 (x, y)
[0021] 其中，V1U, y) = V(x，y，I\(x，y))，V2(x，y) = V(x，y，T2(x，y))，示第一深度， Z2表示第二深度，T i表示第一時深，T 2表示第二時深，V i表示第一疊加速度，V 2表示第二疊 加速度，V表示所述疊加剖面的疊加速度場，X表示疊加剖面的縱測線方向坐標，y表示疊加 剖面的橫測線方向坐標。
[0022] 在一個實施方式中，按照以下公式計算所述第一界面多次波的地震道的能量和所 述第二界面多次波的地震道的能量：
【主權(quán)項】
1. 一種兩界面之間的多次波的疊后預(yù)測方法，其特征在于，包括： 獲取疊加剖面和與所述疊加剖面對應(yīng)的疊加速度場； 獲取給定的第一時深和第二時深，其中，第二時深大于第一時深，第一時深對應(yīng)的地下 反射界面為第一界面，第二時深對應(yīng)的地下反射界面為第二界面； 以所述第一界面將所述疊加剖面分為第一部分數(shù)據(jù)和第二部分數(shù)據(jù)，其中，所述第一 部分數(shù)據(jù)為所述第一界面以上部分的數(shù)據(jù)，所述第二部分數(shù)據(jù)為所述第一界面以下部分的 數(shù)據(jù)； 從所述疊加速度場中查找出所述第一時深對應(yīng)的第一疊加速度和所述第二時深對應(yīng) 的第二疊加速度； 根據(jù)所述第一時深和所述第一疊加速度計算得到所述第一界面的第一深度，根據(jù)所述 第二時深和所述第二疊加速度計算得到所述第二界面的第二深度； 根據(jù)所述第一深度和所述第一疊加速度，采用下行波方程將所述疊加剖面從基準面延 拓到所述第一界面，得到第一波場，再采用上行波方程，將所述第一波場從所述第一界面延 拓到所述基準面，得到所述第一界面的全程和微屈多次波； 根據(jù)所述第二深度和所述第二疊加速度，采用下行波方程將所述疊加剖面從基準面延 拓到所述第二界面，得到第二波場，再采用上行波方程，將所述第二波場從所述第二界面延 拓到所述基準面，得到所述第二界面的全程和微屈多次波； 根據(jù)所述第一界面的全程和微屈多次波計算得到所述第一界面多次波的地震道的能 量，根據(jù)所述第二界面的全程和微屈多次波計算得到所述第二界面多次波的地震道的能 量； 根據(jù)所述第一界面的全程和微屈多次波和所述第二界面的全程和微屈多次波，計算得 到所述第一界面和所述第二界面的多次波之間的互相關(guān)剖面； 根據(jù)所述第一界面多次波的地震道的能量、所述第二界面多次波的地震道的能量和所 述互相關(guān)剖面，計算得到歸一化的互相關(guān)剖面； 根據(jù)所述歸一化的相關(guān)剖面和所述第二部分數(shù)據(jù)，計算得到所述第一界面和所述第二 界面之間的層間多次波。
2. 如權(quán)利要求1所述的方法，其特征在于，按照以下公式計算所述第一深度和第二深 度： Z1 (x, y) = T^x, y)V1(x, y) Z2 (x, y) = T2 (x, y)V2(x, y) 其中，V1 (x, y) = V (x, y, T1 (x, y))，V2 (x, y) = V (x, y, T2 (x, y)