In order to more accurately and effectively understand the intermuscular coupling of different temporal and spatial levels from the perspective of complex networks, a new multi-scale intermuscular coupling network analysis method was proposed in this paper. The multivariate variational modal decomposition (MVMD) and Copula mutual information (Copula MI) were combined to construct an intermuscular coupling network model based on MVMD-Copula MI, and the characteristics of intermuscular coupling of multiple muscles of upper limbs in different time-frequency scales during reaching exercise in healthy subjects were analyzed by using the network parameters such as node strength and clustering coefficient. The experimental results showed that there are obvious differences in the characteristics of intermuscular coupling in the six time-frequency scales. Specifically, the triceps brachii (TB) had relatively high coupling strength with the middle deltoid (MD) and posterior deltoid (PD), and the intermuscular function was closely connected. However, the biceps brachii (BB) was independent of other muscles. The intermuscular coupling network had scale differences. MVMD-Copula MI can quantitatively describe the relationship of multi-scale intermuscular coupling strength, which has good application prospects.