Nowadays, the electric power system and natural gas network are becoming increasingly coupled and interdependent. A harmonized integration of natural gas and electricity network with bi-directional energy conversion is expected to accommodate high penetration levels of renewables in terms of system flexibility. This work focuses on the steady-state analysis of the integrated natural gas and electric power system with bi-directional energy conversion. A unified energy flow formulation is developed to describe the nodal balance and branch flow in both systems and it is solved with the Newton-Raphson method. Both the unification of units and the per-unit system are proposed to simplify the system description and to enhance the computation efficiency. The applicability of the proposed method is demonstrated by analyzing an IEEE-9 test system integrated with a 7-node natural gas network. Later, time series of wind power and power load are used to investigate the mitigation effect of the integrated energy system. At last, the effect of wind power and power demand on the output of Power to Gas (P2G) and gas-fired power generation (GPG) has also been investigated. (C) 2016 Elsevier Ltd. All rights reserved.