Robust optimization of subsurface flow using polynomial chaos and response surface surrogates