Abstract: Considering the influence of in-situ stresses on the mechanical properties of filled jointed rock mass and the attenuation of stress wave amplitude and time delay caused by the viscoelasticity of rock mass, and using the time-domain recursive analysis method, the P-wave propagation equations in viscoelastic filled jointed rock mass under in-situ stresses are established by introducing the quality factors of stress wave. In the analysis, an improved three-element model is adopted for filled joints, and the deformation of the joint contact interfaces is assumed to satisfy the nonlinear hyperbolic deformation model in normal direction and the linear elastic model in tangential direction. The results show that with the increase of in-situ stress, the transmission coefficients of stress wave increase and the energy dissipation of stress wave decreases. The change of the filled joint angle will change the impinge angle and stress distribution of the joint, which has obvious effect on the wave propagation characteristics. The increase of joint filled thickness increases the energy dissipation and decreases the transmission coefficients. The energy dissipation of stress wave decreases obviously with the increases of wave frequency, which reflects the obstruction and filtering effect of high frequency waves. With the increases of initial stiffness of joint contact interface, the transmission coefficients of stress wave increase gradually and the reflection coefficients decrease gradually, while the energy dissipation does not change significantly.