In this study, the partial least squares regression (PLSR) models were developed using no pre-processing, traditional preprocessing, multi-preprocessing 5 range, multi-preprocessing 3 range, genetic algorithm (GA), and successive projection algorithm (SPA) to assess the higher heating value (HHV) and ultimate analysis of grounded biomass for energy usage employing near-infrared (NIR) spectroscopy. A novel approach was utilized based on the assumption that using multiple pretreatment methods across different sections in the entire NIR wavenumber range would enhance the performance of the model. The performance of the model obtained from 200 biomass samples for HHV and 120 samples for ultimate analysis was compared, and the best model was selected based on the coefficient of determination of validation set, root mean square error of prediction, and the ratio of prediction to deviation values. Based on model performance results, the proposed HHV model from GA-PLSR, and the N and O models from the mul-ti-preprocessing PLSR 5 range method could be used for most applications, including research, whereas the C and H models from GA-PLSR performance is fair and applicable only for rough screening. The overall findings highlight that the multi-preprocessing 5-range method, which was attempted as a novel approach in this study to develop the PLSR model, demonstrated better accuracy for HHV, C, N, and O, improving by 4.1839%, 8.1842%, 3.7587%, and 35.9404%, respec-tively. Therefore, it can be considered a reliable and non-destructive alternative method for rap-idly assessing biomass properties for energy usage and can also be used effectively in biomass trading. However, due to the smaller number of samples used in the model development, more samples are needed to update the model for a robust application.