This system is well studied. You can find full description of solutions in 
"Introduction to the theory of numbers" by Leonard E. Dickson. (See Theorem 47). 

If all the variables are between $1$ and $P$ then the number of solutions is
$$\frac{18}{\pi^2}P^3\log P+O(P^3),$$ 
see "An asymptotic formula for the number of solutions of a system of equations"
N.N. Rogovskaya - Diophantine approximations, Part II (Russian), 1986. She used more natural parametrization: if $a_1=x-u$, $a_2=y-v$, $a_3=z-w$, $b_1=x+u$, $b_2=y+v$, $b_3=z+w$, then the system is equivalent to
$$a_1+a_2+a_3=0,\qquad a_1b_1+a_2b_2+a_3b_3=0.$$ All solutions of the equation
$$kb_1+lb_2-(k+l)b_3=0$$
for $(k,l)=1$ are
$$b_1=b_3-lt,\quad b_2=b_3+kt\quad(b_3,t\in\mathbb{Z}).$$

Your system can be also studied with circle method. And this is unique case (among Vinogradov's type systems) when trigonometric integal is known explicitly (from Rogovskaya's result). 

**UPD:** This arguments are valid for symmetrical system of equations
$$x^2+y^2+z^2=u^2+v^2+z^2,\qquad x+y+z=u+v+z.$$
The original system will take this form after changing variables $z\to -z$, $v\to -v$.