The quick progress in quantum entanglement research allows us not one to study quantum systems down to N-bodies but also to take a new look at these systems in different branches of physics; particularly the statistical thermodynamics where the application of the thermo-field dynamics (<I>TFD</I>) method to the investigation of entanglement is fruitful. Because the traditional methods based on the identification of a specific parameter show their limit. The process using (<I>TFD</I>) facilitates the understanding of entanglement because it focuses directly on the eigenstate of the system and it is useful in the equilibrium and the non-equilibrium states also. In this context, the (<I>TFD</I>) method is used in this paper to analyze entanglement of an electron interacting with a two-mode electromagnetic field assimilated to an electron with two harmonic oscillators. Entanglement entropies are derived between concerned, not concerned harmonic oscillator and electron compute when the system is in a thermodynamic equilibrium and non-equilibrium state. For the equilibrium case, an increase in the number of particles per unit volume increases the quantum entanglement consequently entanglement appears more important for the couple oscillator-electron than the one electron, this trend is reversed for the non-equilibrium case. By respecting the entanglement parameters, such results allow us to know the relative equilibrium state of the overall system.