Abstract: Neumann Expansions Monte-Carlo stochastic finite element method(NMCSFEM) for ground water models is discussed.The computing methods of the mean value、the variance and probability in some interval of water head at each node are given.Neumann expansions is introduced into the research field of groudwater to improve computational efficiency.When number of input stochastic variables is large,and the variances of input stochastic variables are big,NMCSFEM is best choice.The 2-D confined groundwater flow example (with an analytical solution) is used for stochastic numerical test.