Fracture identification and evaluation requires data from various resources, such as image logs, core samples, seismic data, and conventional well logs for a meaningful interpretation. However, several wells have some missing data; for instance, expensive cost run for image logs, cost concern for core samples, and occasionally unsuccessful core retrieving process. Thus, a majority of the current research is focused on predicting fracture based on conventional well log data. Interpreting fractures information is very important especially to develop reservoir model and to plan for drilling and field development. This study employed statistical methods such as multiple linear regression (MLR), principal component analysis (PCA), and gene expression programming (GEP) to predict fracture density from conventional well log data. This study explored three wells from a basement metamorphic rock with ten conventional logs of gamma rays, thorium, potassium, uranium, deep resistivity, flushed zone resistivity, bulk density, neutron porosity, sonic porosity, and photoelectric effect. Four different methods were used to predict the fracture density, and the results show that predicting fracture density is possible using MLR, PCA, and GEP. However, GEP predicted the best fracture density with R2 > 0.86 for all investigated wells, although it had limited use in predicting fracture density. All methods used highlighted that flushed zone resistivity and uranium content are the two most significant well log parameters to predict fracture density. GEP was efficient for use in metamorphic rocks as it works well for conventional well log data as the data is nonlinear, and GEP uses nonlinear algorithms. © The Author(s) 2024.