The soliton propagation of the fractional-in-space nonlinear Schrodinger equation (NLSE) is much more complicated than that of the corresponding integer NLSE. The aim of this paper is to discover some novel fractal soliton propagation behaviors (FSPBs) of this fractional-in-space NLSE. Firstly, the exact solution is compared with the present numerical solution, and the validity and accuracy of the present numerical method are verified. Secondly, the effect of fractional derivatives on soliton propagation is explored through the present numerical simulation results. At the same time, the present method is extended to the three-dimensional fractional-order NLSE. Finally, some novel FSPBs of the fractional-in-space NLSE are given.