A double multiple-relaxation-time lattice Boltzmann (MRT LB) model for the axisymmetric thermal flow in porous media at the representative element volume, which uses the Darcy-Brinkmann-Forchheimer equation to describe the flow in porous region, is presented. The present model utilizes non-orthogonal transformation matrixes instead of the classical orthogonal ones. Validations are carried out in term of two numerical tests, including natural convection in a vertical annulus without and with a porous medium, respectively, and the results show good agreement with the previous studies. (C) 2017 Elsevier Ltd. All rights reserved.