According to the blueprint idea which was presented by Freeze and Harlan, the TOPMODEL was improved by using Lattice Boltzmann method(LBM). The method of LBM was adopted to build a numerical model of watershed concentration which was used to solve the movement process of Richards in unsaturated regions. And the Darcy’s equation was used to set up a numerical model of soil water movement in saturation regions. Then a grid distributed model based on LBMGTOPMODEL was built. In the grid distributed model, the runoff processes was divided into precipitation module, evapotranspiration module, unsaturated water vertical movement module, slope concentration module, river/channel module, and statured soil water/unground water module. The storedfull runoff generation pattern and the runoff generation pattern over infiltration were integrated together with considering runoff yield from soil water profiles, soil anisotropy and land slope. The concentration method was taken into ideas from geomorphologic hydrology to find the quantity relationships between the hydrological process and catchment geomorphology. The slope water movement was described through the integration of geomorphology factors and hydrodynamic diffusion, which solved the distribution of water quantity of slope flow. The Zhongtang Basin was taken as the object to apply hydrological simulation. The results showed that the certainty factor was between 0.547~0.883 and the average value was 0.725. This result could provide spatial distribution situation of water data for calculation of slope flow, channel flow and water movement in root area and unsaturated area.