A novel mixed enriched finite element model is developed for coupled nonâlinear thermoâhydroâmechanical simulation of fractured porous media with threeâphase flow and thermal coupling. Simulation of induced acoustic emission (AE) and microseismic emission (ME) due to tensile fracturing and shear slip instability of preâexisting fracture interfaces is carried out and the numerical results of the emitted signals are analysed. The mathematical model is based on the generalized Biot's theory for coupled interaction of solid and fluid phases. A computationally robust nonâlinear solver is developed to handle the severe nonâlinearities arising from fluid saturations, relative permeabilities of fluids, constitutive models of interfaces and convective thermal coupling. To model preâexisting natural fractures and faults, discrete fracture propagation and nucleation of cracks (microâcracking) independently of the original mesh topology, a local PartitionâofâUnity (PU) finite element method, namely, the Phantom Node Method (PNM) is implemented. The cohesive fracture modelling scheme is implemented to account for the nonâlinear behaviour of fracturing and localization, and to rectify the nonâphysical stress singularity condition at the fracture tip. Effects of different system parameters on fracturing, shearâslip instability and the associated induced AEs and MEs are investigated through various numerical results.