Virtually every econometrics textbook describes the OLS estimation procedure by referring to the orthogonality between the OLS fitted residuals and vectors in the column space of the explanatory variables.
