[发明专利]基于反稳相滤波的绕射波成像方法

权利要求书

1.基于反稳相滤波的绕射波成像方法，其特征在于包括:

1）采用常规的地震资料；

2）对地震资料做Kirchhoff叠前时间偏移得到全波场成像剖面；

3）应用平面波解构滤波方式从步骤2得到的全波场成像结果中提取每个成像点处的反射界面倾角信息；

4）利用步骤3的倾角信息构建核函数；

5）利用步骤4得到的核函数构建如下反稳相滤波器；

6）将步骤5中的反稳相滤波器作为Kirchhoff偏移的权函数对原始炮记录进行二次偏移，实现最终的绕射波单独成像。

2.根据权利要求1所述的基于反稳相滤波的绕射波成像方法，其特征在于进一步包括：

1）采用常规的叠前炮记录地震资料；

2）对叠前炮记录做Kirchhoff叠前时间偏移得到全波场成像剖面，基于绕射叠加理论，Kirchhoff偏移的数学表达式如下：

（式1）

式中：表示全波场叠前数据；表示由双平方根方程计算的双程旅行时；是Kirchhoff偏移权函数；当=1时，得到包括绕射波和反射波的全波场成像结果；

3）应用平面波解构滤波器从步骤2得到的全波场成像结果中提取每个成像点处的反射界面倾角信息，其数学表达式如下：

（式2）

式中：表示局部倾角预测误差算子，表示全波场成像剖面，求解式2的最小二乘问题即能实现全波场成像剖面的局部斜率估计；

4）利用步骤3估计的倾角信息构建如下核函数：

（式3）

式中：表示射线旅行时关于成像点的偏导数；表示每个成像点处的界面法线，由步骤3求得的界面倾角信息得到；,对于镜面反射，式3的值接近1，对于镜面反射之外的绕射，式3的值小于1，因此利用该核函数实现反射与绕射波场的分离；

5）利用步骤4得到的核函数构建如下反稳相滤波器：

（式4）

式中：为（3）式所示核函数，和为两个阈值，与地震频带有关，用于控制菲涅尔带的大小，为小于但接近于1的值，代表反射切除下界；为绕射保留上界；

6）将步骤5中的反稳相滤波器作为Kirchhoff偏移（式1）的权函数对原始炮记录进行二次偏移，实现最终的绕射波单独成像。

3.根据权利要求2所述的基于反稳相滤波的绕射波成像方法，其特征在于，步骤4中核函数的取值范围是，步骤5中的两个阈值和的遍历范围均为。