In this research paper, comprehensive thermodynamic modeling of an integrated energy system consisting of a multi-effect desalination system, geothermal energy system, and hydrogen production unit is considered and the system performance is investigated. The system's primary fuel is a geothermal two-phase flow. The system consists of a single flash steam-based power system, ORC, double effect water–lithium bromide absorption cooling system, PEM electrolyzer, and MED with six effects. The effect of numerous design parameters such as geothermal temperature and pressure on the net power of steam turbine and ORC cycle, the cooling capacity of an absorption chiller, the amount of produced hydrogen in PEM electrolyzer, the mass flow rate of distillate water from MED and the total cost rate of the system are studied. The simulation is carried out by both EES and Matlab software. The results indicate the key role of geothermal temperature and show that both total exergy efficiency and total cost rate of the system elevate with increasing geothermal temperature. Also, the impact of changing absorption chiller parameters like evaporator and absorber temperatures on the COP and GOR of the system is investigated. Since some of these parameters have various effects on cost and efficiency as objective functions, a multi-objective optimization is applied based on a Genetic algorithm for this system and a Pareto-Frontier diagram is presented. The results show that geothermal main temperature has a significant effect on both system exergy efficiency and cost of the system. An increase in this temperature from 260 C to 300 C can increase the exergy efficiency of the system for an average of 12% at various working pressure and also increase the cost of the system by 13%.