We consider the open XXZ spin chain with diagonal boundary terms in the framework of algebraic Bethe Ansatz. Using the explicit computation of the scalar products of Bethe states and a revisited version of the bulk inverse problem, we calculate the correlation functions. In the limit of half-infinite chain, they are obtained as multiple integrals, like in the case of periodic boundary conditions.