Electricity, gas and heating networks are increasingly coupled to form an integrated energy system (IES), which is beneficial in improving the operational performance of the system during power outages. Generally, IES involves multiple energy operators, making the centralized framework impracticable. To this end, a novel distributed restoration method for IES is proposed in this paper. The IES restoration problem is formulated as a mixed-integer second-order cone programming (MISOCP) model to minimize the total load losses, which considers the coupling effects among the electricity, gas and heating networks. Then, considering the independent management and information privacy of different energy networks, the distributed restoration model is formulated by decomposing shared boundary information using the alternating direction method of multipliers (ADMM). For the sake of the non-convexity caused by binary load state variables, an iterative loop scheme is introduced in ADMM procedure. This allows us to solve a relaxed distributed model in the first stage, and perform iterative loops in the second stage, which can improve the convergence of ADMM and make the distributed solution feasible. Finally, numerical studies are conducted on 33-7-6-node IES and 123-24-32-node IES to validate the effectiveness and feasibility of the proposed method.