In this work, a robust numerical tool is developed for the thermo-hydromechanical (THM) modeling of crack propagation in porous media.The flow of fluid in porous media is simulated using the Darcys law by Nonconforming Finite Element Method.To consider the influence of crack on heat transfer, the energy conservation equation is solved for the fluid phase and the solid phase separately and sequentially using the time splitting scheme.In the fluid phase, a combination of Discontinuous Galerkin (DG) and Multi-Point Flux Approximation (MPFA) methods is used to solve the advection-diffusion heat transfer equation.The DG method is applied to solve the advective term of the equation, while the diffusive term is discretized with symmetric MPFA method.This combined numerical method (DG and MPFA) is able to reduce the numerical diffusion and guarantee the stability of the computation.Sequentially, the conductive heat transfer equation in the solid phase is solved using those results obtained from the fluid phase as initial temperature field.The extended Finite Element Method (XFEM) is introduced to solve the cracked media by enriching the fields near the discontinuity (crack surface) and near the singularity (at the crack tip).Then the Stress Intensity Factor (SIF) is computed in the post processing stage using the J-integral technique.Moreover, in order to reduce the computational cost and meanwhile maintain the accuracy of the model, the non iterative scheme with time stepping based on local error control is applied to solve this THM system.Several numerical examples are presented to demonstrate the utility of the proposed model in porous media.