Hi! Your hints really helped me: I solved it, now it's clear!
If you wanna check up my calculation, here they are:
$$
r_1>r_2
$$
$$
\begin{aligned}
J(r_1,r_2) &= \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 (\frac{1}{r_{12}}-\frac{1}{r_2}) \\
&= \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \frac{1}{r_{12}} - \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \frac{1}{r_2} \\
&= \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \frac{1}{r_{12}} - \frac{1}{r_2} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \\
&= \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \frac{1}{r_{12}} - \frac{1}{r_2} \int d\Omega_1 \hspace{5pt} Y^2_{00}(\Omega_1) \int d\Omega_2 \hspace{5pt} |Y_{lm}(\Omega_2)|^2 \\
&= \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \frac{1}{r_{12}} - \frac{1}{r_2} \\
&= I(r_1,r_2) - \frac{1}{r_2} \\
\end{aligned}
$$
$$
\begin{aligned}
I(r_1,r_2) &= \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \frac{1}{r_{12}} \\
&= \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \frac{1}{r_1} \sum_{l=0}^{\infty} (\frac{r_2}{r_1})^l \frac{4\pi}{2l+1} \sum_{m=-l}^l Y_{lm}(\Omega_1) Y^*_{lm}(\Omega_2) \\
&= \frac{4\pi}{2\sqrt{\pi}} \frac{1}{r_1} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y_{00}(\Omega_1) |Y_{lm}(\Omega_2)|^2 \sum_{l=0}^{\infty} (\frac{r_2}{r_1})^l \frac{1}{2l+1} \sum_{m=-l}^l Y_{lm}(\Omega_1) Y^*_{lm}(\Omega_2) \\
&= \frac{2\sqrt{\pi}}{r_1} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} |Y_{lm}(\Omega_2)|^2 \sum_{l=0}^{\infty} (\frac{r_2}{r_1})^l \frac{1}{2l+1} \sum_{m=-l}^l Y_{00}(\Omega_1) Y_{lm}(\Omega_1) Y^*_{lm}(\Omega_2) \\
&= \frac{2\sqrt{\pi}}{r_1} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} |Y_{lm}(\Omega_2)|^2 \sum_{l=0}^{\infty} (\frac{r_2}{r_1})^l \frac{1}{2l+1} \sum_{m=-l}^l Y^*_{00}(\Omega_1) Y_{lm}(\Omega_1) Y^*_{lm}(\Omega_2) \\
&= \frac{2\sqrt{\pi}}{r_1} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} |Y_{lm}(\Omega_2)|^2 \sum_{l=0}^{\infty} (\frac{r_2}{r_1})^l \frac{1}{2l+1} \sum_{m=-l}^l Y^2_{00}(\Omega_1) \delta_{0}^{l} \delta_{0}^{m} Y^*_{lm}(\Omega_2) \\
&= \frac{2\sqrt{\pi}}{r_1} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} |Y_{00}(\Omega_2)|^2 (\frac{r_2}{r_1})^0 \frac{1}{2 \cdot 0+1} Y^2_{00}(\Omega_1) Y^*_{00}(\Omega_2) \\
&= \frac{2\sqrt{\pi}}{r_1} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_2) Y^2_{00}(\Omega_1) Y_{00}(\Omega_2) \\
&= \frac{2\sqrt{\pi}}{r_1} \frac{1}{2\sqrt{\pi}} \int\int\ d\Omega_1 d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_2) Y^2_{00}(\Omega_1) \\
&= \frac{1}{r_1} \int d\Omega_1 \hspace{5pt} Y^2_{00}(\Omega_1) \int d\Omega_2 \hspace{5pt} Y^2_{00}(\Omega_2) \\
&= \frac{1}{r_1} \\
\end{aligned}
$$
$$
\begin{aligned}
J(r_1,r_2) &= I(r_1,r_2) - \frac{1}{r_2} \\
J(r_1,r_2) &= \frac{1}{r_1} - \frac{1}{r_2} \hspace{0.5cm} \blacksquare
\end{aligned}
$$
Thanks.
P.s.
I'm quite new to the forum, can you explain me how to mark the post as "Solved"?
Bye!