Dynamics of vertical-cavity surfac-emitting injectio lasers (VCSEL) was simulated. The model of VCSEL is based on the set of self-consistent differential equations in partition derivatives and material relations for quantum well structure. An optical field in active resonator was numerically obtained by effective frequency method in cylindrical coordinates. The variables separation method was used for 3D wave equation resolving. An optical field outside laser mirror was calculated by Fresnel integrals. Two-mode self-pulsing regime at approximately 2 GHz repetition rate was observed at bias current value of about twice threshold. The obtained near- and far-field dynamical behavior of LP01 and LP02 lateral modes represented the processes of spatial energy exchange between the two modes as well as between the parts of particular mode. In the vicinity about of 20 microns near output VCSEL surface (diffraction zone) the sharp spatial oscillations of optical field in the center of beam were observed. It was shown that the above dynamical processes are caused by non-uniform carriers rearrangement in active layer.