The method of invariant grid (MIG) is an iterative procedure for model reduction in chemical kinetics which is based on the notion of Slow Invariant Manifold (SIM) [A.N. Gorban, I.V. Karlin, Method of invariant manifold for chemical kinetics, Chem. Eng. Sci. 58 (2003) 4751-4768; E. Chiavazzo, A.N. Gorban, I.V. Karlin, Comparison of invariant manifolds for model reduction in chemical kinetics, Commun. ]. Important role, in that method, is played by the initial grid which, once refined, gives a description of the invariant manifold: the invariant grid. A convenient way to get a first approximation of the SIM is given by the spectral quasi-equilibrium manifold (SQEM) [A.N. Gorban, I.V. Karlin, Method of invariant manifold for chemical kinetics, Chem. Eng. Sci. 58 (2003) 4751-4768; E. Chiavazzo, A.N. Gorban, I.V. Karlin, Comparison of invariant manifolds for model reduction in chemical kinetics, Commun. Comput. Phys. 2(5) (2007) 964-992].In the present paper, a flexible numerical method to construct the discrete analog of a quasi-equilibrium manifold, in any dimension, is presented. That object is named quasi-equilibrium grid (QEG), while the procedure quasi-equilibrium grid algorithm (QEGA). Extensions of the QEM notion are also suggested. The QEGA is a numerical tool which can be used to find a grid-based approximation for the locus of minima of a convex function under some linear constraints. The method is validated by construction of one and two-dimensional grids for a model of hydrogen oxidation reaction.