Ensemble-Based Uncertainty Quantification of Geothermal Reservoir Models Using Randomized Low-Rank Matrix Algorithms and Adjoint Code