[發(fā)明專利]一種基于稀疏反演的層間多次波壓制方法

無效

• 申請?zhí)枺?/td>: 201310494680.2
• 申請日: 2013-10-21
• 公開（公告）號: CN103558633A
• 公開（公告）日: 2014-02-05
• 發(fā)明（設(shè)計）人: 趙偉;陳小宏;張云鵬;劉永江;翁斌;王小六;郝振江;糜芳;朱振宇;薛東川;彭利麗;孫文博;李景葉
• 申請（專利權(quán)）人: 中國海洋石油總公司;中海油研究總院;中國石油大學(xué)（北京）
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/36
• 代理公司: 北京紀(jì)凱知識產(chǎn)權(quán)代理有限公司 11245
• 代理人: 徐寧;關(guān)暢

1.一種基于稀疏反演的層間多次波壓制方法，其包括以下步驟：

1）在疊加剖面上選取介于強(qiáng)反射界面之間的基準(zhǔn)面，在炮集上分離來自基準(zhǔn)面上地層的波場數(shù)據(jù)δP和下地層的波場數(shù)據(jù)

2）拾取δP初至一次波振幅峰值對應(yīng)的脈沖估計基準(zhǔn)面以上地層脈沖響應(yīng)的初始值，并通過匹配濾波估計震源子波初始值，并設(shè)定反演脈沖響應(yīng)與震源子波的迭代次數(shù)；

3）迭代更新基準(zhǔn)面以上地層的脈沖響應(yīng)、基準(zhǔn)面以下地層的脈沖響應(yīng)以及震源子波；

4）匹配層間多次波預(yù)測結(jié)果與原始地震數(shù)據(jù)得到振幅校正因子，利用振幅校正因子更新脈沖響應(yīng)和震源子波；

5）以步驟4）的參數(shù)更新結(jié)果作為新的初始值，繼續(xù)重復(fù)步驟3）～4）直至達(dá)到給定迭代次數(shù)，停止迭代；

6）從地震數(shù)據(jù)P中減去預(yù)測得到的層間多次波M得到一次波P₀。

2.如權(quán)利要求1所述的一種基于稀疏反演的層間多次波壓制方法，其特征在于：所述步驟3）迭代更新基準(zhǔn)面以上地層的脈沖響應(yīng)、基準(zhǔn)面以下地層的脈沖響應(yīng)以及震源子波的具體計算過程為：

3.1）更新基準(zhǔn)面以上地層的脈沖響應(yīng)δX：

利用最速下降法得到δX在i+1次迭代時的最速下降方向：

ΔδX=VSiH-P‾VHX‾i]]>

變換ΔδX到時間域，拾取每道振幅峰值對應(yīng)的脈沖后變換至頻率域，得到脈沖響應(yīng)的更新方向ΔδX，乘以不依賴于頻率的步長α，更新脈沖響應(yīng)與擬合殘差：

δX_i+1＝δX_i+αΔδX

V＝V+αΔV₁

其中，

α=-Σωvec(ΔV1)Hvec(V)Σωvec(ΔV1)Hvec(ΔV1)]]>

ΔV1=-ΔδX‾Si+X‾iΔδX‾HP‾;]]>

3.2）更新基準(zhǔn)面以下地層的脈沖響應(yīng)

利用最速下降法計算在i+1次迭代時的最速下降方向：

ΔX‾=VSiH-VP‾HδXi+1]]>

同樣的，變換到時間域，拾取每道最大振幅值對應(yīng)的脈沖后變回頻率域得到的更新方向乘以不依賴于頻率的更新步長β，更新脈沖響應(yīng)與擬合殘差：

X‾i+1=X‾i+βΔX‾‾]]>

V＝V+βΔV₂

其中，

β=-Σωvec(ΔV2)Hvec(V)Σωvec(ΔV2)Hvec(ΔV2)]]>

ΔV2=-ΔX‾‾Si+ΔX‾‾δXi+1HP‾;]]>

3.3）匹配脈沖響應(yīng)與更新震源子波。