Optimal operational strategy and planning of a raw natural gas refining complex (RNGRC) is very challenging since it involves highly nonlinear processes, complex thermodynamics, blending, and utility systems. In this paper, we first propose a superstructure integrating a utility system for the RNGRC, involving multiple gas feedstocks and different product specifications. Then, we develop a large-scale nonconvex mixed-integer nonlinear programming (MINLP) optimization model. The This article is protected by copyright. All rights reserved. Finally, the reformulated nonconvex MINLP model is solved to global optimality using state of the art deterministic global optimization approaches. The computational results demonstrate that a significant profit increase is achieved using the proposed approach compared to that from the real operation.