###### Abstract  
Groundwater pollution poses a significant threat to ecosystems and public safety. Traditional remediation methods have limitations, necessitating innovative approaches. This study integrates numerical modeling and bioremediation to address groundwater contamination in an industrial site. It explores the potential of chemotactic bacteria to enhance remediation efficiency. The research establishes groundwater pollutant transport models, analyzes flow fields, and assesses the distribution of various pollutants. The results demonstrate the effectiveness of chemotactic bacteria, particularly chemotactic bacteria that can rapidly adapt as the pollutant concentration decreases, the concentration of chemotactic bacteria in the low-permeability area has increased by 112%. This study provides insights into the practical application of bioremediation and the promising role of bacterial chemotaxis in treating contaminated groundwater. chemotactic bacteria; bioremediation; groundwater pollutant; numerical simulation 202420202430043004300430043004300430043004300430043004300430043004300430043004300430043004300430043004300430043004300430043004304300430043043004300430430  
viable and widely applied in groundwater remediation projects internationally [9]. However, overcoming aquifer heterogeneity, with variations in porosity and permeability across different geological layers [10], remains a critical challenge in improving the efficiency of in situ bioremediation. Chemotaxis is a capability evolved by many bacteria in the natural world, enabling them to sense concentration gradients of specific compounds (beneficial or harmful) and swim towards regions conducive to their survival and proliferation [11]. Several strains of naturally occurring harmless bacteria have been found not only to degrade petroleum pollutants but also to exhibit positive chemotactic responses to pollutants [12]. Studies at the population level have also been conducted. [PERSON] et al. [13] employed capillaries filled with micrometer-sized glass beads as a porous medium and observed the enrichment of the bacterium _Pseudomonas putida_ (_P. putida_) G7 within the capillaries when naphthalene was present in the aqueous phase, a phenomenon not observed in the control group without naphthalene, thus providing early proof of the role of chemotactic bacteria in homogeneous porous media. [PERSON] et al. [14] used chemotaxis to assess the chemotactic behavior of the complete degrading bacterium _Pseudomonas_ sp. ZWL73 and its conjugative transposon ZP8 toward 4-chloronitrobenzene (4 CNB) and 12 compounds potentially related to 4 CNB degradation, most of which could serve as substrates for ZWL73 growth. [PERSON] et al. [15] studied the chemotactic characteristics of the actinobacterium _Arthobacter_ sp. CN2, usingnitrophenol as a substrate, and found a positive correlation between the bacterium's chemotaxis andnitrophenol concentration. [PERSON] et al. [16] created sandbox simulations to mimic cross-sectional homogeneous soil, making the experimental system closer to the actual environment, and observed the chemotactic migration of _P. putida_ F1 toward the contaminated layer. Although these population-level studies have provided ample experimental and theoretical foundations for the application of chemotactic bacteria in porous media, their experimental systems remain overly simplistic, and they share a common limitation of being unable to replicate heterogeneous soil environments. Due to the high cost of groundwater pollution site remediation and the potential for improper remediation methods to cause secondary pollution at the site, many scholars have used numerical simulation methods in recent years to predict the transport behavior of pollutants at the site. [PERSON] et al. focused on the migration of manganese in groundwater using Visual MODFLOW and MT3 DMS software in Xiangtan, China [17]. [PERSON] et al.