Numerical integration based on trivariate $C^2$ quartic spline quasi-interpolants