A Reduced Basis Method for a PDE-constrained optimization formulation in Discrete Fracture Network flow simulations