It very much depends on what you consider a "derivation" and that you consider an "easy" one. You can start assuming that the energy density depends on $(\kappa_1^2+\kappa_2^2)g^{\frac 12}$, where $\kappa_i$ are the principal curvatures, which can be considered quite natural from a physical point of view (in one dimension this corresponds to assuming that the energy density of a 1D elastic continuum depends quadratically on the curvature, and it leads to the non-linear theory of the Euler Elastica). Then you can deduce your integral functional as a linearized form of the problem, as is done in most textbooks (and also sketched here: https://en.wikipedia.org/wiki/Thin_plate_energy_functional).
If instead you want more something like a "microscopic" derivation, then you'll need more sophisticated tools to make sense of the limit functional of a suitable sequence of 3D Cauchy-continuum energy functionals, when one of the dimensions goes to zero (and the coefficients are rescaled accordingly). A standard reference (in my view a very good one, as it covers the mathematical aspects in full rigor) is: Ciarlet, P. G. (1997). Mathematical Elasticity: Volume II: Theory of Plates. Elsevier.
Hope this helps.
Btw: I would not say "Given a thin plate $z=z(x,y)$", as $z$ is in fact the placement function.