In this paper, the hydraulic characteristics of the fracture-karst aquifer and the distribution patterns of petrochemical contaminants are studied. Then, a numerical model using the mixed Eulerian-Lagrangian approach is constructed to predict the distribution and transport of petrochemical compounds in groundwater. The results of numerical modelling and sensitivity analysis show that it may be a workable way for aquifer remediation to combine contamination sources control and capture zone establishment.