Convection initiation (CI) associated with a cold pool gust front (GF) of a mesoscale convective system (MCS) occurred in the Bohai Bay region, North China, is investigated using a real‐data, high‐resolution Weather Research and Forecasting model (WRF) simulation. Observations revealed that many CIs occurred successively near the GF, most of which developed significantly in both size and intensity and eventually merged with the MCS. The WRF simulation captured the general features of the CIs, GF, and MCS. Lagrangian vertical momentum budgets were conducted along the backward trajectories of air parcels within a convective cell initiated near the GF. The total vertical acceleration was decomposed into dynamic and buoyant components. The results showed that the dynamic acceleration played a decisive role for the CI, which can be further decomposed into five terms. The horizontal curvature term was rather small, with the height variation of air density almost zero. Although both vertical and horizontal extension terms played a positive role during the CI process, the vertical twisting term had the greatest contribution to the CI. This reflects the vertical tilting (to the forward of the GF) of an elevated updraft which was generated by the convergence of a forward‐tilting updraft at the GF leading edge and another airflow coming from ahead of the GF.