Coalbed methane (CBM) pilot wells typically exhibit a short production period, necessitating evaluation of their estimated ultimate recovery (EUR) through numerical simulation. Utilizing limited geological data from the pilot areas to finish history matching and subsequent production forecasting presents substantial challenges. This paper introduces a comprehensive numerical simulation workflow for CBM pilot wells, encompassing the following steps. Initially, geological parameters are categorized into two groups based on their statistical distribution trends: trend parameters (i.e., gas content, permeability, Langmuir volume, and Langmuir pressure) and non-trend parameters (i.e., fracture porosity, gas–water relative permeability, and rock compressibility). The probability method is employed to ascertain the probable high and low limits for trend parameter distributions, while empirical or analogous methods are applied to define the boundaries for non-trend parameters. Subsequently, the parameter sensitivity analysis is conducted to understand the influence of varying parameters on cumulative gas and water production. Conclusively, experimental design algorithms generate over 100 simulation cases using the identified sensitive parameters, from which the top ten optimal cases are chosen for EUR prediction. This workflow features two technological innovations: (1) considering the most comprehensive set of reservoir parameters for uncertainty and sensitivity analysis, and (2) considering the matching accuracy of both cumulative production and dynamic production trends when selecting optimal matching cases. This approach was successfully implemented in the C pilot area of the Bowen Basin, Australia. In addition, it offers valuable insights for numerical simulation of unconventional natural gases, such as shale gas.