[发明专利]一种大斜井炮检移动VSP自由表面多次波成像方法和装置有效

申请号：	202010670672.9	申请日：	2020-07-13
公开（公告）号：	CN111781647B	公开（公告）日：	2022-05-20
发明（设计）人：	李建国	申请（专利权）人：	中油奥博（成都）科技有限公司
主分类号：	G01V1/36	分类号：	G01V1/36
代理公司：	成都方圆聿联专利代理事务所(普通合伙) 51241	代理人：	李鹏
地址：	611730 四川省成都市***	国省代码：	四川;51
权利要求书：	查看更多	说明书：	查看更多
摘要：
搜索关键词：	一种斜井移动 vsp 自由表面多次成像方法装置
钻瓜网技术展会专利词库专利权人专利榜在售专利公布日期热门专利

【权利要求书】：

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

S1.输入大斜井炮检移动VSP全波场共检波点道集、2维地质模型、偏移参数；

S2.根据输入数据，建立炮点、检波点在2维地质模型中的坐标；

S3.对输入的共检波点道集做快速傅里叶变换到频率波数域，设置频率域震源；

S4.将频率域震源从检波点深度向上延拓至自由表面，频率波数域炮检移动VSP共检波点全波场、震源波场从自由表面向下延拓，实现单检波点的自由表面多次波单程波叠前深度偏移成像；

所述步骤S4包括：

S401.计算频率域震源向上延拓网格：

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是波数，dzmig是成像深度步长，i是虚数单位，fft是傅里叶变换函数，ifft是傅里叶逆变换函数，k_m是最大倾角dip对应的最大波数，mutes是波数切除因子，kk、g_s0-g_s4、t_s0-t_s4是中间变量，是向上延拓一个深度步长的波场；

S403.循环执行步骤S401～S402实现频率域震源向上延拓至自由表面；

S404.计算向下延拓成像网格：

kzmig＝(k-1)·dzmig

其中，zmig2是成像终止深度，dzmig是成像深度步长，kzmig是第k个成像深度；

S405.将频率波数域炮检移动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是波数，dzmig是成像深度步长，i是虚数单位，fft是傅里叶变换函数，ifft是傅里叶逆变换函数，k_m是最大倾角dip对应的最大波数，mute是波数切除因子，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⁸