This article analyzes the numerical modeling of the stress-strain state of a gravity structure sand base during grouting works using the controlled compensation grouting technique. The modeling is performed using the finite element method and physical nonlinear soil models that take into account the Coulomb-Mohr strength condition. The modeling is done in the computational software JuliS as well as the universal software complex PC ZSoil. The goal of the research is to study the impacts of compensation grouting parameters on the elevation value of base slab. A model was created for the grout body as an expanding spheroid whose volume changes as cement grout is fed in the soil. The authors developed a grouting sequence for grouting solutions based on the location of grout body in the grouting pillar. It is more efficient to feed the grout starting from the lower grout body. The grouting solution feed mode is selected depending on the soil strength condition and accounting for the prevention of hydraulic fracture. As a result of the research, the initial subsidence of the soil on the daytime surface above the injection zone was obtained, which requires experimental confirmation and is a matter for discussion.