The analysis of seismic reflection amplitude variations with incidence angle and azimuth (AVAz) is a powerful tool for studying and characterizing fractured formations. There is evidence to suggest that wave-induced fluid pressure diffusion (FPD) within connected fractures can significantly affect the AVAz response. This, in turn, indicates that this seismic attribute may contain information on the connectivity degree of the fracture network, which represents a key hydraulic property of fractured formations. The assessment of the sensitivity of seismic reflectivity with respect to fracture connectivity accounting for FPD effects to date is limited to 2D models, which is mainly due to the complexity and the associated computational cost of modeling the effective properties of the corresponding generically anisotropic 3D media. To overcome these issues, we use a numerical upscaling procedure based on Biot’s theory of poroelasticity that allows us to obtain effective anisotropic representations of 3D fractured formations accounting for the FPD effects. Using these effective properties, we then perform a plane-wave analysis to compute the reflection coefficients at the top of the corresponding fractured formation. The results of our numerical analyses extend the previous 2D studies and demonstrate that seismic reflection coefficients are highly sensitive to the degree of connectivity of the fracture networks. The reflectivity responses of formations comprising connected and unconnected fractures indicate pronounced differences not only with respect to incidence angle but also with respect to azimuth. Quite importantly, our results also demonstrate that crossplots of commonly used AVAz coefficients permit us to identify regions with higher or lower fracture connectivity.