Abstract Vector boson scattering at the Large Hadron Collider (LHC) is sensitive to anomalous quartic gauge couplings (aQGCs). In this study, we investigate the aQGC contribution to Wγjj production at the LHC with TeV in the context of an effective field theory (EFT). The unitarity bound is applied as a cut on the energy scale of this production process, which is found to have significant suppressive effects on signals. To enhance the statistical significance, we analyze the kinematic and polarization features of the aQGC signals in detail. We find that the polarization effects induced by aQGCs are unique and can discriminate the signals from the SM backgrounds well. With the proposed event selection strategy, we obtain the constraints on the coefficients of dimension-8 operators with current luminosity. The results indicate that the process is powerful for searching for the and operators.