The increasing penetration of renewable energy resources in integrated energy system makes it necessary to consider uncertainty, appropriate investment and scientific operation for energy storage systems, including battery energy storage systems and thermal energy storage systems. In this paper, we propose a two-stage robust planning-operation co-optimization method for Energy Hub considering uncertainties from renewable energy resources as well as multi-load demands, the sizing problem and precise economic model of energy storage systems, in which lifetime loss cost of battery energy storage systems and static loss of thermal energy storage systems are mainly included. Integer variables exist in both stages, thus in the proposed method, dual norms are used to transform uncertain constraints into robust counterparts with which the worst cases are calculated directly, then the problem is merged into a single stage mixed integer linear program. In case studies, optimal planning strategies including multiple storages in Energy Hub are explored under the consideration of sizing problem and uncertainties, and it’s found that simultaneously constructing multiple energy storages can have mutually beneficial relationships and propose more economic planning/operation strategies. Then it’s validated that precise energy storage model can better improve economics and flexibility of Energy Hub both in longer term and short term planning optimization. Meanwhile, convergence speed of the model and solution algorithm in this paper is acceptable for practical planning problem.