[發(fā)明專利]一種基于薄層反射理論的AVA反演方法

有效

• 申請?zhí)枺?/td>: 201910547674.6
• 申請日: 2019-06-24
• 公開（公告）號: CN110187389B
• 公開（公告）日: 2020-07-24
• 發(fā)明（設(shè)計(jì)）人: 楊春;王赟
• 申請（專利權(quán)）人: 中國地質(zhì)大學(xué)（北京）
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京國昊天誠知識(shí)產(chǎn)權(quán)代理有限公司 11315
• 代理人: 劉昕;南霆

1.一種基于薄層反射理論的AVA反演方法，其特征在于，包括如下步驟：

S1、根據(jù)實(shí)際地震記錄，獲得時(shí)間域的薄層實(shí)AVA道集及地震子波；

S2、利用Hilbert變換，由薄層實(shí)AVA道集求得薄層虛AVA道集，并結(jié)合薄層實(shí)AVA道集和薄層虛AVA道構(gòu)造時(shí)間域的薄層復(fù)AVA道集；

S3、針對時(shí)間域的薄層復(fù)AVA道集及地震子波，利用FFT變換，分別獲得頻率域的薄層復(fù)AVA道集及地震子波，兩者通過除法運(yùn)算，即可獲得頻率域中實(shí)際的薄層PP波復(fù)反射系數(shù)矩陣R^obs；

S4、建立初始的薄層模型m；所述薄層模型m包含薄層厚度以及薄層三層介質(zhì)的縱波速度、橫波速度、密度；

S5、利用薄層PP波反射系數(shù)近似式獲得對應(yīng)薄層模型m的PP波反射系數(shù)矩陣R；

S6、求解目標(biāo)函數(shù)Q(m)＝‖R^obs-R‖²,并判斷Q(m)是否滿足精度要求；

S7、當(dāng)Q(m)小于閾值ε時(shí)，Q(m)滿足精度要求，輸出反演結(jié)果m；當(dāng)Q(m)不小于ε時(shí)，Q(m)不滿足精度要求，利用高斯-牛頓算法，計(jì)算更新量Δm，令m＝m+Δm作為新的薄層模型返回步驟S5；

其中，

步驟S2中，常規(guī)薄層AVA地震數(shù)據(jù)s^real(θ,t)為薄層復(fù)AVA道集的實(shí)部，利用Hilbert變換獲得相應(yīng)的薄層虛AVA道集如下：

t為時(shí)間，θ為入射角；

由此，時(shí)間域的薄層復(fù)AVA道集為：

s(θ,t)＝s^real(θ,t)+js^im(θ,t) (2)。

2.根據(jù)權(quán)利要求1所述的基于薄層反射理論的AVA反演方法，其特征在于，步驟S3中，

利用FFT變換，可以獲得頻率域中薄層復(fù)AVA道集：

S(θ,f)＝FFT[s(θ,t)] (3)；

s(θ,t)為時(shí)間域的薄層復(fù)AVA道集，f為頻率；

利用FFT變換，同樣可以獲得頻率域的地震子波：

W(θ,f)＝FFT[w(θ,t)] (4)；

[w(θ,t)]為時(shí)間域的地震子波；

由此，獲得頻率域中實(shí)際的薄層PP波復(fù)反射系數(shù)矩陣R^obs：

R^obs(θ,f)＝S(θ,f)/W(θ,f) (5)；

其中R^obs(θ,f)為不同頻率處實(shí)際的薄層PP波復(fù)反射系數(shù)矩陣。