In logarithmic conformal field theory, primary fields come together with logarithmic partner fields on which the stress-energy tensor acts non-diagonally. Exploiting this fact and global conformal invariance of two-and three-point functions, operator product expansions of logarithmic operators in arbitrary rank logarithmic conformal field theory are derived.