Surface operators in supersymmetric gauge theories and their operator product expansion