Modelling Neuron-Glial Network Interactions at the Whole-Brain Scale for Human Neuroimaging Applications