We investigate the nonperturbative renormalization of composite operators i
n lattice QCD restricting ourselves to operators that are bilinear in the q
uark fields. These include operators which are relevant to the calculation
of moments of hadronic structure functions. The computations are based on M
onte Carlo simulations using quenched Wilson fermions.