For a different viewpoint from the excellent treatments by Scott and Thurston of 3-dimensional geometries, if you are trying to get a feel for the homogeneous Riemannian $3$-manifolds (which, as noted, were first classified by Bianchi), you might want to try looking at them from the point of view of their most basic invariants, their curvatures. This provides a natural classification and made it easy for me to organize the information when I was learning the subject. From this point of view, it doesn't quite divide the way that Bianchi did it, but you shouldn't have any trouble making the comparison. In outline, it goes like this:
Let $(M,g)$ be a connected Riemannian $3$-manifold that is homogeneous, i.e., the group $G$ of isometries of $M$ acts transitively on $M$. Consider the Ricci curvature $\mathrm{Ric}(g)$. This is also a quadratic form on $M$ and, as such, relative to $g$, it has $3$ real eigenvalues $\lambda_1$, $\lambda_2$, and $\lambda_3$, which are constant functions on $M$ because of the homogeneity hypothesis.
Case 1: $\lambda_1=\lambda_2=\lambda_3=\lambda$. In this case, the sectional curvature of $g$ is constant, and so $M$ is a space form, of elliptic, parabolic, or hyperbolic type, depending on whether $\lambda$ is positive, negative, or zero. These are easily visualized, as models for them are straightforward to construct.
Case 2: $\lambda_1\not=\lambda_2=\lambda_3$. In this case, there is a well-defined $G$-invariant line field $L$ on $M$ which is the eigendirection of multiplicity $1$ for $\mathrm{Ric}(g)$. After passing to a double cover, if necessary, there is a unit vector field $Y$ tangent to $L$ and one can compute the Lie derivative of $g$ in the direction $Y$. There are two subcases:
Case 2a: $L_Yg=0$ (i.e., $Y$ is a Killing field). In this case, the isometry group of $g$ has dimension $4$, the integral curves of $Y$ are geodesics and the set of these geodesics is a surface $S$. The projection $M\to S$ with fibers tangent to $L$ is a Riemannian submersion, and the induced metric on $S$ has constant Gauss curvature. Then $M$ has the structure of a principal bundle over $S$ (with group action generated by $Y$) with a connection whose curvature $2$-form is a constant multiple of the area form on $S$. There is a $2$-parameter family of such metrics.
Case 2b: $L_Yg$ is nonzero, in which case, it turns out that there is a coframing $\omega_1$, $\omega_2$, and $\omega_3$ defined up to sign on $M$ so that $g = {\omega_1}^2+{\omega_2}^2+{\omega_3}^2$, while $\mathrm{Ric}(g) = \lambda_1\ {\omega_1}^2$ (i.e., $\lambda_2=0$) and $L_Yg = \mu\ ({\omega_2}^2-{\omega_3}^2)$ for some constant $\mu>0$. This coframe field is invariant under $G$ (at least, up to sign), so $G$ has dimension $3$ and, essentially, the metric $g$ is a left-invariant metric on $G$, which is a covering of $M$. Thus, this subcase, which consists of a $2$-parameter family of metrics (parametrized by $\lambda_1\not=0$ and $\mu>0$) can be covered as a limit of Case 3, which consists of left invariant metrics on $3$-dimensional Lie groups.
Case 3: The $\lambda_i$ are all distinct, in which case, after passing to a cover and restricting $G$ to its identity component if necessary, one sees that $G$ leaves fixed a homogeneous coframing on $M$, so $g$ is essentially a left-invariant metric on $G$. By the usual algebra that classifies the $3$-dimensional Lie groups, one knows that there will be a $g$-orthogonal $G$-invariant coframing $\omega_1$, $\omega_2$, $\omega_3$ so that $g = {\omega_1}^2+{\omega_2}^2+{\omega_3}^2$ and so that
\begin{aligned}
d\omega_1 &= a_1\ \omega_2\wedge\omega_3\\\\
d\omega_2 &= a_2\ \omega_3\wedge\omega_1 + b\ \omega_1\wedge\omega_2\\\\
d\omega_3 &= a_3\ \omega_1\wedge\omega_2 - b\ \omega_3\wedge\omega_1
\end{aligned}
where $a_1$, $a_2$, $a_3$, and $b$ are constants satisfying either $a_1=0$ or $b=0$ (these latter alternatives are necessary in order for the Jacobi identity to be satisfied). If $b=0$, then the condition that the metric $g$ so defined on the corresponding Lie group $G$ have distinct eigenvalues of its Ricci curvature is just that the $a_i$ be distinct and that $a_1{+}a_2{+}a_3\not=2a_i$ for any $i$. If $a_1=0$, then the condition that the metric so defined have all Ricci eigenvalues be distinct is that $a_2\not=a_3$ and $b^2+a_2a_3\not=0$. Thus, there are two $3$-parameter families of these metrics in Case 3. Meanwhile, when one removes some of the inequalities, one gets limiting cases that cover Case 2b.
Thus, in all the homogeneous cases, the curvature (together with, sometimes, derivatives of curvature eigendirections) actually produces a natural coframing and/or foliation of the metric $g$, and, because these invariants are so natural, they are a good place to start to study the local geometry of the homogeneous Riemannian $3$-manifolds.
