Precise computation of thermodynamic properties of natural gas requires applying an accurate Equation of State (EOS) along with component analysis of natural gas. This procedure is a time consuming process which demands expensive apparatuses. In this study, a rigorous model by firstly applying a gross standard EOS called AGA8-GCM is used to predict compressibility factor (Z) and density (ρ). The required input data for this model are operating temperature and pressure, specific gravity at reference conditions, and merely the amount of nitrogen and carbon dioxide in the natural gas. Then, Genetic Algorithm (GA) is employed to estimate component analysis to increase the accuracy of following calculations. Finally, the thermodynamic properties equations, in which partial derivatives were obtained by AGA8-GCM equation is used to calculate other thermodynamic properties: the speed of sound (usound), heat capacity at constant volume and pressure (CV, CP), Joule-Thomson coefficient (μJT), isentropic exponent (κ), internal energy (U), enthalpy (H), and entropy (S). In order to validate the model, the calculated data were compared to experimental ones collected from literature by the Average Absolute Deviation percentage (AAD%). The AAD% amounts for Z, ρ, usound, CV, CP, μJT, U, H, and S were obtained 0.025%, 0.063%, 0.51%, 0.94%, 1.22%, 3.06%, 0.064%, 0.62% and 2.46%, respectively. The pure methane data is utilized for CV, U, H, and S evaluation as for the lack of empirical data. Furthermore, comparing calculated data with that of AGA8-DCM (Detail Characterization Method), requiring detailed composition of natural gas, shows reliability of this model in the custody transfer region. Investigating the model uncertainty showed a figure of about ±0.1% for both the compressibility factor and density. However, as for the other thermodynamic properties, this figure was usually higher (<0.8%).