The subsurface of Groningen gas field composed of several faults. The continuous production has been resulting in several micro-seismicity activities, particularly for the past decades. One of the reasons for the production-induced-seismicity is fault reactivation at depth. The fault that initially in the non-active state becomes active and starts slipping due to pore fluid extraction. The modelling of fault reactivation induced by production can be simulated by coupling of flow and geomechanics. The project studies the coupling of flow and geomechanics for an idealistic subsurface model with fault. The flow occurs due to pore pressure depletion and is predicted by solving mass conservation. The corresponding deformation is estimated by momentum balance equation. The observation will be limited to rock deformation that is quantified by effective stress and displacement on the fault due to a pressure change in the subsurface. The simulation is performed using Stanford's Automatic Differentiation General Purpose Reservoir Simulator (ADGPRS). In the proposed model, the flow equation is discretised using finite volume method, and poromechanics is discretised by finite-element Galerkin's approach. Both problems share the same grid model, and there is no error associated with the information exchange. The coupled problem is being solved fully implicitly where flow and mechanics formulation are solved simultaneously. The case study for the simulation is a simplified Groningen subsurface model with a fault zone. The impact of variability in production dynamic on fault reactivation is studied. The variations in production rate and production strategy are applied in different formation configurations. It is found that the correlation between the production rates and the stress state in the subsurface depends on the formation offset. However, the stress field does not depend on continuity of the production scheme. Several aspects including the presence of Gas-Water contact in the model was studied in this project as well.