In this paper, we consider, from a numerical point of view, a two-temperature poro-thermoelastic problem. The model is written as a coupled linear system of hyperbolic and elliptic partial differential equations. An existence result is proved and energy decay properties are recalled. Then we introduce a fully discrete approximation by using the finite element method and the implicit Euler scheme. Some