$\iint f(\mathbf x) |dS|$
Note that $\mathbf x \in \mathbb R^3$ or $\mathbf x = (x,y,z)$
$f(x) = \sqrt {(x-x_0)^2 + (y-y_0)^2 + (z-z_0)^2}$
Lets do this in spherical
$x = \cos\theta\sin\phi\\
y = \sin\theta\sin\phi\\
z = \cos\phi$
$|dS|= \|(\frac {\partial x}{\partial\phi},\frac {\partial y}{\partial\phi},\frac {\partial z}{\partial\phi})\times(\frac {\partial x}{\partial\theta},\frac {\partial y}{\partial\theta},\frac {\partial z}{\partial\theta})\| = \sin\phi$
$f(x) = \sqrt {x^2 + y^2 + z^2 - 2xx_0 - 2yy_0 -2zz_0 + x_0^2 +y_0^2 + z_0^2}$
$f(x) = \sqrt {1 - 2\cos\theta sin\phi x_0 - 2\sin\theta\sin\phi y_0 -2\cos\phi z_0 + x_0^2 +y_0^2 + z_0^2}$
$\iint \sqrt {1 - 2\cos\theta \sin\phi x_0 - 2\sin\theta\sin\phi y_0 -2\cos\phi z_0 + x_0^2 +y_0^2 + z_0^2} \sin\phi \;d\phi\;d\theta$
The geometry of this is symmetric. We don't need to use the point $(x_0,y_0,z_0)$ the result will be the same regardless of the orientation of the point. All we need is the distance from the origin. $x_0, y_0 = 0, z_0=d$
$\iint \sqrt {1 - 2\cos\phi d + d^2} \sin\phi \;d\phi\;d\theta$
$u^2 = 1 - 2d\cos\phi + d^2\\
2u\; du = 2 d\sin\phi\;d\phi$
$\int_0^{2\pi}\int_{d-1}^{d+1} \frac {u^2}{d} \;du\;d\theta$
Assuming $d>1$. I suppose I could say:
$\int_0^{2\pi}\int_{|d-1|}^{|d+1|} \frac {u^2}{d} \;du\;d\theta$
$\frac 23 \pi (6d + \frac {2}{d})$ if $d>1$
Wich is the surface area times the distance to the center plus a little bit, and that little bit gets increasingly trivial as $d$ gets large.
$\frac 23 \pi (6 + 2d^2)$ if $d<1$
Again equals the surface area plus a little bit and tends to the surface area exactly if $d=0$