Ground Penetrating Radar (GPR) is a shallow geophysical method for detecting and locating subsurface targets. The GPR image echo characteristics of complex underground structures can be obtained by carrying out GPR forward modeling research. The traditional finite-difference time-domain (FDTD) method has low efficiency and accuracy. The alternating direction implicit FDTD (ADI-FDTD) algorithm surmounts the stability limitations of the traditional FDTD method, making it possible to select a larger time step for higher computational efficiency. For circular underground structures, a pseudowave produced by the ladder approximation method can be corrected using the surface conformal technique. This paper proposes a high-efficiency and high-accuracy GPR forward modeling method that combines the ADI-FDTD algorithm and surface conformal technology. The performance of the conformal ADI-FDTD algorithm is verified by a simple two-layer model. Based on the proposed algorithm, the GPR image features of three complex underground structure models are obtained. Finally, a field experiment is used to support the accuracy and usefulness of the conformal ADI-FDTD algorithm. The numerical simulation results and experimental results show that the conformal ADI-FDTD algorithm reduces the pseudodiffraction wave caused by the ladder approximation method and can significantly improve the computing efficiency for complex underground structure models.