To improve the wind resistance performance of large-span coal-storage structures， an aerodynamic shape optimization method for free-form reticulated structures is proposed by utilizing the parametric configuration method based on non-uniform rational B-splines （NURBS） theory and surrogate modeling technique. A secondary development of the program is carried out based on Grasshopper and Fortran to achieve an integrated process for design parameter adjustment， automatic model updating， and aerodynamic response calculation， which significantly enhances the automation level of optimization design. The control points of NURBS curves serve as optimization variables， and the maximum displacement is set as the objective. An optimization model for free-form reticulated structures is established accordingly. On this basis， aerodynamic shape optimizations are performed on tri-cylindrical and spherical coal-storage reticulated structures. The results show that， the proposed method， which combines parametric modeling and a Kriging surrogate model， effectively identifies rational aerodynamic shapes and achieves effective optimization. The maximum displacements of optimal cylindrical shells under wind directions of 0° and 30°are reduced by 32% and 18%， respectively， while the extreme value for spherical shells is reduced by 12%. The wind resistance performance of both types of reticulated shells is improved， and the goal of reducing construction costs is also achieved. Thus， a new approach for wind-resistant design and form selection of large-span reticulated structures can be finally provided.