Activated sludge process is an important stage in Wastewater Treatment Plant (WWTP). In this study, model of the activated sludge process from Bunus Regional Sewage Treatment Plant Kuala Lumpur, Malaysia is developed. This paper focuses on modelling and model reduction of the WWTP system. The model with best fits of higher than 80% and the order of less than 10 is selected. For modelling purposes, data obtained is stimulated and modelled using System Identification technique which employ linear model ARX. For model reduction purposes, the high order model is reduced using model order reduction (MOR) of a combination of Singular Perturbation Approximation (SPA) and Frequency Domain Gramian based Model Reduction (FDIG) method. From the modelling results obtained, the ARX model with best fit of 85.11% is selected. Meanwhile, for the MOR FD-SPA technique, a 9th order model is selected with  reduction error between frequencies 0.05 rad/s and 1.4 rad/s.