速度域的海底多次波预测与追踪衰减方法技术

本发明专利技术涉及一种速度域的海底多次波预测与追踪衰减方法，属于地震资料处理与分析领域。本发明专利技术具体的实施过程主要分为以下四步：1)高分辨率速度域记录的创建；2)基于初始速度域记录预测海底多次波速度谱；3)双曲型多次波同相轴的追踪；4)利用FK扇形滤波的多次波同相轴衰减。与常规多次波压制方法相比，本方法不依赖多次波的周期，不须对输入的CMP道集进行动校正处理；其在时空域中进行多次波压制，不受截断效应的影响；同时可自由设定所截取记录的时窗长度，当时窗较小时能够衰减叠加速度与一次波较为接近的多次波，从而避免了对一次波的损伤；相较于自由界面多次波衰减方法，其可有效压制远偏移距道中的多次波，而且具有明显的效率优势。的效率优势。的效率优势。

【技术实现步骤摘要】
速度域的海底多次波预测与追踪衰减方法


[0001]本专利技术涉及地震资料处理与分析领域，是一种应用于地震数据的速度域海底多次波预测与追踪衰减方法。

技术介绍

[0002]目前，应用于油气勘探的地震勘探方法主要为反射波法，其通常认为一次反射是有效的，而直达波、浅层折射波和多次波等均为干扰波。对于直达波及浅层折射波等干扰，一般可直接在记录上予以切除，而多次波的识别与压制问题则是长期以来人们一直重点研究的热点问题之一。对于海上地震勘探来说，由于海面和海底一般均为强波阻抗界面，相应的地震记录中通常包含大量振幅较强的海底多次波。若不予以消除，其将严重影响地震成像的真实性与可靠性，进而误导后续的地震地质解释。
[0003]当前的多次波剔除方法主要有预测反褶积技术、基于视速度差异的多次波压制方法与基于波动理论的多次波压制方法等。预测反褶积具有高效率的优点，但多次波的周期呈时变、空变特征，因此无法给定准确的预测步长，从而导致多次波残余或损伤有效信号，而且其难以有效衰减长周期多次波。基于视速度差异的多次波压制方法在实际地震数据处理中得到了广泛应用，但由于近偏移距道中的时差较小，常导致强振幅多次波同相轴的严重残余，而且引入的动校正、反动校正等额外处理步骤往往会使处理结果中产生波形畸变。在基于波动理论的多次波压制方法中，自由界面多次波衰减方法(SRME)已成为地震数据处理的首选方法之一,但其难以有效压制远偏移距道中的多次波，而且运算量较大，这影响了SRME在实际生产中的应用效果。

技术实现思路

[0004]本专利技术要解决的技术问题是提供一种速度域的海底多次波预测与追踪衰减方法。其首先将炮集记录选排到CMP域，利用同相加权技术创建高分辨率的叠加速度谱；然后根据海底模型确定当前速度谱位置的海底深度与海水速度，于速度域预测得到仅包含海底多次波能量的叠加速度谱；最后基于同相轴追踪技术确定时空域中的多次波同相轴，通过记录重排的手段将多次波同相轴校正为水平，并应用视速度滤波方法予以压制。该方法实现的关键有两点：即海底多次波速度谱的预测和多次波同相轴的追踪，其通过前者保证追踪的均为多次波同相轴，而通过后者获得所追踪同相轴精确的旅行时信息。
[0005]本专利技术是通过如下技术方案来实现的：
[0006]速度域的海底多次波预测与追踪衰减方法，所述方法的具体步骤如下：
[0007]1)高分辨率速度域记录的创建在海上地震勘探中，通过人工震源连续走航式的激发地震波，接收地震信号，获得包含多次波的系列炮记录；将炮集记录选排到CMP域，对获得的CMP记录d(x,t)作速度叠加变换，将得到横向坐标为叠加速度v、纵坐标为零偏移距时τ的初始速度域记录，其计算公式表示为
[0008][0009]式中：x、t、v和τ分别为偏移距、旅行时、叠加速度和零偏移距时，x1和x2分别表示最小偏移距和最大偏移距，d和u分别表示偏移距域记录和初始速度域记录；
[0010]为使初始速度域记录u(v,τ)具有更高的分辨率，计算各速度值的同相加权因子w(v,τ)，即当地震记录中存在该速度的同相轴时此加权因子的值较大，否则较小或者为零
[0011][0012]式中，l表示时窗长度；N
x
为偏移距域记录的道数；x、t、v和τ分别为偏移距、旅行时、叠加速度和零偏移距时，x1和x2分别表示最小偏移距和最大偏移距，d表示偏移距域记录，C为保证分母不为零的常数，取平均振幅的0.01～0.001；
[0013]利用式(2)计算的同相加权因子w(v,τ)对初始速度域记录u(v,τ)进行同相加权，即可获得高分辨率的速度域记录s(v,τ)
[0014]s(v,τ)＝u(v,τ)w(v,τ)
ꢀꢀ
(3)
[0015]式中，v和τ分别为叠加速度和零偏移距时；
[0016]2)基于公式(3)获得的初始速度域记录预测海底多次波速度谱输入包含海底深度和海水速度的海底模型，根据当前速度域记录的空间位置确定海底深度z和海水速度v
w
；若海底较为平坦，对于高分辨率速度域记录s(v,τ)中任意v、τ处的反射能量，其在水层中沿垂向再往返传播一次，将形成海底多次波的反射能量，因此求出相应海底多次波能量的速度参数v
m
和τ
m
，进而预测出仅包含海底多次波的速度域记录s
m
(v
m
,τ
m
)，即
[0017][0018]式中，m标志与多次波相关的数据或参量；v1和v2表示速度域记录的最小速度值和最大速度值；τ1和τ2表示速度域记录的最小零偏移距时和最大零偏移距时，τ1通常为0；
[0019]对矩阵s
m
取绝对值即可获得仅含有海底多次波能量的叠加速度谱E(v,τ)
[0020]E(v,τ)＝|s
m
(v,τ)|
ꢀꢀ
(5)
[0021]式中，s
m
为仅包含海底多次波的速度域记录，v和τ分别为叠加速度和零偏移距时；
[0022]3)基于公式(5)获得的仅含有海底多次波能量的叠加速度谱E(v,τ)进行双曲型同相轴的追踪对E(v,τ)进行适度平滑后，d(x,t)中叠加速度值v
i
、零偏移距时τ
i
的双曲线同相轴，将在速度谱E(v,τ)中形成以(v
i
,τ
i
)为中心极值的团状结构能量；应用等值线追踪方法求出各能量团的分布范围，并搜索出其极值点位置，则根据该点的坐标(v
i
,τ
i
)拟合出时空域中的相应同相轴，即其所经各道的旅行时为
[0023][0024]式中，i(i≥1)表示追踪同相轴的序号，n表示CMP记录d(x,t)中各道数据的道号，x
n
表示第n道记录的偏移距，v
i
和τ
i
分别代表所追踪同相轴的叠加速度值和零偏移距时；
[0025]4)利用FK扇形滤波进行多次波同相轴衰减基于同相轴追踪过程确定出多个多次波同相轴，将其准确参数τ
i
及v
i
代入公式(6)计算所经地震道的旅行时则通过FK扇形滤波法进行消除；针对每个同相轴的滤波过程为：
①
在各地震道中以为中心截取给定的一个短时窗长度的记录段，使各记录段沿起点位置对齐，从而将目标同相轴校正为水平；
②
以截取的多道记录段作为输入，通过FK扇形滤波滤出已被校正为水平的同相轴；
③
将滤波后的记录反重排，并放回各地震道的原时窗位置。
[0026]本专利技术与现有技术相比的有益效果：
[0027]与常规多次波压制方法相比，本方法不依赖多次波的周期，不须对输入的CMP道集进行动校正处理；其在时空域中进行多次波压制，不受截断效应的影响；同时可自由设定所截取记录的时窗长度，当时窗较小时能够衰减叠加速度与一次波较为接近的多次本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种速度域的海底多次波预测与追踪衰减方法，其特征在于所述方法的具体步骤如下：1)高分辨率速度域记录的创建在海上地震勘探中，通过人工震源连续走航式的激发地震波，接收地震信号，获得包含多次波的系列炮记录；将炮集记录选排到CMP域，对获得的CMP记录d(x,t)作速度叠加变换，将得到横向坐标为叠加速度v、纵坐标为零偏移距时τ的初始速度域记录，其计算公式表示为式中：x、t、v和τ分别为偏移距、旅行时、叠加速度和零偏移距时，x1和x2分别表示最小偏移距和最大偏移距，d和u分别表示偏移距域记录和初始速度域记录；为使初始速度域记录u(v,τ)具有更高的分辨率，计算各速度值的同相加权因子w(v,τ)，即当地震记录中存在该速度的同相轴时此加权因子的值较大，否则较小或者为零式中，l表示时窗长度；N
x
为偏移距域记录的道数；x、t、v和τ分别为偏移距、旅行时、叠加速度和零偏移距时，x1和x2分别表示最小偏移距和最大偏移距，d表示偏移距域记录，C为保证分母不为零的常数，取平均振幅的0.01～0.001；利用式(2)计算的同相加权因子w(v,τ)对初始速度域记录u(v,τ)进行同相加权，即获得高分辨率的速度域记录s(v,τ)s(v,τ)＝u(v,τ)w(v,τ)
ꢀꢀꢀꢀꢀꢀꢀꢀ
(3)式中，v和τ分别为叠加速度和零偏移距时；2)基于公式(3)获得的初始速度域记录预测海底多次波速度谱输入包含海底深度和海水速度的海底模型，根据当前速度域记录的空间位置确定海底深度z和海水速度v
w
；若海底较为平坦，对于高分辨率速度域记录s(v,τ)中任意v、τ处的反射能量，其在水层中沿垂向再往返传播一次，将形成海底多次波的反射能量，因此求出相应海底多次波能量的速度参数v
m
和τ
m
，进而预测出仅包含海底多次波的速度域记录s
m
(v
m
,τ
m
)，即式中，m标志与多次波相关的数据或参量；v1和v2表示速度域记录的最小速度...

【专利技术属性】
技术研发人员：谭军，姜秀萍，宋鹏，李金山，夏冬明，赵波，
申请(专利权)人：中国海洋大学，
类型：发明
国别省市：