Since its introduction 30 years ago, the lattice Boltzmann (LB) method has achieved great success in simulating fluid flows and modeling physics in fluids. Owing to its kinetic nature, the LB method has the capability to incorporate the essential microscopic or mesoscopic physics, and it is particularly successful in modeling transport phenomena involving complex boundaries and interfacial dynamics.The LB method can be considered to be an efficient numerical tool for fluid flow and heat transfer in porous media. Moreover, since the LB method is inherently transient, it is especially useful for investigating transient solid-liquid phase-change processes wherein the interfacial behaviors are very important. In this article, a comprehensive review of the LB methods for single-phase and solid-liquid phase-change heat transfer in porous media at both the pore scale and representative elementary volume (REV) scale. The review first introduces the fundamentals of the LB method for fluid flow and heat transfer. Then the REV-scale LB method for fluid flow and single-phase heat transfer in porous media, and the LB method for solid-liquid phase-change heat transfer, are described. Some applications 2 of the LB methods for single-phase and solid-liquid phase-change heat transfer in porous media are provided. In addition, applications of the LB method to predict effective thermal conductivity of porous materials are also provided. Finally, further developments of the LB method in the related areas are discussed. Keyword: Lattice Boltzmann method; Single-phase heat transfer; Solid-liquid phase change; Porous media proposed the RLBE model with an enhanced collision operator, which was shown to be linearly stable. Based on the Bhatnagar-Gross-Krook (BGK) collision operator [104], the LB model using a single relaxation time, referred to as the BGK-LB model, has been independently proposed by Chen et al. [14], Koelman [15], and Qian et al. [16]. In the BGK-LB model, the equilibrium distribution function is chosen to recover the incompressible N-S equations in the incompressible limit. Due to its high efficiency and extreme simplicity, the BGK-LB model has become the most popular LB model, in spite of some well-known deficiencies (e.g., numerical instability, viscosity dependence of boundary 9 locations). Almost at the same time when the BGK-LB method was developed, d'Humières [17] proposed the multiple-relaxation-time (MRT) LB method. The MRT-LB equation (also referred to as generalized lattice Boltzmann equation (GLBE)) is an important extension of the RLBE proposed by Higuera et al. [12,13]. This work is significant because, in addition to it overcomes some deficiencies of the BGK-LB model, it facilitates the extension of the LB method, expanding the applications to a much wider range of phenomena and processes.In 1997, a direct connection between the LB equation and the continuous Boltzmann equation in kinetic theory was rigorously established by He and Luo [19,20]. In particular, it has been demonstrated that the LB equation c...