[發(fā)明專利]一種復(fù)雜條件下基于角度域共成像點(diǎn)道集的層析速度反演方法

無效

• 申請?zhí)枺?/td>: 201210326923.7
• 申請日: 2012-09-06
• 公開（公告）號: CN102841375A
• 公開（公告）日: 2012-12-26
• 發(fā)明（設(shè)計(jì)）人: 李振春;張凱;秦寧
• 申請（專利權(quán)）人: 中國石油大學(xué)（華東）
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 山東濟(jì)南齊魯科技專利事務(wù)所有限公司 37108
• 代理人: 宋永麗

1.一種復(fù)雜條件下基于角度域共成像點(diǎn)道集的層析速度反演方法，其特征在于反演步驟為：

(a)由常規(guī)速度分析得到的速度模型作為初始速度模型，基于雙平方根延拓疊前深度偏移方法提取角度域共成像點(diǎn)道集；

(b)在提取的角度域共成像點(diǎn)道集拾取剩余深度差，求取旅行時殘差；

(c)在偏移剖面上拾取層位界面，建立射線追蹤正演所需要的參考速度模型，然后利用起伏地表射線追蹤方法求取靈敏度矩陣；

(d)將旅行時殘差、靈敏度矩陣，建立加入正則化因子的反演方程組；利用LSQR方法求解方程組并更新速度模型；

(e)更新得到的速度模型用于雙平方根延拓疊前深度偏移方法和提取角度域共成像點(diǎn)道集；根據(jù)角度域共成像點(diǎn)道集是否拉平以及對速度的精度要求，判斷是否還需再次迭代；如果還未拉平，則重復(fù)上面的流程，直到滿足精度要求；否則進(jìn)行下一層的速度反演；由上到下逐層進(jìn)行分析，直到所有層分析完畢。

2.根據(jù)權(quán)利要求1所述的一種復(fù)雜條件下基于角度域共成像點(diǎn)道集的層析速度反演方法，其特征在于：所述步驟（a）中基于雙平方根疊前深度偏移提取角度域共成像點(diǎn)道集的計(jì)算公式為：

對于常速介質(zhì)情況，地層傾角表示為????????????????????????????????????????????????，入射角，是介質(zhì)速度，炮檢半偏移距，波場延拓過程中表示成像點(diǎn)附近的局部半偏移距，波場傳播總的旅行時，偏移深度，旅行時表示為炮點(diǎn)、接收點(diǎn)、偏移深度的函數(shù)：，這里用中心點(diǎn)和半偏移距來表示炮點(diǎn)和接收點(diǎn)坐標(biāo)，

，????????????????????????（1）

由Snell’s定律，上下行射線參數(shù)表示為：

，?????????????（2）

旅行時相對于深度的導(dǎo)數(shù)用入射射線和反射射線表示為：

???????????????????????（3）

將方程（3）到中心點(diǎn)-半偏移距域，局部偏移距射線參數(shù)表示為：

????????（4）

???????????????????????????（5）

根據(jù)隱函數(shù)求導(dǎo)法則，可以把關(guān)于的導(dǎo)數(shù)變換到關(guān)于：

????????????????????????（6）

在頻率-波數(shù)域，方程（6）式表示為：

??????????????????????????????（7）

其中，分別是局部偏移距波數(shù)和深度波數(shù)；方程（7）是由成像結(jié)果提取角度域共成像點(diǎn)道集的映射公式，描述的是成像點(diǎn)局部入射角度、偏移深度與局部偏移距之間的關(guān)系；此時，是波場延拓過程中的成像點(diǎn)局部半偏移距，是一個成像參數(shù)，不再表示地表炮檢距；該式是假設(shè)常速介質(zhì)情況下推導(dǎo)出來的，只要波場向下延拓時在局部成像點(diǎn)附近滿足的條件，上式仍然成立。

3.根據(jù)權(quán)利要求1所述的一種復(fù)雜條件下基于角度域共成像點(diǎn)道集的層析速度反演方法，其特征在于：所述步驟（a）中基于雙平方根疊前深度偏移提取角度域共成像點(diǎn)道集的方法是：

1）雙平方根方程疊前深度偏移波場延拓過程；

2）部分成像得到炮檢距波數(shù)域的像；

3）相對于的FFT變換到域；

4）根據(jù)（7）式將成像值映射到域；最后反FFT回到域，得到角度域共成像點(diǎn)道集。