The level of detail on modern geological models requires higher resolution grids that may render the simulation of multiphase flow in porous media intractable. Moreover, these models may comprise highly heterogeneous media with phenomena taking place in different scales. Scale transferring techniques allow for the solution of such problems in a lower resolution scale at reduced computational cost. Among these techniques, the Multiscale Finite Volume (MsFV) method constructs a set of numerical operators in order to map quantities from the fine-scale mesh to a coarser one and vice-versa while maintaining flux conservation on both scales. However, the MsFV formulation, as originally stated, is only consistent on k-orthogonal grids since it uses a linear Two-point Flux Approximation (TPFA) method and may struggle to generate consistent primal-dual coarse grids pairs on unstructured grids. The Multiscale Restriction Smoothed-Basis method (MsRSB) improves on the MsFV by introducing a new iterative procedure to find the multiscale operators and the concept of support regions which reduces the method's complexity when applied to unstructured fine and coarse grids. The original version was only consistent on k-orthogonal fine grids due to the TPFA discretization, but filtering methods have been developed to also enable consistent multipoint schemes on the fine scale. Meanwhile, the Multiscale Control Volume method (MsCV) replaces the TPFA by the Multipoint Flux Approximation with a Diamond stencil (MPFA-D) scheme on the fine-scale while further enhancing the generation of the geometric entities to allow unstructured grids on the fine and coarse scales for two-dimensional simulations. In this work we propose an extension to three-dimensional geometries of both the MsCV and the algorithm to obtain the multiscale geometric entities based on the concept of a background grid, a coarser grid used as a proxy for the primal coarse grid. We modify the original MPFA-D method to use the very robust Global Least Squares (GLS) interpolation technique to obtain the required auxiliary nodal unknowns. We also introduce an enhanced version of the 3-D MsCV with the incorporation of the enhanced MsRSB (E-MsRSB) to enforce M-matrix properties and improve convergence. Finally, we employ the MsCV operators in a two-stage smoother to show how it can be used as a good iterative procedure to recover the fine-scale solution. We show that the 3-D MsCV method produces good results employing unstructured grids on both scales to handle the simulation of the single-phase flow in anisotropic, and heterogeneous porous media and that the smoothing procedure is able to accurately retrieve the fine-scale solution within a few iterations.