Monte Carlo techniques have been widely applied in polarized light simulation. Based on different preconditions, there are two main types of sampling strategies for scattering direction: one is the scalar sampling method; the others are polarized sampling approaches, including the one- and two-point rejection methods. The polarized simulation of oceanic lidar involves a variety of mediums, and an efficient scattering sampling method is the basis for the coupling simulation of the atmosphere and ocean. To determine the optimal scattering sampling method for oceanic lidar simulation, we developed a polarized Monte Carlo model and simulated Mie scattering, Rayleigh scattering, and Petzold average-particle scattering experiments. This simulation model has been validated by comparison with Ramella-Roman's program [Opt. Express13, 4420 (2005)OPEXFF1094-408710.1364/OPEX.13.004420], with differences in reflectance and transmittance Stokes less than 1% in Mie scattering. The simulation results show these scattering sampling methods differ in runtime, scattering angle distributions, and reflectance and transmittance Stokes. Considering the current simulation accuracy of oceanic lidar, the differences in reflectance and transmittance Stokes are acceptable; thus, the runtime becomes the main evaluation factor. The one-point rejection method and scalar sampling method are preferable for the oceanic lidar polarized simulation. Under complex atmosphere-ocean coupling systems, scalar sampling methods may be a better choice since the calculation process of the sampling is independent of the incident Stokes vector.