In the present paper, we utilize the ensemble Kalman filter (EnKF) method to optimize parameters of the turbulence model Spalart-Allramas (S-A) to simulate flow fields around three wind turbine airfoils (NACA63415, S809, and DU97W300) at the same stall degree. The optimized parameters of the medium-thickness airfoil S809 are further applied to the other two airfoils to investigate their performance on airfoils with various thicknesses. For each airfoil, comparisons to experimental data regarding the pressure coefficient distribution and separation location show that using the optimized parameters can reduce numerical errors significantly. Even though applying the optimized parameters of S809 to the other two can also considerably reduce the numerical errors, the errors are slightly larger than those obtained by the optimized parameters of their own. In addition, the model parameters with significant changes after assimilation are $C_{\mathrm{b}1}$, $C_{\upsilon \text{1}}$, and σ, of which $C_{\mathrm{b}1}$ undergoes the most prominent change and decreases posterior of assimilation for all three airfoils. It is speculated that $C_{\mathrm{b}1}$ is crucial for the model parameters' universality, at least for the present three airfoils, of the model parameters.