In practical applications, rapid prediction and optimization of heat transfer performance are essential for premixed methane impinging flame jets (PMIFJs). This study uses computational fluid dynamics (CFD) combined with a methane detailed chemical reaction mechanism (GRI–Mech 3.0) to study the equivalence ratio (ϕ), Reynolds number (Re) of the mixture, and the normalized nozzle–to–plate distance (H/d) on the heat transfer performance of PMIFJs. Moreover, the Kriging model (KM) was used to construct a prediction model of PMIFJ heat transfer performance. A genetic algorithm (GA) was used to determine the maximum likelihood function (MLE) of the model parameters for constructing KM and identify the points with the maximum root mean square error (RMSE) as the new infilled points for surrogate–based optimization (SBO). Combining these methods to analyze the simulation results, the results show that the global heat transfer performance of PMIFJs is enhanced with the increase in ϕ, the increase in Re, and the decrease in H/d. Sensitivity analysis points out that Re and ϕ significantly affect enhanced heat transfer, while H/d has a relatively small effect. In addition, GA was also used to search for the optimal heat transfer performance, and the global heat transfer performance at specific conditions was significantly enhanced. This study deepens the understanding of the heat transfer mechanism of impinging flame jets and provides an efficient method framework for practical applications.