首页|期刊导航|物探化探计算技术|三角形网格有限元法逆时偏移的面法向量导向滤波方法

三角形网格有限元法逆时偏移的面法向量导向滤波方法OA

Face normals guided filtering for FEM-based reverse time migration on triangular meshes

中文摘要英文摘要

有限元法基于三角形/四面体网格离散,对起伏地表与复杂构造边界具有天然拟合优势.逆时偏移虽能对复杂地质体实现高精度成像,但其结果常受低频噪声干扰,导致剖面分辨率下降.融合有限元法的几何适应性与逆时偏移的高精度成像能力可提升复杂区块成像可靠性,但成像滤波环节面临网格邻域不规则引起的空间滤波偏差,有限差分法框架下的滤波算子难以直接移用.为此,笔者提出一种面向三角形网格的面法向量导向滤波方法.该方法先对单元面法向量实施双边滤波,在质心位置保留法向量特征;再通过邻域三角形面积加权,将滤波后的法向量信息传递至节点以更新成像值.在波场延拓阶段,空间离散采用优化矩阵带宽方法以降低计算消耗,时间推进则基于 Newmark格式以保证数值稳定.经塔河工区实际地震数据测试,本文方法所得普拉特品质因数较均值滤波方法提升 6.51%.模型试验进一步验证,该方法在有效压制逆时偏移低频噪声的同时,能清晰保持地层边界特征,为强起伏地表及复杂构造区地震勘探提供了一种可行的成像增强技术路径.

The finite element method(FEM)employs triangular/tetrahedral mesh discretization,offering inherent advantages in fitting rugged topography and complex structural boundaries.While reverse time migration(RTM)achieves high-precision imaging of intricate geological bodies,its results are frequently contaminated by low-frequency noise,thereby degrading profile resolution.Integrating the geometric adaptability of the FEM with the imaging fidelity of RTM enhances reliability in complex regions;however,the imaging filtering stage suffers from spatial filtering deviations induced by irregular mesh neighborhoods,rendering it difficult to directly transfer filtering operators developed under the finite difference method framework.To address this issue,this study proposes a face normals-guided filtering approach tailored for triangular meshes.The method first applies bilateral filtering to the face normals of individual elements,preserving normal features at the element centroids.Subsequently,the filtered normals information is propagated to mesh nodes via area-weighted averaging of adjacent triangles to update the imaging amplitudes.During wavefield extrapolation,spatial discretization employs an optimized matrix bandwidth technique to reduce computational overhead,while temporal integration uses the Newmark scheme to ensure numerical stability.Tests on field seismic data from the Tahe area demonstrate that the proposed method achieves a Pratt's figure of merit(FOM)improvement of 6.51%compared to conventional mean filtering.Model experiments further confirm that this approach effectively suppresses RTM low-frequency noise while distinctly preserving stratigraphic boundary characteristics,thereby providing a viable technical pathway for enhancing seismic imaging in rugged topography and structurally complex zones.

李垒;文晓涛;陈明春;雷扬;吴闻静;张陈强

成都理工大学 地球物理学院,成都 610059成都理工大学 地球物理学院,成都 610059中石化石油工程地球物理有限公司南方分公司,成都 610041中石化石油工程地球物理有限公司南方分公司,成都 610041中石化石油工程地球物理有限公司南方分公司,成都 610041中石化石油工程地球物理有限公司南方分公司,成都 610041

天文与地球科学

有限元法逆时偏移面法向量三角形网格

finite element methodreverse time migrationface normalstriangular meshes

《物探化探计算技术》 2026 (4)

595-602,8

国家科技重大专项(2025ZD1004600)

10.12474/wthtjs.20260113-0001

评论