Geological faults are highly developed in the eastern Liaoning Province in China, where Mesozoic granitic intrusions and Archean and Paleoproterozoic metamorphic rocks are widely distributed. Although the heat flow value in eastern Liaoning Province is generally low, the hot springs are very developed. It is obvious that the faults have significant control over the distribution of hot springs, and traditional methods of spatial data analysis such as WofE (weight of evidence) usually do not take into account the direction of the distribution of geothermal resources in the geothermal forecast process, which seriously affects the accuracy of the prediction results. To overcome the deficiency of the traditional evidence weight method, wherein it does not take the direction of evidence factor into account, this study put forward a combination of the Fry and WofE methods, Fry-WofE, based on geological observation, gravity, remote sensing, and DEM (digital elevation model) multivariate data. This study takes eastern Liaoning Province in China as an example, and the geothermal prospect was predicted respectively by the Fry-WofE and WofE methods from the statistical data on the spatial distribution of the exposed space of geothermal anomalies the surface. The result shows that the Fry-WofE method can achieve better prediction results when comparing the accuracy of these two methods. Based on the results of Fry-WofE prediction and water system extraction, 13 favorable geothermal prospect areas are delineated in eastern Liaoning Province. The Fry-WofE method is effective in study areas where the geothermal distribution area is obviously controlled by the fault. We provide not only a new method for solving the similar issue of geothermal exploration, but also a new insight into the distribution of geothermal resources in Liaoning Province.