Multi-scale approach for modeling the transversely isotropic elastic properties of shale considering multi-inclusions and interfacial transition zone