本发明涉及一种压制地震资料深层多次波的方法,属于地震数据处理,尤其涉及地震资料处理时保真提高信噪比处理阶段中压制多次波的方法。
背景技术:
1、复杂地质构造的深层勘探开发展示了深层油气藏勘探的良好前景,然而在深层地质勘探过程中,由于上覆地层存在强反射界面,造成深层多次波发育,严重影响地震资料的信噪比和保真性,制约着油气勘探的进程。
2、现有技术中,为解决地震资料处理中深层多次波的问题,通常利用多次波与一次有效波在频率、速度等方面的差异进行压制,例如,公开号为cn112925023a的中国专利文件,公开了一种全波场反演地震数据多次波压制方法,该方法基于τ-p变换,利用双平面波域对地震数据的压缩特性,通过多次迭代反演进行多次波压制和一次波估计。然而由于层间多次波在频率和速度方面与一次有效波重叠,难以识别,导致对深层多次波的压制效果有限,难以满足当前精细勘探的需求。
技术实现思路
1、本发明的目的在于提供一种压制地震资料深层多次波的方法,用于解决对地震资料处理中对深层多次波压制时难以对层间多次波进行识别和压制的问题。
2、为了实现上述目的,本发明提供了一种压制地震资料深层多次波的方法,包括如下步骤:s1、对待压制地震资料进行τ-p变换,得到τ-p域地震资料;
3、s2、通过预先得到的若干多次波算子,分别对τ-p域地震资料中的层间多次波进行压制,得到对应各多次波算子的一次有效波地震资料;
4、s3、计算各一次有效波地震资料的能量并进行比较,确定出与能量最小的一次有效波地震资料对应的多次波算子作为最优多次波算子;
5、s4、将通过最优多次波算子对τ-p域地震资料中的层间多次波进行压制得到的一次有效波地震资料进行反τ-p变换,得到多次波压制后地震资料。
6、进一步地,在上述方法中,步骤s2中所述的多次波算子根据地层反射系数进行确定,将地层反射系数与待压制地震资料做不同预测步长的褶积运算得到不同多次波算子。
7、进一步地,在上述方法中,通过如下公式得到一次有效波地震资料:
8、p-a(ω)p2+a2(ω)p3+…=p0
9、式中,p为τ-p域地震资料的总波场响应,a(ω)为多次波算子,p0为一次有效波地震资料的一次有效波响应。
10、进一步地,在上述方法中,步骤s3中,通过对如下公式进行求解得到一次有效波地震资料能量:
11、
12、式中,为第i炮的能量,i为炮数,表示第i炮总波场响应,t为地震波在介质中的传播时间,xr表示检波点位置,表示第i个炮点的位置,表示滤波算子,*表示褶积,n为迭代次数,为多次波响应。
13、进一步地,在上述方法中,还通过沿着p轴改变预测步长的方法对τ-p域地震资料进行多次波压制。
14、进一步地,在上述方法中,通过求解如下公式得到所述预测步长δt(p):
15、
16、式中,为n阶多次波旅行时,为n-1阶多次波旅行时,z为介质深度,δt(0)为地面接收点到地下反射点的垂直旅行时,p为地震资料的水平波慢度,v为多次波在介质中的传播速度。
1.一种压制地震资料深层多次波的方法,其特征在于,包括如下步骤:
2.根据权利要求1所述的压制地震资料深层多次波的方法,其特征在于,步骤s2中所述的多次波算子根据地层反射系数进行确定,将地层反射系数与待压制地震资料做不同预测步长的褶积运算得到不同多次波算子。
3.根据权利要求2所述的压制地震资料深层多次波的方法,其特征在于,通过如下公式得到一次有效波地震资料:
4.根据权利要求3所述的压制地震资料深层多次波的方法,其特征在于,步骤s3中,通过对如下公式进行求解得到一次有效波地震资料能量:
5.根据权利要求1所述的压制地震资料深层多次波的方法,其特征在于,还通过沿着p轴改变预测步长的方法对τ-p域地震资料进行多次波压制。
6.根据权利要求2或5所述的压制地震资料深层多次波的方法,其特征在于,通过求解如下公式得到所述预测步长δt(p):