Control Variates for Stochastic Simulation of Chemical Reaction Networks