The spatial power spectrum function of oceanic turbulence is used to calculated the variations of turbulence perturbation item with turbulence parameters in extended Huygens-Fresnel integral by numerical methods. Then, the change of propagation quality factor of partically coherence Hermite-Gaussian (PCHG) beam traveling in oceanic turbulence is calculated by Wigner distribution function and the second moment definition of the beam. The affection of beam parameters on M2 factor is also analyzed. At last, the matrix of M2 factor is used to calculate the track of M2 factor of nonsymmetrical two-dimensional PCHG beam. It is found that the track plots of two dimensional PCHG beams′ spots on the receive plane are different from each other because of different models on x direction and y direction, but the difference decreases with the increase of propagation distance, and finally tends to same circular track.