A mixed integer linear programming (MILP) model is proposed for the reformation of natural gas pipelines. The model is based on the topology of existing pipelines, the load and pressure at each node and the design factors of the region and minimizes the annual substitution depreciation cost of pipelines, the annual construction depreciation cost of compressor stations and the operating cost of existing compressor stations. Considering the nonlinear pressure drop equations, the model is linearized by a piecewise method and solved by the Gurobi optimizer. Two cases of natural gas pipeline networks with hydrogen injection are presented. Several adjustments are applied to the original natural gas pipeline network to ensure that our design scheme can satisfy the safety and economic requirements of gas transportation. Thus, this work is likely to serve as a decision-support tool for the reformation of pipeline networks with hydrogen injection. (C) 2018 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.