  To improve accuracy in calculating QCD effects, we propose a method for
renormalon subtraction in the context of the operator-product expansion. The
method enables subtracting renormalons of various powers in $\Lambda_{\rm QCD}$
efficiently and simultaneously from single-scale observables. We apply it to
different observables and examine consistency with theoretical expectations.
