Adaptive finite element-discrete element method for numerical analysis of the multistage hydrofracturing of horizontal wells in tight reservoirs considering pre-existing fractures, hydromechanical coupling, and leak-off effects