Quantification of roughness plays an important role in modeling strength deformability and fluid flow behaviors of rock joints. A procedure was suggested to simulate joint roughness, and characteristics of the roughness was investigated in this study. Stationary fractional Brownian profiles with known input values of the fractal parameter and other profile properties were generated based on random midpoint displacement method. Also, a procedure to simulate three dimensional roughness surface was suggested using the random midpoint displacement method. Selected statistical roughness parameters were calculated for the generated self-affine profiles to investigate the attribute of roughness. Obtained results show that statistical parameters applied in this study were able to consider correlation structure and amplitude of the profiles. However, effect of data density should be tackled to use statistical parameters for roughness quantification.