Ok, so your nmos is working in weak inversion (subthreshold), as your Vth will probably be around 400-500mV.
I guess that you want to compare your results with the formula Id=W/L I0' *exp[(vgs-vth)/nUT]. Vth is strongly dependent on the length and in case of very small widths also width dependent.
Plot also Vth as function of W and L. you'll see that the Vth is dependent on L. For large L the Vth remains constant, resulting in a linear relation between the length and current. For small W and L, Vth can change a lot and I think that this is the source of the effect you simulated. To know for sure that the non-constant vth it really is the dominant source, you have to check it yourself with the Weak-inversion formula, as i'm not a device physicist.