The stochastic flood process model can generate a large number of flood processes with shapes varying but statistically similar to historical flood data. The combined flood control dispatching of reservoir groups requires coordinated flood resource scheduling based on flood processes from various hydrological stations in the basin, aiming to achieve scientific allocation of flood resources. Therefore, there is an urgent need to develop a multi-station stochastic flood process simulation method. In view of this, a multi-station stochastic flood process simulation method based on the Copula function is proposed. This method first uses a two-dimensional Copula function to simulate the flood process of the main station, and then uses a conditional mixing method to construct a three-dimensional Copula function to sequentially simulate the flood processes of the remaining stations. The results of studying eight hydrological stations in the middle and upper reaches of the Yangtze River indicate that the relative error between the simulated flood data and historical flood statistical parameters is small. Except for the maximum values at Beibei and Wulong stations, the relative errors of the statistical parameters at the hydrological stations are about 5%. This reflects the statistical characteristics of the historical flood process and can support work related to flood control planning and design for hydraulic projects, combined dispatching of reservoir groups, and risk analysis of flood control dispatching.