Vector boson scattering at the Large Hadron Collider (LHC) is sensitive to anomalous quartic gauge coup-lings (aQGCs).In this study,we investigate the aQGC contribution to Wyjj production at the LHC with √s =13 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 polar-ization 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 cur-rent luminosity.The results indicate that the process pp → Wγjj is powerful for searching for the OM2,3,4,5 and OT5,6,7 operators.