Based on the Preissmann implicit scheme for the one-dimensional Saint-Venant equation, the mathematical model for one-dimensional river networks and canal networks was developed and the key issues on the model were expatiated particularly in this article. This model applies the method of three-steps solution for channel-junction-channel to simulate the river networks, and the Gauss elimination method was used to calculate the sparse matrix. This model was applied to simulate the tree-type irrigation canal networks, complex looped channel networks and the Lower Columbia Slough networks. The results of water level and discharge agree with the data from the Adlul and field data. The model is proved to be robust for simulating unsteady flows in river networks with various degrees of complex structure. The calculated results show that this model is useful for engineering applications in complicated river networks. Future research was recommended to focus on setting up ecological numerical model of water quality in river networks and canal networks.