Many recent studies have revealed that spatial interactions of functional brain networks derived from fMRI data can well model functional connectomes of the human brain. However, it has been rarely explored what the energy consumption characteristics are for such spatial interactions of macro-scale functional networks, which remains crucial for the understanding of brain organization, behavior, and dynamics. To explore this unanswered question, this article presents a novel framework for quantitative assessment of energy consumptions of macroscale functional brain network's spatial interactions via two main effective computational methodologies. First, we designed a novel scheme combining dictionary learning and hierarchical clustering to derive macro-scale consistent brain network templates that can be used to define a common reference space for brain network interactions and energy assessments. Second, the control energy consumption for driving the brain networks during their spatial interactions is computed from the viewpoint of the linear network control theory. Especially, the energetically favorable brain networks were identified and their energy characteristics were comprehensively analyzed. Experimental results on the Human Connectome Project (HCP) task-based fMRI (tfMRI) data showed that the proposed methods can reveal meaningful, diverse energy consumption patterns of macro-scale network interactions. In particular, those networks present remarkable differences in energy consumption. The energetically least favorable brain networks are stable and consistent across HCP tasks such as motor, language, social, and working memory tasks. In general, our framework provides a new perspective to characterize human brain functional connectomes by quantitative assessment for the energy consumption of spatial interactions of macro-scale brain networks.