With the energy crisis and ecological environment deterioration, porous thermal-insulating materials become an advanced research hotspot, and the influence of pore distribution cannot be ignored. The mathematical model is established basing on the heat transfer theory, regarding the minimum heat flux density as the objective function, the constant total porosity as a constraint condition, using the BFGS method to optimize the pore distribution. The results show that when the heat flux is the minimum, in the case of the fixed total porosity, the high temperature zone has high porosity, the low temperature zone has low porosity; the maximal fluctuating amplitude of porosity between the adjacent discrete points has great impact on the thermal insulating performance, the greater the fluctuating amplitude, the better the thermal insulating ability. After calculating the temperature field of the corresponding physical model, it can be found that the temperature gradient is non-uniform, the temperature gradient of the high temperature zone is steep, and that of the low temperature zone is gentle. These results have guiding significance for preparation of porous thermal-insulating materials.