This study investigated the point-wise superconvergence of block finite elements for the variable coefficient elliptic equation in a regular family of rectangular partitions of the domain in three-dimensional space. Initially, the estimates for the three-dimensional discrete Greens function and discrete derivative Greens function were presented. Subsequently, employing an interpolation operator of projection type, two essential weak estimates were derived, which were crucial for superconvergence analysis. Ultimately, by combining the aforementioned estimates, we achieved superconvergence estimates for the derivatives and function values of the finite element approximation in the point-wise sense of the $L^{\mathrm{\infty }}$-norm. A numerical example illustrated the theoretical results.