The larger spatial variability of soil properties is due to the difference of soil properties for different soils. So it is necessary to establish a model to show the actual mechanism of the variability. The multi-indicator geology statistical model based on the Markov field is introduced. And the analytical expression for the statistics of the soil parameters is derived and the impact of multi-soil on K-L （Karhunen-Loeve） expansion is analyzed. Then the stochastic seepage simulation is conducted with spectral method. The results show the correctness of Gaussian stationary field-Markov field and the validity of spectral model.