The authors consider a $d$-dimensional random field $u = /{u(t, x)/}$ that solves a non-linear system of stochastic wave equations in spatial dimensions $k /in /{1, 2, 3/}$, driven by a spatially homogeneous Gaussian noise that is white in time. They mainly consider the case where the spatial covariance is given by a Riesz kernel with exponent $/beta$. Using Malliavin calculus, they establish upper and lower bounds on the probabilities that the random field visits a deterministic subset of $/mathbb{R}^d$, in terms, respectively, of Hausdorff measure and Newtonian capacity of this set. The dimension that appears in the Hausdorff measure is close to optimal, and shows that when $d(2-/beta) > 2(k+1)$, points are polar for $u$. Conversely, in low dimensions $d$, points are not polar. There is, however, an interval in which the question of polarity of points remains open.