We study the boundary effects in invasion percolation with and without trapping. We find that the presence of boundaries introduces a new set of surface critical exponents, as in the case of standard percolation. Numerical simulations show a fractal dimension, for the region of the percolating cluster near the boundary, remarkably different from the bulk one. We find a logarithmic cross-over from surface to bulk fractal properties, as one would expect from the finite-size theory of critical systems. The distribution of the quenched variables on the growing interface near the boundary self-organises into an asymptotic shape characterized by a discontinuity at a value $x_c=0.5$, which coincides with the bulk critical threshold. The exponent $tau^{sur}$ of the boundary avalanche distribution for IP without trapping is $tau^{sur}=1.56pm0.05$; this value is very near to the bulk one. Then we conclude that only the geometrical properties (fractal dimension) of the model are affected by the presence of a boundary, while other statistical and dynamical properties are unchanged. Furthermore, we are able to present a theoretical computation of the relevant critical exponents near the boundary. This analysis combines two recently introduced theoretical tools, the Fixed Scale Transformation (FST) and the Run Time Statistics (RTS), which are particularly suited for the study of irreversible self-organised growth models with quenched disorder. Our theoretical results are in rather good agreement with numerical data.