Summary: Migration plays a crucial role in epidemic spreading, and its dynamic can be studied by metapopulation model. Instead of the uniform mixing hypothesis, we adopt networked metapopulation to build the model of the epidemic spreading and the individuals’ migration. In these populations, individuals are connected by contact network and populations are coupled by individuals migration. With the network mean-field and the gravity law of migration, we establish the N-seat intertwined SIR model and obtain its basic reproduction number \(\mathscr{R}_0\). Meanwhile, we devise a non-Markov node-search algorithm for model statistical simulations. Through the static network migration ansatz and \(\mathscr{R}_0\) formula, we discover that migration will not directly increase the epidemic replication capacity. But when \(\mathscr{R}_0>1\), the migration will make the susceptive population evolve from metastable state (disease-free equilibrium) to stable state (endemic equilibrium), and then increase the influence area of epidemic. Re-evoluting the epidemic outbreak in Wuhan, top 94 cities empirical data validate the above mechanism. In addition, we estimate that the positive anti-epidemic measures taken by the Chinese government may have reduced 4 million cases at least during the first wave of COVID-19, which means those measures, such as the epidemiological investigation, nucleic acid detection in medium-high risk areas and isolation of confirmed cases, also play a significant role in preventing epidemic spreading after travel restriction between cities.