The optimal design of low-pollution emission burners plays an important role in controlling pollutant emissions of industrial equipment, and is crucial for the sustainable development of the national economy and environmental protection. However, many uncertain factors challenge the optimal design of low-pollution emission burners. The Latin hypercube sampling (LHS) method was used to obtain sampling data representing the distribution of the uncertain variable. The training dataset was obtained using the turbulent combustion coupling model. A high-precision sparse polynomial chaos expansion (PCE) model was constructed by the degree-adaptive scheme and least angle regression (LAR) algorithm. Furthermore, the Legendre polynomial is introduced to establish a continuous robust optimization model. The model is carried out by the non-dominated sorting genetic algorithm II (NSGA-II). The results show that the excess air coefficient of 1.227 is optimal. Compared with the excess air coefficient of 1.20 under the discrete robust optimization, the optimal coefficient can further reduce pollutant emissions and bring strong robustness to the ethylene cracking furnace. It has also been proven that the continuous robust optimization scheme improves the optimization granularity. Compared with discrete robust optimization, this method reduces the number of samples by 66.7 %.