In this work, numerical optimization based on stochastic gradient methods is used to assist geothermal operators in finding improved field development strategies that are robust to accounted geological uncertainties. Well types, production rate targets and well locations are optimized to maximize the economics of low-enthalpy heat recovery in a real-life case with stacked reservoir formations. Significant improvements are obtained with respect to the strategy designed by engineers. Imposing fault stability constraints impacts significantly the optimal configurations, with coordinated well rates and placement playing a key role to boost efficiency of geothermal production while keeping stress change effects to acceptable limits.