## 1 Introduction
[31] placed vibration monitoring points on mine slopes and established a vibration attenuation model based on the relationship between peak vibration velocity and scaled distance, thereby determining the delay time that meets slope stability requirements. [PERSON] et al. [32] precisely measured and collected data from 112 blasts to protect nearby residences and other structures in a limestone mine, constructing an artificial neural network (ANN) model to predict PPV and identify the optimal delay time. [PERSON] et al. [33] proposed a blast vibration spectrum control scheme suitable for high-bench slope excavation in hydropower projects, achieving vibration prediction and delay optimization based on single-hole vibration data monitored at different height positions. In summary, the aforementioned methods can achieve effective blast vibration control at specific locations. In essence, controlling multi-hole blast vibrations by adjusting delay times relies on the periodicity of vibration waves, utilizing the principle of superimposing wave peaks and troughs to achieve cancellation [34; 35]. However, two critical problems remain unresolved in bench blasting vibration control based on this principle. Firstly, variations in initiation times cause changes in the arrival times of vibration waves from each blast hole to the same monitoring point, and the arrival times of these waves at different monitoring points also vary. These factors significantly impact the superposition and cancellation effects of wave peaks and troughs, thereby affecting vibration reduction. Secondly, the delay times obtained through the misalignment superposition technique for vibration reduction can only achieve vibration control at a single monitoring point. Under this premise, as the initiation times are predetermined, monitoring points at other locations around the blasting area cannot further reduce vibrations using the misalignment superposition technique. Consequently, vibration monitoring at a single point is insufficient to meet the requirements of comprehensive vibration control, potentially leading to misinterpretation of the vibration control situation by technicians and resulting in inadequate corrective measures. On the basis of meeting the vibration control requirements of a single monitoring point, there is an urgent need for a vibration control method that can take into account the positions of multiple monitoring points. To solve the above problems, this paper is conducted under the premise of a predefined initiation sequence, which can be determined using vibration wave misalignment superposition techniques or other established methods. Building upon this foundation, localized adjustments are made to vibration waves arriving at monitoring points within the same time window. By systematically shifting the initiation times of subsequent blast holes, this approach achieves vibration reduction while effectively preserving the fragmentation performance of the original delay scheme. Specifically, a computational method is established to determine the arrival times of vibration waves from each hole to the monitoring points based on the spatiotemporal relationship between blast holes and monitoring points. This enables the identification of the number of vibration waves arriving at the monitoring points within the same time window and the corresponding hole indices. The vibration waves in the same time window with excessive concurrent arrivals (clustering arrival period) are adjusted by modifying the initiation delays of the relevant holes. Subsequently, the reconstructed single-hole vibration curves are superimposed according to the revised arrival sequence to obtain the predicted vibration curve for multi-hole blasting. Furthermore, to investigate vibration control across different spatial locations, a dense array of monitoring points is deployed around the blasting area, and the arrival times of vibration waves from all blast holes to each monitoring point are exhaustively calculated.