We investigate the low-order modeling of collective dynamics in a can-annular combustor consisting of four ring-coupled turbulent lean-premixed combustors. Each combustor is treated as an individual thermoacoustic oscillator, and the entire combustion system is modeled using four Van der Pol oscillators ring-coupled with dissipative, time-delay, and reactive coupling terms. We show that this model, despite its simplicity, can reproduce many collective dynamics observed in experiments under various combinations of equivalence ratios and combustor lengths, such as 2-can anti-phase synchronization, alternating anti-phase synchronization, pairwise anti-phase synchronization, spinning azimuthal mode, and 4 steady thermoacoustic oscillators. The phase relationship in the majority of cases can be quantitatively modeled. Moreover, by incorporating a reactive coupling term, the model is able to reproduce the frequency shift observed experimentally. This study demonstrates the feasibility of using a simple low-order model to reproduce collective dynamics in complex turbulent combustion systems. This suggests that this model could be used (i) to facilitate the interpretation of experimental data within the synchronization framework, (ii) to identify potential parameter regimes leading to amplitude death, and (iii) to serve as a basis for modeling the collective dynamics observed in more complicated multi-combustors.