Abstract: The layered rock mass and its interlayer shearing fracture zone have complex influences on the magnitude and direction of local stress field in underground engineering areas. Obtaining the initial in-situ stress field under complex geological conditions is the prerequisite for analyzing the stability of surrounding rock of underground caverns in layered rock mass. Firstly, aiming at the influences of evolution process, topography and tectonic action of a complex valley, a lateral pressure coefficient method based on the stratum denudation simulation is proposed for the first-stage inversion and the initial ground stress field of big model is obtained. Secondly, based on the influences of layered rock mass on local stress field, a second-stage inverse model for layered rock mass is established. The equivalent tectonic loads are obtained by the stress field calculated by interpolation from the big model, and the second-stage inverse analysis method based on the equivalent tectonic loads is proposed. Combined with the first and the second inverse methods, the optimal inverse analysis method for the three-dimensional initial in-situ stress field of layered rock mass is formed. Finally, the in-situ stress field of underground cavern of Guiyang pumped storage hydroelectric plant is inverted on account of the measured in-situ stress data. The secondary inverse results indicate that the inverse values of the in-situ stress field can meet the requirements of point coincidence at the measuring points and field coincidence reflecting the evolution process of the valley and influences of the layered rock mass. The local stress field in underground engineering areas is significantly affected by the interlayer shearing fracture zone, which is mainly manifested in the following aspects: (1) The local stress value increases slightly near the interface and releases in the soft rock layer. (2) The direction of the local stress field deflects differently due to change of the occurrence of rock strata. Using the proposed method, the influences of anisotropic mechanical properties of layered rock mass on the magnitude, direction and disturbance range of local in-situ stress field are deeply analyzed.