CN109242366A - A kind of multi-period tide optimization method of electric-gas interconnection integrated energy system - Google Patents
A kind of multi-period tide optimization method of electric-gas interconnection integrated energy system Download PDFInfo
- Publication number
- CN109242366A CN109242366A CN201811310889.8A CN201811310889A CN109242366A CN 109242366 A CN109242366 A CN 109242366A CN 201811310889 A CN201811310889 A CN 201811310889A CN 109242366 A CN109242366 A CN 109242366A
- Authority
- CN
- China
- Prior art keywords
- gas
- electric
- integrated energy
- energy system
- period
- 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.)
- Granted
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 28
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 108
- 239000007789 gas Substances 0.000 claims abstract description 79
- 239000003345 natural gas Substances 0.000 claims abstract description 54
- 230000002708 enhancing effect Effects 0.000 claims abstract description 13
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 230000005611 electricity Effects 0.000 claims description 7
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 230000009194 climbing Effects 0.000 claims description 3
- 239000000567 combustion gas Substances 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000000750 progressive effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of multi-period tide optimization methods of electric-gas interconnection integrated energy system, this method comprises: (1), which obtains electric-gas, interconnects integrated energy system information；(2) electric-gas is constructed according to system information and interconnects the multi-period scheduling model of integrated energy system；(3) the equations turned natural gas tide model for enhancing second order cone constraint type of the Nonlinear Nonconvex for the electric-gas being interconnected into the natural gas line flow in the multi-period scheduling model of integrated energy system and pressure；(4) the electric-gas interconnection multi-period scheduling model of integrated energy system after conversion is solved to obtain optimal solution；(5) using the optimal solution as initial value, and linearizing iteration is carried out to the multi-period scheduling model of electric-gas interconnection integrated energy system after conversion using DCP method, until natural gas system it is stringent meet trend constraint；(6) it is exported the last solution at the end of iteration as the optimal load flow solution in future time period.Multi-period tide optimization can be effectively performed in the present invention.
Description
Technical field
The present invention relates to electrical interconnection technology fields more particularly to a kind of electric-gas to interconnect the multi-period of integrated energy system
Tide optimization method.
Background technique
Turn gas technology in the promotion of Generation Side specific gravity and electricity in view of gas turbine to apply in the power system, makes power train
Bidirectional energy stream between uniting natural gas system is possibly realized, the development of natural gas so that electric system with natural gas system by phase
Mutually independently it is changed into intercouple (gradually development is close coupling).It is therefore desirable to break between existing energy resource system independently planning,
The mode of operation, and construct the integrated energy system of heterogeneous energy source interconnections unify, a variety of.Furthermore, energy internet
It can be regarded as on the basis of polymorphic type energy source interconnection (i.e. integrated energy system), the depth of internet thinking and technology incorporates,
Thus the building of integrated energy system will also become the important link of China's energy Internet Strategy.Compared to existing energy system
The advantage of system, electric-gas interconnection integrated energy system is: 1) higher efficiency of energy utilization, bigger economic interests；2) promote
Into renewable energy scale exploitation with it is grid-connected；3) flexibility between increase system and the energy are complementary.
There are slow motion step responses for natural gas system, thus in short-term time scale scheduling, need to consider natural gas system pipe
Line-pack storage characteristics in road.Meanwhile natural gas system tide model is essentially Nonlinear Nonconvex equation, for non-convex
For optimization, locally optimal solution can only be often obtained, and the convergence solved is vulnerable to initial value affecting.The running optimizatin of electric system
Non-convex problem is equally faced, but the linear tide model of direct current can substitute the non-linear trend mould of exchange in engineering practice
Type, therefore efficient linear tide model is retrievable in electric power system optimization.However for natural gas system, existing trend
Model linearization model uses subsection-linear method, need to introduce a large amount of integer variables, greatly increase computation complexity, if considering
A small amount of segmentation integer variable, piecewise linearity precision are difficult to meet engineering practice requirement.The thus day of efficient convex optimization form
Right gas system load flow is particularly important.
Summary of the invention
Goal of the invention: in view of the problems of the existing technology the present invention, provides a kind of electric-gas interconnection integrated energy system
Multi-period tide optimization method guarantees the high efficiency and optimality that understand using second order cone optimization in method, using DCP
(difference-of-convex) method guarantees that the feasibility understood (can meet the stringent physics of natural gas system about
Beam).
Technical solution: the multi-period tide optimization method of electric-gas of the present invention interconnection integrated energy system includes:
(1) power system information and natural gas system information of electric-gas interconnection integrated energy system are obtained respectively；
(2) according to the power system information and natural gas system information, when building electric-gas interconnection integrated energy system is more
Section scheduling model；
(3) electric-gas is interconnected to the natural gas line flow and pressure in the multi-period scheduling model of integrated energy system
Nonlinear Nonconvex it is equations turned for enhancing second order cone constraint type natural gas tide model；
(4) the electric-gas interconnection multi-period scheduling model of integrated energy system after conversion is solved to obtain optimal solution；
(5) using the optimal solution as initial value, and integrated energy system is interconnected to the electric-gas after conversion using DCP method
Multi-period scheduling model carries out linearizing iteration, until natural gas system it is stringent meet trend constraint；
(6) it is exported the last solution at the end of iteration as the optimal load flow solution in future time period.
Further, the power system information obtained in step (1) are as follows: power network topology, branch parameters information, generator ginseng
Count information, the electric load information in future time period, the prediction value information of wind-powered electricity generation；The parameter information of natural gas system are as follows: natural gas
Net topology, pipe parameter information, when the line-pack amount of storage of preceding pipeline, the parameter information of gas source, the gas in future time period is negative
Lotus information.
Further, the electric-gas established in step (2) interconnects the multi-period scheduling model of integrated energy system specifically:
In formula, subscript 0 indicates that benchmark Run-time scenario, subscript t indicate t moment, and i, j, m, n indicate the section in energy resource system
Point；Subscript max indicates upper limit value, and subscript min indicates lower limit value；f0For optimization object function, NGFor generator collection, NgFor combustion
Gas-turbine set, NsFor gas source set, NWFor wind power plant set, T0For when discontinuity surface number, CG,iFor generator cost coefficient, CS,m
For gas source cost coefficient, CW,iFor abandonment cost coefficient,For abandonment percentage,For generator output,
For generator output lower and upper limit,For abandonment ratio,It is expected to contribute for wind-powered electricity generation, PL,i,tFor burden with power,For line
Road i-j active power, EN (i) are and node i connected node set, bijFor route i-j susceptance, θ is node phase angle vector,For route i-j active power lower and upper limit；For the amount of natural gas of gas turbine consumption, η is combustion gas wheel
Unit transformation efficiency,For gas source power output, FD,m,tFor natural gas load, GC (m), GP (m), GN (m) are respectively that node m connects
Pressurizing point, gas turbine and the pipeline set connect,For the absorption flow of pressurizing point k,For the flow for flowing through pressurizing point k；AndRespectively pipeline m-n head end, end and average flow rate, CmnFor pipeline m-n pressure drop constant, πmWith
πnRespectively node m, n pressure,Respectively node m low pressure limit and the upper limit；GLmnFor the line- of pipeline m-n
Pack gas-storing capacity, KmnFor the line-pack parameter of pipeline m-n；For gas driven pressurizing point dissipative coefficient,For power generation
The maximum active climbing of machine,Respectively pressurizing point first and end pressure,WithFor pressurizing point step-up ratio
Upper and lower bound,For gas source contribute lower and upper limit,Climb for gas source maximum,For pipeline m-
N pipeline amount, GLminFor the pipeline amount lower limit of pipeline, GB is pipeline set.
Further, step (3) specifically includes:
By the Nonlinear Nonconvex equation of natural gas line flow and pressure
It is converted into the natural gas tide model of following enhancing second order cone constraint type:
In formula, subscript 0 indicates that benchmark Run-time scenario, subscript t indicate t moment,
Subscript max indicates corresponding upper limit value, and subscript min indicates corresponding lower limit value,<>TIndicate the convex envelope function of quadratic term,Indicate the convex envelope function of bilinear terms, κmnIndicate quadratic term convex closure network variable, λmnIndicate bilinear terms convex closure network
Variable.
Further, step (5) specifically includes:
(5.1) the electric-gas interconnection multi-period scheduling model of integrated energy system is solved to obtain optimal solution x0；
(5.2) convex optimization problem is established:
s.t.smn≥0,x∈X
In formula, f0(x) optimization object function of the multi-period scheduling model of integrated energy system is interconnected for electric-gas, x is state
Variable, X are the feasible zone of x, xrFor the state variable optimal solution solved when the r times iteration, smnFor non-negative slack variable, βrTo punish
Weight coefficient is penalized, r is current iteration number,
(5.3) by x0As the initial value of convex optimization problem, DCP iterative solution, progressive updating x are carried outrNumerical value, until natural
Gas constraint violation index GapcLess than preset value, terminate iteration.
Further, natural gas constraint violation index Gap in step (5.3)cCalculation formula are as follows:
In formula: x*For the state variable value after previous iteration,Respectively state variable x*
Middle respective value.
The utility model has the advantages that compared with prior art, the present invention its remarkable advantage is: the present invention ensure that using second order cone optimization
The high efficiency and optimality of solution guarantee that the feasibility understood (can meet the stringent physics of natural gas system using DCP method
Constraint).
Detailed description of the invention
Fig. 1 is the flow diagram of one embodiment of the present of invention；
Fig. 2 is the integrated energy system figure that IEEE-39 node system and Belgian 20 node systems are constituted.
Specific embodiment
The multi-period tide optimization method for present embodiments providing a kind of electric-gas interconnection integrated energy system, such as Fig. 1 institute
Show, includes the following steps:
S1, the power system information and natural gas system information for obtaining electric-gas interconnection integrated energy system respectively.
Wherein, power system information are as follows: power network topology, branch parameters information, generator parameter information, in future time period
Electric load information, the prediction value information of wind-powered electricity generation；The parameter information of natural gas system are as follows: natural gas grid topology, pipe parameter information,
When the line-pack amount of storage of preceding pipeline, the parameter information of gas source, the gas information on load in future time period.
S2, according to the power system information and natural gas system information, construct electric-gas interconnection integrated energy system it is more when
Section scheduling model:
In formula: subscript 0 indicates that benchmark Run-time scenario, subscript t indicate t moment, and i, j, m, n indicate the section in energy resource system
Point；Subscript max indicates upper limit value, and subscript min indicates lower limit value；f0For optimization object function, NGFor generator collection, NgFor combustion
Gas-turbine set, NsFor gas source set, NWFor wind power plant set, T0For abandonment ratio, CG,iFor generator cost coefficient, CS,mFor gas
Source cost coefficient, CW,iFor abandonment cost coefficient,For abandonment percentage,For generator output,For
Generator output lower and upper limit,For when discontinuity surface number,It is expected to contribute for wind-powered electricity generation, PL,i,tFor burden with power,For
Route i-j active power, EN (i) are and node i connected node set, bijFor route i-j susceptance, θ is node phase angle vector,For route i-j active power lower and upper limit；For the amount of natural gas of gas turbine consumption, η is combustion gas wheel
Unit transformation efficiency,For gas source power output, FD,m,tFor natural gas load, GC (m), GP (m), GN (m) are respectively that node m connects
Pressurizing point, gas turbine and the pipeline set connect,For the absorption flow of pressurizing point k,For the flow for flowing through pressurizing point k；AndRespectively pipeline m-n head end, end and average flow rate, CmnFor pipeline m-n pressure drop constant, πmWith πn
Respectively node m, n pressure,Respectively node m low pressure limit and the upper limit；GLmnFor the line- of pipeline m-n
Pack gas-storing capacity, KmnFor the line-pack parameter of pipeline m-n；For gas driven pressurizing point dissipative coefficient,For power generation
The maximum active climbing of machine,Respectively pressurizing point first and end pressure,WithFor pressurizing point step-up ratio
Upper and lower bound,For gas source contribute lower and upper limit,Climb for gas source maximum,For pipeline m-
N pipeline amount, GLminFor the pipeline amount lower limit of pipeline, GB is pipeline set.
In me, formula (1) is multi-period optimization object function, including non-Gas Generator Set cost of electricity-generating, gas supply cost and abandoning
Eolian.It should be noted that gas supply cost contains the cost of electricity-generating of gas turbine indirectly, thus cost of electricity-generating in (1)
Only meter and non-Gas Generator Set.Formula (2)-(7) are Operation of Electric Systems constraint.Formula (2) is node power Constraints of Equilibrium, and formula (3) is retouched
The linear relationship in DC flow model between line power and first and last end node phase angle difference is stated；Formula (4) and formula (5) are respectively
The constraint of generator bound and Climing constant；Formula (6) is power transmission capacity of pow constraint.Formula (8)-(19) are natural gas system dynamic
Operation constraint.Formula (8) is node flow Constraints of Equilibrium, and formula (9) and formula (10) describe pipeline average flow rate and pipeline first and last end
Non-linear relation between node pressure；Formula (11) indicates that the difference of first and last end flow is equal to adjacent two sections Guan Cunbo in pipeline
It is dynamic；Formula (12) indicates that pipeline pipe is deposited and is proportional to first and last end average pressure；Formula (13) describe pressurizing point absorption flow with flow through
Linear relationship between pressurizing point flow；Formula (14) is the constraint of pressurizing point step-up ratio；Formula (15) is the constraint of pressurizing point transmission capacity；
Formula (16) and formula (17) are gas source supplied capacity and Climing constant；Formula (18) is the constraint of node pressure bound；Formula (19) is T0
Moment natural gas system general pipeline leaves limit constraint.
S3, by the electric-gas interconnect the multi-period scheduling model of integrated energy system in natural gas line flow and pressure
Nonlinear Nonconvex it is equations turned for enhancing second order cone constraint type natural gas tide model.
In the multibreak face traffic control model of integrated energy system being made of formula (1)-(19), formula (9) is Nonlinear Nonconvex
Equation, corresponding Non-linear Optimal Model can inevitably encounter the problems such as sensitive, numerical stability is bad to initial value.Formula (9) first
It can relax as formula (20), further, shown in the standard second order tapered such as formula (21) of formula (20).
Formula (9) is relaxed as formula (21) second order tapered, numerical stability problem is can effectively avoid, however is run in optimal solution
Point, formula (21) may not be with formula (9), i.e. second order cone relaxation is not necessarily stringent.Based on this, the present invention proposes a kind of enhancing second order cone
The natural gas tide model of form.
Enhance the natural gas tide model of second order tapered on the basis of formula (21), deeply considers formula (22).Herein for
Formula (22), has two o'clock to need to illustrate: 1) formula (21) and formula (22) combine stringent of equal value with formula (9)；2) it is different from formula (21), formula
(22) it is still non-convex.
Present invention further propose that using the bilinearity in method relaxation formula (22) of convex closure network (Convex envelope)
Item (the non-convex item of essentially nonlinear).Then left side bilinear terms in formula (22)Formula (23) replacement can be used, and for formula
(22) the non-convex item on right side in, definitionThen formula (A-7) right part can adopt
It is replaced with formula (24).Finally, formula (23)-(25) can replace formula (22) using the method for convex closure network.
So far, formula (21) and formula (23)-(25) constitute the natural gas tide model of enhancing second order tapered.
S4, the electric-gas interconnection multi-period scheduling model of integrated energy system after conversion is solved to obtain optimal solution.
S5, using the optimal solution as initial value, and integrated energy system is interconnected to the electric-gas after conversion using DCP method
Multi-period scheduling model carries out linearizing iteration, until natural gas system it is stringent meet trend constraint.
It is noted that enhancing second order Based On The Conic Model is capable of providing more stringent compared to second order cone natural gas tide model
Optimal solution, however optimal solution still may not meet formula (9), that is, relax non-critical establishment, thus present invention further propose that using DCP's
The feasible solution of method recovery natural gas trend.Definition:
Then formula (22) can be stated are as follows:
gmn(x)-hmn(x)≤0 (26)
Based on current optimal solution xr, DCP method linearisation formula (22) recess portion (i.e. hmn(x)) it, then converts (22) to
Following form:
Based on formula (27), DCP solves following convex optimization problem:
In formula: x is state variable, and X is the feasible zone of x；smnFor non-negative slack variable, βrTo punish that weight coefficient, r are repeatedly
Generation number.
In formula (28), slack variable s is introducedmnIt can guarantee the solvability of formula (28).DCP iteratively solves formula (28), gradually more
New xrNumerical value, until Gap in formula (29)cSufficiently small (i.e. primary nonlinear equation (9) is approximate sets up), terminates iteration.
In formula: GapcFor constraint violation index.
S6, it is exported the last solution at the end of iteration as the optimal load flow solution in future time period.
Emulation testing is carried out to the present invention below.
The example that the present invention tests is as shown in Fig. 2, the synthesis energy being made of IEEE-39 node and Belgian 20 node systems
Source system.Table 1 gives second order cone compared with enhancing second order cone model optimization result, by the table it is found that optimizing in the first stage
In, compared to second order Based On The Conic Model, enhance the duality gap Gap of second order Based On The Conic ModeloSmaller (0.43%VS 0.91%), and constrain
Violate index GapcIt is smaller, illustrate that enhancing second order cone model optimization result is more close with primary nonlinear optimum results.Further
, based on the first stage as a result, second order Based On The Conic Model and enhancing second order Based On The Conic Model can restore feasible solution (Gap in second stagecFoot
It is enough small), but enhance second order Based On The Conic Model and primary nonlinear model more close to (its GapoIt is smaller), thus 1 result verification of table
Enhance the validity of strong second order Based On The Conic Model.
1 second order cone of table is compared with enhancing second order cone model optimization result
* Gap hereinoFor the relative error between second order Based On The Conic Model and nonlinear model optimization target values
Above disclosed is only a preferred embodiment of the present invention, and the right model of the present invention cannot be limited with this
It encloses, therefore equivalent changes made in accordance with the claims of the present invention, is still within the scope of the present invention.
Claims (6)
1. a kind of multi-period tide optimization method of electric-gas interconnection integrated energy system, it is characterised in that this method comprises:
(1) power system information and natural gas system information of electric-gas interconnection integrated energy system are obtained respectively；
(2) it according to the power system information and natural gas system information, constructs electric-gas and interconnects the multi-period tune of integrated energy system
Spend model；
(3) by the electric-gas interconnect the multi-period scheduling model of integrated energy system in natural gas line flow and pressure it is non-
The linear non-convex equations turned natural gas tide model for enhancing second order cone constraint type；
(4) the electric-gas interconnection multi-period scheduling model of integrated energy system after conversion is solved to obtain optimal solution；
(5) using the optimal solution as initial value, and it is more to the electric-gas interconnection integrated energy system after conversion using DCP method when
Section scheduling model carries out linearizing iteration, until natural gas system it is stringent meet trend constraint；
(6) it is exported the last solution at the end of iteration as the optimal load flow solution in future time period.
2. the multi-period tide optimization method of electric-gas interconnection integrated energy system according to claim 1, feature exist
In: it is obtained in step (1)
Power system information are as follows: power network topology, branch parameters information, generator parameter information, the electric load letter in future time period
Breath, the prediction value information of wind-powered electricity generation；
The parameter information of natural gas system are as follows: natural gas grid topology, pipe parameter information, when the line-pack of preceding pipeline is stored
It measures, the parameter information of gas source, the gas information on load in future time period.
3. the multi-period tide optimization method of electric-gas interconnection integrated energy system according to claim 1, feature exist
In: the electric-gas established in step (2) interconnects the multi-period scheduling model of integrated energy system specifically:
In formula, subscript 0 indicates that benchmark Run-time scenario, subscript t indicate t moment, and i, j, m, n indicate the node in energy resource system；On
Marking max indicates upper limit value, and subscript min indicates lower limit value；f0For optimization object function, NGFor generator collection, NgFor gas turbine
Set, NsFor gas source set, NWFor wind power plant set, T0For when discontinuity surface number, CG,iFor generator cost coefficient, CS,mFor gas source
Cost coefficient, CW,iFor abandonment cost coefficient,For abandonment percentage,For generator output,For hair
Motor power output lower and upper limit,For abandonment ratio,It is expected to contribute for wind-powered electricity generation, PL,i,tFor burden with power,For route i-
J active power, EN (i) are and node i connected node set, bijFor route i-j susceptance, θ is node phase angle vector,For route i-j active power lower and upper limit；For the amount of natural gas of gas turbine consumption, η is combustion gas wheel
Unit transformation efficiency,For gas source power output, FD,m,tFor natural gas load, GC (m), GP (m), GN (m) are respectively and node m
Pressurizing point, gas turbine and the pipeline set of connection,For the absorption flow of pressurizing point k,For the stream for flowing through pressurizing point k
Amount；AndRespectively pipeline m-n head end, end and average flow rate, CmnFor pipeline m-n pressure drop constant, πm
With πnRespectively node m, n pressure,Respectively node m low pressure limit and the upper limit；GLmnFor the line- of pipeline m-n
Pack gas-storing capacity, KmnFor the line-pack parameter of pipeline m-n；θkFor gas driven pressurizing point dissipative coefficient,For power generation
The maximum active climbing of machine,Respectively pressurizing point first and end pressure,WithFor pressurizing point step-up ratio
Upper and lower bound,For gas source contribute lower and upper limit,Climb for gas source maximum,For pipeline m-
N pipeline amount, GLminFor the pipeline amount lower limit of pipeline, GB is pipeline set.
4. the multi-period tide optimization method of electric-gas interconnection integrated energy system according to claim 3, feature exist
In: step (3) specifically includes:
By the Nonlinear Nonconvex equation of natural gas line flow and pressureConversion
For the natural gas tide model of following enhancing second order cone constraint type:
In formula, subscript 0 indicates that benchmark Run-time scenario, subscript t indicate t moment,
Subscript max indicates corresponding upper limit value, and subscript min indicates corresponding lower limit value,<>TIndicate the convex envelope function of quadratic term,Indicate the convex envelope function of bilinear terms, κmnIndicate quadratic term convex closure network variable, λmnIndicate bilinear terms convex closure network
Variable.
5. the multi-period tide optimization method of electric-gas interconnection integrated energy system according to claim 4, feature exist
In: step (5) specifically includes:
(5.1) the electric-gas interconnection multi-period scheduling model of integrated energy system is solved to obtain optimal solution x0；
(5.2) convex optimization problem is established:
s.t.smn≥0,x∈X
In formula, f0(x) optimization object function of the multi-period scheduling model of integrated energy system is interconnected for electric-gas, x is state variable,
X is the feasible zone of x, xrFor the state variable optimal solution solved when the r times iteration, smnFor non-negative slack variable, βrFor power of punishment
Weight coefficient, r are current iteration number,
(5.3) by x0As the initial value of convex optimization problem, DCP iterative solution, progressive updating x are carried outrNumerical value, until natural gas is about
Beam violates index GapcLess than preset value, terminate iteration.
6. the multi-period tide optimization method of electric-gas interconnection integrated energy system according to claim 5, feature exist
In: natural gas constraint violation index Gap in step (5.3)cCalculation formula are as follows:
In formula: x*For the state variable value after previous iteration,Respectively state variable x*In it is right
It should be worth.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811310889.8A CN109242366B (en) | 2018-11-06 | 2018-11-06 | Multi-period power flow optimization method of electricity-gas interconnection comprehensive energy system |
| PCT/CN2018/114472 WO2020093295A1 (en) | 2018-11-06 | 2018-11-08 | Multi-period power flow optimization method for electricity-gas interconnection integrated energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811310889.8A CN109242366B (en) | 2018-11-06 | 2018-11-06 | Multi-period power flow optimization method of electricity-gas interconnection comprehensive energy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109242366A true CN109242366A (en) | 2019-01-18 |
| CN109242366B CN109242366B (en) | 2020-08-07 |
Family
ID=65076914
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811310889.8A Active CN109242366B (en) | 2018-11-06 | 2018-11-06 | Multi-period power flow optimization method of electricity-gas interconnection comprehensive energy system |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN109242366B (en) |
| WO (1) | WO2020093295A1 (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109687456A (en) * | 2019-01-25 | 2019-04-26 | 燕山大学 | A kind of dispatching method and system of electric power natural gas system |
| CN109949001A (en) * | 2019-02-22 | 2019-06-28 | 上海市建设工程监理咨询有限公司 | A kind of Project Supervision organization management system |
| CN110070213A (en) * | 2019-03-28 | 2019-07-30 | 广东工业大学 | A kind of dispatching method a few days ago of electric-gas integrated energy system |
| CN110210747A (en) * | 2019-05-28 | 2019-09-06 | 河海大学 | A kind of electric heating gas interconnection energy resource system flexibility dispatching method |
| CN110210104A (en) * | 2019-05-28 | 2019-09-06 | 国电南瑞科技股份有限公司 | A kind of multi-energy system traffic control method |
| CN110322051A (en) * | 2019-06-06 | 2019-10-11 | 国网浙江省电力有限公司经济技术研究院 | Consider the integrated energy system Optimal Configuration Method of N-1 security constraint |
| CN110502859A (en) * | 2019-08-28 | 2019-11-26 | 南方电网科学研究院有限责任公司 | Multi tate dynamic emulation method for electrical couplings garden integrated energy system |
| CN110707705A (en) * | 2019-10-22 | 2020-01-17 | 太原理工大学 | Power flow sequence analysis model of electric-gas coupling comprehensive energy system |
| CN110796295A (en) * | 2019-10-15 | 2020-02-14 | 西安交通大学 | Energy internet air network transmission optimization method |
| CN111815111A (en) * | 2020-06-02 | 2020-10-23 | 天津大学 | Regional comprehensive energy expansion planning method considering pipeline risk level |
| CN112713615A (en) * | 2020-12-23 | 2021-04-27 | 山东大学 | Quick coordination scheduling method and system for electricity-gas integrated energy system |
| CN112990606A (en) * | 2021-04-25 | 2021-06-18 | 国网江西省电力有限公司电力科学研究院 | Comprehensive energy system autonomous regulation and control method and device considering regulation and control cost |
| CN113570117A (en) * | 2021-07-02 | 2021-10-29 | 浙江华云电力工程设计咨询有限公司 | Electricity-gas comprehensive energy microgrid optimal scheduling method based on second-order cone relaxation conversion method |
| CN115296345A (en) * | 2022-06-09 | 2022-11-04 | 南方电网科学研究院有限责任公司 | Start-stop-output-standby combined optimization method and device for gas generator set |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111768036B (en) * | 2020-06-29 | 2023-11-03 | 国网上海市电力公司 | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108649580A (en) * | 2018-05-21 | 2018-10-12 | 武汉大学 | A kind of AC/DC mixed power system Security corrective method based on second order cone |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108023364B (en) * | 2017-11-24 | 2019-07-26 | 天津大学 | Power distribution network distributed generation resource maximum access capability calculation method based on convex difference planning |
| CN108667007B (en) * | 2018-04-16 | 2019-12-13 | 清华大学 | Voltage stability margin calculation method considering constraint of electric-gas coupling system |
-
2018
- 2018-11-06 CN CN201811310889.8A patent/CN109242366B/en active Active
- 2018-11-08 WO PCT/CN2018/114472 patent/WO2020093295A1/en active Application Filing
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108649580A (en) * | 2018-05-21 | 2018-10-12 | 武汉大学 | A kind of AC/DC mixed power system Security corrective method based on second order cone |
Non-Patent Citations (2)
| Title |
|---|
| 原凯等: "基于凸差规划的有源配电网电压无功协调控制方法", 《电力系统及其自动化学报》 * |
| 张伊宁等: "计及需求响应与动态气潮流的电-气综合能源系统优化调度", 《电力系统自动化》 * |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109687456A (en) * | 2019-01-25 | 2019-04-26 | 燕山大学 | A kind of dispatching method and system of electric power natural gas system |
| CN109949001A (en) * | 2019-02-22 | 2019-06-28 | 上海市建设工程监理咨询有限公司 | A kind of Project Supervision organization management system |
| CN109949001B (en) * | 2019-02-22 | 2021-05-07 | 上海市建设工程监理咨询有限公司 | Construction engineering supervision organization management system |
| CN110070213A (en) * | 2019-03-28 | 2019-07-30 | 广东工业大学 | A kind of dispatching method a few days ago of electric-gas integrated energy system |
| WO2020237700A1 (en) * | 2019-05-28 | 2020-12-03 | 国电南瑞科技股份有限公司 | Operation scheduling method for multiple energy systems |
| CN110210747A (en) * | 2019-05-28 | 2019-09-06 | 河海大学 | A kind of electric heating gas interconnection energy resource system flexibility dispatching method |
| CN110210104A (en) * | 2019-05-28 | 2019-09-06 | 国电南瑞科技股份有限公司 | A kind of multi-energy system traffic control method |
| CN110210747B (en) * | 2019-05-28 | 2022-07-29 | 河海大学 | Flexible scheduling method for electric-heating-gas interconnection energy system |
| CN110322051A (en) * | 2019-06-06 | 2019-10-11 | 国网浙江省电力有限公司经济技术研究院 | Consider the integrated energy system Optimal Configuration Method of N-1 security constraint |
| CN110502859A (en) * | 2019-08-28 | 2019-11-26 | 南方电网科学研究院有限责任公司 | Multi tate dynamic emulation method for electrical couplings garden integrated energy system |
| CN110796295A (en) * | 2019-10-15 | 2020-02-14 | 西安交通大学 | Energy internet air network transmission optimization method |
| CN110796295B (en) * | 2019-10-15 | 2022-10-25 | 西安交通大学 | Energy Internet air network transmission optimization method |
| CN110707705A (en) * | 2019-10-22 | 2020-01-17 | 太原理工大学 | Power flow sequence analysis model of electric-gas coupling comprehensive energy system |
| CN111815111A (en) * | 2020-06-02 | 2020-10-23 | 天津大学 | Regional comprehensive energy expansion planning method considering pipeline risk level |
| CN111815111B (en) * | 2020-06-02 | 2022-05-13 | 天津大学 | Regional comprehensive energy expansion planning method considering pipeline risk level |
| CN112713615A (en) * | 2020-12-23 | 2021-04-27 | 山东大学 | Quick coordination scheduling method and system for electricity-gas integrated energy system |
| CN112990606A (en) * | 2021-04-25 | 2021-06-18 | 国网江西省电力有限公司电力科学研究院 | Comprehensive energy system autonomous regulation and control method and device considering regulation and control cost |
| CN113570117A (en) * | 2021-07-02 | 2021-10-29 | 浙江华云电力工程设计咨询有限公司 | Electricity-gas comprehensive energy microgrid optimal scheduling method based on second-order cone relaxation conversion method |
| CN113570117B (en) * | 2021-07-02 | 2024-02-09 | 浙江华云电力工程设计咨询有限公司 | Electric-gas comprehensive energy microgrid optimal scheduling method based on second order cone relaxation conversion method |
| CN115296345A (en) * | 2022-06-09 | 2022-11-04 | 南方电网科学研究院有限责任公司 | Start-stop-output-standby combined optimization method and device for gas generator set |
| CN115296345B (en) * | 2022-06-09 | 2023-08-22 | 南方电网科学研究院有限责任公司 | Start-stop-output-standby combined optimization method and device for gas generator set |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020093295A1 (en) | 2020-05-14 |
| CN109242366B (en) | 2020-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109242366A (en) | A kind of multi-period tide optimization method of electric-gas interconnection integrated energy system | |
| CN103490410B (en) | Micro-grid planning and capacity allocation method based on multi-objective optimization | |
| Martinez-Mares et al. | A robust optimization approach for the interdependency analysis of integrated energy systems considering wind power uncertainty | |
| CN108599154B (en) | Three-phase unbalanced distribution network robust dynamic reconstruction method considering uncertainty budget | |
| CN111652441B (en) | Distribution network optimization method of gas-electricity integrated energy system considering gas-electricity combined demand response | |
| CN109508499A (en) | Multi-period more optimal on-positions of scene distribution formula power supply and capacity research method | |
| CN109978362A (en) | A kind of modeling of gas power grid joint multizone integrated energy system and systems organization method | |
| CN106169108A (en) | Active distribution network short-term active optimization method containing battery energy storage system | |
| CN112701687B (en) | Robust optimization operation method of gas-electricity distribution network system considering price type combined demand response | |
| CN111030110B (en) | Robust cooperative scheduling method for electric power-natural gas coupling system considering electric power conversion gas consumption wind power | |
| CN109818347B (en) | Assessment method for wind power consumption capability of electric power system | |
| CN106786799A (en) | A kind of DC link power step elelctrochemical power generation plan optimization method | |
| CN113705962B (en) | Virtual power plant day-ahead scheduling method based on distribution robust optimization | |
| CN107947245A (en) | Consider the equivalent optimal load flow model building method of natural gas system constraint | |
| CN108493998A (en) | Consider the robust Transmission Expansion Planning in Electric method of demand response and N-1 forecast failures | |
| CN109767127A (en) | Electric-gas association system reliability judgment method based on electrical combined optimization trend | |
| CN108594658A (en) | A kind of electric-gas coupled system maximum probability load margin Model for Multi-Objective Optimization and its method for solving | |
| CN115330021A (en) | Comprehensive energy operation optimization system and method considering methane electric heat utilization ratio | |
| CN107769266A (en) | A kind of Multiple Time Scales generate electricity and standby combined optimization method | |
| CN109066695A (en) | A kind of electrical optimal energy flux computation method of two stages linearisation | |
| CN111768036B (en) | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid | |
| Liang et al. | Capacity configuration optimization of wind-solar combined power generation system based on improved grasshopper algorithm | |
| CN110380447A (en) | A kind of lower electric-gas interconnection energy resource system drop Risk Scheduling method of blower failure | |
| CN111342453B (en) | Electrical comprehensive energy system standby decision method considering various types of standby resources | |
| CN114465226A (en) | Method for establishing multi-level standby acquisition joint optimization model of power 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 |