Modeling of the shock cone formed around a stationary, hairy Horndeski black hole with Bondi-Hoyle-Lyttleton (BHL) accretion has been conducted. We model the dynamical changes of the shock cone resulting from the interaction of matter with the Horndeski black hole, where the scalar field and spacetime have a strong interaction. The effects of the scalar hair, the black hole rotation parameter, and the impacts of the asymptotic speed have been examined, revealing the influence of these parameters on the shock cone and the trapped QPO modes within the cone. Numerical calculations have shown that the hair parameter significantly affects the formation of the shock cone. As the absolute value of the hair parameter increases, the matter in the region of the shock cone is observed to move away from the black hole horizon. The rate of matter expulsion increases as h/M changes. After h/M < -0.6, a visible change in the physical structure of the shock cone occurs, ultimately leading to the complete removal out of the shock cone. On the other hand, it has been revealed that the asymptotic speed significantly affects the formation of the shock cone. As h/M increases in the negative direction and the asymptotic speed increases, the stagnation point moves closer to the black hole horizon. When the value of the hair parameter changes, the rest-mass density of the matter inside the cone decreases, whereas the opposite is observed with the asymptotic speed. Additionally, the formed shock cone has excited QPO modes. The deformation of the cone due to the hair parameter has led to a change or complete disappearance of the QPOs. Meanwhile, at asymptotic speeds of V ∞/c < 0.4, all fundamental frequency modes are formed, while at V ∞/c = 0.4, only the azimuthal mode is excited, and 1:2:3:4:… resonance conditions occur. No QPOs have formed at V ∞/c = 0.6. The results obtained from numerical calculations have been compared with theoretical studies for M87*, and it has been observed that the possible values of h/M found in the numerical simulations are consistent with the theory. Additionally, the results have been compared with those for the GRS 1915+105 black hole, and the hair parameters corresponding to the observed frequencies have been determined.