In this paper, we propose a realistic model for gas distribution of an advanced municipal wastewater treatment works and through minimisation of the total cost of gas distribution we perform retrospective optimisation (RO) using historical plant data. This site is the first in the UK with a mixed operational strategy for biomethane produced on site: to burn in combined heat and power (CHP) engines to create electricity, burn in steam boilers for onsite steam use or inject the biomethane into the National Grid. In addition, natural gas can be imported to make up shortfalls in biomethane if required. Implemented using a novel mixed integer linear programming (MILP) approach, to ensure a fast and robust solution, our results indicate the plant operated optimally within accepted tolerance 98% of the time. However, improving plant robustness (such as reducing unexpected breakdown incidents) could yield a significant increase in gas revenue of 7.8%.