The modeling and multi-energy flow calculation of an integrated energy system (IES) are the bases of its operation and planning. This paper establishes the models of various energy sub-systems and the coupling equipment for an electricity-gas-thermal IES, and an integrated multi-energy flow calculation model of the IES is constructed. A simplified calculation method for the compressor model in a natural gas network, one which is not included in a loop and works in constant compression ratio mode, is also proposed based on the concept of model reduction. In addition, a numerical conversion method for dealing with the conflict between nominal value and per unit value in the multi-energy flow calculation of IES is described. A case study is given to verify the correctness and speed of the proposed method, and the electricity-gas-thermal coupling interaction characteristics among sub-systems are studied.