[發(fā)明專利]大規(guī)模并行的克希霍夫疊前深度偏移方法及裝置

有效

• 申請(qǐng)?zhí)枺?/td>: 201410367547.5
• 申請(qǐng)日: 2014-07-29
• 公開（公告）號(hào): CN104133240B
• 公開（公告）日: 2017-02-01
• 發(fā)明（設(shè)計(jì)）人: 張建磊;李超;趙長(zhǎng)海;王成祥;王獅虎;晏海華
• 申請(qǐng)（專利權(quán)）人: 中國石油天然氣集團(tuán)公司;中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28;G01V1/30;G01V1/34
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司11127
• 代理人: 王天堯

技術(shù)領(lǐng)域

本發(fā)明涉及石油地球物理勘探地震數(shù)據(jù)的處理技術(shù)領(lǐng)域，特別涉及一種大規(guī)模并行的克希霍夫疊前深度偏移方法及裝置。

背景技術(shù)

克希霍夫疊前深度偏移技術(shù)是地球物理勘探地震數(shù)據(jù)處理過程中最重要的一種深度域成像方法，該方法相對(duì)于時(shí)間偏移具有更高的成像質(zhì)量，而且具有不受野外數(shù)據(jù)觀測(cè)方式的限制等優(yōu)點(diǎn)。另外它還能輸出基于炮檢距的共成像點(diǎn)道集，這類道集又是進(jìn)一步速度分析的重要信息。

在水平地表情況下，Kirchhoff積分法偏移公式為：

I(ξ)=∫ΩξW(ξ,m,h)D[t=tD(ξ,m,h),m,h]dmdh---(1)]]>

其中，成像點(diǎn)ξ＝(x_ξ,y_ξ,z_ξ)，對(duì)于深度偏移為ξ＝(x_ξ,y_ξ,z_ξ)，I(ξ)是成像結(jié)果，D[t,m,h]是野外觀測(cè)的地震數(shù)據(jù)，m是共中心點(diǎn)，h是半炮檢距，Ω_ξ是偏移孔徑。疊前偏移過程就是對(duì)一系列觀測(cè)數(shù)據(jù)進(jìn)行加權(quán)求和，W(ξ,m,h)是加權(quán)因子，t_D(ξ,m,h)是由炮點(diǎn)到成像點(diǎn)再到接收點(diǎn)的旅行時(shí)間。

克希霍夫疊前深度偏移方法的數(shù)學(xué)表達(dá)雖然較為抽象，但其物理實(shí)現(xiàn)過程可簡(jiǎn)單描述為：地震道是從待偏移數(shù)據(jù)空間到偏移結(jié)果數(shù)據(jù)空間多對(duì)多的映射，映射關(guān)系為程函方程旅行時(shí)計(jì)算公式，在不考慮時(shí)變孔徑的條件下，克希霍夫疊前深度偏移就是把每一道待偏移地震數(shù)據(jù)(炮檢距為h)映射到偏移結(jié)果數(shù)據(jù)集中的一個(gè)橢圓柱體中，橢圓柱體的軸心是當(dāng)前待偏移地震道按大地坐標(biāo)在偏移結(jié)果數(shù)據(jù)空間中的投影位置。若偏移孔徑采用圓形且孔徑值為300，則一個(gè)待偏移地震采樣點(diǎn)將要映射到282600(3.14*300*300)個(gè)偏移結(jié)果地震采樣點(diǎn)上。

在水平地表情況下，克希霍夫積分法偏移用公式為前面所述公式(1)。(1)可以進(jìn)行深度偏移，當(dāng)進(jìn)行深度偏移時(shí)采用圖1所示的方式計(jì)算旅行時(shí)。其計(jì)算公式采用(2)式和(3)式聯(lián)合計(jì)算。

T＝T_s+T_r??????(2)