A survey of the literature is presented regarding the simulation of port wine stain (PWS) skin color. Knowledge of PWS features, such as the depths and diameters of affected vessels, is essential for informing laser treatment. These may be determined through the inverse application of a skin model. The techniques which have been applied to achieve this are analyzed in detail. Radiative transfer (RT) is found to be the preferred method of simulation. By far the most common approximations to RT are the diffusion approximations, which have been applied successfully in the past and Monte Carlo techniques, which are now the methods of choice. As the requirements for improvement of laser treatment on an individual basis continues, the needs for further work towards accurate estimations of individual optical coefficients and robust, flexible simulation techniques are identified.