From a visual point of view, volumetric information about reservoir portioning and communication such as sweep, flow patterns, and drainage zones are longer better interpreted and pictured when presented by an average volumetric flux calculation. To this hand, finite volume discretization can be used to substitute streamline simulation-based finite difference to assess flow diagnostic information. Herein, we use finite volume-based flow diagnostics to optimize waterflooding. In particular, we discretize in finite volume the flow equation from Darcy's law single-phase incompressible flow steady state combining with two auxiliary flow equations, time of flight and stationary tracers using the two-point flux approximation to describe fluid particles motion and flow lines. In addition, with the estimation of dynamic heterogeneity, we compute the Lorenz coefficient to highlight the reservoir flow and storage capacity characterization. To optimize waterflooding rates, we first, use an objective function the equalized Lorenz coefficient got through the evaluation of average travel time in cells to increase sweep efficiency and decrease the dynamic heterogeneity coefficient. Second, following the same target, we use the flow diagnostic interactive tools to study the volumetric sweep displacement front and harmonize the flooding breakthrough. In this work, our conceptual approach is to see the reservoir initially filled with oil; then, optimizing the Lorenz coefficient leads us to an oil recovery improvement. To be pragmatic, we apply our waterflooding performance optimization model on two case studies, the ninth SPE comparative solution project, a reexamination of black-oil (synthetic case) and ZHNBA Chinese oilfield (real field dataset).