We study the quasi-normal modes for scalar, electromagnetic, and gravitational axial perturbations in the Hayward regular black hole surrounded by quintessence (HBH-ωq). Using the third-order WKB approximation we can determine the dependence of the quasi-normal modes on the parameters of the regular black hole and the parameters on the test fields. We also determine the greybody factor, giving transmission and reflection coefficients of the scattered wave through the effective potentials in the WKB approximation using numerical analysis.