This article explores the dynamic behavior of a fractional-order Schnakenberg chemical reaction model. Specifically, we conduct an analysis of codimension-two bifurcation associated with 1:2, 1:3, and 1:4 resonances. To achieve these results, we utilize the normal form method and bifurcation theory. The findings are illustrated through detailed numerical simulations, including visualizations like two-parameter bifurcation diagrams and maximum Lyapunov exponent plots. These simulations effectively explore the system's behavior under the influence of two varying parameters within a three-dimensional space. Additionally, the simulations vividly demonstrate the theoretical results and offer valuable insights into the underlying dynamics.