This paper develops a robust and efficient PDE-constrained optimal control model for multicomponent pollutions in porous media, which takes into account nonlinear multi-component contamination flows of groundwater. The objective of the pollution optimal control is to identify the optimal injection rates from the top part boundary of domain, which can minimize the least squares error between the concentrations that are simulated and the allowable observed concentrations at observation sites, combining with the affection of costs associated with reducing emissions at injection locations. To discrete the constrained governing system of nonlinear multi-component flows, the splitting improved upwind finite difference scheme is developed for multicomponent PDEs system involving nonlinear chemical reactions of multicomponent pollutants and the pollutant injected rates on the upper part boundary. We employ the differential evolution (DE) optimization algorithm to solve the optimization. We numerically demonstrate the effectiveness of our model by analyzing the flow simulation on a simple geometric aquifer and identifying the optimal injection rates by minimizing the concentration derivation and the abatement costs. We also investigate the simulation of the contamination flow in a more realistic-shaped aquifer, which further validates our model's robustness and efficacy. The developed PDE-constrained control model and algorithm can be applied to applications of groundwater pollution control.