Monte Carlo simulation of oxygen disordering in the YBa2 (Cu1-xMx)3Oy (M = Al, Fe, Co, Ga) system