CN108964041B - Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system - Google Patents

Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system Download PDF

Info

Publication number
CN108964041B
CN108964041B CN201810812846.3A CN201810812846A CN108964041B CN 108964041 B CN108964041 B CN 108964041B CN 201810812846 A CN201810812846 A CN 201810812846A CN 108964041 B CN108964041 B CN 108964041B
Authority
CN
China
Prior art keywords
natural gas
node
flow
constraint
gas
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810812846.3A
Other languages
Chinese (zh)
Other versions
CN108964041A (en
Inventor
卢志刚
顾媛媛
潘丽娜
杨英杰
杨宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yanshan University
Original Assignee
Yanshan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yanshan University filed Critical Yanshan University
Priority to CN201810812846.3A priority Critical patent/CN108964041B/en
Publication of CN108964041A publication Critical patent/CN108964041A/en
Application granted granted Critical
Publication of CN108964041B publication Critical patent/CN108964041B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a control method of the maximum load supply capacity of an electricity-gas interconnection comprehensive energy system, which comprises the following steps: for an electricity-gas interconnection comprehensive energy system comprising an electric power system and a natural gas system, defining the maximum load supply capacity under the operation constraint condition; establishing an electric-gas interconnection comprehensive energy system maximum load supply capacity model considering N-1 safety constraint; establishing a safe operation constraint condition of the electric-gas interconnection comprehensive energy system; linearizing a nonlinear part in the maximum load supply capacity model of the electric-gas interconnection comprehensive energy system; and sequentially inputting parameters of the power system and the natural gas system by using a YALMIP tool box, compiling power system constraints, natural gas system constraints and target function programs, and calling a CPLEX tool box to solve. The method has more accurate result and can be used for optimizing large-scale large systems in actual engineering. The invention is suitable for the failure of the power system or the natural gas system.

Description

Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system
Technical Field
The invention relates to the technical field of planning and operation of interconnected comprehensive energy systems, in particular to a control method for the maximum load supply capacity of an electricity-gas interconnected comprehensive energy system.
Background
In recent years, energy and environmental problems are increasingly highlighted, and the comprehensive energy system becomes a main trend and an important carrier for the development of the future energy system. The comprehensive energy system not only comprises the optimization of new energy and fossil energy, but also comprises the optimization and design of electric power flow, thermal power flow, energy flow, data flow and the like, so that the highest utilization rate of the energy such as electric power, thermal power, natural gas and the like is achieved.
Compared with the traditional primary energy sources such as coal, petroleum and the like, the natural gas has the advantages of cleanness, easiness in storage, safety, reliability and the like. Compared with fossil energy power generation, natural gas power generation also has the advantages of small environmental pollution, obvious energy-saving and carbon-reducing effects, corresponding quick start and stop, low consumption, high profit and the like. As shown in FIG. 1, the gas turbine and the electric gas conversion technology realize the bidirectional flow of the energy flow of the electric-gas interconnected comprehensive energy system. The development of the electricity-to-gas technology converts the surplus output of renewable energy into methane, and then the methane is injected into a natural gas network for transportation or storage, so that the large-scale storage of electric energy becomes possible.
In summary, it is necessary to invent a new control method for maximum load supply capacity for an electrical-gas interconnected integrated energy system, which helps to make system upgrade, modification and extension during planning, and evaluate the margin of the current system operating state during operation, so as to perform better mode arrangement.
Disclosure of Invention
In view of the above technical problems, the present invention provides a method for controlling the maximum load supply capacity of an electrical-pneumatic interconnection energy system.
In order to solve the technical problems, the method of the invention is realized by the following technical scheme:
a control method for the maximum load supply capacity of an electric-gas interconnection comprehensive energy system is characterized by comprising the following steps:
step S1: for an electricity-gas interconnection comprehensive energy system comprising an electric power system and a natural gas system, defining the maximum load supply capacity under the operation constraint condition, wherein the maximum load supply capacity is defined as the maximum load which can be supplied by the electricity-gas interconnection comprehensive energy system under the condition of meeting the operation constraint conditions of the electric power system and the natural gas system, and the operation constraint condition comprises actual operation constraint and N-1 static safety constraint;
step S2: establishing an electric-gas interconnection comprehensive energy system maximum load supply capacity model considering N-1 safety constraint;
step S3: establishing safe operation constraint conditions of the electric-gas interconnected energy system, wherein the safe operation constraint conditions comprise electric power system safe operation constraint, natural gas system safe operation constraint and coupling element constraint;
step S4: linearizing a nonlinear part in the maximum load supply capacity model of the electric-gas interconnection comprehensive energy system;
step S5: using a YALMIP tool box of MATLAB to input parameters of an electric power system and parameters of a natural gas system in sequence; then compiling a linear form power system constraint condition and natural gas system constraint condition program; finally, adding a target function, calling a CPLEX tool box to solve, wherein the flow unit m of the natural gas system is required to be solved when solving3H conversion to power system power unit MW:
Figure GDA0002843514350000021
in the formula, SgasThe pipeline flow rate, measured in MW; ggasIn m for pipe flow velocity3H is measured; LHV is the fixed lower heating value.
In the above technical solution, the model of the maximum load supply capacity of the electrical-electrical interconnection energy system in step S2 is:
TLC＝minTLCk (1)
Figure GDA0002843514350000022
in the formula, TLC is the maximum load supply capacity, and the minimum value of the objective function in the non-fault state and the N-1 state is taken; TLCkThe maximum load supply capacity under the k fault; k is the fault state of the electric-gas interconnected comprehensive energy system, and k is 0,1,2c；ncThe number of accidents is counted; when k is 0, the electric-gas interconnection comprehensive energy system is in a normal state, namely a no-fault state; 1,2, ncIndicating that the interconnected integrated energy systems are respectively positioned at 1 st, 2cA fault condition; g represents the common load of a natural gas system, and g is 1, 2. m is the number of loads of a natural gas system, including a common load and a gas turbine load; z is the number of gas turbine loads;
Figure GDA0002843514350000031
for k failures, TLCkTaking the value of the common load g of the natural gas system at the maximum value; e represents the power system load, e 1, 2. n is the number of loads of the power system;
Figure GDA0002843514350000032
for k failures, TLCkAnd taking the load value of the electric power system load e at the maximum value.
In the above technical solution, the power system safe operation constraint in step S3 considers a dc power flow constraint, and specifically includes the following steps:
step S3101: establishing a power system power balance constraint as follows:
Figure GDA0002843514350000033
in the formula: e is a node-branch incidence matrix of the power system; pl (k)The branch power of the power system is I power under k faults; c is a power system node-generator incidence matrix; pi (k)The output of a generator i of the power system under k faults; d is a power grid node-load incidence matrix;
Figure GDA0002843514350000034
for k failures, TLCkTaking the load value of the electric power system load e at the maximum value;
if f nodes and h branches are set in the power system, E is an f × h matrix, and the structure is as follows:
Figure GDA0002843514350000035
if the p node is at the beginning of the branch l, Epl-1; if the p node is at the end of branch l, E pl1 is ═ 1; when node p has no connection relation with branch l, Epl＝0；
If u generators are set in the power system, C is an f × u matrix, and the structure is as follows:
Figure GDA0002843514350000041
if generator i is connected to node p, then C pi1 is ═ 1; otherwise, Cpi＝0；
Step S3102: calculating the branch load flow of the power system:
Figure GDA0002843514350000042
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000043
the phase angles of nodes p and q at two ends of a branch circuit l of the power system are k times of faults; x is the number ofpqThe reactance of a branch I with nodes at two ends being p and q;
step S3103: establishing upper and lower limit constraints of branch tidal current:
Pl min≤Pl (k)≤Pl max (7)
in the formula, Pl (k)The branch power of the power system is I power under k faults; pl max、Pl minThe upper and lower power limits of the branch circuit I are k times of faults;
step S3104: establishing a balanced node phase angle constraint:
θref＝0 (8)
step S3105: establishing unit output constraint:
Pi min≤Pi (k)≤Pi max (9)
in the formula, Pi (k)The output of a generator i of the power system under k faults; pi max、Pi minThe upper and lower limits of the unit i output when k faults occur;
step S3106: establishing load power constraint:
Figure GDA0002843514350000044
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000045
for k failures, TLCkTaking the load value of the electric power system load e at the maximum value;
Figure GDA0002843514350000046
the upper and lower power limits of the load e of the power system.
In the above technical solution, the constraint on safe operation of the natural gas system in step S3 includes the following steps:
step S3201: and (3) flow balance constraint of a natural gas system:
Figure GDA0002843514350000051
in the formula, a and b are natural gas system nodes;
Figure GDA0002843514350000052
injecting flow into the node a when k faults occur; b e a represents all natural gas nodes connected to node a;
Figure GDA0002843514350000053
the natural gas flow between the nodes a and b under k faults;
Figure GDA0002843514350000054
the natural gas flow through compressor c for k failures; default compressors do not consume natural gas; sgnc(a, b) reflecting the direction of flow through the compressor:
Figure GDA0002843514350000055
in the formula, sgnc(a, b) ═ 1 indicates that node a is the pressurizing station inlet; sgnc(a, b) — 1 indicates that node a is the pressurizing station outlet;
step S3202: pipeline flow equation under natural gas steady state:
Figure GDA0002843514350000056
in the formula, CabConstants related to pipe length, diameter, characteristics, etc.;
Figure GDA0002843514350000057
the node pressures of the natural gas system node a and the node b in k faults are respectively;
Figure GDA0002843514350000058
reflecting the ab flow direction of the k-time fault pipeline, and taking the following values:
Figure GDA0002843514350000059
Figure GDA00028435143500000510
the node pressures of the natural gas system node a and the node b in k faults are respectively;
step S3203: setting the upper and lower limits of the natural gas load:
Figure GDA00028435143500000511
in the formula (I), the compound is shown in the specification,
Figure GDA00028435143500000512
for k failures, TLCkTaking the value of the common load g of the natural gas system at the maximum value;
Figure GDA0002843514350000061
the upper and lower load limits of the common natural gas load g;
step S3204: establishing node pressure constraint:
Figure GDA0002843514350000062
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000063
the node pressure of the natural gas system node a is obtained when k faults occur;
Figure GDA0002843514350000064
the upper and lower pressure limits of the natural gas system node a are set;
step S3205: establishing a natural gas system pipeline flow limit constraint:
Figure GDA0002843514350000065
in the formula, a pipeline w is arranged between nodes a and b of the natural gas system;
Figure GDA0002843514350000066
the flow value of the natural gas pipeline w is the k times of faults;
Figure GDA0002843514350000067
the lower default flow limit is equal to the inverse number of the upper flow limit for the upper w flow limit of the natural gas pipeline;
step S3206: upper and lower gas supply limits:
Figure GDA0002843514350000068
in the formula, Fs (k)The flow is supplied to a natural gas source s when k faults occur;
Figure GDA0002843514350000069
supplying an upper and a lower flow limit for a natural gas source s;
step S3207: establishing compressor constraints including a compressor compression ratio constraint and a compressor flow constraint, wherein the compressor compression ratio constraint is as follows:
Figure GDA00028435143500000610
in the formula (I), the compound is shown in the specification,
Figure GDA00028435143500000611
node pressures of a natural gas system node a and a node b during k faults respectively; tau iscIs the compression ratio of compressor c;
the compressor flow constraint is:
Figure GDA00028435143500000612
in the formula (I), the compound is shown in the specification,
Figure GDA00028435143500000613
the natural gas flow through compressor c for k failures; sgnc(a, b) ═ 1 indicates that node a is the pressurizing station inlet; sgnc(a, b) — 1 indicates that node a is the pressurizing station outlet; cabIs a constant related to the length, diameter, characteristics of the pipe;
Figure GDA0002843514350000071
the node pressures of the natural gas system node a and the node b at k faults are respectively.
In the above solution, the coupling element in step S3 includes a gas turbine, and the constraint of the gas turbine includes:
Figure GDA0002843514350000072
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000073
inputting a heat value for the gas turbine;
Figure GDA0002843514350000074
outputting power for a gas turbine unit i; alpha is alphag,i、βg,i、γg,iDetermined by the gas turbine heat rate curve;
converting heat supplied to the gas turbine unit into a natural gas load:
Figure GDA0002843514350000075
Figure GDA0002843514350000076
equivalent gas load of a natural gas system node a; GHV is a fixed high heating value.
In the above technical solution, step S4 specifically includes the following steps:
step S401: according to the characteristics of the model and the accuracy requirement, determining the appropriate number of the linear stages NPL-1, and segmenting the flow of the natural gas pipeline in the upper and lower limit ranges;
step S402: determining discrete point x of piecewise linearization in independent variable value range1,x2,...,xNPL；
Step S403: calculating f (x) values corresponding to the discrete points, namely calculating square values of the discrete points of the pipeline flow corresponding to the discrete points of the pipeline flow;
step S404: the nonlinear model is represented by linearization according to the following formula:
Figure GDA0002843514350000077
Figure GDA0002843514350000078
δr+1≤ηr,ηr≤δr,r＝1,2,...,NPL-2 (25)
0≤δr≤1,r＝1,2,...,NPL-1 (26)
in the formula, x1,...,xr,xr+1I.e. the flow of the pipeline is above and below the pipelineTaking values of discrete points within a limited range; f (x)1),...,f(xr),f(xr+1) The value is the square value of the corresponding pipeline flow; x, f (x) is the linearized pipeline flow and the pipeline flow squared; deltarThe value range of (a) is 0 to 1, which indicates the position on the r-th segment interval; η is a binary variable.
Due to the adoption of the technical scheme, compared with the prior art, the invention has the following beneficial effects:
1. the coupling of the power system and the natural gas system is considered, the definition of the maximum load supply capacity of the comprehensive energy system is provided, and the method has great significance in planning and operation;
2. considering the N-1 safe operation constraint, the method is suitable for the failure of the power system or the natural gas system;
3. the problem is a complex and nonlinear optimization problem, after linearization processing, the YALMIP toolbox based on MATLAB is used for solving, the result is more accurate than that of the traditional evolutionary algorithm, and the method can be used for the optimization problem of a large-scale large system in actual engineering.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a schematic diagram of a conventional electric-gas interconnected integrated energy system;
FIG. 2 is a block diagram of an integrated energy system utilizing the method of the present invention;
fig. 3 is a flow chart of the method of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention.
FIG. 2 is a block diagram of an integrated energy system utilizing the method of the present invention; in the figure, a 6-node power system is arranged on the left side, and a 7-node natural gas system is arranged on the right side; the left side 1-6 and the right side 1-7 in the figure are respectively the node numbers of the power system and the natural gas system; G1-G5 in the power system are five units, wherein G1 and G2 are gas turbine units; PL1-PL3 are 3 power system loads. In the right natural gas system, GW1 and GW2 are two gas sources of the natural gas system, and GL1-GL3 are common loads of 3 natural gas systems.
Fig. 3 is a flow chart of the method of the present invention. As shown in fig. 2 and 3, the method for controlling the maximum load supply capacity of the electrical-electrical interconnection energy system according to the present invention comprises the following steps:
step S1: for an electricity-gas interconnection comprehensive energy system comprising an electric power system and a natural gas system, defining the maximum load supply capacity under the operation constraint condition, wherein the maximum load supply capacity is defined as the maximum load which can be supplied by the electricity-gas interconnection comprehensive energy system under the condition of meeting the operation constraint conditions of the electric power system and the natural gas system, and the operation constraint condition comprises actual operation constraint and N-1 static safety constraint;
step S2: establishing an electric-gas interconnection comprehensive energy system maximum load supply capacity model considering N-1 safety constraint;
step S3: establishing safe operation constraint conditions of the electric-gas interconnected energy system, wherein the safe operation constraint conditions comprise electric power system safe operation constraint, natural gas system safe operation constraint and coupling element constraint;
step S4: linearizing a nonlinear part in the maximum load supply capacity model of the electric-gas interconnection comprehensive energy system;
step S5: using a YALMIP tool box of MATLAB to input parameters of an electric power system and parameters of a natural gas system in sequence; then compiling the power system constraint conditions and the natural gas system constraint programs in a linear form after arrangement; finally, adding a target function, calling a CPLEX tool box to solve, wherein the power unit of the electric power system is not unified with the flow rate unit of the natural gas system, and solvingThe flow unit m of the natural gas system is required to be adjusted3H conversion to power system power unit MW:
Figure GDA0002843514350000101
in the formula, SgasThe pipeline flow rate, measured in MW; ggasIn m for pipe flow velocity3H is measured; LHV is the fixed lower heating value.
The model of the maximum load supply capacity of the electrical-electrical interconnection comprehensive energy system in the step S2 is as follows:
TLC＝minTLCk (1)
Figure GDA0002843514350000102
in the formula, TLC is the maximum load supply capacity, and the minimum value of the objective function in the non-fault state and the N-1 state is taken; TLCkThe maximum load supply capacity under the k fault; k is the fault state of the electric-gas interconnected comprehensive energy system, and k is 0,1,2c；ncThe number of accidents is counted; when k is 0, the electric-gas interconnection comprehensive energy system is in a normal state, namely a no-fault state; 1,2, ncIndicating that the interconnected integrated energy systems are respectively positioned at 1 st, 2cA fault condition; g represents the common load of a natural gas system, and g is 1, 2. m is the number of loads of a natural gas system, including a common load and a gas turbine load; z is the number of gas turbine loads;
Figure GDA0002843514350000103
for k failures, TLCkTaking the value of the common load g of the natural gas system at the maximum value; e represents the power system load, e 1, 2. n is the number of loads of the power system;
Figure GDA0002843514350000104
for k failures, TLCkAnd taking the load value of the electric power system load e at the maximum value.
The power system safe operation constraint in the step S3 considers a direct current power flow constraint, and specifically includes the following steps:
step S3101: establishing a power system power balance constraint as follows:
Figure GDA0002843514350000105
in the formula: e is a node-branch incidence matrix of the power system; pl (k)The branch power of the power system is I power under k faults; c is a power system node-generator incidence matrix; pi (k)The output of a generator i of the power system under k faults; d is a power grid node-load incidence matrix;
Figure GDA0002843514350000111
for k failures, TLCkTaking the load value of the electric power system load e at the maximum value;
if f nodes and h branches are set in the power system, E is an f × h matrix, and the structure is as follows:
Figure GDA0002843514350000112
if the p node is at the beginning of the branch l, Epl-1; if the p node is at the end of branch l, E pl1 is ═ 1; when node p has no connection relation with branch l, Epl＝0；
If u generators are set in the power system, C is an f × u matrix, and the structure is as follows:
Figure GDA0002843514350000113
if generator i is connected to node p, then C pi1 is ═ 1; otherwise, Cpi＝0；
Step S3102: calculating the branch load flow of the power system:
Figure GDA0002843514350000114
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000115
the phase angles of nodes p and q at two ends of a branch circuit l of the power system are k times of faults; x is the number ofpqThe reactance of a branch circuit l with nodes p and q at two ends of the power system is shown;
step S3103: establishing upper and lower limit constraints of branch tidal current:
Pl min≤Pl (k)≤Pl max (7)
in the formula, Pl (k)The branch power of the power system is I power under k faults; pl max、Pl minThe upper and lower power limits of the branch circuit I are k times of faults;
step S3104: establishing a balanced node phase angle constraint:
θref＝0 (8)
step S3105: establishing unit output constraint:
Pi min≤Pi (k)≤Pi max (9)
in the formula, Pi (k)The output of a generator i of the power system under k faults; pi max、Pi minThe upper and lower limits of the unit i output when k faults occur;
step S3106: establishing load power constraint:
Figure GDA0002843514350000121
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000122
for k failures, TLCkTaking the load value of the electric power system load e at the maximum value, TLCkIs the maximum negative at the k-th faultA load supply capacity;
Figure GDA0002843514350000123
the upper and lower power limits of the load e of the power system.
The safe operation constraint of the natural gas system in the step S3 includes the following steps:
step S3201: and (3) flow balance constraint of a natural gas system:
Figure GDA0002843514350000124
in the formula, a and b are natural gas system nodes;
Figure GDA0002843514350000125
injecting flow into the node a when k faults occur; b e a represents all natural gas nodes connected to node a;
Figure GDA0002843514350000126
the natural gas flow between the nodes a and b under k faults;
Figure GDA0002843514350000127
the natural gas flow through compressor c for k failures; default compressors do not consume natural gas; sgnc(a, b) reflecting the direction of flow through the compressor:
Figure GDA0002843514350000128
in the formula, sgnc(a, b) ═ 1 indicates that node a is the pressurizing station inlet; sgncAnd (a, b) — 1 indicates that the node a is a pressurizing station outlet.
Step S3202: pipeline flow equation under natural gas steady state:
Figure GDA0002843514350000129
in the formula, CabIs related to the length and diameter of the pipelineConstants associated with the characteristics, etc.;
Figure GDA0002843514350000131
the node pressures of the natural gas system node a and the node b in k faults are respectively;
Figure GDA0002843514350000132
reflecting the ab flow direction of the k-time fault pipeline, and taking the following values:
Figure GDA0002843514350000133
step S3203: setting the upper and lower limits of the natural gas load:
Figure GDA0002843514350000134
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000135
for k failures, TLCkTaking the value of the common load g of the natural gas system at the maximum value;
Figure GDA0002843514350000136
the upper and lower load limits of the common natural gas load g;
step S3204: establishing node pressure constraint:
Figure GDA0002843514350000137
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000138
the node pressure of the natural gas system node a is obtained when k faults occur;
Figure GDA0002843514350000139
the upper and lower pressure limits of the natural gas system node a are set;
step S3205: establishing a natural gas system pipeline flow limit constraint:
Figure GDA00028435143500001310
in the formula, a pipeline w is arranged between nodes a and b of the natural gas system;
Figure GDA00028435143500001311
the flow value of the natural gas pipeline w is the k times of faults;
Figure GDA00028435143500001312
the lower default flow limit is equal to the inverse number of the upper flow limit for the upper w flow limit of the natural gas pipeline;
step S3206: upper and lower gas supply limits:
Figure GDA00028435143500001313
in the formula, Fs (k)The flow is supplied to a natural gas source s when k faults occur;
Figure GDA00028435143500001314
the natural gas source s is supplied with upper and lower flow limits.
Step S3207: establishing compressor constraints including a compressor compression ratio constraint and a compressor flow constraint, wherein the compressor compression ratio constraint is as follows:
Figure GDA0002843514350000141
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000142
node pressures of a natural gas system node a and a node b during k faults respectively; tau iscIs the compression ratio of compressor c;
the compressor flow constraint is:
Figure GDA0002843514350000143
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000144
the natural gas flow through compressor c for k failures; sgnc(a, b) ═ 1 indicates that node a is the pressurizing station inlet; sgnc(a, b) — 1 indicates that node a is the pressurizing station outlet; cabIs a constant related to the length, diameter, characteristics of the pipe;
Figure GDA0002843514350000145
the node pressures of the natural gas system node a and the node b at k faults are respectively.
The coupling elements of the natural gas system and the electric power system mainly comprise a gas turbine, an electric gas conversion device and an energy hub, and the coupling constraint considering the gas turbine comprises the following steps:
Figure GDA0002843514350000146
in the formula (I), the compound is shown in the specification,
Figure GDA0002843514350000147
inputting a heat value for the gas turbine;
Figure GDA0002843514350000148
outputting power for a gas turbine unit i; alpha is alphag,i、βg,i、γg,iDetermined by the gas turbine heat rate curve;
converting heat supplied to the gas turbine unit into a natural gas load:
Figure GDA0002843514350000149
Figure GDA00028435143500001410
natural gas systemEquivalent gas load of the node a; GHV is a fixed high heating value.
Step S4 specifically includes the following steps:
step S401: according to the characteristics of the model and the accuracy requirement, determining the appropriate number of the linear stages NPL-1, namely segmenting the flow of the natural gas pipeline in the upper and lower limit ranges;
step S402: determining discrete point x of piecewise linearization in independent variable value range1,x2,...,xNPL；
Step S403: calculating f (x) values corresponding to the discrete points, namely calculating the square value of the discrete points of the pipeline flow corresponding to the discrete points of the pipeline flow in the invention;
step S404: the nonlinear model is represented by linearization according to the following formula:
Figure GDA0002843514350000151
Figure GDA0002843514350000152
δr+1≤ηr,ηr≤δr,r＝1,2,...,NPL-2 (25)
0≤δr≤1,r＝1,2,...,NPL-1 (26)
in the formula, x1,...,xr,xr+1Namely, the value of the discrete point of the pipeline flow in the upper and lower limit ranges is obtained; f (x)1),...,f(xr),f(xr+1) The value is the square value of the corresponding pipeline flow; x, f (x) is the linearized pipeline flow and the pipeline flow squared; deltarThe value range of (a) is 0 to 1, which indicates the position on the r-th segment interval; η is a binary variable.
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes or modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.

Claims (6)

1. A control method for the maximum load supply capacity of an electric-gas interconnection comprehensive energy system is characterized by comprising the following steps:
step S1: for an electricity-gas interconnection comprehensive energy system comprising an electric power system and a natural gas system, defining the maximum load supply capacity under the operation constraint condition, wherein the maximum load supply capacity is defined as the maximum load which can be supplied by the electricity-gas interconnection comprehensive energy system under the condition of meeting the operation constraint conditions of the electric power system and the natural gas system, and the operation constraint condition comprises actual operation constraint and N-1 static safety constraint;
step S2: establishing an electric-gas interconnection comprehensive energy system maximum load supply capacity model considering N-1 safety constraint;
step S3: establishing safe operation constraint conditions of the electric-gas interconnected energy system, wherein the safe operation constraint conditions comprise electric power system safe operation constraint, natural gas system safe operation constraint and coupling element constraint;
step S4: linearizing a nonlinear part in the maximum load supply capacity model of the electric-gas interconnection comprehensive energy system;
step S5: using a YALMIP tool box of MATLAB to input parameters of an electric power system and parameters of a natural gas system in sequence; then compiling a linear form power system constraint condition and natural gas system constraint condition program; finally, adding a target function, calling a CPLEX tool box to solve, wherein the flow unit m of the natural gas system is required to be solved when solving3H conversion to power system power unit MW:
Figure FDA0002843514340000011
in the formula, SgasThe pipeline flow rate, measured in MW; ggasIn m for pipe flow velocity3H is measured; LHV is the fixed lower heating value.
2. The method for controlling the maximum load supply capacity of the electric-pneumatic interconnection energy system according to claim 1, wherein the model of the maximum load supply capacity of the electric-pneumatic interconnection energy system in step S2 is:
TLC＝min TLCk (1)
Figure FDA0002843514340000021
in the formula, TLC is the maximum load supply capacity, and the minimum value of the objective function in the non-fault state and the N-1 state is taken; TLCkThe maximum load supply capacity under the k fault; k is 0,1,2c；ncThe number of accidents is counted; when k is 0, the electric-gas interconnection comprehensive energy system is in a normal state, namely a no-fault state; 1,2, ncIndicating that the interconnected integrated energy systems are respectively positioned at 1 st, 2cA fault condition; g represents the common load of a natural gas system, and g is 1, 2. m is the number of loads of a natural gas system, including a common load and a gas turbine load; z is the number of gas turbine loads;
Figure FDA0002843514340000022
for k failures, TLCkTaking the value of the common load g of the natural gas system at the maximum value; e represents the power system load, e 1, 2. n is the number of loads of the power system;
Figure FDA0002843514340000023
for k failures, TLCkAnd taking the load value of the electric power system load e at the maximum value.
3. The method for controlling the maximum load supply capacity of the electric-gas interconnected integrated energy system according to claim 1, wherein the safety operation constraint of the power system in the step S3 considers a direct current power flow constraint, and specifically comprises the following steps:
step S3101: establishing a power system power balance constraint as follows:
Figure FDA0002843514340000024
in the formula: e is a node-branch incidence matrix of the power system; pl (k)The branch power of the power system is I power under k faults; c is a power system node-generator incidence matrix; pi (k)The output of a generator i of the power system under k faults; d is a power grid node-load incidence matrix;
Figure FDA0002843514340000025
e power of the load of the power system under k faults;
if f nodes and h branches are set in the power system, E is an f × h matrix, and the structure is as follows:
Figure FDA0002843514340000031
if the p node is at the beginning of the branch l, Epl-1; if the p node is at the end of branch l, Epl1 is ═ 1; when node p has no connection relation with branch l, Epl＝0；
If u generators are set in the power system, C is an f × u matrix, and the structure is as follows:
Figure FDA0002843514340000032
if generator i is connected to node p, then Cpi1 is ═ 1; otherwise, Cpi＝0；
Step S3102: calculating the branch load flow of the power system:
Figure FDA0002843514340000033
in the formula (I), the compound is shown in the specification,
Figure FDA0002843514340000034
the phase angles of nodes p and q at two ends of a branch circuit l of the power system are k times of faults; x is the number ofpqThe reactance of a branch I with nodes at two ends being p and q;
step S3103: establishing upper and lower limit constraints of branch tidal current:
Pl min≤Pl (k)≤Pl max (7)
in the formula, Pl (k)The branch power of the power system is I power under k faults; pl max、Pl minThe upper and lower power limits of the branch circuit I are k times of faults;
step S3104: establishing a balanced node phase angle constraint:
θref＝0 (8)
step S3105: establishing unit output constraint:
Pi min≤Pi (k)≤Pi max (9)
in the formula, Pi (k)The output of a generator i of the power system under k faults; pi max、Pi minThe upper and lower limits of the unit i output when k faults occur;
step S3106: establishing load power constraint:
Figure FDA0002843514340000041
in the formula (I), the compound is shown in the specification,
Figure FDA0002843514340000042
for k failures, TLCkTaking the load value of the electric power system load e at the maximum value, TLCkThe maximum load supply capacity under the k fault;
Figure FDA0002843514340000043
the upper and lower power limits of the load e of the power system.
4. The method for controlling the maximum load supply capacity of an electric-gas interconnected integrated energy system as claimed in claim 1, wherein the safety operation constraint of the natural gas system in the step S3 includes the following steps:
step S3201: and (3) flow balance constraint of a natural gas system:
Figure FDA0002843514340000044
in the formula, a and b are natural gas system nodes;
Figure FDA0002843514340000045
injecting flow into the node a when k faults occur; b e a represents all natural gas nodes connected to node a;
Figure FDA0002843514340000046
the natural gas flow between the nodes a and b under k faults;
Figure FDA0002843514340000047
the natural gas flow through compressor c for k failures; default compressors do not consume natural gas; sgnc(a, b) reflecting the direction of flow through the compressor:
Figure FDA0002843514340000048
in the formula, sgnc(a, b) ═ 1 indicates that node a is the pressurizing station inlet; sgnc(a, b) — 1 indicates that node a is the pressurizing station outlet;
step S3202: pipeline flow equation under natural gas steady state:
Figure FDA0002843514340000049
in the formula, CabIs a constant related to the length, diameter, characteristics of the pipe;
Figure FDA00028435143400000410
the node pressures of the natural gas system node a and the node b in k faults are respectively;
Figure FDA00028435143400000411
reflecting the ab flow direction of the k-time fault pipeline, and taking the following values:
Figure FDA0002843514340000051
Figure FDA0002843514340000052
the node pressures of the natural gas system node a and the node b in k faults are respectively;
step S3203: setting the upper and lower limits of the natural gas load:
Figure FDA0002843514340000053
in the formula (I), the compound is shown in the specification,
Figure FDA0002843514340000054
for k failures, TLCkTaking the value of the common load g of the natural gas system at the maximum value;
Figure FDA0002843514340000055
the upper and lower load limits of the common natural gas load g;
step S3204: establishing node pressure constraint:
Figure FDA0002843514340000056
in the formula (I), the compound is shown in the specification,
Figure FDA0002843514340000057
the node pressure of the natural gas system node a is obtained when k faults occur;
Figure FDA0002843514340000058
the upper and lower pressure limits of the natural gas system node a are set;
step S3205: establishing a natural gas system pipeline flow limit constraint:
Figure FDA0002843514340000059
in the formula, a pipeline w is arranged between nodes a and b of the natural gas system;
Figure FDA00028435143400000510
the flow value of the natural gas pipeline w is the k times of faults;
Figure FDA00028435143400000511
the lower default flow limit is equal to the inverse number of the upper flow limit for the upper w flow limit of the natural gas pipeline;
step S3206: upper and lower gas supply limits:
Figure FDA00028435143400000512
in the formula (I), the compound is shown in the specification,
Figure FDA00028435143400000513
the flow is supplied to a natural gas source s when k faults occur;
Figure FDA00028435143400000514
supplying an upper and a lower flow limit for a natural gas source s;
step S3207: establishing compressor constraints including a compressor compression ratio constraint and a compressor flow constraint, wherein the compressor compression ratio constraint is as follows:
Figure FDA0002843514340000061
in the formula (I), the compound is shown in the specification,
Figure FDA0002843514340000062
node pressures of a natural gas system node a and a node b during k faults respectively; tau iscIs the compression ratio of compressor c;
the compressor flow constraint is:
Figure FDA0002843514340000063
in the formula (I), the compound is shown in the specification,
Figure FDA0002843514340000064
the natural gas flow through compressor c for k failures; sgnc(a, b) ═ 1 indicates that node a is the pressurizing station inlet; sgnc(a, b) — 1 indicates that node a is the pressurizing station outlet; cabIs a constant related to the length, diameter, characteristics of the pipe;
Figure FDA0002843514340000065
the node pressures of the natural gas system node a and the node b at k faults are respectively.
5. The method for controlling maximum load supplying capacity of an electric-pneumatic interconnected integrated energy system as claimed in claim 1, wherein the coupling element in step S3 includes a gas turbine, and the constraint of the gas turbine includes:
Figure FDA0002843514340000066
in the formula (I), the compound is shown in the specification,
Figure FDA0002843514340000067
inputting a heat value for the gas turbine;
Figure FDA0002843514340000068
outputting power for a gas turbine unit i; alpha is alphag,i、βg,i、γg,iDetermined by the gas turbine heat rate curve;
converting heat supplied to the gas turbine unit into a natural gas load:
Figure FDA0002843514340000069
Figure FDA00028435143400000610
equivalent gas load of a natural gas system node a; GHV is a fixed high heating value.
6. The method for controlling the maximum load supply capacity of the electric-gas interconnected integrated energy system according to claim 1, wherein the step S4 specifically comprises the following steps:
step S401: according to the characteristics of the model and the accuracy requirement, determining the appropriate number of the linear stages NPL-1, and segmenting the flow of the natural gas pipeline in the upper and lower limit ranges;
step S402: determining discrete point x of piecewise linearization in independent variable value range1,x2,...,xNPL；
Step S403: calculating f (x) values corresponding to the discrete points, namely calculating square values of the discrete points of the pipeline flow corresponding to the discrete points of the pipeline flow;
step S404: the nonlinear model is represented by linearization according to the following formula:
Figure FDA0002843514340000071
Figure FDA0002843514340000072
δr+1≤ηr,ηr≤δr,r＝1,2,...,NPL-2 (25)
0≤δr≤1,r＝1,2,...,NPL-1 (26)
in the formula, x1,...,xr,xr+1Namely, the value of the discrete point of the pipeline flow in the upper and lower limit ranges is obtained; f (x)1),...,f(xr),f(xr+1) The value is the square value of the corresponding pipeline flow; x, f (x) is the linearized pipeline flow and the pipeline flow squared; deltarThe value range of (a) is 0 to 1, which indicates the position on the r-th segment interval; η is a binary variable.
CN201810812846.3A 2018-07-23 2018-07-23 Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system Active CN108964041B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810812846.3A CN108964041B (en) 2018-07-23 2018-07-23 Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810812846.3A CN108964041B (en) 2018-07-23 2018-07-23 Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system

Publications (2)

Publication Number Publication Date
CN108964041A CN108964041A (en) 2018-12-07
CN108964041B true CN108964041B (en) 2021-02-23

Family

ID=64464602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810812846.3A Active CN108964041B (en) 2018-07-23 2018-07-23 Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system

Country Status (1)

Country Link
CN (1) CN108964041B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109687456B (en) * 2019-01-25 2020-07-03 燕山大学 Scheduling method and system of electric power natural gas system
CN109919462B (en) * 2019-02-22 2023-05-12 华南理工大学 Safety constraint optimal energy flow calculation method for electric-gas hybrid system
CN110222970B (en) * 2019-05-30 2022-12-06 天津大学 Elastic scheduling method of gas-electricity coupling comprehensive energy system considering energy storage reserve
CN110504714B (en) * 2019-09-11 2021-01-26 中国南方电网有限责任公司 Control method, device, system, equipment and storage medium of power distribution system
CN111768036B (en) * 2020-06-29 2023-11-03 国网上海市电力公司 Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid
CN112541242B (en) * 2020-12-23 2024-01-12 燕山大学 Expansion planning method for elastically-lifting urban electricity and gas comprehensive energy system

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107808216A (en) * 2017-10-24 2018-03-16 重庆大学 Electrical heat interacted system abandons wind and abandons light and the comprehensive minimum Optimized model construction method of electric thermic load reduction

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104734155B (en) * 2015-04-14 2017-03-01 河海大学 A kind of method obtaining electrical interconnection energy resource system available transmission capacity
JP6402075B2 (en) * 2015-07-13 2018-10-10 株式会社日立製作所 Coal gasification combined power plant
CN105356447B (en) * 2015-10-27 2018-01-19 河海大学 A kind of electrical interconnection integrated energy system Analysis of Steady-state Security Region method
CN205753443U (en) * 2016-05-11 2016-11-30 大连派思新能源发展有限公司 A kind of distributed busbar protection being applied to LNG liquefaction plant
CN206160341U (en) * 2016-10-12 2017-05-10 北京燃气能源发展有限公司 Gas combustion engine drive ice storage air conditioning system
CN106950840B (en) * 2017-05-11 2020-05-19 山东理工大学 Power grid peak clipping-oriented hierarchical distributed coordination control method for comprehensive energy system
CN107257132B (en) * 2017-05-17 2020-04-21 重庆大学 Method for constructing comprehensive load reduction model of electricity-gas interconnection system considering wind power abandonment
CN107451680A (en) * 2017-07-12 2017-12-08 河海大学 A kind of electrical interconnection system random optimum trend heuristic calculation methods

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107808216A (en) * 2017-10-24 2018-03-16 重庆大学 Electrical heat interacted system abandons wind and abandons light and the comprehensive minimum Optimized model construction method of electric thermic load reduction

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A Line Flow Granular Computing Approach for Economic Dispatch With Line Constraints;Xueping Li et al;《IEEE Transactions on Power Systems》;20171130;第32卷(第6期);4832-4842 *
基于供电能力的主动配电系统安全等级划分;卢志刚等;《中国电机工程学报》;20170505;第37卷(第9期);2539-2550 *

Also Published As

Publication number Publication date
CN108964041A (en) 2018-12-07

Similar Documents

Publication Publication Date Title
CN108964041B (en) Control method for maximum load supply capacity of electricity-gas interconnection comprehensive energy system
CN107808216B (en) Method for constructing comprehensive minimum optimization model for wind curtailment and electrical thermal load reduction of electricity-gas-heat interconnection system
CN108173282A (en) A kind of consideration electricity turns gas operating cost integrated energy system Optimization Scheduling
CN109919399B (en) Day-ahead economic dispatching method and system for comprehensive energy system
CN111799777A (en) Comprehensive energy planning method considering coupling of natural gas and electric power
CN112016195B (en) Flexible planning method of electrical coupling system considering energy supply reliability
Camblong et al. Design of a SOFC/GT/SCs hybrid power system to supply a rural isolated microgrid
CN112671046B (en) Coordination optimization configuration method and system for wind, light and fire storage delivery capacity
CN110601203A (en) Piecewise linearization optimal power flow calculation method for electric-gas coupling system
CN110729766A (en) Distributed robust optimization method for multi-region electricity-gas integrated energy system
CN112736923A (en) Natural gas network cascading failure evaluation control method considering power grid failure influence
CN112332461A (en) Multi-energy micro-grid group control method based on two-stage robust model
CN115017723A (en) Power and natural gas system combined optimization load flow calculation method considering hydrogen energy injection
Yang et al. Collaborative planning of integrated natural gas and power supply system considering P2G technique
CN115099063A (en) Operation optimization method for electricity-mixed hydrogen natural gas coupling comprehensive energy system
CN107482664B (en) Island alternating current grid system strength evaluation method and system
CN110829484B (en) Space-time decomposition-based global energy interconnection power balance optimization method
CN110380447B (en) Risk reduction scheduling method for electricity-gas interconnection energy system under failure of fan
Bakken et al. Technical and economic aspects of using a long submarine HVDC connection for frequency control
CN109492809B (en) Wind power plant value evaluation method based on node electricity price
CN113221358B (en) Standby output optimization method of electric-gas coupling system based on reliability parameters
Majidi et al. Recapturing wasted energy in water pressure reducing valves via in-conduit hydropower generators
CN110704935B (en) Automatic debugging and outputting method of SPS software
Zhou et al. Optimal Operation of Integrated Energy Systems Based on Multi-Energy Coupling With Power-to-Gas Technology
CN108985579B (en) Power supply configuration planning method and system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant