One fast simulation method using Markov chains was introduced to simulate angular, energy, and temporal characteristics of pulsed laser beam propagation underwater. Angular dispersion of photons with a different number of collisions was calculated based on scattering function and the state transition matrix of Markov chains. Temporal distribution and energy on the receiving plane were obtained, respectively, by use of a novel successive layering model and receiving ratio. The validity of this method was verified by comparing it with the Monte Carlo ray tracing (MCRT) method. The simulation results were close to those obtained by MCRT but were less time consuming and had smoother curves.