In recent years, the lattice Boltzmann method has achieved tremendous success as a novel numerical method in scientific computing. Compared with traditional computational fluid dynamics methods, the lattice Boltzmann method possesses unique advantages in parallel efficiency and program implementation. In this paper, the lattice Boltzmann method is employed to numerically simulate reactive flow in porous media. Specifically, the lattice Boltzmann method for solving the flow field is derived from incompressible macroscopic equations and formulated for porous media flow, while the lattice Boltzmann method for solving the species transport equation is based on a hybrid regularization approach and exhibits high numerical stability. In this study, the correctness of the program implementation is first verified by comparing the results with those obtained from other computational fluid dynamics methods reported in the literature and from commercial software. Subsequently, a porous media catalytic reactor model is simulated, and the computational results demonstrate that the lattice Boltzmann method is suitable for simulating such scenarios.