This work focuses on the determination of the design and operation parameters of a thermal system depending on the optimization objective set. Its main objective and contribution concern the proposal of a generalized methodological structure involving multiobjective optimization techniques aimed at providing a solution to a practical problem, such as the design and dimensioning of a solar thermal system. The analysis is based on system operation data provided by a dynamic simulation model, leading to the development of multiple surrogate models of the thermal system. The thermal system surrogate models correlate the desired optimization objectives with thermal system design and operation parameters while additional surrogate models of the Pareto frontiers are generated. The implementation of the methodology is demonstrated through the optimal design and operation parameter dimensioning of a solar-assisted geothermal heat pump that provides domestic hot water loads of an office building. Essentially, energy consumption is optimized for a desired domestic hot water thermal load coverage. Implementation of reverse-engineering methods allows the determination of the system parameters representing the optimized criteria.