<p>In this paper, we propose a method for quantitatively analysing spatial distribution faults based on an improved fractal dimension and a 3-D geological model. The method utilizes the Delaunay algorithm to improve the fractal dimension calculation. The 3-D geological model is used to calculate the curvature of the fault plane and the mean curvature of the fault plane is used to correct the number of parameters in the calculation of 3-D fractal dimension. The calculated indexes can reflect the spatial distribution characteristics of the faults in the study area as well as their own structural characteristics. In addition, a comparison analysis was carried out using existing geological reports from the study area to verify the validity of the method. Finally, we calculated the correlation coefficients between the calculated results and the distribution of water inrush points by linear regression analysis. The regression coefficients obtained were &gt; 0.7, which proves that the spatial distribution of faults was strongly correlated with the water inrush conditions and further verifies the validity of the method proposed in this paper.</p>