This paper presents a new boundary integral formulation of the heat energy calculation of steady state heat conduction in heterogeneous media using two methods, i.e. one which is based on some formulations from Hatta and Taya (1986) who investigated equivalent inclusion method, and another one which is closely similar to the method from Christensen (1979) who studied elastic inclusion problem. For the generalized self-consistent scheme (GSCS) model, the corresponding heat energy formulation is also given. Based on the obtained interface heat energy formulation, the boundary element method together with the Maxwell homogeneous scheme are used to calculate the effective heat conduction of 2D and 3D heterogeneous steady state heat conduction media. (C) 2015 Elsevier Ltd. All rights reserved.