Planning is one of the most important components of gas industry production. Most big gas companies usually look for effective planning approaches to accomplish the organizational objectives including cost and time reduction as well as enhancing quality and efficiency. For planning gas refinery production, important parameters including production time, production volume, production cost and storage need to be considered. Planning different units ought to be integrated and coordinated with other departments. This article presents an intensive arithmetic model to determine the production of gas derivatives. The proposed model of this paper is formulated as a mixed integer programming and the resulted problem is solved using NSGA-II algorithm and a hybrid method called BBO/NSGA-II. The problem is also applied for a real-world case study and the results are discussed.