Riverbank filtration (RBF) is increasingly being used as a relatively cheap and sustainable means to improve the quality of surface water. Due to the obvious differences in physical, chemical, and biological characteristics between river water and groundwater, there are strong and complex physical, chemical, and biogeochemical effects in the process of bank filtration. In this paper, multivariate statistical analysis was used to identify the spatial variation of hydrogeochemical groundwater in the process of bank filtration. Firstly, the evolution process of groundwater hydrochemistry during the filtration process was identified through factor analysis. According to the results, the evolution of groundwater hydrochemistry in this area is attributable to four main types of reactions: (1) Leaching; (2) Regional groundwater influence; (3) Aerobic respiration and denitrification; and (4) Mn (IV)/Fe (III)/SO42− reduction. According to the similarity of the geochemistry, the flow path could be divided into four different hydrochemical zones through cluster analysis, revealing the evolution law of groundwater hydrochemistry and its main influencing factors during riverbank infiltration. Large hydraulic gradient in The Zone Strongly Influenced by River Water (The first group) resulted in a weak effect of leaching on groundwater chemistry. Reoxygenation and microorganism respiration occurred in The Zone Moderately Influenced by River Water (The second group), The Zone Weakly Influenced by River Water (The third group), and The Zone Strongly Influenced by Regional Groundwater (The fourth group), resulting in fluctuations in Eh and pH values of groundwater. As a result, sulfate reduction and Mn (IV) and Fe (III) reduction alternated along the flow path. The Zone Strongly Influenced by Regional Groundwater (The fourth group) groundwater chemistry was mainly affected by regional groundwater.