【摘要】 为衰减地质雷达数值模拟时因截断区域厚度或者参数影响产生的倏逝波，提出了一种利用复频移完美匹配层(complex frequency shifted perfectly matched layer,CFS-PML)作为吸收边界，并结合频域高阶有限元算法（higher order finite element method,HO-FEM）求解电磁总场分量分布情况的新方法.该方法改善了截断域介质本构张量矩阵的特性，利用矩阵频移参量α吸收倏逝波，提高掠射角处吸收性能.数值实验结果表明，与传统PML相比，结合CFS-PML的HO-FEM方法精度得到明显提高，反射误差额外降低10～30 dB，节省了17%～30%的计算时间.此外，该方法能应用于复杂地质结构模型的电磁总场计算，为地质雷达数值模拟提供了一种新的方法.更多还原

【Abstract】 In the numerical simulation of ground penetrating radar(GPR), to attenuate the evanescent wave generated by the thickness of the truncated area or the influence of parameters, a new method is proposed to use the complex frequency shifted perfectly matched layer(CFS-PML) as the absorption boundary and combined with the higher-order finite element method(HO-FEM) in the frequency domain to solve the distribution of the total electromagnetic field components. The method improves the properties of the dielectric intrinsic tensor matrix in the truncated domain and uses the matrix frequency shift parameter α to absorb the evanescent wave and improve the absorption performance at the grazing angle. Numerical experiments are conducted on the method. Compared with the traditional PML, the accuracy of the HO-FEM method combined with CFS-PML is significantly improved, the reflection error is reduced by 10-30 dB extra, and the calculation time is saved by 17%-30%. In addition, the method can be applied to the calculation of the total electromagnetic field of a complex geological structure model, which provides a new method for the numerical simulation of geological lidars.更多还原