To localize the source of a contaminant in the subsurface, based on concentration observations at some wells, we propose to test different possible locations and minimize the misfit between observed and simulated concentrations. We use a global optimization technique that relies on an expected improvement criterion, which allows a good exploration of the parameter space, avoids the trapping of local minima and quickly localizes the source of the contaminant on the presented synthetic cases.