Spin correlation functions (up to the 3-site one) of disordered Ising model with the nearest neighbour interaction are calculated and investigated within a two-site cluster approximation for both quenched and annealed cases. The approach yields the exact results for the one-dimensional system. The long-range interaction is taken into account in the mean field approximation.