[发明专利]基于岩相阀值的混阶各向异性扩散地震图像去噪方法

权利要求书

1.基于岩相阀值的混阶各向异性扩散地震图像去噪方法，其特征在于，包括如下步骤：

步骤1，基于岩相信息，对待处理的地震图像进行分割；将地震图像划分为n个子图像S_i(i＝1,2,3,...,n)，每个子图像对应不同岩性的地层，地震子图像的个数等于层位数n；

步骤2，设置最大迭代次数t_max、初始扩散系数和初始岩相阀值；对子图像分别设置初始扩散系数Q＝[Q_i](i＝1,2,...,n)和初始岩相阀值w＝[w_i](i＝1,2,...,n)，即每个分割地震子图像S_i对应一个初始扩散系数和初始岩相阀值，并设置初始岩相阀值均为0；

步骤3，使用混阶各向异性扩散方程对地震子图像S_i分别进行处理；混阶各向异性扩散方程表达式为：

其中

其中，S为地震图像，t为迭代次数，w为岩相阀值，Difu₂(·)和Difu₄(·)分别表示二阶和四阶各向异性扩散项，Gds(·)为梯度终止函数，Q为扩散系数，div表示散度，▽表示梯度，▽²表示拉普拉斯算子；

步骤4，使用梯度数值统计法计算第t次迭代的地震子图像S_i的扩散系数和岩相阀值；扩散系数的表达式为：

其中，S_i^t表示第t次迭代的第i个地震图像，Q_2i和Q_4i分别表示对应Difu₂(·)和Difu₄(·)的扩散系数；Histc₁(·)和Histc₂(·)为数值统计函数，分别表示对图像的一阶和二阶梯度绝对值进行直方图统计，取90％分位数对应的梯度绝对值；

岩相阀值的表达式为：

其中，abs(·)表示取绝对值；

在首次迭代，由于岩相阀值为0，四阶各向异性扩散项的系数为0，混阶各向异性扩散方程即为二阶各向异性扩散方程；在迭代后期，由于岩相阀值为0.5，二阶各向异性扩散项的系数为0，混阶各向异性扩散方程即为四阶各向异性扩散方程；在迭代中期，二阶各向异性扩散方程和四阶各向异性扩散方程同时发挥作用；

步骤5，进入第t+1次迭代，重复步骤3，直至达到最大迭代次数，输出处理后的地震图像。