To comprehensively explore flow control method of profiled end wall for turbine stage at low Reynolds numbers, a surrogate model optimization platform including non-uniform rational B-spline surface parameterization method, support vector regression, and improved chaos particle swarm optimization algorithm is integrated. Optimization designs have been carried out for stator profiled end walls, rotor profiled end wall, and combined end walls, respectively. The results indicate that under the constraint of the output power, the application of various profiled end wall design cases all can effectively improve the aerodynamic performance of the turbine stage. By organizing the flow field of downstream rotor, the profiled end wall of stator can significantly affect the stage efficiency. The flow control benefits of the profiled end wall of the rotor is from the obstruction of the cross migration of the pressure side leg of the horseshoe vortex. The application of profiled end wall on stator has the most practical engineering value. Self-organizing maps and Shapley methods are used to explore potential correlation information of aerodynamic parameters and summarize design experience. The sensitive design variables of profiled end walls are extracted. Based on the local controllability of NURBS surfaces, the regions that affect the stage efficiency are mainly concentrated in the middle of the stator passage, near the stator trailing edge and near the rotor leading edge. The regions with a significant impact on the output power of the turbine stage are near the trailing edge of the rotor and stator. The corresponding design rules of end walls modeling are summarized.