We investigate the stochastic homogenization of a class of turbulent diffusions generated by non-local symmetric Lévy operators with divergence-free drift fields in ergodic random environments, where neither the drift fields nor their associated stream functions are assumed to be bounded. A pivotal step in our proof is the establishment of $W_{loc}^{1,q}$ estimates with $q\in (1,2)$ for the corresponding correctors, under mild prior regularity conditions imposed on the Lévy measure and the stream function.