Due to the symmetry of the paraboloid about the z-axis, you can treat this is a surface of revolution. Consider the curve
, with
, and revolve it about the y-axis. The area of the resulting surface is then

But perhaps you'd like the surface integral treatment. Parameterize the surface by

with
and
, where the third component follows from

Take the normal vector to the surface to be

The precise order of the partial derivatives doesn't matter, because we're ultimately interested in the magnitude of the cross product:

Then the area of the surface is

which reduces to the integral used in the surface-of-revolution setup.