Yes, there is a more intuitive geometric explanation, though it is kind of difficult to see if you just get to see this formula. In differential (Riemannian) geometry one looks at curved (in contrast to flat, like Euclidean space) surfaces or higher dimensional manifolds. The metric you are looking at (more precisely: the pointwise norm associated with the pointwise scalar product it defines) is kind of an infinitesimal measure for distances in these surfaces.
It turns out that in this context you can set up differential calculus, the basic operation of which, when applied to vector fields, is the so called covariant differentiation. If you are looking at surfaces embedded in Euclidean 3 space this is basically a differentiation in that surrounding space followed by an orthogonal projection onto the tangent plane to the surface, but one can define this in an abstract manner, too (without an ambient manifold). If you then take a function and it's gradient (a concept which is also to be defined and depends on the metric) and take the covariant derivative of this object, the trace of this object (wrt the metric) is the Laplacian of the function (as is in Euclidean space, the Laplacian is the trace of the Hessian). In local coordinates this happens to look like the object you wrote down (when $f=1$). While in this form it looks a bit arbitrary it turns out to have some interesting properties. In particular it is invariant under coordinate changes, that is, it is well defined as a differential operator on the surface.
If you want to learn more about this you should fetch some basic textbook on Riemannian Geometry, like do Carmo.