We introduce for SU(2) quantum spin systems the valence bond entanglement entropy as a counting of valence bond spin singlets shared by two subsystems. For a large class of antiferromagnetic systems, it can be calculated in all dimensions with quantum Monte Carlo simulations in the valence bond basis. We show numerically that this quantity displays all features of the von Neumann entanglement entropy for several one-dimensional systems. For two-dimensional Heisenberg models, we find a strict area law for a valence bond solid state and multiplicative logarithmic corrections for the Néel phase.