An oil and gas field requires careful operational planning and management via production optimisation for increased recovery and long-term project profitability. This paper addresses the challenge of production optimisation in a field undergoing secondary recovery by water flooding. The field operates with limited processing capacity at the surface separators, pipeline pressure constraints and water injection constraints; an economic indicator (Net Present Value -NPV) is used as the objective function. The formulated optimisation framework adequately integrates slow-paced subsurface dynamics using reservoir simulation and the fast-paced surface dynamics using sophisticated multiphase flow simulation in the upstream facilities. Optimisation of this holistic longterm model is made possible by developing accurate second order polynomial proxy models at each time step. The resulting formulation is solved as a Nonlinear Program (NLP) using commercially available solvers. A comparative analysis is performed using MATLAB's fmincon solver and the IPOPT solver for their robustness, speed and convergence stability in solving the proposed problem.By implementing 2 synthetic case studies, our mathematical programming approach determines the optimal production and injection rates of all wells and further demonstrates considerable improvement to the NPV obtained by simultaneously applying the tools of streamline, reservoir and surface facility simulation for well rate allocation.