[发明专利]保地质边界的地震图像去噪方法及计算机可读存储介质

权利要求书

1.一种保地质边界的地震图像去噪方法，包括：

针对原始地震图像中任一点，构建结构张量；

针对所述原始地震图像中任一点，基于所述结构张量，构建扩散张量场；

基于所述扩散张量场，构建保地质边界去噪的非线性偏微分方程；

将所述非线性偏微分方程转化为线性方程组；

对所述线性方程组进行求解，获得去噪后的地震图像；

所述保地质边界去噪的非线性偏微分方程为：

其中，X＝(x,z)，表示原始地震图像中的任意一点，f(X)表示原始地震图像，s(X)表示去噪后的地震图像，和分别表示梯度算子和散度算子，α表示去噪因子，D(X)表示扩散张量场，x和z分别表示该点在原始地震图像中沿水平方向和垂直方向的坐标。

2.根据权利要求1所述的保地质边界的地震图像去噪方法，其中，所述结构张量表示为：

其中，G表示结构张量，g_x与g_z分别表示地震图像沿水平方向和垂直方向的梯度，＜·＞表示二维高斯光滑滤波。

3.根据权利要求2所述的保地质边界的地震图像去噪方法，其中，所述扩散张量场为：

其中，D表示扩散张量场，λ₁表示结构张量G的最大特征值，λ₂表示结构张量G的最小特征值，v₁表示局部图像法向方向单位向量，其正交于地震图像的主结构方向，v₂表示局部图像切向方向单位向量，其平行于地震图像的主结构方向。

4.根据权利要求1所述的保地质边界的地震图像去噪方法，其中，所述去噪因子α的取值范围在1至20之间。

5.根据权利要求1所述的保地质边界的地震图像去噪方法，其中，所述线性方程组为：

其中，f表示原始地震图像，s表示去噪后的地震图像，I表示单位阵，和分别表示梯度算子和梯度算子的转置，D表示扩散张量场。

6.根据权利要求5所述的保地质边界的地震图像去噪方法，其中，利用共轭梯度算法对所述线性方程组进行求解。

7.一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：

针对原始地震图像中任一点，构建结构张量；

针对所述原始地震图像中任一点，基于所述结构张量，构建扩散张量场；

基于所述扩散张量场，构建保地质边界去噪的非线性偏微分方程；

将所述非线性偏微分方程转化为线性方程组；

对所述线性方程组进行求解，获得去噪后的地震图像；

所述保地质边界去噪的非线性偏微分方程为：

其中，X＝(x,z)，表示原始地震图像中的任意一点，f(X)表示原始地震图像，s(X)表示去噪后的地震图像，和分别表示梯度算子和散度算子，α表示去噪因子，D(X)表示扩散张量场，x和z分别表示该点在原始地震图像中沿水平方向和垂直方向的坐标。