Abstract: A method for calculating an optimal energy flow of an integrated electricity-gas system comprises: a convex optimization part of an optimal energy flow model is established according to a fuel cost of each thermal power unit node and a gas supply cost of each nature gas source node, and a convex relaxation formula of a quadratic constraint of the optimal energy flow model is established, and a convex function is expanded by a first order Taylor expansion at a relaxation solution to form an expansion formula; a precision requirement threshold value of a non-convex constraint is given, the threshold value is compared with an unbalance magnitude of a non-convex constraint in the expansion formula; when the unbalance magnitude is greater than the threshold value, the expansion is iterated into a solution model of the integrated electricity-gas hybrid energy system until the unbalance magnitude is not greater than the threshold value.