This study took the Tianmu Mountain National Natural Preservation Area in western Zhejiang Province as an example 1) to quantify the species diversity with selected species indices including Shannon-Wiener index (H), Margalef index (K) and Evenness index (E) at different spatial scales, 2) to analyze the spatial distribution patterns of the species diversity by multifractal parameters such as the singularity index α, its fractal dimension f(α), the f(α)–α spectrum range (SR) and its symmetry (Dist) using the multifractal theory, and 3) to determine their relationships. Results of nonlinear regression analysis with power functions showed that increasing spatial scale resulted in increasing H but decreasing E and K, indicating that the scale dependence of species diversity existed. By using the multifractal method, it was indicated that species spatial distribution had multifractal features. Moreover, strong linear relationships of the diversity indices H, E and K with α_min and clear nonlinear associations of the diversity indices H, E and K with power functions for SR and Dist were found. Since interactions of the species diversity and the spatial characteristics are very complex, the above mentioned relationships need further validation along with precise explanations of any correlations among their ecological processes.