Oil sand reservoirs and caprock undergo deformations triggered by pore pressure increases and thermal induced stresses during the steam-assisted-gravity-drainage (SAGD) processes. Geomechanical assessments are mandated by energy regulators to evaluate the caprock integrity and ensure the safe SAGD operations. Commercial reservoir simulation packages started to incorporate geomechanical effects when predicting flow response; however, these geomechanical modules are not able to correctly model the plastic deformations caused by thermal-hydraulic-mechanical (THM) interactions, which has a first order effect on predicting steam chamber propagation and evaluating caprock integrity.
 An integrated coupled THM modeling methodology is proposed here to improve the modeling of reservoir deformations and caprock integrity in a heterogeneous oil sand reservoir with interbedded shale barriers. The pressure and temperature front are found to propagate at different speed and that dominate the elastic and plastic deformations caused by changes of shear and mean effective stress. Therefore, four stages are divided in the SAGD process that can be interpretations of changes in stress paths including buildup of pore pressure, generation and dissipation of thermal induced stresses. The response surfaces of minimum factor of safety (FOS) are introduced and computed to provide a conservative estimate for caprock integrity during SAGD of a heterogeneous reservoir with multiple layers of caprocks in Athabasca oil sands. 
Read full abstract