Abstract: In view of the phenomenon that the surface loess gullies have great influence on the mine pressure appearance of the working face, based on the engineering geological reality of working face 12309 in Wangjialing coal mine, the distribution pattern of the surface loess gullies in the working face is summarized by the method of field investigation, the authors construct the three-dimensional shape database of the loess gully surface in the working face by Get Data software, combine with Rhinoceros, Midas GTS and other auxiliary software, establishe the three-dimensional and full-scale FLAC^3D numerical calculation model of 12309 working face with loess gully landform, use the numerical calculation method to analyze the distribution law of the in-situ rock stress under the loess gully landform, put forward the loess gully landform shadow the response coefficient η, combined with the numerical calculation results, has fitted out the non-linear expression with the coal seam buried depth h as the independent variable η.With the advance bearing pressure of the working face as the early warning threshold, the influence area of the surface loess gully during the mining of the working face is divided into “the area with strong influence on the loess gully landform” “the area with obvious influence on the loess gully landform” “the area with weak influence on the loess gully landform”.The early warning division of 12309 working face is carried out under the influence of the surface loess gullies, which provides a theoretical basis for ensuring the safe and efficient mining of the working face.