[發(fā)明專利]一種大斜井炮檢移動VSP自由表面多次波成像方法和裝置

有效

• 申請?zhí)枺?/td>: 202010670672.9
• 申請日: 2020-07-13
• 公開（公告）號: CN111781647B
• 公開（公告）日: 2022-05-20
• 發(fā)明（設計）人: 李建國
• 申請（專利權）人: 中油奧博（成都）科技有限公司
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 成都方圓聿聯(lián)專利代理事務所(普通合伙) 51241
• 代理人: 李鵬

1.一種大斜井炮檢移動VSP自由表面多次波成像方法，其特征在于：包括以下步驟：

S1.輸入大斜井炮檢移動VSP全波場共檢波點道集、2維地質(zhì)模型、偏移參數(shù)；

S2.根據(jù)輸入數(shù)據(jù)，建立炮點、檢波點在2維地質(zhì)模型中的坐標；

S3.對輸入的共檢波點道集做快速傅里葉變換到頻率波數(shù)域，設置頻率域震源；

S4.將頻率域震源從檢波點深度向上延拓至自由表面，頻率波數(shù)域炮檢移動VSP共檢波點全波場、震源波場從自由表面向下延拓，實現(xiàn)單檢波點的自由表面多次波單程波疊前深度偏移成像；

所述步驟S4包括：

S401.計算頻率域震源向上延拓網(wǎng)格：

kzmig＝(k-1)·dzmig

其中，MRZ_j是第j個檢波點的模型深度，dzmig是成像深度步長，kzmig是第k個向上延拓的深度；

S402.頻率域震源向上延拓一個深度步長：

f＝[f_l:df:f_h]

kk＝f/V

g_s0＝1

g_s1＝-i·π·dzmig/kk

g_s2＝(-i·π·dzmig·kk/4-π²·dzmig²·kk²/2)/kk⁴

g_s3＝(-i·π·dzmig·kk/8+π²·dzmig²·kk²/4-i·π³·dzmig³·kk³/6)/kk⁶

g_s4＝(-i·π·dzmig·kk·5/64-π²·dzmig²·kk²·5/32+i·π³·dzmig³·kk³/8+π⁴·dzmig⁴·kk⁴/24)/kk⁸

其中，dip是最大傾角，fl是最小頻率，fh是最大頻率，f＝[f_l:df:f_h]是頻率向量，表示頻率從f_l到f_h間隔為df，V是速度，k_x是波數(shù)，dzmig是成像深度步長，i是虛數(shù)單位，fft是傅里葉變換函數(shù)，ifft是傅里葉逆變換函數(shù)，k_m是最大傾角dip對應的最大波數(shù)，mutes是波數(shù)切除因子，kk、g_s0-g_s4、t_s0-t_s4是中間變量，是向上延拓一個深度步長的波場；

S403.循環(huán)執(zhí)行步驟S401～S402實現(xiàn)頻率域震源向上延拓至自由表面；

S404.計算向下延拓成像網(wǎng)格：

kzmig＝(k-1)·dzmig

其中，zmig2是成像終止深度，dzmig是成像深度步長，kzmig是第k個成像深度；

S405.將頻率波數(shù)域炮檢移動VSP共檢波點全波場向下延拓一個深度步長：

f＝[f_l:df:f_h]

kk＝f/V

g₀＝1

g₁＝-i·π·dzmig/kk

g₂＝(-i·π·dzmig·kk/4-π²·dzmig²·kk²/2)/kk⁴

g₃＝(-i·π·dzmig·kk/8+π²·dzmig²·kk²/4-i·π³·dzmig³·kk³/6)/kk⁶

g₄＝(-i·π·dzmig·kk·5/64-π²·dzmig²·kk²·5/32+i·π³·dzmig³·kk³/8+π⁴·dzmig⁴·kk⁴/24)/kk⁸

其中，dip是最大傾角，fl是最小頻率，fh是最大頻率，f＝[f_l:df:f_h]是頻率向量，V是速度，k_x是波數(shù)，dzmig是成像深度步長，i是虛數(shù)單位，fft是傅里葉變換函數(shù)，ifft是傅里葉逆變換函數(shù)，k_m是最大傾角dip對應的最大波數(shù)，mute是波數(shù)切除因子，kk、g₀-g₄、t₀-t₄是中間變量，是向下延拓一個深度步長的波場；

S406.將頻率域震源向下延拓一個深度步長：

f＝[f_l:df:f_h]

kk＝f/V

g_s0＝1

g_s1＝i·π·dzmig/kk

g_s2＝(i·π·dzmig·kk/4-π²·dzmig²·kk²/2)/kk⁴

g_s3＝(i·π·dzmig·kk/8+π²·dzmig²·kk²/4+i·π³·dzmig³·kk³/6)/kk⁶

g_s4＝(i·π·dzmig·kk·5/64-π²·dzmig²·kk²·5/32-i·π³·dzmig³·kk³/8+π⁴·dzmig⁴·kk⁴/24)/kk⁸

其中，dip是最大傾角，fl是最小頻率，fh是最大頻率，f＝[f_l:df:f_h]是頻率向量，V是速度，k_x是波數(shù)，dzmig是成像深度步長，i是虛數(shù)單位，fft是傅里葉變換函數(shù)，ifft是傅里葉逆變換函數(shù)，是向下延拓一個深度步長的波場；

S407.相關成像條件提取成像值：

其中，是向下延拓一個深度步長的波場，是向下延拓一個深度步長的波場，是提取的kzmig+dzmig深度的成像值，conj是復共軛函數(shù)；

S408.循環(huán)執(zhí)行步驟S404～S407，直至延拓到成像終止深度zmig2，完成第j個檢波點的炮檢移動VSP自由表面多次波疊前深度偏移成像；

S5.循環(huán)執(zhí)行步驟S4，完成所有檢波點的炮檢移動VSP自由表面多次波疊前深度偏移成像；

S6.將步驟S5的成像重排成共成像道集，共成像道集疊加，得到炮檢移動VSP自由表面多次波疊前深度偏移疊加成像。