Numerical simulation method is proposed for Brillouin optical correlation domain reflectometry in order to verify the characteristics in strain or temperature distributed sensing. The simulation results for the Brilluion gain spectrum obtained by the system agree well with previous theoretical analysis and experimental results, which make us understand more this reflectometry for further performance improvements.