Let $t > 0$ and $\mu$ be a generic Radon measure on $\mathbb{R}^n$ and write
$$I := \int_{\mathbb{R}^n}\int_{\mathbb{R}^n}|x-y|^{-t}d\mu(x)d\mu(y)$$
Mattilas's Geometry of Sets and Measures in Euclidean Spaceas claims in chapte r8 that if the it holds true that
1.) $\mu(\mathbb{R}^n)<\infty$
2.) for some $s > t$ we have $\mu(B(x,r))\leq cr^s$ for some $c>0$ independent of $x,r$
then necessarily $I < \infty$ without any further comments on the matter. Mattila provides the following integral equality prior the claiming that $I < \infty$: For $x\in\mathbb{R}^n$, we have
$$\int_{\mathbb{R}^n}|x-y|^{-t}d\mu(y) = t\int_0^\infty r^{-t-1}\mu(B(x,r))dr$$
Well, I naive person, like myself, just wants to insert that equality brutely just to $I$ and estimate:
$$I = \int_{\mathbb{R}^n}t\int_0^\infty r^{-t-1}\mu(B(x,r))drd\mu(x)\leq \int_{\mathbb{R}^n}t\int_0^\infty r^{-t-1}cr^sdrd\mu(x)$$
whence
$$I\leq ct\int_{\mathbb{R}^n}\int_0^\infty r^{s-t-1}drd\mu(x) \leq ct\mu(\mathbb{R}^n)\left[\frac{1}{s-t}r^{s-t}\right]_0^\infty$$
where
$$\left[\frac{1}{s-t}r^{s-t}\right]_0^\infty = \infty$$
So clearly this is not the way. But then, what is? What should the sequence of equalities/inequalities be?
Combining 1) and 2), we have a better bound on $\mu(B(x,r))$:
$$ \mu(B(x, r)) \leq \min\{ cr^s, \mu(\mathbb{R}^n) \} $$
Plugging this to $I$,
\begin{align*} I &= \int_{\mathbb{R}^n} t \int_{0}^{\infty} r^{-t-1} \mu(B(x,r)) \, \mathrm{d}r\mathrm{d}\mu(x) \\ &\leq \int_{\mathbb{R}^n} t \int_{0}^{\infty} r^{-t-1} \min\{ cr^s, \mu(\mathbb{R}^n) \} \, \mathrm{d}r\mathrm{d}\mu(x) \\ &= \mu(\mathbb{R}^n) t \int_{0}^{\infty} \min\{ cr^{s-t-1}, \mu(\mathbb{R}^n) r^{-t-1} \} \, \mathrm{d}r. \end{align*}
Now denoting the integrand of the last integral by $f(r) := \min\{ cr^{s-t-1}, \mu(\mathbb{R}^n) r^{-t-1} \}$, we know:
$ f(r) = cr^{s-t-1} $ whenever $r$ is sufficiently close to $0$, and
$ f(r) = \mu(\mathbb{R}^n) r^{-t-1} $ whenever $r$ is sufficiently large.
This is enough to show that $I < \infty$.