This paper presents a multi-objective economic-environmental dispatch (MOEED) model for integrated thermal, natural gas, and renewable energy systems considering both pollutant emission levels and total fuel or generation cost aspects. Two cases are carried out with the IEEE 30-bus system by replacing thermal generation units into natural gas units to minimize the amount of toxin emission and fuel cost. Equality, inequality like active, reactive powers, prohibited operating zones (POZs) which represents poor operation in the generation cost function, and security constraints are considered as system constraints. Natural gas units (NGUs) are modeled in detail. Therefore, the flow velocity of gas and pressure pipelines are also considered as system constraints. Multi-objective optimization algorithms, namely multi-objective Harris hawks optimization (MOHHO) and multi-objective flower pollination algorithm (MOFPA) are employed to find Pareto optimal solutions of fuel or generation cost and emission together. Furthermore, the technique for order preference by similarity to ideal solution (TOPSIS) is proposed to obtain the best value of Pareto optimal solutions. Three scenarios are investigated to validate the effectiveness of the proposed model applied to the IEEE 30-bus system with the integration of variable renewable energy sources (VRESs) and natural gas units. The results obtained from Scenario III with NGUs installed instead of two thermal units reveal that the economic dispatching approach presented in this work can greatly minimize emission levels as 0.421 t/h and achieve lower fuel cost as 796.35 $/h. Finally, the results obtained show that the MOHHO outperforms the MOFPA in solving the MOEED problem.