A coupling ground source heat pump system (CGSHP) is established in areas where groundwater is shallow but the seepage velocity is weak, which sets up pumping and injection wells on both sides of borehole heat exchangers (BHEs). A convection-dispersion analytical model of excess temperature in aquifer that considers groundwater forced seepage and axial effects and thermal dispersion effects is proposed. A controllable forced seepage sandbox is built by equation analysis method and similarity criteria. Through indoor test and the proposed analytical model, the correctness and accuracy of the numerical simulation software FEFLOW7.1 is verified. The influence of different pumping-injection flow rate on the heat transfer characteristic of BHEs is studied by numerical simulation. The results show that the average heat efficiency coefficient of BHEs increases and the heat influence range of downstream BHEs expands with the increasing of pumping-injection flow rate. The relation curve between Pe and the increment of heat transfer rate per unit depth of BHEs (Δ`q) is distributed as Gaussian function. The pumping-injection flow rate that makes Darcy velocity reaches 0.6 × 10− 6~1.4 × 10− 6 m∙s− 1 in the aquifer is the best reference range for CGSHP system,so 400 ~ 600 m3∙d− 1 is taken as the best pumping-injection flow rate in this paper.