Abstract: A comprehensive dataset from 594 fracturing wells throughout the Duvernay Formation near Fox Creek, Alberta, is collected to quantify the influences of geological, geomechanical, and operational features on the distribution and magnitude of hydraulic fracturing-induced seismicity. An integrated machine learning-based investigation is conducted to systematically evaluate multiple factors that contribute to induced seismicity. Feature importance indicates that a distance to fault, a distance to basement, minimum principal stress, cumulative fluid injection, initial formation pressure, and the number of fracturing stages are among significant model predictors. Our seismicity prediction map matches the observed spatial seismicity, and the prediction model successfully guides the fracturing job size of a new well to reduce seismicity risks. This study can apply to mitigating potential seismicity risks in other seismicity-frequent regions.