Mathematical modelling and simulation of selected stratigrapy within the province of Saskatchewan, Canada for geothermal energy prospecting