We are devoted to the study of the following sub-Laplacian system with Hardy-type potentials and critical nonlinearities

$$\{\begin{array}{r}-{\mathrm{\Delta }}_{\mathbb{G}}u-{\mu }_1\frac{{\psi }^{2}u}{\text{d}(z{)}^2}={\lambda }_1\frac{{\psi }^{\alpha }|u{|}^{2^{\ast }(\alpha )-2}u}{\text{d}(z{)}^{\alpha }}+\beta p_{1}f(z)\frac{{\psi }^{\gamma }|u{|}^{p_1-2}u|v{|}^{p_2}}{\text{d}(z{)}^{\gamma }}\text{in}\mathbb{G},\\ -{\mathrm{\Delta }}_{\mathbb{G}}v-{\mu }_2\frac{{\psi }^{2}v}{\text{d}(z{)}^2}={\lambda }_2\frac{{\psi }^{\alpha }|v{|}^{2^{\ast }(\alpha )-2}v}{\text{d}(z{)}^{\alpha }}+\beta p_{2}f(z)\frac{{\psi }^{\gamma }|u{|}^{p_1}|v{|}^{p_2-2}v}{\text{d}(z{)}^{\gamma }}\text{in}\mathbb{G},\end{array}$$

where $-{\mathrm{\Delta }}_{\mathbb{G}}$ is the sub-Laplacian on Carnot group $\mathbb{G}$, ${\mu }_1$, ${\mu }_2\in [0,{\mu }_{\mathbb{G}})$, $\alpha ,\gamma \in (0,2)$, ${\lambda }_1$, ${\lambda }_2$, $\beta$, $p_1$, $p_2>0$ with $1<p_1+p_2<2$, $\text{d}(z)$ is the ${\mathrm{\Delta }}_{\mathbb{G}}$-gauge, $\psi =|{\mathrm{\nabla }}_{\mathbb{G}}\text{d}(z)|$, $2^{\ast }(\alpha ):=\frac{2(Q-\alpha )}{Q-2}$ is the critical Sobolev-Hardy exponents, and ${\mu }_{\mathbb{G}}=(\frac{Q-2}{2}{)}^2$ is the best Hardy constant on $\mathbb{G}$. By combining a variant of the symmetric mountain pass theorem with the genus theory, we prove the existence of infinitely many weak solutions whose energy tends to zero when $\beta$ or ${\lambda }_1$, ${\lambda }_2$ belong to a suitable range.