重力异常二维正演中的无网格方法

来源期刊：煤田地质与勘探2018年第6期

论文作者：李俊杰 严家斌

文章页码：181 - 186

关键词：无网格法;点插值法;径向基点插值法;无单元Galerkin法;重力异常;

摘    要：无网格法是一类新型数值算法,具有精度高、高阶形函数构造与物性加载便利等特点,在计算力学领域应用广泛。将无网格方法（PIM、RPIM及EFGM）用于重力异常场二维正演计算:首先从重力异常二维变分问题出发,利用Galerkin法结合高斯积分公式推导了对应的无网格离散系统矩阵表达式;其次通过数值试验得出了RPIM-MQ、RPIM-exp及EFGM-exp形状参数的建议值,最后比较分析了最优形状参数下不同无网格法的计算效果。结果表明:无网格法适用于介质物性分布变化较大的重力异常二维正演,exp函数形状参数c?最优取值区间为[1.5,1.7],?建议值为0.6,MQ函数q取值区间为–4.1～1.9;EFGM较PIM及RPIM具有更高的计算精度。