Abstract The research of viscoelastic media is currently a hot topic in the interpretation and processing of seismic data. To accurately simulate the propagation of seismic waves in viscoelastic media, the fractional viscoelastic equation has emerged as an indispensable method. However, solving this equation numerically has proven to be challenging due to the complexity introduced by its fractional Laplacian operators. Recently, deep learning, especially Fourier neural operators (FNO), has shown excellent performance in learning to fast solve partial differential equations. Traditional FNO methods may face crosstalk problem and make it difficult to achieve satisfactory accuracy, when solving the multi-component fractional order viscoelastic equation. To solve this problem, we introduce a novel approach based on U-net Fourier neural operator (U-FNO). As an enhanced learning method to the traditional FNO-based method, the U-FNO-based method integrates a U-Fourier layer following the standard Fourier layer as a form of regularization, thereby achieving superior prediction accuracy for multi-component equations. Specifically, both the Fourier layers and U-Fourier layers in U-FNO are trained with the solutions of the equation from previous time steps as inputs. This training process enables the U-FNO to efficiently produce more accurate solutions for subsequent wavefield. Numerical simulations reveal that the U-FNO-based method efficiently learns to solve the fractional viscoelastic wave equation independent of fractional Laplacian operators. Additionally, U-FNO-based method offers superior prediction accuracy in comparison with the traditional FNO-based method.