The long-distance multi-stage pressurized pump station water delivery system involves numerous valve closure parameters, complicating the rapid identification of an optimal valve closure scheme that satisfies multiple transient flow oscillation protection requirements. A hydraulic transient model was established based on transient flow calculation theory to address this challenge. Decision biases were identified using the Analytic Hierarchy Process and the Entropy Weight Method. A multi-objective optimization model, incorporating Support Vector Regression (SVR) and the Beluga Whale Optimization (BWO) algorithm, iteratively searches for optimal schemes under different biases. The results indicate that Support Vector Regression exhibits optimal performance, while Beluga Whale Optimization demonstrates excellent performance. The optimal schemes obtained from the multi-objective optimization model meet the transient flow protection requirements of the water delivery system. The study demonstrates that this model effectively solves the multi-objective optimization problem for water hammer protection in multi-stage pressurized pump station water delivery systems.