As far as I could understand, it seems that you want to know whether timelike geodesics can reach the conformal boundary of AdS. If that's the case (please do confirm), the answer is no - no timelike geodesic can reach conformal infinity, it rather gets constantly refocused back into the bulk in a periodic fashion. You need timelike curves which have some acceleration in order to avoid this. Maximally extended null geodesics (i.e. light rays), on the other hand, always reach conformal infinity, both in the past and in the future. An illustration of these facts using Penrose diagrams can be found, for instance, in Section 5.2, pp. 131-134 of the book by S. W. Hawking and G. F. R. Ellis, "The Large Scale Structure of Space-Time" (Cambridge, 1973).
The detailed reasoning behind the above paragraph can be seen in a global, geometric way. In what follows, I'll largely follow the argument presented in the book by B. O'Neill, "Semi-Riemannian Geometry - With Applications to Relativity" (Academic Press, 1983), specially Proposition 4.28 and subsequent remarks, pp. 112-113. For the benefit of those with no access to O'Neill's book, I'll present the self-contained argument in full detail. I'll make use of the fact that $AdS_4$ is the universal covering of the embedded hyperboloid $H_m$ ($m>0$) in $\mathbb{R}^{2,3}=(\mathbb{R}^5,\eta)$
$$ H_m=\{x\in\mathbb{R}^5\ |\ \eta(x,x)\doteq -x_0^2+x_1^2+x_2^2+x_3^2-x_4^2=-m\}\ . $$
The covering map $\Phi:AdS_4\ni(t,r,\theta,\phi)\mapsto (x_0,x_1,x_2,x_3,x_4)\in H_m\subset\mathbb{R}^{2,3}$ through the global coordinates $(t\in\mathbb{R},r\geq 0,0\leq\theta\leq\pi,0\leq\phi<2\pi)$ is given by
$$ x_0=\sqrt{m(1+r^2)}\sin t\ ;$$ $$ x_1=\sqrt{m}r\sin\theta\cos\phi\ ;$$ $$ x_2=\sqrt{m}r\sin\theta\sin\phi\ ;$$ $$ x_3=\sqrt{m}r\cos\theta\ ;$$ $$ x_4=\sqrt{m(1+r^2)}\cos t\ .$$
The pullback of the ambient, flat pseudo-Riemannian metric $\eta$ defined above (with signature $(-+++-)$) by $\Phi$ after restriction to $H_m$ yields the $AdS_4$ metric in the form appearing in the question and in Pedro Figueroa's nice answer up to a constant, positive factor:
$$ds^2= m\left[-(m+r^2)dt^2+(m+r^2)^{-1}dr^2+r^2(d\theta^2+\sin^2\theta d\phi^2)\right]\ .$$
The conformal completion of $AdS_4$, on its turn, is obtained by means of the change of radial variable $u=\sqrt{m+r^2}-r$, so that $r=\frac{m-u^2}{2u}$, $dr=-\frac{1}{u}(\frac{m+u^2}{2u})du$ and $m+r^2=(\frac{m+u^2}{2u})^2$, yielding
$$ds^2=\frac{m}{u^2}\left[-\left(\frac{m+u^2}{2}\right)^2dt^2+du^2+\left(\frac{m-u^2}{2}\right)^2(d\theta^2+\sin^2\theta d\phi^2)\right]\ .$$
Conformal infinity is reached by taking $r\rightarrow+\infty$, which is the same as $u\searrow 0$. The rescaled metric $\Omega^2 ds^2$, $\Omega=m^{-\frac{1}{2}}u$ yields the three-dimensional Einstein static universe as the conformal boundary (i.e. $u=0$).
It's clear that $H_m$ is a level set of the function $f:\mathbb{R}^5\rightarrow\mathbb{R}$ given by $f(x)=\eta(x,x)$. Therefore, the vector field $X_x=\frac{1}{2}\mathrm{grad}_\eta f(x)=x$ (where $\mathrm{grad}_\eta$ is the gradient operator defined with respect to $\eta$) is everywhere normal to $H_m$ - that is, any tangent vector $X_x\in T_x H_m$ satisfies $\eta(X_x,T_x)=0$. Given two vector fields $T,S$ tangent to $H_m$, the intrinsic covariant derivative $\nabla_T S$ on $H_m$ is simply given by the tangential component of the ambient (flat) covariant derivative $(\partial_T S)^a=T^b\partial_b S^a$:
$$ \nabla_T S=\partial_T S-\frac{\eta(X,\partial_T S)}{\eta(X,X)}X=\partial_T S+\frac{\eta(X,\partial_T S)}{m}X\ .$$
The normal component of $\partial_T S$, on its turn, has a special form due to the nature of $H_m$ (notice that $\partial_a X^b=\partial_a x^b=\delta^b_a$):
$$ \eta(X,\partial_T S)=\underbrace{\partial_T(\eta(X,S))}_{=0\ ;}-\eta(S,\partial_T X)=-\eta(S,T)\ \Rightarrow\ \frac{\eta(X,\partial_T S)}{\eta(X,X)}X=\frac{\eta(S,T)}{m}X\ .$$
As such, we conclude that a curve $\gamma:I\ni\lambda\mapsto\gamma(\lambda)\in H_m$ ($I\subset\mathbb{R}$ is an interval with nonvoid interior) is a geodesic of $H_m$ if and only if $\frac{d^2\gamma(\lambda)}{d\lambda^2}(\lambda)\doteq\ddot{\gamma}(\lambda)$ is everywhere normal to $H_m$, that is,
$$\ddot{\gamma}(\lambda)=-\frac{1}{m}\eta(\ddot{\gamma}(\lambda),X_{\gamma(\lambda)})X_{\gamma(\lambda)}=\frac{1}{m}\eta(\dot{\gamma}(\lambda),\dot{\gamma}(\lambda))X_{\gamma(\lambda)}=\frac{1}{m}\eta(\dot{\gamma}(\lambda),\dot{\gamma}(\lambda))\gamma(\lambda)\ .$$
In particular, if $\eta(\dot{\gamma}(\lambda),\dot{\gamma}(\lambda))=0$, then $\gamma$ is also a (null) geodesic in the ambient space $\mathbb{R}^{2,3}$.
Given $x\in H_m$, the linear span of $X_x=x$ and any tangent vector $T_x\neq 0$ to $H_m$ at $x$ defines a 2-plane $P(T_x)$ through the origin of $\mathbb{R}^5$ and containing $x$. In other words,
$$ P(T_x)=\{\alpha X_x +\beta T_x\ |\ \alpha,\beta\in\mathbb{R}\}\ , $$
and therefore
$$ P(T_x)\cap H_m=\{y=\alpha X_x+\beta T_x\ |\ \eta(y,y)=-\alpha^2 m+\beta^2\eta(T_x,T_x)=-m\}\ .$$
This allows us already to classify $P(T_x)\cap H_m$ according to the causal character of $T_x$:
- $T_x$ timelike (i.e. $-k=\eta(T_x,T_x)<0$): we have that $m\alpha^2+k\beta^2=m$ with $k,m>0$, hence $P(T_m)\cap H_m$ is an ellipse;
- $T_x$ spacelike (i.e. $k=\eta(T_x,T_x)>0$): we have that $m\alpha^2-k\beta^2=m$ with $k,m>0$, hence $P(T_m)\cap H_m$ is a pair of hyperbolae, one with $\alpha>0$ and the other with $\alpha<0$. The point $x=X_x$ belongs to the first hyperbola;
- $T_x$ lightlike (i.e. $\eta(T_x,T_x)=0$): we have that $\alpha^2=1$ with $\beta$ arbitrary, hence $P(T_m)\cap H_m$ is a pair of straight lines, one given by $\alpha=1$ and the other by $\alpha=-1$. The point $x=X_x$ belongs to the first line. Notice that each of these lines is a null geodesic both in $H_m$ and in $\mathbb{R}^{2,3}$!
Moreover, $x=\gamma(0)$ and $T_x=\dot{\gamma}(0)$ define a general initial condition for a geodesic $\gamma$ starting at $x$. It remains to show that any curve that stays in $P(T_x)\cap H_m$ is a geodesic in $H_m$. This is clearly true for $T_x$ lightlike, since in this case we have already concluded that $\gamma(\lambda)=x+\lambda T_x$ for all $\lambda\in\mathbb{R}$. For the remaining cases (i.e. $\eta(T_x,T_x)\neq 0$), consider a $\mathscr{C}^2$ curve $\gamma$ in $P(T_x)\cap H_m$ beginning at $\gamma(0)=x$ with $\dot{\gamma}(0)=\dot{\beta}(0)T_x$ (we assume that $\dot{\gamma}(\lambda)\neq 0$ for all $\lambda$). Writing $\gamma(\lambda)=\alpha(\lambda)X_x+\beta(\lambda)T_x$, we conclude from the above classification of $P(T_x)\cap H_m$ that we can choose the parameter $\lambda$ so that
- $T_x$ timelike: $\alpha(\lambda)=\cos\lambda$, $\beta(\lambda)=\sqrt{-\frac{m}{\eta(T_x,T_x)}}\sin\lambda$, so that $\eta(\dot{\gamma}(\lambda),\dot{\gamma}(\lambda))=-m$ with $\dot{\beta}(0)=\sqrt{-\frac{m}{\eta(T_x,T_x)}}$;
- $T_x$ spacelike: $\alpha(\lambda)=\cosh\lambda$, $\beta(\lambda)=\sqrt{\frac{m}{\eta(T_x,T_x)}}\sinh\lambda$, so that $\eta(\dot{\gamma}(\lambda),\dot{\gamma}(\lambda))=+m$ with $\dot{\beta}(0)=\sqrt{\frac{m}{\eta(T_x,T_x)}}$.
In both cases, we conclude that
$$ \ddot{\gamma}(\lambda)=\frac{\eta(\dot{\gamma}(\lambda),\dot{\gamma}(\lambda))}{m}\gamma(\lambda)\ ,$$
i.e. $\gamma$ must satisfy the geodesic equation in $H_m$ with the chosen parametrization, as wished. Since any pair of initial conditions for a geodesic determines a 2-plane through the origin in the above fashion, we conclude that the resulting geodesic in $H_m$ will remain forever in that 2-plane. For later use, I remark that all geodesics of $H_m$ cross at least once the 2-plane $P_0=\{x\in\mathbb{R}^5\ |\ x_1=x_2=x_3=0\}$ - this can be easily seen from the classification of the sets $P(T_x)\cap H_m$. This allows us to prescribe initial conditions in $P_0$ for all geodesics in $H_m$.
Now we have complete knowledge of the geodesics in the fundamental domain $H_m$ of $AdS_4$. What happens when we go back to the universal covering? What happens is that the lifts of spacelike and lightlike geodesics stay confined to a single copy of the fundamental domain, whereas the lifts of timelike geodesics do not. To see this, we exploit the fact that translations in the time coordinate $t$ are isometries and the remark at the end of the previous paragraph to set $$\gamma(0)=X_x=x=(0,0,0,0,\sqrt{m})$$ in $H_m$ (i.e. $\gamma$ is made to start at $P_0$ with $t=0$), so that $$\dot{\gamma}(0)=T_x=(y_0,y_1,y_2,y_3,0)\ .$$ We also normalize $\eta(T_x,T_x)$ to $-m$, $+m$ or zero depending on whether $T_x$ is respectively timelike, spacelike or lightlike. Writing once more $\gamma(\lambda)=\alpha(\lambda)X_x+\beta(\lambda)T_x$, we use the classification of geodesics in $H_m$ by their causal character to write explicit formulae for $\gamma$:
- $T_x$ timelike $\Rightarrow$ $\gamma(\lambda)=(\cos\lambda)X_x+(\sin\lambda)T_x$;
- $T_x$ spacelike $\Rightarrow$ $\gamma(\lambda)=(\cosh\lambda)X_x+(\sinh\lambda)T_x$;
- $T_x$ lightlike $\Rightarrow$ $\gamma(\lambda)=X_x+\lambda T_x$.
The above expressions show that, in the spacelike and lightlike cases, the last component $\gamma(\lambda)_4$ of $\gamma(\lambda)$ never goes to zero, which implies by continuity that the time coordinate $t$ stays within the interval $(-\frac{\pi}{2},\frac{\pi}{2})$, hence the lift of $\gamma$ to $AdS_4$ stays within a single copy of its fundamental domain. One also sees that the spatial components (1,2,3) of $\gamma(\lambda)$ go to infinity as $\lambda\rightarrow\pm\infty$, hence $u\rightarrow 0$ along these geodesics as $\lambda\rightarrow\pm\infty$. In the timelike case, the whole time interval $[0,2\pi]$ is spanned by $\gamma(\lambda)$ as $\lambda$ spans the interval $[0,2\pi]$. Since the curve is closed, its lift to $AdS_4$ spans the whole time line $\mathbb{R}$ as $\lambda$ does so. On the other hand, it's clear that in this case the spatial components of $\gamma(\lambda)$ just keep oscillating within a bounded interval of the coordinate $r$ - hence, the coordinate $u$ stays bounded away from zero. Therefore, a timelike geodesic $\gamma$ never escapes to conformal infinity.
