One of the important sources of renewable energy is geothermal heat. Its special feature of being independent 24/7 ensures the stability and security of the system, either for electricity or heat production. Geothermal energy has a local character and is limited by the geological characteristics of each state. In the Republic of Croatia, the development of geothermal energy is closely related to the development of the oil industry, as geothermal deposits were discovered during oil and gas exploration. Considering the established temperature gradients in Croatia, there is a greater possibility of using geothermal energy, and for this, it is necessary to evaluate its full potential and possibilities of use. The aim of this research is to determine the heat potential of the Croatian part of the Pannonian Basin System (CPBS), a part of Croatia with exceptional geothermal potential, based on the analysis of a large amount of well data with confirmed water inflow. In order to estimate the heat in place, the available data on the presence of inflow, temperature, and porosity, as well as permeability and volume for each well/reservoir included in the assessment, were considered. In geothermal reservoirs, one of the most important pieces of data besides petrophysical and thermodynamic data is the potential of the well, i.e., the maximum flow under certain permeability and porosity conditions. To define this, the productivity index was made dependent on the permeability of each well, and the inflow in each well was risked using Monte Carlo for three main geological phases in CPBS, which subsequently influenced inflow and spacing between production and injection wells. The beta-PERT distribution for permeability is used in Monte Carlo simulation to determine the most likely values and produce a distribution that resembles the real probability distribution. As a result, geothermal potential was mapped according to the obtained values of heat in place for part of the CPBS covered with analysed wells.