Based on N×N-type sampling law of subharmonic method ,incremental sampling is introduced ,and an incremental N×N-type sampling law of subharmonic method is proposed for numerical simulation of ocean turbulence phase screen. This method uses the refractive index fluctuation spectrum proposed by Nikishov to simulate the phase of ocean turbulence ,and verifies the accuracy of the phase screen generated by the phase structure function and compen- sation effect. By changing the number of incremental regions ,sampling levels and subharmonic orders ,the influence of phase screen parameters on the accuracy of phase screen is compared and analyzed. In addition ,the proposed meth- od is compared with the subharmonic method of N× N-type sampling law ,the modified subharmonic method and the improved subharmonic method. The results show that when the subharmonic order of the incremental N×N-type sam- pling law of subharmonic method is 11 ,the compensation effect can reach 0. 999 5 ,and the simulation speed is faster with the same accuracy.