摘要
二维波动方程谱元法通常采用四边形网格,网格质量和属性建模方式对波场模拟精度和效率有较大影响.常规属性建模需要先建立几何模型再对其进行贴体网格剖分.但存在已知稠密属性控制点却难以建立复杂几何模型的问题.针对此问题,提出对模型进行矩形网格剖分并计算GLL点,利用属性控制点对所有GLL点进行插值以提高属性建模的效率.将属性建模方式归纳为单元属性建模方式和节点属性建模方式.在节点属性建模方式中分析了双线性插值、快速Gauss径向基函数插值法的计算效率,将两种插值算法集成于SPECFEM2D的属性建模程序中.实例表明两种属性建模方式对谱元法波场模拟都是有效的.
The spectral element method of 2 D wave equation is usually based on the quadrilateral mesh,and the mesh quality and attribute modeling methods have a great influence on the accuracy and efficiency of wave field simulation. The conventional attribute modeling needs to establish the geometric model and its body fitted mesh. However,the problem is that it is difficult to build a complex geometric model with dense attribute control points. In view of this problem,the model is divided into rectangular meshes and GLL( Gauss-Labotto-Legendre) points ' coordinates are computed.Properties of all GLL points are interpolated by attribute control points to improve the efficiency of attribute modeling. The attribute modeling method is summarized as the element attribute modeling method and node attribute modeling method. The efficiency of bilinear interpolation and fast Gauss radial basis function interpolation method are analyzed in node property modeling,and then these two interpolation algorithms are integrated into the attribute modeling program of SPECFEM2 D. Examples show that both modeling methods are valid for the spectral element wave field simulation.
作者
魏亦文
王有学
张智
WEI Yi-wen;WANG You-xue;ZHANG Zhi(College of Earth Sciences, Guilin University of Technology, Guilin 541004, China)
出处
《地球物理学进展》
CSCD
北大核心
2017年第5期2091-2096,共6页
Progress in Geophysics
基金
国家自然科学基金(41574039
41574078)
广西自然科学基金(2015GXNSFBA139184)
广西高校科研项目(KY2015LX113)联合资助
关键词
波场模拟
谱元法
GLL点
Gauss径向基函数
四边形网格
wave field simulation
spectral element method
GLL points
Gauss radial basis function
quadrilateral mesh