The suggestions by Issam are good. Some simplifications:
* When ∂²w/∂x∂z are eliminated from the last equation, by Issam's Step (1), all references to u are also eliminated.
* So you end up with a first order differential equation (in z) for the variable V == ∂v/∂x. Hence v can be "found" by integrations.
* With this solution inserted into the middle equation, you end up with a first order differential equation (in x) for the variable u, provided a2 - a1 is non-zero (otherwise a consistency condition emerges).
But I am not convinced that the formal solution thus obtained will be immediately useful, since the integrals will involve unknown integration "constants" that can depend on the remaining variables. The nature of the underlying problem must provide further information.
Since we have no information about what system you're trying to model, this suggestion may not apply perfectly, but I'd recommend you to do a nondimensionalization of the system, if possible. It would reduce the number of parameters and make it easier to find exact solutions.