LAUSR

This research presents a 2D conductive thermal representation for geothermal systems through the development of a numerical model based on heat transfer equations using Fortran language and applying the Tri-Diagonal Matrix Algorithm as a numerical solver of the equations. This numerical model simulates the heat transfer between two geothermal wellbores, used as boundary conditions, and calculates the temperature distribution. In addition, the geometry and composition of the model are based on the stratigraphic columns of the wellbores. The effect of the scarcity of geothermal gradient measurements in geothermal prospects is reduced by using a 2D model that requires only information in two wellbores. After applying the numerical model to a Hot Dry Rock zone and a sedimentary basin, it was found that the validation and accuracy estimation of the temperature profiles showed a good correlation and coherence between the measured and the simulated values. These results suggested the use of this numerical model for the reliable description of the thermal state of the upper section of the geothermal prospect with a high correlation with the observed measurements. This model would be very useful for defining the isotherms shape when planning the exploitation of the geothermal resource. Finally, the practical use of this model was highlighted because it only requires the thermal logs and thermal conductivity configuration of only two exploration wellbores as the main input data.