Intermittent operation can improve the coefficient of performance (COP) of a ground source heat pump (GSHP) system. In this paper, an analytical solution to analyze the geo-temperature restoration performance under intermittent operation of borehole heat exchanger (BHE) fields is established. For this purpose, the moving finite line source model is combined with the g-function and the superposition principle. The model takes into account the heat transfer along the borehole, thermal interference between BHEs, and the influence of groundwater flow. The accuracy of the model is validated through comparison with an experiment carried out under intermittent operation. The model makes it possible to analyze the geo-temperature restoration performance and its influencing factors, such as BHE spacing, heat flow rate, operation mode, and groundwater flow. The main conclusions of this work are as follows. The heat transfer along the borehole should be considered when analyzing the geo-temperature restoration performance. When the BHE spacing increases, the soil temperature change decreases and the heat recovery improves. Therefore, adequate borehole separation distance is essential in the case of a multiple BHE system with unbalanced load. The presence of groundwater flow is associated with interference between the BHEs, which should not be ignored. In the case of long-term operation, the groundwater flow is beneficial to the geo-temperature recovery process, even for downstream BHEs. Finally, a greater groundwater flux leads to a better geo-temperature recovery.