The multi-phase flow characteristics and preferential paths in the heterogeneous pore structures of reservoir rocks at high pressures and temperatures are much significant to evaluate the exploitation efficiency of subsurface hydrocarbon resources and carbon sequestration effects. Therefore, these characteristics and paths should be accurately described and understood. However, only a few existing evaluation models consider the effect of thermo-hydro (T-H) coupling on these characteristics and paths. In this study, based on the commercial programs COMSOL with MATLAB, we developed a T-H coupling flow model to characterize the immiscible two-phase interface dynamics and flow paths in heterogeneous pore structures at various two-phase temperature differences. Ten numerical cases were carried to evaluate the effects of the reservoir and injection temperatures on the preferential paths of two-phase flow in the same two-phase temperature difference range. The results indicate that during CO2-water non-isothermal displacement process, the reservoir temperature affects the preferential paths of two-phase flow and injection temperature has a larger effect on branch flow in small pore regions. The injection rate considerably affects the two-phase flow paths as a result of temperature effects. This study provides a reference for quantitatively analyzing multi-phase flow in reservoir rocks at formation temperatures.