This paper presents a comprehensive steady-state gas load flow analysis considering all elements and fundamentals of gas networks. The reason is the increasing usage of natural gas carrier in the industry because they have a higher efficiency. In this paper, interdependencies of the mentioned infrastructures are considered in detail including a nonlinear part-load efficiency performance for units. A Jacobian matrix-based framework is proposed to solve the energy flow problem according to the Newton–Raphson technique. This method omits the use of any extra variable and it is able to model any equipment of natural gas network in the Jacobian matrix. A typical 20-node natural gas network is used to evaluate the outputs of the proposed method. In addition, sensitivity analysis is also provided.