In the storm surge model, the wind drag coefficient Cd is a critical parameter that has a great influence on the forecast of the storm surge level. In the present study, the effect of various wind drag coefficient parameterizations on the storm surge level is investigated in the Bohai Sea, Yellow Sea and East China Sea for Typhoons 7203 and 7303. Firstly, the impacts of initial values of a and b in the linear expression Cd = (a + b × U10) × 10−3 on the pure storm surge model are evaluated based on the data assimilation method. Results indicate that when a and b (i.e., the wind drag coefficients given by Smith, Wu, Geernaert et al. and Mel et al.) are non-zeros, the performance of the model has little difference, and the result from Wu is slightly better. However, they are superior to the performance of the model adopting zero initial values. Then, we discuss the influences of diverse ways of calculating wind drag coefficients, which are inverted by data assimilation method (including both linear and constant Cd) and given in the form of linear formulas, on simulating pure storm surge level. They show that the data assimilation-based coefficients greatly exceed those of the ordinary coefficient formulas. Moreover, the wind drag coefficient in the linear form is a little better than that in constant form when the data assimilation method is used. Finally, the assessment of the impact of astronomical tides on the storm surge level is conducted, and the simulation demonstrates that the storm surge model, which has the combination of four constituents (M2, S2, K1, O1) and wind drag coefficient inverted by the data assimilation method with the linear Cd, exhibits the best performance.