Root architecture systems (RAS) reflect the spatial structure of roots in soil. To clarify the structure and distribution of rice roots and investigate the coupling between roots and soil, wetland rice was selected as the experimental object, and a three-dimensional (3D) growth model of rice root environment-roots (ERoots) based on the parameter Lindenmayer system (L-system) was proposed. ERoots combines a root morphological structure model with a growth model and defines L-system grammar iteration rules with the unit time and unit step length as parameters. At the same time, the basic growth parameters of rice roots were obtained via destructive detection, and 3D growth visualisation of roots was realised via MATLAB. In the soil coupling process, a soil nutrient simulation map was constructed based on the spatial soil characteristics per unit volume, and an adjustment strategy for roots reaching the growth boundary was designed. The flexibility of the model coupled with soil was reflected in the tropisms of root growth, growth rate and root branching strategy. Finally, combined with soil spatial characteristic simulation, geometric growth boundary and 3D root growth model, the ability of 3D growth visualisation of rice roots was verified under three soil conditions: (1) unconfined root growth, (2) confined spatial root growth, and (3) root growth with tropisms. The results indicated that the ERoots root model basically realised coupling with soil and achieved a satisfactory simulation effect in regard to the rice morphological structure. This study provides a reference for 3D growth modelling and visualisation of other crop roots.