The integral
$$\int^{2\pi}_{0}\ln{\sqrt{R^2+r^2-2 R\, r\cos{(\theta-\theta_0)}}}\,\mathrm{d}\theta \tag1$$
does not depend on $\theta_0$, as you can find by changing the variable to $\theta-\theta_0$ and using the periodicity of $\cos$. What's more, it does not depend on $r$ either (as long as $r<R$). Indeed, the function under the integral is harmonic in the disk $|x|<R$, being the fundamental solution of $\Delta$ with the pole at $(R,\theta_0)$. Therefore, its average over circle of radius $r$ is equal to the value at $0$.
$$\int^{2\pi}_{0}\ln{\sqrt{R^2+r^2-2 R\, r\cos{(\theta-\theta_0)}}}\,\mathrm{d}\theta = 2\pi \ln R \tag2$$
Just for completeness: there are other ways to get the findamental solution for the bi-Laplacian.
- Use the fact that every biharmonic function is of the form $r^2 g + h $ where $g,h$ are harmonic.
- Or, write down $\Delta^2$ in polar coordinates and solve the ODE.