A MULTI-SLACK BUS MODEL FOR BI-DIRECTIONAL ENERGY FLOW ANALYSIS OF INTEGRATED POWER-GAS SYSTEMS