The aim of this work was to find an appropriate reaction-scheme to model both pure H2S and CH4 + H2S mixture oxidation systems. The H2S decomposition rate expression constants were adjusted, in which the mean absolute percentage error (MAPE) of this model was 5.91%. Then, three reaction schemes for pure H2S oxidation and two reaction schemes for the CH4 + H2S co-oxidation system were selected, and their kinetic constant were adjusted. The selected reaction schemes with optimized kinetic parameters were used to model H2S and CH4 partial oxidation, in which the corresponding MAPEs varied between 7.90 and 10.32% Finally, a reaction scheme was chosen for complete description of the pure H2S oxidation and H2S + CH4 co-oxidation systems. The proposed kinetic reaction scheme was applied in the kinetic modeling of an industrial Claus reaction furnace (MAPEs of 5.59 and 15.02% for temperature and composition, respectively) and provided a satisfactory representation of the experimental data. Performances of the proposed reaction scheme were validated with comparison to the detailed reaction mechanism.