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