基于有限元法和有限差分法的探地雷达正演模拟需要对计算区域进行网格剖分,剖分形式与网格大小对计算精度有直接影响且计算开销大。本文研究了一种新的探地雷达波动方程正演模拟算法-无单元法。无单元法只需节点,无需单元,摆脱了单元的束缚；采用节点信息及其局部支撑域上的权函数实现局部精确逼近；利用滑动最小二乘法(Moving Least Squares, MLS)构造形函数,形函数及导数具有很好的连续性；使得无单元法具有计算精度高、自适应能力强、灵活度大的特点。论文研究得到了国家自然科学基金项目“复杂GPR模型无网格法正演及多参数混合智能优化反演(41074085)”的支持。本文回顾了探地雷达正演模拟的发展历史、研究现状、无单元法在工程电磁场正演模拟中的应用；详细推导了探地雷达波波动方程及探地雷达波在地下介质中的传播特征；详细介绍了无单元法的基本原理,重点阐述了无单元法中形函数的构造、权函数的选取、影响域、高斯积分网格的划分及其实现等方法和措施,通过与有限元方法作对比,分析了无单元法用于正演模拟的优势。在此基础上,从探地雷达满足的波动方程出发,将无单元法用于探地雷达正演模拟当中,详细推导了探地雷达无单元法波动方程,对探地雷达无单元法正演模拟中的基本边界和吸收边界的加载问题进行了详细研究,重点研究了透射吸收边界、Sarma吸收边界的改进方法。开展了一系列的典型地电模型的探地雷达正演模拟,并与有限元法的计算结果进行对比,分析了无单元法正演模拟中的优越性,为探地雷达的反演奠定了坚实的基础；同时加深了电磁波在地下介质中的传播规律的理解和认识,有助于更好地指导探地雷达资料的地质解释。