We develop a novel large-scale Network Autoregressive model with balance Constraint (NAC) to predict hour-ahead gas flows in the gas transmission network, where the total in- and out-flows of the network are balanced over time. By integrating recent advances in optimization and statistical...