The results of core displacement experiments show that increasing the water drive velocity when it is bigger than the limit value can effectively reduce the residual oil saturation and improve the oil displacement efficiency under the same PV. However, the existing commercial simulators (Eclipse, CMG et al.) cannot simulate the effect of water velocity on the relative permeability curve in the process of numerical simulation.In this article, capillary number (Ca), defined as the dimensionless ratio of viscous force to capillary force, is used to characterize the relationship between water drive velocity and residual oil. Second, a new Boltzmann (BG) equation is proposed to match the nonlinear relationship between Ca and residual oil. The BG equation is a continuous function, which is very beneficial to the stability of numerical calculation. Finally, a new reservoir numerical simulator is established which captures the dynamic variation of residual oil saturation with water drive velocity in a water flooding reservoir based on the black oil model. The new simulator was verified by comparing it with the commercial reservoir simulator ECLIPSE and experimental data. The simulation results show that compared with the common model, the model considering the dynamic variation of residual oil saturation with water drive velocity reduced the residual oil saturation near the main flow line after enhanced injection rate. The oil phase flow capacity in the model is enhanced, the water cut is decreased, and the oil recovery rate is higher. The history matching of the S oilfield in Bohai Bay is achieved with the new simulator, and the history matching accuracy is obviously higher than that of Eclipse. The findings of this study can help with a better understanding of the distribution law and flow law of remaining oil in the high water cut stage of the reservoir and have good theory and application value for water flooding offshore oilfields.