Quasi-Wilson nonconforming finite element approximation for a class of nonlinear Sobolev equations is discussed on rectangular meshes. We first prove that this element has two special characters by novel approaches. One is that (del(h)(u - I(h)u), del(h)v(h))(h) may be estimated as order O(h(2)) when u is an element of H-3(Omega), where I(h)u denotes the bilinear interpolation of u, v(h) is a polynomial belongs to quasi-Wilson finite element space and del(h) denotes the piecewise defined gradient operator, h is the mesh size tending to zero. The other is that the consistency error of this element is of order O(h(2))/O(h(3)) in broken H-1-norm, which is one/two order higher than its interpolation error when u is an element of H-3(Omega)/H-4(Omega). Then we derive the optimal order error estimate and superclose property via mean-value method and the known high accuracy result of bilinear element. Furthermore, we deduce the global superconvergence through interpolation post processing technique. At last, an extrapolation result of order O(h(3)), two order higher than traditional error estimate, is obtained by constructing a new suitable extrapolation scheme.