Three different mixtures of 1-butyl-3-methyl-imidazolium tetrafluoroborate with water have been studied by means of molecular dynamics simulations. Based on the classical Lopes-Padua force field trajectories of approximately 60 ns were computed. This is the third part of a series concerning the collective network of 1-butyl-3-methyl-imidazolium tetrafluoroborate/water mixtures. The first part [C. Schröder et al., J. Chem. Phys. 127, 234503 (2007)] dealt with the orientational structure and static dielectric constants. The second part [C. Schröder et al., J. Chem. Phys. 129, 184501 (2008)] was focused on the decomposition of the dielectric spectrum of these mixtures. In this work the focus lies on the characterization of the neighborhood of ionic liquids by means of the Voronoi decomposition. The Voronoi algorithm is a rational tool to uniquely decompose the space around a reference molecule without using any empirical parameters. Thus, neighborhood relations, direct and indirect ones, can be extracted and were used in combination with g-coefficients. These coefficients represent the generalization of the traditional radial distribution function in order to include the mutual positioning and orientation of anisotropic molecules. Furthermore, the Voronoi method provides, as a by-product, the mutual coordination numbers of molecular species.