This paper proposes the modeling and analysis of multi-energy systems as multilayer networks. The aim is to assess the interdependence between different energy infrastructures. Multilayer network modeling enhances the one-dimensional graph-based approach employed to study the vulnerability and the topological characteristics of power grids. The centrality indices defined for one-dimensional networks are extended to multilayer networks, highlighting the differences emerging when interdependence is considered. The proposed multilayer framework is demonstrated on the 24-bus IEEE power grid coupled with a gas network with 7 nodes.