Solar-induced chlorophyll fluorescence (SIF) is a subtle but informative probe of plant photosynthesis. Quantifying the three-dimensional (3D) distribution of SIF benefits a better understanding of photosynthesis variations over heterogeneous canopies. Although radiative transfer models (RTMs) provide a solid theoretical basis for simulating the 3D SIF distribution, most RTMs use virtual scenes with complex reconstruction processes. This study aims to develop a 3D SIF simulation model (FluorLiDAR) directly driven by terrestrial light detection and ranging (LiDAR) data using leaf and canopy RTMs, including a 3D PAR (photosynthetically active radiation) simulation model, the Fluspect model, the atmosphere radiative transfer module in SCOPE, and the multiple scattering coefficients of sunlit and shaded leaves from the 4-scale model. The results show that (1) the simulated and measured SIF patterns were consistent, with R2 (RMSE) values of 0.73 (0.17 mW/nm/m2/sr) and 0.76 (0.12 mW/nm/m2/sr) for 1-min sampling and 10-min averages, respectively. Moreover, the R2 between FluorLiDAR and the DART simulation reached 0.94. The R2 of FluorLiDAR and DART were both higher than 1D mSCOPE using every 10-min sampling data. (2) Point density, denoted by average NPD (nearest point distance), influenced the performance of our model mainly when it was smaller than 0.1 m. Chlorophyll content had less influence on the model accuracy (R2 and rRMSE), and the bias between simulation and measurement decreased as chlorophyll content increased. (3) The simulated 3D SIF distribution pattern closely resembled PAR within the canopy. Besides, FluorLiDAR can simulate the hot spot effect like DART and mSCOPE though the effect was not as obvious as the other two models. This study highlights the potential of a LiDAR-driven SIF model for 3D SIF simulation over a heterogeneous canopy, which may benefit the understanding of the structural impacts on forest photosynthesis in real forest scenes.