In the Athabasca Basin of Alberta, Canada, the Lower Cretaceous McMurray Formation reservoir is overlain by the Clearwater Formation, a regionally continuous layer composed predominantly of shales with interbedded mudstones. The shales form the reservoir caprock and have the role of blocking the vertical migration of the steam from thermal oil production by confining the stresses and the deformations. We have developed a new method for the 3D seismic reconstruction of the anisotropic stress field, which accounts for the formation pore pressure and the effective stress. We integrated anisotropy estimated from dipole sonic logs with formation pore pressure data from piezometers and elastic properties obtained from multicomponent seismic inversion. The method combines the Terzaghi effective stress with the Schoenberg and Sayers elastic stiffness matrix for horizontal transversely isotropic (HTI) fractured materials. The key points of this method are the estimation of the formation pore pressure in the abnormal regime of the Clearwater Formation, the normal fracture weakness parameter (based on constraints on the compressional velocity of the intact rock under the HTI assumption), and the 3D seismic anisotropic stress field. We expressed the total vertical stress as the weight of the overlying formations, and the total minimum and maximum horizontal stresses as a combination of the total vertical stress, normal fracture weakness, formation pore pressure, Biot-Willis coefficient, and Lamé elastic constants. The effective principal stresses are estimated from the equivalent total principal stresses and the formation pore pressure multiplied by the Biot-Willis coefficient. We observed excellent consistency between the calculated total minimum horizontal stress and mini-frac values. This new method for the 3D seismic reconstruction of the anisotropic stress field allows for the assessment of the caprock integrity and for operational savings based on a reduced number of mini-frac measurements, and it can be used for time-lapse stress estimation within the thermal production reservoir.