In this study, optimization of oil agglomeration parameters of Beyşehir-Bayavşar lignite using Box–Behnken design (BBD) was investigated. BBD model predicted responses, the coefficient of determination R2 for the percentage of ash% and combustible recovery% (CR) of the agglomerates was obtained as 0.97 and 0.90, respectively. These results showed a good correlation between the calculated and observed values. The estimated the percentage of ash content and CR% values were determined as 15.81% and 73.73%, respectively. The best results were obtained in test 22, the zeta potential value was − 11.61 mV, the contact angle was 107°, and the calorific value was 4702 kcal/kg at pH 3.
