US 12,658,697 B2Grant

Method for calculating optimal energy flow of integrated electricity-gas system based on sequential convex programming

Issue Date:2026-06-16

•7 Claims

•2 Drawing Sheets

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.

Metadata

Assignees

SOUTHEAST UNIVERSITY
ELECTRIC POWER RESEARCH INSTITUTE OF YUNNAN POWER GRID CO., LTD.
YUNNAN POWER GRID CO., LTD.

Inventors

Pengxiang Liu
Zhi Wu
Hai Lu
Wei Gu
Yuping Lu
Enbo Luo
Dada Wang
Qirun Sun
Hao Zhang

Application Information

Application Number:US 18/029,973

Filing Date:2022-04-09

Priority Date:2021-08-11

Art Unit:2115

Classifications

IPC:

H02J3/00H02J103/30

Field of Search:

H02J 3/00H02J 2203/20G06F 30/20G06F 2111/04G06F 2113/14G06Q 50/06

Patent Drawings (2 sheets)

Description

Technical Field

[0001] The present disclosure relates to the technical field of an optimal energy flow of an integrated electricity-gas system, and in particular to a method for calculating an optimal energy flow of an integrated electricity-gas system based on sequential convex programming.

Background

[0002] With the rapid development of energy revolution, the hybrid energy system constructed by the coupling of the electrical power system and the nature gas system has become an important development direction on the optimization of energy structure in China. The multi-energy interlinking hybrid energy system has broken the operation barrier of the traditional energy system, has realized the mutual coupling, substitution and supplement of the multi-energy form, and has promoted the diversified utilization of energy. The calculation of optimal energy flow of the integrated electricity-gas system is one of the important theoretical bases for researching the planning and operation of the system.

[0003] At present, the greatest difficulty in solving this problem is the non-convex constraint in the optimal energy flow model of the integrated electricity-gas system. Although the traditional method for calculating the optimal energy flow based on a convex relaxation strategy has relatively high calculation speed, it frequently fails to satisfy the requirements of calculation accuracy, and the obtained solutions cannot guarantee the feasibility requirements for non-convex constraints. Therefore, it is extremely important for the development of the integrated electricity-gas system to provide an algorithm that can efficiently solve the non-convex optimization problem of the optimal energy flow of the integrated electricity-gas system.

Summary

[0004] The objectives of the present disclosure are to provide a method for calculating an optimal energy flow of an integrated electricity-gas system based on sequential convex programming, by adding the unbalance magnitude of the non-convex constraints into the optimization objective as a penalty factor, thereby guaranteeing the feasibility and tightness of the solutions and guaranteeing the feasibility of the non-convex constraints.

[0005] In order to solve the above problems, the present disclosure is achieved by the following technical solutions.

[0006] Provided in the present disclosure is a method for calculating an optimal energy flow of an integrated electricity-gas system based on sequential convex programming. The method comprises the following steps.

[0007] A convex relaxation formula of a quadratic constraint of an optimal energy flow model is established in an integrated electricity-gas hybrid energy system, and a convex function is expanded by a first order Taylor expansion at a relaxation solution to form an expansion.

[0008] A precision requirement threshold value of a non-convex constraint is given, the threshold value is compared with an unbalance magnitude of the non-convex constraint in the expansion formula.

[0009] When the unbalance magnitude is greater than the threshold value, the expansion formula 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, to obtain a relaxation solution in the solution model.

[0010] Further, iterating the expansion formula into a solution model of the integrated electricity-gas hybrid energy system includes as follows.

[0011] The precision requirement threshold value of the non-convex constraint is given, an unbalance magnitude obtained from a difference of the convex relaxation formula of the quadratic constraint is compared with the threshold value.

[0012] When the unbalance magnitude is not greater than a set precision requirement threshold value of the non-convex constraint, the expansion formula is solved according to an energy flow model.

[0013] When the unbalance magnitude is greater than the set precision requirement threshold value of the non-convex constraint, the first order Taylor expansion at the relaxation solution of the convex function is introduced into the model as a penalty term, a relaxation solution is recalculated with the penalty term of the model;

[0014] Before the relaxation solution with the penalty term of the model is calculated, the unbalance magnitude of an solved expansion formula of the relaxation formula after the penalty term is introduced is compared with the threshold value, when the unbalance magnitude is greater than the set precision requirement threshold value of the non-convex constraint, the penalty term is iterated until the unbalance magnitude is not greater than the set precision requirement threshold value of the non-convex constraint.

[0015] Further, before establishing the convex relaxation formula of the quadratic constraint of the optimal energy flow model, a convex optimization part of the 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 in the integrated electricity-gas hybrid energy system;

[0016] The constraints in the electricity-gas interlinking hybrid energy system include a power flow constraint of an electrical power system and a power flow constraint of a gas system, a power flow model of the electrical power system and a power flow model of the gas system are both quadratic nonlinear models.

[0017] Further, an objective function of the convex optimization part of the optimal energy flow model in the integrated electricity-gas system is established as:

[0018]

\min \sum_{i \in Ω_{c}} [{C_{i}^{a} (p_{i}^{c})}^{2} + C_{i}^{b} p_{i}^{c} + C_{i}^{c}] + \sum_{i \in Ω_{s}} C_{m}^{f} f_{m}^{s},

where sets Ω_cand Ω_srepresent a thermal power unit node set and a nature gas source node set respectively, a variable

[0019]

p_{i}^{c}

represents an active power of the thermal power unit, a variable

[0020]

f_{m}^{s}

represents a gas supply rate of the nature gas source, parameters

[0021]

C_{i}^{a}, C_{i}^{b} and C_{i}^{c}

represent a second-order term coefficient, a first-order term coefficient and a zero-order term coefficient of the fuel cost in the thermal power unit respectively, and

[0022]

C_{m}^{f}

represents a gas supply cost coefficient of the nature gas source.

[0023] Further, a linear part of the power flow constraints of the gas system is as follows:

[0024]

\sum_{{ji}^{\in β_{i}}} (p_{ji} - l_{ji} R_{ji}) - \sum_{ik \in α_{i}} p_{ik} + p_{i}^{c} + p_{i}^{g} = D_{i}^{p}, \forall i \in Ω_{b}

\sum_{ji \in β_{i}} (q_{ji} - l_{ji} X_{ji}) - \sum_{ik \in α_{i}} q_{ik} + q_{i}^{c} + q_{i}^{g} = D_{i}^{q}, \forall i \in Ω_{b}

v_{i} - 2 (p_{ij} R_{ij} + q_{ij} X_{ij}) + l_{ij}, (R_{ij}^{2} + X_{ij}^{2}) = v_{j}, \forall ij \in Ω_{l}

{(V_{i}^{\min})}^{2} \leq v_{i} \leq {(V_{i}^{\max})}^{2}, \forall i \in Ω_{b}

P_{c, i}^{\min} \leq p_{i}^{c} \leq P_{c, i}^{\max}, Q_{c, i}^{\min} \leq q_{i}^{c} \leq Q_{c, i}^{\max}, \forall i \in Ω_{c}

P_{g, i}^{\min} \leq p_{i}^{g} \leq P_{g, i}^{\max}, Q_{g, i}^{\min} \leq q_{i}^{g} \leq Q_{g, i}^{\max}, \forall i \in Ω_{g}

P_{ij}^{\min} \leq p_{ij} \leq P_{ij}^{\max}, \forall i j \in Ω_{l}

Q_{ij}^{\min} \leq q_{ij} \leq Q_{ij}^{\max}, \forall i j \in Ω_{l}

0 \leq l_{ij} \leq {(I_{ij}^{\max})}^{2}, \forall i j \in Ω_{l},

where, sets Ω_b, Ω_land Ω_grepresent an electrical grid node set, an electrical transmission lines set and a gas generator node set respectively, sets α_iand β_irepresent an electrical transmission line set with a node i as a head-terminal node and an electrical transmission line set with the node i as a tail-terminal node respectively, variables p_ij, p_jiand p_ikrepresent active powers flowing over lines ij, ji and ik, q_ij, q_jiand q_ikrepresent reactive powers flowing over the lines ij, ji and ik, l_ijand l_jirepresent squares of currents flowing over the lines ij and ji,

[0025]

p_{i}^{c} and q_{i}^{c}

represent an active power and a reactive power output by the thermal power unit respectively,

[0026]

p_{i}^{g} and q_{i}^{g}

represent an active power and a reactive power output by the gas generator respectively, v_iand v_jrepresent squares of voltage amplitudes of the node i and a node j, parameters R_ijand R_jirepresent resistances on the line ij and the line ji respectively, X_ijand X_jirepresent reactances on the line ij and the line ji respectively,

[0027]

D_{i}^{p} and D_{i}^{q}

represent an active load and a reactive load of the node i respectively,

[0028]

V_{i}^{\min} and V_{i}^{\max}

represent a lower limit and an upper limit of the voltage amplitude of the node i respectively,

[0029]

P_{c, i}^{\min} and P_{c, i}^{\max}

represent a lower limit and an upper limit of the active power output by a thermal power unit i respectively,

[0030]

Q_{c, i}^{\min} and Q_{c, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the thermal power unit i respectively,

[0031]

P_{g, i}^{\min} and P_{g, i}^{\max}

represent a lower limit and an upper limit of the active power output by a gas power unit i respectively,

[0032]

Q_{g, i}^{\min} and Q_{g, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the gas power unit i respectively,

[0033]

P_{ij}^{\min} and P_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the active power on an electrical transmission line ij respectively,

[0034]

Q_{ij}^{\min} and Q_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the reactive power on the electrical transmission line ij respectively,

[0035]

I_{ij}^{\max}

represents a heat stable current value for the line ij.

[0036] Further, a linear part of the power flow constraints of the gas system is as follows:

[0037]

\sum_{im \in δ_{m}} (e_{im} - τ_{im}) - \sum_{mn \in γ_{m}} e_{mn} - \sum_{mn \in Ξ_{m}} f_{mn} + f_{m}^{s} - f_{m}^{g} = D_{m}^{g}, \forall m \in Ω_{n}

τ_{mn} = K_{mn} e_{mn}, \forall mn \in Ω_{k}

{(H_{mn}^{\min})}^{2} π_{m} \leq π_{n} \leq {(H_{mn}^{\max})}^{2} π_{m}, \forall mn \in Ω_{k}

p_{i}^{g} = T_{m} f_{m}^{g}, \forall i, m \in Ω_{g}

{(\prod_{m}^{\min})}^{2} \leq π_{m} \leq {(\prod_{m}^{\max})}^{2}, \forall m \in Ω_{n}

- F_{mn}^{\max} \leq f_{mn} \leq F_{mn}^{\max}, \forall mn \in Ω_{p}

0 \leq e_{mn} \leq E_{mn}^{\max}, \forall mn \in Ω_{k}

0 \leq f_{m}^{g} \leq F_{g, m}^{\max}, \forall m \in Ω_{g}

F_{s, m}^{\min} \leq f_{m}^{s} \leq F_{s, m}^{\max}, \forall m \in Ω_{s},

where sets Ω_n, Ω_pand Ω_krepresent a gas node set, a gas pipeline set and a nature gas compressor set respectively, sets δ_mand γ_mrepresent a gas pipeline set with a node m as head-terminal node and a gas pipeline set with the node m as a tail-terminal node respectively, Ξ_mrepresents a nature gas compressor set with node m as a gas-intake node, e_lmand e_mnrepresent a gas amount flowing through a nature gas compressor im and a gas amount flowing through a nature gas compressor mn, τ_imand τ_mnrepresent a gas amount consumed by the nature gas compressor im and a gas amount consumed by the nature gas compressor mn respectively, f_mnrepresents a gas amount flowing through a gas pipeline mn,

[0038]

f_{m}^{s}

represents a nature gas amount injected into the node m by the nature gas source per unit time,

[0039]

f_{m}^{g}

represents a gas consumption amount by a gas generator connected to the node m per unit time, π_mand π_nrepresent squares of gas pressure values at the node m and a node n respectively, a parameter

[0040]

D_{m}^{g}

represents a gas load at the node m, W_mnrepresents a Weymouth coefficient of a gas transmission pipeline mn, K_mnrepresents a ratio coefficient of a gas compression amount and a gas consumption amount per unit time of the nature gas compressor,

[0041]

H_{mn}^{\min} and H_{mn}^{\max}

represent a lower limit and an upper limit of a compression ratio of the nature gas compressor, T_mrepresents a ratio coefficient of the gas consumption amount and a power generation amount of the gas generator,

[0042]

\prod_{m}^{\min} and \prod_{m}^{\max}

represent a lower limit and an upper limit of a nodal pressure,

[0043]

F_{mn}^{\min} and F_{mn}^{\max}

represent a lower limit and an upper limit of a gas transmission amount of the gas pipeline mn per unit time,

[0044]

E_{mn}^{\max}

represents an upper limit of a compression rate of the nature gas compressor,

[0045]

F_{g, m}^{\max}

represents an upper limit of a gas consumption rate of the gas generator,

[0046]

F_{s, m}^{\min} and F_{s, m}^{\max}

represent a lower limit and an upper limit of a gas supply amount of the nature gas source per unit time.

[0047] Further, the model is expressed in a matrix form:

[0048]

\min x^{T} Q x + c x + d

s . t . Ax \leq b,

where matrices Q, c, d represent a second-order term coefficient, a first-order term coefficient and a constant term matrix in the objective function respectively, a matrix A represents a coefficient matrix in the linear constraint, and b represents a constant term coefficient matrix in the linear constraint.

[0049] Further, the convex relaxation formula of the quadratic constraint of the electricity-gas interlinking hybrid energy system is:

[0050]

{(p_{ij})}^{2} + {(q_{ij})}^{2} \leq v_{i} l_{ij}, \forall ij \in Ω_{l} {W_{mn} (f_{mn})}^{2} \leq π_{m} - π_{n}, \forall mn \in Ω_{p},

the first order Taylor expansion is:

[0051]

{(l_{ij} + v_{i})}^{2} - 8 p_{ij}^{⋆} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} - v_{i}^{*})}^{2} - 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l}, \forall ij \in Ω_{l} (π_{m} - π_{n}) - 2 W_{mn} f_{mn}^{*} f_{mn} + {W_{mn} (f_{mn}^{*})}^{2} \leq μ_{mn}^{p}, \forall mn \in Ω_{p}, where p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{mn}^{*}

respectively represents a given value for a corresponding variable, namely, a value for a line active power, a value for a line reactive power, a square value for a line current, a square value for a node voltage, and a value for a pipeline gas flow optimized and obtained in a previous iteration, variables

[0052]

μ_{ij}^{l} and μ_{mn}^{p}

respectively represent an unbalance magnitude of a non-convex constraint of a corresponding line ij and the gas transmission pipeline mn.

[0053] The present disclosure has the following beneficial effects.

[0054] The existing researches on the method for calculating an optimal energy flow of an integrated electricity-gas system lacks the effective means to deal with non-convex constraints. The traditional convex relaxation techniques are usually difficult to ensure the relaxation tightness, which leads to infeasible solutions. The present disclosure adopts the solution idea of sequential convex programming, and adds the unbalance of non-convex constraint into the optimization objective as a penalty factor to ensure the feasibility and tightness of the solution.

[0055] In order to ensure the feasibility of the solution, the existing researches on the method for calculating an optimal energy flow of an integrated electricity-gas system frequently need to adopt a more complex global optimization algorithm to solve, which leads to a low efficiency of the solution and cannot be applied to the power flow calculation of super large-scale energy system. The power flow calculation method based on sequential convex programming provided in the present disclosure is to solve based on the convex optimization problem in each iteration, therefore on the basis of ensuring the feasibility of the solution, the high efficiency of the solution is taken into account.

[0056] Certainly, the implementation of any product of the present disclosure is not necessary to achieve all the advantages described above at the same time.

Brief Description of the Drawings

[0057] In order to illustrate the technical solutions of the embodiments in the present disclosure more clearly, the following will give a brief introduction to the drawings needed for the description of the embodiments. It will be apparent that the drawings described below are only some embodiments of the present disclosure, for ordinary who skilled in the art, other drawings can be obtained according to these drawings without any creative effort.

[0058] The FIGURE illustrates a schematic diagram of solving the relaxation solution of the model.

Detailed Description of the Embodiments

[0059] The following will clearly and completely describe the technical solutions in the embodiments of the present disclosure in combination with the drawings of the embodiments in the present disclosure. It will be apparent that the described embodiments are only some of the embodiments in the present disclosure, but not all of the embodiments. Based on the embodiments in the present disclosure, all other embodiments obtained by those skilled in the art without creative effort are within the protection scope of the present disclosure.

[0060] With reference to the FIGURE, provided in the present disclosure is a method for calculating an optimal energy flow of an integrated electricity-gas system based on sequential convex programming. The method is implemented through the following steps.

[0061] In Step 1, a convex optimization part of an optimal energy flow model in the integrated electricity-gas system is established.

[0062] In Step 101, an objective function is:

[0063]

\min \sum_{i \in Ω_{c}} [{C_{i}^{a} (p_{i}^{c})}^{2} + C_{i}^{b} p_{i}^{c} + C_{i}^{c}] + \sum_{i \in Ω_{s}} C_{m}^{f} f_{m}^{s}

where sets Ω_cand Ω_srepresent a thermal power unit node set and a nature gas source node set respectively, a variable

[0064]

p_{i}^{c}

represents an active power of the thermal power unit, a variable

[0065]

f_{m}^{s}

represents a gas supply rate of the nature gas source, parameters

[0066]

C_{i}^{a}, C_{i}^{b} and C_{i}^{c}

represent a second-order term coefficient, a first-order term coefficient and a zero-order term coefficient of the fuel cost in the thermal power unit respectively, and

[0067]

C_{m}^{f}

represents a gas supply cost coefficient of the nature gas source.

[0068] In Step 102, a linear part of the power flow constraints of the electrical power system is as follows:

[0069]

\sum_{{ji}^{\in β_{i}}} (p_{ji} - l_{ji} R_{ji}) - \sum_{ik \in α_{i}} p_{ik} + p_{i}^{c} + p_{i}^{g} = D_{i}^{p}, \forall i \in Ω_{b}

\sum_{ji \in β_{i}} (q_{ji} - l_{ji} X_{ji}) - \sum_{ik \in α_{i}} q_{ik} + q_{i}^{c} + q_{i}^{g} = D_{i}^{q}, \forall i \in Ω_{b}

v_{i} - 2 (p_{ij} R_{ij} + q_{ij} X_{ij}) + l_{ij}, (R_{ij}^{2} + X_{ij}^{2}) = v_{j}, \forall ij \in Ω_{l}

{(V_{i}^{\min})}^{2} \leq v_{i} \leq {(V_{i}^{\max})}^{2}, \forall i \in Ω_{b}

P_{c, i}^{\min} \leq p_{i}^{c} \leq P_{c, i}^{\max}, Q_{c, i}^{\min} \leq q_{i}^{c} \leq Q_{c, i}^{\max}, \forall i \in Ω_{c}

P_{g, i}^{\min} \leq p_{i}^{g} \leq P_{g, i}^{\max}, Q_{g, i}^{\min} \leq q_{i}^{g} \leq Q_{g, i}^{\max}, \forall i \in Ω_{g}

P_{ij}^{\min} \leq p_{ij} \leq P_{ij}^{\max}, \forall i j \in Ω_{l}

Q_{ij}^{\min} \leq q_{ij} \leq Q_{ij}^{\max}, \forall i j \in Ω_{l}

0 \leq l_{ij} \leq {(I_{ij}^{\max})}^{2}, \forall i j \in Ω_{l},

where, sets Ω_b, Ω_land Ω_grepresent an electrical grid node set, an electrical transmission lines set and a gas generator node set respectively, sets α_iand β_irepresent an electrical transmission line set with a node i as a head-terminal node and an electrical transmission lines set with the node i as a tail-terminal node respectively, variables p_ij, p_jiand p_ikrepresent active powers flowing over lines ij, ji and ik, q_ij, q_jiand q_ikrepresent reactive powers flowing over the lines ij, ji and ik, l_ijand l_jirepresent squares of currents flowing over the lines ij and ji,

[0070]

p_{i}^{c} and q_{i}^{c}

represent an active power and a reactive power output by the thermal power unit respectively,

[0071]

p_{i}^{g} and q_{i}^{g}

represent an active power and a reactive power output by the gas generator respectively, v_iand v_jrepresent squares of voltage amplitudes of the node i and a node j, parameters R_ijand R_jirepresent resistances on the lines ij and the line ji respectively, X_ijand X_jirepresent reactances on the line ij and the line ji respectively,

[0072]

D_{i}^{p} and D_{i}^{q}

represent an active load and a reactive load of the node i respectively,

[0073]

V_{i}^{\min} and V_{i}^{\max}

represent a lower limit and an upper limit of the voltage amplitude of the node i respectively,

[0074]

P_{c, i}^{\min} and P_{c, i}^{\max}

represent a lower limit and an upper limit of the active power output by a thermal power unit i respectively,

[0075]

Q_{c, i}^{\min} and Q_{c, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the thermal power unit i respectively,

[0076]

P_{g, i}^{\min} and P_{g, i}^{\max}

represent a lower limit and an upper limit of the active power output by a gas power unit i respectively,

[0077]

Q_{g, i}^{\min} and Q_{g, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the gas power unit i respectively,

[0078]

P_{ij}^{\min} and P_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the active power on an electrical transmission line ij respectively,

[0079]

Q_{ij}^{\min} and Q_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the reactive power on the electrical transmission line ij respectively,

[0080]

I_{ij}^{\max}

represents a heat stable current value for the line ij.

[0081] In Step 103, a linear part of the power flow constraints of the gas system is as follows:

[0082]

\sum_{im \in δ_{m}} (e_{im} - τ_{im}) - \sum_{mn \in γ_{m}} e_{mn} - \sum_{mn \in Ξ_{m}} f_{mn} + f_{m}^{s} - f_{m}^{g} = D_{m}^{g}, \forall m \in Ω_{n}

τ_{mn} = K_{mn} e_{mn}, \forall mn \in Ω_{k}

{(H_{mn}^{\min})}^{2} π_{m} \leq π_{n} \leq {(H_{mn}^{\max})}^{2} π_{m}, \forall mn \in Ω_{k}

p_{i}^{g} = T_{m} f_{m}^{g}, \forall i, m \in Ω_{g}

{(\prod_{m}^{\min})}^{2} \leq π_{m} \leq {(\prod_{m}^{\max})}^{2}, \forall m \in Ω_{n}

- F_{mn}^{\max} \leq f_{mn} \leq F_{mn}^{\max}, \forall mn \in Ω_{p}

0 \leq e_{mn} \leq E_{mn}^{\max}, \forall mn \in Ω_{k}

0 \leq f_{m}^{g} \leq F_{g, m}^{\max}, \forall m \in Ω_{g}

F_{s, m}^{\min} \leq f_{m}^{s} \leq F_{s, m}^{\max}, \forall m \in Ω_{s},

where sets Ω_n, Ω_pand Ω_krepresent a gas node set, a gas pipeline set and a nature gas compressor set respectively, sets δ_mand γ_mrepresent a gas pipeline set with a node m as a head-terminal node and a gas pipeline set with the node m as a tail-terminal node respectively, Ξ_mrepresents a nature gas compressor set with node m as a gas-intake node, e_lmand e_mnrepresent a gas amount flowing through a nature gas compressor im and a gas amount flowing through a nature gas compressor mn, τ_imand τ_mnrepresent a gas amount consumed by the nature gas compressor im and a gas amount consumed by the nature gas compressor mn respectively, f_mnrepresents a gas amount flowing through a gas pipeline mn,

[0083]

f_{m}^{s}

represents a nature gas amount injected into the node m by the nature gas source per unit time,

[0084]

f_{m}^{g}

represents a gas consumption amount by a gas generator connected to the node m per unit time, π_mand π_nrepresent squares of gas pressure values at the node m and a node n respectively, a parameter

[0085]

D_{m}^{g}

[0086]

H_{m n}^{\min} and H_{m n}^{\max}

[0087]

\prod_{m}^{\min} and \prod_{m}^{\max}

represent a lower limit and an upper limit of a nodal pressure,

[0088]

F_{m n}^{\min} and F_{m n}^{\max}

represent a lower limit and an upper limit of a gas transmission amount of the gas pipeline mn per unit time,

[0089]

E_{m n}^{\max}

represents an upper limit of a compression rate of the nature gas compressor,

[0090]

F_{g, m}^{\max}

represents an upper limit of a gas consumption rate of the gas generator,

[0091]

F_{s, m}^{\min} and F_{s, m}^{\max}

represent a lower limit and an upper limit of a gas supply amount of the nature gas source per unit time.

[0092] In Step 104, the above model is expressed in a matrix form as follows:

[0093]

\min x^{T} Q x + c x + d

s . t . Ax \leq b,

[0094] In Step 2, a convex relaxation formula of a quadratic constraint of the optimal energy flow model and its first order Taylor expansion are established.

[0095] Specifically, the convex relaxation formula of the quadratic constraint of the optimal energy flow model is established in an integrated electricity-gas hybrid energy system, and a convex function is expanded by the first order Taylor expansion at a relaxation solution to form expansion formulas, the expansion formulas described below are Taylor expansion formulas. In other words. That is, linear constraints are established by the optimal energy flow model based on the first-order Taylor expansion instead of the quadratic constraints, and solutions of the model at this time are approximate solutions.

[0096] A precision requirement threshold value of a non-convex constraint is given, then the set threshold value is compared with an unbalance magnitude of the non-convex constraint in the expansion formula, and whether to solve the model is determined according to the comparison results.

[0097] When the unbalance magnitude is greater than the threshold value, the expansion formula 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, to obtain a relaxation solution in the solution model.

[0098] For the unbalance magnitude mentioned above, the approximate solutions may not satisfy an original quadratic constraint, and the shortest distance between the approximate solutions and a feasible domain is the unbalance magnitude.

[0099] Specifically, iterating the expansion formula into a solution model of the integrated electricity-gas hybrid energy system includes as follows.

[0100] The precision requirement threshold value of the non-convex constraint is given, an unbalance magnitude obtained from a difference of the convex relaxation formula of the quadratic constraint is compared with the threshold value.

[0101] When the unbalance magnitude is not greater than a set precision requirement threshold value of the non-convex constraint, the expansion formula is solved according to an energy flow model.

[0102] When the unbalance magnitude is greater than the set precision requirement threshold value of the non-convex constraint, the first order Taylor expansion at the relaxation solution of the convex function is introduced into the model as a penalty term, a relaxation solution is recalculated with the penalty term of the model.

[0103] Before the relaxation solution with the penalty term of the model is calculated, the unbalance magnitude of an solved expansion formula of the relaxation formula after the penalty term is introduced is compared with the threshold value, when the unbalance magnitude is greater than the set precision requirement threshold value of the non-convex constraint, the penalty term is iterated until the unbalance magnitude is not greater than the set precision requirement threshold value of the non-convex constraint.

[0104] The actual calculation processes are as follows.

[0105] In Step 201, the convex relaxation formula of the quadratic constraint of the electricity-gas interconnected hybrid energy system is:

[0106]

{(p_{ij})}^{2} + {(q_{ij})}^{2} \leq v_{i} l_{ij}, \forall ij \in Ω_{l} {W_{mn} (f_{mn})}^{2} \leq π_{m} - π_{n}, \forall mn \in Ω_{p} .

[0107] In Step 202, a first order Taylor expansion formula of the above constraint is:

[0108]

{(l_{ij} + v_{i})}^{2} - 8 p_{ij}^{*} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} + v_{i}^{*})}^{2}

- 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l}, \forall ij \in Ω_{l}

(π_{m} - π_{n}) - 2 W_{m n} f_{m n}^{*} f_{m n} + {W_{m n} (f_{m n}^{*})}^{2} \leq μ_{m n}^{p}, \forall m n \in Ω_{p}, where

p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{m n}^{*}

[0109]

μ_{ij}^{l} and μ_{m n}^{p}

respectively represents an unbalance magnitude of a non-convex constraint of a corresponding line ij and the gas transmission pipeline mn.

[0110] In Step 3, the sequential convex programming strategy is used to solve the optimal energy flow model iteratively.

[0111] In Step 301, the precision requirement threshold value ε of the non-convex constraint is set, and the counter is set to be k=0.

[0112] In Step 302, the model is solved

[0113]

\min x^{T} Q x + c x + d

s . t . Ax \leq b

to obtain the current solutions

[0114]

p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{m n}^{*}

of the variables

[0115]

p_{ij}, q_{ij}, l_{ij}, v_{i} and f_{m n} .

[0116] In Step 303, whether the unbalance magnitude of the non-convex constraint is satisfied is determined

[0117]

v_{i}^{*} l_{ij}^{*} - {(p_{ij}^{*})}^{2} - {(q_{ij}^{*})}^{2} \leq ε, \forall i j \in Ω_{l}

π_{m}^{*} - π_{n}^{*} - {W_{m n} (f_{m n}^{*})}^{2} \leq ε, \forall m n \in Ω_{p}

when the accuracy requirement is satisfied, the calculation results are output; if not, Step 304 is executed.

[0118] In Step 304, the model is solved

[0119]

\min x^{T} Q x + c x + d + 2^{k} (\sum_{ij \in Ω_{l}} μ_{ij}^{l} + \sum_{mn \in Ω_{p}} μ_{mn}^{p})

s . t . Ax \leq b

{(l_{ij} + v_{i})}^{2} - 8 p_{ij}^{*} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} - v_{i}^{*})}^{2} - 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l}, \forall ij \in Ω_{l}

(π_{m} - π_{n}) - 2 W_{mn} f_{mn}^{*} f_{mn} + {W_{mn} (f_{mn}^{*})}^{2} \leq μ_{mn}^{p}, \forall mn \in Ω_{p} .

[0120] In Step 305, whether the unbalance magnitude of the non-convex constraint is satisfied is determined

[0121]

v_{i}^{*} l_{ij}^{*} - {(p_{ij}^{*})}^{2} - {(q_{ij}^{*})}^{2} \leq ε, \forall i j \in Ω_{l}

π_{m}^{*} - π_{n}^{*} - {W_{m n} (f_{m n}^{*})}^{2} \leq ε, \forall m n \in Ω_{p}

when accuracy requirement is satisfied, the calculation results are output; if not, the current solutions

[0122]

p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{m n}^{*}

of the variables

[0123]

p_{ij}, q_{ij}, l_{ij}, v_{i} and f_{m n}

are obtained, and Step 304 is executed.

[0124] The following is an example of testing the integrated electricity-gas hybrid energy system, whose system parameters are as shown in Table 1 to Table 9.

[0125]

TABLE 1

Node parameters for an electrical power system
			Lower	Upper		Load
	Active	Reactive	voltage	voltage	Voltage	rejection
Serial	load	load	limit	limit	reference	cost
number	(MW)	(MVAr)	(p.u.)	(p.u.)	(kV)	($/MWh)

0	97.6	44.2	0.94	1.06	345	1000
1	0	0	0.94	1.06	345	1000
2	322	2.4	0.94	1.06	345	1000
3	500	184	0.94	1.06	345	1000
4	0	0	0.94	1.06	345	1000
5	0	0	0.94	1.06	345	1000
6	233.8	84	0.94	1.06	345	1000
7	522	176.6	0.94	1.06	345	1000
8	6.5	−66.6	0.94	1.06	345	1000
9	0	0	0.94	1.06	345	1000
10	0	0	0.94	1.06	345	1000
11	8.53	88	0.94	1.06	345	1000
12	0	0	0.94	1.06	345	1000
13	0	0	0.94	1.06	345	1000
14	320	153	0.94	1.06	345	1000
15	329	32.3	0.94	1.06	345	1000
16	0	0	0.94	1.06	345	1000
17	158	30	0.94	1.06	345	1000
18	0	0	0.94	1.06	345	1000
19	680	103	0.94	1.06	345	1000
20	274	115	0.94	1.06	345	1000
21	0	0	0.94	1.06	345	1000
22	247.5	84.6	0.94	1.06	345	1000
23	308.6	−92.2	0.94	1.06	345	1000
24	224	47.2	0.94	1.06	345	1000
25	139	17	0.94	1.06	345	1000
26	281	75.5	0.94	1.06	345	1000
27	206	27.6	0.94	1.06	345	1000
28	283.5	26.9	0.94	1.06	345	1000
29	0	0	0.94	1.06	345	1000
30	9.2	4.6	0.94	1.06	345	1000
31	0	0	0.94	1.06	345	1000
32	0	0	0.94	1.06	345	1000
33	0	0	0.94	1.06	345	1000
34	0	0	0.94	1.06	345	1000
35	0	0	0.94	1.06	345	1000
36	0	0	0.94	1.06	345	1000
37	0	0	0.94	1.06	345	1000
38	1104	250	0.94	1.06	345	1000

[0126]

TABLE 2

Node parameters for a gas system
		Lower	Upper
Serial	Gas load	pressure	pressure	Load rejection
number	(1000 m3/h)	limit (bar)	limit (bar)	cost ($/1000 m3)

0	0	40	70	1000
1	0	40	70	1000
2	0	40	70	1000
3	100	40	70	1000
4	120	40	60	1000
5	80	40	60	1000
6	0	40	70	1000
7	0	40	70	1000
8	0	40	70	1000
9	0	40	70	1000
10	0	40	70	1000

[0127]

TABLE 3

Line parameters for an electrical power system
					Suscep-	Transmission
Serial	Head	Tail	Resistance	Reactance	tance	power
number	node	node	(p.u.)	(p.u.)	(p.u.)	(MVA)

0	0	1	0.0035	0.0411	0.6987	600
1	0	38	0.001	0.025	0.75	1000
2	1	2	0.0013	0.0151	0.2572	500
3	1	24	0.007	0.0086	0.146	500
4	1	29	0	0.0181	0	900
5	2	3	0.0013	0.0213	0.2214	500
6	2	17	0.0011	0.0133	0.2138	500
7	3	4	0.0008	0.0128	0.1342	600
8	3	13	0.0008	0.0129	0.1382	500
9	4	5	0.0002	0.0026	0.0434	1200
10	4	7	0.0008	0.0112	0.1476	900
11	5	6	0.0006	0.0092	0.113	900
12	5	10	0.0007	0.0082	0.1389	480
13	5	30	0	0.025	0	1800
14	6	7	0.0004	0.0046	0.078	900
15	7	8	0.0023	0.0363	0.3804	900
16	8	38	0.001	0.025	1.2	900
17	9	10	0.0004	0.0043	0.0729	600
18	9	12	0.0004	0.0043	0.0729	600
19	9	31	0	0.02	0	900
20	11	10	0.0016	0.0435	0	500
21	11	12	0.0016	0.0435	0	500
22	12	13	0.0009	0.0101	0.1723	600
23	13	14	0.0018	0.0217	0.366	600
24	14	15	0.0009	0.0094	0.171	600
25	15	16	0.0007	0.0089	0.1342	600
26	15	18	0.0016	0.0195	0.304	600
27	15	20	0.0008	0.0135	0.2548	600
28	15	23	0.0003	0.0059	0.068	600
29	16	17	0.0007	0.0082	0.1319	600
30	16	26	0.0013	0.0173	0.3216	600
31	18	19	0.0007	0.0138	0	900
32	18	32	0.0007	0.0142	0	900
33	19	33	0.0009	0.018	0	900
34	20	21	0.0008	0.014	0.2565	900
35	21	22	0.0006	0.0096	0.1846	600
36	21	34	0	0.0143	0	900
37	22	23	0.0022	0.035	0.361	600
38	22	35	0.0005	0.0272	0	900
39	24	25	0.0032	0.0323	0.531	600
40	24	36	0.0006	0.0232	0	900
41	25	26	0.0014	0.0147	0.2396	600
42	25	27	0.0043	0.0474	0.7802	600
43	25	28	0.0057	0.0625	1.029	600
44	27	28	0.0014	0.0151	0.249	600
45	28	37	0.0008	0.0156	0	1200

[0128]

TABLE 4

Pipeline parameters for a gas system
			Minimum	Maximum	Weymouth
			transmission	transmission	coefficient
Serial	Head	Tail	rate	rate	((bar/(1000
number	node	node	(1000 m3/h)	(1000 m3/h)	m3/h)){circumflex over ( )}2)

0	0	1	−1100	1100	0.058913
1	6	7	−1100	1100	0.058913
2	2	8	−1100	1100	0.061036
3	7	3	−1100	1100	0.066929
4	7	9	−1100	1100	0.077628
5	8	9	−1100	1100	0.061036
6	10	4	−1100	1100	0.063635
7	10	5	−1100	1100	0.071318

[0129]

TABLE 5

Parameters for a gas valve
			Upper limit	Minimum	Maximum
			of pressure	transmission	transmission
Serial	Head	Tail	difference	rate	rate
number	node	node	(bar)	(1000 m3/h)	(1000 m3/h)

0	6	8	120	−1100	1100

[0130]

TABLE 6

Parameters for a nature gas compressor
					Minimum	Maximum
			Minimum	Maximum	transmission	transmission	Gas-
Serial	Head	Tail	compression	compression	rate	rate	driven
number	node	node	ratio	ratio	(1000 m3/h)	(1000 m3/h)	ratio

0	1	6	1.2	2.8	0	1100	0.035
1	9	10	1.2	2.8	0	1100	0.035

[0131]

TABLE 7

Parameters for a thermal power unit
		Minimum	Maximum	Minimum	Maximum
		active	active	reactive	reactive
Serial		power	power	power	power	c2	c1	c0
number	Node	(MW)	(MW)	(MVAr)	(MVAr)	($/MWh2)	($/MWh)	($/h)

0	29	0	1040	140	400	0.01	0.3	0.2
1	30	0	646	−100	300	0.01	0.3	0.2
2	32	0	652	0	250	0.01	0.3	0.2
3	33	0	508	0	167	0.01	0.3	0.2
4	34	0	687	−100	300	0.01	0.3	0.2
5	35	0	580	0	240	0.01	0.3	0.2
6	36	0	564	0	250	0.01	0.3	0.2
7	38	0	1100	−100	300	0.01	0.3	0.2

[0132]

TABLE 8

Parameters for nature gas source
		Minimum gas	Maximum gas
Serial		supply rate	supply rate	Gas supply cost
number	Node	(1000 m3/h)	(1000 m3/h)	($/1000 m3)

0	0	50	750	88.80555
1	1	0	500	88.80555
2	2	100	500	88.80555

[0133]

TABLE 9

Parameters for a gas generator
			Minimum	Maximum	Minimum	Maximum
	Electric		active	active	reactive	reactive	Conversion
Serial	power	Gas	power	power	power	power	ratio
number	node	node	(MW)	(MW)	(MVAr)	(MVAr)	(1000 m3/MW)

0	31	3	0	725	150	300	0.25
1	37	1	0	865	−150	300	0.25

[0134] Taking the above electricity-gas interlinking hybrid energy system as an example, a calculation for optimal energy flow is carried out for a unit time section, a time scale is set to one hour, and a threshold value for the unbalance magnitude of the non-convex constraint is set to 0.1%. The main technical indexes of the algorithm are as shown in Table 10.

[0135]

TABLE 10

The main technical indexes
Objective function value	Unbalance magnitude	Calculation time

17647.89 RMB	0.06945%	0.194031 seconds

[0136] The results show that the method for calculating an optimal energy flow of an integrated electricity-gas system based on sequential convex programming provided by the present disclosure can effectively obtain the objective function value of the problem, and the calculation time is only 0.194 seconds on the premise of ensuring that the maximum value for the unbalance magnitude of the non-convex constraint does not exceed 0.1% of the threshold value, and the feasibility and efficiency of the optimal energy flow algorithm are considered.

[0137] In the description of this specification, the reference terms “an embodiment”, “example”, “specific example”, etc., are used to refer to the specific features, structures, materials or characteristics described in conjunction with the embodiments or examples, which are included in at least one embodiment or example of the present disclosure. In this specification, schematic representations of the above terms do not necessarily refer to identical embodiments or examples. Furthermore, the specific features, structures, materials or characteristics described may be combined in an appropriate manner in any one or more embodiments or examples.

[0138] The above disclosed preferred embodiments of the present disclosure are only intended to help explain the present disclosure. The preferred embodiments do not describe all of the details and do not limit the present disclosure to the specific embodiments described. It will be apparent that many modifications and changes can be made can be made according to the contents of this specification. These embodiments are selected and described in details in this specification to better explain the principle and practical application of the present disclosure, so that those skilled in the art would better understand and utilize the present disclosure. The present disclosure is only limited by the claims and its full scope and equivalents.

Claims

What is claimed is:

1. A method for calculating an optimal energy flow of an integrated electricity-gas system based on sequential convex programming, wherein integrated electricity-gas system includes an electrical power system having a plurality of thermal power unit nodes and a gas system having nature gas source nodes, and the method comprises following steps:

establishing an electrical power flow model of the electrical power system and a gas power flow model of the gas system;

determining a predetermined parameter of each thermal power unit node of the electrical power system and a predetermined parameter of each nature gas source node of the gas system in the integrated electricity-gas hybrid energy system;

establishing a convex optimization part of the optimal energy flow model according to the predetermined parameter of each thermal power unit node and the predetermined parameter of each nature gas source node;

establishing, in an integrated electricity-gas hybrid energy system, a convex relaxation formula of a quadratic constraint of the optimal energy flow model, and expanding, by a first order Taylor expansion, a convex function at a relaxation solution to form an expansion formula;

giving a precision requirement threshold value of a non-convex constraint, comparing the threshold value with an unbalance magnitude of the non-convex constraint in the expansion formula;

iterating, in response to the unbalance magnitude being greater than the threshold value, the expansion formula into a solution model of the integrated electricity-gas hybrid energy system until the unbalance magnitude is not greater than the threshold value, to obtain a relaxation solution in the solution model,

wherein, constraints in the electricity-gas interlinking hybrid energy system include a power flow constraint of the electrical power system and a power flow constraint of the gas system, the electrical power flow model of the electrical power system and the gas power flow model of the gas system are both quadratic nonlinear models,

wherein the optimal energy flow of the integrated electricity-gas system is efficiently determined based on the obtained relaxation solution, and wherein the integrated electricity-gas system is controlled according to the determined optimal energy flow.

2. The method for calculating the optimal energy flow of the integrated electricity-gas system based on sequential convex programming according to claim 1, wherein the iterating the expansion formula into a solution model of the integrated electricity-gas hybrid energy system includes:

giving, the precision requirement threshold value of the non-convex constraint, comparing an unbalance magnitude obtained from a difference of the convex relaxation formula of the quadratic constraint with the threshold value;

solving, in response to the unbalance magnitude being not greater than a set precision requirement threshold value of the non-convex constraint, the expansion formula according to an energy flow model;

introducing, in response to the unbalance magnitude being greater than the set precision requirement threshold value of the non-convex constraint, the first order Taylor expansion at the relaxation solution of the convex function into the model as a penalty term, recalculating, the model relaxation solution with the penalty term;

comparing, before calculating the model relaxation solution with the penalty term, an unbalance magnitude of an solved expansion formula of the relaxation formula with the penalty term with the threshold value, iterating, in response to the unbalance magnitude being greater than the set precision requirement threshold value of the non-convex constraint, the penalty term until the unbalance magnitude is not greater than the set precision requirement threshold value of the non-convex constraint.

3. The method for calculating the optimal energy flow of the integrated electricity-gas system based on sequential convex programming according to claim 1, wherein establishing, an objective function of the convex optimization part of the optimal energy flow model in the integrated electricity-gas system:

\min \sum_{i \in Ω_{c}} [{C_{i}^{a} (p_{i}^{c})}^{2} + C_{i}^{b} p_{i}^{c} + C_{i}^{c}] + \sum_{i \in Ω_{s}} C_{m}^{f} f_{m}^{s},

where sets Ω_cand Ω_srepresent a thermal power unit node set and a nature gas source node set respectively, a variable

p_{i}^{c}

represents an active power of the thermal power unit, a variable

f_{m}^{s}

represents a gas supply rate of the nature gas source, parameters

C_{i}^{a}, C_{i}^{b} and C_{i}^{c}

represent a second-order term coefficient, a first-order term coefficient and a zero-order term coefficient of the fuel cost in the thermal power unit respectively, and

C_{m}^{f}

represents a gas supply cost coefficient of the nature gas source.

4. The method for calculating the optimal energy flow of the integrated electricity-gas system based on sequential convex programming according to claim 3, wherein a linear part of the power flow constraint of the electrical power system is as follows:

\sum_{ji \in β_{i}} (p_{ji} - l_{ji} R_{ji}) - \sum_{ik \in α_{i}} p_{ik} + p_{i}^{c} + p_{i}^{g} = D_{i}^{p}, \forall i \in Ω_{b}

\sum_{ji \in β_{i}} (q_{ji} - l_{ji} X_{ji}) - \sum_{ik \in α_{i}} q_{ik} + q_{i}^{c} + q_{i}^{g} = D_{i}^{q}, \forall i \in Ω_{b}

v_{i} - 2 (p_{ij} R_{ij} + q_{ij} X_{ij}) + l_{ij} (R_{ij}^{2} + X_{ij}^{2}) = v_{j}, \forall ij \in Ω_{l}

{(V_{i}^{\min})}^{2} \leq v_{i} \leq {(V_{i}^{\max})}^{2}, \forall i \in Ω_{b}

P_{c, i}^{\min} \leq p_{i}^{c} \leq P_{c, i}^{\max}, Q_{c, i}^{\min} \leq q_{i}^{c} \leq Q_{c, i}^{\max}, \forall i \in Ω_{c}

P_{g, i}^{\min} \leq p_{i}^{g} \leq P_{g, i}^{\max}, Q_{g, i}^{\min} \leq q_{i}^{g} \leq Q_{g, i}^{\max}, \forall i \in Ω_{g}

P_{ij}^{\min} \leq p_{ij} \leq P_{ij}^{\max}, \forall ij \in Ω_{l}

Q_{ij}^{\min} \leq q_{ij} \leq Q_{ij}^{\max}, \forall ij \in Ω_{l}

0 \leq l_{ij} \leq {(I_{ij}^{\max})}^{2}, \forall ij \in Ω_{l},

where sets Ω_b, Ω_l, and Ω_grepresent an electrical grid node set, an electrical transmission line set and a gas generator node set respectively, sets α_iand β_irepresent an electrical transmission line set with a node i as a head-terminal node and an electrical transmission line set with the node i as a tail-terminal node respectively, variables p_ij, p_jiand p_ikrepresent active powers flowing over lines ij, ji and ik, q_ij, q_jiand q_ikrepresent reactive powers flowing over the lines ij, ji and ik, l_ijand l_jirepresent squares of currents flowing over the lines ij and ji,

p_{i}^{c} and q_{i}^{c}

represent an active power and a reactive power output by the thermal power unit respectively,

p_{i}^{g} and q_{i}^{g}

D_{i}^{p} and D_{i}^{q}

represent an active load and a reactive load of the node i respectively,

V_{i}^{\min} and V_{i}^{\max}

represent a lower limit and an upper limit of the voltage amplitude of the node i respectively,

P_{c, i}^{\min} and P_{c, i}^{\max}

represent a lower limit and an upper limit of the active power output by a thermal power unit i respectively,

Q_{c, i}^{\min} and Q_{c, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the thermal power unit i respectively,

P_{g, i}^{\min} and P_{g, i}^{\max}

represent a lower limit and an upper limit of the active power output by a gas power unit i respectively,

Q_{g, i}^{\min} and Q_{g, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the gas power unit i respectively,

P_{ij}^{\min} and P_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the active power on an electrical transmission line ij respectively,

Q_{ij}^{\min} and Q_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the reactive power on the electrical transmission line ij respectively,

I_{ij}^{\max}

represents a heat stable current value for the line ij.

5. The method for calculating the optimal energy flow of the integrated electricity-gas system based on sequential convex programming according to claim 4, wherein a linear part of the power flow constraint of the gas system is as follows:

\sum_{im \in δ_{m}} (e_{im} - τ_{im}) - \sum_{mn \in γ_{m}} e_{mn} - \sum_{mn \in Ξ_{m}} f_{mn} + f_{m}^{s} - f_{m}^{g} = D_{m}^{g}, \forall m \in Ω_{n}

τ_{mn} = K_{mn} e_{mn}, \forall mn \in Ω_{k}

{(H_{mn}^{\min})}^{2} π_{m} \leq π_{n} \leq {(H_{mn}^{\max})}^{2} π_{m}, \forall mn \in Ω_{k}

p_{i}^{g} = T_{m} f_{m}^{g}, \forall i, m \in Ω_{g}

{(\prod_{m}^{\min})}^{2} \leq π_{m} \leq {(\prod_{m}^{\max})}^{2}, \forall m \in Ω_{n}

- F_{mn}^{\max} \leq f_{mn} \leq F_{mn}^{\max}, \forall mn \in Ω_{p}

0 \leq e_{mn} \leq E_{mn}^{\max}, \forall mn \in Ω_{k}

0 \leq f_{m}^{g} \leq F_{g, m}^{\max}, \forall m \in Ω_{g}

F_{s, m}^{\min} \leq f_{m}^{s} \leq F_{s, m}^{\max}, \forall m \in Ω_{s},

where sets Ωⁿ, Ω_pand Ω_krepresent a gas node set, a gas pipeline set and a nature gas compressor set respectively, sets δ_mand γ_mrepresent a gas pipeline set with a node m as a head-terminal node and a gas pipeline set with the node m as a tail-terminal node respectively, Ξ_mrepresents a nature gas compressor set with the node m as a gas-intake node, e_lmand e_mnrepresent a gas amount flowing through a nature gas compressor im and a gas amount flowing through a nature gas compressor mn, τ_imand τ_mnrepresent a gas amount consumed by the nature gas compressor im and a gas amount consumed by the nature gas compressor mn respectively, f_mnrepresents a gas amount flowing through a gas pipeline mn,

f_{m}^{s}

represents a nature gas amount injected into the node m by the nature gas source per unit time,

f_{m}^{g}

represents a gas consumption amount by a gas generator connected to the node m per unit time, π_mand π_nrepresent squares of gas pressure values at the node m and a node n respectively, a parameter

D_{m}^{g}

H_{mn}^{\min} and H_{mn}^{\max}

\prod_{m}^{\min} and \prod_{m}^{\max}

represent a lower limit and an upper limit of a nodal pressure, generator,

F_{mn}^{\min} and F_{mn}^{\max}

represent a lower limit and an upper limit of a gas transmission amount of the gas pipeline mn per unit time,

E_{mn}^{\max}

represents an upper limit of a compression rate of the nature gas compressor,

F_{g, m}^{\max}

represents an upper limit of a gas consumption rate of the gas generator,

F_{s, m}^{\min} and F_{s, m}^{\max}

represent a lower limit and an upper limit of a gas supply amount of the nature gas source per unit time.

6. The method for calculating the optimal energy flow of the integrated electricity-gas system based on sequential convex programming according to claim 5, wherein the model is expressed in a matrix form:

min x^TQx+cx+d

s.t. Ax≤b,

7. The method for calculating the optimal energy flow of the integrated electricity-gas system based on sequential convex programming according to claim 3, wherein the convex relaxation formula of the quadratic constraint of the electricity-gas interlinking hybrid energy system is:

(p_ij)²+(q_ij)²≤v_il_ij, ∀ij∈Ω_l

W_mn(f_mn)²≤π_m−π_n, ∀mN∈Ω_p,

the first order Taylor expansion is:

{(I_{ij} + v_{i})}^{2} - 8 p_{ij}^{*} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} - v_{i}^{*})}^{2} - 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l},

\forall ij \in Ω_{l} (π_{m} - π_{n}) - 2 W_{mn} f_{mn}^{*} f_{mn} + W_{mn} {(f_{mn}^{*})}^{2} \leq μ_{mn}^{p}, ⁠ \forall mn \in Ω_{p},

w here p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{mn}^{*}

respectively represent a given value for a corresponding variable, namely, a value for a line active power, a value for a line reactive power, a square value for a line current, a square value for a node voltage, and a value for a pipeline gas flow optimized and obtained in a previous iteration, variables

μ_{ij}^{l} and μ_{mn}^{p}

respectively represent an unbalance magnitude of a non-convex constraint of a corresponding line ij and the gas transmission pipeline mn.

Patent Citations (5)

Patent	Date	Inventor	Cited By
US2017/0184640(A1)	2017-06-01	Pileggi et al.	Applicant
CN108846507(A)	2018-11-01		Examiner
CN110729766(A)	2020-01-01		Examiner
CN112713615(A)	2021-04-01		Examiner
CN113627021(A)	2021-11-01		Applicant

Abstract

Description

Technical Field

Background

Summary

[0005] In order to solve the above problems, the present disclosure is achieved by the following technical solutions.

[0008] A precision requirement threshold value of a non-convex constraint is given, the threshold value is compared with an unbalance magnitude of the non-convex constraint in the expansion formula.

[0010] Further, iterating the expansion formula into a solution model of the integrated electricity-gas hybrid energy system includes as follows.

[0012] When the unbalance magnitude is not greater than a set precision requirement threshold value of the non-convex constraint, the expansion formula is solved according to an energy flow model.

[0017] Further, an objective function of the convex optimization part of the optimal energy flow model in the integrated electricity-gas system is established as:

[0018]

\min \sum_{i \in Ω_{c}} [{C_{i}^{a} (p_{i}^{c})}^{2} + C_{i}^{b} p_{i}^{c} + C_{i}^{c}] + \sum_{i \in Ω_{s}} C_{m}^{f} f_{m}^{s},

where sets Ω_cand Ω_srepresent a thermal power unit node set and a nature gas source node set respectively, a variable

[0019]

p_{i}^{c}

represents an active power of the thermal power unit, a variable

[0020]

f_{m}^{s}

represents a gas supply rate of the nature gas source, parameters

[0021]

C_{i}^{a}, C_{i}^{b} and C_{i}^{c}

represent a second-order term coefficient, a first-order term coefficient and a zero-order term coefficient of the fuel cost in the thermal power unit respectively, and

[0022]

C_{m}^{f}

represents a gas supply cost coefficient of the nature gas source.

[0023] Further, a linear part of the power flow constraints of the gas system is as follows:

[0024]

\sum_{{ji}^{\in β_{i}}} (p_{ji} - l_{ji} R_{ji}) - \sum_{ik \in α_{i}} p_{ik} + p_{i}^{c} + p_{i}^{g} = D_{i}^{p}, \forall i \in Ω_{b}

\sum_{ji \in β_{i}} (q_{ji} - l_{ji} X_{ji}) - \sum_{ik \in α_{i}} q_{ik} + q_{i}^{c} + q_{i}^{g} = D_{i}^{q}, \forall i \in Ω_{b}

v_{i} - 2 (p_{ij} R_{ij} + q_{ij} X_{ij}) + l_{ij}, (R_{ij}^{2} + X_{ij}^{2}) = v_{j}, \forall ij \in Ω_{l}

{(V_{i}^{\min})}^{2} \leq v_{i} \leq {(V_{i}^{\max})}^{2}, \forall i \in Ω_{b}

P_{c, i}^{\min} \leq p_{i}^{c} \leq P_{c, i}^{\max}, Q_{c, i}^{\min} \leq q_{i}^{c} \leq Q_{c, i}^{\max}, \forall i \in Ω_{c}

P_{g, i}^{\min} \leq p_{i}^{g} \leq P_{g, i}^{\max}, Q_{g, i}^{\min} \leq q_{i}^{g} \leq Q_{g, i}^{\max}, \forall i \in Ω_{g}

P_{ij}^{\min} \leq p_{ij} \leq P_{ij}^{\max}, \forall i j \in Ω_{l}

Q_{ij}^{\min} \leq q_{ij} \leq Q_{ij}^{\max}, \forall i j \in Ω_{l}

0 \leq l_{ij} \leq {(I_{ij}^{\max})}^{2}, \forall i j \in Ω_{l},

[0025]

p_{i}^{c} and q_{i}^{c}

represent an active power and a reactive power output by the thermal power unit respectively,

[0026]

p_{i}^{g} and q_{i}^{g}

[0027]

D_{i}^{p} and D_{i}^{q}

represent an active load and a reactive load of the node i respectively,

[0028]

V_{i}^{\min} and V_{i}^{\max}

represent a lower limit and an upper limit of the voltage amplitude of the node i respectively,

[0029]

P_{c, i}^{\min} and P_{c, i}^{\max}

represent a lower limit and an upper limit of the active power output by a thermal power unit i respectively,

[0030]

Q_{c, i}^{\min} and Q_{c, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the thermal power unit i respectively,

[0031]

P_{g, i}^{\min} and P_{g, i}^{\max}

represent a lower limit and an upper limit of the active power output by a gas power unit i respectively,

[0032]

Q_{g, i}^{\min} and Q_{g, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the gas power unit i respectively,

[0033]

P_{ij}^{\min} and P_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the active power on an electrical transmission line ij respectively,

[0034]

Q_{ij}^{\min} and Q_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the reactive power on the electrical transmission line ij respectively,

[0035]

I_{ij}^{\max}

represents a heat stable current value for the line ij.

[0036] Further, a linear part of the power flow constraints of the gas system is as follows:

[0037]

\sum_{im \in δ_{m}} (e_{im} - τ_{im}) - \sum_{mn \in γ_{m}} e_{mn} - \sum_{mn \in Ξ_{m}} f_{mn} + f_{m}^{s} - f_{m}^{g} = D_{m}^{g}, \forall m \in Ω_{n}

τ_{mn} = K_{mn} e_{mn}, \forall mn \in Ω_{k}

{(H_{mn}^{\min})}^{2} π_{m} \leq π_{n} \leq {(H_{mn}^{\max})}^{2} π_{m}, \forall mn \in Ω_{k}

p_{i}^{g} = T_{m} f_{m}^{g}, \forall i, m \in Ω_{g}

{(\prod_{m}^{\min})}^{2} \leq π_{m} \leq {(\prod_{m}^{\max})}^{2}, \forall m \in Ω_{n}

- F_{mn}^{\max} \leq f_{mn} \leq F_{mn}^{\max}, \forall mn \in Ω_{p}

0 \leq e_{mn} \leq E_{mn}^{\max}, \forall mn \in Ω_{k}

0 \leq f_{m}^{g} \leq F_{g, m}^{\max}, \forall m \in Ω_{g}

F_{s, m}^{\min} \leq f_{m}^{s} \leq F_{s, m}^{\max}, \forall m \in Ω_{s},

[0038]

f_{m}^{s}

represents a nature gas amount injected into the node m by the nature gas source per unit time,

[0039]

f_{m}^{g}

represents a gas consumption amount by a gas generator connected to the node m per unit time, π_mand π_nrepresent squares of gas pressure values at the node m and a node n respectively, a parameter

[0040]

D_{m}^{g}

[0041]

H_{mn}^{\min} and H_{mn}^{\max}

[0042]

\prod_{m}^{\min} and \prod_{m}^{\max}

represent a lower limit and an upper limit of a nodal pressure,

[0043]

F_{mn}^{\min} and F_{mn}^{\max}

represent a lower limit and an upper limit of a gas transmission amount of the gas pipeline mn per unit time,

[0044]

E_{mn}^{\max}

represents an upper limit of a compression rate of the nature gas compressor,

[0045]

F_{g, m}^{\max}

represents an upper limit of a gas consumption rate of the gas generator,

[0046]

F_{s, m}^{\min} and F_{s, m}^{\max}

represent a lower limit and an upper limit of a gas supply amount of the nature gas source per unit time.

[0047] Further, the model is expressed in a matrix form:

[0048]

\min x^{T} Q x + c x + d

s . t . Ax \leq b,

[0049] Further, the convex relaxation formula of the quadratic constraint of the electricity-gas interlinking hybrid energy system is:

[0050]

{(p_{ij})}^{2} + {(q_{ij})}^{2} \leq v_{i} l_{ij}, \forall ij \in Ω_{l} {W_{mn} (f_{mn})}^{2} \leq π_{m} - π_{n}, \forall mn \in Ω_{p},

the first order Taylor expansion is:

[0051]

{(l_{ij} + v_{i})}^{2} - 8 p_{ij}^{⋆} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} - v_{i}^{*})}^{2} - 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l}, \forall ij \in Ω_{l} (π_{m} - π_{n}) - 2 W_{mn} f_{mn}^{*} f_{mn} + {W_{mn} (f_{mn}^{*})}^{2} \leq μ_{mn}^{p}, \forall mn \in Ω_{p}, where p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{mn}^{*}

[0052]

μ_{ij}^{l} and μ_{mn}^{p}

respectively represent an unbalance magnitude of a non-convex constraint of a corresponding line ij and the gas transmission pipeline mn.

[0053] The present disclosure has the following beneficial effects.

[0056] Certainly, the implementation of any product of the present disclosure is not necessary to achieve all the advantages described above at the same time.

Brief Description of the Drawings

[0058] The FIGURE illustrates a schematic diagram of solving the relaxation solution of the model.

Detailed Description of the Embodiments

[0061] In Step 1, a convex optimization part of an optimal energy flow model in the integrated electricity-gas system is established.

[0062] In Step 101, an objective function is:

[0063]

\min \sum_{i \in Ω_{c}} [{C_{i}^{a} (p_{i}^{c})}^{2} + C_{i}^{b} p_{i}^{c} + C_{i}^{c}] + \sum_{i \in Ω_{s}} C_{m}^{f} f_{m}^{s}

where sets Ω_cand Ω_srepresent a thermal power unit node set and a nature gas source node set respectively, a variable

[0064]

p_{i}^{c}

represents an active power of the thermal power unit, a variable

[0065]

f_{m}^{s}

represents a gas supply rate of the nature gas source, parameters

[0066]

C_{i}^{a}, C_{i}^{b} and C_{i}^{c}

represent a second-order term coefficient, a first-order term coefficient and a zero-order term coefficient of the fuel cost in the thermal power unit respectively, and

[0067]

C_{m}^{f}

represents a gas supply cost coefficient of the nature gas source.

[0068] In Step 102, a linear part of the power flow constraints of the electrical power system is as follows:

[0069]

\sum_{{ji}^{\in β_{i}}} (p_{ji} - l_{ji} R_{ji}) - \sum_{ik \in α_{i}} p_{ik} + p_{i}^{c} + p_{i}^{g} = D_{i}^{p}, \forall i \in Ω_{b}

\sum_{ji \in β_{i}} (q_{ji} - l_{ji} X_{ji}) - \sum_{ik \in α_{i}} q_{ik} + q_{i}^{c} + q_{i}^{g} = D_{i}^{q}, \forall i \in Ω_{b}

v_{i} - 2 (p_{ij} R_{ij} + q_{ij} X_{ij}) + l_{ij}, (R_{ij}^{2} + X_{ij}^{2}) = v_{j}, \forall ij \in Ω_{l}

{(V_{i}^{\min})}^{2} \leq v_{i} \leq {(V_{i}^{\max})}^{2}, \forall i \in Ω_{b}

P_{c, i}^{\min} \leq p_{i}^{c} \leq P_{c, i}^{\max}, Q_{c, i}^{\min} \leq q_{i}^{c} \leq Q_{c, i}^{\max}, \forall i \in Ω_{c}

P_{g, i}^{\min} \leq p_{i}^{g} \leq P_{g, i}^{\max}, Q_{g, i}^{\min} \leq q_{i}^{g} \leq Q_{g, i}^{\max}, \forall i \in Ω_{g}

P_{ij}^{\min} \leq p_{ij} \leq P_{ij}^{\max}, \forall i j \in Ω_{l}

Q_{ij}^{\min} \leq q_{ij} \leq Q_{ij}^{\max}, \forall i j \in Ω_{l}

0 \leq l_{ij} \leq {(I_{ij}^{\max})}^{2}, \forall i j \in Ω_{l},

where, sets Ω_b, Ω_land Ω_grepresent an electrical grid node set, an electrical transmission lines set and a gas generator node set respectively, sets α_iand β_irepresent an electrical transmission line set with a node i as a head-terminal node and an electrical transmission lines set with the node i as a tail-terminal node respectively, variables p_ij, p_jiand p_ikrepresent active powers flowing over lines ij, ji and ik, q_ij, q_jiand q_ikrepresent reactive powers flowing over the lines ij, ji and ik, l_ijand l_jirepresent squares of currents flowing over the lines ij and ji,

[0070]

p_{i}^{c} and q_{i}^{c}

represent an active power and a reactive power output by the thermal power unit respectively,

[0071]

p_{i}^{g} and q_{i}^{g}

represent an active power and a reactive power output by the gas generator respectively, v_iand v_jrepresent squares of voltage amplitudes of the node i and a node j, parameters R_ijand R_jirepresent resistances on the lines ij and the line ji respectively, X_ijand X_jirepresent reactances on the line ij and the line ji respectively,

[0072]

D_{i}^{p} and D_{i}^{q}

represent an active load and a reactive load of the node i respectively,

[0073]

V_{i}^{\min} and V_{i}^{\max}

represent a lower limit and an upper limit of the voltage amplitude of the node i respectively,

[0074]

P_{c, i}^{\min} and P_{c, i}^{\max}

represent a lower limit and an upper limit of the active power output by a thermal power unit i respectively,

[0075]

Q_{c, i}^{\min} and Q_{c, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the thermal power unit i respectively,

[0076]

P_{g, i}^{\min} and P_{g, i}^{\max}

represent a lower limit and an upper limit of the active power output by a gas power unit i respectively,

[0077]

Q_{g, i}^{\min} and Q_{g, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the gas power unit i respectively,

[0078]

P_{ij}^{\min} and P_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the active power on an electrical transmission line ij respectively,

[0079]

Q_{ij}^{\min} and Q_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the reactive power on the electrical transmission line ij respectively,

[0080]

I_{ij}^{\max}

represents a heat stable current value for the line ij.

[0081] In Step 103, a linear part of the power flow constraints of the gas system is as follows:

[0082]

\sum_{im \in δ_{m}} (e_{im} - τ_{im}) - \sum_{mn \in γ_{m}} e_{mn} - \sum_{mn \in Ξ_{m}} f_{mn} + f_{m}^{s} - f_{m}^{g} = D_{m}^{g}, \forall m \in Ω_{n}

τ_{mn} = K_{mn} e_{mn}, \forall mn \in Ω_{k}

{(H_{mn}^{\min})}^{2} π_{m} \leq π_{n} \leq {(H_{mn}^{\max})}^{2} π_{m}, \forall mn \in Ω_{k}

p_{i}^{g} = T_{m} f_{m}^{g}, \forall i, m \in Ω_{g}

{(\prod_{m}^{\min})}^{2} \leq π_{m} \leq {(\prod_{m}^{\max})}^{2}, \forall m \in Ω_{n}

- F_{mn}^{\max} \leq f_{mn} \leq F_{mn}^{\max}, \forall mn \in Ω_{p}

0 \leq e_{mn} \leq E_{mn}^{\max}, \forall mn \in Ω_{k}

0 \leq f_{m}^{g} \leq F_{g, m}^{\max}, \forall m \in Ω_{g}

F_{s, m}^{\min} \leq f_{m}^{s} \leq F_{s, m}^{\max}, \forall m \in Ω_{s},

[0083]

f_{m}^{s}

represents a nature gas amount injected into the node m by the nature gas source per unit time,

[0084]

f_{m}^{g}

represents a gas consumption amount by a gas generator connected to the node m per unit time, π_mand π_nrepresent squares of gas pressure values at the node m and a node n respectively, a parameter

[0085]

D_{m}^{g}

[0086]

H_{m n}^{\min} and H_{m n}^{\max}

[0087]

\prod_{m}^{\min} and \prod_{m}^{\max}

represent a lower limit and an upper limit of a nodal pressure,

[0088]

F_{m n}^{\min} and F_{m n}^{\max}

represent a lower limit and an upper limit of a gas transmission amount of the gas pipeline mn per unit time,

[0089]

E_{m n}^{\max}

represents an upper limit of a compression rate of the nature gas compressor,

[0090]

F_{g, m}^{\max}

represents an upper limit of a gas consumption rate of the gas generator,

[0091]

F_{s, m}^{\min} and F_{s, m}^{\max}

represent a lower limit and an upper limit of a gas supply amount of the nature gas source per unit time.

[0092] In Step 104, the above model is expressed in a matrix form as follows:

[0093]

\min x^{T} Q x + c x + d

s . t . Ax \leq b,

[0094] In Step 2, a convex relaxation formula of a quadratic constraint of the optimal energy flow model and its first order Taylor expansion are established.

[0099] Specifically, iterating the expansion formula into a solution model of the integrated electricity-gas hybrid energy system includes as follows.

[0101] When the unbalance magnitude is not greater than a set precision requirement threshold value of the non-convex constraint, the expansion formula is solved according to an energy flow model.

[0104] The actual calculation processes are as follows.

[0105] In Step 201, the convex relaxation formula of the quadratic constraint of the electricity-gas interconnected hybrid energy system is:

[0106]

{(p_{ij})}^{2} + {(q_{ij})}^{2} \leq v_{i} l_{ij}, \forall ij \in Ω_{l} {W_{mn} (f_{mn})}^{2} \leq π_{m} - π_{n}, \forall mn \in Ω_{p} .

[0107] In Step 202, a first order Taylor expansion formula of the above constraint is:

[0108]

{(l_{ij} + v_{i})}^{2} - 8 p_{ij}^{*} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} + v_{i}^{*})}^{2}

- 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l}, \forall ij \in Ω_{l}

(π_{m} - π_{n}) - 2 W_{m n} f_{m n}^{*} f_{m n} + {W_{m n} (f_{m n}^{*})}^{2} \leq μ_{m n}^{p}, \forall m n \in Ω_{p}, where

p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{m n}^{*}

[0109]

μ_{ij}^{l} and μ_{m n}^{p}

respectively represents an unbalance magnitude of a non-convex constraint of a corresponding line ij and the gas transmission pipeline mn.

[0110] In Step 3, the sequential convex programming strategy is used to solve the optimal energy flow model iteratively.

[0111] In Step 301, the precision requirement threshold value ε of the non-convex constraint is set, and the counter is set to be k=0.

[0112] In Step 302, the model is solved

[0113]

\min x^{T} Q x + c x + d

s . t . Ax \leq b

to obtain the current solutions

[0114]

p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{m n}^{*}

of the variables

[0115]

p_{ij}, q_{ij}, l_{ij}, v_{i} and f_{m n} .

[0116] In Step 303, whether the unbalance magnitude of the non-convex constraint is satisfied is determined

[0117]

v_{i}^{*} l_{ij}^{*} - {(p_{ij}^{*})}^{2} - {(q_{ij}^{*})}^{2} \leq ε, \forall i j \in Ω_{l}

π_{m}^{*} - π_{n}^{*} - {W_{m n} (f_{m n}^{*})}^{2} \leq ε, \forall m n \in Ω_{p}

when the accuracy requirement is satisfied, the calculation results are output; if not, Step 304 is executed.

[0118] In Step 304, the model is solved

[0119]

\min x^{T} Q x + c x + d + 2^{k} (\sum_{ij \in Ω_{l}} μ_{ij}^{l} + \sum_{mn \in Ω_{p}} μ_{mn}^{p})

s . t . Ax \leq b

{(l_{ij} + v_{i})}^{2} - 8 p_{ij}^{*} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} - v_{i}^{*})}^{2} - 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l}, \forall ij \in Ω_{l}

(π_{m} - π_{n}) - 2 W_{mn} f_{mn}^{*} f_{mn} + {W_{mn} (f_{mn}^{*})}^{2} \leq μ_{mn}^{p}, \forall mn \in Ω_{p} .

[0120] In Step 305, whether the unbalance magnitude of the non-convex constraint is satisfied is determined

[0121]

v_{i}^{*} l_{ij}^{*} - {(p_{ij}^{*})}^{2} - {(q_{ij}^{*})}^{2} \leq ε, \forall i j \in Ω_{l}

π_{m}^{*} - π_{n}^{*} - {W_{m n} (f_{m n}^{*})}^{2} \leq ε, \forall m n \in Ω_{p}

when accuracy requirement is satisfied, the calculation results are output; if not, the current solutions

[0122]

p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{m n}^{*}

of the variables

[0123]

p_{ij}, q_{ij}, l_{ij}, v_{i} and f_{m n}

are obtained, and Step 304 is executed.

[0124] The following is an example of testing the integrated electricity-gas hybrid energy system, whose system parameters are as shown in Table 1 to Table 9.

[0125]

TABLE 1

Node parameters for an electrical power system
			Lower	Upper		Load
	Active	Reactive	voltage	voltage	Voltage	rejection
Serial	load	load	limit	limit	reference	cost
number	(MW)	(MVAr)	(p.u.)	(p.u.)	(kV)	($/MWh)

0	97.6	44.2	0.94	1.06	345	1000
1	0	0	0.94	1.06	345	1000
2	322	2.4	0.94	1.06	345	1000
3	500	184	0.94	1.06	345	1000
4	0	0	0.94	1.06	345	1000
5	0	0	0.94	1.06	345	1000
6	233.8	84	0.94	1.06	345	1000
7	522	176.6	0.94	1.06	345	1000
8	6.5	−66.6	0.94	1.06	345	1000
9	0	0	0.94	1.06	345	1000
10	0	0	0.94	1.06	345	1000
11	8.53	88	0.94	1.06	345	1000
12	0	0	0.94	1.06	345	1000
13	0	0	0.94	1.06	345	1000
14	320	153	0.94	1.06	345	1000
15	329	32.3	0.94	1.06	345	1000
16	0	0	0.94	1.06	345	1000
17	158	30	0.94	1.06	345	1000
18	0	0	0.94	1.06	345	1000
19	680	103	0.94	1.06	345	1000
20	274	115	0.94	1.06	345	1000
21	0	0	0.94	1.06	345	1000
22	247.5	84.6	0.94	1.06	345	1000
23	308.6	−92.2	0.94	1.06	345	1000
24	224	47.2	0.94	1.06	345	1000
25	139	17	0.94	1.06	345	1000
26	281	75.5	0.94	1.06	345	1000
27	206	27.6	0.94	1.06	345	1000
28	283.5	26.9	0.94	1.06	345	1000
29	0	0	0.94	1.06	345	1000
30	9.2	4.6	0.94	1.06	345	1000
31	0	0	0.94	1.06	345	1000
32	0	0	0.94	1.06	345	1000
33	0	0	0.94	1.06	345	1000
34	0	0	0.94	1.06	345	1000
35	0	0	0.94	1.06	345	1000
36	0	0	0.94	1.06	345	1000
37	0	0	0.94	1.06	345	1000
38	1104	250	0.94	1.06	345	1000

[0126]

TABLE 2

Node parameters for a gas system
		Lower	Upper
Serial	Gas load	pressure	pressure	Load rejection
number	(1000 m3/h)	limit (bar)	limit (bar)	cost ($/1000 m3)

0	0	40	70	1000
1	0	40	70	1000
2	0	40	70	1000
3	100	40	70	1000
4	120	40	60	1000
5	80	40	60	1000
6	0	40	70	1000
7	0	40	70	1000
8	0	40	70	1000
9	0	40	70	1000
10	0	40	70	1000

[0127]

TABLE 3

Line parameters for an electrical power system
					Suscep-	Transmission
Serial	Head	Tail	Resistance	Reactance	tance	power
number	node	node	(p.u.)	(p.u.)	(p.u.)	(MVA)

0	0	1	0.0035	0.0411	0.6987	600
1	0	38	0.001	0.025	0.75	1000
2	1	2	0.0013	0.0151	0.2572	500
3	1	24	0.007	0.0086	0.146	500
4	1	29	0	0.0181	0	900
5	2	3	0.0013	0.0213	0.2214	500
6	2	17	0.0011	0.0133	0.2138	500
7	3	4	0.0008	0.0128	0.1342	600
8	3	13	0.0008	0.0129	0.1382	500
9	4	5	0.0002	0.0026	0.0434	1200
10	4	7	0.0008	0.0112	0.1476	900
11	5	6	0.0006	0.0092	0.113	900
12	5	10	0.0007	0.0082	0.1389	480
13	5	30	0	0.025	0	1800
14	6	7	0.0004	0.0046	0.078	900
15	7	8	0.0023	0.0363	0.3804	900
16	8	38	0.001	0.025	1.2	900
17	9	10	0.0004	0.0043	0.0729	600
18	9	12	0.0004	0.0043	0.0729	600
19	9	31	0	0.02	0	900
20	11	10	0.0016	0.0435	0	500
21	11	12	0.0016	0.0435	0	500
22	12	13	0.0009	0.0101	0.1723	600
23	13	14	0.0018	0.0217	0.366	600
24	14	15	0.0009	0.0094	0.171	600
25	15	16	0.0007	0.0089	0.1342	600
26	15	18	0.0016	0.0195	0.304	600
27	15	20	0.0008	0.0135	0.2548	600
28	15	23	0.0003	0.0059	0.068	600
29	16	17	0.0007	0.0082	0.1319	600
30	16	26	0.0013	0.0173	0.3216	600
31	18	19	0.0007	0.0138	0	900
32	18	32	0.0007	0.0142	0	900
33	19	33	0.0009	0.018	0	900
34	20	21	0.0008	0.014	0.2565	900
35	21	22	0.0006	0.0096	0.1846	600
36	21	34	0	0.0143	0	900
37	22	23	0.0022	0.035	0.361	600
38	22	35	0.0005	0.0272	0	900
39	24	25	0.0032	0.0323	0.531	600
40	24	36	0.0006	0.0232	0	900
41	25	26	0.0014	0.0147	0.2396	600
42	25	27	0.0043	0.0474	0.7802	600
43	25	28	0.0057	0.0625	1.029	600
44	27	28	0.0014	0.0151	0.249	600
45	28	37	0.0008	0.0156	0	1200

[0128]

TABLE 4

Pipeline parameters for a gas system
			Minimum	Maximum	Weymouth
			transmission	transmission	coefficient
Serial	Head	Tail	rate	rate	((bar/(1000
number	node	node	(1000 m3/h)	(1000 m3/h)	m3/h)){circumflex over ( )}2)

0	0	1	−1100	1100	0.058913
1	6	7	−1100	1100	0.058913
2	2	8	−1100	1100	0.061036
3	7	3	−1100	1100	0.066929
4	7	9	−1100	1100	0.077628
5	8	9	−1100	1100	0.061036
6	10	4	−1100	1100	0.063635
7	10	5	−1100	1100	0.071318

[0129]

TABLE 5

Parameters for a gas valve
			Upper limit	Minimum	Maximum
			of pressure	transmission	transmission
Serial	Head	Tail	difference	rate	rate
number	node	node	(bar)	(1000 m3/h)	(1000 m3/h)

0	6	8	120	−1100	1100

[0130]

TABLE 6

Parameters for a nature gas compressor
					Minimum	Maximum
			Minimum	Maximum	transmission	transmission	Gas-
Serial	Head	Tail	compression	compression	rate	rate	driven
number	node	node	ratio	ratio	(1000 m3/h)	(1000 m3/h)	ratio

0	1	6	1.2	2.8	0	1100	0.035
1	9	10	1.2	2.8	0	1100	0.035

[0131]

TABLE 7

Parameters for a thermal power unit
		Minimum	Maximum	Minimum	Maximum
		active	active	reactive	reactive
Serial		power	power	power	power	c2	c1	c0
number	Node	(MW)	(MW)	(MVAr)	(MVAr)	($/MWh2)	($/MWh)	($/h)

0	29	0	1040	140	400	0.01	0.3	0.2
1	30	0	646	−100	300	0.01	0.3	0.2
2	32	0	652	0	250	0.01	0.3	0.2
3	33	0	508	0	167	0.01	0.3	0.2
4	34	0	687	−100	300	0.01	0.3	0.2
5	35	0	580	0	240	0.01	0.3	0.2
6	36	0	564	0	250	0.01	0.3	0.2
7	38	0	1100	−100	300	0.01	0.3	0.2

[0132]

TABLE 8

Parameters for nature gas source
		Minimum gas	Maximum gas
Serial		supply rate	supply rate	Gas supply cost
number	Node	(1000 m3/h)	(1000 m3/h)	($/1000 m3)

0	0	50	750	88.80555
1	1	0	500	88.80555
2	2	100	500	88.80555

[0133]

TABLE 9

Parameters for a gas generator
			Minimum	Maximum	Minimum	Maximum
	Electric		active	active	reactive	reactive	Conversion
Serial	power	Gas	power	power	power	power	ratio
number	node	node	(MW)	(MW)	(MVAr)	(MVAr)	(1000 m3/MW)

0	31	3	0	725	150	300	0.25
1	37	1	0	865	−150	300	0.25

[0135]

TABLE 10

The main technical indexes
Objective function value	Unbalance magnitude	Calculation time

17647.89 RMB	0.06945%	0.194031 seconds

Claims

What is claimed is:

establishing an electrical power flow model of the electrical power system and a gas power flow model of the gas system;

giving a precision requirement threshold value of a non-convex constraint, comparing the threshold value with an unbalance magnitude of the non-convex constraint in the expansion formula;

\min \sum_{i \in Ω_{c}} [{C_{i}^{a} (p_{i}^{c})}^{2} + C_{i}^{b} p_{i}^{c} + C_{i}^{c}] + \sum_{i \in Ω_{s}} C_{m}^{f} f_{m}^{s},

where sets Ω_cand Ω_srepresent a thermal power unit node set and a nature gas source node set respectively, a variable

p_{i}^{c}

represents an active power of the thermal power unit, a variable

f_{m}^{s}

represents a gas supply rate of the nature gas source, parameters

C_{i}^{a}, C_{i}^{b} and C_{i}^{c}

represent a second-order term coefficient, a first-order term coefficient and a zero-order term coefficient of the fuel cost in the thermal power unit respectively, and

C_{m}^{f}

represents a gas supply cost coefficient of the nature gas source.

\sum_{ji \in β_{i}} (p_{ji} - l_{ji} R_{ji}) - \sum_{ik \in α_{i}} p_{ik} + p_{i}^{c} + p_{i}^{g} = D_{i}^{p}, \forall i \in Ω_{b}

\sum_{ji \in β_{i}} (q_{ji} - l_{ji} X_{ji}) - \sum_{ik \in α_{i}} q_{ik} + q_{i}^{c} + q_{i}^{g} = D_{i}^{q}, \forall i \in Ω_{b}

v_{i} - 2 (p_{ij} R_{ij} + q_{ij} X_{ij}) + l_{ij} (R_{ij}^{2} + X_{ij}^{2}) = v_{j}, \forall ij \in Ω_{l}

{(V_{i}^{\min})}^{2} \leq v_{i} \leq {(V_{i}^{\max})}^{2}, \forall i \in Ω_{b}

P_{c, i}^{\min} \leq p_{i}^{c} \leq P_{c, i}^{\max}, Q_{c, i}^{\min} \leq q_{i}^{c} \leq Q_{c, i}^{\max}, \forall i \in Ω_{c}

P_{g, i}^{\min} \leq p_{i}^{g} \leq P_{g, i}^{\max}, Q_{g, i}^{\min} \leq q_{i}^{g} \leq Q_{g, i}^{\max}, \forall i \in Ω_{g}

P_{ij}^{\min} \leq p_{ij} \leq P_{ij}^{\max}, \forall ij \in Ω_{l}

Q_{ij}^{\min} \leq q_{ij} \leq Q_{ij}^{\max}, \forall ij \in Ω_{l}

0 \leq l_{ij} \leq {(I_{ij}^{\max})}^{2}, \forall ij \in Ω_{l},

p_{i}^{c} and q_{i}^{c}

represent an active power and a reactive power output by the thermal power unit respectively,

p_{i}^{g} and q_{i}^{g}

D_{i}^{p} and D_{i}^{q}

represent an active load and a reactive load of the node i respectively,

V_{i}^{\min} and V_{i}^{\max}

represent a lower limit and an upper limit of the voltage amplitude of the node i respectively,

P_{c, i}^{\min} and P_{c, i}^{\max}

represent a lower limit and an upper limit of the active power output by a thermal power unit i respectively,

Q_{c, i}^{\min} and Q_{c, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the thermal power unit i respectively,

P_{g, i}^{\min} and P_{g, i}^{\max}

represent a lower limit and an upper limit of the active power output by a gas power unit i respectively,

Q_{g, i}^{\min} and Q_{g, i}^{\max}

represent a lower limit and an upper limit of the reactive power output by the gas power unit i respectively,

P_{ij}^{\min} and P_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the active power on an electrical transmission line ij respectively,

Q_{ij}^{\min} and Q_{ij}^{\max}

represent a transmission lower limit and a transmission upper limit of the reactive power on the electrical transmission line ij respectively,

I_{ij}^{\max}

represents a heat stable current value for the line ij.

\sum_{im \in δ_{m}} (e_{im} - τ_{im}) - \sum_{mn \in γ_{m}} e_{mn} - \sum_{mn \in Ξ_{m}} f_{mn} + f_{m}^{s} - f_{m}^{g} = D_{m}^{g}, \forall m \in Ω_{n}

τ_{mn} = K_{mn} e_{mn}, \forall mn \in Ω_{k}

{(H_{mn}^{\min})}^{2} π_{m} \leq π_{n} \leq {(H_{mn}^{\max})}^{2} π_{m}, \forall mn \in Ω_{k}

p_{i}^{g} = T_{m} f_{m}^{g}, \forall i, m \in Ω_{g}

{(\prod_{m}^{\min})}^{2} \leq π_{m} \leq {(\prod_{m}^{\max})}^{2}, \forall m \in Ω_{n}

- F_{mn}^{\max} \leq f_{mn} \leq F_{mn}^{\max}, \forall mn \in Ω_{p}

0 \leq e_{mn} \leq E_{mn}^{\max}, \forall mn \in Ω_{k}

0 \leq f_{m}^{g} \leq F_{g, m}^{\max}, \forall m \in Ω_{g}

F_{s, m}^{\min} \leq f_{m}^{s} \leq F_{s, m}^{\max}, \forall m \in Ω_{s},

f_{m}^{s}

represents a nature gas amount injected into the node m by the nature gas source per unit time,

f_{m}^{g}

represents a gas consumption amount by a gas generator connected to the node m per unit time, π_mand π_nrepresent squares of gas pressure values at the node m and a node n respectively, a parameter

D_{m}^{g}

H_{mn}^{\min} and H_{mn}^{\max}

\prod_{m}^{\min} and \prod_{m}^{\max}

represent a lower limit and an upper limit of a nodal pressure, generator,

F_{mn}^{\min} and F_{mn}^{\max}

represent a lower limit and an upper limit of a gas transmission amount of the gas pipeline mn per unit time,

E_{mn}^{\max}

represents an upper limit of a compression rate of the nature gas compressor,

F_{g, m}^{\max}

represents an upper limit of a gas consumption rate of the gas generator,

F_{s, m}^{\min} and F_{s, m}^{\max}

represent a lower limit and an upper limit of a gas supply amount of the nature gas source per unit time.

min x^TQx+cx+d

s.t. Ax≤b,

(p_ij)²+(q_ij)²≤v_il_ij, ∀ij∈Ω_l

W_mn(f_mn)²≤π_m−π_n, ∀mN∈Ω_p,

the first order Taylor expansion is:

{(I_{ij} + v_{i})}^{2} - 8 p_{ij}^{*} p_{ij} - 8 q_{ij}^{*} q_{ij} + 4 {(p_{ij}^{*})}^{2} + 4 {(q_{ij}^{*})}^{2} + {(l_{ij}^{*} - v_{i}^{*})}^{2} - 2 (l_{ij}^{*} - v_{i}^{*}) (l_{ij} - v_{i}) \leq μ_{ij}^{l},

\forall ij \in Ω_{l} (π_{m} - π_{n}) - 2 W_{mn} f_{mn}^{*} f_{mn} + W_{mn} {(f_{mn}^{*})}^{2} \leq μ_{mn}^{p}, ⁠ \forall mn \in Ω_{p},

w here p_{ij}^{*}, q_{ij}^{*}, l_{ij}^{*}, v_{i}^{*} and f_{mn}^{*}

μ_{ij}^{l} and μ_{mn}^{p}

respectively represent an unbalance magnitude of a non-convex constraint of a corresponding line ij and the gas transmission pipeline mn.

Patent

Date

Inventor

Cited By

US2017/0184640(A1)

2017-06-01

Pileggi et al.

Applicant

CN108846507(A)

2018-11-01

Examiner

CN110729766(A)

2020-01-01

Examiner

CN112713615(A)

2021-04-01

Examiner

CN113627021(A)

2021-11-01

Applicant