Determination of optical parameters of turbid media is quite useful in the photodynamic therapy and optical noninvasive diagnostics. A kind of real coded genetic algorithm incorporated with inverse Monte Carlo method and graphics processing unit based acceleration technology is proposed, which can determine optical parameters from the spatially resolved diffuse reflectance of turbid media by Monte Carlo simulation. Fitness function of accumulated square differences, random tournament selection operator, uniform random crossover operator with extended radius, uniform mutation operator, champion mutation operator are designed to guaranty the algorithm to converge with good population diversity. In the range 0≤μa≤100 cm-1 and 0≤μs≤1000 cm-1, the average relative errors are 0.25% and 0.58%, and the root mean-square errors (RMSEs) are 0.32 cm-1 and 1.68 cm-1 for the absorption coefficient and for the scattering coefficient, respectively, which means that this algorithm is feasible and accurate enough for determination of optical parameters of turbid media.