We numerically investigate the formation and ordering of staggered oblate particle pairs in three-dimensional straight ducts with a square cross section. The lattice Boltzmann method is employed to simulate rigid particle pairs in a Newtonian liquid. The effects of initial axial spacing, Reynolds number, blockage ratio, and particle aspect ratio on the formation process, migration behavior, and interparticle spacing are explored in detail. Current results indicate that the process from initial to final steady state can be divided into two stages. The first stage is rapid migration from initial positions toward equilibrium positions under shear-induced lift force and wall-induced repulsive force. The second stage is the slow self-assembly of stable particle pairs in the axial direction due to the interparticle interaction. Interestingly, initial axial spacing significantly affects the formation process of particle pairs but does not affect the final steady state. It is found that the equilibrium positions of staggered particle pairs move slightly toward the duct walls, and the axial spacing increases with increasing Reynolds number or particle aspect ratio, or decreasing blockage ratio. For a staggered particle pair, the second particle will occupy the eddy center induced by the first focusing particle. Based on the existing data, a correlation is put forward to predict the axial interparticle spacing of staggered oblate particle pairs in duct flows. The present results may give insights into manipulating and comprehending non-spherical particle dynamics in microfluidic applications.