[發(fā)明專利]基于克希霍夫積分法的繞射波場(chǎng)分離方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201010527864.0
• 申請(qǐng)日: 2010-11-02
• 公開（公告）號(hào): CN102455439A
• 公開（公告）日: 2012-05-16
• 發(fā)明（設(shè)計(jì)）人: 白英哲;孫贊東
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油大學(xué)(北京)
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 暫無(wú)信息
• 代理人: 暫無(wú)信息

基于克希霍夫積分法的繞射波場(chǎng)分離方法

1.本發(fā)明意在提供一種通過(guò)高信噪比地震資料，通過(guò)克希霍夫偏移輸出傾角道集，根據(jù)在傾角道集中繞射波和反射波的明顯區(qū)別來(lái)獲取繞射波的單獨(dú)成像結(jié)果，從而對(duì)地下的地質(zhì)突變點(diǎn)進(jìn)行精細(xì)成像，以便于更好的識(shí)別碳酸鹽巖洞縫儲(chǔ)層。

本專利使用處理過(guò)的高信噪比地震資料以及高精度速度分析得到的層速度場(chǎng)實(shí)現(xiàn)上述目的采取的技術(shù)方案如下：

步驟1：使用輸入的地層速度模型進(jìn)行射線追蹤，得到地面各炮檢點(diǎn)的旅行時(shí)列表。射線追蹤可以選擇Vadile(1988)提出的快速有限差分射線追蹤法。

步驟2：克希霍夫法輸出傾角道集。

步驟3：在步驟2獲得的傾角道集中，反射波表現(xiàn)為一條向上開口的曲線，在繞射點(diǎn)正上方觀測(cè)的繞射波表現(xiàn)為一條水平直線。因此，可以通過(guò)中值濾波的方法將兩者分離。

步驟4：將步驟3得到的濾波輸出的傾角道集疊加，就可以生成繞射波的單獨(dú)成像剖面。

2.根據(jù)權(quán)利要求1所述的基于克希霍夫積分法的繞射波場(chǎng)分離方法，其特征在于，所述的步驟2中的輸出傾角道集的方法為：

首先，讀入一道輸入地震數(shù)據(jù)。并分別讀入步驟1得到的這一道對(duì)應(yīng)炮檢對(duì)的旅行時(shí)列表；

其次，通過(guò)以下公式計(jì)算地下各點(diǎn)處入射線和反射線同豎直方向的夾角：

tanαs=∂τs/∂x∂τs/∂z]]>

tanαr=∂τr/∂x∂τr/∂z]]>

式中：τ_s為炮點(diǎn)旅行時(shí)，τ_r為檢波點(diǎn)旅行時(shí)，α_s為入射線與豎直方向的夾角α_r為反射線與豎直方向的夾角，(x，z)為空間坐標(biāo)。

根據(jù)幾何關(guān)系，可以求出傾角θ、反射角與α_s、α_r之間的關(guān)系

θ=αs+αr2]]>

由以上分析可知通過(guò)射線追蹤，可以建立起炮點(diǎn)坐標(biāo)和計(jì)算點(diǎn)坐標(biāo)以及傾角θ和檢波點(diǎn)坐標(biāo)r的關(guān)系

s＝s(x，z，r，θ)

進(jìn)而克希霍夫偏移公式可以寫成

R(x,z,θ)=∫Γu(s(x,z,r,θ),r,τs+τr)dr]]>

這樣就實(shí)現(xiàn)了Kirchhoff偏移輸出傾角道集。