In this paper, we propose an integral-averaged interpolation operator I τ in a bounded domain Ω ⊂ R n by using Q 1 -element. The interpolation coefficient is defined by the average integral value of the interpolation function u on the interval formed by the midpoints of the neighboring elements. The operator I τ reduces the regularity requirement for the function u while maintaining standard convergence. Moreover, it possesses an important property of I τ u 0,Ω ≤ u 0,Ω . We conduct stability analysis and error estimation for the operator I τ . Finally, we present several numerical examples to test the efficiency and high accuracy of the operator.