Geological structures in the subsurface ranging from fractures to reservoirs can be simplified as ellipsoidal inhomogeneities. For instance, one can model a reservoir as an inclusion by considering possibly different material properties and different fluid pressure in comparison to the surrounding rock. Hence, the stresses and displacements associated with the fluid withdrawal from or fluid injection into the formations can be determined by assuming no hydraulic communication between the inclusion and the surrounding medium. The lack of hydraulic communication could be the result of a cap rock or an impermeable seal/fault. On other hand, in the case of fractures, this assumption in not valid anymore and the hydraulic communication between the inclusion and medium should be considered in the calculations. This paper provides analytical solutions for deformation and stress distribution inside and outside of poroelastic ellipsoidal inclusions. Eshelby theory with the Biot theory of poroelasticity are combined to model the change in stress caused by changes in pore pressure or temperature inside the inclusions. Using the provided analytical solutions, we explore the effect of inclusion size, material properties and pressure/temperature condition. The results confirm that neglecting hydraulic communication between the inclusion and the sounding matrix may result in lower inclusion volume change ratio associated with mode (1) and higher inclusion volume change ratio associated with mode (2).