中文摘要：

应用迭代插值方法构造了插值小波尺度函数,并将该尺度函数的导数用于离散Maxwell方程组的空间微分,使用四阶Runge Kutta（four order Runge Kutta,RK4）算法计算时间导数,导出了插值小波尺度法的探地雷达（ground penetrating radar,GPR）正演公式,与常规的基于中心差分的时域有限差分算法（finite difference time domain,FDTD）相比,插值小波尺度算法提高了GPR波动方程的空间与时间离散精度.首先,采用具有解析解的层状模型,分别将FDTD算法及插值小波尺度法应用于层状模型正演,单道雷达数据与解析解拟合表明：相同的网格剖分方式,插值小波尺度法比FDTD具有更高的精度.然后,将辅助微分方程完全匹配层（auxiliary differential equation perfecting matched layer,ADE-PML）边界条件应用到插值小波尺度法GPR正演中,在均匀介质模型中对比了FDTD-CPML（坐标伸缩完全匹配层）,FDTD-RK4ADE-PML、插值小波尺度RK4ADE-PML的反射误差,结果表明：插值小波尺度RK4ADE-PML吸收效果优于另外两种条件下的吸收边界.最后,应用加载UPML（各向异性完全匹配层）的FDTD和RK4ADE-PML的插值小波尺度法开展了二维GPR模型的正演,展示了RK4ADE-PML对倏逝波的良好吸收效果.

英文摘要：

Ground penetrating radar（GPR） forward is one of the geophysical research directions. Through the forward of geological model, the database of radar model can be enriched and the characteristics of typical geological radar echo images can be understood, which in turn can guide the data interpretation of GPR measured profile, thereby improving the GPR data interpretation level. In this article, the interpolating wavelet scale function by using iterative interpolation method is presented, and the derivative of scale function is used in spatial differentiation of discrete Maxwell equations.The forward modeling formula of GPR based on the interpolation wavelet scale method is derived by using fourth-order Runge-Kutta method（RK4） for calculating the higher time derivative. Compared with the conventional finite difference time domain（FDTD） algorithm based on the central difference method, the interpolation wavelet scale algorithm improves the accuracy of GPR wave equation in both space and time discretization. Firstly, the FDTD algorithm and the interpolation wavelet scale method are applied to the forward modeling of a layered model with analytic solution.Single channel radar data and analytical solution fitting indicate that the interpolation wavelet scale method has higher accuracy than FDTD, with the same mesh generation used. Therefore, auxiliary differential equation perfectly matching layer（ADE-PML） boundary condition is used on an interpolation wavelet scale, and the comparisons between reflection errors obtained using CPML（FDTD）, RK4ADE-PML（FDTD）, and RK4ADE-PML（interpolating wavelet scales） in a homogeneous medium model show that the absorption effect of RK4ADE-PML（interpolating wavelet scales） is better than the other two absorbing boundaries. Finally, interpolation wavelet scale method, with both UPML, FDTD and RK4ADE-PML loaded, is used for two-dimensional GPR forward modeling, showing good absorption effect for evanescent wave. From all the experimental results, the followin