Summary The Central Cordillera of Colombia volcano-tectonic province has suitable geological conditions for geothermal renewable and sustainable energy sources. The use of 3D geologic modeling and simulations are an effective tool for better understanding the subsurface geology and have a better estimation of potential resource assessment. This work aims to characterize and simulate the geothermal and hydrothermal processes in a volcano-hosted system, using a 3D crustal-basin modelling software, to estimate the in-place geothermal energy potential in a high-enthalpy deep geothermal system, establishing spatial relationships between heat sources, flow patterns, lithologies and fault networks. The results confirm the large convective-conductive geothermal-hydrothermal system activity in the metamorphic Cajamarca complex and the andesitic-dacitic lavas within the upflow zone, and the hydrothermal flow pattern in the outflow zone structurally controlled by the fault network. Furthermore, the simulated heat energy in-place results served as input data for a reservoir pre-feasibility energy recovery simulation with closed-loop systems, to investigate the ultimate recoverable heat energy potential exploitable as an EGS.