Think this question in a geometry problem. So we know that $x^2+y^2+z^2=1$ is a sphere in 3D space, and $x+y+z=1$ is a plane. The solutions is the crosssection region of those two objects, which is a curve. So there are infinite solutions.
Maybe the above argument is not that simple to visualise (sorry if you cannot). In order to prove the above argument, here we do it in a more direct way.
Now we solve the system equation. Easily we know $y=1-z-x$, so we can have $x^2+(1-z-x)^2+z^2+1=0$. Wiht little bit algebra, we can know that
$x^2-(1-z)x+z^2-1=0$, solve $x$ out treating $z$ as a free variable, then we can know $x=\frac{1-z\pm\sqrt{(z-1)(-3z-1)}}{2}$, then in order to get the rational solution, we have to have $(z-1)(-3z-1)=p^2$ where $p\in Q$ and also $z^2<=1$, it is easy to see that the there is a specific case meeting this requirement where $(z-1)=1/q^2(-3z-1)$, $q\in Z$. For instance, so $q=1$ will give you $(1,0,0)$, $q=2$ will give you $(3/7,-2/7,6/7)$, $q=3$ will give you $(-1/3,2/3,2/3)$. Following this logic, keep going on, you can have infinite solutions.