According to the Fraunhofer diffraction theory, a theoretical model of far-field intensity distribution of two-dimensional M×N laser coherent array based on phase modulation is established. Numerical simulation on far-field intensity distribution of a 5×5 laser coherent array aiming at the application practice is made and the influence of different phases on intensity distribution is also analyzed. The results show that the positions of the main maximum or minimum and the value of the intensity distribution change with the phase modulation, and the intensity distribution varies with different array structures. More symmetric main maxima and secondary maxima appear in the simulation results when the phase difference loaded in each column and each row is π. All of these provide a helpful reference for the control of far-field intensity distribution by phase modulation.