In this paper, an orthogonalized isoparametric local average model for 2-D continuative stationary random fields is proposed. This model is used in Neumann stochastic finite element method. Numerical example shows the accuracy and the efficiency of the method.